Loading [Contrib]/a11y/accessibility-menu.js
Mini review Special Issues

RodZ: a key-player in cell elongation and cell division in Escherichia coli

  • RodZ is required for determination of cell shape in rod-shaped bacterium, such as Escherichia coli. RodZ is a transmembrane protein and forms a supramolecular complex called the Rod complex with other proteins, such as MreB-actin and peptidoglycan synthesis enzymes (for e.g., PBP2). Deletion of the rodZ gene changes the cell shape from rod to round or ovoid. Another supramolecular complex called divisome that controls cell division mainly consists of FtsZ-tubulin. MreB directly interacts with FtsZ and this interaction is critical to trigger a transition from cell elongation to cell division. Recently, we found that RodZ also directly interacts with FtsZ, and RodZ recruits MreB to the divisome. Formation of the division ring, called Z ring, is delayed if RodZ does not interact with FtsZ, indicating that RodZ might facilitate the formation of the Z ring during the cell division process. In this mini-review, we have summarized the roles of RodZ in cell elongation and cell division, especially based on our recent study.

    Citation: Risa Ago, Daisuke Shiomi. RodZ: a key-player in cell elongation and cell division in Escherichia coli[J]. AIMS Microbiology, 2019, 5(4): 358-367. doi: 10.3934/microbiol.2019.4.358

    Related Papers:

    [1] Abdennasser Chekroun, Mohammed Nor Frioui, Toshikazu Kuniya, Tarik Mohammed Touaoula . Global stability of an age-structured epidemic model with general Lyapunov functional. Mathematical Biosciences and Engineering, 2019, 16(3): 1525-1553. doi: 10.3934/mbe.2019073
    [2] Shaoli Wang, Jianhong Wu, Libin Rong . A note on the global properties of an age-structured viral dynamic model with multiple target cell populations. Mathematical Biosciences and Engineering, 2017, 14(3): 805-820. doi: 10.3934/mbe.2017044
    [3] Yan-Xia Dang, Zhi-Peng Qiu, Xue-Zhi Li, Maia Martcheva . Global dynamics of a vector-host epidemic model with age of infection. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1159-1186. doi: 10.3934/mbe.2017060
    [4] Kento Okuwa, Hisashi Inaba, Toshikazu Kuniya . Mathematical analysis for an age-structured SIRS epidemic model. Mathematical Biosciences and Engineering, 2019, 16(5): 6071-6102. doi: 10.3934/mbe.2019304
    [5] Xia Wang, Shengqiang Liu . Global properties of a delayed SIR epidemic model with multiple parallel infectious stages. Mathematical Biosciences and Engineering, 2012, 9(3): 685-695. doi: 10.3934/mbe.2012.9.685
    [6] Yoichi Enatsu, Yukihiko Nakata, Yoshiaki Muroya . Global stability for a class of discrete SIR epidemic models. Mathematical Biosciences and Engineering, 2010, 7(2): 347-361. doi: 10.3934/mbe.2010.7.347
    [7] Abdelheq Mezouaghi, Salih Djillali, Anwar Zeb, Kottakkaran Sooppy Nisar . Global proprieties of a delayed epidemic model with partial susceptible protection. Mathematical Biosciences and Engineering, 2022, 19(1): 209-224. doi: 10.3934/mbe.2022011
    [8] Jianxin Yang, Zhipeng Qiu, Xue-Zhi Li . Global stability of an age-structured cholera model. Mathematical Biosciences and Engineering, 2014, 11(3): 641-665. doi: 10.3934/mbe.2014.11.641
    [9] Yicang Zhou, Zhien Ma . Global stability of a class of discrete age-structured SIS models with immigration. Mathematical Biosciences and Engineering, 2009, 6(2): 409-425. doi: 10.3934/mbe.2009.6.409
    [10] Xi-Chao Duan, Xue-Zhi Li, Maia Martcheva . Dynamics of an age-structured heroin transmission model with vaccination and treatment. Mathematical Biosciences and Engineering, 2019, 16(1): 397-420. doi: 10.3934/mbe.2019019
  • RodZ is required for determination of cell shape in rod-shaped bacterium, such as Escherichia coli. RodZ is a transmembrane protein and forms a supramolecular complex called the Rod complex with other proteins, such as MreB-actin and peptidoglycan synthesis enzymes (for e.g., PBP2). Deletion of the rodZ gene changes the cell shape from rod to round or ovoid. Another supramolecular complex called divisome that controls cell division mainly consists of FtsZ-tubulin. MreB directly interacts with FtsZ and this interaction is critical to trigger a transition from cell elongation to cell division. Recently, we found that RodZ also directly interacts with FtsZ, and RodZ recruits MreB to the divisome. Formation of the division ring, called Z ring, is delayed if RodZ does not interact with FtsZ, indicating that RodZ might facilitate the formation of the Z ring during the cell division process. In this mini-review, we have summarized the roles of RodZ in cell elongation and cell division, especially based on our recent study.


    Nomenclature

    acrack depth
    a1, a2, athe depths of crack 1, crack 2 and the combined crack, respectively
    a0initial crack depth
    a1,0, a2,0initial depths of cracks 1 and 2, respectively
    chalf length of a crack
    c1, c2, chalf lengths of crack 1, crack 2 and the combined crack, respectively
    c0initial half length of a crack
    c1,0, c2,0initial half lengths of cracks 1 and 2, respectively
    Cmaterial constant in Eqn. (1)
    da/dNfatigue crack growth rate
    Eelastic modulus
    Hhalf height of a plate
    Jenergy release rate
    KIstress intensity factor for mode I loading
    K0normalized stress intensity factor
    KIrange of stress intensity factor
    mmaterial constant in Eqn. (1)
    Rstress ratio
    sdistance between two adjacent flaws
    s0initial minimum distance between two cracks
    tthickness of a plate
    Whalf width of a plate
    vPoisson’s ratio
    σapplied stress


    1. Introduction

    It is not uncommon to detect multiple flaws in structural components [1,2,3]. As reported by Nuclear and Industry Safety Agency, about two thirds of the stress corrosion cracks were existing as multiple flaws in a welded austenitic stainless steel pipe [4]. How to properly assess the structure containing multiple cracks is of great importance to preventing the potential catastrophic failure. For two coplanar surface cracks as shown in Figure 1, essentially based on the enhancement of stress intensity factor (e.g., [5,6]) and some experimental observations (e.g., [7,8]), guidance for the characterization of multiple flaws as displayed in Table 1, have been provided by the current fitness-for-service (FFS) procedures [9,10,11,12]. Although their basic concepts (i.e., when the distance between two adjacent flaws is equal to or less than a specific value, the two flaws will be combined into a virtual larger flaw) are similar, apparent discrepancy can be seen between various combination rules. Which of them is more reasonable and whether there are better options remain open questions [13].

    Figure 1. Finite thickness plate containing two coplanar surface flaws: (a) Geometric model; (b) Definition for location and shape of two flaws and a combined flaw. The surface flaws were idealized as semi-elliptical cracks. W, H and t are the half width, half height and thickness of the plate, respectively; s is the distance between two adjacent flaws; a1, a2 and a’ denote the depths of crack 1, crack 2 and the combined crack, respectively; c1, c2 and c’ denote the half lengths of crack 1, crack 2 and the combined crack, respectively.
    Table 1. Combination rules for coplanar surface flaws in existing FFS codes.
    SourcesCombination rules
    ASME [9]Fatigue or stress corrosion: s = 0
    Others: s ≤ max(0.5a1,0.5a2)
    BS7910 [10]Fatigue: s = 0
    Others: s ≤ min(2c1,2c2) for a1/c1 or a2/c2 > 1
    s ≤ max(0.5a1,0.5a2) for a1/c1 ≤ 1 and a2/c2 ≤ 1
    API579-1 [11]sc1 + c2
    GB/T 19624 [12]s ≤ min(2c1,2c2)
     | Show Table
    DownLoad: CSV

    In the past four decades, a number of researchers have examined the interacting effect of adjacent cracks by fatigue crack growth tests. The pioneer experimental work of Iida et al. [8] in 1976 established an estimation procedure for fatigue life of a plate with multiple surface flaws. Soboyejo et al. [14] investigated the effects of stress ratio and initial crack separation on the interaction and coalescence of coplanar cracks under cyclic bending. By conducting complicated fatigue tests on specimens containing irregularly oriented surface cracks, Tu et al. [15] suggested a criterion for defect re-characterization. Kamaya [16] carried out fatigue crack growth tests for both coplanar and non-coplanar notches and argued that the growth of interacting surface cracks may be dependent on the area of crack face. Considerable efforts have also been devoted into finite element (FE) analyses which are obviously more cost-effective. It should be pointed that most numerical studies [17,18,19,20] only investigated the interacting effect on the crack tip parameters for stationary cracks, which are not necessarily relevant to the fatigue life of components. The attempts made by Lin et al. [21] and Kamaya [16] are notable exceptions that have simulated the crack growth behavior during the coalescence process for limited cases. However, how to rationalize a criterion for the re-characterization of multiple flaws from a series of numerical simulations remains a challenging task.

    The purposes of this study is to assess the conservatism contained in different combination rules given by current FFS codes, and to develop a more rational one, based on the crack growth life of cracked structures. It should be noted that various failure mechanisms (cleavage, plastic collapse, fatigue, creep, corrosion, etc.) may significantly influence the flaw interaction criteria. To help clarify the description and obtain a clear understanding, our attention herein focuses specifically on the applicability of combination rules over the fatigue failure mode. In Section 2, the numerical method employed to simulate the fatigue crack growth is introduced, and then is validated by classic equations for stress intensity factors and by a fatigue test for a cracked plate. In Section 3, fatigue growth analysis of two identical cracks is performed. Furthermore, the conservatism contained in different combination rules provided by fitness-for-service codes and by the present study have been quantitatively assessed. Then, a numerical analysis for two dissimilar cracks is conducted in Section 4. Finally, several concluding remarks are drawn in Section 5.


    2. Materials and Method

    As depicted in Figure 1, a plate containing two coplanar surface flaws under cyclic loading was considered in this study. The surface flaws were idealized as semi-elliptical cracks. To eliminate the effect of the boundary on the calculation, the width and height of the plate were selected to be much larger than the crack dimensions. In this study, the number of cycles required for the re-characterized flaw to penetrate the plate was regarded as the life of the plate. The fatigue crack growth rate (mm×cycle-1) can be described by the well-known Paris equation

    $\frac{{{\text{d}}a}}{{{\text{d}}N}} = C{(\Delta {K_I})^m}$ (1)

    where C (mm×cycle-1 (MPa×mm0.5)-m) and m are material constants. △KI (MPa×mm0.5) denotes the range of stress intensity factor.

    To trace the change of crack shape, local crack advance at an arbitrary point of crack front along the normal direction was computed. Two distinct methods were employed to predict the fatigue crack growth life the plate containing multiple cracks:

    (I) Interaction Method:

    A detailed step-by-step FE analysis was conducted to simulate the propagation, interaction and coalescence of multiple cracks. The interaction in the period of growth up to cracks touching, as well as the saturation of touching cracks to the achievement of a bounding crack, was taken into account. The routine of numerical simulations is described as follows:

    (a) Create the FE model. A fracture mechanics code, ZENCRACK [22], was adopted to create the FE model. In order to obtain accurate numerical results, extremely refined meshes are generated in the crack tip zone.

    (b) Calculate KI. The values of energy release rate, J-integral, at a set of points that constitute the crack front was given by a FE code ABAQUS [23] using the domain integral method. For mode I loading, stress intensity factor, KI, with plane strain assumption can be obtained by

    ${K_I} = \sqrt {\frac{{EJ}}{{1 - {v^2}}}} $ (2)

    where E and v are elastic modulus and Poisson’s ratio, respectively.

    (c) Calculate the increment of crack size. The fatigue crack growth at each point along the crack front in each step can be given by Eqn. (1).

    (d) Establish a new crack profile. Based on the crack advancement at each node, the coordinates of a set of new nodes can be determined. The new crack profile was automatically created by ZENCRACK until the contact of two cracks. Note that at this moment the analysis would be suspended by the distorted mesh. A manual intervention was needed to slightly adjust the crack shape [24]. After that, ZENCRACK would continue to work.

    (e) Return to stage (a). The cracks propagate through the stages from (a) to (d) until they penetrate the plate.

    (II) Non-Interaction Method:

    A similar step-by-step FE analysis was also performed but each crack was assumed to grow independently before the coalescence. Once a combination rule listed in Table 1 has been satisfied, adjacent cracks would immediately join together into a larger crack, with a depth of the deeper of the two pre-existing cracks and a length being equal to the sum of their lengths and minimum distance, as illustrated in Figure 1(b).

    To check the validity of the FE model, stress intensity factors at a set of points that constitute the crack front calculated by ABAQUS have been compared with the Newman-Raju solutions [25], as shown in Figure 2, for semi-elliptical cracks (with aspect ratio a/c = 0.25 and depth ratio a/t = 0.2, 0.4, 0.6 and 0.8) in a plate (with thickness t = 15mm, half width W = 6.67t and half height H = 6.67t) subjected to a tension load (with stress σ = 200MPa). A typical FE mesh of a plate with an isolated crack is depicted in Figure 3(a). Elastic modulus and Poisson’s ratio are given in Table 2. The FE model contains 6517 eight-node isoparametric elements. Stress intensity factors have been normalized by ${K_0} = \sigma \sqrt {\pi t} $. One can see that the present FE solutions agree very well with the Newman-Raju solutions for a/t = 0.2, 0.4 and 0.6, with errors less than 1.2%. For the case with a/t = 0.8 when $2\phi /\pi > 0.6$, stress intensity factors from the Newman-Raju equations are slightly larger (within 7.0%) than the FE results. Similar difference between FE results and the Newman-Raju solutions for a/t = 0.8 has also been reported by Lei. It is believed that this discrepancy is due to the fitting error introduced by Newman-Raju equations.

    Figure 2. Comparison of normalized stress intensity factors at points along surface cracks in a plate subjected to a tension load between present FE solutions and Raju-Newman solutions.
    Figure 3. Typical finite element meshes of plates with: (a) An isolated crack (1/4 model); (b) Two identical cracks (1/4 model); (c) Two dissimilar cracks (1/2 model).
    Table 2. Material constants for a low carbon alloyed steel 16MnR [26].
    Elastic modulusPoisson’s ratioCoefficient in Eqn. (1)Exponent in Eqn. (1)
    E/MPaνC/mm×cycle–1 (MPa×mm0.5)mm
    2090000.287.30 × 10–153.42
     | Show Table
    DownLoad: CSV

    In order to validate the step-by-step numerical method employed, comparison between predictions and experimental results [26] has also been made for a plate (with half width W = 44.75 mm, half height H = 175 mm, and thickness t = 7 mm) containing a single crack (with initial depth a0 = 3.19 mm and half length c0 = 3.47 mm) under a cyclic loading (with maximum and minimum load, 140 kN and 0 kN, respectively). The initial surface defect was manufactured by electric discharge machining. The specimen was tested on a closed loop, servo-hydraulic testing system MTS 880. Sinusoidal tension loading at a frequency of 22 Hz was employed. To determine the crack shape evolution during the propagation, a crack front marking procedure was used at each interval of crack extension of about 0.5-1 mm, in which the maximum load was kept constant and the load range was reduced by one half. After the fatigue test was done, the specimen was opened and then heated to improve the clarity of the beach marks. The material of the plate is 16MnR, a widely-used low carbon alloyed steel for pressure vessels. Material constants used in this study are listed in Table 2 [26]. Due to the symmetry of both geometry and loading, only one quarter of the specimen has been created. As shown in Figure 3(a), the FE mesh contained a total of 7836 C3D8 elements. Figure 4 shows the comparison of predicted crack profiles with those on the fracture surface after test [26]. Clearly, the crack shapes predicted by the current step-by-step FE analysis are extremely similar to the experimental evidence, even after the plate thickness penetration of crack. A further comparison can be made on the variation of aspect ratio with depth growth throughout the plate thickness, as illustrated in Figure 5. The maximum error of aspect ratios, a/c, between prediction and test data for a given depth ratio, a/t, is only about 1.5%. These encouraging results prove the validity and predictive capability of the numerical analysis used in this study.

    Figure 4. Comparison of tested specimen photo [26] with predicted crack profile evolution for a low carbon alloyed steel 16MnR.
    Figure 5. Variation of aspect ratio with depth growth throughout the plate thickness for a low carbon alloyed steel 16MnR.

    3. Fatigue Growth Analysis of Two Identical Cracks

    We started the investigation on multiple surface cracks with identical sizes. These cases may be regarded as probably more dangerous since the interaction effect is strong [5]. As shown in Figure 1, a finite thickness plate (with t = 15 mm, W = 6.67t and H = 6.67t) subjected to a cyclic load (with maximum stress σ = 200 MPa and stress ratio R = 0.1) was considered. The same material constants as listed in Table 2 were chosen. Two cracks (with a2,0/t = 0.20, a1,0/c1,0 = a2,0/c2,0 = 0.50, a1,0/a2,0 = 1.00 and s0/c2,0 = 2.20, where a1,0 and a2,0 are the initial depths of cracks 1 and 2, respectively; c1,0 and c2,0 are the initial half lengths of cracks 1 and 2, respectively; s0 is the initial minimum distance between two cracks; a0 = a1,0 = a2,0 and c0 = c1,0 = c2,0 where a0 and c0 are the initial crack depth and initial crack half length, respectively, when two cracks are identical to each other.) situated in the middle of the plate. With symmetry conditions imposed only one quarter of the plate model needed to be created using 13046 eight-node brick elements, as shown in Figure 3(b). Note that meshes were rebuilt in each step according to the updated crack profile. Meshes of selected finite element models (1/2 model) in various steps are shown in Figure 6. Figure 7 provides the shape development of two identical cracks predicted by the Interaction Method. It can be observed that the individual cracks propagate almost independently until the minimum distance between them becomes comparable with the crack depth. Then, crack grows at a higher rate in the neighboring zone than in other regions. When the coalescence begins, dented portions are formed, leading to a crack profile deviating from a semi-ellipse. In a short time, the concave points along the crack front disappear and the shape of newly-combined crack becomes smooth ultimately. The predicted shape evolution can be found to be in good agreement with many experimental observations (e.g., [16,26]).

    Figure 6. Meshes of selected finite element models (1/2 model) in various steps: (a) Initial step; (b) Intermediate step; (c) Final step.
    Figure 7. Predicted shape evolution of two identical cracks by using the Interaction Method (for two initial cracks with a2,0/t = 0.20, a1,0/c1,0 = a2,0/c2,0 = 0.50, a1,0/a2,0 = 1.00 and s0/c2,0 = 2.20).

    The phenomenon seen in the Figure 7 can be explained by investigating the developments of stress intensity factor at depth point and surface point in the crack growth process, as shown in Figure 8. Note that stress intensity factors have been normalized by ${K_0} = \sigma \sqrt {\pi t} $. Locations of A and B are marked in the imbedded figures on different stages. It can be seen from Figure 8 that the stress intensity factors at A and B are identical initially, whereas their discrepancy gradually become significant as the two cracks get closer to each other. The stress intensity factor of the merging point, A, reaches a maximum value of approximately triple that at the surface point, B, soon after the contacting moment, and decreases until the crack shape become smooth again. On the final stage, the stress intensity factor at A rises due to the reduction of ligament area.

    Figure 8. Variations of normalized stress intensity factor at depth and surface points with the depth growth throughout the plate thickness (for two initial cracks with a2,0/t = 0.20, a1,0/c1,0 = a2,0/c2,0 = 0.50, a1,0/a2,0 = 1.00 and s0/c2,0 = 2.20). In the embedded figure, the period “during coalescence” is the time between the touch of two cracks and the saturation of the newly combined crack. “After coalescence” means the period from the saturation moment onward.

    Figure 9 shows the comparison of crack growth curves predicted by the Interaction Method and the Non-Interaction Methods with different combination rules listed in Table 1, corresponding to the case of Figure 7. The fatigue cycles, N, has been normalized by the fatigue life, N0, taken for crack to penetrate the plate, given by the Interaction Method. One can find that the fatigue life predicted by the Non-Interaction Method with the combination rule, s = 0, in ASME [9] and BS7910 [10] is overestimated by approximate 15%. This indicates the interaction between two adjacent cracks cannot be ignored, otherwise the assessment may be unsafe. On the other hand, for the same case the criterion, s≤2c, provided by API579 [11] and GB/T19624 [12], leads to an excessively conservative estimation (only one half of the life given by the Interaction Method). None of the available combination rules is satisfying, at least in this case. Therefore, partly based on experience gained from the numerical solutions here and those to be presented later, we propose a new combination rule.

    $s \leq ({a_1} + {a_2})/2$ (3)

    The result predicted by the proposed combination rule (degraded to sa for two identical cracks) is also shown in Figure 9. Clearly, the Non-Interaction Method with proposed combination rule gives a best estimation of crack growth curve.

    Figure 9. Comparison of crack growth curves predicted by the Interaction Method and the Non-Interaction Methods with different combination rules (for two initial identical cracks with a2,0/t = 0.20, a1,0/c1,0 = a2,0/c2,0 = 0.50, a1,0/a2,0 = 1.00 and s0/c2,0 = 2.20).

    Numerical analyses were carried out for more cases with various initial crack shapes (a0/c0 = a1,0/c1,0 = a2,0/c2,0 = 0.25,0.50,0.75 and 1.00) and initial crack depth ratios (a0/t = a1,0/t = a2,0/t = 0.1,0.2,0.3 and 0.4). Note that the two cracks are initially separated at a relatively long distance (s0/c2,0 = 2.20). If the initial crack depth ratio, a0/t, is large enough, the cracks will individually penetrate the plate without coalescence. For this reason, the initial crack depth ratio was selected to be no larger than 0.4 in the present study. The normalized fatigue lives, N/N0, due to different combination rules are shown in Figure 10. As expected, the fatigue life obtained by the proposed combination rule is closer to the numerical results considering the interaction effect, whereas estimations given by ASME and BS7910 are essentially not conservative and those according to API579 and GB/T19624 are of excessive conservatism. It is remarkable that, with the increase of initial crack depth ratio, both non-conservatism due to s = 0 and over-conservatism from s≤2c rises, while the proposed combination rule always keeps a reasonable prediction.

    Figure 10. Variations of the normalized fatigue life with the initial crack depth ratio for cases with two identical surface cracks (a1,0/a2,0 = 1.00 and s0/c2,0 = 2.20).

    A parametric study was also performed for different materials with the exponent in the Paris equation, m, ranging from 2.5 to 4.0, which is the approximate range for most metals in the atmospheric environment [27]. As shown in Figure 11, for m = 4, the numerical results from the Non-Interaction Method with combination rule in ASME and BS7910 can lead to a maximum unsafe deviation of 40%, while API579 and GB/T19624 may underestimate the fatigue life by approximate 60%. Again, the results given by proposed criterion coincide with those considering interaction, and the ratio of the former to the latter is in a reasonable range of 80-100%.

    Figure 11. Effect of material constant, m, on the normalized fatigue life (for two initial identical cracks with a2,0/t = 0.40, a1,0/c1,0 = a2,0/c2,0 = 1, a1,0/a2,0 = 1.00 and s0/c2,0 = 2.20) with other material constants fixed.

    4. Fatigue Growth Analysis of Two Dissimilar Cracks

    The same finite thickness plate and loading condition, as described in Section 3, but with two dissimilar coplanar surface cracks were considered in this section. Due to symmetry one half of the plate was taken into account in the FE models using 13460 eight-node brick elements, as shown in Figure 3(c). Material constants of 16MnR listed in Table 2 were used again. Figure 12 illustrates the shape evolution of two dissimilar cracks (with initial sizes a2,0/t = 0.20, a1,0/c1,0 = a2,0/c2,0 = 0.50, a1,0/a2,0 = 0.60 and s0/c2,0 = 2.20) predicted by the Interaction Method. Crack propagation behavior similar to that shown in Figure 7 can be observed. What differs from the former is that the individual dissimilar cracks grow more independently as isolated flaws. It can be seen from Figure 7 that the spacing between the initial and touching crack shapes on the interacting side is significantly larger than that on the non-interacting side due to the interaction of two cracks. For the larger crack in Figure 12, by way of contrast, the intervals between crack profiles are almost the same on both interacting and non-interacting sides. That means the propagation of the larger crack seems less affected by the smaller crack. Obviously, the larger crack advance much more rapidly since the stress intensity factors at the crack front are always greater than those for smaller crack during the growth process. Similar shape development of two cracks with different sizes can be found in Ref. [14,28].

    Figure 12. Predicted shape evolution of two dissimilar cracks by using the Interaction Method (for two initial cracks with a2,0/t = 0.20, a1,0/c1,0 = a2,0/c2,0 = 0.50, a1,0/a2,0 = 0.60 and s0/c2,0 = 2.20).

    Figure 13 displays the comparison of crack growth curves, corresponding to the case of Figure 12, predicted by the Interaction Method and the Non-Interaction Methods with different combination rules. It demonstrates that, for this case, all combination rules lead to conservative estimations when compared with the results by the Interaction Method. Nonetheless, API579’s criterion (sc1 + c2) underestimates the fatigue life by approximate 30%.

    Figure 13. Comparison of crack growth curves predicted by the Interaction Method and the Non-Interaction Methods with different combination rules (for two initial dissimilar cracks with a2,0/t = 0.20, a1,0/c1,0 = a2,0/c2,0 = 0.50, a1,0/a2,0 = 0.60 and s0/c2,0 = 2.20).

    More cases with various crack sizes ratios (a1,0/a2,0 = 0.40, 0.60, 0.80 and 1.00) were also considered in the study. As can be seen in Figure 14, when the sizes of two cracks are considerably dissimilar (a1,0/a2,0≤0.60), the combination rule suggested by ASME and BS7910 (s = 0) will give estimations very close to the numerical results that considers crack interaction. This is because for these cases the two cracks can be regarded as independent to some extent during the propagation process up to cracks touching, as demonstrated in Figure 12. However, when the crack sizes are somewhat similar (a1,0/a2,0≥0.80), the effect of interaction between two adjacent cracks cannot be neglected, as illustrated in Section 3. In this circumstance, ASME and BS7910 will overestimate the fatigue life. By contrast, API 579 (sc1 + c2) always underestimate the results by at least 25%. Interestingly, conservatism of GB/T 19624 (s≤min(2c1, 2c2) or s = 2c1) is highly dependent on the crack size ratio. The noteworthy result is that, the proposed combination rule (s≤(a1 + a2)/2) always yields a reasonable estimation with necessary conservatism, no matter whether the crack sizes are similar or not.

    Figure 14. Variations of the normalized fatigue life with the initial crack size ratio for cases with two dissimilar surface cracks (a2,0/t = 0.20, a1,0/c1,0 = a2,0/c2,0 = 0.50 and s0/c2,0 = 2.20).

    5. Conclusions

    The combination rules for multiple coplanar fatigue cracks provided by different existing FFS codes (ASME, BS7910, API579 and GB/T19624) have been examined by detailed step-by-step FE analyses. In addition, a new combination rule (s≤(a1 + a2)/2) based on the crack growth life of structure has been proposed. Conclusions can be drawn as follows:

    (1) The amount of conservatism contained in the combination rules given by existing FFS procedures is dependent on the initial crack depth ratio, a0/t, and initial crack size ratio, a1,0/a2,0.

    (2) Predictions using the re-characterization guideline (s = 0) provided by ASME and BS7910 can be non-conservative when initial sizes of two cracks are similar.

    (3) Combination rules suggested by API579 (sc1 + c2) and GB/T19624 (s≤2c1) lead to excessively pessimistic predictions of crack growth for almost all the cases considered in the present study.

    (4) The proposed combination rule (s≤(a1 + a2)/2) always yields a reasonable estimation with necessary conservatism, for various initial crack depths, material constants and relative sizes of two cracks.

    Be aware of that the findings of this research may be only valid for fatigue failure mode. It would be interesting to undertake further studies for other failure mechanisms such as creep and corrosion, based on the same concept and similar numerical method.


    Acknowledgements

    Supports from National Natural Science Foundation of China (Grant Nos. 11472105 and 51505149) and the “111 project” are gratefully acknowledged. J.-F also wishes to thank the financial support provided by Shanghai Sailing Program (Grant No. 15YF1402900) and China Postdoctoral Science Foundation (Grant No. 2015M581543). Helpful discussion on the numerical technique with Mr. Chun-Lei Ma at Consys Ltd. And Dr. Ramesh Chandwani at Zentech Int. Ltd. is also greatly appreciated.


    Conflict of Interest

    The authors declare that there is no conflict of interest regarding the publication of this manuscript.



    Acknowledgments



    The work on RodZ was partly carried out in Dr. Hironori Niki's lab (National Institute of Genetics, Japan). The work in our lab was supported by the JSPS KAKENHI (Grant Numbers 24770191, 25117528, 15H01333, and 15H04731) and by the NIG-JOINT (2016-A1[40], 61A2017, 57A2018, and 60A2019).

    Conflict of interest



    The authors declare no conflict of interest in this review.

    [1] Szwedziak P, Löwe J (2013) Do the divisome and elongasome share a common evolutionary past? Curr Opin Microbiol 16: 745–751. doi: 10.1016/j.mib.2013.09.003
    [2] Blaauwen den T, de Pedro MA, Nguyen-Distèche M, et al. (2008) Morphogenesis of rod-shaped sacculi. FEMS Microbiol Rev 32: 321–344. doi: 10.1111/j.1574-6976.2007.00090.x
    [3] Figge RM, Divakaruni AV, Gober JW (2004) MreB, the cell shape-determining bacterial actin homologue, co-ordinates cell wall morphogenesis in Caulobacter crescentus. Mol Microbiol 51: 1321–1332. doi: 10.1111/j.1365-2958.2003.03936.x
    [4] Vats P, Rothfield L (2007) Duplication and segregation of the actin (MreB) cytoskeleton during the prokaryotic cell cycle. Proc Natl Acad Sci USA 104: 17795–17800. doi: 10.1073/pnas.0708739104
    [5] Vats P, Shih YL, Rothfield L (2009) Assembly of the MreB-associated cytoskeletal ring of Escherichia coli. Mol Microbiol 72: 170–182. doi: 10.1111/j.1365-2958.2009.06632.x
    [6] van der Ploeg R, Verheul J, Vischer NOE, et al. (2013) Colocalization and interaction between elongasome and divisome during a preparative cell division phase in Escherichia coli. Mol Microbiol 87: 1074–1087. doi: 10.1111/mmi.12150
    [7] Fenton AK, Gerdes K (2013) Direct interaction of FtsZ and MreB is required for septum synthesis and cell division in Escherichia coli. EMBO J 32: 1953–1965. doi: 10.1038/emboj.2013.129
    [8] Yoshii Y, Niki H, Shiomi D (2019) Division-site localization of RodZ is required for efficient Z ring formation in Escherichia coli. Mol Microbiol 111: 1229–1244. doi: 10.1111/mmi.14217
    [9] Shiomi D, Sakai M, Niki H (2008) Determination of bacterial rod shape by a novel cytoskeletal membrane protein. EMBO J 27: 3081–3091. doi: 10.1038/emboj.2008.234
    [10] Bendezú FO, Hale CA, Bernhardt TG, et al. (2009) RodZ (YfgA) is required for proper assembly of the MreB actin cytoskeleton and cell shape in E. coli. EMBO J 28: 193–204. doi: 10.1038/emboj.2008.264
    [11] Alyahya SA, Alexander R, Costa T, et al. (2009) RodZ, a component of the bacterial core morphogenic apparatus. Proc Natl Acad Sci USa 106: 1239–1244. doi: 10.1073/pnas.0810794106
    [12] Bendezú FO, de Boer PAJ (2008) Conditional lethality, division defects, membrane involution, and endocytosis in mre and mrd shape mutants of Escherichia coli. J Bacteriol 190: 1792–811. doi: 10.1128/JB.01322-07
    [13] Muchová K, Chromiková Z, Barák I (2013) Control of Bacillus subtilis cell shape by RodZ. Environ Microbiol 15: 3259–3271. doi: 10.1111/1462-2920.12200
    [14] Okumura M, Narumi I, Nishida H (2019) Sensitivity of Deinococcus grandis rodZ deletion mutant to calcium ions results in enhanced spheroplast size. AIMS Microbiology 5: 176–185. doi: 10.3934/microbiol.2019.2.176
    [15] van den Ent F, Johnson CM, Persons L, et al. (2010) Bacterial actin MreB assembles in complex with cell shape protein RodZ. EMBO J 29: 1081–1090. doi: 10.1038/emboj.2010.9
    [16] Mitobe J, Yanagihara I, Ohnishi K, et al. (2011) RodZ regulates the post-transcriptional processing of the Shigella sonnei type III secretion system. EMBO Rep 12: 911–916. doi: 10.1038/embor.2011.132
    [17] Pereira AC, Paiva A, Saraiva IH, et al. (2015) Chemical shift assignments and secondary structure determination of the ectodomain of Bacillus subtilis morphogenic protein RodZ. Biomol NMR Assign 9: 285–288. doi: 10.1007/s12104-014-9593-8
    [18] Ikebe R, Kuwabara Y, Chikada T, et al. (2018) The periplasmic disordered domain of RodZ promotes its self-interaction in Escherichia coli. Genes Cells 23: 307–317. doi: 10.1111/gtc.12572
    [19] Daniel RA, Errington J (2003) Control of cell morphogenesis in bacteria: two distinct ways to make a rod-shaped cell. Cell 113: 767–776. doi: 10.1016/S0092-8674(03)00421-5
    [20] Philippe J, Vernet T, Zapun A (2014) The elongation of ovococci. Microb Drug Resist 20: 215–221. doi: 10.1089/mdr.2014.0032
    [21] Jones LJ, Carballido-López R, Errington J (2001) Control of cell shape in bacteria: helical, actin-like filaments in Bacillus subtilis. Cell 104: 913–922. doi: 10.1016/S0092-8674(01)00287-2
    [22] Shih YL, Le T, Rothfield L (2003) Division site selection in Escherichia coli involves dynamic redistribution of Min proteins within coiled structures that extend between the two cell poles. Proc Natl Acad Sci USA 100: 7865–7870. doi: 10.1073/pnas.1232225100
    [23] Bean GJ, Flickinger ST, Westler WM, et al. (2009) A22 disrupts the bacterial actin cytoskeleton by directly binding and inducing a low-affinity state in MreB. Biochemistry 48: 4852–4857. doi: 10.1021/bi900014d
    [24] Iwai N, Nagai K, Wachi M (2002) Novel S-benzylisothiourea compound that induces spherical cells in Escherichia coli probably by acting on a rod-shape-determining protein(s) other than penicillin-binding protein 2. Biosci Biotechnol Biochem 66: 2658–2662. doi: 10.1271/bbb.66.2658
    [25] Swulius MT, Jensen GJ (2012) The helical MreB cytoskeleton in Escherichia coli MC1000/pLE7 is an artifact of the N-Terminal yellow fluorescent protein tag. J Bacteriol 194: 6382–6386. doi: 10.1128/JB.00505-12
    [26] Bratton BP, Shaevitz JW, Gitai Z, et al. (2018) MreB polymers and curvature localization are enhanced by RodZ and predict E. coli's cylindrical uniformity. Nat Commun 9: 2797.
    [27] Morgenstein RM, Bratton BP, Nguyen JP, et al. (2015) RodZ links MreB to cell wall synthesis to mediate MreB rotation and robust morphogenesis. Proc Natl Acad Sci USA 112: 12510–12515. doi: 10.1073/pnas.1509610112
    [28] van der Ploeg R, Goudelis ST, Blaauwen den T (2015) Validation of FRET assay for the screening of growth inhibitors of Escherichia coli reveals elongasome assembly dynamics. Int J Mol Sci 16: 17637–17654. doi: 10.3390/ijms160817637
    [29] Ursell TS, Nguyen J, Monds RD, et al. (2014) Rod-like bacterial shape is maintained by feedback between cell curvature and cytoskeletal localization. Proc Natl Acad Sci USA 111: E1025–E1034. doi: 10.1073/pnas.1317174111
    [30] Kawazura T, Matsumoto K, Kojima K, et al. (2017) Exclusion of assembled MreB by anionic phospholipids at cell poles confers cell polarity for bidirectional growth. Mol Microbiol 104: 472–486. doi: 10.1111/mmi.13639
    [31] Colavin A, Shi H, Huang KC (2018) RodZ modulates geometric localization of the bacterial actin MreB to regulate cell shape. Nat Commun 9: 1280. doi: 10.1038/s41467-018-03633-x
    [32] Dion MF, Kapoor M, Sun Y, et al. (2019) Bacillus subtilis cell diameter is determined by the opposing actions of two distinct cell wall synthetic systems. Nat Microbiol 24: 96.
    [33] Garner EC, Bernard R, Wang W, et al. (2011) Coupled, circumferential motions of the cell wall synthesis machinery and MreB filaments in B. subtilis. Science 333: 222–225.
    [34] Domínguez-Escobar J, Chastanet A, Crevenna AH, et al. (2011) Processive movement of MreB-associated cell wall biosynthetic complexes in bacteria. Science 333: 225–228. doi: 10.1126/science.1203466
    [35] van Teeffelen S, Wang S, Furchtgott L, et al. (2011) The bacterial actin MreB rotates, and rotation depends on cell-wall assembly. Proc Natl Acad Sci USA 108: 15822–15827. doi: 10.1073/pnas.1108999108
    [36] Shiomi D, Toyoda A, Aizu T, et al. (2013) Mutations in cell elongation genes mreB, mrdA and mrdB suppress the shape defect of RodZ-deficient cells. Mol Microbiol 87: 1029–1044. doi: 10.1111/mmi.12148
    [37] van den Ent F, Izoré T, Bharat TA, et al. (2014) Bacterial actin MreB forms antiparallel double filaments. Elife 2014: e02634.
    [38] Pinho MG, Kjos M, Veening JW (2013) How to get (a)round: mechanisms controlling growth and division of coccoid bacteria. Nat Rev Microbiol 11: 601–614. doi: 10.1038/nrmicro3088
    [39] Land AD, Winkler ME (2011) The requirement for pneumococcal MreC and MreD is relieved by inactivation of the gene encoding PBP1a. J Bacteriol 193: 4166–4179. doi: 10.1128/JB.05245-11
    [40] Tavares AC, Fernandes PB, Carballido-Lopez R, et al. (2015) MreC and MreD Proteins are not required for growth of Staphylococcus aureus. PLoS One 10: e0140523. doi: 10.1371/journal.pone.0140523
    [41] Ouellette SP, Karimova G, Subtil A, et al. (2012) Chlamydia co-opts the rod shape-determining proteins MreB and Pbp2 for cell division. Mol Microbiol 85: 164–178. doi: 10.1111/j.1365-2958.2012.08100.x
    [42] Ouellette SP, Rueden KJ, Gauliard E, et al. (2014) Analysis of MreB interactors in Chlamydia reveals a RodZ homolog but fails to detect an interaction with MraY. Front Microbiol 5: 279.
    [43] Kemege KE, Hickey JM, Barta ML, et al. (2015) Chlamydia trachomatis protein CT009 is a structural and functional homolog to the key morphogenesis component RodZ and interacts with division septal plane localized MreB. Mol Microbiol 95: 365–382. doi: 10.1111/mmi.12855
    [44] Niba ETE, Li G, Aoki K, et al. (2010) Characterization of rodZ mutants: RodZ is not absolutely required for the cell shape and motility. FEMS Microbiol Lett 309: 35–42.
    [45] Niba ETE, Naka Y, Nagase M, et al. (2007) A genome-wide approach to identify the genes involved in biofilm formation in E. coli. DNA Res 14: 237–246.
    [46] Jovanovic G, Mehta P, Ying L, et al. (2014) Anionic lipids and the cytoskeletal proteins MreB and RodZ define the spatio-temporal distribution and function of membrane stress controller PspA in Escherichia coli. Microbiology 160: 2374–2386. doi: 10.1099/mic.0.078527-0
    [47] van Beilen J, Blohmke CJ, Folkerts H, et al. (2016) RodZ and PgsA play intertwined roles in membrane homeostasis of Bacillus subtilis and resistance to weak organic acid stress. Front Microbiol 7: 1015–1012.
  • This article has been cited by:

    1. Nawal Kherbouche, Mohamed Helal, Abdennasser Chekroun, Abdelkader Lakmeche, Mathematical analysis and global dynamics for a time-delayed Chronic Myeloid Leukemia model with treatment, 2020, 15, 0973-5348, 68, 10.1051/mmnp/2020038
    2. Mostafa Adimy, Abdennasser Chekroun, Toshikazu Kuniya, Traveling waves of a differential-difference diffusive Kermack-McKendrick epidemic model with age-structured protection phase, 2022, 505, 0022247X, 125464, 10.1016/j.jmaa.2021.125464
    3. Mostafa Adimy, Abdennasser Chekroun, Toshikazu Kuniya, Global asymptotic stability for a distributed delay differential–difference system of a Kermack–McKendrick SIR model, 2022, 0003-6811, 1, 10.1080/00036811.2022.2075352
    4. Mostafa Adimy, Julien Molina, Laurent Pujo-Menjouet, Grégoire Ranson, Jianhong Wu, Forecasting the Effect of Pre-Exposure Prophylaxis (PrEP) on HIV Propagation with a System of Differential–Difference Equations with Delay, 2022, 10, 2227-7390, 4093, 10.3390/math10214093
    5. Fatima Zohra Hathout, Tarik Mohammed Touaoula, Salih Djilali, Mathematical analysis of a triple age dependent epidemiological model with including a protection strategy, 2022, 27, 1531-3492, 7409, 10.3934/dcdsb.2022048
    6. MOHAMMED AMINE MENOUER, NADIA GUL, SALIH DJILALI, ANWAR ZEB, ZAREEN A. KHAN, EFFECT OF TREATMENT AND PROTECTION MEASURES ON THE OUTBREAK OF INFECTIOUS DISEASE USING AN SIR EPIDEMIC MODEL WITH TWO DELAYS, DISCRETE AND DISTRIBUTED, 2022, 30, 0218-348X, 10.1142/S0218348X2240223X
    7. Abdelheq Mezouaghi, Salih Djillali, Anwar Zeb, Kottakkaran Sooppy Nisar, Global proprieties of a delayed epidemic model with partial susceptible protection, 2021, 19, 1551-0018, 209, 10.3934/mbe.2022011
    8. Mostafa Adimy, Abdennasser Chekroun, Bogdan Kazmierczak, Traveling waves for reaction-diffusion PDE coupled to difference equation with nonlocal dispersal term and time delay, 2022, 17, 0973-5348, 17, 10.1051/mmnp/2022021
    9. Tarik Mohammed Touaoula, Long time behaviour for a mixed reaction-diffusion-difference problem with distributed delay and non-local term, 2023, 0022247X, 127264, 10.1016/j.jmaa.2023.127264
    10. Fatima Zohra Hathout, Tarik Mohammed Touaoula, Salih Djilali, Efficiency of Protection in the Presence of Immigration Process for an Age-Structured Epidemiological Model, 2023, 185, 0167-8019, 10.1007/s10440-023-00572-8
    11. Mostafa Adimy, Abdennasser Chekroun, Grégoire Ranson, Laurent Pujo-Menjouet, Stability Analysis of a New Differential-Difference Model Applied to the Pre-exposure Prophylaxis (PrEP) Effect on the Spread of HIV, 2024, 23, 1575-5460, 10.1007/s12346-024-01093-x
    12. Mostafa Adimy, Abdennasser Chekroun, Toshikazu Kuniya, Hanene Meghelli, Global stability of a SEIR discrete delay differential‐difference system with protection phase, 2023, 46, 0170-4214, 17818, 10.1002/mma.9533
    13. Mostafa Adimy, Abdennasser Chekroun, Charlotte Dugourd-Camus, Hanene Meghelli, Global Asymptotic Stability of a Hybrid Differential–Difference System Describing SIR and SIS Epidemic Models with a Protection Phase and a Nonlinear Force of Infection, 2024, 23, 1575-5460, 10.1007/s12346-023-00891-z
    14. Mostafa Adimy, Abdennasser Chekroun, Charlotte Dugourd-Camus, Is it more effective to protect susceptible individuals or isolate infected people to prevent the spread of a SIR-type epidemic?, 2024, 0, 1531-3492, 0, 10.3934/dcdsb.2024172
    15. Hanene Meghelli, Abdennasser Chekroun, Mostafa Adimy, Traveling Waves in a Hybrid Reaction‐Diffusion‐Difference SEIR Epidemic Model With Nonlocal Transmission and Protection Phase With Delay, 2025, 0170-4214, 10.1002/mma.11110
  • Reader Comments
  • © 2019 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(6684) PDF downloads(674) Cited by(11)

Article outline

Figures and Tables

Figures(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog