iMM1865: A New Reconstruction of Mouse Genome-Scale Metabolic Model

Machine learning for the advancement of genome-scale metabolic modeling

Constraint-based modeling (CBM) has evolved as the core systems biology tool to map the interrelations between genotype, phenotype, and external environment. The recent advancement of high-throughput experimental approaches and multi-omics strategies has generated a plethora of new and precise information from wide-ranging biological domains. On the other hand, the continuously growing field of machine learning (ML) and its specialized branch of deep learning (DL) provide essential computational architectures for decoding complex and heterogeneous biological data. In recent years, both multi-omics and ML have assisted in the escalation of CBM. Condition-specific omics data, such as transcriptomics and proteomics, helped contextualize the model prediction while analyzing a particular phenotypic signature. At the same time, the advanced ML tools have eased the model reconstruction and analysis to increase the accuracy and prediction power. However, the development of these multi-disciplinary methodological frameworks mainly occurs independently, which limits the concatenation of biological knowledge from different domains. Hence, we have reviewed the potential of integrating multi-disciplinary tools and strategies from various fields, such as synthetic biology, CBM, omics, and ML, to explore the biochemical phenomenon beyond the conventional biological dogma. How the integrative knowledge of these intersected domains has improved bioengineering and biomedical applications has also been highlighted. We categorically explained the conventional genome-scale metabolic model (GEM) reconstruction tools and their improvement strategies through ML paradigms. Further, the crucial role of ML and DL in omics data restructuring for GEM development has also been briefly discussed. Finally, the case-study-based assessment of the state-of-the-art method for improving biomedical and metabolic engineering strategies has been elaborated. Therefore, this review demonstrates how integrating experimental and in silico strategies can help map the ever-expanding knowledge of biological systems driven by condition-specific cellular information. This multiview approach will elevate the application of ML-based CBM in the biomedical and bioengineering fields for the betterment of society and the environment.

Genome-scale metabolic models in translational medicine: the current status and potential of machine learning in improving the effectiveness of the models

The genome-scale metabolic model (GEM) has emerged as one of the leading modeling approaches for systems-level metabolic studies and has been widely explored for a broad range of organisms and applications. Owing to the development of genome sequencing technologies and available biochemical data, it is possible to reconstruct GEMs for model and non-model microorganisms as well as for multicellular organisms such as humans and animal models. GEMs will evolve in parallel with the availability of biological data, new mathematical modeling techniques and the development of automated GEM reconstruction tools. The use of high-quality, context-specific GEMs, a subset of the original GEM in which inactive reactions are removed while maintaining metabolic functions in the extracted model, for model organisms along with machine learning (ML) techniques could increase their applications and effectiveness in translational research in the near future. Here, we briefly review the current state of GEMs, discuss the potential contributions of ML approaches for more efficient and frequent application of these models in translational research, and explore the extension of GEMs to integrative cellular models.

Modeling Red Blood Cell Metabolism in the Omics Era

Red blood cells (RBCs) are abundant (more than 80% of the total cells in the human body), yet relatively simple, as they lack nuclei and organelles, including mitochondria. Since the earliest days of biochemistry, the accessibility of blood and RBCs made them an ideal matrix for the characterization of metabolism. Because of this, investigations into RBC metabolism are of extreme relevance for research and diagnostic purposes in scientific and clinical endeavors. The relative simplicity of RBCs has made them an eligible model for the development of reconstruction maps of eukaryotic cell metabolism since the early days of systems biology. Computational models hold the potential to deepen knowledge of RBC metabolism, but also and foremost to predict in silico RBC metabolic behaviors in response to environmental stimuli. Here, we review now classic concepts on RBC metabolism, prior work in systems biology of unicellular organisms, and how this work paved the way for the development of reconstruction models of RBC metabolism. Translationally, we discuss how the fields of metabolomics and systems biology have generated evidence to advance our understanding of the RBC storage lesion, a process of decline in storage quality that impacts over a hundred million blood units transfused every year.

Mouse Memory CD8 T Cell Subsets Defined by Tissue-Resident Memory Integrin Expression Exhibit Distinct Metabolic Profiles

Tissue-resident memory CD8 T cells (TRM) principally reside in peripheral nonlymphoid tissues, such as lung and skin, and confer protection against a variety of illnesses ranging from infections to cancers. The functions of different memory CD8 T cell subsets have been linked with distinct metabolic pathways and differ from other CD8 T cell subsets. For example, skin-derived memory T cells undergo fatty acid oxidation and oxidative phosphorylation to a greater degree than circulating memory and naive cells. Lung TRMs defined by the cell-surface expression of integrins exist as distinct subsets that differ in gene expression and function. We hypothesize that TRM subsets with different integrin profiles will use unique metabolic programs. To test this, differential expression and pathway analysis were conducted on RNA sequencing datasets from mouse lung TRMs yielding significant differences related to metabolism. Next, metabolic models were constructed, and the predictions were interrogated using functional metabolite uptake assays. The levels of oxidative phosphorylation, mitochondrial mass, and neutral lipids were measured. Furthermore, to investigate the potential relationships to TRM development, T cell differentiation studies were conducted in vitro with varying concentrations of metabolites. These demonstrated that lipid conditions impact T cell survival, and that glucose concentration impacts the expression of canonical TRM marker CD49a, with no effect on central memory-like T cell marker CCR7. In summary, it is demonstrated that mouse resident memory T cell subsets defined by integrin expression in the lung have unique metabolic profiles, and that nutrient abundance can alter differentiation.

COMO: a pipeline for multi-omics data integration in metabolic modeling and drug discovery

Identifying potential drug targets using metabolic modeling requires integrating multiple modeling methods and heterogeneous biological datasets, which can be challenging without efficient tools. We developed Constraint-based Optimization of Metabolic Objectives (COMO), a user-friendly pipeline that integrates multi-omics data processing, context-specific metabolic model development, simulations, drug databases and disease data to aid drug discovery. COMO can be installed as a Docker Image or with Conda and includes intuitive instructions within a Jupyter Lab environment. It provides a comprehensive solution for the integration of bulk and single-cell RNA-seq, microarrays and proteomics outputs to develop context-specific metabolic models. Using public databases, open-source solutions for model construction and a streamlined approach for predicting repurposable drugs, COMO enables researchers to investigate low-cost alternatives and novel disease treatments. As a case study, we used the pipeline to construct metabolic models of B cells, which simulate and analyze them to predict metabolic drug targets for rheumatoid arthritis and systemic lupus erythematosus, respectively. COMO can be used to construct models for any cell or tissue type and identify drugs for any human disease where metabolic inhibition is relevant. The pipeline has the potential to improve the health of the global community cost-effectively by providing high-confidence targets to pursue in preclinical and clinical studies. The source code of the COMO pipeline is available at https://github.com/HelikarLab/COMO. The Docker image can be pulled at https://github.com/HelikarLab/COMO/pkgs/container/como.

A new metabolic model of Drosophila melanogaster and the integrative analysis of Parkinson's disease

High conservation of the disease-associated genes between flies and humans facilitates the common use of Drosophila melanogaster to study metabolic disorders under controlled laboratory conditions. However, metabolic modeling studies are highly limited for this organism. We here report a comprehensively curated genome-scale metabolic network model of Drosophila using an orthology-based approach. The gene coverage and metabolic information of the draft model derived from a reference human model were expanded via Drosophila-specific KEGG and MetaCyc databases, with several curation steps to avoid metabolic redundancy and stoichiometric inconsistency. Furthermore, we performed literature-based curations to improve gene-reaction associations, subcellular metabolite locations, and various metabolic pathways. The performance of the resulting Drosophila model (8,230 reactions, 6,990 metabolites, and 2,388 genes), iDrosophila1 (https://github.com/SysBioGTU/iDrosophila), was assessed using flux balance analysis in comparison with the other currently available fly models leading to superior or comparable results. We also evaluated the transcriptome-based prediction capacity of iDrosophila1, where differential metabolic pathways during Parkinson's disease could be successfully elucidated. Overall, iDrosophila1 is promising to investigate system-level metabolic alterations in response to genetic and environmental perturbations.

Lifespan-extending interventions induce consistent patterns of fatty acid oxidation in mouse livers

Aging manifests as progressive deteriorations in homeostasis, requiring systems-level perspectives to investigate the gradual molecular dysregulation of underlying biological processes. Here, we report systemic changes in the molecular regulation of biological processes under multiple lifespan-extending interventions. Differential Rank Conservation (DIRAC) analyses of mouse liver proteomics and transcriptomics data show that mechanistically distinct lifespan-extending interventions (acarbose, 17α-estradiol, rapamycin, and calorie restriction) generally tighten the regulation of biological modules. These tightening patterns are similar across the interventions, particularly in processes such as fatty acid oxidation, immune response, and stress response. Differences in DIRAC patterns between proteins and transcripts highlight specific modules which may be tightened via augmented cap-independent translation. Moreover, the systemic shifts in fatty acid metabolism are supported through integrated analysis of liver transcriptomics data with a mouse genome-scale metabolic model. Our findings highlight the power of systems-level approaches for identifying and characterizing the biological processes involved in aging and longevity.

The NAD salvage pathway in mesenchymal cells is indispensable for skeletal development in mice

NAD is an essential co-factor for cellular energy metabolism and multiple other processes. Systemic NAD+ deficiency has been implicated in skeletal deformities during development in both humans and mice. NAD levels are maintained by multiple synthetic pathways but which ones are important in bone forming cells is unknown. Here, we generate mice with deletion of Nicotinamide Phosphoribosyltransferase (Nampt), a critical enzyme in the NAD salvage pathway, in all mesenchymal lineage cells of the limbs. At birth, NamptΔPrx1 exhibit dramatic limb shortening due to death of growth plate chondrocytes. Administration of the NAD precursor nicotinamide riboside during pregnancy prevents the majority of in utero defects. Depletion of NAD post-birth also promotes chondrocyte death, preventing further endochondral ossification and joint development. In contrast, osteoblast formation still occurs in knockout mice, in line with distinctly different microenvironments and reliance on redox reactions between chondrocytes and osteoblasts. These findings define a critical role for cell-autonomous NAD homeostasis during endochondral bone formation.

Systematic diet composition swap in a mouse genome-scale metabolic model reveals determinants of obesogenic diet metabolism in liver cancer

Dietary nutrient availability and gene expression, together, influence tissue metabolic activity. Here, we explore whether altering dietary nutrient composition in the context of mouse liver cancer suffices to overcome chronic gene expression changes that arise from tumorigenesis and western-style diet (WD). We construct a mouse genome-scale metabolic model and estimate metabolic fluxes in liver tumors and non-tumoral tissue after computationally varying the composition of input diet. This approach, called Systematic Diet Composition Swap (SyDiCoS), revealed that, compared to a control diet, WD increases production of glycerol and succinate irrespective of specific tissue gene expression patterns. Conversely, differences in fatty acid utilization pathways between tumor and non-tumor liver are amplified with WD by both dietary carbohydrates and lipids together. Our data suggest that combined dietary component modifications may be required to normalize the distinctive metabolic patterns that underlie selective targeting of tumor metabolism.

Genome-scale models as a vehicle for knowledge transfer from microbial to mammalian cell systems

With the plethora of omics data becoming available for mammalian cell and, increasingly, human cell systems, Genome-scale metabolic models (GEMs) have emerged as a useful tool for their organisation and analysis. The systems biology community has developed an array of tools for the solution, interrogation and customisation of GEMs as well as algorithms that enable the design of cells with desired phenotypes based on the multi-omics information contained in these models. However, these tools have largely found application in microbial cells systems, which benefit from smaller model size and ease of experimentation. Herein, we discuss the major outstanding challenges in the use of GEMs as a vehicle for accurately analysing data for mammalian cell systems and transferring methodologies that would enable their use to design strains and processes. We provide insights on the opportunities and limitations of applying GEMs to human cell systems for advancing our understanding of health and disease. We further propose their integration with data-driven tools and their enrichment with cellular functions beyond metabolism, which would, in theory, more accurately describe how resources are allocated intracellularly.

Compartmentalization of metabolism between cell types in multicellular organisms: a computational perspective

In multicellular organisms, metabolism is compartmentalized at many levels, including tissues and organs, different cell types, and subcellular compartments. Compartmentalization creates a coordinated homeostatic system where each compartment contributes to the production of energy and biomolecules the organism needs to carrying out specific metabolic tasks. Experimentally studying metabolic compartmentalization and metabolic interactions between cells and tissues in multicellular organisms is challenging at a systems level. However, recent progress in computational modeling provides an alternative approach to this problem. Here we discuss how integrating metabolic network modeling with omics data offers an opportunity to reveal metabolic states at the level of organs, tissues and, ultimately, individual cells. We review the current status of genome-scale metabolic network models in multicellular organisms, methods to study metabolic compartmentalization in silico, and insights gained from computational analyses. We also discuss outstanding challenges and provide perspectives for the future directions of the field.

Genome-scale metabolic network reconstruction of model animals as a platform for translational research

Genome-scale metabolic models (GEMs) are used extensively for analysis of mechanisms underlying human diseases and metabolic malfunctions. However, the lack of comprehensive and high-quality GEMs for model organisms restricts translational utilization of omics data accumulating from the use of various disease models. Here we present a unified platform of GEMs that covers five major model animals, including Mouse1 (Mus musculus), Rat1 (Rattus norvegicus), Zebrafish1 (Danio rerio), Fruitfly1 (Drosophila melanogaster), and Worm1 (Caenorhabditis elegans). These GEMs represent the most comprehensive coverage of the metabolic network by considering both orthology-based pathways and species-specific reactions. All GEMs can be interactively queried via the accompanying web portal Metabolic Atlas. Specifically, through integrative analysis of Mouse1 with RNA-sequencing data from brain tissues of transgenic mice we identified a coordinated up-regulation of lysosomal GM2 ganglioside and peptide degradation pathways which appears to be a signature metabolic alteration in Alzheimer's disease (AD) mouse models with a phenotype of amyloid precursor protein overexpression. This metabolic shift was further validated with proteomics data from transgenic mice and cerebrospinal fluid samples from human patients. The elevated lysosomal enzymes thus hold potential to be used as a biomarker for early diagnosis of AD. Taken together, we foresee that this evolving open-source platform will serve as an important resource to facilitate the development of systems medicines and translational biomedical applications.

Guided extraction of genome-scale metabolic models for the integration and analysis of omics data

Omics data can be integrated into a reference model using various model extraction methods (MEMs) to yield context-specific genome-scale metabolic models (GEMs). How to chose the appropriate MEM, thresholding rule and threshold remains a challenge. We integrated mouse transcriptomic data from a Cyp51 knockout mice diet experiment (GSE58271) using five MEMs (GIMME, iMAT, FASTCORE, INIT an tINIT) in a combination with a recently published mouse GEM iMM1865. Except for INIT and tINIT, the size of extracted models varied with the MEM used (t-test: p-value < 0.001). The Jaccard index of iMAT models ranged from 0.27 to 1.0. Out of the three factors under study in the experiment (diet, gender and genotype), gender explained most of the variability ( > 90%) in PC1 for FASTCORE. In iMAT, each of the three factors explained less than 40% of the variability within PC1, PC2 and PC3. Among all the MEMs, FASTCORE captured the most of the true variability in the data by clustering samples by gender. Our results show that for the efficient use of MEMs in the context of omics data integration and analysis, one should apply various MEMs, thresholding rules, and thresholding values to select the MEM and its configuration that best captures the true variability in the data. This selection can be guided by the methodology as proposed and used in this paper. Moreover, we describe certain approaches that can be used to analyse the results obtained with the selected MEM and to put these results in a biological context.

Systematic investigation of mouse models of Parkinson's disease by transcriptome mapping on a brain-specific genome-scale metabolic network

Genome-scale metabolic networks enable systemic investigation of metabolic alterations caused by diseases by providing interpretation of omics data. Although Mus musculus (mouse) is one of the most commonly used model organisms for neurodegenerative diseases, a brain-specific metabolic network model of mice has not yet been reconstructed. Here we reconstructed the first brain-specific metabolic network model of mice, iBrain674-Mm, by a homology-based approach, which consisted of 992 reactions controlled by 674 genes and distributed over 48 pathways. We validated the newly reconstructed network model by showing that it predicts healthy resting-state metabolic phenotypes of mouse brain compatible with the literature. We later used iBrain674-Mm to interpret various experimental mouse models of Parkinson's Disease (PD) at the transcriptome level. To this end, we applied a constraint-based modelling based biomarker prediction method called TIMBR (Transcriptionally Inferred Metabolic Biomarker Response) to predict altered metabolite production from transcriptomic data. Systemic analysis of seven different PD mouse models by TIMBR showed that the neuronal levels of glutamate, lactate, creatine phosphate, neuronal acetylcholine, bilirubin and formate increased in most of the PD mouse models, whereas the levels of melatonin, epinephrine, astrocytic formate and astrocytic bilirubin decreased. Although most of the predictions were consistent with the literature, there were some inconsistencies among different PD mouse models, signifying that there is no perfect experimental model to reflect PD metabolism. The newly reconstructed brain-specific genome-scale metabolic network model of mice can make important contributions to the interpretation and development of experimental mouse models of PD and other neurodegenerative diseases.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

The era of big data: Genome-scale modelling meets machine learning

With omics data being generated at an unprecedented rate, genome-scale modelling has become pivotal in its organisation and analysis. However, machine learning methods have been gaining ground in cases where knowledge is insufficient to represent the mechanisms underlying such data or as a means for data curation prior to attempting mechanistic modelling. We discuss the latest advances in genome-scale modelling and the development of optimisation algorithms for network and error reduction, intracellular constraining and applications to strain design. We further review applications of supervised and unsupervised machine learning methods to omics datasets from microbial and mammalian cell systems and present efforts to harness the potential of both modelling approaches through hybrid modelling.