Stochastic diffusion characterises early colony formation in Mediterranean coral Corallium rubrum.

The colony formation in Mediterranean coral Corallium rubrum is initiated by a larva that metamorphoses into the first polyp of the emerging colony approximately two weeks after settlement. The primary polyp then sets up a slow process that eventually, at least after a few years, gives rise to a tree-like rigid colony structure on which other polyps flourish. For a mature colony, this axial skeleton provides support for new polyps. However, the first emergence of the characteristic axial skeleton can take two years or more from the larva stage. The early colony morphology, instead, is shaped exclusively by the polyps' abundant deposition of sclerites, a magnesian calcite biomineral that has a different granularity from the distinctive red-coloured skeleton. With the appearance of the first polyp, a growing sclerite heap in a mesoglea layer provides a base for the emerging colony. In this paper, to elucidate the mechanical processes of early skeleton development in C. rubrum colonies, we present a computational model whereby the mesoglea layer provides a diffusion medium for the sclerites that the polyps deposit. We show that our stochastic model with three parameters captures the dynamic variability observed in measurements on living colonies. Our simulation results provide evidence for a diffusion process whereby the interplay between polyp budding and sclerite deposition are the main determinants of structure in early colony formation. Our model demonstrates that the frequency of budding events in an early colony can be described as a function of the available mesoglea surface whereas the number of polyps on the colony plays a secondary role in determining this frequency. We show that these model predictions are confirmed by direct observations on the colonies in our sample. Moreover, our results indicate that diffusion is a prevalent mechanism of colony development also at later stages of a colony's life span.

The inverse correlation between robustness and sensitivity to autoregulation in two-component systems.

Two-component systems (TCS) are signal transduction systems in bacteria and many other organisms that relay the sensory signal to genetic components. TCS consist of two proteins: a histidine kinase and a response regulator that the histidine kinase activates. This seemingly simple machinery can generate complex regulatory dynamics that enables the level of gene expression that matches the input signal: many TCS response regulators act on their own genes as transcription factors, resulting in a positive autoregulation mechanism. This regulation, in return, modulates the transcription factor activity as a function of the input signal. Positive autoregulation does not necessarily result in positive feedback. Sensitivity to autoregulation is quantified as the output level amplification resulting from the positive autoregulation mechanism. Another structural property of these systems is formally characterized as "robustness": in a robust TCS, the output of the system is solely a function of the input signal. Thus, a robust TCS remains insensitive to fluctuations in the concentrations of its protein components and, this way, maintains the precision in the output transcription factor activity in response to input stimulus. In this paper, we show with a formal model that TCS operate on a spectrum of inverse correlation between robustness and sensitivity to autoregulation. Our model predicts that the modulation by positive autoregulation is a function of loss in TCS robustness, for example, by spontaneous dephosphorylation of the histidine kinase. Consequently, the loss in robustness provides a proportional modulation by positive autoregulation to widen the response range with a scaled amplification of the output. At the other end of the spectrum, in the presence of a strictly robust TCS machinery, amplification of the transcription factor activity by autoregulation is diminished. We show that our results are in agreement with published experimental results. Our results suggest that these TCS evolve to converge at a trade-off between robustness and positive autoregulation.

Stochastic Simulations as a Tool for Assessing Signal Fidelity in Gene Expression in Synthetic Promoter Design.

The design and development of synthetic biology applications in a workflow often involve connecting modular components. Whereas computer-aided design tools are picking up in synthetic biology as in other areas of engineering, the methods for verifying the correct functioning of living technologies are still in their infancy. Especially, fine-tuning for the right promoter strength to match the design specifications is often a lengthy and expensive experimental process. In particular, the relationship between signal fidelity and noise in synthetic promoter design can be a key parameter that can affect the healthy functioning of the engineered organism. To this end, based on our previous work on synthetic promoters for the E. coli PhoBR two-component system, we make a case for using chemical reaction network models for computational verification of various promoter designs before a lab implementation. We provide an analysis of this system with extensive stochastic simulations at a single-cell level to assess the signal fidelity and noise relationship. We then show how quasi-steady-state analysis via ordinary differential equations can be used to navigate between models with different levels of detail. We compare stochastic simulations with our full and reduced models by using various metrics for assessing noise. Our analysis suggests that strong promoters with low unbinding rates can act as control tools for filtering out intrinsic noise in the PhoBR context. Our results confirm that even simpler models can be used to determine promoters with specific signal to noise characteristics.

Quantifying dynamic mechanisms of auto-regulation in Escherichia coli with synthetic promoter in response to varying external phosphate levels.

Escherichia coli have developed one of the most efficient regulatory response mechanisms to phosphate starvation. The machinery involves a cascade with a two-component system (TCS) that relays the external signal to the genetic circuit, resulting in a feedback response. Achieving a quantitative understanding of this system has implications in synthetic biology and biotechnology, for example, in applications for wastewater treatment. To this aim, we present a computational model and experimental results with a detailed description of the TCS, consisting of PhoR and PhoB, together with the mechanisms of gene expression. The model is parameterised within the feasible range, and fitted to the dynamic response of our experimental data on PhoB as well as PhoA, the product of this network that is used in alkaline phosphatase production. Deterministic and stochastic simulations with our model predict the regulation dynamics in higher external phosphate concentrations while reproducing the experimental observations. In a cycle of simulations and experimental verification, our model predicts and explores phenotypes with various synthetic promoter designs that can optimise the inorganic phosphate intake in E. coli. Sensitivity analysis demonstrates that the Pho-controlled genes have a significant influence over the phosphate response. Together with experimental findings, our model should thus provide insights for the investigations on engineering new sensors and regulators for living technologies.

Mechanistic interplay between ceramide and insulin resistance.

Recent research adds to a growing body of literature on the essential role of ceramides in glucose homeostasis and insulin signaling, while the mechanistic interplay between various components of ceramide metabolism remains to be quantified. We present an extended model of C16:0 ceramide production through both the de novo synthesis and the salvage pathways. We verify our model with a combination of published models and independent experimental data. In silico experiments of the behavior of ceramide and related bioactive lipids in accordance with the observed transcriptomic changes in obese/diabetic murine macrophages at 5 and 16 weeks support the observation of insulin resistance only at the later phase. Our analysis suggests the pivotal role of ceramide synthase, serine palmitoyltransferase and dihydroceramide desaturase involved in the de novo synthesis and the salvage pathways in influencing insulin resistance versus its regulation.

Gener: a minimal programming module for chemical controllers based on DNA strand displacement.

UNLABELLED: : Gener is a development module for programming chemical controllers based on DNA strand displacement. Gener is developed with the aim of providing a simple interface that minimizes the opportunities for programming errors: Gener allows the user to test the computations of the DNA programs based on a simple two-domain strand displacement algebra, the minimal available so far. The tool allows the user to perform stepwise computations with respect to the rules of the algebra as well as exhaustive search of the computation space with different options for exploration and visualization. Gener can be used in combination with existing tools, and in particular, its programs can be exported to Microsoft Research's DSD tool as well as to LaTeX. AVAILABILITY AND IMPLEMENTATION: Gener is available for download at the Cosbi website at http://www.cosbi.eu/research/prototypes/gener as a windows executable that can be run on Mac OS X and Linux by using Mono. CONTACT: ozan@cosbi.eu.

Stochastic flux analysis of chemical reaction networks.

BACKGROUND: Chemical reaction networks provide an abstraction scheme for a broad range of models in biology and ecology. The two common means for simulating these networks are the deterministic and the stochastic approaches. The traditional deterministic approach, based on differential equations, enjoys a rich set of analysis techniques, including a treatment of reaction fluxes. However, the discrete stochastic simulations, which provide advantages in some cases, lack a quantitative treatment of network fluxes. RESULTS: We describe a method for flux analysis of chemical reaction networks, where flux is given by the flow of species between reactions in stochastic simulations of the network. Extending discrete event simulation algorithms, our method constructs several data structures, and thereby reveals a variety of statistics about resource creation and consumption during the simulation. We use these structures to quantify the causal interdependence and relative importance of the reactions at arbitrary time intervals with respect to the network fluxes. This allows us to construct reduced networks that have the same flux-behavior, and compare these networks, also with respect to their time series. We demonstrate our approach on an extended example based on a published ODE model of the same network, that is, Rho GTP-binding proteins, and on other models from biology and ecology. CONCLUSIONS: We provide a fully stochastic treatment of flux analysis. As in deterministic analysis, our method delivers the network behavior in terms of species transformations. Moreover, our stochastic analysis can be applied, not only at steady state, but at arbitrary time intervals, and used to identify the flow of specific species between specific reactions. Our cases study of Rho GTP-binding proteins reveals the role played by the cyclic reverse fluxes in tuning the behavior of this network.

Algorithmic modeling quantifies the complementary contribution of metabolic inhibitions to gemcitabine efficacy.

Gemcitabine (2,2-difluorodeoxycytidine, dFdC) is a prodrug widely used for treating various carcinomas. Gemcitabine exerts its clinical effect by depleting the deoxyribonucleotide pools, and incorporating its triphosphate metabolite (dFdC-TP) into DNA, thereby inhibiting DNA synthesis. This process blocks the cell cycle in the early S phase, eventually resulting in apoptosis. The incorporation of gemcitabine into DNA takes place in competition with the natural nucleoside dCTP. The mechanisms of indirect competition between these cascades for common resources are given with the race for DNA incorporation; in clinical studies dedicated to singling out mechanisms of resistance, ribonucleotide reductase (RR) and deoxycytidine kinase (dCK) and human equilibrative nucleoside transporter1 (hENT1) have been associated to efficacy of gemcitabine with respect to their roles in the synthesis cascades of dFdC-TP and dCTP. However, the direct competition, which manifests itself in terms of inhibitions between these cascades, remains to be quantified. We propose an algorithmic model of gemcitabine mechanism of action, verified with respect to independent experimental data. We performed in silico experiments in different virtual conditions, otherwise difficult in vivo, to evaluate the contribution of the inhibitory mechanisms to gemcitabine efficacy. In agreement with the experimental data, our model indicates that the inhibitions due to the association of dCTP with dCK and the association of gemcitabine diphosphate metabolite (dFdC-DP) with RR play a key role in adjusting the efficacy. While the former tunes the catalysis of the rate-limiting first phosphorylation of dFdC, the latter is responsible for depletion of dCTP pools, thereby contributing to gemcitabine efficacy with a dependency on nucleoside transport efficiency. Our simulations predict the existence of a continuum of non-efficacy to high-efficacy regimes, where the levels of dFdC-TP and dCTP are coupled in a complementary manner, which can explain the resistance to this drug in some patients.

Biophysical mechanism for ras-nanocluster formation and signaling in plasma membrane.

Ras GTPases are lipid-anchored G proteins, which play a fundamental role in cell signaling processes. Electron micrographs of immunogold-labeled Ras have shown that membrane-bound Ras molecules segregate into nanocluster domains. Several models have been developed in attempts to obtain quantitative descriptions of nanocluster formation, but all have relied on assumptions such as a constant, expression-level independent ratio of Ras in clusters to Ras monomers (cluster/monomer ratio). However, this assumption is inconsistent with the law of mass action. Here, we present a biophysical model of Ras clustering based on short-range attraction and long-range repulsion between Ras molecules in the membrane. To test this model, we performed Monte Carlo simulations and compared statistical clustering properties with experimental data. We find that we can recover the experimentally-observed clustering across a range of Ras expression levels, without assuming a constant cluster/monomer ratio or the existence of lipid rafts. In addition, our model makes predictions about the signaling properties of Ras nanoclusters in support of the idea that Ras nanoclusters act as an analog-digital-analog converter for high fidelity signaling.