Towards kinetic modeling of genome-scale metabolic networks without sacrificing stoichiometric,thermodynamic and physiological constraints

Mathematical modeling is an essential tool for the comprehensive understanding of cell metabolism and its interactions with the environmental and process conditions. Recent developments in the construction and analysis of stoichiometric models made it possible to define limits on steady-state metabolic behavior using flux balance analysis. However, detailed information on enzyme kinetics and enzyme regulation is needed to formulate kinetic models that can accurately capture the dynamic metabolic responses. The use of mechanistic enzyme kinetics is a difficult task due to uncertainty in the kinetic properties of enzymes. Therefore, the majority of recent works considered only mass action kinetics for reactions in metabolic networks. Herein, we applied the optimization and risk analysis of complex living entities (ORACLE) framework and constructed a large-scale mechanistic kinetic model of optimally grown Escherichia coli. We investigated the complex interplay between stoichiometry, thermodynamics, and kinetics in determining the flexibility and capabilities of metabolism. Our results indicate that enzyme saturation is a necessary consideration in modeling metabolic networks and it extends the feasible ranges of metabolic fluxes and metabolite concentrations. Our results further suggest that enzymes in metabolic networks have evolved to function at different saturation states to ensure greater flexibility and robustness of cellular metabolism.     

Dynamics and Control of the Central Carbon Metabolism in Hepatoma Cells

Background  

The liver plays a major role in metabolism and performs a number of vital functions in the body. Therefore, the determination of hepatic metabolite dynamics and the analysis of the control of the respective biochemical pathways are of great pharmacological and medical importance. Extra- and intracellular time-series data from stimulus-response experiments are gaining in importance in the identification of in vivo metabolite dynamics, while dynamic network models are excellent tools for analyzing complex metabolic control patterns. This is the first study that has been undertaken on the data-driven identification of a dynamic liver central carbon metabolism model and its application in the analysis of the distribution of metabolic control in hepatoma cells.     

Hydrogen peroxide production and myo‐inositol metabolism as important traits for virulence of Mycoplasma hyopneumoniae

Mycoplasma hyopneumoniae is the causative agent of enzootic pneumonia. In our previous work, we reconstructed the metabolic models of this species along with two other mycoplasmas from the respiratory tract of swine: Mycoplasma hyorhinis, considered less pathogenic but which nonetheless causes disease and Mycoplasma flocculare, a commensal bacterium. We identified metabolic differences that partially explained their different levels of pathogenicity. One important trait was the production of hydrogen peroxide from the glycerol metabolism only in the pathogenic species. Another important feature was a pathway for the metabolism of myo‐inositol in M. hyopneumoniae. Here, we tested these traits to understand their relation to the different levels of pathogenicity, comparing not only the species but also pathogenic and attenuated strains of M. hyopneumoniae. Regarding the myo‐inositol metabolism, we show that only M. hyopneumoniae assimilated this carbohydrate and remained viable when myo‐inositol was the primary energy source. Strikingly, only the two pathogenic strains of M. hyopneumoniae produced hydrogen peroxide in complex medium. We also show that this production was dependent on the presence of glycerol. Although further functional tests are needed, we present in this work two interesting metabolic traits of M. hyopneumoniae that might be directly related to its enhanced virulence.     

Effects of Acetyl-l-Carnitine and Myo-Inositol on High-Energy Phosphate and Membrane Phospholipid Metabolism in Zebra Fish: A 31P-NMR-Spectroscopy Study

Acetyl-l-carnitine (ALCAR) and myo-inositol are reported to enhance motor activity in animal models; modulate membrane phospholipid metabolism (ALCAR and myo-inositol) and high-energy phosphate metabolism (ALCAR) back to normal; and be effective treatments of major depression in humans. Fish in general and zebra fish in particular present unique animal models for the in vivo study of high-energy phosphate and membrane phospholipid metabolism by noninvasive in vivo ^31P NMR. This ^31P NMR study of free-swimming zebra fish showed that both ALCAR and myo-inositol decreased levels of phosphodiesters and inorganic orthophosphate and increased levels of PCr in the fish. These findings demonstrate both ALCAR and myo-inositol modulate membrane phospholipid and high-energy phosphate metabolism in free-swimming zebra fish.     

Statistical epistasis between candidate gene alleles for complex tuber traits in an association mapping population of tetraploid potato

Association mapping using DNA-based markers is a novel tool in plant genetics for the analysis of complex traits. Potato tuber yield, starch content, starch yield and chip color are complex traits of agronomic relevance, for which carbohydrate metabolism plays an important role. At the functional level, the genes and biochemical pathways involved in carbohydrate metabolism are among the best studied in plants. Quantitative traits such as tuber starch and sugar content are therefore models for association genetics in potato based on candidate genes. In an association mapping experiment conducted with a population of 243 tetraploid potato varieties and breeding clones, we previously identified associations between individual candidate gene alleles and tuber starch content, starch yield and chip quality. In the present paper, we tested 190 DNA markers at 36 loci scored in the same association mapping population for pairwise statistical epistatic interactions. Fifty marker pairs were associated mainly with tuber starch content and/or starch yield, at a cut-off value of q ≤ 0.20 for the experiment-wide false discovery rate (FDR). Thirteen marker pairs had an FDR of q ≤ 0.10. Alleles at loci encoding ribulose-bisphosphate carboxylase/oxygenase activase (Rca), sucrose phosphate synthase (Sps) and vacuolar invertase (Pain1) were most frequently involved in statistical epistatic interactions. The largest effect on tuber starch content and starch yield was observed for the paired alleles Pain1-8c and Rca-1a, explaining 9 and 10% of the total variance, respectively. The combination of these two alleles increased the means of tuber starch content and starch yield. Biological models to explain the observed statistical epistatic interactions are discussed.     

Biochemical characteristization of propionyl-coenzyme a carboxylase complex of <Emphasis Type="Italic">Streptomyces toxytricini</Emphasis>

Acyl-CoA carboxylases (ACC) are involved in important primary or secondary metabolic pathways such as fatty acid and/or polyketides synthesis. In the 62 kb fragment of pccB gene locus of Streptomyces toxytricini producing a pancreatic inhibitor lipstatin, 3 distinct subunit genes of presumable propionyl-CoA carboxylase (PCCase) complex, assumed to be one of ACC responsible for the secondary metabolism, were identified along with gene for a biotin protein ligase (Bpl). The subunits of PCCase complex were a subunit (AccA3), P subunit (PccB), and auxiliary ɛ subunit (PccE). In order to disclose the involvement of the PCCase complex in secondary metabolism, some biochemical characteristics of each subunit as well as their complex were examined. In the test of substrate specificity of the PCCase complex, it was confirmed that this complex showed much higher conversion of propionyl-CoA rather than acetyl-CoA. It implies the enzyme complex could play a main role in the production of methylmalonyl-CoA from propionyl-CoA, which is a precursor of secondary polyketide biosynthesis.     

Physiological functions of pyruvate:NADP+ oxidoreductase and 2-oxoglutarate decarboxylase in Euglena gracilis under aerobic and anaerobic conditions

In Euglena gracilis, pyruvate:NADP⁺ oxidoreductase, in addition to the pyruvate dehydrogenase complex, functions for the oxidative decarboxylation of pyruvate in the mitochondria. Furthermore, the 2-oxoglutarate dehydrogenase complex is absent, and instead 2-oxoglutarate decarboxylase is found in the mitochondria. To elucidate the central carbon and energy metabolisms in Euglena under aerobic and anaerobic conditions, physiological significances of these enzymes involved in 2-oxoacid metabolism were examined by gene silencing experiments. The pyruvate dehydrogenase complex was indispensable for aerobic cell growth in a glucose medium, although its activity was less than 1% of that of pyruvate:NADP⁺ oxidoreductase. In contrast, pyruvate:NADP⁺ oxidoreductase was only involved in the anaerobic energy metabolism (wax ester fermentation). Aerobic cell growth was almost completely suppressed when the 2-oxoglutarate decarboxylase gene was silenced, suggesting that the tricarboxylic acid cycle is modified in Euglena and 2-oxoglutarate decarboxylase takes the place of the 2-oxoglutarate dehydrogenase complex in the aerobic respiratory metabolism.     

New ENU-induced semidominant mutation, Ali18, causes inflammatory arthritis, dermatitis, and osteoporosis in the mouse

Inflammation is a complex cellular and humoral response against trauma and infection, and its presence leads to destruction of tissue in humans. The mechanisms that initiate inflammatory diseases remain largely unknown because of complex interactions between multiple genetic and environmental factors during pathogenesis. Animal models for human diseases offer dissection of complex pathogenesis by inbred genetic backgrounds and controlled circumstances. In this article we report a chemically induced new mutation, Ali18 (Abnormal limb), as a mouse model for inflammatory arthritis and dermatitis. Ali18/+ mice exhibit rubor and swelling of footpads in hindlimbs in adults. In Ali18/Ali18 mice, the digits in forelimbs and hindlimbs and tails were necrotic and/or deformed by severe swelling. Histologic analysis revealed infiltration of mixed populations of inflammatory cells into bone marrow, peripheral joints, and skin in the affected areas of Ali18/Ali18 mice. In addition, generalized osteoporosis-like phenotypes were confirmed by dual energy X-ray absorptiometry (DXA), microcomputed tomography (μCT), and peripheral quantitative computed tomography (pQCT) in homozygous animals. Whereas the Ali18 mutation was mapped to a single locus, the phenotype presentation was altered by complex modifier effects from other inbred genetic backgrounds. Detailed analysis of the Ali18 phenotype and identification of the mutation and its modifier genes may provide molecular insights into the complex nature of inflammatory diseases and the relationship between inflammation and bone metabolism.     

Effects of calorie restriction on life span of microorganisms

Calorie restriction (CR) in microorganisms such as budding and fission yeasts has a robust and well-documented impact on longevity. In order to efficiently utilize the limited energy during CR, these organisms shift from primarily fermentative metabolism to mitochondrial respiration. Respiration activates certain conserved longevity factors such as sirtuins and is associated with widespread physiological changes that contribute to increased survival. However, the importance of respiration during CR-mediated longevity has remained controversial. The emergence of several novel metabolically distinct microbial models for longevity has enabled CR to be studied from new perspectives. The majority of CR and life span studies have been conducted in the primarily fermentative Crabtree-positive yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe, but studies in primarily respiratory Crabtree-negative yeast and obligate aerobes can offer complementary insight into the more complex mammalian response to CR. Not only are microorganisms helping characterize a conserved cellular mechanism for CR-mediated longevity, but they can also directly impact mammalian metabolism as part of the natural gut flora. Here, we discuss the contributions of microorganisms to our knowledge of CR and longevity at the level of both the cell and the organism.     

Intracellular killing of bacteria: is Dictyostelium a model macrophage or an alien?

Predation of bacteria by phagocytic cells was first developed during evolution by environmental amoebae. Many of the core mechanisms used by amoebae to sense, ingest and kill bacteria have also been conserved in specialized phagocytic cells in mammalian organisms. Here we focus on recent results revealing how Dictyostelium discoideum senses and kills non‐pathogenic bacteria. In this model, genetic analysis of intracellular killing of bacteria has revealed a surprisingly complex array of specialized mechanisms. These results raise new questions on these processes, and challenge current models based largely on studies in mammalian phagocytes. In addition, recent studies suggest one additional level on complexity by revealing how Dictyostelium recognizes specifically various bacterial species and strains, and adapts its metabolism to process them. It remains to be seen to what extent mechanisms uncovered in Dictyostelium are also used in mammalian phagocytic cells.     

Modelling nonlinear dynamics of Crassulacean acid metabolism productivity and water use for global predictions

Crassulacean acid metabolism (CAM) crops are important agricultural commodities in water‐limited environments across the globe, yet modelling of CAM productivity lacks the sophistication of widely used C3 and C4 crop models, in part due to the complex responses of the CAM cycle to environmental conditions. This work builds on recent advances in CAM modelling to provide a framework for estimating CAM biomass yield and water use efficiency from basic principles. These advances, which integrate the CAM circadian rhythm with established models of carbon fixation, stomatal conductance and the soil–plant‐atmosphere continuum, are coupled to models of light attenuation, plant respiration and biomass partitioning. Resulting biomass yield and transpiration for Opuntia ficus‐indica and Agave tequilana are validated against field data and compared with predictions of CAM productivity obtained using the empirically based environmental productivity index. By representing regulation of the circadian state as a nonlinear oscillator, the modelling approach captures the diurnal dynamics of CAM stomatal conductance, allowing the prediction of CAM transpiration and water use efficiency for the first time at the plot scale. This approach may improve estimates of CAM productivity under light‐limiting conditions when compared with previous methods.     

Network analysis of metabolic enzyme evolution in <Emphasis Type="Italic">Escherichia coli</Emphasis>

Background  

The two most common models for the evolution of metabolism are the patchwork evolution model, where enzymes are thought to diverge from broad to narrow substrate specificity, and the retrograde evolution model, according to which enzymes evolve in response to substrate depletion. Analysis of the distribution of homologous enzyme pairs in the metabolic network can shed light on the respective importance of the two models. We here investigate the evolution of the metabolism in E. coli viewed as a single network using EcoCyc.     

Modeling and identification of endocrine-metabolic systems. theoretical aspects and their importance in practice

Some guidelines for modeling, identification, and validation in endocrinology and metabolism are presented. Fundamentals based on recent theoretical developments are given, and their essential role for solving physiological and clinical problems is shown. Concepts and issues discussed include a priori and a posteriori identifiability; optimal experiment design; assessment of model validity; simple and complex models; compartmental and noncompartmental models. Recent studies on the glucose-insulin system and on ketone body metabolism are used to illustrate the methodological line of thought discussed in the paper.     

Impact of MtrA on phosphate metabolism genes and the response to altered phosphate conditions in Streptomyces

Phosphate metabolism is known to be regulated by the PhoPR regulatory system in Streptomyces and some other bacteria. In this study, we report that MtrA also regulates phosphate metabolism in Streptomyces. Our data showed that, in Streptomyces coelicolor, MtrA regulates not only phosphate metabolism genes such as phoA but also phoP under different phosphate conditions, including growth on rich complex media without added inorganic phosphate and on defined media with low or high concentrations of inorganic phosphate. Cross-regulation was also observed among mtrA, phoP and glnR under these conditions. We demonstrated both in vitro and in vivo binding of MtrA to the promoter regions of genes associated with phosphate metabolism and to the intergenic region between phoR and phoU, indicating that these phosphate metabolism genes are targets of MtrA. We further showed that MtrA in S. lividans and S. venezuelae has detectable regulatory effects on expression of phosphate metabolism genes. Additionally, the MtrA homologue from Corynebacterium glutamicum bound predicted MtrA sites of multiple phosphate metabolism genes, implying its potential for regulating phosphate metabolism in this species. Overall, our findings support MtrA as a major regulator for phosphate metabolism in Streptomyces and also potentially in other actinobacteria.     

The Mre11 complex and the metabolism of chromosome breaks: the importance of communicating and holding things together

The conserved Mre11 complex (Mre11, Rad50, and Nbs1) plays a role in each aspect of chromosome break metabolism. The complex acts as a break sensor and functions in the activation and propagation of signaling pathways that govern cell cycle checkpoint functions in response to DNA damage. In addition, the Mre11 complex influences recombinational DNA repair through promoting recombination between sister chromatids. The Mre11 complex is required for mammalian cell viability but hypomorphic mutants of Mre11 and Nbs1 have been identified in human genetic instability disorders. These hypomorphic mutations, as well as those identified in yeast, have provided a benchmark for establishing mouse models of Mre11 complex deficiency. In addition to consideration of Mre11 complex functions in human cells and yeast, this review will discuss the characterization of mouse models and insight gleaned from those models regarding the metabolism of chromosome breaks. The current picture of break metabolism supports a central role for the Mre11 complex at the interface of chromosome stability and the regulation of cell growth. Further genetic analysis of the Mre11 complex will be an invaluable tool for dissecting its function on an organismal level and determining its role in the prevention of malignancy.