Ectodomain orientation, conformational plasticity and oligomerization of ErbB1 receptors investigated by molecular dynamics

A perspective of fluorescence microscopy for cellular structural biology with EGFR as witness

The epidermal growth factor receptor (EGFR) is a poster child for the understanding of receptor behaviour, and of paramount importance to cell function and human health. Cloned almost forty years ago, the interest in EGFR's structure/function relationships remains unabated, not least because changes in oncogenic EGFR mutants are key drivers of the formation of lung and brain tumours. The structure of the assemblies formed by EGFR have been comprehensibly investigated by techniques such as high-resolution X-ray crystallography, NMR and all-atom molecular dynamics (MD) simulations. However, the complexity embedded in the portfolio of EGFR states that are only possible in the physiological environment of cells has often proved refractory to cell-free structural methods. Conversely, some key inroads made by quantitative fluorescence microscopy and super-resolution have depended on exploiting the wealth of structures available. Here, a brief personal perspective is provided on how quantitative fluorescence microscopy and super-resolution methods have cross-fertilised with cell-free-derived EGFR structural information. I primarily discuss areas in which my research group has made a contribution to fill gaps in EGFR's cellular structural biology and towards developing new tools to investigate macromolecular assemblies in cells.

Structure-based approach: molecular insight of pyranocumarins against α-glucosidase through computational studies

α-glucosidase is an enzyme that catalyzes the release of α-glucose molecules through hydrolysis reactions. Regulation of this enzyme can increase sugar levels in type-2 diabetes mellitus (DM) patients. Pyranocoumarin derivatives have been identified as α-glucosidase inhibitors. Through an in silico approach, this work studied the inhibition of three pyranocoumarin compounds against the α-glucosidase at the molecular level. Molecular docking and molecular dynamics simulation were performed to understand the dynamics behavior of pyranocoumarin derivatives against α-glucosidase. The prediction of free binding energy (ΔG _bind) using the Quantum Mechanics/Molecular Mechanics-Generalized Born (QM/MM-GBSA) approach for each system had the following results, PC1-α-Glu: -13.97 kcal mol^-1, PC2-α-Glu: -3.69 kcal mol^-1, and PC3-α-Glu: -13.68 kcal mol^-1. The interaction energy of each system shows that the grid score, ΔG _bind, and ΔG _exp values had a similar correlation, that was PC1-α-Glu > PC3-α-Glu > PC2-α-Glu. Additionally, the decomposition energy analysis (ΔG ^residue _bind) was carried out to find out the contribution of the key binding residue. The results showed that there were 15 key binding residues responsible for stabilizing pyranocumarin binding with criteria of ΔG ^residue _bind < -1.00 kcal mol^-1. The evaluation presented in this work could provide information on the molecular level about the inhibitory efficiency of pyranocoumarin derivatives against a-glucosidase enzyme based on computational studies.

Computationally-guided design and affinity improvement of a protein binder targeting a specific site on HER2

A protein binder with a desired epitope and binding affinity is critical to the development of therapeutic agents. Here we present computationally-guided design and affinity improvement of a protein binder recognizing a specific site on domain IV of human epidermal growth factor receptor 2 (HER2). As a model, a protein scaffold composed of Leucine-rich repeat (LRR) modules was used. We designed protein binders which appear to bind a target site on domain IV using a computational method. Top 10 designs were expressed and tested with binding assays, and a lead with a low micro-molar binding affinity was selected. Binding affinity of the selected lead was further increased by two-orders of magnitude through mutual feedback between computational and experimental methods. The utility and potential of our approach was demonstrated by determining the binding interface of the developed protein binder through its crystal structure in complex with the HER2 domain IV.

In silico and in vitro identification of candidate SIRT1 activators from Indonesian medicinal plants compounds database

Sirtuin 1 (SIRT1) is a class III family of protein histone deacetylases involved in NAD⁺-dependent deacetylation reactions. It has been suggested that SIRT1 activators may have a protective role against type 2 diabetes, the aging process, and inflammation. This study aimed to explore and identify medicinal plant compounds from Indonesian Herbal Database (HerbalDB) that might potentially become a candidate for SIRT1 activators through a combination of in silico and in vitro methods. Two pharmacophore models were developed using co-crystalized ligands that allosterically bind with SIRT1 similar to the putative ligands used by SIRT1 activators. Then, these were used for the virtual screening of HerbalDB. The identified compounds were subjected to molecular docking and 50 ns molecular dynamics simulation. Molecular dynamics simulation was analyzed using MM-GB(PB)SA methods. The compounds identified by these methods were tested in an in vitro study using a SIRT-Glo™ luminescence assay. Virtual screening using structure-based pharmacophores predicted that mulberrin as the best candidate SIRT1 activator. Virtual screening using ligand-based pharmacophores predicted that gartanin, quinidine, and quinine to be the best candidates as SIRT1 activators. The molecular docking studies showed the important residues involved were Ile223 and Ile227 at the allosteric region. The MM-GB(PB)SA calculations confirmed that mulberrin, gartanin, quinidine, quinine showed activity at allosteric region and their EC₅₀ in vitro values are 2.10; 1.79; 1.71; 1.14 μM, respectively. Based on in silico and in vitro study results, mulberin, gartanin, quinidine, and quinine had good activity as SIRT1 activators.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Interactions of the EphA2 Kinase Domain with PIPs in Membranes: Implications for Receptor Function

EphA2 is a member of the receptor tyrosine kinase family. Interactions of the cytoplasmic region of EphA2 with the cell membrane are functionally important and yet remain incompletely characterized. Molecular dynamics simulations combined with biochemical studies reveal the interactions of the transmembrane, juxtamembrane (JM), and kinase domains with the membrane. We describe how the kinase domain is oriented relative to the membrane and how the JM region can modulate this interaction. We highlight the role of phosphatidylinositol phosphates (PIPs) in mediating the interaction of the kinase domain with the membrane and, conversely, how positively charged patches at the kinase surface and in the JM region induce the formation of nanoclusters of PIP molecules in the membrane. Integration of these results with those from previous studies enable computational reconstitution of a near complete EphA2 receptor within a membrane, suggesting a role for receptor-lipid interactions in modulation of EphA2.

Transient Acceleration of Epidermal Growth Factor Receptor Dynamics Produces Higher-Order Signaling Clusters

Cell signaling depends on spatiotemporally regulated molecular interactions. Although the movements of signaling proteins have been analyzed with various technologies, how spatial dynamics influence the molecular interactions that transduce signals is unclear. Here, we developed a single-molecule method to analyze the spatiotemporal coupling between motility, clustering, and signaling. The analysis was performed with the epidermal growth factor receptor (EGFR), which triggers signaling through its dimerization and phosphorylation after association with EGF. Our results show that the few EGFRs isolated in membrane subdomains were released by an EGF-dependent increase in their diffusion area, facilitating molecular associations and producing immobile clusters. Using a two-color single-molecule analysis, we found that the EGF-induced state transition alters the properties of the immobile clusters, allowing them to interact for extended periods with the cytoplasmic protein, GRB2. Our study reveals a novel correlation between this molecular interaction and its mesoscale dynamics, providing the initial signaling node.

Conformational stability of the epidermal growth factor (EGF) receptor as influenced by glycosylation, dimerization and EGF hormone binding

The epidermal growth factor receptor (EGFR) is an important transmembrane glycoprotein kinase involved the initiation or perpetuation of signal transduction cascades within cells. These processes occur after EGFR binds to a ligand [epidermal growth factor (EGF)], thus inducing its dimerization and tyrosine autophosphorylation. Previous publications have highlighted the importance of glycosylation and dimerization for promoting proper function of the receptor and conformation in membranes; however, the effects of these associations on the protein conformational stability have not yet been described. Molecular dynamics simulations were performed to characterize the conformational preferences of the monomeric and dimeric forms of the EGFR extracellular domain upon binding to EGF in the presence and absence of N-glycan moieties. Structural stability analyses revealed that EGF provides the most conformational stability to EGFR, followed by glycosylation and dimerization, respectively. The findings also support that EGF-EGFR binding takes place through a large-scale induced-fitting mechanism. Proteins 2017; 85:561-570. © 2016 Wiley Periodicals, Inc.

Conformational Changes in the Epidermal Growth Factor Receptor: Role of the Transmembrane Domain Investigated by Coarse-Grained MetaDynamics Free Energy Calculations

The epidermal growth factor receptor (EGFR) is a dimeric membrane protein that regulates key aspects of cellular function. Activation of the EGFR is linked to changes in the conformation of the transmembrane (TM) domain, brought about by changes in interactions of the TM helices of the membrane lipid bilayer. Using an advanced computational approach that combines Coarse-Grained molecular dynamics and well-tempered MetaDynamics (CG-MetaD), we characterize the large-scale motions of the TM helices, simulating multiple association and dissociation events between the helices in membrane, thus leading to a free energy landscape of the dimerization process. The lowest energy state of the TM domain is a right-handed dimer structure in which the TM helices interact through the N-terminal small-X₃-small sequence motif. In addition to this state, which is thought to correspond to the active form of the receptor, we have identified further low-energy states that allow us to integrate with a high level of detail a range of previous experimental observations. These conformations may lead to the active state via two possible activation pathways, which involve pivoting and rotational motions of the helices, respectively. Molecular dynamics also reveals correlation between the conformational changes of the TM domains and of the intracellular juxtamembrane domains, paving the way for a comprehensive understanding of EGFR signaling at the cell membrane.

Determining the geometry of oligomers of the human epidermal growth factor family on cells with <10 nm resolution

There is a limited range of methods available to characterize macromolecular organization in cells on length scales from 5-50 nm. We review methods currently available and show the latest results from a new single-molecule localization-based method, fluorophore localization imaging with photobleaching (FLImP), using the epidermal growth factor (EGF) receptor (EGFR) as an example system. Our measurements show that FLImP is capable of achieving spatial resolution in the order of 6 nm.

Structures of the EphA2 Receptor at the Membrane: Role of Lipid Interactions

Ephs are transmembrane receptors that mediate cell-cell signaling. The N-terminal ectodomain binds ligands and enables receptor clustering, which activates the intracellular kinase. Relatively little is known about the function of the membrane-proximal fibronectin domain 2 (FN2) of the ectodomain. Multiscale molecular dynamics simulations reveal that FN2 interacts with lipid bilayers via a site comprising K441, R443, R465, Q462, S464, S491, W467, F490, and P459-461. FN2 preferentially binds anionic lipids, a preference that is reduced in the mutant K441E + R443E. We confirm these results by measuring the binding of wild-type and mutant FN2 domains to lipid vesicles. In simulations of the complete EphA2 ectodomain plus the transmembrane region, we show that FN2 anchors the otherwise flexible ectodomain at the surface of the bilayer. Altogether, our data suggest that FN2 serves a dual function of interacting with anionic lipids and constraining the structure of the EphA2 ectodomain to adopt membrane-proximal configurations.

A tale of the epidermal growth factor receptor: The quest for structural resolution on cells

The challenge of determining the architecture and geometry of oligomers of the epidermal growth factor receptor (EGFR) on the cell surface has been approached using a variety of biochemical and biophysical methods. This review is intended to provide a narrative of how key concepts in the field of EGFR research have evolved over the years, from the origins of the prevalent EGFR signalling dimer hypothesis through to the development and implementation of methods that are now challenging the conventional view. The synergy between X-ray crystallography and cellular fluorescence microscopy has become particularly important, precisely because the results from these two methods diverged and highlighted the complexity of the challenge. We illustrate how developments in super-resolution microscopy are now bridging this gap. Exciting times lie ahead where knowledge of the nature of the complexes can assist with the development of a new generation of anti-cancer drugs.

Enhanced dimerization drives ligand-independent activity of mutant epidermal growth factor receptor in lung cancer

Epidermal growth factor receptor kinase mutations drive oncogenesis, but the molecular mechanism of pathological signal initiation is poorly understood. Using high-resolution microscopy methods, the authors reveal that these kinase mutations induce structural changes in the receptor ectodomain that lead to enhanced, ligand-independent dimerization.

Mutations within the epidermal growth factor receptor (EGFR/erbB1/Her1) are often associated with tumorigenesis. In particular, a number of EGFR mutants that demonstrate ligand-independent signaling are common in non–small cell lung cancer (NSCLC), including kinase domain mutations L858R (also called L834R) and exon 19 deletions (e.g., ΔL747-P753insS), which collectively make up nearly 90% of mutations in NSCLC. The molecular mechanisms by which these mutations confer constitutive activity remain unresolved. Using multiple subdiffraction-limit imaging modalities, we reveal the altered receptor structure and interaction kinetics of NSCLC-associated EGFR mutants. We applied two-color single quantum dot tracking to quantify receptor dimerization kinetics on living cells and show that, in contrast to wild-type EGFR, mutants are capable of forming stable, ligand-independent dimers. Two-color superresolution localization microscopy confirmed ligand-independent aggregation of EGFR mutants. Live-cell Förster resonance energy transfer measurements revealed that the L858R kinase mutation alters ectodomain structure such that unliganded mutant EGFR adopts an extended, dimerization-competent conformation. Finally, mutation of the putative dimerization arm confirmed a critical role for ectodomain engagement in ligand-independent signaling. These data support a model in which dysregulated activity of NSCLC-associated kinase mutants is driven by coordinated interactions involving both the kinase and extracellular domains that lead to enhanced dimerization.

Compounds identified by virtual docking to a tetrameric EGFR extracellular domain can modulate Grb2 internalization

Background

Overexpression or mutation of the epidermal growth factor receptor (EGFR) potently enhances the growth of many solid tumors. Tumor cells frequently display resistance to mechanistically-distinct EGFR-directed therapeutic agents, making it valuable to develop therapeutics that work by additional mechanisms. Current EGFR-targeting therapeutics include antibodies targeting the extracellular domains, and small molecules inhibiting the intracellular kinase domain. Recent studies have identified a novel prone extracellular tetrameric EGFR configuration, which we identify as a potential target for drug discovery.

Methods

Our focus is on the prone EGFR tetramer, which contains a novel protein-protein interface involving extracellular domain III. This EGFR tetramer is computationally targeted for stabilization by small molecule ligand binding. This study performed virtual screening of a Life Chemicals, Inc. small molecule library of 345,232 drug-like compounds against a molecular dynamics simulation of protein-protein interfaces distinct to the novel tetramer. One hundred nine chemically diverse candidate molecules were selected and evaluated using a cell-based high-content imaging screen that directly assessed induced internalization of the EGFR effector protein Grb2. Positive hits were further evaluated for influence on phosphorylation of EGFR and its effector ERK1/2.

Results

Fourteen hit compounds affected internalization of Grb2, an adaptor responsive to EGFR activation. Most hits had limited effect on cell viability, and minimally influenced EGFR and ERK1/2 phosphorylation. Docked hit compound poses generally include Arg270 or neighboring residues, which are also involved in binding the effective therapeutic cetuximab, guiding further chemical optimization.

Conclusions

These data suggest that the EGFR tetrameric configuration offers a novel cancer drug target.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1415-6) contains supplementary material, which is available to authorized users.

Nanometric molecular separation measurements by single molecule photobleaching

Although considerable progress has been made in imaging distances in cells below the diffraction limit using FRET and super-resolution microscopy, methods for determining the separation of macromolecules in the 10-50 nm range have been elusive. We have developed fluorophore localisation imaging with photobleaching (FLImP), based on the quantised bleaching of individual protein-bound dye molecules, to quantitate the molecular separations in oligomers and nanoscale clusters. We demonstrate the benefits of using our method in studying the nanometric organisation of the epidermal growth factor receptor in cells.

Structure-function relationships and supramolecular organization of the EGFR (epidermal growth factor receptor) on the cell surface

Dimerization and higher-order oligomerization are believed to play an important role in the activation of the EGFR (epidermal growth factor receptor). Understanding of the process has been limited by the lack of availability of suitable methods for the measurement in cells of distances in the range 10-100 nm, too short for imaging methods and too long for spectroscopic methods such as FRET. In the present article, we review the current state of our knowledge of EGFR oligomerization, and describe results from a new single-molecule localization method that has allowed the quantitative characterization of the distribution of EGFR-EGFR distances in cells. Recent data suggest the involvement of cortical actin in regulating the formation of EGFR complexes.

Hotspot mutations in KIT receptor differentially modulate its allosterically coupled conformational dynamics: impact on activation and drug sensitivity

Receptor tyrosine kinase KIT controls many signal transduction pathways and represents a typical allosterically regulated protein. The mutation-induced deregulation of KIT activity impairs cellular physiological functions and causes serious human diseases. The impact of hotspots mutations (D816H/Y/N/V and V560G/D) localized in crucial regulatory segments, the juxtamembrane region (JMR) and the activation (A-) loop, on KIT internal dynamics was systematically studied by molecular dynamics simulations. The mutational outcomes predicted in silico were correlated with in vitro and in vivo activation rates and drug sensitivities of KIT mutants. The allosteric regulation of KIT in the native and mutated forms is described in terms of communication between the two remote segments, JMR and A-loop. A strong correlation between the communication profile and the structural and dynamical features of KIT in the native and mutated forms was established. Our results provide new insight on the determinants of receptor KIT constitutive activation by mutations and resistance of KIT mutants to inhibitors. Depiction of an intra-molecular component of the communication network constitutes a first step towards an integrated description of vast communication pathways established by KIT in physiopathological contexts.

Receptor tyrosine kinase KIT plays a crucial role in the regulation of cell signaling. This allosterically controlled activity may be affected by gain-of-function mutations that promote the development of several cancers. Identification of the molecular basis of KIT constitutive activation and allosteric regulation has inspired computational study of KIT hotspot mutations. In the present contribution, we investigated the mutation-induced effects on KIT conformational dynamics and intra-protein communication conditionally on the mutation location and the nature of the substituting amino acid. Our data elucidate that all studied mutations stabilize an inactive non-autoinhibited state of KIT over the inactive auto-inhibited state prevalent for the native protein. This shift in the protein conformational landscape promotes KIT constitutive activation. Our in silico analysis established correlations between the structural and dynamical effects induced by oncogenic mutations and the mutants auto-activation rates and drug sensitivities measured in vitro and in vivo. Particularly, the A-loop mutations stabilize the drug-resistant forms, while the JMR mutations may facilitate inhibitors binding to the active site. Cross-correlations established between local and long-range structural and dynamical effects demonstrate the allosteric character of the gain-of-function mutations mode of action.

Membrane interaction of bound ligands contributes to the negative binding cooperativity of the EGF receptor

The epidermal growth factor receptor (EGFR) plays a key role in regulating cell proliferation, migration, and differentiation, and aberrant EGFR signaling is implicated in a variety of cancers. EGFR signaling is triggered by extracellular ligand binding, which promotes EGFR dimerization and activation. Ligand-binding measurements are consistent with a negatively cooperative model in which the ligand-binding affinity at either binding site in an EGFR dimer is weaker when the other site is occupied by a ligand. This cooperativity is widely believed to be central to the effects of ligand concentration on EGFR-mediated intracellular signaling. Although the extracellular portion of the human EGFR dimer has been resolved crystallographically, the crystal structures do not reveal the structural origin of this negative cooperativity, which has remained unclear. Here we report the results of molecular dynamics simulations suggesting that asymmetrical interactions of the two binding sites with the membrane may be responsible (perhaps along with other factors) for this negative cooperativity. In particular, in our simulations the extracellular domains of an EGFR dimer spontaneously lay down on the membrane in an orientation in which favorable membrane contacts were made with one of the bound ligands, but could not be made with the other. Similar interactions were observed when EGFR was glycosylated, as it is in vivo.

The ErbB/HER family of protein-tyrosine kinases and cancer

The human epidermal growth factor receptor (EGFR) family consists of four members that belong to the ErbB lineage of proteins (ErbB1-4). These receptors consist of a glycosylated extracellular domain, a single hydrophobic transmembrane segment, and an intracellular portion with a juxtamembrane segment, a protein kinase domain, and a carboxyterminal tail. Seven ligands bind to EGFR including epidermal growth factor and transforming growth factor α, none bind to ErbB2, two bind to ErbB3, and seven ligands bind to ErbB4. The ErbB proteins function as homo and heterodimers. The heterodimer consisting of ErbB2, which lacks a ligand, and ErbB3, which is kinase impaired, is surprisingly the most robust signaling complex of the ErbB family. Growth factor binding to EGFR induces a large conformational change in the extracellular domain, which leads to the exposure of a dimerization arm in domain II of the extracellular segment. Two ligand-EGFR complexes unite to form a back-to-back dimer in which the ligands are on opposite sides of the aggregate. Following ligand binding, EGFR intracellular kinase domains form an asymmetric homodimer that resembles the heterodimer formed by cyclin and cyclin-dependent kinase. The carboxyterminal lobe of the activator kinase of the dimer interacts with the amino-terminal lobe of the receiver kinase thereby leading to its allosteric stimulation. Downstream ErbB signaling modules include the phosphatidylinositol 3-kinase/Akt (PKB) pathway, the Ras/Raf/MEK/ERK1/2 pathway, and the phospholipase C (PLCγ) pathway. Several malignancies are associated with the mutation or increased expression of members of the ErbB family including lung, breast, stomach, colorectal, head and neck, and pancreatic carcinomas and glioblastoma (a brain tumor). Gefitinib, erlotinib, and afatinib are orally effective protein-kinase targeted quinazoline derivatives that are used in the treatment of ERBB1-mutant lung cancer. Lapatinib is an orally effective quinazoline derivative used in the treatment of ErbB2-overexpressing breast cancer. Trastuzumab, pertuzumab, and ado-trastuzumab emtansine, which are given intravenously, are monoclonal antibodies that target the extracellular domain and are used for the treatment of ErbB2-positive breast cancer; ado-trastuzumab emtansine is an antibody-drug conjugate that delivers a cytotoxic drug to cells overexpressing ErbB2. Cetuximab and panitumumab are monoclonal antibodies that target ErbB1 and are used in the treatment of colorectal cancer. Cancers treated with these targeted drugs eventually become resistant to them. The role of combinations of targeted drugs or targeted drugs with cytotoxic therapies is being explored in an effort to prevent or delay drug resistance in the treatment of these malignancies.

Ligand binding and dynamics of the monomeric epidermal growth factor receptor ectodomain

The ectodomain of the human epidermal growth factor receptor (hEGFR) controls input to several cell signalling networks via binding with extracellular growth factors. To gain insight into the dynamics and ligand binding of the ectodomain, the hEGFR monomer was subjected to molecular dynamics simulation. The monomer was found to be substantially more flexible than the ectodomain dimer studied previously. Simulations where the endogeneous ligand EGF binds to either Subdomain I or Subdomain III, or where hEGFR is unbound, show significant differences in dynamics. The molecular mechanics Poisson-Boltzmann surface area method has been used to derive relative free energies of ligand binding, and we find that the ligand is capable of binding either subdomain with a slight preference for III. Alanine-scanning calculations for the effect of selected ligand mutants on binding reproduce the trends of affinity measurements. Taken together, these results emphasize the possible role of the ectodomain monomer in the initial step of ligand binding, and add details to the static picture obtained from crystal structures.

Exploring higher-order EGFR oligomerisation and phosphorylation--a combined experimental and theoretical approach

The epidermal growth factor receptor (EGFR) kinase is generally considered to be activated by either ligand-induced dimerisation or a ligand-induced conformational change within pre-formed dimers. Ligand-induced higher-order EGFR oligomerisation or clustering has been reported but it is not clear how EGFR oligomers, as distinct from EGFR dimers, influence signaling outputs. To address this question, we combined measures of receptor clustering (microscopy; image correlation spectroscopy) and phosphorylation (Western blots) with modelling of mass-action chemical kinetics. A stable BaF/3 cell-line that contains a high proportion (>90%) of inactive dimers of EGFR-eGFP but no secreted ligand and no other detectable ErbB receptors was used as the model cell system. EGF at concentrations of greater than 1 nM was found to cluster EGFR-eGFP dimers into higher-order complexes and cause parallel increases in EGFR phosphorylation. The kinetics of EGFR clustering and phosphorylation were both rapid, plateauing within 2 minutes after stimulation with 30 nM EGF. A rule-based model was formulated to interpret the data. This model took into account ligand binding, ligand-induced conformational changes in the cytosolic tail, monomer-dimer-trimer-tetramer transitions via ectodomain- and kinase-mediated interactions, and phosphorylation. The model predicts that cyclic EGFR tetramers are the predominant phosphorylated species, in which activated receptor dimers adopt a cyclic side-by-side orientation, and that receptor kinase activation is stabilised by the intramolecular interactions responsible for cyclic tetramerization.

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.

Structural modelling and dynamics of proteins for insights into drug interactions

Proteins are the workhorses of biomolecules and their function is affected by their structure and their structural rearrangements during ligand entry, ligand binding and protein-protein interactions. Hence, the knowledge of protein structure and, importantly, the dynamic behaviour of the structure are critical for understanding how the protein performs its function. The predictions of the structure and the dynamic behaviour can be performed by combinations of structure modelling and molecular dynamics simulations. The simulations also need to be sensitive to the constraints of the environment in which the protein resides. Standard computational methods now exist in this field to support the experimental effort of solving protein structures. This review presents a comprehensive overview of the basis of the calculations and the well-established computational methods used to generate and understand protein structure and function and the study of their dynamic behaviour with the reference to lung-related targets.

Conformational transition and energy landscape of ErbB4 activated by neuregulin1β: one microsecond molecular dynamics simulations

ErbB4, a receptor tyrosine kinase of the ErbB family, plays crucial roles in cell growth and differentiation, especially in the development of the heart and nervous system. Ligand binding to its extracellular region could modulate the activation process. To understand the mechanism of ErbB4 activation induced by ligand binding, we performed one microsecond molecular dynamics (MD) simulations on the ErbB4 extracellular region (ECR) with and without its endogenous ligand neuregulin1β (NRG1β). The conformational transition of the ECR-ErbB4/NRG1β complex from a tethered inactive conformation to an extended active-like form has been observed, while such large and function-related conformational change has not been seen in the simulation on the ECR-ErbB4, suggesting that ligand binding is indeed the active inducing force for the conformational transition and further dimerization. On the basis of MD simulations and principal component analysis, we constructed a rough energy landscape for the conformational transition of ECR-ErbB4/NRG1β complex, suggesting that the conformational change from the inactive state to active-like state involves a stable conformation. The energy barrier for the tether opening was estimated as ~2.7 kcal/mol, which is very close to the experimental value (1-2 kcal/mol) reported for ErbB1. On the basis of the simulation results, an atomic mechanism for the ligand-induced activation of ErbB4 was postulated. The present MD simulations provide a new insight into the conformational changes underlying the activation of ErbB4.

Human epidermal growth factor receptor (HER1) aligned on the plasma membrane adopts key features of Drosophila EGFR asymmetry

Current models suggest that ligand-binding heterogeneity in HER1 [human EGFR (epidermal growth factor receptor] arises from negative co-operativity in signalling HER1 dimers, for which the asymmetry of the extracellular region of the Drosophila EGFR has recently provided a structural basis. However, no asymmetry is apparent in the current crystal structure of the isolated extracellular region of HER1. This receptor also differs from the Drosophila EGFR in that negative co-operativity is found only in full-length receptors in cells. Structural insights into HER1 in epithelial cells, derived from FLIM (fluorescence lifetime imaging microscopy) and two-dimensional FRET (Förster resonance energy transfer) combined with Monte Carlo and molecular dynamics simulations, have demonstrated a high-affinity ligand-binding HER1 conformation consistent with the extracellular region aligned flat on the plasma membrane. This conformation shares key features with that of the Drosophila EGFR, suggesting that the structural basis for negative co-operativity is conserved from invertebrates to humans, but that, in HER1, the extracellular region asymmetry requires interactions with the plasma membrane.

Investigating extracellular in situ EGFR structure and conformational changes using FRET microscopy

The crystallographic structures of functional fragments of ErbBs have provided excellent insights into the geometry of growth factor binding and receptor dimerization. By placing together receptor fragments to build structural models of entire receptors, we expect to understand how these enzymes are allosterically regulated; however, several predictions from these models are inconsistent with experimental evidence from cells. The opening of this gap underlines the need to investigate intact ErbBs by combining cellular and structural studies into a full picture.

Evidence for extended YFP-EGFR dimers in the absence of ligand on the surface of living cells

The epidermal growth factor receptor (EGFR) is a member of the erbB tyrosine kinase family of receptors. Structural studies have revealed two distinct conformations of the ectodomain of the EGFR: a compact, tethered, conformation and an untethered extended conformation. In the context of a monomer-dimer transition model, ligand binding is thought to untether the monomeric receptor leading to exposure of a dimerization arm which then facilitates receptor dimerization, kinase activation and signaling. For receptors directed orthogonal to the local plane of the membrane surface, this would lead to a large change in the distance of the receptor N-terminus from the membrane surface. To investigate this experimentally, we produced stable BaF/3 cell lines expressing a biochemically functional yellow fluorescent protein (YFP)-EGFR chimera and determined the vertical separation of the N-terminal YFP tag from the membrane using fluorescence resonance energy transfer (FRET) techniques. Homo-FRET/rFLIM was employed to determine the presence of unliganded dimers and to measure the average distance between the N-terminal tags in those dimers. The results suggest that EGF-induced activation occurs within or between pre-formed and extended dimers with very little change in the extension of the N-terminii from the membrane surface. These results provide constraints on possible models for EGFR activation.

Optics clustered to output unique solutions: a multi-laser facility for combined single molecule and ensemble microscopy

Optics clustered to output unique solutions (OCTOPUS) is a microscopy platform that combines single molecule and ensemble imaging methodologies. A novel aspect of OCTOPUS is its laser excitation system, which consists of a central core of interlocked continuous wave and pulsed laser sources, launched into optical fibres and linked via laser combiners. Fibres are plugged into wall-mounted patch panels that reach microscopy end-stations in adjacent rooms. This allows multiple tailor-made combinations of laser colours and time characteristics to be shared by different end-stations minimising the need for laser duplications. This setup brings significant benefits in terms of cost effectiveness, ease of operation, and user safety. The modular nature of OCTOPUS also facilitates the addition of new techniques as required, allowing the use of existing lasers in new microscopes while retaining the ability to run the established parts of the facility. To date, techniques interlinked are multi-photon/multicolour confocal fluorescence lifetime imaging for several modalities of fluorescence resonance energy transfer (FRET) and time-resolved anisotropy, total internal reflection fluorescence, single molecule imaging of single pair FRET, single molecule fluorescence polarisation, particle tracking, and optical tweezers. Here, we use a well-studied system, the epidermal growth factor receptor network, to illustrate how OCTOPUS can aid in the investigation of complex biological phenomena.

Mutation D816V alters the internal structure and dynamics of c-KIT receptor cytoplasmic region: implications for dimerization and activation mechanisms

The type III receptor tyrosine kinase (RTK) KIT plays a crucial role in the transmission of cellular signals through phosphorylation events that are associated with a switching of the protein conformation between inactive and active states. D816V KIT mutation is associated with various pathologies including mastocytosis and cancers. D816V-mutated KIT is constitutively active, and resistant to treatment with the anti-cancer drug Imatinib. To elucidate the activating molecular mechanism of this mutation, we applied a multi-approach procedure combining molecular dynamics (MD) simulations, normal modes analysis (NMA) and binding site prediction. Multiple 50-ns MD simulations of wild-type KIT and its mutant D816V were recorded using the inactive auto-inhibited structure of the protein, characteristic of type III RTKs. Computed free energy differences enabled us to quantify the impact of D816V on protein stability in the inactive state. We evidenced a local structural alteration of the activation loop (A-loop) upon mutation, and a long-range structural re-organization of the juxta-membrane region (JMR) followed by a weakening of the interaction network with the kinase domain. A thorough normal mode analysis of several MD conformations led to a plausible molecular rationale to propose that JMR is able to depart its auto-inhibitory position more easily in the mutant than in wild-type KIT and is thus able to promote kinase mutant dimerization without the need for extra-cellular ligand binding. Pocket detection at the surface of NMA-displaced conformations finally revealed that detachment of JMR from the kinase domain in the mutant was sufficient to open an access to the catalytic and substrate binding sites.

Protein kinases are involved in a huge amount of cellular processes through phosphorylation, a crucial mechanism in cell signaling, and their misregulation often results in disease. The deactivation of protein tyrosine kinases (PTKs) or their oncogenic activation arises from mutations which affect the protein primary structure and the configuration of the enzymatic site apparently by stabilizing the activation loop (A-loop) extended conformation. Particularly, mutation D816V of receptor tyrosine kinase (RTK) KIT, found in patients with pediatric mastocytosis, acute leukemia or germ cell tumors, can be considered as the archetype of mutation inducing a displacement of the population equilibrium toward the active conformation. We present a comprehensive computational study of the activating mechanism(s) of this mutation. Our multi-approach in silico procedure evidenced a local alteration of the A-loop structure, and a long-range structural re-organization of the juxta-membrane region (JMR) followed by a weakening of the interaction network with the kinase domain. Our results provided a plausible conception of how the observed departure of JMR from kinase domain in the mutant promotes kinase mutant dimerization without requiring extra-cellular ligand binding. The pocket profiles we obtained suggested putative allosteric binding sites that could be targeted by ligands/modulators that trap the mutated enzyme.

Synergy between trastuzumab and pertuzumab for human epidermal growth factor 2 (Her2) from colocalization: an in silico based mechanism

Introduction

Human epidermal growth factor 2 (Her2), a receptor tyrosine kinase, is overexpressed in breast cancers. It has been successfully targeted by small molecule kinase inhibitors and by antibodies. Recent clinical data show a synergistic response in patients when two antibodies, trastuzumab and pertuzumab, are given in combination.

Methods

This unexpected effect is rationalized through computer models and molecular dynamic simulations by hypothesizing that the two antibodies can co-localize on the same molecule of the Her2 extracellular domain.

Results

Simulations suggest that the clinical synergism observed for the two antibodies arises partly from enhanced affinity that originates in cooperative interactions between these two antibodies when they are co-localized on Her2 and "clamp" it; this may inhibit dimerization and possibly higher oligomerizations with neighboring receptors. In the presence of trastuzumab, the receptor becomes highly plastic, especially domains I to III, and this appears to promote increased association with pertuzumab. Further, the presence of pertuzumab evokes novel interactions between the receptor and trastuzumab. Indeed, splicing out of this region in silico results in a big reduction in the interactions of the antibody with the receptor.

Conclusions

If validated, these findings will bring about a new direction in the design of antibodies whereby different epitopes on the same antibody may be targeted to lead to synergistic/cooperative inhibition and contribute to generate more potent therapeutics and to increase clinical efficacy.

Human epidermal growth factor receptor (EGFR) aligned on the plasma membrane adopts key features of Drosophila EGFR asymmetry

The ability of epidermal growth factor receptor (EGFR) to control cell fate is defined by its affinity for ligand. Current models suggest that ligand-binding heterogeneity arises from negative cooperativity in signaling receptor dimers, for which the asymmetry of the extracellular region of the Drosophila EGFR has recently provided a structural basis. However, no asymmetry is apparent in the isolated extracellular region of the human EGFR. Human EGFR also differs from the Drosophila EGFR in that negative cooperativity is found only in full-length receptors in cells. To gain structural insights into the human EGFR in situ, we developed an approach based on quantitative Förster resonance energy transfer (FRET) imaging, combined with Monte Carlo and molecular dynamics simulations, to probe receptor conformation in epithelial cells. We experimentally demonstrate a high-affinity ligand-binding human EGFR conformation consistent with the extracellular region aligned flat on the plasma membrane. We explored the relevance of this conformation to ligand-binding heterogeneity and found that the asymmetry of this structure shares key features with that of the Drosophila EGFR, suggesting that the structural basis for negative cooperativity is conserved from invertebrates to humans but that in human EGFR the extracellular region asymmetry requires interactions with the plasma membrane.

The tethering arm of the EGF receptor is required for negative cooperativity and signal transduction

The EGF receptor is a classical receptor-tyrosine kinase. In the absence of ligand, the receptor adopts a closed conformation in which the dimerization arm of subdomain II interacts with the tethering arm in subdomain IV. Following the binding of EGF, the receptor opens to form a symmetric, back-to-back dimer. Although it is clear that the dimerization arm of subdomain II is central to the formation of receptor dimers, the role of the tethering arm of subdomain IV (residues 561-585) in this configuration is not known. Here we use (125)I-EGF binding studies to assess the functional role of the tethering arm in the EGF receptor dimer. Mutation of the three major residues that contribute to tethering (D563A,H566A,K585A-EGF receptor) did not significantly alter either the ligand binding properties or the signaling properties of the EGF receptor. By contrast, breaking the Cys(558)-Cys(567) disulfide bond through double alanine replacements or deleting the loop entirely led to a decrease in the negative cooperativity in EGF binding and was associated with small changes in downstream signaling. Deletion of the Cys(571)-Cys(593) disulfide bond abrogated cooperativity, resulting in a high affinity receptor and increased sensitivity of downstream signaling pathways to EGF. Releasing the Cys(571)-Cys(593) disulfide bond resulted in extreme negative cooperativity, ligand-independent kinase activity, and impaired downstream signaling. These data demonstrate that the tethering arm plays an important role in supporting cooperativity in ligand binding. Because cooperativity implies subunit-subunit interactions, these results also suggest that the tethering arm contributes to intersubunit interactions within the EGF receptor dimer.

Large conformational changes in proteins: signaling and other functions

Guanine and adenine nucleotide triphosphatases, such as Ras proteins and protein kinases, undergo large conformational changes upon ligand binding in the course of their functions. New computer simulation methods have combined with experimental studies to deepen our understanding of these phenomena. In particular, a 'conformational selection' picture is emerging, where alterations in the relative populations of pre-existing conformations can best describe the conformational switching activity of these important proteins.