Heinrich J. Huber

Mitochondrial and Plasma Membrane Potential of Cultured Cerebellar Neurons during Glutamate-Induced Necrosis, Apoptosis, and Tolerance

A failure of mitochondrial bioenergetics has been shown to be closely associated with the onset of apoptotic and necrotic neuronal injury. Here, we developed an automated computational model that interprets the single-cell fluorescence for tetramethylrhodamine methyl ester (TMRM) as a consequence of changes in either ΔΨm or ΔΨp, thus allowing for the characterization of responses for populations of single cells and subsequent statistical analysis. Necrotic injury triggered by prolonged glutamate excitation resulted in a rapid monophasic or biphasic loss of ΔΨm that was closely associated with a loss of ΔΨp and a rapid decrease in neuronal NADPH and ATP levels. Delayed apoptotic injury, induced by transient glutamate excitation, resulted in a small, reversible decrease in TMRM fluorescence, followed by a sustained hyperpolarization of ΔΨm as confirmed using the ΔΨp-sensitive anionic probe DiBAC2(3). This hyperpolarization of ΔΨm was closely associated with a significant increase in neuronal glucose uptake, NADPH availability, and ATP levels. Statistical analysis of the changes in ΔΨm or ΔΨp at a single-cell level revealed two major correlations; those neurons displaying a more pronounced depolarization of ΔΨp during the initial phase of glutamate excitation entered apoptosis more rapidly, and neurons that displayed a more pronounced hyperpolarization of ΔΨm after glutamate excitation survived longer. Indeed, those neurons that were tolerant to transient glutamate excitation (18%) showed the most significant increases in ΔΨm. Our results indicate that a hyperpolarization of ΔΨm is associated with increased glucose uptake, NADPH availability, and survival responses during excitotoxic injury.

Dynamics of outer mitochondrial membrane permeabilization during apoptosis.

Individual cells within a population undergo apoptosis at distinct, apparently random time points. By analyzing cellular mitotic history, we identified that sibling HeLa cell pairs, in contrast to random cell pairs, underwent apoptosis synchronously. This allowed us to use high-speed cellular imaging to investigate mitochondrial outer membrane permeabilization (MOMP), a highly coordinated, rapid process during apoptosis, at a temporal resolution approximately 100 times higher than possible previously. We obtained new functional and mechanistic insight into the process of MOMP: We were able to determine the kinetics of pore formation in the outer mitochondrial membrane from the initiation phase of cytochrome-c-GFP redistribution, and showed differential pore formation kinetics in response to intrinsic or extrinsic apoptotic stimuli (staurosporine, tumor necrosis factor-related apoptosis-inducing ligand (TRAIL)). We also detected that the onset of mitochondrial permeabilization frequently proceeded as a wave through the cytosol, and that the frequency of wave occurrence in response to TRAIL was reduced by inhibition of protein kinase CK2. Computational analysis by a partial differential equation model suggested that the spread of permeabilization signals could sufficiently be explained by diffusion–adsorption velocities of locally generated permeabilization inducers. Taken together, our study yielded the first comprehensive analysis of clonal cell-to-cell variability in apoptosis execution and allowed to visualize and explain the dynamics of MOMP in cells undergoing apoptosis.

Single-cell quantification of Bax activation and mathematical modelling suggest pore formation on minimal mitochondrial Bax accumulation

Mitochondrial outer membrane permeabilisation (MOMP) during apoptosis is triggered by the activation and oligomerisation of Bax and Bak, but a quantification of these processes in individual cells has not yet been performed. Single-cell imaging of Bax translocation and oligomerisation in Bax-deficient DU-145 cells expressing CFP-Bax and YFP-Bax revealed that both processes started only minutes before or concomitantly with MOMP, with the majority of Bax translocation and oligomerisation occurring downstream of MOMP. Quantification of YFP-Bax concentrations at mitochondria revealed an increase of only 1.8±1.5% at MOMP onset. This was increased to 11.2±3.6% in bak-silenced cells. These data suggested that Bax activation exceeded by far the quantities required for MOMP induction, and that minimal Bax or Bak activation may be sufficient to trigger rapid pore formation. In a cellular automaton modelling approach that incorporated the quantities and movement probabilities of Bax and its inhibitors, activators and enablers in the mitochondrial membrane, we could re-model rapid pore formation kinetics at submaximal Bax activation.

Systems Analysis of BCL2 Protein Family Interactions Establishes a Model to Predict Responses to Chemotherapy

Apoptotic desensitization is a hallmark of cancer cells, but present knowledge of molecular systems controlling apoptosis has yet to provide significant prognostic insights. Here, we report findings from a systems study of the intrinsic pathway of apoptosis by BCL2 family proteins and clinical translation of its findings into a model with applications in colorectal cancer (CRC). By determining absolute protein quantifications in CRC cells and patient tumor samples, we found that BAK and BAX were expressed more highly than their antiapoptotic inhibitors. This counterintuitive finding suggested that sole inhibition of effector BAX and BAK could not be sufficient for systems stability in nonstressed cells. Assuming a model of direct effector activation by BH3-only proteins, we calculated that the amount of stress-induced BH3-only proteins required to activate mitochondrial apoptosis could predict individual death responses of CRC cells to 5-fluorouracil/oxaliplatin. Applying this model predictor to protein profiles in tumor and matched normal tissue samples from 26 patients with CRCs, we found that differences in protein quantities were sufficient to model the increased tumor sensitivity to chemotherapy compared with normal tissue. In addition, these differences were sufficient to differentiate clinical responders from nonresponders with high confidence. Applications of our model, termed DR_MOMP, were used to assess the impact of apoptosis-sensitizing drugs in lowering the necessary dose of state-of-the-art chemotherapy in individual patients. Together, our findings offer a ready clinical tool with the potential to tailor chemotherapy to individual patients.

Glucose metabolism determines resistance of cancer cells to bioenergetic crisis after cytochrome-c release

Many anticancer drugs activate caspases via the mitochondrial apoptosis pathway. Activation of this pathway triggers a concomitant bioenergetic crisis caused by the release of cytochrome‐c (cyt‐c). Cancer cells are able to evade these processes by altering metabolic and caspase activation pathways. In this study, we provide the first integrated system study of mitochondrial bioenergetics and apoptosis signalling and examine the role of mitochondrial cyt‐c release in these events. In accordance with single‐cell experiments, our model showed that loss of cyt‐c decreased mitochondrial respiration by 95% and depolarised mitochondrial membrane potential ΔΨm from −142 to −88 mV, with active caspase‐3 potentiating this decrease. ATP synthase was reversed under such conditions, consuming ATP and stabilising ΔΨm. However, the direction and level of ATP synthase activity showed significant heterogeneity in individual cancer cells, which the model explained by variations in (i) accessible cyt‐c after release and (ii) the cells glycolytic capacity. Our results provide a quantitative and mechanistic explanation for the protective role of enhanced glucose utilisation for cancer cells to avert the otherwise lethal bioenergetic crisis associated with apoptosis initiation.

Predicting pK(a) values of substituted phenols from atomic charges: comparison of different quantum mechanical methods and charge distribution schemes.

The acid dissociation (ionization) constant pK(a) is one of the fundamental properties of organic molecules. We have evaluated different computational strategies and models to predict the pK(a) values of substituted phenols using partial atomic charges. Partial atomic charges for 124 phenol molecules were calculated using 83 approaches containing seven theory levels (MP2, HF, B3LYP, BLYP, BP86, AM1, and PM3), three basis sets (6-31G*, 6-311G, STO-3G), and five population analyses (MPA, NPA, Hirshfeld, MK, and Löwdin). The correlations between pK(a) and various atomic charge descriptors were examined, and the best descriptors were selected for preparing the quantitative structure-property relationship (QSPR) models. One QSPR model was created for each of the 83 approaches to charge calculation, and then the accuracy of all these models was analyzed and compared. The pK(a)s predicted by most of the models correlate strongly with experimental pK(a) values. For example, more than 25% of the models have correlation coefficients (R²) greater than 0.95 and root-mean-square errors smaller than 0.49. All seven examined theory levels are applicable for pK(a) prediction from charges. The best results were obtained for the MP2 and HF level of theory. The most suitable basis set was found to be 6-31G*. The 6-311G basis set provided slightly weaker correlations, and unexpectedly also, the STO-3G basis set is applicable for the QSPR modeling of pK(a). The Mulliken, natural, and Löwdin population analyses provide accurate models for all tested theory levels and basis sets. The results provided by the Hirshfeld population analysis were also acceptable, but the QSPR models based on MK charges show only weak correlations.

Clinical application of a systems model of apoptosis execution for the prediction of colorectal cancer therapy responses and personalisation of therapy

Objective Key to the clinical management of colorectal cancer is identifying tools which aid in assessing patient prognosis and determining more effective and personalised treatment strategies. We evaluated whether an experimental systems biology strategy which analyses the susceptibility of cancer cells to undergo caspase activation can be exploited to predict patient responses to 5-fluorouracil-based chemotherapy and to case-specifically identify potential alternative targeted treatments to reactivate apoptosis. Design We quantified five essential apoptosis-regulating proteins (Pro-Caspases 3 and 9, APAF-1, SMAC and XIAP) in samples of Stage II (n=13) and III (n=17) tumour and normal colonic (n=8) tissue using absolute quantitative immunoblotting and employed systems simulations of apoptosis signalling to predict the susceptibility of tumour cells to execute apoptosis. Additional systems analyses assessed the efficacy of novel apoptosis-inducing therapeutics such as XIAP antagonists, proteasome inhibitors and Pro-Caspase-3-activating compounds in restoring apoptosis execution in apoptosis-incompetent tumours. Results Comparisons of caspase activity profiles demonstrated that the likelihood of colorectal tumours to undergo apoptosis decreases with advancing disease stage. Systems-level analysis correctly predicted positive or negative outcome in 85% (p=0.004) of colorectal cancer patients receiving 5-fluorouracil based chemotherapy and significantly outperformed common uni- and multi-variate statistical approaches. Modelling of individual patient responses to novel apoptosis-inducing therapeutics revealed markedly different inter-individual responses. Conclusions Our study represents the first proof-of-concept example demonstrating the significant clinical potential of systems biology-based approaches for predicting patient outcome and responsiveness to novel targeted treatment paradigms.

Two-step activation of FOXO3 by AMPK generates a coherent feed-forward loop determining excitotoxic cell fate

Cerebral ischemia and excitotoxic injury induce transient or permanent bioenergetic failure, and may result in neuronal apoptosis or necrosis. We have previously shown that ATP depletion and activation of AMP-activated protein kinase (AMPK) during excitotoxic injury induces neuronal apoptosis by transcription of the pro-apoptotic BH3-only protein, Bim. AMPK, however, also exerts pro-survival functions in neurons. The molecular switches that determine these differential outcomes are not well understood. Using an approach combining biochemistry, single-cell imaging and computational modeling, we here demonstrate that excitotoxic injury activated the bim promoter in a FOXO3-dependent manner. The activation of AMPK reduced AKT activation, and led to dephosphorylation and nuclear translocation of FOXO3. Subsequent mutation studies indicated that bim gene activation during excitotoxic injury required direct FOXO3 phosphorylation by AMPK in the nucleus as a second activation step. Inhibition of this phosphorylation prevented Bim expression and protected neurons against excitotoxic and oxygen/glucose deprivation-induced injury. Systems analysis and computational modeling revealed that these two activation steps defined a coherent feed-forward loop; a network motif capable of filtering any effects of short-term AMPK activation on bim gene induction. This may prevent unwanted AMPK-mediated Bim expression and apoptosis during transient or physiological bioenergetic stress.

APOPTO-CELL---a simulation tool and interactive database for analyzing cellular susceptibility to apoptosis

UNLABELLED We have developed a web service that provides a comprehensive analysis of the susceptibility of cells to undergo apoptosis in response to an activation of the mitochondrial apoptotic pathway. Based on ordinary differential equations, (pre-determined) protein concentrations and release kinetics of mitochondrial pro-apoptotic factors, a network of 52 reactions and 19 reaction partners can be employed as a tool to display temporal protein profiles, to identify key regulatory proteins and to determine critical threshold concentrations required for the execution of apoptosis in HeLa cancer cells or other cell types. The web service also provides an interactive database function for the deposition of cell-type-specific quantitative data. In addition, the web service provides an output that can be compared directly to experimental results obtained from real-time single-cell experiments, making this a widely applicable systems biology tool for apoptosis and cancer researchers. AVAILABILITY http://systemsbiology.rcsi.ie/apopto-cell.html. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Activation of executioner caspases is a predictor of progression-free survival in glioblastoma patients: a systems medicine approach

Glioblastoma (GBM) is the most common and aggressive primary brain tumor in adults. GBM cells are highly resistant to apoptosis induced by antitumor drugs and radiotherapy resulting in cancer progression. We assessed whether a systems medicine approach, analysing the ability of tumor cells to execute apoptosis could be utilized to predict the response of GBM patients to treatment. Concentrations of the key proapoptotic proteins procaspase-3, procaspase-9, Smac and Apaf-1 and the antiapopotic protein XIAP were determined in a panel of GBM cell lines and GBM patient tumor resections. These values were used as input for APOPTO-CELL, a systems biological based mathematical model built to predict cellular susceptibility to undergo caspase activation. The modeling was capable of accurately distinguishing between GBM cells that die or survive in response to treatment with temozolomide in 10 of the 11 lines analysed. Importantly the results obtained using GBM patient samples show that APOPTO-CELL was capable of stratifying patients according to their progression-free survival times and predicted the ability of tumor cells to support caspase activation in 16 of the 21 GBM patients analysed. Calculating the susceptibility to apoptosis execution may be a potent tool in predicting GBM patient therapy responsiveness and may allow for the use of APOPTO-CELL in a clinical setting.