Comparative analysis of kinetic realizations of insulin signaling.

Several studies have developed dynamical models to understand the underlying mechanisms of insulin signaling, a signaling cascade that leads to the translocation of glucose, the human body's main source of energy. Fortunately, reaction network analysis allows us to extract properties of dynamical systems without depending on their model parameter values. This study focuses on the comparison of insulin signaling in healthy state (INSMS or INSulin Metabolic Signaling) and in type 2 diabetes (INRES or INsulin RESistance) using reaction network analysis. The analysis uses network decomposition to identify the different subsystems involved in insulin signaling (e.g., insulin receptor binding and recycling, GLUT4 translocation, and ERK signaling pathway, among others). Furthermore, results show that INSMS and INRES are similar with respect to some network, structo-kinetic, and kinetic properties. Their differences, however, provide insights into what happens when insulin resistance occurs. First, the variation in the number of species involved in INSMS and INRES suggests that when irregularities occur in the insulin signaling pathway, other complexes (and, hence, other processes) get involved, characterizing insulin resistance. Second, the loss of concordance exhibited by INRES suggests less restrictive interplay between the species involved in insulin signaling, leading to unusual activities in the signaling cascade. Lastly, GLUT4 losing its absolute concentration robustness in INRES may signify that the transporter has lost its reliability in shuttling glucose to the cell, inhibiting efficient cellular energy production. This study also suggests possible applications of the equilibria parametrization and network decomposition, resulting from the analysis, to potentially establish absolute concentration robustness in a species.

Reaction Network Analysis of Metabolic Insulin Signaling.

Absolute concentration robustness (ACR) and concordance are novel concepts in the theory of robustness and stability within Chemical Reaction Network Theory. In this paper, we have extended Shinar and Feinberg's reaction network analysis approach to the insulin signaling system based on recent advances in decomposing reaction networks. We have shown that the network with 20 species, 35 complexes, and 35 reactions is concordant, implying at most one positive equilibrium in each of its stoichiometric compatibility class. We have obtained the system's finest independent decomposition consisting of 10 subnetworks, a coarsening of which reveals three subnetworks which are not only functionally but also structurally important. Utilizing the network's deficiency-oriented coarsening, we have developed a method to determine positive equilibria for the entire network. Our analysis has also shown that the system has ACR in 8 species all coming from a deficiency zero subnetwork. Interestingly, we have shown that, for a set of rate constants, the insulin-regulated glucose transporter GLUT4 (important in glucose energy metabolism), has stable ACR.

Linear conjugacy of chemical kinetic systems.

Two networks are said to be linearly conjugate if the solution of their dynamic equations can be transformed into each other by a positive linear transformation. The study on dynamical equivalence in chemical kinetic systems was initiated by Craciun and Pantea in 2008 and eventually led to the Johnston-Siegel Criterion for linear conjugacy (JSC). Several studies have applied Mixed Integer Linear Programming (MILP) approach to generate linear conjugates of MAK (mass action kinetic) systems, Bio-CRNs (which is a subset of Hill-type kinetic systems when the network is restricted to digraphs), and PL-RDK (complex factorizable power law kinetic) systems. In this study, we present a general computational solution to construct linear conjugates of any "rate constant-interaction function decomposable" (RID) chemical kinetic systems, wherein each of its rate function is the product of a rate constant and an interaction function. We generate an extension of the JSC to the complex factorizable (CF) subset of RID kinetic systems and show that any non-complex factorizable (NF) RID kinetic system can be dynamically equivalent to a CF system via transformation. We show that linear conjugacy can be generated for any RID kinetic systems by applying the JSC to any NF kinetic system that are transformed to CF kinetic system.

13C based proteinogenic amino acid (PAA) and metabolic flux ratio analysis of Lactococcus lactis reveals changes in pentose phosphate (PP) pathway in response to agitation and temperature related stresses.

Lactococcus lactis subsp. cremoris MG1363 is an important starter culture for dairy fermentation. During industrial fermentations, L. lactis is constantly exposed to stresses that affect the growth and performance of the bacterium. Although the response of L. lactis to several stresses has been described, the adaptation mechanisms at the level of in vivo fluxes have seldom been described. To gain insights into cellular metabolism, 13C metabolic flux analysis and gas chromatography mass spectrometry (GC-MS) were used to measure the flux ratios of active pathways in the central metabolism of L. lactis when subjected to three conditions varying in temperature (30°C, 37°C) and agitation (with and without agitation at 150 rpm). Collectively, the concentrations of proteinogenic amino acids (PAAs) and free fatty acids (FAAs) were compared, and Pearson correlation analysis (r) was calculated to measure the pairwise relationship between PAAs. Branched chain and aromatic amino acids, threonine, serine, lysine and histidine were correlated strongly, suggesting changes in flux regulation in glycolysis, the pentose phosphate (PP) pathway, malic enzyme and anaplerotic reaction catalysed by pyruvate carboxylase (pycA). Flux ratio analysis revealed that glucose was mainly converted by glycolysis, highlighting the stability of L. lactis' central carbon metabolism despite different conditions. Higher flux ratios through oxaloacetate (OAA) from pyruvate (PYR) reaction in all conditions suggested the activation of pyruvate carboxylate (pycA) in L. lactis, in response to acid stress during exponential phase. Subsequently, more significant flux ratio differences were seen through the oxidative and non-oxidative pentose phosphate (PP) pathways, malic enzyme, and serine and C1 metabolism, suggesting NADPH requirements in response to environmental stimuli. These reactions could play an important role in optimization strategies for metabolic engineering in L. lactis. Overall, the integration of systematic analysis of amino acids and flux ratio analysis provides a systems-level understanding of how L. lactis regulates central metabolism under various conditions.

Reaction networks and kinetics of biochemical systems.

This paper further develops the connection between Chemical Reaction Network Theory (CRNT) and Biochemical Systems Theory (BST) that we recently introduced [1]. We first use algebraic properties of kinetic sets to study the set of complex factorizable kinetics CFK(N) on a CRN, which shares many characteristics with its subset of mass action kinetics. In particular, we extend the Theorem of Feinberg-Horn [9] on the coincidence of the kinetic and stoichiometric subsets of a mass action system to CF kinetics, using the concept of span surjectivity. We also introduce the branching type of a network, which determines the availability of kinetics on it and allows us to characterize the networks for which all kinetics are complex factorizable: A "Kinetics Landscape" provides an overview of kinetics sets, their algebraic properties and containment relationships. We then apply our results and those (of other CRNT researchers) reviewed in [1] to fifteen BST models of complex biological systems and discover novel network and kinetic properties that so far have not been widely studied in CRNT. In our view, these findings show an important benefit of connecting CRNT and BST modeling efforts.

Model Construction and Analysis of Respiration in Halobacterium salinarum.

The archaeon Halobacterium salinarum can produce energy using three different processes, namely photosynthesis, oxidative phosphorylation and fermentation of arginine, and is thus a model organism in bioenergetics. Compared to its bacteriorhodopsin-driven photosynthesis, less attention has been devoted to modeling its respiratory pathway. We created a system of ordinary differential equations that models its oxidative phosphorylation. The model consists of the electron transport chain, the ATP synthase, the potassium uniport and the sodium-proton antiport. By fitting the model parameters to experimental data, we show that the model can explain data on proton motive force generation, ATP production, and the charge balancing of ions between the sodium-proton antiporter and the potassium uniport. We performed sensitivity analysis of the model parameters to determine how the model will respond to perturbations in parameter values. The model and the parameters we derived provide a resource that can be used for analytical studies of the bioenergetics of H. salinarum.

Chemical reaction network approaches to Biochemical Systems Theory.

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.

Flux duality in nonlinear GMA systems: implications for metabolic engineering.

Pathway models in biotechnology are customarily designed with metabolite concentrations as the primary, dependent variables, whereas fluxes are derived quantities that are secondarily computed from the primary variables. In other fields of mathematics, such as graph theory and linear systems analysis, it has proven useful to complement primal network model designs in terms of vertices (pools) and edges (processes) with dual designs, where the roles of the primary and secondary quantities are interchanged. The conversion from primal to dual systems is fairly easy in linear systems, but it is unclear to what degree it is possible in nonlinear systems. In this article, we present a method to transform nonlinear primal models of pathways or other biotechnological processes that conform to the Generalized Mass Action (GMA) structure within the formalism of Biochemical Systems Theory (BST) into dual models that focus on fluxes, rather than pools. Interestingly, this transformation is relatively straightforward, once some notational issues are streamlined, and the resulting dual system is again in the format of a GMA system. The transformation is illustrated with the example of glycolysis in Saccharomyces cerevisiae, a well known organism in the food industry. The results suggest how rewriting a model in terms of its fluxes helps bridge the gap between flux balance and dynamic models and also offers a view of the investigated system that is complementary to that of the original model. This dual view is important, because fluxes are sometimes more relevant for the behavior of a biotechnological system than (metabolite) pools. The dual system furthermore offers a systematic approach toward understanding the dynamical constraints under which the system operates.

Metabolic Engineering with power-law and linear-logarithmic systems.

Metabolic Engineering aims to improve the performance of biotechnological processes through rational manipulation rather than random mutagenesis of the organisms involved. Such a strategy can only succeed when a mathematical model of the target process is available. Simplifying assumptions are often needed to cope with the complexity of such models in an efficient way, and the choice of such assumptions often leads to models that fall within a certain structural template or formalism. The most popular formalisms can be grouped in two categories: power-law and linear-logarithmic. As optimization and analysis of a model strongly depends on its structure, most methods in Metabolic Engineering have been defined within a given formalism and never used in any other. In this work, the four most commonly used formalisms (two power-law and two linear-logarithmic) are placed in a common framework defined within Biochemical Systems Theory. This framework defines every model as matrix equations in terms of the same parameters, enabling the formulation of a common steady state analysis and providing means for translating models and methods from one formalism to another. Several Metabolic Engineering methods are analysed here and shown to be variants of a single equation. Particularly, two problem solving philosophies are compared: the application of the design equation and the solution of constrained optimization problems. Generalizing the design equation to all the formalisms shows it to be interchangeable with the direct solution of the rate law in matrix form. Furthermore, optimization approaches are concluded to be preferable since they speed the exploration of the feasible space, implement a better specification of the problem and exclude unrealistic results. Beyond consolidating existing knowledge and enabling comparison, the systematic approach adopted here can fill the gaps between the different methods and combine their strengths.