The landscape of tiered regulation of breast cancer cell metabolism

Control analysis in the identification of key enzymes driving metabolic adaptations: Towards drug target discovery

Metabolic Control Analysis (MCA) marked a turning point in understanding the design principles of metabolic network control by establishing control coefficients as a means to quantify the degree of control that an enzyme exerts on flux or metabolite concentrations. MCA has demonstrated that control of metabolic pathways is distributed among many enzymes rather than depending on a single rate-limiting step. MCA also proved that this distribution depends not only on the stoichiometric structure of the network but also on other kinetic determinants, such as the degree of saturation of the enzyme active site, the distance to thermodynamic equilibrium, and metabolite feedback regulatory loops. Consequently, predicting the alterations that occur during metabolic adaptation in response to strong changes involving a redistribution in such control distribution can be challenging. Here, using the framework provided by MCA, we illustrate how control distribution in a metabolic pathway/network depends on enzyme kinetic determinants and to what extent the redistribution of control affects our predictions on candidate enzymes suitable as targets for small molecule inhibition in the drug discovery process. Our results uncover that kinetic determinants can lead to unexpected control distribution and outcomes that cannot be predicted solely from stoichiometric determinants. We also unveil that the inference of key enzyme-drivers of an observed metabolic adaptation can be dramatically improved using mean control coefficients and ruling out those enzyme activities that are associated with low control coefficients. As the use of constraint-based stoichiometric genome-scale metabolic models (GSMMs) becomes increasingly prevalent for identifying genes/enzymes that could be potential drug targets, we anticipate that incorporating kinetic determinants and ruling out enzymes with low control coefficients into GSMM workflows will facilitate more accurate predictions and reveal novel therapeutic targets.

Comparing mechanism-based and machine learning models for predicting the effects of glucose accessibility on tumor cell proliferation

Glucose plays a central role in tumor metabolism and development and is a target for novel therapeutics. To characterize the response of cancer cells to blockade of glucose uptake, we collected time-resolved microscopy data to track the growth of MDA-MB-231 breast cancer cells. We then developed a mechanism-based, mathematical model to predict how a glucose transporter (GLUT1) inhibitor (Cytochalasin B) influences the growth of the MDA-MB-231 cells by limiting access to glucose. The model includes a parameter describing dose dependent inhibition to quantify both the total glucose level in the system and the glucose level accessible to the tumor cells. Four common machine learning models were also used to predict tumor cell growth. Both the mechanism-based and machine learning models were trained and validated, and the prediction error was evaluated by the coefficient of determination (R²). The random forest model provided the highest accuracy predicting cell dynamics (R² = 0.92), followed by the decision tree (R² = 0.89), k-nearest-neighbor regression (R² = 0.84), mechanism-based (R² = 0.77), and linear regression model (R² = 0.69). Thus, the mechanism-based model has a predictive capability comparable to machine learning models with the added benefit of elucidating biological mechanisms.

A pipeline for the reconstruction and evaluation of context-specific human metabolic models at a large-scale

Constraint-based (CB) metabolic models provide a mathematical framework and scaffold for in silico cell metabolism analysis and manipulation. In the past decade, significant efforts have been done to model human metabolism, enabled by the increased availability of multi-omics datasets and curated genome-scale reconstructions, as well as the development of several algorithms for context-specific model (CSM) reconstruction. Although CSM reconstruction has revealed insights on the deregulated metabolism of several pathologies, the process of reconstructing representative models of human tissues still lacks benchmarks and appropriate integrated software frameworks, since many tools required for this process are still disperse across various software platforms, some of which are proprietary. In this work, we address this challenge by assembling a scalable CSM reconstruction pipeline capable of integrating transcriptomics data in CB models. We combined omics preprocessing methods inspired by previous efforts with in-house implementations of existing CSM algorithms and new model refinement and validation routines, all implemented in the Troppo Python-based open-source framework. The pipeline was validated with multi-omics datasets from the Cancer Cell Line Encyclopedia (CCLE), also including reference fluxomics measurements for the MCF7 cell line. We reconstructed over 6000 models based on the Human-GEM template model for 733 cell lines featured in the CCLE, using MCF7 models as reference to find the best parameter combinations. These reference models outperform earlier studies using the same template by comparing gene essentiality and fluxomics experiments. We also analysed the heterogeneity of breast cancer cell lines, identifying key changes in metabolism related to cancer aggressiveness. Despite the many challenges in CB modelling, we demonstrate using our pipeline that combining transcriptomics data in metabolic models can be used to investigate key metabolic shifts. Significant limitations were found on these models ability for reliable quantitative flux prediction, thus motivating further work in genome-wide phenotype prediction.

Genome-scale models of human metabolism are promising tools capable of contextualising large omics datasets within a framework that enables analysis and manipulation of metabolic phenotypes. Despite various successes in applying these methods to provide mechanistic hypotheses for deregulated metabolism in disease, there is no standardized workflow to extract these models using existing methods and the tools required to do so are mostly implemented using proprietary software. We have assembled a generic pipeline to extract and validate context-specific metabolic models using multi-omics datasets and implemented it using the troppo framework. We first validate our pipeline using MCF7 cell line models and assess their ability to predict lethal gene knockouts as well as flux activity using multi-omics data. We also demonstrate how this approach can be generalized for large-scale transcriptomics datasets and used to generate insights on the metabolic heterogeneity of cancer and relevant features for other data mining approaches. The pipeline is available as part of an open-source framework that is generic for a variety of applications.

Exploring the metabolic landscape of pancreatic ductal adenocarcinoma cells using genome-scale metabolic modeling

Pancreatic ductal adenocarcinoma (PDAC) is a major research focus because of its poor therapy response and dismal prognosis. PDAC cells adapt their metabolism to the surrounding environment, often relying on diverse nutrient sources. Because traditional experimental techniques appear exhaustive to find a viable therapeutic strategy, a highly curated and omics-informed PDAC genome-scale metabolic model was reconstructed using patient-specific transcriptomics data. From the model-predictions, several new metabolic functions were explored as potential therapeutic targets in addition to the known metabolic hallmarks of PDAC. Significant downregulation in the peroxisomal beta oxidation pathway, flux modulation in the carnitine shuttle system, and upregulation in the reactive oxygen species detoxification pathway reactions were observed. These unique metabolic traits of PDAC were correlated with potential drug combinations targeting genes with poor prognosis in PDAC. Overall, this study provides a better understanding of the metabolic vulnerabilities in PDAC and will lead to novel effective therapeutic strategies.

Cardio-onco-metabolism: metabolic remodelling in cardiovascular disease and cancer

Cardiovascular disease and cancer are the two leading causes of morbidity and mortality in the world. The emerging field of cardio-oncology has revealed that these seemingly disparate disease processes are intertwined, owing to the cardiovascular sequelae of anticancer therapies, shared risk factors that predispose individuals to both cardiovascular disease and cancer, as well the possible potentiation of cancer growth by cardiac dysfunction. As a result, interest has increased in understanding the fundamental biological mechanisms that are central to the relationship between cardiovascular disease and cancer. Metabolism, appropriate regulation of energy, energy substrate utilization, and macromolecular synthesis and breakdown are fundamental processes for cellular and organismal survival. In this Review, we explore the emerging data identifying metabolic dysregulation as an important theme in cardio-oncology. We discuss the growing recognition of metabolic reprogramming in cardiovascular disease and cancer and view the novel area of cardio-oncology through the lens of metabolism.

INTEGRATE: Model-based multi-omics data integration to characterize multi-level metabolic regulation

Metabolism is directly and indirectly fine-tuned by a complex web of interacting regulatory mechanisms that fall into two major classes. On the one hand, the expression level of the catalyzing enzyme sets the maximal theoretical flux level (i.e., the net rate of the reaction) for each enzyme-controlled reaction. On the other hand, metabolic regulation controls the metabolic flux through the interactions of metabolites (substrates, cofactors, allosteric modulators) with the responsible enzyme. High-throughput data, such as metabolomics and transcriptomics data, if analyzed separately, do not accurately characterize the hierarchical regulation of metabolism outlined above. They must be integrated to disassemble the interdependence between different regulatory layers controlling metabolism. To this aim, we propose INTEGRATE, a computational pipeline that integrates metabolomics and transcriptomics data, using constraint-based stoichiometric metabolic models as a scaffold. We compute differential reaction expression from transcriptomics data and use constraint-based modeling to predict if the differential expression of metabolic enzymes directly originates differences in metabolic fluxes. In parallel, we use metabolomics to predict how differences in substrate availability translate into differences in metabolic fluxes. We discriminate fluxes regulated at the metabolic and/or gene expression level by intersecting these two output datasets. We demonstrate the pipeline using a set of immortalized normal and cancer breast cell lines. In a clinical setting, knowing the regulatory level at which a given metabolic reaction is controlled will be valuable to inform targeted, truly personalized therapies in cancer patients.

Argininosuccinate lyase is a metabolic vulnerability in breast development and cancer

Epithelial-to-mesenchymal transition (EMT) is fundamental to both normal tissue development and cancer progression. We hypothesized that EMT plasticity defines a range of metabolic phenotypes and that individual breast epithelial metabolic phenotypes are likely to fall within this phenotypic landscape. To determine EMT metabolic phenotypes, the metabolism of EMT was described within genome-scale metabolic models (GSMMs) using either transcriptomic or proteomic data from the breast epithelial EMT cell culture model D492. The ability of the different data types to describe breast epithelial metabolism was assessed using constraint-based modeling which was subsequently verified using ¹³C isotope tracer analysis. The application of proteomic data to GSMMs provided relatively higher accuracy in flux predictions compared to the transcriptomic data. Furthermore, the proteomic GSMMs predicted altered cholesterol metabolism and increased dependency on argininosuccinate lyase (ASL) following EMT which were confirmed in vitro using drug assays and siRNA knockdown experiments. The successful verification of the proteomic GSMMs afforded iBreast2886, a breast GSMM that encompasses the metabolic plasticity of EMT as defined by the D492 EMT cell culture model. Analysis of breast tumor proteomic data using iBreast2886 identified vulnerabilities within arginine metabolism that allowed prognostic discrimination of breast cancer patients on a subtype-specific level. Taken together, we demonstrate that the metabolic reconstruction iBreast2886 formalizes the metabolism of breast epithelial cell development and can be utilized as a tool for the functional interpretation of high throughput clinical data.

Glutamine Modulates Expression and Function of Glucose 6-Phosphate Dehydrogenase via NRF2 in Colon Cancer Cells

Nucleotide pools need to be constantly replenished in cancer cells to support cell proliferation. The synthesis of nucleotides requires glutamine and 5-phosphoribosyl-1-pyrophosphate produced from ribose-5-phosphate via the oxidative branch of the pentose phosphate pathway (ox-PPP). Both PPP and glutamine also play a key role in maintaining the redox status of cancer cells. Enhanced glutamine metabolism and increased glucose 6-phosphate dehydrogenase (G6PD) expression have been related to a malignant phenotype in tumors. However, the association between G6PD overexpression and glutamine consumption in cancer cell proliferation is still incompletely understood. In this study, we demonstrated that both inhibition of G6PD and glutamine deprivation decrease the proliferation of colon cancer cells and induce cell cycle arrest and apoptosis. Moreover, we unveiled that glutamine deprivation induce an increase of G6PD expression that is mediated through the activation of the nuclear factor (erythroid-derived 2)-like 2 (NRF2). This crosstalk between G6PD and glutamine points out the potential of combined therapies targeting oxidative PPP enzymes and glutamine catabolism to combat colon cancer.

Genome-scale integration of transcriptome and metabolome unveils squalene synthase and dihydrofolate reductase as targets against AML cells resistant to chemotherapy

The development of resistance to chemotherapeutic agents, such as Doxorubicin (DOX) and cytarabine (AraC), is one of the greatest challenges to the successful treatment of Acute Myeloid Leukemia (AML). Such acquisition is often underlined by a metabolic reprogramming that can provide a therapeutic opportunity, as it can lead to the emergence of vulnerabilities and dependencies to be exploited as targets against the resistant cells. In this regard, genome-scale metabolic models (GSMMs) have emerged as powerful tools to integrate multiple layers of data to build cancer-specific models and identify putative metabolic vulnerabilities. Here, we use genome-scale metabolic modelling to reconstruct a GSMM of the THP1 AML cell line and two derivative cell lines, one with acquired resistance to AraC and the second with acquired resistance to DOX. We also explore how, adding to the transcriptomic layer, the metabolomic layer enhances the selectivity of the resulting condition specific reconstructions. The resulting models enabled us to identify and experimentally validate that drug-resistant THP1 cells are sensitive to the FDA-approved antifolate methotrexate. Moreover, we discovered and validated that the resistant cell lines could be selectively targeted by inhibiting squalene synthase, providing a new and promising strategy to directly inhibit cholesterol synthesis in AML drug resistant cells.

Incorporating Machine Learning into Established Bioinformatics Frameworks

The exponential growth of biomedical data in recent years has urged the application of numerous machine learning techniques to address emerging problems in biology and clinical research. By enabling the automatic feature extraction, selection, and generation of predictive models, these methods can be used to efficiently study complex biological systems. Machine learning techniques are frequently integrated with bioinformatic methods, as well as curated databases and biological networks, to enhance training and validation, identify the best interpretable features, and enable feature and model investigation. Here, we review recently developed methods that incorporate machine learning within the same framework with techniques from molecular evolution, protein structure analysis, systems biology, and disease genomics. We outline the challenges posed for machine learning, and, in particular, deep learning in biomedicine, and suggest unique opportunities for machine learning techniques integrated with established bioinformatics approaches to overcome some of these challenges.

Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches

Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.

Causal integration of multi-omics data with prior knowledge to generate mechanistic hypotheses

Multi-omics datasets can provide molecular insights beyond the sum of individual omics. Various tools have been recently developed to integrate such datasets, but there are limited strategies to systematically extract mechanistic hypotheses from them. Here, we present COSMOS (Causal Oriented Search of Multi-Omics Space), a method that integrates phosphoproteomics, transcriptomics, and metabolomics datasets. COSMOS combines extensive prior knowledge of signaling, metabolic, and gene regulatory networks with computational methods to estimate activities of transcription factors and kinases as well as network-level causal reasoning. COSMOS provides mechanistic hypotheses for experimental observations across multi-omics datasets. We applied COSMOS to a dataset comprising transcriptomics, phosphoproteomics, and metabolomics data from healthy and cancerous tissue from eleven clear cell renal cell carcinoma (ccRCC) patients. COSMOS was able to capture relevant crosstalks within and between multiple omics layers, such as known ccRCC drug targets. We expect that our freely available method will be broadly useful to extract mechanistic insights from multi-omics studies.

Systems metabolomics: from metabolomic snapshots to design principles

Metabolomics is a rapidly expanding technology that finds increasing application in a variety of fields, form metabolic disorders to cancer, from nutrition and wellness to design and optimization of cell factories. The integration of metabolic snapshots with metabolic fluxes, physiological readouts, metabolic models, and knowledge-informed Artificial Intelligence tools, is required to obtain a system-level understanding of metabolism. The emerging power of multi-omic approaches and the development of integrated experimental and computational tools, able to dissect metabolic features at cellular and subcellular resolution, provide unprecedented opportunities for understanding design principles of metabolic (dis)regulation and for the development of precision therapies in multifactorial diseases, such as cancer and neurodegenerative diseases.