BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matthew D. Johnston

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

ResearcherUnclaimed source record
math.DSMolecular Networksmath.NAphysics.chem-phmath.CAmath.OCphysics.soc-phPopulations and Evolution

Catalog footprint

What is connected

13works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

13 published item(s)

preprint2020arXiv

A Dynamical Framework for Modeling Fear of Infection and Frustration with Social Distancing in COVID-19 Spread

In this paper, we introduce a novel modeling framework for incorporating fear of infection and frustration with social distancing into disease dynamics. We show that the resulting SEIR behavior-perception model has three principal modes of qualitative behavior---no outbreak, controlled outbreak, and uncontrolled outbreak. We also demonstrate that the model can produce transient and sustained waves of infection consistent with secondary outbreaks. We fit the model to cumulative COVID-19 case and mortality data from several regions. Our analysis suggests that regions which experience a significant decline after the first wave of infection, such as Canada and Israel, are more likely to contain secondary waves of infection, whereas regions which only achieve moderate success in mitigating the disease's spread initially, such as the United States, are likely to experience substantial secondary waves or uncontrolled outbreaks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On classes of reaction networks and their associated polynomial dynamical systems

In the study of reaction networks and the polynomial dynamical systems that they generate, special classes of networks with important properties have been identified. These include reversible, weakly reversible}, and, more recently, endotactic networks. While some inclusions between these network types are clear, such as the fact that all reversible networks are weakly reversible, other relationships are more complicated. Adding to this complexity is the possibility that inclusions be at the level of the dynamical systems generated by the networks rather than at the level of the networks themselves. We completely characterize the inclusions between reversible, weakly reversible, endotactic, and strongly endotactic network, as well as other less well studied network types. In particular, we show that every strongly endotactic network in two dimensions can be generated by an extremally weakly reversible network. We also introduce a new class of source-only networks, which is a computationally convenient property for networks to have, and show how this class relates to the above mentioned network types.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Linear Programming Approach to Dynamical Equivalence, Linear Conjugacy, and the Deficiency One Theorem

The well-known Deficiency One Theorem gives structural conditions on a chemical reaction network under which, for any set of parameter values, the steady states of the corresponding mass action system may be easily characterized. It is also known, however, that mass action systems are not uniquely associated with reaction networks and that some representations may satisfy the Deficiency One Theorem while others may not. In this paper we present a mixed-integer linear programming framework capable of determining whether a given mass action system has a dynamically equivalent or linearly conjugate representation which has an underlying network satisfying the Deficiency One Theorem. This extends recent computational work determining linearly conjugate systems which are weakly reversible and have a deficiency of zero.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A computational approach to persistence, permanence, and endotacticity of biochemical reaction systems

We introduce a mixed-integer linear programming (MILP) framework capable of determining whether a chemical reaction network possesses the property of being endotactic or strongly endotactic. The network property of being strongly endotactic is known to lead to persistence and permanence of chemical species under genetic kinetic assumptions, while the same result is conjectured but as yet unproved for general endotactic networks. The algorithms we present are the first capable of verifying endotacticity of chemical reaction networks for systems with greater than two constituent species. We implement the algorithms in the open-source online package CoNtRol and apply them to several well-studied biochemical examples, including the general $n$-site phosphorylation / dephosphorylation networks and a circadian clock mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Computational Approach to Steady State Correspondence of Regular and Generalized Mass Action Systems

It has been recently observed that the dynamical properties of mass action systems arising from many models of biochemical reaction networks can be derived by considering the corresponding properties of a related generalized mass action system. The correspondence process known as network translation in particular has been shown to be useful in characterizing a system's steady states. In this paper, we further develop the theory of network translation with particular focus on a subclass of translations known as improper translations. For these translations, we derive conditions on the network topology of the translated network which are sufficient to guarantee the original and translated systems share the same steady states. We then present a mixed-integer linear programming (MILP) algorithm capable of determining whether a mass action system can be corresponded to a generalized system through the process of network translation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A note on "MAPK networks and their capacity for multistationarity due to toric steady states"

We provide a short supplement to the paper "MAPK networks and their capacity for multistationarity due to toric steady states" by Pérez Millán and Turjanski. We show that the capacity for toric steady states in the three networks analyzed in that paper can be derived using the process of network translation, which corresponds the original mass action system to a generalized mass action system with the same steady states. In all three cases, the translated chemical reaction network is proper, weakly reversible, and has both a structural and kinetic deficiency of zero. This is sufficient to guarantee toric steady states by previously established work on network translations. A basis of the steady state ideal is then derived by consideration of the linkage classes of the translated chemical reaction network.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Translated Chemical Reaction Networks

Many biochemical and industrial applications involve complicated networks of simultaneously occurring chemical reactions. Under the assumption of mass action kinetics, the dynamics of these chemical reaction networks are governed by systems of polynomial ordinary differential equations. The steady states of these mass action systems have been analysed via a variety of techniques, including elementary flux mode analysis, algebraic techniques (e.g. Groebner bases), and deficiency theory. In this paper, we present a novel method for characterizing the steady states of mass action systems. Our method explicitly links a network's capacity to permit a particular class of steady states, called toric steady states, to topological properties of a related network called a translated chemical reaction network. These networks share their reaction stoichiometries with their source network but are permitted to have different complex stoichiometries and different network topologies. We apply the results to examples drawn from the biochemical literature.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Computing Weakly Reversible Linearly Conjugate Chemical Reaction Networks with Minimal Deficiency

Mass-action kinetics is frequently used in systems biology to model the behaviour of interacting chemical species. Many important dynamical properties are known to hold for such systems if they are weakly reversible and have a low deficiency. In particular, the Deficiency Zero and Deficiency One Theorems guarantee strong regularity with regards to the number and stability of positive equilibrium states. It is also known that chemical reaction networks with disparate reaction structure can exhibit the same qualitative dynamics. The theory of linear conjugacy encapsulates the cases where this relationship is captured by a linear transformation. In this paper, we propose a mixed-integer linear programming algorithm capable of determining weakly reversible reaction networks with a minimal deficiency which are linearly conjugate to a given reaction network.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Linear Programming Approach to Weak Reversibility and Linear Conjugacy of Chemical Reaction Networks

A numerically effective procedure for determining weakly reversible chemical reaction networks that are linearly conjugate to a known reaction network is proposed in this paper. The method is based on translating the structural and algebraic characteristics of weak reversibility to logical statements and solving the obtained set of linear (in)equalities in the framework of mixed integer linear programming. The unknowns in the problem are the reaction rate coefficients and the parameters of the linear conjugacy transformation. The efficacy of the approach is shown through numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical Equivalence and Linear Conjugacy of Chemical Reaction Networks: New Results and Methods

In the first part of this paper, we propose new optimization-based methods for the computation of preferred (dense, sparse, reversible, detailed and complex balanced) linearly conjugate reaction network structures with mass action dynamics. The developed methods are extensions of previously published results on dynamically equivalent reaction networks and are based on mixed-integer linear programming. As related theoretical contributions we show that (i) dense linearly conjugate networks define a unique super-structure for any positive diagonal state transformation if the set of chemical complexes is given, and (ii) the existence of linearly conjugate detailed balanced and complex balanced networks do not depend on the selection of equilibrium points. In the second part of the paper it is shown that determining dynamically equivalent realizations to a network that is structurally fixed but parametrically not can also be written and solved as a mixed-integer linear programming problem. Several examples illustrate the presented computation methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Linear Conjugacy of Chemical Reaction Networks

Under suitable assumptions, the dynamic behaviour of a chemical reaction network is governed by an autonomous set of polynomial ordinary differential equations over continuous variables representing the concentrations of the reactant species. It is known that two networks may possess the same governing mass-action dynamics despite disparate network structure. To date, however, there has only been limited work exploiting this phenomenon even for the cases where one network possesses known dynamics while the other does not. In this paper, we bring these known results into a broader unified theory which we call conjugate chemical reaction network theory. We present a theorem which gives conditions under which two networks with different governing mass-action dynamics may exhibit the same qualitative dynamics and use it to extend the scope of the well-known theory of weakly reversible systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Weak Dynamic Non-Emptiability and Persistence of Chemical Kinetics Systems

A frequently desirable characteristic of chemical kinetics systems is that of persistence, the property that if all the species are initially present then none of them may tend toward extinction. It is known that solutions of deterministically modelled mass-action systems may only approach portions of the boundary of the positive orthant which correspond to semi-locking sets (alternatively called siphons). Consequently, most recent work on persistence of these systems has been focused on these sets. In this paper, we focus on a result which states that, for a conservative mass-action system, persistence holds if every critical semi-locking set is dynamically non-emptiable and the system contains no nested locking sets. We will generalize this result by introducing the notion of a weakly dynamically non-emptiable semi-locking set and making novel use of the well-known Farkas' Lemma. We will also connect this result to known results regarding complex balanced systems and systems with facets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Stratum Approach to Global Stability of Complex Balanced Systems

It has long been known that complex balanced mass-action systems exhibit a restrictive form of behaviour known as locally stable dynamics. This means that within each compatibility class $\mathcal{C}_{\mathbf{x}_0}$---the forward invariant space where solutions lies---there is exactly one equilibrium concentration and that this concentration is locally asymptotically stable. It has also been conjectured that this stability extends globally to $\mathcal{C}_{\mathbf{x}_0}$. That is to say, all solutions originating in $\mathcal{C}_{\mathbf{x}_0}$ approach the unique positive equilibrium concentration rather than $\partial \mathcal{C}_{\mathbf{x}_0}$ or $\infty$. To date, however, no general proof of this conjecture has been found. In this paper, we approach the problem of global stability for complex balanced systems through the methodology of dividing the positive orthant into regions called strata. This methodology has been previously applied to detailed balanced systems---a proper subset of complex balanced systems---to show that, within a stratum, trajectories are repelled from any face of $\mathbb{R}_{\geq 0}^m$ adjacent to the stratum. Several known global stability results for detailed balanced systems are generalized to complex balanced systems.

0 citations0 reviewsNo code linkedOpen access