Dataset Information

Bioelectric gene and reaction networks: computational modelling of genetic, biochemical and bioelectrical dynamics in pattern regulation.

ABSTRACT: Gene regulatory networks (GRNs) describe interactions between gene products and transcription factors that control gene expression. In combination with reaction-diffusion models, GRNs have enhanced comprehension of biological pattern formation. However, although it is well known that biological systems exploit an interplay of genetic and physical mechanisms, instructive factors such as transmembrane potential (Vmem) have not been integrated into full GRN models. Here we extend regulatory networks to include bioelectric signalling, developing a novel synthesis: the bioelectricity-integrated gene and reaction (BIGR) network. Using in silico simulations, we highlight the capacity for Vmem to alter steady-state concentrations of key signalling molecules inside and out of cells. We characterize fundamental feedbacks where Vmem both controls, and is in turn regulated by, biochemical signals and thereby demonstrate Vmem homeostatic control, Vmem memory and Vmem controlled state switching. BIGR networks demonstrating hysteresis are identified as a mechanisms through which more complex patterns of stable Vmem spots and stripes, along with correlated concentration patterns, can spontaneously emerge. As further proof of principle, we present and analyse a BIGR network model that mechanistically explains key aspects of the remarkable regenerative powers of creatures such as planarian flatworms. The functional properties of BIGR networks generate the first testable, quantitative hypotheses for biophysical mechanisms underlying the stability and adaptive regulation of anatomical bioelectric pattern.

SUBMITTER: Pietak A

PROVIDER: S-EPMC5636277 | biostudies-literature | 2017 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Bioelectric gene and reaction networks: computational modelling of genetic, biochemical and bioelectrical dynamics in pattern regulation.

Pietak Alexis A Levin Michael M

Journal of the Royal Society, Interface 20170901 134

Gene regulatory networks (GRNs) describe interactions between gene products and transcription factors that control gene expression. In combination with reaction-diffusion models, GRNs have enhanced comprehension of biological pattern formation. However, although it is well known that biological systems exploit an interplay of genetic and physical mechanisms, instructive factors such as transmembrane potential (<i>V</i><sub>mem</sub>) have not been integrated into full GRN models. Here we extend ...[more]

PMID: 28954851

Similar Datasets

Project description:Complex molecular biological processes such as transcription and translation, signal transduction, post-translational modification cascades, and metabolic pathways can be described in principle by biochemical reactions that explicitly take into account the sophisticated network of chemical interactions regulating cell life. The ability to deduce the possible qualitative behaviors of such networks from a set of reactions is a central objective and an ongoing challenge in the field of systems biology. Unfortunately, the construction of complete mathematical models is often hindered by a pervasive problem: despite the wealth of qualitative graphical knowledge about network interactions, the form of the governing nonlinearities and/or the values of kinetic constants are hard to uncover experimentally. The kinetics can also change with environmental variations. This work addresses the following question: given a set of reactions and without assuming a particular form for the kinetics, what can we say about the asymptotic behavior of the network? Specifically, it introduces a class of networks that are "structurally (mono) attractive" meaning that they are incapable of exhibiting multiple steady states, oscillation, or chaos by virtue of their reaction graphs. These networks are characterized by the existence of a universal energy-like function called a Robust Lyapunov function (RLF). To find such functions, a finite set of rank-one linear systems is introduced, which form the extremals of a linear convex cone. The problem is then reduced to that of finding a common Lyapunov function for this set of extremals. Based on this characterization, a computational package, Lyapunov-Enabled Analysis of Reaction Networks (LEARN), is provided that constructs such functions or rules out their existence. An extensive study of biochemical networks demonstrates that LEARN offers a new unified framework. Basic motifs, three-body binding, and genetic networks are studied first. The work then focuses on cellular signalling networks including various post-translational modification cascades, phosphotransfer and phosphorelay networks, T-cell kinetic proofreading, and ERK signalling. The Ribosome Flow Model is also studied.

Project description:BACKGROUND: In this paper we propose a model reduction method for biochemical reaction networks governed by a variety of reversible and irreversible enzyme kinetic rate laws, including reversible Michaelis-Menten and Hill kinetics. The method proceeds by a stepwise reduction in the number of complexes, defined as the left and right-hand sides of the reactions in the network. It is based on the Kron reduction of the weighted Laplacian matrix, which describes the graph structure of the complexes and reactions in the network. It does not rely on prior knowledge of the dynamic behaviour of the network and hence can be automated, as we demonstrate. The reduced network has fewer complexes, reactions, variables and parameters as compared to the original network, and yet the behaviour of a preselected set of significant metabolites in the reduced network resembles that of the original network. Moreover the reduced network largely retains the structure and kinetics of the original model. RESULTS: We apply our method to a yeast glycolysis model and a rat liver fatty acid beta-oxidation model. When the number of state variables in the yeast model is reduced from 12 to 7, the difference between metabolite concentrations in the reduced and the full model, averaged over time and species, is only 8%. Likewise, when the number of state variables in the rat-liver beta-oxidation model is reduced from 42 to 29, the difference between the reduced model and the full model is 7.5%. CONCLUSIONS: The method has improved our understanding of the dynamics of the two networks. We found that, contrary to the general disposition, the first few metabolites which were deleted from the network during our stepwise reduction approach, are not those with the shortest convergence times. It shows that our reduction approach performs differently from other approaches that are based on time-scale separation. The method can be used to facilitate fitting of the parameters or to embed a detailed model of interest in a more coarse-grained yet realistic environment.

Dataset Information

Bioelectric gene and reaction networks: computational modelling of genetic, biochemical and bioelectrical dynamics in pattern regulation.

Publications

Bioelectric gene and reaction networks: computational modelling of genetic, biochemical and bioelectrical dynamics in pattern regulation.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure