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

Transmembrane domain dimerization induces cholesterol rafts in curved lipid bilayers

Are the dimerization of transmembrane (TM) domains and the reorganization of the lipid bilayer two independent events?

GxxxG motifs, phenylalanine, and cholesterol guide the self-association of transmembrane domains of ErbB2 receptors

GxxxG motifs are common in transmembrane domains of membrane proteins and are often introduced to artificial peptides to inhibit or promote association to stable structures. The transmembrane domain of ErbB2 presents two separate such motifs that are proposed to be connected to stability and activity of the dimer. Using molecular simulations, we show that these sequences play a critical role during the recognition stage, forming transient complexes that lead to stable dimers. In pure phospholipid bilayers association occurs by contacts formed at the C-terminus promoted by the presence of phenylalanine residues. Helices subsequently rotate to eventually pack at short separations favored by lipid entropic contributions. In contrast, at intermediate cholesterol concentrations, a different pathway is followed that involves dimers with a weaker interface toward the N-terminus. However, at high cholesterol content, a switch toward the C-terminus is observed with an overall nonmonotonic change of the dimerization affinity. This conformational switch modulated by cholesterol has important implications on the thermodynamic, structural, and kinetic characteristics of helix-helix association in lipid membranes.

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.

Systematic molecular dynamics searching in a lipid bilayer: Application to the glycophorin A and oncogenic ErbB-2 transmembrane domains

Molecular dynamics (MD) simulations of proteins in a lipid bilayer environment are usually undertaken with one or a few starting structures. Here we report a search protocol for systematically exploring the possible interactions in helical bundle transmembrane proteins, a frequently occurring structural motif. The search protocol correctly identifies the experimentally known structure of the dimeric human glycophorin A transmembrane domain as the lowest energy structure among five different models without any prior assumptions, whilst an identical in vacuo search fails to identify the correct structure. The lowest energy structure from the search in a lipid bilayer has a root mean square deviation of 1.1A to the experimental structure. We have applied the same search protocol to the unknown transmembrane structure of the oncogenic mutant ErbB-2 protein, a member of the family of epidermal growth factor receptors. Resulting structures show the role of glutamic acid hydrogen bonding and close helical packing. Water molecules may also play a key role in stabilisation of the transmembrane helix association.

Characterization of the archaeal thermophile Sulfolobus turreted icosahedral virus validates an evolutionary link among double-stranded DNA viruses from all domains of life

Icosahedral nontailed double-stranded DNA (dsDNA) viruses are present in all three domains of life, leading to speculation about a common viral ancestor that predates the divergence of Eukarya, Bacteria, and Archaea. This suggestion is supported by the shared general architecture of this group of viruses and the common fold of their major capsid protein. However, limited information on the diversity and replication of archaeal viruses, in general, has hampered further analysis. Sulfolobus turreted icosahedral virus (STIV), isolated from a hot spring in Yellowstone National Park, was the first icosahedral virus with an archaeal host to be described. Here we present a detailed characterization of the components forming this unusual virus. Using a proteomics-based approach, we identified nine viral and two host proteins from purified STIV particles. Interestingly, one of the viral proteins originates from a reading frame lacking a consensus start site. The major capsid protein (B345) was found to be glycosylated, implying a strong similarity to proteins from other dsDNA viruses. Sequence analysis and structural predication of virion-associated viral proteins suggest that they may have roles in DNA packaging, penton formation, and protein-protein interaction. The presence of an internal lipid layer containing acidic tetraether lipids has also been confirmed. The previously presented structural models in conjunction with the protein, lipid, and carbohydrate information reported here reveal that STIV is strikingly similar to viruses associated with the Bacteria and Eukarya domains of life, further strengthening the hypothesis for a common ancestor of this group of dsDNA viruses from all domains of life.

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.

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.

Molecular dynamics (MD) investigations of preformed structures of the transmembrane domain of the oncogenic Neu receptor dimer in a DMPC bilayer

The critical Val/Glu mutation in the membrane spanning domain of the rat Neu receptor confers the ability for ligand-independent signaling and leads to increased dimerization and transforming ability. There is evidence that the two transmembrane interacting helices play a role in receptor activation by imposing orientation constraints to the intracellular tyrosine kinase domains. By using MD simulations we have attempted to discriminate between correct and improper helix-helix packing by examining the structural and energetic properties of preformed left-handed and right-handed structures in a fully hydrated DMPC bilayer. The best energetic balance between the residues at the helix-helix interface and the residues exposed to the lipids is obtained for helices in symmetrical left-handed interactions packed together via Glu side chain/Ala backbone interhelical hydrogen bonds. Analyses demonstrate the importance of the ATVEG motif in helix-helix packing and point to additional contacting residues necessary for association. Our findings, all consistent with experimental data, suggest that a symmetrical left-handed structure of the helices could be the transmembrane domain configuration that promotes receptor activation and transformation. The present study may provide further insight into signal transduction mechanisms of the ErbB/Neu receptors.

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.