Stability of open pathways

An in silico approach in identification of drug targets in Leishmania: A subtractive genomic and metabolic simulation analysis

Leishmaniasis is one of the major health issue in developing countries. The current therapeutic regimen for this disease is less effective with lot of adverse effects thereby warranting an urgent need to develop not only new and selective drug candidates but also identification of effective drug targets. Here we present subtractive genomics procedure for identification of putative drug targets in Leishmania. Comprehensive druggability analysis has been carried out in the current work for identified metabolic pathways and drug targets. We also demonstrate effective metabolic simulation methodology to pinpoint putative drug targets in threonine biosynthesis pathway. Metabolic simulation data from the current study indicate that decreasing flux through homoserine kinase reaction can be considered as a good therapeutic opportunity. The data from current study is expected to show new avenue for designing experimental strategies in search of anti-leishmanial agents.

Topological and kinetic determinants of the modal matrices of dynamic models of metabolism

Large-scale kinetic models of metabolism are becoming increasingly comprehensive and accurate. A key challenge is to understand the biochemical basis of the dynamic properties of these models. Linear analysis methods are well-established as useful tools for characterizing the dynamic response of metabolic networks. Central to linear analysis methods are two key matrices: the Jacobian matrix (J) and the modal matrix (M-1) arising from its eigendecomposition. The modal matrix M-1 contains dynamically independent motions of the kinetic model near a reference state, and it is sparse in practice for metabolic networks. However, connecting the structure of M-1 to the kinetic properties of the underlying reactions is non-trivial. In this study, we analyze the relationship between J, M-1, and the kinetic properties of the underlying network for kinetic models of metabolism. Specifically, we describe the origin of mode sparsity structure based on features of the network stoichiometric matrix S and the reaction kinetic gradient matrix G. First, we show that due to the scaling of kinetic parameters in real networks, diagonal dominance occurs in a substantial fraction of the rows of J, resulting in simple modal structures with clear biological interpretations. Then, we show that more complicated modes originate from topologically-connected reactions that have similar reaction elasticities in G. These elasticities represent dynamic equilibrium balances within reactions and are key determinants of modal structure. The work presented should prove useful towards obtaining an understanding of the dynamics of kinetic models of metabolism, which are rooted in the network structure and the kinetic properties of reactions.

Qualitative Analysis of an ODE Model of a Class of Enzymatic Reactions : Some Results on Global Stability of Messenger RNA-MicroRNA Interaction

The present paper analyzes an ODE model of a certain class of (open) enzymatic reactions. This type of model is used, for instance, to describe the interactions between messenger RNAs and microRNAs. It is shown that solutions defined by positive initial conditions are well defined and bounded on [Formula: see text] and that the positive octant of [Formula: see text] is a positively invariant set. We prove further that in this positive octant there exists a unique equilibrium point, which is asymptotically stable and a global attractor for any initial state with positive components; a controllability property is emphasized. We also investigate the qualitative behavior of the QSSA system in the phase plane [Formula: see text]. For this planar system we obtain similar results regarding global stability by using Lyapunov theory, invariant regions and controllability.

Steady states and stability in metabolic networks without regulation

Metabolic networks are often extremely complex. Despite intensive efforts many details of these networks, e.g., exact kinetic rates and parameters of metabolic reactions, are not known, making it difficult to derive their properties. Considerable effort has been made to develop theory about properties of steady states in metabolic networks that are valid for any values of parameters. General results on uniqueness of steady states and their stability have been derived with specific assumptions on reaction kinetics, stoichiometry and network topology. For example, deep results have been obtained under the assumptions of mass-action reaction kinetics, continuous flow stirred tank reactors (CFSTR), concordant reaction networks and others. Nevertheless, a general theory about properties of steady states in metabolic networks is still missing. Here we make a step further in the quest for such a theory. Specifically, we study properties of steady states in metabolic networks with monotonic kinetics in relation to their stoichiometry (simple and general) and the number of metabolites participating in every reaction (single or many). Our approach is based on the investigation of properties of the Jacobian matrix. We show that stoichiometry, network topology, and the number of metabolites that participate in every reaction have a large influence on the number of steady states and their stability in metabolic networks. Specifically, metabolic networks with single-substrate-single-product reactions have disconnected steady states, whereas in metabolic networks with multiple-substrates-multiple-product reactions manifolds of steady states arise. Metabolic networks with simple stoichiometry have either a unique globally asymptotically stable steady state or asymptotically stable manifolds of steady states. In metabolic networks with general stoichiometry the steady states are not always stable and we provide conditions for their stability. In order to demonstrate the biological relevance we illustrate the results on the examples of the TCA cycle, the mevalonate pathway and the Calvin cycle.

Global stability of enzymatic chains of full reversible Michaelis-Menten reactions

We consider a chain of metabolic reactions catalyzed by enzymes, of reversible Michaelis-Menten type with full dynamics, i.e. not reduced with any quasi-steady state approximations. We study the corresponding dynamical system and show its global stability if the equilibrium exists. If the system is open, the equilibrium may not exist. The main tool is monotone systems theory. Finally we study the implications of these results for the study of coupled genetic-metabolic systems.

Global stability of reversible enzymatic metabolic chains

We consider metabolic networks with reversible enzymatic reactions. The model is written as a system of ordinary differential equations, possibly with inputs and outputs. We prove the global stability of the equilibrium (if it exists), using techniques of monotone systems and compartmental matrices. We show that the equilibrium does not always exist. Finally, we consider a metabolic system coupled with a genetic network, and we study the dependence of the metabolic equilibrium (if it exists) with respect to concentrations of enzymes. We give some conclusions concerning the dynamical behavior of coupled genetic/metabolic systems.