Integrated Model of Chemical Perturbations of a Biological Pathway Using 18 In Vitro High-Throughput Screening Assays for the Estrogen Receptor.

We demonstrate a computational network model that integrates 18 in vitro, high-throughput screening assays measuring estrogen receptor (ER) binding, dimerization, chromatin binding, transcriptional activation, and ER-dependent cell proliferation. The network model uses activity patterns across the in vitro assays to predict whether a chemical is an ER agonist or antagonist, or is otherwise influencing the assays through a manner dependent on the physics and chemistry of the technology platform ("assay interference"). The method is applied to a library of 1812 commercial and environmental chemicals, including 45 ER positive and negative reference chemicals. Among the reference chemicals, the network model correctly identified the agonists and antagonists with the exception of very weak compounds whose activity was outside the concentration range tested. The model agonist score also correlated with the expected potency class of the active reference chemicals. Of the 1812 chemicals evaluated, 111 (6.1%) were predicted to be strongly ER active in agonist or antagonist mode. This dataset and model were also used to begin a systematic investigation of assay interference. The most prominent cause of false-positive activity (activity in an assay that is likely not due to interaction of the chemical with ER) is cytotoxicity. The model provides the ability to prioritize a large set of important environmental chemicals with human exposure potential for additional in vivo endocrine testing. Finally, this model is generalizable to any molecular pathway for which there are multiple upstream and downstream assays available.

Computational modelling of T-cell formation kinetics: output regulated by initial proliferation-linked deferral of developmental competence.

Bone-marrow-derived progenitors must continually enter the thymus of an adult mouse to sustain T-cell homeostasis, yet only a few input cells per day are sufficient to support a yield of 5 × 10(7) immature T-cells per day and an eventual output of 1-2 × 10(6) mature cells per day. While substantial progress has been made to delineate the developmental pathway of T-cell lineage commitment, still little is known about the relationship between differentiation competence and the remarkable expansion of the earliest (DN1 stage) T-cell progenitors. To address this question, we developed computational models where the probability to progress to the next stage (DN2) is related to division number. To satisfy differentiation kinetics and overall cell yield data, our models require that adult DN1 cells divide multiple times before becoming competent to progress into DN2 stage. Our findings were subsequently tested by in vitro experiments, where putative early and later-stage DN1 progenitors from the thymus were purified and their progression into DN2 was measured. These experiments showed that the two DN1 sub-populations divided with similar rates, but progressed to the DN2 stage with different rates, thus providing experimental evidence that DN1 cells increase their commitment probability in a cell-intrinsic manner as they undergo cell division. Proliferation-linked shifts in eligibility of DN1 cells to undergo specification thus control kinetics of T-cell generation.

Computational multiscale modeling of embryo development.

Recent advances in live imaging and genetics of mammalian development which integrate observations of biochemical activity, cell-cell signaling and mechanical interactions between cells pave the way for predictive mathematical multi-scale modeling. In early mammalian embryo development, two of the most critical events which lead to tissue patterning involve changes in gene expression as well as mechanical interactions between cells. We discuss the relevance of mathematical modeling of multi-cellular systems and in particular in simulating these patterns and describe some of the technical challenges one encounters. Many of these issues are not unique for the embryonic system but are shared by other multi-cellular modeling areas.

Probing the role of stochasticity in a model of the embryonic stem cell: heterogeneous gene expression and reprogramming efficiency.

BACKGROUND: Embryonic stem cells (ESC) have the capacity to self-renew and remain pluripotent, while continuously providing a source of a variety of differentiated cell types. Understanding what governs these properties at the molecular level is crucial for stem cell biology and its application to regenerative medicine. Of particular relevance is to elucidate those molecular interactions which govern the reprogramming of somatic cells into ESC. A computational approach can be used as a framework to explore the dynamics of a simplified network of the ESC with the aim to understand how stem cells differentiate and also how they can be reprogrammed from somatic cells. RESULTS: We propose a computational model of the embryonic stem cell network, in which a core set of transcription factors (TFs) interact with each other and are induced by external factors. A stochastic treatment of the network dynamics suggests that NANOG heterogeneity is the deciding factor for the stem cell fate. In particular, our results show that the decision of staying in the ground state or commitment to a differentiated state is fundamentally stochastic, and can be modulated by the addition of external factors (2i/3i media), which have the effect of reducing fluctuations in NANOG expression. Our model also hosts reprogramming of a committed cell into an ESC by over-expressing OCT4. In this context, we recapitulate the important experimental result that reprogramming efficiency peaks when OCT4 is over-expressed within a specific range of values. CONCLUSIONS: We have demonstrated how a stochastic computational model based upon a simplified network of TFs in ESCs can elucidate several key observed dynamical features. It accounts for (i) the observed heterogeneity of key regulators, (ii) characterizes the ESC under certain external stimuli conditions and (iii) describes the occurrence of transitions from the ESC to the differentiated state. Furthermore, the model (iv) provides a framework for reprogramming from somatic cells and conveys an understanding of reprogramming efficiency as a function of OCT4 over-expression.

Cytokinin signaling as a positional cue for patterning the apical-basal axis of the growing Arabidopsis shoot meristem.

The transcription factor WUSCHEL (WUS) acts from a well-defined domain within the Arabidopsis thaliana shoot apical meristem (SAM) to maintain a stem cell niche. A negative-feedback loop involving the CLAVATA (CLV) signaling pathway regulates the number of WUS-expressing cells and provides the current paradigm for the homeostatic maintenance of stem cell numbers. Despite the continual turnover of cells in the SAM during development, the WUS domain remains patterned at a fixed distance below the shoot apex. Recent work has uncovered a positive-feedback loop between WUS function and the plant hormone cytokinin. Furthermore, loss of function of the cytokinin biosynthetic gene, LONELY GUY (LOG), results in a wus-like phenotype in rice. Herein, we find the Arabidopsis LOG4 gene is expressed in the SAM epidermis. We use this to develop a computational model representing a growing SAM to suggest the plausibility that apically derived cytokinin and CLV signaling, together, act as positional cues for patterning the WUS domain within the stem cell niche. Furthermore, model simulations backed by experimental data suggest a previously unknown negative feedback between WUS function and cytokinin biosynthesis in the Arabidopsis SAM epidermis. These results suggest a plausible dynamic feedback principle by which the SAM stem cell niche is patterned.

Simulating the mammalian blastocyst--molecular and mechanical interactions pattern the embryo.

Mammalian embryogenesis is a dynamic process involving gene expression and mechanical forces between proliferating cells. The exact nature of these interactions, which determine the lineage patterning of the trophectoderm and endoderm tissues occurring in a highly regulated manner at precise periods during the embryonic development, is an area of debate. We have developed a computational modeling framework for studying this process, by which the combined effects of mechanical and genetic interactions are analyzed within the context of proliferating cells. At a purely mechanical level, we demonstrate that the perpendicular alignment of the animal-vegetal (a-v) and embryonic-abembryonic (eb-ab) axes is a result of minimizing the total elastic conformational energy of the entire collection of cells, which are constrained by the zona pellucida. The coupling of gene expression with the mechanics of cell movement is important for formation of both the trophectoderm and the endoderm. In studying the formation of the trophectoderm, we contrast and compare quantitatively two hypotheses: (1) The position determines gene expression, and (2) the gene expression determines the position. Our model, which couples gene expression with mechanics, suggests that differential adhesion between different cell types is a critical determinant in the robust endoderm formation. In addition to differential adhesion, two different testable hypotheses emerge when considering endoderm formation: (1) A directional force acts on certain cells and moves them into forming the endoderm layer, which separates the blastocoel and the cells of the inner cell mass (ICM). In this case the blastocoel simply acts as a static boundary. (2) The blastocoel dynamically applies pressure upon the cells in contact with it, such that cell segregation in the presence of differential adhesion leads to the endoderm formation. To our knowledge, this is the first attempt to combine cell-based spatial mechanical simulations with genetic networks to explain mammalian embryogenesis. Such a framework provides the means to test hypotheses in a controlled in silico environment.

Computational morphodynamics of plants: integrating development over space and time.

The emerging field of computational morphodynamics aims to understand the changes that occur in space and time during development by combining three technical strategies: live imaging to observe development as it happens; image processing and analysis to extract quantitative information; and computational modelling to express and test time-dependent hypotheses. The strength of the field comes from the iterative and combined use of these techniques, which has provided important insights into plant development.

Variability in the control of cell division underlies sepal epidermal patterning in Arabidopsis thaliana.

How growth and proliferation are precisely controlled in organs during development and how the regulation of cell division contributes to the formation of complex cell type patterns are important questions in developmental biology. Such a pattern of diverse cell sizes is characteristic of the sepals, the outermost floral organs, of the plant Arabidopsis thaliana. To determine how the cell size pattern is formed in the sepal epidermis, we iterate between generating predictions from a computational model and testing these predictions through time-lapse imaging. We show that the cell size diversity is due to the variability in decisions of individual cells about when to divide and when to stop dividing and enter the specialized endoreduplication cell cycle. We further show that altering the activity of cell cycle inhibitors biases the timing and changes the cell size pattern as our model predicts. Models and observations together demonstrate that variability in the time of cell division is a major determinant in the formation of a characteristic pattern.

Computational morphodynamics: a modeling framework to understand plant growth.

Computational morphodynamics utilizes computer modeling to understand the development of living organisms over space and time. Results from biological experiments are used to construct accurate and predictive models of growth. These models are then used to make novel predictions that provide further insight into the processes involved, which can be tested experimentally to either confirm or rule out the validity of the computational models. This review highlights two fundamental challenges: (a) to understand the feedback between mechanics of growth and chemical or molecular signaling, and (b) to design models that span and integrate single cell behavior with tissue development. We review different approaches to model plant growth and discuss a variety of model types that can be implemented to demonstrate how the interplay between computational modeling and experimentation can be used to explore the morphodynamics of plant development.

Multiple feedback loops through cytokinin signaling control stem cell number within the Arabidopsis shoot meristem.

A central unanswered question in stem cell biology, both in plants and in animals, is how the spatial organization of stem cell niches are maintained as cells move through them. We address this question for the shoot apical meristem (SAM) which harbors pluripotent stem cells responsible for growth of above-ground tissues in flowering plants. We find that localized perception of the plant hormone cytokinin establishes a spatial domain in which cell fate is respecified through induction of the master regulator WUSCHEL as cells are displaced during growth. Cytokinin-induced WUSCHEL expression occurs through both CLAVATA-dependent and CLAVATA-independent pathways. Computational analysis shows that feedback between cytokinin response and genetic regulators predicts their relative patterning, which we confirm experimentally. Our results also may explain how increasing cytokinin concentration leads to the first steps in reestablishing the shoot stem cell niche in vitro.

Computational modeling of the hematopoietic erythroid-myeloid switch reveals insights into cooperativity, priming, and irreversibility.

Hematopoietic stem cell lineage choices are decided by genetic networks that are turned ON/OFF in a switch-like manner. However, prior to lineage commitment, genes are primed at low expression levels. Understanding the underlying molecular circuitry in terms of how it governs both a primed state and, at the other extreme, a committed state is of relevance not only to hematopoiesis but also to developmental systems in general. We develop a computational model for the hematopoietic erythroid-myeloid lineage decision, which is determined by a genetic switch involving the genes PU.1 and GATA-1. Dynamical models based upon known interactions between these master genes, such as mutual antagonism and autoregulation, fail to make the system bistable, a desired feature for robust lineage determination. We therefore suggest a new mechanism involving a cofactor that is regulated as well as recruited by one of the master genes to bind to the antagonistic partner that is necessary for bistability and hence switch-like behavior. An interesting fallout from this architecture is that suppression of the cofactor through external means can lead to a loss of cooperativity, and hence to a primed state for PU.1 and GATA-1. The PU.1-GATA-1 switch also interacts with another mutually antagonistic pair, C/EBPalpha-FOG-1. The latter pair inherits the state of its upstream master genes and further reinforces the decision due to several feedback loops, thereby leading to irreversible commitment. The genetic switch, which handles the erythroid-myeloid lineage decision, is an example of a network that implements both a primed and a committed state by regulating cooperativity through recruitment of cofactors. Perturbing the feedback between the master regulators and downstream targets suggests potential reprogramming strategies. The approach points to a framework for lineage commitment studies in general and could aid the search for lineage-determining genes.

A computational model for understanding stem cell, trophectoderm and endoderm lineage determination.

BACKGROUND: Recent studies have associated the transcription factors, Oct4, Sox2 and Nanog as parts of a self-regulating network which is responsible for maintaining embryonic stem cell properties: self renewal and pluripotency. In addition, mutual antagonism between two of these and other master regulators have been shown to regulate lineage determination. In particular, an excess of Cdx2 over Oct4 determines the trophectoderm lineage whereas an excess of Gata-6 over Nanog determines differentiation into the endoderm lineage. Also, under/over-expression studies of the master regulator Oct4 have revealed that some self-renewal/pluripotency as well as differentiation genes are expressed in a biphasic manner with respect to the concentration of Oct4. METHODOLOGY/PRINCIPAL FINDINGS: We construct a dynamical model of a minimalistic network, extracted from ChIP-on-chip and microarray data as well as literature studies. The model is based upon differential equations and makes two plausible assumptions; activation of Gata-6 by Oct4 and repression of Nanog by an Oct4-Gata-6 heterodimer. With these assumptions, the results of simulations successfully describe the biphasic behavior as well as lineage commitment. The model also predicts that reprogramming the network from a differentiated state, in particular the endoderm state, into a stem cell state, is best achieved by over-expressing Nanog, rather than by suppression of differentiation genes such as Gata-6. CONCLUSIONS: The computational model provides a mechanistic understanding of how different lineages arise from the dynamics of the underlying regulatory network. It provides a framework to explore strategies of reprogramming a cell from a differentiated state to a stem cell state through directed perturbations. Such an approach is highly relevant to regenerative medicine since it allows for a rapid search over the host of possibilities for reprogramming to a stem cell state.

Oscillatory dynamics arising from competitive inhibition and multisite phosphorylation.

There have been a growing number of observations of oscillating protein levels (p53 and NFkB) in eukaryotic signalling pathways. This has resulted in a renewed interest in the mechanism by which such oscillations might occur. Recent computational work has shown that a multisite phosphorylation mechanism such as that found in the MAPK cascade can theoretically exhibit bistability. The bistable behavior was shown to arise from sequestration and saturation mechanisms for the enzymes that catalyse the multisite phosphorylation cycle. These effects generate the positive feedback necessary for bistability. In this paper we describe two kinds of oscillatory dynamics which can occur in a network by which, both use such bistable multisite phosphorylated cycles. In the first example, the fully phosphorylated form of the phosphorylated cycle represses the production of the kinase, which carries out the phosphorylation of the unphosphorylated states of the cycle. The dynamics of this system leads to a relaxation oscillator. In the second example, we consider a cascade of two cycles, in which the fully phosphorylated form of the kinase, in the first cycle, phosphorylates the unphosphorylated forms in the second cycle. A feedback loop, by which the fully phosphorylated form of the second cycle inhibits the kinase step in the first cycle is also present. In this case we obtain a ring oscillator. Both these networks illustrate the versatility of the multisite bistable network.

Transcriptional dynamics of the embryonic stem cell switch.

Recent ChIP experiments of human and mouse embryonic stem cells have elucidated the architecture of the transcriptional regulatory circuitry responsible for cell determination, which involves the transcription factors OCT4, SOX2, and NANOG. In addition to regulating each other through feedback loops, these genes also regulate downstream target genes involved in the maintenance and differentiation of embryonic stem cells. A search for the OCT4-SOX2-NANOG network motif in other species reveals that it is unique to mammals. With a kinetic modeling approach, we ascribe function to the observed OCT4-SOX2-NANOG network by making plausible assumptions about the interactions between the transcription factors at the gene promoter binding sites and RNA polymerase (RNAP), at each of the three genes as well as at the target genes. We identify a bistable switch in the network, which arises due to several positive feedback loops, and is switched on/off by input environmental signals. The switch stabilizes the expression levels of the three genes, and through their regulatory roles on the downstream target genes, leads to a binary decision: when OCT4, SOX2, and NANOG are expressed and the switch is on, the self-renewal genes are on and the differentiation genes are off. The opposite holds when the switch is off. The model is extremely robust to parameter changes. In addition to providing a self-consistent picture of the transcriptional circuit, the model generates several predictions. Increasing the binding strength of NANOG to OCT4 and SOX2, or increasing its basal transcriptional rate, leads to an irreversible bistable switch: the switch remains on even when the activating signal is removed. Hence, the stem cell can be manipulated to be self-renewing without the requirement of input signals. We also suggest tests that could discriminate between a variety of feedforward regulation architectures of the target genes by OCT4, SOX2, and NANOG.

Conservation analysis of large biochemical networks.

MOTIVATION: Large biochemical networks pose a unique challenge from the point of view of evaluating conservation laws. The computational problem in most cases exceeds the capability of available software tools, often resulting in inaccurate computation of the number and form of conserved cycles. Such errors have profound effects on subsequent calculations, particularly in the evaluation of the Jacobian which is a critical quantity in many other calculations. The goal of this paper is to outline a new algorithm that is computationally efficient and robust at extracting the correct conservation laws for very large biochemical networks. RESULTS: We show that our algorithm can perform the conservation analysis of large biochemical networks, and can evaluate the correct conserved cycles when compared with other similar software tools. Biochemical simulators such as Jarnac and COPASI are successful at extracting only a subset of the conservation laws that our algorithm can. This is illustrated with examples for some large networks which show the advantages of our method.

Bifurcation discovery tool.

MOTIVATION: Biochemical networks often yield interesting behavior such as switching, oscillation and chaotic dynamics. This article describes a tool that is capable of searching for bifurcation points in arbitrary ODE-based reaction networks by directing the user to regions in the parameter space, where such interesting dynamical behavior can be observed. RESULTS: We have implemented a genetic algorithm that searches for Hopf bifurcations, turning points and bistable switches. The software is implemented as a Systems Biology Workbench (SBW) enabled module and accepts the standard SBML model format. The interface permits a user to choose the parameters to be searched, admissible parameter ranges, and the nature of the bifurcation to be sought. The tool will return the parameter values for the model for which the particular behavior is observed. AVAILABILITY: The software, tutorial manual and test models are available for download at the following website: http:/www.sys-bio.org/ under the bifurcation link. The software is an open source and licensed under BSD.

The SBW-MATLAB interface.

SUMMARY: The SBW-MATLAB Interface allows MATLAB users to take advantage of the wide variety of tools available through SBW, the Systems Biology Workbench (Sauro et al. (2003) OMICS, 7, 355-372). It also enables MATLAB users to themselves create SBW-enabled tools which can be freely distributed.