A minimal cascade model for the mitotic oscillator involving cyclin and cdc2 kinase

Newton's cradle: Cell cycle regulation by two mutually inhibitory oscillators

The cell division cycle is a fundamental physiological process displaying a great degree of plasticity during the course of multicellular development. This plasticity is evident in the transition from rapid and stringently-timed divisions of the early embryo to subsequent size-controlled mitotic cycles. Later in development, cells may pause and restart proliferation in response to myriads of internal or external signals, or permanently exit the cell cycle following terminal differentiation or senescence. Beyond this, cells can undergo modified cell division variants, such as endoreplication, which increases their ploidy, or meiosis, which reduces their ploidy. This wealth of behaviours has led to numerous conceptual analogies intended as frameworks for understanding the proliferative program. Here, we aim to unify these mechanisms under one dynamical paradigm. To this end, we take a control theoretical approach to frame the cell cycle as a pair of arrestable and mutually-inhibiting, doubly amplified, negative feedback oscillators controlling chromosome replication and segregation events, respectively. Under appropriate conditions, this framework can reproduce fixed-period oscillations, checkpoint arrests of variable duration, and endocycles. Subsequently, we use phase plane and bifurcation analysis to explain the dynamical basis of these properties. Then, using a physiologically realistic, biochemical model, we show that the very same regulatory structure underpins the diverse functions of the cell cycle control network. We conclude that Newton's cradle may be a suitable mechanical analogy of how the cell cycle is regulated.

Curating models from BioModels: Developing a workflow for creating OMEX files

The reproducibility of computational biology models can be greatly facilitated by widely adopted standards and public repositories. We examined 50 models from the BioModels Database and attempted to validate the original curation and correct some of them if necessary. For each model, we reproduced these published results using Tellurium. Once reproduced we manually created a new set of files, with the model information stored by the Systems Biology Markup Language (SBML), and simulation instructions stored by the Simulation Experiment Description Markup Language (SED-ML), and everything included in an Open Modeling EXchange (OMEX) file, which could be used with a variety of simulators to reproduce the same results. On the one hand, the reproducibility procedure of 50 models developed a manual workflow that we would use to build an automatic platform to help users more easily curate and verify models in the future. On the other hand, these exercises allowed us to find the limitations and possible enhancement of the current curation and tooling to verify and curate models.

Deep learning-driven imaging of cell division and cell growth across an entire eukaryotic life cycle

The life cycle of biomedical and agriculturally relevant eukaryotic microorganisms involves complex transitions between proliferative and non-proliferative states such as dormancy, mating, meiosis, and cell division. New drugs, pesticides, and vaccines can be created by targeting specific life cycle stages of parasites and pathogens. However, defining the structure of a microbial life cycle often relies on partial observations that are theoretically assembled in an ideal life cycle path. To create a more quantitative approach to studying complete eukaryotic life cycles, we generated a deep learning-driven imaging framework to track microorganisms across sexually reproducing generations. Our approach combines microfluidic culturing, life cycle stage-specific segmentation of microscopy images using convolutional neural networks, and a novel cell tracking algorithm, FIEST, based on enhancing the overlap of single cell masks in consecutive images through deep learning video frame interpolation. As proof of principle, we used this approach to quantitatively image and compare cell growth and cell cycle regulation across the sexual life cycle of Saccharomyces cerevisiae. We developed a fluorescent reporter system based on a fluorescently labeled Whi5 protein, the yeast analog of mammalian Rb, and a new High-Cdk1 activity sensor, LiCHI, designed to report during DNA replication, mitosis, meiotic homologous recombination, meiosis I, and meiosis II. We found that cell growth preceded the exit from non-proliferative states such as mitotic G1, pre-meiotic G1, and the G0 spore state during germination. A decrease in the total cell concentration of Whi5 characterized the exit from non-proliferative states, which is consistent with a Whi5 dilution model. The nuclear accumulation of Whi5 was developmentally regulated, being at its highest during meiotic exit and spore formation. The temporal coordination of cell division and growth was not significantly different across three sexually reproducing generations. Our framework could be used to quantitatively characterize other single-cell eukaryotic life cycles that remain incompletely described. An off-the-shelf user interface Yeastvision provides free access to our image processing and single-cell tracking algorithms.

A yeast cell cycle pulse generator model shows consistency with multiple oscillatory and checkpoint mutant datasets

Modeling biological systems holds great promise for speeding up the rate of discovery in systems biology by predicting experimental outcomes and suggesting targeted interventions. However, this process is dogged by an identifiability issue, in which network models and their parameters are not sufficiently constrained by coarse and noisy data to ensure unique solutions. In this work, we evaluated the capability of a simplified yeast cell-cycle network model to reproduce multiple observed transcriptomic behaviors under genomic mutations. We matched time-series data from both cycling and checkpoint arrested cells to model predictions using an asynchronous multi-level Boolean approach. We showed that this single network model, despite its simplicity, is capable of exhibiting dynamical behavior similar to the datasets in most cases, and we demonstrated the drop in severity of the identifiability issue that results from matching multiple datasets.

An oscillating reaction network with an exact closed form solution in the time domain

BACKGROUND

Oscillatory behavior is critical to many life sustaining processes such as cell cycles, circadian rhythms, and notch signaling. Important biological functions depend on the characteristics of these oscillations (hereafter, oscillation characteristics or OCs): frequency (e.g., event timings), amplitude (e.g., signal strength), and phase (e.g., event sequencing). Numerous oscillating reaction networks have been documented or proposed. Some investigators claim that oscillations in reaction networks require nonlinear dynamics in that at least one rate law is a nonlinear function of species concentrations. No one has shown that oscillations can be produced for a reaction network with linear dynamics. Further, no one has obtained closed form solutions for the frequency, amplitude and phase of any oscillating reaction network. Finally, no one has published an algorithm for constructing oscillating reaction networks with desired OCs.

RESULTS

This is a theoretical study that analyzes reaction networks in terms of their representation as systems of ordinary differential equations. Our contributions are: (a) construction of an oscillating, two species reaction network [two species harmonic oscillator (2SHO)] that has no nonlinearity; (b) obtaining closed form formulas that calculate frequency, amplitude, and phase in terms of the parameters of the 2SHO reaction network, something that has not been done for any published oscillating reaction network; and (c) development of an algorithm that parameterizes the 2SHO to achieve desired oscillation, a capability that has not been produced for any published oscillating reaction network.

CONCLUSIONS

Our 2SHO demonstrates the feasibility of creating an oscillating reaction network whose dynamics are described by a system of linear differential equations. Because it is a linear system, we can derive closed form expressions for the frequency, amplitude, and phase of oscillations, something that has not been done for other published reaction networks. With these formulas, we can design 2SHO reaction networks to have desired oscillation characteristics. Finally, our sensitivity analysis suggests an approach to constructing a 2SHO for a biochemical system.

The ups and downs of biological oscillators: a comparison of time-delayed negative feedback mechanisms

Many biochemical oscillators are driven by the periodic rise and fall of protein concentrations or activities. A negative feedback loop underlies such oscillations. The feedback can act on different parts of the biochemical network. Here, we mathematically compare time-delay models where the feedback affects production and degradation. We show a mathematical connection between the linear stability of the two models, and derive how both mechanisms impose different constraints on the production and degradation rates that allow oscillations. We show how oscillations are affected by the inclusion of a distributed delay, of double regulation (acting on production and degradation) and of enzymatic degradation.

Sloppiness: Fundamental study, new formalism and its application in model assessment

Computational modelling of biological processes poses multiple challenges in each stage of the modelling exercise. Some significant challenges include identifiability, precisely estimating parameters from limited data, informative experiments and anisotropic sensitivity in the parameter space. One of these challenges' crucial but inconspicuous sources is the possible presence of large regions in the parameter space over which model predictions are nearly identical. This property, known as sloppiness, has been reasonably well-addressed in the past decade, studying its possible impacts and remedies. However, certain critical unanswered questions concerning sloppiness, particularly related to its quantification and practical implications in various stages of system identification, still prevail. In this work, we systematically examine sloppiness at a fundamental level and formalise two new theoretical definitions of sloppiness. Using the proposed definitions, we establish a mathematical relationship between the parameter estimates' precision and sloppiness in linear predictors. Further, we develop a novel computational method and a visual tool to assess the goodness of a model around a point in parameter space by identifying local structural identifiability and sloppiness and finding the most sensitive and least sensitive parameters for non-infinitesimal perturbations. We demonstrate the working of our method in benchmark systems biology models of various complexities. The pharmacokinetic HIV infection model analysis identified a new set of biologically relevant parameters that can be used to control the free virus in an active HIV infection.

Learning in Transcriptional Network Models: Computational Discovery of Pathway-Level Memory and Effective Interventions

Trainability, in any substrate, refers to the ability to change future behavior based on past experiences. An understanding of such capacity within biological cells and tissues would enable a particularly powerful set of methods for prediction and control of their behavior through specific patterns of stimuli. This top-down mode of control (as an alternative to bottom-up modification of hardware) has been extensively exploited by computer science and the behavioral sciences; in biology however, it is usually reserved for organism-level behavior in animals with brains, such as training animals towards a desired response. Exciting work in the field of basal cognition has begun to reveal degrees and forms of unconventional memory in non-neural tissues and even in subcellular biochemical dynamics. Here, we characterize biological gene regulatory circuit models and protein pathways and find them capable of several different kinds of memory. We extend prior results on learning in binary transcriptional networks to continuous models and identify specific interventions (regimes of stimulation, as opposed to network rewiring) that abolish undesirable network behavior such as drug pharmacoresistance and drug sensitization. We also explore the stability of created memories by assessing their long-term behavior and find that most memories do not decay over long time periods. Additionally, we find that the memory properties are quite robust to noise; surprisingly, in many cases noise actually increases memory potential. We examine various network properties associated with these behaviors and find that no one network property is indicative of memory. Random networks do not show similar memory behavior as models of biological processes, indicating that generic network dynamics are not solely responsible for trainability. Rational control of dynamic pathway function using stimuli derived from computational models opens the door to empirical studies of proto-cognitive capacities in unconventional embodiments and suggests numerous possible applications in biomedicine, where behavior shaping of pathway responses stand as a potential alternative to gene therapy.

Modeling signaling pathways in biology with MaBoSS: From one single cell to a dynamic population of heterogeneous interacting cells

As a result of the development of experimental technologies and the accumulation of data, biological and molecular processes can be described as complex networks of signaling pathways. These networks are often directed and signed, where nodes represent entities (genes/proteins) and arrows interactions. They are translated into mathematical models by adding a dynamic layer onto them. Such mathematical models help to understand and interpret non-intuitive experimental observations and to anticipate the response to external interventions such as drug effects on phenotypes. Several frameworks for modeling signaling pathways exist. The choice of the appropriate framework is often driven by the experimental context. In this review, we present MaBoSS, a tool based on Boolean modeling using a continuous time approach, which predicts time-dependent probabilities of entities in different biological contexts. MaBoSS was initially built to model the intracellular signaling in non-interacting homogeneous cell populations. MaBoSS was then adapted to model heterogeneous cell populations (EnsembleMaBoSS) by considering families of models rather than a unique model. To account for more complex questions, MaBoSS was extended to simulate dynamical interacting populations (UPMaBoSS), with a precise spatial distribution (PhysiBoSS). To illustrate all these levels of description, we show how each of these tools can be used with a running example of a simple model of cell fate decisions. Finally, we present practical applications to cancer biology and studies of the immune response.

Computational modeling of protracted HCMV replication using genome substrates and protein temporal profiles

Human cytomegalovirus (HCMV) is a major cause of illness in immunocompromised individuals. The HCMV lytic cycle contributes to the clinical manifestations of infection. The lytic cycle occurs over ∼96 h in diverse cell types and consists of viral DNA (vDNA) genome replication and temporally distinct expression of hundreds of viral proteins. Given its complexity, understanding this elaborate system can be facilitated by the introduction of mechanistic computational modeling of temporal relationships. Therefore, we developed a multiplicity of infection (MOI)-dependent mechanistic computational model that simulates vDNA kinetics and late lytic replication based on in-house experimental data. The predictive capabilities were established by comparison to post hoc experimental data. Computational analysis of combinatorial regulatory mechanisms suggests increasing rates of protein degradation in association with increasing vDNA levels. The model framework also allows expansion to account for additional mechanisms regulating the processes. Simulating vDNA kinetics and the late lytic cycle for a wide range of MOIs yielded several unique observations. These include the presence of saturation behavior at high MOIs, inefficient replication at low MOIs, and a precise range of MOIs in which virus is maximized within a cell type, being 0.382 IU to 0.688 IU per fibroblast. The predicted saturation kinetics at high MOIs are likely related to the physical limitations of cellular machinery, while inefficient replication at low MOIs may indicate a minimum input material required to facilitate infection. In summary, we have developed and demonstrated the utility of a data-driven and expandable computational model simulating lytic HCMV infection.

A Model for the Proliferation–Quiescence Transition in Human Cells

The process of revitalising quiescent cells in order for them to proliferate plays a pivotal role in the repair of worn-out tissues as well as for tissue homeostasis. This process is also crucial in the growth, development and well-being of higher multi-cellular organisms such as mammals. Deregulation of proliferation-quiescence transition is related to many diseases, such as cancer. Recent studies have revealed that this proliferation–quiescence process is regulated tightly by the Rb−E2F bistable switch mechanism. Based on experimental observations, in this study, we formulate a mathematical model to examine the effect of the growth factor concentration on the proliferation–quiescence transition in human cells. Working with a non-dimensionalised model, we prove the positivity, boundedness and uniqueness of solutions. To understand model solution behaviour close to bifurcation points, we carry out bifurcation analysis, which is further illustrated by the use of numerical bifurcation analysis, sensitivity analysis and numerical simulations. Indeed, bifurcation and numerical analysis of the model predicted a transition between bistable and stable states, which are dependent on the growth factor concentration parameter (GF). The derived predictions confirm experimental observations.

Time-keeping and decision-making in the cell cycle

Cell growth, DNA replication, mitosis and division are the fundamental processes by which life is passed on from one generation of eukaryotic cells to the next. The eukaryotic cell cycle is intrinsically a periodic process but not so much a ‘clock’ as a ‘copy machine’, making new daughter cells as warranted. Cells growing under ideal conditions divide with clock-like regularity; however, if they are challenged with DNA-damaging agents or mitotic spindle disrupters, they will not progress to the next stage of the cycle until the damage is repaired. These ‘decisions’ (to exit and re-enter the cell cycle) are essential to maintain the integrity of the genome from generation to generation. A crucial challenge for molecular cell biologists in the 1990s was to unravel the genetic and biochemical mechanisms of cell cycle control in eukaryotes. Central to this effort were biochemical studies of the clock-like regulation of ‘mitosis promoting factor’ during synchronous mitotic cycles of fertilized frog eggs and genetic studies of the switch-like regulation of ‘cyclin-dependent kinases' in yeast cells. In this review, we uncover some secrets of cell cycle regulation by mathematical modelling of increasingly more complex molecular regulatory networks of cell cycle ‘clocks’ and ‘switches’.

A yeast cell cycle model integrating stress, signaling, and physiology

The cell division cycle in eukaryotic cells is a series of highly coordinated molecular interactions that ensure that cell growth, duplication of genetic material, and actual cell division are precisely orchestrated to give rise to two viable progeny cells. Moreover, the cell cycle machinery is responsible for incorporating information about external cues or internal processes that the cell must keep track of to ensure a coordinated, timely progression of all related processes. This is most pronounced in multicellular organisms, but also a cardinal feature in model organisms such as baker's yeast. The complex and integrative behavior is difficult to grasp and requires mathematical modeling to fully understand the quantitative interplay of the single components within the entire system. Here, we present a self-oscillating mathematical model of the yeast cell cycle that comprises all major cyclins and their main regulators. Furthermore, it accounts for the regulation of the cell cycle machinery by a series of external stimuli such as mating pheromones and changes in osmotic pressure or nutrient quality. We demonstrate how the external perturbations modify the dynamics of cell cycle components and how the cell cycle resumes after adaptation to or relief from stress.

Application of Lie Symmetry to a Mathematical Model that Describes a Cancer Sub-Network

In this paper, a mathematical model of a cancer sub-network is analysed from the view point of Lie symmetry methods. This model discusses a human cancer cell which is developed due to the dysfunction of some genes at the R-checkpoint during the cell cycle. The primary purpose of this paper is to apply the techniques of Lie symmetry to the model and present some approximated solutions for the three-dimensional system of first-order ordinary differential equations describing a cancer sub-network. The result shows that the phosphatase gene (Cdc25A) regulates the cyclin-dependent kinases inhibitor (P27Kip1). Furthermore, this research discovered that the activity that reverses the inhibitory effects on cell cycle progression at the R-checkpoint initiates a pathway.

Cyclin/Forkhead-mediated coordination of cyclin waves: an autonomous oscillator rationalizing the quantitative model of Cdk control for budding yeast

Networks of interacting molecules organize topology, amount, and timing of biological functions. Systems biology concepts required to pin down 'network motifs' or 'design principles' for time-dependent processes have been developed for the cell division cycle, through integration of predictive computer modeling with quantitative experimentation. A dynamic coordination of sequential waves of cyclin-dependent kinases (cyclin/Cdk) with the transcription factors network offers insights to investigate how incompatible processes are kept separate in time during the eukaryotic cell cycle. Here this coordination is discussed for the Forkhead transcription factors in light of missing gaps in the current knowledge of cell cycle control in budding yeast. An emergent design principle is proposed where cyclin waves are synchronized by a cyclin/Cdk-mediated feed-forward regulation through the Forkhead as a transcriptional timer. This design is rationalized by the bidirectional interaction between mitotic cyclins and the Forkhead transcriptional timer, resulting in an autonomous oscillator that may be instrumental for a well-timed progression throughout the cell cycle. The regulation centered around the cyclin/Cdk-Forkhead axis can be pivotal to timely coordinate cell cycle dynamics, thereby to actuate the quantitative model of Cdk control.

Crosstalk between Plk1, p53, cell cycle, and G2/M DNA damage checkpoint regulation in cancer: computational modeling and analysis

Different cancer cell lines can have varying responses to the same perturbations or stressful conditions. Cancer cells that have DNA damage checkpoint-related mutations are often more sensitive to gene perturbations including altered Plk1 and p53 activities than cancer cells without these mutations. The perturbations often induce a cell cycle arrest in the former cancer, whereas they only delay the cell cycle progression in the latter cancer. To study crosstalk between Plk1, p53, and G2/M DNA damage checkpoint leading to differential cell cycle regulations, we developed a computational model by extending our recently developed model of mitotic cell cycle and including these key interactions. We have used the model to analyze the cancer cell cycle progression under various gene perturbations including Plk1-depletion conditions. We also analyzed mutations and perturbations in approximately 1800 different cell lines available in the Cancer Dependency Map and grouped lines by genes that are represented in our model. Our model successfully explained phenotypes of various cancer cell lines under different gene perturbations. Several sensitivity analysis approaches were used to identify the range of key parameter values that lead to the cell cycle arrest in cancer cells. Our resulting model can be used to predict the effect of potential treatments targeting key mitotic and DNA damage checkpoint regulators on cell cycle progression of different types of cancer cells.

Resonating neurons stabilize heterogeneous grid-cell networks

A central theme that governs the functional design of biological networks is their ability to sustain stable function despite widespread parametric variability. Here, we investigated the impact of distinct forms of biological heterogeneities on the stability of a two-dimensional continuous attractor network (CAN) implicated in grid-patterned activity generation. We show that increasing degrees of biological heterogeneities progressively disrupted the emergence of grid-patterned activity and resulted in progressively large perturbations in low-frequency neural activity. We postulated that targeted suppression of low-frequency perturbations could ameliorate heterogeneity-induced disruptions of grid-patterned activity. To test this, we introduced intrinsic resonance, a physiological mechanism to suppress low-frequency activity, either by adding an additional high-pass filter (phenomenological) or by incorporating a slow negative feedback loop (mechanistic) into our model neurons. Strikingly, CAN models with resonating neurons were resilient to the incorporation of heterogeneities and exhibited stable grid-patterned firing. We found CAN models with mechanistic resonators to be more effective in targeted suppression of low-frequency activity, with the slow kinetics of the negative feedback loop essential in stabilizing these networks. As low-frequency perturbations (1/f noise) are pervasive across biological systems, our analyses suggest a universal role for mechanisms that suppress low-frequency activity in stabilizing heterogeneous biological networks.

Mitotic waves in an import-diffusion model with multiple nuclei in a shared cytoplasm

Nuclei import and export proteins, including cell cycle regulators. These import-export processes are modulated periodically by the cell cycle, for example due to the periodic assembly and breakdown of the nuclear envelope. As such, replicated DNA can be segregated between the two daughter cells and the proteins that were localized in the nucleus are free to diffuse throughout the cytoplasm. Here, we study a mathematical import-diffusion model to show how proteins, i.e. cell cycle regulators, could be redistributed in the cytoplasm by nuclei that periodically toggle between interphase and mitosis. We show that when the cell cycle period depends on the local concentration of regulators, the model exhibits mitotic waves. We discuss how the velocity and spatial origin of these mitotic waves depend on the different model parameters. This work is motivated by recent in vitro experiments reporting on mitotic waves in cycling cell-free extracts made with Xenopus laevis frog eggs, where multiple nuclei share the same cytoplasm. Such experiments have shown that nuclei act as pacemakers for the cell cycle and thus play an important role in collectively defining the spatial origin of mitotic waves.

Quantitative model of eukaryotic Cdk control through the Forkhead CONTROLLER

In budding yeast, synchronization of waves of mitotic cyclins that activate the Cdk1 kinase occur through Forkhead transcription factors. These molecules act as controllers of their sequential order and may account for the separation in time of incompatible processes. Here, a Forkhead-mediated design principle underlying the quantitative model of Cdk control is proposed for budding yeast. This design rationalizes timing of cell division, through progressive and coordinated cyclin/Cdk-mediated phosphorylation of Forkhead, and autonomous cyclin/Cdk oscillations. A "clock unit" incorporating this design that regulates timing of cell division is proposed for both yeast and mammals, and has a DRIVER operating the incompatible processes that is instructed by multiple CLOCKS. TIMERS determine whether the clocks are active, whereas CONTROLLERS determine how quickly the clocks shall function depending on external MODULATORS. This "clock unit" may coordinate temporal waves of cyclin/Cdk concentration/activity in the eukaryotic cell cycle making the driver operate the incompatible processes, at separate times.

Finding new edges: systems approaches to MTOR signaling

Cells have evolved highly intertwined kinase networks to finely tune cellular homeostasis to the environment. The network converging on the mechanistic target of rapamycin (MTOR) kinase constitutes a central hub that integrates metabolic signals and adapts cellular metabolism and functions to nutritional changes and stress. Feedforward and feedback loops, crosstalks and a plethora of modulators finely balance MTOR-driven anabolic and catabolic processes. This complexity renders it difficult — if not impossible — to intuitively decipher signaling dynamics and network topology. Over the last two decades, systems approaches have emerged as powerful tools to simulate signaling network dynamics and responses. In this review, we discuss the contribution of systems studies to the discovery of novel edges and modulators in the MTOR network in healthy cells and in disease.

Dynamic bistable switches enhance robustness and accuracy of cell cycle transitions

Bistability is a common mechanism to ensure robust and irreversible cell cycle transitions. Whenever biological parameters or external conditions change such that a threshold is crossed, the system abruptly switches between different cell cycle states. Experimental studies have uncovered mechanisms that can make the shape of the bistable response curve change dynamically in time. Here, we show how such a dynamically changing bistable switch can provide a cell with better control over the timing of cell cycle transitions. Moreover, cell cycle oscillations built on bistable switches are more robust when the bistability is modulated in time. Our results are not specific to cell cycle models and may apply to other bistable systems in which the bistable response curve is time-dependent.

Many systems in nature show bistability, which means they can evolve to one of two stable steady states under exactly the same conditions. Which state they evolve to depends on where the system comes from. Such bistability underlies the switching behavior that is essential for cells to progress in the cell division cycle. A quick switch happens when the cell jumps from one steady state to another steady state. Typical of this switching behavior is its robustness and irreversibility. In this paper, we expand this viewpoint of the dynamics of the cell cycle by considering bistable switches which themselves are changing in time. This gives the cell an extra layer of control over transitions both in time and in space, and can make those transitions more robust. Such dynamically changing bistability can appear very naturally. We show this in a model of mitotic entry, in which we include a nuclear and cytoplasmic compartment. The activity of a crucial cell cycle protein follows a bistable switch in each compartment, but the shape of its response is changing in time as proteins are imported into and exported from the nucleus.

Computer-Aided Whole-Cell Design: Taking a Holistic Approach by Integrating Synthetic With Systems Biology

Computer-aided design (CAD) for synthetic biology promises to accelerate the rational and robust engineering of biological systems. It requires both detailed and quantitative mathematical and experimental models of the processes to (re)design biology, and software and tools for genetic engineering and DNA assembly. Ultimately, the increased precision in the design phase will have a dramatic impact on the production of designer cells and organisms with bespoke functions and increased modularity. CAD strategies require quantitative models of cells that can capture multiscale processes and link genotypes to phenotypes. Here, we present a perspective on how whole-cell, multiscale models could transform design-build-test-learn cycles in synthetic biology. We show how these models could significantly aid in the design and learn phases while reducing experimental testing by presenting case studies spanning from genome minimization to cell-free systems. We also discuss several challenges for the realization of our vision. The possibility to describe and build whole-cells in silico offers an opportunity to develop increasingly automatized, precise and accessible CAD tools and strategies.

Towards Design of Self-Organizing Biomimetic Systems

The ability of designing biosynthetic systems with well-defined functional biomodules from scratch is an ambitious and revolutionary goal to deliver innovative, engineered solutions to future challenges in biotechnology and process systems engineering. In this work, several key challenges including modularization, functional biomodule identification, and assembly are discussed. In addition, an in silico protocell modeling approach is presented as a foundation for a computational model-based toolkit for rational analysis and modular design of biomimetic systems.

Network mechanisms and dysfunction within an integrated computational model of progression through mitosis in the human cell cycle

The cellular protein-protein interaction network that governs cellular proliferation (cell cycle) is highly complex. Here, we have developed a novel computational model of human mitotic cell cycle, integrating diverse cellular mechanisms, for the purpose of generating new hypotheses and predicting new experiments designed to help understand complex diseases. The pathogenic state investigated is infection by a human herpesvirus. The model starts at mitotic entry initiated by the activities of Cyclin-dependent kinase 1 (CDK1) and Polo-like kinase 1 (PLK1), transitions through Anaphase-promoting complex (APC/C) bound to Cell division cycle protein 20 (CDC20), and ends upon mitotic exit mediated by APC/C bound to CDC20 homolog 1 (CDH1). It includes syntheses and multiple mechanisms of degradations of the mitotic proteins. Prior to this work, no such comprehensive model of the human mitotic cell cycle existed. The new model is based on a hybrid framework combining Michaelis-Menten and mass action kinetics for the mitotic interacting reactions. It simulates temporal changes in 12 different mitotic proteins and associated protein complexes in multiple states using 15 interacting reactions and 26 ordinary differential equations. We have defined model parameter values using both quantitative and qualitative data and using parameter values from relevant published models, and we have tested the model to reproduce the cardinal features of human mitosis determined experimentally by numerous laboratories. Like cancer, viruses create dysfunction to support infection. By simulating infection of the human herpesvirus, cytomegalovirus, we hypothesize that virus-mediated disruption of APC/C is necessary to establish a unique mitotic collapse with sustained CDK1 activity, consistent with known mechanisms of virus egress. With the rapid discovery of cellular protein-protein interaction networks and regulatory mechanisms, we anticipate that this model will be highly valuable in helping us to understand the network dynamics and identify potential points of therapeutic interventions.

EasyModel: user-friendly tool for building and analysis of simple mathematical models in systems biology

SUMMARY

EasyModel is a new user-friendly web application that contains ready-for-simulation versions of the BioModels Database, and allows for the intuitive creation of new models. Its main target audience is the experimental biologist and students of bioinformatics or systems biology without programming skills. Expert users can also benefit from it by implementing basic models quickly and downloading the code for further tailoring.

AVAILABILITY AND IMPLEMENTATION

Freely available on the web at https://easymodel.udl.cat. Implementation is described in its own section.

Inferring Leading Interactions in the p53/Mdm2/Mdmx Circuit through Live-Cell Imaging and Modeling

The tumor-suppressive transcription factor p53 is a master regulator of stress responses. In non-stressed conditions, p53 is maintained at low levels by the ubiquitin ligase Mdm2 and its binding partner Mdmx. Mdmx depletion leads to a biphasic p53 response, with an initial post-mitotic pulse followed by oscillations. The mechanism underlying this dynamical behavior is unknown. Two different roles for Mdmx have been proposed: enhancing p53 ubiquitination by Mdm2 and inhibiting p53 activity. Here, we developed a mathematical model of the p53/Mdm2/Mdmx network to investigate which Mdmx functions quantitatively affect specific features of p53 dynamics under various conditions. We found that enhancement of Mdm2 activity was the most critical role of Mdmx under unstressed conditions. The model also accurately predicted p53 dynamics in Mdmx-depleted cells following DNA damage. This work outlines a strategy for rapidly testing possible network interactions to reveal those most impactful in regulating the dynamics of key proteins.

Mathematical Model Explaining the Role of CDC6 in the Diauxic Growth of CDK1 Activity during the M-Phase of the Cell Cycle

In this paper we propose a role for the CDC 6 protein in the entry of cells into mitosis. This has not been considered in the literature so far. Recent experiments suggest that CDC 6 , upon entry into mitosis, inhibits the appearance of active CDK 1 and cyclin B complexes. This paper proposes a mathematical model which incorporates the dynamics of kinase CDK 1 , its regulatory protein cyclin B, the regulatory phosphatase CDC 25 and the inhibitor CDC 6 known to be involved in the regulation of active CDK 1 and cyclin B complexes. The experimental data lead us to formulate a new hypothesis that CDC 6 slows down the activation of inactive complexes of CDK 1 and cyclin B upon mitotic entry. Our mathematical model, based on mass action kinetics, provides a possible explanation for the experimental data. We claim that the dynamics of active complexes CDK 1 and cyclin B have a similar nature to diauxic dynamics introduced by Monod in 1949. In mathematical terms we state it as the existence of more than one inflection point of the curve defining the dynamics of the complexes.

Quantitative Studies for Cell-Division Cycle Control

The cell-division cycle (CDC) is driven by cyclin-dependent kinases (CDKs). Mathematical models based on molecular networks, as revealed by molecular and genetic studies, have reproduced the oscillatory behavior of CDK activity. Thus, one basic system for representing the CDC is a biochemical oscillator (CDK oscillator). However, genetically clonal cells divide with marked variability in their total duration of a single CDC round, exhibiting non-Gaussian statistical distributions. Therefore, the CDK oscillator model does not account for the statistical nature of cell-cycle control. Herein, we review quantitative studies of the statistical properties of the CDC. Over the past 70 years, studies have shown that the CDC is driven by a cluster of molecular oscillators. The CDK oscillator is coupled to transcriptional and mitochondrial metabolic oscillators, which cause deterministic chaotic dynamics for the CDC. Recent studies in animal embryos have raised the possibility that the dynamics of molecular oscillators underlying CDC control are affected by allometric volume scaling among the cellular compartments. Considering these studies, we discuss the idea that a cluster of molecular oscillators embedded in different cellular compartments coordinates cellular physiology and geometry for successful cell divisions.

Signaling pathways in context

The last decade has seen a rise in the development of methods and models to analyze cellular networks on all levels. The applications of this knowledge are, however, often confined to specifics of the network in concrete conditions and leveraging it is hampered by the lack of information about this context and its implications on the system. While not all cellular networks have been deciphered yet, even for well-studied networks their versatility in different contexts is barely considered. Here, we focus on challenges and potentials when integrating signaling networks into their encompassing structures. We highlight three different consequences of this process: a) its fundamental importance for whole-cell and large-scale models, b) significant changes in contextual behavior imposed on entire systems by genetic variations, and c) species-specific conservation or divergence of signaling motifs can give important clues on how to handle cellular context. While important studies have been conducted on these topics to some extent, an increased focus on developing and exploiting solutions for integrative contextualization should turn out as a fruitful path for both theoretical and experimental research.

A comprehensive, mechanistically detailed, and executable model of the cell division cycle in Saccharomyces cerevisiae

Understanding how cellular functions emerge from the underlying molecular mechanisms is a key challenge in biology. This will require computational models, whose predictive power is expected to increase with coverage and precision of formulation. Genome-scale models revolutionised the metabolic field and made the first whole-cell model possible. However, the lack of genome-scale models of signalling networks blocks the development of eukaryotic whole-cell models. Here, we present a comprehensive mechanistic model of the molecular network that controls the cell division cycle in Saccharomyces cerevisiae. We use rxncon, the reaction-contingency language, to neutralise the scalability issues preventing formulation, visualisation and simulation of signalling networks at the genome-scale. We use parameter-free modelling to validate the network and to predict genotype-to-phenotype relationships down to residue resolution. This mechanistic genome-scale model offers a new perspective on eukaryotic cell cycle control, and opens up for similar models-and eventually whole-cell models-of human cells.

Application of Parameter Optimization to Search for Oscillatory Mass-Action Networks Using Python

Biological systems can be described mathematically to model the dynamics of metabolic, protein, or gene-regulatory networks, but locating parameter regimes that induce a particular dynamic behavior can be challenging due to the vast parameter landscape, particularly in large models. In the current work, a Pythonic implementation of existing bifurcation objective functions, which reward systems that achieve a desired bifurcation behavior, is implemented to search for parameter regimes that permit oscillations or bistability. A differential evolution algorithm progressively approximates the specified bifurcation type while performing a global search of parameter space for a candidate with the best fitness. The user-friendly format facilitates integration with systems biology tools, as Python is a ubiquitous programming language. The bifurcation–evolution software is validated on published models from the BioModels Database and used to search populations of randomly-generated mass-action networks for oscillatory dynamics. Results of this search demonstrate the importance of reaction enrichment to provide flexibility and enable complex dynamic behaviors, and illustrate the role of negative feedback and time delays in generating oscillatory dynamics.

A Pycellerator Tutorial

We present a tutorial on using Pycellerator for biomolecular simulations. Models are described in human readable (and editable) text files (UTF8 or ASCII) containing collections of reactions, assignments, initial conditions, function definitions, and rate constants. These models are then converted into a Python program that can optionally solve the system, e.g., as a system of differential equations using ODEINT, or be run by another program. The input language implements an extended version of the Cellerator arrow notation, including mass action, Hill functions, S-Systems, MWC, and reactions with user-defined kinetic laws. Simple flux balance analysis is also implemented. We will demonstrate the implementation and analysis of progressively more complex models, starting from simple mass action through indexed cascades. Pycellerator can be used as a library that is integrated into other programs, run as a command line program, or in iPython notebooks. It is implemented in Python 2.7 and available under an open source GPL license.

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.

Dissipative structures in biological systems: bistability, oscillations, spatial patterns and waves

The goal of this review article is to assess how relevant is the concept of dissipative structure for understanding the dynamical bases of non-equilibrium self-organization in biological systems, and to see where it has been applied in the five decades since it was initially proposed by Ilya Prigogine. Dissipative structures can be classified into four types, which will be considered, in turn, and illustrated by biological examples: (i) multistability, in the form of bistability and tristability, which involve the coexistence of two or three stable steady states, or in the form of birhythmicity, which involves the coexistence between two stable rhythms; (ii) temporal dissipative structures in the form of sustained oscillations, illustrated by biological rhythms; (iii) spatial dissipative structures, known as Turing patterns; and (iv) spatio-temporal structures in the form of propagating waves. Rhythms occur with widely different periods at all levels of biological organization, from neural, cardiac and metabolic oscillations to circadian clocks and the cell cycle; they play key roles in physiology and in many disorders. New rhythms are being uncovered while artificial ones are produced by synthetic biology. Rhythms provide the richest source of examples of dissipative structures in biological systems. Bistability has been observed experimentally, but has primarily been investigated in theoretical models in an increasingly wide range of biological contexts, from the genetic to the cell and animal population levels, both in physiological conditions and in disease. Bistable transitions have been implicated in the progression between the different phases of the cell cycle and, more generally, in the process of cell fate specification in the developing embryo. Turing patterns are exemplified by the formation of some periodic structures in the course of development and by skin stripe patterns in animals. Spatio-temporal patterns in the form of propagating waves are observed within cells as well as in intercellular communication. This review illustrates how dissipative structures of all sorts abound in biological systems.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 1)'.

A minimal "push-pull" bistability model explains oscillations between quiescent and proliferative cell states

A minimal model for oscillating between quiescent and growth/proliferation states, dependent on the availability of a central metabolic resource, is presented. From the yeast metabolic cycles, metabolic oscillations in oxygen consumption are represented as transitions between quiescent and growth states. We consider metabolic resource availability, growth rates, and switching rates (between states) to model a relaxation oscillator explaining transitions between these states. This frustrated bistability model reveals a required communication between the metabolic resource that determines oscillations and the quiescent and growth state cells. Cells in each state reflect memory, or hysteresis of their current state, and “push–pull” cells from the other state. Finally, a parsimonious argument is made for a specific central metabolite as the controller of switching between quiescence and growth states. We discuss how an oscillator built around the availability of such a metabolic resource is sufficient to generally regulate oscillations between growth and quiescence through committed transitions.

Biological Oscillators in Nanonetworks-Opportunities and Challenges

One of the major issues in molecular communication-based nanonetworks is the provision and maintenance of a common time knowledge. To stay true to the definition of molecular communication, biological oscillators are the potential solutions to achieve that goal as they generate oscillations through periodic fluctuations in the concentrations of molecules. Through the lens of a communication systems engineer, the scope of this survey is to explicitly classify, for the first time, existing biological oscillators based on whether they are found in nature or not, to discuss, in a tutorial fashion, the main principles that govern the oscillations in each oscillator, and to analyze oscillator parameters that are most relevant to communication engineer researchers. In addition, the survey highlights and addresses the key open research issues pertaining to several physical aspects of the oscillators and the adoption and implementation of the oscillators to nanonetworks. Moreover, key research directions are discussed.

Modeling-Based Investigation of the Effect of Noise in Cellular Systems

Noise is pervasive in cellular biology and inevitably affects the dynamics of cellular processes. Biological systems have developed regulatory mechanisms to ensure robustness with respect to noise or to take advantage of stochasticity. We review here, through a couple of selected examples, some insights on possible robustness factors and constructive roles of noise provided by computational modeling. In particular, we focus on (1) factors that likely contribute to the robustness of oscillatory processes such as the circadian clocks and the cell cycle, (2) how reliable coding/decoding of calcium-mediated signaling could be achieved in presence of noise and, in some cases, enhanced through stochastic resonance, and (3) how embryonic cell differentiation processes can exploit stochasticity to create heterogeneity in a population of identical cells.

Delay models for the early embryonic cell cycle oscillator

Time delays are known to play a crucial role in generating biological oscillations. The early embryonic cell cycle in the frog Xenopus laevis is one such example. Although various mathematical models of this oscillating system exist, it is not clear how to best model the required time delay. Here, we study a simple cell cycle model that produces oscillations due to the presence of an ultrasensitive, time-delayed negative feedback loop. We implement the time delay in three qualitatively different ways, using a fixed time delay, a distribution of time delays, and a delay that is state-dependent. We analyze the dynamics in all cases, and we use experimental observations to interpret our results and put constraints on unknown parameters. In doing so, we find that different implementations of the time delay can have a large impact on the resulting oscillations.

Parameter-robustness analysis for a biochemical oscillator model describing the social-behaviour transition phase of myxobacteria

We develop a tool based on bifurcation analysis for parameter-robustness analysis for a class of oscillators and, in particular, examine a biochemical oscillator that describes the transition phase between social behaviours of myxobacteria. Myxobacteria are a particular group of soil bacteria that have two dogmatically different types of social behaviour: when food is abundant they live fairly isolated forming swarms, but when food is scarce, they aggregate into a multicellular organism. In the transition between the two types of behaviours, spatial wave patterns are produced, which is generally believed to be regulated by a certain biochemical clock that controls the direction of myxobacteria's motion. We provide a detailed analysis of such a clock and show that, for the proposed model, there exists some interval in parameter space where the behaviour is robust, i.e. the system behaves similarly for all parameter values. In more mathematical terms, we show the existence and convergence of trajectories to a limit cycle, and provide estimates of the parameter under which such a behaviour occurs. In addition, we show that the reported convergence result is robust, in the sense that any small change in the parameters leads to the same qualitative behaviour of the solution.

Exploring the Underlying Mechanisms of the Xenopus laevis Embryonic Cell Cycle

The cell cycle is an indispensable process in proliferation and development. Despite significant efforts, global quantification and physical understanding are still challenging. In this study, we explored the mechanisms of the Xenopus laevis embryonic cell cycle by quantifying the underlying landscape and flux. We uncovered the Mexican hat landscape of the Xenopus laevis embryonic cell cycle with several local basins and barriers on the oscillation path. The local basins characterize the different phases of the Xenopus laevis embryonic cell cycle, and the local barriers represent the checkpoints. The checkpoint mechanism of the cell cycle is revealed by the landscape basins and barriers. While landscape shape determines the stabilities of the states on the oscillation path, the curl flux force determines the stability of the cell cycle flow. Replication is fundamental for biology of living cells. We quantify the input energy (through the entropy production) as the thermodynamic requirement for initiation and sustainability of single cell life (cell cycle). Furthermore, we also quantify curl flux originated from the input energy as the dynamical requirement for the emergence of a new stable phase (cell cycle). This can provide a new quantitative insight for the origin of single cell life. In fact, the curl flux originated from the energy input or nutrition supply determines the speed and guarantees the progression of the cell cycle. The speed of the cell cycle is a hallmark of cancer. We characterized the quality of the cell cycle by the coherence time and found it is supported by the flux and energy cost. We are also able to quantify the degree of time irreversibility by the cross correlation function forward and backward in time from the stochastic traces in the simulation or experiments, providing a way for the quantification of the time irreversibility and the flux. Through global sensitivity analysis upon landscape and flux, we can identify the key elements for controlling the cell cycle speed. This can help to design an effective strategy for drug discovery against cancer.

The Importance of Kinase-Phosphatase Integration: Lessons from Mitosis

Kinases and phosphatases work antagonistically to control the behaviour of individual substrate molecules. This can be incorrectly extrapolated to imply that they also work antagonistically on the signals or processes that these molecules control. In fact, in many situations kinases and phosphatases work together to positively drive signal responses. We explain how this 'cooperativity' is critical for setting the amplitude, localisation, timing, and shape of phosphorylation signals. We use mitosis to illustrate why these properties are important for controlling mitotic entry, sister chromatid cohesion, kinetochore-microtubule attachments, the spindle assembly checkpoint, mitotic spindle elongation, and mitotic exit. These examples provide a rationale to explain how complex signalling behaviour could rely on similar types of integration within many other biological processes.

Biophysical clocks face a trade-off between internal and external noise resistance

Many organisms use free running circadian clocks to anticipate the day night cycle. However, others organisms use simple stimulus-response strategies ('hourglass clocks') and it is not clear when such strategies are sufficient or even preferable to free running clocks. Here, we find that free running clocks, such as those found in the cyanobacterium Synechococcus elongatus and humans, can efficiently project out light intensity fluctuations due to weather patterns ('external noise') by exploiting their limit cycle attractor. However, such limit cycles are necessarily vulnerable to 'internal noise'. Hence, at sufficiently high internal noise, point attractor-based 'hourglass' clocks, such as those found in a smaller cyanobacterium with low protein copy number, Prochlorococcus marinus, can outperform free running clocks. By interpolating between these two regimes in a diverse range of oscillators drawn from across biology, we demonstrate biochemical clock architectures that are best suited to different relative strengths of external and internal noise.

Notions of similarity for systems biology models

Systems biology models are rapidly increasing in complexity, size and numbers. When building large models, researchers rely on software tools for the retrieval, comparison, combination and merging of models, as well as for version control. These tools need to be able to quantify the differences and similarities between computational models. However, depending on the specific application, the notion of 'similarity' may greatly vary. A general notion of model similarity, applicable to various types of models, is still missing. Here we survey existing methods for the comparison of models, introduce quantitative measures for model similarity, and discuss potential applications of combined similarity measures. To frame model comparison as a general problem, we describe a theoretical approach to defining and computing similarities based on a combination of different model aspects. The six aspects that we define as potentially relevant for similarity are underlying encoding, references to biological entities, quantitative behaviour, qualitative behaviour, mathematical equations and parameters and network structure. We argue that future similarity measures will benefit from combining these model aspects in flexible, problem-specific ways to mimic users' intuition about model similarity, and to support complex model searches in databases.

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.

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.

PharmML in Action: an Interoperable Language for Modeling and Simulation

PharmML is an XML-based exchange format created with a focus on nonlinear mixed-effect (NLME) models used in pharmacometrics, but providing a very general framework that also allows describing mathematical and statistical models such as single-subject or nonlinear and multivariate regression models. This tutorial provides an overview of the structure of this language, brief suggestions on how to work with it, and use cases demonstrating its power and flexibility.

The Design Space of the Embryonic Cell Cycle Oscillator

One of the main tasks in the analysis of models of biomolecular networks is to characterize the domain of the parameter space that corresponds to a specific behavior. Given the large number of parameters in most models, this is no trivial task. We use a model of the embryonic cell cycle to illustrate the approaches that can be used to characterize the domain of parameter space corresponding to limit cycle oscillations, a regime that coordinates periodic entry into and exit from mitosis. Our approach relies on geometric construction of bifurcation sets, numerical continuation, and random sampling of parameters. We delineate the multidimensional oscillatory domain and use it to quantify the robustness of periodic trajectories. Although some of our techniques explore the specific features of the chosen system, the general approach can be extended to other models of the cell cycle engine and other biomolecular networks.

Modeling and simulation of biological systems using SPICE language

The article deals with BB-SPICE (SPICE for Biochemical and Biological Systems), an extension of the famous Simulation Program with Integrated Circuit Emphasis (SPICE). BB-SPICE environment is composed of three modules: a new textual and compact description formalism for biological systems, a converter that handles this description and generates the SPICE netlist of the equivalent electronic circuit and NGSPICE which is an open-source SPICE simulator. In addition, the environment provides back and forth interfaces with SBML (System Biology Markup Language), a very common description language used in systems biology. BB-SPICE has been developed in order to bridge the gap between the simulation of biological systems on the one hand and electronics circuits on the other hand. Thus, it is suitable for applications at the interface between both domains, such as development of design tools for synthetic biology and for the virtual prototyping of biosensors and lab-on-chip. Simulation results obtained with BB-SPICE and COPASI (an open-source software used for the simulation of biochemical systems) have been compared on a benchmark of models commonly used in systems biology. Results are in accordance from a quantitative viewpoint but BB-SPICE outclasses COPASI by 1 to 3 orders of magnitude regarding the computation time. Moreover, as our software is based on NGSPICE, it could take profit of incoming updates such as the GPU implementation, of the coupling with powerful analysis and verification tools or of the integration in design automation tools (synthetic biology).

Ca2+-dependent hyperpolarization hypothesis for mammalian sleep

The detailed molecular mechanisms underlying the regulation of sleep/wake cycles in mammals are elusive. In this regulation, at least two mechanisms with fast and slow time scales are involved. In the faster time scale, a state of non-rapid-eye-movement (NREM) sleep can be microscopically characterized by the millisecond-to-second-order electrical behavior of neurons, namely slow-wave oscillations described by electrophysiology. In the slower time scale, the total duration of NREM sleep is homeostatically regulated by sleep pressure (the need for sleep), which is usually sustained for hours or even days and can be macroscopically described by electroencephalogram (EEG). The longer dynamics of sleep regulation are often explained by the accumulation of sleep-inducing substances (SISs). However, we still do not have a concrete model to connect fast, microscopic dynamics and slow, macroscopic dynamics. In this review, we introduce a recent Ca²⁺-dependent hyperpolarization hypothesis, in which the Ca²⁺-dependent hyperpolarization of cortical-membrane potential induces slow-wave oscillation. Slow dynamics of the Ca²⁺-dependent hyperpolarization pathway might be regulated by recently identified sleep-promoting kinases as well as classical SISs. Therefore, cortical Ca²⁺-dependent hyperpolarization may be a fundamental mechanism connecting fast neural activity to the slow dynamics of sleep pressure.