Identification and analysis of mutational hotspots in oncogenes and tumour suppressors

Mutational hotspots lead to robust but suboptimal adaptive outcomes in certain environments

The observed mutational spectrum of adaptive outcomes can be constrained by many factors. For example, mutational biases can narrow the observed spectrum by increasing the rate of mutation at isolated sites in the genome. In contrast, complex environments can shift the observed spectrum by defining fitness consequences of mutational routes. We investigate the impact of different nutrient environments on the evolution of motility in Pseudomonas fluorescens Pf0-2x (an engineered non-motile derivative of Pf0-1) in the presence and absence of a strong mutational hotspot. Previous work has shown that this mutational hotspot can be built and broken via six silent mutations, which provide rapid access to a mutation that rescues swimming motility and confers the strongest swimming phenotype in specific environments. Here, we evolved a hotspot and non-hotspot variant strain of Pf0-2x for motility under nutrient-rich (LB) and nutrient-limiting (M9) environmental conditions. We observed the hotspot strain consistently evolved faster across all environmental conditions and its mutational spectrum was robust to environmental differences. However, the non-hotspot strain had a distinct mutational spectrum that changed depending on the nutrient environment. Interestingly, while alternative adaptive mutations in nutrient-rich environments were equal to, or less effective than, the hotspot mutation, the majority of these mutations in nutrient-limited conditions produced superior swimmers. Our competition experiments mirrored these findings, underscoring the role of environment in defining both the mutational spectrum and the associated phenotype strength. This indicates that while mutational hotspots working in concert with natural selection can speed up access to robust adaptive mutations (which can provide a competitive advantage in evolving populations), they can limit exploration of the mutational landscape, restricting access to potentially stronger phenotypes in specific environments.

Detection and Quantification of ctDNA for Longitudinal Monitoring of Treatment in Non-Small Cell Lung Cancer Patients Using a Universal Mutant Detection Assay by Denaturing Capillary Electrophoresis

Background: Observation of anticancer therapy effect by monitoring of minimal residual disease (MRD) is becoming an important tool in management of non-small cell lung cancer (NSCLC). The approach is based on periodic detection and quantification of tumor-specific somatic DNA mutation in circulating tumor DNA (ctDNA) extracted from patient plasma. For such repetitive testing, complex liquid-biopsy techniques relying on ultra-deep NGS sequencing are impractical. There are other, cost-effective, methods for ctDNA analysis, typically based on quantitative PCR or digital PCR, which are applicable for detecting specific individual mutations in hotspots. While such methods are routinely used in NSCLC therapy prediction, however, extension to cover broader spectrum of mutations (e.g., in tumor suppressor genes) is required for universal longitudinal MRD monitoring.

Methods: For a set of tissue samples from 81 NSCLC patients we have applied a denaturing capillary electrophoresis (DCE) for initial detection of somatic mutations within 8 predesigned PCR amplicons covering oncogenes and tumor suppressor genes. Mutation-negative samples were then subjected to a large panel NGS sequencing. For each patient mutation found in tissue was then traced over time in ctDNA by DCE.

Results: In total we have detected a somatic mutation in tissue of 63 patients. For those we have then prospectively analyzed ctDNA from collected plasma samples over a period of up to 2 years. The dynamics of ctDNA during the initial chemotherapy therapy cycles as well as in the long-term follow-up matched the clinically observed response.

Conclusion: Detection and quantification of tumor-specific mutations in ctDNA represents a viable complement to MRD monitoring during therapy of NSCLC patients. The presented approach relying on initial tissue mutation detection by DCE combined with NGS and a subsequent ctDNA mutation testing by DCE only represents a cost-effective approach for its routine implementation.

Exploiting protein family and protein network data to identify novel drug targets for bladder cancer

Bladder cancer remains one of the most common forms of cancer and yet there are limited small molecule targeted therapies. Here, we present a computational platform to identify new potential targets for bladder cancer therapy. Our method initially exploited a set of known driver genes for bladder cancer combined with predicted bladder cancer genes from mutationally enriched protein domain families. We enriched this initial set of genes using protein network data to identify a comprehensive set of 323 putative bladder cancer targets. Pathway and cancer hallmarks analyses highlighted putative mechanisms in agreement with those previously reported for this cancer and revealed protein network modules highly enriched in potential drivers likely to be good targets for targeted therapies. 21 of our potential drug targets are targeted by FDA approved drugs for other diseases — some of them are known drivers or are already being targeted for bladder cancer (FGFR3, ERBB3, HDAC3, EGFR). A further 4 potential drug targets were identified by inheriting drug mappings across our in-house CATH domain functional families (FunFams). Our FunFam data also allowed us to identify drug targets in families that are less prone to side effects i.e., where structurally similar protein domain relatives are less dispersed across the human protein network. We provide information on our novel potential cancer driver genes, together with information on pathways, network modules and hallmarks associated with the predicted and known bladder cancer drivers and we highlight those drivers we predict to be likely drug targets.

Computational methods for detecting cancer hotspots

Cancer mutations that are recurrently observed among patients are known as hotspots. Hotspots are highly relevant because they are, presumably, likely functional. Known hotspots in BRAF, PIK3CA, TP53, KRAS, IDH1 support this idea. However, hundreds of hotspots have never been validated experimentally. The detection of hotspots nevertheless is challenging because background mutations obscure their statistical and computational identification. Although several algorithms have been applied to identify hotspots, they have not been reviewed before. Thus, in this mini-review, we summarize more than 40 computational methods applied to detect cancer hotspots in coding and non-coding DNA. We first organize the methods in cluster-based, 3D, position-specific, and miscellaneous to provide a general overview. Then, we describe their embed procedures, implementations, variations, and differences. Finally, we discuss some advantages, provide some ideas for future developments, and mention opportunities such as application to viral integrations, translocations, and epigenetics.

A missense variant remote from the active site impairs stability of human phosphoglucomutase 1

Missense variants of human phosphoglucomutase 1 (PGM1) cause the inherited metabolic disease known as PGM1 deficiency. This condition is categorised as both a glycogen storage disease and a congenital disorder of glycosylation. Approximately 20 missense variants of PGM1 are linked to PGM1 deficiency, and biochemical studies have suggested that they fall into two general categories: those affecting the active site and catalytic efficiency, and those that appear to impair protein folding and/or stability. In this study, we characterise a novel variant of Arg422, a residue distal from the active site of PGM1 and the site of a previously identified disease-related variant (Arg422Trp). In prior studies, the R422W variant was found to produce insoluble protein in a recombinant expression system, precluding further in vitro characterisation. Here we investigate an alternative variant of this residue, Arg422Gln, which is amenable to experimental characterisation presumably due to its more conservative physicochemical substitution. Biochemical, crystallographic, and computational studies of R422Q establish that this variant causes only minor changes in catalytic efficiency and 3D structure, but is nonetheless dramatically reduced in stability. Unexpectedly, binding of a substrate analog is found to further destabilise the protein, in contrast to its stabilising effect on wild-type PGM1 and several other missense variants. This work establishes Arg422 as a lynchpin residue for the stability of PGM1 and supports the impairment of protein stability as a pathomechanism for variants that cause PGM1 deficiency. SYNOPSIS: Biochemical and structural studies of a missense variant far from the active site of human PGM1 identify a residue with a key role in enzyme stability.

Modeling and analysis of site-specific mutations in cancer identifies known plus putative novel hotspots and bias due to contextual sequences

In cancer, recurrently mutated sites in DNA and proteins, called hotspots, are thought to be raised by positive selection and therefore important due to its potential functional impact. Although recent evidence for APOBEC enzymatic activity have shown that specific types of sequences are likely to be false, the identification of putative hotspots is important to confirm either its functional role or its mechanistic bias. In this work, an algorithm and a statistical model is presented to detect hotspots. The model consists of a beta-binomial component plus fixed effects that efficiently fits the distribution of mutated sites. The algorithm employs an optimal stepwise approach to find the model parameters. Simulations show that the proposed algorithmic model is highly accurate for common hotspots. The approach has been applied to TCGA mutational data from 33 cancer types. The results show that well-known cancer hotspots are easily detected. Besides, novel hotspots are also detected. An analysis of the sequence context of detected hotspots show a preference for TCG sites that may be related to APOBEC or other unknown mechanistic biases. The detected hotspots are available online in http://bioinformatica.mty.itesm.mx/HotSpotsAnnotations.

Integrating context of tumor biology and vaccine design to shape multidimensional immunotherapies

Advances in cancer therapy have offered great promise but only modest clinical benefits as monotherapies to date. Patients usually respond well to therapies targeted at specific mutations, but only for a short time. Conversely, immunotherapies help fewer patients, but increase survival. Combination therapies, which could offer the best of both worlds, are currently limited by substantial toxicity. While recent advances in genomics and proteomics have yielded an unprecedented depth of enabling datasets, it has also shifted the focus toward  in silico predictions. Designing the next wave of multidimensional immunotherapies will require leveraging this knowledge while providing a renewed emphasis on tumor biology and vaccine design. This includes careful selection of tumor clinical stage in the context of pre-existing tumor microenvironments, target antigen and technology platform selections to maximize their effect, and treatment staging. Here, we review strategies on how to approach an increasingly complex landscape of immunotherapeutic agents for use in combination therapies.

Passenger Hotspot Mutations in Cancer

Current statistical models for assessing hotspot significance do not properly account for variation in site-specific mutability, thereby yielding many false-positives. We thus (i) detail a Log-normal-Poisson (LNP) background model that accounts for this variability in a manner consistent with models of mutagenesis; (ii) use it to show that passenger hotspots arise from all common mutational processes; and (iii) apply it to a ∼10,000-patient cohort to nominate driver hotspots with far fewer false-positives compared with conventional methods. Overall, we show that many cancer hotspot mutations recurring at the same genomic site across multiple tumors are actually passenger events, recurring at inherently mutable genomic sites under no positive selection.

Towards quantitative and multiplexed in vivo functional cancer genomics

Large-scale sequencing of human tumours has uncovered a vast array of genomic alterations. Genetically engineered mouse models recapitulate many features of human cancer and have been instrumental in assigning biological meaning to specific cancer-associated alterations. However, their time, cost and labour-intensive nature limits their broad utility; thus, the functional importance of the majority of genomic aberrations in cancer remains unknown. Recent advances have accelerated the functional interrogation of cancer-associated alterations within in vivo models. Specifically, the past few years have seen the emergence of CRISPR-Cas9-based strategies to rapidly generate increasingly complex somatic alterations and the development of multiplexed and quantitative approaches to ascertain gene function in vivo.

Mitochondria as playmakers of apoptosis, autophagy and senescence

Mitochondria are the key energy-producing organelles and cellular source of reactive species. They are responsible for managing cell life and death by a balanced homeostasis passing through a network of structures, regulated principally via fission and fusion. Herein we discuss about the most advanced findings considering mitochondria as dynamic biophysical systems playing compelling roles in the regulation of energy metabolism in both physiologic and pathologic processes controlling cell death and survival. Precisely, we focus on the mitochondrial commitment to the onset, maintenance and counteraction of apoptosis, autophagy and senescence in the bioenergetic reprogramming of cancer cells. In this context, looking for a pharmacological manipulation of cell death processes as a successful route for future targeted therapies, there is major biotechnological challenge in underlining the location, function and molecular mechanism of mitochondrial proteins. Based on the critical role of mitochondrial functions for cellular health, a better knowledge of the main molecular players in mitochondria disfunction could be decisive for the therapeutical control of degenerative diseases, including cancer.

Predicting synthetic lethal interactions using conserved patterns in protein interaction networks

In response to a need for improved treatments, a number of promising novel targeted cancer therapies are being developed that exploit human synthetic lethal interactions. This is facilitating personalised medicine strategies in cancers where specific tumour suppressors have become inactivated. Mainly due to the constraints of the experimental procedures, relatively few human synthetic lethal interactions have been identified. Here we describe SLant (Synthetic Lethal analysis via Network topology), a computational systems approach to predicting human synthetic lethal interactions that works by identifying and exploiting conserved patterns in protein interaction network topology both within and across species. SLant out-performs previous attempts to classify human SSL interactions and experimental validation of the models predictions suggests it may provide useful guidance for future SSL screenings and ultimately aid targeted cancer therapy development.

A CATH domain functional family based approach to identify putative cancer driver genes and driver mutations

Tumour sequencing identifies highly recurrent point mutations in cancer driver genes, but rare functional mutations are hard to distinguish from large numbers of passengers. We developed a novel computational platform applying a multi-modal approach to filter out passengers and more robustly identify putative driver genes. The primary filter identifies enrichment of cancer mutations in CATH functional families (CATH-FunFams) – structurally and functionally coherent sets of evolutionary related domains. Using structural representatives from CATH-FunFams, we subsequently seek enrichment of mutations in 3D and show that these mutation clusters have a very significant tendency to lie close to known functional sites or conserved sites predicted using CATH-FunFams. Our third filter identifies enrichment of putative driver genes in functionally coherent protein network modules confirmed by literature analysis to be cancer associated. Our approach is complementary to other domain enrichment approaches exploiting Pfam families, but benefits from more functionally coherent groupings of domains. Using a set of mutations from 22 cancers we detect 151 putative cancer drivers, of which 79 are not listed in cancer resources and include recently validated cancer associated genes EPHA7, DCC netrin-1 receptor and zinc-finger protein ZNF479.

Biophysical simulations and structure-based modeling of residue interaction networks in the tumor suppressor proteins reveal functional role of cancer mutation hotspots in molecular communication

In the current study, we have combined molecular simulations and energetic analysis with dynamics-based network modeling and perturbation response scanning to determine molecular signatures of mutational hotspot residues in the p53, PTEN, and SMAD4 tumor suppressor proteins. By examining structure, energetics and dynamics of these proteins, we have shown that inactivating mutations preferentially target a group of structurally stable residues that play a fundamental role in global propagation of dynamic fluctuations and mediating allosteric interaction networks. Through integration of long-range perturbation dynamics and network-based approaches, we have quantified allosteric potential of residues in the studied proteins. The results have revealed that mutational hotspot sites often correspond to high centrality mediating centers of the residue interaction networks that are responsible for coordination of global dynamic changes and allosteric signaling. Our findings have also suggested that structurally stable mutational hotpots can act as major effectors of allosteric interactions and mutations in these positions are typically associated with severe phenotype. Modeling of shortest inter-residue pathways has shown that mutational hotspot sites can also serve as key mediating bridges of allosteric communication in the p53 and PTEN protein structures. Multiple regression models have indicated that functional significance of mutational hotspots can be strongly associated with the network signatures serving as robust predictors of critical regulatory positions responsible for loss-of-function phenotype. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of mutational hotspots, providing a plausible rationale for explaining and localizing disease-causing mutations in tumor suppressor genes.

Mapping genetic variations to three-dimensional protein structures to enhance variant interpretation: a proposed framework

The translation of personal genomics to precision medicine depends on the accurate interpretation of the multitude of genetic variants observed for each individual. However, even when genetic variants are predicted to modify a protein, their functional implications may be unclear. Many diseases are caused by genetic variants affecting important protein features, such as enzyme active sites or interaction interfaces. The scientific community has catalogued millions of genetic variants in genomic databases and thousands of protein structures in the Protein Data Bank. Mapping mutations onto three-dimensional (3D) structures enables atomic-level analyses of protein positions that may be important for the stability or formation of interactions; these may explain the effect of mutations and in some cases even open a path for targeted drug development. To accelerate progress in the integration of these data types, we held a two-day Gene Variation to 3D (GVto3D) workshop to report on the latest advances and to discuss unmet needs. The overarching goal of the workshop was to address the question: what can be done together as a community to advance the integration of genetic variants and 3D protein structures that could not be done by a single investigator or laboratory? Here we describe the workshop outcomes, review the state of the field, and propose the development of a framework with which to promote progress in this arena. The framework will include a set of standard formats, common ontologies, a common application programming interface to enable interoperation of the resources, and a Tool Registry to make it easy to find and apply the tools to specific analysis problems. Interoperability will enable integration of diverse data sources and tools and collaborative development of variant effect prediction methods.

Bioinformatics in translational drug discovery

Bioinformatics approaches are becoming ever more essential in translational drug discovery both in academia and within the pharmaceutical industry. Computational exploitation of the increasing volumes of data generated during all phases of drug discovery is enabling key challenges of the process to be addressed. Here, we highlight some of the areas in which bioinformatics resources and methods are being developed to support the drug discovery pipeline. These include the creation of large data warehouses, bioinformatics algorithms to analyse 'big data' that identify novel drug targets and/or biomarkers, programs to assess the tractability of targets, and prediction of repositioning opportunities that use licensed drugs to treat additional indications.