Deciphering the Structural Effects of Activating EGFR Somatic Mutations with Molecular Dynamics Simulation

Numerous somatic mutations occurring in the epidermal growth factor receptor (EGFR) family (ErbB) of receptor tyrosine kinases (RTK) have been reported from cancer patients, although relatively few have been tested and shown to cause functional changes in ErbBs. The ErbB receptors are dimerized and activated upon ligand binding, and dynamic conformational changes of the receptors are inherent for induction of downstream signaling. For two mutations shown experimentally to alter EGFR function, A702V and the Delta(746)ELREA(750) deletion mutation, we illustrate in the following protocol how molecular dynamics (MD) simulations can probe the (1) conformational stability of the mutant tyrosine kinase structure in comparison with wild-type EGFR; (2) structural consequences and conformational transitions and their relationship to observed functional changes; (3) effects of mutations on the strength of binding ATP as well as for binding between the kinase domains in the activated asymmetric dimer; and (4) effects of the mutations on key interactions within the EGFR binding site associated with the activated enzyme. The protocol provides a detailed step-by-step procedure as well as guidance that can be more generally useful for investigation of protein structures using MD simulations as a means to probe structural dynamics and the relationship to biological function

Structural Basis for the Functional Changes by EGFR Exon 20 Insertion Mutations

Simple SummaryNon-small cell lung cancer (NSCLC) is the most common type of lung cancer that claims the lives of many worldwide. Activating mutations occurring on the epidermal growth factor receptor (EGFR) protein have been associated with the pathogenesis of NSCLC, among which exon 20 insertion mutations play a significant role. The objective of this study is to examine the dynamic structural changes occurring on the EGFR protein as a result of two common EGFR exon 20 insertion mutations, V769insASV and D770insNPG. The study further aims to uncover the mechanisms by which the insertion mutations increase kinase activity. Our results suggest that the insertion mutations stabilize structural elements key to maintaining the active EGFR conformation. Furthermore, the insertions disrupt an interaction essential in stabilizing the inactive conformation, which could drive the kinase from an inactive to an active EGFR state. AbstractActivating somatic mutations of the epidermal growth factor receptor (EGFR) are frequently implicated in non-small cell lung cancer (NSCLC). While L858R and exon 19 deletion mutations are most prevalent, exon 20 insertions are often observed in NSCLC. Here, we investigated the structural implications of two common EGFR exon 20 insertions in NSCLC, V769insASV and D770insNPG. The active and inactive conformations of wild-type, D770insNPG and V769insASV EGFRs were probed with molecular dynamics simulations to identify local and global alterations that the mutations exert on the EGFR kinase domain, highlighting mechanisms for increased enzymatic activity. In the active conformation, the mutations increase interactions that stabilize the alpha C helix that is essential for EGFR activity. Moreover, the key Lys745-Glu762 salt bridge was more conserved in the insertion mutations. The mutants also preserved the state of the structurally critical aspartate-phenylalanine-glycine (DFG)-motif and regulatory spine (R-spine), which were altered in wild-type EGFR. The insertions altered the structure near the ATP-binding pocket, e.g., the P-loop, which may be a factor for the clinically observed tyrosine kinase inhibitor (TKI) insensitivity by the insertion mutants. The inactive state simulations also showed that the insertions disrupt the Ala767-Arg776 interaction that is key for maintaining the "alpha C-out" inactive conformation, which could consequently fuel the transition from the inactive towards the active EGFR state.</div

Identification of Predictive ERBB Mutations by Leveraging Publicly Available Cell Line Databases

While targeted therapies can be effective for a subgroup of patients, identification of individuals who benefit from the treatments is challenging. At the same time, the predictive significance of the vast majority of the thousands of mutations observed in the cancer tissues remains unknown. Here, we describe the identification of novel predictive biomarkers for ERBB-targeted tyrosine kinase inhibitors (TKI) by leveraging the genetic and drug screening data available in the public cell line databases: Cancer Cell Line Encyclopedia (CCLE), Genomics of Drug Sensitivity in Cancer (GDSC), and Cancer Therapeutics Response Portal (CTRP). We assessed the potential of 412 ERBB mutations in 296 cell lines to predict responses to 10 different ERBB-targeted TKIs. Seventy-six ERBB mutations were identified that were associated with ERBB TKI sensitivity comparable to non-small cell lung cancer cell lines harboring the well-established predictive EGFR L858R mutation or exon 19 deletions. Fourteen (18.4 %) of these mutations were classified as oncogenic by the cBioPortal database, whereas 62 (81.6 %) were regarded as novel potentially predictive mutations. Out of nine functionally validated novel mutations, EGFR Y1069C and ERBB2 E936K were transforming in Ba/F3 cells and demonstrated enhanced signaling activity. Mechanistically, the EGFR Y1069C mutation disrupted the binding of the ubiquitin ligase c-CBL to EGFR, whereas the ERBB2 E936K mutation selectively enhanced the activity of ERBB heterodimers. These findings indicate that integrating data from publicly available cell line databases can be used to identify novel, predictive non-hotspot mutations, potentially expanding the patient population benefiting from existing cancer therapies

Immunogenic SARS-CoV-2 Epitopes: In Silico Study Towards Better Understanding of COVID-19 Disease-Paving the Way for Vaccine Development

The emergence of the COVID-19 outbreak at the end of 2019, caused by the novel coronavirus SARS-CoV-2, has, to date, led to over 13.6 million infections and nearly 600,000 deaths. Consequently, there is an urgent need to better understand the molecular factors triggering immune defense against the virus and to develop countermeasures to hinder its spread. Using in silico analyses, we showed that human major histocompatibility complex (MHC) class I cell-surface molecules vary in their capacity for binding different SARS-CoV-2-derived epitopes, i.e., short sequences of 8-11 amino acids, and pinpointed five specific SARS-CoV-2 epitopes that are likely to be presented to cytotoxic T-cells and hence activate immune responses. The identified epitopes, each one of nine amino acids, have high sequence similarity to the equivalent epitopes of SARS-CoV virus, which are known to elicit an effective T cell response in vitro. Moreover, we give a structural explanation for the binding of SARS-CoV-2-epitopes to MHC molecules. Our data can help us to better understand the differences in outcomes of COVID-19 patients and may aid the development of vaccines against SARS-CoV-2 and possible future outbreaks of novel coronaviruses

Structural characterization of EGFR exon 19 deletion mutation using molecular dynamics simulation.

Epidermal growth factor receptor (EGFR) is a tyrosine kinase receptor important in diverse biological processes including cell proliferation and survival. Upregulation of EGFR activity due to over-expression or mutation is widely implicated in cancer. Activating somatic mutations of the EGFR kinase are postulated to affect the conformation and/or stability of the protein, shifting the EGFR inactive-active state equilibrium towards the activated state. Here, we examined a common EGFR deletion mutation, Δ746ELREA750, which is frequently observed in non-small cell lung cancer patients. By using molecular dynamics simulation, we investigated the structural effects of the mutation that lead to the experimentally reported increases in kinase activity. Simulations of the active form wild-type and ΔELREA EGFRs revealed the deletion stabilizes the αC helix of the kinase domain, which is located adjacent to the deletion site, by rigidifying the flexible β3-αC loop that accommodates the ELREA sequence. Consequently, the αC helix is stabilized in the "αC-in" active conformation that would prolong the time of the activated state. Moreover, in the mutant kinase, a salt bridge between E762 and K745, which is key for EGFR activity, was also stabilized during the simulation. Additionally, the interaction between EGFR and ATP was favored by ΔELREA EGFR over wild-type EGFR, as reflected by the number of hydrogen bonds formed and the free energy of binding. Simulation of inactive EGFR suggested the deletion would promote a shift from the inactive conformation towards active EGFR, which is supported by the inward movement of the αC helix. The MDS results also align with the effects of tyrosine kinase inhibitors on ΔELREA and wild-type EGFR lung cancer cell lines, where more pronounced inhibition was observed against ΔELREA than for wild-type EGFR by inhibitors recognizing the active kinase conformation

Discovery of Compounds Inhibiting the ADP-Ribosyltransferase Activity of Pertussis Toxin

The targeted pathogen-selective approach to drug development holds promise to minimize collateral damage to the beneficial microbiome. The AB$_{5}$-topology pertussis toxin (PtxS1-S5) is a major virulence factor of Bordetella pertussis, the causative agent of the highly contagious respiratory disease whooping cough. Once internalized into the host cell, PtxS1 ADP-ribosylates α-subunits of the heterotrimeric Gαi-superfamily, thereby disrupting G-protein-coupled receptor signaling. Here, we report the discovery of the first small molecules inhibiting the ADP-ribosyltransferase activity of pertussis toxin. We developed protocols to purify milligram-levels of active recombinant B. pertussis PtxS1 from Escherichia coli and an in vitro high throughput-compatible assay to quantify NAD$^{+}$ consumption during PtxS1-catalyzed ADP-ribosylation of Gαi. Two inhibitory compounds (NSC228155 and NSC29193) with low micromolar IC$_{50}$-values (3.0 μM and 6.8 μM) were identified in the in vitro NAD$^{+}$ consumption assay that also were potent in an independent in vitro assay monitoring conjugation of ADP-ribose to Gαi. Docking and molecular dynamics simulations identified plausible binding poses of NSC228155 and in particular of NSC29193, most likely owing to the rigidity of the latter ligand, at the NAD$^{+}$-binding pocket of PtxS1. NSC228155 inhibited the pertussis AB$_{5}$ holotoxin-catalyzed ADP-ribosylation of Gαi in living human cells with a low micromolar IC$_{50}$-value (2.4 μM). NSC228155 and NSC29193 might prove to be useful hit compounds in targeted B. pertussis-selective drug development

Discovery of compounds inhibiting the ADP-ribosyltransferase activity of pertussis toxin

Abstract The targeted pathogen-selective approach to drug development holds promise to minimize collateral damage to the beneficial microbiome. The AB₅-topology pertussis toxin (PtxS1-S5) is a major virulence factor of Bordetella pertussis, the causative agent of the highly contagious respiratory disease whooping cough. Once internalized into the host cell, PtxS1 ADP-ribosylates α-subunits of the heterotrimeric Gαi-superfamily, thereby disrupting G-protein-coupled receptor signaling. Here, we report the discovery of the first small molecules inhibiting the ADP-ribosyltransferase activity of pertussis toxin. We developed protocols to purify milligram-levels of active recombinant B. pertussis PtxS1 from Escherichia coli and an in vitro high throughput-compatible assay to quantify NAD⁺ consumption during PtxS1-catalyzed ADP-ribosylation of Gαi. Two inhibitory compounds (NSC228155 and NSC29193) with low micromolar IC₅₀-values (3.0 μM and 6.8 μM) were identified in the in vitro NAD⁺ consumption assay that also were potent in an independent in vitro assay monitoring conjugation of ADP-ribose to Gαi. Docking and molecular dynamics simulations identified plausible binding poses of NSC228155 and in particular of NSC29193, most likely owing to the rigidity of the latter ligand, at the NAD⁺-binding pocket of PtxS1. NSC228155 inhibited the pertussis AB₅ holotoxin-catalyzed ADP-ribosylation of Gαi in living human cells with a low micromolar IC₅₀-value (2.4 μM). NSC228155 and NSC29193 might prove to be useful hit compounds in targeted B. pertussis-selective drug development

An unbiased in vitro screen for activating epidermal growth factor receptor mutations

Abstract Cancer tissues harbor thousands of mutations, and a given oncogene may be mutated at hundreds of sites, yet only a few of these mutations have been functionally tested. Here, we describe an unbiased platform for the functional characterization of thousands of variants of a single receptor tyrosine kinase (RTK) gene in a single assay. Our in vitro screen for activating mutations (iSCREAM) platform enabled rapid analysis of mutations conferring gain-of-function RTK activity promoting clonal growth. The screening strategy included a somatic model of cancer evolution and utilized a library of 7,216 randomly mutated epidermal growth factor receptor (EGFR) single-nucleotide variants that were tested in murine lymphoid Ba/F3 cells. These cells depend on exogenous interleukin-3 (IL-3) for growth, but this dependence can be compensated by ectopic EGFR overexpression, enabling selection for gain-of-function EGFR mutants. Analysis of the enriched mutants revealed EGFR A702V, a novel activating variant that structurally stabilized the EGFR kinase dimer interface and conferred sensitivity to kinase inhibition by afatinib. As proof of concept for our approach, we recapitulated clinical observations and identified the EGFR L858R as the major enriched EGFR variant. Altogether, iSCREAM enabled robust enrichment of 21 variants from a total of 7,216 EGFR mutations. These findings indicate the power of this screening platform for unbiased identification of activating RTK variants that are enriched under selection pressure in a model of cancer heterogeneity and evolution