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

The main sources of molecular organization in the cell. Atlas of self-organized and self-regulated dynamic biostructures

One of the most important goals of contemporary biology is to understand the principles of the molecular order underlying the complex dynamic architecture of cells. Here, we present an overview of the main driving forces involved in the cellular molecular complexity and in the emergent functional dynamic structures, spanning from the most basic molecular organization levels to the complex emergent integrative systemic behaviors. First, we address the molecular information processing which is essential in many complex fundamental mechanisms such as the epigenetic memory, alternative splicing, regulation of transcriptional system, and the adequate self-regulatory adaptation to the extracellular environment. Next, we approach the biochemical self-organization, which is central to understand the emergency of metabolic rhythms, circadian oscillations, and spatial traveling waves. Such a complex behavior is also fundamental to understand the temporal compartmentalization of the cellular metabolism and the dynamic regulation of many physiological activities. Numerous examples of biochemical self-organization are considered here, which show that practically all the main physiological processes in the cell exhibit this type of dynamic molecular organization. Finally, we focus on the biochemical self-assembly which, at a primary level of organization, is a basic but important mechanism for the order in the cell allowing biomolecules in a disorganized state to form complex aggregates necessary for a plethora of essential structures and physiological functions. In total, more than 500 references have been compiled in this review. Due to these main sources of order, systemic functional structures emerge in the cell, driving the metabolic functionality towards the biological complexity.

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.

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.

A Stochastic Binary Model for the Regulation of Gene Expression to Investigate Responses to Gene Therapy

In this manuscript, we use an exactly solvable stochastic binary model for the regulation of gene expression to analyze the dynamics of response to a treatment aiming to modulate the number of transcripts of a master regulatory switching gene. The challenge is to combine multiple processes with different time scales to control the treatment response by a switching gene in an unavoidable noisy environment. To establish biologically relevant timescales for the parameters of the model, we select the RKIP gene and two non-specific drugs already known for changing RKIP levels in cancer cells. We demonstrate the usefulness of our method simulating three treatment scenarios aiming to reestablish RKIP gene expression dynamics toward a pre-cancerous state: (1) to increase the promoter's ON state duration; (2) to increase the mRNAs' synthesis rate; and (3) to increase both rates. We show that the pre-treatment kinetic rates of ON and OFF promoter switching speeds and mRNA synthesis and degradation will affect the heterogeneity and time for treatment response. Hence, we present a strategy for reaching increased average mRNA levels with diminished heterogeneity while reducing drug dosage by simultaneously targeting multiple kinetic rates that effectively represent the chemical processes underlying the regulation of gene expression. The decrease in heterogeneity of treatment response by a target gene helps to lower the chances of emergence of resistance. Our approach may be useful for inferring kinetic constants related to the expression of antimetastatic genes or oncogenes and for the design of multi-drug therapeutic strategies targeting the processes underpinning the expression of master regulatory genes.

A computational model of mutual antagonism in the mechano-signaling network of RhoA and nitric oxide

BACKGROUND

RhoA is a master regulator of cytoskeletal contractility, while nitric oxide (NO) is a master regulator of relaxation, e.g., vasodilation. There are multiple forms of cross-talk between the RhoA/ROCK pathway and the eNOS/NO/cGMP pathway, but previous work has not studied their interplay at a systems level. Literature review suggests that the majority of their cross-talk interactions are antagonistic, which motivates us to ask whether the RhoA and NO pathways exhibit mutual antagonism in vitro, and if so, to seek the theoretical implications of their mutual antagonism.

RESULTS

Experiments found mutual antagonism between RhoA and NO in epithelial cells. Since mutual antagonism is a common motif for bistability, we sought to explore through theoretical simulations whether the RhoA-NO network is capable of bistability. Qualitative modeling showed that there are parameters that can cause bistable switching in the RhoA-NO network, and that the robustness of the bistability would be increased by positive feedback between RhoA and mechanical tension.

CONCLUSIONS

We conclude that the RhoA-NO bistability is robust enough in silico to warrant the investment of further experimental testing. Tension-dependent bistability has the potential to create sharp concentration gradients, which could contribute to the localization and self-organization of signaling domains during cytoskeletal remodeling and cell migration.

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.

An Inactivation Switch Enables Rhythms in a Neurospora Clock Model

Autonomous endogenous time-keeping is ubiquitous across many living organisms, known as the circadian clock when it has a period of about 24 h. Interestingly, the fundamental design principle with a network of interconnected negative and positive feedback loops is conserved through evolution, although the molecular components differ. Filamentous fungus Neurospora crassa is a well-established chrono-genetics model organism to investigate the underlying mechanisms. The core negative feedback loop of the clock of Neurospora is composed of the transcription activator White Collar Complex (WCC) (heterodimer of WC1 and WC2) and the inhibitory element called FFC complex, which is made of FRQ (Frequency protein), FRH (Frequency interacting RNA Helicase) and CK1a (Casein kinase 1a). While exploring their temporal dynamics, we investigate how limit cycle oscillations arise and how molecular switches support self-sustained rhythms. We develop a mathematical model of 10 variables with 26 parameters to understand the interactions and feedback among WC1 and FFC elements in nuclear and cytoplasmic compartments. We performed control and bifurcation analysis to show that our novel model produces robust oscillations with a wild-type period of 22.5 h. Our model reveals a switch between WC1-induced transcription and FFC-assisted inactivation of WC1. Using the new model, we also study the possible mechanisms of glucose compensation. A fairly simple model with just three nonlinearities helps to elucidate clock dynamics, revealing a mechanism of rhythms' production. The model can further be utilized to study entrainment and temperature compensation.

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.