Integrating systemic and molecular levels to infer key drivers sustaining metabolic adaptations.

Metabolic adaptations to complex perturbations, like the response to pharmacological treatments in multifactorial diseases such as cancer, can be described through measurements of part of the fluxes and concentrations at the systemic level and individual transporter and enzyme activities at the molecular level. In the framework of Metabolic Control Analysis (MCA), ensembles of linear constraints can be built integrating these measurements at both systemic and molecular levels, which are expressed as relative differences or changes produced in the metabolic adaptation. Here, combining MCA with Linear Programming, an efficient computational strategy is developed to infer additional non-measured changes at the molecular level that are required to satisfy these constraints. An application of this strategy is illustrated by using a set of fluxes, concentrations, and differentially expressed genes that characterize the response to cyclin-dependent kinases 4 and 6 inhibition in colon cancer cells. Decreases and increases in transporter and enzyme individual activities required to reprogram the measured changes in fluxes and concentrations are compared with down-regulated and up-regulated metabolic genes to unveil those that are key molecular drivers of the metabolic response.

Metabolomics: The Stethoscope for the Twenty-First Century.

Metabolomics encompasses the systematic identification and quantification of all metabolic products in the human body. This field could provide clinicians with novel sets of diagnostic biomarkers for disease states in addition to quantifying treatment response to medications at an individualized level. This literature review aims to highlight the technology underpinning metabolic profiling, identify potential applications of metabolomics in clinical practice, and discuss the translational challenges that the field faces. We searched PubMed, MEDLINE, and EMBASE for primary and secondary research articles regarding clinical applications of metabolomics. Metabolic profiling can be performed using mass spectrometry and nuclear magnetic resonance-based techniques using a variety of biological samples. This is carried out in vivo or in vitro following careful sample collection, preparation, and analysis. The potential clinical applications constitute disruptive innovations in their respective specialities, particularly oncology and metabolic medicine. Outstanding issues currently preventing widespread clinical use are scalability of data interpretation, standardization of sample handling practice, and e-infrastructure. Routine utilization of metabolomics at a patient and population level will constitute an integral part of future healthcare provision.

Interoperable and scalable data analysis with microservices: applications in metabolomics.

MOTIVATION: Developing a robust and performant data analysis workflow that integrates all necessary components whilst still being able to scale over multiple compute nodes is a challenging task. We introduce a generic method based on the microservice architecture, where software tools are encapsulated as Docker containers that can be connected into scientific workflows and executed using the Kubernetes container orchestrator. RESULTS: We developed a Virtual Research Environment (VRE) which facilitates rapid integration of new tools and developing scalable and interoperable workflows for performing metabolomics data analysis. The environment can be launched on-demand on cloud resources and desktop computers. IT-expertise requirements on the user side are kept to a minimum, and workflows can be re-used effortlessly by any novice user. We validate our method in the field of metabolomics on two mass spectrometry, one nuclear magnetic resonance spectroscopy and one fluxomics study. We showed that the method scales dynamically with increasing availability of computational resources. We demonstrated that the method facilitates interoperability using integration of the major software suites resulting in a turn-key workflow encompassing all steps for mass-spectrometry-based metabolomics including preprocessing, statistics and identification. Microservices is a generic methodology that can serve any scientific discipline and opens up for new types of large-scale integrative science. AVAILABILITY AND IMPLEMENTATION: The PhenoMeNal consortium maintains a web portal (https://portal.phenomenal-h2020.eu) providing a GUI for launching the Virtual Research Environment. The GitHub repository https://github.com/phnmnl/ hosts the source code of all projects. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

p13CMFA: Parsimonious 13C metabolic flux analysis.

Deciphering the mechanisms of regulation of metabolic networks subjected to perturbations, including disease states and drug-induced stress, relies on tracing metabolic fluxes. One of the most informative data to predict metabolic fluxes are 13C based metabolomics, which provide information about how carbons are redistributed along central carbon metabolism. Such data can be integrated using 13C Metabolic Flux Analysis (13C MFA) to provide quantitative metabolic maps of flux distributions. However, 13C MFA might be unable to reduce the solution space towards a unique solution either in large metabolic networks or when small sets of measurements are integrated. Here we present parsimonious 13C MFA (p13CMFA), an approach that runs a secondary optimization in the 13C MFA solution space to identify the solution that minimizes the total reaction flux. Furthermore, flux minimization can be weighted by gene expression measurements allowing seamless integration of gene expression data with 13C data. As proof of concept, we demonstrate how p13CMFA can be used to estimate intracellular flux distributions from 13C measurements and transcriptomics data. We have implemented p13CMFA in Iso2Flux, our in-house developed isotopic steady-state 13C MFA software. The source code is freely available on GitHub (https://github.com/cfoguet/iso2flux/releases/tag/0.7.2).

Network modules uncover mechanisms of skeletal muscle dysfunction in COPD patients.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) patients often show skeletal muscle dysfunction that has a prominent negative impact on prognosis. The study aims to further explore underlying mechanisms of skeletal muscle dysfunction as a characteristic systemic effect of COPD, potentially modifiable with preventive interventions (i.e. muscle training). The research analyzes network module associated pathways and evaluates the findings using independent measurements. METHODS: We characterized the transcriptionally active network modules of interacting proteins in the vastus lateralis of COPD patients (n = 15, FEV1 46 ± 12% pred, age 68 ± 7 years) and healthy sedentary controls (n = 12, age 65 ± 9  years), at rest and after an 8-week endurance training program. Network modules were functionally evaluated using experimental data derived from the same study groups. RESULTS: At baseline, we identified four COPD specific network modules indicating abnormalities in creatinine metabolism, calcium homeostasis, oxidative stress and inflammatory responses, showing statistically significant associations with exercise capacity (VO2 peak, Watts peak, BODE index and blood lactate levels) (P < 0.05 each), but not with lung function (FEV1). Training-induced network modules displayed marked differences between COPD and controls. Healthy subjects specific training adaptations were significantly associated with cell bioenergetics (P < 0.05) which, in turn, showed strong relationships with training-induced plasma metabolomic changes; whereas, effects of training in COPD were constrained to muscle remodeling. CONCLUSION: In summary, altered muscle bioenergetics appears as the most striking finding, potentially driving other abnormal skeletal muscle responses. Trial registration The study was based on a retrospectively registered trial (May 2017), ClinicalTrials.gov identifier: NCT03169270.

De novo MYC addiction as an adaptive response of cancer cells to CDK4/6 inhibition.

Cyclin-dependent kinases (CDK) are rational cancer therapeutic targets fraught with the development of acquired resistance by tumor cells. Through metabolic and transcriptomic analyses, we show that the inhibition of CDK4/6 leads to a metabolic reprogramming associated with gene networks orchestrated by the MYC transcription factor. Upon inhibition of CDK4/6, an accumulation of MYC protein ensues which explains an increased glutamine metabolism, activation of the mTOR pathway and blunting of HIF-1α-mediated responses to hypoxia. These MYC-driven adaptations to CDK4/6 inhibition render cancer cells highly sensitive to inhibitors of MYC, glutaminase or mTOR and to hypoxia, demonstrating that metabolic adaptations to antiproliferative drugs unveil new vulnerabilities that can be exploited to overcome acquired drug tolerance and resistance by cancer cells.

Metabolic Reprogramming and Dependencies Associated with Epithelial Cancer Stem Cells Independent of the Epithelial-Mesenchymal Transition Program.

In solid tumors, cancer stem cells (CSCs) can arise independently of epithelial-mesenchymal transition (EMT). In spite of recent efforts, the metabolic reprogramming associated with CSC phenotypes uncoupled from EMT is poorly understood. Here, by using metabolomic and fluxomic approaches, we identify major metabolic profiles that differentiate metastatic prostate epithelial CSCs (e-CSCs) from non-CSCs expressing a stable EMT. We have found that the e-CSC program in our cellular model is characterized by a high plasticity in energy substrate metabolism, including an enhanced Warburg effect, a greater carbon and energy source flexibility driven by fatty acids and amino acid metabolism and an essential reliance on the proton buffering capacity conferred by glutamine metabolism. An analysis of transcriptomic data yielded a metabolic gene signature for our e-CSCs consistent with the metabolomics and fluxomics analyses that correlated with tumor progression and metastasis in prostate cancer and in 11 additional cancer types. Interestingly, an integrated metabolomics, fluxomics, and transcriptomics analysis allowed us to identify key metabolic players regulated at the post-transcriptional level, suggesting potential biomarkers and therapeutic targets to effectively forestall metastasis. Stem Cells 2016;34:1163-1176.

HepatoDyn: A Dynamic Model of Hepatocyte Metabolism That Integrates 13C Isotopomer Data.

The liver performs many essential metabolic functions, which can be studied using computational models of hepatocytes. Here we present HepatoDyn, a highly detailed dynamic model of hepatocyte metabolism. HepatoDyn includes a large metabolic network, highly detailed kinetic laws, and is capable of dynamically simulating the redox and energy metabolism of hepatocytes. Furthermore, the model was coupled to the module for isotopic label propagation of the software package IsoDyn, allowing HepatoDyn to integrate data derived from 13C based experiments. As an example of dynamical simulations applied to hepatocytes, we studied the effects of high fructose concentrations on hepatocyte metabolism by integrating data from experiments in which rat hepatocytes were incubated with 20 mM glucose supplemented with either 3 mM or 20 mM fructose. These experiments showed that glycogen accumulation was significantly lower in hepatocytes incubated with medium supplemented with 20 mM fructose than in hepatocytes incubated with medium supplemented with 3 mM fructose. Through the integration of extracellular fluxes and 13C enrichment measurements, HepatoDyn predicted that this phenomenon can be attributed to a depletion of cytosolic ATP and phosphate induced by high fructose concentrations in the medium.

ChainRank, a chain prioritisation method for contextualisation of biological networks.

BACKGROUND: Advances in high throughput technologies and growth of biomedical knowledge have contributed to an exponential increase in associative data. These data can be represented in the form of complex networks of biological associations, which are suitable for systems analyses. However, these networks usually lack both, context specificity in time and space as well as the distinctive borders, which are usually assigned in the classical pathway view of molecular events (e.g. signal transduction). This complexity and high interconnectedness call for automated techniques that can identify smaller targeted subnetworks specific to a given research context (e.g. a disease scenario). RESULTS: Our method, named ChainRank, finds relevant subnetworks by identifying and scoring chains of interactions that link specific network components. Scores can be generated from integrating multiple general and context specific measures (e.g. experimental molecular data from expression to proteomics and metabolomics, literature evidence, network topology). The performance of the novel ChainRank method was evaluated on recreating selected signalling pathways from a human protein interaction network. Specifically, we recreated skeletal muscle specific signaling networks in healthy and chronic obstructive pulmonary disease (COPD) contexts. The analysis showed that ChainRank can identify main mediators of context specific molecular signalling. An improvement of up to factor 2.5 was shown in the precision of finding proteins of the recreated pathways compared to random simulation. CONCLUSIONS: ChainRank provides a framework, which can integrate several user-defined scores and evaluate their combined effect on ranking interaction chains linking input data sets. It can be used to contextualise networks, identify signaling and regulatory path amongst targeted genes or to analyse synthetic lethality in the context of anticancer therapy. ChainRank is implemented in R programming language and freely available at https://github.com/atenyi/ChainRank.

Sampling with poling-based flux balance analysis: optimal versus sub-optimal flux space analysis of Actinobacillus succinogenes.

BACKGROUND: Flux balance analysis is traditionally implemented to identify the maximum theoretical flux for some specified reaction and a single distribution of flux values for all the reactions present which achieve this maximum value. However it is well known that the uncertainty in reaction networks due to branches, cycles and experimental errors results in a large number of combinations of internal reaction fluxes which can achieve the same optimal flux value. RESULTS: In this work, we have modified the applied linear objective of flux balance analysis to include a poling penalty function, which pushes each new set of reaction fluxes away from previous solutions generated. Repeated poling-based flux balance analysis generates a sample of different solutions (a characteristic set), which represents all the possible functionality of the reaction network. Compared to existing sampling methods, for the purpose of generating a relatively "small" characteristic set, our new method is shown to obtain a higher coverage than competing methods under most conditions. The influence of the linear objective function on the sampling (the linear bias) constrains optimisation results to a subspace of optimal solutions all producing the same maximal fluxes. Visualisation of reaction fluxes plotted against each other in 2 dimensions with and without the linear bias indicates the existence of correlations between fluxes. This method of sampling is applied to the organism Actinobacillus succinogenes for the production of succinic acid from glycerol. CONCLUSIONS: A new method of sampling for the generation of different flux distributions (sets of individual fluxes satisfying constraints on the steady-state mass balances of intermediates) has been developed using a relatively simple modification of flux balance analysis to include a poling penalty function inside the resulting optimisation objective function. This new methodology can achieve a high coverage of the possible flux space and can be used with and without linear bias to show optimal versus sub-optimal solution spaces. Basic analysis of the Actinobacillus succinogenes system using sampling shows that in order to achieve the maximal succinic acid production CO2 must be taken into the system. Solutions involving release of CO2 all give sub-optimal succinic acid production.

Dual feedback loops in the GAL regulon suppress cellular heterogeneity in yeast.

Transcriptional noise is known to be an important cause of cellular heterogeneity and phenotypic variation. The extent to which molecular interaction networks may have evolved to either filter or exploit transcriptional noise is a much debated question. The yeast genetic network regulating galactose metabolism involves two proteins, Gal3p and Gal80p, that feed back positively and negatively, respectively, on GAL gene expression. Using kinetic modeling and experimental validation, we demonstrate that these feedback interactions together are important for (i) controlling the cell-to-cell variability of GAL gene expression and (ii) ensuring that cells rapidly switch to an induced state for galactose uptake.

Workforce preparation: the Biohealth computing model for Master and PhD students.

The article addresses the strategic role of workforce preparation in the process of adoption of Systems Medicine as a driver of biomedical research in the new health paradigm. It reports on relevant initiatives, like CASyM, fostering Systems Medicine at EU level. The chapter focuses on the BioHealth Computing Program as a reference for multidisciplinary training of future systems-oriented researchers describing the productive interactions with the Synergy-COPD project.

Relevance of the MEK/ERK signaling pathway in the metabolism of activated macrophages: a metabolomic approach.

The activation of immune cells in response to a pathogen involves a succession of signaling events leading to gene and protein expression, which requires metabolic changes to match the energy demands. The metabolic profile associated with the MAPK cascade (ERK1/2, p38, and JNK) in macrophages was studied, and the effect of its inhibition on the specific metabolic pattern of LPS stimulation was characterized. A [1,2-[(13)C](2)]glucose tracer-based metabolomic approach was used to examine the metabolic flux distribution in these cells after MEK/ERK inhibition. Bioinformatic tools were used to analyze changes in mass isotopomer distribution and changes in glucose and glutamine consumption and lactate production in basal and LPS-stimulated conditions in the presence and absence of the selective inhibitor of the MEK/ERK cascade, PD325901. Results showed that PD325901-mediated ERK1/2 inhibition significantly decreased glucose consumption and lactate production but did not affect glutamine consumption. These changes were accompanied by a decrease in the glycolytic flux, consistent with the observed decrease in fructose-2,6-bisphosphate concentration. The oxidative and nonoxidative pentose phosphate pathways and the ratio between them also decreased. However, tricarboxylic acid cycle flux did not change significantly. LPS activation led to the opposite responses, although all of these were suppressed by PD325901. However, LPS also induced a small decrease in pentose phosphate pathway fluxes and an increase in glutamine consumption that were not affected by PD325901. We concluded that inhibition of the MEK/ERK cascade interferes with central metabolism, and this cross-talk between signal transduction and metabolism also occurs in the presence of LPS.

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.

Carbon metabolism and the sign of control coefficients in metabolic adaptations underlying K-ras transformation.

Metabolic adaptations are associated with changes in enzyme activities. These adaptations are characterized by patterns of positive and negative changes in metabolic fluxes and concentrations of intermediate metabolites. Knowledge of the mechanism and parameters governing enzyme kinetics is rarely available. However, the signs-increases or decreases-of many of these changes can be predicted using the signs of metabolic control coefficients. These signs require the only knowledge of the structure of the metabolic network and a limited qualitative knowledge of the regulatory dependences, which is widely available for carbon metabolism. Here, as a case study, we identified control coefficients with fixed signs in order to predict the pattern of changes in key enzyme activities which can explain the observed changes in fluxes and concentrations underlying the metabolic adaptations in oncogenic K-ras transformation in NIH-3T3 cells. The fixed signs of control coefficients indicate that metabolic changes following the oncogenic transformation-increased glycolysis and oxidative branch of the pentose-phosphate pathway, and decreased concentration in sugar-phosphates-could be associated with increases in activity for glucose-6-phosphate dehydrogenase, pyruvate kinase and lactate dehydrogenase, and decrease for transketolase. These predictions were validated experimentally by measuring specific activities. We conclude that predictions based on fixed signs of control coefficients are a very robust tool for the identification of changes in enzyme activities that can explain observed metabolic adaptations in carbon metabolism.

Enzymatic and metabolic characterization of the phosphoglycerate kinase deficiency associated with chronic hemolytic anemia caused by the PGK-Barcelona mutation.

Recently, we reported a new mutation of phosphoglycerate kinase (PGK), called PGK-Barcelona, which causes chronic hemolytic anemia associated with progressive neurological impairment. We found a 140T→A substitution that produces an Ile46Asn change located at the N-domain of the enzyme and we suggested that the decrease of the PGK activity is probably related to a loss of enzyme stability. In this paper, by analyzing whole hemolysates and cloned enzymes, we show that both enzymes possess similar kinetic properties (although some differences are observed in the Km values) and the same electrophoretic mobility. However, PGK-Barcelona has higher thermal instability. Therefore, we confirm that the decrease of the red blood cell (RBC) PGK activity caused by the PGK-Barcelona mutation is more closely related to a loss of enzyme stability than to a decrease of enzyme catalytic function. Furthermore, we have measured the levels of glycolytic metabolites and adenine nucleotides in the RBC from controls and from the patient. The increase of 2,3-bisphosphoglycerate and the decrease of ATP RBC levels are the only detected metabolic changes that could cause hemolytic anemia.

Metabolic network adaptations in cancer as targets for novel therapies.

Metabolite concentrations and fluxes are the system variables that characterize metabolism. The systematic study of metabolite profiles is known as metabolomics; however, knowledge of the complete set of metabolites may not be enough to predict distinct phenotypes. A complete understanding of metabolic processes requires detailed knowledge of enzyme-controlled intracellular fluxes. These can be estimated through quantitative measurements of metabolites at different times or by analysing the stable isotope patterns obtained after incubation with labelled substrates. We have identified distinct intracellular fluxes associated with metabolic adaptations accompanying cancer. The maintenance of an imbalance between fluxes for the oxidative and non-oxidative PPP (pentose phosphate pathway) has been shown to be critical for angiogenesis and cancer cell survival. Mouse NIH 3T3 cells transformed by different mutated K-ras oncogenes have differential routing of glucose to anaerobic glycolysis, the PPP and the Krebs cycle. These results indicate that knowledge of metabolic fingerprints associated with an altered genetic profile could be exploited in the rational design of new therapies. We conclude that the understanding of the multifactorial nature of metabolic adaptations in cancer may open new ways to develop novel multi-hit antitumoral therapies.

Software Supporting a Workflow of Quantitative Dynamic Flux Maps Estimation in Central Metabolism from SIRM Experimental Data.

Stable isotope-resolved metabolomics (SIRM), based on the analysis of biological samples from living cells incubated with artificial isotope enriched substrates, enables mapping the rates of biochemical reactions (metabolic fluxes). We developed software supporting a workflow of analysis of SIRM data obtained with mass spectrometry (MS). The evaluation of fluxes starting from raw MS recordings requires at least three steps of computer support: first, extraction of mass spectra of metabolites of interest, then correction of the spectra for natural isotope abundance, and finally, evaluation of fluxes by simulation of the corrected spectra using a corresponding mathematical model. A kinetic model based on ordinary differential equations (ODEs) for isotopomers of metabolites of the corresponding biochemical network supports the final part of the analysis, which provides a dynamic flux map.

In silico strategy to rationally engineer metabolite production: A case study for threonine in Escherichia coli.

Genetic engineering of metabolic pathways is a standard strategy to increase the production of metabolites of economic interest. However, such flux increases could very likely lead to undesirable changes in metabolite concentrations, producing deleterious perturbations on other cellular processes. These negative effects could be avoided by implementing a balanced increase of enzyme concentrations according to the Universal Method [Kacser and Acerenza (1993) Eur J Biochem 216:361-367]. Exact application of the method usually requires modification of many reactions, which is difficult to achieve in practice. Here, improvement of threonine production via pyruvate kinase deletion in Escherichia coli is used as a case study to demonstrate a partial application of the Universal Method, which includes performing sensitivity analysis. Our analysis predicts that manipulating a few reactions is sufficient to obtain an important increase in threonine production without major perturbations of metabolite concentrations.

PhenoMeNal: processing and analysis of metabolomics data in the cloud.

BACKGROUND: Metabolomics is the comprehensive study of a multitude of small molecules to gain insight into an organism's metabolism. The research field is dynamic and expanding with applications across biomedical, biotechnological, and many other applied biological domains. Its computationally intensive nature has driven requirements for open data formats, data repositories, and data analysis tools. However, the rapid progress has resulted in a mosaic of independent, and sometimes incompatible, analysis methods that are difficult to connect into a useful and complete data analysis solution. FINDINGS: PhenoMeNal (Phenome and Metabolome aNalysis) is an advanced and complete solution to set up Infrastructure-as-a-Service (IaaS) that brings workflow-oriented, interoperable metabolomics data analysis platforms into the cloud. PhenoMeNal seamlessly integrates a wide array of existing open-source tools that are tested and packaged as Docker containers through the project's continuous integration process and deployed based on a kubernetes orchestration framework. It also provides a number of standardized, automated, and published analysis workflows in the user interfaces Galaxy, Jupyter, Luigi, and Pachyderm. CONCLUSIONS: PhenoMeNal constitutes a keystone solution in cloud e-infrastructures available for metabolomics. PhenoMeNal is a unique and complete solution for setting up cloud e-infrastructures through easy-to-use web interfaces that can be scaled to any custom public and private cloud environment. By harmonizing and automating software installation and configuration and through ready-to-use scientific workflow user interfaces, PhenoMeNal has succeeded in providing scientists with workflow-driven, reproducible, and shareable metabolomics data analysis platforms that are interfaced through standard data formats, representative datasets, versioned, and have been tested for reproducibility and interoperability. The elastic implementation of PhenoMeNal further allows easy adaptation of the infrastructure to other application areas and 'omics research domains.