Modelling Cell Growth and its Modulation of the G1/S Transition

Systems Level Modeling of the Cell Cycle Using Budding Yeast

Proteins involved in the regulation of the cell cycle are highly conserved across all eukaryotes, and so a relatively simple eukaryote such as yeast can provide insight into a variety of cell cycle perturbations including those that occur in human cancer. To date, the budding yeast Saccharomyces cerevisiae has provided the largest amount of experimental and modeling data on the progression of the cell cycle, making it a logical choice for in-depth studies of this process. Moreover, the advent of methods for collection of high-throughput genome, transcriptome, and proteome data has provided a means to collect and precisely quantify simultaneous cell cycle gene transcript and protein levels, permitting modeling of the cell cycle on the systems level. With the appropriate mathematical framework and sufficient and accurate data on cell cycle components, it should be possible to create a model of the cell cycle that not only effectively describes its operation, but can also predict responses to perturbations such as variation in protein levels and responses to external stimuli including targeted inhibition by drugs. In this review, we summarize existing data on the yeast cell cycle, proteomics technologies for quantifying cell cycle proteins, and the mathematical frameworks that can integrate this data into representative and effective models. Systems level modeling of the cell cycle will require the integration of high-quality data with the appropriate mathematical framework, which can currently be attained through the combination of dynamic modeling based on proteomics data and using yeast as a model organism.

A quantitative model of the initiation of DNA replication in Saccharomyces cerevisiae predicts the effects of system perturbations

Background

Eukaryotic cell proliferation involves DNA replication, a tightly regulated process mediated by a multitude of protein factors. In budding yeast, the initiation of replication is facilitated by the heterohexameric origin recognition complex (ORC). ORC binds to specific origins of replication and then serves as a scaffold for the recruitment of other factors such as Cdt1, Cdc6, the Mcm2-7 complex, Cdc45 and the Dbf4-Cdc7 kinase complex. While many of the mechanisms controlling these associations are well documented, mathematical models are needed to explore the network’s dynamic behaviour. We have developed an ordinary differential equation-based model of the protein-protein interaction network describing replication initiation.

Results

The model was validated against quantified levels of protein factors over a range of cell cycle timepoints. Using chromatin extracts from synchronized Saccharomyces cerevisiae cell cultures, we were able to monitor the in vivo fluctuations of several of the aforementioned proteins, with additional data obtained from the literature. The model behaviour conforms to perturbation trials previously reported in the literature, and accurately predicts the results of our own knockdown experiments. Furthermore, we successfully incorporated our replication initiation model into an established model of the entire yeast cell cycle, thus providing a comprehensive description of these processes.

Conclusions

This study establishes a robust model of the processes driving DNA replication initiation. The model was validated against observed cell concentrations of the driving factors, and characterizes the interactions between factors implicated in eukaryotic DNA replication. Finally, this model can serve as a guide in efforts to generate a comprehensive model of the mammalian cell cycle in order to explore cancer-related phenotypes.

Molecular systems biology of Sic1 in yeast cell cycle regulation through multiscale modeling

Cell cycle control is highly regulated to guarantee the precise timing of events essential for cell growth, i.e., DNA replication onset and cell division. Failure of this control plays a role in cancer and molecules called cyclin-dependent kinase (Cdk) inhibitors (Ckis) exploit a critical function in cell cycle timing. Here we present a multiscale modeling where experimental and computational studies have been employed to investigate structure, function and temporal dynamics of the Cki Sic1 that regulates cell cycle progression in Saccharomyces cerevisiae. Structural analyses reveal molecular details of the interaction between Sic1 and Cdk/cyclin complexes, and biochemical investigation reveals Sic1 function in analogy to its human counterpart p27(Kip1), whose deregulation leads to failure in timing of kinase activation and, therefore, to cancer. Following these findings, a bottom-up systems biology approach has been developed to characterize modular networks addressing Sic1 regulatory function. Through complementary experimentation and modeling, we suggest a mechanism that underlies Sic1 function in controlling temporal waves of cyclins to ensure correct timing of the phase-specific Cdk activities.

A low number of SIC1 mRNA molecules ensures a low noise level in cell cycle progression of budding yeast

The budding yeast genome comprises roughly 6000 genes generating a number of about 10 000 mRNA copies, which gives a general estimation of 1-2 mRNA copies generated per gene. What does this observation implicate for cellular processes and their regulation? Whether the number of mRNA molecules produced is important for setting the amount of proteins implicated in a particular function is at present unknown. In this context, we studied cell cycle control as one of the highly fine tuned processes that guarantee the precise timing of events essential for cell growth. Here, we developed a stochastic model that addresses the effect of varying the mRNA amount of Sic1, inhibitor of the Cdk1-Clb5 kinase activity, and the resulting noise on Sic1/Clb5 balance at the G1/S transition. We considered a range of SIC1 transcripts number according to our experimental data derived from the MS2 mRNA tagging system. Computational simulation revealed that an increased amount of SIC1 mRNAs lead to an amplified dispersion of Sic1 protein levels, suggesting mRNA control being critical to set timing of Sic1 downregulation and, therefore, S phase onset. Moreover, Sic1/Clb5 balance is strongly influenced by Clb5 production in both daughter and mother cells in order to maintain the characteristic time of S phase entry overall the population. Furthermore, CLB5 mRNA molecules calculated to reproduce temporal dynamics of Sic1 and Clb5 for daughter and mother cells agree with recent data obtained from more complex networks. Thus, the results presented here provide novel insights into the influence that the mRNA amount and, indirectly, the transcription process exploit on cell cycle progression.