Overcoming the stress of starvation is one of an organism’s most challenging phenotypic responses. Those organisms that frequently survive the challenge, by virtue of their fitness, will have evolved genomes that are shaped by their specific environments. Understanding this genotype–environment–phenotype relationship at a deep level will require quantitative predictive models of the complex molecular systems that link these aspects of an organism’s existence. Here, we treat one of the most fundamental molecular systems, protein synthesis, and the amino acid biosynthetic pathways involved in the stringent response to starvation. These systems face an inherent logical dilemma: Building an amino acid biosynthetic pathway to synthesize its product—the cognate amino acid of the pathway—may require that very amino acid when it is no longer available. To study this potential “catch-22,” we have created a generic model of amino acid biosynthesis in response to sudden starvation. Our mathematical analysis and computational results indicate that there are two distinctly different outcomes: Partial recovery to a new steady state, or full system failure. Moreover, the cell’s fate is dictated by the cognate bias, the number of cognate amino acids in the corresponding biosynthetic pathway relative to the average number of that amino acid in the proteome. We test these implications by analyzing the proteomes of over 1,800 sequenced microbes, which reveals statistically significant evidence of low cognate bias, a genetic trait that would avoid the biosynthetic quandary. Furthermore, these results suggest that the pattern of cognate bias, which is readily derived by genome sequencing, may provide evolutionary clues to an organism’s natural environment.

Proteome

Model Organism

10.1093/molbev/msu225

Cite

Citations (2)

The tricarboxylic acid cycle in Dictyostelium discoideum. Systemic effects of including protein turnover in the current model

Journal of Biological Chemistry (1993)

Fumihidé Shiraishi Michael A. Savageau

The current model for the tricarboxylic acid cycle in Dictyostelium discoideum is based on extensive experimental studies of enzyme kinetics in vitro and of metabolite fluxes measured in vivo. In the previous papers (Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22912-22918; 22919-22925; 22926-22933; 22934-22943) of this series we have carried out extensive analyses of the current model within the framework of biochemical systems theory with a view toward understanding the behavior of the integrated system. The model was found to be ill determined with respect to at least three of its features. In this paper we propose a minimal modification in the model that is consistent with previous experimental data but also includes recycling of amino acids for protein synthesis, one of the neglected features identified as important in the previous analysis. We again perform an analysis within the framework of biochemical systems theory to determine the systemic consequences of this change. The results show that the robustness of the modified model, as determined by the parameter sensitivities, is improved by 2 orders of magnitude over that of the previous model. Analysis of the dynamics shows that the turnover times for the pools of alanine, glutamate, and aspartate are reduced by 2 orders of magnitude and made more physiologically realistic. The distribution of flux is no longer rigidly fixed, and problems previously centered on the metabolism of pyruvate have been partially alleviated. Continued discrepancies lead us to question the degree to which kinetic data obtained with purified enzymes in vitro faithfully reflect the kinetic behavior of the integrated enzyme system in vivo. We must continue to re-examine the manner in which the kinetics of reactions in vivo are represented and to reassess the physical conditions that prevail in vitro and in vivo. Results in this paper direct our attention toward specific aspects of the system where these efforts should be focused. Thus, a minimal modification of the previous model has led to several improvements that make it more representative of the tricarboxylic acid cycle in D. discoideum, and the analysis in this paper leads to further predictions for improving the model.

Tricarboxylic acid

Protein turnover

10.1016/s0021-9258(19)85283-3

Cite

Citations (39)

Recasting nonlinear differential equations as S-systems: a canonical nonlinear form

Mathematical Biosciences (1987)

Michael A. Savageau Eberhard O. Voit

Canonical form

10.1016/0025-5564(87)90035-6

Cite

Citations (306)

Elucidating the genotype–phenotype map by automatic enumeration and analysis of the phenotypic repertoire

npj Systems Biology and Applications (2015)

Jason G. Lomnitz Michael A. Savageau

Abstract Background: The gap between genotype and phenotype is filled by complex biochemical systems most of which are poorly understood. Because these systems are complex, it is widely appreciated that quantitative understanding can only be achieved with the aid of mathematical models. However, formulating models and measuring or estimating their numerous rate constants and binding constants is daunting. Here we present a strategy for automating difficult aspects of the process. Methods: The strategy, based on a system design space methodology, is applied to a class of 16 designs for a synthetic gene oscillator that includes seven designs previously formulated on the basis of experimentally measured and estimated parameters. Results: Our strategy provides four important innovations by automating: (1) enumeration of the repertoire of qualitatively distinct phenotypes for a system; (2) generation of parameter values for any particular phenotype; (3) simultaneous realization of parameter values for several phenotypes to aid visualization of transitions from one phenotype to another, in critical cases from functional to dysfunctional; and (4) identification of ensembles of phenotypes whose expression can be phased to achieve a specific sequence of functions for rationally engineering synthetic constructs. Our strategy, applied to the 16 designs, reproduced previous results and identified two additional designs capable of sustained oscillations that were previously missed. Conclusions: Starting with a system’s relatively fixed aspects, its architectural features, our method enables automated analysis of nonlinear biochemical systems from a global perspective, without first specifying parameter values. The examples presented demonstrate the efficiency and power of this automated strategy.

Enumeration

Identification

Realization (probability)

10.1038/npjsba.2015.3

Cite

Citations (18)

Effects of alternative connectivity on behavior of randomly constructed Boolean networks

Physica D Nonlinear Phenomena (2002)

Chikoo Oosawa Michael A. Savageau

Boolean network

Rank (graph theory)

Boolean model

Gene regulatory network

10.1016/s0167-2789(02)00530-4

Cite

Citations (75)

Transport of Biosynthetic Intermediates: Regulation of Homoserine and Threonine Uptake in Escherichia coli

Journal of Bacteriology (1974)

Bonnie Templeton Michael A. Savageau

Homoserine is transported by a single system that it shares with alanine, isoleucine, leucine, phenylalanine, threonine, valine and perhaps cysteine, methionine, serine, and tyrosine. We investigated the regulation of this transport system and found that alanine, isoleucine, leucine, methionine, and valine each repress the homoserine-transporting system. From the concentration resulting in 50% repression of this transport system and the maximal amount of repression, we ranked the amino acids according to their effectiveness in repressing homoserine transport (in decreasing order): leucine>methionine>alanine>valine>isoleucine. The exponential rate of decrease in transport capacity after leucine addition equals the exponential growth rate of the culture, and protein synthesis is necessary for the derepression seen when leucine is removed. Threonine, in addition to using the above system, is transported by a second system shared with serine. We present further evidence for this serine-threonine transport system and show that it is not regulated like the homoserine-transporting system.

Homoserine

Isoleucine

Alanine

Derepression

Amino acid synthesis

10.1128/jb.120.1.114-120.1974

Cite

Citations (26)