Cyclin-dependent kinases CDK4 and CDK6 are essential for the control of the cell cycle through the G1 phase. Aberrant expression of CDK4 and CDK6 is a hallmark of cancer, which would suggest that CDK4 and CDK6 are attractive targets for cancer therapy. Herein, we report that calcein AM (the calcein acetoxymethyl-ester) is a potent specific inhibitor of CDK4 and CDK6 in HCT116 human colon adenocarcinoma cells, inhibiting retinoblastoma protein (pRb) phosphorylation and inducing cell cycle arrest in the G1 phase. The metabolic effects of calcein AM on HCT116 cells were also evaluated and the flux between the oxidative and non-oxidative branches of the pentose phosphate pathway was significantly altered. To elucidate whether these metabolic changes were due to the inhibition of CDK4 and CDK6, we also characterized the metabolic profile of a CDK4, CDK6 and CDK2 triple knockout of mouse embryonic fibroblasts. The results show that the metabolic profile associated with the depletion of CDK4, CDK6 and CDK2 coincides with the metabolic changes induced by calcein AM on HCT116 cells, thus confirming that the inhibition of CDK4 and CDK6 disrupts the balance between the oxidative and non-oxidative branches of the pentose phosphate pathway. Taken together, these results indicate that low doses of calcein can halt cell division and kill tumor cells. Thus, selective inhibition of CDK4 and CDK6 may be of greater pharmacological interest, since inhibitors of these kinases affect both cell cycle progression and the robust metabolic profile of tumors.
Cyclin-dependent kinases; CDK-inhibitor; Tracer-based metabolomics; Pentose phosphate pathway; Phase-plane analysis
Chemoprevention is a pragmatic approach to reduce the risk of colorectal cancer, one of the leading causes of cancer-related death in western countries. In this regard, maslinic acid (MA), a pentacyclic triterpene extracted from wax-like coatings of olives, is known to inhibit proliferation and induce apoptosis in colon cancer cell lines without affecting normal intestinal cells. The present study evaluated the chemopreventive efficacy and associated mechanisms of maslinic acid treatment on spontaneous intestinal tumorigenesis in ApcMin/+ mice. Twenty-two mice were randomized into 2 groups: control group and MA group, fed with a maslinic acid–supplemented diet for six weeks. MA treatment reduced total intestinal polyp formation by 45% (P<0.01). Putative molecular mechanisms associated with suppressing intestinal polyposis in ApcMin/+ mice were investigated by comparing microarray expression profiles of MA-treated and control mice and by analyzing the serum metabolic profile using NMR techniques. The different expression phenotype induced by MA suggested that it exerts its chemopreventive action mainly by inhibiting cell-survival signaling and inflammation. These changes eventually induce G1-phase cell cycle arrest and apoptosis. Moreover, the metabolic changes induced by MA treatment were associated with a protective profile against intestinal tumorigenesis. These results show the efficacy and underlying mechanisms of MA against intestinal tumor development in the ApcMin/+ mice model, suggesting its chemopreventive potential against colorectal cancer.
The mitochondrial electron transport chain transforms energy satisfying cellular demand and generates reactive oxygen species (ROS) that act as metabolic signals or destructive factors. Therefore, knowledge of the possible modes and bifurcations of electron transport that affect ROS signaling provides insight into the interrelationship of mitochondrial respiration with cellular metabolism. Here, a bifurcation analysis of a sequence of the electron transport chain models of increasing complexity was used to analyze the contribution of individual components to the modes of respiratory chain behavior. Our algorithm constructed models as large systems of ordinary differential equations describing the time evolution of the distribution of redox states of the respiratory complexes. The most complete model of the respiratory chain and linked metabolic reactions predicted that condensed mitochondria produce more ROS at low succinate concentration and less ROS at high succinate levels than swelled mitochondria. This prediction was validated by measuring ROS production under various swelling conditions. A numerical bifurcation analysis revealed qualitatively different types of multistationary behavior and sustained oscillations in the parameter space near a region that was previously found to describe the behavior of isolated mitochondria. The oscillations in transmembrane potential and ROS generation, observed in living cells were reproduced in the model that includes interaction of respiratory complexes with the reactions of TCA cycle. Whereas multistationarity is an internal characteristic of the respiratory chain, the functional link of respiration with central metabolism creates oscillations, which can be understood as a means of auto-regulation of cell metabolism.
The mitochondrial respiratory chain shows a variety of modes of behavior. In living cells, flashes of ROS production and oscillations accompanied by a decrease of transmembrane potential can be registered. The mechanisms of such complex behavior are difficult to rationalize without a mathematical formalization of mitochondrial respiration. Our most complete model of mitochondrial respiration accounts for the details of electron transport, reproducing the observed types of behavior, which includes the existence of multiple steady states and periodic oscillations. This most detailed model contains hundreds of differential equations, and such complexity makes it difficult to grasp the main determinants of its behavior. Therefore the full model was reduced to a simplified description of complex III only, and numerical bifurcation analysis was used to study its behavior. Then the evolution of its behavior was traced in a sequence of models with increasing complexity leading back to the full model. This analysis revealed the mechanism of switching between the modes of behavior and the conditions for persistence in a given state, which defines ATP production, ROS signaling and destructive effects. This is important for understanding the biochemical basics of many systemic diseases.
Transketolase is an enzyme involved in a critical step of the non-oxidative branch of the pentose phosphate pathway whose inhibition could lead to new anticancer drugs. Here, we report new human transketolase inhibitors, based on the phenyl urea scaffold, found by applying structure-based virtual screening. These inhibitors are designed to cover a hot spot in the dimerization interface of the homodimer of the enzyme, providing for the first time compounds with a suggested novel binding mode not based on mimicking the thiamine pyrophosphate cofactor.
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.
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 13C 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.
The analysis of compartmentation of metabolic networks based on the measured 13C 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 13C-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.
Transketolase-like 1 (TKTL1) induces glucose degradation through anaerobic pathways, even in presence of oxygen, favoring the malignant aerobic glycolytic phenotype characteristic of tumor cells. As TKTL1 appears to be a valid biomarker for cancer prognosis, the aim of the current study was to correlate its expression with tumor stage, probability of tumor recurrence and survival, in a series of colorectal cancer patients.
Tumor tissues from 63 patients diagnosed with colorectal cancer at different stages of progression were analyzed for TKTL1 by immunohistochemistry. Staining was quantified by computational image analysis, and correlations between enzyme expression, local growth, lymph-node involvement and metastasis were assessed. The highest values for TKTL1 expression were detected in the group of stage III tumors, which showed significant differences from the other groups (Kruskal-Wallis test, P = 0.000008). Deeper analyses of T, N and M classifications revealed a weak correlation between local tumor growth and enzyme expression (Mann-Whitney test, P = 0.029), a significant association of the enzyme expression with lymph-node involvement (Mann-Whitney test, P = 0.0014) and a significant decrease in TKTL1 expression associated with metastasis (Mann-Whitney test, P = 0.0004).
To our knowledge, few studies have explored the association between variations in TKTL1 expression in the primary tumor and metastasis formation. Here we report downregulation of enzyme expression when metastasis appears, and a correlation between enzyme expression and regional lymph-node involvement in colon cancer. This finding may improve our understanding of metastasis and lead to new and more efficient therapies against cancer.
Chronic Obstructive Pulmonary Disease (COPD) is an inflammatory process of the lung inducing persistent airflow limitation. Extensive systemic effects, such as skeletal muscle dysfunction, often characterize these patients and severely limit life expectancy. Despite considerable research efforts, the molecular basis of muscle degeneration in COPD is still a matter of intense debate. In this study, we have applied a network biology approach to model the relationship between muscle molecular and physiological response to training and systemic inflammatory mediators. Our model shows that failure to co-ordinately activate expression of several tissue remodelling and bioenergetics pathways is a specific landmark of COPD diseased muscles. Our findings also suggest that this phenomenon may be linked to an abnormal expression of a number of histone modifiers, which we discovered correlate with oxygen utilization. These observations raised the interesting possibility that cell hypoxia may be a key factor driving skeletal muscle degeneration in COPD patients.
Chronic Obstructive Pulmonary Disease (COPD) is a major life threatening disease of the lungs, characterized by airflow limitation and chronic inflammation. Progressive reduction of the body muscle mass is a condition linked to COPD that significantly decreases quality of life and survival. Physical exercise has been proposed as a therapeutic option but its utility is still a matter of debate. The mechanisms underlying muscle wasting are also still largely unknown. The results presented in this paper show that diseased muscles are largely unable to coordinate the expression of muscle remodelling and bioenergetics pathways and that the cause of this phenomena may be tissue hypoxia. These findings contrast with current hypotheses based on the role of chronic inflammation and show that a mechanism based on an oxygen driven, epigenetic control of these two important functions may be an important disease mechanism.
Phenazine 5,10-dioxides (PDOs) are a new class of bioreductive cytotoxins, which could act towards tumours containing hypoxic regions. The PDOs selective-hypoxic bioreduction was probed in vitro; however, the mechanism of action has not been completely explained. Besides, PDOs in vivo antitumour activities have not been demonstrated hitherto. We study the mechanism of hypoxic/normoxic cytotoxicity of PDO representative members. Electron spin resonance is used to confirm •OH production, alkaline comet assay to determine genotoxicity, and gel electrophoresis and flow cytometry to analyze DNA fragmentation and cell cycle distribution. Chemically induced rat breast tumours are employed to evaluate in vivo activities. For the most selective cytotoxin, 7(8)-bromo-2-hydroxyphenazine 5,10-dioxide (PDO1), exclusive hypoxic •OH production is evidenced, while for the unselective ones, •OH is produced in both conditions (normoxia and simulated hypoxia). In normoxia (Caco-2 cells), PDO1 induces cell-cycle arrest and DNA fragmentation but does not significantly induce apoptosis neither at IC50 nor IC80. No difference in the comet-assay scores are observed in normoxia and simulated hypoxia being the unselective 2-amino-7(8)-bromophenazine 5,10-dioxide (PDO2) the most genotoxic. The in vivo efficacy with the absence of systemic toxicity of PDO1 and PDO2 is checked out. Results from this study highlight the potential of PDOs as new therapeutics for cancer.
Maslinic acid, a pentacyclic triterpene found in the protective wax-like coating of the leaves and fruit of Olea europaea L., is a promising agent for the prevention of colon cancer. We have shown elsewhere that maslinic acid inhibits cell proliferation to a significant extent and activates mitochondrial apoptosis in colon cancer cells. In our latest work we have investigated further this compound's apoptotic molecular mechanism.
We used HT29 adenocarcinoma cells. Changes genotoxicity were analyzed by single-cell gel electrophoresis (comet assay). The cell cycle was determined by flow cytometry. Finally, changes in protein expression were examined by western blotting. Student's t-test was used for statistical comparison.
HT29 cells treated with maslinic acid showed significant increases in genotoxicity and cell-cycle arrest during the G0/G1 phase after 72 hours' treatment and an apoptotic sub-G0/G1 peak after 96 hours. Nevertheless, the molecular mechanism for this cytotoxic effect of maslinic acid has never been properly explored. We show here that the anti-tumoral activity of maslinic acid might proceed via p53-mediated apoptosis by acting upon the main signaling components that lead to an increase in p53 activity and the induction of the rest of the factors that participate in the apoptotic pathway. We found that in HT29 cells maslinic acid activated the expression of c-Jun NH2-terminal kinase (JNK), thus inducing p53. Treatment of tumor cells with maslinic acid also resulted in an increase in the expression of Bid and Bax, repression of Bcl-2, release of cytochrome-c and an increase in the expression of caspases -9, -3, and -7. Moreover, maslinic acid produced belated caspase-8 activity, thus amplifying the initial mitochondrial apoptotic signaling.
All these results suggest that maslinic acid induces apoptosis in human HT29 colon-cancer cells through the JNK-Bid-mediated mitochondrial apoptotic pathway via the activation of p53. Thus we propose a plausible sequential molecular mechanism for the expression of the different proteins responsible for the intrinsic mitochondrial apoptotic pathway. Further studies with other cell lines will be needed to confirm the general nature of these findings.
colon-cancer cells; JNK; maslinic acid; mitochondrial apoptotic pathway; p53-mediated apoptosis; triterpenes
Reactive oxygen species (ROS) produced in the mitochondrial respiratory chain (RC) are primary signals that modulate cellular adaptation to environment, and are also destructive factors that damage cells under the conditions of hypoxia/reoxygenation relevant for various systemic diseases or transplantation. The important role of ROS in cell survival requires detailed investigation of mechanism and determinants of ROS production. To perform such an investigation we extended our rule-based model of complex III in order to account for electron transport in the whole RC coupled to proton translocation, transmembrane electrochemical potential generation, TCA cycle reactions, and substrate transport to mitochondria. It fits respiratory electron fluxes measured in rat brain mitochondria fueled by succinate or pyruvate and malate, and the dynamics of NAD+ reduction by reverse electron transport from succinate through complex I. The fitting of measured characteristics gave an insight into the mechanism of underlying processes governing the formation of free radicals that can transfer an unpaired electron to oxygen-producing superoxide and thus can initiate the generation of ROS. Our analysis revealed an association of ROS production with levels of specific radicals of individual electron transporters and their combinations in species of complexes I and III. It was found that the phenomenon of bistability, revealed previously as a property of complex III, remains valid for the whole RC. The conditions for switching to a state with a high content of free radicals in complex III were predicted based on theoretical analysis and were confirmed experimentally. These findings provide a new insight into the mechanisms of ROS production in RC.
Respiration at the level of mitochondria is considered as delivery of electrons and protons from NADH or succinate to oxygen through a set of transporters constituting the respiratory chain (RC). Mitochondrial respiration, dealing with transfer of unpaired electrons, may produce reactive oxygen species (ROS) such as O2− and subsequently H2O2 as side products. ROS are chemically very active and can cause oxidative damage to cellular components. The production of ROS, normally low, can increase under stress to the levels incompatible with cell survival; thus, understanding the ways of ROS production in the RC represents a vital task in research. We used mathematical modeling to analyze experiments with isolated brain mitochondria aimed to study relations between electron transport and ROS production. Elsewhere we reported that mitochondrial complex III can operate in two distinct steady states at the same microenvironmental conditions, producing either low or high levels of ROS. Here, this property of bistability was confirmed for the whole RC. The associations between measured ROS production and computed individual free radical levels in complexes I and III were established. The discovered phenomenon of bistability is important as a basis for new strategies in organ transplantation and therapy.
To enhance our understanding of complex biological systems like diseases we need to put all of the available data into context and use this to detect relations, pattern and rules which allow predictive hypotheses to be defined. Life science has become a data rich science with information about the behaviour of millions of entities like genes, chemical compounds, diseases, cell types and organs, which are organised in many different databases and/or spread throughout the literature. Existing knowledge such as genotype - phenotype relations or signal transduction pathways must be semantically integrated and dynamically organised into structured networks that are connected with clinical and experimental data. Different approaches to this challenge exist but so far none has proven entirely satisfactory.
To address this challenge we previously developed a generic knowledge management framework, BioXM™, which allows the dynamic, graphic generation of domain specific knowledge representation models based on specific objects and their relations supporting annotations and ontologies. Here we demonstrate the utility of BioXM for knowledge management in systems biology as part of the EU FP6 BioBridge project on translational approaches to chronic diseases. From clinical and experimental data, text-mining results and public databases we generate a chronic obstructive pulmonary disease (COPD) knowledge base and demonstrate its use by mining specific molecular networks together with integrated clinical and experimental data.
We generate the first semantically integrated COPD specific public knowledge base and find that for the integration of clinical and experimental data with pre-existing knowledge the configuration based set-up enabled by BioXM reduced implementation time and effort for the knowledge base compared to similar systems implemented as classical software development projects. The knowledgebase enables the retrieval of sub-networks including protein-protein interaction, pathway, gene - disease and gene - compound data which are used for subsequent data analysis, modelling and simulation. Pre-structured queries and reports enhance usability; establishing their use in everyday clinical settings requires further simplification with a browser based interface which is currently under development.
Metabolic flux profiling based on the analysis of distribution of stable isotope tracer in metabolites is an important method widely used in cancer research to understand the regulation of cell metabolism and elaborate new therapeutic strategies. Recently, we developed software Isodyn, which extends the methodology of kinetic modeling to the analysis of isotopic isomer distribution for the evaluation of cellular metabolic flux profile under relevant conditions. This tool can be applied to reveal the metabolic effect of proapoptotic drug edelfosine in leukemia Jurkat cell line, uncovering the mechanisms of induction of apoptosis in cancer cells.
The study of 13C distribution of Jukat cells exposed to low edelfosine concentration, which induces apoptosis in ≤5% of cells, revealed metabolic changes previous to the development of apoptotic program. Specifically, it was found that low dose of edelfosine stimulates the TCA cycle. These metabolic perturbations were coupled with an increase of nucleic acid synthesis de novo, which indicates acceleration of biosynthetic and reparative processes. The further increase of the TCA cycle fluxes, when higher doses of drug applied, eventually enhance reactive oxygen species (ROS) production and trigger apoptotic program.
The application of Isodyn to the analysis of mechanism of edelfosine-induced apoptosis revealed primary drug-induced metabolic changes, which are important for the subsequent initiation of apoptotic program. Initiation of such metabolic changes could be exploited in anticancer therapy.
Cell differentiation is an orderly process that begins with modifications in gene expression. This process is regulated by the acetylation state of histones. Removal of the acetyl groups of histones by specific enzymes (histone deacetylases, HDAC) usually downregulates expression of genes that can cause cells to differentiate, and pharmacological inhibitors of these enzymes have been shown to induce differentiation in several colon cancer cell lines. Butyrate at high (mM) concentration is both a precursor for acetyl-CoA and a known HDAC inhibitor that induces cell differentiation in colon cells. The dual role of butyrate raises the question whether its effects on HT29 cell differentiation are due to butyrate metabolism or to its HDAC inhibitor activity. To distinguish between these two possibilities, we used a tracer-based metabolomics approach to compare the metabolic changes induced by two different types of HDAC inhibitors (butyrate and the non-metabolic agent trichostatin A) and those induced by other acetyl-CoA precursors that do not inhibit HDAC (caprylic and capric acids). [1,2-13C2]-d-glucose was used as a tracer and its redistribution among metabolic intermediates was measured to estimate the contribution of glycolysis, the pentose phosphate pathway and the Krebs cycle to the metabolic profile of HT29 cells under the different treatments. The results demonstrate that both HDAC inhibitors (trichostatin A and butyrate) induce a common metabolic profile that is associated with histone deacetylase inhibition and differentiation of HT29 cells whereas the metabolic effects of acetyl-CoA precursors are different from those of butyrate. The experimental findings support the concept of crosstalk between metabolic and cell signalling events, and provide an experimental approach for the rational design of new combined therapies that exploit the potential synergism between metabolic adaptation and cell differentiation processes through modification of HDAC activity.
Butyrate; Cell differentiation; Histone deacetylase; Medium chain fatty acids; Metabolism; Trichostatin A
Increased production of reactive oxygen species (ROS) in mitochondria underlies major systemic diseases, and this clinical problem stimulates a great scientific interest in the mechanism of ROS generation. However, the mechanism of hypoxia-induced change in ROS production is not fully understood. To mathematically analyze this mechanism in details, taking into consideration all the possible redox states formed in the process of electron transport, even for respiratory complex III, a system of hundreds of differential equations must be constructed. Aimed to facilitate such tasks, we developed a new methodology of modeling, which resides in the automated construction of large sets of differential equations. The detailed modeling of electron transport in mitochondria allowed for the identification of two steady state modes of operation (bistability) of respiratory complex III at the same microenvironmental conditions. Various perturbations could induce the transition of respiratory chain from one steady state to another. While normally complex III is in a low ROS producing mode, temporal anoxia could switch it to a high ROS producing state, which persists after the return to normal oxygen supply. This prediction, which we qualitatively validated experimentally, explains the mechanism of anoxia-induced cell damage. Recognition of bistability of complex III operation may enable novel therapeutic strategies for oxidative stress and our method of modeling could be widely used in systems biology studies.
The levels of reactive oxygen species (ROS) that are generated as a side product of mitochondrial respiratory electron transport largely define the extent of oxidative stress in living cells. Free radicals formed in electron transport, such as ubisemiquinone, could pass their non-paired electron directly to oxygen, thus producing superoxide radical that gives rise to a variety of ROS. It is well known in clinical practice that upon recommencing oxygen supply after anoxia a tissue produces much more ROS than before the anoxia, and the state of high ROS production is stable. The mechanism of switching from low to high ROS production by temporal anoxia was unknown, in part because of the lack of detailed mathematical description of hundreds of redox states of respiratory complexes, which are formed in the process of electron transport. A new methodology of automated construction of large systems of differential equations allowed us to describe the system in detail and predicts that the mechanism of paradoxical effect of anoxia-reoxygenation could be defined by the properties of complex III of mitochondrial respiratory chain. Our experiments confirmed that the effect of hypoxia-reoxygenation is confined by intramitochondrial processes since it is observed in isolated mitochondria.
Reactive oxygen species (ROS) generation in mitochondria as a side product
of electron and proton transport through the inner membrane is important for
normal cell operation as well as development of pathology. Matrix and cytosol
alkalization stabilizes semiquinone radical, a potential superoxide producer,
and we hypothesized that proton deficiency under the excess of electron donors
enhances reactive oxygen species generation. We tested this hypothesis by
measuring pH dependence of reactive oxygen species released by mitochondria.
The experiments were performed in the media with pH varying from 6 to 8 in the
presence of complex II substrate succinate or under more physiological
conditions with complex I substrates glutamate and malate. Matrix pH was
manipulated by inorganic phosphate, nigericine, and low concentrations of
uncoupler or valinomycin. We found that high pH strongly increased the rate of
free radical generation in all of the conditions studied, even when ΔpH
= 0 in the presence of nigericin. In the absence of inorganic phosphate, when
the matrix was the most alkaline, pH shift in the medium above 7 induced
permeability transition accompanied by the decrease of ROS production. ROS
production increase induced by the alkalization of medium was observed with
intact respiring mitochondria as well as in the presence of complex I
inhibitor rotenone, which enhanced reactive oxygen species release. The
phenomena revealed in this report are important for understanding mechanisms
governing mitochondrial production of reactive oxygen species, in particular
that related with uncoupling proteins.
A current trend in neuroscience research is the use of stable isotope tracers in order to address metabolic processes in vivo. The tracers produce a huge number of metabolite forms that differ according to the number and position of labeled isotopes in the carbon skeleton (isotopomers) and such a large variety makes the analysis of isotopomer data highly complex. On the other hand, this multiplicity of forms does provide sufficient information to address cell operation in vivo. By the end of last millennium, a number of tools have been developed for estimation of metabolic flux profile from any possible isotopomer distribution data. However, although well elaborated, these tools were limited to steady state analysis, and the obtained set of fluxes remained disconnected from their biochemical context. In this review we focus on a new numerical analytical approach that integrates kinetic and metabolic flux analysis. The related computational algorithm estimates the dynamic flux based on the time-dependent distribution of all possible isotopomers of metabolic pathway intermediates that are generated from a labeled substrate. The new algorithm connects specific tracer data with enzyme kinetic characteristics, thereby extending the amount of data available for analysis: it uses enzyme kinetic data to estimate the flux profile, and vice versa, for the kinetic analysis it uses in vivo tracer data to reveal the biochemical basis of the estimated metabolic fluxes.