Causal structure of oscillations in gene regulatory networks: Boolean analysis of ordinary differential equation attractors

On the dynamical aspects of local translation at the activated synapse

Background

The key role in the dynamic regulation of synaptic protein turnover belongs to the Fragile X Mental Retardation Protein, which regulates the efficiency of dendritic mRNA translation in response to stimulation of metabotropic glutamate receptors at excitatory synapses of the hippocampal pyramidal cells. Its activity is regulated via positive and negative regulatory loops that function in different time ranges, which is an absolute factor for the formation of chaotic regimes that lead to disrupted proteome stability. The indicated condition may cause a number of neuropsychiatric diseases, including autism and epilepsy. The present study is devoted to a theoretical analysis of the local translation system dynamic properties and identification of parameters affecting the chaotic potential of the system.

Results

A mathematical model that describes the maintenance of a specific pool of active receptors on the postsynaptic membrane via two mechanisms – de novo synthesis of receptor proteins and restoration of protein function during the recycling process – has been developed. Analysis of the model revealed that an increase in the values of the parameters describing the impact of protein recycling on the maintenance of a pool of active receptors in the membrane, duration of the signal transduction via the mammalian target of rapamycin pathway, influence of receptors on the translation activation, as well as reduction of the rate of synthesis and integration of de novo synthesized proteins into the postsynaptic membrane – contribute to the reduced complexity of the local translation system dynamic state. Formation of these patterns significantly depends on the complexity and non-linearity of the mechanisms of exposure of de novo synthesized receptors to the postsynaptic membrane, the correct evaluation of which is currently problematic.

Conclusions

The model predicts that an increase of “receptor recycling” and reduction of the rate of synthesis and integration of de novo synthesized proteins into the postsynaptic membrane contribute to the reduced complexity of the local translation system dynamic state. Herewith, stable stationary states occur much less frequently than cyclic states. It is possible that cyclical nature of functioning of the local translation system is its “normal” dynamic state.

Limit cycles in models of circular gene networks regulated by negative feedback loops

BACKGROUND

The regulatory feedback loops that present in structural and functional organization of molecular-genetic systems and the phenomenon of the regulatory signal delay, a time period between the moment of signal reception and its implementation, provide natural conditions for complicated dynamic regimes in these systems. The delay phenomenon at the intracellular level is a consequence of the matrix principle of data transmission, implemented through the rather complex processes of transcription and translation.However, the rules of the influence of system structure on system dynamics are not clearly understood. Knowledge of these rules is particularly important for construction of synthetic gene networks with predetermined properties.

RESULTS

We study dynamical properties of models of simplest circular gene networks regulated by negative feedback mechanisms. We have shown existence and stability of oscillating trajectories (cycles) in these models. Two algorithms of construction and localization of these cycles have been proposed. For one of these models, we have solved an inverse problem of parameters identification.

CONCLUSIONS

The modeling results demonstrate that non-stationary dynamics in the models of circular gene networks with negative feedback loops is achieved by a high degree of non-linearity of the mechanism of the autorepressor influence on its own expression, by the presence of regulatory signal delay, the value of which must exceed a certain critical value, and transcription/translation should be initiated from a sufficiently strong promoter/Shine-Dalgarno site. We believe that the identified patterns are key elements of the oscillating construction design.

Ensemble Method of Feature Selection and Reverse Construction of Gene Logical Network Based on Information Entropy

In this paper, we propose a novel ensemble gene selection method to obtain a gene subset. Then we provide a reverse construction method of gene network derived from expression profile data of the gene subset. The uncertainty coefficient based on information entropy are used to define the existence of logical relations among these genes. If the uncertainty coefficient between some genes exceeds predefined thresholds, the gene nodes will be connected by directed edges. Thus, a gene network is generated, which we define as gene logical network. This method is applied to the breast cancer data including control group and experimental group, with comparisons of the 2nd-order logic type distribution, average degree as well as average path length of the networks. It is found that these structures with different networks are quite distinct. By the comparison of the degree difference between control group and experimental group, the key genes are picked up. By defining the dynamics evolution rules of state transition based on the logical regulation among the key genes in the network, the dynamic behaviors for normal breast cells and cells with cancer of different stages are simulated numerically. Some of them are highly related to the development of breast cancer through literature inquiry. The study may provide a useful revelation to the biological mechanism in the formation and development of cancer.

Super-transient scaling in time-delay autonomous Boolean network motifs

Autonomous Boolean networks are commonly used to model the dynamics of gene regulatory networks and allow for the prediction of stable dynamical attractors. However, most models do not account for time delays along the network links and noise, which are crucial features of real biological systems. Concentrating on two paradigmatic motifs, the toggle switch and the repressilator, we develop an experimental testbed that explicitly includes both inter-node time delays and noise using digital logic elements on field-programmable gate arrays. We observe transients that last millions to billions of characteristic time scales and scale exponentially with the amount of time delays between nodes, a phenomenon known as super-transient scaling. We develop a hybrid model that includes time delays along network links and allows for stochastic variation in the delays. Using this model, we explain the observed super-transient scaling of both motifs and recreate the experimentally measured transient distributions.

Boolean modeling: a logic-based dynamic approach for understanding signaling and regulatory networks and for making useful predictions

The biomolecules inside or near cells form a complex interacting system. Cellular phenotypes and behaviors arise from the totality of interactions among the components of this system. A fruitful way of modeling interacting biomolecular systems is by network-based dynamic models that characterize each component by a state variable, and describe the change in the state variables due to the interactions in the system. Dynamic models can capture the stable state patterns of this interacting system and can connect them to different cell fates or behaviors. A Boolean or logic model characterizes each biomolecule by a binary state variable that relates the abundance of that molecule to a threshold abundance necessary for downstream processes. The regulation of this state variable is described in a parameter free manner, making Boolean modeling a practical choice for systems whose kinetic parameters have not been determined. Boolean models integrate the body of knowledge regarding the components and interactions of biomolecular systems, and capture the system's dynamic repertoire, for example the existence of multiple cell fates. These models were used for a variety of systems and led to important insights and predictions. Boolean models serve as an efficient exploratory model, a guide for follow-up experiments, and as a foundation for more quantitative models.

Regulatory logic and pattern formation in the early sea urchin embryo

We model the endomesoderm tissue specification process in the vegetal half of the early sea urchin embryo using Boolean models with continuous-time updating to represent the regulatory network that controls gene expression. Our models assume that the network interaction rules remain constant over time and the dynamics plays out on a predetermined program of cell divisions. An exhaustive search of two-node models, in which each node may represent a module of several genes in the real regulatory network, yields a unique network architecture that can accomplish the pattern formation task at hand--the formation of three latitudinal tissue bands from an initial state with only two distinct cell types. Analysis of an eight-gene model constructed from available experimental data reveals that it has a modular structure equivalent to the successful two-node case. Our results support the hypothesis that the gene regulatory network provides sufficient instructions for producing the correct pattern of tissue specification at this stage of development (between the fourth and tenth cleavages in the urchin embryo).

Piecewise linear and Boolean models of chemical reaction networks

Models of biochemical networks are frequently complex and high-dimensional. Reduction methods that preserve important dynamical properties are therefore essential for their study. Interactions in biochemical networks are frequently modeled using Hill functions ([Formula: see text]). Reduced ODEs and Boolean approximations of such model networks have been studied extensively when the exponent [Formula: see text] is large. However, while the case of small constant [Formula: see text] appears in practice, it is not well understood. We provide a mathematical analysis of this limit and show that a reduction to a set of piecewise linear ODEs and Boolean networks can be mathematically justified. The piecewise linear systems have closed-form solutions that closely track those of the fully nonlinear model. The simpler, Boolean network can be used to study the qualitative behavior of the original system. We justify the reduction using geometric singular perturbation theory and compact convergence, and illustrate the results in network models of a toggle switch and an oscillator.

Controllability of asynchronous Boolean multiplex control networks

In this article, the controllability of asynchronous Boolean multiplex control networks (ABMCNs) is studied. First, the model of Boolean multiplex control networks under Harvey' asynchronous update is presented. By means of semi-tensor product approach, the logical dynamics is converted into linear representation, and a generalized formula of control-depending network transition matrices is achieved. Second, a necessary and sufficient condition is proposed to verify that only control-depending fixed points of ABMCNs can be controlled with probability one. Third, using two types of controls, the controllability of system is studied and formulae are given to show: (a) when an initial state is given, the reachable set at time s under a group of specified controls; (b) the reachable set at time s under arbitrary controls; (c) the specific probability values from a given initial state to destination states. Based on the above formulae, an algorithm to calculate overall reachable states from a specified initial state is presented. Moreover, we also discuss an approach to find the particular control sequence which steers the system between two states with maximum probability. Examples are shown to illustrate the feasibility of the proposed scheme.

FOUNDATIONS FOR MODELING THE DYNAMICS OF GENE REGULATORY NETWORKS: A MULTILEVEL-PERSPECTIVE REVIEW

A promising alternative for unraveling the principles under which the dynamic interactions among genes lead to cellular phenotypes relies on mathematical and computational models at different levels of abstraction, from the molecular level of protein-DNA interactions to the system level of functional relationships among genes. This review article presents, under a bottom–up perspective, a hierarchy of approaches to modeling gene regulatory network dynamics, from microscopic descriptions at the single-molecule level in the spatial context of an individual cell to macroscopic models providing phenomenological descriptions at the population-average level. The reviewed modeling approaches include Molecular Dynamics, Particle-Based Brownian Dynamics, the Master Equation approach, Ordinary Differential Equations, and the Boolean logic abstraction. Each of these frameworks is motivated by a particular biological context and the nature of the insight being pursued. The setting of gene network dynamic models from such frameworks involves assumptions and mathematical artifacts often ignored by the non-specialist. This article aims at providing an entry point for biologists new to the field and computer scientists not acquainted with some recent biophysically-inspired models of gene regulation. The connections promoting intuition between different abstraction levels and the role that approximations play in the modeling process are highlighted throughout the paper.

Introduction to focus issue: quantitative approaches to genetic networks

All cells of living organisms contain similar genetic instructions encoded in the organism's DNA. In any particular cell, the control of the expression of each different gene is regulated, in part, by binding of molecular complexes to specific regions of the DNA. The molecular complexes are composed of protein molecules, called transcription factors, combined with various other molecules such as hormones and drugs. Since transcription factors are coded by genes, cellular function is partially determined by genetic networks. Recent research is making large strides to understand both the structure and the function of these networks. Further, the emerging discipline of synthetic biology is engineering novel gene circuits with specific dynamic properties to advance both basic science and potential practical applications. Although there is not yet a universally accepted mathematical framework for studying the properties of genetic networks, the strong analogies between the activation and inhibition of gene expression and electric circuits suggest frameworks based on logical switching circuits. This focus issue provides a selection of papers reflecting current research directions in the quantitative analysis of genetic networks. The work extends from molecular models for the binding of proteins, to realistic detailed models of cellular metabolism. Between these extremes are simplified models in which genetic dynamics are modeled using classical methods of systems engineering, Boolean switching networks, differential equations that are continuous analogues of Boolean switching networks, and differential equations in which control is based on power law functions. The mathematical techniques are applied to study: (i) naturally occurring gene networks in living organisms including: cyanobacteria, Mycoplasma genitalium, fruit flies, immune cells in mammals; (ii) synthetic gene circuits in Escherichia coli and yeast; and (iii) electronic circuits modeling genetic networks using field-programmable gate arrays. Mathematical analyses will be essential for understanding naturally occurring genetic networks in diverse organisms and for providing a foundation for the improved development of synthetic genetic networks.