A discrete model with density dependent fast migration

The aim of this work is to develop an approximate aggregation method for certain non-linear discrete models. Approximate aggregation consists in describing the dynamics of a general system involving many coupled variables by means of the dynamics of a reduced system with a few global variables. We present discrete models with two different time scales, the slow one considered to be linear and the fast one non-linear because of its transition matrix depends on the global variables. In our discrete model the time unit is chosen to be the one associated to the slow dynamics, and then we approximate the effect of fast dynamics by using a sufficiently large power of its corresponding transition matrix. In a previous work the same system is treated in the case of fast dynamics considered to be linear, conservative in the global variables and inducing a stable frequency distribution of the state variables. A similar non-linear model has also been studied which uses as time unit the one associated to the fast dynamics and has the non-linearity in the slow part of the system. In the present work we transform the system to make the global variables explicit, and we justify the quick derivation of the aggregated system. The local asymptotic behaviour of the aggregated system entails that of the general system under certain conditions, for instance, if the aggregated system has a stable hyperbolic fixed point then the general system has one too. The method is applied to aggregate a multiregional Leslie model with density dependent migration rates.     

The Validity of Quasi-Steady-State Approximations in Discrete Stochastic Simulations

In biochemical networks, reactions often occur on disparate timescales and can be characterized as either fast or slow. The quasi-steady-state approximation (QSSA) utilizes timescale separation to project models of biochemical networks onto lower-dimensional slow manifolds. As a result, fast elementary reactions are not modeled explicitly, and their effect is captured by nonelementary reaction-rate functions (e.g., Hill functions). The accuracy of the QSSA applied to deterministic systems depends on how well timescales are separated. Recently, it has been proposed to use the nonelementary rate functions obtained via the deterministic QSSA to define propensity functions in stochastic simulations of biochemical networks. In this approach, termed the stochastic QSSA, fast reactions that are part of nonelementary reactions are not simulated, greatly reducing computation time. However, it is unclear when the stochastic QSSA provides an accurate approximation of the original stochastic simulation. We show that, unlike the deterministic QSSA, the validity of the stochastic QSSA does not follow from timescale separation alone, but also depends on the sensitivity of the nonelementary reaction rate functions to changes in the slow species. The stochastic QSSA becomes more accurate when this sensitivity is small. Different types of QSSAs result in nonelementary functions with different sensitivities, and the total QSSA results in less sensitive functions than the standard or the prefactor QSSA. We prove that, as a result, the stochastic QSSA becomes more accurate when nonelementary reaction functions are obtained using the total QSSA. Our work provides an apparently novel condition for the validity of the QSSA in stochastic simulations of biochemical reaction networks with disparate timescales.     

The Validity of Quasi-Steady-State Approximations in Discrete Stochastic Simulations

In biochemical networks, reactions often occur on disparate timescales and can be characterized as either fast or slow. The quasi-steady-state approximation (QSSA) utilizes timescale separation to project models of biochemical networks onto lower-dimensional slow manifolds. As a result, fast elementary reactions are not modeled explicitly, and their effect is captured by nonelementary reaction-rate functions (e.g., Hill functions). The accuracy of the QSSA applied to deterministic systems depends on how well timescales are separated. Recently, it has been proposed to use the nonelementary rate functions obtained via the deterministic QSSA to define propensity functions in stochastic simulations of biochemical networks. In this approach, termed the stochastic QSSA, fast reactions that are part of nonelementary reactions are not simulated, greatly reducing computation time. However, it is unclear when the stochastic QSSA provides an accurate approximation of the original stochastic simulation. We show that, unlike the deterministic QSSA, the validity of the stochastic QSSA does not follow from timescale separation alone, but also depends on the sensitivity of the nonelementary reaction rate functions to changes in the slow species. The stochastic QSSA becomes more accurate when this sensitivity is small. Different types of QSSAs result in nonelementary functions with different sensitivities, and the total QSSA results in less sensitive functions than the standard or the prefactor QSSA. We prove that, as a result, the stochastic QSSA becomes more accurate when nonelementary reaction functions are obtained using the total QSSA. Our work provides an apparently novel condition for the validity of the QSSA in stochastic simulations of biochemical reaction networks with disparate timescales.     

Application of Quasi-Steady-State Methods to Nonlinear Models of Intracellular Transport by Molecular Motors

Molecular motors such as kinesin and dynein are responsible for transporting material along microtubule networks in cells. In many contexts, motor dynamics can be modelled by a system of reaction–advection–diffusion partial differential equations (PDEs). Recently, quasi-steady-state (QSS) methods have been applied to models with linear reactions to approximate the behaviour of the full PDE system. Here, we extend this QSS reduction methodology to certain nonlinear reaction models. The QSS method relies on the assumption that the nonlinear binding and unbinding interactions of the cellular motors occur on a faster timescale than the spatial diffusion and advection processes. The full system dynamics are shown to be well approximated by the dynamics on the slow manifold. The slow manifold is parametrized by a single scalar quantity that satisfies a scalar nonlinear PDE, called the QSS PDE. We apply the QSS method to several specific nonlinear models for the binding and unbinding of molecular motors, and we use the resulting approximations to draw conclusions regarding the parameter dependence of the spatial distribution of motors for these models.     

Time hierarchy, equilibrium and non-equilibrium in metabolic systems.

A metabolic system consists of cooperating biochemical reactions. The motion is described by differential equations in the metabolites. The right-hand sides of these equations are linear combinations of the velocities of the individual reactions. These velocities depend in a non-linear manner on the metabolite concentrations (according to the law of mass action). A characteristic "metabolic" time may be defined for the motion of the whole system. It scales the essential metabolic events whose evolution time is comparable to this metabolite time unit. The constituent reactions of the metabolic system have an individual characteristic time which need not coincide with the general metabolic time. The individual time characterises the approach to the individual equilibrium of the isolated undisturbed reaction. According to the ratio of these two time scales, a single reaction may be fast, or slow, or essential, as compared with the metabolic events. Characteristic time of a single reaction and its steady-state deviation from equilibrium are closely related. It can be shown that the relative deviation from equilibrium of a reaction within the metabolic network is of the same numerical order as the ratio between individual time to metabolic time. The interaction of many reactions with different characteristic times introduces a time hierarchy into the system. This can be made transparent by appropriate scaling and by linear transformation of the system. The subsystem of fast cooperating reactions (dehydrogenases, phosphotransferases) attains a state which is near to the individual equilibrium and reestablishes this state after perturbation. The equilibration is fast; an ultrarapid phase of cofactor equilibrium can be distinguished from the fast phase of substrate equilibrium (exchange of metabolic material between different pathways). During the slower metabolic phase these near-equilibria manifest themselves as stoichiometric linkage between unrelated metabolites. The latter cease to be independent variables and combine to metabolic pools. It can be strictly shown that the essential variables at the metabolic time scale are carrier pools and the degree of occupancy of these carriers by metabolic groups. Chemically different types of carrier pools may be functionally linked together by fast reactions. A consequence of such an arrangement of reactions are distance effects: Changes at one end of a metabolic map may be directly conveyed to other pathways via stoichiometric linkage brought about by fast equilibration of cofactor reactions.     

Role of the autonomic nervous system in generating non-linear dynamics in short-term heart period variability.

We evaluated the role played by the autonomic nervous system in producing non-linear dynamics in short heart period variability (HPV) series recorded in healthy young humans. Non-linear dynamics are detected using an index of predictability based on a local non-linear predictor and a surrogate data approach. Different types of surrogates are utilized: (i) phase-randomized Fourier-transform based (FT) data; (ii) amplitude-adjusted FT (AAFT) data; and (iii) iteratively refined AAFT (IAAFT) data of two types (IAAFT-1 and IAAFT-2). The approach was applied to experimental protocols activating or blocking the sympathetic or parasympathetic branches of the autonomic nervous system or periodically perturbing cardiovascular control via paced respiration at different breathing rates. We found that short-term HPV was mostly linear at rest. Experimental protocols activating the sympathetic or parasympathetic nervous system did not produce non-linear dynamics. In contrast, paced respiration, especially at slow breathing rates, elicited significantly non-linear dynamics. Therefore, in short-term HPV ( approximately 300 beats) the use of non-linear models is not supported by the data, except under conditions whereby the subject is constrained to a slow respiratory rate.     

What Can Causal Networks Tell Us about Metabolic Pathways?

Graphical models describe the linear correlation structure of data and have been used to establish causal relationships among phenotypes in genetic mapping populations. Data are typically collected at a single point in time. Biological processes on the other hand are often non-linear and display time varying dynamics. The extent to which graphical models can recapitulate the architecture of an underlying biological processes is not well understood. We consider metabolic networks with known stoichiometry to address the fundamental question: “What can causal networks tell us about metabolic pathways?”. Using data from an Arabidopsis BaySha population and simulated data from dynamic models of pathway motifs, we assess our ability to reconstruct metabolic pathways using graphical models. Our results highlight the necessity of non-genetic residual biological variation for reliable inference. Recovery of the ordering within a pathway is possible, but should not be expected. Causal inference is sensitive to subtle patterns in the correlation structure that may be driven by a variety of factors, which may not emphasize the substrate-product relationship. We illustrate the effects of metabolic pathway architecture, epistasis and stochastic variation on correlation structure and graphical model-derived networks. We conclude that graphical models should be interpreted cautiously, especially if the implied causal relationships are to be used in the design of intervention strategies.     

Reduction of slow-fast discrete models coupling migration and demography

This work deals with a general class of two-time scales discrete nonlinear dynamical systems which are susceptible of being studied by means of a reduced system that is obtained using the so-called aggregation of variables method. This reduction process is applied to several models of population dynamics driven by demographic and migratory processes which take place at two different time scales: slow and fast. An analysis of these models exchanging the role of the slow and fast dynamics is provided: when a Leslie type demography is faster than migrations, a multi-attractor scenario appears for the reduced dynamics; on the other hand, when the migratory process is faster than demography, the reduction process gives rise to new interpretations of well known discrete models, including some Allee effect scenarios.     

Simplified modelling of metabolic pathways for flux prediction and optimization: lessons from an in vitro reconstruction of the upper part of glycolysis

Explicit modelling of metabolic networks relies on well-known mathematical tools and specialized computer programs. However, identifying and estimating the values of the very numerous enzyme parameters inherent to the models remain a tedious and difficult task, and the rate equations of the reactions are usually not known in sufficient detail. A way to circumvent this problem is to use 'non-mechanistic' models, which may account for the behaviour of the systems with a limited number of parameters. Working on the first part of glycolysis reconstituted in vitro, we showed how to derive, from titration experiments, values of effective enzyme activity parameters that do not include explicitly any of the classical kinetic constants. With a maximum of only two parameters per enzyme, this approach produced very good estimates for the flux values, and enabled us to determine the optimization conditions of the system, i.e. to calculate the set of enzyme concentrations that maximizes the flux. This fast and easy method should be valuable in the context of integrative biology or for metabolic engineering, where the challenge is to deal with the dramatic increase in the number of parameters when the systems become complex.     

Non-local concepts and models in biology

In this paper, we consider local and non-local spatially explicit mathematical models for biological phenomena. We show that, when rate differences between fast and slow local dynamics are great enough, non-local models are adequate simplifications of local models. Non-local models thus avoid describing fast processes in mechanistic detail, instead describing the effects of fast processes on slower ones. As a consequence, non-local models are helpful to biologists because they describe biological systems on scales that are convenient to observation, data collection, and insight. We illustrate these arguments by comparing local and non-local models for the aggregation of hypothetical organisms, and we support theoretical ideas with concrete examples from cell biology and animal behavior.     

Neural rate equations for bursting dynamics derived from conductance-based equations

A method of obtaining rate equations from conductance-based equations is developed and applied to fast-spiking and bursting neocortical neurons. It involves splitting systems of conductance-based equations into fast and slow subsystems, and averaging the effects of fast terms that drive the slowly varying quantities by showing that their average is closely proportional to the firing rate. The dependence of the firing rate on the injected current is then approximated in the analysis. The resulting behavior of the slow variables is then substituted back into the fast equations, with the further approximation of replacing the fast voltages in these terms by effective values. For bursting neurons the method yields two coupled limit-cycle oscillators: a self-exciting oscillator for the slow variables that commences limit-cycle oscillations at a critical current and modulates a fast spike-generating oscillator, thereby leading to slowly modulated bursts with a group of spikes in each burst. The dynamics of these coupled oscillators are then verified against those of the conductance-based equations. Finally, it is shown how to place the results in a form suitable for use in mean-field equations for neural population dynamics.     

Global Self-Regulation of the Cellular Metabolic Structure

Background

Different studies have shown that cellular enzymatic activities are able to self-organize spontaneously, forming a metabolic core of reactive processes that remain active under different growth conditions while the rest of the molecular catalytic reactions exhibit structural plasticity. This global cellular metabolic structure appears to be an intrinsic characteristic common to all cellular organisms. Recent work performed with dissipative metabolic networks has shown that the fundamental element for the spontaneous emergence of this global self-organized enzymatic structure could be the number of catalytic elements in the metabolic networks.

Methodology/Principal Findings

In order to investigate the factors that may affect the catalytic dynamics under a global metabolic structure characterized by the presence of metabolic cores we have studied different transitions in catalytic patterns belonging to a dissipative metabolic network. The data were analyzed using non-linear dynamics tools: power spectra, reconstructed attractors, long-term correlations, maximum Lyapunov exponent and Approximate Entropy; and we have found the emergence of self-regulation phenomena during the transitions in the metabolic activities.

Conclusions/Significance

The analysis has also shown that the chaotic numerical series analyzed correspond to the fractional Brownian motion and they exhibit long-term correlations and low Approximate Entropy indicating a high level of predictability and information during the self-regulation of the metabolic transitions. The results illustrate some aspects of the mechanisms behind the emergence of the metabolic self-regulation processes, which may constitute an important property of the global structure of the cellular metabolism.     

Model reduction of genome-scale metabolic models as a basis for targeted kinetic models

Constraint-based, genome-scale metabolic models are an essential tool to guide metabolic engineering. However, they lack the detail and time dimension that kinetic models with enzyme dynamics offer. Model reduction can be used to bridge the gap between the two methods and allow for the integration of kinetic models into the Design-Built-Test-Learn cycle. Here we show that these reduced size models can be representative of the dynamics of the original model and demonstrate the automated generation and parameterisation of such models. Using these minimal models of metabolism could allow for further exploration of dynamic responses in metabolic networks.     

A Bayesian Approach to the Evolution of Metabolic Networks on a Phylogeny

The availability of genomes of many closely related bacteria with diverse metabolic capabilities offers the possibility of tracing metabolic evolution on a phylogeny relating the genomes to understand the evolutionary processes and constraints that affect the evolution of metabolic networks. Using simple (independent loss/gain of reactions) or complex (incorporating dependencies among reactions) stochastic models of metabolic evolution, it is possible to study how metabolic networks evolve over time. Here, we describe a model that takes the reaction neighborhood into account when modeling metabolic evolution. The model also allows estimation of the strength of the neighborhood effect during the course of evolution. We present Gibbs samplers for sampling networks at the internal node of a phylogeny and for estimating the parameters of evolution over a phylogeny without exploring the whole search space by iteratively sampling from the conditional distributions of the internal networks and parameters. The samplers are used to estimate the parameters of evolution of metabolic networks of bacteria in the genus Pseudomonas and to infer the metabolic networks of the ancestral pseudomonads. The results suggest that pathway maps that are conserved across the Pseudomonas phylogeny have a stronger neighborhood structure than those which have a variable distribution of reactions across the phylogeny, and that some Pseudomonas lineages are going through genome reduction resulting in the loss of a number of reactions from their metabolic networks.     

Functional Roles of Slow Enzyme Conformational Changes in Network Dynamics

Extensive studies from different fields reveal that many macromolecules, especially enzymes, show slow transitions among different conformations. This phenomenon is named such things as dynamic disorder, heterogeneity, hysteretic or mnemonic enzymes across these different fields, and has been directly demonstrated by single molecule enzymology and NMR studies recently. We analyzed enzyme slow conformational changes in the context of regulatory networks. A single enzymatic reaction with slow conformational changes can filter upstream network noises, and can either resonantly respond to the system stimulus at certain frequencies or respond adaptively for sustained input signals of the network fluctuations. It thus can serve as a basic functional motif with properties that are normally for larger intermolecular networks in the field of systems biology. We further analyzed examples including enzymes functioning against pH fluctuations, metabolic state change of Artemia embryos, and kinetic insulation of fluctuations in metabolic networks. The study also suggests that hysteretic enzymes may be building blocks of synthetic networks with various properties such as narrow-banded filtering. The work fills the missing gap between studies on enzyme biophysics and network level dynamics, and reveals that the coupling between the two is functionally important; it also suggests that the conformational dynamics of some enzymes may be evolutionally selected.