Circadian rhythms in a wide range of organisms are mediated by molecular mechanisms based on transcription-translation feedback. In this paper, we use bifurcation theory to explore mathematical models of genetic oscillators, based on Kim & Forger's interpretation of the circadian clock in mammals. At the core of their models is a negative feedback loop whereby PER proteins (PER1 and PER2) bind to and inhibit their transcriptional activator, BMAL1. For oscillations to occur, the dissociation constant of the PER:BMAL1 complex, [Formula: see text], must be ≤ 0.04 nM, which is orders of magnitude smaller than a reasonable expectation of 1-10 nM for this protein complex. We relax this constraint by two modifications to Kim & Forger's 'single negative feedback' (SNF) model: first, by introducing a multistep reaction chain for posttranscriptional modifications of Per mRNA and posttranslational phosphorylations of PER, and second, by replacing the first-order rate law for degradation of PER in the nucleus by a Michaelis-Menten rate law. These modifications increase the maximum allowable [Formula: see text] to ~2 nM. In a third modification, we consider an alternative rate law for gene transcription to resolve an unrealistically large rate of Per2 transcription at very low concentrations of BMAL1. Additionally, we studied extensions of the SNF model to include a second negative feedback loop (involving REV-ERB) and a supplementary positive feedback loop (involving ROR). Contrary to Kim & Forger's observations of these extended models, we find that, with our modifications, the supplementary positive feedback loop makes the oscillations more robust than observed in the models with one or two negative feedback loops. However, all three models are similarly robust when accounting for circadian rhythms (~24 h period) with [Formula: see text] ≥ 1 nM. Our results provide testable predictions for future experimental studies.

Abstract Background Second messengers, c-di-GMP and (p)ppGpp, are vital regulatory molecules in bacteria, influencing cellular processes such as biofilm formation, transcription, virulence, quorum sensing, and proliferation. While c-di-GMP and (p)ppGpp are both synthesized from GTP molecules, they play antagonistic roles in regulating the cell cycle. In C. crescentus, c-di-GMP works as a major regulator of pole morphogenesis and cell development. It inhibits cell motility and promotes S-phase entry by inhibiting the activity of the master regulator, CtrA. Intracellular (p)ppGpp accumulates under starvation, which helps bacteria to survive under stressful conditions through regulating nucleotide levels and halting proliferation. (p)ppGpp responds to nitrogen levels through RelA-SpoT homolog enzymes, detecting glutamine concentration using a nitrogen phosphotransferase system (PTS Ntr). This work relates the guanine nucleotide-based second messenger regulatory network with the bacterial PTS Ntr system and investigates how bacteria respond to nutrient availability. Results We propose a mathematical model for the dynamics of c-di-GMP and (p)ppGpp in C. crescentus and analyze how the guanine nucleotide-based second messenger system responds to certain environmental changes communicated through the PTS Ntr system. Our mathematical model consists of seven ODEs describing the dynamics of nucleotides and PTS Ntr enzymes. Our simulations are consistent with experimental observations and suggest, among other predictions, that SpoT can effectively decrease c-di-GMP levels in response to nitrogen starvation just as well as it increases (p)ppGpp levels. Thus, the activity of SpoT (or its homologues in other bacterial species) can likely influence the cell cycle by influencing both c-di-GMP and (p)ppGpp. Conclusions In this work, we integrate current knowledge and experimental observations from the literature to formulate a novel mathematical model. We analyze the model and demonstrate how the PTS Ntr system influences (p)ppGpp, c-di-GMP, GMP and GTP concentrations. While this model does not consider all aspects of PTS Ntr signaling, such as cross-talk with the carbon PTS system, here we present our first effort to develop a model of nutrient signaling in C. crescentus.

For over a century, the Michaelis-Menten (MM) rate law has been used to describe the rates of enzyme-catalyzed reactions and gene expression. Despite the ubiquity of the MM rate law, it accurately captures the dynamics of underlying biochemical reactions only so long as it is applied under the right condition, namely, that the substrate is in large excess over the enzyme-substrate complex. Unfortunately, in circumstances where its validity condition is not satisfied, especially so in protein interaction networks, the MM rate law has frequently been misused. In this review, we illustrate how inappropriate use of the MM rate law distorts the dynamics of the system, provides mistaken estimates of parameter values, and makes false predictions of dynamical features such as ultrasensitivity, bistability, and oscillations. We describe how these problems can be resolved with a slightly modified form of the MM rate law, based on the total quasi-steady state approximation (tQSSA). Furthermore, we show that the tQSSA can be used for accurate stochastic simulations at a lower computational cost than using the full set of mass-action rate laws. This review describes how to use quasi-steady state approximations in the right context, to prevent drawing erroneous conclusions from in silico simulations.

BACKGROUND: Cell size is a key characteristic that significantly affects many aspects of cellular physiology. There are specific control mechanisms during cell cycle that maintain the cell size within a range from generation to generation. Such control mechanisms introduce substantial variabilities to important properties of the cell cycle such as growth and division. To quantitatively study the effect of such variability in progression through cell cycle, detailed stochastic models are required. RESULTS: In this paper, a new hybrid stochastic model is proposed to study the effect of molecular noise and size control mechanism on the variabilities in cell cycle of the budding yeast Saccharomyces cerevisiae. The proposed model provides an accurate, yet computationally efficient approach for simulation of an intricate system by integrating the deterministic and stochastic simulation schemes. The developed hybrid stochastic model can successfully capture several key features of the cell cycle observed in experimental data. In particular, the proposed model: 1) confirms that the majority of noise in size control stems from low copy numbers of transcripts in the G1 phase, 2) identifies the size and time regulation modules in the size control mechanism, and 3) conforms with phenotypes of early G1 mutants in exquisite detail. CONCLUSIONS: Hybrid stochastic modeling approach can be used to provide quantitative descriptions for stochastic properties of the cell cycle within a computationally efficient framework.

Parameter estimation in discrete or continuous deterministic cell cycle models is challenging for several reasons, including the nature of what can be observed, and the accuracy and quantity of those observations. The challenge is even greater for stochastic models, where the number of simulations and amount of empirical data must be even larger to obtain statistically valid parameter estimates. The two main contributions of this work are (1) stochastic model parameter estimation based on directly matching multivariate probability distributions, and (2) a new quasi-Newton algorithm class QNSTOP for stochastic optimization problems. QNSTOP directly uses the random objective function value samples rather than creating ensemble statistics. QNSTOP is used here to directly match empirical and simulated joint probability distributions rather than matching summary statistics. Results are given for a current state-of-the-art stochastic cell cycle model of budding yeast, whose predictions match well some summary statistics and one-dimensional distributions from empirical data, but do not match well the empirical joint distributions. The nature of the mismatch provides insight into the weakness in the stochastic model.

Motivation: Mathematical models of cellular processes can systematically predict the phenotypes of novel combinations of multi-gene mutations. Searching for informative predictions and prioritizing them for experimental validation is challenging since the number of possible combinations grows exponentially in the number of mutations. Moreover, keeping track of the crosses needed to make new mutants and planning sequences of experiments is unmanageable when the experimenter is deluged by hundreds of potentially informative predictions to test. Results: We present CrossPlan, a novel methodology for systematically planning genetic crosses to make a set of target mutants from a set of source mutants. We base our approach on a generic experimental workflow used in performing genetic crosses in budding yeast. We prove that the CrossPlan problem is NP-complete. We develop an integer-linear-program (ILP) to maximize the number of target mutants that we can make under certain experimental constraints. We apply our method to a comprehensive mathematical model of the protein regulatory network controlling cell division in budding yeast. We also extend our solution to incorporate other experimental conditions such as a delay factor that decides the availability of a mutant and genetic markers to confirm gene deletions. The experimental flow that underlies our work is quite generic and our ILP-based algorithm is easy to modify. Hence, our framework should be relevant in plant and animal systems as well. Availability and implementation: CrossPlan code is freely available under GNU General Public Licence v3.0 at https://github.com/Murali-group/crossplan. Supplementary information: Supplementary data are available at Bioinformatics online.

In recent years, it has become increasingly apparent that antisense transcription plays an important role in the regulation of gene expression. The circadian clock is no exception: an antisense transcript of the mammalian core-clock gene PERIOD2 (PER2), which we shall refer to as Per2AS RNA, oscillates with a circadian period and a nearly 12 h phase shift from the peak expression of Per2 mRNA. In this paper, we ask whether Per2AS plays a regulatory role in the mammalian circadian clock by studying in silico the potential effects of interactions between Per2 and Per2AS RNAs on circadian rhythms. Based on the antiphasic expression pattern, we consider two hypotheses about how Per2 and Per2AS mutually interfere with each other's expression. In our pre-transcriptional model, the transcription of Per2AS RNA from the non-coding strand represses the transcription of Per2 mRNA from the coding strand and vice versa. In our post-transcriptional model, Per2 and Per2AS transcripts form a double-stranded RNA duplex, which is rapidly degraded. To study these two possible mechanisms, we have added terms describing our alternative hypotheses to a published mathematical model of the molecular regulatory network of the mammalian circadian clock. Our pre-transcriptional model predicts that transcriptional interference between Per2 and Per2AS can generate alternative modes of circadian oscillations, which we characterize in terms of the amplitude and phase of oscillation of core clock genes. In our post-transcriptional model, Per2/Per2AS duplex formation dampens the circadian rhythm. In a model that combines pre- and post-transcriptional controls, the period, amplitude and phase of circadian proteins exhibit non-monotonic dependencies on the rate of expression of Per2AS. All three models provide potential explanations of the observed antiphasic, circadian oscillations of Per2 and Per2AS RNAs. They make discordant predictions that can be tested experimentally in order to distinguish among these alternative hypotheses.

The size of a cell sets the scale for all biochemical processes within it, thereby affecting cellular fitness and survival. Hence, cell size needs to be kept within certain limits and relatively constant over multiple generations. However, how cells measure their size and use this information to regulate growth and division remains controversial. Here, we present two mechanistic mathematical models of the budding yeast (S. cerevisiae) cell cycle to investigate competing hypotheses on size control: inhibitor dilution and titration of nuclear sites. Our results suggest that an inhibitor-dilution mechanism, in which cell growth dilutes the transcriptional inhibitor Whi5 against the constant activator Cln3, can facilitate size homeostasis. This is achieved by utilising a positive feedback loop to establish a fixed size threshold for the Start transition, which efficiently couples cell growth to cell cycle progression. Yet, we show that inhibitor dilution cannot reproduce the size of mutants that alter the cell's overall ploidy and WHI5 gene copy number. By contrast, size control through titration of Cln3 against a constant number of genomic binding sites for the transcription factor SBF recapitulates both size homeostasis and the size of these mutant strains. Moreover, this model produces an imperfect 'sizer' behaviour in G1 and a 'timer' in S/G2/M, which combine to yield an 'adder' over the whole cell cycle; an observation recently made in experiments. Hence, our model connects these phenomenological data with the molecular details of the cell cycle, providing a systems-level perspective of budding yeast size control.

SUMMARY: Networks have become ubiquitous in systems biology. Visualization is a crucial component in their analysis. However, collaborations within research teams in network biology are hampered by software systems that are either specific to a computational algorithm, create visualizations that are not biologically meaningful, or have limited features for sharing networks and visualizations. We present GraphSpace, a web-based platform that fosters team science by allowing collaborating research groups to easily store, interact with, layout and share networks. AVAILABILITY AND IMPLEMENTATION: Anyone can upload and share networks at http://graphspace.org. In addition, the GraphSpace code is available at http://github.com/Murali-group/graphspace if a user wants to run his or her own server. CONTACT: [email protected]. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The cell division cycle of eukaryotes is governed by a complex network of cyclin-dependent protein kinases (CDKs) and auxiliary proteins that govern CDK activities. The control system must function reliably in the context of molecular noise that is inevitable in tiny yeast cells, because mistakes in sequencing cell cycle events are detrimental or fatal to the cell or its progeny. To assess the effects of noise on cell cycle progression requires not only extensive, quantitative, experimental measurements of cellular heterogeneity but also comprehensive, accurate, mathematical models of stochastic fluctuations in the CDK control system. In this paper we provide a stochastic model of the budding yeast cell cycle that accurately accounts for the variable phenotypes of wild-type cells and more than 20 mutant yeast strains simulated in different growth conditions. We specifically tested the role of feedback regulations mediated by G1- and SG2M-phase cyclins to minimize the noise in cell cycle progression. Details of the model are informed and tested by quantitative measurements (by fluorescence in situ hybridization) of the joint distributions of mRNA populations in yeast cells. We use the model to predict the phenotypes of ~30 mutant yeast strains that have not yet been characterized experimentally.

In present-day eukaryotes, the cell division cycle is controlled by a complex network of interacting proteins, including members of the cyclin and cyclin-dependent protein kinase (Cdk) families, and the Anaphase Promoting Complex (APC). Successful progression through the cell cycle depends on precise, temporally ordered regulation of the functions of these proteins. In light of this complexity, it is surprising that in fission yeast, a minimal Cdk network consisting of a single cyclin-Cdk fusion protein can control DNA synthesis and mitosis in a manner that is indistinguishable from wild type. To improve our understanding of the cell cycle regulatory network, we built and analysed a mathematical model of the molecular interactions controlling the G1/S and G2/M transitions in these minimal cells. The model accounts for all observed properties of yeast strains operating with the fusion protein. Importantly, coupling the model's predictions with experimental analysis of alternative minimal cells, we uncover an explanation for the unexpected fact that elimination of inhibitory phosphorylation of Cdk is benign in these strains while it strongly affects normal cells. Furthermore, in the strain without inhibitory phosphorylation of the fusion protein, the distribution of cell size at division is unusually broad, an observation that is accounted for by stochastic simulations of the model. Our approach provides novel insights into the organization and quantitative regulation of wild type cell cycle progression. In particular, it leads us to propose a new mechanistic model for the phenomenon of mitotic catastrophe, relying on a combination of unregulated, multi-cyclin-dependent Cdk activities.

Cell-fate asymmetry in the predivisional cell of Caulobacter crescentus requires that the regulatory protein DivL localizes to the new pole of the cell where it up-regulates CckA kinase, resulting in a gradient of CtrA~P across the cell. In the preceding stage of the cell cycle (the "stalked" cell), DivL is localized uniformly along the cell membrane and maintained in an inactive form by DivK~P. It is unclear how DivL overcomes inhibition by DivK~P in the predivisional cell simply by changing its location to the new pole. It has been suggested that co-localization of DivL with PleC phosphatase at the new pole is essential to DivL's activity there. However, there are contrasting views on whether the bifunctional enzyme, PleC, acts as a kinase or phosphatase at the new pole. To explore these ambiguities, we formulated a mathematical model of the spatiotemporal distributions of DivL, PleC and associated proteins (DivJ, DivK, CckA, and CtrA) during the asymmetric division cycle of a Caulobacter cell. By varying localization profiles of DivL and PleC in our model, we show how the physiologically observed spatial distributions of these proteins are essential for the transition from a stalked cell to a predivisional cell. Our simulations suggest that PleC is a kinase in predivisional cells, and that, by sequestering DivK~P, the kinase form of PleC enables DivL to be reactivated at the new pole. Hence, co-localization of PleC kinase and DivL is essential to establishing cellular asymmetry. Our simulations reproduce the experimentally observed spatial distribution and phosphorylation status of CtrA in wild-type and mutant cells. Based on the model, we explore novel combinations of mutant alleles, making predictions that can be tested experimentally.

The free-living aquatic bacterium, Caulobacter crescentus, exhibits two different morphologies during its life cycle. The morphological change from swarmer cell to stalked cell is a result of changes of function of two bi-functional histidine kinases, PleC and CckA. Here, we describe a detailed molecular mechanism by which the function of PleC changes between phosphatase and kinase state. By mathematical modeling of our proposed molecular interactions, we derive conditions under which PleC, CckA and its response regulators exhibit bistable behavior, thus providing a scenario for robust switching between swarmer and stalked states. Our simulations are in reasonable agreement with in vitro and in vivo experimental observations of wild type and mutant phenotypes. According to our model, the kinase form of PleC is essential for the swarmer-to-stalked transition and to prevent premature development of the swarmer pole. Based on our results, we reconcile some published experimental observations and suggest novel mutants to test our predictions.

The innate immune system, acting as the first line of host defense, senses and adapts to foreign challenges through complex intracellular and intercellular signaling networks. Endotoxin tolerance and priming elicited by macrophages are classic examples of the complex adaptation of innate immune cells. Upon repetitive exposures to different doses of bacterial endotoxin (lipopolysaccharide) or other stimulants, macrophages show either suppressed or augmented inflammatory responses compared to a single exposure to the stimulant. Endotoxin tolerance and priming are critically involved in both immune homeostasis and the pathogenesis of diverse inflammatory diseases. However, the underlying molecular mechanisms are not well understood. By means of a computational search through the parameter space of a coarse-grained three-node network with a two-stage Metropolis sampling approach, we enumerated all the network topologies that can generate priming or tolerance. We discovered three major mechanisms for priming (pathway synergy, suppressor deactivation, activator induction) and one for tolerance (inhibitor persistence). These results not only explain existing experimental observations, but also reveal intriguing test scenarios for future experimental studies to clarify mechanisms of endotoxin priming and tolerance.

The reciprocal differentiation of T helper 17 (T(H)17) cells and induced regulatory T (iT(reg)) cells plays a critical role in both the pathogenesis and resolution of diverse human inflammatory diseases. Although initial studies suggested a stable commitment to either the T(H)17 or the iT(reg) lineage, recent results reveal remarkable plasticity and heterogeneity, reflected in the capacity of differentiated effectors cells to be reprogrammed among T(H)17 and iT(reg) lineages and the intriguing phenomenon that a group of naïve precursor CD4(+) T cells can be programmed into phenotypically diverse populations by the same differentiation signal, transforming growth factor beta. To reconcile these observations, we have built a mathematical model of T(H)17/iT(reg) differentiation that exhibits four different stable steady states, governed by pitchfork bifurcations with certain degrees of broken symmetry. According to the model, a group of precursor cells with some small cell-to-cell variability can differentiate into phenotypically distinct subsets of cells, which exhibit distinct levels of the master transcription-factor regulators for the two T cell lineages. A dynamical control system with these properties is flexible enough to be steered down alternative pathways by polarizing signals, such as interleukin-6 and retinoic acid and it may be used by the immune system to generate functionally distinct effector cells in desired fractions in response to a range of differentiation signals. Additionally, the model suggests a quantitative explanation for the phenotype with high expression levels of both master regulators. This phenotype corresponds to a re-stabilized co-expressing state, appearing at a late stage of differentiation, rather than a bipotent precursor state observed under some other circumstances. Our simulations reconcile most published experimental observations and predict novel differentiation states as well as transitions among different phenotypes that have not yet been observed experimentally.

The timing of DNA synthesis, mitosis and cell division is regulated by a complex network of biochemical reactions that control the activities of a family of cyclin-dependent kinases. The temporal dynamics of this reaction network is typically modeled by nonlinear differential equations describing the rates of the component reactions. This approach provides exquisite details about molecular regulatory processes but is hampered by the need to estimate realistic values for the many kinetic constants that determine the reaction rates. It is difficult to estimate these kinetic constants from available experimental data. To avoid this problem, modelers often resort to 'qualitative' modeling strategies, such as Boolean switching networks, but these models describe only the coarsest features of cell cycle regulation. In this paper we describe a hybrid approach that combines the best features of continuous differential equations and discrete Boolean networks. Cyclin abundances are tracked by piecewise linear differential equations for cyclin synthesis and degradation. Cyclin synthesis is regulated by transcription factors whose activities are represented by discrete variables (0 or 1) and likewise for the activities of the ubiquitin-ligating enzyme complexes that govern cyclin degradation. The discrete variables change according to a predetermined sequence, with the times between transitions determined in part by cyclin accumulation and degradation and as well by exponentially distributed random variables. The model is evaluated in terms of flow cytometry measurements of cyclin proteins in asynchronous populations of human cell lines. The few kinetic constants in the model are easily estimated from the experimental data. Using this hybrid approach, modelers can quickly create quantitatively accurate, computational models of protein regulatory networks in cells.

Models of regulatory networks become more difficult to construct and understand as they grow in size and complexity. Large models are usually built up from smaller models, representing subsets of reactions within the larger network. To assist modelers in this composition process, we present a formal approach for model composition, a wizard-style program for implementing the approach, and suggested language extensions to the Systems Biology Markup Language to support model composition. To illustrate the features of our approach and how to use the JigCell Composition Wizard, we build up a model of the eukaryotic cell cycle "engine" from smaller pieces.

MOTIVATION: Models of regulatory networks become more difficult to construct and understand as they grow in size and complexity. Modelers naturally build large models from smaller components that each represent subsets of reactions within the larger network. To assist modelers in this process, we present model aggregation, which defines models in terms of components that are designed for the purpose of being combined. RESULTS: We have implemented a model editor that incorporates model aggregation, and we suggest supporting extensions to the Systems Biology Markup Language (SBML) Level 3. We illustrate aggregation with a model of the eukaryotic cell cycle 'engine' created from smaller pieces. AVAILABILITY: Java implementations are available in the JigCell Aggregation Connector. See http://jigcell.biol.vt.edu. CONTACT: [email protected]

The asymmetric cell division cycle of Caulobacter crescentus is orchestrated by an elaborate gene-protein regulatory network, centered on three major control proteins, DnaA, GcrA and CtrA. The regulatory network is cast into a quantitative computational model to investigate in a systematic fashion how these three proteins control the relevant genetic, biochemical and physiological properties of proliferating bacteria. Different controls for both swarmer and stalked cell cycles are represented in the mathematical scheme. The model is validated against observed phenotypes of wild-type cells and relevant mutants, and it predicts the phenotypes of novel mutants and of known mutants under novel experimental conditions. Because the cell cycle control proteins of Caulobacter are conserved across many species of alpha-proteobacteria, the model we are proposing here may be applicable to other genera of importance to agriculture and medicine (e.g., Rhizobium, Brucella).

Progression of a cell through the division cycle is tightly controlled at different steps to ensure the integrity of genome replication and partitioning to daughter cells. From published experimental evidence, we propose a molecular mechanism for control of the cell division cycle in Caulobacter crescentus. The mechanism, which is based on the synthesis and degradation of three "master regulator" proteins (CtrA, GcrA, and DnaA), is converted into a quantitative model, in order to study the temporal dynamics of these and other cell cycle proteins. The model accounts for important details of the physiology, biochemistry, and genetics of cell cycle control in stalked C. crescentus cell. It reproduces protein time courses in wild-type cells, mimics correctly the phenotypes of many mutant strains, and predicts the phenotypes of currently uncharacterized mutants. Since many of the proteins involved in regulating the cell cycle of C. crescentus are conserved among many genera of alpha-proteobacteria, the proposed mechanism may be applicable to other species of importance in agriculture and medicine.

In metabolic networks, metabolites are usually present in great excess over the enzymes that catalyze their interconversion, and describing the rates of these reactions by using the Michaelis-Menten rate law is perfectly valid. This rate law assumes that the concentration of enzyme-substrate complex (C) is much less than the free substrate concentration (S0). However, in protein interaction networks, the enzymes and substrates are all proteins in comparable concentrations, and neglecting C with respect to S0 is not valid. Borghans, DeBoer, and Segel developed an alternative description of enzyme kinetics that is valid when C is comparable to S0. We extend this description, which Borghans et al. call the total quasi-steady state approximation, to networks of coupled enzymatic reactions. First, we analyze an isolated Goldbeter-Koshland switch when enzymes and substrates are present in comparable concentrations. Then, on the basis of a real example of the molecular network governing cell cycle progression, we couple two and three Goldbeter-Koshland switches together to study the effects of feedback in networks of protein kinases and phosphatases. Our analysis shows that the total quasi-steady state approximation provides an excellent kinetic formalism for protein interaction networks, because (1) it unveils the modular structure of the enzymatic reactions, (2) it suggests a simple algorithm to formulate correct kinetic equations, and (3) contrary to classical Michaelis-Menten kinetics, it succeeds in faithfully reproducing the dynamics of the network both qualitatively and quantitatively.

Converting a biochemical reaction network to a set of kinetic rate equations is tedious and error prone. We describe known interface paradigms for inputing models of intracellular regulatory networks: graphical layout (diagrams), wizards, scripting languages, and direct entry of chemical equations. We present the JigCell Model Builder, which allows users to define models as a set of reaction equations using a spreadsheet (an example of direct entry of equations) and outputs model definitions in the Systems Biology Markup Language, Level 2. We present the results of two usability studies. The spreadsheet paradigm demonstrated its effectiveness in reducing the number of errors made by modelers when compared to hand conversion of a wiring diagram to differential equations. A comparison of representatives of the four interface paradigms for a simple model of the cell cycle was conducted which measured time, mouse clicks, and keystrokes to enter the model, and the number of screens needed to view the contents of the model. All four paradigms had similar data entry times. The spreadsheet and scripting language approaches require significantly fewer screens to view the models than do the wizard or graphical layout approaches.

Summary form only given. We present a sophisticated and efficient parallel scheme for the DIRECT global optimization algorithm of Jones et al. (1993). Although several sequential implementations for this algorithm have been successfully applied to large scale MDO problems, few parallel versions of the DIRECT algorithm have addressed well algorithm characteristics such as a single starting point, an unpredictable workload, and a strong data dependency. These challenges engender many interesting design issues including domain decomposition, data access and management, and workload balancing. A hierarchical parallel scheme has been developed to address these challenges at three levels. Each level is supported by parallel and distributed data structures to access shared data sets, distribute workload, or exchange messages. Parameter estimation problems in systems biology provide an ideal application context for the present work. Global nonlinear parameter estimation results obtained on a 200 node Linux cluster are given for a cell cycle model for frog eggs.

SUMMARY: We describe the JigCell Model Builder (JCMB), a tool for creating biochemical reaction network models. JCMB is designed for ease of use and its interface uses the standard spreadsheet metaphor. The JigCell Run Manager (JCRM) is a tool for organizing the large collections of simulation runs typically required by reaction network modeling activities. AVAILABILITY: JCMB and JCRM are part of the JigCell suite available at http://jigcell.biol.vt.edu.