Structure prediction of the dimeric neu/ErbB-2 transmembrane domain from multi-nanosecond molecular dynamics simulations

Dimerization of the neu/ErbB-2 receptor tyrosine kinase is a necessary but not a sufficient step for signaling. Despite the efforts expended to identify the molecular interactions responsible for receptor-receptor contacts and particularly those involving the transmembrane domain, structural details are still unknown. In this work, molecular dynamics simulations of the helical transmembrane domain (TM) of neu and ErbB-2 receptors are used to predict their dimer structure both in the wild and oncogenic forms. A global conformational search method, applied to define the best orientations of parallel helices, showed an energetically favorable configuration with the specific mutation site within the interface, common for both the nontransforming and the transforming neu/ErbB-2 TM dimers. Starting from this configuration, a total of 10 simulations, about 1.4 ns each, performed in vacuum, without any constraints, show that the two helices preferentially wrap in left-handed interactions with a packing angle at about 20°. The resulting structures are nonsymmetric and the hydrogen bond network analysis shows that helices experience π local distortions that facilitate inter-helix hydrogen bond interactions and may result in a change in the helix packing, leading to a symmetric interface. For the mutated sequences, we show that the Glu side chain interacts directly with its cognate or with carbonyl groups of the facing backbone. We show that the connectivity between interfacial residues conforms to the knobs-into-holes packing mode of transmembrane helices. The dimeric interface described in our models is discussed with respect to mutagenesis studies. Received: 12 March 1999 / Revised version: 23 August 1999 / Accepted: 23 August 1999     

Molecular dynamics simulations of the E1/E2 transmembrane domain of the Semliki Forest virus

Transmembrane (TM) helix-helix interactions are important for virus budding and fusion. We have developed a simulation strategy that reveals the main features of the helical packing between the TM domains of the two glycoproteins E1 and E2 of the alpha-virus Semliki Forest virus and that can be extrapolated to sketch TM helical packing in other alpha-viruses. Molecular dynamics simulations were performed in wild-type and mutant peptides, both isolated and forming E1/E2 complexes. The simulations revealed that the isolated wild-type E1 peptide formed a more flexible helix than the rest of peptides and that the wild-type E1/E2 complex consists of two helices that intimately pack their N-terminals. The residues located at the interhelical interface displayed the typical motif of the left-handed coiled-coils. These were small and medium residues as Gly, Ala, Ser, and Leu, which also had the possibility to form interhelical Calpha-H...O hydrogen bonds. Results from the mutant complexes suggested that correct packing is a compromise between these residues at both E1 and E2 interhelical interfaces. This compromise allowed prediction of E1-E2 contact residues in the TM spanning domain of other alphaviruses even though the sequence identity of E2 peptides is low in this domain.     

Molecular dynamics simulations of the transmembrane domain of the oncogenic ErbB2 receptor dimer in a DMPC bilayer

Molecular dynamics simulations of an atomic model of the transmembrane domain of the oncogenic ErbB2 receptor dimer embedded in an explicit dimyristoylphosphatidylcholine (DMPC) bilayer were performed for more than 4 ns. The oncogenic Glu mutation in the membrane spanning segment plays a major role in tyrosine kinase activity and receptor dimerization, and is thought to be partly responsible for the structure of the transmembrane domain of the active receptor. MD results show that the interactions between the two transmembrane helices are characteristic of a left-handed packing as previously demonstrated from in vacuo simulations. Moreover, MD results reveal the absence of persistent hydrogen bonds between the Glu side chains in a membrane environment, which raise the question of the ability for Glu alone to stabilize the TM domain of the ErbB2 receptor. Interestingly the formation of the alpha-pi motif in the two ErbB2 transmembrane helices confirms the concept of intrinsic sequence-induced conformational flexibility. From a careful analysis of our MD results, we suggest that the left-handed helix-helix packing could be the key to correctly orient the intracellular domain of the activated receptor dimer. The prediction of such interactions from computer simulations represents a new step towards the understanding of signaling mechanisms.     

Molecular dynamics simulations of the NGF-TrkA domain 5 complex and comparison with biological data

The nerve growth factor (NGF) is an important pharmacological target for Alzheimer's and other neurodegenerative diseases. Its action derives partly from its binding to the tyrosine kinase A receptor (TrkA). Here we study energetics and dynamics of the NGF-TrkA complex by carrying out multinanosecond molecular dynamics simulations, accompanied by electrostatic calculations based on the Poisson-Boltzmann equation. Our calculations, which are based on the x-ray structure of the complex, suggest that some of the mutations affecting dramatically the affinity of the complex involve residues that form highly favorable, direct or water-mediated hydrogen bond interactions at the ligand-receptor interface and, in some cases, that also critically participate to the large-scale motions of the complex. Furthermore, our calculations offer a rationale for the small effect on binding affinity observed upon specific mutations involving large changes in electrostatics (i.e., the charged-to-neutral mutations). Finally, these calculations, used along with the mutagenesis data, provide a basis for designing new peptides that mimic NGF in TrkA binding function.     

Molecular dynamics simulations of E. coli MsbA transmembrane domain: formation of a semipore structure

The human P-glycoprotein (MDR1/P-gp) is an ATP-binding cassette (ABC) transporter involved in cellular response to chemical stress and failures of anticancer chemotherapy. In the absence of a high-resolution structure for P-gp, we were interested in the closest P-gp homolog for which a crystal structure is available: the bacterial ABC transporter MsbA. Here we present the molecular dynamics simulations performed on the transmembrane domain of the open-state MsbA in a bilayer composed of palmitoyl oleoyl phosphatidylethanolamine lipids. The system studied contained more than 90,000 atoms and was simulated for 50 ns. This simulation shows that the open-state structure of MsbA can be stable in a membrane environment and provides invaluable insights into the structural relationships between the protein and its surrounding lipids. This study reveals the formation of a semipore-like structure stabilized by two key phospholipids which interact with the hinge region of the protein during the entire simulation. Multiple sequence alignments of ABC transporters reveal that one of the residues involved in the interaction with these two phospholipids are under a strong selection pressure specifically applied on the bacterial homologs of MsbA. Hence, comparison of molecular dynamics simulation and phylogenetic data appears as a powerful approach to investigate the functional relevance of molecular events occurring during simulations.     

Molecular dynamics simulations of the ligand-binding domain of an N-methyl-D-aspartate receptor

The mechanism of partial agonism at N-methyl-D-aspartate receptors is an unresolved issue, especially with respect to the role of protein dynamics. We have performed multiple molecular dynamics simulations (7 x 20 ns) to examine the behavior of the ligand-binding core of the NR1 subunit with a series of ligands. Our results show that water plays an important role in stabilizing different conformations of the core and how a closed cleft conformation of the protein might be stabilized in the absence of ligands. In the case of ligand-bound simulations with both full and partial agonists, we observed that ligands within the binding cleft may undergo distinct conformational changes, without grossly influencing the degree of cleft closure within the ligand-binding domain. In agreement with recently published crystallographic data, we also observe similar changes in backbone torsions corresponding to the hinge region between the two lobes for the partial agonist, D-cycloserine. This observation rationalizes the classification of D-cycloserine as a partial agonist and should provide a basis with which to predict partial agonism in this class of receptor by analyzing the behavior of these torsions with other potential ligands.     

Molecular dynamics simulation of the M2 helices within the nicotinic acetylcholine receptor transmembrane domain: structure and collective motions

Multiple nanosecond duration molecular dynamics simulations were performed on the transmembrane region of the Torpedo nicotinic acetylcholine receptor embedded within a bilayer mimetic octane slab. The M2 helices and M2-M3 loop regions were free to move, whereas the outer (M1, M3, M4) helix bundle was backbone restrained. The M2 helices largely retain their hydrogen-bonding pattern throughout the simulation, with some distortions in the helical end and loop regions. All of the M2 helices exhibit bending motions, with the hinge point in the vicinity of the central hydrophobic gate region (corresponding to residues alphaL251 and alphaV255). The bending motions of the M2 helices lead to a degree of dynamic narrowing of the pore in the region of the proposed hydrophobic gate. Calculations of Born energy profiles for various structures along the simulation trajectory suggest that the conformations of the M2 bundle sampled correspond to a closed conformation of the channel. Principal components analyses of each of the M2 helices, and of the five-helix M2 bundle, reveal concerted motions that may be relevant to channel function. Normal mode analyses using the anisotropic network model reveal collective motions similar to those identified by principal components analyses.     

Molecular dynamics simulations on the interaction of the transmembrane NavAb channel with cholesterol and lipids in the membrane

Increased cholesterol levels are associated with multiple pathological conditions. In this work, molecular dynamics simulations were applied to observe the influence of membrane cholesterol levels on a voltage-gated sodium channel. Different lipid compositions are modeled around the channel to obtain information about the possible effects by which cholesterol influences Na_vAb channels. Cholesterol was normally not directly interacting with either the closed or inactivated conformation. Cholesterol increased lipid packing implying that it plays a crucial role in restricting lipid movement in the region around 1 nm of the channel in a 1-palmitoyl-2-oeleoyl phosphatidylcholine matrix. Our results provide the first computational indication of an indirect modulation of Na_vAb channels by membrane cholesterol.     

Molecular dynamics simulations of membrane proteins

Membrane proteins control the traffic across cell membranes and thereby play an essential role in cell function from transport of various solutes to immune response via molecular recognition. Because it is very difficult to determine the structures of membrane proteins experimentally, computational methods have been increasingly used to study their structure and function. Here we focus on two classes of membrane proteins—ion channels and transporters—which are responsible for the generation of action potentials in nerves, muscles, and other excitable cells. We describe how computational methods have been used to construct models for these proteins and to study the transport mechanism. The main computational tool is the molecular dynamics (MD) simulation, which can be used for everything from refinement of protein structures to free energy calculations of transport processes. We illustrate with specific examples from gramicidin and potassium channels and aspartate transporters how the function of these membrane proteins can be investigated using MD simulations.     

The transmembrane domain of the oncogenic mutant ErbB-2 receptor: a structure obtained from site-specific infrared dichroism and molecular dynamics

ErbB-2 is a member of the family of epidermal growth factor receptors, which shows an oncogenic mutation in the rat gene neu, Val664Glu in the transmembrane domain that causes permanent dimerisation and subsequently leads to uncontrollable cell division and tumour formation. We have obtained the alpha-helical structure of the mutant transmembrane domain dimer experimentally with site-specific infrared dichroism (SSID) based on six transmembrane peptides with 13C18O carbonyl group-labelled residues. The derived orientational data indicate a local helix tilt ranging from 28(+/-6) degrees to 22(+/-4) degrees. Altogether using orientational constraints from SSID and experimental alpha-helical constraints while performing a systematic conformational search including molecular dynamics simulation in a lipid bilayer, we have obtained a unique experimentally defined atomic structure. The resulting structure consists of a right handed alpha-helical bundle with the residues Ile659, Val663, Leu667, Ile671, Val674 and Leu679 in the dimerisation interface. The right-handed bundle is in contrast to the left-handed structures obtained in previous modelling efforts. In order to facilitate tight helical packing, the spacious Glu664 residues do not interact directly but with water molecules that enter the bilayer.     

Molecular dynamics simulations of the TrkH membrane protein

TrkH is a transmembrane protein that mediates uptake of K(+) through the cell membrane. Despite the recent determination of its crystallographic structure, the nature of the permeation mechanism is still unknown, that is, whether K(+) ions move across TrkH by active transport or passive diffusion. Here, molecular dynamics simulations and the umbrella sampling technique have been employed to shed light on this question. The existence of binding site S3 and two alternative binding sites have been characterized. Analysis of the coordination number renders values that are almost constant, with a full contribution from the carbonyls of the protein only at S3. This observation contrasts with observations of K(+) channels, where the contribution of the protein to the coordination number is roughly constant in all four binding sites. An intramembrane loop is found immediately after the selectivity filter at the intracellular side of the protein, which obstructs the permeation pathway, and this is reflected in the magnitude of the energy barriers.     

Molecular dynamics simulations of isolated transmembrane helices of potassium channels

In the middle of the S6 helix in voltage-gated potassium channels there is a highly conserved Pro-Val-Pro motif, while the equivalent M2 helix of inward rectifier potassium channels contains a conserved glycine residue in a comparable position. The structural implications of these conserved motifs are of interest given the evidence that S6 and M2 are components of the lining of their respective pores. Multiple sequence alignment and TM helix prediction methods were used to define consensus regions for S6 and M2. Ensembles of 50 structures for each helix were generated by simulated annealing and restrained molecular dynamics. Time-dependent fluctuations of S6 and M2 were investigated by long time scale molecular dynamics simulations on representative members of each ensemble carried out in vacuo in the presence and absence of a hydrophobic potential that mimics a lipid bilayer. The results are discussed in terms of the structural basis of the kink in S6 and M2 and of a putative functional role for flexible helices as “molecular swivels.” © 1996 John Wiley & Sons, Inc.     

Molecular dynamics simulations of ethanol binding to the transmembrane domain of the glycine receptor: implications for the channel potentiation mechanism

The glycine receptor (GlyR) is potentiated by ethanol and other anesthetics. The potentiation mechanism at the molecular level is unknown and remains elusive, but mutagenic studies have shown that ethanol and other volatile anesthetics bind to a pocket between TM1, TM2, and TM3. The present study extends previous studies (Cheng et al., Proteins 2007;68:581-593) wherein we conducted homology modeling and molecular dynamics (MD) simulations to construct models of the homopentameric alpha1 subunits of the GlyR transmembrane domain in open and closed states. To understand the potentiation of GlyR by ethanol we compare the binding of ethanol molecules to the channel in these different states. We observe that ethanol stably resides inside solvent-accessible cavities found in the open state of GlyR that are formed by I229 (of TM1) in one subunit and S267 and A288 (of TM2 and TM3, respectively) in the adjacent subunit. The volume of these putative binding pockets is state-dependent. Selective binding to the open states of receptors has been proposed to explain the potentiating actions of this class of anesthetics. In accordance with this model, our MD simulations suggest that the potentiation of ethanol on GlyR may be effected through preferential binding of ethanol molecules to an inter-subunit binding pocket in the open state.     

The transmembrane domain of Neu in a lipid bilayer: molecular dynamics simulations

The results of full-atom molecular dynamics simulations of the transmembrane domains (TMDs) of both native, and Glu⁶⁶⁴-mutant (either protonated or unprotonated) Neu in an explicit fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer are presented. For the native TMD peptide, a 10.05 ns trajectory was collected, while for the mutant TMD peptides 5.05 ns trajectories were collected for each. The peptides in all three simulations display stable predominantly -helical hydrogen bonding throughout the trajectories. The only significant exception occurs near the C-terminal end of the native and unprotonated mutant TMDs just outside the level of the lipid headgroups, where -helical hydrogen bonding develops, introducing a kink in the backbone structure. However, there is no indication of the formation of a bulge within the hydrophobic region of either native or mutant peptides. Over the course of the simulation of the mutant peptide, it is found that a significant number of water molecules penetrate the hydrophobic region of the surrounding lipid molecules, effectively hydrating Glu⁶⁶⁴. If the energy cost of such water penetration is significant enough, this may be a factor in the enhanced dimerization affinity of Glu⁶⁶⁴-mutant Neu.     

Molecular dynamics simulations on the free and complexed N-terminal SH2 domain of SHP-2

Signal transduction events are often mediated by small protein domains such as SH2 (Src homology 2) domains that recognize phosphotyrosines (pY) and flanking sequences. In case of the SHP-2 receptor tyrosine phosphatase an N-terminal SH2 domain binds and inactivates the phosphatase (PTP) domain. The pY-peptide-binding site on the N-terminal SH2 domain does not overlap with the PTP binding region. Nevertheless, pY-peptide binding causes domain dissociation and phosphatase activation. Comparative multi-nanosecond molecular dynamics simulations on the N-SH2 domain in ligand-bound and free states have been performed to study the allosteric mechanism that leads to domain dissociation upon pY-peptide binding. Significant ligand-dependent differences in the conformational flexibility of regions that are involved in SH2-PTP domain association have been observed. The results support a mechanism of signal transduction where SH2-peptide binding modulates the domain flexibility and reduces its capacity to fit into the entrance of the PTP catalytic domain of SHP-2.     

Dynamics of the Hck-SH3 domain: comparison of experiment with multiple molecular dynamics simulations

Molecular dynamics calculations provide a method by which the dynamic properties of molecules can be explored over timescales and at a level of detail that cannot be obtained experimentally from NMR or X-ray analyses. Recent work (Philippopoulos M, Mandel AM, Palmer AG III, Lim C, 1997, Proteins 28:481-493) has indicated that the accuracy of these simulations is high, as measured by the correspondence of parameters extracted from these calculations to those determined through experimental means. Here, we investigate the dynamic behavior of the Src homology 3 (SH3) domain of hematopoietic cell kinase (Hck) via 5N backbone relaxation NMR studies and a set of four independent 4 ns solvated molecular dynamics calculations. We also find that molecular dynamics simulations accurately reproduce fast motion dynamics as estimated from generalized order parameter (S2) analysis for regions of the protein that have experimentally well-defined coordinates (i.e., stable secondary structural elements). However, for regions where the coordinates are not well defined, as indicated by high local root-mean-square deviations among NMR-determined structural family members or high B-factors/low electron density in X-ray crystallography determined structures, the parameters calculated from a short to moderate length (less than 5-10 ns) molecular dynamics trajectory are dependent on the particular coordinates chosen as a starting point for the simulation.     

Molecular dynamics simulations of the ligand-binding domain of the ionotropic glutamate receptor GluR2.

Ionotropic glutamate receptors are essential for fast synaptic nerve transmission. Recent x-ray structures for the ligand-binding (S1S2) region of the GluR2 alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA)-sensitive receptor have suggested how differences in protein/ligand interactions may determine whether a ligand will behave as a full agonist. We have used multiple molecular dynamics simulations of 2-5 ns duration to explore the structural dynamics of GluR2 S1S2 in the presence and absence of glutamate and in a complex with kainate. Our studies indicate that not only is the degree of domain closure dependent upon interactions with the ligand, but also that protein/ligand interactions influence the motion of the S2 domain with respect to S1. Differences in domain mobility between the three states (apo-S1S2, glutamate-bound, and kainate-bound) are surprisingly clear-cut. We discuss how these changes in dynamics may provide an explanation relating the mechanism of transmission of the agonist-binding event to channel opening. We also show here how the glutamate may adopt an alternative mode of binding not seen in the x-ray structure, which involves a key threonine (T480) side chain flipping into a new conformation. This new conformation results in an altered pattern of hydrogen bonding at the agonist-binding site.     

Molecular dynamics simulations of LysRS: an asymmetric state

We report molecular dynamics simulations of the Escherichia coli Lysyl-tRNA synthetase LysU isoform carried out as a benchmark for mutant simulations in in silico protein engineering efforts. Unlike previous studies of aminoacyl-tRNA synthetases, LysU is modelled in its full dimeric form with explicit solvent. While developing a suitable simulation protocol, we observed an asymmetry that persists despite improvements to the model. This prediction has directly led to experiments that establish a functional asymmetry in nucleotide binding by LysU. The development of a simulation protocol and validation of the model are presented here. The observed asymmetry is described and the role of protein flexibility in developing the asymmetry is discussed.     

Stability and membrane interactions of an autotransport protein: MD simulations of the Hia translocator domain in a complex membrane environment

Hia is a trimeric autotransporter found in the outer membrane of Haemphilus influenzae. The X-ray structure of Hia translocator domain revealed each monomer to consist of an α-helix connected via a loop to a 4-stranded β-sheet, thus the topology of the trimeric translocator domain is a 12-stranded β-barrel containing 3 α-helices that protrude from the mouth of the β-barrel into the extracellular medium. Molecular dynamics simulations of the Hia monomer and trimer have been employed to explore the interactions between the helices, β-barrel and connecting loops that may contribute to the stability of the trimer. In simulations of the Hia monomer we show that the central α-helix may stabilise the fold of the 4-stranded β-sheet. In simulations of the Hia trimer, a H-bond network involving residues in the β-barrel, α-helices and loops has been identified as providing stability for the trimeric arrangement of the monomers. Glutamine residues located in the loops connecting the α-helices to the β-barrel are orientated in a triangular arrangement such that each forms 2 hydrogen bonds to each of the corresponding glutamines in the other loops. In the absence of the loops, the β‐barrel becomes distorted. Simulations show that while the trimeric translocator domain β-barrel is inherently flexible, it is unlikely to accommodate the passenger domain in a folded conformation. Simulations of Hia in an asymmetric model of the outer membrane have revealed membrane–protein interactions that anchor the protein within its native membrane environment.