Quantifying the propagation of parametric uncertainty on flux balance analysis

Flux balance analysis (FBA) and associated techniques operating on stoichiometric genome-scale metabolic models play a central role in quantifying metabolic flows and constraining feasible phenotypes. At the heart of these methods lie two important assumptions: (i) the biomass precursors and energy requirements neither change in response to growth conditions nor environmental/genetic perturbations, and (ii) metabolite production and consumption rates are equal at all times (i.e., steady-state). Despite the stringency of these two assumptions, FBA has been shown to be surprisingly robust at predicting cellular phenotypes. In this paper, we formally assess the impact of these two assumptions on FBA results by quantifying how uncertainty in biomass reaction coefficients, and departures from steady-state due to temporal fluctuations could propagate to FBA results. In the first case, conditional sampling of parameter space is required to re-weigh the biomass reaction so as the molecular weight remains equal to 1 g mmol⁻¹, and in the second case, metabolite (and elemental) pool conservation must be imposed under temporally varying conditions. Results confirm the importance of enforcing the aforementioned constraints and explain the robustness of FBA biomass yield predictions.     

Predictive Potential of Flux Balance Analysis of Saccharomyces cerevisiae Using as Optimization Function Combinations of Cell Compartmental Objectives

Background

The main objective of flux balance analysis (FBA) is to obtain quantitative predictions of metabolic fluxes of an organism, and it is necessary to use an appropriate objective function to guarantee a good estimation of those fluxes.

Methodology

In this study, the predictive performance of FBA was evaluated, using objective functions arising from the linear combination of different cellular objectives. This approach is most suitable for eukaryotic cells, owing to their multiplicity of cellular compartments. For this reason, Saccharomyces cerevisiae was used as model organism, and its metabolic network was represented using the genome-scale metabolic model iMM904. As the objective was to evaluate the predictive performance from the FBA using the kind of objective function previously described, substrate uptake and oxygen consumption were the only input data used for the FBA. Experimental information about microbial growth and exchange of metabolites with the environment was used to assess the quality of the predictions.

Conclusions

The quality of the predictions obtained with the FBA depends greatly on the knowledge of the oxygen uptake rate. For the most of studied classifications, the best predictions were obtained with “maximization of growth”, and with some combinations that include this objective. However, in the case of exponential growth with unknown oxygen exchange flux, the objective function “maximization of growth, plus minimization of NADH production in cytosol, plus minimization of NAD(P)H consumption in mitochondrion” gave much more accurate estimations of fluxes than the obtained with any other objective function explored in this study.     

Growth of an anaerobic sulfate-reducing bacterium sustained by oxygen respiratory energy conservation after O2-driven experimental evolution

Desulfovibrio species are representatives of microorganisms at the boundary between anaerobic and aerobic lifestyles, since they contain the enzymatic systems required for both sulfate and oxygen reduction. However, the latter has been shown to be solely a protective mechanism. By implementing the oxygen-driven experimental evolution of Desulfovibrio vulgaris Hildenborough, we have obtained strains that have evolved to grow with energy derived from oxidative phosphorylation linked to oxygen reduction. We show that a few mutations are sufficient for the emergence of this phenotype and reveal two routes of evolution primarily involving either inactivation or overexpression of the gene encoding heterodisulfide reductase. We propose that the oxygen respiration for energy conservation that sustains the growth of the O₂-evolved strains is associated with a rearrangement of metabolite fluxes, especially NAD⁺/NADH, leading to an optimized O₂ reduction. These evolved strains are the first sulfate-reducing bacteria that exhibit a demonstrated oxygen respiratory process that enables growth.     

Effects of sodium and potassium ions on oxidative phosphorylation in relation to respiratory control by a cell-membrane adenosine triphosphatase

1. A study has been made of the oxygen consumption of kidney homogenates in relation to the ADP concentration as regulated by the cell-membrane adenosine triphosphatase. Stimulation of this enzymic activity by Na(+) and K(+) caused parallel increases in oxygen consumption and ADP concentration. Similarly, inhibition with ouabain caused a parallel fall. The membrane adenosine triphosphatase concerned in active transport therefore appears to regulate respiration through its control of ADP concentration. 2. The respiration of homogenates and mitochondria was also stimulated by K(+) in a way independent of adenosine-triphosphatase activity. It was shown that K(+) facilitates oxidative phosphorylation and the respiratory response to ADP. A K(+) concentration of 25-50mm was needed for maximum oxidative phosphorylation in the presence of physiological concentration of Na(+). Na(+) counteracted K(+) in the effects on mitochondria. It is concluded that K(+) regulates cellular respiration at two structures, one directly in mitochondria, and the second indirectly through control of ADP production at the cell membrane.     

In silico design and adaptive evolution of Escherichia coli for production of lactic acid

The development and validation of new methods to help direct rational strain design for metabolite overproduction remains an important problem in metabolic engineering. Here we show that computationally predicted E. coli strain designs, calculated from a genome-scale metabolic model, can lead to successful production strains and that adaptive evolution of the engineered strains can lead to improved production capabilities. Three strain designs for lactate production were implemented yielding a total of 11 evolved production strains that were used to demonstrate the utility of this integrated approach. Strains grown on 2 g/L glucose at 37 degrees C showed lactate titers ranging from 0.87 to 1.75 g/L and secretion rates that were directly coupled to growth rates.     

Improving the accuracy of flux balance analysis through the implementation of carbon availability constraints for intracellular reactions

Constraint-based modeling methods, such as Flux Balance Analysis (FBA), have been extensively used to decipher complex, information rich -omics datasets to elicit system-wide behavioral patterns of cellular metabolism. FBA has been successfully used to gain insight in a wide range of applications, such as range of substrate utilization, product yields and to design metabolic engineering strategies to improve bioprocess performance. A well-known challenge associated with large genome-scale metabolic networks is that they result in underdetermined problem formulations. Consequently, rather than unique solutions, FBA and related methods examine ranges of reaction flux values that are consistent with the studied physiological conditions. The wider the reported flux ranges, the higher the uncertainty in the determination of basic reaction properties, limiting interpretability of and confidence in the results. Herein, we propose a new, computationally efficient approach that refines flux range predictions by constraining reaction fluxes on the basis of the elemental balance of carbon. We compared carbon constraint FBA (ccFBA) against experimentally-measured intracellular fluxes using the latest CHO GEM (iCHO1766) and were able to substantially improve the accuracy of predicted flux values compared with FBA. ccFBA can be used as a stand-alone method but is also compatible with and complimentary to other constraint-based approaches.     

Molecular cloning, primary structure and disruption of the structural gene of aldolase from Saccharomyces cerevisiae

A yeast cDNA genetic library in a bacteriophage expression vector was screened using an antiserum reacting with fructose 1,6-bisphosphate aldolase from Saccharomyces cerevisiae. Radio-labelled probes of selected immunopositive clones were used for screening of a yeast genomic library. From the genomic clones a yeast/Escherichia coli shuttle plasmid was constructed containing on a 1990-base-pair fragment the entire structural gene FBA1 coding for yeast aldolase. The primary structure of the FBA1 gene was determined. An open reading frame comprises 1077 base pairs coding for a protein of 359 amino acids with a predicted molecular mass of 39,608 Da. As observed for other strongly expressed yeast genes, codon usage is extremely biased. The 810 base pairs at the 5' end and the 90 base pairs at the 3' end of the coding region of the cloned FBA1 gene are sufficient for normal expression and show characteristic elements present in the noncoding sequences of other yeast genes. Aldolase is the major protein in yeast cells transformed with a high-copy-number plasmid containing the FBA1 gene. The aldolase gene was disrupted by insertion of the yeast URA3 gene into the coding region of one FBA1 allele in a homozygous diploid ura3 strain. The haploid offsprings with the defective aldolase allele fba1::URA3 lack aldolase enzymatic activity and fail to grow in media containing as a carbon source metabolites of only one side of the aldolase reaction.     

Ontogeny of the Alimentary Canal and Respiratory Physiology of Larval Hoplosternum littorale (Hancock, 1828): an Intestinal Air-Breathing Teleost

Synopsis Newly hatched larval Hoplosternum littorale do not exhibit bimodal respiration upon hatching but depend initially on branchial respiration. As the respiratory intestine develops, its increasing length leads to an increasing dependence on aerial respiration. The respiratory intestine first appears when the fish is 9 days old but is microscopic at this age and becomes detectable by gross morphological examination at 13 days old. Development of the gut for the remainder of the juvenile period (up to 32 days old) consists of increasing length, mucosal folding, concentration of taste buds, gastric glands, mucous glands and goblet cells. Stage 1 commenced with newly hatched larvae until the larval fish were 11 days old. During this stage the fish did not have the capacity to air-breathe and thus behaved as oxygen conformers. Stage 2 extended from day 12 to day 23 and the fish had the capacity to air-breathe although the respiratory intestine was still developing. Stage 3 began from day 24 to day 32 and at this point the respiratory intestine was well-developed and functioned efficiently. The rate of oxygen consumption of the larval fish was found to be related to dissolved oxygen concentration, age of the larval fish and various interactions of these two factors.     

The Ability of Flux Balance Analysis to Predict Evolution of Central Metabolism Scales with the Initial Distance to the Optimum

The most powerful genome-scale framework to model metabolism, flux balance analysis (FBA), is an evolutionary optimality model. It hypothesizes selection upon a proposed optimality criterion in order to predict the set of internal fluxes that would maximize fitness. Here we present a direct test of the optimality assumption underlying FBA by comparing the central metabolic fluxes predicted by multiple criteria to changes measurable by a ¹³C-labeling method for experimentally-evolved strains. We considered datasets for three Escherichia coli evolution experiments that varied in their length, consistency of environment, and initial optimality. For ten populations that were evolved for 50,000 generations in glucose minimal medium, we observed modest changes in relative fluxes that led to small, but significant decreases in optimality and increased the distance to the predicted optimal flux distribution. In contrast, seven populations evolved on the poor substrate lactate for 900 generations collectively became more optimal and had flux distributions that moved toward predictions. For three pairs of central metabolic knockouts evolved on glucose for 600–800 generations, there was a balance between cases where optimality and flux patterns moved toward or away from FBA predictions. Despite this variation in predictability of changes in central metabolism, two generalities emerged. First, improved growth largely derived from evolved increases in the rate of substrate use. Second, FBA predictions bore out well for the two experiments initiated with ancestors with relatively sub-optimal yield, whereas those begun already quite optimal tended to move somewhat away from predictions. These findings suggest that the tradeoff between rate and yield is surprisingly modest. The observed positive correlation between rate and yield when adaptation initiated further from the optimum resulted in the ability of FBA to use stoichiometric constraints to predict the evolution of metabolism despite selection for rate.     

Respiration-Dependent Utilization of Sugars in Yeasts: a Determinant Role for Sugar Transporters

In many yeast species, including Kluyveromyces lactis, growth on certain sugars (such as galactose, raffinose, and maltose) occurs only under respiratory conditions. If respiration is blocked by inhibitors, mutation, or anaerobiosis, growth does not take place. This apparent dependence on respiration for the utilization of certain sugars has often been suspected to be associated with the mechanism of the sugar uptake step. We hypothesized that in many yeast species, the permease activities for these sugars are not sufficient to ensure the high substrate flow that is necessary for fermentative growth. By introducing additional sugar permease genes, we have obtained K. lactis strains that were capable of growing on galactose and raffinose in the absence of respiration. High dosages of both the permease and maltase genes were indeed necessary for K. lactis cells to grow on maltose in the absence of respiration. These results strongly suggest that the sugar uptake step is the major bottleneck in the fermentative assimilation of certain sugars in K. lactis and probably in many other yeasts.     

A unified stoichiometric model for oxidative and oxidoreductive growth of yeasts

Yeasts degrade glucose through different metabolic pathways, where the choice of the pathway is dependent on the nature of the limitation in the various substrates. When oxygen is limiting in addition to glucose, yeasts often grow according to a mixture of oxidative and reductive metabolism. Oxygen may be limiting either by supply or by inherent biological restrictions such as the respiratory bottleneck in Saccharomyces cerevisiae or by both. A unified model incorporating both supply and biological limitations is proposed for the quantitative prediction of growth rates, consumption and production rates, as well as key metabolite concentrations during mixed oxidoreductive metabolism occuring as a result of such oxygen limitations. This simple unstructured model can be applied to different yeast strains while at the same time requiring a minimum number of measured parameters. "Estimators" are utilized in order to predict the presence of supply-side or biological limitations. The values of these estimators also characterize the relative importance of oxidative to total metabolism. Results from the aerobic and oxygen-limited chemostat cultures were used to corroborate the model predictions. During these experiments, the heat released by the yeast cultures was also monitored on-line. The model correctly predicted the overall stoichiometry, steady-state concentrations, and rates including heat dissipation rates measured in the various situations of oxygen limitations. Direct continuous measurements such as heat can be used in conjunction with the unified model for on-line proces control. (c) 1992 John Wiley & Sons, Inc.     

Suppression of mycelia formation by NO produced endogenously in Candida tropicalis

In the experiments reported here we found that enzymatic NO synthesis in the yeast Candida tropicalis resembles the one in animal tissues with respect to the substrate arginine as well as its sensitivity to potential competitive inhibitors. Both, NO produced by the yeast's nitric oxide synthase and NO derived from an artificial donor, suppressed the formation of pseudomycelia. These results suggest to make use of NO as a tool in elucidating the mechanism controlling mycelia generation in this yeast. The apparent K(m) towards oxygen of the yeast's nitric oxide synthase (about 50 microM) was found to be high as compared to the apparent K(m) value of the yeast's respiratory chain (about 170 nM). From this observation it may be concluded that under conditions of little oxygen supply the nitric oxide synthase will unsuccessfully compete for oxygen with respiration. Therefore, the formation of mycelia spontaneously occurring in yeast cultures grown in sealed chambers can be attributed to a reduced internal NO level rather than limited respiratory activity.     

Mitochondrial Morphology and Fundamental Parameters of the Mitochondrial Respiratory Chain Are Altered in Caenorhabditis elegans Strains Deficient in Mitochondrial Dynamics and Homeostasis Processes

Mitochondrial dysfunction has been linked to myriad human diseases and toxicant exposures, highlighting the need for assays capable of rapidly assessing mitochondrial health in vivo. Here, using the Seahorse XF^e24 Analyzer and the pharmacological inhibitors dicyclohexylcarbodiimide and oligomycin (ATP-synthase inhibitors), carbonyl cyanide 4-(trifluoromethoxy) phenylhydrazone (mitochondrial uncoupler) and sodium azide (cytochrome c oxidase inhibitor), we measured the fundamental parameters of mitochondrial respiratory chain function: basal oxygen consumption, ATP-linked respiration, maximal respiratory capacity, spare respiratory capacity and proton leak in the model organism Caenhorhabditis elegans. Since mutations in mitochondrial homeostasis genes cause mitochondrial dysfunction and have been linked to human disease, we measured mitochondrial respiratory function in mitochondrial fission (drp-1)-, fusion (fzo-1)-, mitophagy (pdr-1, pink-1)-, and electron transport chain complex III (isp-1)-deficient C. elegans. All showed altered function, but the nature of the alterations varied between the tested strains. We report increased basal oxygen consumption in drp-1; reduced maximal respiration in drp-1, fzo-1, and isp-1; reduced spare respiratory capacity in drp-1 and fzo-1; reduced proton leak in fzo-1 and isp-1; and increased proton leak in pink-1 nematodes. As mitochondrial morphology can play a role in mitochondrial energetics, we also quantified the mitochondrial aspect ratio for each mutant strain using a novel method, and for the first time report increased aspect ratios in pdr-1- and pink-1-deficient nematodes.     

Protamine inhibition of the oxidative phosphorylation in intact, cytochrome c-depleted and restored mitochondria.

On the basis of polarographic data it is shown that protamine has a biphasic effect on the respiration of intact mitochondria. At lower protamine concentrations respiration is stimulated and this combined with a decrease of the respiratory control index; at higher ones respiration is inhibited and respiratory control is lost. In cytochrome c-depleted and restored mitochondria protamine effect on oxidative phosphorylation is only inhibitory. Increasing cytochrome c concentrations restore respiration in protamine-treated cytochrome c depleted mitochondria but not the respiratory control. Binding of cytochrome c to mitochondria is studied by determining from Scatchard plots the number of high affinity binding sites (n) and their stability constants (K). In absence of protamine in intact mitochondria n = 2.7 and K = 4.67-10(6) M-1; in cotochrome c depleted mitochondria n = 4.7 and K = 5.16-10(6) M-1. In both types of mitochondria protamine decreases significantly n as well as K. These data show that protamine may affect oxidative phosphorylation by causing desorption of cytochrome c from the inner mitochondrial membrane.     

Manipulating respiratory levels in Escherichia coli for aerobic formation of reduced chemical products

Optimizing the productivity of bioengineered strains requires balancing ATP generation and carbon atom conservation through fine-tuning cell respiration and metabolism. Traditional approaches manipulate cell respiration by altering air feeding, which are technically difficult especially in large bioreactors. An approach based on genetic regulation may better serve this purpose. With excess oxygen supply to the culture, we efficiently manipulated Escherichia coli cell respiration by adding different amount of coenzyme Q1 to strains lacking the ubiCA genes, which encode two critical enzymes for ubiquinone synthesis. As a proof-of-concept, the metabolic effect of the ubiCA gene knockout and coenzyme Q1 supplementation were characterized, and the metabolic profiles of the experimental strains showed clear correlations with coenzyme Q1 concentrations. Further proof-of-principle experiments were performed to illustrate that the approach can be used to optimize cell respiration for the production of chemicals of interest such as ethanol. This study showed that controlled respiration through genetic manipulation can be exploited to allow much larger operating windows for reduced product formation even under fully aerobic conditions.     

Recognition and analysis of protein-coding genes in severe acute respiratory syndrome associated coronavirus

MOTIVATION: The recent outbreak of severe acute respiratory syndrome (SARS) caused by SARS coronavirus (SARS-CoV) has necessitated an in-depth molecular understanding of the virus to identify new drug targets. The availability of complete genome sequence of several strains of SARS virus provides the possibility of identification of protein-coding genes and defining their functions. Computational approach to identify protein-coding genes and their putative functions will help in designing experimental protocols. RESULTS: In this paper, a novel analysis of SARS genome using gene prediction method GeneDecipher developed in our laboratory has been presented. Each of the 18 newly sequenced SARS-CoV genomes has been analyzed using GeneDecipher. In addition to polyprotein 1ab(1), polyprotein 1a and the four genes coding for major structural proteins spike (S), small envelope (E), membrane (M) and nucleocapsid (N), six to eight additional proteins have been predicted depending upon the strain analyzed. Their lengths range between 61 and 274 amino acids. Our method also suggests that polyprotein 1ab, polyprotein 1a, S, M and N are proteins of viral origin and others are of prokaryotic. Putative functions of all predicted protein-coding genes have been suggested using conserved peptides present in their open reading frames. AVAILABILITY: Detailed results of GeneDecipher analysis of all the 18 strains of SARS-CoV genomes are available at http://www.igib.res.in/sarsanalysis.html