Oncogenic Mutations Counteract Intrinsic Disorder in the EGFR Kinase and Promote Receptor Dimerization

OpenStructure: an integrated software framework for computational structural biology

Research projects in structural biology increasingly rely on combinations of heterogeneous sources of information, e.g. evolutionary information from multiple sequence alignments, experimental evidence in the form of density maps and proximity constraints from proteomics experiments. The OpenStructure software framework, which allows the seamless integration of information of different origin, has previously been introduced. The software consists of C++ libraries which are fully accessible from the Python programming language. Additionally, the framework provides a sophisticated graphics module that interactively displays molecular structures and density maps in three dimensions. In this work, the latest developments in the OpenStructure framework are outlined. The extensive capabilities of the framework will be illustrated using short code examples that show how information from molecular-structure coordinates can be combined with sequence data and/or density maps. The framework has been released under the LGPL version 3 license and is available for download from http://www.openstructure.org.

Transitions to catalytically inactive conformations in EGFR kinase

The epidermal growth factor receptor (EGFR) is a key protein in cellular signaling, and its kinase domain (EGFR kinase) is an intensely pursued target of small-molecule drugs. Although both catalytically active and inactive conformations of EGFR kinase have been resolved crystallographically, experimental characterization of the transitions between these conformations remains difficult. Using unbiased, all-atom molecular dynamics simulations, we observed EGFR kinase spontaneously transition from the active to the so-called "Src-like inactive" conformation by way of two sets of intermediate conformations: One corresponds to a previously identified locally disordered state and the other to previously undescribed "extended" conformations, marked by the opening of the ATP-binding site between the two lobes of the kinase domain. We also simulated the protonation-dependent transition of EGFR kinase to another ["Asp-Phe-Gly-out" ("DFG-out")] inactive conformation and observed similar intermediate conformations. A key element observed in the simulated transitions is local unfolding, or "cracking," which supports a prediction of energy landscape theory. We used hydrogen-deuterium (H/D) exchange measurements to corroborate our simulations and found that the simulated intermediate conformations correlate better with the H/D exchange data than existing active or inactive EGFR kinase crystal structures. The intermediate conformations revealed by our simulations of the transition process differ significantly from the existing crystal structures and may provide unique possibilities for structure-based drug discovery.

Architecture and membrane interactions of the EGF receptor

Dimerization-driven activation of the intracellular kinase domains of the epidermal growth factor receptor (EGFR) upon extracellular ligand binding is crucial to cellular pathways regulating proliferation, migration, and differentiation. Inactive EGFR can exist as both monomers and dimers, suggesting that the mechanism regulating EGFR activity may be subtle. The membrane itself may play a role but creates substantial difficulties for structural studies. Our molecular dynamics simulations of membrane-embedded EGFR suggest that, in ligand-bound dimers, the extracellular domains assume conformations favoring dimerization of the transmembrane helices near their N termini, dimerization of the juxtamembrane segments, and formation of asymmetric (active) kinase dimers. In ligand-free dimers, by holding apart the N termini of the transmembrane helices, the extracellular domains instead favor C-terminal dimerization of the transmembrane helices, juxtamembrane segment dissociation and membrane burial, and formation of symmetric (inactive) kinase dimers. Electrostatic interactions of EGFR's intracellular module with the membrane are critical in maintaining this coupling.

From systems to structure: bridging networks and mechanism

There is a wide gap between the generation of large-scale biological data sets and more-detailed, structural and mechanistic studies. However, recent studies that explicitly combine data from systems and structural biological approaches are having a profound effect on our ability to predict how mutations and small molecules affect atomic-level mechanisms, disrupt systems-level networks, and ultimately lead to changes in organismal fitness. In fact, we argue that a shared framework for analysis of nonadditive genetic and thermodynamic responses to perturbations will accelerate the integration of reductionist and global approaches. A stronger bridge between these two areas will allow for a deeper and more-complete understanding of complex biological phenomenon and ultimately provide needed breakthroughs in biomedical research.

The contribution of intrinsic disorder prediction to the elucidation of protein function

Current knowledge of protein function stems primarily from the detailed characterization of some protein domain families and from extensive homology-based annotation transfers. Intrinsically disordered proteins exploit complementary molecular recognition mechanisms and their study has potential to make a significant contribution to the field. Building on established computational biology techniques, efforts are now aimed at unveiling their biological roles at the molecular and system level. Recent developments towards identifying functional sites in disordered regions, and analysing how the occurrence of protein disorder in a gene can provide insight into its biochemical and cellular function are discussed.

The Arabidopsis thaliana SERK1 kinase domain spontaneously refolds to an active state in vitro

Auto-phosphorylating kinase activity of plant leucine-rich-repeat receptor-like kinases (LRR-RLK's) needs to be under tight negative control to avoid unscheduled activation. One way to achieve this would be to keep these kinase domains as intrinsically disordered protein (IDP) during synthesis and transport to its final location. Subsequent folding, which may depend on chaperone activity or presence of interaction partners, is then required for full activation of the kinase domain. Bacterially produced SERK1 kinase domain was previously shown to be an active Ser/Thr kinase. SERK1 is predicted to contain a disordered region in kinase domains X and XI. Here, we show that loss of structure of the SERK1 kinase domain during unfolding is intimately linked to loss of activity. Phosphorylation of the SERK1 kinase domain neither changes its structure nor its stability. Unfolded SERK1 kinase has no autophosphorylation activity and upon removal of denaturant about one half of the protein population spontaneously refolds to an active protein in vitro. Thus, neither chaperones nor interaction partners are required during folding of this protein to its catalytically active state.

Molecular dynamics simulation of PNPLA3 I148M polymorphism reveals reduced substrate access to the catalytic cavity

A missense mutation I148M in PNPLA3 (patatin-like phospholipase domain-containing 3 protein) is significantly correlated with nonalcoholic fatty liver disease (NAFLD). To glean insights into mutation's effect on enzymatic activity, we performed molecular dynamics simulation and flexible docking studies. Our data show that the size of the substrate-access entry site is significantly reduced in mutants, which limits the access of palmitic acid to the catalytic dyad. Besides, the binding free energy calculations suggest low affinity for substrate to mutant enzyme. The substrate-bound system simulations reveal that the spatial arrangement of palmitic acid is distinct in wild-type from that in mutant. The substrate recognition specificity is lost due to the loop where the I148M mutation was located. Our results provide strong evidence for the mechanism by which I148M affects the enzyme activity and suggest that mediating the dynamics may offer a potential avenue for NAFLD.

Tissue-specific splicing of disordered segments that embed binding motifs rewires protein interaction networks

Alternative inclusion of exons increases the functional diversity of proteins. Among alternatively spliced exons, tissue-specific exons play a critical role in maintaining tissue identity. This raises the question of how tissue-specific protein-coding exons influence protein function. Here we investigate the structural, functional, interaction, and evolutionary properties of constitutive, tissue-specific, and other alternative exons in human. We find that tissue-specific protein segments often contain disordered regions, are enriched in posttranslational modification sites, and frequently embed conserved binding motifs. Furthermore, genes containing tissue-specific exons tend to occupy central positions in interaction networks and display distinct interaction partners in the respective tissues, and are enriched in signaling, development, and disease genes. Based on these findings, we propose that tissue-specific inclusion of disordered segments that contain binding motifs rewires interaction networks and signaling pathways. In this way, tissue-specific splicing may contribute to functional versatility of proteins and increases the diversity of interaction networks across tissues.

Probing the diverse landscape of protein flexibility and binding

Protein flexibility spans a broad spectrum, from highly stable folded to intrinsically disordered states. In this review, we discuss how various techniques, including X-ray crystallography, nuclear magnetic resonance spectroscopy and ensemble-modeling strategies employing various experimental measurements, have enabled detailed structural and dynamic characterizations of proteins in their free and bound states. This has revealed a variety of possible binding scenarios in which flexibility can either decrease or increase upon binding. Furthermore, dynamic free-state ensembles have repeatedly been observed to contain transiently formed conformations that partially or completely resemble bound states. These results demonstrate an intimate connection between protein flexibility and protein interactions and illustrate the huge diversity of structure and dynamics in both free proteins and protein complexes.

EGFR in limbo

The epidermal growth factor receptor (EGFR) is normally activated by ligand-induced dimerization. Oncogenic mutations in EGFR promote activation in a largely ligand-independent manner. Shan et al. uncover a partially disordered state of EGFR kinase, providing evidence that oncogenic mutations counteract this intrinsic structural instability to promote dimerization and aberrant activation.

Temporal resolution of autophosphorylation for normal and oncogenic forms of EGFR and differential effects of gefitinib

Epidermal growth factor receptor (EGFR) is a member of the ErbB family of receptor tyrosine kinases (RTK). EGFR overexpression or mutation in many different forms of cancers has highlighted its role as an important therapeutic target. Gefitinib, the first small molecule inhibitor of EGFR kinase function to be approved for the treatment of nonsmall cell lung cancer (NSCLC) by the FDA, demonstrates clinical activity primarily in patients with tumors that harbor somatic kinase domain mutations in EGFR. Here, we compare wild-type EGFR autophosphorylation kinetics to the L834R (also called L858R) EGFR form, one of the most common mutations in lung cancer patients. Using rapid chemical quench, time-resolved electrospray mass spectrometry (ESI-MS), and Western blot analyses, we examined the order of autophosphorylation in wild-type (WT) and L834R EGFR and the effect of gefitinib (Iressa) on the phosphorylation of individual tyrosines. These studies establish that there is a temporal order of autophosphorylation of key tyrosines involved in downstream signaling for WT EGFR and a loss of order for the oncogenic L834R mutant. These studies also reveal unique signature patterns of drug sensitivity for inhibition of tyrosine autophosphorylation by gefitinib: distinct for WT and oncogenic L834R mutant forms of EGFR. Fluorescence studies show that for WT EGFR the binding affinity for gefitinib is weaker for the phosphorylated protein while for the oncogenic mutant, L834R EGFR, the binding affinity of gefitinib is substantially enhanced and likely contributes to the efficacy observed clinically. This mechanistic information is important in understanding the molecular details underpinning clinical observations as well as to aid in the design of more potent and selective EGFR inhibitors.