A skeleton model for the network of cyclin-dependent kinases driving the mammalian cell cycle

A comprehensive review of computational cell cycle models in guiding cancer treatment strategies

This article reviews the current knowledge and recent advancements in computational modeling of the cell cycle. It offers a comparative analysis of various modeling paradigms, highlighting their unique strengths, limitations, and applications. Specifically, the article compares deterministic and stochastic models, single-cell versus population models, and mechanistic versus abstract models. This detailed analysis helps determine the most suitable modeling framework for various research needs. Additionally, the discussion extends to the utilization of these computational models to illuminate cell cycle dynamics, with a particular focus on cell cycle viability, crosstalk with signaling pathways, tumor microenvironment, DNA replication, and repair mechanisms, underscoring their critical roles in tumor progression and the optimization of cancer therapies. By applying these models to crucial aspects of cancer therapy planning for better outcomes, including drug efficacy quantification, drug discovery, drug resistance analysis, and dose optimization, the review highlights the significant potential of computational insights in enhancing the precision and effectiveness of cancer treatments. This emphasis on the intricate relationship between computational modeling and therapeutic strategy development underscores the pivotal role of advanced modeling techniques in navigating the complexities of cell cycle dynamics and their implications for cancer therapy.

Generalized Michaelis-Menten rate law with time-varying molecular concentrations

The Michaelis-Menten (MM) rate law has been the dominant paradigm of modeling biochemical rate processes for over a century with applications in biochemistry, biophysics, cell biology, systems biology, and chemical engineering. The MM rate law and its remedied form stand on the assumption that the concentration of the complex of interacting molecules, at each moment, approaches an equilibrium (quasi-steady state) much faster than the molecular concentrations change. Yet, this assumption is not always justified. Here, we relax this quasi-steady state requirement and propose the generalized MM rate law for the interactions of molecules with active concentration changes over time. Our approach for time-varying molecular concentrations, termed the effective time-delay scheme (ETS), is based on rigorously estimated time-delay effects in molecular complex formation. With particularly marked improvements in protein-protein and protein-DNA interaction modeling, the ETS provides an analytical framework to interpret and predict rich transient or rhythmic dynamics (such as autogenously-regulated cellular adaptation and circadian protein turnover), which goes beyond the quasi-steady state assumption.

An autonomous mathematical model for the mammalian cell cycle

A mathematical model for the mammalian cell cycle is developed as a system of 13 coupled nonlinear ordinary differential equations. The variables and interactions included in the model are based on detailed consideration of available experimental data. A novel feature of the model is inclusion of cycle tasks such as origin licensing and initiation, nuclear envelope breakdown and kinetochore attachment, and their interactions with controllers (molecular complexes involved in cycle control). Other key features are that the model is autonomous, except for a dependence on external growth factors; the variables are continuous in time, without instantaneous resets at phase boundaries; mechanisms to prevent rereplication are included; and cycle progression is independent of cell size. Eight variables represent cell cycle controllers: the Cyclin D1-Cdk4/6 complex, APC^Cdh1, SCF^βTrCP, Cdc25A, MPF, NuMA, the securin-separase complex, and separase. Five variables represent task completion, with four for the status of origins and one for kinetochore attachment. The model predicts distinct behaviors corresponding to the main phases of the cell cycle, showing that the principal features of the mammalian cell cycle, including restriction point behavior, can be accounted for in a quantitative mechanistic way based on known interactions among cycle controllers and their coupling to tasks. The model is robust to parameter changes, in that cycling is maintained over at least a five-fold range of each parameter when varied individually. The model is suitable for exploring how extracellular factors affect cell cycle progression, including responses to metabolic conditions and to anti-cancer therapies.

Modeling the Circadian Control of the Cell Cycle and Its Consequences for Cancer Chronotherapy

The mammalian cell cycle is governed by a network of cyclin/Cdk complexes which signal the progression into the successive phases of the cell division cycle. Once coupled to the circadian clock, this network produces oscillations with a 24 h period such that the progression into each phase of the cell cycle is synchronized to the day-night cycle. Here, we use a computational model for the circadian clock control of the cell cycle to investigate the entrainment in a population of cells characterized by some variability in the kinetic parameters. Our numerical simulations showed that successful entrainment and synchronization are only possible with a sufficient circadian amplitude and an autonomous period close to 24 h. Cellular heterogeneity, however, introduces some variability in the entrainment phase of the cells. Many cancer cells have a disrupted clock or compromised clock control. In these conditions, the cell cycle runs independently of the circadian clock, leading to a lack of synchronization of cancer cells. When the coupling is weak, entrainment is largely impacted, but cells maintain a tendency to divide at specific times of day. These differential entrainment features between healthy and cancer cells can be exploited to optimize the timing of anti-cancer drug administration in order to minimize their toxicity and to maximize their efficacy. We then used our model to simulate such chronotherapeutic treatments and to predict the optimal timing for anti-cancer drugs targeting specific phases of the cell cycle. Although qualitative, the model highlights the need to better characterize cellular heterogeneity and synchronization in cell populations as well as their consequences for circadian entrainment in order to design successful chronopharmacological protocols.

Stability analysis of a multiscale model of cell cycle dynamics coupled with quiescent and proliferating cell populations

In this paper, we perform a mathematical analysis of our proposed nonlinear, multiscale mathematical model of physiologically structured quiescent and proliferating cell populations at the macroscale and cell-cycle proteins at the microscale. Cell cycle dynamics (microscale) are driven by growth factors derived from the total cell population of quiescent and proliferating cells. Cell-cycle protein concentrations, on the other hand, determine the rates of transition between the two subpopulations. Our model demonstrates the underlying impact of cell cycle dynamics on the evolution of cell population in a tissue. We study the model's well-posedness, derive steady-state solutions, and find sufficient conditions for the stability of steady-state solutions using semigroup and spectral theory. Finally, we performed numerical simulations to see how the parameters affect the model's nonlinear dynamics.

Machine learning alternative to systems biology should not solely depend on data

In recent years, artificial intelligence (AI)/machine learning has emerged as a plausible alternative to systems biology for the elucidation of biological phenomena and in attaining specified design objective in synthetic biology. Although considered highly disruptive with numerous notable successes so far, we seek to bring attention to both the fundamental and practical pitfalls of their usage, especially in illuminating emergent behaviors from chaotic or stochastic systems in biology. Without deliberating on their suitability and the required data qualities and pre-processing approaches beforehand, the research and development community could experience similar 'AI winters' that had plagued other fields. Instead, we anticipate the integration or combination of the two approaches, where appropriate, moving forward.

Multi-synchronization and other patterns of multi-rhythmicity in oscillatory biological systems

While experimental and theoretical studies have established the prevalence of rhythmic behaviour at all levels of biological organization, less common is the coexistence between multiple oscillatory regimes (multi-rhythmicity), which has been predicted by a variety of models for biological oscillators. The phenomenon of multi-rhythmicity involves, most commonly, the coexistence between two (birhythmicity) or three (trirhythmicity) distinct regimes of self-sustained oscillations. Birhythmicity has been observed experimentally in a few chemical reactions and in biological examples pertaining to cardiac cell physiology, neurobiology, human voice patterns and ecology. The present study consists of two parts. We first review the mechanisms underlying multi-rhythmicity in models for biochemical and cellular oscillations in which the phenomenon was investigated over the years. In the second part, we focus on the coupling of the cell cycle and the circadian clock and show how an additional source of multi-rhythmicity arises from the bidirectional coupling of these two cellular oscillators. Upon bidirectional coupling, the two oscillatory networks generally synchronize in a unique manner characterized by a single, common period. In some conditions, however, the two oscillators may synchronize in two or three different ways characterized by distinct waveforms and periods. We refer to this type of multi-rhythmicity as ‘multi-synchronization’.

Teaching computational systems biology with an eye on quantitative systems pharmacology at the undergraduate level: Why do it, who would take it, and what should we teach?

Computational systems biology (CSB) is a field that emerged primarily as the product of research activities. As such, it grew in several directions in a distributed and uncoordinated manner making the area appealing and fascinating. The idea of not having to follow a specific path but instead creating one fueled innovation. As the field matured, several interdisciplinary graduate programs emerged attempting to educate future generations of computational systems biologists. These educational initiatives coordinated the dissemination of information across student populations that had already decided to specialize in this field. However, we are now entering an era where CSB, having established itself as a valuable research discipline, is attempting the next major step: Entering undergraduate curricula. As interesting as this endeavor may sound, it has several difficulties, mainly because the field is not uniformly defined. In this manuscript, we argue that this diversity is a significant advantage and that several incarnations of an undergraduate-level CSB biology course could, and should, be developed tailored to programmatic needs. In this manuscript, we share our experiences creating a course as part of a Biomedical Engineering program.

A modular approach for modeling the cell cycle based on functional response curves

Modeling biochemical reactions by means of differential equations often results in systems with a large number of variables and parameters. As this might complicate the interpretation and generalization of the obtained results, it is often desirable to reduce the complexity of the model. One way to accomplish this is by replacing the detailed reaction mechanisms of certain modules in the model by a mathematical expression that qualitatively describes the dynamical behavior of these modules. Such an approach has been widely adopted for ultrasensitive responses, for which underlying reaction mechanisms are often replaced by a single Hill function. Also time delays are usually accounted for by using an explicit delay in delay differential equations. In contrast, however, S-shaped response curves, which by definition have multiple output values for certain input values and are often encountered in bistable systems, are not easily modeled in such an explicit way. Here, we extend the classical Hill function into a mathematical expression that can be used to describe both ultrasensitive and S-shaped responses. We show how three ubiquitous modules (ultrasensitive responses, S-shaped responses and time delays) can be combined in different configurations and explore the dynamics of these systems. As an example, we apply our strategy to set up a model of the cell cycle consisting of multiple bistable switches, which can incorporate events such as DNA damage and coupling to the circadian clock in a phenomenological way.

Bistability plays an important role in many biochemical processes and typically emerges from complex interaction patterns such as positive and double negative feedback loops. Here, we propose to theoretically study the effect of bistability in a larger interaction network. We explicitly incorporate a functional expression describing an S-shaped input-output curve in the model equations, without the need for considering the underlying biochemical events. This expression can be converted into a functional module for an ultrasensitive response, and a time delay is easily included as well. Exploiting the fact that several of these modules can easily be combined in larger networks, we construct a cell cycle model consisting of multiple bistable switches and show how this approach can account for a number of known properties of the cell cycle.

Suppression of synuclein gamma inhibits the movability of endometrial carcinoma cells by PI3K/AKT/ERK signaling pathway

BACKGROUND

Although overexpression of synuclein gamma (SNCG) has been reported in several cancers, few studies have been performed onSNCG in endometrial carcinomas.

OBJECTIVE

This study aimed to investigate the role of SNCG in the progression of endometrial carcinoma.

METHODS

The expression pattern and function ofSNCG gene were analyzed using the Gene Expression Omnibus (GEO) and Gene Set Enrichment Analysis (GSEA) datasets. Two vector types, containing either SNCG or negative control shRNAs, were used to evaluate cell proliferation, apoptosis, and metastasis using Cell Counting Kit 8, colony formation, flow cytometry, wound-healing, transwell, and invasion assays. The relative protein levels of N-cadherin, E-cadherin, vimentin, p-PI3K, PI3K, p-AKT, AKT, p-ERK, and ERK were determined by western bloting.

RESULTS

Our results revealed thatSNCG mRNA expression and SNCG protein levels in shRNA-treated SPEC2 cells were lower than in the negative control cells. Furthermore, cell proliferation, migration, and invasion were significantly inhibited in SNCG shRNA-treated cells, but apoptosis was increased. The results of western blot analysis indicated that SNCG silencing reduced the protein levels of N-cadherin, vimentin, p-PI3K, p-AKT, and p-ERK, but not those of total PI3K, AKT, and ERK.

CONCLUSIONS

Therefore, shRNA-mediated suppression of SNCG inhibited SPEC2 cell proliferation, migration, and invasion, and promoted SPEC2 cell apoptosis, which was presumably accomplished via regulation of the PI3K/AKT/ERK signaling pathway.

The transcriptome dynamics of single cells during the cell cycle

The cell cycle is among the most basic phenomena in biology. Despite advances in single-cell analysis, dynamics and topology of the cell cycle in high-dimensional gene expression space remain largely unknown. We developed a linear analysis of transcriptome data which reveals that cells move along a planar circular trajectory in transcriptome space during the cycle. Non-cycling gene expression adds a third dimension causing helical motion on a cylinder. We find in immortalized cell lines that cell cycle transcriptome dynamics occur largely independently from other cellular processes. We offer a simple method ("Revelio") to order unsynchronized cells in time. Precise removal of cell cycle effects from the data becomes a straightforward operation. The shape of the trajectory implies that each gene is upregulated only once during the cycle, and only two dynamic components represented by groups of genes drive transcriptome dynamics. It indicates that the cell cycle has evolved to minimize changes of transcriptional activity and the related regulatory effort. This design principle of the cell cycle may be of relevance to many other cellular differentiation processes.

Quantifying the Landscape and Transition Paths for Proliferation-Quiescence Fate Decisions

The cell cycle, essential for biological functions, experiences delicate spatiotemporal regulation. The transition between G1 and S phase, which is called the proliferation–quiescence decision, is critical to the cell cycle. However, the stability and underlying stochastic dynamical mechanisms of the proliferation–quiescence decision have not been fully understood. To quantify the process of the proliferation–quiescence decision, we constructed its underlying landscape based on the relevant gene regulatory network. We identified three attractors on the landscape corresponding to the G0, G1, and S phases, individually, which are supported by single-cell data. By calculating the transition path, which quantifies the potential barrier, we built expression profiles in temporal order for key regulators in different transitions. We propose that the two saddle points on the landscape characterize restriction point (RP) and G1/S checkpoint, respectively, which provides quantitative and physical explanations for the mechanisms of Rb governing the RP while p21 controlling the G1/S checkpoint. We found that Emi1 inhibits the transition from G0 to G1, while Emi1 in a suitable range facilitates the transition from G1 to S. These results are partially consistent with previous studies, which also suggested new roles of Emi1 in the cell cycle. By global sensitivity analysis, we identified some critical regulatory factors influencing the proliferation–quiescence decision. Our work provides a global view of the stochasticity and dynamics in the proliferation–quiescence decision of the cell cycle.

Sample-based modeling reveals bidirectional interplay between cell cycle progression and extrinsic apoptosis

Apoptotic cell death can be initiated through the extrinsic and intrinsic signaling pathways. While cell cycle progression promotes the responsiveness to intrinsic apoptosis induced by genotoxic stress or spindle poisons, this has not yet been studied conclusively for extrinsic apoptosis. Here, we combined fluorescence-based time-lapse monitoring of cell cycle progression and cell death execution by long-term time-lapse microscopy with sampling-based mathematical modeling to study cell cycle dependency of TRAIL-induced extrinsic apoptosis in NCI-H460/geminin cells. In particular, we investigated the interaction of cell death timing and progression of cell cycle states. We not only found that TRAIL prolongs cycle progression, but in reverse also that cell cycle progression affects the kinetics of TRAIL-induced apoptosis: Cells exposed to TRAIL in G1 died significantly faster than cells stimulated in S/G2/M. The connection between cell cycle state and apoptosis progression was captured by developing a mathematical model, for which parameter estimation revealed that apoptosis progression decelerates in the second half of the cell cycle. Similar results were also obtained when studying HCT-116 cells. Our results therefore reject the null hypothesis of independence between cell cycle progression and extrinsic apoptosis and, supported by simulations and experiments of synchronized cell populations, suggest that unwanted escape from TRAIL-induced apoptosis can be reduced by enriching the fraction of cells in G1 phase. Besides novel insight into the interrelation of cell cycle progression and extrinsic apoptosis signaling kinetics, our findings are therefore also relevant for optimizing future TRAIL-based treatment strategies.

Critical slowing down and attractive manifold: A mechanism for dynamic robustness in the yeast cell-cycle process

Biological processes that execute complex multiple functions, such as the cell cycle, must ensure the order of sequential events and maintain dynamic robustness against various fluctuations. Here, we examine the mechanisms and fundamental structure that achieve these properties in the cell cycle of the budding yeast Saccharomyces cerevisiae. We show that this process behaves like an excitable system containing three well-decoupled saddle-node bifurcations to execute DNA replication and mitosis events. The yeast cell-cycle regulatory network can be divided into three modules-the G1/S phase, early M phase, and late M phase-wherein both positive feedback loops in each module and interactions among modules play important roles. Specifically, when the cell-cycle process operates near the critical points of the saddle-node bifurcations, a critical slowing down effect takes place. Such interregnum then allows for an attractive manifold and sufficient duration for cell-cycle events, within which to assess the completion of DNA replication and mitosis, e.g., spindle assembly. Moreover, such arrangement ensures that any fluctuation in an early module or event will not transmit to a later module or event. Thus, our results suggest a possible dynamical mechanism of the cell-cycle process to ensure event order and dynamic robustness and give insight into the evolution of eukaryotic cell-cycle processes.

Robust synchronization of the cell cycle and the circadian clock through bidirectional coupling

The cell cycle and the circadian clock represent major cellular rhythms, which appear to be coupled. Thus the circadian factor BMAL1 controls the level of cell cycle proteins such as Cyclin E and WEE1, the latter of which inhibits the kinase CDK1 that governs the G2/M transition. In reverse the cell cycle impinges on the circadian clock through direct control by CDK1 of REV-ERBα, which negatively regulates BMAL1. These observations provide evidence for bidirectional coupling of the cell cycle and the circadian clock. By merging detailed models for the two networks in mammalian cells, we previously showed that unidirectional coupling to the circadian clock can entrain the cell cycle to 24 or 48 h, depending on the cell cycle autonomous period, while complex oscillations occur when entrainment fails. Here we show that the reverse unidirectional coupling via phosphorylation of REV-ERBα or via mitotic inhibition of transcription, both controlled by CDK1, can elicit entrainment of the circadian clock by the cell cycle. We then determine the effect of bidirectional coupling of the cell cycle and circadian clock as a function of their relative coupling strengths. In contrast to unidirectional coupling, bidirectional coupling markedly reduces the likelihood of complex oscillations. While the two rhythms oscillate independently as long as both couplings are weak, one rhythm entrains the other if one of the couplings dominates. If the couplings in both directions become stronger and of comparable magnitude, the two rhythms synchronize, generally at an intermediate period within the range defined by the two autonomous periods prior to coupling. More surprisingly, synchronization may also occur at a period slightly below or above this range, while in some conditions the synchronization period can even be much longer. Two or even three modes of synchronization may sometimes coexist, yielding examples of birhythmicity or trirhythmicity. Because synchronization readily occurs in the form of simple periodic oscillations over a wide range of coupling strengths and in the presence of multiple connections between the two oscillatory networks, the results indicate that bidirectional coupling favours the robust synchronization of the cell cycle and the circadian clock.

Modeling the Dynamics of Let-7-Coupled Gene Regulatory Networks Linking Cell Proliferation to Malignant Transformation

Let-7 microRNA controls the expression of proteins that belong to two distinct gene regulatory networks, namely, a cyclin-dependent kinase (Cdk) network driving the cell cycle and a cell transformation network that can undergo an epigenetic switch between a non-transformed and a malignant transformed cell state. Using mathematical modeling and transcriptomic data analysis, we here investigate how Let-7 controls the Cdk-dependent cell cycle network, and how it couples the latter with the transformation network. We also assess the consequence of this coupling on cancer progression. Our analysis shows that the switch from a quiescent to a proliferative state depends on the relative levels of Let-7 and several cell cycle activators. Numerical simulations further indicate that the Let-7-coupled cell cycle and transformation networks mutually control each other, and our model identifies key players for this mutual control. Transcriptomic data analysis from The Cancer Genome Atlas (TCGA) suggests that the two networks are activated in cancer, in particular in gastrointestinal cancers, and that the activation levels vary significantly among patients affected by a same cancer type. Our mathematical model, when applied to a heterogeneous cell population, suggests that heterogeneity among tumors may in part result from stochastic switches between a non-transformed cell state with low proliferative capability and a transformed cell state with high proliferative property. The model further predicts that Let-7 may reduce tumor heterogeneity by decreasing the occurrence of stochastic switches toward a transformed, proliferative cell state. In conclusion, we identified the key components responsible for the qualitative dynamics of two networks interconnected by Let-7. The two networks are heterogeneously activated in several cancers, thereby stressing the need to consider patient’s specific characteristics to optimize therapeutic strategies.

Multi-rhythmicity generated by coupling two cellular rhythms

The cell cycle and the circadian clock represent two major cellular rhythms, which are coupled because the circadian clock governs the synthesis of several proteins of the network that drives the mammalian cell cycle. Analysis of a detailed model for these coupled cellular rhythms previously showed that the cell cycle can be entrained at the circadian period of 24 h, or at a period of 48 h, depending on the autonomous period of the cell cycle and on the coupling strength. We show by means of numerical simulations that multiple stable periodic regimes, i.e. multi-rhythmicity, may originate from the coupling of the two cellular rhythms. In these conditions, the cell cycle can evolve to any one of two (birhythmicity) or three stable periodic regimes (trirhythmicity). When applied at the right phase, transient perturbations of appropriate duration and magnitude can induce the switch between the different oscillatory states. Such switching is characterized by final state sensitivity, which originates from the complex structure of the attraction basins. By providing a novel instance of multi-rhythmicity in a realistic model for the coupling of two major cellular rhythms, the results throw light on the conditions in which multiple stable periodic regimes may coexist in biological systems.

Revisiting a skeleton model for the mammalian cell cycle: From bistability to Cdk oscillations and cellular heterogeneity

A network of cyclin-dependent kinases (Cdks) regulated by multiple negative and positive feedback loops controls progression in the mammalian cell cycle. We previously proposed a detailed computational model for this network, which consists of four coupled Cdk modules. Both this detailed model and a reduced, skeleton version show that the Cdk network is capable of temporal self-organization in the form of sustained Cdk oscillations, which correspond to the orderly progression along the different cell cycle phases G1, S (DNA replication), G2 and M (mitosis). We use the skeleton model to revisit the role of positive feedback (PF) loops on the dynamics of the mammalian cell cycle by showing that the multiplicity of PF loops extends the range of bistability in the isolated Cdk modules controlling the G1/S and G2/M transitions. Resorting to stochastic simulations we show that, through their effect on the range of bistability, multiple PF loops enhance the robustness of Cdk oscillations with respect to molecular noise. The model predicts that a rise in the total level of Cdk1 also enlarges the domain of bistability in the isolated Cdk modules as well as the range of oscillations in the full Cdk network. Surprisingly, stochastic simulations indicate that Cdk1 overexpression reduces the robustness of Cdk oscillations towards molecular noise; this result is due to the increased distance between the two branches of the bistable switch at higher levels of Cdk1. At intermediate levels of growth factor stochastic simulations show that cells may randomly switch between cell cycle arrest and cell proliferation, as a consequence of fluctuations. In the presence of Cdk1 overexpression, these transitions occur even at low levels of growth factor. Extending stochastic simulations from single cells to cell populations suggests that stochastic switches between cell cycle arrest and proliferation may provide a source of heterogeneity in a cell population, as observed in cancer cells characterized by Cdk1 overexpression.

Modeling the Influence of Seasonal Differences in the HPA Axis on Synchronization of the Circadian Clock and Cell Cycle

Synchronization of biological functions to environmental signals enables organisms to anticipate and appropriately respond to daily external fluctuations and is critical to the maintenance of homeostasis. Misalignment of circadian rhythms with environmental cues is associated with adverse health outcomes. Cortisol, the downstream effector of hypothalamic-pituitary-adrenal (HPA) activity, facilitates synchronization of peripheral biological processes to the environment. Cortisol levels exhibit substantial seasonal rhythmicity, with peak levels occurring during the short-photoperiod winter months and reduced levels occurring in the long-photoperiod summer season. Seasonal changes in cortisol secretion could therefore alter its entraining capabilities, resulting in a season-dependent modification in the alignment of biological activities with the environment. We develop a mathematical model to investigate the influence of photoperiod-induced seasonal differences in the circadian rhythmicity of the HPA axis on the synchronization of the peripheral circadian clock and cell cycle in a heterogeneous cell population. Model simulations predict that the high-amplitude cortisol rhythms in winter result in the greatest entrainment of peripheral oscillators. Furthermore, simulations predict a circadian gating of the cell cycle with respect to the expression of peripheral clock genes. Seasonal differences in cortisol rhythmicity are also predicted to influence mitotic synchrony, with a high-amplitude winter rhythm resulting in the greatest synchrony and a shift in timing of the cell cycle phases, relative to summer. Our results highlight the primary interactions among the HPA axis, the peripheral circadian clock, and the cell cycle and thereby provide an improved understanding of the implications of circadian misalignment on the synchronization of peripheral regulatory processes.

Minimum Action Path Theory Reveals the Details of Stochastic Transitions Out of Oscillatory States

Cell state determination is the outcome of intrinsically stochastic biochemical reactions. Transitions between such states are studied as noise-driven escape problems in the chemical species space. Escape can occur via multiple possible multidimensional paths, with probabilities depending nonlocally on the noise. Here we characterize the escape from an oscillatory biochemical state by minimizing the Freidlin-Wentzell action, deriving from it the stochastic spiral exit path from the limit cycle. We also use the minimized action to infer the escape time probability density function.

Coarse-graining and hybrid methods for efficient simulation of stochastic multi-scale models of tumour growth

The development of hybrid methodologies is of current interest in both multi-scale modelling and stochastic reaction-diffusion systems regarding their applications to biology. We formulate a hybrid method for stochastic multi-scale models of cells populations that extends the remit of existing hybrid methods for reaction-diffusion systems. Such method is developed for a stochastic multi-scale model of tumour growth, i.e. population-dynamical models which account for the effects of intrinsic noise affecting both the number of cells and the intracellular dynamics. In order to formulate this method, we develop a coarse-grained approximation for both the full stochastic model and its mean-field limit. Such approximation involves averaging out the age-structure (which accounts for the multi-scale nature of the model) by assuming that the age distribution of the population settles onto equilibrium very fast. We then couple the coarse-grained mean-field model to the full stochastic multi-scale model. By doing so, within the mean-field region, we are neglecting noise in both cell numbers (population) and their birth rates (structure). This implies that, in addition to the issues that arise in stochastic-reaction diffusion systems, we need to account for the age-structure of the population when attempting to couple both descriptions. We exploit our coarse-graining model so that, within the mean-field region, the age-distribution is in equilibrium and we know its explicit form. This allows us to couple both domains consistently, as upon transference of cells from the mean-field to the stochastic region, we sample the equilibrium age distribution. Furthermore, our method allows us to investigate the effects of intracellular noise, i.e. fluctuations of the birth rate, on collective properties such as travelling wave velocity. We show that the combination of population and birth-rate noise gives rise to large fluctuations of the birth rate in the region at the leading edge of front, which cannot be accounted for by the coarse-grained model. Such fluctuations have non-trivial effects on the wave velocity. Beyond the development of a new hybrid method, we thus conclude that birth-rate fluctuations are central to a quantitatively accurate description of invasive phenomena such as tumour growth.

Dissipative structures and biological rhythms

Sustained oscillations abound in biological systems. They occur at all levels of biological organization over a wide range of periods, from a fraction of a second to years, and with a variety of underlying mechanisms. They control major physiological functions, and their dysfunction is associated with a variety of physiological disorders. The goal of this review is (i) to give an overview of the main rhythms observed at the cellular and supracellular levels, (ii) to briefly describe how the study of biological rhythms unfolded in the course of time, in parallel with studies on chemical oscillations, (iii) to present the major roles of biological rhythms in the control of physiological functions, and (iv) the pathologies associated with the alteration, disappearance, or spurious occurrence of biological rhythms. Two tables present the main examples of cellular and supracellular rhythms ordered according to their period, and their role in physiology and pathophysiology. Among the rhythms discussed are neural and cardiac rhythms, metabolic oscillations such as those occurring in glycolysis in yeast, intracellular Ca⁺⁺ oscillations, cyclic AMP oscillations in Dictyostelium amoebae, the segmentation clock that controls somitogenesis, pulsatile hormone secretion, circadian rhythms which occur in all eukaryotes and some bacteria with a period close to 24 h, the oscillatory dynamics of the enzymatic network driving the cell cycle, and oscillations in transcription factors such as NF-ΚB and tumor suppressors such as p53. Ilya Prigogine's concept of dissipative structures applies to temporal oscillations and allows us to unify within a common framework the various rhythms observed at different levels of biological organization, regardless of their period and underlying mechanism.

Steady-State-Preserving Simulation of Genetic Regulatory Systems

A novel family of exponential Runge-Kutta (expRK) methods are designed incorporating the stable steady-state structure of genetic regulatory systems. A natural and convenient approach to constructing new expRK methods on the base of traditional RK methods is provided. In the numerical integration of the one-gene, two-gene, and p53-mdm2 regulatory systems, the new expRK methods are shown to be more accurate than their prototype RK methods. Moreover, for nonstiff genetic regulatory systems, the expRK methods are more efficient than some traditional exponential RK integrators in the scientific literature.

Stochastic multi-scale models of competition within heterogeneous cellular populations: Simulation methods and mean-field analysis

We propose a modelling framework to analyse the stochastic behaviour of heterogeneous, multi-scale cellular populations. We illustrate our methodology with a particular example in which we study a population with an oxygen-regulated proliferation rate. Our formulation is based on an age-dependent stochastic process. Cells within the population are characterised by their age (i.e. time elapsed since they were born). The age-dependent (oxygen-regulated) birth rate is given by a stochastic model of oxygen-dependent cell cycle progression. Once the birth rate is determined, we formulate an age-dependent birth-and-death process, which dictates the time evolution of the cell population. The population is under a feedback loop which controls its steady state size (carrying capacity): cells consume oxygen which in turn fuels cell proliferation. We show that our stochastic model of cell cycle progression allows for heterogeneity within the cell population induced by stochastic effects. Such heterogeneous behaviour is reflected in variations in the proliferation rate. Within this set-up, we have established three main results. First, we have shown that the age to the G1/S transition, which essentially determines the birth rate, exhibits a remarkably simple scaling behaviour. Besides the fact that this simple behaviour emerges from a rather complex model, this allows for a huge simplification of our numerical methodology. A further result is the observation that heterogeneous populations undergo an internal process of quasi-neutral competition. Finally, we investigated the effects of cell-cycle-phase dependent therapies (such as radiation therapy) on heterogeneous populations. In particular, we have studied the case in which the population contains a quiescent sub-population. Our mean-field analysis and numerical simulations confirm that, if the survival fraction of the therapy is too high, rescue of the quiescent population occurs. This gives rise to emergence of resistance to therapy since the rescued population is less sensitive to therapy.

Exponentially Fitted Two-Derivative Runge-Kutta Methods for Simulation of Oscillatory Genetic Regulatory Systems

Oscillation is one of the most important phenomena in the chemical reaction systems in living cells. The general purpose simulation algorithms fail to take into account this special character and produce unsatisfying results. In order to enhance the accuracy of the integrator, the second-order derivative is incorporated in the scheme. The oscillatory feature of the solution is captured by the integrators with an exponential fitting property. Three practical exponentially fitted TDRK (EFTDRK) methods are derived. To test the effectiveness of the new EFTDRK methods, the two-gene system with cross-regulation and the circadian oscillation of the period protein in Drosophila are simulated. Each EFTDRK method has the best fitting frequency which minimizes the global error. The numerical results show that the new EFTDRK methods are more accurate and more efficient than their prototype TDRK methods or RK methods of the same order and the traditional exponentially fitted RK method in the literature.

Dynamics of the mammalian cell cycle in physiological and pathological conditions

A network of cyclin-dependent kinases (Cdks) controls progression along the successive phases G1, S, G2, and M of the mammalian cell cycle. Deregulations in the expression of molecular components in this network often lead to abusive cell proliferation and cancer. Given the complex nature of the Cdk network, it is fruitful to resort to computational models to grasp its dynamical properties. Investigated by means of bifurcation diagrams, a detailed computational model for the Cdk network shows how the balance between quiescence and proliferation is affected by activators (oncogenes) and inhibitors (tumor suppressors) of cell cycle progression, as well as by growth factors and other external factors such as the extracellular matrix (ECM) and cell contact inhibition. Suprathreshold changes in all these factors can trigger a switch in the dynamical behavior of the network corresponding to a bifurcation between a stable steady state, associated with cell cycle arrest, and sustained oscillations of the various cyclin/Cdk complexes, corresponding to cell proliferation. The model for the Cdk network accounts for the dependence or independence of cell proliferation on serum and/or cell anchorage to the ECM. Such computational approach provides an integrated view of the control of cell proliferation in physiological or pathological conditions. Whether the balance is tilted toward cell cycle arrest or cell proliferation depends on the direction in which the threshold associated with the bifurcation is passed once the cell integrates the multiple signals, internal or external to the Cdk network, that promote or impede progression in the cell cycle.

A hybrid model of cell cycle in mammals

Time plays an essential role in many biological systems, especially in cell cycle. Many models of biological systems rely on differential equations, but parameter identification is an obstacle to use differential frameworks. In this paper, we present a new hybrid modeling framework that extends René Thomas' discrete modeling. The core idea is to associate with each qualitative state "celerities" allowing us to compute the time spent in each state. This hybrid framework is illustrated by building a 5-variable model of the mammalian cell cycle. Its parameters are determined by applying formal methods on the underlying discrete model and by constraining parameters using timing observations on the cell cycle. This first hybrid model presents the most important known behaviors of the cell cycle, including quiescent phase and endoreplication.

A 1.26 μW Cytomimetic IC Emulating Complex Nonlinear Mammalian Cell Cycle Dynamics: Synthesis, Simulation and Proof-of-Concept Measured Results

Cytomimetic circuits represent a novel, ultra low-power, continuous-time, continuous-value class of circuits, capable of mapping on silicon cellular and molecular dynamics modelled by means of nonlinear ordinary differential equations (ODEs). Such monolithic circuits are in principle able to emulate on chip, single or multiple cell operations in a highly parallel fashion. Cytomimetic topologies can be synthesized by adopting the Nonlinear Bernoulli Cell Formalism (NBCF), a mathematical framework that exploits the striking similarities between the equations describing weakly-inverted Metal-Oxide Semiconductor (MOS) devices and coupled nonlinear ODEs, typically appearing in models of naturally encountered biochemical systems. The NBCF maps biological state variables onto strictly positive subthreshold MOS circuit currents. This paper presents the synthesis, the simulation and proof-of-concept chip results corresponding to the emulation of a complex cellular network mechanism, the skeleton model for the network of Cyclin-dependent Kinases (CdKs) driving the mammalian cell cycle. This five variable nonlinear biological model, when appropriate model parameter values are assigned, can exhibit multiple oscillatory behaviors, varying from simple periodic oscillations, to complex oscillations such as quasi-periodicity and chaos. The validity of our approach is verified by simulated results with realistic process parameters from the commercially available AMS 0.35 μm technology and by chip measurements. The fabricated chip occupies an area of 2.27 mm2 and consumes a power of 1.26 μW from a power supply of 3 V. The presented cytomimetic topology follows closely the behavior of its biological counterpart, exhibiting similar time-dependent solutions of the Cdk complexes, the transcription factors and the proteins.

Geometric analysis of the Goldbeter minimal model for the embryonic cell cycle

A minimal model describing the embryonic cell division cycle at the molecular level in eukaryotes is analyzed mathematically. It is known from numerical simulations that the corresponding three-dimensional system of ODEs has periodic solutions in certain parameter regimes. We prove the existence of a stable limit cycle and provide a detailed description on how the limit cycle is generated. The limit cycle corresponds to a relaxation oscillation of an auxiliary system, which is singularly perturbed and has the same orbits as the original model. The singular perturbation character of the auxiliary problem is caused by the occurrence of small Michaelis constants in the model. Essential pieces of the limit cycle of the auxiliary problem consist of segments of slow motion close to several branches of a two dimensional critical manifold which are connected by fast jumps. In addition, a new phenomenon of exchange of stability occurs at lines, where the branches of the two-dimensional critical manifold intersect. This novel type of relaxation oscillations is studied by combining standard results from geometric singular perturbation with several suitable blow-up transformations.

Modelling cell cycle synchronisation in networks of coupled radial glial cells

Radial glial cells play a crucial role in the embryonic mammalian brain. Their proliferation is thought to be controlled, in part, by ATP mediated calcium signals. It has been hypothesised that these signals act to locally synchronise cell cycles, so that clusters of cells proliferate together, shedding daughter cells in uniform sheets. In this paper we investigate this cell cycle synchronisation by taking an ordinary differential equation model that couples the dynamics of intracellular calcium and the cell cycle and extend it to populations of cells coupled via extracellular ATP signals. Through bifurcation analysis we show that although ATP mediated calcium release can lead to cell cycle synchronisation, a number of other asynchronous oscillatory solutions including torus solutions dominate the parameter space and cell cycle synchronisation is far from guaranteed. Despite this, numerical results indicate that the transient and not the asymptotic behaviour of the system is important in accounting for cell cycle synchronisation. In particular, quiescent cells can be entrained on to the cell cycle via ATP mediated calcium signals initiated by a driving cell and crucially will cycle in near synchrony with the driving cell for the duration of neurogenesis. This behaviour is highly sensitive to the timing of ATP release, with release at the G1/S phase transition of the cell cycle far more likely to lead to near synchrony than release during mid G1 phase. This result, which suggests that ATP release timing is critical to radial glia cell cycle synchronisation, may help us to understand normal and pathological brain development.

The balance between cell cycle arrest and cell proliferation: control by the extracellular matrix and by contact inhibition

To understand the dynamics of the cell cycle, we need to characterize the balance between cell cycle arrest and cell proliferation, which is often deregulated in cancers. We address this issue by means of a detailed computational model for the network of cyclin-dependent kinases (Cdks) driving the mammalian cell cycle. Previous analysis of the model focused on how this balance is controlled by growth factors (GFs) or the levels of activators (oncogenes) and inhibitors (tumour suppressors) of cell cycle progression. Supra-threshold changes in the level of any of these factors can trigger a switch in the dynamical behaviour of the Cdk network corresponding to a bifurcation between a stable steady state, associated with cell cycle arrest, and sustained oscillations of the various cyclin/Cdk complexes, corresponding to cell proliferation. Here, we focus on the regulation of cell proliferation by cellular environmental factors external to the Cdk network, such as the extracellular matrix (ECM), and contact inhibition, which increases with cell density. We extend the model for the Cdk network by including the phenomenological effect of both the ECM, which controls the activation of the focal adhesion kinase (FAK) that promotes cell cycle progression, and cell density, which inhibits cell proliferation via the Hippo/YAP pathway. The model shows that GFs and FAK activation are capable of triggering in a similar dynamical manner the transition to cell proliferation, while the Hippo/YAP pathway can arrest proliferation once cell density passes a critical threshold. The results account for the dependence or independence of cell proliferation on serum and/or cell anchorage to ECM. Whether the balance in the Cdk network is tilted towards cell cycle arrest or proliferation depends on the direction in which the threshold associated with the bifurcation is passed once the cell integrates the multiple, internal or external signals that promote or impede progression in the cell cycle.

Dynamical and topological robustness of the mammalian cell cycle network: a reverse engineering approach

A common gene regulatory network model is the threshold Boolean network, used for example to model the Arabidopsis thaliana floral morphogenesis network or the fission yeast cell cycle network. In this paper, we analyze a logical model of the mammalian cell cycle network and its threshold Boolean network equivalent. Firstly, the robustness of the network was explored with respect to update perturbations, in particular, what happened to the attractors for all the deterministic updating schemes. Results on the number of different limit cycles, limit cycle lengths, basin of attraction size, for all the deterministic updating schemes were obtained through mathematical and computational tools. Secondly, we analyzed the topology robustness of the network, by reconstructing synthetic networks that contained exactly the same attractors as the original model by means of a swarm intelligence approach. Our results indicate that networks may not be very robust given the great variety of limit cycles that a network can obtain depending on the updating scheme. In addition, we identified an omnipresent network with interactions that match with the original model as well as the discovery of new interactions. The techniques presented in this paper are general, and can be used to analyze other logical or threshold Boolean network models of gene regulatory networks.

Oscillatory enzyme reactions and Michaelis-Menten kinetics

Oscillations occur in a number of enzymatic systems as a result of feedback regulation. How Michaelis-Menten kinetics influences oscillatory behavior in enzyme systems is investigated in models for oscillations in the activity of phosphofructokinase (PFK) in glycolysis and of cyclin-dependent kinases in the cell cycle. The model for the PFK reaction is based on a product-activated allosteric enzyme reaction coupled to enzymatic degradation of the reaction product. The Michaelian nature of the product decay term markedly influences the period, amplitude and waveform of the oscillations. Likewise, a model for oscillations of Cdc2 kinase in embryonic cell cycles based on Michaelis-Menten phosphorylation-dephosphorylation kinetics shows that the occurrence and amplitude of the oscillations strongly depend on the ultrasensitivity of the enzymatic cascade that controls the activity of the cyclin-dependent kinase.

MicroRNA-mediated regulation in biological systems with oscillatory behavior

As a class of small noncoding RNAs, microRNAs (miRNAs) regulate stability or translation of mRNA transcripts. Some reports bring new insights into possible roles of microRNAs in modulating cell cycle. In this paper, we focus on the mechanism and effectiveness of microRNA-mediated regulation in the cell cycle. We first describe two specific regulatory circuits that incorporate base-pairing microRNAs and show their fine-tuning roles in the modulation of periodic behavior. Furthermore, we analyze the effects of miR369-3 on the modulation of the cell cycle, confirming that miR369-3 plays a role in shortening the period of the cell cycle. These results are consistent with experimental observations.

Targeting cell cycle regulation in cancer therapy

Cell proliferation is an essential mechanism for growth, development and regeneration of eukaryotic organisms; however, it is also the cause of one of the most devastating diseases of our era: cancer. Given the relevance of the processes in which cell proliferation is involved, its regulation is of paramount importance for multicellular organisms. Cell division is orchestrated by a complex network of interactions between proteins, metabolism and microenvironment including several signaling pathways and mechanisms of control aiming to enable cell proliferation only in response to specific stimuli and under adequate conditions. Three main players have been identified in the coordinated variation of the many molecules that play a role in cell cycle: i) The cell cycle protein machinery including cyclin-dependent kinases (CDK)-cyclin complexes and related kinases, ii) The metabolic enzymes and related metabolites and iii) The reactive-oxygen species (ROS) and cellular redox status. The role of these key players and the interaction between oscillatory and non-oscillatory species have proved essential for driving the cell cycle. Moreover, cancer development has been associated to defects in all of them. Here, we provide an overview on the role of CDK-cyclin complexes, metabolic adaptations and oxidative stress in regulating progression through each cell cycle phase and transitions between them. Thus, new approaches for the design of innovative cancer therapies targeting crosstalk between cell cycle simultaneous events are proposed.

From quiescence to proliferation: Cdk oscillations drive the mammalian cell cycle

We recently proposed a detailed model describing the dynamics of the network of cyclin-dependent kinases (Cdks) driving the mammalian cell cycle (Gérard and Goldbeter, 2009). The model contains four modules, each centered around one cyclin/Cdk complex. Cyclin D/Cdk4–6 and cyclin E/Cdk2 promote progression in G1 and elicit the G1/S transition, respectively; cyclin A/Cdk2 ensures progression in S and the transition S/G2, while the activity of cyclin B/Cdk1 brings about the G2/M transition. This model shows that in the presence of sufficient amounts of growth factor the Cdk network is capable of temporal self-organization in the form of sustained oscillations, which correspond to the ordered, sequential activation of the various cyclin/Cdk complexes that control the successive phases of the cell cycle. The results suggest that the switch from cellular quiescence to cell proliferation corresponds to the transition from a stable steady state to sustained oscillations in the Cdk network. The transition depends on a finely tuned balance between factors that promote or hinder progression in the cell cycle. We show that the transition from quiescence to proliferation can occur in multiple ways that alter this balance. By resorting to bifurcation diagrams, we analyze the mechanism of oscillations in the Cdk network. Finally, we show that the complexity of the detailed model can be greatly reduced, without losing its key dynamical properties, by considering a skeleton model for the Cdk network. Using such a skeleton model for the mammalian cell cycle we show that positive feedback (PF) loops enhance the amplitude and the robustness of Cdk oscillations with respect to molecular noise. We compare the relative merits of the detailed and skeleton versions of the model for the Cdk network driving the mammalian cell cycle.

"The Octet": Eight Protein Kinases that Control Mammalian DNA Replication

Development of a fertilized human egg into an average sized adult requires about 29 trillion cell divisions, thereby producing enough DNA to stretch to the Sun and back 200 times (DePamphilis and Bell, 2011)! Even more amazing is the fact that throughout these mitotic cell cycles, the human genome is duplicated once and only once each time a cell divides. If a cell accidentally begins to re-replicate its nuclear DNA prior to cell division, checkpoint pathways trigger apoptosis. And yet, some cells are developmentally programmed to respond to environmental cues by switching from mitotic cell cycles to endocycles, a process in which multiple S phases occur in the absence of either mitosis or cytokinesis. Endocycles allow production of viable, differentiated, polyploid cells that no longer proliferate. What is surprising is that among the 516 (Manning et al., 2002) to 557 (BioMart web site) protein kinases encoded by the human genome, only eight regulate nuclear DNA replication directly. These are Cdk1, Cdk2, Cdk4, Cdk6, Cdk7, Cdc7, Checkpoint kinase-1 (Chk1), and Checkpoint kinase-2. Even more remarkable is the fact that only four of these enzymes (Cdk1, Cdk7, Cdc7, and Chk1) are essential for mammalian development. Here we describe how these protein kinases determine when DNA replication occurs during mitotic cell cycles, how mammalian cells switch from mitotic cell cycles to endocycles, and how cancer cells can be selectively targeted for destruction by inducing them to begin a second S phase before mitosis is complete.

Systems biology of cellular rhythms

Rhythms abound in biological systems, particularly at the cellular level where they originate from the feedback loops present in regulatory networks. Cellular rhythms can be investigated both by experimental and modeling approaches, and thus represent a prototypic field of research for systems biology. They have also become a major topic in synthetic biology. We review advances in the study of cellular rhythms of biochemical rather than electrical origin by considering a variety of oscillatory processes such as Ca++ oscillations, circadian rhythms, the segmentation clock, oscillations in p53 and NF-κB, synthetic oscillators, and the oscillatory dynamics of cyclin-dependent kinases driving the cell cycle. Finally we discuss the coupling between cellular rhythms and their robustness with respect to molecular noise.

Entrainment of the mammalian cell cycle by the circadian clock: modeling two coupled cellular rhythms

The cell division cycle and the circadian clock represent two major cellular rhythms. These two periodic processes are coupled in multiple ways, given that several molecular components of the cell cycle network are controlled in a circadian manner. For example, in the network of cyclin-dependent kinases (Cdks) that governs progression along the successive phases of the cell cycle, the synthesis of the kinase Wee1, which inhibits the G2/M transition, is enhanced by the complex CLOCK-BMAL1 that plays a central role in the circadian clock network. Another component of the latter network, REV-ERBα, inhibits the synthesis of the Cdk inhibitor p21. Moreover, the synthesis of the oncogene c-Myc, which promotes G1 cyclin synthesis, is repressed by CLOCK-BMAL1. Using detailed computational models for the two networks we investigate the conditions in which the mammalian cell cycle can be entrained by the circadian clock. We show that the cell cycle can be brought to oscillate at a period of 24 h or 48 h when its autonomous period prior to coupling is in an appropriate range. The model indicates that the combination of multiple modes of coupling does not necessarily facilitate entrainment of the cell cycle by the circadian clock. Entrainment can also occur as a result of circadian variations in the level of a growth factor controlling entry into G1. Outside the range of entrainment, the coupling to the circadian clock may lead to disconnected oscillations in the cell cycle and the circadian system, or to complex oscillatory dynamics of the cell cycle in the form of endoreplication, complex periodic oscillations or chaos. The model predicts that the transition from entrainment to 24 h or 48 h might occur when the strength of coupling to the circadian clock or the level of growth factor decrease below critical values.

The cell cycle and the circadian clock are two major cellular rhythms. These two periodic processes are tightly coupled through multiple regulatory interactions; several components of the cell cycle machinery are indeed controlled by the circadian network. By using detailed computational models for the cell cycle and circadian networks we investigate the conditions in which the mammalian cell cycle can be entrained by the circadian clock. We show that entrainment to a circadian period can occur when the period of the cell cycle prior to coupling is either smaller or larger than 24 h. Entrainment to 48 h can also be observed. The presence of multiple modes of coupling does not enlarge the domain of entrainment. Coupling to the circadian clock may also lead to complex oscillatory dynamics of the cell cycle in the form of endoreplication, complex periodic oscillations, or chaotic oscillations. The model predicts that entrainment of the cell cycle could also result from the circadian variation of a growth factor gating entry into G1, and that the transition from an entrained period of 24 h to 48 h might result from a decrease in coupling strength or in the level of growth factor.

Optimizing cancer pharmacotherapeutics using mathematical modeling and a systems biology approach

Research in mathematics and in mathematical biology on cancer and its treatments has been soaring in the past 10 years at an unprecedented speed. Such thriving is likely due as much to new findings in fundamental biology as to an emerging general interest from mathematicians and engineers towards applications in biology and medicine and to their subsequently designed representations and predictions of tumor processes that are now allowed by modern means of computation and simulation. This article, which does not claim the status of an extended review paper on mathematical models of cancer and its treatment, is focused on modeling in a systems biology perspective. I will list here the most necessary mathematical methods, in my opinion, which, while enforcing already existing methods, should be further developed towards designing and applying optimized individualized treatments of cancer in the clinic.

An automaton model for the cell cycle

WE CONSIDER AN AUTOMATON MODEL THAT PROGRESSES SPONTANEOUSLY THROUGH THE FOUR SUCCESSIVE PHASES OF THE CELL CYCLE: G1, S (DNA replication), G2 and M (mitosis). Each phase is characterized by a mean duration τ and a variability V. As soon as the prescribed duration of a given phase has passed, the transition to the next phase of the cell cycle occurs. The time at which the transition takes place varies in a random manner according to a distribution of durations of the cell cycle phases. Upon completion of the M phase, the cell divides into two cells, which immediately enter a new cycle in G1. The duration of each phase is reinitialized for the two newborn cells. At each time step in any phase of the cycle, the cell has a certain probability to be marked for exiting the cycle and dying at the nearest G1/S or G2/M transition. To allow for homeostasis, which corresponds to maintenance of the total cell number, we assume that cell death counterbalances cell replication at mitosis. In studying the dynamics of this automaton model, we examine the effect of factors such as the mean durations of the cell cycle phases and their variability, the type of distribution of the durations, the number of cells, the regulation of the cell population size and the independence of steady-state proportions of cells in each phase with respect to initial conditions. We apply the stochastic automaton model for the cell cycle to the progressive desynchronization of cell populations and to their entrainment by the circadian clock. A simple deterministic model leads to the same steady-state proportions of cells in the four phases of the cell cycle.