Review Topical Sections

Targeting biofilms in medical devices using probiotic cells: a systematic review

  • Biofilm-related infections are becoming a major clinical problem due to the increasingly widespread ability of pathogens to develop persistent biofilms in medical devices. The inadequate response of conventional antimicrobial strategies to counteract biofilm development demands urgent alternatives. An increasing interest in promoting a natural approach to health has intensified the research in the field of probiotics to battle pathogenic biofilms. This study aims to summarize the recent evidence supporting the effects of probiotic cells on the control and prevention of medical device-associated biofilms using a PRISMA-oriented (Preferred Reporting Items for Systematic reviews and Meta-Analyses) systematic search. This review demonstrated that probiotic cells have the potential to be used as biocontrol agents against biofilm formation by a broad spectrum of microorganisms. The restriction of biofilm growth caused by probiotics seemed to be strain-specific and independent of the antibiofilm strategy used (displacement, exclusion or competition). Lactobacillus, Lactococcus and Streptococcus were the most studied genus of probiotics and those with the higher capability to hinder biofilm formation, causing reductions up to 99.9%. These findings will pave the way to further experiments on the topic so that probiotic cells may become part of the clinical arsenal for the prevention and treatment of medical device-associated infections.

    Citation: Fábio M. Carvalho, Rita Teixeira-Santos, Filipe J. M. Mergulhão, Luciana C. Gomes. Targeting biofilms in medical devices using probiotic cells: a systematic review[J]. AIMS Materials Science, 2021, 8(4): 501-523. doi: 10.3934/matersci.2021031

    Related Papers:

    [1] Nattawan Chuchalerm, Wannika Sawangtong, Benchawan Wiwatanapataphee, Thanongchai Siriapisith . Study of Non-Newtonian blood flow - heat transfer characteristics in the human coronary system with an external magnetic field. Mathematical Biosciences and Engineering, 2022, 19(9): 9550-9570. doi: 10.3934/mbe.2022444
    [2] B. Wiwatanapataphee, D. Poltem, Yong Hong Wu, Y. Lenbury . Simulation of Pulsatile Flow of Blood in Stenosed Coronary Artery Bypass with Graft. Mathematical Biosciences and Engineering, 2006, 3(2): 371-383. doi: 10.3934/mbe.2006.3.371
    [3] Lorenzo Civilla, Agnese Sbrollini, Laura Burattini, Micaela Morettini . An integrated lumped-parameter model of the cardiovascular system for the simulation of acute ischemic stroke: description of instantaneous changes in hemodynamics. Mathematical Biosciences and Engineering, 2021, 18(4): 3993-4010. doi: 10.3934/mbe.2021200
    [4] Benchawan Wiwatanapataphee, Yong Hong Wu, Thanongchai Siriapisith, Buraskorn Nuntadilok . Effect of branchings on blood flow in the system of human coronary arteries. Mathematical Biosciences and Engineering, 2012, 9(1): 199-214. doi: 10.3934/mbe.2012.9.199
    [5] Fan He, Minru Li, Xinyu Wang, Lu Hua, Tingting Guo . Numerical investigation of quantitative pulmonary pressure ratio in different degrees of stenosis. Mathematical Biosciences and Engineering, 2024, 21(2): 1806-1818. doi: 10.3934/mbe.2024078
    [6] Scott R. Pope, Laura M. Ellwein, Cheryl L. Zapata, Vera Novak, C. T. Kelley, Mette S. Olufsen . Estimation and identification of parameters in a lumped cerebrovascular model. Mathematical Biosciences and Engineering, 2009, 6(1): 93-115. doi: 10.3934/mbe.2009.6.93
    [7] Honghui Zhang, Jun Xia, Yinlong Yang, Qingqing Yang, Hongfang Song, Jinjie Xie, Yue Ma, Yang Hou, Aike Qiao . Branch flow distribution approach and its application in the calculation of fractional flow reserve in stenotic coronary artery. Mathematical Biosciences and Engineering, 2021, 18(5): 5978-5994. doi: 10.3934/mbe.2021299
    [8] Xintong Wu, Yingyi Geng, Xinhong Wang, Jucheng Zhang, Ling Xia . Continuous extraction of coronary artery centerline from cardiac CTA images using a regression-based method. Mathematical Biosciences and Engineering, 2023, 20(3): 4988-5003. doi: 10.3934/mbe.2023231
    [9] Ziyu Jin, Ning Li . Diagnosis of each main coronary artery stenosis based on whale optimization algorithm and stacking model. Mathematical Biosciences and Engineering, 2022, 19(5): 4568-4591. doi: 10.3934/mbe.2022211
    [10] Hongtao Liu, Shuqin Liu, Xiaoxu Ma, Yunpeng Zhang . A numerical model applied to the simulation of cardiovascular hemodynamics and operating condition of continuous-flow left ventricular assist device. Mathematical Biosciences and Engineering, 2020, 17(6): 7519-7543. doi: 10.3934/mbe.2020384
  • Biofilm-related infections are becoming a major clinical problem due to the increasingly widespread ability of pathogens to develop persistent biofilms in medical devices. The inadequate response of conventional antimicrobial strategies to counteract biofilm development demands urgent alternatives. An increasing interest in promoting a natural approach to health has intensified the research in the field of probiotics to battle pathogenic biofilms. This study aims to summarize the recent evidence supporting the effects of probiotic cells on the control and prevention of medical device-associated biofilms using a PRISMA-oriented (Preferred Reporting Items for Systematic reviews and Meta-Analyses) systematic search. This review demonstrated that probiotic cells have the potential to be used as biocontrol agents against biofilm formation by a broad spectrum of microorganisms. The restriction of biofilm growth caused by probiotics seemed to be strain-specific and independent of the antibiofilm strategy used (displacement, exclusion or competition). Lactobacillus, Lactococcus and Streptococcus were the most studied genus of probiotics and those with the higher capability to hinder biofilm formation, causing reductions up to 99.9%. These findings will pave the way to further experiments on the topic so that probiotic cells may become part of the clinical arsenal for the prevention and treatment of medical device-associated infections.



    The report released by the American Heart Association states that cardiovascular disease maintains its status as the primary global cause of death, claiming over 17.6 million lives annually, with forecasts suggesting this figure will rise to exceed 23.6 million by 2030 [1]. The coronary artery, the vascular system dedicated to providing nutrients and energy to the heart, has been the main research focus in recent years.

    The lumped parameter model, also known as the 0D model, was first developed by Otto Frank in 1898. It is the simplest mathematical model to describe the cardiovascular system and it is frequently employed to simulate the systemic characteristics of the cardiovascular system. The 0D model is known for its simple structure, easy-to-understand features, low computational cost, and efficient computing power. It has a wide range of simulations and can cover many fields and application scenarios. Compared with the coronary artery multi-scale model, its accuracy is relatively low and difficult to handle in the simulation of some complex physiological conditions, but it is still an effective modeling method in many fields [2]. Rideout et al. simulated the human pulmonary circulation system [3], accounting for the interactive effects of left atrial and right ventricular contractions on pressures in the pulmonary veins. They also developed an analog computer model to research the human cardiovascular system[4,5], with a particular emphasis on pressure and flow dynamics within veins, encompassing gravitational effects, respiration, venous collapse, and valve function. Liang et al. [6] developed a human cardiovascular system lumped parameter model featuring a detailed compartmental representation of both the main vascular circulations and the heart.

    Wang et al. [7] employed a closed-loop circuit model to simulate the left coronary arterial system and studied the impacts of stenosis, heart rate, and other factors on coronary hemodynamics. Mantero et al. [8] studied the role of the coronary artery bed in the cardiovascular system. Pietrabissa et al. [9] improved the cardiovascular system closed-loop circuit model, incorporating a more intricate depiction of the left coronary artery, which was also utilized to assess the coronary artery hemodynamics of surgical reconstruction. Duanmu et al. [10] established a new closed-loop lumped parameter model for the whole coronary network (including the right coronary artery). They first proposed the head loss model, which describes the physiological phenomenon of a sudden pressure drop at the coronary artery's entrance.

    Since considering the entire coronary artery becomes prohibitively expensive computationally, the model needs an outflow boundary condition to simulate the truncated tree. The three-element Windkessel model, put forth by Westerhof et al. [11], is the most classical lumped parameter model and is widely employed as the outflow boundary condition. At present, there is no strict theoretical support for the selection of parameters in the Windkessel model, which is usually obtained through experiments. Olufsen et al. [12,13,14] proposed the structured tree model that can also serve as a boundary condition, with its evaluation based on the root radius of each blood vessel.

    Arterial system stenosis can cause a variety of cardiovascular diseases. Zhang et al. [15] studied carotid artery stenosis. Pagiatakis et al. [16,17] studied the influence of stenosis eccentricity and stenosis severity on coronary bifurcation lesions. The fractional flow reserve (FFR), a nondimensional index, can combine the physiological effects of the coronary artery with the severity of stenosis to evaluate the hemodynamic significance of the coronary artery disease noninvasively. FFR is referred to as the ratio of the maximum hyperemic blood flow in the coronary artery with stenosis to the maximum hyperemic blood flow that occurs in the same coronary artery under normal conditions [18,19,20]. Clinical research has demonstrated that FFR between 0.75 and 0.8 (that is, under maximum congestion conditions, when the pressure of the distal coronary artery is 75–80% of the proximal) is usually considered to be the threshold for ischemic attacks. Lesions below this threshold are considered to be the cause of ischemia [21,22].

    Our purpose is to establish a coronary artery lumped parameter model that can effectively simulate its blood supply capacity. The aortic valve fluid-structure interaction (AV FSI) simulator can simulate the interaction between the heart and the arterial system and provide an in-depth understanding of the function of the cardiovascular system. The 0D model simplifies the description of the cardiovascular system and provides a series of algebraic equations to simulate changes in blood flow and pressure. In this paper, combining the AV FSI simulator and the 0D model can make full use of the advantages of both to achieve more accurate simulation and prediction. The head loss model is incorporated at the entrance of the coronary artery. The outlet connects with the symmetric structured tree model proposed by Olufsen et al. The resistance-inductance-capacitance (RLC) parameters in the 0D coronary artery model are optimized in a special way. In addition, the coronary artery stenosis model is established by changing the local RLC parameters, and the blood supply capacity of the coronary artery under different degrees of stenosis is analyzed using FFR and relative flow.

    The AV FSI simulator is built using a porcine pericardial valve, as demonstrated in our previous research [23,24]. The thickness of the valve leaflet is 0.04 cm, and the outer wall connecting with it has a thickness of about 0.15 cm and an inner diameter of 1.3 cm. The length of the whole blood vessel is 13 cm. The three-element Windkessel model, frequently used to simulate the afterload [25], is coupled in the outlet. Figure 1 shows the schematic illustration of the AV FSI simulator.

    Figure 1.  A schematic illustration of the AV FSI simulator.

    Aortic valve leaflets are considered to be nonlinear, incompressible, anisotropic, and hyperelastic. Specifically, the strain energy function of aortic valve leaflets [26] is

    W=C10[eC01(I13)1]+k12k2[ek2(I41)21], (2.1)

    in which k1, k2, C10, and C01 are material parameters; k1, k2 are the material constants related to the stress and the dimensionless parameter, respectively; C10, C01 are utilized to characterize the properties of the matrix material. There are two strain invariants, I1=tr(C) and I4=f0(Cf0), along with the reference fiber direction (f0). C is the right Cauchy-Green deformation tensor. These parameter values are presented in Table 1.

    Table 1.  The parameters of the AV FSI simulator, adapted from [23] and [27].
    C10(kPa) C01 k1(kPa) k2 C(ml/mmHg) Rp(mmHgs/ml) Rc (mmHgs/ml)
    1.21 7.99 24.23 57.62 1.75 0.79 0.033

     | Show Table
    DownLoad: CSV

    The Windkessel model [28] is

    CdPWKdt+PWKRp=QAO,PAOPWKRc=QAO, (2.2)

    in which C represents the vascular compliance, PWK represents the pressure stored within the Windkessel model, Rp represents the peripheral resistance, and Rc represents the characteristic resistance. These parameter values are presented in Table 1.

    In the Immersion Boundary/Finite Element (IB/FE) method [29], the fluid domain (Ωf(0)) uses the Eulerian grids and finite difference discretization, while the structure domain (Ωs(0)) employs the Lagrangian grids and finite element discretization [23,30]. The influence of the structure on the fluid is represented by incorporating a force source term into the Navier-Stokes equation [31]. The Dirac delta function δ(x) is utilized to establish the coupling between the structure described in Lagrangian form and the fluid represented in Eulerian form. The IB/FE method can be described as follows [31,32]:

    ρ(u(x,t)t+u(x,t)u(x,t))=p(x,t)+μΔu(x,t)+fe(x,t),u(x,t)=0,fe(x,t) = Ωs(0)Fe(X,t)δ(xχ(X,t))dX,χ(X,t)t = Ωu(x,t)δ(xχ(X,t))dx,Ωs(0)Fe(X,t)V(X)dX = Ωs(0)Pe(X,t):XV(X)dX, (2.3)

    where x=(x1,x2,x3)Ω represents the Euler coordinate system, and Ω represents the physiil region occupied by the aortic valve model. X=(X1,X2,X3)Ωs(0) represents the Lagrange coordinate system. χ(X,t)Ω represents the Eulerian coordinate of the position that a Lagrangian coordinate point X occupies at time t. u(x,t) represents the Eulerian velocity, and p(x,t) represents the Eulerian pressure. fe(x,t) and Fe(X,t) are the Eulerian elastic force and the Lagrangian elastic force, respectively. Pe(X,t) represents the first Piola-Kirchhoff stress tensor, which can be derived from the strain energy function. In general, the density of blood (ρ) and its dynamic viscosity (μ) are assumed to remain constant, e.g., ρ=1.0 g/cm3 and μ=4.0×103 kg/ms. For more details of the AV FSI simulator and the IB/FE method, see [24,29,30].

    The blood regarded as a continuous, incompressible Newtonian fluid, and the blood vessel has a cylindrical wall that is rigid [33,34]. According to the analogy between the hemodynamic system and the circuit system, we can establish an equivalent circuit model of a blood vessel; see Figure 2. The blood flow satisfies the Newton equation and the continuity equation for any minuscule artery length [7]. The governing equations [7,35] for flow rate and pressure of each blood vessel is

    PinPout=LQint+RQin,QinQout=CPoutt, (2.4)

    in which Pin and Pout are the average inlet and outlet pressures that vary with time. Similarly, Qin and Qout represent the average inlet and outlet flow rates changing with time, respectively. R, L, and C are the blood flow viscous resistance, blood flow inertia, and vascular compliance, respectively, analogous to the resistance, inductance, and capacitance in the circuit system.

    Figure 2.  Based on the simplified structure model of arterial vessel (a), the corresponding equivalent circuit model (b) is obtained.

    The coronary artery emerges from the aortic sinus at the root of the aorta and is divided into the left and right branches that contain many vascular branches. Based on the CT of a patient-specific coronary [10] and literature [9], we consider a coronary artery model composed of 20 segments: 17 segments attributed to the left coronary artery (three additional septal branch (SPT) vascular segments are added) and 3 segments pertaining to the right coronary artery. The schematic diagram of the coronary artery simplified structure model is shown in Figure 3, and the 0D model of the whole coronary artery is detailed in Figure 4. Each segment of the coronary artery can be modeled with Eq (2.4). The complete system can be mathematically modeled by the ordinary differential equations (ODEs) as

    yt=Ay+b. (2.5)

    In this system, y denotes the vector consisting of unknown vessel flow rates and pressures, A represents a square matrix composed of parameter values, and b represents a column vector derived from the specified boundary conditions.

    Figure 3.  The simplified structure model of coronary artery, adapted from [10].
    Figure 4.  The whole coronary artery lumped parameter model, adapted from [7]. In the model, the left main coronary artery (LMCA) includes left anterior descending branch (LAD), left circumflex branch (LCX), SPT, diagonal branch (DIAG), and marginal branch (MARG), and the right coronary artery (RCA) includes posterior descending artery (PDA) and posterior lateral artery (PLA). Rh is a variable resistance related to flow rate [10]. Z1-Z13 are the root impedances for the symmetric structured tree models. The ventricular pressure generator (VPG) simulates the feedback pressure in the myocardium.

    Initial conditions, inlet and outlet boundary conditions, bifurcation boundary conditions, and so on are needed to solve the hemodynamics of the coronary artery described by Eq (2.5). In this study, we give the following conditions:

    (1) After several numerical simulations, it is observed that the initial conditions mainly have a significant effect on the speed of reaching the stable solution but have a small effect on the final numerical results of the steady state. According to the Wiggers diagram [36], the initial pressure and flow rate at each segment can be set to 80 mmHg and 0 ml/s.

    (2) In order to obtain the physiologically realistic blood pressure and flow rate in the coronary artery, the aortic root pressure calculated by the AV FSI simulator (see Section 2.1) is selected as the inlet boundary condition.

    (3) At the entrance of LMCA and RCA, the sudden decrease in blood vessel radius will lead to an abrupt decrease in blood pressure. The pressure drop can be described by the following empirical formula:

    Δp=ζρv22, (2.6)

    where ζ=2.88 is the local energy loss coefficient [37] and v is the average velocity. Therefore, the flow resistance Rh can be calculated using the following Eq [10]:

    Rh=ΔpQ=ζρ2A2Q, (2.7)

    in which A represents the cross-sectional area of the blood vessel. It can be seen that Rh is a variable flow resistance proportional to the flow rate.

    (4) Each vessel's terminal is connected to the capillaries downstream, making it challenging to assess the boundary conditions at the outlet. An outflow boundary condition, the structured tree model presented by Olufsen et al. [12], has been used in literature, in which each vessel is modeled as a compliant linear segment, and the root impedance (Z) depends largely on the minimum vessel diameter at the end truncation (rmin). The structured tree has nothing to do with the actual geometry of the vessel but a topological structure based on the literature data [14]. In order to simplify the model, we choose the symmetric structured tree model and set rmin = 0.01 mm, indicating that the structured tree terminates when the vessel radius is less than 0.01 mm. The root impedance calculation formula [13] of the N-generation symmetric structured tree model is

    Z=8μλπr302α3(2α3)N2α31, (2.8)

    in which λ=50 defines the ratio of vessel length to radius [38], r0 represents the root radius of the initial vessel, and the ratio of the radius of the daughter vessel (rd) to the radius of the parent vessel (rp) is α=0.778. The terminal impedance of each coronary artery outlet is shown in Table 2.

    Table 2.  Parameter values for the coronary artery lumped parameter model.
    Segment Results of Eq (2.9) Optimization results Z
    R L C R L C
    LMCA 0.139 0.034 0.001 2.0159 0.01 0.00218
    LAD 0.473 0.069 0.001 0.3685 0.001 0.0012
    LAD1 3.190 0.238 0.001 2.9761 0.2064 0.00163 145.6
    LAD2 0.922 0.100 0.001 0.747 0.0017 0.00107
    SPT1 3.184 0.235 0.001 2.689 0.2031 0.0012 155.6
    SPT2 2.906 0.287 0.002 3.8576 0.0017 0.0001 407.5
    SPT3 2.819 0.281 0.002 3.6832 0.0014 0.0001 327
    LAD3 6.803 0.238 0.0002 3.6025 0.009 0.00018 155
    LAD4 4.775 0.353 0.001 3.731 0.1658 0.00188 80.2
    DIAG 1.941 0.184 0.001 2.0509 0.094 0.00194 55.1
    LCX 0.336 0.050 0.001 0.4362 0.001 0.00104
    MARG1 1.477 0.140 0.001 1.5214 0.184 0.00146 155
    LCX1 0.401 0.048 0.001 0.2699 0.0264 0.00075
    MARG2 2.292 0.161 0.001 2.014 0.150 0.001 168.4
    LCX2 0.856 0.082 0.001 0.9891 0.0082 0.00049
    MARG3 2.463 0.186 0.001 2.7604 0.0794 0.00083 72
    LCX3 2.807 0.203 0.001 3.7397 0.1812 0.00095 170.1
    RCA 1.602 0.257 0.049 2.5284 0.5 0.02061
    PLA 1.195 0.138 0.013 5.0 0.1896 0.01017 59.81
    PDA 2.078 0.170 0.008 4.6671 0.001 0.00913 139.72
    F 1.228 0.345

     | Show Table
    DownLoad: CSV

    (5) Blood flow in the coronary artery is distinct from that in other parts of the arterial system, as it is significantly affected by the relaxation and contraction of the heart. In order to simulate coronary blood flow veritably, it is necessary to include the cardiac pressure acting on the coronary artery vessel (i.e., ventricular pressure) [39]. At the end of the model, the VPG is connected to simulate the intramyocardial pressure, which selected left ventricular pressure according to the location of the coronary artery.

    (6) The principles of flow conservation and pressure continuity are ensured at the points where blood vessels bifurcate. Specifically, the blood flow in the parent vessel is equivalent to the aggregate sum of the flows in all its daughter vessels, and the blood pressure of the parent vessel is equal to the blood pressure of each daughter vessel.

    In this paper, a fourth-order Runge-Kutta method [40] is employed for the solution of the ODEs describing the hemodynamics of the coronary artery. In order to obtain stable numerical results, each simulation performs three cycles.

    Generally, there are two methods to determine the RLC parameter values in the 0D coronary artery model. One is the manual estimate method, which depends on experimental measurements, empirical data, or iterative adjustment of parameters. This method takes a lot of time. The other is calculated through the lumped parameter modeling formula, but the formula is not unique, there is [7],

    R=128μl/πD4,L=4ρl/πD2,C=πD3l/4Gh, (2.9)

    in which the Young's modulus of blood vessel G=2×105 Pa [41], l is the length of the blood vessel, D represents the diameter of the blood vessel, and the thickness of the blood vessel wall h=0.08D [42]. These parameter values related to the anatomical structure of the patient's coronary artery can be obtained through a series of medical examinations, such as coronary angiography and echocardiography.

    In order to obtain realistic coronary artery blood pressure and flow rate, the value range of RLC parameters in the 0D coronary artery model is preliminarily estimated by Eq (2.9), and then optimized by the intelligent optimization algorithm, common artificial intelligence [43], genetic algorithm [44], and particle swarm optimization (PSO) algorithm [45]. The PSO algorithm was first developed by Russell Eberhart and James Kennedy in 1995, inspired by the foraging mode of birds (the birds allow groups to find optimal purposes through information sharing). In the PSO algorithm, each optimization solution is simulated as a 'particle'. Each particle navigates through the search space, continually updating its position and velocity information to iteratively converge toward the optimal solution. The algorithm is simple, easy to understand, and easy to program. It has fast convergence speed and can find the approximate optimal solution to the problem in a short time. The number of parameters is small, and the adjustment of parameters is relatively simple [46]. It is suitable for solving high-dimensional optimization problems. The summary of the PSO algorithm can be seen in Appendix A.

    The optimization objective function is defined as the root mean square error between the reference solution (yn) and the numerical result (yn), expressed as:

    F=1KKn=1|ynyn|2, (2.10)

    in which K represents the number of sample points. The optimal results are obtained when the objective function attains its minimum value. The termination condition is to reach the maximum number of iterations or obtain an acceptable satisfactory solution.

    Coronary stenosis leads to diminished blood flow reaching the heart, causing myocardial ischemia, which in turn leads to various cardiovascular diseases. It is necessary to research the influence of coronary artery stenosis on local blood flow. To establish the lumped parameter model of coronary artery stenosis, we adjust the RLC parameter values of the corresponding vascular segments in our 0D coronary artery model. The reduction rate of the vascular cross-sectional area (β) is defined as the ratio between the stenosed vessel's cross-sectional area (As) and the healthy cross-sectional area (A0), namely,

    β=AsA0. (2.11)

    According to the Eq (2.9), the relationship [47] between the RLC parameter values in the stenosis state and the healthy state can be derived as follows:

    Rs=R0β2 ,Ls=L0β1 ,Cs=C0β3/2, (2.12)

    where Rs, Ls, and Cs represent the resistance, inductance, and capacitance in the stenosis state, respectively. R0, L0, and C0 represent the ones in the healthy state. The value (1β) represents the vascular stenosis rate, usually expressed as a percentage. It is assumed that other parameters remain constant in the stenosis state.

    FFR quantifies the percentage of normal, hyperemic myocardial blood flow retained when the supplying artery has stenosis [17], has become an important clinical indicator for judging the degree of arterial stenosis that prevents blood flow from reaching the myocardium [20], and is also the gold standard for evaluating whether coronary artery stenosis causes ischemia [48]. In clinical practice, FFR can be expressed as the ratio between the time-averaged proximal pressure Pp and the time-averaged distal pressure Pd in coronary artery stenosis (under conditions of maximal myocardial hyperemia) [20], i.e.,

    FFR=PdPp. (2.13)

    For the flow rates predicted by the 0D model, the average flow is calculated by integrating the flow rate waveform over one heartbeat. This average flow is subsequently normalized by dividing it with the average flow rate that would occur in the absence of stenosis, thus obtaining the relative flow [7]. Relative flow represents the change of flow under different stenosis degrees based on normal conditions. It can be employed to compare the blood flow under different conditions and is an effective index for gauging the coronary artery's ability to supply blood to the myocardium.

    The cross-sectional area of the coronary artery is decreased by 20%, 40%, 60%, and 80% to simulate different degrees of stenosis. The pressure-related FFR and flow rate-related relative flow are analyzed to study the blood supply capacity under coronary artery stenosis.

    In all numerical simulations, the heart rate is set to 75 beats/min, so one cardiac cycle is 60/75 = 0.8 s. Using the IBAMR open source software platform (http://github.com/IBAMR/IBAMR) to implement the AV FSI simulator introduced in Section 2.1, the aortic root pressure and the flow rate through the aortic valve are in line with the real physiological conditions, as shown in Figure 5. The systolic/diastolic blood pressure of the aorta is maintained within the normal physiological range (120/80 mmHg). With the aortic valve opening, blood flows into the aorta, and the flow rate increases to a peak of 572.43 ml/s. When the aortic valve closes, regurgitation occurrs and the backflow rate can reach up to 170.70 ml/s. Over the course of a single cardiac cycle, the volume of blood that flows into the aorta is 89.73 ml. The numerical results align well with the experimental data reported in articles [49] and [50].

    Figure 5.  Aortic and left ventricular pressure (a), aortic valve flow rate (b) in a cardiac cycle.

    It should be noted that there are oscillations in the pressure and flow rate curves in Figure 5, which are mainly caused by numerical effects and also involve physiological effects. The oscillations occur at about 0.4 s. This period is the stage when the aortic valve is about to close completely, and the hemodynamic behavior is particularly complex. The AV FSI simulator cannot accurately simulate this phenomenon.

    In order to obtain the parameters that can better describe the properties, the RLC parameters in the coronary artery lumped parameter model are optimized as introduced in Section 2.3. For the right coronary artery model, the numerical results of Duanmu et al. [10] (the flow rate and pressure of RCA and PLA) are selected as the reference solution, and their left ventricular pressure and aortic pressure are used as the boundary conditions to complete the parameter optimization. The comparison results are shown in Figure 6. The black circular data are the results of Duanmu et al. [10], and the red triangle data are numerical results after parameter optimization. It can be seen that the numerical results of the flow rate and pressure after parameter optimization are in good agreement with the reference solution. The fitness values of pressure and flow rate waveforms are calculated separately. Specifically, the pressure-related fitness value is 0.195, while the flow rate-related fitness value is 0.239. Therefore, the optimization result of the pressure waveform seems to be more ideal than that of the flow rate waveform.

    Figure 6.  Comparison of right coronary artery flow rate and pressure (RCA and PLA).

    For the left coronary artery model, the numerical results of Pietrabissa et al. [9] (the flow rate of LMCA, LAD, and SPT2) are selected as the reference solution, and their left ventricular pressure and aortic pressure are used as the boundary conditions to obtain the parameter optimization values. Figure 7 presents the comparison of the flow rate numerical results with those obtained by Pietrabissa et al. [9]. It can be seen that the numerical results of the flow rate after parameter optimization are similar to the waveforms of the reference solution. After parameter optimization, the objective function values of the left and right coronary arteries calculated by Eq (2.10) are 0.526 and 0.209, respectively. This shows that optimizing the right coronary artery model is easier than optimizing the left coronary artery model. On the one hand, the complexity of the structure of the left and right coronary arteries is different; the left coronary artery has more branches, its structure is more complex, and the parameter optimization is more difficult. On the other hand, the choice of the reference solution is different. The left coronary artery lacks the pressure reference solution, and the number of reference solutions is relatively scarce. To address this issue, we employed two strategies. First, the RLC parameters with an ideal optimization effect are fixed, reducing the number of variables while dynamically adjusting them based on optimization outcomes. Second, the numerical results with an ideal optimization effect can be approximated as a reference solution, increasing the available data. Table 2 lists the optimized parameter values of the 0D coronary artery model.

    Figure 7.  Comparison of left coronary artery flow rate (LMCA, LAD, and SPT2).

    We first set the aortic root pressure calculated by the AV FSI simulator as the inlet boundary condition, and the corresponding flow rate of each segment of the left coronary artery is shown in Figure 8. The blood flow dynamics in the coronary artery differ from those in other sections of the arterial system. During the systolic period (the time is 0–0.39 s), the left ventricle contracts, the intramyocardial pressure increases, the extravascular pressure of the left coronary artery increases, and the blood flow is restricted, resulting in a decrease in the flow rate, which is not synchronized with aortic flow. Then, it increases with the decrease in intramyocardial pressure until entering the diastolic period (the time is 0.39–0.8 s). The trend is well-captured in Figure 8. The relative position of the coronary artery in the heart also has an effect on the flow rate. For example, the flow rate of the epicardial coronary artery (e.g., LMCA, LAD, and LCX) in the systolic period is not as pronounced as that of the endocardial coronary artery (e.g., SPT2). Specifically, the flow rate of the endocardial coronary artery decreases rapidly at the beginning of the systolic period and maintains a low level for a period of time. The flow rate of the epicardial coronary artery decreases slowly during systole, and there is a slight increase. During diastole, the flow rate of the endocardial coronary artery and epicardial coronary artery increases rapidly.

    Figure 8.  Flow rate curve of each segment of the left coronary artery.

    Figure 9(a) depicts the flow rate waveforms at various spatial locations along the main trunk, starting from the entrance of the LMCA trunk (LMCA, LAD, LAD2, and LAD4). It can be observed that the magnitude of the flow rate wave decreases significantly, but the waveforms remain basically unchanged. The flow rate waveforms in the LCX trunk (LCX, LCX1, LCX2, LCX3) shown in Figure 9(b) also present the same trend. In contrast, the pressure waveforms change little at different positions in the same trunk, as illustrated in Figure 9(c).

    Figure 9.  Numerical results of LMCA trunk and LCX trunk.

    Figure 10 shows the flow rate and pressure of each segment of the right coronary artery. In the systolic period, the flow rate in the right coronary artery does not significantly decrease compared with the left coronary artery. This is mainly due to the lower pressure environment in which the right ventricle functions, with its peak pressure being only about one-fifth of the mean aortic pressure [51]. Based on the numerical results obtained from the 0D coronary artery model, it is estimated that during a single cardiac cycle, the blood flow entering the coronary artery distributes as follows: 3.96 ml into the overall coronary artery, with 3.11 ml flowing to the left coronary artery and 0.85 ml directed toward the right coronary artery. In total, the coronary artery blood flow accounted for about 4.41% of the cardiac output.

    Figure 10.  Flow rate and pressure of each segment of the right coronary artery.

    Coronary artery stenosis is very important for the study of coronary artery blood supply capacity, which leads to a decrease in the cross-sectional area of blood vessels, hindering blood flow. In severe cases, it can cause myocardial ischemia. The coronary artery stenosis model is established by changing the RLC parameters of a specific segment; see Eqs (2.11) and (2.12). In these simulations, we observe the changes in flow rate and pressure at different degrees of stenosis (20%, 40%, 60%, and 80%) to study how the stenosis rate affects flow rate and pressure. Figure 11 illustrates the relationship between flow rate, pressure, and stenosis rate in the PDA of the right coronary artery. We find that when the stenosis rate rises, the flow rate and the pressure decrease, and when the stenosis rate is higher than 60%, the blood pressure decreases significantly.

    Figure 11.  The flow rate (a) and pressure (b) of PDA under different stenosis rates.

    In the left coronary artery, six typical blood vessels (LAD, LAD2, SPT1, LCX, LCX2, and MARG2) are selected for stenosis simulation according to different locations (different branches and different generations). The relationship between flow rate, pressure, and stenosis rate of the left coronary artery is similar to that of the right coronary artery. Figure 12 displays the variations in flow rate and pressure across six different segments at varying stenosis rates. The relationship between flow rate, pressure, and stenosis rate in the left coronary artery is analogous to that in the right coronary artery. Compared with flow rate, pressure decreases gradually with an increase in stenosis rate, which indicates that blood pressure is a more stable diagnostic criterion than flow rate in clinical medicine.

    Figure 12.  The flow rate and pressure of some segments of the left coronary artery under different stenosis rates.

    According to the flow rate and pressure at different stenosis rates, the relative flow and FFR can be calculated. Figure 13(a) shows the relative flow of the right coronary artery PDA when PDA is stenosed, and its FFR is shown in Figure 13(b), which is calculated by Eq (2.13). The FFR of PDA in healthy state is 0.98. With the increase in stenosis rate, the relative flow and FFR both decrease. The rate of change in relative flow and FFR increases significantly when the stenosis rate exceeds 60%, consistent with Figure 11. When the stenosis rate is higher than 72%, the FFR is lower than 0.8, reaching the threshold of myocardial ischemia.

    Figure 13.  The relative flow (a) and FFR (b) of PDA.

    In the left coronary artery, the relative flow and FFR of LMCA trunk (LAD, LAD2, and LAD4), LCX trunk (LCX1, LCX2, and LCX3), and MARGs at different locations (MARG1, MARG2, and MARG3) are calculated and shown in Figure 14. Similar to PDA, the relative flow and FFR of the left coronary artery decrease with increased stenosis. The FFR of each vessel segment of the LMCA trunk in the healthy state is 0.990, 0.989, and 0.967, respectively. The FFR of the LCX trunk and MARGs are also above 0.97. When the stenosis rate is between 50% and 60%, the relative flow and FFR begin to decrease significantly. In addition, for the LMCA trunk, LCX trunk, and MARGs, the decrease in relative flow and FFR becomes more and more significant with the increase in generation.

    Figure 14.  The relative flow and FFR of some segments of the left coronary artery.

    In a cardiac cycle, the movement of the aortic valve occurs mainly in the systolic period. The strong contraction of the ventricular myocardium leads to a dramatic rise in the intraventricular pressure. Once the left ventricular pressure surpasses the aortic pressure, the aortic valve opens, and the blood starts to flow into the aorta. When the left ventricular pressure reaches its peak, the aortic valve is completely opened, and blood rushes into the aorta. The maximum flow rate was 572.43 ml/s. As the left ventricular pressure decreases below that of the aorta, the aortic valve begins to close. The AV FSI simulator can well simulate the hemodynamics of the aortic valve in a cardiac cycle and is consistent with clinical measurement data and published literature results [49,50]. Therefore, the AV FSI simulator can provide the boundary condition in a healthy state for the 0D model to achieve a more realistic simulation. In the future, we plan to introduce more boundary conditions in different states to simulate the performance of the cardiovascular system in various complex situations.

    Taking the results of the coronary artery published by Pietrabissa et al. [9] and Duanmu et al. [10] as the reference solution, the PSO algorithm is utilized to optimize the RLC parameters in the 0D model. Because the right coronary artery has a simple structure (including only three vessels) and gives more reference solution, its optimization outcome has superior performance compared to that of the left coronary artery. Table 2 presents the RLC parameter values for each segment of the coronary artery. These optimized parameters are then implemented in subsequent studies.

    The 0D coronary artery model established in this paper can simulate the flow rate and pressure of all vessel segments in a cardiac cycle and evaluate the blood supply capacity of the coronary artery. The predictions align well with previous research findings and clinical experimental data. However, there is oscillation in the numerical result, which may be due to the influence of discretization error, mesh accuracy, boundary condition processing, and other factors. The flow rate in the coronary artery is highly correlated with the heart and arterial system. During the systolic period, the driving force for perfusing the left myocardial is the difference between the intramyocardial pressure and the aortic pressure, and the flow rate decreases with myocardial contraction [10,17,51] and then increases with the reduction of intramyocardial pressure until the beginning of diastole. However, the flow rate in the right coronary artery does not significantly decrease during the systolic period because the lower pressure environment within the right ventricle is compared to the left ventricle [9,39]. In addition, we also studied the flow rate and pressure waveforms along the different spatial positions of the trunk. Since the different cross-sectional areas of vessels along the trunk, the flow rate shows different amplitudes. Concretely, the flow rate waveforms tend to scale (normalized relative to their average) into a curve except for a small phase difference, and the pressure waveforms change little. Huo et al. [52] discovered that the attenuation of the pressure waveforms' amplitude is particularly pronounced in blood vessels with a diameter less than 10 μm. Under normal circumstances, the average coronary blood flow constitutes 4% to 5% of the total cardiac output [51], and our numerical result is 4.41%.

    To investigate the impact of coronary artery stenosis on both flow rate and pressure dynamics and evaluate the blood supply capacity of the stenotic coronary artery, the stenosis model is established to simulate the flow rate and pressure of each segment of the coronary artery when the stenosis rate varied between 20% and 80%. As the stenosis rate increases, the flow rate and pressure decrease to varying degrees. In clinical medicine, pressure is a very stable diagnostic criterion because the reduction of pressure is more regular [10]. The flow rate-related relative flow and pressure-related FFR show good concurrence with the previously published experimental data [7,53]. With the increase in stenosis rate, the relative flow and FFR have the same trend, which is decreasing. Hamilos et al. [54] showed that the FFR value of the stenosis rate between 60% and 70% is between 0.68 and 0.97. Koo et al. [55] found that when the stenosis rate is 79 ± 11%, the average FFR is 0.81 ± 0.12, and only one-third of lesions exhibit an FFR value below 0.75. For intermediate blood vessel segments, the FFR in healthy status is not less than 0.97, while for the terminal vessels, the range is between 0.9 and 0.97 [16].

    However, the coronary artery model still has some limitations. This model ignored the interplay between the heart, the arterial system, and the coronary artery, and it could not predict wave propagation. We assume that the RLC parameters are always constant, which is approximately true in the diastolic period, but in the systolic period, the parameters will change due to ventricular contraction [7]. For a solitary lesion, the occurrence of ischemia is affected by many factors, such as the shape and location of the lesion, the size and function of the myocardium being perfused, and the degree of collateral circulation [56,57,58]. These factors are not taken into account in the coronary artery stenosis model.

    In this paper, we establish a complete coronary artery lumped parameter model based on the similarity between the circuit system and the cardiovascular system, including the left and right coronary arteries. The model uses RLC parameters to describe each vessel segment and introduces the AV FSI simulator to calculate the inlet boundary condition. It can predict the pressure and flow rate distribution throughout the coronary artery tree up to the terminal blood vessel and evaluate the blood supply capacity under health status and different degrees of stenosis. The relative flow and FFR align well with the published literature. The model can effectively simulate the common characteristics of coronary artery hemodynamics and be easily transformed into clinical practice for real-time clinical diagnosis.

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

    This work is supported by National Natural Science Foundation of China (Grant No. 12271440), which are gratefully acknowledged.

    The authors declare there is no conflict of interest.

    The PSO algorithm is programmed on the MATLAB platform without relying on any ready-made libraries or toolkits. It is constructed according to the core idea of the PSO algorithm. Its pseudocode is as follows:

    Algorithm 1 PSO algorithm.
    for each particle i do
     Initialize velocity vi and position xi for particle i
     Evaluate particle i and set optimal position of individual history pBesti=xi
    end for
    Optimal position of group history gBest=min{pBesti}
    while no stop do
     for i=1 to N do
      Update the velocity and position of particle i using Eq (9.1)
      Evaluate particle i
      if Fmin(xi)<Fmin(pBesti) then
       pBesti=xi
      end if
      if Fmin(pBesti)<Fmin(gBest) then
       gBest=pBesti
      end if
     end for
    end while

    vk+1id=ωvkid+c1rand1(pBestkidxkid)+c2rand2(gBestkdxkid),xk+1id=xkid+vkid, (A.1)

    where N is the particle swarm size set to 200. k is the current number of iterations and cannot exceed the maximum number of iterations of 400. d represents the particle dimension, that is, the number of independent variables, and there are 60 RLC parameters to be optimized in this paper. The inertia weight ω = 0.7, the individual learning factor (c1), and group learning factor (c2) are set to 0.5. rand1 and rand2 are two random numbers over [0, 1] intervals.



    [1] Srivastava A, Chandra N, Kumar S (2019) The role of biofilms in medical devices and implants, Biofilms in Human Diseases: Treatment and Control, Springer, Cham, 151-165.
    [2] Sun F, Qu F, Ling Y, et al. (2013) Biofilm-associated infections: Antibiotic resistance and novel therapeutic strategies. Future Microbiol 8: 877-886. doi: 10.2217/fmb.13.58
    [3] Khatoon Z, McTiernan CD, Suuronen EJ, et al. (2018) Bacterial biofilm formation on implantable devices and approaches to its treatment and prevention. Heliyon 4: e01067. doi: 10.1016/j.heliyon.2018.e01067
    [4] Percival SL, Suleman L, Vuotto C, et al. (2015) Healthcare-associated infections, medical devices and biofilms: Risk, tolerance and control. J Med Microbiol 64: 323-334. doi: 10.1099/jmm.0.000032
    [5] Van Epps JS, Younger JG (2016) Implantable device-related infection. Shock 46: 597-608. doi: 10.1097/SHK.0000000000000692
    [6] Tunney MM, Gorman SP, Patrick S (2002) Infection associated with medical devices. Int J Gen Syst 31: 195-205.
    [7] Donlan RM (2001) Biofilms and device-associated infections. Emerg Infect Dis 7: 277-281. doi: 10.3201/eid0702.010226
    [8] Vertes A, Hitchins V, Phillips KS (2012) Analytical challenges of microbial biofilms on medical devices. Anal Chem 84: 3858-3866. doi: 10.1021/ac2029997
    [9] Veerachamy S, Yarlagadda T, Manivasagam G, et al. (2014) Bacterial adherence and biofilm formation on medical implants: A review. P I Mech Eng H 228: 1083-1099. doi: 10.1177/0954411914556137
    [10] Ramstedt M, Ribeiro IAC, Bujdakova H, et al. (2019) Evaluating efficacy of antimicrobial and antifouling materials for urinary tract medical devices: Challenges and recommendations. Macromol Biosci 19: e1800384. doi: 10.1002/mabi.201800384
    [11] Reid G (1999) Biofilms in infectious disease and on medical devices. Int J Antimicrob Ag 11: 223-226. doi: 10.1016/S0924-8579(99)00020-5
    [12] Bryers JD (2008) Medical biofilms. Biotechnol Bioeng 100: 1-18. doi: 10.1002/bit.21838
    [13] Darouiche RO (2004) Treatment of infections associated with surgical implants. N Engl J Med 350: 1422-1429. doi: 10.1056/NEJMra035415
    [14] Darouiche RO (2001) Device-associated infections: A macroproblem that starts with microadherence. Clin Infect Dis 33: 1567-1572. doi: 10.1086/323130
    [15] Campoccia D, Montanaro L, Arciola CR (2006) The significance of infection related to orthopedic devices and issues of antibiotic resistance. Biomaterials 27: 2331-2339. doi: 10.1016/j.biomaterials.2005.11.044
    [16] Azevedo AS, Almeida C, Melo LF, et al. (2017) Impact of polymicrobial biofilms in catheter-associated urinary tract infections. Crit Rev Microbiol 43: 423-439. doi: 10.1080/1040841X.2016.1240656
    [17] World Health Organization, Report on the Burden of Endemic Health Care-Associated Infection Worldwide, 2011. Available from: https://apps.who.int/iris/bitstream/handle/10665/80135/9789241501507_eng.pdf;jsessionid=C2A2CD7A9910C3B4AC1FC3B18F1DBD11?sequence=1.
    [18] Feneley RCL, Hopley IB, Wells PNT (2015) Urinary catheters: History, current status, adverse events and research agenda. J Med Eng Technol 39: 459-470. doi: 10.3109/03091902.2015.1085600
    [19] Flemming HC, Wingender J (2010) The biofilm matrix. Nat Rev Microbiol 8: 623-633. doi: 10.1038/nrmicro2415
    [20] Rabin N, Zheng Y, Opoku-Temeng C, et al. (2015) Biofilm formation mechanisms and targets for developing antibiofilm agents. Future Med Chem 7: 493-512. doi: 10.4155/fmc.15.6
    [21] Deva AK, Adams WP, Vickery K (2013) The role of bacterial biofilms in device-associated infection. Plast Reconstr Surg 132: 1319-1328. doi: 10.1097/PRS.0b013e3182a3c105
    [22] Garrett TR, Bhakoo M, Zhang Z (2008) Bacterial adhesion and biofilms on surfaces. Prog Nat Sci 18: 1049-1056. doi: 10.1016/j.pnsc.2008.04.001
    [23] Lindsay D, von Holy A (2006) Bacterial biofilms within the clinical setting: What healthcare professionals should know. J Hosp Infect 64: 313-325. doi: 10.1016/j.jhin.2006.06.028
    [24] Corte L, Pierantoni DC, Tascini C, et al. (2019) Biofilm specific activity: A measure to quantify microbial biofilm. Microorganisms 7: 73. doi: 10.3390/microorganisms7030073
    [25] Marić S, Vraneš J (2007) Characteristics and significance of microbial biofilm formation. Period Biol 109: 115-121.
    [26] Donlan RM (2002) Biofilms: Microbial life on surfaces. Emerg Infect Dis 8: 881-890. doi: 10.3201/eid0809.020063
    [27] Vasudevan R (2014) Biofilms: Microbial cities of scientific significance. J Microbiol Exp 1: 84-98.
    [28] Floyd KA, Eberly AR, Hadjifrangiskou M (2017) Adhesion of bacteria to surfaces and biofilm formation on medical devices, Biofilms and Implantable Medical Devices: Infection and Control, Woodhead Publishing, 47-95.
    [29] Chen M, Yu Q, Sun H (2013) Novel strategies for the prevention and treatment of biofilm related infections. Int J Mol Sci 14: 18488-18501. doi: 10.3390/ijms140918488
    [30] Zhu Z, Wang Z, Li S, et al. (2019) Antimicrobial strategies for urinary catheters. J Biomed Mater Res A 107: 445-467. doi: 10.1002/jbm.a.36561
    [31] Dwivedi P, Narvi SS, Tewari RP (2013) Application of polymer nanocomposites in the nanomedicine landscape: Envisaging strategies to combat implant associated infections. J Appl Biomater Funct Mater 11: 129-142.
    [32] Singha P, Locklin J, Handa H (2017) A review of the recent advances in antimicrobial coatings for urinary catheters. Acta Biomater 50: 20-40. doi: 10.1016/j.actbio.2016.11.070
    [33] Prabhurajeshwar C, Chandrakanth RK (2017) Probiotic potential of lactobacilli with antagonistic activity against pathogenic strains: An in vitro validation for the production of inhibitory substances. Biomed J 40: 270-283. doi: 10.1016/j.bj.2017.06.008
    [34] de Melo Pereira GV, de Oliveira Coelho B, Magalhães Júnior AI, et al. (2018) How to select a probiotic? A review and update of methods and criteria. Biotechnol Adv 36: 2060-2076. doi: 10.1016/j.biotechadv.2018.09.003
    [35] Radaic A, Ye C, Parks B, et al. (2020) Modulation of pathogenic oral biofilms towards health with nisin probiotic. J Oral Microbiol 12: 1809302. doi: 10.1080/20002297.2020.1809302
    [36] FAO (2006) Probiotics in food: Health and nutritional properties and guidelines for evaluation, 2006. Available from: http://www.fao.org/3/a-a0512e.pdf.
    [37] Carvalho FM, Teixeira-Santos R, Mergulhão FJM, et al. (2021) The use of probiotics to fight biofilms in medical devices: A systematic review and meta-analysis. Microorganisms 9: 1-26.
    [38] Barzegari A, Kheyrolahzadeh K, Mahdi S, et al. (2020) The battle of probiotics and their derivatives against biofilms. Infect Drug Resist 13: 659-672. doi: 10.2147/IDR.S232982
    [39] Reid G (1999) The scientific basis for probiotic strains of Lactobacillus. Appl Environ Microb 65: 3763-3766. doi: 10.1128/AEM.65.9.3763-3766.1999
    [40] Markowiak P, Ślizewska K (2017) Effects of probiotics, prebiotics, and synbiotics on human health. Nutrients 9: 1021. doi: 10.3390/nu9091021
    [41] Suez J, Zmora N, Segal E, et al. (2019) The pros, cons, and many unknowns of probiotics. Nat Med 25: 716-729. doi: 10.1038/s41591-019-0439-x
    [42] Salas-Jara M, Ilabaca A, Vega M, et al. (2016) Biofilm forming Lactobacillus: New challenges for the development of probiotics. Microorganisms 4: 35. doi: 10.3390/microorganisms4030035
    [43] Saxelin M, Tynkkynen S, Mattila-Sandholm T, et al. (2005) Probiotic and other functional microbes: From markets to mechanisms. Curr Opin Biotech 16: 204-211. doi: 10.1016/j.copbio.2005.02.003
    [44] Fioramonti J, Theodorou V, Bueno L (2003) Probiotics: What are they? What are their effects on gut physiology? Best Pract Res Cl Ga 17: 711-724.
    [45] Gogineni VK, Morrow LE (2013) Probiotics: Mechanisms of action and clinical applications. J Prob Health 1: 1-11. doi: 10.4172/2329-8901.1000101
    [46] Hemaiswarya S, Raja R, Ravikumar R, et al. (2013) Mechanism of action of probiotics. Braz Arch Biol Techn 56: 113-119. doi: 10.1590/S1516-89132013000100015
    [47] Khalighi A, Behdani R, Kouhestani S (2016) Probiotics: A comprehensive review of their classification, mode of action and role in human nutrition, In: Rao V, Rao L, Probiotics and Prebiotics in Human Nutrition and Health, Croatia: InTech, 10: 63646.
    [48] Moher D, Liberati A, Tetzlaff J, et al. (2009) Preferred reporting items for systematic reviews and meta-analyses: The PRISMA statement. Plos Med 6: e1000097. doi: 10.1371/journal.pmed.1000097
    [49] Nemcová R (1997) Criteria for selection of lactobacilli for probiotic use. Vet Med 42: 19-27.
    [50] Kim AR, Ahn KB, Yun CH, et al. (2019) Lactobacillus plantarum lipoteichoic acid inhibits oral multispecies biofilm. J Endod 45: 310-315. doi: 10.1016/j.joen.2018.12.007
    [51] Bermudez-Brito M, Plaza-Díaz J, Muñoz-Quezada S, et al. (2012) Probiotic mechanisms of action. Ann Nutr Metab 61: 160-174. doi: 10.1159/000342079
    [52] Muñoz M, Mosquera A, Alméciga-Díaz CJ, et al. (2012) Fructooligosaccharides metabolism and effect on bacteriocin production in Lactobacillus strains isolated from ensiled corn and molasses. Anaerobe 18: 321-330. doi: 10.1016/j.anaerobe.2012.01.007
    [53] Kechagia M, Basoulis D, Konstantopoulou S, et al. (2013) Health benefits of probiotics: A review. ISRN Nutr 2013: 481651. doi: 10.5402/2013/481651
    [54] Holzapfel WH, Haberer P, Geisen R, et al. (2001) Taxonomy and important features of probiotic microorganisms in food and nutrition. Am J Clin Nutr 73: 365-373. doi: 10.1093/ajcn/73.2.365s
    [55] Carr FJ, Chill D, Maida N (2002) The lactic acid bacteria: A literature survey. Crit Rev Microbiol 28: 281-370. doi: 10.1080/1040-840291046759
    [56] Leroy F, De Vuyst L (2004) Lactic acid bacteria as functional starter cultures for the food fermentation industry. Trends Food Sci Tech 15: 67-78. doi: 10.1016/j.tifs.2003.09.004
    [57] Ng SC, Hart AL, Kamm MA, et al. (2009) Mechanisms of action of probiotics: Recent advances. Inflamm Bowel Dis 15: 300-310. doi: 10.1002/ibd.20602
    [58] Chen Q, Zhu Z, Wang J, et al. (2017) Probiotic E. coli Nissle 1917 biofilms on silicone substrates for bacterial interference against pathogen colonization. Acta Biomater 50: 353-360.
    [59] Maldonado NC, Silva De Ruiz C, Cecilia M, et al. (2007) A simple technique to detect Klebsiella biofilm-forming-strains. Inhibitory potential of Lactobacillus fermentum CRL 1058 whole cells and products, In: Méndez-Vilas A, Communicating Current Research and Educational Topics and Trends in Applied Microbiology, Formatex, 52-59.
    [60] McMillan A, Dell M, Zellar MP, et al. (2011) Disruption of urogenital biofilms by lactobacilli. Colloid Surface B 86: 58-64. doi: 10.1016/j.colsurfb.2011.03.016
    [61] Velraeds MMC, van Belt-Gritter B De, Busscher HJ, et al. (2000) Inhibition of uropathogenic biofilm growth on silicone rubber in human urine by lactobacilli—A teleologic approach. World J Urol 18: 422-426. doi: 10.1007/PL00007084
    [62] Reid G, Tieszer C (1994) Use of lactobacilli to reduce the adhesion of Staphylococcus aureus to catheters. Int Biodeter Biodegr 34: 73-83. doi: 10.1016/0964-8305(95)00011-9
    [63] Millsap K, Reid G, van der Mei HC, et al. (1994) Displacement of Enterococcus faecalis from hydrophobic and hydrophilic substrata by Lactobacillus and Streptococcus spp. as studied in a parallel plate flow chamber. Appl Environ Microb 60: 1867-1874.
    [64] Ifeoma ME, Jennifer UK (2016) Inhibition of biofilms on urinary catheters using immobilized Lactobacillus cells. African J Microbiol Res 10: 920-929. doi: 10.5897/AJMR2016.8056
    [65] Derakhshandeh S, Shahrokhi N, Khalaj V, et al. (2020) Surface display of uropathogenic Escherichia coli FimH in Lactococcus lactis: In vitro characterization of recombinant bacteria and its protectivity in animal model. Microb Pathogensis 141: 103974. doi: 10.1016/j.micpath.2020.103974
    [66] Busscher HJ, van Hoogmoed CG, Geertsema-Doornbusch GI, et al. (1997) Streptococcus thermophilus and its biosurfactants inhibit adhesion by Candida spp. on silicone rubber. Appl Environ Microb 63: 3810-3817.
    [67] van der Mei HC, van de Belt-Gritter B, van Weissenbruch R, et al. (1999) Effect of consumption of dairy products with probiotic bacteria on biofilm formation on silicone rubber implant surfaces in an artificial throat. Food Bioprod Process 77: 156-158. doi: 10.1205/096030899532303
    [68] van der Mei HC, Free RH, Elving GJ, et al. (2000) Effect of probiotic bacteria on prevalence of yeasts in oropharyngeal biofilms on silicone rubber voice prostheses in vitro. J Med Microbiol 49: 713-718. doi: 10.1099/0022-1317-49-8-713
    [69] van der Mei HC, Buijssen KJDA, van der Laan BFAM, et al. (2014) Voice prosthetic biofilm formation and Candida morphogenic conversions in absence and presence of different bacterial strains and species on silicone-rubber. Plos One 9: e104508. doi: 10.1371/journal.pone.0104508
    [70] Jiang Q, Stamatova I, Kainulainen V, et al. (2016) Interactions between Lactobacillus rhamnosus GG and oral micro-organisms in an in vitro biofilm model. BMC Microbiol 16: 149. doi: 10.1186/s12866-016-0759-7
    [71] Kang MS, Oh JS, Lee HC, et al. (2011) Inhibitory effect of Lactobacillus reuteri on periodontopathic and cariogenic bacteria. J Microbiol 49: 193-199. doi: 10.1007/s12275-011-0252-9
    [72] Teanpaisan R, Piwat S, Dahlén G (2011) Inhibitory effect of oral Lactobacillus against oral pathogens. Lett Appl Microbiol 53: 452-459. doi: 10.1111/j.1472-765X.2011.03132.x
    [73] Marttinen AM, Haukioja AL, Keskin M, et al. (2013) Effects of Lactobacillus reuteri PTA 5289 and L. paracasei DSMZ16671 on the adhesion and biofilm formation of Streptococcus mutans. Curr Microbiol 67: 193-199.
    [74] Söderling EM, Marttinen AM, Haukioja AL (2011) Probiotic lactobacilli interfere with Streptococcus mutans biofilm formation in vitro. Curr Microbiol 62: 618-622. doi: 10.1007/s00284-010-9752-9
    [75] Lee SH, Kim YJ (2014) A comparative study of the effect of probiotics on cariogenic biofilm model for preventing dental caries. Arch Microbiol 196: 601-609. doi: 10.1007/s00203-014-0998-7
    [76] Jaffar N, Ishikawa Y, Mizuno K, et al. (2016) Mature biofilm degradation by potential probiotics: Aggregatibacter actinomycetemcomitans versus Lactobacillus spp. Plos One 11: e0159466. doi: 10.1371/journal.pone.0159466
    [77] Vilela SFG, Barbosa JO, Rossoni RD, et al. (2015) Lactobacillus acidophilus ATCC 4356 inhibits biofilm formation by C. albicans and attenuates the experimental candidiasis in Galleria mellonella. Virulence 6: 29-39.
    [78] Krzyściak W, Kościelniak D, Papież M, et al. (2017) Effect of a Lactobacillus salivarius probiotic on a double-species Streptococcus mutans and Candida albicans caries biofilm. Nutrients 9: 1242. doi: 10.3390/nu9111242
    [79] Rossoni RD, de Barros PP, de Alvarenga JA, et al. (2018) Antifungal activity of clinical Lactobacillus strains against Candida albicans biofilms: Identification of potential probiotic candidates to prevent oral candidiasis. Biofouling 34: 212-225. doi: 10.1080/08927014.2018.1425402
    [80] Lin X, Chen X, Tu Y, et al. (2017) Effect of probiotic lactobacilli on the growth of Streptococcus mutans and multispecies biofilms isolated from children with active caries. Med Sci Monit 23: 4175-4181. doi: 10.12659/MSM.902237
    [81] Song YG, Lee SH (2017) Inhibitory effects of Lactobacillus rhamnosus and Lactobacillus casei on Candida biofilm of denture surface. Arch Oral Biol 76: 1-6. doi: 10.1016/j.archoralbio.2016.12.014
    [82] Ciandrini E, Campana R, Baffone W (2017) Live and heat-killed Lactobacillus spp. interfere with Streptococcus mutans and Streptococcus oralis during biofilm development on titanium surface. Arch Oral Biol 78: 48-57.
    [83] Comelli EM, Guggenheim B, Stingele F, et al. (2002) Selection of dairy bacterial strains as probiotics for oral health. Eur J Oral Sci 110: 218-224. doi: 10.1034/j.1600-0447.2002.21216.x
    [84] James KM, MacDonald KW, Chanyi RM, et al. (2016) Inhibition of Candida albicans biofilm formation and modulation of gene expression by probiotic cells and supernatant. J Med Microbiol 65: 328-336. doi: 10.1099/jmm.0.000226
    [85] Fernández CE, Giacaman RA, Tenuta LM, et al. (2015) Effect of the probiotic Lactobacillus rhamnosus LB21 on the cariogenicity of Streptococcus mutans UA159 in a dual-species biofilm model. Caries Res 49: 583-590. doi: 10.1159/000439315
    [86] Wu CC, Lin CT, Wu CY, et al. (2015) Inhibitory effect of Lactobacillus salivarius on Streptococcus mutans biofilm formation. Mol Oral Microbiol 30: 16-26. doi: 10.1111/omi.12063
    [87] Tahmourespour A, Kermanshahi RK (2011) The effect of a probiotic strain (Lactobacillus acidophilus) on the plaque formation of oral streptococci. Bosn J Basic Med Sci 11: 37-40. doi: 10.17305/bjbms.2011.2621
    [88] Vacca C, Contu MP, Rossi C, et al. (2020) In vitro interactions between Streptococcus intermedius and Streptococcus salivarius K12 on a titanium cylindrical surface. Pathogens 9: 1069. doi: 10.3390/pathogens9121069
    [89] Fang K, Jin X, Hong SH (2018) Probiotic Escherichia coli inhibits biofilm formation of pathogenic E. coli via extracellular activity of DegP. Sci Rep 8: 4939.
    [90] Hancock V, Dahl M, Klemm P (2010) Probiotic Escherichia coli strain Nissle 1917 outcompetes intestinal pathogens during biofilm formation. J Med Microbiol 59: 392-399. doi: 10.1099/jmm.0.008672-0
    [91] Kıvanç M, Er S (2020) Biofilm formation of Candida spp. isolated from the vagina and antibiofilm activities of lactic acid bacteria on the these Candida isolates. Afr Health Sci 20: 641-648.
    [92] Matsubara VH, Wang Y, Bandara HMHN, et al. (2016) Probiotic lactobacilli inhibit early stages of Candida albicans biofilm development by reducing their growth, cell adhesion, and filamentation. Appl Microbiol Biot 100: 6415-6426. doi: 10.1007/s00253-016-7527-3
    [93] Alexandre Y, Le Berre R, Barbier G, et al. (2014) Screening of Lactobacillus spp. for the prevention of Pseudomonas aeruginosa pulmonary infections. BMC Microbiol 14: 107.
    [94] Song H, Zhang J, Qu J, et al. (2019) Lactobacillus rhamnosus GG microcapsules inhibit Escherichia coli biofilm formation in coculture. Biotechnol Lett 41: 1007-1014. doi: 10.1007/s10529-019-02694-2
    [95] Sambanthamoorthy K, Feng X, Patel R, et al. (2014) Antimicrobial and antibiofilm potential of biosurfactants isolated from lactobacilli against multi-drug-resistant pathogens. BMC Microbiol 14: 197. doi: 10.1186/1471-2180-14-197
    [96] Kaur S, Sharma P, Kalia N, et al. (2018) Anti-biofilm properties of the fecal probiotic lactobacilli against Vibrio spp. Front Cell Infect Mi 8: 120. doi: 10.3389/fcimb.2018.00120
    [97] Otero MC, Nader-Macías ME (2006) Inhibition of Staphylococcus aureus by H2O2-producing Lactobacillus gasseri isolated from the vaginal tract of cattle. Anim Reprod Sci 96: 35-46. doi: 10.1016/j.anireprosci.2005.11.004
    [98] Azeredo J, Azevedo NF, Briandet R, et al. (2017) Critical review on biofilm methods. Crit Rev Microbiol 43: 313-351. doi: 10.1080/1040841X.2016.1208146
    [99] Gomes M, Gomes LC, Santos RT, et al. (2020) PDMS in urinary tract devices: Applications, problems and potential solutions, In: Carlsen PN, 1st Ed, Polydimethylsiloxane: Structure and Applications, Nova Science Publishers.
    [100] Khan F, Tabassum N, Kim YM (2020) A strategy to control colonization of pathogens: Embedding of lactic acid bacteria on the surface of urinary catheter. Appl Microbiol Biot 104: 9053-9066. doi: 10.1007/s00253-020-10903-6
    [101] Zhu Z, Yu F, Chen H, et al. (2017) Coating of silicone with mannoside-PAMAM dendrimers to enhance formation of non-pathogenic Escherichia coli biofilms against colonization of uropathogens. Acta Biomater 64: 200-210. doi: 10.1016/j.actbio.2017.10.008
    [102] Lopez AI, Kumar A, Planas MR, et al. (2011) Biomaterials biofunctionalization of silicone polymers using poly(amidoamine) dendrimers and a mannose derivative for prolonged interference against pathogen colonization. Biomaterials 32: 4336-4346. doi: 10.1016/j.biomaterials.2011.02.056
    [103] Gross G, van Der Meulen J, Snel J, et al. (2008) Mannose-specific interaction of Lactobacillus plantarum with porcine jejunal epithelium. FEMS Immunol Med Mic 54: 215-223. doi: 10.1111/j.1574-695X.2008.00466.x
    [104] Barth KA, Coullerez G, Nilsson LM, et al. (2008) An engineered mannoside presenting platform: Escherichia coli adhesion under static and dynamic conditions. Adv Funct Mater 18: 1459-1469. doi: 10.1002/adfm.200701246
    [105] Tan L, Fu J, Feng F, et al. (2020) Engineered probiotics biofilm enhances osseointegration via immunoregulation and anti-infection. Sci Adv 6: eaba5723.
    [106] Piqué N, Berlanga M, Miñana-Galbis D (2019) Health benefits of heat-killed (Tyndallized) probiotics: An overview. Int J Mol Sci 20: 2534. doi: 10.3390/ijms20102534
  • Reader Comments
  • © 2021 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(6412) PDF downloads(390) Cited by(11)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog