Modelling of biochemical reaction networks
Gloppen Jørgensen, Arne Gunnar
2011-01-01
This report investigates signalling in reaction kinetic networks. The main topic is signalling between a substance being controlled by another substance and how this can be related to control theory. Different types of so-called natural controllers are compared and certain properties are investigated. Natural controllers are models on how a catalyst enzyme controls, for example the concentration, of a substance. There are sixteen different combinations of signalling between these substanc...
Mean field interaction in biochemical reaction networks
Tembine, Hamidou
2011-09-01
In this paper we establish a relationship between chemical dynamics and mean field game dynamics. We show that chemical reaction networks can be studied using noisy mean field limits. We provide deterministic, noisy and switching mean field limits and illustrate them with numerical examples. © 2011 IEEE.
Modeling stochasticity in biochemical reaction networks
Constantino, P. H.; Vlysidis, M.; Smadbeck, P.; Kaznessis, Y. N.
2016-03-01
Small biomolecular systems are inherently stochastic. Indeed, fluctuations of molecular species are substantial in living organisms and may result in significant variation in cellular phenotypes. The chemical master equation (CME) is the most detailed mathematical model that can describe stochastic behaviors. However, because of its complexity the CME has been solved for only few, very small reaction networks. As a result, the contribution of CME-based approaches to biology has been very limited. In this review we discuss the approach of solving CME by a set of differential equations of probability moments, called moment equations. We present different approaches to produce and to solve these equations, emphasizing the use of factorial moments and the zero information entropy closure scheme. We also provide information on the stability analysis of stochastic systems. Finally, we speculate on the utility of CME-based modeling formalisms, especially in the context of synthetic biology efforts.
Chemical reaction network approaches to Biochemical Systems Theory.
Arceo, Carlene Perpetua P; Jose, Editha C; Marin-Sanguino, Alberto; Mendoza, Eduardo R
2015-11-01
This paper provides a framework to represent a Biochemical Systems Theory (BST) model (in either GMA or S-system form) as a chemical reaction network with power law kinetics. Using this representation, some basic properties and the application of recent results of Chemical Reaction Network Theory regarding steady states of such systems are shown. In particular, Injectivity Theory, including network concordance [36] and the Jacobian Determinant Criterion [43], a "Lifting Theorem" for steady states [26] and the comprehensive results of Müller and Regensburger [31] on complex balanced equilibria are discussed. A partial extension of a recent Emulation Theorem of Cardelli for mass action systems [3] is derived for a subclass of power law kinetic systems. However, it is also shown that the GMA and S-system models of human purine metabolism [10] do not display the reactant-determined kinetics assumed by Müller and Regensburger and hence only a subset of BST models can be handled with their approach. Moreover, since the reaction networks underlying many BST models are not weakly reversible, results for non-complex balanced equilibria are also needed. PMID:26363083
Reduction of dynamical biochemical reactions networks in computational biology
Radulescu, O.; Gorban, A.N.; Zinovyev, A.; Noel, V.
2012-01-01
Biochemical networks are used in computational biology, to model mechanistic details of systems involved in cell signaling, metabolism, and regulation of gene expression. Parametric and structural uncertainty, as well as combinatorial explosion are strong obstacles against analyzing the dynamics of large models of this type. Multiscaleness, an important property of these networks, can be used to get past some of these obstacles. Networks with many well separated time scales, can be reduced to...
Exact results for noise power spectra in linear biochemical reaction networks
Warren, Patrick B.; Tanase-Nicola, Sorin; Wolde, Pieter Rein ten
2005-01-01
We present a simple method for determining the exact noise power spectra in linear chemical reaction networks. We apply the method to networks which are representative of biochemical processes such as gene expression and signal detection. Our results clarify how noise is transmitted by signal detection motifs, and indicate how to coarse-grain networks by the elimination of fast reactions.
Turing-Hopf instability in biochemical reaction networks arising from pairs of subnetworks.
Mincheva, Maya; Roussel, Marc R
2012-11-01
Network conditions for Turing instability in biochemical systems with two biochemical species are well known and involve autocatalysis or self-activation. On the other hand general network conditions for potential Turing instabilities in large biochemical reaction networks are not well developed. A biochemical reaction network with any number of species where only one species moves is represented by a simple digraph and is modeled by a reaction-diffusion system with non-mass action kinetics. A graph-theoretic condition for potential Turing-Hopf instability that arises when a spatially homogeneous equilibrium loses its stability via a single pair of complex eigenvalues is obtained. This novel graph-theoretic condition is closely related to the negative cycle condition for oscillations in ordinary differential equation models and its generalizations, and requires the existence of a pair of subnetworks, each containing an even number of positive cycles. The technique is illustrated with a double-cycle Goodwin type model. PMID:22698892
Furusawa, Chikara; Kaneko, Kunihiko
2005-01-01
The evolutionary origin of universal statistics in biochemical reaction network is studied, to explain the power-law distribution of reaction links and the power-law distributions of chemical abundances. Using cell models with catalytic reaction network, we find evidence that the power-law distribution in abundances of chemicals emerges by the selection of cells with higher growth speeds. Through the further evolution, this inhomogeneity in chemical abundances is shown to be embedded in the d...
Marchetti, Luca; Priami, Corrado; Thanh, Vo Hong
2016-07-01
This paper introduces HRSSA (Hybrid Rejection-based Stochastic Simulation Algorithm), a new efficient hybrid stochastic simulation algorithm for spatially homogeneous biochemical reaction networks. HRSSA is built on top of RSSA, an exact stochastic simulation algorithm which relies on propensity bounds to select next reaction firings and to reduce the average number of reaction propensity updates needed during the simulation. HRSSA exploits the computational advantage of propensity bounds to manage time-varying transition propensities and to apply dynamic partitioning of reactions, which constitute the two most significant bottlenecks of hybrid simulation. A comprehensive set of simulation benchmarks is provided for evaluating performance and accuracy of HRSSA against other state of the art algorithms.
A moment-convergence method for stochastic analysis of biochemical reaction networks
Zhang, Jiajun; Nie, Qing; Zhou, Tianshou
2016-05-01
Traditional moment-closure methods need to assume that high-order cumulants of a probability distribution approximate to zero. However, this strong assumption is not satisfied for many biochemical reaction networks. Here, we introduce convergent moments (defined in mathematics as the coefficients in the Taylor expansion of the probability-generating function at some point) to overcome this drawback of the moment-closure methods. As such, we develop a new analysis method for stochastic chemical kinetics. This method provides an accurate approximation for the master probability equation (MPE). In particular, the connection between low-order convergent moments and rate constants can be more easily derived in terms of explicit and analytical forms, allowing insights that would be difficult to obtain through direct simulation or manipulation of the MPE. In addition, it provides an accurate and efficient way to compute steady-state or transient probability distribution, avoiding the algorithmic difficulty associated with stiffness of the MPE due to large differences in sizes of rate constants. Applications of the method to several systems reveal nontrivial stochastic mechanisms of gene expression dynamics, e.g., intrinsic fluctuations can induce transient bimodality and amplify transient signals, and slow switching between promoter states can increase fluctuations in spatially heterogeneous signals. The overall approach has broad applications in modeling, analysis, and computation of complex biochemical networks with intrinsic noise.
Lisa M. Bishop
2010-09-01
Full Text Available We develop the stochastic, chemical master equation as a unifying approach to the dynamics of biochemical reaction systems in a mesoscopic volume under a living environment. A living environment provides a continuous chemical energy input that sustains the reaction system in a nonequilibrium steady state with concentration fluctuations. We discuss the linear, unimolecular single-molecule enzyme kinetics, phosphorylation-dephosphorylation cycle (PdPC with bistability, and network exhibiting oscillations. Emphasis is paid to the comparison between the stochastic dynamics and the prediction based on the traditional approach based on the Law of Mass Action. We introduce the difference between nonlinear bistability and stochastic bistability, the latter has no deterministic counterpart. For systems with nonlinear bistability, there are three different time scales: (a individual biochemical reactions, (b nonlinear network dynamics approaching to attractors, and (c cellular evolution. For mesoscopic systems with size of a living cell, dynamics in (a and (c are stochastic while that with (b is dominantly deterministic. Both (b and (c are emergent properties of a dynamic biochemical network; We suggest that the (c is most relevant to major cellular biochemical processes such as epi-genetic regulation, apoptosis, and cancer immunoediting. The cellular evolution proceeds with transitions among the attractors of (b in a “punctuated equilibrium” manner.
Qian, Hong; Bishop, Lisa M
2010-01-01
We develop the stochastic, chemical master equation as a unifying approach to the dynamics of biochemical reaction systems in a mesoscopic volume under a living environment. A living environment provides a continuous chemical energy input that sustains the reaction system in a nonequilibrium steady state with concentration fluctuations. We discuss the linear, unimolecular single-molecule enzyme kinetics, phosphorylation-dephosphorylation cycle (PdPC) with bistability, and network exhibiting oscillations. Emphasis is paid to the comparison between the stochastic dynamics and the prediction based on the traditional approach based on the Law of Mass Action. We introduce the difference between nonlinear bistability and stochastic bistability, the latter has no deterministic counterpart. For systems with nonlinear bistability, there are three different time scales: (a) individual biochemical reactions, (b) nonlinear network dynamics approaching to attractors, and (c) cellular evolution. For mesoscopic systems with size of a living cell, dynamics in (a) and (c) are stochastic while that with (b) is dominantly deterministic. Both (b) and (c) are emergent properties of a dynamic biochemical network; We suggest that the (c) is most relevant to major cellular biochemical processes such as epi-genetic regulation, apoptosis, and cancer immunoediting. The cellular evolution proceeds with transitions among the attractors of (b) in a "punctuated equilibrium" manner. PMID:20957107
Ruess, Jakob
2015-12-01
Many stochastic models of biochemical reaction networks contain some chemical species for which the number of molecules that are present in the system can only be finite (for instance due to conservation laws), but also other species that can be present in arbitrarily large amounts. The prime example of such networks are models of gene expression, which typically contain a small and finite number of possible states for the promoter but an infinite number of possible states for the amount of mRNA and protein. One of the main approaches to analyze such models is through the use of equations for the time evolution of moments of the chemical species. Recently, a new approach based on conditional moments of the species with infinite state space given all the different possible states of the finite species has been proposed. It was argued that this approach allows one to capture more details about the full underlying probability distribution with a smaller number of equations. Here, I show that the result that less moments provide more information can only stem from an unnecessarily complicated description of the system in the classical formulation. The foundation of this argument will be the derivation of moment equations that describe the complete probability distribution over the finite state space but only low-order moments over the infinite state space. I will show that the number of equations that is needed is always less than what was previously claimed and always less than the number of conditional moment equations up to the same order. To support these arguments, a symbolic algorithm is provided that can be used to derive minimal systems of unconditional moment equations for models with partially finite state space.
Dynamic analysis of biochemical network using complex network method
Wang Shuqiang; Shen Yanyan; Hu Jinxing; Li Ning; Zeng Dewei
2015-01-01
In this study, the stochastic biochemical reaction model is proposed based on the law of mass action and complex network theory. The dynamics of biochemical reaction system is presented as a set of non-linear differential equations and analyzed at the molecular-scale. Given the initial state and the evolution rules of the biochemical reaction system, the system can achieve homeostasis. Compared with random graph, the biochemical reaction network has larger ...
Milias-Argeitis, Andreas; Engblom, Stefan; Bauer, Pavol; Khammash, Mustafa
2015-12-01
Nature presents multiple intriguing examples of processes that proceed with high precision and regularity. This remarkable stability is frequently counter to modellers' experience with the inherent stochasticity of chemical reactions in the regime of low-copy numbers. Moreover, the effects of noise and nonlinearities can lead to 'counterintuitive' behaviour, as demonstrated for a basic enzymatic reaction scheme that can display stochastic focusing (SF). Under the assumption of rapid signal fluctuations, SF has been shown to convert a graded response into a threshold mechanism, thus attenuating the detrimental effects of signal noise. However, when the rapid fluctuation assumption is violated, this gain in sensitivity is generally obtained at the cost of very large product variance, and this unpredictable behaviour may be one possible explanation of why, more than a decade after its introduction, SF has still not been observed in real biochemical systems. In this work, we explore the noise properties of a simple enzymatic reaction mechanism with a small and fluctuating number of active enzymes that behaves as a high-gain, noisy amplifier due to SF caused by slow enzyme fluctuations. We then show that the inclusion of a plausible negative feedback mechanism turns the system from a noisy signal detector to a strong homeostatic mechanism by exchanging high gain with strong attenuation in output noise and robustness to parameter variations. Moreover, we observe that the discrepancy between deterministic and stochastic descriptions of stochastically focused systems in the evolution of the means almost completely disappears, despite very low molecule counts and the additional nonlinearity due to feedback. The reaction mechanism considered here can provide a possible resolution to the apparent conflict between intrinsic noise and high precision in critical intracellular processes. PMID:26609065
Dynamic analysis of biochemical network using complex network method
Wang Shuqiang
2015-01-01
Full Text Available In this study, the stochastic biochemical reaction model is proposed based on the law of mass action and complex network theory. The dynamics of biochemical reaction system is presented as a set of non-linear differential equations and analyzed at the molecular-scale. Given the initial state and the evolution rules of the biochemical reaction system, the system can achieve homeostasis. Compared with random graph, the biochemical reaction network has larger information capacity and is more efficient in information transmission. This is consistent with theory of evolution.
RMBNToolbox: random models for biochemical networks
Niemi Jari
2007-05-01
Full Text Available Abstract Background There is an increasing interest to model biochemical and cell biological networks, as well as to the computational analysis of these models. The development of analysis methodologies and related software is rapid in the field. However, the number of available models is still relatively small and the model sizes remain limited. The lack of kinetic information is usually the limiting factor for the construction of detailed simulation models. Results We present a computational toolbox for generating random biochemical network models which mimic real biochemical networks. The toolbox is called Random Models for Biochemical Networks. The toolbox works in the Matlab environment, and it makes it possible to generate various network structures, stoichiometries, kinetic laws for reactions, and parameters therein. The generation can be based on statistical rules and distributions, and more detailed information of real biochemical networks can be used in situations where it is known. The toolbox can be easily extended. The resulting network models can be exported in the format of Systems Biology Markup Language. Conclusion While more information is accumulating on biochemical networks, random networks can be used as an intermediate step towards their better understanding. Random networks make it possible to study the effects of various network characteristics to the overall behavior of the network. Moreover, the construction of artificial network models provides the ground truth data needed in the validation of various computational methods in the fields of parameter estimation and data analysis.
St Laurent Georges
2010-03-01
Full Text Available Abstract Background Signal transduction networks represent the information processing systems that dictate which dynamical regimes of biochemical activity can be accessible to a cell under certain circumstances. One of the major concerns in molecular systems biology is centered on the elucidation of the robustness properties and information processing capabilities of signal transduction networks. Achieving this goal requires the establishment of causal relations between the design principle of biochemical reaction systems and their emergent dynamical behaviors. Methods In this study, efforts were focused in the construction of a relatively well informed, deterministic, non-linear dynamic model, accounting for reaction mechanisms grounded on standard mass action and Hill saturation kinetics, of the canonical reaction topology underlying Toll-like receptor 4 (TLR4-mediated signaling events. This signaling mechanism has been shown to be deployed in macrophages during a relatively short time window in response to lypopolysaccharyde (LPS stimulation, which leads to a rapidly mounted innate immune response. An extensive computational exploration of the biochemical reaction space inhabited by this signal transduction network was performed via local and global perturbation strategies. Importantly, a broad spectrum of biologically plausible dynamical regimes accessible to the network in widely scattered regions of parameter space was reconstructed computationally. Additionally, experimentally reported transcriptional readouts of target pro-inflammatory genes, which are actively modulated by the network in response to LPS stimulation, were also simulated. This was done with the main goal of carrying out an unbiased statistical assessment of the intrinsic robustness properties of this canonical reaction topology. Results Our simulation results provide convincing numerical evidence supporting the idea that a canonical reaction mechanism of the TLR4
Elston Timothy C
2004-03-01
Full Text Available Abstract Background Intrinsic fluctuations due to the stochastic nature of biochemical reactions can have large effects on the response of biochemical networks. This is particularly true for pathways that involve transcriptional regulation, where generally there are two copies of each gene and the number of messenger RNA (mRNA molecules can be small. Therefore, there is a need for computational tools for developing and investigating stochastic models of biochemical networks. Results We have developed the software package Biochemical Network Stochastic Simulator (BioNetS for efficientlyand accurately simulating stochastic models of biochemical networks. BioNetS has a graphical user interface that allows models to be entered in a straightforward manner, and allows the user to specify the type of random variable (discrete or continuous for each chemical species in the network. The discrete variables are simulated using an efficient implementation of the Gillespie algorithm. For the continuous random variables, BioNetS constructs and numerically solvesthe appropriate chemical Langevin equations. The software package has been developed to scale efficiently with network size, thereby allowing large systems to be studied. BioNetS runs as a BioSpice agent and can be downloaded from http://www.biospice.org. BioNetS also can be run as a stand alone package. All the required files are accessible from http://x.amath.unc.edu/BioNetS. Conclusions We have developed BioNetS to be a reliable tool for studying the stochastic dynamics of large biochemical networks. Important features of BioNetS are its ability to handle hybrid models that consist of both continuous and discrete random variables and its ability to model cell growth and division. We have verified the accuracy and efficiency of the numerical methods by considering several test systems.
Programmability of Chemical Reaction Networks
Cook, Matthew; Soloveichik, David; Winfree, Erik; Bruck, Jehoshua
Motivated by the intriguing complexity of biochemical circuitry within individual cells we study Stochastic Chemical Reaction Networks (SCRNs), a formal model that considers a set of chemical reactions acting on a finite number of molecules in a well-stirred solution according to standard chemical kinetics equations. SCRNs have been widely used for describing naturally occurring (bio)chemical systems, and with the advent of synthetic biology they become a promising language for the design of artificial biochemical circuits. Our interest here is the computational power of SCRNs and how they relate to more conventional models of computation. We survey known connections and give new connections between SCRNs and Boolean Logic Circuits, Vector Addition Systems, Petri nets, Gate Implementability, Primitive Recursive Functions, Register Machines, Fractran, and Turing Machines. A theme to these investigations is the thin line between decidable and undecidable questions about SCRN behavior.
Robust simplifications of multiscale biochemical networks
Zinovyev Andrei
2008-10-01
Full Text Available Abstract Background Cellular processes such as metabolism, decision making in development and differentiation, signalling, etc., can be modeled as large networks of biochemical reactions. In order to understand the functioning of these systems, there is a strong need for general model reduction techniques allowing to simplify models without loosing their main properties. In systems biology we also need to compare models or to couple them as parts of larger models. In these situations reduction to a common level of complexity is needed. Results We propose a systematic treatment of model reduction of multiscale biochemical networks. First, we consider linear kinetic models, which appear as "pseudo-monomolecular" subsystems of multiscale nonlinear reaction networks. For such linear models, we propose a reduction algorithm which is based on a generalized theory of the limiting step that we have developed in 1. Second, for non-linear systems we develop an algorithm based on dominant solutions of quasi-stationarity equations. For oscillating systems, quasi-stationarity and averaging are combined to eliminate time scales much faster and much slower than the period of the oscillations. In all cases, we obtain robust simplifications and also identify the critical parameters of the model. The methods are demonstrated for simple examples and for a more complex model of NF-κB pathway. Conclusion Our approach allows critical parameter identification and produces hierarchies of models. Hierarchical modeling is important in "middle-out" approaches when there is need to zoom in and out several levels of complexity. Critical parameter identification is an important issue in systems biology with potential applications to biological control and therapeutics. Our approach also deals naturally with the presence of multiple time scales, which is a general property of systems biology models.
Probabilistic sensitivity analysis of biochemical reaction systems
Zhang, Hong-Xuan; Dempsey, William P.; Goutsias, John
2009-01-01
Sensitivity analysis is an indispensable tool for studying the robustness and fragility properties of biochemical reaction systems as well as for designing optimal approaches for selective perturbation and intervention. Deterministic sensitivity analysis techniques, using derivatives of the system response, have been extensively used in the literature. However, these techniques suffer from several drawbacks, which must be carefully considered before using them in problems of systems biology. ...
Tanase-Nicola, Sorin; Warren, Patrick B.; Wolde, Pieter Rein ten
2005-01-01
Understanding cell function requires an accurate description of how noise is transmitted through biochemical networks. We present an analytical result for the power spectrum of the output signal of a biochemical network that takes into account the correlations between the noise in the input signal (the extrinsic noise) and the noise in the reactions that constitute the network (the intrinsic noise). These correlations arise from the fact that the reactions by which biochemical signals are det...
Autocatalysis in reaction networks.
Deshpande, Abhishek; Gopalkrishnan, Manoj
2014-10-01
The persistence conjecture is a long-standing open problem in chemical reaction network theory. It concerns the behavior of solutions to coupled ODE systems that arise from applying mass-action kinetics to a network of chemical reactions. The idea is that if all reactions are reversible in a weak sense, then no species can go extinct. A notion that has been found useful in thinking about persistence is that of "critical siphon." We explore the combinatorics of critical siphons, with a view toward the persistence conjecture. We introduce the notions of "drainable" and "self-replicable" (or autocatalytic) siphons. We show that: Every minimal critical siphon is either drainable or self-replicable; reaction networks without drainable siphons are persistent; and nonautocatalytic weakly reversible networks are persistent. Our results clarify that the difficulties in proving the persistence conjecture are essentially due to competition between drainable and self-replicable siphons. PMID:25245394
SIMULATING BIOCHEMICAL SIGNALING NETWORKS IN COMPLEX MOVING GEOMETRIES.
Strychalski, Wanda; Adalsteinsson, David; Elston, Timothy C
2010-01-01
Signaling networks regulate cellular responses to environmental stimuli through cascades of protein interactions. External signals can trigger cells to polarize and move in a specific direction. During migration, spatially localized activity of proteins is maintained. To investigate the effects of morphological changes on intracellular signaling, we developed a numerical scheme consisting of a cut cell finite volume spatial discretization coupled with level set methods to simulate the resulting advection-reaction-diffusion system. We then apply the method to several biochemical reaction networks in changing geometries. We found that a Turing instability can develop exclusively by cell deformations that maintain constant area. For a Turing system with a geometry-dependent single or double peak solution, simulations in a dynamically changing geometry suggest that a single peak solution is the only stable one, independent of the oscillation frequency. The method is also applied to a model of a signaling network in a migrating fibroblast. PMID:24086102
Transient amplification limits noise suppression in biochemical networks
Dixon, John; Lindemann, Anika; McCoy, Jonathan H.
2016-01-01
Cell physiology is orchestrated, on a molecular level, through complex networks of biochemical reactions. The propagation of random fluctuations through these networks can significantly impact cell behavior, raising challenging questions about how network design shapes the cell's ability to suppress or exploit these fluctuations. Here, drawing on insights from statistical physics, fluid dynamics, and systems biology, we explore how transient amplification phenomena arising from network connectivity naturally limit a biochemical system's ability to suppress small fluctuations around steady-state behaviors. We find that even a simple system consisting of two variables linked by a single interaction is capable of amplifying small fluctuations orders of magnitude beyond the levels predicted by linear stability theory. We also find that adding additional interactions can promote further amplification, even when these interactions implement classic design strategies known to suppress fluctuations. These results establish that transient amplification is an essential factor determining baseline noise levels in stable intracellular networks. Significantly, our analysis is not bound to specific systems or interaction mechanisms: we find that noise amplification is an emergent phenomenon found near steady states in any network containing sufficiently strong interactions, regardless of its form or function.
BNDB – The Biochemical Network Database
Kaufmann Michael
2007-10-01
Full Text Available Abstract Background Technological advances in high-throughput techniques and efficient data acquisition methods have resulted in a massive amount of life science data. The data is stored in numerous databases that have been established over the last decades and are essential resources for scientists nowadays. However, the diversity of the databases and the underlying data models make it difficult to combine this information for solving complex problems in systems biology. Currently, researchers typically have to browse several, often highly focused, databases to obtain the required information. Hence, there is a pressing need for more efficient systems for integrating, analyzing, and interpreting these data. The standardization and virtual consolidation of the databases is a major challenge resulting in a unified access to a variety of data sources. Description We present the Biochemical Network Database (BNDB, a powerful relational database platform, allowing a complete semantic integration of an extensive collection of external databases. BNDB is built upon a comprehensive and extensible object model called BioCore, which is powerful enough to model most known biochemical processes and at the same time easily extensible to be adapted to new biological concepts. Besides a web interface for the search and curation of the data, a Java-based viewer (BiNA provides a powerful platform-independent visualization and navigation of the data. BiNA uses sophisticated graph layout algorithms for an interactive visualization and navigation of BNDB. Conclusion BNDB allows a simple, unified access to a variety of external data sources. Its tight integration with the biochemical network library BN++ offers the possibility for import, integration, analysis, and visualization of the data. BNDB is freely accessible at http://www.bndb.org.
Markovian Dynamics on Complex Reaction Networks
Goutsias, John
2012-01-01
Complex networks, comprised of individual elements that interact with each other through reaction channels, are ubiquitous across many scientific and engineering disciplines. Examples include biochemical, pharmacokinetic, epidemiological, ecological, social, neural, and multi-agent networks. A common approach to modeling such networks is by a master equation that governs the dynamic evolution of the joint probability mass function of the underling population process and naturally leads to Markovian dynamics for such process. Due however to the nonlinear nature of most reactions, the computation and analysis of the resulting stochastic population dynamics is a difficult task. This review article provides a coherent and comprehensive coverage of recently developed approaches and methods to tackle this problem. After reviewing a general framework for modeling Markovian reaction networks and giving specific examples, the authors present numerical and computational techniques capable of evaluating or approximating...
Markovian dynamics on complex reaction networks
Goutsias, J., E-mail: goutsias@jhu.edu; Jenkinson, G., E-mail: jenkinson@jhu.edu
2013-08-10
Complex networks, comprised of individual elements that interact with each other through reaction channels, are ubiquitous across many scientific and engineering disciplines. Examples include biochemical, pharmacokinetic, epidemiological, ecological, social, neural, and multi-agent networks. A common approach to modeling such networks is by a master equation that governs the dynamic evolution of the joint probability mass function of the underlying population process and naturally leads to Markovian dynamics for such process. Due however to the nonlinear nature of most reactions and the large size of the underlying state-spaces, computation and analysis of the resulting stochastic population dynamics is a difficult task. This review article provides a coherent and comprehensive coverage of recently developed approaches and methods to tackle this problem. After reviewing a general framework for modeling Markovian reaction networks and giving specific examples, the authors present numerical and computational techniques capable of evaluating or approximating the solution of the master equation, discuss a recently developed approach for studying the stationary behavior of Markovian reaction networks using a potential energy landscape perspective, and provide an introduction to the emerging theory of thermodynamic analysis of such networks. Three representative problems of opinion formation, transcription regulation, and neural network dynamics are used as illustrative examples.
Markovian dynamics on complex reaction networks
Goutsias, J.; Jenkinson, G.
2013-08-01
Complex networks, comprised of individual elements that interact with each other through reaction channels, are ubiquitous across many scientific and engineering disciplines. Examples include biochemical, pharmacokinetic, epidemiological, ecological, social, neural, and multi-agent networks. A common approach to modeling such networks is by a master equation that governs the dynamic evolution of the joint probability mass function of the underlying population process and naturally leads to Markovian dynamics for such process. Due however to the nonlinear nature of most reactions and the large size of the underlying state-spaces, computation and analysis of the resulting stochastic population dynamics is a difficult task. This review article provides a coherent and comprehensive coverage of recently developed approaches and methods to tackle this problem. After reviewing a general framework for modeling Markovian reaction networks and giving specific examples, the authors present numerical and computational techniques capable of evaluating or approximating the solution of the master equation, discuss a recently developed approach for studying the stationary behavior of Markovian reaction networks using a potential energy landscape perspective, and provide an introduction to the emerging theory of thermodynamic analysis of such networks. Three representative problems of opinion formation, transcription regulation, and neural network dynamics are used as illustrative examples.
Coarse-graining stochastic biochemical networks: adiabaticity and fast simulations
Nemenman, Ilya [Los Alamos National Laboratory; Sinitsyn, Nikolai [Los Alamos National Laboratory; Hengartner, Nick [Los Alamos National Laboratory
2008-01-01
We propose a universal approach for analysis and fast simulations of stiff stochastic biochemical kinetics networks, which rests on elimination of fast chemical species without a loss of information about mesoscoplc, non-Poissonian fluctuations of the slow ones. Our approach, which is similar to the Born-Oppenhelmer approximation in quantum mechanics, follows from the stochastic path Integral representation of the cumulant generating function of reaction events. In applications with a small number of chemIcal reactions, It produces analytical expressions for cumulants of chemical fluxes between the slow variables. This allows for a low-dimensional, Interpretable representation and can be used for coarse-grained numerical simulation schemes with a small computational complexity and yet high accuracy. As an example, we derive the coarse-grained description for a chain of biochemical reactions, and show that the coarse-grained and the microscopic simulations are in an agreement, but the coarse-gralned simulations are three orders of magnitude faster.
Thermodynamically consistent Bayesian analysis of closed biochemical reaction systems
Goutsias John
2010-11-01
Full Text Available Abstract Background Estimating the rate constants of a biochemical reaction system with known stoichiometry from noisy time series measurements of molecular concentrations is an important step for building predictive models of cellular function. Inference techniques currently available in the literature may produce rate constant values that defy necessary constraints imposed by the fundamental laws of thermodynamics. As a result, these techniques may lead to biochemical reaction systems whose concentration dynamics could not possibly occur in nature. Therefore, development of a thermodynamically consistent approach for estimating the rate constants of a biochemical reaction system is highly desirable. Results We introduce a Bayesian analysis approach for computing thermodynamically consistent estimates of the rate constants of a closed biochemical reaction system with known stoichiometry given experimental data. Our method employs an appropriately designed prior probability density function that effectively integrates fundamental biophysical and thermodynamic knowledge into the inference problem. Moreover, it takes into account experimental strategies for collecting informative observations of molecular concentrations through perturbations. The proposed method employs a maximization-expectation-maximization algorithm that provides thermodynamically feasible estimates of the rate constant values and computes appropriate measures of estimation accuracy. We demonstrate various aspects of the proposed method on synthetic data obtained by simulating a subset of a well-known model of the EGF/ERK signaling pathway, and examine its robustness under conditions that violate key assumptions. Software, coded in MATLAB®, which implements all Bayesian analysis techniques discussed in this paper, is available free of charge at http://www.cis.jhu.edu/~goutsias/CSS%20lab/software.html. Conclusions Our approach provides an attractive statistical methodology for
Limits for Stochastic Reaction Networks
Cappelletti, Daniele
Reaction systems have been introduced in the 70s to model biochemical systems. Nowadays their range of applications has increased and they are fruitfully used in dierent elds. The concept is simple: some chemical species react, the set of chemical reactions form a graph and a rate function is...... associated with each reaction. Such functions describe the speed of the dierent reactions, or their propensities. Two modelling regimes are then available: the evolution of the dierent species concentrations can be deterministically modelled through a system of ODE, while the counts of the dierent species at...... a certain time are stochastically modelled by means of a continuous-time Markov chain. Our work concerns primarily stochastic reaction systems, and their asymptotic properties. In Paper I, we consider a reaction system with intermediate species, i.e. species that are produced and fast degraded along...
Modeling of fluctuating reaction networks
Full Text:Various dynamical systems are organized as reaction networks, where the population size of one component affects the populations of all its neighbors. Such networks can be found in interstellar surface chemistry, cell biology, thin film growth and other systems. I cases where the populations of reactive species are large, the network can be modeled by rate equations which provide all reaction rates within mean field approximation. However, in small systems that are partitioned into sub-micron size, these populations strongly fluctuate. Under these conditions rate equations fail and the master equation is needed for modeling these reactions. However, the number of equations in the master equation grows exponentially with the number of reactive species, severely limiting its feasibility for complex networks. Here we present a method which dramatically reduces the number of equations, thus enabling the incorporation of the master equation in complex reaction networks. The method is examplified in the context of reaction network on dust grains. Its applicability for genetic networks will be discussed. 1. Efficient simulations of gas-grain chemistry in interstellar clouds. Azi Lipshtat and Ofer Biham, Phys. Rev. Lett. 93 (2004), 170601. 2. Modeling of negative autoregulated genetic networks in single cells. Azi Lipshtat, Hagai B. Perets, Nathalie Q. Balaban and Ofer Biham, Gene: evolutionary genomics (2004), In press
Modelling biochemical reaction systems by stochastic differential equations with reflection.
Niu, Yuanling; Burrage, Kevin; Chen, Luonan
2016-05-01
In this paper, we gave a new framework for modelling and simulating biochemical reaction systems by stochastic differential equations with reflection not in a heuristic way but in a mathematical way. The model is computationally efficient compared with the discrete-state Markov chain approach, and it ensures that both analytic and numerical solutions remain in a biologically plausible region. Specifically, our model mathematically ensures that species numbers lie in the domain D, which is a physical constraint for biochemical reactions, in contrast to the previous models. The domain D is actually obtained according to the structure of the corresponding chemical Langevin equations, i.e., the boundary is inherent in the biochemical reaction system. A variant of projection method was employed to solve the reflected stochastic differential equation model, and it includes three simple steps, i.e., Euler-Maruyama method was applied to the equations first, and then check whether or not the point lies within the domain D, and if not perform an orthogonal projection. It is found that the projection onto the closure D¯ is the solution to a convex quadratic programming problem. Thus, existing methods for the convex quadratic programming problem can be employed for the orthogonal projection map. Numerical tests on several important problems in biological systems confirmed the efficiency and accuracy of this approach. PMID:26920245
An effective method for computing the noise in biochemical networks
Zhang, Jiajun; Nie, Qing; He, Miao; Zhou, Tianshou
2013-02-01
We present a simple yet effective method, which is based on power series expansion, for computing exact binomial moments that can be in turn used to compute steady-state probability distributions as well as the noise in linear or nonlinear biochemical reaction networks. When the method is applied to representative reaction networks such as the ON-OFF models of gene expression, gene models of promoter progression, gene auto-regulatory models, and common signaling motifs, the exact formulae for computing the intensities of noise in the species of interest or steady-state distributions are analytically given. Interestingly, we find that positive (negative) feedback does not enlarge (reduce) noise as claimed in previous works but has a counter-intuitive effect and that the multi-OFF (or ON) mechanism always attenuates the noise in contrast to the common ON-OFF mechanism and can modulate the noise to the lowest level independently of the mRNA mean. Except for its power in deriving analytical expressions for distributions and noise, our method is programmable and has apparent advantages in reducing computational cost.
Catalytic reaction dynamics in inhomogeneous networks.
Watanabe, Akitomo; Yakubo, Kousuke
2014-05-01
Biochemical reactions in a cell can be modeled by a catalytic reaction network (CRN). It has been reported that catalytic chain reactions occur intermittently in the CRN with a homogeneous random-graph topology and its avalanche-size distribution obeys a power law with the exponent 4/3 [A. Awazu and K. Kaneko, Phys. Rev. E 80, 010902(R) (2009)]. This fact indicates that the catalytic reaction dynamics in homogeneous CRNs exhibits self-organized criticality (SOC). Structures of actual CRNs are, however, known to be highly inhomogeneous. We study the influence of various types of inhomogeneities found in real-world metabolic networks on the universality class of SOC. Our numerical results clarify that SOC keeps its universality class even for networks possessing structural inhomogeneities such as the scale-free property, community structures, and degree correlations. In contrast, if the CRN has inhomogeneous catalytic functionality, the universality class of SOC depends on how widely distributed the number of reaction paths catalyzed by a single chemical species is. PMID:25353843
Simulation of biochemical reactions with time-dependent rates by the rejection-based algorithm
We address the problem of simulating biochemical reaction networks with time-dependent rates and propose a new algorithm based on our rejection-based stochastic simulation algorithm (RSSA) [Thanh et al., J. Chem. Phys. 141(13), 134116 (2014)]. The computation for selecting next reaction firings by our time-dependent RSSA (tRSSA) is computationally efficient. Furthermore, the generated trajectory is exact by exploiting the rejection-based mechanism. We benchmark tRSSA on different biological systems with varying forms of reaction rates to demonstrate its applicability and efficiency. We reveal that for nontrivial cases, the selection of reaction firings in existing algorithms introduces approximations because the integration of reaction rates is very computationally demanding and simplifying assumptions are introduced. The selection of the next reaction firing by our approach is easier while preserving the exactness
Ubiquitous ``glassy'' relaxation in catalytic reaction networks
Awazu, Akinori; Kaneko, Kunihiko
2009-10-01
Study of reversible catalytic reaction networks is important not only as an issue for chemical thermodynamics but also for protocells. From extensive numerical simulations and theoretical analysis, slow relaxation dynamics to sustain nonequlibrium states are commonly observed. These dynamics show two types of salient behaviors that are reminiscent of glassy behavior: slow relaxation along with the logarithmic time dependence of the correlation function and the emergence of plateaus in the relaxation-time course. The former behavior is explained by the eigenvalue distribution of a Jacobian matrix around the equilibrium state that depends on the distribution of kinetic coefficients of reactions. The latter behavior is associated with kinetic constraints rather than metastable states and is due to the absence of catalysts for chemicals in excess and the negative correlation between two chemical species. Examples are given and generality is discussed with relevance to bottleneck-type dynamics in biochemical reactions as well.
Beard, Daniel A.; Liang, Shou-Dan; Qian, Hong; Biegel, Bryan (Technical Monitor)
2001-01-01
Predicting behavior of large-scale biochemical metabolic networks represents one of the greatest challenges of bioinformatics and computational biology. Approaches, such as flux balance analysis (FBA), that account for the known stoichiometry of the reaction network while avoiding implementation of detailed reaction kinetics are perhaps the most promising tools for the analysis of large complex networks. As a step towards building a complete theory of biochemical circuit analysis, we introduce energy balance analysis (EBA), which compliments the FBA approach by introducing fundamental constraints based on the first and second laws of thermodynamics. Fluxes obtained with EBA are thermodynamically feasible and provide valuable insight into the activation and suppression of biochemical pathways.
In silico evolution of biochemical networks
Francois, Paul
2010-03-01
We use computational evolution to select models of genetic networks that can be built from a predefined set of parts to achieve a certain behavior. Selection is made with the help of a fitness defining biological functions in a quantitative way. This fitness has to be specific to a process, but general enough to find processes common to many species. Computational evolution favors models that can be built by incremental improvements in fitness rather than via multiple neutral steps or transitions through less fit intermediates. With the help of these simulations, we propose a kinetic view of evolution, where networks are rapidly selected along a fitness gradient. This mathematics recapitulates Darwin's original insight that small changes in fitness can rapidly lead to the evolution of complex structures such as the eye, and explain the phenomenon of convergent/parallel evolution of similar structures in independent lineages. We will illustrate these ideas with networks implicated in embryonic development and patterning of vertebrates and primitive insects.
LucidDraw: Efficiently visualizing complex biochemical networks within MATLAB
Shi Guiyang
2010-01-01
Full Text Available Abstract Background Biochemical networks play an essential role in systems biology. Rapidly growing network data and versatile research activities call for convenient visualization tools to aid intuitively perceiving abstract structures of networks and gaining insights into the functional implications of networks. There are various kinds of network visualization software, but they are usually not adequate for visual analysis of complex biological networks mainly because of the two reasons: 1 most existing drawing methods suitable for biochemical networks have high computation loads and can hardly achieve near real-time visualization; 2 available network visualization tools are designed for working in certain network modeling platforms, so they are not convenient for general analyses due to lack of broader range of readily accessible numerical utilities. Results We present LucidDraw as a visual analysis tool, which features (a speed: typical biological networks with several hundreds of nodes can be drawn in a few seconds through a new layout algorithm; (b ease of use: working within MATLAB makes it convenient to manipulate and analyze the network data using a broad spectrum of sophisticated numerical functions; (c flexibility: layout styles and incorporation of other available information about functional modules can be controlled by users with little effort, and the output drawings are interactively modifiable. Conclusions Equipped with a new grid layout algorithm proposed here, LucidDraw serves as an auxiliary network analysis tool capable of visualizing complex biological networks in near real-time with controllable layout styles and drawing details. The framework of the algorithm enables easy incorporation of extra biological information, if available, to influence the output layouts with predefined node grouping features.
Finite time thermodynamic coupling in a biochemical network
Dasgupta, Anjan Kr
2014-01-01
The paper describes some thermodynamic constrains and relations in biochemical or metabolic network and provides a basis for entropy enthalpy compensation. Conventional definition of macroscopic forces and fluxes leads to a paradox namely, non-existence of positive efficiency of a chemically driven process. This paradox is resolved by deriving an appropriate definition of macroscopic force using the local balance equations. Entropy enthalpy compensation, whose thermodynamic basis is so far un...
Feedback Regulation and Its Efficiency in Biochemical Networks
Kobayashi, Tetsuya J.; Yokota, Ryo; Aihara, Kazuyuki
2016-03-01
Intracellular biochemical networks fluctuate dynamically due to various internal and external sources of fluctuation. Dissecting the fluctuation into biologically relevant components is important for understanding how a cell controls and harnesses noise and how information is transferred over apparently noisy intracellular networks. While substantial theoretical and experimental advancement on the decomposition of fluctuation was achieved for feedforward networks without any loop, we still lack a theoretical basis that can consistently extend such advancement to feedback networks. The main obstacle that hampers is the circulative propagation of fluctuation by feedback loops. In order to define the relevant quantity for the impact of feedback loops for fluctuation, disentanglement of the causally interlocked influences between the components is required. In addition, we also lack an approach that enables us to infer non-perturbatively the influence of the feedback to fluctuation in the same way as the dual reporter system does in the feedforward networks. In this work, we address these problems by extending the work on the fluctuation decomposition and the dual reporter system. For a single-loop feedback network with two components, we define feedback loop gain as the feedback efficiency that is consistent with the fluctuation decomposition for feedforward networks. Then, we clarify the relation of the feedback efficiency with the fluctuation propagation in an open-looped FF network. Finally, by extending the dual reporter system, we propose a conjugate feedback and feedforward system for estimating the feedback efficiency non-perturbatively only from the statistics of the system.
Graphical reduction of reaction networks by linear elimination of species
Sáez, Meritxell; Wiuf, Carsten; Feliu, Elisenda
2015-01-01
The quasi-steady state approximation and time-scale separation are commonly applied methods to simplify models of biochemical reaction networks based on ordinary differential equations (ODEs). The concentrations of the "fast" species are assumed effectively to be at steady state with respect to the "slow" species. Under this assumption the steady state equations can be used to eliminate the "fast" variables and a new ODE system with only the slow species can be obtained. We interpret a reduce...
Modeling the spatiotemporal dynamics of biochemical reaction systems at single-molecule resolution has become feasible with the increase of computing power and is applied especially to cellular signal transduction. For an association reaction the two molecules have to be in contact. Hence, a physically faithful model of the molecular interaction assumes non-overlapping molecules that interact at their surfaces (boundary scheme). For performance reasons, this model can be replaced by particles that can overlap and react when they are closer than a certain distance with a reaction probability (volume scheme). Here we present an analytical approximation for the reaction probability in the volume scheme and compare the volume- with the boundary scheme. A dissociation reaction, in contrast, creates two molecules next to each other. If the reaction is reversible, these two products can directly re-bind again, leading to an overestimation of the dimerized state in the simulation. We show how the correct recombination rate can be achieved if the products of the dissociation are placed at identical positions, but cannot react for a certain timespan. This refractory time corresponds to the completion of the diffusion-controlled dissociation of the two molecules to their contact distance ri + rj at t = τ-hat ×(ri+rj)2/(Di+Dj) with τ-hat = 1/10 for molecules with radii ri and rj and diffusion coefficients Di and Dj, respectively. (paper)
Molecular codes in biological and chemical reaction networks.
Dennis Görlich
Full Text Available Shannon's theory of communication has been very successfully applied for the analysis of biological information. However, the theory neglects semantic and pragmatic aspects and thus cannot directly be applied to distinguish between (bio- chemical systems able to process "meaningful" information from those that do not. Here, we present a formal method to assess a system's semantic capacity by analyzing a reaction network's capability to implement molecular codes. We analyzed models of chemical systems (martian atmosphere chemistry and various combustion chemistries, biochemical systems (gene expression, gene translation, and phosphorylation signaling cascades, an artificial chemistry, and random reaction networks. Our study suggests that different chemical systems possess different semantic capacities. No semantic capacity was found in the model of the martian atmosphere chemistry, the studied combustion chemistries, and highly connected random networks, i.e. with these chemistries molecular codes cannot be implemented. High semantic capacity was found in the studied biochemical systems and in random reaction networks where the number of second order reactions is twice the number of species. We conclude that our approach can be applied to evaluate the information processing capabilities of a chemical system and may thus be a useful tool to understand the origin and evolution of meaningful information, e.g. in the context of the origin of life.
Genetic Recombination as a Chemical Reaction Network
Müller, Stefan; Hofbauer, Josef
2015-01-01
The process of genetic recombination can be seen as a chemical reaction network with mass-action kinetics. We review the known results on existence, uniqueness, and global stability of an equilibrium in every compatibility class and for all rate constants, from both the population genetics and the reaction networks point of view.
DNA reaction networks: Providing a panoramic view
Wang, Fei; Fan, Chunhai
2016-08-01
A quantitative understanding of the functional landscape of a biochemical circuit can reveal the design rules required to optimize the circuit. Now, a high-throughput droplet-based microfluidic platform has been developed which enables high-resolution mapping of bifurcation diagrams for two nonlinear DNA networks.
Stochastic thermodynamics of chemical reaction networks
Schmiedl, Tim; Seifert, Udo
2006-01-01
For chemical reaction networks described by a master equation, we define energy and entropy on a stochastic trajectory and develop a consistent nonequilibrium thermodynamic description along a single stochastic trajectory of reaction events. A first-law like energy balance relates internal energy, applied (chemical) work and dissipated heat for every single reaction. Entropy production along a single trajectory involves a sum over changes in the entropy of the network itself and the entropy o...
Antibiogram, Biochemical Reactions and Biotyping of Biofield Treated Providencia rettgeri
Trivedi, Mahendra Kumar
2015-01-01
Providencia rettgeri (P. rettgeri) is the key organism for gastrointestinal tract infections due to its high virulence properties. The current study was designed to investigate the effect of Mr. Trivedi’s biofield energy treatment on P. rettgeri in lyophilized as well as revived state for antimicrobial susceptibility pattern, biochemical characteristics, and biotype number. The lyophilized strain of P. rettgeri (ATCC 9250) was divided into two parts, Group (Gr.)...
A network dynamics approach to chemical reaction networks
van der Schaft, Abraham; Rao, S.; Jayawardhana, B.
2016-01-01
A treatment of chemical reaction network theory is given from the perspective of nonlinear network dynamics, in particular of consensus dynamics. By starting from the complex-balanced assumption the reaction dynamics governed by mass action kinetics can be rewritten into a form which allows for a ve
Stochastic Modeling and Simulation of Reaction-Diffusion Biochemical Systems
LI Fei
2016-01-01
Reaction Diffusion Master Equation (RDME) framework, characterized by the discretization of the spatial domain, is one of the most widely used methods in the stochastic simulation of reaction-diffusion systems. Discretization sizes for RDME have to be appropriately chosen such that each discrete compartment is "well-stirred" and the computational cost is not too expensive. An efficient discretization size based on the reaction-diffusion dynamics of each species is derived in this disserta...
Qian Hong
2008-05-01
Full Text Available Abstract Background: Several approaches, including metabolic control analysis (MCA, flux balance analysis (FBA, correlation metric construction (CMC, and biochemical circuit theory (BCT, have been developed for the quantitative analysis of complex biochemical networks. Here, we present a comprehensive theory of linear analysis for nonequilibrium steady-state (NESS biochemical reaction networks that unites these disparate approaches in a common mathematical framework and thermodynamic basis. Results: In this theory a number of relationships between key matrices are introduced: the matrix A obtained in the standard, linear-dynamic-stability analysis of the steady-state can be decomposed as A = SRT where R and S are directly related to the elasticity-coefficient matrix for the fluxes and chemical potentials in MCA, respectively; the control-coefficients for the fluxes and chemical potentials can be written in terms of RT BS and ST BS respectively where matrix B is the inverse of A; the matrix S is precisely the stoichiometric matrix in FBA; and the matrix eAt plays a central role in CMC. Conclusion: One key finding that emerges from this analysis is that the well-known summation theorems in MCA take different forms depending on whether metabolic steady-state is maintained by flux injection or concentration clamping. We demonstrate that if rate-limiting steps exist in a biochemical pathway, they are the steps with smallest biochemical conductances and largest flux control-coefficients. We hypothesize that biochemical networks for cellular signaling have a different strategy for minimizing energy waste and being efficient than do biochemical networks for biosynthesis. We also discuss the intimate relationship between MCA and biochemical systems analysis (BSA.
Modeling isotopomer distributions in biochemical networks using isotopomer mapping matrices
Schmidt, Karsten; Carlsen, Morten; Nielsen, Jens Bredal;
1997-01-01
Within the last decades NMR spectroscopy has undergone tremendous development and has become a powerful analytical tool for the investigation of intracellular flux distributions in biochemical networks using C-13-labeled substrates. Not only are the experiments much easier to conduct than...... experiments employing radioactive tracer elements, but NMR spectroscopy also provides additional information on the labeling pattern of the metabolites. Whereas the maximum amount of information obtainable with C-14-labeled substrates is the fractional enrichment in the individual carbon atom positions, NMR...... spectroscopy can also provide information on the degree of labeling at neighboring carbon atom positions by analyzing multiplet patterns in NMR spectra or using 2-dimensional NMR spectra. It is possible to quantify the mole fractions of molecules that show a specific labeling pattern, i.e., information of the...
Multilayer Network Analysis of Nuclear Reactions.
Zhu, Liang; Ma, Yu-Gang; Chen, Qu; Han, Ding-Ding
2016-01-01
The nuclear reaction network is usually studied via precise calculation of differential equation sets, and much research interest has been focused on the characteristics of nuclides, such as half-life and size limit. In this paper, however, we adopt the methods from both multilayer and reaction networks, and obtain a distinctive view by mapping all the nuclear reactions in JINA REACLIB database into a directed network with 4 layers: neutron, proton, (4)He and the remainder. The layer names correspond to reaction types decided by the currency particles consumed. This combined approach reveals that, in the remainder layer, the β-stability has high correlation with node degree difference and overlapping coefficient. Moreover, when reaction rates are considered as node strength, we find that, at lower temperatures, nuclide half-life scales reciprocally with its out-strength. The connection between physical properties and topological characteristics may help to explore the boundary of the nuclide chart. PMID:27558995
Multilayer Network Analysis of Nuclear Reactions
Zhu, Liang; Ma, Yu-Gang; Chen, Qu; Han, Ding-Ding
2016-01-01
The nuclear reaction network is usually studied via precise calculation of differential equation sets, and much research interest has been focused on the characteristics of nuclides, such as half-life and size limit. In this paper, however, we adopt the methods from both multilayer and reaction networks, and obtain a distinctive view by mapping all the nuclear reactions in JINA REACLIB database into a directed network with 4 layers: neutron, proton, 4He and the remainder. The layer names correspond to reaction types decided by the currency particles consumed. This combined approach reveals that, in the remainder layer, the β-stability has high correlation with node degree difference and overlapping coefficient. Moreover, when reaction rates are considered as node strength, we find that, at lower temperatures, nuclide half-life scales reciprocally with its out-strength. The connection between physical properties and topological characteristics may help to explore the boundary of the nuclide chart. PMID:27558995
Multilayer network analysis of nuclear reactions
Zhu, Liang; Chen, Qu; Han, Ding-Ding
2016-01-01
The nuclear reaction network is usually studied via precise calculation of differential equation sets, and much research interest has been focused on the characteristics of nuclides, such as half-life and size limit. In this paper, however, we adopt the methods from both multilayer and reaction networks, and obtain a distinctive view by mapping all the nuclear reactions in JINA REACLIB database into a directed network with 4 layers: neutron, proton, $^4$He and the remainder. The layer names correspond to reaction types decided by the currency particles consumed. This combined approach reveals that, in the remainder layer, the $\\beta$-stability has high correlation with node degree difference and overlapping coefficient. Moreover, when reaction rates are considered as node strength, we find that, at lower temperatures, nuclide half-life scales reciprocally with its out-strength. The connection between physical properties and topological characteristics may help to explore the boundary of the nuclide chart.
Unveiling the hidden structure of complex stochastic biochemical networks
Valleriani, Angelo; Li, Xin; Kolomeisky, Anatoly B.
2014-02-01
Complex Markov models are widely used and powerful predictive tools to analyze stochastic biochemical processes. However, when the network of states is unknown, it is necessary to extract information from the data to partially build the network and estimate the values of the rates. The short-time behavior of the first-passage time distributions between two states in linear chains has been shown recently to behave as a power of time with an exponent equal to the number of intermediate states. For a general Markov model we derive the complete Taylor expansion of the first-passage time distribution between two arbitrary states. By combining algebraic methods and graph theory approaches it is shown that the first term of the Taylor expansion is determined by the shortest path from the initial state to the final state. When this path is unique, we prove that the coefficient of the first term can be written in terms of the product of the transition rates along the path. It is argued that the application of our results to first-return times may be used to estimate the dependence of rates on external parameters in experimentally measured time distributions.
Trade-Offs in Delayed Information Transmission in Biochemical Networks
Mancini, F.; Marsili, M.; Walczak, A. M.
2016-03-01
In order to transmit biochemical signals, biological regulatory systems dissipate energy with concomitant entropy production. Additionally, signaling often takes place in challenging environmental conditions. In a simple model regulatory circuit given by an input and a delayed output, we explore the trade-offs between information transmission and the system's energetic efficiency. We determine the maximally informative network, given a fixed amount of entropy production and a delayed response, exploring both the case with and without feedback. We find that feedback allows the circuit to overcome energy constraints and transmit close to the maximum available information even in the dissipationless limit. Negative feedback loops, characteristic of shock responses, are optimal at high dissipation. Close to equilibrium positive feedback loops, known for their stability, become more informative. Asking how the signaling network should be constructed to best function in the worst possible environment, rather than an optimally tuned one or in steady state, we discover that at large dissipation the same universal motif is optimal in all of these conditions.
Unveiling the hidden structure of complex stochastic biochemical networks
Complex Markov models are widely used and powerful predictive tools to analyze stochastic biochemical processes. However, when the network of states is unknown, it is necessary to extract information from the data to partially build the network and estimate the values of the rates. The short-time behavior of the first-passage time distributions between two states in linear chains has been shown recently to behave as a power of time with an exponent equal to the number of intermediate states. For a general Markov model we derive the complete Taylor expansion of the first-passage time distribution between two arbitrary states. By combining algebraic methods and graph theory approaches it is shown that the first term of the Taylor expansion is determined by the shortest path from the initial state to the final state. When this path is unique, we prove that the coefficient of the first term can be written in terms of the product of the transition rates along the path. It is argued that the application of our results to first-return times may be used to estimate the dependence of rates on external parameters in experimentally measured time distributions
Law of Localization in Chemical Reaction Networks
Okada, Takashi; Mochizuki, Atsushi
2016-07-01
In living cells, chemical reactions are connected by sharing their products and substrates, and form complex networks, e.g., metabolic pathways. Here we developed a theory to predict the sensitivity, i.e., the responses of concentrations and fluxes to perturbations of enzymes, from network structure alone. Nonzero response patterns turn out to exhibit two characteristic features, localization and hierarchy. We present a general theorem connecting sensitivity with network topology that explains these characteristic patterns. Our results imply that network topology is an origin of biological robustness. Finally, we suggest a strategy to determine real networks from experimental measurements.
Law of localization in chemical reaction networks
Okada, Takashi
2016-01-01
In living cells, chemical reactions are connected by sharing their products and substrates, and form complex networks, e.g. metabolic pathways. Here we developed a theory to predict the sensitivity, i.e. the responses of concentrations and fluxes to perturbations of enzymes, from network structure alone. Responses turn out to exhibit two characteristic patterns, $localization$ and $hierarchy$. We present a general theorem connecting sensitivity with network topology that explains these characteristic patterns. Our results imply that network topology is an origin of biological robustness. Finally, we suggest a strategy to determine real networks from experimental measurements.
Temperature Control System for Biochemical Reactions in Microchip-Based Devices
荆高山; 张坚; 朱小山; 冯继宏; 谭智敏; 刘理天; 程京
2001-01-01
A silicon-glass chip based microreactor has been designed and fabricated for biochemical reactions such as polymerase chain reactions (PCR). The chip based microreactor has integrated resistive heating elements. The computer-controlled temperature control system is highly reliable with precise temperature control, excellent temperature uniformity, and rapid heating and cooling capabilities. The development of the microreaction system is an important step towards the construction of a lab-on-a-chip system.
Radical-ion-pair reactions are the biochemical equivalent of the optical double slit experiment
Kominis, I. K.
2010-01-01
Radical-ion-pair reactions were recently shown to represent a rich biophysical laboratory for the application of quantum measurement theory methods and concepts. We here show that radical-ion-pair reactions essentially form a non-linear biochemical double slit interferometer. Quantum coherence effects are visible when "which-path" information is limited, and the incoherent limit is approached when measurement-induced decoherence sets in. Based on this analogy with the optical double slit expe...
Gökalp, Hüsnü Yusuf; Hüdayi ERCOŞKUN; Ahmet Hilmi ÇON
1998-01-01
Glycolysis, lipolysis and proteolysis are the main biochemical reactions which effect characteristic flavour, aroma, color and texture development of fermented meat products. These reactions are results of the microorganisms or endogenous meat enzymes. The products which are results of glycolysis are firstly responsible from decreasing of pH. Volatile compounds from glycolysis effect aroma. Compounds formed by lipolysis are the main reason that effect aroma. 60 % volatile compounds formed...
田允; 黄继云; 王锐; 陶蓉蓉; 卢应梅; 廖美华; 陆楠楠; 李静; 芦博; 韩峰
2012-01-01
Autism is a highly heritable neurodevelopmental condition characterized by impaired social interaction and communication. However, the role of synaptic dysfunction during development of autism remains unclear. In the present study, we address the alterations of biochemical signaling in hippocampal network following induction of the autism in experimental animals. Here, the an- imal disease model and DNA array being used to investigate the differences in transcriptome or- ganization between autistic and normal brain by gene co--expression network analysis.
Reduced models of networks of coupled enzymatic reactions
Kumar, Ajit
2011-01-01
The Michaelis-Menten equation has played a central role in our understanding of biochemical processes. It has long been understood how this equation approximates the dynamics of irreversible enzymatic reactions. However, a similar approximation in the case of networks, where the product of one reaction can act as an enzyme in another, has not been fully developed. Here we rigorously derive such an approximation in a class of coupled enzymatic networks where the individual interactions are of Michaelis-Menten type. We show that the sufficient conditions for the validity of the total quasi steady state assumption (tQSSA), obtained in a single protein case by Borghans, de Boer and Segel can be extended to sufficient conditions for the validity of the tQSSA in a large class of enzymatic networks. Secondly, we derive reduced equations that approximate the network's dynamics and involve only protein concentrations. This significantly reduces the number of equations necessary to model such systems. We prove the vali...
A network dynamics approach to chemical reaction networks
van der Schaft, A. J.; Rao, S.; Jayawardhana, B.
2016-04-01
A treatment of a chemical reaction network theory is given from the perspective of nonlinear network dynamics, in particular of consensus dynamics. By starting from the complex-balanced assumption, the reaction dynamics governed by mass action kinetics can be rewritten into a form which allows for a very simple derivation of a number of key results in the chemical reaction network theory, and which directly relates to the thermodynamics and port-Hamiltonian formulation of the system. Central in this formulation is the definition of a balanced Laplacian matrix on the graph of chemical complexes together with a resulting fundamental inequality. This immediately leads to the characterisation of the set of equilibria and their stability. Furthermore, the assumption of complex balancedness is revisited from the point of view of Kirchhoff's matrix tree theorem. Both the form of the dynamics and the deduced behaviour are very similar to consensus dynamics, and provide additional perspectives to the latter. Finally, using the classical idea of extending the graph of chemical complexes by a 'zero' complex, a complete steady-state stability analysis of mass action kinetics reaction networks with constant inflows and mass action kinetics outflows is given, and a unified framework is provided for structure-preserving model reduction of this important class of open reaction networks.
Deciphering Time Scale Hierarchy in Reaction Networks.
Nagahata, Yutaka; Maeda, Satoshi; Teramoto, Hiroshi; Horiyama, Takashi; Taketsugu, Tetsuya; Komatsuzaki, Tamiki
2016-03-01
Markovian dynamics on complex reaction networks are one of the most intriguing subjects in a wide range of research fields including chemical reactions, biological physics, and ecology. To represent the global kinetics from one node (corresponding to a basin on an energy landscape) to another requires information on multiple pathways that directly or indirectly connect these two nodes through the entire network. In this paper we present a scheme to extract a hierarchical set of global transition states (TSs) over a discrete-time Markov chain derived from first-order rate equations. The TSs can naturally take into account the multiple pathways connecting any pair of nodes. We also propose a new type of disconnectivity graph (DG) to capture the hierarchical organization of different time scales of reactions that can capture changes in the network due to changes in the time scale of observation. The crux is the introduction of the minimum conductance cut (MCC) in graph clustering, corresponding to the dividing surface across the network having the "smallest" transition probability between two disjoint subnetworks (superbasins on the energy landscape) in the network. We present a new combinatorial search algorithm for finding this MCC. We apply our method to a reaction network of Claisen rearrangement of allyl vinyl ether that consists of 23 nodes and 66 links (saddles on the energy landscape) connecting them. We compare the kinetic properties of our DG to those of the transition matrix of the rate equations and show that our graph can properly reveal the hierarchical organization of time scales in a network. PMID:26641663
Thomas Philipp
2012-05-01
Full Text Available Abstract Background It is well known that the deterministic dynamics of biochemical reaction networks can be more easily studied if timescale separation conditions are invoked (the quasi-steady-state assumption. In this case the deterministic dynamics of a large network of elementary reactions are well described by the dynamics of a smaller network of effective reactions. Each of the latter represents a group of elementary reactions in the large network and has associated with it an effective macroscopic rate law. A popular method to achieve model reduction in the presence of intrinsic noise consists of using the effective macroscopic rate laws to heuristically deduce effective probabilities for the effective reactions which then enables simulation via the stochastic simulation algorithm (SSA. The validity of this heuristic SSA method is a priori doubtful because the reaction probabilities for the SSA have only been rigorously derived from microscopic physics arguments for elementary reactions. Results We here obtain, by rigorous means and in closed-form, a reduced linear Langevin equation description of the stochastic dynamics of monostable biochemical networks in conditions characterized by small intrinsic noise and timescale separation. The slow-scale linear noise approximation (ssLNA, as the new method is called, is used to calculate the intrinsic noise statistics of enzyme and gene networks. The results agree very well with SSA simulations of the non-reduced network of elementary reactions. In contrast the conventional heuristic SSA is shown to overestimate the size of noise for Michaelis-Menten kinetics, considerably under-estimate the size of noise for Hill-type kinetics and in some cases even miss the prediction of noise-induced oscillations. Conclusions A new general method, the ssLNA, is derived and shown to correctly describe the statistics of intrinsic noise about the macroscopic concentrations under timescale separation conditions
An updated nuclear reaction network for BBN
Serpico, P. D.
2004-01-01
The key Standard-Physics inputs of the Big Bang Nucleosynthesis (BBN) are the light nuclei reaction rates. Both the network and the nuclear rates have been recently reanalyzed and updated, and cosmological and New-Physics constraints (taking into account the WMAP Cosmic Microwave Background anisotropies measurement) obtained with a new code are presented.
A feedback approach to bifurcation analysis in biochemical networks with many parameters
Waldherr, Steffen; Allgower, Frank
2007-01-01
Feedback circuits in biochemical networks which underly cellular signaling pathways are important elements in creating complex behavior. A specific aspect thereof is how stability of equilibrium points depends on model parameters. For biochemical networks, which are modelled using many parameters, it is typically very difficult to estimate the influence of parameters on stability. Finding parameters which result in a change in stability is a key step for a meaningful bifurcation analysis. We ...
Radical-ion-pair reactions are the biochemical equivalent of the optical double-slit experiment.
Kominis, Iannis K
2011-05-01
Radical-ion-pair reactions were recently shown to represent a rich biophysical laboratory for the application of quantum measurement theory methods and concepts. Here we show that radical-ion-pair reactions essentially form a nonlinear biochemical double-slit interferometer. Quantum coherence effects are visible when "which-path" information is limited, and the incoherent limit is approached when measurement-induced decoherence sets in. Based on this analogy with the optical double-slit experiment we derive and elaborate on the fundamental master equation of spin-selective radical-ion-pair reactions that covers the continuous range from complete incoherence to maximum singlet-triplet coherence. PMID:21728616
Insights into the organization of biochemical regulatory networks using graph theory analyses.
Ma'ayan, Avi
2009-02-27
Graph theory has been a valuable mathematical modeling tool to gain insights into the topological organization of biochemical networks. There are two types of insights that may be obtained by graph theory analyses. The first provides an overview of the global organization of biochemical networks; the second uses prior knowledge to place results from multivariate experiments, such as microarray data sets, in the context of known pathways and networks to infer regulation. Using graph analyses, biochemical networks are found to be scale-free and small-world, indicating that these networks contain hubs, which are proteins that interact with many other molecules. These hubs may interact with many different types of proteins at the same time and location or at different times and locations, resulting in diverse biological responses. Groups of components in networks are organized in recurring patterns termed network motifs such as feedback and feed-forward loops. Graph analysis revealed that negative feedback loops are less common and are present mostly in proximity to the membrane, whereas positive feedback loops are highly nested in an architecture that promotes dynamical stability. Cell signaling networks have multiple pathways from some input receptors and few from others. Such topology is reminiscent of a classification system. Signaling networks display a bow-tie structure indicative of funneling information from extracellular signals and then dispatching information from a few specific central intracellular signaling nexuses. These insights show that graph theory is a valuable tool for gaining an understanding of global regulatory features of biochemical networks. PMID:18940806
Combining flux and energy balance analysis to model large-scale biochemical networks.
Heuett, William J; Qian, Hong
2006-12-01
Stoichiometric Network Theory is a constraints-based, optimization approach for quantitative analysis of the phenotypes of large-scale biochemical networks that avoids the use of detailed kinetics. This approach uses the reaction stoichiometric matrix in conjunction with constraints provided by flux balance and energy balance to guarantee mass conserved and thermodynamically allowable predictions. However, the flux and energy balance constraints have not been effectively applied simultaneously on the genome scale because optimization under the combined constraints is non-linear. In this paper, a sequential quadratic programming algorithm that solves the non-linear optimization problem is introduced. A simple example and the system of fermentation in Saccharomyces cerevisiae are used to illustrate the new method. The algorithm allows the use of non-linear objective functions. As a result, we suggest a novel optimization with respect to the heat dissipation rate of a system. We also emphasize the importance of incorporating interactions between a model network and its surroundings. PMID:17245812
Modified Step Variational Iteration Method for Solving Fractional Biochemical Reaction Model
R. Yulita Molliq
2011-01-01
Full Text Available A new method called the modification of step variational iteration method (MoSVIM is introduced and used to solve the fractional biochemical reaction model. The MoSVIM uses general Lagrange multipliers for construction of the correction functional for the problems, and it runs by step approach, which is to divide the interval into subintervals with time step, and the solutions are obtained at each subinterval as well adopting a nonzero auxiliary parameter ℏ to control the convergence region of series' solutions. The MoSVIM yields an analytical solution of a rapidly convergent infinite power series with easily computable terms and produces a good approximate solution on enlarged intervals for solving the fractional biochemical reaction model. The accuracy of the results obtained is in a excellent agreement with the Adam Bashforth Moulton method (ABMM.
A new computational method to split large biochemical networks into coherent subnets
Verwoerd Wynand S
2011-02-01
Full Text Available Abstract Background Compared to more general networks, biochemical networks have some special features: while generally sparse, there are a small number of highly connected metabolite nodes; and metabolite nodes can also be divided into two classes: internal nodes with associated mass balance constraints and external ones without. Based on these features, reclassifying selected internal nodes (separators to external ones can be used to divide a large complex metabolic network into simpler subnetworks. Selection of separators based on node connectivity is commonly used but affords little detailed control and tends to produce excessive fragmentation. The method proposed here (Netsplitter allows the user to control separator selection. It combines local connection degree partitioning with global connectivity derived from random walks on the network, to produce a more even distribution of subnetwork sizes. Partitioning is performed progressively and the interactive visual matrix presentation used allows the user considerable control over the process, while incorporating special strategies to maintain the network integrity and minimise the information loss due to partitioning. Results Partitioning of a genome scale network of 1348 metabolites and 1468 reactions for Arabidopsis thaliana encapsulates 66% of the network into 10 medium sized subnets. Applied to the flavonoid subnetwork extracted in this way, it is shown that Netsplitter separates this naturally into four subnets with recognisable functionality, namely synthesis of lignin precursors, flavonoids, coumarin and benzenoids. A quantitative quality measure called efficacy is constructed and shows that the new method gives improved partitioning for several metabolic networks, including bacterial, plant and mammal species. Conclusions For the examples studied the Netsplitter method is a considerable improvement on the performance of connection degree partitioning, giving a better balance of
Switching dynamics in reaction networks induced by molecular discreteness
To study the fluctuations and dynamics in chemical reaction processes, stochastic differential equations based on the rate equation involving chemical concentrations are often adopted. When the number of molecules is very small, however, the discreteness in the number of molecules cannot be neglected since the number of molecules must be an integer. This discreteness can be important in biochemical reactions, where the total number of molecules is not significantly larger than the number of chemical species. To elucidate the effects of such discreteness, we study autocatalytic reaction systems comprising several chemical species through stochastic particle simulations. The generation of novel states is observed; it is caused by the extinction of some molecular species due to the discreteness in their number. We demonstrate that the reaction dynamics are switched by a single molecule, which leads to the reconstruction of the acting network structure. We also show the strong dependence of the chemical concentrations on the system size, which is caused by transitions to discreteness-induced novel states
Performance Improvements for Nuclear Reaction Network Integration
Longland, Richard; José, Jordi
2014-01-01
Aims: The aim of this work is to compare the performance of three reaction network integration methods used in stellar nucleosynthesis calculations. These are the Gear's backward differentiation method, Wagoner's method (a 2nd-order Runge-Kutta method), and the Bader-Deuflehard semi-implicit multi-step method. Methods: To investigate the efficiency of each of the integration methods considered here, a test suite of temperature and density versus time profiles is used. This suite provides a range of situations ranging from constant temperature and density to the dramatically varying conditions present in white dwarf mergers, novae, and x-ray bursts. Some of these profiles are obtained separately from full hydrodynamic calculations. The integration efficiencies are investigated with respect to input parameters that constrain the desired accuracy and precision. Results: Gear's backward differentiation method is found to improve accuracy, performance, and stability in integrating nuclear reaction networks. For te...
Nielsen, Jens Bredal
1997-01-01
of the thermokinetic description of reaction rates to include the influence of effecters. Here the reaction rate is written as a linear function of the logarithm of the metabolite concentrations. With this type of rate function it is shown that the approach of Delgado and Liao [Biochem. J. (1992) 282......, 919-927] can be much more widely applied, although it was originally based on linearized kinetics. The methodology of determining elasticity coefficients directly from pool levels is illustrated with an analysis of the first two steps of the biosynthetic pathway of penicillin. The results compare well...
Neural Networks in Chemical Reaction Dynamics
Raff, Lionel; Hagan, Martin
2011-01-01
This monograph presents recent advances in neural network (NN) approaches and applications to chemical reaction dynamics. Topics covered include: (i) the development of ab initio potential-energy surfaces (PES) for complex multichannel systems using modified novelty sampling and feedforward NNs; (ii) methods for sampling the configuration space of critical importance, such as trajectory and novelty sampling methods and gradient fitting methods; (iii) parametrization of interatomic potential functions using a genetic algorithm accelerated with a NN; (iv) parametrization of analytic interatomic
Ubiquitous "glassy" relaxation in catalytic reaction networks
Awazu, Akinori; Kaneko, Kunihiko
2009-01-01
Study of reversible catalytic reaction networks is important not only as an issue for chemical thermodynamics but also for protocells. From extensive numerical simulations and theoretical analysis, slow relaxation dynamics to sustain nonequlibrium states are commonly observed. These dynamics show two types of salient behaviors that are reminiscent of glassy behavior: slow relaxation along with the logarithmic time dependence of the correlation function and the emergence of plateaus in the rel...
A New Approach for Modeling Procaryotic Biochemical Networks With Differential Equations
Gebert, Jutta; Radde, Nicole
2006-06-01
One major challenge in Computational Biology is the simulation of the processes in a biological cell, which makes it necessary to understand the interactions between cell components. It is convenient to model the entirety of such interactions as biochemical networks. In this paper we present our novel approach to describe these biochemical networks with piecewise linear differential equations and analyze it theoretically. Then we will discuss methods for the parameter estimation from time series measurements including inference of the network topology. Finally we show an application of our model for the bacterium Corynebacterium glutamicum.
From Catalytic Reaction Networks to Protocells
Kaneko, Kunihiko
2013-12-01
In spite of recent advances, there still remains a large gape between a set of chemical reactions and a biological cell. Here we discuss several theoretical efforts to fill in the gap. The topics cover (i) slow relaxation to equilibrium due to glassy behavior in catalytic reaction networks (ii) consistency between molecule replication and cell growth, as well as energy metabolism (iii) control of a system by minority molecules in mutually catalytic system, which work as a carrier of genetic information, and leading to evolvability (iv) generation of a compartmentalized structure as a cluster of molecules centered around the minority molecule, and division of the cluster accompanied by the replication of minority molecule (v) sequential, logical process over several states from concurrent reaction dynamics, by taking advantage of discreteness in molecule number.
Graphical reduction of reaction networks by linear elimination of species
Sáez, Meritxell; Wiuf, Carsten; Feliu, Elisenda
2015-01-01
We give a graphically based procedure to reduce a reaction network to a smaller reaction network with fewer species after linear elimination of a set of noninteracting species. We give a description of the reduced reaction network, its kinetics and conservations laws, and explore properties of the...
Johnston, Matthew D; Pantea, Casian; Donnell, Pete
2016-01-01
We introduce a mixed-integer linear programming (MILP) framework capable of determining whether a chemical reaction network possesses the property of being endotactic or strongly endotactic. The network property of being strongly endotactic is known to lead to persistence and permanence of chemical species under genetic kinetic assumptions, while the same result is conjectured but as yet unproved for general endotactic networks. The algorithms we present are the first capable of verifying endotacticity of chemical reaction networks for systems with greater than two constituent species. We implement the algorithms in the open-source online package CoNtRol and apply them to a large sample of networks from the European Bioinformatics Institute's BioModels Database. We use strong endotacticity to establish for the first time the permanence of a well-studied circadian clock mechanism. PMID:25986743
Breitling, Rainer; Gilbert, David; Heiner, Monika; Orton, Richard
2008-01-01
Quantitative models of biochemical networks (signal transduction cascades, metabolic pathways, gene regulatory circuits) are a central component of modern systems biology. Building and managing these complex models is a major challenge that can benefit from the application of formal methods adopted from theoretical computing science. Here we provide a general introduction to the field of formal modelling, which emphasizes the intuitive biochemical basis of the modelling process, but is also a...
Discreteness-Induced Criticality in Random Catalytic Reaction Networks
Awazu, Akinori; Kaneko, Kunihiko
2009-01-01
Universal intermittent dynamics in a random catalytic reaction network, induced by smallness in the molecule number is reported. Stochastic simulations for a random catalytic reaction network subject to a flow of chemicals show that the system undergoes a transition from a stationary to an intermittent reaction phase when the flow rate is decreased. In the intermittent reaction phase, two temporal regimes with active and halted reactions alternate. The number frequency of reaction events at e...
Hahnel, G.B.; Gould, R.W.
1982-01-01
Incubation temperatures of 11°, 18° and 28° did not substantially affect biochemical reactions of either virulent or avirulent forms of Aeromonas salmonicida subspecies salmonicida. The only change observed, amygdalin fermentation, was positive at 11° and 18° but negative at 28°C. Several isolates utilized sucrose, a characteristic not normally recognized for A. salmonicida subspecies salmonicida.Antimicrobial susceptibility screening indicated resistance to novobiocin increased at the higher incubation temperatures. Standardized drug sensitivity testing procedures and precise zone diameter interpretive standards for bacterial fish pathogens are needed.
Introduction to dynamic spin chemistry magnetic field effects on chemical and biochemical reactions
Hayashi, Hisaharu
2004-01-01
This book presents a detailed account of one of the most mysterious problems in science - whether ordinary magnetic fields can exert an appreciable influence on chemical and biochemical reactions. The first aim of the book is to introduce this research, through theoretical and dynamic spin chemistry, to graduate students and researchers, by means of detailed theoretical and experimental descriptions. The second aim is to review typical recent investigations, which will stimulate new interest and applications in the 21st century. Because dynamic spin chemistry is based on established science, i
Computational analysis of the roles of biochemical reactions in anomalous diffusion dynamics
Naruemon, Rueangkham; Charin, Modchang
2016-04-01
Most biochemical processes in cells are usually modeled by reaction–diffusion (RD) equations. In these RD models, the diffusive process is assumed to be Gaussian. However, a growing number of studies have noted that intracellular diffusion is anomalous at some or all times, which may result from a crowded environment and chemical kinetics. This work aims to computationally study the effects of chemical reactions on the diffusive dynamics of RD systems by using both stochastic and deterministic algorithms. Numerical method to estimate the mean-square displacement (MSD) from a deterministic algorithm is also investigated. Our computational results show that anomalous diffusion can be solely due to chemical reactions. The chemical reactions alone can cause anomalous sub-diffusion in the RD system at some or all times. The time-dependent anomalous diffusion exponent is found to depend on many parameters, including chemical reaction rates, reaction orders, and chemical concentrations. Project supported by the Thailand Research Fund and Mahidol University (Grant No. TRG5880157), the Thailand Center of Excellence in Physics (ThEP), CHE, Thailand, and the Development Promotion of Science and Technology.
Proton mediated control of biochemical reactions with bioelectronic pH modulation
Deng, Yingxin; Miyake, Takeo; Keene, Scott; Josberger, Erik E.; Rolandi, Marco
2016-04-01
In Nature, protons (H+) can mediate metabolic process through enzymatic reactions. Examples include glucose oxidation with glucose dehydrogenase to regulate blood glucose level, alcohol dissolution into carboxylic acid through alcohol dehydrogenase, and voltage-regulated H+ channels activating bioluminescence in firefly and jellyfish. Artificial devices that control H+ currents and H+ concentration (pH) are able to actively influence biochemical processes. Here, we demonstrate a biotransducer that monitors and actively regulates pH-responsive enzymatic reactions by monitoring and controlling the flow of H+ between PdHx contacts and solution. The present transducer records bistable pH modulation from an “enzymatic flip-flop” circuit that comprises glucose dehydrogenase and alcohol dehydrogenase. The transducer also controls bioluminescence from firefly luciferase by affecting solution pH.
Proton mediated control of biochemical reactions with bioelectronic pH modulation
Deng, Yingxin; Miyake, Takeo; Keene, Scott; Josberger, Erik E.; Rolandi, Marco
2016-01-01
In Nature, protons (H+) can mediate metabolic process through enzymatic reactions. Examples include glucose oxidation with glucose dehydrogenase to regulate blood glucose level, alcohol dissolution into carboxylic acid through alcohol dehydrogenase, and voltage-regulated H+ channels activating bioluminescence in firefly and jellyfish. Artificial devices that control H+ currents and H+ concentration (pH) are able to actively influence biochemical processes. Here, we demonstrate a biotransducer that monitors and actively regulates pH-responsive enzymatic reactions by monitoring and controlling the flow of H+ between PdHx contacts and solution. The present transducer records bistable pH modulation from an “enzymatic flip-flop” circuit that comprises glucose dehydrogenase and alcohol dehydrogenase. The transducer also controls bioluminescence from firefly luciferase by affecting solution pH. PMID:27052724
Open complex-balanced mass action chemical reaction networks
Rao, Shodhan; van der Schaft, Arjan; Jayawardhana, Bayu
2014-01-01
We consider open chemical reaction networks, i.e. ones with inflows and outflows. We assume that all the inflows to the network are constant and all outflows obey the mass action kinetics rate law. We define a complex-balanced open reaction network as one that admits a complex-balanced steady state.
SBMLsqueezer: A CellDesigner plug-in to generate kinetic rate equations for biochemical networks
Schröder Adrian
2008-04-01
Full Text Available Abstract Background The development of complex biochemical models has been facilitated through the standardization of machine-readable representations like SBML (Systems Biology Markup Language. This effort is accompanied by the ongoing development of the human-readable diagrammatic representation SBGN (Systems Biology Graphical Notation. The graphical SBML editor CellDesigner allows direct translation of SBGN into SBML, and vice versa. For the assignment of kinetic rate laws, however, this process is not straightforward, as it often requires manual assembly and specific knowledge of kinetic equations. Results SBMLsqueezer facilitates exactly this modeling step via automated equation generation, overcoming the highly error-prone and cumbersome process of manually assigning kinetic equations. For each reaction the kinetic equation is derived from the stoichiometry, the participating species (e.g., proteins, mRNA or simple molecules as well as the regulatory relations (activation, inhibition or other modulations of the SBGN diagram. Such information allows distinctions between, for example, translation, phosphorylation or state transitions. The types of kinetics considered are numerous, for instance generalized mass-action, Hill, convenience and several Michaelis-Menten-based kinetics, each including activation and inhibition. These kinetics allow SBMLsqueezer to cover metabolic, gene regulatory, signal transduction and mixed networks. Whenever multiple kinetics are applicable to one reaction, parameter settings allow for user-defined specifications. After invoking SBMLsqueezer, the kinetic formulas are generated and assigned to the model, which can then be simulated in CellDesigner or with external ODE solvers. Furthermore, the equations can be exported to SBML, LaTeX or plain text format. Conclusion SBMLsqueezer considers the annotation of all participating reactants, products and regulators when generating rate laws for reactions. Thus, for
On the rejection-based algorithm for simulation and analysis of large-scale reaction networks
Stochastic simulation for in silico studies of large biochemical networks requires a great amount of computational time. We recently proposed a new exact simulation algorithm, called the rejection-based stochastic simulation algorithm (RSSA) [Thanh et al., J. Chem. Phys. 141(13), 134116 (2014)], to improve simulation performance by postponing and collapsing as much as possible the propensity updates. In this paper, we analyze the performance of this algorithm in detail, and improve it for simulating large-scale biochemical reaction networks. We also present a new algorithm, called simultaneous RSSA (SRSSA), which generates many independent trajectories simultaneously for the analysis of the biochemical behavior. SRSSA improves simulation performance by utilizing a single data structure across simulations to select reaction firings and forming trajectories. The memory requirement for building and storing the data structure is thus independent of the number of trajectories. The updating of the data structure when needed is performed collectively in a single operation across the simulations. The trajectories generated by SRSSA are exact and independent of each other by exploiting the rejection-based mechanism. We test our new improvement on real biological systems with a wide range of reaction networks to demonstrate its applicability and efficiency
On the rejection-based algorithm for simulation and analysis of large-scale reaction networks
Thanh, Vo Hong; Zunino, Roberto; Priami, Corrado
2015-06-01
Stochastic simulation for in silico studies of large biochemical networks requires a great amount of computational time. We recently proposed a new exact simulation algorithm, called the rejection-based stochastic simulation algorithm (RSSA) [Thanh et al., J. Chem. Phys. 141(13), 134116 (2014)], to improve simulation performance by postponing and collapsing as much as possible the propensity updates. In this paper, we analyze the performance of this algorithm in detail, and improve it for simulating large-scale biochemical reaction networks. We also present a new algorithm, called simultaneous RSSA (SRSSA), which generates many independent trajectories simultaneously for the analysis of the biochemical behavior. SRSSA improves simulation performance by utilizing a single data structure across simulations to select reaction firings and forming trajectories. The memory requirement for building and storing the data structure is thus independent of the number of trajectories. The updating of the data structure when needed is performed collectively in a single operation across the simulations. The trajectories generated by SRSSA are exact and independent of each other by exploiting the rejection-based mechanism. We test our new improvement on real biological systems with a wide range of reaction networks to demonstrate its applicability and efficiency.
Chang Yu-Te
2008-11-01
Full Text Available Abstract Background Gene networks in nanoscale are of nonlinear stochastic process. Time delays are common and substantial in these biochemical processes due to gene transcription, translation, posttranslation protein modification and diffusion. Molecular noises in gene networks come from intrinsic fluctuations, transmitted noise from upstream genes, and the global noise affecting all genes. Knowledge of molecular noise filtering and biochemical process delay compensation in gene networks is crucial to understand the signal processing in gene networks and the design of noise-tolerant and delay-robust gene circuits for synthetic biology. Results A nonlinear stochastic dynamic model with multiple time delays is proposed for describing a gene network under process delays, intrinsic molecular fluctuations, and extrinsic molecular noises. Then, the stochastic biochemical processing scheme of gene regulatory networks for attenuating these molecular noises and compensating process delays is investigated from the nonlinear signal processing perspective. In order to improve the robust stability for delay toleration and noise filtering, a robust gene circuit for nonlinear stochastic time-delay gene networks is engineered based on the nonlinear robust H∞ stochastic filtering scheme. Further, in order to avoid solving these complicated noise-tolerant and delay-robust design problems, based on Takagi-Sugeno (T-S fuzzy time-delay model and linear matrix inequalities (LMIs technique, a systematic gene circuit design method is proposed to simplify the design procedure. Conclusion The proposed gene circuit design method has much potential for application to systems biology, synthetic biology and drug design when a gene regulatory network has to be designed for improving its robust stability and filtering ability of disease-perturbed gene network or when a synthetic gene network needs to perform robustly under process delays and molecular noises.
New techniques for integrating nuclear reaction networks
The development of high resolution stellar models in multiple dimensions has moved us one step closer to answering fundamental questions of how stars burn their fuel. These models rely on solving a computationally expensive set of conservation, transportation, and energy production equations, the latter of which is determined by a complex network of nuclear reactions. The size of this network is constrained by the nature of the stellar environment in question and whether detailed nucleosynthesis yields are required. Alternatively, the nucleosynthesis yields can be computed by decoupling the calculation into a post-processing nucleosynthesis part, which relies on the output of a simplified hydrodynamics-focused calculation. Here, we investigate improving post-processing calculation performance by using more advanced integration methods than those typically adopted in nuclear astrophysics research. These are the Bader-Deuflehard multi-step method and the Gear backward differentiation method, which are compared with the traditional Wagoner two-step method. Integration method performance is quantified for a range of stellar environment models including Novae, X-ray bursts, hydrostatic core helium burning, and explosive hydrogen-helium burning in white dwarf merger events. We show that by using these methods, nucleosynthesis integration accuracy can be greatly improved and computation time can be significantly reduced by up to two orders of magnitude. Applicability of the integration methods to full hydrodynamical models will also be outlined. (author)
Efficient Characterization of Parametric Uncertainty of Complex (Biochemical Networks.
Claudia Schillings
2015-08-01
Full Text Available Parametric uncertainty is a particularly challenging and relevant aspect of systems analysis in domains such as systems biology where, both for inference and for assessing prediction uncertainties, it is essential to characterize the system behavior globally in the parameter space. However, current methods based on local approximations or on Monte-Carlo sampling cope only insufficiently with high-dimensional parameter spaces associated with complex network models. Here, we propose an alternative deterministic methodology that relies on sparse polynomial approximations. We propose a deterministic computational interpolation scheme which identifies most significant expansion coefficients adaptively. We present its performance in kinetic model equations from computational systems biology with several hundred parameters and state variables, leading to numerical approximations of the parametric solution on the entire parameter space. The scheme is based on adaptive Smolyak interpolation of the parametric solution at judiciously and adaptively chosen points in parameter space. As Monte-Carlo sampling, it is "non-intrusive" and well-suited for massively parallel implementation, but affords higher convergence rates. This opens up new avenues for large-scale dynamic network analysis by enabling scaling for many applications, including parameter estimation, uncertainty quantification, and systems design.
Hysteresis-driven structure formation in biochemical networks
Klein
1998-09-21
A mechanism of structure formation, based on hysteresis behaviour is presented. A bisubstrate kinetic system with substrate inhibition, discussed previously in the context of Turing structure formation, may show hysteresis behaviour, when embedded in a metabolic network: the system may possess multiple steady states and may be switched from one stable fixpoint to the other. When cells containing this type of system are diffusively coupled, under certain conditions patterns result, which, as is demonstrated, are not of the Turing type. The main difference to diffusion-driven (Turing) structures is the fact that the hysteresis-driven patterns emerge under diffusive conditions, under which both the homogeneous and the asymmetrical steady state is stable. The resulting special properties and biological implications are discussed.Copyright 1998 Academic Press Limited PMID:9778438
Inoue, Kentaro; Maeda, Kazuhiro; Miyabe, Takaaki; Matsuoka, Yu; Kurata, Hiroyuki
2014-09-01
Mathematical modeling has become a standard technique to understand the dynamics of complex biochemical systems. To promote the modeling, we had developed the CADLIVE dynamic simulator that automatically converted a biochemical map into its associated mathematical model, simulated its dynamic behaviors and analyzed its robustness. To enhance the feasibility by CADLIVE and extend its functions, we propose the CADLIVE toolbox available for MATLAB, which implements not only the existing functions of the CADLIVE dynamic simulator, but also the latest tools including global parameter search methods with robustness analysis. The seamless, bottom-up processes consisting of biochemical network construction, automatic construction of its dynamic model, simulation, optimization, and S-system analysis greatly facilitate dynamic modeling, contributing to the research of systems biology and synthetic biology. This application can be freely downloaded from http://www.cadlive.jp/CADLIVE_MATLAB/ together with an instruction. PMID:24623466
Open chemical reaction networks, steady-state loads and Braess-like paradox
Banerjee, Kinshuk
2014-01-01
Open chemical reaction systems involve matter-exchange with the surroundings. As a result, species can accumulate inside a system during the course of the reaction. We study the role of network topology in governing the concentration build-up inside a fixed reaction volume at steady state, particularly focusing on the effect of additional paths. The problem is akin to that in traffic networks where an extra route, surprisingly, can increase the overall travel time. This is known as the Braess' paradox. Here, we report chemical analogues of such a paradox in suitably chosen reaction networks, where extra reaction step(s) can inflate the total concentration, denoted as `load', at steady state. It is shown that, such counter-intuitive behavior emerges in a qualitatively similar pattern in networks of varying complexities. We then explore how such extra routes affect the load in a biochemical scheme of uric acid degradation. From a thorough analysis of this network, we propose a functional role of some decomposit...
Deterministic Function Computation with Chemical Reaction Networks
Chen, Ho-Lin; Soloveichik, David
2012-01-01
We study the deterministic computation of functions on tuples of natural numbers by chemical reaction networks (CRNs). CRNs have been shown to be efficiently Turing-universal when allowing for a small probability of error. CRNs that are guaranteed to converge on a correct answer, on the other hand, have been shown to decide only the semilinear predicates. We introduce the notion of function, rather than predicate, computation by representing the output of a function f:N^k --> N^l by a count of some molecular species, i.e., if the CRN starts with n_1,...,n_k molecules of some "input" species X_1,...,X_k, the CRN is guaranteed to converge to having f(n_1,...,n_k) molecules of the "output" species Y_1,...,Y_l. We show that a function f:N^k --> N^l is deterministically computed by a CRN if and only if its graph {(x,y) \\in N^k x N^l | f(x) = y} is a semilinear set. Finally, we show that each semilinear function f can be computed on input x in expected time O(polylog |x|).
Accelerated Sensitivity Analysis in High-Dimensional Stochastic Reaction Networks
Arampatzis, Georgios; Katsoulakis, Markos A.; Pantazis, Yannis
2015-01-01
Existing sensitivity analysis approaches are not able to handle efficiently stochastic reaction networks with a large number of parameters and species, which are typical in the modeling and simulation of complex biochemical phenomena. In this paper, a two-step strategy for parametric sensitivity analysis for such systems is proposed, exploiting advantages and synergies between two recently proposed sensitivity analysis methodologies for stochastic dynamics. The first method performs sensitivity analysis of the stochastic dynamics by means of the Fisher Information Matrix on the underlying distribution of the trajectories; the second method is a reduced-variance, finite-difference, gradient-type sensitivity approach relying on stochastic coupling techniques for variance reduction. Here we demonstrate that these two methods can be combined and deployed together by means of a new sensitivity bound which incorporates the variance of the quantity of interest as well as the Fisher Information Matrix estimated from the first method. The first step of the proposed strategy labels sensitivities using the bound and screens out the insensitive parameters in a controlled manner. In the second step of the proposed strategy, a finite-difference method is applied only for the sensitivity estimation of the (potentially) sensitive parameters that have not been screened out in the first step. Results on an epidermal growth factor network with fifty parameters and on a protein homeostasis with eighty parameters demonstrate that the proposed strategy is able to quickly discover and discard the insensitive parameters and in the remaining potentially sensitive parameters it accurately estimates the sensitivities. The new sensitivity strategy can be several times faster than current state-of-the-art approaches that test all parameters, especially in “sloppy” systems. In particular, the computational acceleration is quantified by the ratio between the total number of parameters over
Mélykúti, Bence
2010-01-01
The Chemical Langevin Equation (CLE), which is a stochastic differential equation driven by a multidimensional Wiener process, acts as a bridge between the discrete stochastic simulation algorithm and the deterministic reaction rate equation when simulating (bio)chemical kinetics. The CLE model is valid in the regime where molecular populations are abundant enough to assume their concentrations change continuously, but stochastic fluctuations still play a major role. The contribution of this work is that we observe and explore that the CLE is not a single equation, but a parametric family of equations, all of which give the same finite-dimensional distribution of the variables. On the theoretical side, we prove that as many Wiener processes are sufficient to formulate the CLE as there are independent variables in the equation, which is just the rank of the stoichiometric matrix. On the practical side, we show that in the case where there are m1 pairs of reversible reactions and m2 irreversible reactions there is another, simple formulation of the CLE with only m1 + m2 Wiener processes, whereas the standard approach uses 2 m1 + m2. We demonstrate that there are considerable computational savings when using this latter formulation. Such transformations of the CLE do not cause a loss of accuracy and are therefore distinct from model reduction techniques. We illustrate our findings by considering alternative formulations of the CLE for a human ether a-go-go related gene ion channel model and the Goldbeter-Koshland switch. © 2010 American Institute of Physics.
Mutual information in time-varying biochemical systems
Tostevin, Filipe; Wolde, Pieter Rein ten
2010-01-01
Cells must continuously sense and respond to time-varying environmental stimuli. These signals are transmitted and processed by biochemical signalling networks. However, the biochemical reactions making up these networks are intrinsically noisy, which limits the reliability of intracellular signalling. Here we use information theory to characterise the reliability of transmission of time-varying signals through elementary biochemical reactions in the presence of noise. We calculate the mutual...
On the Complexity of Reconstructing Chemical Reaction Networks
Fagerberg, Rolf; Flamm, Christoph; Merkle, Daniel;
2013-01-01
The analysis of the structure of chemical reaction networks is crucial for a better understanding of chemical processes. Such networks are well described as hypergraphs. However, due to the available methods, analyses regarding network properties are typically made on standard graphs derived from...... the full hypergraph description, e.g. on the so-called species and reaction graphs. However, a reconstruction of the underlying hypergraph from these graphs is not necessarily unique. In this paper, we address the problem of reconstructing a hypergraph from its species and reaction graph and show NP...
Paugam, A; Benchetrit, M; Fiacre, A; Tourte-Schaefer, C; Dupouy-Camet, J
1999-02-01
We compared the ability of four commercially available yeast identification systems for routine laboratory hospital use: Auxacolor (AUX) (Sanofi Diagnostics Pasteur, Marne-la-Coquette), Fungichrom I (FUC) and Fungifast I Twin (FUF) (International Microbio, Toulon), Api Candida (API) (bioMérieux, Lyon). These systems are based on obtaining a biochemical profile easily defined by colorimetric reactions. We tested 202 yeasts belonging to 19 species which were included or were not included in the manufacturer's data base of the identification systems. Without extra tests, for all the organisms tested, after 24 h of incubation, the percentage of organisms correctly identified was 48% for AUX, 75% for FUC, 77% for FUF and 81% for API. However, if we consider the ratio of the number of correct identifications without extra tests with the number of yeasts included in the manufacturers' data bases (sensitivity) the results increased to 61% for AUX, 81% for FUC, 91% for FUF and 83% for API. These systems are particularly well adapted to medical use, they are simple to set up, interpret, and have very good efficiency for the yeasts most commonly isolated in clinical specimens. The findings reported here indicate that the most favourable results were obtained with FUF and API systems. PMID:10200929
The nonlinear dynamics of biochemical reactions in a small-sized system on the order of a cell are stochastic. Assuming spatial homogeneity, the populations of n molecular species follow a multi-dimensional birth-and-death process on Zn. We introduce the Delbrück–Gillespie process, a continuous-time Markov jump process, whose Kolmogorov forward equation has been known as the chemical master equation, and whose stochastic trajectories can be computed via the Gillespie algorithm. Using simple models, we illustrate that a system of nonlinear ordinary differential equations on Rn emerges in the infinite system size limit. For finite system size, transitions among multiple attractors of the nonlinear dynamical system are rare events with exponentially long transit times. There is a separation of time scales between the deterministic ODEs and the stochastic Markov jumps between attractors. No diffusion process can provide a global representation that is accurate on both short and long time scales for the nonlinear, stochastic population dynamics. On the short time scale and near deterministic stable fixed points, Ornstein–Uhlenbeck Gaussian processes give linear stochastic dynamics that exhibit time-irreversible circular motion for open, driven chemical systems. Extending this individual stochastic behaviour-based nonlinear population theory of molecular species to other biological systems is discussed. (invited article)
Weber's Law in Autocatalytic Reaction Networks
Inoue, Masayo; Kaneko, Kunihiko
2011-01-01
Biological responses often obey Weber's law, according to which the magnitude of the response depends only on the fold change in the external input. In this study, we demonstrate that a system involving a simple autocatalytic reaction shows such response when a chemical is slowly synthesized by the reaction from a faster influx process. We also show that an autocatalytic reaction process occurring in series or in parallel can obey Weber's law with an oscillatory adaptive response. Considering...
S-system-based analysis of the robust properties common to many biochemical network models.
Matsuoka, Yu; Jahan, Nusrat; Kurata, Hiroyuki
2016-05-01
Robustness is a key feature to characterize the adaptation of organisms to changes in their internal and external environments. A broad range of kinetic or dynamic models of biochemical systems have been developed. Robustness analyses are attractive for exploring some common properties of many biochemical models. To reveal such features, we transform different types of mathematical equations into a standard or intelligible formula and use the multiple parameter sensitivity (MPS) to identify some factors critically responsible for the total robustness to many perturbations. The MPS would be determined by the top quarter of the highly sensitive parameters rather than the single parameter with the maximum sensitivity. The MPS did not show any correlation to the network size. The MPS is closely related to the standard deviation of the sensitivity profile. A decrease in the standard deviation enhanced the total robustness, which shows the hallmark of distributed robustness that many factors (pathways) involve the total robustness. PMID:26861555
Neural networks for the prediction organic chemistry reactions
Wei, Jennifer N; Aspuru-Guzik, Alán
2016-01-01
Reaction prediction remains one of the great challenges for organic chemistry. Solving this problem computationally requires the programming of a vast amount of knowledge and intuition of the rules of organic chemistry and the development of algorithms for their application. It is desirable to develop algorithms that, like humans, "learn" from being exposed to examples of the application of the rules of organic chemistry. In this work, we introduce a novel algorithm for predicting the products of organic chemistry reactions using machine learning to first identify the reaction type. In particular, we trained deep convolutional neural networks to predict the outcome of reactions based example reactions, using a new reaction fingerprint model. Due to the flexibility of neural networks, the system can attempt to predict reactions outside the domain where it was trained. We test this capability on problems from a popular organic chemistry textbook.
Jao, Tun; Schröter, Manuel; Chen, Chao-Long; Cheng, Yu-Fan; Lo, Chun-Yi Zac; Chou, Kun-Hsien; Patel, Ameera X; Lin, Wei-Che; Lin, Ching-Po; Bullmore, Edward T
2015-11-15
Functional properties of the brain may be associated with changes in complex brain networks. However, little is known about how properties of large-scale functional brain networks may be altered stepwise in patients with disturbance of consciousness, e.g., an encephalopathy. We used resting-state fMRI data on patients suffering from various degrees of hepatic encephalopathy (HE) to explore how topological and spatial network properties of functional brain networks changed at different cognitive and consciousness states. Severity of HE was measured clinically and by neuropsychological tests. Fifty-eight non-alcoholic liver cirrhosis patients and 62 normal controls were studied. Patients were subdivided into liver cirrhosis with no outstanding HE (NoHE, n=23), minimal HE with cognitive impairment only detectable by neuropsychological tests (MHE, n=28), and clinically overt HE (OHE, n=7). From the earliest stage, the NoHE, functional brain networks were progressively more random, less clustered, and less modular. Since the intermediate stage (MHE), increased ammonia level was accompanied by concomitant exponential decay of mean connectivity strength, especially in the primary cortical areas and midline brain structures. Finally, at the OHE stage, there were radical reorganization of the topological centrality-i.e., the relative importance-of the hubs and reorientation of functional connections between nodes. In summary, this study illustrated progressively greater abnormalities in functional brain network organization in patients with clinical and biochemical evidence of more severe hepatic encephalopathy. The early-than-expected brain network dysfunction in cirrhotic patients suggests that brain functional connectivity and network analysis may provide useful and complementary biomarkers for more aggressive and earlier intervention of hepatic encephalopathy. Moreover, the stepwise deterioration of functional brain networks in HE patients may suggest that hierarchical
Modelling biochemical networks with intrinsic time delays: a hybrid semi-parametric approach
Oliveira Rui
2010-09-01
Full Text Available Abstract Background This paper presents a method for modelling dynamical biochemical networks with intrinsic time delays. Since the fundamental mechanisms leading to such delays are many times unknown, non conventional modelling approaches become necessary. Herein, a hybrid semi-parametric identification methodology is proposed in which discrete time series are incorporated into fundamental material balance models. This integration results in hybrid delay differential equations which can be applied to identify unknown cellular dynamics. Results The proposed hybrid modelling methodology was evaluated using two case studies. The first of these deals with dynamic modelling of transcriptional factor A in mammalian cells. The protein transport from the cytosol to the nucleus introduced a delay that was accounted for by discrete time series formulation. The second case study focused on a simple network with distributed time delays that demonstrated that the discrete time delay formalism has broad applicability to both discrete and distributed delay problems. Conclusions Significantly better prediction qualities of the novel hybrid model were obtained when compared to dynamical structures without time delays, being the more distinctive the more significant the underlying system delay is. The identification of the system delays by studies of different discrete modelling delays was enabled by the proposed structure. Further, it was shown that the hybrid discrete delay methodology is not limited to discrete delay systems. The proposed method is a powerful tool to identify time delays in ill-defined biochemical networks.
A microfluidic platform for controlled biochemical stimulation of twin neuronal networks.
Biffi, Emilia; Piraino, Francesco; Pedrocchi, Alessandra; Fiore, Gianfranco B; Ferrigno, Giancarlo; Redaelli, Alberto; Menegon, Andrea; Rasponi, Marco
2012-06-01
Spatially and temporally resolved delivery of soluble factors is a key feature for pharmacological applications. In this framework, microfluidics coupled to multisite electrophysiology offers great advantages in neuropharmacology and toxicology. In this work, a microfluidic device for biochemical stimulation of neuronal networks was developed. A micro-chamber for cell culturing, previously developed and tested for long term neuronal growth by our group, was provided with a thin wall, which partially divided the cell culture region in two sub-compartments. The device was reversibly coupled to a flat micro electrode array and used to culture primary neurons in the same microenvironment. We demonstrated that the two fluidically connected compartments were able to originate two parallel neuronal networks with similar electrophysiological activity but functionally independent. Furthermore, the device allowed to connect the outlet port to a syringe pump and to transform the static culture chamber in a perfused one. At 14 days invitro, sub-networks were independently stimulated with a test molecule, tetrodotoxin, a neurotoxin known to block action potentials, by means of continuous delivery. Electrical activity recordings proved the ability of the device configuration to selectively stimulate each neuronal network individually. The proposed microfluidic approach represents an innovative methodology to perform biological, pharmacological, and electrophysiological experiments on neuronal networks. Indeed, it allows for controlled delivery of substances to cells, and it overcomes the limitations due to standard drug stimulation techniques. Finally, the twin network configuration reduces biological variability, which has important outcomes on pharmacological and drug screening. PMID:22655017
A computational framework for the automated construction of glycosylation reaction networks.
Gang Liu
Full Text Available Glycosylation is among the most common and complex post-translational modifications identified to date. It proceeds through the catalytic action of multiple enzyme families that include the glycosyltransferases that add monosaccharides to growing glycans, and glycosidases which remove sugar residues to trim glycans. The expression level and specificity of these enzymes, in part, regulate the glycan distribution or glycome of specific cell/tissue systems. Currently, there is no systematic method to describe the enzymes and cellular reaction networks that catalyze glycosylation. To address this limitation, we present a streamlined machine-readable definition for the glycosylating enzymes and additional methodologies to construct and analyze glycosylation reaction networks. In this computational framework, the enzyme class is systematically designed to store detailed specificity data such as enzymatic functional group, linkage and substrate specificity. The new classes and their associated functions enable both single-reaction inference and automated full network reconstruction, when given a list of reactants and/or products along with the enzymes present in the system. In addition, graph theory is used to support functions that map the connectivity between two or more species in a network, and that generate subset models to identify rate-limiting steps regulating glycan biosynthesis. Finally, this framework allows the synthesis of biochemical reaction networks using mass spectrometry (MS data. The features described above are illustrated using three case studies that examine: i O-linked glycan biosynthesis during the construction of functional selectin-ligands; ii automated N-linked glycosylation pathway construction; and iii the handling and analysis of glycomics based MS data. Overall, the new computational framework enables automated glycosylation network model construction and analysis by integrating knowledge of glycan structure and enzyme
Cheali, Peam; Gernaey, Krist; Sin, Gürkan
2013-01-01
This study presents the development of an expanded biorefinery processing network for producing biofuels that combines biochemical and thermochemical conversion platforms. The expanded network is coupled to a framework that uses a superstructure based optimization approach to generate and compare...... for identifying at early stage optimal biorefinery concept with respect to technical, economic and environmental criteria....
Cellular metabolic network analysis: discovering important reactions in Treponema pallidum.
Chen, Xueying; Zhao, Min; Qu, Hong
2015-01-01
T. pallidum, the syphilis-causing pathogen, performs very differently in metabolism compared with other bacterial pathogens. The desire for safe and effective vaccine of syphilis requests identification of important steps in T. pallidum's metabolism. Here, we apply Flux Balance Analysis to represent the reactions quantitatively. Thus, it is possible to cluster all reactions in T. pallidum. By calculating minimal cut sets and analyzing topological structure for the metabolic network of T. pallidum, critical reactions are identified. As a comparison, we also apply the analytical approaches to the metabolic network of H. pylori to find coregulated drug targets and unique drug targets for different microorganisms. Based on the clustering results, all reactions are further classified into various roles. Therefore, the general picture of their metabolic network is obtained and two types of reactions, both of which are involved in nucleic acid metabolism, are found to be essential for T. pallidum. It is also discovered that both hubs of reactions and the isolated reactions in purine and pyrimidine metabolisms play important roles in T. pallidum. These reactions could be potential drug targets for treating syphilis. PMID:26495292
Ye, Jianxiong; Feng, Enmin; Wang, Lei; Xiu, Zhilong; Sun, Yaqin
Glycerol bioconversion to 1,3-propanediol (1,3-PD) by Klebsiella pneumoniae (K. pneumoniae) can be characterized by an intricate network of interactions among biochemical fluxes, metabolic compounds, key enzymes and genetic regulatory. To date, there still exist some uncertain factors in this complex network because of the limitation in bio-techniques, especially in measuring techniques for intracellular substances. In this paper, among these uncertain factors, we aim to infer the transport mechanisms of glycerol and 1,3-PD across the cell membrane, which have received intensive interest in recent years. On the basis of different inferences of the transport mechanisms, we reconstruct various metabolic networks correspondingly and subsequently develop their dynamical systems (S-systems). To determine the most reasonable metabolic network from all possible ones, we establish a quantitative definition of biological robustness and undertake parameter identification and robustness analysis for each system. Numerical results show that it is most possible that both glycerol and 1,3-PD pass the cell membrane by active transport and passive diffusion.
SBMLsqueezer 2: context-sensitive creation of kinetic equations in biochemical networks
Draeger, Andreas; Zielinski, Daniel C.; Keller, Roland;
2015-01-01
simplified the network reconstruction process, but building kinetic models for these systems is still a manually intensive task. Appropriate kinetic equations, based upon reaction rate laws, must be constructed and parameterized for each reaction. The complex test-and-evaluation cycles that can be involved...... during kinetic model construction would thus benefit from automated methods for rate law assignment. Results: We present a high-throughput algorithm to automatically suggest and create suitable rate laws based upon reaction type according to several criteria. The criteria for choices made by the...... algorithm can be influenced in order to assign the desired type of rate law to each reaction. This algorithm is implemented in the software package SBMLsqueezer 2. In addition, this program contains an integrated connection to the kinetics database SABIO-RK to obtain experimentally-derived rate laws when...
Turing instability in reaction-diffusion models on complex networks
Ide, Yusuke; Izuhara, Hirofumi; Machida, Takuya
2016-09-01
In this paper, the Turing instability in reaction-diffusion models defined on complex networks is studied. Here, we focus on three types of models which generate complex networks, i.e. the Erdős-Rényi, the Watts-Strogatz, and the threshold network models. From analysis of the Laplacian matrices of graphs generated by these models, we numerically reveal that stable and unstable regions of a homogeneous steady state on the parameter space of two diffusion coefficients completely differ, depending on the network architecture. In addition, we theoretically discuss the stable and unstable regions in the cases of regular enhanced ring lattices which include regular circles, and networks generated by the threshold network model when the number of vertices is large enough.
Integration of metabolome data with metabolic networks reveals reporter reactions
Çakir, Tunahan; Patil, Kiran Raosaheb; Önsan, Zeynep Ilsen; Ülgen, Kutlu Özergin; Kirdar, Betül; Nielsen, Jens
2006-01-01
Interpreting quantitative metabolome data is a difficult task owing to the high connectivity in metabolic networks and inherent interdependency between enzymatic regulation, metabolite levels and fluxes. Here we present a hypothesis-driven algorithm for the integration of such data with metabolic...... network topology. The algorithm thus enables identification of reporter reactions, which are reactions where there are significant coordinated changes in the level of surrounding metabolites following environmental/genetic perturbations. Applicability of the algorithm is demonstrated by using data from...... is measured. By combining the results with transcriptome data, we further show that it is possible to infer whether the reactions are hierarchically or metabolically regulated. Hereby, the reported approach represents an attempt to map different layers of regulation within metabolic networks through...
生化反应的五行归属%Anfive elements classification of bio-chemical reactions
徐天成
2015-01-01
五行学说提供了对事物进行系统分类的合理方法。人类体内充满了复杂的生物化学反应体系，能否用五行学说的基本规律归纳生物化学反应体系值得探讨。在以往相关研究的基础上，本文利用五行思想对生化反应的物质和反应类型进行分类，用简洁的中医学规律研究复杂的生化反应过程，这种新思路对药物研究、中西医结合学科的发展等具有挖掘价值。%Thefive elements theory provides us a system classification of all things so we can do the research in a much reasonable way. Further more, our bodies are full of complicated bio-chemical reactions and thefive elements theory should also be applicable to those reaction systems. This paper will try to make use of thefive elements theory to carry on an innovative classification of different types of bio-chemical reactions with the help of the former foundation of related researches, to make the complicated bio-chemical reactions easy to understand, and point out that this kind of lately academic thought is of great significance to the medicine researches as well as the integrated traditional Chinese and western medicine.
Cellular Metabolic Network Analysis: Discovering Important Reactions in Treponema pallidum
Xueying Chen; Min Zhao; Hong Qu
2015-01-01
T. pallidum, the syphilis-causing pathogen, performs very differently in metabolism compared with other bacterial pathogens. The desire for safe and effective vaccine of syphilis requests identification of important steps in T. pallidum's metabolism. Here, we apply Flux Balance Analysis to represent the reactions quantitatively. Thus, it is possible to cluster all reactions in T. pallidum. By calculating minimal cut sets and analyzing topological structure for the metabolic network of T. pall...
Genetic programming-based approach to elucidate biochemical interaction networks from data.
Kandpal, Manoj; Kalyan, Chakravarthy Mynampati; Samavedham, Lakshminarayanan
2013-02-01
Biochemical systems are characterised by cyclic/reversible reciprocal actions, non-linear interactions and a mixed relationship structures (linear and non-linear; static and dynamic). Deciphering the architecture of such systems using measured data to provide quantitative information regarding the nature of relationships that exist between the measured variables is a challenging proposition. Causality detection is one of the methodologies that are applied to elucidate biochemical networks from such data. Autoregressive-based modelling approach such as granger causality, partial directed coherence, directed transfer function and canonical variate analysis have been applied on different systems for deciphering such interactions, but with limited success. In this study, the authors propose a genetic programming-based causality detection (GPCD) methodology which blends evolutionary computation-based procedures along with parameter estimation methods to derive a mathematical model of the system. Application of the GPCD methodology on five data sets that contained the different challenges mentioned above indicated that GPCD performs better than the other methods in uncovering the exact structure with less false positives. On a glycolysis data set, GPCD was able to fill the 'interaction gaps' which were missed by other methods. PMID:23848052
Casey, Fergal
2011-08-22
Abstract Background Gene and protein interactions are commonly represented as networks, with the genes or proteins comprising the nodes and the relationship between them as edges. Motifs, or small local configurations of edges and nodes that arise repeatedly, can be used to simplify the interpretation of networks. Results We examined triplet motifs in a network of quantitative epistatic genetic relationships, and found a non-random distribution of particular motif classes. Individual motif classes were found to be associated with different functional properties, suggestive of an underlying biological significance. These associations were apparent not only for motif classes, but for individual positions within the motifs. As expected, NNN (all negative) motifs were strongly associated with previously reported genetic (i.e. synthetic lethal) interactions, while PPP (all positive) motifs were associated with protein complexes. The two other motif classes (NNP: a positive interaction spanned by two negative interactions, and NPP: a negative spanned by two positives) showed very distinct functional associations, with physical interactions dominating for the former but alternative enrichments, typical of biochemical pathways, dominating for the latter. Conclusion We present a model showing how NNP motifs can be used to recognize supportive relationships between protein complexes, while NPP motifs often identify opposing or regulatory behaviour between a gene and an associated pathway. The ability to use motifs to point toward underlying biological organizational themes is likely to be increasingly important as more extensive epistasis mapping projects in higher organisms begin.
Jiang, Wen Zhi; Adamec, Jiri; Weeks, Donald P
2013-11-15
A simple, inexpensive microdistillation device is described for capturing methanol or formaldehyde as end products of biochemical reactions or in environmental samples. We demonstrate that the microdistillation protocol, coupled with the use of alcohol oxidase and the formaldehyde-sensitive reagent Purpald (4-amino-3-hydrazino-5-mercapto-1,2,4-triazole), serves as a quick and inexpensive alternative to chromatographic and mass spectrometer analyses for determining if formaldehyde or methanol is a product of reactions that contain substances that interfere with the Purpald reaction. These techniques were used to affirm formaldehyde as the end product of the dicamba monooxygenase-catalyzed O-demethylation of the herbicide dicamba (2-methoxy-3,6-dichlorobenzoic acid). PMID:23938775
ARWAR: A network approach for predicting Adverse Drug Reactions.
Rahmani, Hossein; Weiss, Gerhard; Méndez-Lucio, Oscar; Bender, Andreas
2016-01-01
Predicting novel drug side-effects, or Adverse Drug Reactions (ADRs), plays an important role in the drug discovery process. Existing methods consider mainly the chemical and biological characteristics of each drug individually, thereby neglecting information hidden in the relationships among drugs. Complementary to the existing individual methods, in this paper, we propose a novel network approach for ADR prediction that is called Augmented Random-WAlk with Restarts (ARWAR). ARWAR, first, applies an existing method to build a network of highly related drugs. Then, it augments the original drug network by adding new nodes and new edges to the network and finally, it applies Random Walks with Restarts to predict novel ADRs. Empirical results show that the ARWAR method presented here outperforms the existing network approach by 20% with respect to average Fmeasure. Furthermore, ARWAR is capable of generating novel hypotheses about drugs with respect to novel and biologically meaningful ADR. PMID:26638149
On the steady states of weakly reversible chemical reaction networks
Deng, Jian; Jones, Christopher; Feinberg, Martin; Nachman, Adrian
2011-01-01
A natural condition on the structure of the underlying chemical reaction network, namely weak reversibility, is shown to guarantee the existence of an equilibrium (steady state) in each positive stoichiometric compatibility class for the associated mass-action system. Furthermore, an index formula is given for the set of equilibria in a given stoichiometric compatibility class.
Injectivity, multiple zeros, and multistationarity in reaction networks
Feliu, Elisenda
2015-01-01
Polynomial dynamical systems are widely used to model and study real phenomena. In biochemistry, they are the preferred choice for modelling the concentration of chemical species in reaction networks with mass-action kinetics. These systems are typically parametrized by many (unknown) parameters....
Complex and detailed balancing of chemical reaction networks revisited
van der Schaft, Abraham; Rao, Shodhan; Jayawardhana, Bayu
2015-01-01
The characterization of the notions of complex and detailed balancing for mass action kinetics chemical reaction networks is revisited from the perspective of algebraic graph theory, in particular Kirchhoff’s Matrix Tree theorem for directed weighted graphs. This yields an elucidation of previously
Kentaro Inoue
Full Text Available BACKGROUND: For visualizing large-scale biochemical network maps, it is important to calculate the coordinates of molecular nodes quickly and to enhance the understanding or traceability of them. The grid layout is effective in drawing compact, orderly, balanced network maps with node label spaces, but existing grid layout algorithms often require a high computational cost because they have to consider complicated positional constraints through the entire optimization process. RESULTS: We propose a hybrid grid layout algorithm that consists of a non-grid, fast layout (preprocessor algorithm and an approximate pattern matching algorithm that distributes the resultant preprocessed nodes on square grid points. To demonstrate the feasibility of the hybrid layout algorithm, it is characterized in terms of the calculation time, numbers of edge-edge and node-edge crossings, relative edge lengths, and F-measures. The proposed algorithm achieves outstanding performances compared with other existing grid layouts. CONCLUSIONS: Use of an approximate pattern matching algorithm quickly redistributes the laid-out nodes by fast, non-grid algorithms on the square grid points, while preserving the topological relationships among the nodes. The proposed algorithm is a novel use of the pattern matching, thereby providing a breakthrough for grid layout. This application program can be freely downloaded from http://www.cadlive.jp/hybridlayout/hybridlayout.html.
Wu, Qianqian; Tian, Tianhai
2016-01-01
To deal with the growing scale of molecular systems, sophisticated modelling techniques have been designed in recent years to reduce the complexity of mathematical models. Among them, a widely used approach is delayed reaction for simplifying multistep reactions. However, recent research results suggest that a delayed reaction with constant time delay is unable to describe multistep reactions accurately. To address this issue, we propose a novel approach using state-dependent time delay to approximate multistep reactions. We first use stochastic simulations to calculate time delay arising from multistep reactions exactly. Then we design algorithms to calculate time delay based on system dynamics precisely. To demonstrate the power of proposed method, two processes of mRNA degradation are used to investigate the function of time delay in determining system dynamics. In addition, a multistep pathway of metabolic synthesis is used to explore the potential of the proposed method to simplify multistep reactions with nonlinear reaction rates. Simulation results suggest that the state-dependent time delay is a promising and accurate approach to reduce model complexity and decrease the number of unknown parameters in the models. PMID:27553753
David eGomez
2015-06-01
Full Text Available Molecular crowding is ubiquitous within cells and affects many biological processes including protein-protein binding, enzyme activities and gene regulation. Here we revisit some generic effects of crowding using a combination of lattice simulations and reaction-diffusion simulations with the program ReaDDy. Specifically, we implement three reactions, simple binding, a diffusion-limited reaction and a reaction with Michaelis-Menten kinetics. Histograms of binding and unbinding times provide a detailed picture how crowding affects these reactions and how the separate effects of crowding on binding equilibrium and on diffusion act together. In addition, we discuss how crowding affects processes related to gene expression such as RNA polymerase-promoter binding and translation elongation.
Gomez, David; Klumpp, Stefan
2015-06-01
Molecular crowding is ubiquitous within cells and affects many biological processes including protein-protein binding, enzyme activities and gene regulation. Here we revisit some generic effects of crowding using a combination of lattice simulations and reaction-diffusion simulations with the program ReaDDy. Specifically, we implement three reactions, simple binding, a diffusion-limited reaction and a reaction with Michaelis-Menten kinetics. Histograms of binding and unbinding times provide a detailed picture how crowding affects these reactions and how the separate effects of crowding on binding equilibrium and on diffusion act together. In addition, we discuss how crowding affects processes related to gene expression such as RNA polymerase-promoter binding and translation elongation.
Stochastic analysis of complex reaction networks using binomial moment equations.
Barzel, Baruch; Biham, Ofer
2012-09-01
The stochastic analysis of complex reaction networks is a difficult problem because the number of microscopic states in such systems increases exponentially with the number of reactive species. Direct integration of the master equation is thus infeasible and is most often replaced by Monte Carlo simulations. While Monte Carlo simulations are a highly effective tool, equation-based formulations are more amenable to analytical treatment and may provide deeper insight into the dynamics of the network. Here, we present a highly efficient equation-based method for the analysis of stochastic reaction networks. The method is based on the recently introduced binomial moment equations [Barzel and Biham, Phys. Rev. Lett. 106, 150602 (2011)]. The binomial moments are linear combinations of the ordinary moments of the probability distribution function of the population sizes of the interacting species. They capture the essential combinatorics of the reaction processes reflecting their stoichiometric structure. This leads to a simple and transparent form of the equations, and allows a highly efficient and surprisingly simple truncation scheme. Unlike ordinary moment equations, in which the inclusion of high order moments is prohibitively complicated, the binomial moment equations can be easily constructed up to any desired order. The result is a set of equations that enables the stochastic analysis of complex reaction networks under a broad range of conditions. The number of equations is dramatically reduced from the exponential proliferation of the master equation to a polynomial (and often quadratic) dependence on the number of reactive species in the binomial moment equations. The aim of this paper is twofold: to present a complete derivation of the binomial moment equations; to demonstrate the applicability of the moment equations for a representative set of example networks, in which stochastic effects play an important role. PMID:23030885
Urszula Małolepsza; Henryk Urbanek; Justyna Polit
2013-01-01
The reactions of strawberry plants to infection with B. cinerea and treatment with salicylic acid has been studied. Infection of leaves with B. cinerea resulted in early increases in active oxygen species generation, superoxide dismutase and peroxidase activities and phenolic compounds content. Some increases of the above reactions were noticed in plants treated with salicylic acid but not in the plants treated with SA and then later infected with B. cinerea.
Urszula Małolepsza
2013-12-01
Full Text Available The reactions of strawberry plants to infection with B. cinerea and treatment with salicylic acid has been studied. Infection of leaves with B. cinerea resulted in early increases in active oxygen species generation, superoxide dismutase and peroxidase activities and phenolic compounds content. Some increases of the above reactions were noticed in plants treated with salicylic acid but not in the plants treated with SA and then later infected with B. cinerea.
Complex Chemical Reaction Networks from Heuristics-Aided Quantum Chemistry
Rappoport, Dmitrij; Galvin, Cooper J.; Zubarev, Dmitry; Aspuru-Guzik, Alan
2014-01-01
While structures and reactivities of many small molecules can be computed efficiently and accurately using quantum chemical methods, heuristic approaches remain essential for modeling complex structures and large-scale chemical systems. Here, we present a heuristics-aided quantum chemical methodology applicable to complex chemical reaction networks such as those arising in cell metabolism and prebiotic chemistry. Chemical heuristics offer an expedient way of traversing high-dimensional reacti...
Variable elimination in chemical reaction networks with mass action kinetics
Feliu, Elisenda; Wiuf, Carsten
2011-01-01
We consider chemical reaction networks taken with mass action kinetics. The steady states of such a system are solutions to a system of polynomial equations. Even for small systems the task of finding the solutions is daunting. We develop an algebraic framework and procedure for linear elimination of variables. The procedure reduces the variables in the system to a set of "core" variables by eliminating variables corresponding to a set of non-interacting species. The steady states are paramet...
A chemical reaction network solver for the astrophysics code NIRVANA
Ziegler, U.
2016-02-01
Context. Chemistry often plays an important role in astrophysical gases. It regulates thermal properties by changing species abundances and via ionization processes. This way, time-dependent cooling mechanisms and other chemistry-related energy sources can have a profound influence on the dynamical evolution of an astrophysical system. Modeling those effects with the underlying chemical kinetics in realistic magneto-gasdynamical simulations provide the basis for a better link to observations. Aims: The present work describes the implementation of a chemical reaction network solver into the magneto-gasdynamical code NIRVANA. For this purpose a multispecies structure is installed, and a new module for evolving the rate equations of chemical kinetics is developed and coupled to the dynamical part of the code. A small chemical network for a hydrogen-helium plasma was constructed including associated thermal processes which is used in test problems. Methods: Evolving a chemical network within time-dependent simulations requires the additional solution of a set of coupled advection-reaction equations for species and gas temperature. Second-order Strang-splitting is used to separate the advection part from the reaction part. The ordinary differential equation (ODE) system representing the reaction part is solved with a fourth-order generalized Runge-Kutta method applicable for stiff systems inherent to astrochemistry. Results: A series of tests was performed in order to check the correctness of numerical and technical implementation. Tests include well-known stiff ODE problems from the mathematical literature in order to confirm accuracy properties of the solver used as well as problems combining gasdynamics and chemistry. Overall, very satisfactory results are achieved. Conclusions: The NIRVANA code is now ready to handle astrochemical processes in time-dependent simulations. An easy-to-use interface allows implementation of complex networks including thermal processes
Spreading out of perturbations in reversible reaction networks.
Maslov, Sergei; Sneppen, Kim; Ispolatov, I
2007-08-17
Using an example of physical interactions between proteins, we study how a perturbation propagates in the equilibrium of a network of reversible reactions governed by the law of mass action. We introduce a matrix formalism to describe the linear response of all equilibrium concentrations to shifts in total abundances of individual reactants, and reveal its heuristic analogy to the flow of electric current in a network of resistors. Our main conclusion is that, on average, the induced changes in equilibrium concentrations decay exponentially as a function of network distance from the source of perturbation. We analyze how this decay is influenced by such factors as the topology of a network, binding strength, and correlations between concentrations of neighboring nodes. We find that the minimal branching of the network, small values of dissociation constants, and low equilibrium free (unbound) concentrations of reacting substances all decrease the decay constant and thus increase the range of propagation. Exact analytic expressions for the decay constant are obtained for the case of equally strong interactions and uniform as well as oscillating concentrations on the Bethe lattice. Our general findings are illustrated using a real network of protein-protein interactions in baker's yeast with experimentally determined protein concentrations. PMID:18046464
Spreading out of perturbations in reversible reaction networks
Maslov, Sergei [Department of Condensed Matter Physics, Brookhaven National Laboratory, Upton, NY 11973 (United States); Sneppen, Kim [Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen Oe (Denmark); Ispolatov, I [Ariadne Genomics Inc., 9430 Key West Ave. 113 Rockville, MD 20850 (United States)
2007-08-15
Using an example of physical interactions between proteins, we study how a perturbation propagates in the equilibrium of a network of reversible reactions governed by the law of mass action. We introduce a matrix formalism to describe the linear response of all equilibrium concentrations to shifts in total abundances of individual reactants, and reveal its heuristic analogy to the flow of electric current in a network of resistors. Our main conclusion is that, on average, the induced changes in equilibrium concentrations decay exponentially as a function of network distance from the source of perturbation. We analyze how this decay is influenced by such factors as the topology of a network, binding strength, and correlations between concentrations of neighboring nodes. We find that the minimal branching of the network, small values of dissociation constants, and low equilibrium free (unbound) concentrations of reacting substances all decrease the decay constant and thus increase the range of propagation. Exact analytic expressions for the decay constant are obtained for the case of equally strong interactions and uniform as well as oscillating concentrations on the Bethe lattice. Our general findings are illustrated using a real network of protein-protein interactions in baker's yeast with experimentally determined protein concentrations.
Patterns of Stochastic Behavior in Dynamically Unstable High-Dimensional Biochemical Networks
Rosenfeld, Simon
2009-01-01
The question of dynamical stability and stochastic behavior of large biochemical networks is discussed. It is argued that stringent conditions of asymptotic stability have very little chance to materialize in a multidimensional system described by the differential equations of chemical kinetics. The reason is that the criteria of asymptotic stability (Routh-Hurwitz, Lyapunov criteria, Feinberg’s Deficiency Zero theorem) would impose the limitations of very high algebraic order on the kinetic rates and stoichiometric coefficients, and there are no natural laws that would guarantee their unconditional validity. Highly nonlinear, dynamically unstable systems, however, are not necessarily doomed to collapse, as a simple Jacobian analysis would suggest. It is possible that their dynamics may assume the form of pseudo-random fluctuations quite similar to a shot noise, and, therefore, their behavior may be described in terms of Langevin and Fokker-Plank equations. We have shown by simulation that the resulting pseudo-stochastic processes obey the heavy-tailed Generalized Pareto Distribution with temporal sequence of pulses forming the set of constituent-specific Poisson processes. Being applied to intracellular dynamics, these properties are naturally associated with burstiness, a well documented phenomenon in the biology of gene expression. PMID:19838330
Chen, Bor-Sen; Chen, Po-Wei
2009-12-01
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. PMID:19788895
Bosonic reaction-diffusion processes on scale-free networks
Baronchelli, Andrea; Catanzaro, Michele; Pastor-Satorras, Romualdo
2008-07-01
Reaction-diffusion processes can be adopted to model a large number of dynamics on complex networks, such as transport processes or epidemic outbreaks. In most cases, however, they have been studied from a fermionic perspective, in which each vertex can be occupied by at most one particle. While still useful, this approach suffers from some drawbacks, the most important probably being the difficulty to implement reactions involving more than two particles simultaneously. Here we develop a general framework for the study of bosonic reaction-diffusion processes on complex networks, in which there is no restriction on the number of interacting particles that a vertex can host. We describe these processes theoretically by means of continuous-time heterogeneous mean-field theory and divide them into two main classes: steady-state and monotonously decaying processes. We analyze specific examples of both behaviors within the class of one-species processes, comparing the results (whenever possible) with the corresponding fermionic counterparts. We find that the time evolution and critical properties of the particle density are independent of the fermionic or bosonic nature of the process, while differences exist in the functional form of the density of occupied vertices in a given degree class k . We implement a continuous-time Monte Carlo algorithm, well suited for general bosonic simulations, which allows us to confirm the analytical predictions formulated within mean-field theory. Our results, at both the theoretical and numerical levels, can be easily generalized to tackle more complex, multispecies, reaction-diffusion processes and open a promising path for a general study and classification of this kind of dynamical systems on complex networks.
Singh, Shalini; Samal, Areejit; Giri, Varun; Krishna, Sandeep; Raghuram, Nandula; Jain, Sanjay
2013-05-01
Unraveling the structure of complex biological networks and relating it to their functional role is an important task in systems biology. Here we attempt to characterize the functional organization of the large-scale metabolic networks of three microorganisms. We apply flux balance analysis to study the optimal growth states of these organisms in different environments. By investigating the differential usage of reactions across flux patterns for different environments, we observe a striking bimodal distribution in the activity of reactions. Motivated by this, we propose a simple algorithm to decompose the metabolic network into three subnetworks. It turns out that our reaction classifier, which is blind to the biochemical role of pathways, leads to three functionally relevant subnetworks that correspond to input, output, and intermediate parts of the metabolic network with distinct structural characteristics. Our decomposition method unveils a functional bow-tie organization of metabolic networks that is different from the bow-tie structure determined by graph-theoretic methods that do not incorporate functionality. PMID:23767567
A new method for studying caffeine's antioxygenic property: Peroxidase-Oxidase biochemical reaction
WANG Jun; CAI Ruxiu; LIN Zhixin; LIU Zhihong
2003-01-01
The effect of Caffeine on Peroxidase-Oxidase (PO) reaction was studied systematically in this paper. We proved that the valley of PO oscillation is the best phase angle which was used to research antioxygenic property by the Analyte Pulse Perturbation Technique (APP), based on investigating the mechanism. Area integral calculus was proposed to use in quantitative analysis for the first time. There is a good linear relationship (R = 0.9950) between the ratio of amplitude changes of PO oscillation and the concentration of caffeine in the range 4.61×10-7 mol/L-1.84×10-5 mol/L. A new method for analysis by PO oscillation was set up. We also investigated two-dimensional projections and Fourier spectrum of nonlinear complicate system--PO reaction which was perturbed by caffeine, in order to provide a theoretical basis for studying effects of kinds of antioxidants on life system.
On multistage homotopy-perturbation method applied to nonlinear biochemical reaction model
The total time evolution of the reactant concentrations in the basic enzyme-substrate reaction is simulated by an adaptation of the standard analytic homotopy-perturbation method (HPM). The standard HPM is converted into a hybrid numeric-analytic method called the multistage HPM (MHPM). The numerical results obtained from the MHPM and the classical fourth-order Runge-Kutta (RK4) method are in complete agreement
Recent developments in research on catalytic reaction networks
Roberto Serra
2013-09-01
Full Text Available Over the last years, analyses performed on a stochastic model of catalytic reaction networks have provided some indications about the reasons why wet-lab experiments hardly ever comply with the phase transition typically predicted by theoretical models with regard to the emergence of collectively self-replicating sets of molecule (also defined as autocatalytic sets, ACSs, a phenomenon that is often observed in nature and that is supposed to have played a major role in the emergence of the primitive forms of life. The model at issue has allowed to reveal that the emerging ACSs are characterized by a general dynamical fragility, which might explain the difficulty to observe them in lab experiments. In this work, the main results of the various analyses are reviewed, with particular regard to the factors able to affect the generic properties of catalytic reactions network, for what concerns, not only the probability of ACSs to be observed, but also the overall activity of the system, in terms of production of new species, reactions and matter.
N. O. Khromykh
2015-07-01
Full Text Available Morphometrical indexes, and spectrophotometrically measured protein and glutathione (GSH, GSSG contents and activity of peroxidase (POD, EC 1.11.1.7, glutathione-reductase (GR, EC 1.6.4.2 and glutathione S-transferase (GST, EС 2.5.1.18 were examined in Hordeum vulgare L. seedlings after 0.01 and 0.1 mg/l AgNPs treatment during 24 h. We tested the hypothesis that the action of nanoparticles has a stressful effect on the physiological and biochemical processes of seedlings. Growth of roots was inhibited and fresh weight decreased by 29% and 21% under low and high concentrations respectively. Conversely, leaf growth was intensified, and leaf length (16% and 18% and fresh weight (35% and 44% increased at low and high concentrations respectively. POD activity in roots increased by 26% and 7%, and decreased in leaves to 57% and 81% of control at low and high concentrations respectively. GSH content changed insignificantly, but GSSG content increased in roots (2 and 2.5-fold and in leaves (13% and 30% at both AgNPs concentrations. GSH/GSSG-ratio decreased in roots (1.9 and 2.6-fold and in leaves (1.1 and 1.3-fold at low and high concentrations respectively. GR activity decreased at a concentration of 0.01 mg/l (7% in roots and 17% in leaves respectively and increased at 0.1 mg/l (52% in roots and 6% in leaves. GST activity increased in leaves (52% and 78% at low and high concentrations but decreased by 17% in roots under high concentration of nanosilver. Thus, the action of AgNPs on barley seedlings had a dose-dependent and organ-specific character. The various directions of changes in growth, metabolic processes and activity of antioxidant defense systems appear to be a stress response of barley seedlings to the impact of AgNPs, which underlines the necessity of detailed study of plant intracellular processes exposed to the action of nanomaterial.
Hierarchical feedback modules and reaction hubs in cell signaling networks.
Jianfeng Xu
Full Text Available Despite much effort, identification of modular structures and study of their organizing and functional roles remain a formidable challenge in molecular systems biology, which, however, is essential in reaching a systematic understanding of large-scale cell regulation networks and hence gaining capacity of exerting effective interference to cell activity. Combining graph theoretic methods with available dynamics information, we successfully retrieved multiple feedback modules of three important signaling networks. These feedbacks are structurally arranged in a hierarchical way and dynamically produce layered temporal profiles of output signals. We found that global and local feedbacks act in very different ways and on distinct features of the information flow conveyed by signal transduction but work highly coordinately to implement specific biological functions. The redundancy embodied with multiple signal-relaying channels and feedback controls bestow great robustness and the reaction hubs seated at junctions of different paths announce their paramount importance through exquisite parameter management. The current investigation reveals intriguing general features of the organization of cell signaling networks and their relevance to biological function, which may find interesting applications in analysis, design and control of bio-networks.
Structural parameter identifiability analysis for dynamic reaction networks
Davidescu, Florin Paul; Jørgensen, Sten Bay
2008-01-01
A fundamental problem in model identification is to investigate whether unknown parameters in a given model structure potentially can be uniquely recovered from experimental data. This issue of global or structural identifiability is essential during nonlinear first principles model development...... where for a given set of measured variables it is desirable to investigate which parameters may be estimated prior to spending computational effort on the actual estimation. This contribution addresses the structural parameter identifiability problem for the typical case of reaction network models. The...
Adaptive hybrid simulations for multiscale stochastic reaction networks
The probability distribution describing the state of a Stochastic Reaction Network (SRN) evolves according to the Chemical Master Equation (CME). It is common to estimate its solution using Monte Carlo methods such as the Stochastic Simulation Algorithm (SSA). In many cases, these simulations can take an impractical amount of computational time. Therefore, many methods have been developed that approximate sample paths of the underlying stochastic process and estimate the solution of the CME. A prominent class of these methods include hybrid methods that partition the set of species and the set of reactions into discrete and continuous subsets. Such a partition separates the dynamics into a discrete and a continuous part. Simulating such a stochastic process can be computationally much easier than simulating the exact discrete stochastic process with SSA. Moreover, the quasi-stationary assumption to approximate the dynamics of fast subnetworks can be applied for certain classes of networks. However, as the dynamics of a SRN evolves, these partitions may have to be adapted during the simulation. We develop a hybrid method that approximates the solution of a CME by automatically partitioning the reactions and species sets into discrete and continuous components and applying the quasi-stationary assumption on identifiable fast subnetworks. Our method does not require any user intervention and it adapts to exploit the changing timescale separation between reactions and/or changing magnitudes of copy-numbers of constituent species. We demonstrate the efficiency of the proposed method by considering examples from systems biology and showing that very good approximations to the exact probability distributions can be achieved in significantly less computational time. This is especially the case for systems with oscillatory dynamics, where the system dynamics change considerably throughout the time-period of interest
Biochemical toxicology studies of methomyl and its kinetic reaction with cholinesterase in rats
The effect of the pesticide methomyl on acetylcholinesterase activities in brain and blood of male rats in vivo and the kinetics involved in their reaction in vitro were studied. Also, its effect on peroxidase action of catalase in vivo was studied using 14C- formate. The results showed that methomyl is a competitive inhibitor for acetylcholinesterase and the concentration levels that caused 50% inhibition of the enzyme activity were 2.1 x 10-2 M and 1.9 x 10-4 for brain and blood- Ache, respectively. The inhibition of acetylcholinesterase activity decreased to 67.5 % at the time of appearance of toxicity symptoms. The radioactivity eliminated in both the expired air and in urine was reduced
Cellulose was produced heterotrophically from different carbon substrates by carrot tissue cultures and Acetobacter xylinum (a cellulose-producing bacterium) and by castor bean seeds germinated in the dark, in each case in the presence of water having known concentration of oxygen-18 (18O). We used the relationship between the amount of 18O in the water and in the cellulose that was synthesized to determine the number and 18O content of the substrate oxygens that exchanged with water during the reactions leading to cellulose synthesis. Our observations support the hypothesis that oxygen isotope ratios of plant cellulose are determined by isotopic exchange occurring during hydration of carbonyl groups of the intermediates of cellulose synthesis. (author)
Crack, Jason C; Hutchings, Matthew I; Thomson, Andrew J; Le Brun, Nick E
2016-03-01
In Paracoccus denitrificans, three CRP/FNR family regulatory proteins, NarR, NnrR and FnrP, control the switch between aerobic and anaerobic (denitrification) respiration. FnrP is a [4Fe-4S] cluster-containing homologue of the archetypal O2 sensor FNR from E. coli and accordingly regulates genes encoding aerobic and anaerobic respiratory enzymes in response to O2, and also NO, availability. Here we show that FnrP undergoes O2-driven [4Fe-4S] to [2Fe-2S] cluster conversion that involves up to 2 O2 per cluster, with significant oxidation of released cluster sulfide to sulfane observed at higher O2 concentrations. The rate of the cluster reaction was found to be ~sixfold lower than that of E. coli FNR, suggesting that FnrP can remain transcriptionally active under microaerobic conditions. This is consistent with a role for FnrP in activating expression of the high O2 affinity cytochrome c oxidase under microaerobic conditions. Cluster conversion resulted in dissociation of the transcriptionally active FnrP dimer into monomers. Therefore, along with E. coli FNR, FnrP belongs to the subset of FNR proteins in which cluster type is correlated with association state. Interestingly, two key charged residues, Arg140 and Asp154, that have been shown to play key roles in the monomer-dimer equilibrium in E. coli FNR are not conserved in FnrP, indicating that different protomer interactions are important for this equilibrium. Finally, the FnrP [4Fe-4S] cluster is shown to undergo reaction with multiple NO molecules, resulting in iron nitrosyl species and dissociation into monomers. PMID:26790880
COEL: A Cloud-based Reaction Network Simulator
Peter eBanda
2016-04-01
Full Text Available Chemical Reaction Networks (CRNs are a formalism to describe the macroscopic behavior of chemical systems. We introduce COEL, a web- and cloud-based CRN simulation framework that does not require a local installation, runs simulations on a large computational grid, provides reliable database storage, and offers a visually pleasing and intuitive user interface. We present an overview of the underlying software, the technologies, and the main architectural approaches employed. Some of COEL's key features include ODE-based simulations of CRNs and multicompartment reaction networks with rich interaction options, a built-in plotting engine, automatic DNA-strand displacement transformation and visualization, SBML/Octave/Matlab export, and a built-in genetic-algorithm-based optimization toolbox for rate constants.COEL is an open-source project hosted on GitHub (http://dx.doi.org/10.5281/zenodo.46544, which allows interested research groups to deploy it on their own sever. Regular users can simply use the web instance at no cost at http://coel-sim.org. The framework is ideally suited for a collaborative use in both research and education.
无
2000-01-01
A new analytical method for the determination of urea-urease system based on biochemical reaction heat induced laser beam deflection is presented in this paper. With the method, the Michaelis constant (Km) of urease and apparent inhibition constant (Ki) of some metal ion inhibitors were measured respectively. This method was also used for the quantitative determination of metal ions with satisfactory result.
Multiplexing oscillatory biochemical signals.
de Ronde, Wiet; ten Wolde, Pieter Rein
2014-04-01
In recent years it has been increasingly recognized that biochemical signals are not necessarily constant in time and that the temporal dynamics of a signal can be the information carrier. Moreover, it is now well established that the protein signaling network of living cells has a bow-tie structure and that components are often shared between different signaling pathways. Here we show by mathematical modeling that living cells can multiplex a constant and an oscillatory signal: they can transmit these two signals simultaneously through a common signaling pathway, and yet respond to them specifically and reliably. We find that information transmission is reduced not only by noise arising from the intrinsic stochasticity of biochemical reactions, but also by crosstalk between the different channels. Yet, under biologically relevant conditions more than 2 bits of information can be transmitted per channel, even when the two signals are transmitted simultaneously. These observations suggest that oscillatory signals are ideal for multiplexing signals. PMID:24685537
Experimental (Network) and Evaluated Nuclear Reaction Data at NDS
Dr Simakov of Nuclear Data Services Unit in the Nuclear Data Section (NDS) gave a brief overview of the data compilation and evaluation activities in the nuclear data community: experimental nuclear reaction data (EXFOR, http://www-nds.iaea.org/exfor/) and evaluated nuclear reaction data (ENDF, http://www-nds.iaea.org/endf). The International Network of Nuclear Reaction Data Centres (NRDC) coordinated by NDS includes 14 Centres in 8 Countries (China, Hungary, India, Japan, Korea, Russian, Ukraine, USA) and 2 International Organizations (NEA, IAEA). It had the first meeting of four core centres (Brookhaven, Saclay, Obninsk, Vienna) in 1966 and the EXFOR was adopted as an official data exchange format. In 2000, IAEA implemented the EXFOR database as a relational multiform database and the EXFOR is a trusted, increasing and living database with 19100 experimental works (as of September 2011) and 141600 data tables. The EXFOR provides a compilation control system for selection of articles and compilation of data and the NRDC home page provides manuals, documents and codes. The nuclear data can be retrieved by the web-retrieval system or distributed on a DVD on request. The EXFOR data play a critical role in the development of evaluated nuclear reaction data. There are several major general purpose libraries: ENDF (US), CENDL (China), JEFF (EU), JENDL (Japan) and RUSFOND (Russia). In addition, there are special libraries for particular applications: EAF (European Activation File), FENDL (Fusion Evaluated Nuclear Data Library for ITER neutronics), IBANDL (Ion Beam Analysis Nuclear Data Library for surface analysis of solids), IRDF, DXS (Dosimetry, radiation damage and gas production data) and Medical portal. Dr V. Zerkin of NDS demonstrated the data retrieval from the EXFOR database and the ENDF library.
Moraes, Alvaro
2015-01-01
Epidemics have shaped, sometimes more than wars and natural disasters, demo- graphic aspects of human populations around the world, their health habits and their economies. Ebola and the Middle East Respiratory Syndrome (MERS) are clear and current examples of potential hazards at planetary scale. During the spread of an epidemic disease, there are phenomena, like the sudden extinction of the epidemic, that can not be captured by deterministic models. As a consequence, stochastic models have been proposed during the last decades. A typical forward problem in the stochastic setting could be the approximation of the expected number of infected individuals found in one month from now. On the other hand, a typical inverse problem could be, given a discretely observed set of epidemiological data, infer the transmission rate of the epidemic or its basic reproduction number. Markovian epidemic models are stochastic models belonging to a wide class of pure jump processes known as Stochastic Reaction Networks (SRNs), that are intended to describe the time evolution of interacting particle systems where one particle interacts with the others through a finite set of reaction channels. SRNs have been mainly developed to model biochemical reactions but they also have applications in neural networks, virus kinetics, and dynamics of social networks, among others. 4 This PhD thesis is focused on novel fast simulation algorithms and statistical inference methods for SRNs. Our novel Multi-level Monte Carlo (MLMC) hybrid simulation algorithms provide accurate estimates of expected values of a given observable of SRNs at a prescribed final time. They are designed to control the global approximation error up to a user-selected accuracy and up to a certain confidence level, and with near optimal computational work. We also present novel dual-weighted residual expansions for fast estimation of weak and strong errors arising from the MLMC methodology. Regarding the statistical inference
Variable elimination in chemical reaction networks with mass action kinetics
Feliu, Elisenda
2011-01-01
We consider chemical reaction networks taken with mass action kinetics. The steady states of such a system are solutions to a system of polynomial equations. Even for small systems the task of finding the solutions is daunting. We develop an algebraic framework and procedure for linear elimination of variables. The procedure reduces the variables in the system to a set of "core" variables by eliminating variables corresponding to a set of non-interacting species. The steady states are parameterized algebraically by the core variables, and a graphical condition is given for when a steady state with positive core variables necessarily have all variables positive. Further, we characterize graphically the sets of eliminated variables that are constrained by a conservation law and show that this conservation law takes a specific form.
Phelps, Michael E.
2009-09-01
Radiotracer techniques are used in environmental sciences, geology, biology and medicine. Radiotracers with Positron Emission Tomography (PET) provided biological examinations of ~3 million patients 2008. Despite the success of positron labeled tracers in many sciences, there is limited access in an affordable and convenient manner to develop and use new tracers. Integrated microfluidic chips are a new technology well matched to the concentrations of tracers. Our goal is to develop microfluidic chips and new synthesis approaches to enable wide dissemination of diverse types of tracers at low cost, and to produce new generations of radiochemists for which there are many unfilled jobs. The program objectives are to: 1. Develop an integrated microfluidic platform technology for synthesizing and 18F-labeling diverse arrays of different classes of molecules. 2. Incorporate microfluidic chips into small PC controlled devices (“Synthesizer”) with a platform interfaced to PC for electronic and fluid input/out control. 3. Establish a de-centralized model with Synthesizers for discovering and producing molecular imaging probes, only requiring delivery of inexpensive [18F]fluoride ion from commercial PET radiopharmacies vs the centralized approach of cyclotron facilities synthesizing and shipping a few different types of 18F-probes. 4. Develop a position sensitive avalanche photo diode (PSAPD) camera for beta particles embedded in a microfluidic chip for imaging and measuring transport and biochemical reaction rates to valid new 18F-labeled probes in an array of cell cultures. These objectives are met within a research and educational program integrating radio-chemistry, synthetic chemistry, biochemistry, engineering and biology in the Crump Institute for Molecular Imaging. The Radiochemistry Training Program exposes PhD and post doctoral students to molecular imaging in vitro in cells and microorganisms in microfluidic chips and in vivo with PET, from new technologies
Biochemical switching device: biomimetic approach and application to neural network study.
Okamoto, M
1992-06-01
There are many examples of enzymes that share substrates or cofactors in a cyclic manner. Techniques have been developed that use cyclic enzyme systems to assay quantitatively small amounts of biochemical substances (cofactor, substrate), however, only a few studies of the control of these systems have been published. The author previously showed with computer simulations that cyclic enzyme systems have the reliability of ON-OFF types of operation (McCulloch-Pitts' neuronic equation) capable of storing short-memory, and the applicability for a switching circuit in a biocomputer. This paper introduces a unique switching mechanism of cyclic enzyme system (basic switching element), and next, building the integrated biochemical switching system being composed of the basic switching element, shows the physiological phenomenon termed 'selective elimination of synapses' generally produced as a result of low-frequency train of electrical stimuli to the synapses (Kuroda, Y. 1989) Neurochem. Int. 14, 309-319). PMID:1368350
Rule-based models, which are typically formulated to represent cell signaling systems, can now be simulated via various network-free simulation methods. In a network-free method, reaction rates are calculated for rules that characterize molecular interactions, and these rule rates, which each correspond to the cumulative rate of all reactions implied by a rule, are used to perform a stochastic simulation of reaction kinetics. Network-free methods, which can be viewed as generalizations of Gillespie's method, are so named because these methods do not require that a list of individual reactions implied by a set of rules be explicitly generated, which is a requirement of other methods for simulating rule-based models. This requirement is impractical for rule sets that imply large reaction networks (i.e. long lists of individual reactions), as reaction network generation is expensive. Here, we compare the network-free simulation methods implemented in RuleMonkey and NFsim, general-purpose software tools for simulating rule-based models encoded in the BioNetGen language. The method implemented in NFsim uses rejection sampling to correct overestimates of rule rates, which introduces null events (i.e. time steps that do not change the state of the system being simulated). The method implemented in RuleMonkey uses iterative updates to track rule rates exactly, which avoids null events. To ensure a fair comparison of the two methods, we developed implementations of the rejection and rejection-free methods specific to a particular class of kinetic models for multivalent ligand–receptor interactions. These implementations were written with the intention of making them as much alike as possible, minimizing the contribution of irrelevant coding differences to efficiency differences. Simulation results show that performance of the rejection method is equal to or better than that of the rejection-free method over wide parameter ranges. However, when parameter values are
Defect reaction network in Si-doped InAs. Numerical predictions.
Schultz, Peter A. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
2015-05-01
This Report characterizes the defects in the def ect reaction network in silicon - doped, n - type InAs predicted with first principles density functional theory. The reaction network is deduced by following exothermic defect reactions starting with the initially mobile interstitial defects reacting with common displacement damage defects in Si - doped InAs , until culminating in immobile reaction p roducts. The defect reactions and reaction energies are tabulated, along with the properties of all the silicon - related defects in the reaction network. This Report serves to extend the results for the properties of intrinsic defects in bulk InAs as colla ted in SAND 2013 - 2477 : Simple intrinsic defects in InAs : Numerical predictions to include Si - containing simple defects likely to be present in a radiation - induced defect reaction sequence . This page intentionally left blank
Comprehensive estimation of input signals and dynamics in biochemical reaction networks
Schelker, M; Raue, A.; Timmer, Jens; Kreutz, C.
2012-01-01
Motivation: Cellular information processing can be described mathematically using differential equations. Often, external stimulation of cells by compounds such as drugs or hormones leading to activation has to be considered. Mathematically, the stimulus is represented by a time-dependent input function. Parameters such as rate constants of the molecular interactions are often unknown and need to be estimated from experimental data, e.g. by maximum likelihood estimation. For this purpose, the...
der Schaft, Arjan van; Rao, Shodhan; Jayawardhana, Bayu
2011-01-01
Motivated by recent progress on the interplay between graph theory, dynamics, and systems theory, we revisit the analysis of chemical reaction networks described by mass action kinetics. For reaction networks possessing a thermodynamic equilibrium we derive a compact formulation exhibiting at the same time the structure of the complex graph and the stoichiometry of the network, and which admits a direct thermodynamical interpretation. This formulation allows us to easily characterize the set ...
Properties of Random Complex Chemical Reaction Networks and Their Relevance to Biological Toy Models
Bigan, Erwan; Steyaert, Jean-Marc; Douady, Stéphane
2013-01-01
We investigate the properties of large random conservative chemical reaction networks composed of elementary reactions endowed with either mass-action or saturating kinetics, assigning kinetic parameters in a thermodynamically-consistent manner. We find that such complex networks exhibit qualitatively similar behavior when fed with external nutrient flux. The nutrient is preferentially transformed into one specific chemical that is an intrinsic property of the network. We propose a self-consi...
Sensitivity of chemical reaction networks: a structural approach 3. Regular multimolecular systems
Brehm, Bernhard
2016-01-01
We present a systematic mathematical analysis of the qualitative steady-state response to rate perturbations in large classes of reaction networks. This includes multimolecular reactions and allows for catalysis, enzymatic reactions, multiple reaction products, nonmonotone rate functions, and non-closed autonomous systems. Our structural sensitivity analysis is based on the stoichiometry of the reaction network, only. It does not require numerical data on reaction rates. Instead, we impose mild and generic nondegeneracy conditions of algebraic type. From the structural data, only, we derive which steady-state concentrations are sensitive to, and hence influenced by, changes of any particular reaction rate - and which are not. We also establish transitivity properties for influences involving rate perturbations. This allows us to derive an influence graph which globally summarizes the influence pattern of the given network. The influence graph allows the computational, but meaningful, automatic identification ...
Metabolome based reaction graphs of M. tuberculosis and M. leprae: a comparative network analysis.
Ketki D Verkhedkar
Full Text Available BACKGROUND: Several types of networks, such as transcriptional, metabolic or protein-protein interaction networks of various organisms have been constructed, that have provided a variety of insights into metabolism and regulation. Here, we seek to exploit the reaction-based networks of three organisms for comparative genomics. We use concepts from spectral graph theory to systematically determine how differences in basic metabolism of organisms are reflected at the systems level and in the overall topological structures of their metabolic networks. METHODOLOGY/PRINCIPAL FINDINGS: Metabolome-based reaction networks of Mycobacterium tuberculosis, Mycobacterium leprae and Escherichia coli have been constructed based on the KEGG LIGAND database, followed by graph spectral analysis of the network to identify hubs as well as the sub-clustering of reactions. The shortest and alternate paths in the reaction networks have also been examined. Sub-cluster profiling demonstrates that reactions of the mycolic acid pathway in mycobacteria form a tightly connected sub-cluster. Identification of hubs reveals reactions involving glutamate to be central to mycobacterial metabolism, and pyruvate to be at the centre of the E. coli metabolome. The analysis of shortest paths between reactions has revealed several paths that are shorter than well established pathways. CONCLUSIONS: We conclude that severe downsizing of the leprae genome has not significantly altered the global structure of its reaction network but has reduced the total number of alternate paths between its reactions while keeping the shortest paths between them intact. The hubs in the mycobacterial networks that are absent in the human metabolome can be explored as potential drug targets. This work demonstrates the usefulness of constructing metabolome based networks of organisms and the feasibility of their analyses through graph spectral methods. The insights obtained from such studies provide a
Analysis of weblike network structures of directed graphs for chemical reactions in methane plasmas
Chemical reactions of molecular gases like methane are so complicated that a chart of decomposed and/or synthesized species originating from molecules in plasma resembles a weblike network in which we write down species and reactions among them. Here we consider properties of the network structures of chemical reactions in methane plasmas. In the network, atoms/molecules/radical species are assumed to form nodes and chemical reactions correspond to directed edges in the terminology of graph theory. Investigation of the centrality index reveals importance of CH3 in the global chemical reaction, and difference of an index for each radical species between cases with and without electrons clarifies that the electrons are at an influential position to tighten the network structure
Efficient Characterization of Parametric Uncertainty of Complex (Bio)chemical Networks.
Claudia Schillings; Mikael Sunnåker; Jörg Stelling; Christoph Schwab
2015-01-01
Parametric uncertainty is a particularly challenging and relevant aspect of systems analysis in domains such as systems biology where, both for inference and for assessing prediction uncertainties, it is essential to characterize the system behavior globally in the parameter space. However, current methods based on local approximations or on Monte-Carlo sampling cope only insufficiently with high-dimensional parameter spaces associated with complex network models. Here, we propose an alternat...
Efficient Characterization of Parametric Uncertainty of Complex (Bio)chemical Networks
Claudia Schillings; Mikael Sunnåker; Jörg Stelling; Christoph Schwab
2015-01-01
Parametric uncertainty is a particularly challenging and relevant aspect of systems analysis in domains such as systems biology where, both for inference and for assessing prediction uncertainties, it is essential to characterize the system behavior globally in the parameter space. However, current methods based on local approximations or on Monte-Carlo sampling cope only insufficiently with high-dimensional parameter spaces associated with complex network models. Here, we propose an alternat...
Unravelling the Maillard reaction network by multiresponse kinetic modelling
Martins, S.I.F.S.
2003-01-01
The Maillard reaction is an important reaction in food industry. It is responsible for the formation of colour and aroma, as well as toxic compounds as the recent discovered acrylamide. The knowledge of kinetic parameters, such as rate constants and activation energy, is necessary to predict its extent and, consequently, to optimise it. Each of the chapters presented in this thesis can be seen as a necessary step to succeed in applying multiresponse kinetic modelling in a complex reaction, su...
Teodorowicz, M.; Fiedorowicz, E.; Kostyra, H.; Wichers, H J; Kostyra, E.
2013-01-01
Purpose The purpose of this study was to determine the influence of Maillard reaction (MR, glycation) on biochemical and biological properties of the major peanut allergen Ara h 1. Methods Three different time/temperature conditions of treatment were applied (37, 60, and 145 °C). The extent of MR was assessed by SDS-PAGE, loss of free amino groups, fluorescence intensity, content of bound sugar and fructosamine. The Caco-2 model system was applied to study effects of hydrolysed and non-hydrol...
Efficient Characterization of Parametric Uncertainty of Complex (Bio)chemical Networks.
Schillings, Claudia; Sunnåker, Mikael; Stelling, Jörg; Schwab, Christoph
2015-08-01
Parametric uncertainty is a particularly challenging and relevant aspect of systems analysis in domains such as systems biology where, both for inference and for assessing prediction uncertainties, it is essential to characterize the system behavior globally in the parameter space. However, current methods based on local approximations or on Monte-Carlo sampling cope only insufficiently with high-dimensional parameter spaces associated with complex network models. Here, we propose an alternative deterministic methodology that relies on sparse polynomial approximations. We propose a deterministic computational interpolation scheme which identifies most significant expansion coefficients adaptively. We present its performance in kinetic model equations from computational systems biology with several hundred parameters and state variables, leading to numerical approximations of the parametric solution on the entire parameter space. The scheme is based on adaptive Smolyak interpolation of the parametric solution at judiciously and adaptively chosen points in parameter space. As Monte-Carlo sampling, it is "non-intrusive" and well-suited for massively parallel implementation, but affords higher convergence rates. This opens up new avenues for large-scale dynamic network analysis by enabling scaling for many applications, including parameter estimation, uncertainty quantification, and systems design. PMID:26317784
Lou Xuyang [Research Center of Control Science and Engineering, Southern Yangtze University, 1800 Lihu Road, Wuxi, Jiangsu 214122 (China); Cui Baotong [Research Center of Control Science and Engineering, Southern Yangtze University, 1800 Lihu Road, Wuxi, Jiangsu 214122 (China)]. E-mail: btcui@sohu.com
2007-07-15
Employing Lyapunov functional method, we analyze the ultimate boundedness and global exponential stability of a class of reaction-diffusion cellular neural networks with time-varying delays. Some new criteria are obtained to ensure ultimate boundedness and global exponential stability of delayed reaction-diffusion cellular neural networks (DRCNNs). Without assuming that the activation functions f {sub ijl}(.) are bounded, the results extend and improve the earlier publications.
Uysal, Ebru; Meral, Yüce; Meral, Yuce; Hasan KURT
2014-01-01
Recently, the Polymerase Chain Reaction technique has begun to benefit from nanotechnology. In this paper, effects of carbon nanotubes in the Polymerase Chain Reaction were investigated by Electrophoresis, Circular Dichroism Spectrometry and Dynamic Light Scattering Techniques. The unique ability to amplify low copy number DNA within minutes has made in vitro Polymerase Chain Reaction (PCR) one of the most essential techniques in modern biology. In order to harness this technique to its full ...
Mehta, Pankaj; Lang, Alex H.; Schwab, David J.
2016-03-01
A central goal of synthetic biology is to design sophisticated synthetic cellular circuits that can perform complex computations and information processing tasks in response to specific inputs. The tremendous advances in our ability to understand and manipulate cellular information processing networks raises several fundamental physics questions: How do the molecular components of cellular circuits exploit energy consumption to improve information processing? Can one utilize ideas from thermodynamics to improve the design of synthetic cellular circuits and modules? Here, we summarize recent theoretical work addressing these questions. Energy consumption in cellular circuits serves five basic purposes: (1) increasing specificity, (2) manipulating dynamics, (3) reducing variability, (4) amplifying signal, and (5) erasing memory. We demonstrate these ideas using several simple examples and discuss the implications of these theoretical ideas for the emerging field of synthetic biology. We conclude by discussing how it may be possible to overcome these limitations using "post-translational" synthetic biology that exploits reversible protein modification.
On the graph and systems analysis of reversible chemical reaction networks with mass action kinetics
Rao, Shodhan; Jayawardhana, Bayu; Schaft, Arjan van der
2012-01-01
Motivated by the recent progresses on the interplay between the graph theory and systems theory, we revisit the analysis of reversible chemical reaction networks described by mass action kinetics by reformulating it using the graph knowledge of the underlying networks. Based on this formulation, we
The Impact of Different Sources of Fluctuations on Mutual Information in Biochemical Networks.
Chevalier, Michael; Venturelli, Ophelia; El-Samad, Hana
2015-10-01
Stochastic fluctuations in signaling and gene expression limit the ability of cells to sense the state of their environment, transfer this information along cellular pathways, and respond to it with high precision. Mutual information is now often used to quantify the fidelity with which information is transmitted along a cellular pathway. Mutual information calculations from experimental data have mostly generated low values, suggesting that cells might have relatively low signal transmission fidelity. In this work, we demonstrate that mutual information calculations might be artificially lowered by cell-to-cell variability in both initial conditions and slowly fluctuating global factors across the population. We carry out our analysis computationally using a simple signaling pathway and demonstrate that in the presence of slow global fluctuations, every cell might have its own high information transmission capacity but that population averaging underestimates this value. We also construct a simple synthetic transcriptional network and demonstrate using experimental measurements coupled to computational modeling that its operation is dominated by slow global variability, and hence that its mutual information is underestimated by a population averaged calculation. PMID:26484538
Drift-Implicit Multi-Level Monte Carlo Tau-Leap Methods for Stochastic Reaction Networks
Ben Hammouda, Chiheb
2015-05-12
In biochemical systems, stochastic e↵ects can be caused by the presence of small numbers of certain reactant molecules. In this setting, discrete state-space and stochastic simulation approaches were proved to be more relevant than continuous state-space and deterministic ones. These stochastic models constitute the theory of stochastic reaction networks (SRNs). Furthermore, in some cases, the dynamics of fast and slow time scales can be well separated and this is characterized by what is called sti↵ness. For such problems, the existing discrete space-state stochastic path simulation methods, such as the stochastic simulation algorithm (SSA) and the explicit tau-leap method, can be very slow. Therefore, implicit tau-leap approxima- tions were developed to improve the numerical stability and provide more e cient simulation algorithms for these systems. One of the interesting tasks for SRNs is to approximate the expected values of some observables of the process at a certain fixed time T. This is can be achieved using Monte Carlo (MC) techniques. However, in a recent work, Anderson and Higham in 2013, proposed a more computationally e cient method which combines multi-level Monte Carlo (MLMC) technique with explicit tau-leap schemes. In this MSc thesis, we propose new fast stochastic algorithm, particularly designed 5 to address sti↵ systems, for approximating the expected values of some observables of SRNs. In fact, we take advantage of the idea of MLMC techniques and drift-implicit tau-leap approximation to construct a drift-implicit MLMC tau-leap estimator. In addition to accurately estimating the expected values of a given observable of SRNs at a final time T , our proposed estimator ensures the numerical stability with a lower cost than the MLMC explicit tau-leap algorithm, for systems including simultane- ously fast and slow species. The key contribution of our work is the coupling of two drift-implicit tau-leap paths, which is the basic brick for
Frank, T. D.
2013-08-01
We derive a nonlinear limit cycle model for oscillatory mood variations as observed in patients with cycling bipolar disorder. To this end, we consider two signaling pathways leading to the activation of two enzymes that play a key role for cellular and neural processes. We model pathway cross-talk in terms of an inhibitory impact of the first pathway on the second and an excitatory impact of the second on the first. The model also involves a negative feedback loop (inhibitory self-regulation) for the first pathway and a positive feedback loop (excitatory self-regulation) for the second pathway. We demonstrate that due to the cross-talk the biochemical dynamics is described by an oscillator equation. Under disease-free conditions the oscillatory system exhibits a stable fixed point. The breakdown of the self-inhibition of the first pathway at higher concentration levels is studied by means of a scalar control parameter ξ, where ξ equal to zero refers to intact self-inhibition at all concentration levels. Under certain conditions, stable limit cycle solutions emerge at critical parameter values of ξ larger than zero. These oscillations mimic pathological cycling mood variations that emerge due to a disease-induced bifurcation. Consequently, our modeling analysis supports the notion of bipolar disorder as a dynamical disease. In addition, our study establishes a connection between mechanistic biochemical modeling of bipolar disorder and phenomenological nonlinear oscillator approaches to bipolar disorder suggested in the literature.
Preclusion of switch behavior in reaction networks with mass-action kinetics
Feliu, Elisenda; Wiuf, C.
2012-01-01
We study networks taken with mass-action kinetics and provide a Jacobian criterion that applies to an arbitrary network to preclude the existence of multiple positive steady states within any stoichiometric class for any choice of rate constants. We are concerned with the characterization of...... injective networks, that is, networks for which the species formation rate function is injective in the interior of the positive orthant within each stoichiometric class. We show that a network is injective if and only if the determinant of the Jacobian of a certain function does not vanish. The function...... precludes the existence of degenerate steady states. Further, we relate injectivity of a network to that of the network obtained by adding outflow, or degradation, reactions for all species....
Automatic Verification of Biochemical Network Using Model Checking Method%基于模型校核的生化网络自动辨别方法
Jinkyung Kim; Younghee Lee; Il Moon
2008-01-01
This study focuses on automatic searching and verifying methods for the reachability, transition logics and hierarchical structure in all possible paths of biological processes using model checking. The automatic search and verification for alternative paths within complex and large networks in biological process can provide a consid-erable amount of solutions, which is difficult to handle manually. Model checking is an automatic method for veri-fying if a circuit or a condition, expressed as a concurrent transition system, satisfies a set of properties expressed ina temporal logic, such as computational tree logic (CTL). This article represents that model checking is feasible in biochemical network verification and it shows certain advantages over simulation for querying and searching of special behavioral properties in biochemical processes.
Reaction-Diffusion Processes on Random and Scale-Free Networks
Banerjee, Subhasis; Mallick, Shrestha Basu; Bose, Indrani
We study the discrete Gierer-Meinhardt model of reaction-diffusion on three different types of networks: regular, random and scale-free. The model dynamics lead to the formation of stationary Turing patterns in the steady state in certain parameter regions. Some general features of the patterns are studied through numerical simulation. The results for the random and scale-free networks show a marked difference from those in the case of the regular network. The difference may be ascribed to the small world character of the first two types of networks.
The US nuclear reaction data network. Summary of the first meeting, March 13 & 14 1996
NONE
1996-03-01
The first meeting of the US Nuclear Reaction Data Network (USNRDN) was held at the Colorado School of Mines, March 13-14, 1996 chaired by F. Edward Cecil. The Agenda of the meeting is attached. The Network, its mission, products and services; related nuclear data and data networks, members, and organization are described in Attachment 1. The following progress reports from the members of the USNRDN were distributed prior to the meeting and are given as Attachment 2. (1) Measurements and Development of Analytic Techniques for Basic Nuclear Physics and Nuclear Applications; (2) Nuclear Reaction Data Activities at the National Nuclear Data Center; (3) Studies of nuclear reactions at very low energies; (4) Nuclear Reaction Data Activities, Nuclear Data Group; (5) Progress in Neutron Physics at Los Alamos - Experiments; (6) Nuclear Reaction Data Activities in Group T2; (7) Progress Report for the US Nuclear Reaction Data Network Meeting; (8) Nuclear Astrophysics Research Group (ORNL); (9) Progress Report from Ohio University; (10) Exciton Model Phenomenology; and (11) Progress Report for Coordination Meeting USNRDN.
The US nuclear reaction data network. Summary of the first meeting, March 13 ampersand 14 1996
The first meeting of the US Nuclear Reaction Data Network (USNRDN) was held at the Colorado School of Mines, March 13-14, 1996 chaired by F. Edward Cecil. The Agenda of the meeting is attached. The Network, its mission, products and services; related nuclear data and data networks, members, and organization are described in Attachment 1. The following progress reports from the members of the USNRDN were distributed prior to the meeting and are given as Attachment 2. (1) Measurements and Development of Analytic Techniques for Basic Nuclear Physics and Nuclear Applications; (2) Nuclear Reaction Data Activities at the National Nuclear Data Center; (3) Studies of nuclear reactions at very low energies; (4) Nuclear Reaction Data Activities, Nuclear Data Group; (5) Progress in Neutron Physics at Los Alamos - Experiments; (6) Nuclear Reaction Data Activities in Group T2; (7) Progress Report for the US Nuclear Reaction Data Network Meeting; (8) Nuclear Astrophysics Research Group (ORNL); (9) Progress Report from Ohio University; (10) Exciton Model Phenomenology; and (11) Progress Report for Coordination Meeting USNRDN
Reconfigurable neuromorphic computation in biochemical systems.
Chiang, Hui-Ju Katherine; Jiang, Jie-Hong R; Fages, Francois
2015-08-01
Implementing application-specific computation and control tasks within a biochemical system has been an important pursuit in synthetic biology. Most synthetic designs to date have focused on realizing systems of fixed functions using specifically engineered components, thus lacking flexibility to adapt to uncertain and dynamically-changing environments. To remedy this limitation, an analog and modularized approach to realize reconfigurable neuromorphic computation with biochemical reactions is presented. We propose a biochemical neural network consisting of neuronal modules and interconnects that are both reconfigurable through external or internal control over the concentrations of certain molecular species. Case studies on classification and machine learning applications using the DNA strain displacement technology demonstrate the effectiveness of our design in both reconfiguration and autonomous adaptation. PMID:26736417
A cascade reaction network mimicking the basic functional steps of adaptive immune response
Han, Da; Wu, Cuichen; You, Mingxu; Zhang, Tao; Wan, Shuo; Chen, Tao; Qiu, Liping; Zheng, Zheng; Liang, Hao; Tan, Weihong
2015-10-01
Biological systems use complex ‘information-processing cores’ composed of molecular networks to coordinate their external environment and internal states. An example of this is the acquired, or adaptive, immune system (AIS), which is composed of both humoral and cell-mediated components. Here we report the step-by-step construction of a prototype mimic of the AIS that we call an adaptive immune response simulator (AIRS). DNA and enzymes are used as simple artificial analogues of the components of the AIS to create a system that responds to specific molecular stimuli in vitro. We show that this network of reactions can function in a manner that is superficially similar to the most basic responses of the vertebrate AIS, including reaction sequences that mimic both humoral and cellular responses. As such, AIRS provides guidelines for the design and engineering of artificial reaction networks and molecular devices.
Nakagawa, Masaki; Togashi, Yuichi
2016-01-01
Cell activities primarily depend on chemical reactions, especially those mediated by enzymes, and this has led to these activities being modeled as catalytic reaction networks. Although deterministic ordinary differential equations of concentrations (rate equations) have been widely used for modeling purposes in the field of systems biology, it has been pointed out that these catalytic reaction networks may behave in a way that is qualitatively different from such deterministic representation when the number of molecules for certain chemical species in the system is small. Apart from this, representing these phenomena by simple binary (on/off) systems that omit the quantities would also not be feasible. As recent experiments have revealed the existence of rare chemical species in cells, the importance of being able to model potential small-number phenomena is being recognized. However, most preceding studies were based on numerical simulations, and theoretical frameworks to analyze these phenomena have not been sufficiently developed. Motivated by the small-number issue, this work aimed to develop an analytical framework for the chemical master equation describing the distributional behavior of catalytic reaction networks. For simplicity, we considered networks consisting of two-body catalytic reactions. We used the probability generating function method to obtain the steady-state solutions of the chemical master equation without specifying the parameters. We obtained the time evolution equations of the first- and second-order moments of concentrations, and the steady-state analytical solution of the chemical master equation under certain conditions. These results led to the rank conservation law, the connecting state to the winner-takes-all state, and analysis of 2-molecules M-species systems. A possible interpretation of the theoretical conclusion for actual biochemical pathways is also discussed. PMID:27047384
Lyapunov Functions, Stationary Distributions, and Non-equilibrium Potential for Reaction Networks
Anderson, David F.; Craciun, Gheorghe; Gopalkrishnan, Manoj; Wiuf, Carsten Henrik
2015-01-01
We consider the relationship between stationary distributions for stochastic models of reaction systems and Lyapunov functions for their deterministic counterparts. Specifically, we derive the well-known Lyapunov function of reaction network theory as a scaling limit of the non-equilibrium potent......We consider the relationship between stationary distributions for stochastic models of reaction systems and Lyapunov functions for their deterministic counterparts. Specifically, we derive the well-known Lyapunov function of reaction network theory as a scaling limit of the non......-equilibrium potential of the stationary distribution of stochastically modeled complex balanced systems. We extend this result to general birth–death models and demonstrate via example that similar scaling limits can yield Lyapunov functions even for models that are not complex or detailed balanced, and may even have...
Report on the IAEA Technical Meeting of the International Network of Nuclear Reaction Data Centres
This report summarizes the IAEA Technical Meeting of the International Network of Nuclear Reaction Data Centres (biennial Data Centre Heads Meeting), held at the Japan Nuclear Reaction Data Centre, Hokkaido University, Sapporo, Japan, from 20 - 23 April 2010. The meeting was attended by 27 participants from 12 cooperating data centres of seven Member States and two International Organizations. The report contains a summary of the meeting, the conclusions and actions, the lists of working papers and presentations presented at the meeting. This report summarizes the IAEA Technical Meeting of the International Network of Nuclear Reaction Data Centres (biennial Data Centre Heads Meeting), held at the Japan Nuclear Reaction Data Centre, Hokkaido University, Sapporo, Japan, from 20 - 23 April 2010. The meeting was attended by 27 participants from 12 cooperating data centres of seven Member States and two International Organizations. The report contains a summary of the meeting, the conclusions and actions, the lists of working papers and presentations presented at the meeting. (author)
The origin of large molecules in primordial autocatalytic reaction networks
Giri, Varun
2011-01-01
Large molecules such as proteins and nucleic acids are crucial for life, yet their primordial origin remains a major puzzle. The production of large molecules, as we know it today, requires good catalysts, and the only good catalysts we know that can accomplish this task consist of large molecules. Thus the origin of large molecules is a chicken and egg problem in chemistry. Here we present a mechanism, based on autocatalytic sets (ACSs), that is a possible solution to this problem. We discuss a mathematical model describing the population dynamics of molecules in a stylized but prebiotically plausible chemistry. Large molecules can be produced in this chemistry by the coalescing of smaller ones, with the smallest molecules, the `food set', being buffered. Some of the reactions can be catalyzed by molecules within the chemistry with varying catalytic strengths. Normally the concentrations of large molecules in such a scenario are very small, diminishing exponentially with their size. ACSs, if present in the c...
Julio Michael Stern; Fabio Nakano
2014-01-01
This article presents a simple derivation of optimization models for reaction networks leading to a generalized form of the mass-action law, and compares the formal structure of Minimum Information Divergence, Quadratic Programming and Kirchhoff type network models. These optimization models are used in related articles to develop and illustrate the operation of ontology alignment algorithms and to discuss closely connected issues concerning the epistemological and statistical significance of...
On the graph and systems analysis of reversible chemical reaction networks with mass action kinetics
Rao, Shodhan; Jayawardhana, Bayu; der Schaft, Arjan van
2012-01-01
Motivated by the recent progresses on the interplay between the graph theory and systems theory, we revisit the analysis of reversible chemical reaction networks described by mass action kinetics by reformulating it using the graph knowledge of the underlying networks. Based on this formulation, we can characterize the space of equilibrium points and provide simple dynamical analysis on the state space modulo the space of equilibrium points.
The cost and capacity of signaling in the Escherichia coli protein reaction network
In systems biology new ways are required to analyze the large amount of existing data on regulation of cellular processes. Recent work can be roughly classified into either dynamical models of well-described subsystems, or coarse-grained descriptions of the topology of the molecular networks at the scale of the whole organism. In order to bridge these two disparate approaches one needs to develop simplified descriptions of dynamics and topological measures which address the propagation of signals in molecular networks. Transmission of a signal across a reaction node depends on the presence of other reactants. It will typically be more demanding to transmit a signal across a reaction node with more input links. Sending signals along a path with several subsequent reaction nodes also increases the constraints on the presence of other proteins in the overall network. Therefore counting in and out links along reactions of a potential pathway can give insight into the signaling properties of a particular molecular network. Here, we consider the directed network of protein regulation in E. coli, characterizing its modularity in terms of its potential to transmit signals. We demonstrate that the simplest measure based on identifying subnetworks of strong components, within which each node could send a signal to every other node, does indeed partition the network into functional modules. We suggest that the total number of reactants needed to send a signal between two nodes in the network can be considered as the cost associated with transmitting this signal. Similarly we define spread as the number of reaction products that could be influenced by transmission of a successful signal. Our considerations open for a new class of network measures that implicitly utilize the constrained repertoire of chemical modifications of any biological molecule. The counting of cost and spread connects the topology of networks to the specificity of signaling across the network. Thereby, we
Al-Radhawi, M. Ali; Angeli, David
2014-01-01
Piecewise-Linear in Rates (PWLR) Lyapunov functions are introduced for a class of Chemical Reaction Networks (CRNs). In addition to their simple structure, these functions are robust with respect to arbitrary monotone reaction rates, of which mass-action is a special case. The existence of such functions ensures the convergence of trajectories towards equilibria, and guarantee their asymptotic stability with respect to the corresponding stoichiometric compatibility class. We give the definiti...
Aikawa, M.; Arnould, M.; Goriely, S.; Jorissen, A.; Takahashi, K.
2005-01-01
Nuclear reaction rates are quantities of fundamental importance in astrophysics. Substantial efforts have been devoted in the last decades to measure or calculate them. The present paper presents for the first time a detailed description of the Brussels nuclear reaction rate library BRUSLIB and of the nuclear network generator NETGEN so as to make these nuclear data packages easily accessible to astrophysicists for a large variety of applications. BRUSLIB is made of two parts. The first one c...
Growth states of catalytic reaction networks exhibiting energy metabolism
Kondo, Yohei; Kaneko, Kunihiko
2011-07-01
All cells derive nutrition by absorbing some chemical and energy resources from the environment; these resources are used by the cells to reproduce the chemicals within them, which in turn leads to an increase in their volume. In this study we introduce a protocell model exhibiting catalytic reaction dynamics, energy metabolism, and cell growth. Results of extensive simulations of this model show the existence of four phases with regard to the rates of both the influx of resources and cell growth. These phases include an active phase with high influx and high growth rates, an inefficient phase with high influx but low growth rates, a quasistatic phase with low influx and low growth rates, and a death phase with negative growth rate. A mean field model well explains the transition among these phases as bifurcations. The statistical distribution of the active phase is characterized by a power law, and that of the inefficient phase is characterized by a nearly equilibrium distribution. We also discuss the relevance of the results of this study to distinct states in the existing cells.
Toward a self-consistent and unitary reaction network for big bang nucleosynthesis
Unitarity, the mathematical expression of the conservation of probability in multichannel reactions, is an essential ingredient in the development of accurate nuclear reaction networks appropriate for nucleosynthesis in a variety of environments. We describe our ongoing program to develop a 'unitary reaction network' for the big-bang nucleosynthesis environment and look at an example of the need and power of unitary parametrizations of nuclear scattering and reaction data. Recent attention has been focused on the possible role of the 9B compound nuclear system in the resonant destruction of 7Li during primordial nucleosynthesis. We have studied reactions in the 9B compound system with a multichannel, two-body unitary R-matrix code (EDA) using the known elastic and reaction data, in a four-channel treatment. The data include elastic 6Li(3He,3He)6Li differential cross sections from 0.7 to 2.0 MeV, integrated reaction cross sections for energies from 0.7 to 5.0 MeV for 6Li(3He,p)8Be* and from 0.4 to 5.0 MeV for the 6Li(3He,γ)7Be reaction. Capture data have been added to the previous analysis with integrated cross section measurements from 0.7 to 0.825 MeV for 6Li(3He,γ)9B. The resulting resonance parameters are compared with tabulated values from TUNL Nuclear Data Group analyses. Previously unidentified resonances are noted and the relevance of this analysis and a unitary reaction network for big-bang nucleosynthesis are emphasized. (author)
Adaptive exponential synchronization of delayed neural networks with reaction-diffusion terms
Sheng Li [School of Communication and Control Engineering, Jiangnan University, 1800 Lihu Road, Wuxi, Jiangsu 214122 (China); Yang Huizhong [School of Communication and Control Engineering, Jiangnan University, 1800 Lihu Road, Wuxi, Jiangsu 214122 (China)], E-mail: victory8209@yahoo.com.cn; Lou Xuyang [School of Communication and Control Engineering, Jiangnan University, 1800 Lihu Road, Wuxi, Jiangsu 214122 (China)
2009-04-30
This paper presents an exponential synchronization scheme for a class of neural networks with time-varying and distributed delays and reaction-diffusion terms. An adaptive synchronization controller is derived to achieve the exponential synchronization of the drive-response structure of neural networks by using the Lyapunov stability theory. At the same time, the update laws of parameters are proposed to guarantee the synchronization of delayed neural networks with all parameters unknown. It is shown that the approaches developed here extend and improve the ideas presented in recent literatures.
Aikawa, M; Goriely, S; Jorissen, A; Takahashi, K
2005-01-01
Nuclear reaction rates are quantities of fundamental importance in astrophysics. Substantial efforts have been devoted in the last decades to measure or calculate them. The present paper presents for the first time a detailed description of the Brussels nuclear reaction rate library BRUSLIB and of the nuclear network generator NETGEN so as to make these nuclear data packages easily accessible to astrophysicists for a large variety of applications. BRUSLIB is made of two parts. The first one contains the 1999 NACRE compilation based on experimental data for 86 reactions with (mainly) stable targets up to Si. The second part of BRUSLIB concerns nuclear reaction rate predictions calculated within a statistical Hauser-Feshbach approximation, which limits the reliability of the rates to reactions producing compound nuclei with a high enough level density. These calculations make use of global and coherent microscopic nuclear models for the quantities entering the rate calculations. The use of such models is utterl...
Variable elimination in chemical reaction networks with mass-action kinetics
Feliu, Elisenda; Wiuf, C.
2012-01-01
We consider chemical reaction networks taken with mass-action kinetics. The steady states of such a system are solutions to a system of polynomial equations. Even for small systems the task of finding the solutions is daunting. We develop an algebraic framework and procedure for linear elimination...
Cui Baotong [Research Center of Control Science and Engineering, Southern Yangtze University, 1800 Lihu Rd., Wuxi, Jiangsu 214122 (China)] e-mail: btcui@sohu.com; Lou Xuyang [Research Center of Control Science and Engineering, Southern Yangtze University, 1800 Lihu Rd., Wuxi, Jiangsu 214122 (China)
2006-03-01
The global asymptotic stability of bi-directional associative memory neural networks with distributed delays and reaction-diffusion terms are studied by using the analysis technique and Lyapunov functional. A sufficient condition is proposed. Two numerical examples are given to show the correctness of our analysis.
MSU SINP CDFE nuclear data activities in the nuclear reaction data centres network
This paper is the progress report of the Centre for Photonuclear Experiments Data, Moscow. It is a short review of the works carried out by the CDFE concerning the IAEA nuclear reaction data centers network activities from May 2001 until May 2002. and the description of the main results obtained. (a.n.)
An IAEA Technical Meeting of the International Network of Nuclear Reaction Data Centres was held at the IAEA Headquarters in Vienna from 25 to 26 May 2009. The meeting was attended by 23 participants from 13 cooperating data centres. A summary of the meeting is given in this report, along with the conclusions, actions, and status report of the participating data centres. (author)
Wang, Jin-Liang; Wu, Huai-Ning; Huang, Tingwen; Ren, Shun-Yan
2016-04-01
Two types of coupled neural networks with reaction-diffusion terms are considered in this paper. In the first one, the nodes are coupled through their states. In the second one, the nodes are coupled through the spatial diffusion terms. For the former, utilizing Lyapunov functional method and pinning control technique, we obtain some sufficient conditions to guarantee that network can realize synchronization. In addition, considering that the theoretical coupling strength required for synchronization may be much larger than the needed value, we propose an adaptive strategy to adjust the coupling strength for achieving a suitable value. For the latter, we establish a criterion for synchronization using the designed pinning controllers. It is found that the coupled reaction-diffusion neural networks with state coupling under the given linear feedback pinning controllers can realize synchronization when the coupling strength is very large, which is contrary to the coupled reaction-diffusion neural networks with spatial diffusion coupling. Moreover, a general criterion for ensuring network synchronization is derived by pinning a small fraction of nodes with adaptive feedback controllers. Finally, two examples with numerical simulations are provided to demonstrate the effectiveness of the theoretical results. PMID:25955852
On a theory of stability for nonlinear stochastic chemical reaction networks
We present elements of a stability theory for small, stochastic, nonlinear chemical reaction networks. Steady state probability distributions are computed with zero-information (ZI) closure, a closure algorithm that solves chemical master equations of small arbitrary nonlinear reactions. Stochastic models can be linearized around the steady state with ZI-closure, and the eigenvalues of the Jacobian matrix can be readily computed. Eigenvalues govern the relaxation of fluctuation autocorrelation functions at steady state. Autocorrelation functions reveal the time scales of phenomena underlying the dynamics of nonlinear reaction networks. In accord with the fluctuation-dissipation theorem, these functions are found to be congruent to response functions to small perturbations. Significant differences are observed in the stability of nonlinear reacting systems between deterministic and stochastic modeling formalisms
Water formation at low temperatures by surface O2 hydrogenation II: the reaction network
Cuppen, H M; Romanzin, C; Linnartz, H; 10.1039/C0CP00251H
2010-01-01
Water is abundantly present in the Universe. It is the main component of interstellar ice mantles and a key ingredient for life. Water in space is mainly formed through surface reactions. Three formation routes have been proposed in the past: hydrogenation of surface O, O2, and O3. In a previous paper [Ioppolo et al., Astrophys. J., 2008, 686, 1474] we discussed an unexpected non-standard zeroth-order H2O2 production behaviour in O2 hydrogenation experiments, which suggests that the proposed reaction network is not complete, and that the reaction channels are probably more interconnected than previously thought. In this paper we aim to derive the full reaction scheme for O2 surface hydrogenation and to constrain the rates of the individual reactions. This is achieved through simultaneous H-atom and O2 deposition under ultra-high vacuum conditions for astronomically relevant temperatures. Different H/O2 ratios are used to trace different stages in the hydrogenation network. The chemical changes in the forming ...
Global exponential stability of Hopfield reaction-diffusion neural networks with time-varying delays
王林山; 徐道义
2003-01-01
The authors analyze the existence of the equilibrium point and global exponential stability for Hopfield reaction-diffusion neural networks with time-varying delays by means of the topological degree theory and generalized Halanay inequality. Since the diffusion phenomena and time delay could not be ignored in neural networks and electric circuits, the model presented here is close to the actual systems, and the sufficient conditions on global exponential stability established in this paper,which are easily verifiable, have a wider adaptive range.
Report on the IAEA technical meeting on network of nuclear reaction data centres
This report summarizes the IAEA Technical Meeting on the Network of Nuclear Reaction Data Centres (biennial Data Centre Heads Meeting), held at the Brookhaven National Laboratory, Upton, NY, USA from 4-7 October 2004. The meeting was attended by 20 participants from 11 co-operating data centres of six Member States and two International Organizations. The report contains a summary of the meeting, the conclusions and actions, status reports of the participating data centres, and a revised technical protocol for the cooperation of the network. (author)
Coupling sample paths to the partial thermodynamic limit in stochastic chemical reaction networks
Levien, Ethan
2016-01-01
We present a new technique for reducing the variance in Monte Carlo estimators of stochastic chemical reaction networks. Our method makes use of the fact that many stochastic reaction networks converge to piecewise deterministic Markov processes in the large system-size limit. The statistics of the piecewise deterministic process can be obtained much more efficiently than those of the exact process. By coupling sample paths of the exact model to the piecewise deterministic process we are able to reduce the variance, and hence the computational complexity of the Monte Carlo estimator. In addition to rigorous results concerning the asymptotic behavior of our method, numerical simulations are performed on some simple biological models suggesting that significant computational gains are made for even moderate system-sizes.
Pd-PEG-PU polymer networks: A convenient catalyst for hydrogenation and Suzuki coupling reactions
Polyethyleneglycol-polyurethane (PEG-PU) polymer network films with sodium tetrachloropalladate were fabricated by solution casting technique. The incorporated palladium ions were reduced using sodium borohydride to obtain zero-valent (metallic) palladium in the films. These films were characterized by X-ray diffraction, Fourier transform-infra red spectroscopy, atomic absorption spectroscopy and X-ray photoelectron spectroscopy for their structure and chemical composition. The films containing palladium metal are found to be highly efficient toward hydrogenation of a number of unsaturated organic compounds and Suzuki coupling reactions of haloarenes. The novelty of this catalyst system is the incorporation of palladium in the polymer network that renders the catalyst system's reusability because there is no leaching of the metal from the polymer matrix and as the catalyst system is in the form of a compact film it can be easily separated from the reaction mixture and can be reused.
Doubly Periodic Traveling Waves in a Cellular Neural Network with Linear Reaction
Lin JianJhong
2009-01-01
Full Text Available Szekeley observed that the dynamic pattern of the locomotion of salamanders can be explained by periodic vector sequences generated by logical neural networks. Such sequences can mathematically be described by "doubly periodic traveling waves" and therefore it is of interest to propose dynamic models that may produce such waves. One such dynamic network model is built here based on reaction-diffusion principles and a complete discussion is given for the existence of doubly periodic waves as outputs. Since there are 2 parameters in our model and 4 a priori unknown parameters involved in our search of solutions, our results are nontrivial. The reaction term in our model is a linear function and hence our results can also be interpreted as existence criteria for solutions of a nontrivial linear problem depending on 6 parameters.
Henri, Christopher V.; Fernàndez-Garcia, Daniel
2015-12-01
The interplay between the spatial variability of the aquifer hydraulic properties, mass transfer due to sub-grid heterogeneity and chemical reactions often complicates reactive transport simulations. It is well documented that hydro-biochemical properties are ubiquitously heterogeneous and that diffusion and slow advection at the sub-grid scale typically leads to the conceptualization of an aquifer as a multi-porosity system. Within this context, chemical reactions taking place in mobile/immobile water regions can be substantially different between each other. This paper presents a particle-based method that can efficiently simulate heterogeneity, network reactions and multi-rate mass transfer. The approach is based on the development of transition probabilities that describe the likelihood that particles belonging to a given species and mobile/immobile domain at a given time will be transformed into another species and mobile/immobile domain afterwards. The joint effect of mass transfer and sequential degradation is shown to be non-trivial. A characteristic rebound of degradation products can be observed. This late rebound of concentrations is not driven by any change in the flow regime (e.g., pumping ceases in the pump-and-treat remediation strategy) but due to the natural interplay between mass transfer and chemical reactions. To illustrate that the method can simultaneously represent mass transfer, spatially varying properties and network reactions without numerical problems, we have simulated the degradation of tetrachloroethylene (PCE) in a three-dimensional fully heterogeneous aquifer subjected to rate-limited mass transfer. Two types of degradation modes were considered to compare the effect of an active biofilm with that of clay pods present in the aquifer. Results of the two scenarios display significantly differences. Biofilms that promote the degradation of compounds in an immobile region are shown to significantly enhance degradation, rapidly producing
INFLUENCE OF NOISE AND DELAY ON REACTION-DIFFUSION RECURRENT NEURAL NETWORKS
Li Wu
2006-01-01
In this paper, the influence of the noise and delay upon the stability property of reaction-diffusion recurrent neural networks (RNNs) with the time-varying delay is discussed. The new and easily verifiable conditions to guarantee the mean value exponential stability of an equilibrium solution are derived. The rate of exponential convergence can be estimated by means of a simple computation based on these criteria.
Report on the IAEA technical meeting of the International Network of Nuclear Reaction Data Centres
An IAEA Technical Meeting of the International Network of Nuclear Reaction Data Centres was held at IAEA Headquarters, Vienna, Austria, from 8 to 10 October 2007. The meeting was attended by 19 participants from 11 cooperating data centres of six Member States and two international organizations. A summary of the meeting is given in this report, along with the conclusions, actions, and status reports of the participating data centres. (author)
Report on the IAEA technical meeting on the network of nuclear reaction data centres
Results of the IAEA Technical meeting on the Network of Nuclear Reaction Data Centres held at the IAEA Headquarters, Vienna, Austria, 12 to 14 October 2005, are summarized in this report. The meeting was attended by 16 participants from 11 co-operating data centres of six Member States and two International Organizations. The report contains a summary of the meeting, the conclusions and actions, and status reports of the participating data centres. (author)
Product-form stationary distributions for deficiency zero chemical reaction networks
Anderson, David F.; Craciun, Gheorghe; Kurtz, Thomas G.
2008-01-01
We consider stochastically modeled chemical reaction systems with mass-action kinetics and prove that a product-form stationary distribution exists for each closed, irreducible subset of the state space if an analogous deterministically modeled system with mass-action kinetics admits a complex balanced equilibrium. Feinberg's deficiency zero theorem then implies that such a distribution exists so long as the corresponding chemical network is weakly reversible and has a deficiency of zero. The...
Report on the IAEA advisory group meeting on network of nuclear reaction data centres
This report summarizes the IAEA Advisory Group Meeting (AGM) on Network of Nuclear Reaction Data Centres, hold at the Institute of Physics and Power Engineering, Obninsk, Russia, 15 to 19 May 2000. The meeting was attended by 28 participants from 13 co-operating data centres from seven Member States and two International Organizations. The report contains a meeting summary, the conclusions and actions, progress and status reports of the participating data centres and working papers considered at the meeting. (author)
Summary Report of the Technical Meeting on International Network of Nuclear Reaction Data Centres
This report summarizes the IAEA Technical Meeting on the International Network of Nuclear Reaction Data Centres, held at the IAEA Headquarters in Vienna, Austria from 23 to 25 April 2013. The meeting was attended by 24 participants representing 13 cooperative centres from 8 Member States and 2 International Organisations. A summary of the meeting is given in this report along with the conclusions and actions. (author)
Report on the IAEA technical meeting on network of nuclear reaction data centres
An IAEA Technical Meeting on the Network of Nuclear Reaction Data Centres (and the biennial Data Centre Heads' Meeting) was held at the OECD Nuclear Energy Agency, Issy-les-Moulineaux (near Paris), France, from 27 to 30 May 2002. The meeting was attended by 21 participants from 12 co-operating data centres of six Member States and two international organizations. This report contains the meeting summary, conclusions and actions, status reports of the participating data centres, and working papers considered. (author)
Report on the IAEA technical meeting on network of nuclear reaction data centres
An IAEA Technical Meeting on the Network of Nuclear Reaction Data Centres (biennial Data Centre Heads Meeting) was held at IAEA Headquarters, Vienna, Austria, from 25 to 28 September 2006. The meeting was attended by 19 participants from 10 cooperating data centres of six Member States and two international organizations. A summary of the meeting is given in this report, along with the conclusions, actions, and status reports of the participating data centres. (author)
A reaction-diffusion model of ROS-induced ROS release in a mitochondrial network.
Lufang Zhou
2010-01-01
Full Text Available Loss of mitochondrial function is a fundamental determinant of cell injury and death. In heart cells under metabolic stress, we have previously described how the abrupt collapse or oscillation of the mitochondrial energy state is synchronized across the mitochondrial network by local interactions dependent upon reactive oxygen species (ROS. Here, we develop a mathematical model of ROS-induced ROS release (RIRR based on reaction-diffusion (RD-RIRR in one- and two-dimensional mitochondrial networks. The nodes of the RD-RIRR network are comprised of models of individual mitochondria that include a mechanism of ROS-dependent oscillation based on the interplay between ROS production, transport, and scavenging; and incorporating the tricarboxylic acid (TCA cycle, oxidative phosphorylation, and Ca(2+ handling. Local mitochondrial interaction is mediated by superoxide (O2.- diffusion and the O2.(--dependent activation of an inner membrane anion channel (IMAC. In a 2D network composed of 500 mitochondria, model simulations reveal DeltaPsi(m depolarization waves similar to those observed when isolated guinea pig cardiomyocytes are subjected to a localized laser-flash or antioxidant depletion. The sensitivity of the propagation rate of the depolarization wave to O(2.- diffusion, production, and scavenging in the reaction-diffusion model is similar to that observed experimentally. In addition, we present novel experimental evidence, obtained in permeabilized cardiomyocytes, confirming that DeltaPsi(m depolarization is mediated specifically by O2.-. The present work demonstrates that the observed emergent macroscopic properties of the mitochondrial network can be reproduced in a reaction-diffusion model of RIRR. Moreover, the findings have uncovered a novel aspect of the synchronization mechanism, which is that clusters of mitochondria that are oscillating can entrain mitochondria that would otherwise display stable dynamics. The work identifies the
Pattern formation on networks with reactions: A continuous-time random-walk approach
Angstmann, C. N.; Donnelly, I. C.; Henry, B. I.
2013-03-01
We derive the generalized master equation for reaction-diffusion on networks from an underlying stochastic process, the continuous time random walk (CTRW). The nontrivial incorporation of the reaction process into the CTRW is achieved by splitting the derivation into two stages. The reactions are treated as birth-death processes and the first stage of the derivation is at the single particle level, taking into account the death process, while the second stage considers an ensemble of these particles including the birth process. Using this model we have investigated different types of pattern formation across the vertices on a range of networks. Importantly, the CTRW defines the Laplacian operator on the network in a non-ad hoc manner and the pattern formation depends on the structure of this Laplacian. Here we focus attention on CTRWs with exponential waiting times for two cases: one in which the rate parameter is constant for all vertices and the other where the rate parameter is proportional to the vertex degree. This results in nonsymmetric and symmetric CTRW Laplacians, respectively. In the case of symmetric Laplacians, pattern formation follows from the Turing instability. However in nonsymmetric Laplacians, pattern formation may be possible with or without a Turing instability.
A unified biological modeling and simulation system for analyzing biological reaction networks
Yu, Seok Jong; Tung, Thai Quang; Park, Junho; Lim, Jongtae; Yoo, Jaesoo
2013-12-01
In order to understand the biological response in a cell, a researcher has to create a biological network and design an experiment to prove it. Although biological knowledge has been accumulated, we still don't have enough biological models to explain complex biological phenomena. If a new biological network is to be created, integrated modeling software supporting various biological models is required. In this research, we design and implement a unified biological modeling and simulation system, called ezBioNet, for analyzing biological reaction networks. ezBioNet designs kinetic and Boolean network models and simulates the biological networks using a server-side simulation system with Object Oriented Parallel Accelerator Library framework. The main advantage of ezBioNet is that a user can create a biological network by using unified modeling canvas of kinetic and Boolean models and perform massive simulations, including Ordinary Differential Equation analyses, sensitivity analyses, parameter estimates and Boolean network analysis. ezBioNet integrates useful biological databases, including the BioModels database, by connecting European Bioinformatics Institute servers through Web services Application Programming Interfaces. In addition, we employ Eclipse Rich Client Platform, which is a powerful modularity framework to allow various functional expansions. ezBioNet is intended to be an easy-to-use modeling tool and a simulation system for understanding the control mechanism by monitoring the change of each component in a biological network. The simulation result can be managed and visualized on ezBioNet, which is available free of charge at http://ezbionet.sourceforge.net or http://ezbionet.cbnu.ac.kr.
FERN – a Java framework for stochastic simulation and evaluation of reaction networks
Erhard, Florian; Friedel, Caroline C; Zimmer, Ralf
2008-01-01
Background Stochastic simulation can be used to illustrate the development of biological systems over time and the stochastic nature of these processes. Currently available programs for stochastic simulation, however, are limited in that they either a) do not provide the most efficient simulation algorithms and are difficult to extend, b) cannot be easily integrated into other applications or c) do not allow to monitor and intervene during the simulation process in an easy and intuitive way. Thus, in order to use stochastic simulation in innovative high-level modeling and analysis approaches more flexible tools are necessary. Results In this article, we present FERN (Framework for Evaluation of Reaction Networks), a Java framework for the efficient simulation of chemical reaction networks. FERN is subdivided into three layers for network representation, simulation and visualization of the simulation results each of which can be easily extended. It provides efficient and accurate state-of-the-art stochastic simulation algorithms for well-mixed chemical systems and a powerful observer system, which makes it possible to track and control the simulation progress on every level. To illustrate how FERN can be easily integrated into other systems biology applications, plugins to Cytoscape and CellDesigner are included. These plugins make it possible to run simulations and to observe the simulation progress in a reaction network in real-time from within the Cytoscape or CellDesigner environment. Conclusion FERN addresses shortcomings of currently available stochastic simulation programs in several ways. First, it provides a broad range of efficient and accurate algorithms both for exact and approximate stochastic simulation and a simple interface for extending to new algorithms. FERN's implementations are considerably faster than the C implementations of gillespie2 or the Java implementations of ISBJava. Second, it can be used in a straightforward way both as a stand
FERN – a Java framework for stochastic simulation and evaluation of reaction networks
Zimmer Ralf
2008-08-01
Full Text Available Abstract Background Stochastic simulation can be used to illustrate the development of biological systems over time and the stochastic nature of these processes. Currently available programs for stochastic simulation, however, are limited in that they either a do not provide the most efficient simulation algorithms and are difficult to extend, b cannot be easily integrated into other applications or c do not allow to monitor and intervene during the simulation process in an easy and intuitive way. Thus, in order to use stochastic simulation in innovative high-level modeling and analysis approaches more flexible tools are necessary. Results In this article, we present FERN (Framework for Evaluation of Reaction Networks, a Java framework for the efficient simulation of chemical reaction networks. FERN is subdivided into three layers for network representation, simulation and visualization of the simulation results each of which can be easily extended. It provides efficient and accurate state-of-the-art stochastic simulation algorithms for well-mixed chemical systems and a powerful observer system, which makes it possible to track and control the simulation progress on every level. To illustrate how FERN can be easily integrated into other systems biology applications, plugins to Cytoscape and CellDesigner are included. These plugins make it possible to run simulations and to observe the simulation progress in a reaction network in real-time from within the Cytoscape or CellDesigner environment. Conclusion FERN addresses shortcomings of currently available stochastic simulation programs in several ways. First, it provides a broad range of efficient and accurate algorithms both for exact and approximate stochastic simulation and a simple interface for extending to new algorithms. FERN's implementations are considerably faster than the C implementations of gillespie2 or the Java implementations of ISBJava. Second, it can be used in a straightforward
Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network
Förster, Jochen; Famili, I.; Fu, P.;
2003-01-01
environment were included. A total of 708 structural open reading frames (ORFs) were accounted for in the reconstructed network, corresponding to 1035 metabolic reactions. Further, 140 reactions were included on the basis of biochemical evidence resulting in a genome-scale reconstructed metabolic network...... Escherichia coli. The reconstructed metabolic network is the first comprehensive network for a eukaryotic organism, and it may be used as the basis for in silico analysis of phenotypic functions....
ezBioNet: A modeling and simulation system for analyzing biological reaction networks
Yu, Seok Jong; Tung, Thai Quang; Park, Junho; Lim, Jongtae; Yoo, Jaesoo
2012-10-01
To achieve robustness against living environments, a living organism is composed of complicated regulatory mechanisms ranging from gene regulations to signal transduction. If such life phenomena are to be understand, an integrated analysis tool that should have modeling and simulation functions for biological reactions, as well as new experimental methods for measuring biological phenomena, is fundamentally required. We have designed and implemented modeling and simulation software (ezBioNet) for analyzing biological reaction networks. The software can simultaneously perform an integrated modeling of various responses occurring in cells, ranging from gene expressions to signaling processes. To support massive analysis of biological networks, we have constructed a server-side simulation system (VCellSim) that can perform ordinary differential equations (ODE) analysis, sensitivity analysis, and parameter estimates. ezBioNet integrates the BioModel database by connecting the european bioinformatics institute (EBI) servers through Web services APIs and supports the handling of systems biology markup language (SBML) files. In addition, we employed eclipse RCP (rich client platform) which is a powerful modularity framework allowing various functional expansions. ezBioNet is intended to be an easy-to-use modeling tool, as well as a simulation system, to understand the control mechanism by monitoring the change of each component in a biological network. A researcher may perform the kinetic modeling and execute the simulation. The simulation result can be managed and visualized on ezBioNet, which is freely available at http://ezbionet.cbnu.ac.kr.
Kaihong Zhao
2010-01-01
Full Text Available The existence of equilibrium solutions to reaction-diffusion recurrent neural networks with Dirichlet boundary conditions on time scales is proved by the topological degree theory and M-matrix method. Under some sufficient conditions, we obtain the uniqueness and global exponential stability of equilibrium solution to reaction-diffusion recurrent neural networks with Dirichlet boundary conditions on time scales by constructing suitable Lyapunov functional and inequality skills. One example is given to illustrate the effectiveness of our results.
Grytskyy, Dmytro; Diesmann, Markus; Helias, Moritz
2016-06-01
Self-organized structures in networks with spike-timing dependent synaptic plasticity (STDP) are likely to play a central role for information processing in the brain. In the present study we derive a reaction-diffusion-like formalism for plastic feed-forward networks of nonlinear rate-based model neurons with a correlation sensitive learning rule inspired by and being qualitatively similar to STDP. After obtaining equations that describe the change of the spatial shape of the signal from layer to layer, we derive a criterion for the nonlinearity necessary to obtain stable dynamics for arbitrary input. We classify the possible scenarios of signal evolution and find that close to the transition to the unstable regime metastable solutions appear. The form of these dissipative solitons is determined analytically and the evolution and interaction of several such coexistent objects is investigated.
Leporace, Gustavo; Batista, Luiz Alberto; Metsavaht, Leonardo; Nadal, Jurandir
2015-08-01
The aim of the study was to analyze and compare the residuals obtained from ground reaction force (GRF) models developed using two different neural network configurations (one network with three outputs; and three networks with one output each), based on accelerometer data. Seventeen healthy subjects walked along a walkway, with a force plate embedded, with a three dimensional accelerometer attached to the shank. Multilayer perceptron networks (MLP) models were developed with the 3D accelerometer data as inputs to predict the GRF. The residuals of these models were evaluated graphically and numerically to verify the fitting. A visual analysis of the simulated signals suggests the model was able to adequately predict the GRF. The errors and correlations found in the MLP models for the 3D GRF is at least similar to other studies, although some of them showed higher errors. There was not difference between the two MLP configurations. However, despite the high correlation coefficient and closeness to a normal probability distribution, the residual analysis still presented a higher kurtosis and skewness, suggesting that the inclusion of other variables and the increase of the validation sample size could increase the fitting of the simulation. PMID:26736876
Baldi, P.; Blanke, Mogens; Castaldi, P.; Mimmo, N.; Simani, S.
2015-01-01
This paper suggests a novel diagnosis scheme for detection, isolation and estimation of faults affecting satellite reaction wheels. Both spin rate measurements and actuation torque defects are dealt with. The proposed system consists of a fault detection and isolation module composed by a bank of...... residual filters organized in a generalized scheme, followed by a fault estimation module consisting of a bank of adaptive estimation filters. The residuals are decoupled from aerodynamic disturbances thanks to the Nonlinear Geometric Approach. The use of Radial Basis Function Neural Networks is shown to...... allow design of generalized fault estimation filters, which do not need a priori information about the faults internal model. Simulation results with a detailed nonlinear spacecraft model, which includes disturbances, show that the proposed diagnosis scheme can deal with faults affecting both reaction...
Elucidation of Diels-Alder Reaction Network of 2,5-Dimethylfuran and Ethylene on HY Zeolite Catalyst
Do, Phuong T. M. [Univ. of Delaware, Newark, DE (United States); McAtee, Jesse R. [Univ. of Delaware, Newark, DE (United States); Watson, Donald A. [Univ. of Delaware, Newark, DE (United States); Lobo, Raul F. [Univ. of Delaware, Newark, DE (United States)
2012-12-12
The reaction of 2,5-dimethylfuran and ethylene to produce p-xylene represents a potentially important route for the conversion of biomass to high-value organic chemicals. Current preparation methods suffer from low selectivity and produce a number of byproducts. Using modern separation and analytical techniques, the structures of many of the byproducts produced in this reaction when HY zeolite is employed as a catalyst have been identified. From these data, a detailed reaction network is proposed, demonstrating that hydrolysis and electrophilic alkylation reactions compete with the desired Diels–Alder/dehydration sequence. This information will allow the rational identification of more selective catalysts and more selective reaction conditions.
Ignat Drozdov
Full Text Available Small intestinal (SI neuroendocrine tumors (NET are increasing in incidence, however little is known about their biology. High throughput techniques such as inference of gene regulatory networks from microarray experiments can objectively define signaling machinery in this disease. Genome-wide co-expression analysis was used to infer gene relevance network in SI-NETs. The network was confirmed to be non-random, scale-free, and highly modular. Functional analysis of gene co-expression modules revealed processes including 'Nervous system development', 'Immune response', and 'Cell-cycle'. Importantly, gene network topology and differential expression analysis identified over-expression of the GPCR signaling regulators, the cAMP synthetase, ADCY2, and the protein kinase A, PRKAR1A. Seven CREB response element (CRE transcripts associated with proliferation and secretion: BEX1, BICD1, CHGB, CPE, GABRB3, SCG2 and SCG3 as well as ADCY2 and PRKAR1A were measured in an independent SI dataset (n = 10 NETs; n = 8 normal preparations. All were up-regulated (p<0.035 with the exception of SCG3 which was not differently expressed. Forskolin (a direct cAMP activator, 10(-5 M significantly stimulated transcription of pCREB and 3/7 CREB targets, isoproterenol (a selective ß-adrenergic receptor agonist and cAMP activator, 10(-5 M stimulated pCREB and 4/7 targets while BIM-53061 (a dopamine D(2 and Serotonin [5-HT(2] receptor agonist, 10(-6 M stimulated 100% of targets as well as pCREB; CRE transcription correlated with the levels of cAMP accumulation and PKA activity; BIM-53061 stimulated the highest levels of cAMP and PKA (2.8-fold and 2.5-fold vs. 1.8-2-fold for isoproterenol and forskolin. Gene network inference and graph topology analysis in SI NETs suggests that SI NETs express neural GPCRs that activate different CRE targets associated with proliferation and secretion. In vitro studies, in a model NET cell system, confirmed that transcriptional
Kim, Daesang
2013-11-01
We have extended reactive flow simulation in pore-network models to include geometric changes in the medium from dissolution effects. These effects include changes in pore volume and reactive surface area, as well as topological changes that open new connections. The computed changes were based upon a mineral map from an X-ray computed tomography image of a sandstone core. We studied the effect of these changes on upscaled (pore-scale to core-scale) reaction rates and compared against the predictions of a continuum model. Specifically, we modeled anorthite and kaolinite reactions under acidic flow conditions during which the anorthite reactions remain far from equilibrium (dissolution only), while the kaolinite reactions can be near-equilibrium. Under dissolution changes, core-scale reaction rates continuously and nonlinearly evolved in time. At higher injection rates, agreement with predictions of the continuum model degraded significantly. For the far-from-equilibrium reaction, our results indicate that the ability to correctly capture the heterogeneity in dissolution changes in the reactive mineral surface area is critical to accurately predict upscaled reaction rates. For the near-equilibrium reaction, the ability to correctly capture the heterogeneity in the saturation state remains critical. Inclusion of a Nernst-Planck term to ensure neutral ionic currents under differential diffusion resulted in at most a 9% correction in upscaled rates.
Das, Biswajit; Gangopadhyay, Gautam, E-mail: gautam@bose.res.in [S. N. Bose National Centre For Basic Sciences, Block-JD, Sector-III, Salt Lake, Kolkata 700 098 (India); Banerjee, Kinshuk [Department of Chemistry, University of Calcutta, 92 A.P.C. Road, Kolkata 700 009 (India)
2013-12-28
In this work, we develop an approach to nonequilibrium thermodynamics of an open chemical reaction network in terms of the elementary reaction propensities. The method is akin to the microscopic formulation of the dissipation function in terms of the Kullback-Leibler distance of phase space trajectories in Hamiltonian system. The formalism is applied to a single oligomeric enzyme kinetics at chemiostatic condition that leads the reaction system to a nonequilibrium steady state, characterized by a positive total entropy production rate. Analytical expressions are derived, relating the individual reaction contributions towards the total entropy production rate with experimentally measurable reaction velocity. Taking a real case of Escherichia coli β-galactosidase enzyme obeying Michaelis-Menten kinetics, we thoroughly analyze the temporal as well as the steady state behavior of various thermodynamic quantities for each elementary reaction. This gives a useful insight in the relative magnitudes of various energy terms and the dissipated heat to sustain a steady state of the reaction system operating far-from-equilibrium. It is also observed that, the reaction is entropy-driven at low substrate concentration and becomes energy-driven as the substrate concentration rises.
Wadén Johan
2009-10-01
Full Text Available Background Cardiovascular disease is the main cause of premature death in patients with type 1 diabetes. Patients with diabetic kidney disease have an increased risk of heart attack or stroke. Accurate knowledge of the complex inter-dependencies between the risk factors is critical for pinpointing the best targets for research and treatment. Therefore, the aim of this study was to describe the association patterns between clinical and biochemical features of diabetic complications. Methods Medical records and serum and urine samples of 4,197 patients with type 1 diabetes were collected from health care centers in Finland. At baseline, the mean diabetes duration was 22 years, 52% were male, 23% had kidney disease (urine albumin excretion over 300 mg/24 h or end-stage renal disease and 8% had a history of macrovascular events. All-cause mortality was evaluated after an average of 6.5 years of follow-up (25,714 patient years. The dataset comprised 28 clinical and 25 biochemical variables that were regarded as the nodes of a network to assess their mutual relationships. Results The networks contained cliques that were densely inter-connected (r > 0.6, including cliques for high-density lipoprotein (HDL markers, for triglycerides and cholesterol, for urinary excretion and for indices of body mass. The links between the cliques showed biologically relevant interactions: an inverse relationship between HDL cholesterol and the triglyceride clique (r P -16, a connection between triglycerides and body mass via C-reactive protein (r > 0.3, P -16 and intermediate-density cholesterol as the connector between lipoprotein metabolism and albuminuria (r > 0.3, P -16. Aging and macrovascular disease were linked to death via working ability and retinopathy. Diabetic kidney disease, serum creatinine and potassium, retinopathy and blood pressure were inter-connected. Blood pressure correlations indicated accelerated vascular aging in individuals with kidney disease
Evolution of Autocatalytic Sets in Computational Models of Chemical Reaction Networks
Hordijk, Wim
2016-06-01
Several computational models of chemical reaction networks have been presented in the literature in the past, showing the appearance and (potential) evolution of autocatalytic sets. However, the notion of autocatalytic sets has been defined differently in different modeling contexts, each one having some shortcoming or limitation. Here, we review four such models and definitions, and then formally describe and analyze them in the context of a mathematical framework for studying autocatalytic sets known as RAF theory. The main results are that: (1) RAF theory can capture the various previous definitions of autocatalytic sets and is therefore more complete and general, (2) the formal framework can be used to efficiently detect and analyze autocatalytic sets in all of these different computational models, (3) autocatalytic (RAF) sets are indeed likely to appear and evolve in such models, and (4) this could have important implications for a possible metabolism-first scenario for the origin of life.
Evolution of Autocatalytic Sets in Computational Models of Chemical Reaction Networks.
Hordijk, Wim
2016-06-01
Several computational models of chemical reaction networks have been presented in the literature in the past, showing the appearance and (potential) evolution of autocatalytic sets. However, the notion of autocatalytic sets has been defined differently in different modeling contexts, each one having some shortcoming or limitation. Here, we review four such models and definitions, and then formally describe and analyze them in the context of a mathematical framework for studying autocatalytic sets known as RAF theory. The main results are that: (1) RAF theory can capture the various previous definitions of autocatalytic sets and is therefore more complete and general, (2) the formal framework can be used to efficiently detect and analyze autocatalytic sets in all of these different computational models, (3) autocatalytic (RAF) sets are indeed likely to appear and evolve in such models, and (4) this could have important implications for a possible metabolism-first scenario for the origin of life. PMID:26499126
A model for lignin alteration - Part I: A kinetic reaction-network model
Payne, D.F.; Ortoleva, P.J.
2001-01-01
A new quantitative model is presented which simulates the maturation of lignin-derived sedimentary organic matter under geologic conditions. In this model, compositionally specific reactants evolve to specific intermediate and mobile products through balanced, nth order processes, by way of a network of sequential and parallel reactions. The chemical kinetic approach is based primarily on published observed structural transformations of naturally matured, lignin-derived, sedimentary organic matter. Assuming that Upper Cretaceous Williams Fork coal in the Piceance Basin is primarily lignin-derived, the model is calibrated for the Multi-Well Experiment(MWX) Site in this basin. This kind of approach may be applied to other selectively preserved chemical components of sedimentary organic matter. ?? 2001 Elsevier Science Ltd. All rights reserved.
Parlett, Christopher M. A.; Isaacs, Mark A.; Beaumont, Simon K.; Bingham, Laura M.; Hondow, Nicole S.; Wilson, Karen; Lee, Adam F.
2016-02-01
The chemical functionality within porous architectures dictates their performance as heterogeneous catalysts; however, synthetic routes to control the spatial distribution of individual functions within porous solids are limited. Here we report the fabrication of spatially orthogonal bifunctional porous catalysts, through the stepwise template removal and chemical functionalization of an interconnected silica framework. Selective removal of polystyrene nanosphere templates from a lyotropic liquid crystal-templated silica sol-gel matrix, followed by extraction of the liquid crystal template, affords a hierarchical macroporous-mesoporous architecture. Decoupling of the individual template extractions allows independent functionalization of macropore and mesopore networks on the basis of chemical and/or size specificity. Spatial compartmentalization of, and directed molecular transport between, chemical functionalities affords control over the reaction sequence in catalytic cascades; herein illustrated by the Pd/Pt-catalysed oxidation of cinnamyl alcohol to cinnamic acid. We anticipate that our methodology will prompt further design of multifunctional materials comprising spatially compartmentalized functions.
Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network
Förster, Jochen; Famili, I.; Fu, P.; Palsson, B.O.; Nielsen, Jens
2003-01-01
The metabolic network in the yeast Saccharomyces cerevisiae was reconstructed using currently available genomic, biochemical, and physiological information. The metabolic reactions were compartmentalized between the cytosol and the mitochondria, and transport steps between the compartments and the...... containing 1175 metabolic reactions and 584 metabolites. The number of gene functions included in the reconstructed network corresponds to similar to16% of all characterized ORFs in S. cerevisiae. Using the reconstructed network, the metabolic capabilities of S. cerevisiae were calculated and compared with...
Singh, Shalini; Samal, Areejit; Giri, Varun; Krishna, Sandeep; Raghuram, Nandula; Jain, Sanjay
2012-01-01
Unraveling the structure of complex biological networks and relating it to their functional role is an important task in systems biology. Here we attempt to characterize the functional organization of the large-scale metabolic networks of three microorganisms. We apply flux balance analysis to study the optimal growth states of these organisms in different environments. By investigating the differential usage of reactions across flux patterns for different environments, we observe a striking ...