Review Special Issues

Chemical structure, properties and potential applications of surfactin, as well as advanced strategies for improving its microbial production

  • #Authors contributed equally to this work
  • Received: 01 December 2022 Revised: 08 February 2023 Accepted: 08 March 2023 Published: 14 March 2023
  • Surfactin, a cyclic lipopeptide produced by microbes belonging to the genus Bacillus, is one of the most effective biosurfactants available in many industrial fields. However, its low production and high cost have intensively constrained its commercial applications. In this review, we first summarize the molecular structure, biological properties, beneficial roles and potential applications of surfactin in the fields of medical care and food safety, highlighting the great medical and commercial values of making its industrial production into reality. Further, genetic regulation for surfactin biosynthesis and advanced strategies for enhancing its microbial production, including optimizing fermentation conditions, rational genetic engineering and synthetic biology combined with metabolic engineering approaches, are elucidated. Finally, prospects for improving surfactin biosynthesis are discussed, and the establishment of suitable chassis hosts for exogenous production of surfactin might serve as an important strategy in future research.

    Citation: Cheng Zhen, Xian-Feng Ge, Yi-Ting Lu, Wen-Zheng Liu. Chemical structure, properties and potential applications of surfactin, as well as advanced strategies for improving its microbial production[J]. AIMS Microbiology, 2023, 9(2): 195-217. doi: 10.3934/microbiol.2023012

    Related Papers:

    [1] San-Xing Wu, Xin-You Meng . Dynamics of a delayed predator-prey system with fear effect, herd behavior and disease in the susceptible prey. AIMS Mathematics, 2021, 6(4): 3654-3685. doi: 10.3934/math.2021218
    [2] Binfeng Xie, Na Zhang . Influence of fear effect on a Holling type III prey-predator system with the prey refuge. AIMS Mathematics, 2022, 7(2): 1811-1830. doi: 10.3934/math.2022104
    [3] Fatao Wang, Ruizhi Yang, Yining Xie, Jing Zhao . Hopf bifurcation in a delayed reaction diffusion predator-prey model with weak Allee effect on prey and fear effect on predator. AIMS Mathematics, 2023, 8(8): 17719-17743. doi: 10.3934/math.2023905
    [4] Jie Liu, Qinglong Wang, Xuyang Cao, Ting Yu . Bifurcation and optimal harvesting analysis of a discrete-time predator–prey model with fear and prey refuge effects. AIMS Mathematics, 2024, 9(10): 26283-26306. doi: 10.3934/math.20241281
    [5] A. Q. Khan, Ibraheem M. Alsulami . Complicate dynamical analysis of a discrete predator-prey model with a prey refuge. AIMS Mathematics, 2023, 8(7): 15035-15057. doi: 10.3934/math.2023768
    [6] Xiaoming Su, Jiahui Wang, Adiya Bao . Stability analysis and chaos control in a discrete predator-prey system with Allee effect, fear effect, and refuge. AIMS Mathematics, 2024, 9(5): 13462-13491. doi: 10.3934/math.2024656
    [7] Kottakkaran Sooppy Nisar, G Ranjith Kumar, K Ramesh . The study on the complex nature of a predator-prey model with fractional-order derivatives incorporating refuge and nonlinear prey harvesting. AIMS Mathematics, 2024, 9(5): 13492-13507. doi: 10.3934/math.2024657
    [8] Weili Kong, Yuanfu Shao . Bifurcations of a Leslie-Gower predator-prey model with fear, strong Allee effect and hunting cooperation. AIMS Mathematics, 2024, 9(11): 31607-31635. doi: 10.3934/math.20241520
    [9] Yaping Wang, Yuanfu Shao, Chuanfu Chai . Dynamics of a predator-prey model with fear effects and gestation delays. AIMS Mathematics, 2023, 8(3): 7535-7559. doi: 10.3934/math.2023378
    [10] Jing Zhang, Shengmao Fu . Hopf bifurcation and Turing pattern of a diffusive Rosenzweig-MacArthur model with fear factor. AIMS Mathematics, 2024, 9(11): 32514-32551. doi: 10.3934/math.20241558
  • Surfactin, a cyclic lipopeptide produced by microbes belonging to the genus Bacillus, is one of the most effective biosurfactants available in many industrial fields. However, its low production and high cost have intensively constrained its commercial applications. In this review, we first summarize the molecular structure, biological properties, beneficial roles and potential applications of surfactin in the fields of medical care and food safety, highlighting the great medical and commercial values of making its industrial production into reality. Further, genetic regulation for surfactin biosynthesis and advanced strategies for enhancing its microbial production, including optimizing fermentation conditions, rational genetic engineering and synthetic biology combined with metabolic engineering approaches, are elucidated. Finally, prospects for improving surfactin biosynthesis are discussed, and the establishment of suitable chassis hosts for exogenous production of surfactin might serve as an important strategy in future research.



    In population ecology, understanding how predators and primary producers influence nutrient flow relative to each other is important. Ecosystem interactions and predator-prey relationships are governed by predation and the delivery of resource processes. The identification of ecological factors that can alter or control dynamic behavior requires theoretical and experimental research. One way to study these questions is by means of experimental control, and another useful way is via mathematical modeling as well as computer simulations. Over decades of theoretical ecology and biomathematics development, mathematical modeling has become an indispensable tool for scientists in related fields to study ecosystems. Since Lotka [1] and Volterra [2], as cornerstones of theoretical ecology, published the first study of predator-prey dynamics, any species in nature can be a predator or prey, and due to its prevalence, it has become one of the most popular topics for researchers to study [3,4,5]. Besides, because biological resources are renewable and have the most unique development mechanisms, the over-utilization of biological resources and the destruction of the environment by humans will directly affect the balance of the ecosystem. Maintaining ecological balance and meeting humans material needs have attracted the most attention from researchers focused on the scientific management of renewable resource development [6,7,8].

    Shelter serves as a defense strategy. It refers broadly to a series of behaviors by prey to avoid predators in order to increase their survival rate. The concept of sanctuary was first developed by Maynard-Smith [9] and Gause et al. [10], and its popularity has been very high, garnering widespread attention from many scholars [11,12,13,14,15]. Sih et al.[16] investigated the effects of prey refuge in a three-species model and concluded that the system's stability is strongly related to the refuge. Also, similar findings can be displayed in [17,18,19,20,21,22]. The two modes of refuge analyzed by Gonzalez-Olivares et al. [17] have diverse stability domains in terms of the parameter space. Qi et al.[21] ensure the stability of the system by varying the strength of the refuge.

    Through reviewing a large amount of literature, we begin to consider [23,24] as a basis for the two prey and one predator species that will be modeled in this article. We assume that at a certain time t, the populations of the two prey and one predator are x1(t), x2(t), and y(t), respectively. Based on the above, we construct the following model:

    {dx1dt=r1x1(1x1K1)a1x1x2c1(1m1)x1yq1E1x1,dx2dt=r2x2(1x2K2)a2x1x2c2(1m2)x2yq2E2x2,dydt=e1(1m1)x1y+e2(1m2)x2ydyq3E3y. (1.1)

    The significance of the full parameters is annotated in Table 1.

    Table 1.  Biological meaning of parameters.
    Parameters Biological meaning
    r1,r2 Growth rates of prey x1 and prey x2
    K1,K2 Carrying capacity of prey x1 and prey x2
    a1,a2 Interspecific competition between prey x1 and prey x2
    c1,c2 Predation coefficients for prey x1 and prey x2
    m1,m2 Refuge rates of prey x1 and prey x2
    e1,e2 Conversion factors for prey x1 and prey x2
    q1,q2,q3 Captureability factors for prey x1, prey x2 and predator y
    E1,E2,E3 Harvesting efforts for prey x1, prey x2 and predator y
    d Predator y mortality rate

     | Show Table
    DownLoad: CSV

    Most species in nature, including humans, are influenced by fear. Fear may cause an abnormal state and behavior to arise. As usual, prey have an innate fear of predators. The ecology of fear is related to combining the optimal behavior of prey and predators with their population densities [25,26]. In view of reality, it is a fact that prey fear predators, which is seen as a psychological effect that can have a lasting impact on prey populations. This psychological influence is often easy to overlook, but it is necessary to consider it in the context of practical ecology [27]. Wang et al. [28] first considered the effect of the fear factor on the model and first proposed the fear of prey F(k,y). Afterwards, some researchers have investigated the effects of the fear effect and predator interferences in some three-dimensional systems as well as explored the generation of Hopf bifurcation conditions in the presence of a fear parameter as a bifurcation parameter [29,30,31,32]. Zanette et al. [33] observed that prey will reduce reproducing because of fear of being killed by predators, thus decreasing the risk of being killed after giving birth, which also leads directly to a decline in prey birth rates. According to the above discussion, our paper considers the different fears ki caused by predators for the two prey species.

    In reality, when prey feel the crisis of being hunted, they will reproduce less and increase their survival rate. These conditions about the fear factor F(ki,y) (i=1,2) are listed as follows:

    1) F(0,y)=1: prey production does not decrease when the prey does not fear the predator;

    2) F(ki,0)=1: even though the prey will develop a fear of predators and there will be no predators, prey production will still not decline;

    3) limkiF(ki,y)=0: when the prey's fear of the predator is very high, this will result in the prey production tending to zero;

    4) limyF(ki,y)=0: prey have a fear of predators, and when predator numbers are too large, this can also lead to prey production tending to zero;

    5) F(ki,y)ki<0: the greater the prey's fear of predators, the less productive it will be;

    6) F(ki,y)y<0: predators are inversely proportional to their prey.

    For ease of analysis, we draw on Wang et al. [28] to consider the fear effect:

    F(ki,y)=11+kiy(i=1,2), (1.2)

    obviously, F(ki,y) (i=1,2) in (1.2) satisfies conditions 1)6). Based on the above conditions, this study will consider the effect of fear on system (1.1) to obtain system (1.3).

    {dx1dt=r1x11+k1y(1x1K1)a1x1x2c1(1m1)x1yq1E1x1,dx2dt=r2x21+k2y(1x2K2)a2x1x2c2(1m2)x2yq2E2x2,dydt=e1(1m1)x1y+e2(1m2)x2ydyq3E3y. (1.3)

    Notably, most biological parameters in much of the literature are fixed constants. However, in reality, the survival of species is full of unknowns, and all data are not always constant, which can lead to deviations from the ideal model with fixed parameters. In order to make the model more relevant and the results more accurate, we cannot just consider fixed parameters. Therefore, to make the study more convincing, it is necessary to target imprecise parameters. Professor Zadeh [34], who first proposed the fuzzy set theory, also argued that the application of fuzzy differential equations is a more accurate method for modeling biological dynamics in the absence of accurate data conditions [35]. Moreover, the first introduction of the idea of fuzzy derivatives came from Chang and Zadeh [36]. Further, Kaleva [37] studied the generalized fuzzy derivatives based on Hukuhara differentiability, the Zadeh extension principle, and the strong generalized differentiability concept. Bede et al. [38] employed the notion of strongly generalized differentiability to investigate fuzzy differential equations. Khastan and Nieto [39] solved the margin problem for fuzzy differential equations in their article. Motivated by the method of Pal [13] and Wang [23], we assume that the imprecise parameters ~r1, ~r2, ~a1, ~a2, ~c1, ~c2, ~e1, ~e1 and ˜d represent all triangular fuzzy numbers (the relevant theories of fuzzy sets are detailed in Appendix A), then the system (1.3) can be written as

    {~dx1dt=~r1x11+k1y(1x1K1)~a1x1x2~c1(1m1)x1yq1E1x1,~dx2dt=~r2x21+k2y(1x2K2)~a2x1x2~c2(1m2)x2yq2E2x2,~dydt=~e1(1m1)x1y+~e2(1m2)x2y˜dyq3E3y, (1.4)

    we cut these imprecise parameters ~r1, ~r2, ~a1, ~a2, ~c1, ~c2, ~e1, ~e1, and ˜d by using α-level. System (1.4) can be expressed as follows:

    {(dx1dt)αL=rα1Lx11+k1yrα1R1+k1yx21K1aα1Rx1x2cα1R(1m1)x1yq1E1x1,(dx1dt)αR=rα1Rx11+k1yrα1L1+k1yx21K1aα1Lx1x2cα1L(1m1)x1yq1E1x1,(dx2dt)αL=rα2Lx21+k2yrα2R1+k2yx22K2aα2Rx1x2cα2R(1m2)x2yq2E2x2,(dx2dt)αR=rα2Rx21+k2yrα2L1+k2yx22K2aα2Lx1x2cα2L(1m2)x2yq2E2x2,(dydt)αL=eα1L(1m1)x1y+eα2L(1m2)x2ydαRyq3E3y,(dydt)αR=eα1R(1m1)x1y+eα2R(1m2)x2ydαLyq3E3y. (1.5)

    Introducing weighted sum, we change (1.5) to (1.6)

    {dx1dt=w1(dx1dt)αL+w2(dx1dt)αR,dx2dt=w1(dx2dt)αL+w2(dx2dt)αR,dydt=w1(dydt)αL+w2(dydt)αR, (1.6)

    where w1 and w2 are satisfied with w1+w2=1, and w1,w20. Simplifying the system (1.6), we obtain

    {dx1dt=A11+k1yx1A21+k1yx21K1A3x1x2A4(1m1)x1yq1E1x1,dx2dt=B11+k2yx2B21+k2yx22K2B3x1x2B4(1m2)x2yq2E2x2,dydt=C1(1m1)x1y+C2(1m2)x2yC3yq3E3y, (1.7)

    where

    A1=w1rα1L+w2rα1R,  A2=w1rα1R+w2rα1L,  A3=w1aα1R+w2aα1L,A4=w1cα1R+w2cα1L,  B1=w1rα2L+w2rα2R,  B2=w1rα2R+w2rα2L,B3=w1aα2R+w2aα2L,  B4=w1cα2R+w2cα2L,  C1=w1eα1L+w2eα1R,C2=w1eα2L+w2eα2R,  C3=w1dαR+w2dαL.

    The rest of the paper is shown below: In Section 2, we first prove the nonnegativity and boundedness of the system (1.7). Sections 3 and 4 discuss all possible equilibria and give conditions for the local asymptotic stability and global asymptotic stability of the equilibria. Immediately after that, in Section 5, we analyze the Hopf bifurcation by using the normal form theory. In Section 6, we numerically simulate the theoretical results of Sections 4 and 5. Finally, the article ends with detailed conclusions.

    In this section, we give the following theorem to ensure the boundedness and nonnegativity of the solutions of the system (1.7).

    Theorem 2.1. Provided that the initial values x1(0)>0, x2(0)>0, and y(0)>0, all solutions of system (1.7) are nonnegative.

    Proof. It is not difficult to find that the right half of the system (1.7) fulfills the local Lipschitzian condition. Integrating both sides of the system (1.7) at the same time yields

    x1(t)=x1(0)[expt0(A11+k1yA21+k1yx1K1A3x2A4(1m1)yq1E1)ds]>0,x2(t)=x2(0)[expt0(B11+k2yB21+k2yx2K2B3x1B4(1m2)yq2E2)ds]>0,y(t)=y(0)[expt0(C1(1m1)x1C2(1m2)x2C3q3E3)ds]>0. (2.1)

    If the solution curve starts at any internal point of R3+={(x1(t),x2(t),y(t))R3:x1(t)0, x2(t)0, y(t)0}, then x1(t), x2(t), and y(t) will always be nonnegative.

    Theorem 2.2. Assume that the initial values x1(0), x2(0), and y(0) are all greater than zero. The feasible region Ω is a positive invariant set of the system (1.7) defined by

    Ω={(x1(t),x2(t),y(t))R3+:C1A4x1(t)+C2B4x2(t)+y(t)ϕμ},
    (i.e  Ω={(x1(t),x2(t),y(t))R3+:w1eα1L+w2eα1Rw1cα1R+w2cα1Lx1(t)+w1eα2L+w2eα2Rw1cα2R+w2cα2Lx2(t)+y(t)ϕμ}),

    where μ=min{q1E1,q2E2,C3+q3E3}.

    Proof. Define a function

    W(t)=C1A4x1(t)+C2B4x2(t)+y(t). (2.2)

    After taking the derivative on both sides of (2.2), we obtain

    dWdt=C1A4dx1dt+C2B4dx2dt+dydt. (2.3)

    Furthermore, we can obtain

    dWdt+μW=C1A4(1+k1y)(A1x1A2x21K1)C1A3A4x1x2C1(1m1)x1ydWdt+μW=+C2B4(1+k2y)(B1x2B2x22K2)C2B3B4x1x2C2(1m2)x2ydWdt+μW=+C1(1m1)x1y+C2(1m2)x2yC3yq3E3y+C1μA4x1dWdt+μW=C1A4q1E1x1C2B4q2E2x2+C2μB4x2+μy,dWdt+μW=C1A4(1+k1y)(A1x1A2x21K1)+C2B4(1+k2y)(B1x2B2x22K2)dWdt+μW=+C1A4x1(μq1E1)+C2B4x2(μq2E2)+y(μC3q3E3)dWdt+μW=(C1A3A4+C2B3B4)x1x2, (2.4)

    where μ=min{q1E1,q2E2,C3+q3E3}. Let ϕ1=A21K14A2,ϕ2=B21K24B2, ϕC1A4ϕ1+C2B4ϕ2, we have

    dWdt+μWC1A4ϕ1+C2B4ϕ2=ϕ. (2.5)

    Therefore, it can be deduced that

    Wϕμ+Neμt, (2.6)

    where N is a positive constant. Then we can further obtain

    lim suptWϕμ, (2.7)

    which indicates that the feasible domain Ω is a positive invariant set.

    In this section, we discuss the existence of all equilibria in the system (1.7). All equilibria for system (1.7) are provided by

    (1) Trivial equilibrium P1=(0,0,0).

    (2) Axial equilibrium P2=(xς1,0,0) exists if A1>q1E1, where xς1=K1(A1q1E1)A2.

    (3) Axial equilibrium P3=(0,xΥ2,0) exists if B1>q2E2, where xΥ2=K2(B1q2E2)B2.

    (4) Axial equilibrium P4=(xΨ1,0,yΨ) exists if A1A2xΨ1K1+q1E1 and Δ1>A4(1m1)+k1q1E1, where

    xΨ1=C3+q3E3C1(1m1),yΨ=Δ1(A4(1m1)+k1q1E1)2k1A4(1m1),Δ1=4(k1A4(1m1))(A1A2xΨ1K1q1E1)+(A4m1A4k1q1E1)2. (3.1)

    (5) Axial equilibrium P5=(0,xν2,yν) exists if B1B2xν2K2+q2E2 and Δ2>B4(1m2)+k2q2E2, where

    xν2=C3+q3E3C2(1m2),yν=Δ2(B4(1m2)+k2q2E2)2k2B4(1m2),Δ2=4(k2B4(1m2))(B1B2xν2K2q2E2)+(B4m2B4k2q2E2)2. (3.2)

    (6) Axial equilibrium P6=(xι1,xι2,0) exists if B1>B3xι1+q2E2, K1K2A3B3>A2B2 and A3B1K2+B2q1E1>A3K2q2E2+A1B2, where

    xι1=A3K1K2(B1q2E2)+K1B2(q1E1A1)K1K2A3B3A2B2,xι2=K2(B1B3xι1q2E2)B2. (3.3)

    (7) Internal equilibrium P7=(x1,x2,y) exists, and its value will be given in the proof of Theorem 3.1.

    Theorem 3.1. When g3>0 and g4g5<0 are met, there is an internal equilibrium P7.

    Proof. We derive that from the second equation of the system (1.7)

    g1y2+g2y+g3=0, (3.4)

    where

    g1=k2B4(1m2),  g2=k2B3x1B4(1m2)q2E2k2,  g3=(B1x2K2B2B3x1q2E2).

    It follows from the Descartes law of signs that Eq (3.4) has one and only one solution y greater than zero if and only if g3>0, i.e., B1>x2K2B2+B3x1+q2E2. Substituting y into the algebra expression on the right side of the first equation of the system (1.7) equals zero; furthermore, we obtain

    x2=A1A3+A3k1yA2x1K1A3+K1A3k1yA4(1m1)yA3q1E1A3. (3.5)

    Introduce (3.5) into the right side of the third equation of the system (1.7), which meets zero, it simplifies to obtain

    g4x1g5=0, (3.6)

    where

    g4=[C1(1m1)A2C2(1m2)K1A3+K1A3k1y],g5=[A1C2(1m2)A3+A3k1yA4C2(1m1)(1m2)yA3C2q3E3(1m2)A3C3q3E3].

    Reusing the Descartes law of signs, we can assert that there exists at least one positive solution x1 of Eq (3.6) if and only if g4g5<0. And then we can deduce that

    x2=A1A3+A3k1yA2x1K1A3+K1A3k1yA4(1m1)yA3q1E1A3.

    then the interior equilibrium P7(x1,x2,y) exists.

    In this section, the Jocabian matrix will be used to prove the local stability of all equilibria. Moreover, we prove the global stability of the internal equilibrium P7 by constructing a Lyapunov function.

    The Jocabian matrix for system (1.7) is given below:

    M=(M11M12M13M21M22M23M31M32M33), (4.1)

    where

    M11=A11+k1y2A21+k1yx1K1A3x2A4(1m1)yq1E1,M12=A3x1,M13=k1A1x1(1+k1y)2+k1A2(1+k1y)2x21K1A4(1m1)x1,M21=B3x2,M22=B11+k2y2B21+k2yx2K2B3x1B4(1m2)yq2E2,M23=k2B1x2(1+k2y)2+k2B2(1+k2y)2x22K2B4(1m2)x2,M31=C1(1m1)y,M32=C2(1m2)y,  M33=C1(1m1)x1+C2(1m2)x2C3q3E3. (4.2)

    Through simple calculation, we directly draw the conclusion that trivial and axial equilibria are locally asymptotically stable:

    (1) P1(0,0,0) is locally asymptotically stable if

    A1q1E1<0andB1q2E2<0. (4.3)

    (2) P2(xς1,0,0) is locally asymptotically stable if

    B1q2E2<B3xς1q2<B3(C3+q3E3)C1(1m1)q2. (4.4)

    (3) P3(0,xΥ2,0) is locally asymptotically stable if

    A1q1E1<A3xΥ2q1<A3(C3+q3E3)C2(1m2)q1. (4.5)

    (4) P4(xΨ1,0,yΨ) is locally asymptotically stable if

    B1q2(1+k2yΨ)E2<B3xΨ1+B4(1m2)yΨq2. (4.6)

    (5) P5(0,xν2,yν) is locally asymptotically stable if

    A1q1(1+k1yν)E1<A3xν2+A4(1m1)yνq1. (4.7)

    (6) P6(xι1,xι2,0) is locally asymptotically stable if

    C1(1m1)xι1+C2(1m2)xι2<C3+q3E3. (4.8)

    We draw the conclusion that the internal equilibrium P7(x1,x2,y) is locally asymptotically stable from the proof of Theorem 4.1.

    Theorem 4.1. The internal equilibrium P7 is locally asymptotically stable if it exists and the following conditions are fulfilled:

    ψ1>0,ψ1ψ2>0,ψ3>0. (4.9)

    Proof. The Jocabian matrix of system (1.7) at (x1,x2,y) is

    (L1L2L3L4L5L6L7L8L9), (4.10)

    where

    L1=A2(1+k1y)x1K1<0,L2=A3x1<0,L3=k1A1x1(1+k1y)2A4(1m1)x1+k1A2(1+k1y)2(x1)2K1,L4=B3x2<0,L5=B2(1+k2y)x2K2<0,L6=k2B1x2(1+k2y)2B4(1m2)x2+k2B2(1+k2y)2(x2)2K2,L7=C1(1m1)y>0,L8=C2(1m2)y>0,L9=0. (4.11)

    Therefore, the characteristic equation at P7 can be expressed as

    η3+ψ1η2+ψ2η+ψ3=0, (4.12)

    where

    ψ1=L1L5,ψ2=L1L5L6L8L3L7L2L4,ψ3=L8(L1L6L3L4)+L7(L3L5+L2L6). (4.13)

    The Routh-Hurwitz criterion shows that the internal equilibrium P7 is locally asymptotically stable; the following conditions need to be met: ψ1>0, ψ1ψ2>0, and ψ3>0.

    This subsection studies the global asymptotic stability of interior equilibrium P7.

    Theorem 4.2. If condition 4Γ1Γ2l1l2A2B2(1+k1y)(1+k2y) > (l1A3+l2B3)2 (i.e. 4Γ1Γ2l1l2(w1rα1R+w2rα1L)(w1rα2R+w2rα2L)(1+k1y)(1+k2y) > (l1(w1aα1R+w2aα1L)+l2(w1aα2R+w2aα2L))2) holds, then P7 is globally asymptotically stable.

    Proof. We construct a Lyapunov function:

    V(x1,x2,y)=l1[x1x1x1ln(x1x1)]+l2[x2x2x2ln(x2x2)]  +yyyln(yy). (4.14)

    Obviously, xixixiln(xixi)0 (i=1,2) and yyyln(yy)0, thus V0. Taking the derivative of V(x1,x2,y) over t, one has

    dVdt=l1(x1x1x1)dx1dt+l2(x2x2x2)dx2dt+yyydydt, (4.15)

    where

    x1x1x1dx1dt=k1A1(1+k1y)(1+k1y)(x1x1)(yy)A2(x1x1)2K1(1+k1y)(1+k1y)x1x1x1dx1dt=A2k1(x1yx1y)K1(1+k1y)(1+k1y)(x1x1)A3(x1x1)(x2x2)x1x1x1dx1dt=A4(1m1)(x1x1)(yy),x2x2x2dx2dt=k2B1(1+k2y)(1+k2y)(x2x2)(yy)B2(x2x2)2K2(1+k2y)(1+k2y)x2x2x2dx2dt=B2k2(x2yx2y)K2(1+k2y)(1+k2y)(x2x2)B3(x1x1)(x2x2)x2x2x2dx2dt=B4(1m2)(x2x2)(yy),yyydydt=C1(1m1)(x1x1)(yy)+C2(1m2)(x2x2)(yy). (4.16)

    To simplify the calculation, let

    x1yx1y=y(x1x1)x1(yy),  x2yx2y=y(x2x2)x2(yy),Γ1=1K1(1+k1y)(1+k1y),  Γ2=1K2(1+k2y)(1+k2y),l1=C1(1m1)Γ1k1(K1A1+A2x1)+A4(1m1),  l2=C2(1m2)Γ2k2(K2B1+B2x2)+A4(1m1). (4.17)

    We obtain

    dVdt={Γ1l1A2(1+k1y)(x1x1)2+(l1A3+l2B3)(x1x1)(x2x2)+Γ2l2B2(1+k2y)(x2x2)2}dVdt  =YGY, (4.18)

    where

    Y=[(x1x1),(x2x2)],  G=(Γ1l1A2(1+k1y)l1A3+l2B32l1A3+l2B32Γ2l2B2(1+k2y)).

    Therefore, dVdt<0 if and only if 4Γ1Γ2l1l2A2B2(1+k1y)(1+k2y)>(l1A3+l2B3)2.

    In this section, we will use the normal form theory introduced by Hassard et al.[40] and the central manifold theory [41] to study the Hopf bifurcation of the system (1.7). When the system (1.7) undergoes Hopf bifurcation, the corresponding characteristic equation must have a pair of conjugate pure imaginary roots, that is,

    η1,2=±iω,i=1. (5.1)

    Consider the parameter k1 as a bifurcation parameter. When the value of parameter k1 changes near the critical point kΞ1 of Hopf bifurcation, the pure imaginary roots ±iω will become a complex eigenvalue η=ρ+i˜ω. Substituting η=ρ+i˜ω into Eq (4.12), we need to separate the imaginary and real parts to get

    ρ3+ψ3+ρψ2+ρ2ψ13ρ˜ω2ψ1˜ω2=0, (5.2)
    3ρ2˜ω+ψ2˜ω+2ρψ1˜ω˜ω3=0. (5.3)

    By simplifying Eqs (5.2) and (5.3), we obtain

    ψ38ρ32ρψ28ρ2ψ1ψ1ψ22ρψ21=0, (5.4)

    at k1=kΞ1, taking the derivative of Eq (5.4) over k1 yields

    dρdk1|k1=kΞ1=12(dψ3dk1ψ1dψ2dk1ψ2dψ1dk1)/(ψ2+ψ21). (5.5)

    If it satisfies dρdk1|k1=kΞ10, the system (1.7) will generate Hopf bifurcation, which indicates that when parameter k1 crosses the bifurcation critical point kΞ1, the population state evolves from stable equilibrium to periodic oscillation over time.

    When the system (1.7) undergoes Hopf bifurcation at k1=kΞ1, the final decision condition is also met. Considering that the characteristic roots of Eq (4.12) are η1,2=±iω and η3=ψ1, in order to obtain this condition, we introduce

    z1=x1x1,z2=x2x2,z3=yy. (5.6)

    Substituting (5.6) into the system (1.7) and separating the linear and nonlinear parts, it can be obtained that

    (˙z1˙z2˙z3)=J(P7)(z1z2z3)+(F1(z1,z2,z3)F2(z1,z2,z3)F3(z1,z2,z3)), (5.7)

    where

    F1(z1,z2,z3)=2j1+j2+j33tj1j2j3zj11zj22zj33+O((|z1|+|z2|+|z3|)4),F2(z1,z2,z3)=2j1+j2+j33nj1j2j3zj11zj22zj33+O((|z1|+|z2|+|z3|)4),F3(z1,z2,z3)=2j1+j2+j33lj1j2j3zj11zj22zj33+O((|z1|+|z2|+|z3|)4), (5.8)

    where O((|z1|+|z2|+|z3|)4) is a fourth-order polynomial function about variables (|z1|,|z2|,|z3|), while tj1j2j3, nj1j2j3, and lj1j2j3 can be obtained through calculation:

    t101=A1k12(1+k1y)2+A2k1(1+k1y)2x1K1A4(1m1)2,  t002=A1x1k21(1+k1y)3A2k21(1+k1y)3(x1)2K1,t102=A1k213(1+k1y)3A2k213(1+k1y)32x1K1,  t003=A1x1k31(1+k1y)4+A2k31(1+k1y)4(x1)2K1,t110=A32,  t200=A21+k1y1K1,  t201=A2k13(1+k1y)21K1,t011=t020=t030=t012=t021=t111=t120=t210=t300=0,n011=B1k22(1+k2y)2+B2k2(1+k2y)2x2K2B4(1m2)2,  n002=B1x2k22(1+k2y)3B2k22(1+k2y)3(x2)2K2,n012=B1k223(1+k2y)3B2k223(1+k2y)32x2K2,  n003=B1x2k32(1+k2y)4+B2k32(1+k2y)4(x2)2K2,n110=B32,  n020=B21+k2y1K2,  n021=B2k23(1+k2y)21K2,n101=n200=n030=n111=n120=n102=n210=n201=n300=0,l101=C1(1m1)2,  l011=C2(1m2)2,l110=l200=l020=l002=l030=l003=l012=l021=l111=l120=l102=l210=l201=l300=0. (5.9)

    By introducing a reversible transformation

    (z1z2z3)=(101q21q22q23q31q32q33)(y1y2y3), (5.10)

    which q21q32q22q31+q22q33q23q320, the expression for the coefficient qij(i=2,3;j=1,2,3) is

    q21=L2L3L4L6+L1L3L5L6L1L2L26+L3L6ω2L23L4L5L23L252L2L3L5L6+L22L26+L23ω2,q22=(L23L4L1L3L6+L3L5L6L2L6)ωL23L252L2L3L5L6+L22L26+L23ω2,q23=L1L6L6η3L3L4L3L5L2L6L3η3,q31=L2L3L4L5L1L3L25L22L4L6+L1L2L5L6L1L3ω2L2L6ω2L23L252L2L3L5L6+L22L26+L23ω2,q32=ω(L1L2L6+L2L5L6L2L3L4L3L25L3ω2)L23L252L2L3L5L6+L22L26+L23ω2,q33=L2L4L1L5+L1η3+L5η3η23L3L5L2L6L3η3, (5.11)

    Li(i=1,2,,9) is defined in (4.11). So the standard type of system (5.7) can be written as

    (˙y1˙y2˙y3)=(0ω0ω0000η3)(y1y2y3)+(~F1(y1,y2,y3)~F2(y1,y2,y3)~F3(y1,y2,y3)), (5.12)

    where

    (~F1(y1,y2,y3)~F2(y1,y2,y3)~F3(y1,y2,y3))=(101q21q22q23q31q32q33)1(F1(z1,z2,z3)F2(z1,z2,z3)F3(z1,z2,z3)). (5.13)

    In Eq (5.13), the coefficients of polynomial ˜F(y1,y2,y3) are ˜tj1j2j3, ˜nj1j2j3 and ˜lj1j2j3.

    Based on the central manifold theory, use the center manifold Wc(0,0,0) existing at the origin to reduce the dimension of the system (5.12), that is

    Wc(0,0,0)={(y1,y2,y3)R3|y3=h(y1,y2),h(0,0)=0,Dh(0,0)=0}, (5.14)

    where

    h(y1,y2)=1v1+v23hv1v2y1y2+O((|y1|+|y2|)4). (5.15)

    Through calculation, we can obtain

    h10=h01=h11=0,h20=˜l200+˜l020ω2η3(1+ω4),h02=˜l020˜l200ω21+ω4,h30=1(1+ω4)(η23+ω6)(˜l200˜n101˜l300η32˜l020˜n200η3ω+˜l020˜n101ω2+˜l011˜l020ω3h30=˜l030ω3+2˜l200˜n200η3ω3˜l300η3ω4˜l011˜l200ω5˜l030ω7),h21=(˜l200+˜l020ω2)(˜l011η3+2˜t200η3ω+˜n110η23ω32˜t110ω4+˜n101η3ω52˜n020η3ω6)η3(1+ω4)(η23+ω6)h21=˜l020(2˜n110η3ω+˜n101ω32˜n020ω4)η23+ω6,h12=(˜l200+˜l020ω2)(2˜t110η3ω˜n101η23ω2+˜l011ω3+2˜n020η23ω3+2˜t200ω4+2˜n110η3ω6)η3(1+ω4)(η23+ω6)h12=˜l020(2˜n020η3ω+2˜n110ω4˜n101η3)η23+ω6,h03=(˜l200+˜l020ω2)(˜l011η23+ω2+˜n101ω3+2˜n200η3ω6)η3(1+ω4)(η23+ω6)h03=+˜l011˜l020η3˜l030η3+˜l300ω3+2˜l020˜n200ω4η23+ω6. (5.16)

    Correspondingly, the dynamic properties of the system are limited to the central flow Wc(0,0,0), and in conjunction with Eq (5.14), system (5.12) can be simplified as

    (˙y1˙y2)=(U(y1,y2)N(y1,y2)), (5.17)

    where

    U(y1,y2)=ωy2+1y21+2y1y2+3y22+4y31+5y21y2+6y1y22+7y32,N(y1,y2)=ωy1+ȷ1y21+ȷ2y1y2+ȷ3y22+ȷ4y31+ȷ5y21y2+ȷ6y1y22+ȷ7y32, (5.18)

    in the formula, we have

    1=˜t200,2=˜t110,3=˜t020,4=h20˜t101+˜t300,5=˜t210,6=h02˜t101,7=0,ȷ1=˜n200,ȷ2=˜n110,ȷ3=˜n020,ȷ4=˜n300,ȷ5=h20˜n011,ȷ6=0,ȷ7=h02˜n011+˜n030. (5.19)

    We introduce the partial derivative sign

    Uy1(ykΞ1)=Uy1,  3Uy21y2(ykΞ1)=Uy1y1y2,  2Ny22(ykΞ1)=Ny2y2,  , (5.20)

    where subscripts y1 and y2 indicate partial derivatives for the first and second variable, respectively. Based on Eq (5.18), it can be obtained that Uy1=0, Uy20, Ny10, Ny2=0, and Uy2Ny10. In addition, it ensures that the system (5.18) has pure virtual feature roots ±i|Uy2Ny1|. Thus, it can be determined that system (1.7) produces Hopf bifurcation; the direction of the bifurcation is determined by the following equation:

    QkΞ1=116ω(3+5+ȷ5+ȷ7)+116ω(13ȷ2ȷ3ȷ1ȷ2ȷ11). (5.21)

    Theorem 5.1. If dρdk1|k1=kΞ10, then system (1.7) will generate Hopf bifurcation at interior equilibrium P7. In addition, when dρdk1|k1=kΞ1<0, if QkΞ1<0 and 0<k1kΞ11, then system (1.7) will generate supercritical Hopf bifurcation and form a stable periodic orbit, or if QkΞ1>0 and 0<k1kΞ11, then system (1.7) will generate subcritical Hopf bifurcation and form a stable periodic orbit.

    In this section, we first discussed equilibria P1 to P7 of system (1.7) with distinct values of α, w1, and w2. Consider the parameter values as follows: ~r1=(2.8,3,3.2), ~r2=(2.8,3,3.2), ~c1=(0.1,0.2,0.3), ~c2=(0.5,0.6,0.7), ~a1=(0.1,0.2,0.3), ~a2=(0.2,0.3,0.4), ~e1=(0.2,0.3,0.4), ~e2=(0.3,0.4,0.5), and ˜d=(0.1,0.2,0.3). Tables 28 showed that the trivial equilibrium P1 retained constant at (0, 0, 0), the values of prey x1, prey x2, and predator y always maintained at 0; the values of prey x1 in P2 and prey x2 in P3 severally decreased with increasing w1 under the same α; the values of prey x1 and predator y in P4 increased with increasing w1, and for P5 the value of prey x2 and predator y rose with growing w1; the values of prey x1 and x2 in P6 decreased with growing w1; and for the same α, considering interior equilibrium P7, the values of prey x1, prey x2, and predator y decreased with growing w1.

    Table 2.  The trivial equilibrium P1 for k1=0.1, k2=0.7, q1=0.7, q2=0.5, q3=0.7, E1=0.3, E2=0.2, E3=0.2, K1=5, K2=5, m1=0.9, m2=0.3.
    w1 w2 P1 at α=0 P1 at α=0.3 P1 at α=0.6 P1 at α=0.9
    0 1 (0,0,0) (0,0,0) (0,0,0) (0,0,0)
    0.2 0.8 (0,0,0) (0,0,0) (0,0,0) (0,0,0)
    0.4 0.6 (0,0,0) (0,0,0) (0,0,0) (0,0,0)
    0.6 0.4 (0,0,0) (0,0,0) (0,0,0) (0,0,0)
    0.8 0.2 (0,0,0) (0,0,0) (0,0,0) (0,0,0)
    1 0 (0,0,0) (0,0,0) (0,0,0) (0,0,0)

     | Show Table
    DownLoad: CSV
    Table 3.  The axial equilibrium P2 for k1=0.1, k2=0.7, q1=0.7, q2=0.5, q3=0.7, E1=0.3, E2=0.2, E3=0.2, K1=5, K2=5, m1=0.9, m2=0.3.
    w1 w2 P2 at α=0 P2 at α=0.3 P2 at α=0.6 P2 at α=0.9
    0 1 (5.3393,0,0) (5.1224,0,0) (4.9144,0,0) (4.7148,0,0)
    0.2 0.8 (5.0521,0,0) (4.9280,0,0) (4.8069,0,0) (4.6888,0,0)
    0.4 0.6 (4.7804,0,0) (4.7409,0,0) (4.7017,0,0) (4.6629,0,0)
    0.6 0.4 (4.5230,0,0) (4.5608,0,0) (4.5988,0,0) (4.6372,0,0)
    0.8 0.2 (4.2788,0,0) (4.3872,0,0) (4.4980,0,0) (4.6116,0,0)
    1 0 (4.0469,0,0) (4.2197,0,0) (4.3994,0,0) (4.5861,0,0)

     | Show Table
    DownLoad: CSV
    Table 4.  The axial equilibrium P3 for k1=0.7, k2=0.1, q1=0.7, q2=0.5, q3=0.7, E1=0.3, E2=0.2, E3=0.2, K1=5, K2=5, m1=0.3, m2=0.9.
    w1 w2 P3 at α=0 P3 at α=0.3 P3 at α=0.6 P3 at α=0.9
    0 1 (0,5.5357,0) (0,5.3147,0) (0,5.1027,0) (0,4.8993,0)
    0.2 0.8 (0,5.2431,0) (0,5.1166,0) (0,4.9932,0) (0,4.8728,0)
    0.4 0.6 (0,4.9662,0) (0,4.9260,0) (0,4.8861,0) (0,4.8465,0)
    0.6 0.4 (0,4.7039,0) (0,4.7424,0) (0,4.7812,0) (0,4.8202,0)
    0.8 0.2 (0,4.4551,0) (0,4.5655,0) (0,4.6785,0) (0,4.7942,0)
    1 0 (0,4.2187,0) (0,4.3949,0) (0,4.5779,0) (0,4.7682,0)

     | Show Table
    DownLoad: CSV
    Table 5.  The axial equilibrium P4 for k1=0.3, k2=0.7, q1=0.7, q2=0.5, q3=0.7, E1=0.3, E2=0.2, E3=0.2, K1=20, K2=20, m1=0.2, m2=0.6.
    w1 w2 P4 at α=0 P4 at α=0.3 P4 at α=0.6 P4 at α=0.9
    0 1 (0.4800,0,0.1496) (0.5838,0,0.1767) (0.7059,0,0.2040) (0.8516,0,0.2317)
    0.2 0.8 (0.6222,0,0.1858) (0.6971,0,0.2022) (0.7802,0,0.2188) (0.8732,0,0.2354)
    0.4 0.6 (0.8000,0,0.2224) (0.8306,0,0.2280) (0.8623,0,0.2335) (0.8954,0,0.2391)
    0.6 0.4 (1.0286,0,0.2595) (0.9902,0,0.2539) (0.9534,0,0.2484) (0.9181,0,0.2428)
    0.8 0.2 (1.3333,0,0.2969) (1.1845,0,0.2801) (1.0551,0,0.2633) (0.9415,0,0.2465)
    1 0 (1.7600,0,0.3343) (1.4261,0,0.3062) (1.1692,0,0.2782) (0.9655,0,0.2502)

     | Show Table
    DownLoad: CSV
    Table 6.  The axial equilibrium P5 for k1=0.3, k2=0.7, q1=0.7, q2=0.5, q3=0.7, E1=0.3, E2=0.2, E3=0.2, K1=20, K2=20, m1=0.8, m2=0.6.
    w1 w2 P5 at α=0 P5 at α=0.3 P5 at α=0.6 P5 at α=0.9
    0 1 (0,0.1920,0.3638) (0,0.2298,0.3833) (0,0.2727,0.4029) (0,0.3220,0.4226)
    0.2 0.8 (0,0.2435,0.3899) (0,0.2697,0.4016) (0,0.2981,0.4134) (0,0.3291,0.4252)
    0.4 0.6 (0,0.3048,0.4160) (0,0.3150,0.4199) (0,0.3255,0.4239) (0,0.3363,0.4278)
    0.6 0.4 (0,0.3789,0.4422) (0,0.3668,0.4383) (0,0.3551,0.4343) (0,0.3437,0.4304)
    0.8 0.2 (0,0.4706,0.4683) (0,0.4268,0.4566) (0,0.3872,0.4448) (0,0.3513,0.4330)
    1 0 (0,0.5867,0.4942) (0,0.4970,0.4748) (0,0.4222,0.4553) (0,0.3590,0.4357)

     | Show Table
    DownLoad: CSV
    Table 7.  The axial equilibrium P6 for k1=0.3, k2=0.4, q1=0.7, q2=0.5, q3=0.7, E1=0.6, E2=0.1, E3=0.2, K1=5, K2=5, m1=0.9, m2=0.4.
    w1 w2 P6 at α=0 P6 at α=0.3 P6 at α=0.6 P6 at α=0.9
    0 1 (4.9826,3.8455,0) (4.6419,3.5356,0) (4.3385,3.2568,0) (4.0671,3.0043,0)
    0.2 0.8 (4.5369,3.4395,0) (4.3577,3.2745,0) (4.1901,3.1191,0) (4.0331,2.9724,0)
    0.4 0.6 (4.1543,3.0858,0) (4.1016,3.0366,0) (4.0500,2.9883,0) (3.9995,2.9408,0)
    0.6 0.4 (3.8232,2.7740,0) (3.8700,2.8184,0) (3.9177,2.8636,0) (3.9665,2.9097,0)
    0.8 0.2 (3.5351,2.4958,0) (3.6599,2.6172,0) (3.7926,2.7447,0) (3.9338,2.8789,0)
    1 0 (3.2834,2.2448,0) (3.4690,2.4307,0) (3.6742,2.6311,0) (3.9017,2.8485,0)

     | Show Table
    DownLoad: CSV
    Table 8.  The interior equilibrium P7 for k1=0.4, k2=0.5, q1=0.6, q2=0.4, q3=0.2, E1=0.2, E2=0.3, E3=0.2, K1=100, K2=100, m1=0.9, m2=0.3.
    w1 w2 P7 at α=0 P7 at α=0.3
    0 1 (1.7240,3.6276,1.1986) (1.5910,3.4597,1.0271)
    0.2 0.8 (1.6168,3.5123,1.0583) (1.5042,3.3223,0.9139)
    0.4 0.6 (1.5417,3.4296,0.9337) (1.4211,3.1793,0.8023)
    0.6 0.4 (1.4309,3.2886,0.7930) (1.3913,3.1413,0.7269)
    0.8 0.2 (1.4291,3.2451,0.6884) (1.3195,3.0467,0.6350)
    1 0 (1.3764,2.9890,0.5245) (1.2797,2.9766,0.5529)
    w1 w2 P7 at α=0.6 P7 at α=0.9
    0 1 (1.7066,3.7444,1.0113) (1.5325,3.4231,0.7999)
    0.2 0.8 (1.5291,3.4203,0.8678) (1.4899,3.3773,0.7762)
    0.4 0.6 (1.5273,3.4215,0.8306) (1.4789,3.3677,0.7638)
    0.6 0.4 (1.4894,3.3737,0.7779) (1.4685,3.3587,0.7516)
    0.8 0.2 (1.4594,3.3320,0.7279) (1.4590,3.3504,0.7397)
    1 0 (1.4328,3.2552,0.6687) (1.4124,3.2841,0.7105)

     | Show Table
    DownLoad: CSV

    Considering four sets of different initial values, it could be seen from Figure 1 that different orbits eventually converged to the same value, which concluded that the interior equilibrium of the system (1.7) fulfills the character of globally asymptotical stability. Figure 2 plotted the bifurcation graph of system (1.7) with the horizontal coordinates k1, and the Hopf bifurcation of the system occurred with k1 taking values in the range of 0.01k10.7. When 0.01k1<0.384, the system oscillates periodically, while it maintains a stable steady-state when 0.384<k10.7. Therefore, based on Figure 2, it could be concluded that the fear of prey x1 for predator y affected the stability of the system. We further observed that as k1 increased, the prey x1 density continued to decrease while the predator y density kept increasing. Thus, the result also suggested that greater fear of predators had a negative impact on prey populations while having a positive impact on predator populations. Correspondingly, Figures 3 and 4 showed the waveform plots and phase diagram at k1=0.1 and k1=0.7, respectively.

    Figure 1.  Global stability of the internal equilibrium P7 = (5.665, 1.668, 2.047) of system (1.7) is given by the following parameter values: α=1, w1+w2=1, A1=2.0, A2=2.0, B1=2.0, B2=2.0, k1=0.2, k2=0.1, q1=0.4, q2=0.4, q3=0.2, E1=0.2, E2=0.2, E3=0.2, A3=0.1, B3=0.1, A4=0.3, B4=0.6, K1=10, K2=10, m1=0.4, m2=0.4, C1=0.1, C2=0.2, C3=0.5.
    Figure 2.  Hopf bifurcation occurs as a bifurcation parameter k1, and the remaining parameters take the following values: α=1, w1+w2=1, A1=3.0, A2=3.0, B1=3.0, B2=3.0, k2=0.4, q1=0.6, q2=0.4, q3=0.2, E1=0.2, E2=0.2, E3=0.2, A3=0.2, B3=0.1, A4=0.3, B4=0.6, K1=10, K2=70, m1=0.9, m2=0.3, C1=0.1, C2=0.2, C3=0.5.
    Figure 3.  Waveform plots and phase diagram of system (1.7) with k1=0.1, and α=1, w1+w2=1.
    Figure 4.  Waveform plots and phase diagram of system (1.7) with k1=0.7, and α=1, w1+w2=1.

    In addition, Figure 5 also plots the bifurcation graph with changing m1. As can be seen in Figure 5, m1 took values from 0.3 to 1, in which the system also underwent a Hopf bifurcation. When the value m1 ranged from 0.3 to 0.657, the system (1.7) was stable; nevertheless, it would become unstable at 0.657<m11. Correspondingly, Figures 6 and 7 showed the waveform plots and phase diagram at m1=0.6 and m1=0.9, respectively.

    Figure 5.  Hopf bifurcation occurs as a bifurcation parameter of system (1.7) parameter m1, and the remaining parameters take the following values: α=1, w1+w2=1, A1=2.0, A2=2.0, B1=2.0, B2=2.0, k1=0.1, k2=0.4, q1=0.7, q2=0.4, q3=0.2, E1=0.2, E2=0.2, E3=0.3, A3=0.1, B3=0.1, A4=0.3, B4=0.6, K1=10, K2=70, m2=0.4, C1=0.1, C2=0.2, C3=0.5.
    Figure 6.  Waveform plots and phase diagram of system (1.7) with m1=0.6 and α=1, w1+w2=1.
    Figure 7.  Waveform plots and phase diagram of system (1.7) with m1=0.9 and α=1, w1+w2=1.

    Further, we find an interesting dynamic phenomenon through some numerical simulations. System (1.7) appears as a chaotic phenomenon, as shown in Figure 8.

    Figure 8.  Waveform plots and phase diagram of chaotic phenomena with the following parameter values: α=1, w1+w2=1, A1=2.0, A2=2.0, B1=3.0, B2=3.0, k1=0.2, k2=0.5, q1=0.6, q2=0.4, q3=0.2, E1=0.2, E2=0.3, E3=0.2, A3=0.2, B3=0.3, A4=0.3, B4=0.6, K1=10, K2=70, m1=0.9, m2=0.3, C1=0.1, C2=0.2, C3=0.5.

    In this work, we develop a model of one-predator and two-prey interactions in a fuzzy environment, examine the effects of fear and prey refuge on the system, and provide insight into the dynamic complexity. The proofs of the theoretical parts of this paper are based on system (1.7). It has been proven that all equilibria in system (1.7) are locally asymptotically stable, and interior equilibrium P7 is also globally asymptotically stable. We have been further concerned about the appearance and direction of Hopf bifurcation. With the support of theoretical research, our numerical simulations have been able to display a wealth of charts and graphs.

    First of all, different equilibria are displayed from Tables 28 with different α,w1,w2, respectively. Throughout Figure 1, we have verified the global asymptotical stability of interior equilibrium P7, and find that the system is from unstable to stable with the increase of fear k1, which demonstrates that the fear effect may be an important factor influencing the stability of the system (see Figures 24). Furthermore, it has also been observed that an increase in prey refuge m1 leads to oscillatory phenomena (see Figures 57). Finally, through studying the Hopf bifurcation, we have discovered some interesting biological phenomena, namely that system (1.7) appears to be in a chaotic state (see Figure 8).

    Xuyang Cao: Conceptualization, Investigation, Methodology, Validation, Writing-original draft, Formal analysis, Software; Qinglong Wang: Conceptualization, Methodology, Formal analysis, Writing-review and editing, Supervision; Jie Liu: Validation, Visualization, Data curation.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors thank the editor and referees for their careful reading and valuable comments.

    The work is supported by the Natural Science Foundation of Hubei Province (No. 2023AFB1095) and the National Natural Science Foundation of China (No. 12101211) and the Program for Innovative Research Team of the Higher Education Institution of Hubei Province (No. T201812) and the Teaching Research Project of Education Department of Hubei Province (No. 2022367) and the Graduate Education Innovation Project of Hubei Minzu University (Nos. MYK2024071, MYK2023042).

    The authors declare that there are no conflicts of interest regarding the publication of this paper.

    Definition 1. [34] Fuzzy set: A fuzzy set ˜ in a universe of discourse S is denoted by the set of pairs

    ˜={(s,μ˜(s)):sS},

    where the mapping μ˜:S[0,1] is the membership function of the fuzzy set ˜ and μ˜ is the membership value or degree of membership of sS in the fuzzy set ˜.

    Definition 2. [42] α-cut of fuzzy set: For any α(0,1], the α-cut of fuzzy set ˜ defined by α={s:μ˜(s))α} is a crisp set. For α=0 the support of ˜ is defined as 0=Supp(˜)=¯{sR,μ˜(s)>0}.

    Definition 3. [43] Fuzzy number: A fuzzy number satisfying the property S=R is called a convex fuzzy set.

    Definition 4. [44] Triangular fuzzy number: A triangular fuzzy number (TFN) ˜(b1,b2,b3) represent fuzzy set of the real line R satisfying the property that the membership function μ˜:R[0,1] can be espressed by

    μ˜={sb1b2b1  if  b1sb2,b3sb3b2  if  b2sb3,0otherwise.

    Hence, the α-cut of triangular fuzzy number meets boundedness and encapsulation on [L(α),R(α)], in which L(α)=infs:μ˜(s)α=b1+α(b2b1) and R(α)=sup{s:μ˜(s)α}=b3+α(b3b2).

    Lemma 1. [45] In weighted sum method, wj stands for the weight of jth objective. wjgj represent a utility function for jth objective, and the total utility function π is represented by

    π=ljwjgj,j=1,2,,l,

    where wj>0 and ljwj=1 are satisfied.


    Acknowledgments



    This work was supported by the National Natural Science Foundation of China (grant number 32202121).

    Conflict of interest



    The authors declare no competing interests.

    [1] Kaper JB, Nataro JP, Mobley H (2004) Pathogenic Escherichia coli. Nat Rev Microbiol 2: 123-140. https://doi.org/10.1038/nrmicro818
    [2] Wang YX, Ye ZZ, Si CY, et al. (2012) Application of aptamer based biosensors for detection of pathogenic microorganisms. Chinese J Anal Chem 40: 634-642. https://doi.org/10.1016/S1872-2040(11)60542-2
    [3] Cunha BA (2001) Antibiotic side effects. Med Clin N Am 85: 149. https://doi.org/10.1016/S0025-7125(05)70309-6
    [4] Rangarajan V, Clarke KG (2015) Process development and intensification for enhanced production of Bacillus lipopeptides. Biotechnol Genet Eng Revs 31: 46-68. https://doi.org/10.1080/02648725.2016.1166335
    [5] Wei YH, Wang LF, Chang JS, et al. (2003) Identification of induced acidification in iron-enriched cultures of Bacillus subtilis during biosurfactant fermentation. J Biosci Bioeng 96: 174-178. https://doi.org/10.1263/jbb.96.174
    [6] Seydlova G, Svobodova J (2008) Review of surfactin chemical properties and the potential biomedical applications. Cent Eur J Med 3: 123-133. https://doi.org/10.2478/s11536-008-0002-5
    [7] Das P, Mukherjee S, Sen R (2008) Antimicrobial potential of a lipopeptide biosurfactant derived from a marine Bacillus circulans. J Appl Microbiol 104: 1675-1684. https://doi.org/10.1111/j.1365-2672.2007.03701.x
    [8] Boettcher C, Kell H, Holzwarth JF, et al. (2010) Flexible loops of thread-like micelles are formed upon interaction of L-alpha-dimyristoyl-phosphatidylcholine with the biosurfactant surfactin as revealed by cryo-electron tomography. Biophys Chem 149: 22-27. https://doi.org/10.1016/j.bpc.2010.03.006
    [9] Sachdev DP, Cameotra SS (2013) Biosurfactants in agriculture. Appl Microbiol Biot 97: 1005-1016. https://doi.org/10.1007/s00253-012-4641-8
    [10] Zhang YY, Liu C, Dong B, et al. (2015) Anti-inflammatory activity and mechanism of surfactin in lipopolysaccharide-activated macrophages. Inflammation 38: 756-764. https://doi.org/10.1007/s10753-014-9986-y
    [11] Gang HZ, Liu JF, Mu BZ (2011) Molecular dynamics study of surfactin monolayer at the air/water interface. J Phys Chem B 115: 12770-12777. https://doi.org/10.1021/jp206350j
    [12] Barros F, de Quadros CP, Marstica MR, et al. (2007) Surfactin: Chemical, technological and functional properties for food applications. Quim Nova 30: 409-414. https://doi.org/10.1590/S0100-40422007000200031
    [13] Zhai SW, Chen XH, Wang MH (2017) Effects of different levels of dietary surfactin supplementation on intestinal morphology, and intestinal microflora of growth retarded marbled eel juveniles (Anguilla marmaorata). Isr J Aquacult-Bamid 69: 1433.
    [14] Xia L, Wen JP (2022) Available strategies for improving the biosynthesis of surfactin: A review. Crit Rev Biotechnol . https://doi.org/10.1080/07388551.2022.2095252
    [15] Arima K, Kakinuma A, Tamura G (1968) Surfactin, a crystalline peptidelipid surfactant produced by Bacillus subtilis: isolation, characterization and its inhibition of fibrin clot formation. Biochem Bioph Res Co 31: 488-494. https://doi.org/10.1016/0006-291X(68)90503-2
    [16] Hue N, Serani L, Laprevote O (2001) Structural investigation of cyclic peptidolipids from Bacillus subtilis by high-energy tandem mass spectrometry. Rapid Commun Mass Sp 15: 203-209. https://doi.org/10.1002/1097-0231(20010215)15:3<203::AID-RCM212>3.0.CO;2-6
    [17] Fei D, Zhou GW, Yu ZQ, et al. (2020) Low-toxic and nonirritant biosurfactant surfactin and its performances in detergent formulations. J Surfactants Deterg 23: 109-118. https://doi.org/10.1002/jsde.12356
    [18] Long XW, He N, He YK, et al. (2017) Biosurfactant surfactin with pH-regulated emulsification activity for efficient oil separation when used as emulsifier. Bioresource Technol 241: 200-206. https://doi.org/10.1016/j.biortech.2017.05.120
    [19] Shakerifard P, Gancel F, Jacques P, et al. (2009) Effect of different Bacillus subtilis lipopeptides on surface hydrophobicity and adhesion of Bacillus cereus 98/4 spores to stainless steel and Teflon. Biofouling 25: 533-541. https://doi.org/10.1080/08927010902977943
    [20] Marcelino L, Puppin-Rontani J, Coutte F, et al. (2019) Surfactin application for a short period (10/20 s) increases the surface wettability of sound dentin. Amino Acids 51: 1233-1240. https://doi.org/10.1007/s00726-019-02750-
    [21] Yang ZY, Zu YQ, Zhu JS, et al. (2020) Application of biosurfactant surfactin as a pH-switchable biodemulsifier for efficient oil recovery from waste crude oil. Chemosphere 240: 124946. https://doi.org/10.1016/j.chemosphere.2019.124946
    [22] Zouari R, Besbes S, Ellouze-Chaabouni S, et al. (2016) Cookies from composite wheat-sesame peels flours: Dough quality and effect of Bacillus subtilis SPB1 biosurfactant addition. Food Chem 194: 758-769. https://doi.org/10.1016/j.foodchem.2015.08.064
    [23] Ohadi M, Shahravan A, Dehghannoudeh N, et al. (2020) Potential use of microbial surfactant in microemulsion drug delivery system: A systematic review. Drug Des Dev Ther 14: 541-550. https://doi.org/10.2147/DDDT.S232325
    [24] Heerklotz H, Seelig J (2001) Detergent-like action of the antibiotic peptide surfactin on lipid membranes. Biophys J 81: 1547-1554. https://doi.org/10.1016/S0006-3495(01)75808-0
    [25] Yeh MS, Wei YH, Chang JS (2005) Enhanced production of surfactin from Bacillus subtilis by addition of solid carriers. Biotechnol Progr 21: 1329-1334. https://doi.org/10.1021/bp050040c
    [26] Liu Q, Lin JZ, Wang WD, et al. (2015) Production of surfactin isoforms by Bacillus subtilis BS-37 and its applicability to enhanced oil recovery under laboratory conditions. Biochem Eng J 93: 31-37. https://doi.org/10.1016/j.bej.2014.08.023
    [27] Das P, Yang XP, Ma L (2014) Analysis of biosurfactants from industrially viable Pseudomonas strain isolated from crude oil suggests how rhamnolipids congeners affect emulsification property and antimicrobial activity. Front Microbiol 5: 696. https://doi.org/10.3389/fmicb.2014.00696
    [28] Bernheimer AW, Avigad LS (1970) Nature and properties of a cytolytic agent produced by Bacillus subtilis. J Gen Appl Microbiol 61: 361-369. https://doi.org/10.1099/00221287-61-3-361
    [29] Bouffioux O, Berquand A, Eeman M, et al. (2007) Molecular organization of surfactin-phospholipid monolayers: Effect of phospholipid chain length and polar head. Bba-Biomembranes 1768: 1758-1768. https://doi.org/10.1016/j.bbamem.2007.04.015
    [30] Maget-Dana R, Ptak M (1995) Interactions of surfactin with membrane models. Biophys J 68: 1937-1943. https://doi.org/10.1016/S0006-3495(95)80370-X
    [31] Carrillo C, Teruel JA, Aranda FJ, et al. (2003) Molecular mechanism of membrane permeabilization by the peptide antibiotic surfactin. Bba-Biomembranes 1611: 91-97. https://doi.org/10.1016/S0005-2736(03)00029-4
    [32] Thimon L, Peypoux F, Maget-Dana R, et al. (1992) Interactions of bioactive lipopeptides, iturin A and surfactin from Bacillus subtilis. Biotechnol Appl Bioc 16: 144-151.
    [33] Hwang YH, Park BK, Lim JH, et al. (2007) Lipopolysaccharide-binding and neutralizing activities of surfactin C in experimental models of septic shock. Eur J Pharmacol 556: 166-171. https://doi.org/10.1016/j.ejphar.2006.10.031
    [34] Hwang MH, Lim JH, Yun HI, et al. (2005) Surfactin C inhibits the lipopolysaccharide-induced transcription of interleukin-1 beta and inducible nitric oxide synthase and nitric oxide production in murine RAW 264.7 cells. Biotechnol Lett 27: 1605-1608. https://doi.org/10.1007/s10529-005-2515-1
    [35] Huang X, Wei Z, Gao X, et al. (2008) Study on antiviral activity of surfactin on Porcine parvovirus in vitro. Acta Agriculturae Zhejiangensis 20: 349-352.
    [36] Vollenbroich D, Ozel M, Vater J, et al. (1997) Mechanism of inactivation of enveloped viruses by the biosurfactant surfactin from Bacillus subtilis. Biologicals 25: 289-297. https://doi.org/10.1006/biol.1997.0099
    [37] Vollenbroich D, Pauli G, Ozel M, et al. (1997) Antimycoplasma properties and application in cell culture of surfactin, a lipopeptide antibiotic from Bacillus subtilis. Appl Environ Microb 63: 44-49. https://doi.org/10.1128/AEM.63.1.44-49.1997
    [38] Park SY, Kim JH, Lee YJ, et al. (2013) Surfactin suppresses TPA-induced breast cancer cell invasion through the inhibition of MMP-9 expression. Int J Oncol 42: 287-296. https://doi.org/10.3892/ijo.2012.1695
    [39] Peterlik M, Grant WB, Cross HS (2009) Calcium, vitamin D and cancer. Anticancer Res 29: 3687-3698.
    [40] Lee JH, Nam SH, Seo WT, et al. (2012) The production of surfactin during the fermentation of cheonggukjang by potential probiotic Bacillus subtilis CSY191 and the resultant growth suppression of MCF-7 human breast cancer cells. Food Chem 131: 1347-1354. https://doi.org/10.1016/j.foodchem.2011.09.133
    [41] Wu YS, Ngai SC, Goh BH, et al. (2017) Anticancer activities of surfactin and potential application of nanotechnology assisted surfactin delivery. Front Pharmacol 8: 761. https://doi.org/10.3389/fphar.2017.00761
    [42] Chen XY, Lu YJ, Shan MY, et al. (2022) A mini-review: mechanism of antimicrobial action and application of surfactin. World J Microb Biot 38: 143. https://doi.org/10.1007/s11274-022-03323-3
    [43] Huang XQ, Gao XP, Zheng LY, et al. (2009) Optimization of sterilization of salmonella enteritidis in meat by surfactin and iturin using a response surface method. Int J Pept Res Ther 15: 61-67. https://doi.org/10.1007/s10989-008-9164-x
    [44] Wang Y, Tian JH, Shi FF, et al. (2021) Protective effect of surfactin on copper sulfate-induced inflammation, oxidative stress, and hepatic injury in zebrafish. Microbiol Immunol 65: 410-421. https://doi.org/10.1111/1348-0421.12924
    [45] Horng YB, Yu YH, Dybus A, et al. (2019) Antibacterial activity of Bacillus species-derived surfactin on Brachyspira hyodysenteriae and Clostridium perfringens. Amb Express 9: 188. https://doi.org/10.1186/s13568-019-0914-2
    [46] Rowland I, Gibson G, Heinken A, et al. (2018) Gut microbiota functions: metabolism of nutrients and other food components. Eur J Nutr 57: 1-24. https://doi.org/10.1007/s00394-017-1445-8
    [47] Jandhyala SM, Talukdar R, Subramanyam C, et al. (2015) Role of the normal gut microbiota. World J Gastroentero 21: 8787-8803. https://doi.org/10.3748/wjg.v21.i29.8787
    [48] Yu YH, Wu CM, Chen WJ, et al. (2021) Effectiveness of Bacillus licheniformis-fermented products and their derived antimicrobial lipopeptides in controlling coccidiosis in broilers. Animals-Basel 11: 3576. https://doi.org/10.3390/ani11123576
    [49] Zhai SW, Shi QC, Chen XH (2016) Effects of dietary surfactin supplementation on growth performance, intestinal digestive enzymes activities and some serum biochemical parameters of tilapia (Oreochromis niloticus) fingerlings. Ital J Anim Sci 15: 318-324. https://doi.org/10.1080/1828051X.2016.1175325
    [50] Piewngam P, Zheng Y, Nguyen TH, et al. (2018) Pathogen elimination by probiotic Bacillus via signalling interference. Nature 562: 532. https://doi.org/10.1038/s41586-018-0616-y
    [51] Bais HP, Fall R, Vivanco JM (2004) Biocontrol of Bacillus subtilis against infection of Arabidopsis roots by Pseudomonas syringae is facilitated by biofilm formation and surfactin production. Plant Physiol 134: 307-319. https://doi.org/10.1104/pp.103.028712
    [52] Huang X, Wei Z, Gao X, et al. (2009) Study of antiviral activity of surfactin on pseudorabies virus in vitro. J Biol 26: 41-43.
    [53] Yuan LF, Zhang SA, Wang YH, et al. (2018) Surfactin inhibits membrane fusion during invasion of epithelial cells by enveloped viruses. J Virol 92. https://doi.org/10.1128/JVI.00809-18
    [54] Zhou SN, Liu G, Wu SM (2020) Marine bacterial surfactin CS30-2 induced necrosis-like cell death in Huh7.5 liver cancer cells. J Oceanol Limnol 38: 826-833. https://doi.org/10.1007/s00343-019-9129-2
    [55] Wang D, Liu Y, Yang Z, et al. (2008) Application of surfactin in microbial enhanced oil recovery. Acta Petrolei Sinica 29: 111-115.
    [56] He Z, Zhao H, Lu Z (2017) Effect of surfactin as surfactant on physical and oxidation stability of O/W DHA-rich algae oil emulsion. Food Sci 38: 146-151.
    [57] Ganesan NG, Rangarajan V (2021) A kinetics study on surfactin production from Bacillus subtilis MTCC 2415 for application in green cosmetics. Biocatal Agr Biotech 33. https://doi.org/10.1016/j.bcab.2021.102001
    [58] Whang LM, Liu P, Ma CC, et al. (2008) Application of biosurfactants, rhamnolipid, and surfactin, for enhanced biodegradation of diesel-contaminated water and soil. J Hazard Mater 151: 155-163. https://doi.org/10.1016/j.jhazmat.2007.05.063
    [59] Mulligan CN (2005) Environmental applications for biosurfactants. Environ Pollut 133: 183-198. https://doi.org/10.1016/j.envpol.2004.06.009
    [60] Chen WC, Juang RS, Wei YH (2015) Applications of a lipopeptide biosurfactant, surfactin, produced by microorganisms. Biochem Eng J 103: 158-169. https://doi.org/10.1016/j.bej.2015.07.009
    [61] Kaneda T (1977) Fatty acids of the genus Bacillus: an example of branched-chain preference. Bacteriol Rev 41: 391-418. https://doi.org/10.1128/MMBR.41.2.391-418.1977
    [62] Youssef NH, Duncan KE, McInerney MJ (2005) Importance of 3-hydroxy fatty acid composition of lipopeptides for biosurfactant activity. Appl Environ Microb 71: 7690-7695. https://doi.org/10.1128/AEM.71.12.7690-7695.2005
    [63] Ongena M, Jacques P (2008) Bacillus lipopeptides: versatile weapons for plant disease biocontrol. Trends Microbiol 16: 115-125. https://doi.org/10.1016/j.tim.2007.12.009
    [64] Kopp F, Marahiel MA (2007) Macrocyclization strategies in polyketide and nonribosomal peptide biosynthesis. Nat Prod Rep 24: 735-749. https://doi.org/10.1039/b613652b
    [65] Hamoen LW, Venema G, Kuipers OP (2003) Controlling competence in Bacillus subtilis: shared use of regulators. Microbiol-Sgm 149: 9-17. https://doi.org/10.1099/mic.0.26003-0
    [66] Nakano MM, Corbell N, Besson J, et al. (1992) Isolation and characterization of sfp: a gene that functions in the production of the lipopeptide biosurfactant, surfactin, in Bacillus subtilis. Mol Gen Genet 232: 313-321. https://doi.org/10.1007/BF00280011
    [67] Lambalot RH, Gehring AM, Flugel RS, et al. (1996) A new enzyme superfamily-the phosphopantetheinyl transferases. Chem Biol 3: 923-936. https://doi.org/10.1016/S1074-5521(96)90181-7
    [68] Hu FX, Liu YY, Li S (2019) Rational strain improvement for surfactin production: enhancing the yield and generating novel structures. Microb Cell Fact 18: 42. https://doi.org/10.1186/s12934-019-1089-x
    [69] Banat IM, Satpute SK, Cameotra SS, et al. (2014) Cost effective technologies and renewable substrates for biosurfactants' production. Front Microbiol 5: 697. https://doi.org/10.3389/fmicb.2014.00697
    [70] Zhi Y, Wu Q, Xu Y (2017) Genome and transcriptome analysis of surfactin biosynthesis in Bacillus amyloliquefaciens MT45. Sci Rep-Uk 7: 40976. https://doi.org/10.1038/srep40976
    [71] Abdel-Mawgoud AM, Aboulwafa MM, Hassouna N (2008) Optimization of surfactin production by Bacillus subtilis isolate BS5. Appl Biochem Biotech 150: 305-325. https://doi.org/10.1007/s12010-008-8155-x
    [72] Amani H, Haghighi M, Keshtkar MJ (2013) Production and optimization of microbial surfactin by Bacillus subtilis for ex situ enhanced oil recovery. Petrol Sci Technol 31: 1249-1258. https://doi.org/10.1080/10916466.2010.542416
    [73] Yang N, Wu Q, Xu Y (2020) Fermentation optimization for the production of surfactin by Bacillus amyloliquefaciens. Chinese Biotechnol 40: 51-58.
    [74] Satpute SK, Banat IM, Dhakephalkar PK, et al. (2010) Biosurfactants, bioemulsifiers and exopolysaccharides from marine microorganisms. Biotechnol Adv 28: 436-450. https://doi.org/10.1016/j.biotechadv.2010.02.006
    [75] Nazareth TC, Zanutto CP, Tripathi L, et al. (2020) The use of low-cost brewery waste product for the production of surfactin as a natural microbial biocide. Biotechnol Rep 28: e537. https://doi.org/10.1016/j.btre.2020.e00537
    [76] Ramirez IM, Vaz DA, Banat IM, et al. (2016) Hydrolysis of olive mill waste to enhance rhamnolipids and surfactin production. Bioresource Technol 205: 1-6. https://doi.org/10.1016/j.biortech.2016.01.016
    [77] Zhi Y, Wu Q, Xu Y (2017) Production of surfactin from waste distillers' grains by co-culture fermentation of two Bacillus amyloliquefaciens strains. Bioresource Technol 235: 96-103. https://doi.org/10.1016/j.biortech.2017.03.090
    [78] Willenbacher J, Mohr T, Henkel M, et al. (2016) Substitution of the native srfA promoter by constitutive P-veg in two B. subtilis strains and evaluation of the effect on Surfactin production. J Biotechnol 224: 14-17. https://doi.org/10.1016/j.jbiotec.2016.03.002
    [79] Dhali D, Coutte F, Arias AA, et al. (2017) Genetic engineering of the branched fatty acid metabolic pathway of Bacillus subtilis for the overproduction of surfactin C-14 isoform. Biotechnol J 12: 1600574. https://doi.org/10.1002/biot.201600574
    [80] Wang MM, Yu HM, Shen ZY (2019) Antisense RNA-based strategy for enhancing surfactin production in Bacillus subtilis TS1726 via overexpression of the unconventional biotin carboxylase II to enhance ACCase activity. Acs Synth Biol 8: 251-256. https://doi.org/10.1021/acssynbio.8b00459
    [81] Wang CY, Cao YX, Wang YP, et al. (2019) Enhancing surfactin production by using systematic CRISPRi repression to screen amino acid biosynthesis genes in Bacillus subtilis. Microb Cell Fact 18: 90. https://doi.org/10.1186/s12934-019-1139-4
    [82] Wang MM, Yu HM, Li X, et al. (2020) Single-gene regulated non-spore-forming Bacillus subtilis: Construction, transcriptome responses, and applications for producing enzymes and surfactin. Metab Eng 62: 235-248. https://doi.org/10.1016/j.ymben.2020.08.008
    [83] Wu Q, Zhi Y, Xu Y (2019) Systematically engineering the biosynthesis of a green biosurfactant surfactin by Bacillus subtilis 168. Metab Eng 52: 87-97. https://doi.org/10.1016/j.ymben.2018.11.004
    [84] Li X, Yang H, Zhang DL, et al. (2015) Overexpression of specific proton motive force-dependent transporters facilitate the export of surfactin in Bacillus subtilis. J Ind Microbiol Biot 42: 93-103. https://doi.org/10.1007/s10295-014-1527-z
    [85] Jiao S, Li X, Yu HM, et al. (2017) In situ enhancement of surfactin biosynthesis in Bacillus subtilis using novel artificial inducible promoters. Biotechnol Bioeng 114: 832-842. https://doi.org/10.1002/bit.26197
    [86] Atwa NA, El-Shatoury E, Elazzazy A, et al. (2013) Enhancement of surfactin production by Bacillus velezensis NRC-1 strain using a modified bench-top bioreactor. J Food Agric Environ 11: 169-174.
    [87] Jung J, Yu KO, Ramzi AB, et al. (2012) Improvement of surfactin production in Bacillus subtilis using synthetic wastewater by overexpression of specific extracellular signaling peptides, comX and phrC. Biotechnol Bioeng 109: 2349-2356. https://doi.org/10.1002/bit.24524
    [88] Zhang F, Huo KY, Song XY, et al. (2020) Engineering of a genome-reduced strain Bacillus amyloliquefaciens for enhancing surfactin production. Microb Cell Fact 19: 223. https://doi.org/10.1186/s12934-020-01485-z
    [89] Cheng JT, Guan CR, Cui WJ, et al. (2016) Enhancement of a high efficient autoinducible expression system in Bacillus subtilis by promoter engineering. Protein Expres Purif 127: 81-87. https://doi.org/10.1016/j.pep.2016.07.008
    [90] Pottathil M, Jung A, Lazazzera BA (2008) CSF, a species-specific extracellular signaling peptide for communication among strains of Bacillus subtilis and Bacillus mojavensis. J Bacteriol 190: 4095-4099. https://doi.org/10.1128/JB.00187-08
    [91] Shank EA, Kolter R (2011) Extracellular signaling and multicellularity in Bacillus subtilis. Curr Opin Microbiol 14: 741-747. https://doi.org/10.1016/j.mib.2011.09.016
    [92] Ohsawa T, Tsukahara K, Sato T, et al. (2006) Superoxide stress decreases expression of srfA through inhibition of transcription of the comQXP quorum-sensing locus in Bacillus subtilis. J Biochem 139: 203-211. https://doi.org/10.1093/jb/mvj023
    [93] Coutte F, Niehren J, Dhali D, et al. (2015) Modeling leucine's metabolic pathway and knockout prediction improving the production of surfactin, a biosurfactant from Bacillus subtilis. Biotechnol J 10: 1216-1234. https://doi.org/10.1002/biot.201400541
    [94] Hayashi K, Ohsawa T, Kobayashi K, et al. (2005) The H2O2 stress-responsive regulator PerR positively regulates srfA expression in Bacillus subtilis. J Bacteriol 187: 6659-6667. https://doi.org/10.1128/JB.187.19.6659-6667.2005
    [95] Tsuge K, Ohata Y, Shoda M (2001) Gene yerP, involved in surfactin self-resistance in Bacillus subtilis. Antimicrob Agents Ch 45: 3566-3573. https://doi.org/10.1128/AAC.45.12.3566-3573.2001
    [96] Nah HJ, Pyeon HR, Kang SH, et al. (2017) Cloning and heterologous expression of a large-sized natural product biosynthetic gene cluster in Streptomyces species. Front Microbiol 8: 394. https://doi.org/10.3389/fmicb.2017.00394
    [97] Weihmann R, Domrose A, Drepper T, et al. (2020) Protocols for yTREX/Tn5-based gene cluster expression in Pseudomonas putida. Microb Biotechnol 13: 250-262. https://doi.org/10.1111/1751-7915.13402
    [98] Yamanaka K, Reynolds KA, Kersten RD, et al. (2014) Direct cloning and refactoring of a silent lipopeptide biosynthetic gene cluster yields the antibiotic taromycin A. P Natl Acad Sci Usa 111: 1957-1962. https://doi.org/10.1073/pnas.1319584111
    [99] Lee N, Larionov V, Kouprina N (2015) Highly efficient CRISPR/Cas9-mediated TAR cloning of genes and chromosomal loci from complex genomes in yeast. Nucleic Acids Res 43: e55. https://doi.org/10.1093/nar/gkv112
    [100] Krispin O, Allmansberger R (1998) The Bacillus subtilis AraE protein displays a broad substrate specificity for several different sugars. J Bacteriol 180: 3250-3252. https://doi.org/10.1128/JB.180.12.3250-3252.1998
    [101] Hu FX, Liu YY, Lin JZ, et al. (2020) Efficient production of surfactin from xylose-rich corncob hydrolysate using genetically modified Bacillus subtilis 168. Appl Microbiol Biot 104: 4017-4026. https://doi.org/10.1007/s00253-020-10528-9
    [102] Naseri G, Koffas M (2020) Application of combinatorial optimization strategies in synthetic biology. Nat Commun 11: 2446. https://doi.org/10.1038/s41467-020-16175-y
    [103] Stephens C, Christen B, Fuchs T, et al. (2007) Genetic analysis of a novel pathway for D-xylose metabolism in Caulobacter crescentus. J Bacteriol 189: 2181-2185. https://doi.org/10.1128/JB.01438-06
    [104] Nitschke M, Pastore GM (2004) Biosurfactant production by Bacillus subtilis using cassava-processing effluent. Appl Biochem Biotech 112: 163-172. https://doi.org/10.1385/ABAB:112:3:163
  • This article has been cited by:

    1. Yuan Tian, Hua Guo, Wenyu Shen, Xinrui Yan, Jie Zheng, Kaibiao Sun, Dynamic analysis and validation of a prey-predator model based on fish harvesting and discontinuous prey refuge effect in uncertain environments, 2025, 33, 2688-1594, 973, 10.3934/era.2025044
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4631) PDF downloads(359) Cited by(40)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog