Accurate reconstruction of metabolic pathways is an important prerequisite for interpreting metabolomics changes and understanding the diverse biological processes in disease models. A tracer-based... Show moreAccurate reconstruction of metabolic pathways is an important prerequisite for interpreting metabolomics changes and understanding the diverse biological processes in disease models. A tracer-based metabolomics strategy utilizes stable isotope-labeled precursors to resolve complex pathways by tracing the labeled atom(s) to downstream metabolites through enzymatic reactions. Isotope enrichment analysis is informative and achieved by counting total labeled atoms and acquiring the mass isotopologue distribution (MID) of the intact metabolite. However, quantitative analysis of labeled metabolite substructures/moieties (MS2 fragments) can offer more valuable insights into the reaction connections through measuring metabolite transformation. In order to acquire the isotopic labeling information at the intact metabolite and moiety level simultaneously, we developed a method that couples hydrophilic interaction liquid chromatography (HILIC) with Zeno trap-enabled high-resolution multiple reaction monitoring (MRMHR). The method enabled accurate and reproducible MID quantification for intact metabolites as well as their fragmented moieties, with notably high sensitivity in the MS2 fragmentation mode based on the measurement of 13C- or 15N-labeled cellular samples. The method was applied to human-induced pluripotent stem cell-derived neurons to trace the fate of 13C/15N atoms from D-13C6-glucose/L-15N2-glutamine added to the media. With the MID analysis of both intact metabolites and fragmented moieties, we validated the pathway reconstruction of de novo glutathione synthesis in mid-brain neurons. We discovered increased glutathione oxidization from both basal and newly synthesized glutathione pools under neuronal oxidative stress. Furthermore, the significantly decreased de novo glutathione synthesis was investigated and associated with altered activities of several key enzymes, as evidenced by suppressed glutamate supply via glucose metabolism and a diminished flux of glutathione synthetic reaction in the neuronal model of rotenone-induced neurodegeneration. Show less
We aimed at investigating host-virus co-metabolism during SARS-CoV-2 infection. Therefore, we extended comprehensive sex-specific, whole-body organ resolved models of human metabolism with the... Show moreWe aimed at investigating host-virus co-metabolism during SARS-CoV-2 infection. Therefore, we extended comprehensive sex-specific, whole-body organ resolved models of human metabolism with the necessary reactions to replicate SARS-CoV-2 in the lung as well as selected peripheral organs. Using this comprehensive host-virus model, we obtained the following key results: 1. The predicted maximal possible virus shedding rate was limited by isoleucine availability. 2. The supported initial viral load depended on the increase in CD4+ T-cells, consistent with the literature. 3. During viral infection, the whole-body metabolism changed including the blood metabolome, which agreed well with metabolomic studies from COVID-19 patients and healthy controls. 4. The virus shedding rate could be reduced by either inhibition of the guanylate kinase 1 or availability of amino acids, e.g., in the diet. 5. The virus variants differed in their maximal possible virus shedding rates, which could be inversely linked to isoleucine occurrences in the sequences. Taken together, this study presents the metabolic crosstalk between host and virus and emphasises the role of amino acid metabolism during SARS-CoV-2 infection, in particular of isoleucine. As such, it provides an example of how computational modelling can complement more canonical approaches to gain insight into host-virus crosstalk and to identify potential therapeutic strategies. Show less
MOTIVATION\nRESULTS\nAVAILABILITY AND IMPLEMENTATION\nSUPPLEMENTARY INFORMATION\nGenome-scale metabolic reconstructions have been assembled for thousands of organisms using a wide range of tools.... Show moreMOTIVATION\nRESULTS\nAVAILABILITY AND IMPLEMENTATION\nSUPPLEMENTARY INFORMATION\nGenome-scale metabolic reconstructions have been assembled for thousands of organisms using a wide range of tools. However, metabolite annotations, required to compare and link metabolites between reconstructions, remain incomplete. Here, we aim to further extend metabolite annotation coverage using various databases and chemoinformatic approaches.\nWe developed a COBRA toolbox extension, deemed MetaboAnnotator, which facilitates the comprehensive annotation of metabolites with database independent and dependent identifiers, obtains molecular structure files, and calculates metabolite formula and charge at pH 7.2. The resulting metabolite annotations allow for subsequent cross-mapping between reconstructions and mapping of, e.g., metabolomic data.\nMetaboAnnotator and tutorials are freely available at https://github.com/opencobra.\nSupplementary data are available at Bioinformatics online. Show less
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... Show moreComprehensive 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. Show less
In constraint-based metabolic modelling, physical and biochemical constraints define a polyhedral convex set of feasible flux vectors. Uniform sampling of this set provides an unbiased... Show moreIn constraint-based metabolic modelling, physical and biochemical constraints define a polyhedral convex set of feasible flux vectors. Uniform sampling of this set provides an unbiased characterization of the metabolic capabilities of a biochemical network. However, reliable uniform sampling of genome-scale biochemical networks is challenging due to their high dimensionality and inherent anisotropy. Here, we present an implementation of a new sampling algorithm, coordinate hit-and-run with rounding (CHRR). This algorithm is based on the provably efficient hit-and-run random walk and crucially uses a preprocessing step to round the anisotropic flux set. CHRR provably converges to a uniform stationary sampling distribution. We apply it to metabolic networks of increasing dimensionality. We show that it converges several times faster than a popular artificial centering hit-and-run algorithm, enabling reliable and tractable sampling of genome-scale biochemical networks. Show less