BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Eduardo R. Mendoza

Eduardo R. Mendoza contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.DSMolecular Networksmath.AG

Trust snapshot

Quick read

Trust 21 - Emerging
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Positive equilibria of power law kinetics on networks with independent linkage classes

Studies about the set of positive equilibria ($E_+$) of kinetic systems have been focused on mass action, and not that much on power law kinetic (PLK) systems, even for PL-RDK systems (PLK systems where two reactions with identical reactant complexes have the same kinetic order vectors). For mass action, reactions with different reactants have different kinetic order rows. A PL-RDK system satisfying this property is called factor span surjective (PL-FSK). In this work, we show that a cycle terminal PL-FSK system with $E_+\ne \varnothing$ and has independent linkage classes (ILC) is a poly-PLP system, i.e., $E_+$ is the disjoint union of log-parametrized sets. The key insight for the extension is that factor span surjectivity induces an isomorphic digraph structure on the kinetic complexes. The result also completes, for ILC networks, the structural analysis of the original complex balanced generalized mass action systems (GMAS) by Müller and Regensburger. We also identify a large set of PL-RDK systems where non-emptiness of $E_+$ is a necessary and sufficient condition for non-emptiness of each set of positive equilibria for each linkage class. These results extend those of Boros on mass action systems with ILC. We conclude this paper with two applications of our results. Firstly, we consider absolute complex balancing (ACB), i.e., the property that each positive equilibrium is complex balanced, in poly-PLP systems. Finally, we use the new results to study absolute concentration robustness (ACR) in these systems. In particular, we obtain a species hyperplane containment criterion to determine ACR in the system species.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Weakly Reversible CF-Decompositions of Chemical Kinetic Systems

This paper studies chemical kinetic systems which decompose into weakly reversible complex factorizable (CF) systems. Among power law kinetic systems, CF systems (denoted as PL-RDK systems) are those where branching reactions of a reactant complex have identical rows in the kinetic order matrix. Mass action and generalized mass action systems (GMAS) are well-known examples. Schmitz's global carbon cycle model is a previously studied non-complex factorizable (NF) power law system (denoted as PL-NDK). We derive novel conditions for the existence of weakly reversible CF-decompositions and present an algorithm for verifying these conditions. We discuss methods for identifying independent decompositions, i.e., those where the stoichiometric subspaces of the subnetworks form a direct sum, as such decompositions relate positive equilibria sets of the subnetworks to that of the whole network. We then use the results to determine the positive equilibria sets of PL-NDK systems which admit an independent weakly reversible decomposition into PL-RDK systems of PLP type, i.e., the positive equilibria are log-parametrized, which is a broad generalization of a Deficiency Zero Theorem of Fortun et al. (2019).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Computational Approach to Multistationarity in Poly-PL Kinetic Systems

One important question that interests those who work in chemical reaction network theory (CRNT) is this: Does the system obtained from a reaction network admit a positive equilibrium and if it does, can there be more than one within a stoichiometric class? The higher deficiency algorithm (HDA) of Ji and Feinberg provided a method of determining the multistationarity capacity of a CRN with mass action kinetics (MAK). An extension of this, called Multistationarity Algorithm (MSA), recently came into the scene tackling CRNs with power law kinetics (PLK), a kinetic system which is more general (having MAK systems as a special case). For this paper, we provide a computational approach to study the multistationarity feature of reaction networks endowed with kinetics which are non-negative linear combinations of power law functions called poly-PL kinetics (PYK). The idea is to use MSA and combine it with a transformation called STAR-MSC (i.e., $S$-invariant Termwise Addition of Reactions via Maximal Stoichiometric Coefficients) producing PLKs that are dynamically equivalent to PYKs. This leads us to being able to determinine the multistationarity capacity of a much larger class of kinetic systems. We show that if the transformed dynamically equivalent PLK system is multistationary for a stoichiometric class for a set of particular rate constants, then so is its original corresponding PYK system. Moreover, the monostationarity property of the transformed PLK system also implies the monostationarity property of the original PYK system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Computational Approach to Multistationarity of Power-Law Kinetic Systems

This paper presents a computational solution to determine if a chemical reaction network endowed with power-law kinetics (PLK system) has the capacity for multistationarity, i.e., whether there exist positive rate constants such that the corresponding differential equations admit multiple positive steady states within a stoichiometric class. The approach, which is called the "Multistationarity Algorithm for PLK systems" (MSA), combines (i) the extension of the "higher deficiency algorithm" of Ji and Feinberg for mass action to PLK systems with reactant-determined interactions, and (ii) a method that transforms any PLK system to a dynamically equivalent one with reactant-determined interactions. Using this algorithm, we obtain two new results: the monostationarity of a popular model of anaerobic yeast fermentation pathway, and the multistationarity of a global carbon cycle model with climate engineering, both in the generalized mass action format of biochemical systems theory. We also provide examples of the broader scope of our approach for deficiency one PLK systems in comparison to the extension of Feinberg's "deficiency one algorithm" to such systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Absolute concentration robustness in power law kinetic systems

Absolute concentration robustness (ACR) is a condition wherein a species in a chemical kinetic system possesses the same value for any positive steady state the network may admit regardless of initial conditions. Thus far, results on ACR center on chemical kinetic systems with deficiency one. In this contribution, we use the idea of dynamic equivalence of chemical reaction networks to derive novel results that guarantee ACR for some classes of power law kinetic systems with deficiency zero. Furthermore, using network decomposition, we identify ACR in higher deficiency networks (i.e. deficiency $\geq$ 2) by considering the presence of a low deficiency subnetwork with ACR. Network decomposition also enabled us to recognize and define a weaker form of concentration robustness than ACR, which we named as `balanced concentration robustness'. Finally, we also discuss and emphasize our view of ACR as a primarily kinetic character rather than a condition that arises from structural sources.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chemical reaction network decompositions and realizations of S-systems

This paper presents novel decomposition classes of chemical reaction networks (CRNs) derived from S-system kinetics. Based on the network decomposition theory initiated by Feinberg in 1987, we introduce the concept of incidence independent decompositions and develop the theory of $\mathscr{C}$- and $\mathscr{C}^*$- decompositions which partition the set of complexes and the set of nonzero complexes respectively, including their structure theorems in terms of linkage classes. Analogous to Feinberg's independent decomposition, we demonstrate the important relationship between sets of complex balance equilibria for an incidence independent decomposition of weakly reversible subnetworks for any kinetics. We show that the $\mathscr{C}^*$-decompositions are also incidence independent. We also introduce in this paper a new realization for an S-system that is analyzed using a newly defined class of species coverable CRNs. This led to the extension of the deficiency formula and characterization of fundamental decompositions of species decomposable reaction networks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fundamental Decompositions and Multistationarity of Power-Law Kinetic Systems

The fundamental decomposition of a chemical reaction network (also called its "$\mathscr{F}$-decomposition") is the set of subnetworks generated by the partition of its set of reactions into the "fundamental classes" introduced by Ji and Feinberg in 2011 as the basis of their "higher deficiency algorithm" for mass action systems. The first part of this paper studies the properties of the $\mathscr{F}$-decomposition, in particular, its independence (i.e., the network's stoichiometric subspace is the direct sum of the subnetworks' stoichiometric subspaces) and its incidence-independence (i.e., the image of the network's incidence map is the direct sum of the incidence maps' images of the subnetworks). We derive necessary and sufficient conditions for these properties and identify network classes where the $\mathscr{F}$-decomposition coincides with other known decompositions. The second part of the paper applies the above-mentioned results to improve the Multistationarity Algorithm for power-law kinetic systems (MSA), a general computational approach that we introduced in previous work. We show that for systems with non-reactant determined interactions but with an independent $\mathscr{F}$-decomposition, the transformation to a dynamically equivalent system with reactant-determined interactions -- required in the original MSA -- is not necessary. We illustrate this improvement with the subnetwork of Schmitz's carbon cycle model recently analyzed by Fortun et al.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Positive Equilibria of Hill-Type Kinetic Systems

This work introduces a novel approach to study properties of positive equilibria of a chemical reaction network $\mathscr{N}$ endowed with Hill-type kinetics $K$, called a Hill-type kinetic (HTK) system $\left(\mathscr{N},K\right)$, including their multiplicity and concentration robustness in a species. We associate a unique positive linear combination of power-law kinetic systems called poly-PL kinetic (PYK) system $\left( {\mathscr{N},{K_\text{PY}}} \right)$ to the given HTK system. The associated system has the key property that its equilibria sets coincide with those of the Hill-type system, i.e., ${E_ + }\left( {\mathscr{N},K} \right) = {E_ + }\left( {\mathscr{N},{K_\text{PY}}} \right)$ and ${Z_ + }\left( {\mathscr{N},K} \right) = {Z_ + }\left( {\mathscr{N},{K_\text{PY}}} \right)$. This allows us to identify two novel subsets of the Hill-type kinetics, called PL-equilibrated and PL-complex balanced kinetics, to which recent results on absolute concentration robustness (ACR) of species and complex balancing at positive equilibria of power-law (PL) kinetic systems can be applied. Our main results also include the Shinar-Feinberg ACR Theorem for PL-equilibrated HT-RDK systems (i.e., subset of complex factorizable HTK systems), which establishes a foundation for the analysis of ACR in HTK systems, and the extension of the results of Müller and Regensburger on generalized mass action systems to PL-complex balanced HT-RDK systems. In addition, we derive the theory of balanced concentration robustness (BCR) in an analogous manner to ACR for PL-equilibrated systems. Finally, we provide further extensions of our results to a more general class of kinetics, which includes quotients of poly-PL functions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Robustness in power law kinetic systems with reactant-determined interactions

Robustness against the presence of environmental disruptions can be observed in many systems of chemical reaction network. However, identifying the underlying components of a system that give rise to robustness is often elusive. The influential work of Shinar and Feinberg established simple yet subtle network-based conditions for absolute concentration robustness (ACR), a phenomena in which a species in a mass-action system has the same concentration for any steady state the network may admit. In this contribution, we extend this result to embrace kinetic systems more general than mass-action systems, namely, power-law kinetic systems with reactant-determined interactions (denoted by "PL-RDK"). In PL-RDK, the kinetic order vectors (which we call "interactions") of reactions with the same reactant complex are identical. As illustration, we considered a scenario in the pre-industrial state of global carbon cycle. A power-law approximation of the dynamical system of this scenario is found to be dynamically equivalent to an ACR-possessing PL-RDK system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structural Properties of an S-system Model of Mycobacterium Tuberculosis Gene Regulation

Magombedze and Mulder in 2013 studied the gene regulatory system of Mycobacterium Tuberculosis (Mtb) by partitioning this into three subsystems based on putative gene function and role in dormancy/latency development. Each subsystem, in the form of S-system, is represented by an embedded chemical reaction network (CRN), defined by a species subset and a reaction subset induced by the set of digraph vertices of the subsystem. For the embedded networks of S-system, we showed interesting structural properties and proved that all S-system CRNs (with at least two species) are discordant. Analyzing the subsystems as subnetworks, where arcs between vertices belonging to different subsystems are retained, we formed a digraph homomorphism from the corresponding subnetworks to the embedded networks. Lastly, we explored the modularity concept of CRN in the context of digraph.

0 citations0 reviewsNo code linkedOpen access