Zhuan Liu and
Kunkun Guo*
College of Materials Science and Engineering, Hunan University, Changsha, 410082, China. E-mail: kunkunguo@hnu.edu.cn
First published on 17th October 2014
FtsZ filaments play a central role in bacterial cell division. A theoretical framework is developed by combining a phase field model for rod-shaped cells with a kinetic description for FtsZ ring maintenance, in order to investigate cell morphodynamics during bacterial cell division. The cell division time and cell shape are collectively determined by the curvature elastic energy of the cell membrane/wall and the constriction force generated by FtsZ rings. The dependence of cell morphodynamics during cell division on different initial states of rod-shaped cells and FtsZ rings, such as the aspect ratio and the FtsZ concentration in cells, ZT, is extensively studied. The obtained results with the measured experimental parameters are found to be well comparable to the observed results physiologically. Likewise, it is found that the quasi-steady state of FtsZ rings accords with the theoretical result derived from the kinetic description of FtsZ rings. In addition, the morphological phase diagram is presented as a function of the FtsZ concentration in cells, ZT, and the aspect ratio of rod-shaped cells, given that rod-shaped cells with high ZTs and/or aspect ratios tend to divide. It will be straightforward to extend the theoretical framework to other complicated systems relevant to filamentous proteins and cells, for example, maintaining structural integrity, serving as a template for cell growth, cell adhesion and motility, and mechanical signal transduction.
The FtsZ ring has a dynamic structure, where filaments assemble into the ring and disassemble from it.13 The mechanism of its dynamic behavior can not be completely explained by the dynamic instability of microtubules or treadmilling of other filamentous proteins.14,15 Of course, the GDP-bound FtsZ monomers that have been depolymerized from the FtsZ ring can be rapidly exchanged with GTP-bound ones in the cytoplasm for subsequent polymerizations, so that the cytosolic filaments consist entirely of GTP-bound monomers. Likewise, FtsZ monomers are observed to be rapidly exchanged between FtsZ rings and the cytoplasmic pool, and a half time of turnover is measured experimentally as about 10 s.16
One mechanism for force generation is proposed in several theoretical models, arising from curvature transitions induced by GTP hydrolysis.5,7–9,17–20 However, the conformation switch is observed experimentally to be independent of the nucleotide binding states of the FstZ monomers. At the same time, crystal structures of single FtsZ monomers reveal the structure of the GTP-bound monomer to be similar to that of the GDP-bound one.17 Subsequently, the effects of nucleotide binding on the structure of a FtsZ dimer are extensively studied by molecular dynamic simulations, supposing that the FtsZ dimer structure depends on its nucleotide binding state.18 As a whole, in the presence of GTP hydrolysis, FtsZ filaments are expected to possess three different curved conformations.5,7–9,19,20 The first is a straight conformation with a curvature radius larger than 0.5 μm when the FtsZ protein binds GTP molecules. The intermediate curved conformation corresponds to a 2.5 degree bend between monomers, producing a curvature radius of about 200 nm. Finally, the highly-curved conformation (miniring) contains 16 monomers bound to GDP molecules with a 23 degree bend at each interface, and its intrinsic curvature radius ranges from 12.5 to 100 nm. Therefore, the coupling of GTP hydrolysis at each of the FtsZ monomers within the FtsZ filament could cause a transition from a straight conformation to a curved one, by which the FtsZ filament would provide the force to constrict the cell wall/membrane.
A completely different mechanism for force generation has been suggested in other theoretical models.21–24 These models propose that the FtsZ ring consists of several short filaments, and these short filaments can interact with each other via lateral and longitudinal bonds. These lateral and longitudinal bonds provide a negative (favorable) free energy. Meanwhile, a constriction force would be generated as the cooperativity of the polymerization, the condensation and bundling of the FtsZ filaments proceed, together with the increased number of lateral and longitudinal bonds.
So far, research focused on the force generation of FtsZ rings has mainly disregarded the shape transformations of cells. The coupling of FtsZ rings and cell shape is expressed in such a way that the constriction force is a free parameter, independent of the dynamic turnover of FtsZ rings.25 The phase field method has a reputation to be very general and applicable to complex microstructural phenomena, and it has been extended to study vesicle shape dynamics recently,26–29 in particular, cell morphodynamics due to phase transitions between cross-linked actin filaments and bundles.30,31 In the present study, one theoretical framework is firstly developed to study the coupling of cells and FtsZ ring systems, and combines a kinetic description of FtsZ ring maintenance with phase field dynamics for cell morphodynamics. The greatest advantage of the theoretical framework is to capture the collective behaviors of shape transformations of cells and the constriction forces generated by dynamic FtsZ rings. Furthermore, the theoretical framework presented here is easily extended to investigate other complicated systems relevant to filamentous proteins and cells.
This paper is organized as follows. Section II contains a detailed description of the theoretical framework by combination of the phase field model for cells and one kinetic model for FtsZ ring turnover, together with the numerical calculation method for the combined theoretical model. In Section III, a typical dynamic process during bacterial cell division is presented, and the steady states of cells are discussed as functions of different initial states of cells and FtsZ rings, such as the total concentration of FtsZ in a cell, the total number of FtsZ monomers within the ring, the total number of FtsZ–GDP monomers within the ring, the initial mean length of FtsZ filaments within the ring, λ, and the aspect ratio of cells. Additionally, a morphological phase diagram is presented as a function of the total concentration of FtsZ in the cells and the aspect ratio of cells. Finally, a brief summary and outlook are given in Section IV.
A typical bacterium, such as E. coli, is approximately assumed as a two-dimensional rod-shaped cell with a fixed surface area A0, which is equivalent to the volume in three dimensions, and the dynamic turnover of short FtsZ filaments is considered to couple with GTP hydrolysis, see Fig. 1. The FtsZ ring that consists of several short FtsZ filaments, as expected, is always distributed in the midcell. In order to distinguish the interior and exterior of the rod-shaped cell, an auxiliary phase field, ϕ, is introduced, where ϕ takes on the value of 1 in the interior of the cell wall but ϕ = 0 in the cell exterior. This field varies abruptly in the diffusive interface between two limited values, ϕ = 1 and ϕ = 0. The width of this interface, ε, is used to describe the thickness of the cell wall.
Therefore, the dynamic evolution equation for the phase field is given by26,27,30
![]() | (1) |
In the present study, the shape of the cell membrane or cell wall can be determined by the cooperation of various forces, including the surface tension, the bending force, and the pressure that constrains the cell area of vesicles, as can the volume in three dimensions. In addition, the perimeter of the cell wall is observed experimentally not to be fixed,33 allowing proteins or lipid molecules to enter or escape during either cell wall expansion or cell division. We also consider the radial and constriction forces generated from FtsZ rings in the midcell, and the effective friction caused by cell division and cell motility. At first, the surface energy that is proportional to the cell’s perimeter L, is able to be implemented in the phase field formulations as follows28,29
![]() | (2) |
![]() | (3) |
Here, this area density can be converted into a line density with F′dr = Fte ε|∇ϕ|2 dr.30 Therefore, the surface tension force with a line density is defined as
![]() | (4) |
The bending energy Hbe of the cell wall32 is written in the phase field formulation as
![]() | (5) |
![]() | (6) |
The cell area in two dimensions, A = ∫ϕdr, which is equivalent to the volume in three dimensions, is observed experimentally33 to be conserved during cell division, indicating that the perimeter is not highly conserved. Therefore, a constraint term is required to guarantee the cell area and can be expressed in the form of phase field.
![]() | (7) |
![]() | (8) |
The coupling of FtsZ rings with cells provides a retraction force to the cell wall/membrane, thereby leading to cell division. In order to simply describe the mechanical properties of the FtsZ ring, a stable and smooth FtsZ ring is assumed to form in the midcell. This assumption of a smooth ring is justified by the mechanical averaging across the thickness of the ring due to the uniform distribution of hydrolyzed and unhydrolyzed monomers. The mechanical energy Hz stored in the FtsZ ring as a function of its radius and the size of FtsZ monomers is assumed to obey linear elasticity, given by34
![]() | (9) |
![]() | (10) |
In order to achieve the mechanical energy, Hz, of the FtsZ-ring given by eqn (9), the total number of FtsZ monomers and GDP-bound FtsZ monomers within the ring have to be obtained in advance. This involves a kinetic description of FtsZ ring turnover, including the incorporation of free short filaments and FtsZ monomers into the ring, hydrolysis of GTP-bound monomers within the ring, disassembly of both GTP-bound and GDP-bound monomers at the tip of the filament within the ring, as well as the mechanical characterization of the FtsZ ring in terms of force generation. Compared to the previous kinetic model proposed by Cytrynbaum,13 the contributions of GTP-bound monomers within the ring to disassembly and to force generation are considered. Two main physical quantities are introduced, and these are the length distribution of filaments in the ring, p(l, t), and the number of hydrolyzed monomers, SD(t) in the ring, where l is the number of monomers within the filament at time t. Subsequently, the total number of FtsZ monomers, S, and filament tips, F, in the ring can be derived by and
, respectively. The cytosolic concentration of FtsZ monomers can be obtained from ZT − S/(NcA), where ZT is the total concentration of FtsZ monomers in the cell, and the factor NcA is used to transform the number of molecules to concentration in micromoles.
The cytosolic filaments are experimentally observed to satisfy an exponential distribution when the cytosolic concentration is above the critical concentration.10 As concluded from polymer dynamics, independent of cooperative versus isodesmic assembly or fragmentation and annealing, filament lengths are also found to obey a quasi-steady exponential distribution with a mean length, λ.13 Therefore, the lengths of cytosolic filaments are reasonably assumed to possess an exponential distribution with a mean length of λ = 30 monomers (typical length), and the amplitude of the distribution can be determined by the cytosolic concentration. The incorporation rate of filaments with a length, l, is proportional to the cytosolic concentration of the same filament length. In addition, the GDP-bound monomers depolymerized from the ring can be rapidly exchanged with GTP-bound monomers, thereby giving cytosolic filaments that consist entirely of GTP-bound monomers. As previously presented,14,15 the depolymerization of one monomer is governed by not only its position within a filament but also its nucleotide state. κToff and κDoff are defined as the depolymerization rates of GTP-bound and GDP-bound FtsZ monomers at the filament tips, respectively. However, the GDP-bound monomers at the filament tips are experimentally found to disassemble more rapidly than GTP-bound ones.4 The fraction of GDP-bound monomers in the ring is fD = SD/S.
Filaments are assumed to be randomly distributed throughout the ring. Changes in the length distribution p(l, t) arise from the incorporation of both FtsZ filaments and free monomers in the cytosol, and the disassociation of monomers from filament tips dependent on their nucleotide states. Therefore, the filament length distribution p(l, t) in the ring obeys
![]() | (11) |
The time evolution of the number of GDP-bound FtsZ monomers in the ring, SD, is given by
![]() | (12) |
The first term on the right side in eqn (12) accounts for the hydrolysis, in that any GTP-bound monomers in the FtsZ ring can hydrolyze in a stochastic manner, independent of the state of adjacent monomers. The second term on the right side in eqn (12) accounts for the disassembly of GDP-bound monomers at all filament tips within the ring F. The total number of filament tips F is obtained by integrating the length distribution of filaments within the ring, as presented in the preceding paragraph. κhy is the rate of GTP-bound monomers hydrolyzed into GDP-bound ones. By combining eqn (11) and (12), the quasi-steady state for this subsystem (FtsZ ring) is given by the following
![]() | (13) |
![]() | (14) |
![]() | (15) |
![]() | (16) |
During cell division, an effective friction between cell wall/membrane and fluid environment is proportional to the local speed, u, obeying Ffr = −τu. At the quasi-steady state, the total force approaches zero, that is Fte + Fbe + Far + Fz + Ffr = 0. Therefore, the evolution equation of the phase field model in the presence of FtsZ rings according to eqn (1)–(10) can be derived as follows
![]() | (17) |
Due to the presence of MinCDE proteins in cells, the spatial distribution of the FtsZ rings is regulated to the middle of cells.1,35 An FtsZ ring, which generally consists laterally of 6–10 filaments, localizes in the midcell along the horizontal axis at a wide range of 25 nm.11 Therefore, the force generated by FtsZ rings in other positions is very small, and can be approximately neglected. The evolution equation of the phase field model in the absence of FtsZ rings is given by
![]() | (18) |
The fourth-order nonlinear partial differential equations (17) and (18) are solved using an alternating direction-implicit scheme and a second order backward differentiation formula, while the partial differential equations (11) and (12) are explicitly solved using a forward Euler differential scheme. These position-dependent equations are solved on a 400 × 200 rectangle with the box size of 40r0 × 20r0, where the unit length r0 is about 0.1 μm. Therefore, the average length of the sides of the box is about 3 μm, which is available to in vivo bacteria.1,11 The simulation time step of these evolution equations is Δt = 4.0 × 10−5 s. A typical simulation starts with a stationary rod-shaped cell with surface area A0 = 2.78 μm2 and an aspect ratio equal to 3.0, which are also comparable with in vivo bacteria.1,11 The diffusive interface that describes the thickness of the cell wall is chosen as one unit length, r0 = 100 nm, which is comparable to the experimentally observed thickness of the cell wall in the range of 20 to 50 nm.11 Compared with the experimental result for the lateral size of the FtsZ rings,11 the dynamic FtsZ ring is assumed to be distributed in the middle of a rod-shaped cell, and has a width of 0.3r0 (about 30 nm). Therefore, the evolution equation of the phase field in the presence of the force generated by FtsZ rings, eqn (17), is only solved on three grid sizes (about 30 nm) in the middle of a rod-shaped cell. The lengths of the FtsZ filaments are presented to satisfy a quasi-steady exponential distribution: 1/λe−l/λ, see eqn (16), where λ is the mean length of the FtsZ filaments. Therefore, the total number of the monomers within the ring at the initial state is obtained by Sinitial = ∫lp(l, 0)dl. Likewise, the filament length distribution within the ring is determined by p(l, 0) = Sinitial/λ2e−l/λ in the initial state. The length distribution and the total number of GDP-bound FtsZ monomers within the ring for the next time step are able to be obtained according to the evolution of eqn (11) and (12). The resulting values of S and SD are inserted into the evolution eqn (17) and (18) of the phase field model to produce a new shape of a cell. These equations are solved iteratively to produce the new cell shape, as well as the updated distribution of FtsZ filaments. The steps are finished when a steady state is eventually reached where the cell has a stationary shape and the FtsZ filaments possess a stable distribution. The parameters used in the numerical calculations, which are all comparable to the experimental values,13,31 are provided in Table 1.
Description | Value | |
---|---|---|
γ | Surface tension | 50 pN μm−1 |
κ | Bending rigidity | 0.5 pN μm |
τ | Friction coefficient | 5 × 103 pN s μm−3 |
B | Bending modulus of FtsZ filament | 1.2 × 10−2 pN μm2 |
ε | Boundary width | 0.1 μm |
δ | Size of FtsZ monomer | 4 nm |
Γ | Langrange multiplier | 0.1 |
κT | Intrinsic curvature of GTP-bound FtsZ | 0.8 μm−1 |
κD | Intrinsic curvature of GDP-bound FtsZ | 80 μm−1 |
Nc | Conversion factor | 600 μM−1 μm−3 |
κin | Membrane association rate | 0.0021 μM−1 μm−1 s−1 |
κToff | Depolymerization rate of GTP-bound FtsZ | 0.1 s−1 |
κDoff | Depolymerization rate of GDP-bound FtsZ | 3.5 s−1 |
κhy | GTP hydrolysis rate | 0.03 s−1 |
MA | Area constraint | 1.0 pN μm−2 |
λ | Average number of monomers in a filament | 30 monomers |
ZT | Total concentration of FtsZ monomers in cell | 13 μM |
![]() | ||
Fig. 3 The effect of the mean length for FtsZ filaments, λ: (a) constriction forces generated by FtsZ rings at the initial and steady states, and contraction rates in the inset; (b) S, SD and fD at the steady state, the red lines are obtained from the theoretical values, eqn (13)–(15); (c) the radius of the constriction ring at the steady state, cell shapes at the steady state given in the insets; (d) length distribution with λ = 10 and λ = 30 at the steady state. |
Cell wall expansion is experimentally observed to be mainly contributed to by cell-wall biosynthesis, such as peptidoglycan synthesis proteins.36,37 Meanwhile, the turgor pressure is found not to be required for cell-wall biosynthesis, but does play a simple role, whereby it stretches recently assembled, unextended peptidoglycan.36,37 Once the FtsZ ring is formed and is anchored to cell wall/membrane, other proteins, such as the septum synthesis proteins, are recruited to the ring site, and cell wall growth and remodeling can proceed. In addition, many bacteria, for example, E. coli, generally elongate with both an exponential growth and a constant diameter. Therefore, the dynamic processes of cell division and steady states of cells are extensively investigated in response to the preset aspect ratio for rod-shaped cells. The results are given in Fig. 4. The cell area of rod-shaped cells, as expected, is altered upon the change in the aspect ratio. However, the initial concentration of FtsZ monomers in cells is chosen as a constant, 13 μM, the fraction of FtsZ monomers within the ring contains 30% of the total number of FtsZ monomers, and the fraction of hydrolyzed monomers in the ring, fD, is 0.4. As shown in Fig. 4(c), when the aspect ratio for rod-shaped cells shifts from 1.6 to 3.8, the radius at the midcell in the steady state, R, at first increases from 0.5 μm to close to 1.0 μm, where the shape is almost spherical, and then decreases gradually until cell division finishes. Typical shapes of cells at the initial and steady states are given in the inset of Fig. 4(c) and their aspect ratios are 2.4 and 3.8. It is noted that one rod-shaped cell would be inclined to divide at the ring site when it has a high aspect ratio in the initial state. That is, a rod-shaped cell that has a higher aspect ratio, would tend to undergo cell division more easily. The result is in good agreement with experiments that show that cell division of prokaryotes happens as a result of cell wall expansion with a constant radius.37 In the present study, the horizontal length of rod-shaped cells is kept constant, while the longitudinal length has to adjust to satisfy the preset aspect ratio. Therefore, the radius of rod-shaped cells, R, would reduce along with the increase in the aspect ratio, as well as the fact that the areas of rod-shaped cells would respond to the transformation of the aspect ratios. If the area of rod-shaped cells with different aspect ratios is fixed, the radius of a rod-shaped cell would decrease along with the increase in the aspect ratio. The initial constriction force generated by FtsZ rings would be augmented with respect to the contractile radius, when S and SD are kept unchanged in the initial state, see eqn (10). As shown in Fig. 4(a), the initial constriction force is confirmed to increase as a function of the aspect ratio. Additionally, it can be inferred from eqn (13) and (14) that the total number of FtsZ and GDP-bound FtsZ monomers within the ring would be decreased along with the reduction of R in the steady state, which is confirmed in Fig. 4(b). However, the constriction force at the steady state would depend on the competition between its radius and the number of FtsZ monomers within the ring. Therefore, the forces at the steady state shown in Fig. 4(a) exhibit evident fluctuations in response to the aspect ratio. As a whole, cell division depends on the aspect ratio for rod-shaped cells, as well as cell wall expansion in vivo.
Fig. 5 illustrates the dependence of the dynamic process in cell division on the total concentration of FtsZ monomers in cells, ZT, where the number of FtsZ monomers in the ring, S, is 30% of the total number of FtsZ monomers in the cell, and the fraction of hydrolyzed monomers in the ring, fD, is 40%. It is demonstrated experimentally10,11 that at saturated FtsZ concentrations, overlapping FtsZ filaments remain, allowing additional rings to form, whereas at low FtsZ concentrations of about 3–5 μM, a broad distribution of membrane-bound FtsZ monomers is observed in the midcell, but no FtsZ rings form. A critical concentration for ring formation is required, given that only a single ring forms when the concentration of FtsZ monomers is about 8–10 μM, and multiple rings form when the concentration is above 10 μM. Therefore, not only the formation of FtsZ rings but also the constriction force are expected to depend on the concentration of FtsZ monomers in cells. Fig. 5(a) presents the dependences of the constriction forces at the initial state and steady state on ZT. It can be seen that the constriction forces, regardless of the initial and final states, Fz,initial and Fz,final, increase along with the increase of ZT. When ZT is less than 12 μM, the constriction force generated by the FtsZ rings is found not to be large enough to overcome the curvature elastic energy for cells, therefore, the rod-shaped cell would eventually become spherical in shape (energy benefit). As ZT further increases, the cell tends to retract at the midcell, and the radius at the midcell becomes less. When ZT increases up to 13 μM, the constriction force continuously maintains a high level, and then cell division proceeds. The radius at the midcell is presented in Fig. 5(c) as a function of ZT, and the variation trend is found to be in accordance with the previous discussions. Fig. 5(b) gives the dependences of the total number of FtsZ and GDP-bound monomers, S and SD, on ZT. It can be seen that at ZT < 13 μM, S and SD increase as a function of ZT, but at ZT > 13 μM, both S and SD gradually decline. Likewise, it is found that at ZT > 13 μM, cell division starts to happen, and the radius at the midcell rapidly reduces, thereby leading to the decrease in both S and SD. In addition, it is found in the inset of Fig. 5(b) that the fraction of hydrolyzed monomers in the ring, fD, keeps fluctuating as a function of ZT, consistent with the theoretical result given in eqn (15). At the same time, as shown in Fig. 5(b), the numerical results of S and SD (symbols) are found to obey the relations of eqn (13) and (14) (red lines).
![]() | ||
Fig. 5 The effect of total concentration of FtsZ monomers in cells, ZT: (a) constriction forces generated by FtsZ rings at the initial and steady states, and contraction rates in the inset; (b) S, SD and fD for FtsZ rings at the steady state, the red lines are obtained from the theoretical values, eqn (13)–(15); (c) the radius of FtsZ rings at the steady state, a cell shape at the steady state given in the inset. |
Fig. 6 presents the dependence of cell morphodynamics on the initial number of FtsZ monomers within the FtsZ rings, Sinitial, where the initial fraction of hydrolyzed monomers in the ring, fD, is still kept at 40%. It can be seen from Fig. 6(a) that the initial constriction force generated by FtsZ rings increases linearly as a function of Sinitial, as concluded from eqn (10). However, we find that at Sinitial < 6500 monomers, the constriction force, S, SD and R in the steady state are independent of Sinitial, where the constriction force of about 8.4 pN at the steady state is far less than the force required for cell division. At Sinitial > 6500 monomers, the initial force (about 28 pN) is larger than the critical force for cell division, the rod-shaped cell gradually shrinks at the midcell, and then cell division begins. In addition, the fraction of hydrolyzed monomers in the ring, fD, is found to be close to 0.162 ± 0.002 and independent of Sinitial, see the inset of Fig. 6(b).
Fig. 7 gives the effect of the initial number of GDP-bound FtsZ monomers within the ring, SD,initial, on cell morphodynamics, where ZT and Sinitial are chosen as 13 μM and 6500 monomers, respectively. In Fig. 7(a), the initial constriction force generated by the FtsZ ring, Fz,initial, is found to increase linearly along with the enhancement of SD,initial. This trend obeys the relation of eqn (10). Likewise, it is also found that at SD,initial < 2600 monomer, the constriction force, S, SD and R in the steady state remain almost unchanged, and are independent of SD,initial, where the constriction force is about 8.4 pN and far less than the critical force (∼28 pN) for cell division. At SD,initial > 2600 monomers, the initial force is about 28 pN, close to the critical force for cell division. The radius at the midcell becomes much less, and then cell division is under way. In addition, the fraction of hydrolyzed monomers in the ring, fD, is found to be approximately 0.162 ± 0.002 regardless of the alteration in SD,initial, see the inset of Fig. 7(b).
As a result, the constriction force generated by the FtsZ rings, whether in the initial or in the steady state, is required to be above the critical force of 28 pN in order to arrive at cell division, regardless of the changes in the total concentration of FtsZ monomers in cells, the mean length of FtsZ filaments, and the aspect ratio for rod-shaped cells. This critical constriction force has the same magnitude as in the previous result, indicating that cell division can succeed for a wide range of FtsZ rings with a force between 8 pN and 80 pN.25 Additionally, it is found that when the initial constriction force generated is less than the critical force, the constriction forces, S, SD, fD and the radius at the midcell in the steady state are independent of the initial number of FtsZ monomers and GDP-bound FtsZ monomers in the FtsZ rings, Sinitial and SD,initial. However, except for the mean length of FtsZ filaments, the initial constriction forces, Fz,initial, are found to increase along with the enhancement of ZT, the aspect ratio for rod-shaped cells, Sinitial and SD,initial. In addition, it is noted that S and SD to some extent start to decrease at the beginning of such a reduction in the radius at the midcell. At the same time, the fraction of hydrolyzed monomers in the steady state, fD is observed to be dependent on the mean length of FtsZ filaments, λ, but independent of ZT, the aspect ratio for rod-shaped cells, Sinitial and SD,initial. Whereas, S and SD in the steady state are observed to depend on λ, ZT and radius at the midcell, R. As a whole, these changes in S, SD and fD obtained in the numerical calculations are found to agree well with the theoretical results derived from eqn (13)–(15), which describe the quasi-steady state of the FtsZ ring subsystem.
![]() | ||
Fig. 8 The phase diagram as a function of aspect ratio for rod-shaped cells and total concentrations of FtsZ in cells, ZT. |
The constriction force generated by FtsZ rings determines the contractile direction, velocity, as well as cell shape. The constriction force presented here are altered along with the dynamic process during cell division, as well as the changes in cell shapes and dynamic turnover of FtsZ rings. We find that with reasonable parameters given in Table 1, a rod-shaped cell retracts at the midcell when the constriction forces are above the critical force of about 28 pN that has the same magnitude as in the previous report.25 Cell shapes and cell division times obtained in our calculations are all comparable to physiologically observed results.11 FtsZ filaments within the ring are also found to obey an exponential distribution that has its characteristic length close to the preset mean length of FtsZ filaments, λ. Likewise, cell morphodynamics during cell division are extensively investigated as functions of different initial states of rod-shaped cells and FtsZ rings, such as the aspect ratio, the mean length of FtsZ filaments, the concentration of FtsZ monomers in cells, ZT, and the total number of FtsZ and GDP-bound FtsZ monomers within the ring. It is found that with the increase in the aspect ratio, the mean length of FtsZ filaments, and ZT in the initial state, rod-shaped cells contract at the midcell and eventually divide when the constriction force is above the force required for cell division. However, when the constriction force is below the critical force, the cell shape eventually becomes close to spherical, and is independent of the total number of FtsZ and GDP-bound FtsZ monomers within the ring in the initial state. It is noted that the total number of FtsZ and GDP-bound FtsZ monomers within the ring, S and SD, start to reduce to some extent at the beginning of the reduction of the radius at the midcell. In addition, the quasi-steady state for the dynamic turnover of FtsZ rings is derived, thereby giving the theoretical expressions for Sqss, SDqss and fDqss, see eqn (13)–(15). The obtained results for S, SD and fD in the numerical calculations are found to be consistent with the theoretical results, and fD approaches the experimental value of about 20%.11
In order to explicitly elucidate the relation of cell division and some parameters, the morphological phase diagram is presented as a function of the FtsZ concentration in cells, ZT, and the aspect ratio of rod-shaped cells, given that cell division tends to be observed in rod-shaped cells with high ZT and/or high aspect ratios. As experimentally observed,10,11 only a single ring forms at ZT < 10 μM, but at this concentration, the constriction force generated by a single ring is not high enough to initiate cell division. Likewise, one rod-shaped cell that possesses a higher aspect ratio would easily undergo cell division, in good agreement with the experimental result that cell division of prokaryotes happens due to cell wall expansion with a constant radius.37 Therefore, the theoretical framework presented here can effectively predict cell morphodynamics during division with the coupling of dynamic FtsZ rings and cell morphology.
In the present work, the dependence of cell morphodynamics during division on different initial states of cells and FtsZ rings is the focus. A more detailed discussion relevant to various kinetic rates within the FtsZ ring and division times will be presented in future work. The basic theoretical framework developed here, which has the evident advantage to integrate cell shape with dynamic FtsZ rings, can be extended to study division in eukaryotic cells, such as budding and fission in yeast and plant cells, which can provide valuable insights into critical cellular functions. In future research, it is also of great interest to develop other space-dependent theoretical frameworks aimed to study the spatial structure and dynamic turnover of FtsZ rings.
This journal is © The Royal Society of Chemistry 2014 |