Molecular dynamics simulations of the ErbB-2 transmembrane domain within an explicit membrane environment: comparison with vacuum simulations

ORF8a of SARS-CoV forms an ion channel: experiments and molecular dynamics simulations

ORF8a protein is 39 residues long and contains a single transmembrane domain. The protein is synthesized using solid phase peptide synthesis and reconstituted into artificial lipid bilayers that forms cation-selective ion channels with a main conductance level of 8.9±0.8pS at elevated temperature (38.5°C). Computational modeling studies including multi nanosecond molecular dynamics simulations in a hydrated POPC lipid bilayer are done with a 22 amino acid transmembrane helix to predict a putative homooligomeric helical bundle model. A structural model of a pentameric bundle is proposed with cysteines, serines and threonines facing the pore.

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.

Transmembrane Helix Packing of ErbB/Neu Receptor in Membrane Environment: A Molecular Dynamics Study

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains. Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the Neu(TM) helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization-competent structure responsible of the proper topology necessary for receptor activation.

Transmembrane Peptides from Tyrosine Kinase Receptor. Mutation-related Behavior in a Lipid Bilayer Investigated by Molecular Dynamics Simulations

Polar mutations in transmembrane alpha helices may alter the structural details of the hydrophobic sequences and control intermolecular contacts. We have performed molecular dynamics simulations on the transmembrane domain of the proto-oncogenic and the oncogenic forms of the Neu receptor in a fluid DMPC bilayer to test whether the Glu mutation which replaces the Val residue at position 664 may alter the helical structure and its insertion in the membrane. The simulations show that the wild and the mutant forms of the transmembrane domain have a different behavior in the bilayer. The native transmembrane sequence is found to be more flexible than in the presence of the Glu mutation, characterized by a tendency to pi deformation to accommodate the helix length to the membrane thickness. The mutant form of this domain does not evidence helical deformation in the present simulation. Hydrophobic matching is achieved both by a larger helix tilt and a vertical shift of the helix towards the membrane interface, favoring the accessibility of the Glu side chain to the membrane environment. A rapid exchange of hydrogen bond interactions with the surrounding water molecules and the lipid headgroups is observed. The difference in the behavior between the two peptides in a membrane environment was also observed experimentally. Both simulation and experimental results agree with the hypothesis that water may act as an intermediate for the formation of cross links between the facing Glu side chains stabilizing the dimer.

A structural model of EmrE, a multi-drug transporter from Escherichia coli

Using a recently reported computational method, we describe an approach to model the structure of EmrE, a proton coupled multi-drug transporter of Escherichia coli. EmrE is the smallest ion-coupled transporter known; it functions as an oligomer and each monomer comprises four transmembrane segments. Because of its size, EmrE provides a unique experimental paradigm. The computational method does not afford a unique solution for the monomer. The experimental constraints available were used to select the most likely structure and to dock two monomers together to yield a dimer. The model is further validated by modeling of Hsmr, an EmrE homolog with a remarkable amino acid composition with over 40% of Ala and Val. The Hsmr model is similar to that of EmrE, with the majority of the Ala or Val residues facing the lipid. In addition, the model of EmrE features a putative substrate-binding site very similar to that observed in BmrR, a transcription activator of multi-drug transporters, with a similar substrate profile. The two crucial residues that couple proton fluxes with substrate binding in the homo-dimer of EmrE, Glu-14, have a spatial arrangement that agrees with proposed molecular mechanisms of transport.

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.

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 Glu664-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 alpha-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 pi-helical hydrogen bonding develops, introducing a kink in the backbone structure. However, there is no indication of the formation of a pi 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 Glu664. If the energy cost of such water penetration is significant enough, this may be a factor in the enhanced dimerization affinity of Glu664-mutant Neu.

Enhanced sampling by multiple molecular dynamics trajectories: carbonmonoxy myoglobin 10 μs A0 → A1–3 transition from ten 400 picosecond simulations

The utility of multiple trajectories to extend the time scale of molecular dynamics simulations is reported for the spectroscopic A-states of carbonmonoxy myoglobin (MbCO). Experimentally, the A0-->A(1-3) transition has been observed to be 10 micros at 300 K, which is beyond the time scale of standard molecular dynamics simulations. To simulate this transition, 10 short (400 ps) and two longer time (1.2 ns) molecular dynamics trajectories, starting from five different crystallographic and solution phase structures with random initial velocities centered in a 37 A radius sphere of water, have been used to sample the native-fold of MbCO. Analysis of the ensemble of structures gathered over the cumulative 5.6 ns reveals two biomolecular motions involving the side chains of His64 and Arg45 to explain the spectroscopic states of MbCO. The 10 micros A0-->A(1-3) transition involves the motion of His64, where distance between His64 and CO is found to vary up to 8.8 +/- 1.0 A during the transition of His64 from the ligand (A(1-3)) to bulk solvent (A0). The His64 motion occurs within a single trajectory only once, however the multiple trajectories populate the spectroscopic A-states fully. Consequently, multiple independent molecular dynamics simulations have been found to extend biomolecular motion from 5 ns of total simulation to experimental phenomena on the microsecond time scale.

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.

Conformational dynamics of helix S6 from Shaker potassium channel: simulation studies

Prolines in transmembrane (TM) alpha-helices are believed to play an important structural and/or functional role in membrane proteins. At a structural level a proline residue distorts alpha-helical structure due to the loss of at least one stabilizing backbone hydrogen bond, and introduces flexibility in the helix that may result in substantial kink and swivel motions about the effective "hinge." At a functional level, for example in Kv channels, it is believed that proline-induced molecular hinges may have a direct role in gating, i.e., the conformational change linked to opening/closing the channel to movement of ions. In this article we study the conformational dynamics of the S6 TM helix from of the Kv channel Shaker, which possesses the motif PVP--a motif that is conserved in Kv channels. We perform multiple molecular dynamics simulations of single S6 helices in a membrane-mimetic environment in order to effectively map the kink-swivel conformational space of the protein, exploiting the ability of multiple simulations to achieve greater sampling. We show that the presence of proline locally perturbs the helix, disrupting local dihedral angles and producing local twist and unwinding in the region of the hinge--an effect that is relaxed with distance from the PVP motif. We furthermore show that motions about the hinge are highly anisotropic, reflecting a preferred region of kink-swivel conformation space that may have implications for the gating process.

Mapping the energy surface of transmembrane helix-helix interactions

Transmembrane helices are no longer believed to be just hydrophobic segments that exist solely to anchor proteins to a lipid bilayer, but rather they appear to have the capacity to specify function and structure. Specific interactions take place between hydrophobic segments within the lipid bilayer whereby subtle mutations that normally would be considered innocuous can result in dramatic structural differences. That such specificity takes place within the lipid bilayer implies that it may be possible to identify the most favorable interaction surface of transmembrane alpha-helices based on computational methods alone, as shown in this study. Herein, an attempt is made to map the energy surface of several transmembrane helix-helix interactions for several homo-oligomerizing proteins, where experimental data regarding their structure exist (glycophorin A, phospholamban, Influenza virus A M2, Influenza virus C CM2, and HIV vpu). It is shown that due to symmetry constraints in homo-oligomers the computational problem can be simplified. The results obtained are mostly consistent with known structural data and may additionally provide a view of possible alternate and intermediate configurations.

Dimer interface of transmembrane domains for neu/erbB-2 receptor dimerization and transforming activation: a model revealed by molecular dynamics simulations

The specific point mutation Val-->Glu664 within the transmembrane domain of the neu/erbB-2 receptor is associated with increased receptor dimerization and increased receptor tyrosine kinase activity resulting in malignant transformation of cells. It is well established that Glu and residues in proximity are necessary for receptor dimerization but many studies suggest that other intramembrane constraints, not yet elucidated, are determinant for transformation. In this work, we investigated dimer models both to understand the structural role of the Glu mutation in the transmembrane domain association and to determine helix-helix contacts required for oncogenic transformation. Different types of helix-helix association based on data resulting from Cys mutational studies of the full wild receptor and spectroscopic data of transmembrane neu peptides have been explored by molecular dynamics simulations. The study leads to propose a model for the dimeric association of the transmembrane domains of the oncogenic neu receptor showing left-handed interactions of the two helices stabilized by symmetrical hydrogen bonding interactions involving the Glu side chain on one helix and the facing carbonyl of Ala661 on the second helix. Contacting residues observed in the symmetric interface explain the transforming activity or the non transforming activity of many neu mutants. Moreover the left-handed coiled coil structure is fully consistent with recent results proving the role of rotational linkage of the transmembrane domain with the kinase domain. Comparison between the predicted dimer model and those presumed from experiments strongly suggests helix flexibility in the extracellular juxtamembrane region.

Molecular dynamics and (2)H-NMR study of the influence of an amphiphilic peptide on membrane order and dynamics

A molecular dynamics simulation of a fully hydrated model membrane consisting of 12 molecules of 1, 2-dimyristoyl-sn-glycero-3-phosphocholine, one amphiphilic peptide with the sequence acetyl-Lys-Lys-Gly-Leu(16)-Lys-Lys-Ala-amide, and 593 water molecules was performed for 1.06 ns (Belohorcova, K., J. H. Davis, T. B. Woolf, and B. Roux. 1997. Biophys. J. 73:3039-3055). The analysis presented here is primarily focused on the phospholipid component and the results are compared with experimental (2)H-NMR studies of the lipid component of mixtures of the same peptide and lipid at a molar ratio of 1:32, and with earlier studies of closely related peptide/lipid mixtures. The phospholipid chain and headgroup isomer populations and isomerization rates compare favorably with previous simulations and experimental measurements. Of particular interest is the effect of the peptide on the phospholipid headgroup and hydrocarbon chain orientational order calculated from the simulation, which also agree well with experimental measurements performed on this and closely related systems. Comparison of the experimental results with the simulations not only shows that there is significant agreement between the two methods, but also provides new insight into the effect of the peptide on the lipid dynamics. In particular, these results confirm that a membrane spanning peptide has little effect on lipid chain order, and bilayer thickness if its hydrophobic length closely matches the lipid hydrocarbon thickness. In addition, we find that the peptide can have a strong ordering effect if it is longer than the lipid hydrophobic thickness.