Identification of alterations in the Jacobian of biochemical reaction networks from steady state covariance data at two conditions

Model building of biochemical reaction networks typically involves experiments in which changes in the behavior due to natural or experimental perturbations are observed. Computational models of reaction networks are also used in a systems biology approach to study how transitions from a healthy to a diseased state result from changes in genetic or environmental conditions. In this paper we consider the nonlinear inverse problem of inferring information about the Jacobian of a Langevin type network model from covariance data of steady state concentrations associated to two different experimental conditions. Under idealized assumptions on the Langevin fluctuation matrices we prove that relative alterations in the network Jacobian can be uniquely identified when comparing the two data sets. Based on this result and the premise that alteration is locally confined to separable parts due to network modularity we suggest a computational approach using hybrid stochastic-deterministic optimization for the detection of perturbations in the network Jacobian using the sparsity promoting effect of $\ell _p$ -penalization. Our approach is illustrated by means of published metabolomic and signaling reaction networks.     

Nonlinear Dynamics Analysis of a Self-Organizing Recurrent Neural Network: Chaos Waning

Self-organization is thought to play an important role in structuring nervous systems. It frequently arises as a consequence of plasticity mechanisms in neural networks: connectivity determines network dynamics which in turn feed back on network structure through various forms of plasticity. Recently, self-organizing recurrent neural network models (SORNs) have been shown to learn non-trivial structure in their inputs and to reproduce the experimentally observed statistics and fluctuations of synaptic connection strengths in cortex and hippocampus. However, the dynamics in these networks and how they change with network evolution are still poorly understood. Here we investigate the degree of chaos in SORNs by studying how the networks'' self-organization changes their response to small perturbations. We study the effect of perturbations to the excitatory-to-excitatory weight matrix on connection strengths and on unit activities. We find that the network dynamics, characterized by an estimate of the maximum Lyapunov exponent, becomes less chaotic during its self-organization, developing into a regime where only few perturbations become amplified. We also find that due to the mixing of discrete and (quasi-)continuous variables in SORNs, small perturbations to the synaptic weights may become amplified only after a substantial delay, a phenomenon we propose to call deferred chaos.     

Kinetic asymmetry as a key source of functional diversity in biochemical networks

From the analysis of the dynamic properties of various symmetric and asymmetric kinetic schemes, the present report demonstrates that all kinetic schemes, which can be hypothetically divided into two equal halves about an axis of mirror symmetry, are endowed with structural metastability under mass-closed conditions. In mass-closed symmetric schemes, absolute symmetry in reaction conditions in two halves is essential for the occurrence of ordered dynamic behaviour. Even an infinitesimal deviation from the symmetry relations instantaneously drives such systems from limit-cycles to turbulence. Reaction schemes with no axes of symmetry may exhibit a large variety of complex, structurally stable temporal order for wide ranges of values of system parameters and variables. Kinetic asymmetry, therefore, may confer to biochemical networks the functional diversity as well as stability against environmental perturbations.     

Coherent coupling of feedback loops: a design principle of cell signaling networks

MOTIVATION: It is widely accepted that cell signaling networks have been evolved to be robust against perturbations. To investigate the topological characteristics resulting in such robustness, we have examined large-scale signaling networks and found that a number of feedback loops are present mostly in coupled structures. In particular, the coupling was made in a coherent way implying that same types of feedback loops are interlinked together. RESULTS: We have investigated the role of such coherently coupled feedback loops through extensive Boolean network simulations and found that a high proportion of coherent couplings can enhance the robustness of a network against its state perturbations. Moreover, we found that the robustness achieved by coherently coupled feedback loops can be kept evolutionarily stable. All these results imply that the coherent coupling of feedback loops might be a design principle of cell signaling networks devised to achieve the robustness.     

Characterizing multistationarity regimes in biochemical reaction networks

Switch like responses appear as common strategies in the regulation of cellular systems. Here we present a method to characterize bistable regimes in biochemical reaction networks that can be of use to both direct and reverse engineering of biological switches. In the design of a synthetic biological switch, it is important to study the capability for bistability of the underlying biochemical network structure. Chemical Reaction Network Theory (CRNT) may help at this level to decide whether a given network has the capacity for multiple positive equilibria, based on their structural properties. However, in order to build a working switch, we also need to ensure that the bistability property is robust, by studying the conditions leading to the existence of two different steady states. In the reverse engineering of biological switches, knowledge collected about the bistable regimes of the underlying potential model structures can contribute at the model identification stage to a drastic reduction of the feasible region in the parameter space of search. In this work, we make use and extend previous results of the CRNT, aiming not only to discriminate whether a biochemical reaction network can exhibit multiple steady states, but also to determine the regions within the whole space of parameters capable of producing multistationarity. To that purpose we present and justify a condition on the parameters of biochemical networks for the appearance of multistationarity, and propose an efficient and reliable computational method to check its satisfaction through the parameter space.     

Bistability and oscillations in chemical reaction networks

Bifurcation theory is one of the most widely used approaches for analysis of dynamical behaviour of chemical and biochemical reaction networks. Some of the interesting qualitative behaviour that are analyzed are oscillations and bistability (a situation where a system has at least two coexisting stable equilibria). Both phenomena have been identified as central features of many biological and biochemical systems. This paper, using the theory of stoichiometric network analysis (SNA) and notions from algebraic geometry, presents sufficient conditions for a reaction network to display bifurcations associated with these phenomena. The advantage of these conditions is that they impose fewer algebraic conditions on model parameters than conditions associated with standard bifurcation theorems. To derive the new conditions, a coordinate transformation will be made that will guarantee the existence of branches of positive equilibria in the system. This is particularly useful in mathematical biology, where only positive variable values are considered to be meaningful. The first part of the paper will be an extended introduction to SNA and algebraic geometry-related methods which are used in the coordinate transformation and set up of the theorems. In the second part of the paper we will focus on the derivation of bifurcation conditions using SNA and algebraic geometry. Conditions will be derived for three bifurcations: the saddle-node bifurcation, a simple branching point, both linked to bistability, and a simple Hopf bifurcation. The latter is linked to oscillatory behaviour. The conditions derived are sufficient and they extend earlier results from stoichiometric network analysis as can be found in (Aguda and Clarke in J Chem Phys 87:3461–3470, 1987; Clarke and Jiang in J Chem Phys 99:4464–4476, 1993; Gatermann et al. in J Symb Comput 40:1361–1382, 2005). In these papers some necessary conditions for two of these bifurcations were given. A set of examples will illustrate that algebraic conditions arising from given sufficient bifurcation conditions are not more difficult to interpret nor harder to calculate than those arising from necessary bifurcation conditions. Hence an increasing amount of information is gained at no extra computational cost. The theory can also be used in a second step for a systematic bifurcation analysis of larger reaction networks. We have added a dedication of the paper to K. Gatermann.     

The phases of small networks of chemical reactors and neurons

We present an experimental study of the phase relationships observed in small reactor networks consisting of two and three continuous flow stirred tank reactors. In the three-reactor network one chemical oscillator is coupled to two other reactors in parallel in analogy to a small neural net. Each reactor contains an identical reaction mixture of the excitable Belousov-Zhabotinsky reaction which is characterized by its bifurcation diagram, where the electrical current is the bifurcation parameter. Coupling between the reactors is electrical via Pt-working electrodes and it can be either repulsive (inhibitory) or attractive (excitatory). An external electrical stimulus is applied to all three reactors in the form of an asymmetric electrical current pulse which sweeps across the bifurcation diagram. As a consequence, all three reactors oscillate with characteristic oscillation patterns or remain silent in analogy to the firing of neurons. The observed phase behavior depends on the type of coupling in a complex way. This situation is analogous to the in vivo measurements on single neurons (local neurons and projection neurons) performed by G. Laurent and co-workers on the olfactory system of the locust. We propose a simple neural network similar to the reactor network using the Hodgkin-Huxley model to simulate the action potentials of the coupled single neurons. Analogies between the reactor network and the neural network are discussed.     

Computational multiple steady states for enzymatic esterification of ethanol and oleic acid in an isothermal CSTR

The capacity of complex biochemical reaction networks (consisting of 11 coupled non-linear ordinary differential equations) to show multiple steady states, was investigated. The system involved esterification of ethanol and oleic acid by lipase in an isothermal continuous stirred tank reactor (CSTR). The Deficiency One Algorithm and the Subnetwork Analysis were applied to determine the steady state multiplicity. A set of rate constants and two corresponding steady states are computed. The phenomena of bistability, hysteresis and bifurcation are discussed. Moreover, the capacity of steady state multiplicity is extended to the family of the studied reaction networks.     

A new measure of the robustness of biochemical networks

MOTIVATION: The robustness of a biochemical network is defined as the tolerance of variations in kinetic parameters with respect to the maintenance of steady state. Robustness also plays an important role in the fail-safe mechanism in the evolutionary process of biochemical networks. The purposes of this paper are to use the synergism and saturation system (S-system) representation to describe a biochemical network and to develop a robustness measure of a biochemical network subject to variations in kinetic parameters. Since most biochemical networks in nature operate close to the steady state, we consider only the robustness measurement of a biochemical network at the steady state. RESULTS: We show that the upper bound of the tolerated parameter variations is related to the system matrix of a biochemical network at the steady state. Using this upper bound, we can calculate the tolerance (robustness) of a biochemical network without testing many parametric perturbations. We find that a biochemical network with a large tolerance can also better attenuate the effects of variations in rate parameters and environments. Compensatory parameter variations and network redundancy are found to be important mechanisms for the robustness of biochemical networks. Finally, four biochemical networks, such as a cascaded biochemical network, the glycolytic-glycogenolytic pathway in a perfused rat liver, the tricarboxylic acid cycle in Dictyostelium discoideum and the cAMP oscillation network in bacterial chemotaxis, are used to illustrate the usefulness of the proposed robustness measure.     

Identification of small scale biochemical networks based on general type system perturbations

New technologies enable acquisition of large data-sets containing genomic, proteomic and metabolic information that describe the state of a cell. These data-sets call for systematic methods enabling relevant information about the inner workings of the cell to be extracted. One important issue at hand is the understanding of the functional interactions between genes, proteins and metabolites. We here present a method for identifying the dynamic interactions between biochemical components within the cell, in the vicinity of a steady-state. Key features of the proposed method are that it can deal with data obtained under perturbations of any system parameter, not only concentrations of specific components, and that the direct effect of the perturbations does not need to be known. This is important as concentration perturbations are often difficult to perform in biochemical systems and the specific effects of general type perturbations are usually highly uncertain, or unknown. The basis of the method is a linear least-squares estimation, using time-series measurements of concentrations and expression profiles, in which system states and parameter perturbations are estimated simultaneously. An important side-effect of also employing estimation of the parameter perturbations is that knowledge of the system's steady-state concentrations, or activities, is not required and that deviations from steady-state prior to the perturbation can be dealt with. Time derivatives are computed using a zero-order hold discretization, shown to yield significant improvements over the widely used Euler approximation. We also show how network interactions with dynamics that are too fast to be captured within the available sampling time can be determined and excluded from the network identification. Known and unknown moiety conservation relationships can be processed in the same manner. The method requires that the number of samples equals at least the number of network components and, hence, is at present restricted to relatively small-scale networks. We demonstrate herein the performance of the method on two small-scale in silico genetic networks.     

Systems biology and the origins of life? Part I. Are biochemical networks possible ancestors of living systems? Reproduction,identity and sensitivity to signals of biochemical networks

The set of these two theoretical papers offers an alternative to the hypothesis of a primordial RNA-world. The basic idea of these papers is to consider that the first prebiotic systems could have been networks of catalysed reactions encapsulated by a membrane. In order to test this hypothesis it was attempted to list the main obligatory features of living systems and see whether encapsulated biochemical networks could possibly display these features. The traits of living systems are the following: the ability they have to reproduce; the fact they possess an identity; the fact that biological events should be considered in the context of a history; the fact that living systems are able to evolve by selection of alterations of their structure and self-organization. The aim of these two papers is precisely to show that encapsulated biochemical networks can possess these properties and can be considered good candidates for the first prebiotic systems. In the present paper it is shown that if the proteinoids are not very specific catalysts and if some of the reactions of the network are autocatalytic whereas others are not, the resulting system does not reach a steady-state and tends to duplicate. In the same line, these biochemical networks possess an identity, viz. an information, defined from the probability of occurrence of these nodes. Moreover interaction of two ligands can increase, or decrease, this information. In the first case, the system is defined as emergent, in the second case it is considered integrated. Another property of living systems is that their behaviour is defined in the context of a time-arrow. For instance, they are able to sense whether the intensity of a signal is reached after an increase, or a decrease. This property can be mimicked by a simple physico-chemical system made up of the diffusion of a ligand followed by its chemical transformation catalysed by a proteinoid displaying inhibition by excess substrate. Under these conditions the system reacts differently depending on whether the same ligand concentration is reached after an increase or a decrease.     

Robust synchronization analysis in nonlinear stochastic cellular networks with time-varying delays, intracellular perturbations and intercellular noise

Naturally, a cellular network consisted of a large amount of interacting cells is complex. These cells have to be synchronized in order to emerge their phenomena for some biological purposes. However, the inherently stochastic intra and intercellular interactions are noisy and delayed from biochemical processes. In this study, a robust synchronization scheme is proposed for a nonlinear stochastic time-delay coupled cellular network (TdCCN) in spite of the time-varying process delay and intracellular parameter perturbations. Furthermore, a nonlinear stochastic noise filtering ability is also investigated for this synchronized TdCCN against stochastic intercellular and environmental disturbances. Since it is very difficult to solve a robust synchronization problem with the Hamilton-Jacobi inequality (HJI) matrix, a linear matrix inequality (LMI) is employed to solve this problem via the help of a global linearization method. Through this robust synchronization analysis, we can gain a more systemic insight into not only the robust synchronizability but also the noise filtering ability of TdCCN under time-varying process delays, intracellular perturbations and intercellular disturbances. The measures of robustness and noise filtering ability of a synchronized TdCCN have potential application to the designs of neuron transmitters, on-time mass production of biochemical molecules, and synthetic biology. Finally, a benchmark of robust synchronization design in Escherichia coli repressilators is given to confirm the effectiveness of the proposed methods.     

Growing actin networks form lamellipodium and lamellum by self-assembly

Many different cell types are able to migrate by formation of a thin actin-based cytoskeletal extension. Recently, it became evident that this extension consists of two distinct substructures, designated lamellipodium and lamellum, which differ significantly in their kinetic and kinematic properties as well as their biochemical composition. We developed a stochastic two-dimensional computer simulation that includes chemical reaction kinetics, G-actin diffusion, and filament transport to investigate the formation of growing actin networks in migrating cells. Model parameters were chosen based on experimental data or theoretical considerations. In this work, we demonstrate the system's ability to form two distinct networks by self-organization. We found a characteristic transition in mean filament length as well as a distinct maximum in depolymerization flux, both within the first 1-2 μm. The separation into two distinct substructures was found to be extremely robust with respect to initial conditions and variation of model parameters. We quantitatively investigated the complex interplay between ADF/cofilin and tropomyosin and propose a plausible mechanism that leads to spatial separation of, respectively, ADF/cofilin- or tropomyosin-dominated compartments. Tropomyosin was found to play an important role in stabilizing the lamellar actin network. Furthermore, the influence of filament severing and annealing on the network properties is explored, and simulation data are compared to existing experimental data.     

Protein fragment complementation strategies for biochemical network mapping

The organization of biochemical networks that make up the living cell can be defined by studying the dynamics of protein-protein interactions. To this end, experimental strategies based on protein fragment complementation assays (PCAs) have been used to map biochemical networks and to identify novel components of these networks. Pharmacological perturbations of the interactions can be observed, and the resulting pharmacological profiles and subcellular locations of interactions allow each gene product to be 'placed' at its relevant point in a network. Network mapping by PCA could be used with, or instead of, traditional target-based drug discovery strategies to increase the quantity and quality of information about the actions of small molecules on living cells and the intricate networks that make up their chemical machinery.     

On the estimation of robustness and filtering ability of dynamic biochemical networks under process delays, internal parametric perturbations and external disturbances

Inherently, biochemical regulatory networks suffer from process delays, internal parametrical perturbations as well as external disturbances. Robustness is the property to maintain the functions of intracellular biochemical regulatory networks despite these perturbations. In this study, system and signal processing theories are employed for measurement of robust stability and filtering ability of linear and nonlinear time-delay biochemical regulatory networks. First, based on Lyapunov stability theory, the robust stability of biochemical network is measured for the tolerance of additional process delays and additive internal parameter fluctuations. Then the filtering ability of attenuating additive external disturbances is estimated for time-delay biochemical regulatory networks. In order to overcome the difficulty of solving the Hamilton Jacobi inequality (HJI), the global linearization technique is employed to simplify the measurement procedure by a simple linear matrix inequality (LMI) method. Finally, an example is given in silico to illustrate how to measure the robust stability and filtering ability of a nonlinear time-delay perturbative biochemical network. This robust stability and filtering ability measurement for biochemical network has potential application to synthetic biology, gene therapy and drug design.     

The thermodynamic meaning of metabolic exchange fluxes

Metabolic flux analysis (MFA) deals with the experimental determination of steady-state fluxes in metabolic networks. An important feature of the ¹³C MFA method is its capability to generate information on both directions of bidirectional reaction steps given by exchange fluxes. The biological interpretation of these exchange fluxes and their relation to thermodynamic properties of the respective reaction steps has never been systematically investigated. As a central result, it is shown here that for a general class of enzyme reaction mechanisms the quotients of net and exchange fluxes measured by ¹³C MFA are coupled to Gibbs energies of the reaction steps. To establish this relation the concept of apparent flux ratios of enzymatic isotope-labeling networks is introduced and some computing rules for these flux ratios are given. Application of these rules reveals a conceptional pitfall of ¹³C MFA, which is the inherent dependency of measured exchange fluxes on the chosen tracer atom. However, it is shown that this effect can be neglected for typical biochemical reaction steps under physiological conditions. In this situation, the central result can be formulated as a two-sided inequality relating fluxes, pool sizes, and standard Gibbs energies. This relation has far-reaching consequences for metabolic flux analysis, quantitative metabolomics, and network thermodynamics.     

Edgetic perturbation models of human inherited disorders

Cellular functions are mediated through complex systems of macromolecules and metabolites linked through biochemical and physical interactions, represented in interactome models as ‘nodes’ and ‘edges’, respectively. Better understanding of genotype‐to‐phenotype relationships in human disease will require modeling of how disease‐causing mutations affect systems or interactome properties. Here we investigate how perturbations of interactome networks may differ between complete loss of gene products (‘node removal’) and interaction‐specific or edge‐specific (‘edgetic’) alterations. Global computational analyses of ～50 000 known causative mutations in human Mendelian disorders revealed clear separations of mutations probably corresponding to those of node removal versus edgetic perturbations. Experimental characterization of mutant alleles in various disorders identified diverse edgetic interaction profiles of mutant proteins, which correlated with distinct structural properties of disease proteins and disease mechanisms. Edgetic perturbations seem to confer distinct functional consequences from node removal because a large fraction of cases in which a single gene is linked to multiple disorders can be modeled by distinguishing edgetic network perturbations. Edgetic network perturbation models might improve both the understanding of dissemination of disease alleles in human populations and the development of molecular therapeutic strategies.     

Self-Organization in Living Cells: Networks of Protein Machines and Nonequilibrium Soft Matter

Microscopic self-organization phenomena inside a living cell should not represent merely a reduced copy of self-organization in macroscopic systems. A cell is populated by active protein machines that communicate via small molecules diffusing through the cytoplasm. Mutual synchronization of machine cycles can spontaneously develop in such networks – an effect which is similar to coherent laser generation. On the other hand, an interplay between reactions, diffusion and phase transitions in biological soft matter may lead to the formation of stationary or traveling nonequilibrium nanoscale structures.