High-resolution mapping of cancer cell networks using co-functional interactions.

Powerful new technologies for perturbing genetic elements have recently expanded the study of genetic interactions in model systems ranging from yeast to human cell lines. However, technical artifacts can confound signal across genetic screens and limit the immense potential of parallel screening approaches. To address this problem, we devised a novel PCA-based method for correcting genome-wide screening data, bolstering the sensitivity and specificity of detection for genetic interactions. Applying this strategy to a set of 436 whole genome CRISPR screens, we report more than 1.5 million pairs of correlated "co-functional" genes that provide finer-scale information about cell compartments, biological pathways, and protein complexes than traditional gene sets. Lastly, we employed a gene community detection approach to implicate core genes for cancer growth and compress signal from functionally related genes in the same community into a single score. This work establishes new algorithms for probing cancer cell networks and motivates the acquisition of further CRISPR screen data across diverse genotypes and cell types to further resolve complex cellular processes

A predictive computational model reveals that GIV/girdin serves as a tunable valve for EGFR-stimulated cyclic AMP signals.

Cellular levels of the versatile second messenger cyclic (c)AMP are regulated by the antagonistic actions of the canonical G protein → adenylyl cyclase pathway that is initiated by G-protein-coupled receptors (GPCRs) and attenuated by phosphodiesterases (PDEs). Dysregulated cAMP signaling drives many diseases; for example, its low levels facilitate numerous sinister properties of cancer cells. Recently, an alternative paradigm for cAMP signaling has emerged in which growth factor-receptor tyrosine kinases (RTKs; e.g., EGFR) access and modulate G proteins via a cytosolic guanine-nucleotide exchange modulator (GEM), GIV/girdin; dysregulation of this pathway is frequently encountered in cancers. In this study, we present a network-based compartmental model for the paradigm of GEM-facilitated cross-talk between RTKs and G proteins and how that impacts cellular cAMP. Our model predicts that cross-talk between GIV, Gαs, and Gαi proteins dampens ligand-stimulated cAMP dynamics. This prediction was experimentally verified by measuring cAMP levels in cells under different conditions. We further predict that the direct proportionality of cAMP concentration as a function of receptor number and the inverse proportionality of cAMP concentration as a function of PDE concentration are both altered by GIV levels. Taking these results together, our model reveals that GIV acts as a tunable control valve that regulates cAMP flux after growth factor stimulation. For a given stimulus, when GIV levels are high, cAMP levels are low, and vice versa. In doing so, GIV modulates cAMP via mechanisms distinct from the two most often targeted classes of cAMP modulators, GPCRs and PDEs

Modeling cancer metabolism on a genome scale

Cancer cells have fundamentally altered cellular metabolism that is associated with their tumorigenicity and malignancy. In addition to the widely studied Warburg effect, several new key metabolic alterations in cancer have been established over the last decade, leading to the recognition that altered tumor metabolism is one of the hallmarks of cancer. Deciphering the full scope and functional implications of the dysregulated metabolism in cancer requires both the advancement of a variety of omics measurements and the advancement of computational approaches for the analysis and contextualization of the accumulated data. Encouragingly, while the metabolic network is highly interconnected and complex, it is at the same time probably the best characterized cellular network. Following, this review discusses the challenges that genome‐scale modeling of cancer metabolism has been facing. We survey several recent studies demonstrating the first strides that have been done, testifying to the value of this approach in portraying a network‐level view of the cancer metabolism and in identifying novel drug targets and biomarkers. Finally, we outline a few new steps that may further advance this field

DRUG-NEM: Optimizing drug combinations using single-cell perturbation response to account for intratumoral heterogeneity.

An individual malignant tumor is composed of a heterogeneous collection of single cells with distinct molecular and phenotypic features, a phenomenon termed intratumoral heterogeneity. Intratumoral heterogeneity poses challenges for cancer treatment, motivating the need for combination therapies. Single-cell technologies are now available to guide effective drug combinations by accounting for intratumoral heterogeneity through the analysis of the signaling perturbations of an individual tumor sample screened by a drug panel. In particular, Mass Cytometry Time-of-Flight (CyTOF) is a high-throughput single-cell technology that enables the simultaneous measurements of multiple ([Formula: see text]40) intracellular and surface markers at the level of single cells for hundreds of thousands of cells in a sample. We developed a computational framework, entitled Drug Nested Effects Models (DRUG-NEM), to analyze CyTOF single-drug perturbation data for the purpose of individualizing drug combinations. DRUG-NEM optimizes drug combinations by choosing the minimum number of drugs that produce the maximal desired intracellular effects based on nested effects modeling. We demonstrate the performance of DRUG-NEM using single-cell drug perturbation data from tumor cell lines and primary leukemia samples

ATM in focus:a damage sensor and cancer target

The ability of a cell to conserve and maintain its native DNA sequence is fundamental for the survival and normal functioning of the whole organism and protection from cancer development. Here we review recently obtained results and current topics concerning the role of the ataxia-telangiectasia mutated (ATM) protein kinase as a damage sensor and its potential as therapeutic target for treating cancer. This monograph discusses DNA repair mechanisms activated after DNA double-strand breaks (DSBs), i.e. non-homologous end joining, homologous recombination and single strand annealing and the role of ATM in the above types of repair. In addition to DNA repair, ATM participates in a diverse set of physiological processes involving metabolic regulation, oxidative stress, transcriptional modulation, protein degradation and cell proliferation. Full understanding of the complexity of ATM functions and the design of therapeutics that modulate its activity to combat diseases such as cancer necessitates parallel theoretical and experimental efforts. This could be best addressed by employing a systems biology approach, involving mathematical modelling of cell signalling pathways

Reverse engineering of drug induced DNA damage response signalling pathway reveals dual outcomes of ATM kinase inhibition

The DNA Damage Response (DDR) pathway represents a signalling mechanism that is activated in eukaryotic cells following DNA damage and comprises of proteins involved in DNA damage detection, DNA repair, cell cycle arrest and apoptosis. This pathway consists of an intricate network of signalling interactions driving the cellular ability to recognise DNA damage and recruit specialised proteins to take decisions between DNA repair or apoptosis. ATM and ATR are central components of the DDR pathway. The activities of these kinases are vital in DNA damage induced phosphorylational induction of DDR substrates. Here, firstly we have experimentally determined DDR signalling network surrounding the ATM/ATR pathway induced following double stranded DNA damage by monitoring and quantifying time dependent inductions of their phosphorylated forms and their key substrates. We next involved an automated inference of unsupervised predictive models of time series data to generate in silico (molecular) interaction maps. We characterized the complex signalling network through system analysis and gradual utilisation of small time series measurements of key substrates through a novel network inference algorithm. Furthermore, we demonstrate an application of an assumption-free reverse engineering of the intricate signalling network of the activated ATM/ATR pathway. We next studied the consequences of such drug induced inductions as well as of time dependent ATM kinase inhibition on cell survival through further biological experiments. Intermediate and temporal modelling outcomes revealed the distinct signaling profile associated with ATM kinase activity and inhibition and explained the underlying signalling mechanism for dual ATM functionality in cytotoxic and cytoprotective pathways

A microRNA molecular modeling extension for prediction of colorectal cancer treatment

Background: Several studies show that the regulatory impact of microRNAs (miRNAs) is an essential contribution to the pathogenesis of colorectal cancer (CRC). The expression levels of diverse miRNAs are associated with specific clinical diagnoses and prognoses of CRC. However, this association reveals very little actionable information with regard to how or whether to treat a CRC patient. To address this problem, we use miRNA expression data along with other molecular information to predict individual response of CRC cell lines and CRC patients. Methods: A strategy has been developed to join four types of information: molecular, kinetic, genetic and treatment data for prediction of individual treatment response of CRC. Results: Information on miRNA regulation, including miRNA target regulation and transcriptional regulation of miRNA, in integrated into an in silico molecular model for colon cancer. This molecular model is applied to study responses of seven CRC cell lines from NCI-60 to ten agents targeting signaling pathways. Predictive results of models without and with implemented miRNA information are compared and advantages are shown for the extended model. Finally, the extended model was applied to the data of 22 CRC patients to predict response to treatments of sirolimus and LY294002. The in silico results can also replicate the oncogenic and tumor suppression roles of miRNA on the therapeutic response as reported in the literature. Conclusions: In summary, the results reveal that detailed molecular events can be combined with individual genetic data, including gene/miRNA expression data, to enhance in silico prediction of therapeutic response of individual CRC tumors. The study demonstrates that miRNA information can be applied as actionable information regarding individual therapeutic response

Understanding resistance to combination chemotherapy

available in PMC 2014 April 04The current clinical application of combination chemotherapy is guided by a historically successful set of practices that were developed by basic and clinical researchers 50–60 years ago. Thus, in order to understand how emerging approaches to drug development might aid the creation of new therapeutic combinations, it is critical to understand the defining principles underlying classic combination therapy and the original experimental rationales behind them. One such principle is that the use of combination therapies with independent mechanisms of action can minimize the evolution of drug resistance. Another is that in order to kill sufficient cancer cells to cure a patient, multiple drugs must be delivered at their maximum tolerated dose – a condition that allows for enhanced cancer cell killing with manageable toxicity. In light of these models, we aim to explore recent genomic evidence underlying the mechanisms of resistance to the combination regimens constructed on these principles. Interestingly, we find that emerging genomic evidence contradicts some of the rationales of early practitioners in developing commonly used drug regimens. However, we also find that the addition of recent targeted therapies has yet to change the current principles underlying the construction of anti-cancer combinatorial regimens, nor have they made substantial inroads into the treatment of most cancers. We suggest that emerging systems/network biology approaches have an immense opportunity to impact the rational development of successful drug regimens. Specifically, by examining drug combinations in multivariate ways, next generation combination therapies can be constructed with a clear understanding of how mechanisms of resistance to multi-drug regimens differ from single agent resistance.Massachusetts Institute of Technology (Eisen and Chang Career Development Associate Professor of Biology)National Cancer Institute (U.S.) (NCI Integrative Cancer Biology Program (ICBP), #U54-CA112967-06)National Institutes of Health (U.S.) (NIH RO1-CA128803-04

Systematic Analysis of Quantitative Logic Model Ensembles Predicts Drug Combination Effects on Cell Signaling Networks

A major challenge in developing anticancer therapies is determining the efficacies of drugs and their combinations in physiologically relevant microenvironments. We describe here our application of “constrained fuzzy logic” (CFL) ensemble modeling of the intracellular signaling network for predicting inhibitor treatments that reduce the phospho-levels of key transcription factors downstream of growth factors and inflammatory cytokines representative of hepatocellular carcinoma (HCC) microenvironments. We observed that the CFL models successfully predicted the effects of several kinase inhibitor combinations. Furthermore, the ensemble predictions revealed ambiguous predictions that could be traced to a specific structural feature of these models, which we resolved with dedicated experiments, finding that IL-1α activates downstream signals through TAK1 and not MEKK1 in HepG2 cells. We conclude that CFL-Q2LM (Querying Quantitative Logic Models) is a promising approach for predicting effective anticancer drug combinations in cancer-relevant microenvironments.United States. Army Research Office (W911NF-09-0001

INTEGRATIVE ANALYSIS OF OMICS DATA IN ADULT GLIOMA AND OTHER TCGA CANCERS TO GUIDE PRECISION MEDICINE

Transcriptomic profiling and gene expression signatures have been widely applied as effective approaches for enhancing the molecular classification, diagnosis, prognosis or prediction of therapeutic response towards personalized therapy for cancer patients. Thanks to modern genome-wide profiling technology, scientists are able to build engines leveraging massive genomic variations and integrating with clinical data to identify “at risk” individuals for the sake of prevention, diagnosis and therapeutic interventions. In my graduate work for my Ph.D. thesis, I have investigated genomic sequencing data mining to comprehensively characterise molecular classifications and aberrant genomic events associated with clinical prognosis and treatment response, through applying high-dimensional omics genomic data to promote the understanding of gene signatures and somatic molecular alterations contributing to cancer progression and clinical outcomes. Following this motivation, my dissertation has been focused on the following three topics in translational genomics. 1) Characterization of transcriptomic plasticity and its association with the tumor microenvironment in glioblastoma (GBM). I have integrated transcriptomic, genomic, protein and clinical data to increase the accuracy of GBM classification, and identify the association between the GBM mesenchymal subtype and reduced tumorpurity, accompanied with increased presence of tumor-associated microglia. Then I have tackled the sole source of microglial as intrinsic tumor bulk but not their corresponding neurosphere cells through both transcriptional and protein level analysis using a panel of sphere-forming glioma cultures and their parent GBM samples.FurthermoreI have demonstrated my hypothesis through longitudinal analysis of paired primary and recurrent GBM samples that the phenotypic alterations of GBM subtypes are not due to intrinsic proneural-to-mesenchymal transition in tumor cells, rather it is intertwined with increased level of microglia upon disease recurrence. Collectively I have elucidated the critical role of tumor microenvironment (Microglia and macrophages from central nervous system) contributing to the intra-tumor heterogeneity and accurate classification of GBM patients based on transcriptomic profiling, which will not only significantly impact on clinical perspective but also pave the way for preclinical cancer research. 2) Identification of prognostic gene signatures that stratify adult diffuse glioma patientsharboring1p/19q co-deletions. I have compared multiple statistical methods and derived a gene signature significantly associated with survival by applying a machine learning algorithm. Then I have identified inflammatory response and acetylation activity that associated with malignant progression of 1p/19q co-deleted glioma. In addition, I showed this signature translates to other types of adult diffuse glioma, suggesting its universality in the pathobiology of other subset gliomas. My efforts on integrative data analysis of this highly curated data set usingoptimizedstatistical models will reflect the pending update to WHO classification system oftumorsin the central nervous system (CNS). 3) Comprehensive characterization of somatic fusion transcripts in Pan-Cancers. I have identified a panel of novel fusion transcripts across all of TCGA cancer types through transcriptomic profiling. Then I have predicted fusion proteins with kinase activity and hub function of pathway network based on the annotation of genetically mobile domains and functional domain architectures. I have evaluated a panel of in -frame gene fusions as potential driver mutations based on network fusion centrality hypothesis. I have also characterised the emerging complexity of genetic architecture in fusion transcripts through integrating genomic structure and somatic variants and delineating the distinct genomic patterns of fusion events across different cancer types. Overall my exploration of the pathogenetic impact and clinical relevance of candidate gene fusions have provided fundamental insights into the management of a subset of cancer patients by predicting the oncogenic signalling and specific drug targets encoded by these fusion genes. Taken together, the translational genomic research I have conducted during my Ph.D. study will shed new light on precision medicine and contribute to the cancer research community. The novel classification concept, gene signature and fusion transcripts I have identified will address several hotly debated issues in translational genomics, such as complex interactions between tumor bulks and their adjacent microenvironments, prognostic markers for clinical diagnostics and personalized therapy, distinct patterns of genomic structure alterations and oncogenic events in different cancer types, therefore facilitating our understanding of genomic alterations and moving us towards the development of precision medicine