Anne Shiu

Toric dynamical systems

Toric dynamical systems are known as complex balancing mass action systems in the mathematical chemistry literature, where many of their remarkable properties have been established. They include as special cases all deficiency zero systems and all detailed balancing systems. One feature is that the steady state locus of a toric dynamical system is a toric variety, which has a unique point within each invariant polyhedron. We develop the basic theory of toric dynamical systems in the context of computational algebraic geometry and show that the associated moduli space is also a toric variety. It is conjectured that the complex balancing state is a global attractor. We prove this for detailed balancing systems whose invariant polyhedron is two-dimensional and bounded.

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 parameterized (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.

Convex Rank Tests and Semigraphoids

Convex rank tests are partitions of the symmetric group which have desirable geometric properties. The statistical tests defined by such partitions involve counting all permutations in the equivalence classes. Each class consists of the linear extensions of a partially ordered set specified by data. Our methods refine existing rank tests of nonparametric statistics, such as the sign test and the runs test, and are useful for exploratory analysis of ordinal data. We establish a bijection between convex rank tests and probabilistic conditional independence structures known as semigraphoids. The subclass of submodular rank tests is derived from faces of the cone of submodular functions or from Minkowski summands of the permutohedron. We enumerate all small instances of such rank tests. Of particular interest are graphical tests, which correspond to both graphical models and to graph associahedra.

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.

Comparison of pattern detection methods in microarray time series of the segmentation clock.

While genome-wide gene expression data are generated at an increasing rate, the repertoire of approaches for pattern discovery in these data is still limited. Identifying subtle patterns of interest in large amounts of data (tens of thousands of profiles) associated with a certain level of noise remains a challenge. A microarray time series was recently generated to study the transcriptional program of the mouse segmentation clock, a biological oscillator associated with the periodic formation of the segments of the body axis. A method related to Fourier analysis, the Lomb-Scargle periodogram, was used to detect periodic profiles in the dataset, leading to the identification of a novel set of cyclic genes associated with the segmentation clock. Here, we applied to the same microarray time series dataset four distinct mathematical methods to identify significant patterns in gene expression profiles. These methods are called: Phase consistency, Address reduction, Cyclohedron test and Stable persistence, and are based on different conceptual frameworks that are either hypothesis- or data-driven. Some of the methods, unlike Fourier transforms, are not dependent on the assumption of periodicity of the pattern of interest. Remarkably, these methods identified blindly the expression profiles of known cyclic genes as the most significant patterns in the dataset. Many candidate genes predicted by more than one approach appeared to be true positive cyclic genes and will be of particular interest for future research. In addition, these methods predicted novel candidate cyclic genes that were consistent with previous biological knowledge and experimental validation in mouse embryos. Our results demonstrate the utility of these novel pattern detection strategies, notably for detection of periodic profiles, and suggest that combining several distinct mathematical approaches to analyze microarray datasets is a valuable strategy for identifying genes that exhibit novel, interesting transcriptional patterns.

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.

THE DYNAMICS OF WEAKLY REVERSIBLE POPULATION PROCESSES NEAR FACETS

This paper concerns the dynamical behavior of weakly reversible, deterministically modeled population processes near the facets (codimension-one faces) of their invariant manifolds and proves that the facets of such systems are “repelling.” It has been conjectured that any population process whose network graph is weakly reversible (has strongly connected components) is persistent. We prove this conjecture to be true for the subclass of weakly reversible systems for which only facets of the invariant manifold are associated with semilocking sets, or siphons. An important application of this work pertains to chemical reaction systems that are complex-balancing. For these systems it is known that within the interior of each invariant manifold there is a unique equilibrium. The global attractor conjecture states that each of these equilibria is globally asymptotically stable relative to the interior of the invariant manifold in which it lies. Our results pertaining to weakly reversible systems imply that thi...

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 moleculari...

Three counter-examples on semi-graphoids

Semi-graphoids are combinatorial structures that arise in statistical learning theory. They are equivalent to convex rank tests and to polyhedral fans that coarsen the reflection arrangement of the symmetric group Sn. In this paper we resolve two problems on semi-graphoids posed in Studenýs book (2005), and we answer a related question of Postnikov, Reiner and Williams on generalized permutohedra. We also study the semigroup and the toric ideal associated with semi-graphoids.

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.