Time scale and dimension analysis of a budding yeast cell cycle model

The progress through the eukaryotic cell division cycle is driven by an underlying molecular regulatory network. Cell cycle progression can be considered as a series of irreversible transitions from one steady state to another in the correct order. Although this view has been put forward some time ago, it has not been quantitatively proven yet. Bifurcation analysis of a model for the budding yeast cell cycle has identified only two different steady states (one for G1 and one for mitosis) using cell mass as a bifurcation parameter. By analyzing the same model, using different methods of dynamical systems theory, we provide evidence for transitions among several different steady states during the budding yeast cell cycle. By calculating the eigenvalues of the Jacobian of kinetic differential equations we have determined the stability of the cell cycle trajectories of the Chen model. Based on the sign of the real part of the eigenvalues, the cell cycle can be divided into excitation and relaxation periods. During an excitation period, the cell cycle control system leaves a formerly stable steady state and, accordingly, excitation periods can be associated with irreversible cell cycle transitions like START, entry into mitosis and exit from mitosis. During relaxation periods, the control system asymptotically approaches the new steady state. We also show that the dynamical dimension of the Chen’s model fluctuates by increasing during excitation periods followed by decrease during relaxation periods. In each relaxation period the dynamical dimension of the model drops to one, indicating a period where kinetic processes are in steady state and all concentration changes are driven by the increase of cytoplasmic growth.We apply two numerical methods, which have not been used to analyze biological control systems. These methods are more sensitive than the bifurcation analysis used before because they identify those transitions between steady states that are not controlled by a bifurcation parameter (e.g. cell mass). Therefore by applying these tools for a cell cycle control model, we provide a deeper understanding of the dynamical transitions in the underlying molecular network

USE OF GROWTH ANALYSIS TO EVALUATE GENETIC MECHANISMS AFFECTING ACHENE YIELD FORMATION OF SUNFLOWER

The main objective of this study was to investigate the process of dry matter accumulation (DMA) in achenes during the grain-filling period of fifty sunflower genotypes by using the functional method of growth analysis in a field trial at Bicsérd, Hungary. The Hunt-formula of lnY= P0 + P1*X + P2*X2 was fitted to data. Maximum yield (Ymax), the average of the absolute growth rate (AGRavg), maximum growth rate (AGRmax), date of the maximum growth rate (Xagrmax), and the average of the relative growth rate (RGRavg) were calculated from growth curves for hybrids and replications. Significant differences among hybrids and their interaction with sampling dates indicate hybrid differences in the intensity of DMA accumulation. The strongest correlation was observed between the parameters of Ymax and AGRmax

Model of For3p-Mediated Actin Cable Assembly in Fission Yeast

Formin For3p nucleates actin cables at the tips of fission yeast cells for polarized cell growth. The results of prior experiments have suggested a possible mechanism for actin cable assembly that involves association of For3p near cell tips, For3p-mediated actin polymerization, retrograde flow of actin cables toward the cell center, For3p dissociation from cell tips, and cable disassembly. We used analytical and computational modeling to test the validity and implications of the proposed coupled For3p/actin mechanism. We compared the model to prior experiments quantitatively and generated predictions for the expected behavior of the actin cable system upon changes of parameter values. We found that the model generates stable steady states with realistic values of rate constants and actin and For3p concentrations. Comparison of our results to previous experiments monitoring the FRAP of For3p-3GFP and the response of actin cables to treatments with actin depolymerizing drugs provided further support for the model. We identified the set of parameter values that produces results in agreement with experimental observations. We discuss future experiments that will help test the model's predictions and eliminate other possible mechanisms. The results of the model suggest that flow of actin cables may establish actin and For3p concentration gradients in the cytoplasm that could be important in global cell patterning

Transcriptional Regulation Is a Major Controller of Cell Cycle Transition Dynamics

DNA replication, mitosis and mitotic exit are critical transitions of the cell cycle which normally occur only once per cycle. A universal control mechanism was proposed for the regulation of mitotic entry in which Cdk helps its own activation through two positive feedback loops. Recent discoveries in various organisms showed the importance of positive feedbacks in other transitions as well. Here we investigate if a universal control system with transcriptional regulation(s) and post-translational positive feedback(s) can be proposed for the regulation of all cell cycle transitions. Through computational modeling, we analyze the transition dynamics in all possible combinations of transcriptional and post-translational regulations. We find that some combinations lead to ‘sloppy’ transitions, while others give very precise control. The periodic transcriptional regulation through the activator or the inhibitor leads to radically different dynamics. Experimental evidence shows that in cell cycle transitions of organisms investigated for cell cycle dependent periodic transcription, only the inhibitor OR the activator is under cyclic control and never both of them. Based on these observations, we propose two transcriptional control modes of cell cycle regulation that either STOP or let the cycle GO in case of a transcriptional failure. We discuss the biological relevance of such differences

Quantifying the Dynamics of Coupled Networks of Switches and Oscillators

Complex network dynamics have been analyzed with models of systems of coupled switches or systems of coupled oscillators. However, many complex systems are composed of components with diverse dynamics whose interactions drive the system's evolution. We, therefore, introduce a new modeling framework that describes the dynamics of networks composed of both oscillators and switches. Both oscillator synchronization and switch stability are preserved in these heterogeneous, coupled networks. Furthermore, this model recapitulates the qualitative dynamics for the yeast cell cycle consistent with the hypothesized dynamics resulting from decomposition of the regulatory network into dynamic motifs. Introducing feedback into the cell-cycle network induces qualitative dynamics analogous to limitless replicative potential that is a hallmark of cancer. As a result, the proposed model of switch and oscillator coupling provides the ability to incorporate mechanisms that underlie the synchronized stimulus response ubiquitous in biochemical systems

An Extended Gene Protein/Products Boolean Network Model Including Post-Transcriptional Regulation

Background: Networks Biology allows the study of complex interactions between biological systems using formal, well structured, and computationally friendly models. Several different network models can be created, depending on the type of interactions that need to be investigated. Gene Regulatory Networks (GRN) are an effective model commonly used to study the complex regulatory mechanisms of a cell. Unfortunately, given their intrinsic complexity and non discrete nature, the computational study of realistic-sized complex GRNs requires some abstractions. Boolean Networks (BNs), for example, are a reliable model that can be used to represent networks where the possible state of a node is a boolean value (0 or 1). Despite this strong simplification, BNs have been used to study both structural and dynamic properties of real as well as randomly generated GRNs. Results: In this paper we show how it is possible to include the post-transcriptional regulation mechanism (a key process mediated by small non-coding RNA molecules like the miRNAs) into the BN model of a GRN. The enhanced BN model is implemented in a software toolkit (EBNT) that allows to analyze boolean GRNs from both a structural and a dynamic point of view. The open-source toolkit is compatible with available visualization tools like Cytoscape and allows to run detailed analysis of the network topology as well as of its attractors, trajectories, and state-space. In the paper, a small GRN built around the mTOR gene is used to demonstrate the main capabilities of the toolkit. Conclusions: The extended model proposed in this paper opens new opportunities in the study of gene regulation. Several of the successful researches done with the support of BN to understand high-level characteristics of regulatory networks, can now be improved to better understand the role of post-transcriptional regulation for example as a network-wide noise-reduction or stabilization mechanism

Syntax-Guided Optimal Synthesis for Chemical Reaction Networks

We study the problem of optimal syntax-guided synthesis of stochastic Chemical Reaction Networks (CRNs) that plays a fundamental role in design automation of molecular devices and in the construction of predictive biochemical models. We propose a sketching language for CRNs that concisely captures syntactic constraints on the network topology and allows its under-specification. Given a sketch, a correctness specification, and a cost function defined over the CRN syntax, our goal is to find a CRN that simultaneously meets the constraints, satisfies the specification and minimizes the cost function. To ensure computational feasibility of the synthesis process, we employ the Linear Noise Approximation allowing us to encode the synthesis problem as a satisfiability modulo theories problem over a set of parametric Ordinary Differential Equations (ODEs). We design and implement a novel algorithm for the optimal synthesis of CRNs that employs almost complete refutation procedure for SMT over reals and ODEs, and exploits a meta-sketching abstraction controlling the search strategy. Through relevant case studies we demonstrate that our approach significantly improves the capability of existing methods for synthesis of biochemical systems and paves the way towards their automated and provably-correct design

A Novel Method to Verify Multilevel Computational Models of Biological Systems Using Multiscale Spatio-Temporal Meta Model Checking

Insights gained from multilevel computational models of biological systems can be translated into real-life applications only if the model correctness has been verified first. One of the most frequently employed in silico techniques for computational model verification is model checking. Traditional model checking approaches only consider the evolution of numeric values, such as concentrations, over time and are appropriate for computational models of small scale systems (e.g. intracellular networks). However for gaining a systems level understanding of how biological organisms function it is essential to consider more complex large scale biological systems (e.g. organs). Verifying computational models of such systems requires capturing both how numeric values and properties of (emergent) spatial structures (e.g. area of multicellular population) change over time and across multiple levels of organization, which are not considered by existing model checking approaches. To address this limitation we have developed a novel approximate probabilistic multiscale spatio-temporal meta model checking methodology for verifying multilevel computational models relative to specifications describing the desired/expected system behaviour. The methodology is generic and supports computational models encoded using various high-level modelling formalisms because it is defined relative to time series data and not the models used to generate it. In addition, the methodology can be automatically adapted to case study specific types of spatial structures and properties using the spatio-temporal meta model checking concept. To automate the computational model verification process we have implemented the model checking approach in the software tool Mule (http://mule.modelchecking.org). Its applicability is illustrated against four systems biology computational models previously published in the literature encoding the rat cardiovascular system dynamics, the uterine contractions of labour, the Xenopus laevis cell cycle and the acute inflammation of the gut and lung. Our methodology and software will enable computational biologists to efficiently develop reliable multilevel computational models of biological systems