Hybrid modeling and simulation of stochastic effects on progression through the eukaryotic cell cycle.
Ontology highlight
ABSTRACT:
SUBMITTER: Liu Z
PROVIDER: S-EPMC3272065 | biostudies-other | 2012 Jan
REPOSITORIES: biostudies-other
Similar Datasets
Project description:The eukaryotic cell cycle is regulated by a complicated chemical reaction network. Although many deterministic models have been proposed, stochastic models are desired to capture noise in the cell resulting from low numbers of critical species. However, converting a deterministic model into one that accurately captures stochastic effects can result in a complex model that is hard to build and expensive to simulate. In this paper, we first apply a hybrid (mixed deterministic and stochastic) simulation method to such a stochastic model. With proper partitioning of reactions between deterministic and stochastic simulation methods, the hybrid method generates the same primary characteristics and the same level of noise as Gillespie's stochastic simulation algorithm, but with better efficiency. By studying the results generated by various partitionings of reactions, we developed a new strategy for hybrid stochastic modeling of the cell cycle. The new approach is not limited to using mass-action rate laws. Numerical experiments demonstrate that our approach is consistent with characteristics of noisy cell cycle progression, and yields cell cycle statistics in accord with experimental observations
2020-04-14 | MODEL2004140002 | BioModels
Project description:The growth and division of eukaryotic cells are regulated by complex, multi-scale networks. In this process, the mechanism of controlling cell-cycle progression has to be robust against inherent noise in the system. In this paper, a hybrid stochastic model is developed to study the effects of noise on the control mechanism of the budding yeast cell cycle. The modeling approach leverages, in a single multi-scale model, the advantages of two regimes: (1) the computational efficiency of a deterministic approach, and (2) the accuracy of stochastic simulations. Our results show that this hybrid stochastic model achieves high computational efficiency while generating simulation results that match very well with published experimental measurements.
| S-EPMC7101447 | biostudies-literature
Project description:The environment of a living cell is vastly different from that of an in vitro reaction system, an issue that presents great challenges to the use of in vitro models, or computer simulations based on them, for understanding biochemistry in vivo. Virus capsids make an excellent model system for such questions because they typically have few distinct components, making them amenable to in vitro and modeling studies, yet their assembly can involve complex networks of possible reactions that cannot be resolved in detail by any current experimental technology. We previously fit kinetic simulation parameters to bulk in vitro assembly data to yield a close match between simulated and real data, and then used the simulations to study features of assembly that cannot be monitored experimentally. The present work seeks to project how assembly in these simulations fit to in vitro data would be altered by computationally adding features of the cellular environment to the system, specifically the presence of nucleic acid about which many capsids assemble. The major challenge of such work is computational: simulating fine-scale assembly pathways on the scale and in the parameter domains of real viruses is far too computationally costly to allow for explicit models of nucleic acid interaction. We bypass that limitation by applying analytical models of nucleic acid effects to adjust kinetic rate parameters learned from in vitro data to see how these adjustments, singly or in combination, might affect fine-scale assembly progress. The resulting simulations exhibit surprising behavioral complexity, with distinct effects often acting synergistically to drive efficient assembly and alter pathways relative to the in vitro model. The work demonstrates how computer simulations can help us understand how assembly might differ between the in vitro and in vivo environments and what features of the cellular environment account for these differences.
| S-EPMC4887116 | biostudies-literature
Project description:Deterministic models of biochemical processes at the subcellular level might become inadequate when a cascade of chemical reactions is induced by a few molecules. Inherent randomness of such phenomena calls for the use of stochastic simulations. However, being computationally intensive, such simulations become infeasible for large and complex reaction networks. To improve their computational efficiency in handling these networks, we present a hybrid approach, in which slow reactions and fluxes are handled through exact stochastic simulation and their fast counterparts are treated partially deterministically through chemical Langevin equation. The classification of reactions as fast or slow is accompanied by the assumption that in the time-scale of fast reactions, slow reactions do not occur and hence do not affect the probability of the state. Our new approach also handles reactions with complex rate expressions such as Michaelis-Menten kinetics. Fluxes which cannot be modeled explicitly through reactions, such as flux of Ca(2+) from endoplasmic reticulum to the cytosol through inositol 1,4,5-trisphosphate receptor channels, are handled deterministically. The proposed hybrid algorithm is used to model the regulation of the dynamics of cytosolic calcium ions in mouse macrophage RAW 264.7 cells. At relatively large number of molecules, the response characteristics obtained with the stochastic and deterministic simulations coincide, which validates our approach in the limit of large numbers. At low doses, the response characteristics of some key chemical species, such as levels of cytosolic calcium, predicted with stochastic simulations, differ quantitatively from their deterministic counterparts. These observations are ubiquitous throughout dose response, sensitivity, and gene-knockdown response analyses. While the relative differences between the peak-heights of the cytosolic [Ca(2+)] time-courses obtained from stochastic (mean of 16 realizations) and deterministic simulations are merely 1%-4% for most perturbations, it is specially sensitive to levels of G(βγ) (relative difference as large as 90% at very low G(βγ)).
| S-EPMC2998048 | biostudies-other
Project description:Casein kinase 2 (CK2) is a ubiquitous eukaryotic Ser/Thr protein kinase that plays an important role in cell cycle progression. Although its function in this process remains unclear, it is known to be required for the G(1) and G(2)/M phase transitions in yeast. Here, we show that CK2 activity changes notably during cell cycle progression and is increased within 3 h of serum stimulation of quiescent cells. During the time period in which it exhibits high enzymatic activity, CK2 associates with and phosphorylates a key molecule for translation initiation, eukaryotic translation initiation factor (eIF) 5. Using MS, we show that Ser-389 and -390 of eIF5 are major sites of phosphorylation by CK2. This is confirmed using eIF5 mutants that lack CK2 sites; the phosphorylation levels of mutant eIF5 proteins are significantly reduced, relative to WT eIF5, both in vitro and in vivo. Expression of these mutants reveals that they have a dominant-negative effect on phosphorylation of endogenous eIF5, and that they perturb synchronous progression of cells through S to M phase, resulting in a significant reduction in growth rate. Furthermore, the formation of mature eIF5/eIF2/eIF3 complex is reduced in these cells, and, in fact, restricted diffusional motion of WT eIF5 was almost abolished in a GFP-tagged eIF5 mutant lacking CK2 phosphorylation sites, as measured by fluorescence correlation spectroscopy. These results suggest that CK2 may be involved in the regulation of cell cycle progression by associating with and phosphorylating a key molecule for translation initiation.
| S-EPMC1266118 | biostudies-literature
Project description:It is now well-established that the host's adaptive immune system plays an important role in identifying and eliminating cancer cells in much the same way that intracellular pathogens are cleared during an adaptive immune response to infection. From a therapeutic standpoint, the adaptive immune system is unique in that it can co-evolve alongside a developing tumor. Tumor acquisition of immune evasive phenotypes, such as class-I MHC down-regulation, remains a major limitation of successful T-cell immunotherapy. Here, we consider a population dynamical model coupling tumor and adaptive immune compartments in order to study the dynamics and survival of an evolving threat when faced with adaptive immune pressure. We demonstrate that predicted optimal growth strategies depend on whether or not the threat may acquire an immune-evasive phenotype as well as the mode of immune detection. We parameterize adaptive immune functioning by T-cell turnover and repertoire diversity and predict that decreases in the latter quantity which occur in advanced age may substantially affect the ability to recognize, and therefore control, an immune evasive threat like cancer. This framework recapitulates general features of age-dependent AML incidence, thereby providing a probable association between cancer frequency and adaptive immune functioning. Lastly, we quantify therapeutic efficacy of adjuvant immunotherapeutic strategies, and predict their benefits and limitations with regard to handling immune evasion. Our model generates survival behavior consistent with known growth-dependent characteristics, and serves as a first attempt at modeling stochastic cancer evolution alongside an adaptive immune compartment.
| S-EPMC6225903 | biostudies-literature
Project description:DNA replication in eukaryotic cells initiates from hundreds of origins along their genomes, leading to complete duplication of genetic information before cell division. The large number of potential origins, coupled with system uncertainty, dictates the need for new analytical tools to capture spatial and temporal patterns of DNA replication genome-wide. We have developed a stochastic hybrid model that reproduces DNA replication throughout a complete genome. The model can capture different modes of DNA replication and is applicable to various organisms. Using genome-wide data on the location and firing efficiencies of origins in the fission yeast, we show how the DNA replication process evolves during S-phase in the presence of stochastic origin firing. Simulations reveal small regions of the genome that extend S-phase to three times its reported duration. The low levels of late replication predicted by the model are below the detection limit of techniques used to measure S-phase length. Parameter sensitivity analysis shows that increased replication fork speeds genome-wide, or additional origins are not sufficient to reduce S-phase to its reported length. We model the redistribution of a limiting initiation factor during S-phase and show that it could shorten S-phase to the reported duration. Alternatively, S-phase may be extended, and what has traditionally been defined as G2 may be occupied by low levels of DNA synthesis with the onset of mitosis delayed by activation of the G2/M checkpoint.
| S-EPMC2527905 | biostudies-other
Project description:BACKGROUND:The glomerular podocyte is a highly specialized cell type with the ability to ultrafilter blood and support glomerular capillary pressure. However, little is known about the genetic programs leading to this functionality or the final phenotype. RESULTS:In the current study, we found that the expression of a myocardin/MKL family member, MKL1, was significantly upregulated during cell cycle arrest induced by a temperature switch in murine podocyte clone 5 (MPC5) cells. Further investigation demonstrated that overexpression of MKL1 led to inhibition of cell proliferation by decreasing the number of cells in S phase of the cell cycle. In contrast, MKL1 knockdown by RNA interference had the opposite effect, highlighting a potential role of MKL1 in blocking G1/S transition of the cell cycle in MPC5 cells. Additionally, using an RT(2) Profiler PCR Array, p21 was identified as a direct target of MKL1. We further revealed that MKL1 activated p21 transcription by recruitment to the CArG element in its promoter, thus resulting in cell cycle arrest. In addition, the expression of MKL1 is positively correlated with that of p21 in podocytes in postnatal mouse kidney and significantly upregulated during the morphological switch of podocytes from proliferation to differentiation. CONCLUSIONS:Our observations demonstrate that MKL1 has physiological roles in the maturation and development of podocytes, and thus its misregulation might lead to glomerular and renal dysfunction.
| S-EPMC4330937 | biostudies-literature
Project description:Histone H2A variants are highly conserved proteins found ubiquitously in nature and thought to perform specialized functions in the cell. Studies in yeast on the histone H2A variant H2A.Z have shown a role for this protein in transcription as well as chromosome segregation. Our studies have focused on understanding the role of H2A.Z during cell cycle progression. We found that htz1delta cells were delayed in DNA replication and progression through the cell cycle. Furthermore, cells lacking H2A.Z required the S-phase checkpoint pathway for survival. We also found that H2A.Z localized to the promoters of cyclin genes, and cells lacking H2A.Z were delayed in the induction of these cyclin genes. Several different models are proposed to explain these observations.
| S-EPMC1346916 | biostudies-literature
Project description:Mathematical cancer models are immensely powerful tools that are based in part on the fractal nature of biological structures, such as the geometry of the lung. Cancers of the lung provide an opportune model to develop and apply algorithms that capture changes and disease phenotypes. We reviewed mathematical models that have been developed for biological sciences and applied them in the context of small cell lung cancer (SCLC) growth, mutational heterogeneity, and mechanisms of metastasis. The ultimate goal is to develop the stochastic and deterministic nature of this disease, to link this comprehensive set of tools back to its fractalness and to provide a platform for accurate biomarker development. These techniques may be particularly useful in the context of drug development research, such as combination with existing omics approaches. The integration of these tools will be important to further understand the biology of SCLC and ultimately develop novel therapeutics.
| S-EPMC5995226 | biostudies-other