![[Feedback]](/icons/banner/feedbackACT.gif){kind=icon}
![[For Subscribers]](/icons/banner/subscriberACT.gif){kind=icon}
![[Archive]](/icons/banner/archiveACT.gif){kind=icon}
![[Search]](/icons/banner/searchACT.gif){kind=icon}
![[Contents]](/icons/banner/tocACT.gif){kind=icon}
|
|
|
|
|
||
Vol. 14, Issue 11, 4695-4706, November 2003
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Center for Biomedical Imaging Technology, Department of Physiology, University of Connecticut Health Center, Farmington, Connecticut 06032
Submitted April 21, 2003;
Revised June 29, 2003;
Accepted July 3, 2003
Monitoring Editor: Mark Solomon
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
; Murray and Hunt, 1993). Cdc2 has an obligatory requirement for cyclin B binding (Cdc2-cyclin B is also known as M phase promoting factor or MPF). Cdc2-cyclin B activity is regulated by the activating phosphatase Cdc25 and the inactivating Wee1 family kinases, which act on two key phosphorylation sites (Solomon, 1993; Figure 1). In cells, Cdc2 activation is initiated either by changing the balance in favor of Cdc25 over Wee1 or by increasing the total Cdc2-cyclin B. The first way is thought to be more common. The second method has been demonstrated in experimental situations where exogenous cyclin B is added and binds to endogenous Cdc2, which is in excess (Swenson et al., 1986; Murray and Kirschner, 1989); it may also occur during rapid early embryonic divisions. A crucial feature of the Cdc2 activation pathway is that active Cdc2-cyclin B activates its activator Cdc25 and inactivates its inactivator Wee1 (Figure 1). This positive feedback is thought to result in all or none activation that does not flicker on and off once it is turned on.
|
Theoretically, the irreversible character of Cdc2 activation can be described in terms of bistability (Tyson et al., 2001). The fact that cells can be arrested in the inactive state and when activated remain in the active state even after the stimulus is removed implies that the system can be in at least two different stable states separated by some kind of a "barrier." Recently, experimental evidence for bistable behavior in Xenopus laevis egg extracts has been reported (Pomerening et al., 2003; Sha et al., 2003). This type of behavior can be elegantly explained by the possible bistability of a biochemical network (Thron, 1997; Ferrell and Xiong, 2001) illustrated qualitatively in Figure 2A for a hypothetical one-variable model. Usually, autoregulation implies a negative slope in the rate dependence (AB and CD portions of the graph in Figure 2A) that leads to a stable steady state, but the positive feedback mentioned above can change the sign of the slope in a certain concentration range (portion BC). As a result, the system may acquire two stable steady states, ss1 and ss3, separated by the unstable steady state ss2 (hence the term bistability). The system behavior in this case will be exactly like that in the cell cycle: if unperturbed, the system can stay infinitely long in the inactive steady state but "thrown over the unstable state" by a stimulus, it will be in the domain of attraction of the active steady state (interval [ss2, ss3] in Figure 2A) and will never return to the inactive state by itself.
|
Novak and Tyson (1993) have developed a detailed model of Cdc2 activation based on bistability of a biochemical pathway. An example of a switch-like behavior of the model triggered by inactivation of Wee1 is shown in Figure 2B. They found that an increase in total Cdc2-cyclin B could also trigger Cdc2 activation. The model thus reproduces key features of Cdc2 activation.
However, there are uncertainties, primarily due to the lack of experimental data necessary to make a complete model. The kinetic rate constants for the reactions shown in Figure 1 are not known, and the total concentrations (inactive plus active) of Cdc25, Wee1 and Cdc2-cyclin B are not established. These quantities, also known as parameter values, determine whether a biochemical network model is bistable.
There is another difficulty associated with parameter values. The total concentration is actually not a constant in cells, because the concentration can vary with time or in different regions. Likewise, the kinetic rate "constants" can vary due to differences in chemical potential. This leads us to the issue of robustness of the system: will the system behavior be retained in the face of biological variability? In the specific case of Cdc2 activation, over what range of total concentrations or rate constants will there be an irreversible activation of Cdc2? The two stable states of a bistable system are usually "antagonistic" in the sense that they can coexist only in a limited parameter range, whereas outside of this range only one of them prevails. This raises a concern whether bistability of a given biochemical pathway is robust enough to be the sole mechanism behind bio-switching.
Robustness of bistability has been discussed earlier in the context of cell-cycle models (Borisuk and Tyson, 1998) as well as in other contexts (Gardner et al., 2000). It has been recognized that mechanisms that result in steeply sigmoidal responses help expand parameter range where the model is bistable (Ferrell and Xiong, 2001). These mechanisms can be enzyme saturation or multistep phosphorylation in the positive feedback reactions [the former is suggested in Novak and Tyson (1993), the latter is assumed in Pomerening et al. (2003)], or cooperative binding in one or more steps of the loop (Gardner et al., 2000). It has also been noted that the presence of a buffer for MPF can, to an extent, improve robustness of bistability (Thron, 1997; see also "Analysis of the Effect of Independent Multiple Binding Sites" below). In this article, we demonstrate that bistability of a system can become essentially independent of kinetic parameters if multistability of a biochemical network couples with a phase equilibrium in the cell, so that the bio-switch involves a phase transition such as dissolution of precipitates.
The work presented here was stimulated by the discovery of aggregates of Cdc2-cyclin B in starfish oocytes (see Terasaki et al., in this issue). After applying the maturation-inducing hormone, the aggregates disappear at the time of MPF activation in a wave-like manner. In addition, experiments with fluorescence recovery after photobleaching (FRAP) indicate turnover of MPF between the aggregates and the cytoplasm.
We show how the Cdc2-cyclin B aggregates could have a central role in providing a decisive and irreversible transition. We first use a simplified model to find limits for bistability in the Cdc2 activation system in the absence of aggregates and then show that the robustness of a bio-switch can enhance significantly if the biochemical instability is coupled with a phase change in which the aggregates dissolve in the cytoplasm. Also, a simple spatial model is used to study the directionality of the Cdc2 activation wave.
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
To facilitate the analysis, we describe the pathway presented in Figure 1 with a simplified compartmental (nonspatial) model that can be solved analytically. Introducing notations as follows: m is the concentration of active MPF, c is the concentration of active Cdc25, and w is the concentration of active Myt1 (Myt1 is a representative of the "Wee1" family of proteins, hence the origin of the notation), and assuming, for simplicity, that the enzymes are unsaturated, the equations corresponding to the scheme of Figure 1 in the absence of aggregates can be written as
 | (1) |
)]. These parameters are rates of activation and inactivation of MPF (k and p), Cdc25 (ka and kb), and Myt1 (ke and kf), respectively. It is convenient for bistability analysis to introduce a nondimensional parameter 
= pwt/(kct), measuring the relative strengths of inhibitor versus activator, and two concentration-type parameters: mc = kb/ka, mw = ke/kf. The latter have the meaning of EC50 values, the equivalent concentrations of active MPF that under steady-state conditions would activate half of the total amounts of Cdc25 and Myt1, respectively. One side effect of the simplification made in system 1, compared with the original Tyson model (see Appendix 1 for comparison), is that the inactive state (m = 0, c = 0, w = wt) always exists, whereas in more realistic approaches, m and c are small but nonzero in the inactive state and this state does not exist beyond the stability threshold. However, our numerical tests show that this artifact does not significantly affect results of stability analysis, so the estimates of bistability constrains obtained below for the system 1 provide an accurate approximation for the systems with small nonzero m and c in the inactive state.1
It can be shown (for details of stability analysis, see Appendix 2), that the system described by Equation 1 may be bistable only under following condition:
 | (2) |
Moreover, given Eq. 2, the bistability exists only for the particular interval of mt:
 | (3) |
, only the inactive steady state is stable, whereas for 
the only stable state is the active one. The bistability constraints 2 and 3 are, in fact, more restrictive than may seem. For example, in the case of the equal inhibitor/activator "strengths", = 1, it follows from Eq. 2 that the system would be bistable only if the effect of MPF on its inhibitor were at least twice as strong as on its activator. If we further assume (for illustration only) that mc = 10 nM and mw = 2.5 nM, then the bistability exists only for mt between 8.7 and 10 nM (see Eqs. 3) (the Web site http://terasaki.uchc.edu/cdc2/provides an opportunity for numerical experimentation with this model). In general, as is shown in Figure 3, A and B, to obtain an "all or none" switch within a relatively large interval of mt, the inequality 2 should be replaced with ( + 1)mw << mc. Indeed, the width of the bistability interval, 
, is a decreasing function of ( +1)mw/mc and drops to less than half-maximum at ( + 1)mw/mc = 0.5 (Figure 3A). Also, tightening the interval 3 is accompanied by deviating from the all or none response (Figure 3B).
|
This analysis demonstrates that there is a relatively narrow range of parameters in which the Cdc2 activation system is bistable with all or none response. One could argue that the evolution process might tune a system to satisfy restrictive constraints; however, given the variability in biochemical kinetic parameters, it seems unlikely that the system would rely on fine-tuning in such a fundamental process as cell cycle regulation.
Total Cdc2-Cyclin B Can Act as an Effective Switch
In mathematical terms, there is an intrinsic, parameter-independent property of the system 1. This is the instability of the inactive state (with respect to activation) above a certain parameter threshold. The stability condition with respect to the total MPF reads as mt = mc (see Eq. 2.5 in Appendix 2). This implies that varying the total amount of Cdc2-cyclin B is an effective switch.
That total Cdc2-cyclin B can be an effective switch was shown by numerical simulations in the original Novak and Tyson (1993) model. In the cell, one way that total Cdc2-cyclin B could exceed a threshold is by synthesis of cyclin B. Cdc2 is usually present in excess, so that synthesis of cyclin B and its tight binding to Cdc2 will result in increasing total Cdc2-cyclin B. The rapid divisions of early embryonic development may be driven by synthesis of cyclin B rather than by regulation of Cdc2 activity by Cdc25 and Wee1 family kinases. There are also classic experiments in which increasing cyclin alone causes entry into M phase. In this work, expression of cyclin B by mRNA injection into frog oocytes caused maturation (Swenson et al., 1986), and in frog extracts depleted of endogenous mRNA, addition of cyclin mRNA can drive the extract through several cell cycles (Murray and Kirschner, 1989).
How Aggregates of Cdc2-Cyclin B Could Cause a Robust Switch
In starfish oocytes, protein synthesis inhibitors do not block maturation (Houk and Epel, 1974); this eliminates the possibility of increasing total Cdc2-cyclin B by synthesis. In this case, however, the aggregates of inactive MPF, being in equilibrium with inactive MPF in solution, can play a role of "stores," where inactive MPF is inaccessible to the activator but ready to be released into the cytosol after the hormone is applied. Because the hormone-triggered inactivation of Myt1 depletes the soluble inactive MPF, the turnover will lead to dissolution of the aggregates, thus increasing the total amount of MPF accessible for activation. To test this idea, we extend the compartmental model discussed in the previous section.
In the presence of the aggregates, the total amount of MPF in the solution, mt, becomes a variable, and our interpretation of the aggregates will determine the form of a governing equation for this variable. If the aggregates are the new insoluble phase of inactive MPF similar to a precipitate in a saturated solution, this equation comes from the kinetic theory of the first order phase transitions (Landau and Lifshitz, 1979). As shown in Appendix 3, the equation for mt in this case takes the form
 | (4) |
mt - m. The second term of the upper line originates from the energy cost of creating an interface between the phases and depends on the surface tension of an aggregate, so 
is a positive constant: = ac/am where ac is a critical size of nuclei at supersaturation of the order of saturation (Landau and Lifshitz, 1979) and am is the average size of the aggregate in the absence of the soluble phase (see Appendix 3). The kinetic constant kaggr is mainly determined by the diffusion coefficient of the soluble MPF and the number of aggregates per unit volume.
The system is now described by Eqs. 1 and 4. It follows from Eq. 4 that there are two families of the system steady states independent of each other: one, with mi 
ms, relates to the aggregates being in equilibrium with the solution, and the other, with mt = M, corresponds to the situation where the aggregates are dissolved. The overall system 1, 4 will be bistable if the inactive steady state is stable in the first family and there is a stable active steady state in the second family, whereas the bistability of the subsystem 1 is no longer required. As a consequence, constraints 2 and 3 are lifted. Instead, we get much more accommodating conditions, ms < mc and M > mc, which require that the concentration of the saturated soluble inactive MPF is below a certain threshold, whereas the concentration of all MPF is above it. Indeed, as shown in Appendix 3, the analysis of bistability of the full system can be reduced to that of stability of the inactive steady state in the subsystem 1, first for mt = mi ms (because m = 0 in the inactive state) and then for mt = M, and the constraints on ms and M shown above then follow from the stability threshold mt = mc. The fact that switching the system from inactive to active is no longer constrained by the requirements 2 and 3 indicates that robustness of bio-switching is enhanced dramatically when instability in the biochemical pathway is accompanied by the dissolution of the aggregates. This is illustrated in Figure 4 by results of numerical simulations based on Eqs. 1 and 4.
|
What will happen in the model developed above if cyclin B is overexpressed? It might seem that the system will remain bistable at arbitrary large M because the bistability constraints 3 are lifted. Indeed, in the pure saturated solutions, the level of saturation, ms, would remain independent of M and all the extra MPF would go to precipitates, so the system would never activate by overexpression.
However, in such a complex system as the cell, the additional MPF will likely be redistributed among its soluble and insoluble forms. In other words, ms may grow with M and at some point overcome the activation threshold. Also, it was implied that MPF in aggregates is inaccessible for the activator, which in reality might be true only to a certain degree. As aggregates grow, more MPF could be activated directly from them and enter the cytoplasm thus also contributing to the increase of the total concentration of the soluble MPF to the point where it exceeds the threshold. Therefore, the enhanced robustness of bistability does not necessarily rule out activation at sufficient overexpression of MPF.
Overall, we have shown that total Cdc2-cyclin B can act as a switch during starfish oocyte maturation. In this scenario, the aggregates behave like a precipitate, and their dissolution increases the total soluble Cdc2-cyclin B beyond a threshold, thereby activating the Cdc2 system. The robustness of the switch is enhanced significantly over the situation where all of the components are soluble.
Analysis of the Effect of Multiple Independent Binding Sites
In the previous section, we showed that if the aggregates behave like a precipitate, then there is robust switching. However, the aggregates could instead consist of a scaffold which has multiple independent binding sites for inactive Cdc2-cyclin B. We investigated how well this type of aggregate would act as a switch.
In this case, because the inactive MPF is partly immobilized due to binding to some kind of existing binding sites (an immobile buffer), Eq. 4 should be replaced with
 | (5) |
(M(1), M(2)) where M(2) = mc(1 + bt/(K + mc)) and 
where 
is the root of the equation 1 + btK/(K + x)2 = mw(mc - mw)/(x - mw)2 lying between mw and mc. Thus, the existing binding sites for inactive MPF do not enhance robustness of switching from the perspective of parameter independence, although because of binding to the buffer, the bistability interval for the full MPF concentration M shifts to larger values and expands. Also, because of that, the amount of active MPF in the activated state can increase significantly. However, in this scenario, the clumps of inactive MPF will never disappear entirely. The simulation results in Figure 5 illustrate activation of the system 1, 5 under various conditions. In particular, switching the system state from inactive to active fails if the bistability conditions are violated.
|
Other Characteristics of the Aggregates
In both of the aggregate models presented above, it is implied that activation occurs in solution, not in the aggregates. One alternative would be that all inactive MPF is immobilized, whereas the active MPF is soluble. We explored whether under these conditions, a simple spatial model described in the subsequent section (see also Appendix 1) can account for the observed Cdc2-cyclin B clumping.
Numerical simulations of the model run under the assumptions specified above immediately show that, to obtain the stable inactive steady state in this model, it is absolutely necessary that Myt1 colocalizes with inactive Cdc2-cyclin B. In fact, Cdc2-cyclin B clumps would be the consequence of Myt1 spatial aggregation. The problem with this explanation is that a turnover observed in FRAP experiments, in this interpretation, is the turnover between inactive and active forms of Cdc2-cyclin B governed by the same processes that activate Cdc2-cyclin B after applying the maturation-inducing hormone. Because the amount of active Cdc2-cyclin B in the inactive steady state should be very low, the clump turnover in the model is much slower compared with the activation time in contrast with the experimental findings that the time scale of turnover is comparable to that of activation (see Terasaki et al., in this issue). This suggests that the turnover, observed in FRAP experiments, is the exchange of inactive Cdc2-cyclin B between its soluble and insoluble forms, the assumption used in the preceding sections.
Another implication of both aggregate models is that of the total cellular Cdc2-cyclin B, a large fraction should be in the aggregates compared with in solution.
The results of FRAP experiments can be used to estimate the ratio of the amounts of Cdc2-cyclin B immobilized in the clumps and soluble in the cytoplasm. In those experiments, the characteristic time for recovery was 30 s. As shown in Appendix 3, this ratio can be evaluated as 3.33NDaav
/R3 where N is the number of the aggregates in the cell, D is the Cdc2-cyclin B diffusion coefficient, aav is the average radius of an aggregate, and R is the radius of the cell. With the parameter values determined from the experiments as N 
2000, D 10 µm2/s, aav 
1 µm, R 90 µm, and as above, the estimate yields the value of 2.7, which means that the amount of Cdc2-cyclin B immobilized in clumps is at least equal to if not greater than that of soluble in the cytoplasm and therefore significant.
Directionality of the MPF Activation Wave
Imaging of cyclin B-GFP in starfish oocytes (see Terasaki et al., in this issue) provides evidence that MPF activation occurs as a wave originating at the animal pole. Because it is well known that the activation of bistable reaction-diffusion systems can occur in the form of a "tidal wave" (Grindrod, 1996), we combine the kinetics of the system MPF-Cdc25-Myt1 (see diagram in Figure 1) with diffusion to analyze the spatial pattern of MPF activation. For this purpose, we use spatially uniform initial distributions for all molecular species involved in the model.
Mathematically, the system is described by a set of partial and ordinary differential equations (see Appendix 1). Using the Virtual Cell modeling framework (Schaff et al., 2000), we set up a model and solved it on the two-dimensional geometry mimicking the morphology of the cell. Because the maturation of the starfish oocytes is initiated by a maturation-inducing hormone that binds to the oocyte plasma membrane and enables a signaling cascade resulting, according to (Okumura et al., 2002), in production of the kinase Akt, we model this production as a flux originating from the membrane, by using an appropriate boundary condition. A set of parameters used in the simulations is close to the one suggested in Novak and Tyson (1993) and provides system bistability (see Appendix 1 for the model description and computation details).
Simulation results are presented in Figure 6. Interestingly, even with all the initial conditions and applied membrane flux of Akt being uniform, the wave of Cdc2-cyclin B activation (disappearance of inactive Cdc2-cyclin B in Figure 6A) initiates at the animal pole and propagates in the direction of the cell axis. Why it is so becomes clear from Figure 6B showing the Akt dynamics. Even though the flux of Akt production is uniform across the plasma membrane, the accumulation of Akt in the region of animal pole is initially much higher because of enhanced "local surface-to-volume ratio." Thus, our results indicate that the asymmetry in positioning of the nucleus inside the cell can be among the factors determining initiation and directionality of the Cdc2-cyclin B activation wave.
|
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
Because the key enzymatic interactions that govern Cdc2-cyclin B activation have an intrinsic threshold depending on the total amount of soluble Cdc2-cyclin B, the most convenient way of switching is to somehow regulate this amount, therefore the inactive MPF caged in aggregates is ideal for this purpose. At the same time, the stored MPF must "know" when it should be released, and the turnover between the aggregates and the soluble inactive MPF being in equilibrium with each other seems to be a perfect sensor for that: when, as a result of the hormone application, the cytoplasmic inactive MPF goes below saturation, the aggregates will dissolve thus increasing the total amount of MPF in the cytoplasm.
Finally, the new state must not be readily reversible, and here lies the difference between aggregation and buffering, which otherwise might look similar. In the scenario with the phase change, the new state with no aggregates is separated from the initial one by an energy barrier and, therefore, is at least metastable (characterized by the local minimum of free energy), the energy barrier being a cost of creating an interface between soluble and insoluble phases. In other words, to reverse the transition in this scenario, "binding sites" (the stable nuclei of the insoluble phase of sufficient size) must be created first, before binding could occur, whereas in the case of buffering the binding sites are readily available. Mathematically, this difference manifests itself in opposite signs of dmt/dt in Eqs. 4 and 5 as mt 
M.
The existence of two states with local minima of free energy (in one state, aggregates are in equilibrium with soluble MPF, in the other, aggregates are dissolved) separated by the energy barrier can also be interpreted in terms of bistability. Therefore, in this scenario, the system can be bistable even if the subsystem 1 is not and the parameter constraints 2 and 3 are not met. This significantly contributes to robustness of switching, because bistability of the system becomes largely insensitive to kinetic parameters as long as the saturating concentration of the soluble MPF, ms, is relatively low, whereas the total amount of MPF (soluble plus insoluble), M, is relatively high. However, in the case of overexpression, the additional MPF will likely be redistributed among its soluble and insoluble forms that eventually would lead to spontaneous activation.
In light of the hypotheses presented in this article, the most interesting issues that can be addressed experimentally concern Cdc2-cyclin B associations in the aggregates: are there other molecules present within the aggregates? What features of Cdc2-cyclin B could promote aggregation when inactive and prevent aggregation when active? Do the aggregates behave like a precipitate or a multiple binding site? On the theoretical side, the next obvious step for future development is the full spatial model that would include Cdc2-cyclin B aggregation.
It is possible that this analysis has more general relevance. There must be complex cellular mechanisms to regulate the amounts of all the different proteins, but it seems likely that these amounts fluctuate significantly over time and in different regions. In the Cdc2 activation pathway, we are proposing that aggregates of Cdc2-cyclin B relieve the cell of controlling amounts of the pathway components precisely. Perhaps aggregates or other mechanisms have the same role in other pathways whose response characteristics have a similar dependence on protein amounts.
|
APPENDIX 1 |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
) is combined with diffusion, it gives rise to a reaction-diffusion system that was used in spatial computer simulations of MPF activation. The model is a set of the parabolic-type partial differential equations,
 | (1.1) |
) and Cdc25, Myt1, and Akt act as unsaturated enzymes:
 | (1.2) |
:
 | (1.3) |
) and provide system bistability. The initial conditions are uniform and correspond to the inactive steady state in the absence of Akt. The values of cytoplasmic diffusion coefficients of MPF, Cdc25, and Akt have been reduced to 3 µm2/s to account for the hindering effect of yolk platelets, which are abundant in the cytoplasm (Terasaki et al., 2001).
|
Using the Virtual Cell computational framework (Schaff et al., 2000), the system 1.1–1.3 has been simulated numerically on two-dimensional geometry mimicking the morphology of a starfish oocyte (Figure 6). The computational domain has been reduced in half to 190 x 95 µm due to axial symmetry and contained 5000 mesh points. The simulated time is 25 min with the integration being performed with 0.05-min time step. The numerical error is estimated to be <1.5%.
|
APPENDIX 2 |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
 | (2.1) |
 | (2.2) |
 | (2.3) |
of the Jacobian of the system 1 evaluated at a steady state are negative. In other words, all of the roots of the equation
 | (2.4) |
 |
 | (2.5) |
1() = 2() with
 |
1 and 2 that two lesser roots, 1,2, of the equation 1() = 2() satisfy inequalities 1 < -ke -kfm < 0, 2 < -kam - kb < 0, whereas 3 is negative if 1(0) > 2(0). Using 2.1 and 2.3, one can show that the latter condition is equivalent to m2 > mw(mc - mt), which, with m from 2.2 and under the condition 2.5, yields
 | (2.6) |
|
APPENDIX 3 |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
). Consider a spherical aggregate with the radius of a, immersed in the solution of inactive MPF with concentration mi. The turnover between the aggregate and the solution is due to diffusion and the concentration of the soluble inactive MPF near the aggregate is that of saturation. This concentration, msa, is however different from the bulk saturation concentration, ms, because of the surface tension effects (Landau and Lifshitz, 1980):
 | (3.1) |
is the aggregate surface tension, v is the specific volume (the volume per molecule) in the aggregate, T is absolute temperature, and kB is the Boltzmann constant. The spatial distribution around the aggregate can be estimated from the diffusion equation under steady-state conditions,
 |
) = mi, mi(a) = msa has the form
 |
 | (3.2) |
 | (3.3) |
), the size distribution right after nucleation is linear and later tends to a Dirac delta-function: all aggregates gradually acquire approximately same size. In any case, 
and
 | (3.4) |
am3n/3 = Mv.
The second term in 3.3 is associated with the energy cost of creating an interface between the aggregate and solution (Landau and Lifshitz, 1980). Introducing the following notation: ac = 2v/(kBT), = ac/am, and kaggr = 4Dnam, we obtain from 3.3 and 3.4
 | (3.5) |
Based on the idea of diffusion-limited turnover, it is possible to estimate from the FRAP experiments the ratio of the amounts of MPF immobilized in the aggregates and soluble MPF in the cytoplasm. This ratio is 1.11N(aav/R)3(ma/mi) where N is the total number of aggregates in the cell, the concentration of MPF in the aggregate, 1/v, is denoted as ma, R is the cell radius, and we took into account that the nucleus takes 10% of the cell volume. The ratio of concentrations can be evaluated from the balance of fluxes in turnover; at steady state, the incoming flux from the solution, 4Daavmi, is balanced by the outgoing flux 
, where is the characteristic time of fluorescence recovery: 
. Substituting this back in the expression above, we obtain the result used in the main text.
Finally, we perform the steady-state stability analysis for the system defined by Eqs. 1 and 4 of the text. The eigenvalue equation for this system,
 |
M, to
 | (3.6) |
 |
ms, mt is stable with -kaggr < 0, and the problem reduces to the stability of the subsystem 1. For the steady states with mt = M, stability of mt follows directly from Eq. 4 rather than from Eq. 3.6 because of discontinuity of dmt/dt at mt = M, and the problem again reduces to the stability of the subsystem 1. |
APPENDIX 4 |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
 | (4.1) |
It is easy to verify that the stability condition of the inactive state of the system 1, 5, m = c = 0, w = wt, mt = (1/2)(M - K - bt + ((M - K - bt)2 + 4KM)1/2, as in the case of aggregation (Appendix 3), reduces to that of stability of the subsystem 1, mt < mc. This gives us the upper bound of the bistability interval for M:
 | (4.2) |
mt - m: m = 
1(mi) = 2(mi) with 1(x) = M - x - btx/(x + K) and 2(x) = mw(mc - x)/(x - mw). Schematic plots of these functions indicate that the solutions for the active state compatible with 4.2 are possible only if mw < mc. With decreasing M, the plot of 1(x) shifts down, and the solutions exist until the functions 1(x) and 2(x) become tangent to each other: d1(x)/dx = d2(x)/dx. This determines the lower bound of the bistability interval for 
,, where is the root of the equation
 |
(mw,mc). |
ACKNOWLEDGMENTS |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
|
Footnotes |
|---|
) and replaced kc with k(c + 
(ct - c)) where << 1. 
* Corresponding author. E-mail address: boris{at}neuron.uchc.edu.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTSDISCUSSIONAPPENDIX 1APPENDIX 2APPENDIX 3APPENDIX 4ACKNOWLEDGMENTSREFERENCES |
|---|
Ferrell, Jr., J.E., and Xiong, W. (2001). Bistability in cell signaling: how to make continuous processes discontinuous, and reversible processes irreversible. Chaos 11, 227-236.[CrossRef][Medline]
Gardner, T.S., Cantor, C.R., and Collins, J.J. (2000). Construction of a genetic toggle switch in Escherichia coli. Nature 403, 339-342.[CrossRef][Medline]
Grindrod, P. (1996). The Theory and Applications of Reaction Diffusion Equations. Patterns and Waves, Oxford, UK: Oxford University Press.
Houk, M.S., and Epel, D. (1974). Protein synthesis during hormonally induced meiotic maturation and fertilization in starfish oocytes. Dev. Biol. 40, 298-310.[CrossRef][Medline]
Landau, L.D., and Lifshitz, E.M. (1979). Physical Kinetics [in Russian], Nauka, Moscow, Secs. 99, 100.
Landau, L.D. and Lifshitz, E.M. (1980). Statistical Physics, Oxford, UK: Pergamon Press.
Marchenko, V.I. (1996). On the theory of fog. JETP Lett. 64, 66-69.[CrossRef]
Murray, A., and Hunt, T. (1993). The Cell Cycle. An Introduction, New York: W. H. Freeman Co.
Murray, A.W., and Kirschner, M.W. (1989). Cyclin synthesis drives the early embryonic cell cycle. Nature 339, 275-280.[CrossRef][Medline]
Nurse, P. (1990). Universal control mechanism regulating onset of M-phase. Nature 344, 503-508.[CrossRef][Medline]
Novak, B., and Tyson, J.J. (1993). Numerical analysis of a comprehensive model of M-phase control in Xenopus and intact embryos. J. Cell Sci. 106, 1153-1168.
Okumura, E., Fukuhara, T., Yoshida, H., Hanada, S., Kozutsumi, R., Mori, M., Tachibana, K., and Kishimoto, T. (2002). Akt inhibits Myt1 in the signaling pathway that leads to meiotic G2/M-phase transition. Nat. Cell Biol. 4, 111-115.[CrossRef][Medline]
Pomerening, J.R., Sontag, E.D., and Ferrell, Jr., J.E. 2003. Building a cell cycle oscillator: hysteresis and bistability in the activation of Cdc2. Nat. Cell Biol. 5, 346-351.[CrossRef][Medline]
Schaff, J.C., Slepchenko, B.M., and Loew, L.M. (2000). Physiological modeling with the Virtual Cell framework. In: Methods in Enzymology, vol. 321 (ed. M. Johnson), San Diego: Academic Press, 1-23.[Medline]
Sha, W., Moore, J., Chen, K., Lassaletta, A.D., Yi, C.-S., Tyson, J.J., and Sible, J.C. (2003). Hysteresis drives cell-cycle transitions in Xenopus laevis egg extracts. Proc. Natl. Acad. Sci. USA 100, 975-980.
Solomon, M.J. (1993). Activation of the various cyclin/cdc2 protein kinases. Curr. Opin. Cell Biol. 5, 180-186.[CrossRef][Medline]
Swenson, K.I., Farrell, K.M., and Ruderman, J.V. (1986). The clam embryo protein cyclin A induces entry into M phase and the resumption of meiosis in Xenopus oocytes. Cell 47, 861-870.[CrossRef][Medline]
Terasaki, M., Campagnola, P., Rolls, M.M., Stein, P. A., Ellenberg, J., Hinkle, B., and Slepchenko, B. (2001). A new model for nuclear envelope breakdown. Mol. Biol. Cell 12, 503-510.
Thron, C.D. (1997). Bistable biochemical switching and the control of the events of cell cycle. Oncogene 25, 317-325.
Tyson, J.J., Chen K., and Novak, B. (2001). Network dynamics and cell physiology. Nat. Rev. Mol. Cell. Biol. 2, 908-916.[CrossRef][Medline]
This article has been cited by other articles:
|
|
 |
|
  J.-Q. Wu, V. Sirotkin, D. R. Kovar, M. Lord, C. C. Beltzner, J. R. Kuhn, and T. D. Pollard Assembly of the cytokinetic contractile ring from a broad band of nodes in fission yeast J. Cell Biol., July 31, 2006; 174(3): 391 - 402. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 P. Marangos and J. Carroll The dynamics of cyclin B1 distribution during meiosis I in mouse oocytes Reproduction, August 1, 2004; 128(2): 153 - 162. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||
