Web-based NGS data analysis using miRMaster: a large-scale meta-analysis of human miRNAs.

The analysis of small RNA NGS data together with the discovery of new small RNAs is among the foremost challenges in life science. For the analysis of raw high-throughput sequencing data we implemented the fast, accurate and comprehensive web-based tool miRMaster. Our toolbox provides a wide range of modules for quantification of miRNAs and other non-coding RNAs, discovering new miRNAs, isomiRs, mutations, exogenous RNAs and motifs. Use-cases comprising hundreds of samples are processed in less than 5 h with an accuracy of 99.4%. An integrative analysis of small RNAs from 1836 data sets (20 billion reads) indicated that context-specific miRNAs (e.g. miRNAs present only in one or few different tissues / cell types) still remain to be discovered while broadly expressed miRNAs appear to be largely known. In total, our analysis of known and novel miRNAs indicated nearly 22 000 candidates of precursors with one or two mature forms. Based on these, we designed a custom microarray comprising 11 872 potential mature miRNAs to assess the quality of our prediction. MiRMaster is a convenient-to-use tool for the comprehensive and fast analysis of miRNA NGS data. In addition, our predicted miRNA candidates provided as custom array will allow researchers to perform in depth validation of candidates interesting to them.

A Stepwise Integrated Approach to Personalized Risk Predictions in Stage III Colorectal Cancer.

Purpose: Apoptosis is essential for chemotherapy responses. In this discovery and validation study, we evaluated the suitability of a mathematical model of apoptosis execution (APOPTO-CELL) as a stand-alone signature and as a constituent of further refined prognostic stratification tools.Experimental Design: Apoptosis competency of primary tumor samples from patients with stage III colorectal cancer (n = 120) was calculated by APOPTO-CELL from measured protein concentrations of Procaspase-3, Procaspase-9, SMAC, and XIAP. An enriched APOPTO-CELL signature (APOPTO-CELL-PC3) was synthesized to capture apoptosome-independent effects of Caspase-3. Furthermore, a machine learning Random Forest approach was applied to APOPTO-CELL-PC3 and available molecular and clinicopathologic data to identify a further enhanced signature. Association of the signature with prognosis was evaluated in an independent colon adenocarcinoma cohort (TCGA COAD, n = 136).Results: We identified 3 prognostic biomarkers (P = 0.04, P = 0.006, and P = 0.0004 for APOPTO-CELL, APOPTO-CELL-PC3, and Random Forest signatures, respectively) with increasing stratification accuracy for patients with stage III colorectal cancer.The APOPTO-CELL-PC3 signature ranked highest among all features. The prognostic value of the signatures was independently validated in stage III TCGA COAD patients (P = 0.01, P = 0.04, and P = 0.02 for APOPTO-CELL, APOPTO-CELL-PC3, and Random Forest signatures, respectively). The signatures provided further stratification for patients with CMS1-3 molecular subtype.Conclusions: The integration of a systems-biology-based biomarker for apoptosis competency with machine learning approaches is an appealing and innovative strategy toward refined patient stratification. The prognostic value of apoptosis competency is independent of other available clinicopathologic and molecular factors, with tangible potential of being introduced in the clinical management of patients with stage III colorectal cancer. Clin Cancer Res; 23(5); 1200-12. ©2016 AACR.

Predicting the cell death responsiveness and sensitization of glioma cells to TRAIL and temozolomide.

Genotoxic chemotherapy with temozolomide (TMZ) is a mainstay of treatment for glioblastoma (GBM); however, at best, TMZ provides only modest survival benefit to a subset of patients. Recent insight into the heterogeneous nature of GBM suggests a more personalized approach to treatment may be necessary to overcome cancer drug resistance and improve patient care. These include novel therapies that can be used both alone and with TMZ to selectively reactivate apoptosis within malignant cells. For this approach to work, reliable molecular signatures that can accurately predict treatment responsiveness need to be identified first. Here, we describe the first proof-of-principle study that merges quantitative protein-based analysis of apoptosis signaling networks with data- and knowledge-driven mathematical systems modeling to predict treatment responsiveness of GBM cell lines to various apoptosis-inducing stimuli. These include monotherapies with TMZ and TRAIL, which activate the intrinsic and extrinsic apoptosis pathways, respectively, as well as combination therapies of TMZ+TRAIL. We also successfully employed this approach to predict whether individual GBM cell lines could be sensitized to TMZ or TRAIL via the selective targeting of Bcl-2/Bcl-xL proteins with ABT-737. Our findings suggest that systems biology-based approaches could assist in personalizing treatment decisions in GBM to optimize cell death induction.

An Analysis of the Truncated Bid- and ROS-dependent Spatial Propagation of Mitochondrial Permeabilization Waves during Apoptosis.

Apoptosis is a form of programmed cell death that is essential for the efficient elimination of surplus, damaged, and transformed cells during metazoan embryonic development and adult tissue homeostasis. Situated at the interface of apoptosis initiation and execution, mitochondrial outer membrane permeabilization (MOMP) represents one of the most fundamental processes during apoptosis signal transduction. It was shown that MOMP can spatiotemporally propagate through cells, in particular in response to extrinsic apoptosis induction. Based on apparently contradictory experimental evidence, two distinct molecular mechanisms have been proposed to underlie the propagation of MOMP signals, namely a reaction-diffusion mechanism governed by anisotropies in the production of the MOMP-inducer truncated Bid (tBid), or a process that drives the spatial propagation of MOMP by sequential bursts of reactive oxygen species. We therefore generated mathematical models for both scenarios and performed in silico simulations of spatiotemporal MOMP signaling to identify which one of the two mechanisms is capable of qualitatively and quantitatively reproducing the existing data. We found that the explanatory power of each model was limited in that only a subset of experimental findings could be replicated. However, the integration of both models into a combined mathematical description of spatiotemporal tBid and reactive oxygen species signaling accurately reproduced all available experimental data and furthermore, provided robustness to spatial MOMP propagation when mitochondria are spatially separated. Our study therefore provides a theoretical framework that is sufficient to describe and mechanistically explain the spatiotemporal propagation of one of the most fundamental processes during apoptotic cell death.

A systems biology analysis of apoptosome formation and apoptosis execution supports allosteric procaspase-9 activation.

The protease caspase-9 is activated on the apoptosome, a multiprotein signal transduction platform that assembles in response to mitochondria-dependent apoptosis initiation. Despite extensive molecular research, the assembly of the holo-apoptosome and the process of caspase-9 activation remain incompletely understood. Here, we therefore integrated quantitative data on the molecular interactions and proteolytic processes during apoptosome formation and apoptosis execution and conducted mathematical simulations to investigate the resulting biochemical signaling, quantitatively and kinetically. Interestingly, when implementing the homodimerization of procaspase-9 as a prerequisite for activation, the calculated kinetics of apoptosis execution and the efficacy of caspase-3 activation failed to replicate experimental data. In contrast, assuming a scenario in which procaspase-9 is activated allosterically upon binding to the apoptosome backbone, the mathematical simulations quantitatively and kinetically reproduced all experimental data. These data included a XIAP threshold concentration at which apoptosis execution is suppressed in HeLa cervical cancer cells, half-times of procaspase-9 processing, as well as the molecular timer function of the apoptosome. Our study therefore provides novel mechanistic insight into apoptosome-dependent apoptosis execution and suggests that caspase-9 is activated allosterically by binding to the apoptosome backbone. Our findings challenge the currently prevailing dogma that all initiator procaspases require homodimerization for activation.

The central role of initiator caspase-9 in apoptosis signal transduction and the regulation of its activation and activity on the apoptosome.

Key structural and catalytic features are conserved across the entire family of cysteine-dependent aspartate-specific proteases (caspases). Of the caspases involved in apoptosis signal transduction, the initiator caspases-2, -8 and -9 are activated at multi-protein activation platforms, and activation is thought to involve homo-dimerisation of the monomeric zymogens. Caspase-9, the essential initiator caspase required for apoptosis signalling through the mitochondrial pathway, is activated on the apoptosome complex, and failure to activate caspase-9 has profound pathophysiological consequences. Here, we review the pertinent literature on which the currently prevalent understanding of caspase-9 activation is based, extend this view by insight obtained from recent structural and kinetic studies on caspase-9 signalling, and describe an emerging model for the regulation of caspase-9 activation and activity that arise from the complexity of multi-protein interactions at the apoptosome. This integrated view allows us to postulate and to discuss functional consequences for caspase-9 activation and apoptosis execution that may take centre stage in future experimental cell research on apoptosis signalling.

The caspase-8 dimerization/dissociation balance is a highly potent regulator of caspase-8, -3, -6 signaling.

Apoptosis is driven by positive feedback activation between aspartate-specific cysteinyl proteases (caspases). These feedback loops ensure the swift and efficient elimination of cells upon initiation of apoptosis execution. At the same time, the signaling network must be insensitive to erroneous, mild caspase activation to avoid unwanted, excessive cell death. Sublethal caspase activation in fact was shown to be a requirement for the differentiation of multiple cell types but might also occur accidentally during short, transient cellular stress conditions. Here we carried out an in silico comparison of the molecular mechanisms that so far have been identified to impair the amplification of caspase activities via the caspase-8, -3, -6 loop. In a systems model resembling HeLa cervical cancer cells, the dimerization/dissociation balance of caspase-8 potently suppressed the amplification of caspase responses, surprisingly outperforming or matching known caspase-8 and -3 inhibitors such as bifunctional apoptosis repressor or x-linked inhibitor of apoptosis protein. These findings were further substantiated in global sensitivity analyses based on combinations of protein concentrations from the sub- to superphysiological range to screen the full spectrum of biological variability that can be expected within cell populations and between distinct cell types. Additional modeling showed that the combined effects of x-linked inhibitor of apoptosis protein and caspase-8 dimerization/dissociation processes can also provide resistance to larger inputs of active caspases. Our study therefore highlights a central and so far underappreciated role of caspase-8 dimerization/dissociation in avoiding unwanted cell death by lethal amplification of caspase responses via the caspase-8, -3, -6 loop.