Quantitative computational models play an increasingly important role in modern biology. Such models typically involve many free parameters, and assigning their values is often a substantial obstacle to model development. Directly measuring in vivo biochemical parameters is difficult, and collectively fitting them to other experimental data often yields large parameter uncertainties. Nevertheless, in earlier work we showed in a growth-factor-signaling model that collective fitting could yield well-constrained predictions, even when it left individual parameters very poorly constrained. We also showed that the model had a ‘‘sloppy’’ spectrum of parameter sensitivities, with eigenvalues roughly evenly distributed over many decades. Here we use a collection of models from the literature to test whether such sloppy spectra are common in systems biology. Strikingly, we find that every model we examine has a sloppy spectrum of sensitivities. We also test several consequences of this sloppiness for building predictive models. In particular, sloppiness suggests that collective fits to even large amounts of ideal time-series data will often leave many parameters poorly constrained. Tests over our model collection are consistent with this suggestion. This difficulty with collective fits may seem to argue for direct parameter measurements, but sloppiness also implies that such measurements must be formidably precise and complete to usefully constrain many model predictions. We confirm this implication in our growth-factor-signaling model. Our results suggest that sloppy sensitivity spectra are universal in systems biology models. The prevalence of sloppiness highlights the power of collective fits and suggests that modelers should focus on predictions rather than on parameters.

https://www.lassp.cornell.edu/sethna/pubPDF/SloppyEverywhere.pdf

Universally sloppy parameter sensitivities in systems biology models.

https://digital.library.unt.edu/ark:/67531/metadc488175/m2/1/high_res_d/2003_phytochemistry_817.pdf

Plant metabolomics: large-scale phytochemistry in the functional genomics era

Complex assemblies of interacting proteins carry out most of the interesting jobs in a cell, such as metabolism, DNA synthesis, movement and information processing. These physiological properties play out as a subtle molecular dance, choreographed by underlying regulatory networks. To understand this dance, a new breed of theoretical molecular biologists reproduces these networks in computers and in the mathematical language of dynamical systems.

/pdf/network-dynamics-and-cell-physiology-4ydb3fmmx2.pdf

Network dynamics and cell physiology

Eukaryotic cells have evolved various response mechanisms to counteract the deleterious consequences of oxidative stress Among these processes, metabolic alterations seem to play an important role We recently discovered that yeast cells with reduced activity of the key glycolytic enzyme triosephosphate isomerase exhibit an increased resistance to the thiol-oxidizing reagent diamide Here we show that this phenotype is conserved in Caenorhabditis elegans and that the underlying mechanism is based on a redirection of the metabolic flux from glycolysis to the pentose phosphate pathway, altering the redox equilibrium of the cytoplasmic NADP(H) pool Remarkably, another key glycolytic enzyme, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), is known to be inactivated in response to various oxidant treatments, and we show that this provokes a similar redirection of the metabolic flux The naturally occurring inactivation of GAPDH functions as a metabolic switch for rerouting the carbohydrate flux to counteract oxidative stress As a consequence, altering the homoeostasis of cytoplasmic metabolites is a fundamental mechanism for balancing the redox state of eukaryotic cells under stress conditions

/pdf/dynamic-rerouting-of-the-carbohydrate-flux-is-key-to-ntdzqa95q1.pdf

Dynamic rerouting of the carbohydrate flux is key to counteracting oxidative stress

Because of its profound effects on the rates of biological processes such as aerobic metabolism, environmental temperature plays an important role in shaping the distribution and abundance of species. As temperature increases, the rate of metabolism increases and then rapidly declines at higher temperatures – a response that can be described using a thermal performance curve (TPC). Although the shape of the TPC for aerobic metabolism is often attributed to the competing effects of thermodynamics, which can be described using the Arrhenius equation, and the effects of temperature on protein stability, this account represents an over-simplification of the factors acting even at the level of single proteins. In addition, it cannot adequately account for the effects of temperature on complex multistep processes, such as aerobic metabolism, that rely on mechanisms acting across multiple levels of biological organization. The purpose of this review is to explore our current understanding of the factors that shape the TPC for aerobic metabolism in response to acute changes in temperature, and to highlight areas where this understanding is weak or insufficient. Developing a more strongly grounded mechanistic model to account for the shape of the TPC for aerobic metabolism is crucial because these TPCs are the foundation of several recent attempts to predict the responses of species to climate change, including the metabolic theory of ecology and the hypothesis of oxygen and capacity-limited thermal tolerance.

The effects of temperature on aerobic metabolism: towards a mechanistic understanding of the responses of ectotherms to a changing environment

This paper examines whether the in vivo behavior of yeast glycolysis can be understood in terms of the in vitro kinetic properties of the constituent enzymes. In nongrowing, anaerobic, compressed Saccharomyces cerevisiae the values of the kinetic parameters of most glycolytic enzymes were determined. For the other enzymes appropriate literature values were collected. By inserting these values into a kinetic model for glycolysis, fluxes and metabolites were calculated. Under the same conditions fluxes and metabolite levels were measured. In our first model, branch reactions were ignored. This model failed to reach the stable steady state that was observed in the experimental flux measurements. Introduction of branches towards trehalose, glycogen, glycerol and succinate did allow such a steady state. The predictions of this branched model were compared with the empirical behavior. Half of the enzymes matched their predicted flux in vivo within a factor of 2. For the other enzymes it was calculated what deviation between in vivo and in vitro kinetic characteristics could explain the discrepancy between in vitro rate and in vivo flux.

Corinne A. Reijenga

Papers

Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry.