Interaction of a two-transmembrane-helix peptide with lipid bilayers and dodecyl sulfate micelles

To probe structural changes that occur when a membrane protein is transferred from lipid bilayers to SDS micelles, a fragment of bacteriorhodopsin containing transmembrane helical segments A and B was studied by fluorescence spectroscopy, molecular dynamics (MD) simulation, and stopped flow kinetics. In lipid bilayers, F?rster resonance energy transfer (FRET) was observed between tyrosine 57 on helix B and tryptophans 10 and 12 on helix A. FRET efficiency decreased substantially when the peptide was transferred to SDS. MD simulation showed no evidence for significant disruption of helix-helix interactions in SDS micelles. However, a cluster of water molecules was observed to form a hydrogen-bonded network with the phenolic hydroxyl group of tyrosine 57, which probably causes the disappearance of tyrosine-to-tryptophan FRET in SDS. The tryptophan quantum yield decreased in SDS, and the change occurred at nearly the same rate as membrane solubilization. The results provide a clear example of the importance of corroborating distance changes inferred from FRET by using complementary methods.     

Molecular Dynamics Simulations of Micelle Formation around Dimeric Glycophorin A Transmembrane Helices

Insertion and formation of membrane proteins involves the interaction of protein helices with one another in lipid environments. Researchers have studied glycophorin A (GpA) transmembrane helices embedded in sodium dodecyl sulfate (SDS) micelles to identify contacts significant for helix dimerization. However, a detailed picture of the conformation and dynamics of the GpA-SDS system cannot be obtained solely through experiment. Molecular dynamics simulations of SDS and a GpA dimer can provide an atomic-level picture of SDS aggregation and helix association. We report 2.5-ns simulations of GpA wild-type and mutants in a preformed micelle as well as a 32-ns simulation showing the formation of a complete micelle around wild-type GpA from an initially random placement of SDS molecules in an aqueous environment. In the latter case, an initial instability of GpA helices in water is reversed after the helices become surrounded by SDS. The properties of the spontaneously formed micelle surrounding the GpA are indistinguishable from those of the preformed micelle surrounding the GpA dimer.     

Defining the transmembrane helix of M2 protein from influenza A by molecular dynamics simulations in a lipid bilayer

Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules.     

Tertiary structure of bacteriorhodopsin. Positions and orientations of helices A and B in the structural map determined by neutron diffraction

Positions and rotations of two helices in the tertiary structure of bacteriorhodopsin have been studied by neutron diffraction using reconstituted, hybrid purple membrane samples. Purple membrane was biosynthetically 2H-labeled at non-exchangeable hydrogen positions of leucine and tryptophan residues. Two chymotryptic fragments were purified, encompassing either the first two or the last five of the seven putative transmembrane segments identified in the amino acid sequence of bacteriorhodopsin. The 2H-labeled fragments, diluted to variable extents with the identical, unlabeled fragment, were mixed with their unlabeled counterpart; bacteriorhodopsin was then renatured and reconstituted. The crystalline purple membrane samples thus obtained contained hybrid bacteriorhodopsin molecules in which certain transmembrane segments had been selectively 2H-labeled to various degrees. Neutron diffraction powder patterns were recorded and analyzed both by calculating difference Fourier maps and by model building. The two analyses yielded consistent results. The first and second transmembrane segments in the sequence correspond to helices 1 and 7 of the three-dimensional structure, respectively. Rotational orientations of these two helices were identified using best fits to the observed diffraction intensities. The data also put restrictions on the position of the third transmembrane segment. These observations are discussed in the context of folding models for bacteriorhodopsin, the environment of the retinal Schiff base, and site-directed mutagenesis experiments.     

Secondary structures comparison of aquaporin-1 and bacteriorhodopsin: a Fourier transform infrared spectroscopy study of two-dimensional membrane crystals.

Aquaporins are integral membrane proteins found in diverse animal and plant tissues that mediate the permeability of plasma membranes to water molecules. Projection maps of two-dimensional crystals of aquaporin-1 (AQP1) reconstituted in lipid membranes suggested the presence of six to eight transmembrane helices in the protein. However, data from other sequence and spectroscopic analyses indicate that this protein may adopt a porin-like beta-barrel fold. In this paper, we use Fourier transform infrared spectroscopy to characterize the secondary structure of highly purified native and proteolyzed AQP1 reconstituted in membrane crystalline arrays and compare it to bacteriorhodopsin. For this analysis the fractional secondary structure contents have been determined by using several different algorithms. In addition, a neural network-based evaluation of the Fourier transform infrared spectra in terms of numbers of secondary structure segments and their interconnections [sij] has been performed. The following conclusions were reached: 1) AQP1 is a highly helical protein (42-48% alpha-helix) with little or no beta-sheet content. 2) The alpha-helices have a transmembrane orientation, but are more tilted (21 degrees or 27 degrees, depending on the considered refractive index) than the bacteriorhodopsin helices. 3) The helices in AQP1 undergo limited hydrogen/deuterium exchange and thus are not readily accessible to solvent. Our data support the AQP1 structural model derived from sequence prediction and epitope insertion experiments: AQP1 is a protein with at least six closely associated alpha-helices that span the lipid membrane.     

Proline residues in transmembrane alpha helices affect the folding of bacteriorhodopsin

Proline residues occur frequently in transmembrane alpha helices, which contrasts with their behaviour as helix-breakers in water-soluble proteins. The three membrane-embedded proline residues of bacteriorhodopsin have been replaced individually by alanine and glycine to give P50A, or P50G on helix B, P91A, or P91G on helix C, and P186A or P186G on helix F, and the effect on the protein folding kinetics has been investigated. The rate-limiting apoprotein folding step, which results in formation of a seven transmembrane, alpha helical state, was slower than wild-type protein for the Pro50 and Pro91 mutants, regardless of whether they were mutated to Ala or Gly. These proline residues give rise to several inter-helix contacts, which are therefore important in folding to the seven transmembrane helix state. No evidence for cis-trans isomerisations of the peptidyl prolyl bonds was found during this rate-limiting apoprotein folding step. Mutations of all three membrane-embedded proline residues affected the subsequent retinal binding and final folding to bacteriorhodopsin, suggesting that these proline residues contribute to formation of the retinal binding pocket within the helix bundle, again via helix/helix interactions. These results point to proline residues in transmembrane alpha helices being important in the folding of integral membrane proteins. The helix/helix interactions and hydrogen bonds that arise from the presence of proline residues in transmembrane alpha helices can affect the formation of transmembrane alpha helix bundles as well as cofactor binding pockets.     

Interactions of alpha-helices with lipid bilayers: a review of simulation studies

Membrane proteins, of which the majority seem to contain one or more alpha-helix, constitute approx. 30% of most genomes. A complete understanding of the nature of helix/bilayer interactions is necessary for an understanding of the structural principles underlying membrane proteins. This review describes computer simulation studies of helix/bilayer interactions. Key experimental studies of the interactions of alpha-helices and lipid bilayers are briefly reviewed. Surface associated helices are found in some membrane-bound enzymes (e.g. prostaglandin synthase), and as stages in the mechanisms of antimicrobial peptides and of pore-forming bacterial toxins. Transmembrane alpha-helices are found in most integral membrane proteins, and also in channels formed by amphipathic peptides or by bacterial toxins. Mean field simulations, in which the lipid bilayer is approximated as a hydrophobic continuum, have been used in studies of membrane-active peptides (e.g. alamethicin, melittin, magainin and dermaseptin) and of simple membrane proteins (e.g. phage Pf1 coat protein). All atom molecular dynamics simulations of fully solvated bilayers with transmembrane helices have been applied to: the constituent helices of bacteriorhodopsin; peptide-16 (a simple model TM helix); and a number of pore-lining helices from ion channels. Surface associated helices (e.g. melittin and dermaseptin) have been simulated, as have alpha-helical bundles such as bacteriorhodopsin and alamethicin. From comparison of the results from the two classes of simulation, it emerges that a major theoretical challenge is to exploit the results of all atom simulations in order to improve the mean field approach.     

Molecular dynamics of individual alpha-helices of bacteriorhodopsin in dimyristol phosphatidylocholine. I. Structure and dynamics.

Understanding the role of the lipid bilayer in membrane protein structure and dynamics is needed for tertiary structure determination methods. However, the molecular details are not well understood. Molecular dynamics computer calculations can provide insight into these molecular details of protein:lipid interactions. This paper reports on 10 simulations of individual alpha-helices in explicit lipid bilayers. The 10 helices were selected from the bacteriorhodopsin structure as representative alpha-helical membrane folding components. The bilayer is constructed of dimyristoyl phosphatidylcholine molecules. The only major difference between simulations is the primary sequence of the alpha-helix. The results show dramatic differences in motional behavior between alpha-helices. For example, helix A has much smaller root-mean-squared deviations than does helix D. This can be understood in terms of the presence of aromatic residues at the interface for helix A that are not present in helix D. Additional motions are possible for the helices that contain proline side chains relative to other amino acids. The results thus provide insight into the types of motion and the average structures possible for helices within the bilayer setting and demonstrate the strength of molecular simulations in providing molecular details that are not directly visualized in experiments.     

Kinetics of an individual transmembrane helix during bacteriorhodopsin folding

The kinetics of an individual helix of bacteriorhodopsin have been monitored during folding of the protein into lipid bilayer vesicles. A fluorescence probe was introduced at individual sites throughout helix D of bacteriorhodopsin and the changes in the fluorescence of the label were time-resolved. Partially denatured, labelled bacteriorhodopsin in SDS was folded directly into phosphatidylcholine lipid vesicles. Stopped-flow mixing of the reactants allowed the folding kinetics to be monitored with millisecond time resolution by time-resolving changes in the label fluorescence, intrinsic protein fluorescence as well as in the absorption of the retinal chromophore. Monitoring specific positions on helix D showed that two kinetic phases were altered compared to those determined by monitoring the average protein behaviour. These two phases, of 6.7 s(-1) and 0.33 s(-1), were previously assigned to formation of a key apoprotein intermediate during bacteriorhodopsin folding. The faster 6.7s(-1) phase was missing when time-resolving fluorescence changes of labels attached to the middle of helix D. The amplitude of the 0.33 s(-1) phase increased along the helix, as single labels were attached in turn from the cytoplasmic to the extracellular side. An interpretation of these results is that the 6.7 s(-1) phase involves partitioning of helix D within the lipid headgroups of the bilayer vesicle, while the 0.33 s(-1) phase could reflect transmembrane insertion of this helix. In addition, a single site on helix G was monitored during folding. The results indicate that, unlike helix D, the insertion of helix G cannot be differentiated from the average protein behaviour. The data show that, while folding of bacteriorhodopsin from SDS into lipids is a co-operative process, it is nevertheless possible to obtain information on specific regions of a membrane protein during folding in vitro.     

Structure of the integrin alphaIIb transmembrane segment

Integrin cell-adhesion receptors transduce signals bidirectionally across the plasma membrane via the single-pass transmembrane segments of each alpha and beta subunit. While the beta3 transmembrane segment consists of a linear 29-residue alpha-helix, the structure of the alphaIIb transmembrane segment reveals a linear 24-residue alpha-helix (Ile-966 -Lys-989) followed by a backbone reversal that packs Phe-992-Phe-993 against the transmembrane helix. The length of the alphaIIb transmembrane helix implies the absence of a significant transmembrane helix tilt in contrast to its partnering beta3 subunit. Sequence alignment shows Gly-991-Phe-993 to be fully conserved among all 18 human integrin alpha subunits, suggesting that their unusual structural motif is prototypical for integrin alpha subunits. The alphaIIb transmembrane structure demonstrates a level of complexity within the membrane that is beyond simple transmembrane helices and forms the structural basis for assessing the extent of structural and topological rearrangements upon alphaIIb-beta3 association, i.e. integrin transmembrane signaling.     

Simulation of the packing of idealized transmembrane α-helix bundles

The aim of this study is to investigate if the packing motifs of native transmembrane helices can be produced by simulations with simple potentials and to develop a method for the rapid generation of initial candidate models for integral membrane proteins composed of bundles of transmembrane helices. Constituent residues are mapped along the helix axis in order to maintain the amino acid sequence-dependent properties of the helix. Helix packing is optimized according to a semi-empirical potential mainly composed of four components: a bilayer potential, a crossing angle potential, a helix dipole potential and a helix-helix distance potential. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. Preliminary testing of the method has been conducted with idealized seven Ala₂₀ helix bundles. The structures generated show a high degree of compactness. It was observed that both bacteriorhodopsin-like and δ-endotoxin-like structures are generated in seven-helix bundle simulations, within which the composition varies dependent upon the cooling rate. The simulation method has also been employed to explore the packing of N = 4 and N = 12 transmembrane helix bundles. The results suggest that seven and 12 transmembrane helix bundles resembling those observed experimentally (e.g., bacteriorhodopsin, rhodopsin and cytochrome c oxidase subunit I) may be generated by simulations using simple potentials. Received: 16 November 1998 / Revised version: 26 March 1999 / Accepted: 8 April 1999     

Proline kinks in transmembrane alpha-helices

Integral membrane proteins often contain proline residues in their presumably alpha-helical transmembrane segments. This is in marked contrast to globular proteins, where proline is rarely found inside alpha-helices. Proline residues cause kinks in helices, and, in addition to leaving the i-4 backbone carbonyl without its normal hydrogen bond donor, also sterically prevent the (i-3)-carbonyl-(i + l)-amide backbone hydrogen bond from forming. Here, some structural aspects of proline kinks in transmembrane helices are discussed on the basis of an analysis of Pro-kinked helices in the photosynthetic reaction center and bacteriorhodopsin, as well as results from an analysis of Pro-containing transmembrane segments identified in the NBRF Protein Sequence Databank.     

Unfolding a transmembrane helix dimer: A FRET study in mixed micelles

The exact nature of membrane protein folding and assembly is not understood in detail yet. Addition of SDS to a membrane protein dissolved in mild, non-polar detergent results in formation of mixed micelles and in subsequent denaturation of higher ordered membrane protein structures. The exact nature of this denaturation event is, however, enigmatic, and separation of an individual helix pair in mixed micelles has also not been reported yet. Here we followed unfolding of the human glycophorin A transmembrane helix dimer in mixed micelles by fluorescence spectroscopy. Energy transfer between differently labelled glycophorin A transmembrane helices decreased with increasing SDS mole fractions albeit without modifying the helicity of the peptides. The energetics and kinetics of the dimer dissociation in mixed micelles is analyzed and discussed, and the experimental data demonstrate that mixed micelles can be used as a general method to investigate unfolding of α-helical membrane proteins.     

Interaction of transmembrane helices by a knobs-into-holes packing characteristic of soluble coiled coils

Membrane-embedded protein domains frequently exist as α-helical bundles, as exemplified by photosynthetic reaction centers, bacteriorhodopsin, and cytochrome C oxidase. The sidechain packing between their transmembrane helices was investigated by a nearest-neighbor analysis which identified sets of interfacial residues for each analyzed helix–helix interface. For the left-handed helix–helix pairs, the interfacial residues almost exclusively occupy positions a, d, e, or g within a heptad motif (abcdefg) which is repeated two to three times for each interacting helical surface. The connectivity between the interfacial residues of adjacent helices conforms to the knobs-into-holes type of sidechain packing known from soluble coiled coils. These results demonstrate on a quantitative basis that the geometry of sidechain packing is similar for left-handed helix–helix pairs embedded in membranes and coiled coils of soluble proteins. The transmembrane helix–helix interfaces studied are somewhat less compact and regular as compared to soluble coiled coils and tolerate all hydrophobic amino acid types to similar degrees. The results are discussed with respect to previous experimental findings which demonstrate that specific interactions between transmembrane helices are important for membrane protein folding and/or oligomerization. Proteins 31:150–159, 1998. © 1998 Wiley-Liss, Inc.     

Identification of Helical Packing Motifs Common to Bacteriorhodopsin and G Protein-Coupled Receptors

A sequence analysis and comparison of transmembrane helices in bacteriorhodopsin (BR) and G protein-coupled receptors (GPCRs) is presented to identify potential regions of homology across protein families. The results show a common pattern of residues is conserved within the interhelical contact regions of BR that fit a knob-into-hole structural motif previously postulated for globular proteins and photosynthetic reaction centers. Based on an alignment of conserved prolines in transmembrane helices, it is inferred that analogous helix packing arrangements are possible in the rhodopsin-like GPCRs. Molecular models of GPCR helices V and VI indicate these interactions occur between aromatic and hydrophobic residues flanking the highly conserved prolines in these sequences. A similar packing arrangement is shown to occur in the X-ray structure of the melittin which also displays a unique pairing of proline-linked helices. The contact pattern identified is further applied to predict the packing of pairs of proline-containing helices in the pheromone-like and cAMP GPCRs. A potential role in stabilizing structure formation is also suggested for the contacts. The results and conclusions are supported by recent biophysical studies of zinc binding to kappa-opioid receptor mutants.     

Identification of Helical Packing Motifs Common to Bacteriorhodopsin and G Protein-Coupled Receptors

A sequence analysis and comparison of transmembrane helices in bacteriorhodopsin (BR) and G protein-coupled receptors (GPCRs) is presented to identify potential regions of homology across protein families. The results show a common pattern of residues is conserved within the interhelical contact regions of BR that fit a knob-into-hole structural motif previously postulated for globular proteins and photosynthetic reaction centers. Based on an alignment of conserved prolines in transmembrane helices, it is inferred that analogous helix packing arrangements are possible in the rhodopsin-like GPCRs. Molecular models of GPCR helices V and VI indicate these interactions occur between aromatic and hydrophobic residues flanking the highly conserved prolines in these sequences. A similar packing arrangement is shown to occur in the X-ray structure of the melittin which also displays a unique pairing of proline-linked helices. The contact pattern identified is further applied to predict the packing of pairs of proline-containing helices in the pheromone-like and cAMP GPCRs. A potential role in stabilizing structure formation is also suggested for the contacts. The results and conclusions are supported by recent biophysical studies of zinc binding to kappa-opioid receptor mutants.     

Assembly of transmembrane helices of simple polytopic membrane proteins from sequence conservation patterns

The transmembrane (TM) domains of most membrane proteins consist of helix bundles. The seemingly simple task of TM helix bundle assembly has turned out to be extremely difficult. This is true even for simple TM helix bundle proteins, i.e., those that have the simple form of compact TM helix bundles. Herein, we present a computational method that is capable of generating native-like structural models for simple TM helix bundle proteins having modest numbers of TM helices based on sequence conservation patterns. Thus, the only requirement for our method is the presence of more than 30 homologous sequences for an accurate extraction of sequence conservation patterns. The prediction method first computes a number of representative well-packed conformations for each pair of contacting TM helices, and then a library of tertiary folds is generated by overlaying overlapping TM helices of the representative conformations. This library is scored using sequence conservation patterns, and a subsequent clustering analysis yields five final models. Assuming that neighboring TM helices in the sequence contact each other (but not that TM helices A and G contact each other), the method produced structural models of Calpha atom root-mean-square deviation (CA RMSD) of 3-5 A from corresponding crystal structures for bacteriorhodopsin, halorhodopsin, sensory rhodopsin II, and rhodopsin. In blind predictions, this type of contact knowledge is not available. Mimicking this, predictions were made for the rotor of the V-type Na(+)-adenosine triphosphatase without such knowledge. The CA RMSD between the best model and its crystal structure is only 3.4 A, and its contact accuracy reaches 55%. Furthermore, the model correctly identifies the binding pocket for sodium ion. These results demonstrate that the method can be readily applied to ab initio structure prediction of simple TM helix bundle proteins having modest numbers of TM helices.     

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     

A Continuum theory for the prediction of lateral and rotational positioning of α-helices in membrane proteins: Bacteriorhodopsin

We have developed a new method for the prediction of the lateral and the rotational positioning of transmembrane helices, based upon the present status of knowledge about the dominant interaction of the tertiary structure formation. The basic assumption about the interaction is that the interhelix binding is due to the polar interactions and that very short extramembrane loop segments restrict the relative position of the helices. Another assumption is made for the simplification of the prediction that a helix may be regarded as a continuum rod having polar interaction fields around it. The polar interaction field is calculated by a probe helix method, using a copolymer of serine and alanine as probe helices. The lateral position of helices is determined by the strength of the interhelix binding estimated from the polar interaction field together with the length of linking loop segments. The rotational positioning is determined by the polar interaction field, assuming the optimum lateral configuration. The structural change due to the binding of a prosthetic group is calculated, fixing the rotational freedom of a helix that is connected to the prosthetic group. Applying this method to bacteriorhodopsin, the optimum lateral and rotational positioning of transmembrane helices that are very similar to the experimental configuration was obtained. This method was implemented by a software system, which was developed for this work, and automatic calculation became possible for membrane proteins comprised of several transmembrane helices. © 1995 Wiley-Liss, Inc.