Cancer Metabolism: A Modeling Perspective

Tumor cells alter their metabolism to maintain unregulated cellular proliferation and survival, but this transformation leaves them reliant on constant supply of nutrients and energy. In addition to the widely studied dysregulated glucose metabolism to fuel tumor cell growth, accumulating evidences suggest that utilization of amino acids and lipids contributes significantly to cancer cell metabolism. Also recent progresses in our understanding of carcinogenesis have revealed that cancer is a complex disease and cannot be understood through simple investigation of genetic mutations of cancerous cells. Cancer cells present in complex tumor tissues communicate with the surrounding microenvironment and develop traits which promote their growth, survival, and metastasis. Decoding the full scope and targeting dysregulated metabolic pathways that support neoplastic transformations and their preservation requires both the advancement of experimental technologies for more comprehensive measurement of omics as well as the advancement of robust computational methods for accurate analysis of the generated data. Here, we review cancer-associated reprogramming of metabolism and highlight the capability of genome-scale metabolic modeling approaches in perceiving a system-level perspective of cancer metabolism and in detecting novel selective drug targets

Personalized Cardiovascular Disease Prediction and Treatment-A Review of Existing Strategies and Novel Systems Medicine Tools

Cardiovascular disease (CVD) continues to constitute the leading cause of death globally. CVD risk stratification is an essential tool to sort through heterogeneous populations and identify individuals at risk of developing CVD. However, applications of current risk scores have recently been shown to result in considerable misclassification of high-risk subjects. In addition, despite long standing beneficial effects in secondary prevention, current CVD medications have in a primary prevention setting shown modest benefit in terms of increasing life expectancy. A systems biology approach to CVD risk stratification may be employed for improving risk-estimating algorithms through addition of high-throughput derived omics biomarkers. In addition, modeling of personalized benefit-of-treatment may help in guiding choice of intervention. In the area of medicine, realizing that CVD involves perturbations of large complex biological networks, future directions in drug development may involve moving away from a reductionist approach toward a system level approach. Here, we review current CVD risk scores and explore how novel algorithms could help to improve the identification of risk and maximize personalized treatment benefit. We also discuss possible future directions in the development of effective treatment strategies for CVD through the use of genome-scale metabolic models (GEMs) as well as other biological network-based approaches

Predicting growth of the healthy infant using a genome scale metabolic model

An estimated 165 million children globally have stunted growth, and extensive growth data are available. Genome scale metabolic models allow the simulation of molecular flux over each metabolic enzyme, and are well adapted to analyze biological systems. We used a human genome scale metabolic model to simulate the mechanisms of growth and integrate data about breast-milk intake and composition with the infant\u27s biomass and energy expenditure of major organs. The model predicted daily metabolic fluxes from birth to age 6 months, and accurately reproduced standard growth curves and changes in body composition. The model corroborates the finding that essential amino and fatty acids do not limit growth, but that energy is the main growth limiting factor. Disruptions to the supply and demand of energy markedly affected the predicted growth, indicating that elevated energy expenditure may be detrimental. The model was used to simulate the metabolic effect of mineral deficiencies, and showed the greatest growth reduction for deficiencies in copper, iron, and magnesium ions which affect energy production through oxidative phosphorylation. The model and simulation method were integrated to a platform and shared with the research community. The growth model constitutes another step towards the complete representation of human metabolism, and may further help improve the understanding of the mechanisms underlying stunting

Identification of Discriminating Metabolic Pathways and Metabolites in Human PBMCs Stimulated by Various Pathogenic Agents

Immunity and cellular metabolism are tightly interconnected but it is not clear whether different pathogens elicit specific metabolic responses. To address this issue, we studied differential metabolic regulation in peripheral blood mononuclear cells (PBMCs) of healthy volunteers challenged by Candida albicans, Borrelia burgdorferi, lipopolysaccharide, and Mycobacterium tuberculosis in vitro. By integrating gene expression data of stimulated PBMCs of healthy individuals with the KEGG pathways, we identified both common and pathogen-specific regulated pathways depending on the time of incubation. At 4 h of incubation, pathogenic agents inhibited expression of genes involved in both the glycolysis and oxidative phosphorylation pathways. In contrast, at 24 h of incubation, particularly glycolysis was enhanced while genes involved in oxidative phosphorylation remained unaltered in the PBMCs. In general, differential gene expression was less pronounced at 4 h compared to 24 h of incubation. KEGG pathway analysis allowed differentiation between effects induced by Candida and bacterial stimuli. Application of genome-scale metabolic model further generated a Candida-specific set of 103 reporter metabolites (e.g., desmosterol) that might serve as biomarkers discriminating Candida stimulated PBMCs from bacteria-stimuated PBMCs. Our analysis also identified a set of 49 metabolites that allowed discrimination between the effects of Borrelia burgdorferi, lipopolysaccharide and Mycobacterium tuberculosis. We conclude that analysis of pathogen-induced effects on PBMCs by a combination of KEGG pathways and genome-scale metabolic model provides deep insight in the metabolic changes coupled to host defense

The gut microbiota modulates host amino acid and glutathione metabolism in mice

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV-R) and germ-free (GF) mice using gene expression data and tissue-specific genome-scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue-specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV-R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N-acetylated AAs in the hepatic portal vein of CONV-R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV-R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice

Multi-Omics biological embeddings for ML-models

Machine learning algorithms have led to the development of numerous vector embeddings for biological entities such as metabolites, proteins, genes, and enzymes. However, these embeddings often lack contextual information due to their specialized focus on individual omics. Disease progression and biosynthesis pathways are increasingly understood through complex, multi-layered networks that integrate diverse omics data and intricate signaling and reaction sequences. Capturing these relationships in a meaningful way requires embeddings that account for both functional and multi-modal dependencies. We propose an embedding approach that unifies these different biological modalities by treating them as directions in a shared space rather than as isolated data types. Similar to how word embeddings in natural language processing reveal meaningful relationships (e.g., Tokyo – Japan + UK = London, indicating a directional representation of capitals), we can model genes and proteins in a way that captures their inherent connections. A gene implies information about the protein it encodes, and vice versa, forming a structured and interpretable representation of biological pathways. Our model, inspired by NLP techniques, breaks down pathway sequences into contextual pairs spanning different omics types. By aligning pathway steps in proximity, the embeddings reflect biologically relevant relationships, enhancing their interpretability and utility. Because these embeddings are generated from pathway sequences, they can be applied to optimize reaction pathways, aiding retrosynthesis in microbiomes, drug development, and even human health interventions

Multi-Omic Data Interpretation to Repurpose Subtype Specific Drug Candidates for Breast Cancer

Triple-negative breast cancer (TNBC), which is largely synonymous with the basal-like molecular subtype, is the 5th leading cause of cancer deaths for women in the United States. The overall prognosis for TNBC patients remains poor given that few treatment options exist; including targeted therapies (not FDA approved), and multi-agent chemotherapy as standard-of-care treatment. TNBC like other complex diseases is governed by the perturbations of the complex interaction networks thereby elucidating the underlying molecular mechanisms of this disease in the context of network principles, which have the potential to identify targets for drug development. Here, we present an integrated “omics” approach based on the use of transcriptome and interactome data to identify dynamic/active protein-protein interaction networks (PPINs) in TNBC patients. We have identified three highly connected modules, EED, DHX9, and AURKA, which are extremely activated in TNBC tumors compared to both normal tissues and other breast cancer subtypes. Based on the functional analyses, we propose that these modules are potential drivers of proliferation and, as such, should be considered candidate molecular targets for drug development or drug repositioning in TNBC. Consistent with this argument, we repurposed steroids, anti-inflammatory agents, anti-infective agents, cardiovascular agents for patients with basal-like breast cancer. Finally, we have performed essential metabolite analysis on personalized genome-scale metabolic models and found that metabolites such as sphingosine-1-phosphate and cholesterol-sulfate have utmost importance in TNBC tumor growth