Dataset Information

Interactions of the EGFR juxtamembrane domain with PIP2-containing lipid bilayers: Insights from multiscale molecular dynamics simulations.

ABSTRACT:

Background

The epidermal growth factor receptor (EGFR) is the best characterised member of the receptor tyrosine kinases, which play an important role in signalling across mammalian cell membranes. The EGFR juxtamembrane (JM) domain is involved in the mechanism of activation of the receptor, interacting with the anionic lipid phosphatidylinositol 4,5-bisphosphate (PIP2) in the intracellular leaflet of the cell membrane.

Methods

Multiscale MD simulations were used to characterize PIP2-JM interactions. Simulations of the transmembrane helix plus JM region (TM-JM) dimer (PDB:2M20) in both PIP2-containing and PIP2-depleted lipid bilayer membranes revealed the interactions of the JM with PIP2 and other lipids.

Results

PIP2 forms strong interactions with the basic residues in the R645-R647 motif of the JM domain resulting in clustering of PIP2 around the protein. This association of PIP2 and the JM domain aids stabilization of JM-A dimer away from the membrane. Mutation (R645N/R646N/R647N) or PIP2-depletion results in deformation of the JM-A dimer and changes in JM-membrane interactions.

Conclusions

These simulations support the proposal that the positively charged residues at the start of the JM-A domain stabilize the JM-A helices in an orientation away from the membrane surface through binding to PIP2, thus promoting a conformation corresponding to an asymmetric (i.e. activated) kinase.

General significance

This study indicates that MD simulations may be used to characterise JM/lipid interactions, thus helping to define their role in the mechanisms of receptor tyrosine kinases. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

SUBMITTER: Abd Halim KB

PROVIDER: S-EPMC4547087 | biostudies-literature | 2015 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Interactions of the EGFR juxtamembrane domain with PIP2-containing lipid bilayers: Insights from multiscale molecular dynamics simulations.

Abd Halim Khairul Bariyyah KB Koldsø Heidi H Sansom Mark S P MSP

Biochimica et biophysica acta 20140916 5

<h4>Background</h4>The epidermal growth factor receptor (EGFR) is the best characterised member of the receptor tyrosine kinases, which play an important role in signalling across mammalian cell membranes. The EGFR juxtamembrane (JM) domain is involved in the mechanism of activation of the receptor, interacting with the anionic lipid phosphatidylinositol 4,5-bisphosphate (PIP2) in the intracellular leaflet of the cell membrane.<h4>Methods</h4>Multiscale MD simulations were used to characterize P ...[more]

PMID: 25219456

Similar Datasets

Project description:We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.

Project description:Diphtheria toxin translocation (T) domain undergoes conformational changes in acidic solution and associates with the lipid membranes, followed by refolding and transmembrane insertion of two nonpolar helices. This process is an essential step in delivery of the toxic catalytic domain of the diphtheria toxin to the infected cell, yet its molecular determinants are poorly characterized and understood. Therefore, an atomistic model of the T-domain-membrane interaction is needed to help characterize factors responsible for such association. In this work, we present atomistic model structures of T-domain membrane-bound conformations and investigate structural factors responsible for T-domain affinity with the lipid bilayer in acidic solution using all-atom molecular dynamics (MD) simulations. The initial models of the protein conformations and protein-membrane association that serve as starting points in the present work were developed using atomistic simulations of partial unfolding of the T-domain in acidic solution (Kurnikov, I. V.; et al. J. Mol. Biol. 2013, 425, 2752-2764), and coarse-grained simulations of the T-domain association with the membranes of various compositions (Flores-Canales, J. C.; et al. J. Membr. Biol. 2015, 248, 529-543). In this work we present atomistic level modeling of two distinct configurations of the T-domain in association with the anionic lipid bilayer. In microsecond-long MD simulations both conformations retain their compact structure and gradually penetrate deeper into the bilayer interface. One membrane-bound conformation is stabilized by the protein contacts with the lipid hydrophobic core. The second modeled conformation is initially inserted less deeply and forms multiple contacts with the lipid at the interface (headgroup) region. Such contacts are formed by the charged and hydrophilic groups of partially unfolded terminal helixes and loops. Neutralization of the acidic residues at the membrane interface allows for deeper insertion of the protein and reorientation of the protein at the membrane interface, which corroborates that acidic residue protonation as well as presence of the anionic lipids may play a role in the membrane association and further membrane insertion of the T-domain as implicated in experiments. All simulations reported in this work were performed using AMBER force-field on Anton supercomputer. To perform these reported simulations, we developed and carefully tested a force-field for the anionic 1-palmitoyl-2-oleoyl-phosphatidyl-glycerol (POPG) lipid, compatible with the Amber 99SB force-field and stable in microsecond-long MD simulations in isothermal-isobaric ensemble.

Project description:N-methyl-d-aspartate receptors (NMDARs) are members of the ionotropic glutamate receptor family that mediate excitatory synaptic transmission in the central nervous system. The channels of NMDARs are permeable to Ca2+ but blocked by Mg2+, distinctive properties that underlie essential brain processes such as induction of synaptic plasticity. However, due to limited structural information about the NMDAR transmembrane ion channel forming domain, the mechanism of divalent cation permeation and block is understood poorly. In this paper we developed an atomistic model of the transmembrane domain (TMD) of NMDARs composed of GluN1 and GluN2A subunits (GluN1/2A receptors). The model was generated using (a) a homology model based on the structure of the NaK channel and a partially resolved structure of an AMPA receptor (AMPAR), and (b) a partially resolved X-ray structure of GluN1/2B NMDARs. Refinement and extensive Molecular Dynamics (MD) simulations of the NMDAR TMD model were performed in explicit lipid bilayer membrane and water. Targeted MD with simulated annealing was introduced to promote structure refinement. Putative positions of the Mg2+ and Ca2+ ions in the ion channel divalent cation binding site are proposed. Differences in the structural and dynamic behavior of the channel protein in the presence of Mg2+ or Ca2+ are analyzed. NMDAR protein conformational flexibility was similar with no ion bound to the divalent cation binding site and with Ca2+ bound, whereas Mg2+ binding reduced protein fluctuations. While bound at the binding site both ions retained their preferred ligand coordination numbers: 6 for Mg2+, and 7-8 for Ca2+. Four asparagine side chain oxygens, a back-bone oxygen, and a water molecule participated in binding a Mg2+ ion. The Ca2+ ion first coordination shell ligands typically included four to five side-chain oxygen atoms of the binding site asparagine residues, two water molecules and zero to two backbone oxygens of the GluN2B subunits. These results demonstrate the importance of high-resolution channel structures for elucidation of mechanisms of NMDAR permeation and block.

Dataset Information

Interactions of the EGFR juxtamembrane domain with PIP2-containing lipid bilayers: Insights from multiscale molecular dynamics simulations.

Background

Methods

Results

Conclusions

General significance

Publications

Interactions of the EGFR juxtamembrane domain with PIP2-containing lipid bilayers: Insights from multiscale molecular dynamics simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets