A new method to model membrane protein structure based on silent amino acid substitutions.

The importance of accurately modeling membrane proteins cannot be overstated, in lieu of the difficulties in solving their structures experimentally. Often, however, modeling procedures (e.g., global searching molecular dynamics) generate several possible candidates rather then pointing to a single model. Herein we present a new approach to select among candidate models based on the general hypothesis that silent amino acid substitutions, present in variants identified from evolutionary conservation data or mutagenesis analysis, do not affect the stability of a native structure but may destabilize the non-native structures also found. The proof of this hypothesis has been tested on the alpha-helical transmembrane domains of two homodimers, human glycophorin A and human CD3-zeta, a component of the T-cell receptor. For both proteins, only one structure was identified using all the variants. For glycophorin A, this structure is virtually identical to the structure determined experimentally by NMR. We present a model for the transmembrane domain of CD3-zeta that is consistent with predictions based on mutagenesis, homology modeling, and the presence of a disulfide bond. Our experiments suggest that this method allows the prediction of transmembrane domain structure based only on widely available evolutionary conservation data.     

Transmembrane four-helix bundle of influenza A M2 protein channel: structural implications from helix tilt and orientation.

The transmembrane portion of the M2 protein from the Influenza A virus has been studied in hydrated dimyristroylphosphotidylcholine lipid bilayers with solid-state NMR. Orientational constraints were obtained from isotopically labeled peptide samples mechanically aligned between thin glass plates. 15N chemical shifts from single site labeled samples constrain the molecular frame with respect to the magnetic field. When these constraints are applied to the peptide, modeled as a uniform alpha-helix, the tilt of the helix with respect to the bilayer normal was determined to be 33 degrees +/- 3 degrees. Furthermore, the orientation about the helix axis was also determined within an error of +/- 30 degrees. These results imply that the packing of this tetrameric protein is in a left-handed four-helix bundle. Only with such a large tilt angle are the hydrophilic residues aligned to the channel axis.     

Hydrophobic matching controls the tilt and stability of the dimeric platelet-derived growth factor receptor (PDGFR) β transmembrane segment

The platelet-derived growth factor receptor β is a member of the cell surface receptor tyrosine kinase family and dimerizes upon activation. We determined the structure of the transmembrane segment in dodecylphosphocholine micelles by liquid-state NMR and found that it forms a stable left-handed helical dimer. Solid-state NMR and oriented circular dichroism were used to measure the tilt angle of the helical segments in macroscopically aligned model membranes with different acyl chain lengths. Both methods showed that decreasing bilayer thickness (DEPC-POPC-DMPC) led to an increase in the helix tilt angle from 10° to 30° with respect to the bilayer normal. At the same time, reconstitution of the comparatively long hydrophobic segment became less effective, eventually resulting in complete protein aggregation in the short-chain lipid DLPC. Unrestrained molecular dynamics simulations of the dimer were carried out in explicit lipid bilayers (DEPC, POPC, DMPC, sphingomyelin), confirming the observed dependence of the helix tilt angle on bilayer thickness. Notably, molecular dynamics revealed that the left-handed dimer gets tilted en bloc, whereas conformational transitions to alternative (e.g. right-handed dimeric) states were not supported. The experimental data along with the simulation results demonstrate a pronounced interplay between the platelet-directed growth factor receptor β transmembrane segment and the bilayer thickness. The effect of hydrophobic mismatch might play a key role in the redistribution and activation of the receptor within different lipid microdomains of the plasma membrane in vivo.     

A conserved tetrameric interaction of cry toxin helix α3 suggests a functional role for toxin oligomerization

Crystal (Cry) toxins are widely used for insect control, but their mechanism of toxicity is still uncertain. These toxins can form lytic pores in vitro, and water soluble tetrameric pre-pore intermediates have been reported. Even the precise oligomeric state of the toxin in membranes, trimeric or tetrameric, is still a debated issue. Based on previous reports, we have assumed that interactions between toxin monomers in solution are at least partly mediated by domain I, and we have analyzed in silico the homo-oligomerization tendencies of the domain I α-helices individually. Using many homologous sequences for each α-helix, our strategy allows selection of evolutionarily conserved interactions. These interactions appeared only in helices α3 and α5, but only α3 produced a suitably oriented or α-helical sample in lipid bilayers, forming homotetramers in C14-betaine, and allowing determination of its rotational orientation in lipid bilayers using site-specific infrared dichroism (SSID). The determined orientation in the tetrameric model is in agreement with only one of the evolutionarily conserved models. In addition mutation R99E, which was found to inhibit oligomerization experimentally, greatly destabilized the tetramer in molecular dynamic simulations. In this model, helix 3 is able to form inter-monomer interactions without significant rearrangements of domain I, which is compatible with the available crystal structure of Cry toxins in solution. The model presented here at least partially explains the reported tetrameric oligomerization of Cry toxins in solution and the inhibition of this oligomerization by a synthetic α3 peptide.     

Interpretation of H-NMR Experiments on the Orientation of the Transmembrane Helix WALP23 by Computer Simulations

Orientation, dynamics, and packing of transmembrane helical peptides are important determinants of membrane protein structure, dynamics, and function. Because it is difficult to investigate these aspects by studying real membrane proteins, model transmembrane helical peptides are widely used. NMR experiments provide information on both orientation and dynamics of peptides, but they require that motional models be interpreted. Different motional models yield different interpretations of quadrupolar splittings (QS) in terms of helix orientation and dynamics. Here, we use coarse-grained (CG) molecular dynamics (MD) simulations to investigate the behavior of a well-known model transmembrane peptide, WALP23, under different hydrophobic matching/mismatching conditions. We compare experimental ²H-NMR QS (directly measured in experiments), as well as helix tilt angle and azimuthal rotation (not directly measured), with CG MD simulation results. For QS, the agreement is significantly better than previously obtained with atomistic simulations, indicating that equilibrium sampling is more important than atomistic details for reproducing experimental QS. Calculations of helix orientation confirm that the interpretation of QS depends on the motional model used. Our simulations suggest that WALP23 can form dimers, which are more stable in an antiparallel arrangement. The origin of the preference for the antiparallel orientation lies not only in electrostatic interactions but also in better surface complementarity. In most cases, a mixture of monomers and antiparallel dimers provides better agreement with NMR data compared to the monomer and the parallel dimer. CG MD simulations allow predictions of helix orientation and dynamics and interpretation of QS data without requiring any assumption about the motional model.     

Structure of the influenza C virus CM2 protein transmembrane domain obtained by site-specific infrared dichroism and global molecular dynamics searching

The 115-residue protein CM2 from Influenza C virus has been recently characterized as a tetrameric integral membrane glycoprotein. Infrared spectroscopy and site-directed infrared dichroism were utilized here to determine its transmembrane structure. The transmembrane domain of CM2 is alpha-helical, and the helices are tilted by beta = (14.6 +/- 3.0) degrees from the membrane normal. The rotational pitch angle about the helix axis omega for the 1-(13)C-labeled residues Gly(59) and Leu(66) is omega = (218 +/- 17) degrees, where omega is defined as zero for a residue pointing in the direction of the helix tilt. A detailed structure was obtained from a global molecular dynamics search utilizing the orientational data as an energy refinement term. The structure consists of a left-handed coiled-coil with a helix crossing angle of Omega = 16 degrees. The putative transmembrane pore is occluded by the residue Met(65). In addition hydrogen/deuterium exchange experiments show that the core is not accessible to water.     

Molecular dynamics simulation of transmembrane polypeptide orientational fluctuations

The orientation and motion of a model lysine-terminated transmembrane polypeptide were investigated by molecular dynamics simulation. Recent 2H NMR studies of synthetic polypeptides with deuterated alanine side chains suggest that such transmembrane polypeptides undergo fast, axially symmetric reorientation about the bilayer normal but have a preferred average azimuthal orientation about the helix axis. In this work, interactions that might contribute to this behavior were investigated in a simulated system consisting of 64 molecules of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and one alpha-helical polypeptide with the sequence acetyl-KK-(LA)11-KK-amide. In one simulation, initiated with the peptide oriented along the bilayer normal, the system was allowed to evolve for 8.5 ns at 1 atm of pressure and a temperature of 55 degrees C. A second simulation was initiated with the peptide orientation chosen to match a set of experimentally observed alanine methyl deuteron quadrupole splittings and allowed to proceed for 10 ns. Simulated alanine methyl group orientations were found to be inequivalent, a result that is consistent with 2H NMR observations of specifically labeled polypeptides in POPC bilayers. Helix tilt varied substantially over the durations of both simulations. In the first simulation, the peptide tended toward an orientation about the helix axis similar to that suggested by experiment. In the second simulation, orientation about the helix axis tended to return to this value after an excursion. These results provide some insight into how interactions at the bilayer surface can constrain reorientation about the helix axis while accommodating large changes in helix tilt.     

Orientation and Dynamics of Peptides in Membranes Calculated from H-NMR Data

Solid-state ²H-NMR is routinely used to determine the alignment of membrane-bound peptides. Here we demonstrate that it can also provide a quantitative measure of the fluctuations around the distinct molecular axes. Using several dynamic models with increasing complexity, we reanalyzed published ²H-NMR data on two representative α-helical peptides: 1), the amphiphilic antimicrobial peptide PGLa, which permeabilizes membranes by going from a monomeric surface-bound to a dimeric tilted state and finally inserting as an oligomeric pore; and 2), the hydrophobic WALP23, which is a typical transmembrane segment, although previous analysis had yielded helix tilt angles much smaller than expected from hydrophobic mismatch and molecular dynamics simulations. Their ²H-NMR data were deconvoluted in terms of the two main helix orientation angles (representing the time-averaged peptide tilt and azimuthal rotation), as well as the amplitudes of fluctuation about the corresponding molecular axes (providing the dynamic picture). The mobility of PGLa is found to be moderate and to correlate well with the respective oligomeric states. WALP23 fluctuates more vigorously, now in better agreement with the molecular dynamics simulations and mismatch predictions. The analysis demonstrates that when ²H-NMR data are fitted to extract peptide orientation angles, an explicit representation of the peptide rigid-body angular fluctuations should be included.     

A structure for the trimeric MHC class II-associated invariant chain transmembrane domain

The major histocompatibility complex (MHC)-associated invariant chain (Ii) contains a single transmembrane domain that forms trimers. Ii is involved in the assembly of the MHC and antigen presentation, and is thus central to the function of the immune system. Here, we show by attenuated total reflectance, Fourier transform infrared (ATR-FTIR) spectroscopy that the transmembrane domain is alpha-helical and we provide a structural model of the transmembrane domain obtained by a combination of site-specific infrared dichroism and molecular dynamics (MD) simulations. This work resolves the backbone structure of a transmembrane peptide by multiple (13)C=(18)O labelling at ten different residues. A second purely computational approach, based on MD simulations of Ii transmembrane homologous sequences, yields a similar structure that is consistent with our experimental results. The structure presented forms a left-handed coiled coil with an average helix tilt of 13(+/-6) degrees; the residue Gln47 implicated in trimer formation forms strong interhelical contacts, Thr50 points to the inside of the trimeric coil and forms a network of hydrogen bonds.     

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of α-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 × 5 μs) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2× rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of ∼16°. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.     

Transmembrane helices that form two opposite homodimeric interactions: an asparagine scan study of alphaM and beta2 integrins

Integrins are alpha/beta heterodimers, but recent in vitro and in vivo experiments also suggest an ability to associate through their transmembrane domains to form homomeric interactions. While the results of some in vitro experiments are consistent with an interaction mediated by a GxxxG-like motif, homo-oligomers observed after in vivo cross-linking are consistent with an almost opposite helix-helix interface. We have shown recently that both models of interaction are compatible with evolutionary conservation data, and we predicted that the alpha-helices in both models would have a similar rotational orientation. Herein, we have tested our prediction using in vitro asparagine scan of five consecutive residues along the GxxxG-like motif of the transmembrane domain of alpha and beta integrins, alphaM and beta2. We show that Asn-mediated dimerization occurs twice for every turn of the helix, consistent with two almost opposite forms of interaction as suggested previously for alphaIIb and beta3 transmembrane domains. The orientational parameters helix tilt and rotational orientation of each of these two Asn-stabilized dimers were measured by site-specific infrared dichroism (SSID) in model lipid bilayers and were found to be consistent with our predicted computational models. Our results highlight an intrinsic tendency for integrin transmembrane alpha-helices to form two opposite types of homomeric interaction in addition to their heteromeric interactions and suggest that integrins may form complex and specific networks at the transmembrane domain during function.     

Competing interactions contributing to alpha-helical stability in aqueous solution.

The stability of a 15-residue peptide has been investigated using CD spectroscopy and molecular simulation techniques. The sequence of the peptide was designed to include key features that are known to stabilize alpha-helices, including ion pairs, helix dipole capping, peptide bond capping, and aromatic interactions. The degree of helicity has been determined experimentally by CD in three solvents (aqueous buffer, methanol, and trifluoroethanol) and at two temperatures. Simulations of the peptide in the aqueous system have been performed over 500 ps at the same two temperatures using a fully explicit solvent model. Consistent with the CD data, the degree of helicity is decreased at the higher temperature. Our analysis of the simulation results has focused on competition between different side-chain/side-chain and side-chain/main-chain interactions, which can, in principle, stabilize the helix. The unfolding in aqueous solution occurs at the amino terminus because the side-chain interactions are insufficient to stabilize both the helix dipole and the peptide hydrogen bonds. Loss of capping of the peptide backbone leads to water insertion within the first peptide hydrogen bond and hence unfolding. In contrast, the carboxy terminus of the alpha-helix is stable in both simulations because the C-terminal lysine residue stabilizes the helix dipole, but at the expense of an ion pair.     

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.     

Recursive use of evolutionary conservation data in molecular modeling of membrane proteins A model of the multidrug H+ antiporter emrE.

Membrane proteins are currently the most biomedically important family of proteins, serving as targets for the majority of pharmaceutical agents. It is also clear that they are invariably abundant in all of the genomes sequence so far, representing up to a third of all open reading frames. Finally, and regrettably, it is clear that they are highly resistant to structural elucidation, representing less than 0.2% of the Protein Data Bank. Recent accomplishments in genome sequencing efforts, however, may help offset this imbalance through the availability of evolutionary conservation data. Herein, we develop a novel approach, utilizing a combination of evolutionary conservation data and global searching molecular dynamics simulations to model membrane proteins, deriving a model for the multidrug H+ antiporter EmrE, a transmembrane four-helix bundle. Structures resulting from an extensive, rotational molecular dynamics search, were evaluated by comparing the residue specific interaction energy and the evolutionary conservation data. Subsequent rounds of molecular dynamics, in which confinement of the search space was undertaken in order to achieve a self consistent result, point to a structure that best satisfies the evolutionary conservation data. As the conservation patterns calculated for each of the helices suggested that the different conservation pattern for helix 3 (as well as being the most conserved) might be due to the oligomeric nature of EmrE, a dodecamer of helices was constructed based on the result of a search of helix 3 as a trimer. The resulting interaction energy per residue in the final model is in reasonable agreement with the evolutionary data and consistent with recent site directed mutagenesis experiments, pointing to the strength of this method as a general tool.     

A test for the structure of the left helix of poly-L-proline type II in peptides

The spectral criterion of a left-handed helix of the poly-L-proline II type was elaborated during the study of a number of synthesized oligopeptides (in a solid state and solution): (Gly-Pro-Pro)1-8, (Gly-Pro)1, (Gly-Pro-Ala)1-4, (Gly-Pro-Gly)1-4, (Gly-Pro-Pro) X (Gly-Pro-Gly)1-2(Gly-Pro-Pro), (Gly-Pro-Pro)n, (Orn3-Gly)n and also rat skin collagen by X-ray diffraction, circular dichroism and infrared spectroscopy methods; the characteristic shape of the left-handed helix CD spectrum was found. The change of spectral characteristics with the change of left-handed helix distortion was established. The linear noncooperative melting process of the left-handed conformation was demonstrated. The data obtained allow to determine qualitatively the presence of the left-handed helix in different polypeptides and proteins.     

Molecular modeling studies on CNG channel from bovine retinal rod: a structural model of the cyclic nucleotide-binding domain

A dimeric model of the cyclic nucleotide-binding domain of the all-alpha homomeric cyclic nucleotide-gated channel from bovine retinal rod is constructed. The model, based on the structure of the fairly homologous catabolite gene activator protein (Weber and Steitz, J Mol Biol 1987;198:311-326), is obtained by use of comparative modeling and molecular dynamics simulations. Our model provides a structural basis for the experimentally measured difference in activity between cAMP and cGMP, as well as the different solvent accessibilities of GLY597 in the complex with cGMP, with cAMP and in the protein in free state. In addition, it provides support for the rearrangement of the domain C helix on ligand binding and releasing proposed by Matulef et al. (Neuron 1999;24:443-452).     

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.     

Coarse-grained molecular dynamics of tetrameric transmembrane peptide bundles within a lipid bilayer

The conformations of model transmembrane peptides are studied to understand the structural and dynamical aspects of tetrameric bundles using a series of coarse grain (CG) molecular dynamics (MD) simulations since membrane proteins play a crucial role in cell function. In this work, two different amphipathic models have been constructed using similar hydrophobic/hydrophilic characteristics with two structurally distinct morphologies to evaluate the effect of roughness and hydrophilic topology on the structure of tetrameric bundles, one class that forms an ion-channel and one class that does not. Free energy calculations of typical amphipathic peptide topologies show that using a relatively smooth surface morphology allows for a stable conformation of the tetramer bundle in a diamond formation. However, the model with side chains attached to the core in order to roughen the surface has a stable square tetramer bundle which is consistent with experimental data and all-atom (AA) MD simulations. Comparisons of the CG simulations with AA MD simulations are in reasonable agreement with the formation of tetrameric homo-oligomers, partitioning within the lipid bilayer and tilt angle with respect to the bilayer normal. We concluded that a square or diamond shape tetrameric homo-oligomers could be stabilized by rational design of the peptide morphology and topology of the surface, thus allowing us to tune the permeability of the bundle or channel.     

Phospholemman transmembrane structure reveals potential interactions with Na+/K+-ATPase

Phospholemman (PLM) is a 72-residue bitopic cardiac transmembrane protein, which acts as a modulator of the Na(+)/K(+)-ATPase and the Na(+)/Ca(2+) exchanger and possibly forms taurine channels in nonheart tissue. This work presents a high resolution structural model obtained from a combination of site-specific infrared spectroscopy and experimentally constrained high throughput molecular dynamics (MD) simulations. Altogether, 37 experimental constraints, including nine long range orientational constraints, have been used during MD simulations in an explicit lipid bilayer/water system. The resulting tetrameric alpha-helical bundle has an average helix tilt of 7.3 degrees and a crossing angle close to 0 degrees . It does not reveal a hydrophilic pore, but instead strong interactions between various residues occlude any pore. The helix-helix packing is unusual, with Gly(19) and Gly(20) pointing to the outside of the helical bundle, facilitating potential interaction with other transmembrane proteins, thus providing a structural basis for the modulatory effect of PLM on the Na(+)/K(+)-ATPase. A two-stage model of interaction between PLM and the Na(+)/K(+)-ATPase is discussed involving PLM-ATPase interaction and subsequent formation of an unstable PLM trimer, which readily interacts with surrounding ATPase molecules. Further unconstrained MD simulations identified other packing models of PLM, one of which could potentially undergo a conformational transition to an open pore.