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.

A framework for deriving analytic steady states of biochemical reaction networks.

The long-term behaviors of biochemical systems are often described by their steady states. Deriving these states directly for complex networks arising from real-world applications, however, is often challenging. Recent work has consequently focused on network-based approaches. Specifically, biochemical reaction networks are transformed into weakly reversible and deficiency zero generalized networks, which allows the derivation of their analytic steady states. Identifying this transformation, however, can be challenging for large and complex networks. In this paper, we address this difficulty by breaking the complex network into smaller independent subnetworks and then transforming the subnetworks to derive the analytic steady states of each subnetwork. We show that stitching these solutions together leads to the analytic steady states of the original network. To facilitate this process, we develop a user-friendly and publicly available package, COMPILES (COMPutIng anaLytic stEady States). With COMPILES, we can easily test the presence of bistability of a CRISPRi toggle switch model, which was previously investigated via tremendous number of numerical simulations and within a limited range of parameters. Furthermore, COMPILES can be used to identify absolute concentration robustness (ACR), the property of a system that maintains the concentration of particular species at a steady state regardless of any initial concentrations. Specifically, our approach completely identifies all the species with and without ACR in a complex insulin model. Our method provides an effective approach to analyzing and understanding complex biochemical systems.

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.