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

Anne Shiu

Anne Shiu 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.AGmath.CONeurons and Cognitionmath.ACmath.CASystems and Control

Catalog footprint

What is connected

18works
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

18 published item(s)

preprint2022arXiv

Absolute concentration robustness in networks with low-dimensional stoichiometric subspace

A reaction system exhibits "absolute concentration robustness" (ACR) in some species if the positive steady-state value of that species does not depend on initial conditions. Mathematically, this means that the positive part of the variety of the steady-state ideal lies entirely in a hyperplane of the form $x_i=c$, for some $c>0$. Deciding whether a given reaction system -- or those arising from some reaction network -- exhibits ACR is difficult in general, but here we show that for many simple networks, assessing ACR is straightforward. Indeed, our criteria for ACR can be performed by simply inspecting a network or its standard embedding into Euclidean space. Our main results pertain to networks with many conservation laws, so that all reactions are parallel to one other. Such "one-dimensional" networks include those networks having only one species. We also consider networks with only two reactions, and show that ACR is characterized by a well-known criterion of Shinar and Feinberg. Finally, up to some natural ACR-preserving operations -- relabeling species, lengthening a reaction, and so on -- only three families of networks with two reactions and two species have ACR. Our results are proven using algebraic and combinatorial techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nondegenerate Neural Codes and Obstructions to Closed-Convexity

Previous work on convexity of neural codes has produced codes that are open-convex but not closed-convex -- or vice-versa. However, why a code is one but not the other, and how to detect such discrepancies are open questions. We tackle these questions in two ways. First, we investigate the concept of degeneracy introduced by Cruz et al., and extend their results to show that nondegeneracy precisely captures the situation when taking closures or interiors of open or closed realizations, respectively, yields another realization of the code. Second, we give the first general criteria for precluding a code from being closed-convex (without ruling out open-convexity), unifying ad-hoc geometric arguments in prior works. One criterion is built on a phenomenon we call a rigid structure, while the other can be stated algebraically, in terms of the neural ideal of the code. These results complement existing criteria having the opposite purpose: precluding open-convexity but not closed-convexity. Finally, we show that a family of codes shown by Jeffs to be not open-convex is in fact closed-convex and realizable in dimension two.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamics of ERK regulation in the processive limit

We consider a model of extracellular signal-regulated kinase (ERK) regulation by dual-site phosphorylation and dephosphorylation, which exhibits bistability and oscillations, but loses these properties in the limit in which the mechanisms underlying phosphorylation and dephosphorylation become processive. Our results suggest that anywhere along the way to becoming processive, the model remains bistable and oscillatory. More precisely, in simplified versions of the model, precursors to bistability and oscillations (specifically, multistationarity and Hopf bifurcations, respectively) exist at all "processivity levels". Finally, we investigate whether bistability and oscillations can exist together.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Identifiability of Linear Compartmental Models: The Effect of Moving Inputs, Outputs, and Leaks

A mathematical model is identifiable if its parameters can be recovered from data. Here we investigate, for linear compartmental models, whether (local, generic) identifiability is preserved when parts of the model -- specifically, inputs, outputs, leaks, and edges -- are moved, added, or deleted. Our results are as follows. First, for certain catenary, cycle, and mammillary models, moving or deleting the leak preserves identifiability. Next, for cycle models with up to one leak, moving inputs or outputs preserves identifiability. Thus, every cycle model with up to one leak (and at least one input and at least one output) is identifiable. Next, we give conditions under which adding leaks renders a cycle model unidentifiable. Finally, for certain cycle models with no leaks, adding specific edges again preserves identifiability. Our proofs, which are algebraic and combinatorial in nature, rely on results on elementary symmetric polynomials and the theory of input-output equations for linear compartmental models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neural Codes With Three Maximal Codewords: Convexity and Minimal Embedding Dimension

Neural codes, represented as collections of binary strings called codewords, are used to encode neural activity. A code is called convex if its codewords are represented as an arrangement of convex open sets in Euclidean space. Previous work has focused on addressing the question: how can we tell when a neural code is convex? Giusti and Itskov identified a local obstruction and proved that convex neural codes have no local obstructions. The converse is true for codes on up to four neurons, but false in general. Nevertheless, we prove this converse holds for codes with up to three maximal codewords, and moreover the minimal embedding dimension of such codes is at most two.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Analyzing Multistationarity in Chemical Reaction Networks using the Determinant Optimization Method

Multistationary chemical reaction networks are of interest to scientists and mathematicians alike. While some criteria for multistationarity exist, obtaining explicit reaction rates and steady states that exhibit multistationarity for a given network -- in order to check nondegeneracy or determine stability of the steady states, for instance -- is nontrivial. Nonetheless, we accomplish this task for a certain family of sequestration networks. Additionally, our results allow us to prove the existence of nondegenerate steady states for some of these sequestration networks, thereby resolving a subcase of a conjecture of Joshi and Shiu. Our work relies on the determinant optimization method, developed by Craciun and Feinberg, for asserting that certain networks are multistationary. More precisely, we implement the construction of reaction rates and multiple steady states which appears in the proofs that underlie their method. Furthermore, we describe in detail the steps of this construction so that other researchers can more easily obtain, as we did, multistationary rates and steady states.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Obstructions to convexity in neural codes

How does the brain encode spatial structure? One way is through hippocampal neurons called place cells, which become associated to convex regions of space known as their receptive fields: each place cell fires at a high rate precisely when the animal is in the receptive field. The firing patterns of multiple place cells form what is known as a convex neural code. How can we tell when a neural code is convex? To address this question, Giusti and Itskov identified a local obstruction, defined via the topology of a code's simplicial complex, and proved that convex neural codes have no local obstructions. Curto et al. proved the converse for all neural codes on at most four neurons. Via a counterexample on five neurons, we show that this converse is false in general. Additionally, we classify all codes on five neurons with no local obstructions. This classification is enabled by our enumeration of connected simplicial complexes on 5 vertices up to isomorphism. Finally, we examine how local obstructions are related to maximal codewords (maximal sets of neurons that co-fire). Curto et al. proved that a code has no local obstructions if and only if it contains certain "mandatory" intersections of maximal codewords. We give a new criterion for an intersection of maximal codewords to be non-mandatory, and prove that it classifies all such non-mandatory codewords for codes on up to 5 neurons.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

What makes a neural code convex?

Neural codes allow the brain to represent, process, and store information about the world. Combinatorial codes, comprised of binary patterns of neural activity, encode information via the collective behavior of populations of neurons. A code is called convex if its codewords correspond to regions defined by an arrangement of convex open sets in Euclidean space. Convex codes have been observed experimentally in many brain areas, including sensory cortices and the hippocampus, where neurons exhibit convex receptive fields. What makes a neural code convex? That is, how can we tell from the intrinsic structure of a code if there exists a corresponding arrangement of convex open sets? In this work, we provide a complete characterization of local obstructions to convexity. This motivates us to define max intersection-complete codes, a family guaranteed to have no local obstructions. We then show how our characterization enables one to use free resolutions of Stanley-Reisner ideals in order to detect violations of convexity. Taken together, these results provide a significant advance in understanding the intrinsic combinatorial properties of convex codes.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A global convergence result for processive multisite phosphorylation systems

Multisite phosphorylation plays an important role in intracellular signaling. There has been much recent work aimed at understanding the dynamics of such systems when the phosphorylation/dephosphorylation mechanism is distributive, that is, when the binding of a substrate and an enzyme molecule results in addition or removal of a single phosphate group and repeated binding therefore is required for multisite phosphorylation. In particular, such systems admit bistability. Here we analyze a different class of multisite systems, in which the binding of a substrate and an enzyme molecule results in addition or removal of phosphate groups at all phosphorylation sites. That is, we consider systems in which the mechanism is processive, rather than distributive. We show that in contrast with distributive systems, processive systems modeled with mass-action kinetics do not admit bistability and, moreover, exhibit rigid dynamics: each invariant set contains a unique equilibrium, which is a global attractor. Additionally, we obtain a monomial parametrization of the steady states. Our proofs rely on a technique of Johnston for using "translated" networks to study systems with "toric steady states", recently given sign conditions for injectivity of polynomial maps, and a result from monotone systems theory due to Angeli and Sontag.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sign conditions for injectivity of generalized polynomial maps with applications to chemical reaction networks and real algebraic geometry

We give necessary and sufficient conditions in terms of sign vectors for the injectivity of families of polynomial maps with arbitrary real exponents defined on the positive orthant. Our work relates and extends existing injectivity conditions expressed in terms of Jacobian matrices and determinants. In the context of chemical reaction networks with power-law kinetics, our results can be used to preclude as well as to guarantee multiple positive steady states. In the context of real algebraic geometry,our work recognizes a prior result of Craciun, Garcia-Puente, and Sottile, together with work of two of the authors, as the first partial multivariate generalization of the classical Descartes' rule, which bounds the number of positive real roots of a univariate real polynomial in terms of the number of sign variations of its coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A geometric approach to the Global Attractor Conjecture

This paper introduces the class of "strongly endotactic networks", a subclass of the endotactic networks introduced by G. Craciun, F. Nazarov, and C. Pantea. The main result states that the global attractor conjecture holds for complex-balanced systems that are strongly endotactic: every trajectory with positive initial condition converges to the unique positive equilibrium allowed by conservation laws. This extends a recent result by D. F. Anderson for systems where the reaction diagram has only one linkage class (connected component). The results here are proved using differential inclusions, a setting that includes power-law systems. The key ideas include a perspective on reaction kinetics in terms of combinatorial geometry of reaction diagrams, a projection argument that enables analysis of a given system in terms of systems with lower dimension, and an extension of Birch's theorem, a well-known result about intersections of affine subspaces with manifolds parameterized by monomials.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Projection Argument for Differential Inclusions, with Applications to Persistence of Mass-Action Kinetics

Motivated by questions in mass-action kinetics, we introduce the notion of vertexical family of differential inclusions. Defined on open hypercubes, these families are characterized by particular good behavior under projection maps. The motivating examples are certain families of reaction networks -- including reversible, weakly reversible, endotactic, and strongly endotactic reaction networks -- that give rise to vertexical families of mass-action differential inclusions. We prove that vertexical families are amenable to structural induction. Consequently, a trajectory of a vertexical family approaches the boundary if and only if either the trajectory approaches a vertex of the hypercube, or a trajectory in a lower-dimensional member of the family approaches the boundary. With this technology, we make progress on the global attractor conjecture, a central open problem concerning mass-action kinetics systems. Additionally, we phrase mass-action kinetics as a functor on reaction networks with variable rates.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chemical reaction systems with toric steady states

Mass-action chemical reaction systems are frequently used in Computational Biology. The corresponding polynomial dynamical systems are often large (consisting of tens or even hundreds of ordinary differential equations) and poorly parametrized (due to noisy measurement data and a small number of data points and repetitions). Therefore, it is often difficult to establish the existence of (positive) steady states or to determine whether more complicated phenomena such as multistationarity exist. If, however, the steady state ideal of the system is a binomial ideal, then we show that these questions can be answered easily. The focus of this work is on systems with this property, and we say that such systems have toric steady states. Our main result gives sufficient conditions for a chemical reaction system to have toric steady states. Furthermore, we analyze the capacity of such a system to exhibit positive steady states and multistationarity. Examples of systems with toric steady states include weakly-reversible zero-deficiency chemical reaction systems. An important application of our work concerns the networks that describe the multisite phosphorylation of a protein by a kinase/phosphatase pair in a sequential and distributive mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Siphons in chemical reaction networks

Siphons in a chemical reaction system are subsets of the species that have the potential of being absent in a steady state. We present a characterization of minimal siphons in terms of primary decomposition of binomial ideals, we explore the underlying geometry, and we demonstrate the effective computation of siphons using computer algebra software. This leads to a new method for determining whether given initial concentrations allow for various boundary steady states.

0 citations0 reviewsNo code linkedOpen access