A structural classification of candidate oscillatory and multistationary biochemical systems

Fibromyalgia pathogenesis explained by a neuroendocrine multistable model

Fibromyalgia (FM) is a central disorder characterized by chronic pain, fatigue, insomnia, depression, and other minor symptoms. Knowledge about pathogenesis is lacking, diagnosis difficult, clinical approach puzzling, and patient management disappointing. We conducted a theoretical study based on literature data and computational analysis, aimed at developing a comprehensive model of FM pathogenesis and addressing suitable therapeutic targets. We started from the evidence that FM must involve a dysregulation of central pain processing, is female prevalent, suggesting a role for the hypothalamus-pituitary-gonadal (HPG) axis, and is stress-related, suggesting a role for the HP-adrenocortical (HPA) axis. Central pathogenesis was supposed to involve a pain processing loop system including the thalamic ventroposterolateral nucleus (VPL), the primary somatosensory cortex (SSC), and the thalamic reticular nucleus (TRN). For decreasing GABAergic and/or increasing glutamatergic transmission, the loop system crosses a bifurcation point, switching from monostable to bistable, and converging on a high-firing-rate steady state supposed to be the pathogenic condition. Thereafter, we showed that GABAergic transmission is positively correlated with gonadal-hormone-derived neurosteroids, notably allopregnanolone, whereas glutamatergic transmission is positively correlated with stress-induced glucocorticoids, notably cortisol. Finally, we built a dynamic model describing a multistable, double-inhibitory loop between HPG and HPA axes. This system has a high-HPA/low-HPG steady state, allegedly reached in females under combined premenstrual/postpartum brain allopregnanolone withdrawal and stress condition, driving the thalamocortical loop to the high-firing-rate steady state, and explaining the connection between endocrine and neural mechanisms in FM pathogenesis. Our model accounts for FM female prevalence and stress correlation, suggesting the use of neurosteroid drugs as a possible solution to currently unsolved problems in the clinical treatment of the disease.

Chemical reaction motifs driving non-equilibrium behaviours in phase separating materials

Chemical reactions that couple to systems that phase separate have been implicated in diverse contexts from biology to materials science. However, how a particular set of chemical reactions (chemical reaction network, CRN) would affect the behaviours of a phase separating system is difficult to fully predict theoretically. In this paper, we analyse a mean field theory coupling CRNs to a combined system of phase separating and non-phase separating materials and analyse how the properties of the CRNs affect different classes of non-equilibrium behaviour: microphase separation or temporally oscillating patterns. We examine the problem of achieving microphase separated condensates by statistical analysis of the Jacobians, of which the most important motifs are negative feedback of the phase separating component and combined inhibition/activation by the non-phase separating components. We then identify CRN motifs that are likely to yield microphase by examining randomly generated networks and parameters. Molecular sequestration of the phase separating motif is shown to be the most robust towards yielding microphase separation. Subsequently, we find that dynamics of the phase separating species is promoted most easily by inducing oscillations in the diffusive components coupled to the phase separating species. Our results provide guidance towards the design of CRNs that manage the formation, dissolution and organization of compartments.

Bayesian Optimization for Design of Multiscale Biological Circuits

Recent advances in synthetic biology have enabled the construction of molecular circuits that operate across multiple scales of cellular organization, such as gene regulation, signaling pathways, and cellular metabolism. Computational optimization can effectively aid the design process, but current methods are generally unsuited for systems with multiple temporal or concentration scales, as these are slow to simulate due to their numerical stiffness. Here, we present a machine learning method for the efficient optimization of biological circuits across scales. The method relies on Bayesian optimization, a technique commonly used to fine-tune deep neural networks, to learn the shape of a performance landscape and iteratively navigate the design space toward an optimal circuit. This strategy allows the joint optimization of both circuit architecture and parameters, and provides a feasible approach to solve a highly nonconvex optimization problem in a mixed-integer input space. We illustrate the applicability of the method on several gene circuits for controlling biosynthetic pathways with strong nonlinearities, multiple interacting scales, and using various performance objectives. The method efficiently handles large multiscale problems and enables parametric sweeps to assess circuit robustness to perturbations, serving as an efficient in silico screening method prior to experimental implementation.

Vertex results for the robust analysis of uncertain biochemical systems

We consider the problem of assessing the sensitivity of uncertain biochemical systems in the presence of input perturbations (either constant or periodic) around a stable steady state. In particular, we propose approaches for the robust sensitivity analysis of systems with uncertain parameters assumed to take values in a hyper-rectangle. We highlight vertex results, which allow us to check whether a property is satisfied for all parameter choices in the hyper-rectangle by simply checking whether it is satisfied for all parameter choices at the vertices of the hyper-rectangle. We show that, for a vast class of systems, including (bio)chemical reaction networks with mass-action kinetics, the system Jacobian has a totally multiaffine structure (namely, all minors of the Jacobian matrix are multiaffine functions of the uncertain parameters), which can be exploited to obtain several vertex results. We consider different problems: robust non-singularity; robust stability of the steady-state; robust steady-state sensitivity analysis, in the case of constant perturbations; robust frequency-response sensitivity analysis, in the presence of periodic perturbations; and robust adaptation analysis. The developed theory is then applied to gain insight into some examples of uncertain biochemical systems, including the incoherent feed-forward loop, the coherent feed-forward loop, the Brusselator oscillator and the Goldbeter oscillator.

A general hypothesis of multistable systems in pathophysiology

Despite intensive investigations numerous diseases remain etiologically puzzling and recalcitrant to treatments. A hypothesis is proposed here assuming that these difficulties are due to an unsuitable approach to the mechanisms of life, which is subjugated by an apparent complexity and fails to grasp the uniformity that lays behind. The stability of metabolism, despite the enormous complex of chemical reactions, suggests that reciprocal control is a prerequisite of life. Negative feedback loops have been known for a long time to maintain homeostasis, while more recently, different life processes involved in transitions or changes have been modeled by positive loops giving rise to bistable switches, also including various diseases. The present hypothesis makes a generalization, by assuming that any functional element of a biological system is involved in a positive or a negative feedback loop. Consequently, the hypothesis holds that the starting mechanism of any disease that affects a healthy human can be conceptually reduced to a bistable or multistationary loop system, thus providing a unifying model leading to the discovery of critical therapeutic targets.

A general theory of multistable systems in pathophysiology

Despite intensive investigations numerous diseases remain etiologically puzzling and recalcitrant to treatments. A theory is proposed here assuming that these difficulties are due to an unsuitable approach to the mechanisms of life, which is subjugated by an apparent complexity and fails to grasp the uniformity that lays behind. The stability of metabolism, despite the enormous complex of chemical reactions, suggests that reciprocal control is a prerequisite of life. Negative feedback loops have been known for a long time to maintain homeostasis, while more recently, different life processes involved in transitions or changes have been modeled by positive loops giving rise to bistable switches, also including various diseases. The present theory makes a generalization, by assuming that any functional element of a biological system is involved in a positive or a negative feedback loop. Consequently, the theory holds that the starting mechanism of any disease that affects a healthy human can be conceptually reduced to a bistable or multistationary loop system, thus providing a unifying model leading to the discovery of critical therapeutic targets.

A general hypothesis of multistable systems in pathophysiology

Despite intensive investigations numerous diseases remain etiologically puzzling and recalcitrant to treatments. A hypothesis is proposed here assuming that these difficulties are due to an unsuitable approach to the mechanisms of life, which is subjugated by an apparent complexity and fails to grasp the uniformity that lays behind. The stability of metabolism, despite the enormous complex of chemical reactions, suggests that reciprocal control is a prerequisite of life. Negative feedback loops have been known for a long time to maintain homeostasis, while more recently, different life processes involved in transitions or changes have been modeled by positive loops giving rise to bistable switches, also including various diseases. The present hypothesis makes a generalization, by assuming that any functional element of a biological system is involved in a positive or a negative feedback loop. Consequently, the hypothesis holds that the starting mechanism of any disease that affects a healthy human can be conceptually reduced to a bistable or multistationary loop system, thus providing a unifying model leading to the discovery of critical therapeutic targets.

Thalamocortical bistable switch as a theoretical model of fibromyalgia pathogenesis inferred from a literature survey

Fibromyalgia (FM) is an unsolved central pain processing disturbance. We aim to provide a unifying model for FM pathogenesis based on a loop network involving thalamocortical regions, i.e., the ventroposterior lateral thalamus (VPL), the somatosensory cortex (SC), and the thalamic reticular nucleus (TRN). The dynamics of the loop have been described by three differential equations having neuron mean firing rates as variables and containing Hill functions to model mutual interactions among the loop elements. A computational analysis conducted with MATLAB has shown a transition from monostability to bistability of the loop behavior for a weakening of GABAergic transmission between TRN and VPL. This involves the appearance of a high-firing-rate steady state, which becomes dominant and is assumed to represent pathogenic pain processing giving rise to chronic pain. Our model is consistent with a bulk of literature evidence, such as neuroimaging and pharmacological data collected on FM patients, and with correlations between FM and immunoendocrine conditions, such as stress, perimenopause, chronic inflammation, obesity, and chronic dizziness. The model suggests that critical targets for FM treatment are to be found among immunoendocrine pathways leading to GABA/glutamate imbalance having an impact on the thalamocortical system.

A New Bistable Switch Model of Alzheimer’s Disease Pathogenesis

We propose a model to explain the pathogenesis of Alzheimer’s disease (AD) based on the theory that any disease affecting a healthy organism originates from a bistable feedback loop that shifts the system from a physiological to a pathological condition. We focused on the known double inhibitory loop involving the cellular prion protein (PrPC) and the enzyme BACE1 that produces amyloid-beta (Aβ) peptides. BACE1 is inhibited by PrPC, but its inhibitory activity is lost when PrPC binds to Aβ oligomers (Aβo). Excessive Aβo formation would switch the loop to a pathogenic condition involving the Aβo-PrPC-mGluR5 complex, Fyn kinase activation, tau, and NMDAR phosphorylation, ultimately leading to neurodegeneration. Based on the emerging role of cyclic nucleotides in Aβ production, and thereby in synaptic plasticity and cognitive processes, cAMP and cGMP can be considered as modulatory factors capable of inducing the transition from a physiological steady state to a pathogenic one. This would imply that critical pharmacological targets for AD treatment lie within pathways that lead to an imbalance of cyclic nucleotides in neurons. If this hypothesis is confirmed, it will provide precise indications for the development of preventive or therapeutic treatments for the disease.

Context-Dependent Stability and Robustness of Genetic Toggle Switches with Leaky Promoters

Multistable switches are ubiquitous building blocks in both systems and synthetic biology. Given their central role, it is thus imperative to understand how their fundamental properties depend not only on the tunable biophysical properties of the switches themselves, but also on their genetic context. To this end, we reveal in this article how these factors shape the essential characteristics of toggle switches implemented using leaky promoters such as their stability and robustness to noise, both at single-cell and population levels. In particular, our results expose the roles that competition for scarce transcriptional and translational resources, promoter leakiness, and cell-to-cell heterogeneity collectively play. For instance, the interplay between protein expression from leaky promoters and the associated cost of relying on shared cellular resources can give rise to tristable dynamics even in the absence of positive feedback. Similarly, we demonstrate that while promoter leakiness always acts against multistability, resource competition can be leveraged to counteract this undesirable phenomenon. Underpinned by a mechanistic model, our results thus enable the context-aware rational design of multistable genetic switches that are directly translatable to experimental considerations, and can be further leveraged during the synthesis of large-scale genetic systems using computer-aided biodesign automation platforms.

A multistationary loop model of ALS unveils critical molecular interactions involving mitochondria and glucose metabolism

Amyotrophic lateral sclerosis (ALS) is a poor-prognosis disease with puzzling pathogenesis and inconclusive treatments. We develop a mathematical model of ALS based on a system of interactive feedback loops, focusing on the mutant SOD1G93A mouse. Misfolded mutant SOD1 aggregates in motor neuron (MN) mitochondria and triggers a first loop characterized by oxidative phosphorylation impairment, AMP kinase over-activation, 6-phosphofructo-2-kinase (PFK3) rise, glucose metabolism shift from pentose phosphate pathway (PPP) to glycolysis, cell redox unbalance, and further worsening of mitochondrial dysfunction. Oxidative stress then triggers a second loop, involving the excitotoxic glutamatergic cascade, with cytosolic Ca2+ overload, increase of PFK3 expression, and further metabolic shift from PPP to glycolysis. Finally, cytosolic Ca2+ rise is also detrimental to mitochondria and oxidative phosphorylation, thus closing a third loop. These three loops are overlapped and positive (including an even number of inhibitory steps), hence they form a candidate multistationary (bistable) system. To describe the system dynamics, we model the interactions among the functional agents with differential equations. The system turns out to admit two stable equilibria: the healthy state, with high oxidative phosphorylation and preferential PPP, and the pathological state, with AMP kinase activation, PFK3 over expression, oxidative stress, excitotoxicity and MN degeneration. We demonstrate that the loop system is monotone: all functional agents consistently act toward the healthy or pathological condition, depending on low or high mutant SOD1 input. We also highlight that molecular interactions involving PFK3 are crucial, as their deletion disrupts the system's bistability leading to a single healthy equilibrium point. Hence, our mathematical model unveils that promising ALS management strategies should be targeted to mechanisms that keep low PFK3 expression and activity within MNs.

Mal de Debarquement Syndrome: A Matter of Loops?

Introduction: Mal de Debarquement Syndrome (MdDS) is a poorly understood neurological disorder affecting mostly perimenopausal women. MdDS has been hypothesized to be a maladaptation of the vestibulo-ocular reflex, a neuroplasticity disorder, and a consequence of neurochemical imbalances and hormonal changes. Our hypothesis considers elements from these theories, but presents a novel approach based on the analysis of functional loops, according to Systems and Control Theory.

Hypothesis: MdDS is characterized by a persistent sensation of self-motion, usually occurring after sea travels. We assume the existence of a neuronal mechanism acting as an oscillator, i.e., an adaptive internal model, that may be able to cancel a sinusoidal disturbance of posture experienced aboard, due to wave motion. Thereafter, we identify this mechanism as a multi-loop neural network that spans between vestibular nuclei and the flocculonodular lobe of the cerebellum. We demonstrate that this loop system has a tendency to oscillate, which increases with increasing strength of neuronal connections. Therefore, we hypothesize that synaptic plasticity, specifically long-term potentiation, may play a role in making these oscillations poorly damped. Finally, we assume that the neuromodulator Calcitonin Gene-Related Peptide, which is modulated in perimenopausal women, exacerbates this process thus rendering the transition irreversible and consequently leading to MdDS.

Conclusion and Validation: The concept of an oscillator that becomes noxiously permanent can be used as a model for MdDS, given a high correlation between patients with MdDS and sea travels involving undulating passive motion, and an alleviation of symptoms when patients are re-exposed to similar passive motion. The mechanism could be further investigated utilizing posturography tests to evaluate if subjective perception of motion matches with objective postural instability. Neurochemical imbalances that would render individuals more susceptible to developing MdDS could be investigated through hormonal profile screening. Alterations in the connections between vestibular nuclei and cerebellum, notably GABAergic fibers, could be explored by neuroimaging techniques as well as transcranial magnetic stimulation. If our hypothesis were tested and verified, optimal targets for MdDS treatment could be found within both the neural networks and biochemical factors that are deemed to play a fundamental role in loop functioning and synaptic plasticity.

Loop analysis of blood pressure/volume homeostasis

We performed a mathematical analysis of the dynamic control loops regulating the vasomotor tone of vascular smooth muscle, blood volume, and mean arterial pressure, which involve the arginine vasopressin (AVP) system, the atrial natriuretic peptide system (ANP), and the renin-angiotensin-aldosterone system (RAAS). Our loop analysis of the AVP-ANP-RAAS system revealed the concurrent presence of two different regulatory mechanisms, which perform the same qualitative function: one affects blood pressure by regulating vasoconstriction, the other by regulating blood volume. Both the systems are candidate oscillators consisting of the negative-feedback loop of a monotone system: they admit a single equilibrium that can either be stable or give rise to oscillatory instability. Also a subsystem, which includes ANP and AVP stimulation of vascular smooth muscle cells, turns out to be a candidate oscillator composed of a monotone system with multiple negative feedback loops, and we show that its oscillatory potential is higher when the delays along all feedback loops are comparable. Our results give insight into the physiological mechanisms ruling long-term homeostasis of blood hydraulic parameters, which operate based on dynamical loops of interactions.

The efficiency and resilience of our body are guaranteed by the presence of myriads of dynamic control loops that regulate fundamental vital functions. In this work, we studied the regulatory mechanisms that govern the interplay of vasoconstriction/vasodilation, blood volume and mean arterial pressure. We analysed the loops in the system and showed the presence of two coexisting mechanisms for blood pressure regulation, which perform the same qualitative function, conferring robustness to the system: one mechanism tunes vasoconstriction, the other blood volume. We showed that both systems are candidate oscillators: either they are stable or they oscillate regularly around their unique equilibrium. We analysed a subsystem that describes the stimulation of vascular smooth muscle cells due to the hormones arginine vasopressin (AVP) and atrial natriuretic peptide (ANP): also this system is a candidate oscillator ruled by multiple negative-feedback loops, and its potential for oscillations is higher when all the loops have similar delay. Our results cast light on the fundamental physiological phenomena that preserve the stable functioning of blood pressure and volume. This could have even wider relevance if other homeostasis and endocrine systems displayed similar features, with repercussions on the management of adverse homeostasis shifts like hypertension.

Homogeneous Time Constants Promote Oscillations in Negative Feedback Loops

Biological oscillators are present in nearly all self-regulating systems, from individual cells to entire organisms. In any oscillator structure, a negative feedback loop is necessary, but not sufficient to guarantee the emergence of periodic behaviors. The likelihood of oscillations can be improved by careful tuning of the system time constants and by increasing the loop gain, yet it is unclear whether there is any general relationship between optimal time constants and loop gain. This issue is particularly relevant in genetic oscillators resulting from a chain of different subsequent biochemical events, each with distinct (and uncertain) kinetics. Using two families of genetic oscillators as model examples, we show that the loop gain required for oscillations is minimum when all elements in the loop have the same time constant. On the contrary, we show that homeostasis is ensured if a single element is considerably slower than the others.

Ceramide-transfer protein-mediated ceramide transfer is a structurally tunable flow-inducing mechanism with structural feed-forward loops

This paper considers two models of ceramide-transfer protein (CERT)-mediated ceramide transfer at the trans-Golgi network proposed in the literature, short distance shuttle and neck swinging, and seeks structural (parameter-free) features of the two models, which rely exclusively on the peculiar interaction network and not on specific parameter values. In particular, it is shown that both models can be seen as flow-inducing systems, where the flows between pairs of species are tuned by the concentrations of other species, and suitable external inputs can structurally regulate ceramide transfer. In the short distance shuttle model, the amount of transferred ceramide is structurally tuned by active protein kinase D (PKD), both directly and indirectly, in a coherent feed-forward loop motif. In the neck-swinging model, the amount of transferred ceramide is structurally tuned by active PI4KIIIβ, while active PKD has an ambivalent effect, due to the presence of an incoherent feed-forward loop motif that directly inhibits ceramide transfer and indirectly promotes it; the structural role of active PKD is to favour CERT mobility in the cytosol. It is also shown that the influences among key variables often have structurally determined steady-state signs, which can help falsify the models against experimental traces.

A Robust Molecular Network Motif for Period-Doubling Devices

Life is sustained by a variety of cyclic processes such as cell division, muscle contraction, and neuron firing. The periodic signals powering these processes often direct a variety of other downstream systems, which operate at different time scales and must have the capacity to divide or multiply the period of the master clock. Period modulation is also an important challenge in synthetic molecular systems, where slow and fast components may have to be coordinated simultaneously by a single oscillator whose frequency is often difficult to tune. Circuits that can multiply the period of a clock signal (frequency dividers), such as binary counters and flip-flops, are commonly encountered in electronic systems, but design principles to obtain similar devices in biological systems are still unclear. We take inspiration from the architecture of electronic flip-flops, and we propose to build biomolecular period-doubling networks by combining a bistable switch with negative feedback modules that preprocess the circuit inputs. We identify a network motif and we show it can be "realized" using different biomolecular components; two of the realizations we propose rely on transcriptional gene networks and one on nucleic acid strand displacement systems. We examine the capacity of each realization to perform period-doubling by studying how bistability of the motif is affected by the presence of the input; for this purpose, we employ mathematical tools from algebraic geometry that provide us with valuable insights on the input/output behavior as a function of the realization parameters. We show that transcriptional network realizations operate correctly also in a stochastic regime when processing oscillations from the repressilator, a canonical synthetic in vivo oscillator. Finally, we compare the performance of different realizations in a range of realistic parameters via numerical sensitivity analysis of the period-doubling region, computed with respect to the input period and amplitude. Our mathematical and computational analysis suggests that the motif we propose is generally robust with respect to specific implementation details: functionally equivalent circuits can be built as long as the species-interaction topology is respected. This indicates that experimental construction of the circuit is possible with a variety of components within the rapidly expanding libraries available in synthetic biology.

Stability analysis of an artificial biomolecular oscillator with non-cooperative regulatory interactions

Oscillators are essential to fuel autonomous behaviours in molecular systems. Artificial oscillators built with programmable biological molecules such as DNA and RNA are generally easy to build and tune, and can serve as timers for biological computation and regulation. We describe a new artificial nucleic acid biochemical reaction network, and we demonstrate its capacity to exhibit oscillatory solutions. This network can be built in vitro using nucleic acids and three bacteriophage enzymes, and has the potential to be implemented in cells. Numerical simulations suggest that oscillations occur in a realistic range of reaction rates and concentrations.

Qualitative and quantitative responses to press perturbations in ecological networks

Predicting the sign of press perturbation responses in ecological networks is challenging, due to the poor knowledge of the strength of the direct interactions among the species, and to the entangled coexistence of direct and indirect effects. We show in this paper that, for a class of networks that includes mutualistic and monotone networks, the sign of press perturbation responses can be qualitatively determined based only on the sign pattern of the community matrix, without any knowledge of parameter values. For other classes of networks, we show that a semi-qualitative approach yields sufficient conditions for community matrices with a given sign pattern to exhibit mutualistic responses to press perturbations; quantitative conditions can be provided as well for community matrices that are eventually nonnegative. We also present a computational test that can be applied to any class of networks so as to check whether the sign of the responses to press perturbations is constant in spite of parameter variations.

Type-dependent stochastic Ising model describing the dynamics of a non-symmetric feedback module

We study an alternative approach to model the dynamical behaviors of biological feedback loop, that is, a type-dependent spin system, this class of stochastic models was introduced by Fernández et. al [13], and are useful since take account to inherent variability of gene expression. We analyze a non-symmetric feedback module being an extension for the repressilator, the first synthetic biological oscillator, invented by Elowitz and Leibler [7]. We consider a mean-field dynamics for a type-dependent Ising model, and then study the empirical-magnetization vector representing concentration of molecules. We apply a convergence result from stochastic jump processes to deterministic trajectories and present a bifurcation analysis for the associated dynamical system. We show that non-symmetric module under study can exhibit very rich behaviours, including the empirical oscillations described by repressilator.

Molecular Titration Promotes Oscillations and Bistability in Minimal Network Models with Monomeric Regulators

Molecular titration is emerging as an important biochemical interaction mechanism within synthetic devices built with nucleic acids and the CRISPR/Cas system. We show that molecular titration in the context of feedback circuits is a suitable mechanism to enhance the emergence of oscillations and bistable behaviors. We consider biomolecular modules that can be inhibited or activated by input monomeric regulators; the regulators compete with constitutive titrating species to determine the activity of their target. By tuning the titration rate and the concentration of titrating species, it is possible to modulate the delay and convergence speed of the transient response, and the steepness and dead zone of the stationary response of the modules. These phenomena favor the occurrence of oscillations when modules are interconnected to create a negative feedback loop; bistability is favored in a positive feedback interconnection. Numerical simulations are supported by mathematical analysis showing that the capacity of the closed loop systems to exhibit oscillations or bistability is structural.

An analytical approach to bistable biological circuit discrimination using real algebraic geometry

Biomolecular circuits with two distinct and stable steady states have been identified as essential components in a wide range of biological networks, with a variety of mechanisms and topologies giving rise to their important bistable property. Understanding the differences between circuit implementations is an important question, particularly for the synthetic biologist faced with determining which bistable circuit design out of many is best for their specific application. In this work we explore the applicability of Sturm's theorem--a tool from nineteenth-century real algebraic geometry--to comparing 'functionally equivalent' bistable circuits without the need for numerical simulation. We first consider two genetic toggle variants and two different positive feedback circuits, and show how specific topological properties present in each type of circuit can serve to increase the size of the regions of parameter space in which they function as switches. We then demonstrate that a single competitive monomeric activator added to a purely monomeric (and otherwise monostable) mutual repressor circuit is sufficient for bistability. Finally, we compare our approach with the Routh-Hurwitz method and derive consistent, yet more powerful, parametric conditions. The predictive power and ease of use of Sturm's theorem demonstrated in this work suggest that algebraic geometric techniques may be underused in biomolecular circuit analysis.

Computing the structural influence matrix for biological systems

We consider the problem of identifying structural influences of external inputs on steady-state outputs in a biological network model. We speak of a structural influence if, upon a perturbation due to a constant input, the ensuing variation of the steady-state output value has the same sign as the input (positive influence), the opposite sign (negative influence), or is zero (perfect adaptation), for any feasible choice of the model parameters. All these signs and zeros can constitute a structural influence matrix, whose (i, j) entry indicates the sign of steady-state influence of the jth system variable on the ith variable (the output caused by an external persistent input applied to the jth variable). Each entry is structurally determinate if the sign does not depend on the choice of the parameters, but is indeterminate otherwise. In principle, determining the influence matrix requires exhaustive testing of the system steady-state behaviour in the widest range of parameter values. Here we show that, in a broad class of biological networks, the influence matrix can be evaluated with an algorithm that tests the system steady-state behaviour only at a finite number of points. This algorithm also allows us to assess the structural effect of any perturbation, such as variations of relevant parameters. Our method is applied to nontrivial models of biochemical reaction networks and population dynamics drawn from the literature, providing a parameter-free insight into the system dynamics.

Using sign patterns to detect the possibility of periodicity in biological systems

Many models in the physical and life sciences, formulated as dynamical systems, exhibit a positive steady state, with its local qualitative behavior determined by the eigenvalues of its Jacobian matrix. Our interest lies in detecting if this steady state is linearly stable or if the system has periodic solutions arising from a Hopf bifurcation. We address this by considering the sign pattern of the Jacobian matrix and its set of allowed refined inertias. The refined inertia of a matrix, which is an extension of the classical matrix inertia, is a property of its eigenvalues. A Hopf bifurcation, leading to periodic solutions, may be possible if the sign pattern of the Jacobian matrix allows a specific set of refined inertias. For most systems, we also need to consider magnitude restrictions on the entries of the Jacobian matrix that are a consequence of the particular biological model. The usefulness of sign pattern analysis to detect linear stability or the possibility of periodicity is illustrated with several biological examples, including metabolic-genetic circuits, biochemical reaction networks, predator-prey and competition systems.