Modeling the segmentation clock as a network of coupled oscillations in the Notch, Wnt and FGF signaling pathways

The formation of somites in the course of vertebrate segmentation is governed by an oscillator known as the segmentation clock, which is characterized by a period ranging from 30 min to a few hours depending on the organism. This oscillator permits the synchronized activation of segmentation genes in successive cohorts of cells in the presomitic mesoderm in response to a periodic signal emitted by the segmentation clock, thereby defining the future segments. Recent microarray experiments [Dequeant, M.L., Glynn, E., Gaudenz, K., Wahl, M., Chen, J., Mushegian, A., Pourquie, O., 2006. A complex oscillating network of signaling genes underlies the mouse segmentation clock. Science 314, 1595-1598] indicate that the Notch, Wnt and Fibroblast Growth Factor (FGF) signaling pathways are involved in the mechanism of the segmentation clock. By means of computational modeling, we investigate the conditions in which sustained oscillations occur in these three signaling pathways. First we show that negative feedback mediated by the Lunatic Fringe protein on intracellular Notch activation can give rise to periodic behavior in the Notch pathway. We then show that negative feedback exerted by Axin2 on the degradation of β-catenin through formation of the Axin2 destruction complex can produce oscillations in the Wnt pathway. Likewise, negative feedback on FGF signaling mediated by the phosphatase product of the gene MKP3/Dusp6 can produce oscillatory gene expression in the FGF pathway. Coupling the Wnt, Notch and FGF oscillators through common intermediates can lead to synchronized oscillations in the three signaling pathways or to complex periodic behavior, depending on the relative periods of oscillations in the three pathways. The phase relationships between cycling genes in the three pathways depend on the nature of the coupling between the pathways and on their relative autonomous periods. The model provides a framework for analyzing the dynamics of the segmentation clock in terms of a network of oscillating modules involving the Wnt, Notch and FGF signaling pathways.     

Circadian clock effects on cellular proliferation: Insights from theory and experiments

Oscillations of the cellular circadian clock have emerged as an important regulator of many physiological processes, both in health and in disease. One such process, cellular proliferation, is being increasingly recognized to be affected by the circadian clock. Here, we review how a combination of experimental and theoretical work has furthered our understanding of the way circadian clocks couple to the cell cycle and play a role in tissue homeostasis and cancer. Finally, we discuss recently introduced methods for modeling coupling of clocks based on techniques from survival analysis and machine learning and highlight their potential importance for future studies.     

A spatio-temporal model of Notch signalling in the zebrafish segmentation clock: conditions for synchronised oscillatory dynamics

In the vertebrate embryo, tissue blocks called somites are laid down in head-to-tail succession, a process known as somitogenesis. Research into somitogenesis has been both experimental and mathematical. For zebrafish, there is experimental evidence for oscillatory gene expression in cells in the presomitic mesoderm (PSM) as well as evidence that Notch signalling synchronises the oscillations in neighbouring PSM cells. A biological mechanism has previously been proposed to explain these phenomena. Here we have converted this mechanism into a mathematical model of partial differential equations in which the nuclear and cytoplasmic diffusion of protein and mRNA molecules is explicitly considered. By performing simulations, we have found ranges of values for the model parameters (such as diffusion and degradation rates) that yield oscillatory dynamics within PSM cells and that enable Notch signalling to synchronise the oscillations in two touching cells. Our model contains a Hill coefficient that measures the co-operativity between two proteins (Her1, Her7) and three genes (her1, her7, deltaC) which they inhibit. This coefficient appears to be bounded below by the requirement for oscillations in individual cells and bounded above by the requirement for synchronisation. Consistent with experimental data and a previous spatially non-explicit mathematical model, we have found that signalling can increase the average level of Her1 protein. Biological pattern formation would be impossible without a certain robustness to variety in cell shape and size; our results possess such robustness. Our spatially-explicit modelling approach, together with new imaging technologies that can measure intracellular protein diffusion rates, is likely to yield significant new insight into somitogenesis and other biological processes.     

simBio: a Java package for the development of detailed cell models

Quantitative dynamic computer models, which integrate a variety of molecular functions into a cell model, provide a powerful tool to create and test working hypotheses. We have developed a new modeling tool, the simBio package (freely available from http://www.sim-bio.org/), which can be used for constructing cell models, such as cardiac cells (the Kyoto model from Matsuoka et al., 2003, 2004a, b, the LRd model from Faber and Rudy, 2000, and the Noble 98 model from Noble et al., 1998), epithelial cells (Strieter et al., 1990) and pancreatic β cells (Magnus and Keizer, 1998). The simBio package is written in Java, uses XML and can solve ordinary differential equations. In an attempt to mimic biological functional structures, a cell model is, in simBio, composed of independent functional modules called Reactors, such as ion channels and the sarcoplasmic reticulum, and dynamic variables called Nodes, such as ion concentrations. The interactions between Reactors and Nodes are described by the graph theory and the resulting graph represents a blueprint of an intricate cellular system. Reactors are prepared in a hierarchical order, in analogy to the biological classification. Each Reactor can be composed or improved independently, and can easily be reused for different models. This way of building models, through the combination of various modules, is enabled through the use of object-oriented programming concepts. Thus, simBio is a straightforward system for the creation of a variety of cell models on a common database of functional modules.     

Targeting stem cell signaling pathways for drug discovery: advances in the Notch and Wnt pathways

Signaling pathways transduce extracellular stimuli into cells through molecular cascades to regulate cellular functions.In stem cells,a small number of pathways,notably those of TGF-?/BMP,Hedgehog,Notch,and Wnt,are responsible for the regulation of pluripotency and differentiation.During embryonic development,these pathways govern cell fate specifications as well as the formation of tissues and organs.In adulthood,their normal functions are important for tissue homeostasis and regeneration,whereas aberrations result in diseases,such as cancer and degenerative disorders.In complex biological systems,stem cell signaling pathways work in concert as a network and exhibit crosstalk,such as the negative crosstalk between Wnt and Notch.Over the past decade,genetic and genomic studies have identified a number of potential drug targets that are involved in stem cell signaling pathways.Indeed,discovery of new targets and drugs for these pathways has become one of the most active areas in both the research community and pharmaceutical industry.Remarkable progress has been made and several promising drug candidates have entered into clinical trials.This review focuses on recent advances in the discovery of novel drugs which target the Notch and Wnt pathways.     

Modeling the mammalian circadian clock: sensitivity analysis and multiplicity of oscillatory mechanisms

We extend the study of a computational model recently proposed for the mammalian circadian clock (Proc. Natl Acad. Sci. USA 100 (2003) 7051). The model, based on the intertwined positive and negative regulatory loops involving the Per, Cry, Bmal1, and Clock genes, can give rise to sustained circadian oscillations in conditions of continuous darkness. These limit cycle oscillations correspond to circadian rhythms autonomously generated by suprachiasmatic nuclei and by some peripheral tissues. By using different sets of parameter values producing circadian oscillations, we compare the effect of the various parameters and show that both the occurrence and the period of the oscillations are generally most sensitive to parameters related to synthesis or degradation of Bmal1 mRNA and BMAL1 protein. The mechanism of circadian oscillations relies on the formation of an inactive complex between PER and CRY and the activators CLOCK and BMAL1 that enhance Per and Cry expression. Bifurcation diagrams and computer simulations nevertheless indicate the possible existence of a second source of oscillatory behavior. Thus, sustained oscillations might arise from the sole negative autoregulation of Bmal1 expression. This second oscillatory mechanism may not be functional in physiological conditions, and its period need not necessarily be circadian. When incorporating the light-induced expression of the Per gene, the model accounts for entrainment of the oscillations by light-dark (LD) cycles. Long-term suppression of circadian oscillations by a single light pulse can occur in the model when a stable steady state coexists with a stable limit cycle. The phase of the oscillations upon entrainment in LD critically depends on the parameters that govern the level of CRY protein. Small changes in the parameters governing CRY levels can shift the peak in Per mRNA from the L to the D phase, or can prevent entrainment. The results are discussed in relation to physiological disorders of the sleep-wake cycle linked to perturbations of the human circadian clock, such as the familial advanced sleep phase syndrome or the non-24h sleep-wake syndrome.     

Interactions between Wnt and Vg1 signalling pathways initiate primitive streak formation in the chick embryo

The posterior marginal zone (PMZ) of the chick embryo has Nieuwkoop centre-like properties: when transplanted to another part of the marginal zone, it induces a complete embryonic axis, without making a cellular contribution to the induced structures. However, when the PMZ is removed, the embryo can initiate axis formation from another part of the remaining marginal zone. Chick Vg1 can mimic the axis-inducing ability of the PMZ, but only when misexpressed somewhere within the marginal zone. We have investigated the properties that define the marginal zone as a distinct region. We show that the competence of the marginal zone to initiate ectopic primitive streak formation in response to cVg1 is dependent on Wnt activity. First, within the Wnt family, only Wnt8C is expressed in the marginal zone, in a gradient decreasing from posterior to anterior. Second, misexpression of Wnt1 in the area pellucida enables this region to form a primitive streak in response to cVg1. Third, the Wnt antagonists Crescent and Dkk-1 block the primitive streak-inducing ability of cVg1 in the marginal zone. These findings suggest that Wnt activity defines the marginal zone and allows cVg1 to induce an axis. We also present data suggesting some additional complexity: first, the Vg1 and Wnt pathways appear to regulate the expression of downstream components of each other's pathway; and second, misexpression of different Wnt antagonists suggests that different classes of Wnts may cooperate with each other to regulate axis formation in the normal embryo.     

The dynamics of single-substrate continuous cultures: the role of transport enzymes

A chemostat limited by a single growth-limiting substrate displays a rich spectrum of dynamics. Depending on the flow rate and feed concentration, the chemostat settles into a steady state or executes sustained oscillations. The transients in response to abrupt increases in the flow rate or the feed concentration are also quite complex. For example, if the increase in the flow rate is small, there is no perceptible change in the substrate concentration. If the increase in the flow rate is large, there is a large increase in the substrate concentration lasting several hours or days before the culture adjusts to a new steady state. In the latter case, the substrate concentration and cell density frequently undergo damped oscillations during their approach to the steady state. In this work, we formulate a simple structured model containing the inducible transport enzyme as the key intracellular variable. The model displays the foregoing dynamics under conditions similar to those employed in the experiments. The model suggests that long recovery times (on the order of several hours to several days) can occur because the initial transport enzyme level is too small to cope with the increased substrate supply. The substrate concentration, therefore, increases until the enzyme level is built up to a sufficiently high level by the slow process of enzyme induction. Damped and sustained oscillations can occur because transport enzyme synthesis is autocatalytic, and hence, destabilizing. At low dilution rates, the response of stabilizing processes, such as enzyme dilution and substrate consumption, becomes very slow, leading to damped and sustained oscillations.     

Extracellular regulation of VEGF: Isoforms,proteolysis, and vascular patterning

The regulation of vascular endothelial growth factor A (VEGF) is critical to neovascularization in numerous tissues under physiological and pathological conditions. VEGF has multiple isoforms, created by alternative splicing or proteolytic cleavage, and characterized by different receptor-binding and matrix-binding properties. These isoforms are known to give rise to a spectrum of angiogenesis patterns marked by differences in branching, which has functional implications for tissues. In this review, we detail the extensive extracellular regulation of VEGF and the ability of VEGF to dictate the vascular phenotype. We explore the role of VEGF-releasing proteases and soluble carrier molecules on VEGF activity. While proteases such as MMP9 can ‘release’ matrix-bound VEGF and promote angiogenesis, for example as a key step in carcinogenesis, proteases can also suppress VEGF's angiogenic effects. We explore what dictates pro- or anti-angiogenic behavior. We also seek to understand the phenomenon of VEGF gradient formation. Strong VEGF gradients are thought to be due to decreased rates of diffusion from reversible matrix binding, however theoretical studies show that this scenario cannot give rise to lasting VEGF gradients in vivo. We propose that gradients are formed through degradation of sequestered VEGF. Finally, we review how different aspects of the VEGF signal, such as its concentration, gradient, matrix-binding, and NRP1-binding can differentially affect angiogenesis. We explore how this allows VEGF to regulate the formation of vascular networks across a spectrum of high to low branching densities, and from normal to pathological angiogenesis. A better understanding of the control of angiogenesis is necessary to improve upon limitations of current angiogenic therapies.     

Cycling expression and cooperative operator interaction in the trp operon of Escherichia coli

Oscillatory behaviour in the tryptophan operon of an Escherichia coli mutant strain lacking the enzyme-inhibition regulatory mechanism has been observed by Bliss et al. but not confirmed by others. This behaviour could be important from the standpoint of synthetic biology, whose goals include the engineering of intracellular genetic oscillators. This work is devoted to investigating, from a mathematical modelling point of view, the possibility that the trp operon of the E. coli inhibition-free strain expresses cyclically. For that we extend a previously introduced model for the regulatory pathway of the tryptophan operon in Escherichia coli to account for the observed multiplicity and cooperativity of repressor binding sites. Thereafter we investigate the model dynamics using deterministic numeric solutions, stochastic simulations, and analytic studies. Our results suggest that a quasi-periodic behaviour could be observed in the trp operon expression level of single bacteria.     

A computational model of how an interaction between the thalamocortical and thalamic reticular neurons transforms the low-frequency oscillations of the globus pallidus

In Parkinson’s disease, neurons of the internal segment of the globus pallidus (GPi) display the low-frequency tremor-related oscillations. These oscillatory activities are transmitted to the thalamic relay nuclei. Computer models of the interacting thalamocortical (TC) and thalamic reticular (RE) neurons were used to explore how the TC-RE network processes the low-frequency oscillations of the GPi neurons. The simulation results show that, by an interaction between the TC and RE neurons, the TC-RE network transforms a low-frequency oscillatory activity of the GPi neurons to a higher frequency of oscillatory activity of the TC neurons (the superharmonic frequency transformation). In addition to the interaction between the TC and RE neurons, the low-threshold calcium current in the RE and TC neurons and the hyperpolarization-activated cation current (I _h) in the TC neurons have significant roles in the superharmonic frequency transformation property of the TC-RE network. The external globus pallidus (GPe) oscillatory activity, which is directly transmitted to the RE nucleus also displays a significant modulatory effect on the superharmonic frequency transformation property of the TC-RE network. Action Editor: John Rinzel     

Interplay between activator-inhibitor coupling and cell-matrix adhesion in a cellular automaton model for chondrogenic patterning

We present a stochastic cellular automaton model for the behavior of limb bud precartilage mesenchymal cells undergoing chondrogenic patterning. This "agent-oriented" model represents cells by points on a lattice that obey rules motivated by experimental findings. The "cells" follow these rules as autonomous agents, interacting with other cells and with the microenvironments cell activities produce. The rules include random cell motion, production and lateral deposition of a substrate adhesion molecule (SAM, corresponding to fibronectin), production and release of a diffusible growth factor ("activator," corresponding to TGF-beta) that stimulates production of the SAM, and another diffusible factor ("inhibitor") that suppresses the activity of the activator. We implemented the cellular automaton on a two-dimensional (2D) square lattice to emulate the quasi-2D micromass culture extensively used to study patterning in avian limb bud precartilage cells. We identified parameters that produce nodular patterns that resemble, in size and distribution, cell condensations in leg-cell cultures, thus establishing a correspondence between in vitro and in silico results. We then studied the in vitro and in silico micromass cultures experimentally. We altered the standard in vitro micromass culture by diluting the initial cell density, transiently exposing it to exogenous activator, suppressing the inhibitor, and constitutively activating fibronectin production. We altered the standard in silico micromass culture in each case by changing the corresponding parameter. In vitro and in silico experiments agreed well. We also used the model to test hypotheses for differences in the in vitro patterns of cells derived from chick embryo forelimb and hindlimb. We discuss the applicability of this model to limb development in vivo and to other organ development.     

Similarities and differences between the Wnt and reelin pathways in the forming brain

One of the key features in development is the reutilization of successful signaling pathways. Here, we emphasize the involvement of the Wnt pathway, one of the five kinds of signal transduction pathway predominating early embryonic development of all animals, in regulating the formation of brain structure. We discuss the interrelationships between the Wnt and reelin pathways in the regulation of cortical layering. We summarize data emphasizing key molecules, which, when mutated, result in abnormal brain development. This integrated view, which is based on conservation of pathways, reveals the relative position of participants in the pathway, points to control mechanisms, and allows raising testable working hypotheses. Nevertheless, although signaling pathways are highly conserved from flies to humans, the overall morphology is not. We propose that future studies directed at understanding of diversification will provide fruitful insights on mammalian brain formation.     

The interaction between urokinase receptor and vitronectin in cell adhesion and signalling

The extracellular matrix (ECM) is a complex structural entity surrounding and supporting cells present in all tissue and organs. Cell-matrix interactions play fundamental roles during embryonic development, morphogenesis, tissue homoeostasis, wound healing, and tumourigenesis. Cell-matrix communication is kept in balance by physical contact and by transmembrane integrin receptors providing the dynamic link between the extracellular and intracellular environments through bi-directional signalling. The urokinase-type plasminogen activator receptor (uPAR) is a plasma membrane receptor overexpressed during inflammation and in almost all human cancers. One of its functions is to endorse ECM remodelling through the activation of plasminogen and downstream proteases, including matrix-metalloproteases (MMPs). Beside its role in ECM degradation, uPAR modulates cell-matrix contact through a direct engagement with the ECM component, vitronectin (Vn), and by regulating the activity state of integrins thus promoting or inhibiting integrin signalling and integrin-mediated cell adhesion to other ECM components, like fibronectin and collagen. In this review we have centred our attention on the non-proteolytic function of uPAR as a mediator of cell adhesion and downstream signalling.     

The interaction between dimethylsulfoxide and nitrate reducing pathways in Rhodobacter capsulatus

Abstract The interaction between nitrate- and dimethyl-sulphoxide (DMSO)-reducing pathways was demonstrated in intact cells of Rhodobacter capsulatus AD2 removed from cultures grown under different conditions. The results provide evidence of competition between the DMSO and nitrate reductases for a common electron pool. Furthermore, strong inhibition was observed of the anaerobic dark DMSO-dependent growth of R. capsulatus by nitrate in the growth medium. This phenomenon is also discussed.