A reactive species reactions module for integration into genome-scale metabolic models for improved insights: Application to cancer.

Reactive species (RS) play significant roles in many disease contexts. Despite their crucial roles in diseases including cancer, the RS are not adequately modeled in the genome-scale metabolic (GSM) models, which are used to understand cell metabolism in disease contexts. We have developed a scalable RS reactions module that can be integrated with any Recon 3D-derived human metabolic model, or after fine-tuning, with any metabolic model. With RS-integration, the GSM models of three cancers (basal-like triple negative breast cancer (TNBC), high grade serous ovarian carcinoma (HGSOC) and colorectal cancer (CRC)) built from Recon 3D, precisely highlighted the increases/decreases in fluxes (dysregulation) occurring in important pathways of these cancers. These dysregulations were not prominent in the standard cancer models without the RS module. Further, the results from these RS-integrated cancer GSM models suggest the following decreasing order in the ease of ferroptosis-targeting to treat the cancers: TNBC > HGSOC > CRC.

A Computational Framework for Studying Gut-Brain Axis in Autism Spectrum Disorder.

Introduction: The integrity of the intestinal epithelium is crucial for human health and is harmed in autism spectrum disorder (ASD). An aberrant gut microbial composition resulting in gut-derived metabolic toxins was found to damage the intestinal epithelium, jeopardizing tissue integrity. These toxins further reach the brain via the gut-brain axis, disrupting the normal function of the brain. A mechanistic understanding of metabolic disturbances in the brain and gut is essential to design effective therapeutics and early intervention to block disease progression. Herein, we present a novel computational framework integrating constraint based tissue specific metabolic (CBM) model and whole-body physiological pharmacokinetics (PBPK) modeling for ASD. Furthermore, the role of gut microbiota, diet, and oxidative stress is analyzed in ASD. Methods: A representative gut model capturing host-bacteria and bacteria-bacteria interaction was developed using CBM techniques and patient data. Simultaneously, a PBPK model of toxin metabolism was assembled, incorporating multi-scale metabolic information. Furthermore, dynamic flux balance analysis was performed to integrate CBM and PBPK. The effectiveness of a probiotic and dietary intervention to improve autism symptoms was tested on the integrated model. Results: The model accurately highlighted critical metabolic pathways of the gut and brain that are associated with ASD. These include central carbon, nucleotide, and vitamin metabolism in the host gut, and mitochondrial energy and amino acid metabolisms in the brain. The proposed dietary intervention revealed that a high-fiber diet is more effective than a western diet in reducing toxins produced inside the gut. The addition of probiotic bacteria Lactobacillus acidophilus, Bifidobacterium longum longum, Akkermansia muciniphila, and Prevotella ruminicola to the diet restores gut microbiota balance, thereby lowering oxidative stress in the gut and brain. Conclusion: The proposed computational framework is novel in its applicability, as demonstrated by the determination of the whole-body distribution of ROS toxins and metabolic association in ASD. In addition, it emphasized the potential for developing novel therapeutic strategies to alleviate autism symptoms. Notably, the presented integrated model validates the importance of combining PBPK modeling with COBRA -specific tissue details for understanding disease pathogenesis.

Metabolic modeling of host-microbe interactions for therapeutics in colorectal cancer.

The onset of colorectal cancer (CRC) is often attributed to gut bacterial dysbiosis, and thus gut microbiota are highly relevant in devising treatment strategies. Certain gut microbes, like Enterococcus spp., exhibit remarkable anti-neoplastic and probiotic properties, which can aid in silver nanoparticle (AgNPs) induced reactive oxygen species (ROS)-based CRC treatment. However, the effects of AgNPs on gut microbial metabolism have not been reported thus far. In this study, a detailed systems-level understanding of ROS metabolism in Enterococcus durans (E. durans), a representative gut microbe, was gained using constraint-based modeling, wherein, the critical association between ROS and folate metabolism was established. Experimental studies involving low AgNP concentration treatment of E. durans cultures confirmed these modeling predictions (an increased extracellular folate concentration by 52%, at the 9th h of microbial growth, was observed). Besides, the computational studies established various metabolic pathways involving amino acids, energy metabolites, nucleotides, and SCFAs as the key players in elevating folate levels on ROS exposure. The anti-cancer potential of E. durans was also studied through MTT analysis of HCT 116 cells treated with microbial culture (AgNP treated) supernatant. A decrease in cell viability by 19% implicated the role of microbial metabolites (primarily folate) in causing cell death. The genome-scale modeling approach was then extended to extensively model CRC metabolism, as well as CRC-E. durans interactions in the context of CRC treatment, using tissue-specific metabolic models of CRC and healthy colon. These findings on further validation can facilitate the development of robust and effective cancer therapy.

Personalized whole-body models integrate metabolism, physiology, and the gut microbiome.

Comprehensive molecular-level models of human metabolism have been generated on a cellular level. However, models of whole-body metabolism have not been established as they require new methodological approaches to integrate molecular and physiological data. We developed a new metabolic network reconstruction approach that used organ-specific information from literature and omics data to generate two sex-specific whole-body metabolic (WBM) reconstructions. These reconstructions capture the metabolism of 26 organs and six blood cell types. Each WBM reconstruction represents whole-body organ-resolved metabolism with over 80,000 biochemical reactions in an anatomically and physiologically consistent manner. We parameterized the WBM reconstructions with physiological, dietary, and metabolomic data. The resulting WBM models could recapitulate known inter-organ metabolic cycles and energy use. We also illustrate that the WBM models can predict known biomarkers of inherited metabolic diseases in different biofluids. Predictions of basal metabolic rates, by WBM models personalized with physiological data, outperformed current phenomenological models. Finally, integrating microbiome data allowed the exploration of host-microbiome co-metabolism. Overall, the WBM reconstructions, and their derived computational models, represent an important step toward virtual physiological humans.

Mechanism of the oxidative stress-mediated increase in lipid accumulation by the bacterium, R. opacus PD630: Experimental analysis and genome-scale metabolic modeling.

Appropriate species of oleaginous bacteria, with their high growth rates and lipid accumulation capabilities, can be good contenders for industrial triacylglycerol (TAG) production, compared to microalgae. Further, oxidative stress (OS) can be used to significantly increase TAG yields in oleaginous microbes, but the mechanism is unexplored. In a first, this study explored the mechanism behind OS-mediated increase in TAG accumulation by the bacterium, Rhodococccus opacus PD630, through experimental analysis and metabolic modelling. Two mechanisms that could increase acetyl-CoA (TAG-precursor) levels were hypothesized based on literature information. One was OS-mediated inactivation of the aconitase (TCA cycle), and another was the inactivation of the triosephosphate isomerase (TPI; glycolysis). The results negated the involvement of aconitase in increased acetyl-CoA levels. Analysis of the metabolic model showed that inactivation of TPI, re-routed the flux through the pentose phosphate pathway (PPP), supplying both NADPH and acetyl-CoA for TAG synthesis. Additionally, inactivation of TPI increased TAG flux by 143%, whereas, inactivating both TPI and aconitase, increased it by 152%. We present experimental evidence for OS-mediated decrease in TPI activity and increase in activity of glucose-6-phosphate dehydrogenase (PPP enzyme). The findings indicate that increased flux through PPP can be explored to improve TAG accumulation on a large-scale.

Metabolite systems profiling identifies exploitable weaknesses in retinoblastoma.

Retinoblastoma (RB) is a childhood eye cancer. Currently, chemotherapy, local therapy, and enucleation are the main ways in which these tumors are managed. The present work is the first study that uses constraint-based reconstruction and analysis approaches to identify and explain RB-specific survival strategies, which are RB tumor specific. Importantly, our model-specific secretion profile is also found in RB1-depleted human retinal cells in vitro and suggests that novel biomarkers involved in lipid metabolism may be important. Finally, RB-specific synthetic lethals have been predicted as lipid and nucleoside transport proteins that can aid in novel drug target development.

The Virtual Metabolic Human database: integrating human and gut microbiome metabolism with nutrition and disease.

A multitude of factors contribute to complex diseases and can be measured with 'omics' methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, www.vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources 'Human metabolism', 'Gut microbiome', 'Disease', 'Nutrition', and 'ReconMaps'. The VMH captures 5180 unique metabolites, 17 730 unique reactions, 3695 human genes, 255 Mendelian diseases, 818 microbes, 632 685 microbial genes and 8790 food items. The VMH's unique features are (i) the hosting of the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; (ii) seven human metabolic maps for data visualization; (iii) a nutrition designer; (iv) a user-friendly webpage and application-programming interface to access its content; (v) user feedback option for community engagement and (vi) the connection of its entities to 57 other web resources. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.

Recon3D enables a three-dimensional view of gene variation in human metabolism.

Genome-scale network reconstructions have helped uncover the molecular basis of metabolism. Here we present Recon3D, a computational resource that includes three-dimensional (3D) metabolite and protein structure data and enables integrated analyses of metabolic functions in humans. We use Recon3D to functionally characterize mutations associated with disease, and identify metabolic response signatures that are caused by exposure to certain drugs. Recon3D represents the most comprehensive human metabolic network model to date, accounting for 3,288 open reading frames (representing 17% of functionally annotated human genes), 13,543 metabolic reactions involving 4,140 unique metabolites, and 12,890 protein structures. These data provide a unique resource for investigating molecular mechanisms of human metabolism. Recon3D is available at http://vmh.life.

Modeling the effects of commonly used drugs on human metabolism.

Metabolism contributes significantly to the pharmacokinetics and pharmacodynamics of a drug. In addition, diet and genetics have a profound effect on cellular metabolism with respect to both health and disease. In the present study, we assembled a comprehensive, literature-based drug metabolic reconstruction of the 18 most highly prescribed drug groups, including statins, anti-hypertensives, immunosuppressants and analgesics. This reconstruction captures in detail our current understanding of their absorption, intracellular distribution, metabolism and elimination. We combined this drug module with the most comprehensive reconstruction of human metabolism, Recon 2, yielding Recon2_DM1796, which accounts for 2803 metabolites and 8161 reactions. By defining 50 specific drug objectives that captured the overall drug metabolism of these compounds, we investigated the effects of dietary composition and inherited metabolic disorders on drug metabolism and drug-drug interactions. Our main findings include: (a) a shift in dietary patterns significantly affects statins and acetaminophen metabolism; (b) disturbed statin metabolism contributes to the clinical phenotype of mitochondrial energy disorders; and (c) the interaction between statins and cyclosporine can be explained by several common metabolic and transport pathways other than the previously established CYP3A4 connection. This work holds the potential for studying adverse drug reactions and designing patient-specific therapies.

Membrane transporters in a human genome-scale metabolic knowledgebase and their implications for disease.

Membrane transporters enable efficient cellular metabolism, aid in nutrient sensing, and have been associated with various diseases, such as obesity and cancer. Genome-scale metabolic network reconstructions capture genomic, physiological, and biochemical knowledge of a target organism, along with a detailed representation of the cellular metabolite transport mechanisms. Since the first reconstruction of human metabolism, Recon 1, published in 2007, progress has been made in the field of metabolite transport. Recently, we published an updated reconstruction, Recon 2, which significantly improved the metabolic coverage and functionality. Human metabolic reconstructions have been used to investigate the role of metabolism in disease and to predict biomarkers and drug targets. Given the importance of cellular transport systems in understanding human metabolism in health and disease, we analyzed the coverage of transport systems for various metabolite classes in Recon 2. We will review the current knowledge on transporters (i.e., their preferred substrates, transport mechanisms, metabolic relevance, and disease association for each metabolite class). We will assess missing coverage and propose modifications and additions through a transport module that is functional when combined with Recon 2. This information will be valuable for further refinements. These data will also provide starting points for further experiments by highlighting areas of incomplete knowledge. This review represents the first comprehensive overview of the transporters involved in central metabolism and their transport mechanisms, thus serving as a compendium of metabolite transporters specific for human metabolic reconstructions.

A community-driven global reconstruction of human metabolism.

Multiple models of human metabolism have been reconstructed, but each represents only a subset of our knowledge. Here we describe Recon 2, a community-driven, consensus 'metabolic reconstruction', which is the most comprehensive representation of human metabolism that is applicable to computational modeling. Compared with its predecessors, the reconstruction has improved topological and functional features, including ∼2× more reactions and ∼1.7× more unique metabolites. Using Recon 2 we predicted changes in metabolite biomarkers for 49 inborn errors of metabolism with 77% accuracy when compared to experimental data. Mapping metabolomic data and drug information onto Recon 2 demonstrates its potential for integrating and analyzing diverse data types. Using protein expression data, we automatically generated a compendium of 65 cell type-specific models, providing a basis for manual curation or investigation of cell-specific metabolic properties. Recon 2 will facilitate many future biomedical studies and is freely available at http://humanmetabolism.org/.

Predicting the impact of diet and enzymopathies on human small intestinal epithelial cells.

Small intestinal epithelial cells (sIECs) have a significant share in whole body metabolism as they perform enzymatic digestion and absorption of nutrients. Furthermore, the diet plays a key role in a number of complex diseases including obesity and diabetes. The impact of diet and altered genetic backgrounds on human metabolism may be studied by using computational modeling. A metabolic reconstruction of human sIECs was manually assembled using the literature. The resulting sIEC model was subjected to two different diets to obtain condition-specific metabolic models. Fifty defined metabolic tasks evaluated the functionalities of these models, along with the respective secretion profiles, which distinguished between impacts of different dietary regimes. Under the average American diet, the sIEC model resulted in higher secretion flux for metabolites implicated in metabolic syndrome. In addition, enzymopathies were analyzed in the context of the sIEC metabolism. Computed results were compared with reported gastrointestinal (GI) pathologies and biochemical defects as well as with biomarker patterns used in their diagnosis. Based on our simulations, we propose that (i) sIEC metabolism is perturbed by numerous enzymopathies, which can be used to study cellular adaptive mechanisms specific for such disorders, and in the identification of novel co-morbidities, (ii) porphyrias are associated with both heme synthesis and degradation and (iii) disturbed intestinal gamma-aminobutyric acid synthesis may be linked to neurological manifestations of various enzymopathies. Taken together, the sIEC model represents a comprehensive, biochemically accurate platform for studying the function of sIEC and their role in whole body metabolism.

Systems-level characterization of a host-microbe metabolic symbiosis in the mammalian gut.

The human gut microbiota consists of ten times more microorganisms than there are cells in our body, processes otherwise indigestible nutrients, and produces important energy precursors, essential amino acids, and vitamins. In this study, we assembled and validated a genome-scale metabolic reconstruction of Bacteroides thetaiotaomicron (iAH991), a prominent representative of the human gut microbiota, consisting of 1488 reactions, 1152 metabolites, and 991 genes. To create a comprehensive metabolic model of host-microbe interactions, we integrated iAH991 with a previously published mouse metabolic reconstruction, which was extended for intestinal transport and absorption reactions. The two metabolic models were linked through a joint compartment, the lumen, allowing metabolite exchange and providing a route for simulating different dietary regimes. The resulting model consists of 7239 reactions, 5164 metabolites, and 2769 genes. We simultaneously modeled growth of mouse and B. thetaiotaomicron on five different diets varying in fat, carbohydrate, and protein content. The integrated model captured mutually beneficial cross-feeding as well as competitive interactions. Furthermore, we identified metabolites that were exchanged between the two organisms, which were compared with published metabolomics data. This analysis resulted for the first time in a comprehensive description of the co-metabolism between a host and its commensal microbe. We also demonstrate in silico that the presence of B. thetaiotaomicron could rescue the growth phenotype of the host with an otherwise lethal enzymopathy and vice versa. This systems approach represents a powerful tool for modeling metabolic interactions between a gut microbe and its host in health and disease.

A compendium of inborn errors of metabolism mapped onto the human metabolic network.

Inborn errors of metabolism (IEMs) are hereditary metabolic defects, which are encountered in almost all major metabolic pathways occurring in man. Many IEMs are screened for in neonates through metabolomic analysis of dried blood spot samples. To enable the mapping of these metabolomic data onto the published human metabolic reconstruction, we added missing reactions and pathways involved in acylcarnitine (AC) and fatty acid oxidation (FAO) metabolism. Using literary data, we reconstructed an AC/FAO module consisting of 352 reactions and 139 metabolites. When this module was combined with the human metabolic reconstruction, the synthesis of 39 acylcarnitines and 22 amino acids, which are routinely measured, was captured and 235 distinct IEMs could be mapped. We collected phenotypic and clinical features for each IEM enabling comprehensive classification. We found that carbohydrate, amino acid, and lipid metabolism were most affected by the IEMs, while the brain was the most commonly affected organ. Furthermore, we analyzed the IEMs in the context of metabolic network topology to gain insight into common features between metabolically connected IEMs. While many known examples were identified, we discovered some surprising IEM pairs that shared reactions as well as clinical features but not necessarily causal genes. Moreover, we could also re-confirm that acetyl-CoA acts as a central metabolite. This network based analysis leads to further insight of hot spots in human metabolism with respect to IEMs. The presented comprehensive knowledge base of IEMs will provide a valuable tool in studying metabolic changes involved in inherited metabolic diseases.