Time-scale separation--Michaelis and Menten's old idea, still bearing fruit

Regulatory properties of transcription factors with diverse mechanistic function

Transcription factors (TFs) regulate the process of transcription through the modulation of different kinetic steps. Although models can often describe the observed transcriptional output of a measured gene, predicting a TFs role on a given promoter requires an understanding of how the TF alters each step of the transcription process. In this work, we use a simple model of transcription to assess the role of promoter identity, and the degree to which TFs alter binding of RNAP (stabilization) and initiation of transcription (acceleration) on three primary characteristics: the range of steady-state regulation, cell-to-cell variability in expression, and the dynamic response time of a regulated gene. We find that steady state regulation and the response time of a gene behave uniquely for TFs that regulate incoherently, i.e that speed up one step but slow the other. We also find that incoherent TFs have dynamic implications, with one type of incoherent mode configuring the promoter to respond more slowly at intermediate TF concentrations. We also demonstrate that the noise of gene expression for these TFs is sensitive to promoter strength, with a distinct non-monotonic profile that is apparent under stronger promoters. Taken together, our work uncovers the coupling between promoters and TF regulatory modes with implications for understanding natural promoters and engineering synthetic gene circuits with desired expression properties.

The Hill function is the universal Hopfield barrier for sharpness of input-output responses

The Hill functions, [Formula: see text], have been widely used in biology for over a century but, with the exception of [Formula: see text], they have had no justification other than as a convenient fit to empirical data. Here, we show that they are the universal limit for the sharpness of any input-output response arising from a Markov process model at thermodynamic equilibrium. Models may represent arbitrary molecular complexity, with multiple ligands, internal states, conformations, coregulators, etc, under core assumptions that are detailed in the paper. The model output may be any linear combination of steady-state probabilities, with components other than the chosen input ligand held constant. This formulation generalizes most of the responses in the literature. We use a coarse-graining method in the graph-theoretic linear framework to show that two sharpness measures for input-output responses fall within an effectively bounded region of the positive quadrant, [Formula: see text], for any equilibrium model with [Formula: see text] input binding sites. [Formula: see text] exhibits a cusp which approaches, but never exceeds, the sharpness of [Formula: see text], but the region and the cusp can be exceeded when models are taken away from thermodynamic equilibrium. Such fundamental thermodynamic limits are called Hopfield barriers, and our results provide a biophysical justification for the Hill functions as the universal Hopfield barriers for sharpness. Our results also introduce an object, [Formula: see text], whose structure may be of mathematical interest, and suggest the importance of characterizing Hopfield barriers for other forms of cellular information processing.

The Hill function is the universal Hopfield barrier for sharpness of input-output responses

The Hill functions, ℋ h ( x ) = x h / 1 + x h , have been widely used in biology for over a century but, with the exception of ℋ 1 , they have had no justification other than as a convenient fit to empirical data. Here, we show that they are the universal limit for the sharpness of any input-output response arising from a Markov process model at thermodynamic equilibrium. Models may represent arbitrary molecular complexity, with multiple ligands, internal states, conformations, co-regulators, etc, under core assumptions that are detailed in the paper. The model output may be any linear combination of steady-state probabilities, with components other than the chosen input ligand held constant. This formulation generalises most of the responses in the literature. We use a coarse-graining method in the graph-theoretic linear framework to show that two sharpness measures for input-output responses fall within an effectively bounded region of the positive quadrant, Ω m ⊂ ℝ + 2 , for any equilibrium model with m input binding sites. Ω m exhibits a cusp which approaches, but never exceeds, the sharpness of ℋ m but the region and the cusp can be exceeded when models are taken away from thermodynamic equilibrium. Such fundamental thermodynamic limits are called Hopfield barriers and our results provide a biophysical justification for the Hill functions as the universal Hopfield barriers for sharpness. Our results also introduce an object, Ω m , whose structure may be of mathematical interest, and suggest the importance of characterising Hopfield barriers for other forms of cellular information processing.

Is it selfish to be filamentous in biofilms? Individual-based modeling links microbial growth strategies with morphology using the new and modular iDynoMiCS 2.0

Microbial communities are found in all habitable environments and often occur in assemblages with self-organized spatial structures developing over time. This complexity can only be understood, predicted, and managed by combining experiments with mathematical modeling. Individual-based models are particularly suited if individual heterogeneity, local interactions, and adaptive behavior are of interest. Here we present the completely overhauled software platform, the individual-based Dynamics of Microbial Communities Simulator, iDynoMiCS 2.0, which enables researchers to specify a range of different models without having to program. Key new features and improvements are: (1) Substantially enhanced ease of use (graphical user interface, editor for model specification, unit conversions, data analysis and visualization and more). (2) Increased performance and scalability enabling simulations of up to 10 million agents in 3D biofilms. (3) Kinetics can be specified with any arithmetic function. (4) Agent properties can be assembled from orthogonal modules for pick and mix flexibility. (5) Force-based mechanical interaction framework enabling attractive forces and non-spherical agent morphologies as an alternative to the shoving algorithm. The new iDynoMiCS 2.0 has undergone intensive testing, from unit tests to a suite of increasingly complex numerical tests and the standard Benchmark 3 based on nitrifying biofilms. A second test case was based on the "biofilms promote altruism" study previously implemented in BacSim because competition outcomes are highly sensitive to the developing spatial structures due to positive feedback between cooperative individuals. We extended this case study by adding morphology to find that (i) filamentous bacteria outcompete spherical bacteria regardless of growth strategy and (ii) non-cooperating filaments outcompete cooperating filaments because filaments can escape the stronger competition between themselves. In conclusion, the new substantially improved iDynoMiCS 2.0 joins a growing number of platforms for individual-based modeling of microbial communities with specific advantages and disadvantages that we discuss, giving users a wider choice.

Generalized Michaelis-Menten rate law with time-varying molecular concentrations

The Michaelis-Menten (MM) rate law has been the dominant paradigm of modeling biochemical rate processes for over a century with applications in biochemistry, biophysics, cell biology, systems biology, and chemical engineering. The MM rate law and its remedied form stand on the assumption that the concentration of the complex of interacting molecules, at each moment, approaches an equilibrium (quasi-steady state) much faster than the molecular concentrations change. Yet, this assumption is not always justified. Here, we relax this quasi-steady state requirement and propose the generalized MM rate law for the interactions of molecules with active concentration changes over time. Our approach for time-varying molecular concentrations, termed the effective time-delay scheme (ETS), is based on rigorously estimated time-delay effects in molecular complex formation. With particularly marked improvements in protein-protein and protein-DNA interaction modeling, the ETS provides an analytical framework to interpret and predict rich transient or rhythmic dynamics (such as autogenously-regulated cellular adaptation and circadian protein turnover), which goes beyond the quasi-steady state assumption.

Design principles of 3D epigenetic memory systems

Cells remember their identities, in part, by using epigenetic marks-chemical modifications placed along the genome. How can mark patterns remain stable over cell generations despite their constant erosion by replication and other processes? We developed a theoretical model that reveals that three-dimensional (3D) genome organization can stabilize epigenetic memory as long as (i) there is a large density difference between chromatin compartments, (ii) modifying "reader-writer" enzymes spread marks in three dimensions, and (iii) the enzymes are limited in abundance relative to their histone substrates. Analogous to an associative memory that encodes memory in neuronal connectivity, mark patterns are encoded in a 3D network of chromosomal contacts. Our model provides a unified account of diverse observations and reveals a key role of 3D genome organization in epigenetic memory.

The linear framework II: using graph theory to analyse the transient regime of Markov processes

The linear framework uses finite, directed graphs with labelled edges to model biomolecular systems. Graph vertices represent chemical species or molecular states, edges represent reactions or transitions and edge labels represent rates that also describe how the system is interacting with its environment. The present paper is a sequel to a recent review of the framework that focussed on how graph-theoretic methods give insight into steady states as rational algebraic functions of the edge labels. Here, we focus on the transient regime for systems that correspond to continuous-time Markov processes. In this case, the graph specifies the infinitesimal generator of the process. We show how the moments of the first-passage time distribution, and related quantities, such as splitting probabilities and conditional first-passage times, can also be expressed as rational algebraic functions of the labels. This capability is timely, as new experimental methods are finally giving access to the transient dynamic regime and revealing the computations and information processing that occur before a steady state is reached. We illustrate the concepts, methods and formulas through examples and show how the results may be used to illuminate previous findings in the literature.

Algorithmic criteria for the validity of quasi-steady state and partial equilibrium models: the Michaelis-Menten reaction mechanism

We present "on the fly" algorithmic criteria for the accuracy and stability (non-stiffness) of reduced models constructed with the quasi-steady state and partial equilibrium approximations. The criteria comprise those introduced in Goussis (Combust Theor Model 16:869-926, 2012) that addressed the case where each fast time scale is due to one reaction and a new one that addresses the case where a fast time scale is due to more than one reactions. The development of these criteria is based on the ability to approximate accurately the fast and slow subspaces of the tangent space. Their validity is assessed on the basis of the Michaelis-Menten reaction mechanism, for which extensive literature is available regarding the validity of the existing various reduced models. The criteria predict correctly the regions in both the parameter and phase spaces where each of these models is valid. The findings are supported by numerical computations at indicative points in the parameter space. Due to their algorithmic character, these criteria can be readily employed for the reduction of large and complex mathematical models.

Thermodynamic bounds on ultrasensitivity in covalent switching

Switch-like motifs are among the basic building blocks of biochemical networks. A common motif that can serve as an ultrasensitive switch consists of two enzymes acting antagonistically on a substrate, one making and the other removing a covalent modification. To work as a switch, such covalent modification cycles must be held out of thermodynamic equilibrium by continuous expenditure of energy. Here, we exploit the linear framework for timescale separation to establish tight bounds on the performance of any covalent-modification switch in terms of the chemical potential difference driving the cycle. The bounds apply to arbitrary enzyme mechanisms, not just Michaelis-Menten, with arbitrary rate constants and thereby reflect fundamental physical constraints on covalent switching.

Protein-protein complexes can undermine ultrasensitivity-dependent biological adaptation

Robust perfect adaptation (RPA) is a ubiquitously observed signalling response across all scales of biological organization. A major class of network architectures that drive RPA in complex networks is the Opposer module-a feedback-regulated network into which specialized integral-computing 'opposer node(s)' are embedded. Although ultrasensitivity-generating chemical reactions have long been considered a possible mechanism for such adaptation-conferring opposer nodes, this hypothesis has relied on simplified Michaelian models, which neglect the presence of protein-protein complexes. Here we develop complex-complete models of interlinked covalent-modification cycles with embedded ultrasensitivity, explicitly capturing all molecular interactions and protein complexes. Strikingly, we demonstrate that the presence of protein-protein complexes thwarts the network's capacity for RPA in any 'free' active protein form, conferring RPA capacity instead on the concentration of a larger protein pool consisting of two distinct forms of a single protein. We further show that the presence of enzyme-substrate complexes, even at comparatively low concentrations, play a crucial and previously unrecognized role in controlling the RPA response-significantly reducing the range of network inputs for which RPA can obtain, and imposing greater parametric requirements on the RPA response. These surprising results raise fundamental new questions as to the biochemical requirements for adaptation-conferring Opposer modules within complex cellular networks.

SBML to bond graphs: From conversion to composition

The Systems Biology Markup Language (SBML) is a popular software-independent XML-based format for describing models of biological phenomena. The BioModels Database is the largest online repository of SBML models. Several tools and platforms are available to support the reuse and composition of SBML models. However, these tools do not explicitly assess whether models are physically plausible or thermodynamically consistent. This often leads to ill-posed models that are physically impossible, impeding the development of realistic complex models in biology. Here, we present a framework that can automatically convert SBML models into bond graphs, which imposes energy conservation laws on these models. The new bond graph models are easily mergeable, resulting in physically plausible coupled models. We illustrate this by automatically converting and coupling a model of pyruvate distribution to a model of the pentose phosphate pathway.

The linear framework: using graph theory to reveal the algebra and thermodynamics of biomolecular systems

The linear framework uses finite, directed graphs with labelled edges to model biomolecular systems. Graph vertices represent biochemical species or molecular states, edges represent reactions or transitions and labels represent rates. The graph yields a linear dynamics for molecular concentrations or state probabilities, with the graph Laplacian as the operator, and the labels encode the nonlinear interactions between system and environment. The labels can be specified by vertices of other graphs or by conservation laws or, when the environment consists of thermodynamic reservoirs, they may be constants. In the latter case, the graphs correspond to infinitesimal generators of Markov processes. The key advantage of the framework has been that steady states are determined as rational algebraic functions of the labels by the Matrix-Tree theorems of graph theory. When the system is at thermodynamic equilibrium, this prescription recovers equilibrium statistical mechanics but it continues to hold for non-equilibrium steady states. The framework goes beyond other graph-based approaches in treating the graph as a mathematical object, for which general theorems can be formulated that accommodate biomolecular complexity. It has been particularly effective at analysing enzyme-catalysed modification systems and input-output responses.

Artificial neural networks enable genome-scale simulations of intracellular signaling

Mammalian cells adapt their functional state in response to external signals in form of ligands that bind receptors on the cell-surface. Mechanistically, this involves signal-processing through a complex network of molecular interactions that govern transcription factor activity patterns. Computer simulations of the information flow through this network could help predict cellular responses in health and disease. Here we develop a recurrent neural network framework constrained by prior knowledge of the signaling network with ligand-concentrations as input and transcription factor-activity as output. Applied to synthetic data, it predicts unseen test-data (Pearson correlation r = 0.98) and the effects of gene knockouts (r = 0.8). We stimulate macrophages with 59 different ligands, with and without the addition of lipopolysaccharide, and collect transcriptomics data. The framework predicts this data under cross-validation (r = 0.8) and knockout simulations suggest a role for RIPK1 in modulating the lipopolysaccharide response. This work demonstrates the feasibility of genome-scale simulations of intracellular signaling.

Designing Biological Circuits: From Principles to Applications

Genetic circuit design is a well-studied problem in synthetic biology. Ever since the first genetic circuits─the repressilator and the toggle switch─were designed and implemented, many advances have been made in this area of research. The current review systematically organizes a number of key works in this domain by employing the versatile framework of generalized morphological analysis. Literature in the area has been mapped on the basis of (a) the design methodologies used, ranging from brute-force searches to control-theoretic approaches, (b) the modeling techniques employed, (c) various circuit functionalities implemented, (d) key design characteristics, and (e) the strategies used for the robust design of genetic circuits. We conclude our review with an outlook on multiple exciting areas for future research, based on the systematic assessment of key research gaps that have been readily unravelled by our analysis framework.

A numerical approach for detecting switch-like bistability in mass action chemical reaction networks with conservation laws

BACKGROUND

Theoretical analysis of signaling pathways can provide a substantial amount of insight into their function. One particular area of research considers signaling pathways capable of assuming two or more stable states given the same amount of signaling ligand. This phenomenon of bistability can give rise to switch-like behavior, a mechanism that governs cellular decision making. Investigation of whether or not a signaling pathway can confer bistability and switch-like behavior, without knowledge of specific kinetic rate constant values, is a mathematically challenging problem. Recently a technique based on optimization has been introduced, which is capable of finding example parameter values that confer switch-like behavior for a given pathway. Although this approach has made it possible to analyze moderately sized pathways, it is limited to reaction networks that presume a uniterminal structure. It is this limited structure we address by developing a general technique that applies to any mass action reaction network with conservation laws.

RESULTS

In this paper we developed a generalized method for detecting switch-like bistable behavior in any mass action reaction network with conservation laws. The method involves (1) construction of a constrained optimization problem using the determinant of the Jacobian of the underlying rate equations, (2) minimization of the objective function to search for conditions resulting in a zero eigenvalue, (3) computation of a confidence level that describes if the global minimum has been found and (4) evaluation of optimization values, using either numerical continuation or directly simulating the ODE system, to verify that a bistability region exists. The generalized method has been tested on three motifs known to be capable of bistability.

CONCLUSIONS

We have developed a variation of an optimization-based method for the discovery of bistability, which is not limited to uniterminal chemical reaction networks. Successful completion of the method provides an S-shaped bifurcation diagram, which indicates that the network acts as a bistable switch for the given optimization parameters.

Beyond the Michaelis-Menten: Bayesian Inference for Enzyme Kinetic Analysis

Although the Michaelis-Menten (MM) rate law has been widely used to estimate enzyme kinetic parameters, it works only under the condition of extremely low enzyme concentration. Furthermore, even when this condition is satisfied, parameter estimation is often imprecise due to the parameter identifiability issue. To overcome these limitations of the canonical approach to enzyme kinetics, we developed a Bayesian approach based on a modified form of the MM rate law, which is derived with the total quasi-steady state approximation. Here, we illustrate how to perform the Bayesian inference for the progress curve assay with our user-friendly computational R package. We also describe an optimal experimental design for the progress curve assay, with which enzyme kinetic parameters can be accurately and precisely estimated from minimal measurements of the progress curves.

When and why direct transmission models can be used for environmentally persistent pathogens

Variants of the susceptible-infected-removed (SIR) model of Kermack & McKendrick (1927) enjoy wide application in epidemiology, offering simple yet powerful inferential and predictive tools in the study of diverse infectious diseases across human, animal and plant populations. Direct transmission models (DTM) are a subset of these that treat the processes of disease transmission as comprising a series of discrete instantaneous events. Infections transmitted indirectly by persistent environmental pathogens, however, are examples where a DTM description might fail and are perhaps better described by models that comprise explicit environmental transmission routes, so-called environmental transmission models (ETM). In this paper we discuss the stochastic susceptible-exposed-infected-removed (SEIR) DTM and susceptible-exposed-infected-removed-pathogen (SEIR-P) ETM and we show that the former is the timescale separation limit of the latter, with ETM host-disease dynamics increasingly resembling those of a DTM when the pathogen’s characteristic timescale is shortened, relative to that of the host population. Using graphical posterior predictive checks (GPPC), we investigate the validity of the SEIR model when fitted to simulated SEIR-P host infection and removal times. Such analyses demonstrate how, in many cases, the SEIR model is robust to departure from direct transmission. Finally, we present a case study of white spot disease (WSD) in penaeid shrimp with rates of environmental transmission and pathogen decay (SEIR-P model parameters) estimated using published results of experiments. Using SEIR and SEIR-P simulations of a hypothetical WSD outbreak management scenario, we demonstrate how relative shortening of the pathogen timescale comes about in practice. With atttempts to remove diseased shrimp from the population every 24h, we see SEIR and SEIR-P model outputs closely conincide. However, when removals are 6-hourly, the two models’ mean outputs diverge, with distinct predictions of outbreak size and duration.

Mathematical models of the spread and progression of communicable disease in populations are important tools in efforts to prevent and control outbreaks. A common class of disease models assume that infection is transmitted directly from infectious to susceptible individuals when they are in close proximity—so called direct transmission models. These are used widely and have proven invaluable as simplified descriptions of a wide array of infectious diseases in diverse populations. However, many pathogens spread through indirect, environmental routes of transmission, for example via contact with contaminated water sources in the case of cholera, or inhalation of infectious airborne droplets for respiratory infections, such as Covid-19.

We show that direct transmission models work well for such pathogens with short environmental lifetimes and where hosts shed pathogens into the environment at high rates. This means that we do not require information about environmental pathogen levels to understand the behaviour of outbreaks caused by these pathogens. When shedding rates are also low, e.g., with macroparasitic infections, or when variable environmental factors play a role in transmissibility, then explicit modelling of both the pathogen and environmental transmission will provide a more accurate picture than a direct transmission approximation.

Universally valid reduction of multiscale stochastic biochemical systems using simple non-elementary propensities

Biochemical systems consist of numerous elementary reactions governed by the law of mass action. However, experimentally characterizing all the elementary reactions is nearly impossible. Thus, over a century, their deterministic models that typically contain rapid reversible bindings have been simplified with non-elementary reaction functions (e.g., Michaelis-Menten and Morrison equations). Although the non-elementary reaction functions are derived by applying the quasi-steady-state approximation (QSSA) to deterministic systems, they have also been widely used to derive propensities for stochastic simulations due to computational efficiency and simplicity. However, the validity condition for this heuristic approach has not been identified even for the reversible binding between molecules, such as protein-DNA, enzyme-substrate, and receptor-ligand, which is the basis for living cells. Here, we find that the non-elementary propensities based on the deterministic total QSSA can accurately capture the stochastic dynamics of the reversible binding in general. However, serious errors occur when reactant molecules with similar levels tightly bind, unlike deterministic systems. In that case, the non-elementary propensities distort the stochastic dynamics of a bistable switch in the cell cycle and an oscillator in the circadian clock. Accordingly, we derive alternative non-elementary propensities with the stochastic low-state QSSA, developed in this study. This provides a universally valid framework for simplifying multiscale stochastic biochemical systems with rapid reversible bindings, critical for efficient stochastic simulations of cell signaling and gene regulation. To facilitate the framework, we provide a user-friendly open-source computational package, ASSISTER, that automatically performs the present framework.

Beyond the adiabatic limit in systems with fast environments: A τ-leaping algorithm

We propose a τ-leaping simulation algorithm for stochastic systems subject to fast environmental changes. Similar to conventional τ-leaping the algorithm proceeds in discrete time steps, but as a principal addition it captures environmental noise beyond the adiabatic limit. The key idea is to treat the input rates for the τ-leaping as (clipped) Gaussian random variables with first and second moments constructed from the environmental process. In this way, each step of the algorithm retains environmental stochasticity to subleading order in the timescale separation between system and environment. We test the algorithm on several toy examples with discrete and continuous environmental states and find good performance in the regime of fast environmental dynamics. At the same time, the algorithm requires significantly less computing time than full simulations of the combined system and environment. In this context we also discuss several methods for the simulation of stochastic population dynamics in time-varying environments with continuous states.

Ultrasensitivity and bistability in covalent-modification cycles with positive autoregulation

Switch-like behaviours in biochemical networks are of fundamental significance in biological signal processing, and exist as two distinct types: ultra-sensitivity and bistability. Here we propose two new models of a reversible covalent-modification cycle with positive autoregulation (PAR), a motif structure that is thought to be capable of both ultrasensitivity and bistability in different parameter regimes. These new models appeal to a modelling framework that we call complex-complete, which accounts fully for the molecular complexities of the underlying signalling mechanisms. Each of the two new models encodes a specific molecular mechanism for PAR. We demonstrate that the modelling simplifications for PAR models that have been used in previous work, which rely on Michaelian approximations, are unable to accurately recapitulate the qualitative signalling responses supported by our detailed models. Strikingly, we show that complex-complete PAR models are capable of new qualitative responses such as one-way switches and a 'prozone' effect, depending on the specific PAR-encoding mechanism, which are not supported by Michaelian simplifications. Our results highlight the critical importance of accurately representing the molecular details of biochemical signalling mechanisms, and strongly suggest that the Michaelian approximation is inadequate for predictive models of enzyme-mediated chemical reactions with added regulations such as PAR.

A guide to the Michaelis-Menten equation: steady state and beyond

The modern definition of enzymology is synonymous with the Michaelis-Menten equation instituted by Leonor Michaelis and Maud Menten. Most textbooks, or chapters within, discussing enzymology start with the derivation of the equation under the assumption of rapid equilibrium (as done by Michaelis-Menten) or steady state (as modified by Briggs and Haldane) conditions to highlight the importance of this equation as the bedrock on which interpretation of enzyme kinetic results is dependent. However, few textbooks or monographs take the effort of placing the equation within its right historical context and discuss the assumptions that have gone into its institution. This guide will dwell on these in substantial detail. Further, this guide will attempt to instil a sense of appreciation for the mathematical curve rectangular hyperbola, its unique attributes and how ubiquitous the curve is in biological systems. To conclude, this guide will discuss the limitations of the equation, and the method it embodies, and trace the journey of how investigators are attempting to move beyond the steady-state approach and the Michaelis-Menten equation into full progress curve, pre-steady state and single-turnover kinetic analysis to obtain greater insights into enzyme kinetics and catalysis.

Multi-scale modeling of intensive macroalgae cultivation and marine nitrogen sequestration

Multi-scale macroalgae growth models are required for the efficient design of sustainable, economically viable, and environmentally safe farms. Here, we develop a multi-scale model for Ulva sp. macroalgae growth and nitrogen sequestration in an intensive cultivation farm, regulated by temperature, light, and nutrients. The model incorporates a range of scales by incorporating spatial effects in two steps: light extinction at the reactor scale (1 m) and nutrient absorption at the farm scale (1 km). The model was validated on real data from an experimental reactor installed in the sea. Biomass production rates, chemical compositions, and nitrogen removal were simulated under different seasons, levels of dilution in the environment and water-exchange rate in the reactor. This multi-scale model provides an important tool for environmental authorities and seaweed farmers who desire to upscale to large bioremediation and/or macroalgae biomass production farms, thus promoting the marine sustainable development and the macroalgae-based bioeconomy.

Zollmann et al. develop a multi-scale model for Ulva sp. macroalgae growth, biochemical composition, and nitrogen sequestration, regulated by temperature, light and nutrients. Their results demonstrate that this model can be used to design environmentally friendly and economically sustainable seaweed farms.

Allosteric conformational ensembles have unlimited capacity for integrating information

Integration of binding information by macromolecular entities is fundamental to cellular functionality. Recent work has shown that such integration cannot be explained by pairwise cooperativities, in which binding is modulated by binding at another site. Higher-order cooperativities (HOCs), in which binding is collectively modulated by multiple other binding events, appear to be necessary but an appropriate mechanism has been lacking. We show here that HOCs arise through allostery, in which effective cooperativity emerges indirectly from an ensemble of dynamically interchanging conformations. Conformational ensembles play important roles in many cellular processes but their integrative capabilities remain poorly understood. We show that sufficiently complex ensembles can implement any form of information integration achievable without energy expenditure, including all patterns of HOCs. Our results provide a rigorous biophysical foundation for analysing the integration of binding information through allostery. We discuss the implications for eukaryotic gene regulation, where complex conformational dynamics accompanies widespread information integration.

On rapid oscillations driving biological processes at disparate timescales

We consider a generic biological process described by a dynamical system, subject to an input signal with a high-frequency periodic component. The rapid oscillations of the input signal induce inherently multiscale dynamics, motivating order-reduction techniques. It is intuitive that the system behaviour is well approximated by its response to the averaged input signal. However, changes to the high-frequency component that preserve the average signal are beyond the reach of such intuitive reasoning. In this study, we explore system response under the influence of such an input signal by exploiting the timescale separation between high-frequency input variations and system response time. Employing the asymptotic method of multiple scales, we establish that, in some circumstances, the intuitive approach is simply the leading-order asymptotic contribution. We focus on higher-order corrections that capture the response to the details of the high-frequency component beyond its average. This approach achieves a reduction in system complexity while providing valuable insight into the structure of the response to the oscillations. We develop the general theory for nonlinear systems, while highlighting the important case of systems affine in the state and the input signal, presenting examples of both discrete and continuum state spaces. Importantly, this class of systems encompasses biochemical reaction networks described by the chemical master equation and its continuum approximations. Finally, we apply the framework to a nonlinear system describing mRNA translation and protein expression previously studied in the literature. The analysis shines new light on several aspects of the system quantification and both extends and simplifies results previously obtained.

An Early Stage Researcher's Primer on Systems Medicine Terminology

Background: Systems Medicine is a novel approach to medicine, that is, an interdisciplinary field that considers the human body as a system, composed of multiple parts and of complex relationships at multiple levels, and further integrated into an environment. Exploring Systems Medicine implies understanding and combining concepts coming from diametral different fields, including medicine, biology, statistics, modeling and simulation, and data science. Such heterogeneity leads to semantic issues, which may slow down implementation and fruitful interaction between these highly diverse fields. Methods: In this review, we collect and explain more than100 terms related to Systems Medicine. These include both modeling and data science terms and basic systems medicine terms, along with some synthetic definitions, examples of applications, and lists of relevant references. Results: This glossary aims at being a first aid kit for the Systems Medicine researcher facing an unfamiliar term, where he/she can get a first understanding of them, and, more importantly, examples and references for digging into the topic.

Stable maintenance of hidden switches as a strategy to increase the gene expression stability

In response to severe genetic and environmental perturbations, wild-type organisms can express hidden alternative phenotypes adaptive to such adverse conditions. While our theoretical understanding of the population-level fitness advantage and evolution of phenotypic switching under variable environments has grown, the mechanism by which these organisms maintain phenotypic switching capabilities under static environments remains to be elucidated. Here, using computational simulations, we analyzed the evolution of gene circuits under natural selection and found that different strategies evolved to increase the gene expression stability near the optimum level. In a population comprising bistable individuals, a strategy of maintaining bistability and raising the potential barrier separating the bistable regimes was consistently taken. Our results serve as evidence that hidden bistable switches can be stably maintained during environmental stasis-an essential property enabling the timely release of adaptive alternatives with small genetic changes in the event of substantial perturbations.

Gene Regulation in and out of Equilibrium

Determining whether and how a gene is transcribed are two of the central processes of life. The conceptual basis for understanding such gene regulation arose from pioneering biophysical studies in eubacteria. However, eukaryotic genomes exhibit vastly greater complexity, which raises questions not addressed by this bacterial paradigm. First, how is information integrated from many widely separated binding sites to determine how a gene is transcribed? Second, does the presence of multiple energy-expending mechanisms, which are absent from eubacterial genomes, indicate that eukaryotes are capable of improved forms of genetic information processing? An updated biophysical foundation is needed to answer such questions. We describe the linear framework, a graph-based approach to Markov processes, and show that it can accommodate many previous studies in the field. Under the assumption of thermodynamic equilibrium, we introduce a language of higher-order cooperativities and show how it can rigorously quantify gene regulatory properties suggested by experiment. We point out that fundamental limits to information processing arise at thermodynamic equilibrium and can only be bypassed through energy expenditure. Finally, we outline some of the mathematical challenges that must be overcome to construct an improved biophysical understanding of gene regulation.

Misuse of the Michaelis-Menten rate law for protein interaction networks and its remedy

For over a century, the Michaelis-Menten (MM) rate law has been used to describe the rates of enzyme-catalyzed reactions and gene expression. Despite the ubiquity of the MM rate law, it accurately captures the dynamics of underlying biochemical reactions only so long as it is applied under the right condition, namely, that the substrate is in large excess over the enzyme-substrate complex. Unfortunately, in circumstances where its validity condition is not satisfied, especially so in protein interaction networks, the MM rate law has frequently been misused. In this review, we illustrate how inappropriate use of the MM rate law distorts the dynamics of the system, provides mistaken estimates of parameter values, and makes false predictions of dynamical features such as ultrasensitivity, bistability, and oscillations. We describe how these problems can be resolved with a slightly modified form of the MM rate law, based on the total quasi-steady state approximation (tQSSA). Furthermore, we show that the tQSSA can be used for accurate stochastic simulations at a lower computational cost than using the full set of mass-action rate laws. This review describes how to use quasi-steady state approximations in the right context, to prevent drawing erroneous conclusions from in silico simulations.

Reversal symmetries for cyclic paths away from thermodynamic equilibrium

If a system is at thermodynamic equilibrium, an observer cannot tell whether a film of it is being played forward or in reverse: any transition will occur with the same frequency in the forward as in the reverse direction. However, if expenditure of energy changes the rate of even a single transition to yield a nonequilibrium steady state, such time-reversal symmetry undergoes a widespread breakdown, far beyond the point at which the energy is expended. An explosion of interdependency also arises, with steady-state probabilities of system states depending in a complicated manner on the rate of every transition in the system. Nevertheless, in the midst of this global nonequilibrium complexity, we find that cyclic paths have reversibility properties that remain local, and which can exhibit symmetry, no matter how far the system is from thermodynamic equilibrium. Specifically, given any cycle of reversible transitions, the ratio of the frequencies with which the cycle occurs in one direction versus the other is determined, in the long-time limit, only by the thermodynamic force on the cycle itself, without requiring knowledge of transition rates elsewhere in the system. In particular, if there is no net energy expenditure on the cycle, then, over long times, the cycle occurrence frequencies are the same in either direction.

Robustness and parameter geography in post-translational modification systems

Biological systems are acknowledged to be robust to perturbations but a rigorous understanding of this has been elusive. In a mathematical model, perturbations often exert their effect through parameters, so sizes and shapes of parametric regions offer an integrated global estimate of robustness. Here, we explore this “parameter geography” for bistability in post-translational modification (PTM) systems. We use the previously developed “linear framework” for timescale separation to describe the steady-states of a two-site PTM system as the solutions of two polynomial equations in two variables, with eight non-dimensional parameters. Importantly, this approach allows us to accommodate enzyme mechanisms of arbitrary complexity beyond the conventional Michaelis-Menten scheme, which unrealistically forbids product rebinding. We further use the numerical algebraic geometry tools Bertini, Paramotopy, and alphaCertified to statistically assess the solutions to these equations at ∼10⁹ parameter points in total. Subject to sampling limitations, we find no bistability when substrate amount is below a threshold relative to enzyme amounts. As substrate increases, the bistable region acquires 8-dimensional volume which increases in an apparently monotonic and sigmoidal manner towards saturation. The region remains connected but not convex, albeit with a high visibility ratio. Surprisingly, the saturating bistable region occupies a much smaller proportion of the sampling domain under mechanistic assumptions more realistic than the Michaelis-Menten scheme. We find that bistability is compromised by product rebinding and that unrealistic assumptions on enzyme mechanisms have obscured its parametric rarity. The apparent monotonic increase in volume of the bistable region remains perplexing because the region itself does not grow monotonically: parameter points can move back and forth between monostability and bistability. We suggest mathematical conjectures and questions arising from these findings. Advances in theory and software now permit insights into parameter geography to be uncovered by high-dimensional, data-centric analysis.

Biological organisms are often said to have robust properties but it is difficult to understand how such robustness arises from molecular interactions. Here, we use a mathematical model to study how the molecular mechanism of protein modification exhibits the property of multiple internal states, which has been suggested to underlie memory and decision making. The robustness of this property is revealed by the size and shape, or “geography,” of the parametric region in which the property holds. We use advances in reducing model complexity and in rapidly solving the underlying equations, to extensively sample parameter points in an 8-dimensional space. We find that under realistic molecular assumptions the size of the region is surprisingly small, suggesting that generating multiple internal states with such a mechanism is much harder than expected. While the shape of the region appears straightforward, we find surprising complexity in how the region grows with increasing amounts of the modified substrate. Our approach uses statistical analysis of data generated from a model, rather than from experiments, but leads to precise mathematical conjectures about parameter geography and biological robustness.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For ¹³C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

Kinetic control of stationary flux ratios for a wide range of biochemical processes

One of the most intriguing features of biological systems is their ability to regulate the steady-state fluxes of the underlying biochemical reactions; however, the regulatory mechanisms and their physicochemical properties are not fully understood. Fundamentally, flux regulation can be explained with a chemical kinetic formalism describing the transitions between discrete states, with the reaction rates defined by an underlying free energy landscape. Which features of the energy landscape affect the flux distribution? Here we prove that the ratios of the steady-state fluxes of quasi-first-order biochemical processes are invariant to energy perturbations of the discrete states and are only affected by the energy barriers. In other words, the nonequilibrium flux distribution is under kinetic and not thermodynamic control. We illustrate the generality of this result for three biological processes. For the network describing protein folding along competing pathways, the probabilities of proceeding via these pathways are shown to be invariant to the stability of the intermediates or to the presence of additional misfolded states. For the network describing protein synthesis, the error rate and the energy expenditure per peptide bond is proven to be independent of the stability of the intermediate states. For molecular motors such as myosin-V, the ratio of forward to backward steps and the number of adenosine 5'-triphosphate (ATP) molecules hydrolyzed per step is demonstrated to be invariant to energy perturbations of the intermediate states. These findings place important constraints on the ability of mutations and drug perturbations to affect the steady-state flux distribution for a wide class of biological processes.

A computational framework for a Lyapunov-enabled analysis of biochemical reaction networks

Complex molecular biological processes such as transcription and translation, signal transduction, post-translational modification cascades, and metabolic pathways can be described in principle by biochemical reactions that explicitly take into account the sophisticated network of chemical interactions regulating cell life. The ability to deduce the possible qualitative behaviors of such networks from a set of reactions is a central objective and an ongoing challenge in the field of systems biology. Unfortunately, the construction of complete mathematical models is often hindered by a pervasive problem: despite the wealth of qualitative graphical knowledge about network interactions, the form of the governing nonlinearities and/or the values of kinetic constants are hard to uncover experimentally. The kinetics can also change with environmental variations. This work addresses the following question: given a set of reactions and without assuming a particular form for the kinetics, what can we say about the asymptotic behavior of the network? Specifically, it introduces a class of networks that are "structurally (mono) attractive" meaning that they are incapable of exhibiting multiple steady states, oscillation, or chaos by virtue of their reaction graphs. These networks are characterized by the existence of a universal energy-like function called a Robust Lyapunov function (RLF). To find such functions, a finite set of rank-one linear systems is introduced, which form the extremals of a linear convex cone. The problem is then reduced to that of finding a common Lyapunov function for this set of extremals. Based on this characterization, a computational package, Lyapunov-Enabled Analysis of Reaction Networks (LEARN), is provided that constructs such functions or rules out their existence. An extensive study of biochemical networks demonstrates that LEARN offers a new unified framework. Basic motifs, three-body binding, and genetic networks are studied first. The work then focuses on cellular signalling networks including various post-translational modification cascades, phosphotransfer and phosphorelay networks, T-cell kinetic proofreading, and ERK signalling. The Ribosome Flow Model is also studied.

Dynamic Crowding Regulates Transcription

The nuclear environment is highly crowded by biological macromolecules, including chromatin and mobile proteins, which alter the kinetics and efficiency of transcriptional machinery. These alterations have been described, both theoretically and experimentally, for steady-state crowding densities; however, temporal changes in crowding density ("dynamic crowding") have yet to be integrated with gene expression. Dynamic crowding is pertinent to nuclear biology because processes such as chromatin translocation and protein diffusion lend to highly mobile biological crowders. Therefore, to capture such dynamic crowding and investigate its influence on transcription, we employ a three-pronged, systems-molecular approach. A system of chemical reactions represents the transcription pathway, the rates of which are determined by molecular-scale simulations; Brownian dynamics and Monte Carlo simulations quantify protein diffusion and DNA-protein binding affinity, dependent on macromolecular density. Altogether, this approach shows that transcription depends critically on dynamic crowding as the gene expression resultant from dynamic crowding can be profoundly different than that of steady-state crowding. In fact, expression levels can display both amplification and suppression and are notably different for genes or gene populations with different chemical and structural properties. These properties can be exploited to impose circadian expression, which is asymmetric and varies in strength, or to explain expression in cells under biomechanical stress. Therefore, this work demonstrates that dynamic crowding nontrivially alters transcription kinetics and presents dynamic crowding within the bulk nuclear nanoenvironment as a novel regulatory framework for gene expression.

On signalling and estimation limits for molecular birth-processes

Understanding and uncovering the mechanisms or motifs that molecular networks employ to regulate noise is a key problem in cell biology. As it is often difficult to obtain direct and detailed insight into these mechanisms, many studies instead focus on assessing the best precision attainable on the signalling pathways that compose these networks. Molecules signal one another over such pathways to solve noise regulating estimation and control problems. Quantifying the maximum precision of these solutions delimits what is achievable and allows hypotheses about underlying motifs to be tested without requiring detailed biological knowledge. The pathway capacity, which defines the maximum rate of transmitting information along it, is a widely used proxy for precision. Here it is shown, for estimation problems involving elementary yet biologically relevant birth-process networks, that capacity can be surprisingly misleading. A time-optimal signalling motif, called birth-following, is derived and proven to better the precision expected from the capacity, provided the maximum signalling rate constraint is large and the mean one above a certain threshold. When the maximum constraint is relaxed, perfect estimation is predicted by the capacity. However, the true achievable precision is found highly variable and sensitive to the mean constraint. Since the same capacity can map to different combinations of rate constraints, it can only equivocally measure precision. Deciphering the rate constraints on a signalling pathway may therefore be more important than computing its capacity.

A new Michaelis-Menten equation valid everywhere multi-scale dynamics prevails

The Michaelis-Menten reaction scheme is among the most influential models in the field of biochemistry, since it led to a very popular expression for the rate of an enzyme-catalysed reaction. After the realisation that this expression is valid in a limited region of the parameter space, two additional expressions were later introduced. The range of validity of these three expressions has been studied thoroughly, since the significance of a reliable rate is not based only on the accuracy of its predictive abilities but also on the physical insight that is acquired in the process of its construction. Here a new expression for the rate is introduced that is valid in practically the full parameter space and reduces to the expressions in the literature, when considering appropriate limits. The new expression is produced by employing algorithmic tools for asymptotic analysis, so that its construction is not hindered by the dimensional form and the complexity of the full model or by a wide parameter range of interest. These tools can be employed for the derivation of enzyme rate expressions of much more complex kinetics mechanisms.

Monostationarity and Multistationarity in Tree Networks of Goldbeter-Koshland Loops

A major challenge in systems biology is to elicit general properties in the face of molecular complexity. Here, we introduce a class of enzyme-catalysed biochemical networks and examine how the existence of a single positive steady state (monostationarity) depends on the network structure, enzyme mechanisms, kinetic rate laws and parameter values. We consider Goldbeter-Koshland (GK) covalent modification loops arranged in a tree network, so that a substrate form in one loop can be an enzyme in another loop. GK loops are a canonical motif in cell signalling and trees offer a generalisation of linear cascades which accommodate network complexity while remaining mathematically tractable. In particular, they permit a modular, recursive proof strategy which may be more widely applicable. We show that if each enzyme follows its own complex reaction mechanism under mass action kinetics, then any network is monostationary for all appropriate parameter values. If the kinetics is non-mass action with a plausible monotonicity requirement, and each enzyme follows the Michaelis-Menten mechanism, then monostationarity is preserved. Surprisingly, a single GK loop with a complex enzyme mechanism under non-mass action monotone kinetics can have more than one positive steady state (multistationarity). The broader interplay between network structure, enzyme mechanism and kinetics remains an intriguing open problem.

Modeling cardiac complexity: Advancements in myocardial models and analytical techniques for physiological investigation and therapeutic development in vitro

Cardiomyopathies, heart failure, and arrhythmias or conduction blockages impact millions of patients worldwide and are associated with marked increases in sudden cardiac death, decline in the quality of life, and the induction of secondary pathologies. These pathologies stem from dysfunction in the contractile or conductive properties of the cardiomyocyte, which as a result is a focus of fundamental investigation, drug discovery and therapeutic development, and tissue engineering. All of these foci require in vitro myocardial models and experimental techniques to probe the physiological functions of the cardiomyocyte. In this review, we provide a detailed exploration of different cell models, disease modeling strategies, and tissue constructs used from basic to translational research. Furthermore, we highlight recent advancements in imaging, electrophysiology, metabolic measurements, and mechanical and contractile characterization modalities that are advancing our understanding of cardiomyocyte physiology. With this review, we aim to both provide a biological framework for engineers contributing to the field and demonstrate the technical basis and limitations underlying physiological measurement modalities for biologists attempting to take advantage of these state-of-the-art techniques.

Negative reciprocity, not ordered assembly, underlies the interaction of Sox2 and Oct4 on DNA

The mode of interaction of transcription factors (TFs) on eukaryotic genomes remains a matter of debate. Single-molecule data in living cells for the TFs Sox2 and Oct4 were previously interpreted as evidence of ordered assembly on DNA. However, the quantity that was calculated does not determine binding order but, rather, energy expenditure away from thermodynamic equilibrium. Here, we undertake a rigorous biophysical analysis which leads to the concept of reciprocity. The single-molecule data imply that Sox2 and Oct4 exhibit negative reciprocity, with expression of Sox2 increasing Oct4’s genomic binding but expression of Oct4 decreasing Sox2’s binding. Models show that negative reciprocity can arise either from energy expenditure or from a mixture of positive and negative cooperativity at distinct genomic loci. Both possibilities imply unexpected complexity in how TFs interact on DNA, for which single-molecule methods provide novel detection capabilities.

The bodies of humans and other animals contain many types of cells that perform different roles in the body. Most cells in the body carry the same DNA, which is arranged into sections known as genes. The marked differences between cell types arise because different sets of genes are switched on or ‘expressed’.

Proteins called transcription factors control which genes are expressed by binding to DNA and recruiting groups of accessory proteins. However, it is not clear how they interact with each other and with the accessory proteins to decide whether to express a gene. For instance, it is thought that some accessory proteins may provide energy for this process, but it is unknown whether the energy is used continuously or only for a short time. Insights from physics suggest that the former could have more powerful effects.

In 2014, a team of researchers reported using a microscopy approach, known as single-molecule imaging, to follow two transcription factors called Sox2 and Oct4 in cells from mice. After analyzing the data, the researchers concluded that Sox2 and Oct4 had a specific order of binding to DNA, with Sox2 often binding first and then assisting Oct4 to bind. Now Biddle et al. report that the claim made in the 2014 study is unsupported because of errors in the way the data were analyzed. In particular, Biddle et al. argue that what the earlier study actually calculated is not an order of binding but a measure of whether energy is being continuously used when Sox2 and Oct4 bind to DNA.

Biddle et al. reanalyzed the data from the 2014 work and concluded that Sox2 increases the extent of Oct4 binding to DNA, while Oct4 decreases the amount of Sox2 binding to DNA. Mathematical models suggest this may be due to the continuous use of energy as the two proteins bind to DNA. Alternatively, it could also happen if Sox2 and Oct4 helped each other to bind at some sites on DNA and hindered each other from binding in other places, even if energy is only used for a short time. These findings reveal that there is unexpected complexity in how transcription factors bind to DNA.

The next step following on from this work is to carry out experiments that test the two possible explanations for how Sox2 and Oct4 interact on DNA. Including physics in the analysis may help describe more accurately the biology of how transcription factors determine gene expression. Understanding this process will shed new light on many important biological questions and may aid the development of gene therapy and other new medical techniques.

Dynamic Network Modeling of Stem Cell Metabolism

Stem cell metabolism is intrinsically tied to stem cell pluripotency and function. Yet, understanding metabolic rewiring in stem cells has been challenging due to the complex and highly interconnected nature of the metabolic network. Genome-scale metabolic network models are increasingly used to holistically model the metabolic behavior of various cells and tissues using transcriptomics data. However, these powerful approaches that model steady-state behavior have limited utility for studying dynamic stem cell state transitions. To address this complexity, we recently developed the dynamic flux activity (DFA) approach; DFA is a genome-scale modeling approach that uses time-course metabolic data to predict metabolic flux rewiring. This protocol outlines the steps for modeling steady-state and dynamic metabolic behavior using transcriptomics and time-course metabolomics data, respectively. Using data from naive and primed pluripotent stem cells, we demonstrate how we can use genome-scale modeling and DFA to comprehensively characterize the metabolic differences between these states.

Quantitative flux coupling analysis

Flux coupling analysis (FCA) aims to describe the functional dependencies among reactions in a metabolic network. Currently studied coupling relations are qualitative in the sense that they identify pairs of reactions for which the activity of one reaction necessitates the activity of the other one, but without giving any numerical bounds relating the possible activity rates. The potential applications of FCA are heavily investigated, however apart from some trivial cases there is no clue of what bottleneck in the metabolic network causes each dependency. In this article, we introduce a quantitative approach to the same flux coupling problem named quantitative flux coupling analysis (QFCA). It generalizes the current concepts as we show that all the qualitative information provided by FCA is readily available in the quantitative flux coupling equations of QFCA, without the need for any additional analysis. Moreover, we design a simple algorithm to efficiently identify these flux coupling equations which scales up to the genome-scale metabolic networks with thousands of reactions and metabolites in an effective way. Furthermore, this framework enables us to quantify the "strength" of the flux coupling relations. We also provide different biologically meaningful interpretations, including one which gives an intuitive certificate of precisely which metabolites in the network enforce each flux coupling relation. Eventually, we conclude by suggesting the probable application of QFCA to the metabolic gap-filling problem, which we only begin to address here and is left for future research to further investigate.

Cooperativity in proton sensing by PIP aquaporins

One of the most intriguing properties of plasma membrane intrinsic protein (PIP) aquaporins (AQPs) is their ability to modulate water transport by sensing different levels of intracellular pH through the assembly of homo- and heterotetrameric molecular species in the plasma membrane. In this work, using a phenomenological modeling approach, we demonstrate that cooperativity in PIP biological response cannot be directly attributed to a cooperative proton binding, as it is usually considered, since it could also be the consequence of a cooperative conformation transition between open and closed states of the channel. Moreover, our results show that, when mixed populations of homo- and heterotetrameric PIP channels are coexpressed in the plasma membrane of the same cell, the observed decrease in the degree of positive cooperativity would result from the simultaneous presence of molecular species with different levels of proton sensing. Indeed, the random mixing between different PIP paralogues as subunits in a single tetramer, plus the possibility of mixed populations of homo- and heterotetrameric PIP channels widen the spectrum of cooperative responses of a cell membrane. Our approach offers a deep understanding of cooperative transport of AQP channels, as members of a multiprotein family where the relevant proton binding sites of each member have not been clearly elucidated yet.

Structural conditions on complex networks for the Michaelis-Menten input-output response

The Michaelis-Menten (MM) fundamental formula describes how the rate of enzyme catalysis depends on substrate concentration. The familiar hyperbolic relationship was derived by timescale separation for a network of three reactions. The same formula has subsequently been found to describe steady-state input-output responses in many biological contexts, including single-molecule enzyme kinetics, gene regulation, transcription, translation, and force generation. Previous attempts to explain its ubiquity have been limited to networks with regular structure or simplifying parametric assumptions. Here, we exploit the graph-based linear framework for timescale separation to derive general structural conditions under which the MM formula arises. The conditions require a partition of the graph into two parts, akin to a "coarse graining" into the original MM graph, and constraints on where and how the input variable occurs. Other features of the graph, including the numerical values of parameters, can remain arbitrary, thereby explaining the formula's ubiquity. For systems at thermodynamic equilibrium, we derive a necessary and sufficient condition. For systems away from thermodynamic equilibrium, especially those with irreversible reactions, distinct structural conditions arise and a general characterization remains open. Nevertheless, our results accommodate, in much greater generality, all examples known to us in the literature.

Model discrimination for Ca2+ -dependent regulation of myosin light chain kinase in smooth muscle contraction

Excitation-contraction coupling in smooth muscle is mediated by the Ca²⁺ - and calmodulin-dependent regulation of myosin light chain kinase. The precise mechanism of this regulation remains controversial, and several mathematical models have been proposed for the interaction of the three species. These models have previously been analyzed at steady state primarily by numerical simulation of differential equations, for which parameter values must be estimated from data. Here, we use the linear framework for timescale separation to demonstrate that models of this general kind can be solved analytically for an equilibrium steady state, without having to determine parameter values. This analysis leads to parameter-independent methods for discriminating between the models, for which we propose experiments that could be performed with existing methods.

Energy-speed-accuracy relation in complex networks for biological discrimination

Discriminating between correct and incorrect substrates is a core process in biology, but how is energy apportioned between the conflicting demands of accuracy (μ), speed (σ), and total entropy production rate (P)? Previous studies have focused on biochemical networks with simple structure or relied on simplifying kinetic assumptions. Here, we use the linear framework for timescale separation to analytically examine steady-state probabilities away from thermodynamic equilibrium for networks of arbitrary complexity. We also introduce a method of scaling parameters that is inspired by Hopfield's treatment of kinetic proofreading. Scaling allows asymptotic exploration of high-dimensional parameter spaces. We identify in this way a broad class of complex networks and scalings for which the quantity σln(μ)/P remains asymptotically finite whenever accuracy improves from equilibrium, so that μ_{eq}/μ→0. Scalings exist, however, even for Hopfield's original network, for which σln(μ)/P is asymptotically infinite, illustrating the parametric complexity. Outside the asymptotic regime, numerical calculations suggest that, under more restrictive parametric assumptions, networks satisfy the bound, σln(μ/μ_{eq})/P<1, and we discuss the biological implications for discrimination by ribosomes and DNA polymerase. The methods introduced here may be more broadly useful for analyzing complex networks that implement other forms of cellular information processing.

Abstracting the dynamics of biological pathways using information theory: a case study of apoptosis pathway

Motivation

Quantitative models are increasingly used in systems biology. Usually, these quantitative models involve many molecular species and their associated reactions. When simulating a tissue with thousands of cells, using these large models becomes computationally and time limiting.

Results

In this paper, we propose to construct abstractions using information theory notions. Entropy is used to discretize the state space and mutual information is used to select a subset of all original variables and their mutual dependencies. We apply our method to an hybrid model of TRAIL-induced apoptosis in HeLa cell. Our abstraction, represented as a Dynamic Bayesian Network (DBN), reduces the number of variables from 92 to 10, and accelerates numerical simulation by an order of magnitude, yet preserving essential features of cell death time distributions.

Availability and Implementation

This approach is implemented in the tool DBNizer, freely available at http://perso.crans.org/genest/DBNizer .

Contact

gregory.batt@inria.fr or bgenest@irisa.fr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Steady-State Differential Dose Response in Biological Systems

In pharmacology and systems biology, it is a fundamental problem to determine how biological systems change their dose-response behavior upon perturbations. In particular, it is unclear how topologies, reactions, and parameters (differentially) affect the dose response. Because parameters are often unknown, systematic approaches should directly relate network structure and function. Here, we outline a procedure to compare general non-monotone dose-response curves and subsequently develop a comprehensive theory for differential dose responses of biochemical networks captured by non-equilibrium steady-state linear framework models. Although these models are amenable to analytical derivations of non-equilibrium steady states in principle, their size frequently increases (super) exponentially with model size. We extract general principles of differential responses based on a model's graph structure and thereby alleviate the combinatorial explosion. This allows us, for example, to determine reactions that affect differential responses, to identify classes of networks with equivalent differential, and to reject hypothetical models reliably without needing to know parameter values. We exemplify such applications for models of insulin signaling.

From word models to executable models of signaling networks using automated assembly

Word models (natural language descriptions of molecular mechanisms) are a common currency in spoken and written communication in biomedicine but are of limited use in predicting the behavior of complex biological networks. We present an approach to building computational models directly from natural language using automated assembly. Molecular mechanisms described in simple English are read by natural language processing algorithms, converted into an intermediate representation, and assembled into executable or network models. We have implemented this approach in the Integrated Network and Dynamical Reasoning Assembler (INDRA), which draws on existing natural language processing systems as well as pathway information in Pathway Commons and other online resources. We demonstrate the use of INDRA and natural language to model three biological processes of increasing scope: (i) p53 dynamics in response to DNA damage, (ii) adaptive drug resistance in BRAF‐V600E‐mutant melanomas, and (iii) the RAS signaling pathway. The use of natural language makes the task of developing a model more efficient and it increases model transparency, thereby promoting collaboration with the broader biology community.