Oncogenic potential is related to activating effect of cancer single and double somatic mutations in receptor tyrosine kinases

Genetic features of patients with MPS type IIIB: Description of five pathogenic gene variations

BACKGROUND

There are four distinct forms of Sanfilippo syndrome (MPS type III), each of which is an autosomal lysosomal storage disorder. These forms are caused by abnormalities in one of four lysosomal enzymes. This study aimed to identify possible genetic variants that contribute to Sanfilippo IIIB in 14 independent families in Southwest Iran.

METHODS

Patients were included if their clinical features and enzyme assay results were suggestive. The patients were subsequently subjected to Sanger Sequencing to screen for Sanfilippo-related genes. Additional investigations have been conducted using various computational analyses to determine the probable functional effects of diagnosed variants.

RESULTS

Five distinct variations were identified in the NAGLU gene. This included two novel variants in two distinct families and three previously reported variants in 12 distinct families. All of these variations were recognized as pathogenic using the MutationTaster web server. In silico analysis showed that all detected variants affected protein structural stability; four destabilized protein structures, and the fifth variation had the opposite effect.

CONCLUSION

In this study, two novel variations in the NAGLU gene were identified. The results of this study positively contribute to the mutation diversity of the NAGLU gene. To identify new disease biomarkers and therapeutic targets, precision medicine must precisely characterize and account for genetic variations. New harmful gene variants are valuable for updating gene databases concerning Sanfilippo disease variations and NGS gene panels. This may also improve genetic counselling for rapid risk examinations and disease surveillance.

Applying Bioinformatic Platforms, In Vitro, and In Vivo Functional Assays in the Characterization of Genetic Variants in the GH/IGF Pathway Affecting Growth and Development

Heritability accounts for over 80% of adult human height, indicating that genetic variability is the main determinant of stature. The rapid technological development of Next-Generation Sequencing (NGS), particularly Whole Exome Sequencing (WES), has resulted in the characterization of several genetic conditions affecting growth and development. The greatest challenge of NGS remains the high number of candidate variants identified. In silico bioinformatic tools represent the first approach for classifying these variants. However, solving the complicated problem of variant interpretation requires the use of experimental approaches such as in vitro and, when needed, in vivo functional assays. In this review, we will discuss a rational approach to apply to the gene variants identified in children with growth and developmental defects including: (i) bioinformatic tools; (ii) in silico modeling tools; (iii) in vitro functional assays; and (iv) the development of in vivo models. While bioinformatic tools are useful for a preliminary selection of potentially pathogenic variants, in vitro—and sometimes also in vivo—functional assays are further required to unequivocally determine the pathogenicity of a novel genetic variant. This long, time-consuming, and expensive process is the only scientifically proven method to determine causality between a genetic variant and a human genetic disease.

PremPS: Predicting the impact of missense mutations on protein stability

Computational methods that predict protein stability changes induced by missense mutations have made a lot of progress over the past decades. Most of the available methods however have very limited accuracy in predicting stabilizing mutations because existing experimental sets are dominated by mutations reducing protein stability. Moreover, few approaches could consistently perform well across different test cases. To address these issues, we developed a new computational method PremPS to more accurately evaluate the effects of missense mutations on protein stability. The PremPS method is composed of only ten evolutionary- and structure-based features and parameterized on a balanced dataset with an equal number of stabilizing and destabilizing mutations. A comprehensive comparison of the predictive performance of PremPS with other available methods on nine benchmark datasets confirms that our approach consistently outperforms other methods and shows considerable improvement in estimating the impacts of stabilizing mutations. A protein could have multiple structures available, and if another structure of the same protein is used, the predicted change in stability for structure-based methods might be different. Thus, we further estimated the impact of using different structures on prediction accuracy, and demonstrate that our method performs well across different types of structures except for low-resolution structures and models built based on templates with low sequence identity. PremPS can be used for finding functionally important variants, revealing the molecular mechanisms of functional influences and protein design. PremPS is freely available at https://lilab.jysw.suda.edu.cn/research/PremPS/, which allows to do large-scale mutational scanning and takes about four minutes to perform calculations for a single mutation per protein with ~ 300 residues and requires ~ 0.4 seconds for each additional mutation.

The development of computational methods to accurately predict the impacts of amino acid substitutions on protein stability is of paramount importance for the field of protein design and understanding the roles of missense mutations in disease. However, most of the available methods have very limited predictive accuracy for mutations increasing stability and few could consistently perform well across different test cases. Here we present a new computational approach PremPS, which is capable of predicting the effects of single point mutations on protein stability. PremPS employs only ten evolutionary- and structure-based features and is trained on a symmetrical dataset consisting of the same number of cases of stabilizing and destabilizing mutations. Our method was tested against numerous blind datasets and shows a considerable improvement especially in evaluating the effects of stabilizing mutations, outperforming previously developed methods. PremPS is freely available as a user-friendly web server at http://lilab.jysw.suda.edu.cn/research/PremPS/, which is fast enough to handle the large number of cases.

Limitations and challenges in protein stability prediction upon genome variations: towards future applications in precision medicine

Protein stability predictions are becoming essential in medicine to develop novel immunotherapeutic agents and for drug discovery. Despite the large number of computational approaches for predicting the protein stability upon mutation, there are still critical unsolved problems: 1) the limited number of thermodynamic measurements for proteins provided by current databases; 2) the large intrinsic variability of ΔΔG values due to different experimental conditions; 3) biases in the development of predictive methods caused by ignoring the anti-symmetry of ΔΔG values between mutant and native protein forms; 4) over-optimistic prediction performance, due to sequence similarity between proteins used in training and test datasets. Here, we review these issues, highlighting new challenges required to improve current tools and to achieve more reliable predictions. In addition, we provide a perspective of how these methods will be beneficial for designing novel precision medicine approaches for several genetic disorders caused by mutations, such as cancer and neurodegenerative diseases.

Mutational signatures and mutable motifs in cancer genomes

Cancer is a genetic disorder, meaning that a plethora of different mutations, whether somatic or germ line, underlie the etiology of the 'Emperor of Maladies'. Point mutations, chromosomal rearrangements and copy number changes, whether they have occurred spontaneously in predisposed individuals or have been induced by intrinsic or extrinsic (environmental) mutagens, lead to the activation of oncogenes and inactivation of tumor suppressor genes, thereby promoting malignancy. This scenario has now been recognized and experimentally confirmed in a wide range of different contexts. Over the past decade, a surge in available sequencing technologies has allowed the sequencing of whole genomes from liquid malignancies and solid tumors belonging to different types and stages of cancer, giving birth to the new field of cancer genomics. One of the most striking discoveries has been that cancer genomes are highly enriched with mutations of specific kinds. It has been suggested that these mutations can be classified into 'families' based on their mutational signatures. A mutational signature may be regarded as a type of base substitution (e.g. C:G to T:A) within a particular context of neighboring nucleotide sequence (the bases upstream and/or downstream of the mutation). These mutational signatures, supplemented by mutable motifs (a wider mutational context), promise to help us to understand the nature of the mutational processes that operate during tumor evolution because they represent the footprints of interactions between DNA, mutagens and the enzymes of the repair/replication/modification pathways.

Effect of double mutations T790M/L858R on conformation and drug-resistant mechanism of epidermal growth factor receptor explored by molecular dynamics simulations

Epidermal growth factor receptor (EGFR) is one of the most promising targets for the treatment of cancers. Double mutations T790M/L858R lead to different degrees of drug resistance toward inhibitors. In this study, molecular dynamics (MD) simulations followed by principal component analysis are performed to study the conformational changes of EGFR induced by T790M/L858R. The results suggest that T790M/L858R cause obvious disturbance of the structural stability of EGFR. Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) and residue-based free energy decomposition methods are integrated to explore the drug-resistant mechanism of T790M/L858R toward inhibitors. The results show that the decrease in van der Waals interactions of inhibitors with the mutated EFGR relative to the wild-type (WT) one is the main force inducing drug resistance of T790M/L858R toward inhibitors TAK-285, while drug resistance toward W2P and HKI-272 is dominated by the decrease in van der Waals interactions and the increase in polar interactions. We expect that the information obtained from this study can aid rational design of effective drugs to relieve drug resistance of EGFR induced by T790M/L858R.

Computational Approaches to Prioritize Cancer Driver Missense Mutations

Cancer is a complex disease that is driven by genetic alterations. There has been a rapid development of genome-wide techniques during the last decade along with a significant lowering of the cost of gene sequencing, which has generated widely available cancer genomic data. However, the interpretation of genomic data and the prediction of the association of genetic variations with cancer and disease phenotypes still requires significant improvement. Missense mutations, which can render proteins non-functional and provide a selective growth advantage to cancer cells, are frequently detected in cancer. Effects caused by missense mutations can be pinpointed by in silico modeling, which makes it more feasible to find a treatment and reverse the effect. Specific human phenotypes are largely determined by stability, activity, and interactions between proteins and other biomolecules that work together to execute specific cellular functions. Therefore, analysis of missense mutations’ effects on proteins and their complexes would provide important clues for identifying functionally important missense mutations, understanding the molecular mechanisms of cancer progression and facilitating treatment and prevention. Herein, we summarize the major computational approaches and tools that provide not only the classification of missense mutations as cancer drivers or passengers but also the molecular mechanisms induced by driver mutations. This review focuses on the discussion of annotation and prediction methods based on structural and biophysical data, analysis of somatic cancer missense mutations in 3D structures of proteins and their complexes, predictions of the effects of missense mutations on protein stability, protein-protein and protein-nucleic acid interactions, and assessment of conformational changes in protein conformations induced by mutations.

Monitoring of clonal evolution of double C-KIT exon 17 mutations by Droplet Digital PCR in patients with core-binding factor acute myeloid leukemia

C-KIT gene mutations result in the constitutive activation of tyrosine kinase activity, and greatly affect the pathogenesis and prognosis of core-binding factor acute myeloid leukemia (CBF-AML). C-KIT mutations are often found as single point mutations. However, the rate of double mutations has recently increased in AML patients. In this study, we detected six cases (18.8%) harboring double C-KIT exon17 mutations in 75 patients with CBF-AML. The clone composition and dynamic evolution were analyzed by sequencing and droplet digital PCR (ddPCR). Results revealed that these double mutations can be occurred in either the same or different clones. Different clones of double mutations may result in different sensitivity to the treatment of CBF-AML. The clones with N822 mutation responded better to treatment as compared to those with D816 mutation. Moreover, D816 clone was readily transformed into a predominant clone at relapse. Meanwhile, the predominant clones in the same patient may change during the progression of disease. The emerging mutation can originate from a small quantity of clones at diagnosis or newly acquired during the course of disease. Furthermore, patients with double mutations had better overall survival (OS) and event-free survival (EFS) than those with single mutation, but the differences did not reach statistical significance (P > 0.05). The ddPCR is an effective method for monitoring clonal evolution in patients with CBF-AML.

Computational and Experimental Characterization of Patient Derived Mutations Reveal an Unusual Mode of Regulatory Spine Assembly and Drug Sensitivity in EGFR Kinase

The catalytic activation of protein kinases requires precise positioning of key conserved catalytic and regulatory motifs in the kinase core. The Regulatory Spine (RS) is one such structural motif that is dynamically assembled upon kinase activation. The RS is also a mutational hotspot in cancers; however, the mechanisms by which cancer mutations impact RS assembly and kinase activity are not fully understood. In this study, through mutational analysis of patient derived mutations in the RS of EGFR kinase, we identify an activating mutation, M766T, at the RS3 position. RS3 is located in the regulatory αC-helix, and a series of mutations at the RS3 position suggest a strong correlation between the amino acid type present at the RS3 position and ligand (EGF) independent EGFR activation. Small polar amino acids increase ligand independent activity, while large aromatic amino acids decrease kinase activity. M766T relies on the canonical asymmetric dimer for full activation. Molecular modeling and molecular dynamics simulations of WT and mutant EGFR suggest a model in which M766T activates the kinase domain by disrupting conserved autoinhibitory interactions between M766 and hydrophobic residues in the activation segment. In addition, a water mediated hydrogen bond network between T766, the conserved K745-E762 salt bridge, and the backbone amide of the DFG motif is identified as a key determinant of M766T-mediated activation. M766T is resistant to FDA approved EGFR inhibitors such as gefitinib and erlotinib, and computational estimation of ligand binding free energy identifies key residues associated with drug sensitivity. In sum, our studies suggest an unusual mode of RS assembly and oncogenic EGFR activation, and provide new clues for the design of allosteric protein kinase inhibitors.

Annotating Mutational Effects on Proteins and Protein Interactions: Designing Novel and Revisiting Existing Protocols

In this review we describe a protocol to annotate the effects of missense mutations on proteins, their functions, stability, and binding. For this purpose we present a collection of the most comprehensive databases which store different types of sequencing data on missense mutations, we discuss their relationships, possible intersections, and unique features. Next, we suggest an annotation workflow using the state-of-the art methods and highlight their usability, advantages, and limitations for different cases. Finally, we address a particularly difficult problem of deciphering the molecular mechanisms of mutations on proteins and protein complexes to understand the origins and mechanisms of diseases.

Balancing Protein Stability and Activity in Cancer: A New Approach for Identifying Driver Mutations Affecting CBL Ubiquitin Ligase Activation

Oncogenic mutations in the monomeric Casitas B-lineage lymphoma (Cbl) gene have been found in many tumors, but their significance remains largely unknown. Several human c-Cbl (CBL) structures have recently been solved, depicting the protein at different stages of its activation cycle and thus providing mechanistic insight underlying how stability-activity tradeoffs in cancer-related proteins-may influence disease onset and progression. In this study, we computationally modeled the effects of missense cancer mutations on structures representing four stages of the CBL activation cycle to identify driver mutations that affect CBL stability, binding, and activity. We found that recurrent, homozygous, and leukemia-specific mutations had greater destabilizing effects on CBL states than random noncancer mutations. We further tested the ability of these computational models, assessing the changes in CBL stability and its binding to ubiquitin-conjugating enzyme E2, by performing blind CBL-mediated EGFR ubiquitination assays in cells. Experimental CBL ubiquitin ligase activity was in agreement with the predicted changes in CBL stability and, to a lesser extent, with CBL-E2 binding affinity. Two thirds of all experimentally tested mutations affected the ubiquitin ligase activity by either destabilizing CBL or disrupting CBL-E2 binding, whereas about one-third of tested mutations were found to be neutral. Collectively, our findings demonstrate that computational methods incorporating multiple protein conformations and stability and binding affinity evaluations can successfully predict the functional consequences of cancer mutations on protein activity, and provide a proof of concept for mutations in CBL.

ProKinO: a unified resource for mining the cancer kinome

Protein kinases represent a large and diverse family of evolutionarily related proteins that are abnormally regulated in human cancers. Although genome sequencing studies have revealed thousands of variants in protein kinases, translating "big" genomic data into biological knowledge remains a challenge. Here, we describe an ontological framework for integrating and conceptualizing diverse forms of information related to kinase activation and regulatory mechanisms in a machine readable, human understandable form. We demonstrate the utility of this framework in analyzing the cancer kinome, and in generating testable hypotheses for experimental studies. Through the iterative process of aggregate ontology querying, hypothesis generation and experimental validation, we identify a novel mutational hotspot in the αC-β4 loop of the kinase domain and demonstrate the functional impact of the identified variants in epidermal growth factor receptor (EGFR) constitutive activity and inhibitor sensitivity. We provide a unified resource for the kinase and cancer community, ProKinO, housed at http://vulcan.cs.uga.edu/prokino.

'Latent drivers' expand the cancer mutational landscape

A major challenge facing the community involves identification of mutations that drive cancer. Analyses of cancer genomes to detect, and distinguish, 'driver' from 'passenger' mutations are daunting tasks. Here we suggest that there is a third 'latent driver' category. 'Latent driver' mutations behave as passengers, and do not confer a cancer hallmark. However, coupled with other emerging mutations, they drive cancer development and drug resistance. 'Latent drivers' emerge prior to and during cancer evolution. These allosteric mutations can work through 'AND' all-or-none or incremental 'Graded' logic gate mechanisms. Current diagnostic platforms generally assume that actionable 'driver' mutations are those appearing most frequently in cancer. We propose that 'latent driver' detection may help forecast cancer progression and modify personalized drug regimes.

Structural signature of the G719S-T790M double mutation in the EGFR kinase domain and its response to inhibitors

Some individuals with non-small-cell lung cancer (NSCLC) benefit from therapies targeting epidermal growth factor receptor (EGFR), and the characterization of a new mechanism of resistance to the EGFR-specific antibody gefitinib will provide valuable insight into how therapeutic strategies might be designed to overcome this particular resistance mechanism. The G719S and T790M mutations and their combination were involved in causing different conformational redistribution of EGFR. In the present computational study, we analyzed the impact and structural influence of G719S/T790M double mutation (DM) in EGFR with ligand (gefitinib) through molecular dynamic simulation (50 ns) and docking analysis. We observed the escalation in distance between the functional loop and activation loop with respect to T790M mutation compared to the G719S mutation. Furthermore, we confirmed that the G719S mutation causes the ligand to move closer to the hinge region, whereas T790M makes the ligand escape from the binding pocket. Obtained results provide with an explanation for the resistance induced by T790M and a vital clue for the design of drugs to combat gefitinib resistance.

Prediction and prioritization of rare oncogenic mutations in the cancer Kinome using novel features and multiple classifiers

Cancer is a genetic disease that develops through a series of somatic mutations, a subset of which drive cancer progression. Although cancer genome sequencing studies are beginning to reveal the mutational patterns of genes in various cancers, identifying the small subset of “causative” mutations from the large subset of “non-causative” mutations, which accumulate as a consequence of the disease, is a challenge. In this article, we present an effective machine learning approach for identifying cancer-associated mutations in human protein kinases, a class of signaling proteins known to be frequently mutated in human cancers. We evaluate the performance of 11 well known supervised learners and show that a multiple-classifier approach, which combines the performances of individual learners, significantly improves the classification of known cancer-associated mutations. We introduce several novel features related specifically to structural and functional characteristics of protein kinases and find that the level of conservation of the mutated residue at specific evolutionary depths is an important predictor of oncogenic effect. We consolidate the novel features and the multiple-classifier approach to prioritize and experimentally test a set of rare unconfirmed mutations in the epidermal growth factor receptor tyrosine kinase (EGFR). Our studies identify T725M and L861R as rare cancer-associated mutations inasmuch as these mutations increase EGFR activity in the absence of the activating EGF ligand in cell-based assays.

Cancer progresses by accumulation of mutations in a subset of genes that confer growth advantage. The 518 protein kinase genes encoded in the human genome, collectively called the kinome, represent one of the largest families of oncogenes. Targeted sequencing studies of many different cancers have shown that the mutational landscape comprises both cancer-causing “driver” mutations and harmless “passenger” mutations. While the frequent recurrence of some driver mutations in human cancers helps distinguish them from the large number of passenger mutations, a significant challenge is to identify the rare “driver” mutations that are less frequently observed in patient samples and yet are causative. Here we combine computational and experimental approaches to identify rare cancer-associated mutations in Epidermal Growth Factor receptor kinase (EGFR), a signaling protein frequently mutated in cancers. Specifically, we evaluate a novel multiple-classifier approach and features specific to the protein kinase super-family in distinguishing known cancer-associated mutations from benign mutations. We then apply the multiple classifier to identify and test the functional impact of rare cancer-associated mutations in EGFR. We report, for the first time, that the EGFR mutations T725M and L861R, which are infrequently observed in cancers, constitutively activate EGFR in a manner analogous to the frequently observed driver mutations.

A graph theoretic approach to utilizing protein structure to identify non-random somatic mutations

Background

It is well known that the development of cancer is caused by the accumulation of somatic mutations within the genome. For oncogenes specifically, current research suggests that there is a small set of "driver" mutations that are primarily responsible for tumorigenesis. Further, due to recent pharmacological successes in treating these driver mutations and their resulting tumors, a variety of approaches have been developed to identify potential driver mutations using methods such as machine learning and mutational clustering. We propose a novel methodology that increases our power to identify mutational clusters by taking into account protein tertiary structure via a graph theoretical approach.

Results

We have designed and implemented GraphPAC (Graph Protein Amino acid Clustering) to identify mutational clustering while considering protein spatial structure. Using GraphPAC, we are able to detect novel clusters in proteins that are known to exhibit mutation clustering as well as identify clusters in proteins without evidence of prior clustering based on current methods. Specifically, by utilizing the spatial information available in the Protein Data Bank (PDB) along with the mutational data in the Catalogue of Somatic Mutations in Cancer (COSMIC), GraphPAC identifies new mutational clusters in well known oncogenes such as EGFR and KRAS. Further, by utilizing graph theory to account for the tertiary structure, GraphPAC discovers clusters in DPP4, NRP1 and other proteins not identified by existing methods. The R package is available at: http://bioconductor.org/packages/release/bioc/html/GraphPAC.html.

Conclusion

GraphPAC provides an alternative to iPAC and an extension to current methodology when identifying potential activating driver mutations by utilizing a graph theoretic approach when considering protein tertiary structure.

Predicting the Impact of Missense Mutations on Protein-Protein Binding Affinity

The crucial prerequisite for proper biological function is the protein's ability to establish highly selective interactions with macromolecular partners. A missense mutation that alters the protein binding affinity may cause significant perturbations or complete abolishment of the function, potentially leading to diseases. The availability of computational methods to evaluate the impact of mutations on protein-protein binding is critical for a wide range of biomedical applications. Here, we report an efficient computational approach for predicting the effect of single and multiple missense mutations on protein-protein binding affinity. It is based on a well-tested simulation protocol for structure minimization, modified MM-PBSA and statistical scoring energy functions with parameters optimized on experimental sets of several thousands of mutations. Our simulation protocol yields very good agreement between predicted and experimental values with Pearson correlation coefficients of 0.69 and 0.63 and root-mean-square errors of 1.20 and 1.90 kcal mol^-1 for single and multiple mutations, respectively. Compared with other available methods, our approach achieves high speed and prediction accuracy and can be applied to large datasets generated by modern genomics initiatives. In addition, we report a crucial role of water model and the polar solvation energy in estimating the changes in binding affinity. Our analysis also reveals that prediction accuracy and effect of mutations on binding strongly depends on the type of mutation and its location in a protein complex.

The structural pathway of interleukin 1 (IL-1) initiated signaling reveals mechanisms of oncogenic mutations and SNPs in inflammation and cancer

Interleukin-1 (IL-1) is a large cytokine family closely related to innate immunity and inflammation. IL-1 proteins are key players in signaling pathways such as apoptosis, TLR, MAPK, NLR and NF-κB. The IL-1 pathway is also associated with cancer, and chronic inflammation increases the risk of tumor development via oncogenic mutations. Here we illustrate that the structures of interfaces between proteins in this pathway bearing the mutations may reveal how. Proteins are frequently regulated via their interactions, which can turn them ON or OFF. We show that oncogenic mutations are significantly at or adjoining interface regions, and can abolish (or enhance) the protein-protein interaction, making the protein constitutively active (or inactive, if it is a repressor). We combine known structures of protein-protein complexes and those that we have predicted for the IL-1 pathway, and integrate them with literature information. In the reconstructed pathway there are 104 interactions between proteins whose three dimensional structures are experimentally identified; only 15 have experimentally-determined structures of the interacting complexes. By predicting the protein-protein complexes throughout the pathway via the PRISM algorithm, the structural coverage increases from 15% to 71%. In silico mutagenesis and comparison of the predicted binding energies reveal the mechanisms of how oncogenic and single nucleotide polymorphism (SNP) mutations can abrogate the interactions or increase the binding affinity of the mutant to the native partner. Computational mapping of mutations on the interface of the predicted complexes may constitute a powerful strategy to explain the mechanisms of activation/inhibition. It can also help explain how an oncogenic mutation or SNP works.

Activating mutations cluster in the "molecular brake" regions of protein kinases and do not associate with conserved or catalytic residues

Mutations leading to activation of proto-oncogenic protein kinases (PKs) are a type of drivers crucial for understanding tumorogenesis and as targets for antitumor drugs. However, bioinformatics tools so far developed to differentiate driver mutations, typically based on conservation considerations, systematically fail to recognize activating mutations in PKs. Here, we present the first comprehensive analysis of the 407 activating mutations described in the literature, which affect 41 PKs. Unexpectedly, we found that these mutations do not associate with conserved positions and do not directly affect ATP binding or catalytic residues. Instead, they cluster around three segments that have been demonstrated to act, in some PKs, as "molecular brakes" of the kinase activity. This finding led us to hypothesize that an auto inhibitory mechanism mediated by such "brakes" is present in all PKs and that the majority of activating mutations act by releasing it. Our results also demonstrate that activating mutations of PKs constitute a distinct group of drivers and that specific bioinformatics tools are needed to identify them in the numerous cancer sequencing projects currently underway. The clustering in three segments should represent the starting point of such tools, a hypothesis that we tested by identifying two somatic mutations in EPHA7 that might be functionally relevant.

Molecular mechanisms of disease-causing missense mutations

Genetic variations resulting in a change of amino acid sequence can have a dramatic effect on stability, hydrogen bond network, conformational dynamics, activity and many other physiologically important properties of proteins. The substitutions of only one residue in a protein sequence, so-called missense mutations, can be related to many pathological conditions and may influence susceptibility to disease and drug treatment. The plausible effects of missense mutations range from affecting the macromolecular stability to perturbing macromolecular interactions and cellular localization. Here we review the individual cases and genome-wide studies that illustrate the association between missense mutations and diseases. In addition, we emphasize that the molecular mechanisms of effects of mutations should be revealed in order to understand the disease origin. Finally, we report the current state-of-the-art methodologies that predict the effects of mutations on protein stability, the hydrogen bond network, pH dependence, conformational dynamics and protein function.

Effects of oncogenic mutations on the conformational free-energy landscape of EGFR kinase

Activating mutations in the epidermal growth factor receptor (EGFR) tyrosine kinase are frequently found in many cancers. It has been suggested that changes in the equilibrium between its active and inactive conformations are linked to its oncogenic potential. Here, we quantify the effects of some of the most common single (L858R and T790M) and double (T790M-L858R) oncogenic mutations on the conformational free-energy landscape of the EGFR kinase domain by using massive molecular dynamics simulations together with parallel tempering, metadynamics, and one of the best force-fields available. Whereas the wild-type EGFR catalytic domain monomer is mostly found in an inactive conformation, our results show a clear shift toward the active conformation for all of the mutants. The L858R mutation stabilizes the active conformation at the expense of the inactive conformation and rigidifies the αC-helix. The T790M gatekeeper mutant favors activation by stabilizing a hydrophobic cluster. Finally, T790M with L858R shows a significant positive epistasis effect. This combination not only stabilizes the active conformation, but in nontrivial ways changes the free-energy landscape lowering the transition barriers.