Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStochastic Simulation of Biomolecular Networks in Dynamic Environments
131
0
0.0
(
0
)
Added by
Margaritis Voliotis
Publication date
2015
fields
Biology
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Simulation of biomolecular networks is now indispensable for studying biological systems, from small reaction networks to large ensembles of cells. Here we present a novel approach for stochastic simulation of networks embedded in the dynamic environment of the cell and its surroundings. We thus sample trajectories of the stochastic process described by the chemical master equation with time-varying propensities. A comparative analysis shows that existing approaches can either fail dramatically, or else can impose impractical computational burdens due to numerical integration of reaction propensities, especially when cell ensembles are studied. Here we introduce the Extrande method which, given a simulated time course of dynamic network inputs, provides a conditionally exact and several orders-of-magnitude faster simulation solution. The new approach makes it feasible to demonstrate, using decision-making by a large population of quorum sensing bacteria, that robustness to fluctuations from upstream signaling places strong constraints on the design of networks determining cell fate. Our approach has the potential to significantly advance both understanding of molecular systems biology and design of synthetic circuits.
rate research
Read More
Infection by many viruses begins with fusion of viral and cellular lipid membranes, followed by entry of viral contents into the target cell and ultimately, after many biochemical steps, integration of viral DNA into that of the host cell. The early steps of membrane fusion and viral capsid entry are mediated by adsorption to the cell surface, and receptor and coreceptor binding. HIV-1 specifically targets CD4+ helper T-cells of the human immune system and binds to the receptor CD4 and coreceptor CCR5 before fusion is initiated. Previous experiments have been performed using a cell line (293-Affinofile) in which the expression of CD4 and CCR5 concentration were independently controlled. After exposure to HIV-1 of various strains, the resulting infectivity was measured through the fraction of infected cells. To design and evaluate the effectiveness of drug therapies that target the inhibition of the entry processes, an accurate functional relationship between the CD4/CCR5 concentrations and infectivity is desired in order to more quantitatively analyze experimental data. We propose three kinetic models describing the possible mechanistic processes involved in HIV entry and fit their predictions to infectivity measurements, contrasting and comparing different outcomes. Our approach allows interpretation of the clustering of infectivity of different strains of HIV-1 in the space of mechanistic kinetic parameters. Our model fitting also allows inference of nontrivial stoichiometries of receptor and coreceptor binding and provides a framework through which to quantitatively investigate the effectiveness of fusion inhibitors and neutralizing antibodies.
Systems biology and whole-cell modelling are demanding increasingly comprehensive mathematical models of cellular biochemistry. These models require the development of simplified models of specific processes which capture essential biophysical features but without unnecessarily complexity. Recently there has been renewed interest in thermodynamically-based modelling of cellular processes. Here we present an approach to developing of simplified yet thermodynamically consistent (hence physically plausible) models which can readily be incorporated into large scale biochemical descriptions but which do not require full mechanistic detail of the underlying processes. We illustrate the approach through development of a simplified, physically plausible model of the mitochondrial electron transport chain and show that the simplified model behaves like the full system.
Biological cells are often found to sense their chemical environment near the single-molecule detection limit. Surprisingly, this precision is higher than simple estimates of the fundamental physical limit, hinting towards active sensing strategies. In this work, we analyse the effect of cell memory, e.g. from slow biochemical processes, on the precision of sensing by cell-surface receptors. We derive analytical formulas, which show that memory significantly improves sensing in weakly fluctuating environments. However, surprisingly when memory is adjusted dynamically, the precision is always improved, even in strongly fluctuating environments. In support of this prediction we quantify the directional biases in chemotactic Dictyostelium discoideum cells in a flow chamber with alternating chemical gradients. The strong similarities between cell sensing and control engineering suggest universal problem-solving strategies of living matter.
Stochasticity is an indispensable aspect of biochemical processes at the cellular level. Studies on how the noise enters and propagates in biochemical systems provided us with nontrivial insights into the origins of stochasticity, in total however they constitute a patchwork of different theoretical analyses. Here we present a flexible and generally applicable noise decomposition tool, that allows us to calculate contributions of individual reactions to the total variability of a systems output. With the package it is therefore possible to quantify how the noise enters and propagates in biochemical systems. We also demonstrate and exemplify using the JAK-STAT signalling pathway that it is possible to infer noise contributions resulting from individual reactions directly from experimental data. This is the first computational tool that allows to decompose noise into contributions resulting from individual reactions.
Systems Biology is a fundamental field and paradigm that introduces a new era in Biology. The crux of its functionality and usefulness relies on metabolic networks that model the reactions occurring inside an organism and provide the means to understand the underlying mechanisms that govern biological systems. Even more, metabolic networks have a broader impact that ranges from resolution of ecosystems to personalized medicine.The analysis of metabolic networks is a computational geometry oriented field as one of the main operations they depend on is sampling uniformly points from polytopes; the latter provides a representation of the steady states of the metabolic networks. However, the polytopes that result from biological data are of very high dimension (to the order of thousands) and in most, if not all, the cases are considerably skinny. Therefore, to perform uniform random sampling efficiently in this setting, we need a novel algorithmic and computational framework specially tailored for the properties of metabolic networks.We present a complete software framework to handle sampling in metabolic networks. Its backbone is a Multiphase Monte Carlo Sampling (MMCS) algorithm that unifies rounding and sampling in one pass, obtaining both upon termination. It exploits an improved variant of the Billiard Walk that enjoys faster arithmetic complexity per step. We demonstrate the efficiency of our approach by performing extensive experiments on various metabolic networks. Notably, sampling on the most complicated human metabolic network accessible today, Recon3D, corresponding to a polytope of dimension 5 335 took less than 30 hours. To our knowledge, that is out of reach for existing software.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
