Source author record

Badal Joshi

Badal Joshi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

11 published item(s)

preprint2022arXiv

Multistationarity in cyclic sequestration-transmutation networks

We consider a natural class of reaction networks which consist of reactions where either two species can inactivate each other (i.e., sequestration), or some species can be transformed into another (i.e., transmutation), in a way that gives rise to a feedback cycle. We completely characterize the capacity of multistationarity of these networks. This is especially interesting because such networks provide simple examples of "atoms of multistationarity", i.e., minimal networks that can give rise to multiple positive steady states

preprint2022arXiv

Reaction Network Motifs for Static and Dynamic Absolute Concentration Robustness

Networks with absolute concentration robustness (ACR) have the property that a translation of a coordinate hyperplane either contains all steady states (static ACR) or attracts all trajectories (dynamic ACR). The implication for the underlying biological system is robustness in the concentration of one of the species independent of the initial conditions as well as independent of the concentration of all other species. Identifying network conditions for dynamic ACR is a challenging problem. We lay the groundwork in this paper by studying small reaction networks, those with 2 reactions and 2 species. We give a complete classification by ACR properties of these minimal reaction networks. The dynamics is rich even within this simple setting. Insights obtained from this work will help illuminate the properties of more complex networks with dynamic ACR.

preprint2022arXiv

Transition graph decomposition for complex balanced reaction networks with non-mass-action kinetics

Reaction networks are widely used models to describe biochemical processes. Stochastic fluctuations in the counts of biological macromolecules have amplified consequences due to their small population sizes. This makes it necessary to favor stochastic, discrete population, continuous time models. The stationary distributions provide snapshots of the model behavior at the stationary regime, and as such finding their expression in terms of the model parameters is of great interest. The aim of the present paper is to describe when the stationary distributions of the original model, whose state space is potentially infinite, coincide exactly with the stationary distributions of the process truncated to finite subsets of states, up to a normalizing constant. The finite subsets of states we identify are called copies and are inspired by the modular topology of reaction network models. With such a choice we prove a novel graphical characterization of the concept of complex balancing for stochastic models of reaction networks. The results of the paper hold for the commonly used mass-action kinetics but are not restricted to it, and are in fact stated for more general setting.

preprint2021arXiv

On reaction network implementations of neural networks

This paper is concerned with the utilization of deterministically modeled chemical reaction networks for the implementation of (feed-forward) neural networks. We develop a general mathematical framework and prove that the ordinary differential equations (ODEs) associated with certain reaction network implementations of neural networks have desirable properties including (i) existence of unique positive fixed points that are smooth in the parameters of the model (necessary for gradient descent), and (ii) fast convergence to the fixed point regardless of initial condition (necessary for efficient implementation). We do so by first making a connection between neural networks and fixed points for systems of ODEs, and then by constructing reaction networks with the correct associated set of ODEs. We demonstrate the theory by constructing a reaction network that implements a neural network with a smoothed ReLU activation function, though we also demonstrate how to generalize the construction to allow for other activation functions (each with the desirable properties listed previously). As there are multiple types of "networks" utilized in this paper, we also give a careful introduction to both reaction networks and neural networks, in order to disambiguate the overlapping vocabulary in the two settings and to clearly highlight the role of each network's properties.

preprint2016arXiv

Which small reaction networks are multistationary?

Reaction networks taken with mass-action kinetics arise in many settings, from epidemiology to population biology to systems of chemical reactions. Bistable reaction networks are posited to underlie biochemical switches, which motivates the following question: which reaction networks have the capacity for multiple steady states? Mathematically, this asks: among certain parametrized families of polynomial systems, which admit multiple positive roots? No complete answer is known. This work analyzes the smallest networks, those with only a few chemical species or reactions. For these "smallest" networks, we completely answer the question of multistationarity and, in some cases, multistability too, thereby extending related work of Boros. Our results highlight the role played by the Newton polytope of a network (the convex hull of the reactant vectors). Also, our work is motivated by recent results that explain how a given network's capacity for multistationarity arises from that of certain related networks which are typically smaller. Hence, we are interested in classifying small multistationary networks, and our work forms a first step in this direction.

preprint2015arXiv

A survey of methods for deciding whether a reaction network is multistationary

Which reaction networks, when taken with mass-action kinetics, have the capacity for multiple steady states? There is no complete answer to this question, but over the last 40 years various criteria have been developed that can answer this question in certain cases. This work surveys these developments, with an emphasis on recent results that connect the capacity for multistationarity of one network to that of another. In this latter setting, we consider a network $N$ that is embedded in a larger network $G$, which means that $N$ is obtained from $G$ by removing some subsets of chemical species and reactions. This embedding relation is a significant generalization of the subnetwork relation. For arbitrary networks, it is not true that if $N$ is embedded in $G$, then the steady states of $N$ lift to $G$. Nonetheless, this does hold for certain classes of networks; one such class is that of fully open networks. This motivates the search for embedding-minimal multistationary networks: those networks which admit multiple steady states but no proper, embedded networks admit multiple steady states. We present results about such minimal networks, including several new constructions of infinite families of these networks.

preprint2014arXiv

A detailed balanced reaction network is sufficient but not necessary for its Markov chain to be detailed balanced

Certain chemical reaction networks (CRNs) when modeled as a deterministic dynamical system taken with mass-action kinetics have the property of reaction network detailed balance (RNDB) which is achieved by imposing network-related constraints on the reaction rate constants. Markov chains (whether arising as models of CRNs or otherwise) have their own notion of detailed balance, imposed by the network structure of the graph of the transition matrix of the Markov chain. When considering Markov chains arising from chemical reaction networks with mass-action kinetics, we will refer to this property as Markov chain detailed balance (MCDB). Finally, we refer to the stochastic analog of RNDB as Whittle stochastic detailed balance (WSDB). It is known that RNDB and WSDB are equivalent. We prove that WSDB and MCDB are also intimately related but are not equivalent. While RNDB implies MCDB, the converse is not true. The conditions on rate constants that result in networks with MCDB but without RNDB are stringent, and thus examples of this phenomenon are rare, a notable exception is a network whose Markov chain is a birth and death process. We give a new algorithm to find conditions on the rate constants that are required for MCDB.

preprint2013arXiv

Complete characterization by multistationarity of fully open networks with one non-flow reaction

This article characterizes certain small multistationary chemical reaction networks. We consider the set of fully open networks, those for which all chemical species participate in inflow and outflow, containing one non-flow (reversible or irreversible) reaction. We show that such a network admits multiple positive mass-action steady states if and only if the stoichiometric coefficients in the non-flow reaction satisfy a certain simple arithmetic relation. The multistationary fully open one-reaction networks are identified with the chemical process of autocatalysis. Using the notion of `embedded network' defined recently by Joshi and Shiu, we provide new sufficient conditions for establishing multistationarity of fully open networks, applicable well beyond the one-reaction setting.

preprint2012arXiv

Atoms of multistationarity in chemical reaction networks

Chemical reaction systems are dynamical systems that arise in chemical engineering and systems biology. In this work, we consider the question of whether the minimal (in a precise sense) multistationary chemical reaction networks, which we propose to call `atoms of multistationarity,' characterize the entire set of multistationary networks. Our main result states that the answer to this question is `yes' in the context of fully open continuous-flow stirred-tank reactors (CFSTRs), which are networks in which all chemical species take part in the inflow and outflow. In order to prove this result, we show that if a subnetwork admits multiple steady states, then these steady states can be lifted to a larger network, provided that the two networks share the same stoichiometric subspace. We also prove an analogous result when a smaller network is obtained from a larger network by `removing species.' Our results provide the mathematical foundation for a technique used by Siegal-Gaskins et al. of establishing bistability by way of `network ancestry.' Additionally, our work provides sufficient conditions for establishing multistationarity by way of atoms and moreover reduces the problem of classifying multistationary CFSTRs to that of cataloging atoms of multistationarity. As an application, we enumerate and classify all 386 bimolecular and reversible two-reaction networks. Of these, exactly 35 admit multiple positive steady states. Moreover, each admits a unique minimal multistationary subnetwork, and these subnetworks form a poset (with respect to the relation of `removing species') which has 11 minimal elements (the atoms of multistationarity).

preprint2012arXiv

Simplifying the Jacobian Criterion for precluding multistationarity in chemical reaction networks

Chemical reaction networks taken with mass-action kinetics are dynamical systems that arise in chemical engineering and systems biology. In general, determining whether a chemical reaction network admits multiple steady states is difficult, as this requires determining existence of multiple positive solutions to a large system of polynomials with unknown coefficients. However, in certain cases, various easy criteria can be applied. One such test is the Jacobian Criterion, due to Craciun and Feinberg, which gives sufficient conditions for ruling out the possibility of multiple steady states. A chemical reaction network is said to pass the Jacobian Criterion if all terms in the determinant expansion of its parametrized Jacobian matrix have the same sign. In this article, we present a procedure which simplifies the application of the Jacobian Criterion, and as a result, we identify a new class of networks for which multiple steady states is precluded: those in which all chemical species have total molecularity of at most two. The total molecularity of a species refers to the sum of all of its stoichiometric coefficients in the network. We illustrate our results by examining enzyme catalysis networks.

preprint2011arXiv

Order of magnitude time-reversible Markov chains and characterization of clustering processes

We introduce the notion of order of magnitude reversibility (OM-reversibility) in Markov chains that are parametrized by a positive parameter $\ep$. OM-reversibility is a weaker condition than reversibility, and requires only the knowledge of order of magnitude of the transition probabilities. For an irreducible, OM-reversible Markov chain on a finite state space, we prove that the stationary distribution satisfies order of magnitude detailed balance (analog of detailed balance in reversible Markov chains). The result characterizes the states with positive probability in the limit of the stationary distribution as $\ep \to 0$, which finds an important application in the case of singularly perturbed Markov chains that are reducible for $\ep=0$. We show that OM-reversibility occurs naturally in macroscopic systems, involving many interacting particles. Clustering is a common phenomenon in biological systems, in which particles or molecules aggregate at one location. We give a simple condition on the transition probabilities in an interacting particle Markov chain that characterizes clustering. We show that such clustering processes are OM-reversible, and we find explicitly the order of magnitude of the stationary distribution. Further, we show that the single pole states, in which all particles are at a single vertex, are the only states with positive probability in the limit of the stationary distribution as the rate of diffusion goes to zero.