Metabolic control analysis of respiration in human cancer tissue

Metabolic control analysis as a strategy to identify therapeutic targets, the case of cancer glycolysis

The use of kinetic modeling and metabolic control analysis (MCA) to identify possible therapeutic targets and to investigate the controlling and regulatory mechanisms in cancer glycolysis is here reviewed. The glycolytic pathway has been considered a target to decrease cancer cell growth; however, its occurrence in normal cells makes it difficult to design therapeutic strategies that target this pathway in pathological cells. Notwithstanding, the over-expression of all enzymes and transporters, as well as the expression of isoenzymes with different kinetic and regulatory properties in cancer cells, suggested a different distribution of the control of glycolytic flux than that observed in normal cells. Kinetic models of glycolysis are constructed with enzyme kinetics experimental data, validated with the steady-state metabolite concentrations and glycolytic fluxes; applying MCA, permitted us to identify the steps with the highest control of glycolysis in cancer cells, but low control in normal cells. The cancer glycolysis main controlling steps under several metabolic conditions were: glucose transport, hexokinase and hexose-6-phosphate isomerase (HPI); whereas in normal cells were: the first two and phosphofructokinase-1. HPI is the best therapeutic target because it exerts high control in cancer glycolytic flux, but not in normal cells. Furthermore, kinetic modeling also contributed to identifying new feed-back and feed-forward regulatory loops in cancer cells glycolysis, and to understanding the mode of metabolic action of glycolytic inhibitors. Thus, MCA and metabolic modeling allowed to propose new strategies for inhibiting glycolysis in cancer cells.

Physiological Cell Culture Media Tune Mitochondrial Bioenergetics and Drug Sensitivity in Cancer Cell Models

Simple Summary

Cell biologists trust in standard media for analyzing cellular functions and for the specification of target-oriented drug efficacies in cell culture settings. Here, we present a general applicable workflow for the constant monitoring of bioenergetic states of cells grown in 2D cell models to accompany tailored drug discovery efforts. Using in-depth high-resolution respirometry analyses (HRR) of mitochondrial function, we unveiled that the human-plasma-like media (HPLM) altered cellular energetic states. In a systematic HRR setup for drug profiling experiments, we revealed an unexpected side effect of an FDA-approved cancer drug on mitochondrial function, exclusively in HPLM. Thus, we believe that both the recordings of bioenergetic states and the use of more physiological media would improve and reshape cell-based drug discovery ventures.

Abstract

Two-dimensional cell cultures are established models in research for studying and perturbing cell-type specific functions. However, many limitations apply to the cell growth in a monolayer using standard cell culture media. Although they have been used for decades, their formulations do not mimic the composition of the human cell environment. In this study, we analyzed the impact of a newly formulated human plasma-like media (HPLM) on cell proliferation, mitochondrial bioenergetics, and alterations of drug efficacies using three distinct cancer cell lines. Using high-resolution respirometry, we observed that cells grown in HPLM displayed significantly altered mitochondrial bioenergetic profiles, particularly related to mitochondrial density and mild uncoupling of respiration. Furthermore, in contrast to standard media, the growth of cells in HPLM unveiled mitochondrial dysfunction upon exposure to the FDA-approved kinase inhibitor sunitinib. This seemingly context-dependent side effect of this drug highlights that the selection of the cell culture medium influences the assessment of cancer drug sensitivities. Thus, we suggest to prioritize media with a more physiological composition for analyzing bioenergetic profiles and to take it into account for assigning drug efficacies in the cell culture model of choice.

Active mitochondrial respiration in cancer: a target for the drug

The relative contribution of mitochondrial respiration and subsequent energy production in malignant cells has remained controversial to date. Enhanced aerobic glycolysis and impaired mitochondrial respiration have gained more attention in the metabolic study of cancer. In contrast to the popular concept, mitochondria of cancer cells oxidize a diverse array of metabolic fuels to generate a majority of the cellular energy by respiration. Several mitochondrial respiratory chain (MRC) subunits' expressions are critical for the growth, metastasis, and cancer cell invasion. Also, the assembly factors, which regulate the integration of individual MRC complexes into native super-complexes, are upregulated in cancer. Moreover, a series of anti-cancer drugs function by inhibiting respiration and ATP production. In this review, we have specified the roles of mitochondrial fuels, MRC subunits, and super-complex assembly factors that promote active respiration across different cancer types and discussed the potential roles of MRC inhibitor drugs in controlling cancer.

Molecular and Supramolecular Structure of the Mitochondrial Oxidative Phosphorylation System: Implications for Pathology

Under aerobic conditions, mitochondrial oxidative phosphorylation (OXPHOS) converts the energy released by nutrient oxidation into ATP, the currency of living organisms. The whole biochemical machinery is hosted by the inner mitochondrial membrane (mtIM) where the protonmotive force built by respiratory complexes, dynamically assembled as super-complexes, allows the F1FO-ATP synthase to make ATP from ADP + Pi. Recently mitochondria emerged not only as cell powerhouses, but also as signaling hubs by way of reactive oxygen species (ROS) production. However, when ROS removal systems and/or OXPHOS constituents are defective, the physiological ROS generation can cause ROS imbalance and oxidative stress, which in turn damages cell components. Moreover, the morphology of mitochondria rules cell fate and the formation of the mitochondrial permeability transition pore in the mtIM, which, most likely with the F1FO-ATP synthase contribution, permeabilizes mitochondria and leads to cell death. As the multiple mitochondrial functions are mutually interconnected, changes in protein composition by mutations or in supercomplex assembly and/or in membrane structures often generate a dysfunctional cascade and lead to life-incompatible diseases or severe syndromes. The known structural/functional changes in mitochondrial proteins and structures, which impact mitochondrial bioenergetics because of an impaired or defective energy transduction system, here reviewed, constitute the main biochemical damage in a variety of genetic and age-related diseases.

Mitochondrial Respiration in KRAS and BRAF Mutated Colorectal Tumors and Polyps

This study aimed to characterize the ATP-synthesis by oxidative phosphorylation in colorectal cancer (CRC) and premalignant colon polyps in relation to molecular biomarkers KRAS and BRAF. This prospective study included 48 patients. Resected colorectal polyps and postoperative CRC tissue with adjacent normal tissue (control) were collected. Patients with polyps and CRC were divided into three molecular groups: KRAS mutated, BRAF mutated and KRAS/BRAF wild-type. Mitochondrial respiration in permeabilized tissue samples was observed using high resolution respirometry. ADP-activated respiration rate (V_max) and an apparent affinity of mitochondria to ADP, which is related to mitochondrial outer membrane (MOM) permeability, were determined. Clear differences were present between molecular groups. KRAS mutated CRC group had lower V_max values compared to wild-type; however, the V_max value was higher than in the control group, while MOM permeability did not change. This suggests that KRAS mutation status might be involved in acquiring oxidative phenotype. KRAS mutated polyps had higher V_max values and elevated MOM permeability as compared to the control. BRAF mutated CRC and polyps had reduced respiration and altered MOM permeability, indicating a glycolytic phenotype. To conclude, prognostic biomarkers KRAS and BRAF are likely related to the metabolic phenotype in CRC and polyps. Assessment of the tumor mitochondrial ATP synthesis could be a potential component of patient risk stratification.

Tubulin βII and βIII Isoforms as the Regulators of VDAC Channel Permeability in Health and Disease

In recent decades, there have been several models describing the relationships between the cytoskeleton and the bioenergetic function of the cell. The main player in these models is the voltage-dependent anion channel (VDAC), located in the mitochondrial outer membrane. Most metabolites including respiratory substrates, ADP, and Pi enter mitochondria only through VDAC. At the same time, high-energy phosphates are channeled out and directed to cellular energy transfer networks. Regulation of these energy fluxes is controlled by β-tubulin, bound to VDAC. It is also thought that β-tubulin‒VDAC interaction modulates cellular energy metabolism in cancer, e.g., switching from oxidative phosphorylation to glycolysis. In this review we focus on the described roles of unpolymerized αβ-tubulin heterodimers in regulating VDAC permeability for adenine nucleotides and cellular bioenergetics. We introduce the Mitochondrial Interactosome model and the function of the βII-tubulin subunit in this model in muscle cells and brain synaptosomes, and also consider the role of βIII-tubulin in cancer cells.

Identification of a co-target for enhancing efficacy of sorafenib in HCC through a quantitative modeling approach

Sorafenib (SFB), a multi-kinase inhibitor, is the only approved drug for treating hepatocellular carcinoma (HCC). However, SFB shows low efficacy in many cases. HCC related mortality therefore remains to be high worldwide. SFB, a multi-kinase inhibitor is also known to modulate the redox homeostasis in cancer cells. To understand the effect of SFB on the redox status, a quantitative understanding of the system is necessary. Kinetic modeling of the relevant pathways is a useful approach for obtaining a quantitative understanding of the pathway dynamics and to rank the individual factors based on the extent of influence they wield on the pathway. Here, we report a comprehensive model of the glutathione reaction network (GSH_net ), consisting of four modules and includes SFB-induced redox stress. We compared GSH_net simulations for HCC of six different etiologies with healthy liver, and correctly identified the expected variations in cancer. Next, we studied alterations induced in the system upon SFB treatment and observed differential H₂ O₂ dynamics in all the conditions. Using metabolic control analysis, we identified glutathione S-transferase (GST) as the enzyme with the highest selective control coefficient, making it an attractive co-target for potentiating the action of SFB across all six etiologies. As a proof-of-concept, we selected ethacrynic acid (EA), a known inhibitor of GST, and verified ex vivo that EA synergistically potentiates the cytotoxic effect of SFB. Being an FDA approved drug, EA is a promising candidate for repurposing as a combination therapy with SFB for HCC treatment.

Comparative analysis of the bioenergetics of human adenocarcinoma Caco-2 cell line and postoperative tissue samples from colorectal cancer patients

The aim of this work was to explore the key bioenergetic properties for mitochondrial respiration in the widely-used Caco-2 cell line and in human colorectal cancer (HCC) postoperational tissue samples. Oxygraphy and metabolic control analysis (MCA) were applied to estimate the function of oxidative phosphorylation in cultured Caco-2 cells and HCC tissue samples. The mitochondria of Caco-2 cells and HCC tissues displayed larger functional activity of respiratory complex (C)II compared with CI, whereas in normal colon tissue an inverse pattern in the ratio of CI to CII activity was observed. MCA showed that the respiration in Caco-2 and HCC tissue cells is regulated by different parts of electron transport chain. In HCC tissues, this control is performed essentially at the level of respiratory chain complexes I-IV, whereas in Caco-2 cells at the level of CIV (cytochrome c oxidase) and the ATP synthasome. The differences we found in the regulation of respiratory chain activity and glycose index could represent an adaptive response to distinct growth conditions; this highlights the importance of proper validation of results obtained from in-vitro models before their extrapolation to the more complex in-vivo systems.

Mitochondrial Respiration in Human Colorectal and Breast Cancer Clinical Material Is Regulated Differently

We conducted quantitative cellular respiration analysis on samples taken from human breast cancer (HBC) and human colorectal cancer (HCC) patients. Respiratory capacity is not lost as a result of tumor formation and even though, functionally, complex I in HCC was found to be suppressed, it was not evident on the protein level. Additionally, metabolic control analysis was used to quantify the role of components of mitochondrial interactosome. The main rate-controlling steps in HBC are complex IV and adenine nucleotide transporter, but in HCC, complexes I and III. Our kinetic measurements confirmed previous studies that respiratory chain complexes I and III in HBC and HCC can be assembled into supercomplexes with a possible partial addition from the complex IV pool. Therefore, the kinetic method can be a useful addition in studying supercomplexes in cell lines or human samples. In addition, when results from culture cells were compared to those from clinical samples, clear differences were present, but we also detected two different types of mitochondria within clinical HBC samples, possibly linked to two-compartment metabolism. Taken together, our data show that mitochondrial respiration and regulation of mitochondrial membrane permeability have substantial differences between these two cancer types when compared to each other to their adjacent healthy tissue or to respective cell cultures.

Decoding Warburg's hypothesis: tumor-related mutations in the mitochondrial respiratory chain

Otto Warburg observed that cancer cells derived their energy from aerobic glycolysis by converting glucose to lactate. This mechanism is in opposition to the higher energy requirements of cancer cells because oxidative phosphorylation (OxPhos) produces more ATP from glucose. Warburg hypothesized that this phenomenon occurs due to the malfunction of mitochondria in cancer cells. The rediscovery of Warburg's hypothesis coincided with the discovery of mitochondrial tumor suppressor genes that may conform to Warburg's hypothesis along with the demonstrated negative impact of HIF-1 on PDH activity and the activation of HIF-1 by oncogenic signals such as activated AKT. This work summarizes the alterations in mitochondrial respiratory chain proteins that have been identified and their involvement in cancer. Also discussed is the fact that most of the mitochondrial mutations have been found in homoplasmy, indicating a positive selection during tumor evolution, thereby supporting their causal role.

Metabolic remodeling in human colorectal cancer and surrounding tissues: alterations in regulation of mitochondrial respiration and metabolic fluxes

The aim of the work was to evaluate whether or not there is glycolytic reprogramming in the neighboring cells of colorectal cancer (CRC). Using postoperative material we have compared the functional capacity of oxidative phosphorylation (OXPHOS) in CRC cells, their glycolytic activity and their inclination to aerobic glycolysis, with those of the surrounding and healthy colon tissue cells. Experiments showed that human CRC cannot be considered a hypoxic tumor, since the malignancy itself and cells surrounding it exhibited even higher rates of OXPHOS than healthy large intestine. The absence of acute hypoxia in colorectal carcinomas was also confirmed by their practically equal glucose-phosphorylating capacity as compared with surrounding non-tumorous tissue and by upregulation of VEGF family and their ligands. Studies indicated that human CRC cells in vivo exert a strong distant effect on the energy metabolism of neighboring cells, so that they acquire the bioenergetic parameters specific to the tumor itself. The growth of colorectal carcinomas was associated with potent downregulation of the creatine kinase system. As compared with healthy colon tissue, the tumor surrounding cells display upregulation of OXPHOS and have high values of basal and ADP activated respiration rates. Strong differences between the normal and CRC cells in the affinity of their mitochondria for ADP were revealed; the corresponding K_m values were measured as 93.6±7.7 µM for CRC cells and 84.9±9.9 µM for nearby tissue; both these apparent K_m (ADP) values were considerably (by almost 3 times) lower in comparison with healthy colon tissue cells (256±34 µM).

•

Human colorectal cancer is not a pure hypoxic tumor of the Warburg phenotype.

•

The total hexokinase activity of CRC cells is close to that in nearby tissues.

•

In the tumor there is overexpression of VEGFs (A, B, and C) and their receptors.

•

CRC has higher rates of OXPHOS as compared with healthy tissue cells.

•

Tumor-surrounding cells cannot fuel via a lactate shunt the growth of CRC cells.

A method aimed at assessing the functional consequences of the supramolecular organization of the respiratory electron transfer chain by time-resolved studies

A steadily increasing number of physiological, biochemical, and structural studies have provided a growing support to the notion that the respiratory electron transfer chain may contain supra-molecular edifices made of the assembly of some, if not all, of its individual links. This structure, usually referred to as the solid state model-in comparison to the liquid state model in which the electron transfer reactions between the membrane bound enzymes are diffusion controlled-is seen as conferring specific kinetic properties to the chain and thus as being highly relevant from a functional point of view. Although the assumption that structural changes are mirrored by functional adjustment is undoubtedly legitimate, experimental evidences supporting it remain scarce. Here we review a recent methodological development aimed at tackling the functional relevance of the supramolecular organization of the respiratory electron transfer chain in intact cells.

The role of tubulin in the mitochondrial metabolism and arrangement in muscle cells

Tubulin, a well-known component of the microtubule in the cytoskeleton, has an important role in the transport and positioning of mitochondria in a cell type dependent manner. This review describes different functional interactions of tubulin with cellular protein complexes and its functional interaction with the mitochondrial outer membrane. Tubulin is present in oxidative as well as glycolytic type muscle cells, but the kinetics of the in vivo regulation of mitochondrial respiration in these muscle types is drastically different. The interaction between VDAC and tubulin is probably influenced by such factors as isoformic patterns of VDAC and tubulin, post-translational modifications of tubulin and phosphorylation of VDAC. Important factor of the selective permeability of VDAC is the mitochondrial creatine kinase pathway which is present in oxidative cells, but is inactive or missing in glycolytic muscle and cancer cells. As the tubulin-VDAC interaction reduces the permeability of the channel by adenine nucleotides, energy transfer can then take place effectively only through the mitochondrial creatine kinase/phosphocreatine pathway. Therefore, closure of VDAC by tubulin may be one of the reasons of apoptosis in cells without the creatine kinase pathway. An important question in tubulin regulated interactions is whether other proteins are interacting with tubulin. The functional interaction may be direct, through other proteins like plectins, or influenced by simultaneous interaction of other complexes with VDAC.