A quantitative model of the effect of unreplicated DNA on cell cycle progression in frog egg extracts.

A critical goal in cell biology is to develop a systems-level perspective of eukaryotic cell cycle controls. Among these controls, a complex signaling network (called 'checkpoints') arrests progression through the cell cycle when there is a threat to genomic integrity such as unreplicated or damaged DNA. Understanding the regulatory principles of cell cycle checkpoints is important because loss of checkpoint regulation may be a requisite step on the roadway to cancer. Mathematical modeling has proved to be a useful guide to cell cycle regulation by revealing the importance of bistability, hysteresis and time lags in governing cell cycle transitions and checkpoint mechanisms. In this report, we propose a mathematical model of the frog egg cell cycle including effects of unreplicated DNA on progression into mitosis. By a stepwise approach utilizing parameter estimation tools, we build a model that is grounded in fundamental behaviors of the cell cycle engine (hysteresis and time lags), includes new elements in the signaling network (Myt1 and Chk1 kinases), and fits a large and diverse body of data from the experimental literature. The model provides a validated framework upon which to build additional aspects of the cell cycle checkpoint signaling network, including those control signals in the mammalian cell cycle that are commonly mutated in cancer.

Modeling molecular regulatory networks with JigCell and PET.

We demonstrate how to model macromolecular regulatory networks with JigCell and the Parameter Estimation Toolkit (PET). These software tools are designed specifically to support the process typically used by systems biologists to model complex regulatory circuits. A detailed example illustrates how a model of the cell cycle in frog eggs is created and then refined through comparison of simulation output with experimental data. We show how parameter estimation tools automatically generate rate constants that fit a model to experimental data.

Parameter estimation for a mathematical model of the cell cycle in frog eggs.

Parameter values for a kinetic model of the nuclear replication-division cycle in frog eggs are estimated by fitting solutions of the kinetic equations (nonlinear ordinary differential equations) to a suite of experimental observations. A set of optimal parameter values is found by minimizing an objective function defined as the orthogonal distance between the data and the model. The differential equations are solved by LSODAR and the objective function is minimized by ODRPACK. The optimal parameter values are close to the "guesstimates" of the modelers who first studied this problem. These tools are sufficiently general to attack more complicated problems, where guesstimation is impractical or unreliable.

Modeling regulatory networks at Virginia Tech.

The life of a cell is governed by the physicochemical properties of a complex network of interacting macromolecules (primarily genes and proteins). Hence, a full scientific understanding of and rational engineering approach to cell physiology require accurate mathematical models of the spatial and temporal dynamics of these macromolecular assemblies, especially the networks involved in integrating signals and regulating cellular responses. The Virginia Tech Consortium is involved in three specific goals of DARPA's computational biology program (Bio-COMP): to create effective software tools for modeling gene-protein-metabolite networks, to employ these tools in creating a new generation of realistic models, and to test and refine these models by well-conceived experimental studies. The special emphasis of this group is to understand the mechanisms of cell cycle control in eukaryotes (yeast cells and frog eggs). The software tools developed at Virginia Tech are designed to meet general requirements of modeling regulatory networks and are collected in a problem-solving environment called JigCell.