Edelfosine-induced metabolic changes in cancer cells that precede the overproduction of reactive oxygen species and apoptosis

Bayesian kinetic modeling for tracer-based metabolomic data

BACKGROUND

Stable Isotope Resolved Metabolomics (SIRM) is a new biological approach that uses stable isotope tracers such as uniformly [Formula: see text]-enriched glucose ([Formula: see text]-Glc) to trace metabolic pathways or networks at the atomic level in complex biological systems. Non-steady-state kinetic modeling based on SIRM data uses sets of simultaneous ordinary differential equations (ODEs) to quantitatively characterize the dynamic behavior of metabolic networks. It has been increasingly used to understand the regulation of normal metabolism and dysregulation in the development of diseases. However, fitting a kinetic model is challenging because there are usually multiple sets of parameter values that fit the data equally well, especially for large-scale kinetic models. In addition, there is a lack of statistically rigorous methods to compare kinetic model parameters between different experimental groups.

RESULTS

We propose a new Bayesian statistical framework to enhance parameter estimation and hypothesis testing for non-steady-state kinetic modeling of SIRM data. For estimating kinetic model parameters, we leverage the prior distribution not only to allow incorporation of experts' knowledge but also to provide robust parameter estimation. We also introduce a shrinkage approach for borrowing information across the ensemble of metabolites to stably estimate the variance of an individual isotopomer. In addition, we use a component-wise adaptive Metropolis algorithm with delayed rejection to perform efficient Monte Carlo sampling of the posterior distribution over high-dimensional parameter space. For comparing kinetic model parameters between experimental groups, we propose a new reparameterization method that converts the complex hypothesis testing problem into a more tractable parameter estimation problem. We also propose an inference procedure based on credible interval and credible value. Our method is freely available for academic use at https://github.com/xuzhang0131/MCMCFlux .

CONCLUSIONS

Our new Bayesian framework provides robust estimation of kinetic model parameters and enables rigorous comparison of model parameters between experimental groups. Simulation studies and application to a lung cancer study demonstrate that our framework performs well for non-steady-state kinetic modeling of SIRM data.

Rafting on the Plasma Membrane: Lipid Rafts in Signaling and Disease

The plasma membrane is not a uniform phospholipid bilayer; it has specialized membrane nano- or microdomains called lipid rafts. Lipid rafts are small cholesterol and sphingolipid-rich plasma membrane islands. Although their existence was long debated, their presence in the plasma membrane of living cells is now well accepted with the advent of super-resolution imaging techniques. It is interesting to note that lipid rafts function to compartmentalize receptors and their regulators and substantially modulate cellular signaling. In this review, we will examine the role of lipid rafts and caveolae-lipid raft-like microdomains with a distinct 3D morphology-in cellular signaling. Moreover, we will investigate how raft compartmentalized signaling regulates diverse physiological processes such as proliferation, apoptosis, immune signaling, and development. Also, the deregulation of lipid raft-mediated signaling during tumorigenesis and metastasis will be explored.

The landscape of tiered regulation of breast cancer cell metabolism

Altered metabolism is a hallmark of cancer, but little is still known about its regulation. In this study, we measure transcriptomic, proteomic, phospho-proteomic and fluxomics data in a breast cancer cell-line (MCF7) across three different growth conditions. Integrating these multiomics data within a genome scale human metabolic model in combination with machine learning, we systematically chart the different layers of metabolic regulation in breast cancer cells, predicting which enzymes and pathways are regulated at which level. We distinguish between two types of reactions, directly and indirectly regulated. Directly-regulated reactions include those whose flux is regulated by transcriptomic alterations (~890) or via proteomic or phospho-proteomics alterations (~140) in the enzymes catalyzing them. We term the reactions that currently lack evidence for direct regulation as (putative) indirectly regulated (~930). Many metabolic pathways are predicted to be regulated at different levels, and those may change at different media conditions. Remarkably, we find that the flux of predicted indirectly regulated reactions is strongly coupled to the flux of the predicted directly regulated ones, uncovering a tiered hierarchical organization of breast cancer cell metabolism. Furthermore, the predicted indirectly regulated reactions are predominantly reversible. Taken together, this architecture may facilitate rapid and efficient metabolic reprogramming in response to the varying environmental conditions incurred by the tumor cells. The approach presented lays a conceptual and computational basis for mapping metabolic regulation in additional cancers.

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.

Primary Metabolism and Medium-Chain Fatty Acid Alterations Precede Long-Chain Fatty Acid Changes Impacting Neutral Lipid Metabolism in Response to an Anticancer Lysophosphatidylcholine Analogue in Yeast

The nonmetabolizable lysophosphatidylcholine (LysoPC) analogue edelfosine is the prototype of a class of compounds being investigated for their potential as selective chemotherapeutic agents. Edelfosine targets membranes, disturbing cellular homeostasis. Is not clear at this point how membrane alterations are communicated between intracellular compartments leading to growth inhibition and eventual cell death. In the present study, a combined metabolomics/lipidomics approach for the unbiased identification of metabolic pathways altered in yeast treated with sublethal concentrations of the LysoPC analogue was employed. Mass spectrometry of polar metabolites, fatty acids, and lipidomic profiling was used to study the effects of edelfosine on yeast metabolism. Amino acid and sugar metabolism, the Krebs cycle, and fatty acid profiles were most disrupted, with polar metabolites and short-medium chain fatty acid changes preceding long and very long-chain fatty acid variations. Initial increases in metabolites such as trehalose, proline, and γ-amino butyric acid with a concomitant decrease in metabolites of the Krebs cycle, citrate and fumarate, are interpreted as a cellular attempt to offset oxidative stress in response to mitochondrial dysfunction induced by the treatment. Notably, alanine, inositol, and myristoleic acid showed a steady increase during the period analyzed (2, 4, and 6 h after treatment). Of importance was the finding that edelfosine induced significant alterations in neutral glycerolipid metabolism resulting in a significant increase in the signaling lipid diacylglycerol.

MIDcor, an R-program for deciphering mass interferences in mass spectra of metabolites enriched in stable isotopes

Background

Tracing stable isotopes, such as ¹³C using various mass spectrometry (MS) methods provides a valuable information necessary for the study of biochemical processes in cells. However, extracting such information requires special care, such as a correction for naturally occurring isotopes, or overlapping mass spectra of various components of the cell culture medium. Developing a method for a correction of overlapping peaks is the primary objective of this study.

Results

Our computer program-MIDcor (free at https://github.com/seliv55/mid_correct) written in the R programming language, corrects the raw MS spectra both for the naturally occurring isotopes and for the overlapping of peaks corresponding to various substances. To this end, the mass spectra of unlabeled metabolites measured in two media are necessary: in a minimal medium containing only derivatized metabolites and chemicals for derivatization, and in a complete cell incubated medium. The MIDcor program calculates the difference (D) between the theoretical and experimentally measured spectra of metabolites containing only the naturally occurring isotopes. The result of comparison of D in the two media determines a way of deciphering the true spectra. (1) If D in the complete medium is greater than that in the minimal medium in at least one peak, then unchanged D is subtracted from the raw spectra of the labeled metabolite. (2) If D does not depend on the medium, then the spectrum probably overlaps with a derivatized fragment of the same metabolite, and D is modified proportionally to the metabolite labeling. The program automatically reaches a decision regarding the way of correction. For some metabolites/fragments in the case (2) D was found to decrease when the tested substance was ¹³C labeled, and this isotopic effect also can be corrected automatically, if the user provides a measured spectrum of the substance in which the ¹³C labeling is known a priori.

Conclusion

Using the developed program improves the reliability of stable isotope tracer data analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1513-3) contains supplementary material, which is available to authorized users.

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.

Applications of NMR spectroscopy to systems biochemistry

The past decades of advancements in NMR have made it a very powerful tool for metabolic research. Despite its limitations in sensitivity relative to mass spectrometric techniques, NMR has a number of unparalleled advantages for metabolic studies, most notably the rigor and versatility in structure elucidation, isotope-filtered selection of molecules, and analysis of positional isotopomer distributions in complex mixtures afforded by multinuclear and multidimensional experiments. In addition, NMR has the capacity for spatially selective in vivo imaging and dynamical analysis of metabolism in tissues of living organisms. In conjunction with the use of stable isotope tracers, NMR is a method of choice for exploring the dynamics and compartmentation of metabolic pathways and networks, for which our current understanding is grossly insufficient. In this review, we describe how various direct and isotope-edited 1D and 2D NMR methods can be employed to profile metabolites and their isotopomer distributions by stable isotope-resolved metabolomic (SIRM) analysis. We also highlight the importance of sample preparation methods including rapid cryoquenching, efficient extraction, and chemoselective derivatization to facilitate robust and reproducible NMR-based metabolomic analysis. We further illustrate how NMR has been applied in vitro, ex vivo, or in vivo in various stable isotope tracer-based metabolic studies, to gain systematic and novel metabolic insights in different biological systems, including human subjects. The pathway and network knowledge generated from NMR- and MS-based tracing of isotopically enriched substrates will be invaluable for directing functional analysis of other 'omics data to achieve understanding of regulation of biochemical systems, as demonstrated in a case study. Future developments in NMR technologies and reagents to enhance both detection sensitivity and resolution should further empower NMR in systems biochemical research.

How to measure metabolic fluxes: a taxonomic guide for (13)C fluxomics

Metabolic reaction rates (fluxes) contribute fundamentally to our understanding of metabolic phenotypes and mechanisms of cellular regulation. Stable isotope-based fluxomics integrates experimental data with biochemical networks and mathematical modeling to 'measure' the in vivo fluxes within an organism that are not directly observable. In recent years, (13)C fluxomics has evolved into a technology with great experimental, analytical, and mathematical diversity. This review aims at establishing a unified taxonomy by means of which the various fluxomics methods can be compared to each other. By linking the developed modeling approaches to recent studies, their challenges and opportunities are put into perspective. The proposed classification serves as a guide for scientific 'travelers' who are striving to resolve research questions with the currently available (13)C fluxomics toolset.

(13)C metabolic flux analysis shows that resistin impairs the metabolic response to insulin in L6E9 myotubes

Background

It has been suggested that the adipokine resistin links obesity and insulin resistance, although how resistin acts on muscle metabolism is controversial. We aimed to quantitatively analyse the effects of resistin on the glucose metabolic flux profile and on insulin response in L6E9 myotubes at the metabolic level using a tracer-based metabolomic approach and our in-house developed software, Isodyn.

Results

Resistin significantly increased glucose uptake and glycolysis, altering pyruvate utilisation by the cell. In the presence of resistin, insulin only slightly increased glucose uptake and glycolysis, and did not alter the flux profile around pyruvate induced by resistin. Resistin prevented the increase in gene expression in pyruvate dehydrogenase-E1 and the sharp decrease in gene expression in cytosolic phosphoenolpyruvate carboxykinase-1 induced by insulin.

Conclusions

These data suggest that resistin impairs the metabolic activation of insulin. This impairment cannot be explained by the activity of a single enzyme, but instead due to reorganisation of the whole metabolic flux distribution.

Edelfosine lipid nanosystems overcome drug resistance in leukemic cell lines

Although current therapies have improved leukemia survival rates, adverse drug effects and relapse are frequent. Encapsulation of edelfosine (ET) in lipid nanoparticles (LNs) improves its oral bioavailability and decreases its toxicity. Here we evaluated the efficacy of ET-LN in myeloid leukemia cell lines. Drug-loaded LN were as effective as free ET in sensitive leukemia cell lines. Moreover, the encapsulated drug overcame the resistance of the K562 cell line to the drug. LN containing ET might be used as a promising drug delivery system in leukemia due to their capacity to overcome the in vivo pitfalls of the free drug and their efficacy in vitro in leukemia cell lines.

Stable isotope-assisted metabolomics to detect metabolic flux changes in mammalian cell cultures

The determination of metabolic fluxes provides detailed information of cellular physiology, and the assessment of metabolic flux changes upon a certain perturbation can help to improve biotechnological and pharmaceutical processes. Stable isotope-assisted metabolomics using tracer-labeled substrates is the method of choice to determine the fluxes. Though well-established for microbial cultures, the application to mammalian cells is generally complex and still limited. However, there have been great achievements in recent years and it is now emerging that stable isotope-assisted metabolic flux analysis in mammalian cell cultures will help improving biotechnological production and will also support drug development and discovery.

Potential metabolomic biomarkers for evaluation of adriamycin efficacy using a urinary 1H-NMR spectroscopy

A metabolomics approach using proton nuclear magnetic resonance (NMR) was applied to investigate metabolic alterations following adriamycin (ADR) treatment for gastric adenocarcinoma. After BALB/c-nu/nu mice were implanted with human gastric adenocarcinoma, ADR (1 or 3 mg kg(-1) per day) was intraperitoneally administered for 5 days. Urine was collected on days 2 and 5 and analyzed by NMR. The levels of trimethylamine oxide (TMAO, ×0.3), hippurate (×0.3) and taurine (×0.6) decreased significantly (P < 0.05), whereas the levels of 3-indoxylsulfate (×12.6), trigonelline (×1.5), citrate (×2.5), trimethylamine (TMA, ×2.0) and 2-oxoglutarate (×2.3) increased significantly (P < 0.05) in the tumor model. After ADR treatment, TMAO, hippuarte and taurine were increased significantly on day 5 compared with those of the tumor model. The levels of 2-oxoglutarate, 3-indoxylsulfate, trigonelline, TMA and citrate, which increased in the tumor model, significantly decreased to those of normal control by ADR treatment. Furthermore, the ratio between TMA and TMAO was dramatically altered in both tumor and ADR-treated groups. Overall, metabolites such as TMAO, TMA, 3-indoxylsulfate, hippurate, trigonelline, citrate and 2-oxoglutarate related to the tricarboxylic acid (TCA) cycle might be considered as therapeutic targets to potentiate the efficacy of ADR. Thus, these results suggest that the metabolomics analysis of tumor response to ADR treatment may be applicable for demonstrating the efficacy of anticancer agent, ADR and treatment adaptation.

Stable isotope-resolved metabolomics and applications for drug development

Advances in analytical methodologies, principally nuclear magnetic resonance spectroscopy (NMR) and mass spectrometry (MS), during the last decade have made large-scale analysis of the human metabolome a reality. This is leading to the reawakening of the importance of metabolism in human diseases, particularly cancer. The metabolome is the functional readout of the genome, functional genome, and proteome; it is also an integral partner in molecular regulations for homeostasis. The interrogation of the metabolome, or metabolomics, is now being applied to numerous diseases, largely by metabolite profiling for biomarker discovery, but also in pharmacology and therapeutics. Recent advances in stable isotope tracer-based metabolomic approaches enable unambiguous tracking of individual atoms through compartmentalized metabolic networks directly in human subjects, which promises to decipher the complexity of the human metabolome at an unprecedented pace. This knowledge will revolutionize our understanding of complex human diseases, clinical diagnostics, as well as individualized therapeutics and drug response. In this review, we focus on the use of stable isotope tracers with metabolomics technologies for understanding metabolic network dynamics in both model systems and in clinical applications. Atom-resolved isotope tracing via the two major analytical platforms, NMR and MS, has the power to determine novel metabolic reprogramming in diseases, discover new drug targets, and facilitates ADME studies. We also illustrate new metabolic tracer-based imaging technologies, which enable direct visualization of metabolic processes in vivo. We further outline current practices and future requirements for biochemoinformatics development, which is an integral part of translating stable isotope-resolved metabolomics into clinical reality.

Compartmentation of glycogen metabolism revealed from 13C isotopologue distributions

Background

Stable isotope tracers are used to assess metabolic flux profiles in living cells. The existing methods of measurement average out the isotopic isomer distribution in metabolites throughout the cell, whereas the knowledge of compartmental organization of analyzed pathways is crucial for the evaluation of true fluxes. That is why we accepted a challenge to create a software tool that allows deciphering the compartmentation of metabolites based on the analysis of average isotopic isomer distribution.

Results

The software Isodyn, which simulates the dynamics of isotopic isomer distribution in central metabolic pathways, was supplemented by algorithms facilitating the transition between various analyzed metabolic schemes, and by the tools for model discrimination. It simulated ¹³C isotope distributions in glucose, lactate, glutamate and glycogen, measured by mass spectrometry after incubation of hepatocytes in the presence of only labeled glucose or glucose and lactate together (with label either in glucose or lactate). The simulations assumed either a single intracellular hexose phosphate pool, or also channeling of hexose phosphates resulting in a different isotopic composition of glycogen. Model discrimination test was applied to check the consistency of both models with experimental data. Metabolic flux profiles, evaluated with the accepted model that assumes channeling, revealed the range of changes in metabolic fluxes in liver cells.

Conclusions

The analysis of compartmentation of metabolic networks based on the measured ¹³C distribution was included in Isodyn as a routine procedure. The advantage of this implementation is that, being a part of evaluation of metabolic fluxes, it does not require additional experiments to study metabolic compartmentation. The analysis of experimental data revealed that the distribution of measured ¹³C-labeled glucose metabolites is inconsistent with the idea of perfect mixing of hexose phosphates in cytosol. In contrast, the observed distribution indicates the presence of a separate pool of hexose phosphates that is channeled towards glycogen synthesis.