How do helix–helix interactions help determine the folds of membrane proteins? Perspectives from the study of homo‐oligomeric helical bundles

The final, structure-determining step in the folding of membrane proteins involves the coalescence of preformed transmembrane helices to form the native tertiary structure. Here, we review recent studies on small peptide and protein systems that are providing quantitative data on the interactions that drive this process. Gel electrophoresis, analytical ultracentrifugation, and fluorescence resonance energy transfer (FRET) are useful methods for examining the assembly of homo-oligomeric transmembrane helical proteins. These methods have been used to study the assembly of the M2 proton channel from influenza A virus, glycophorin, phospholamban, and several designed membrane proteins-all of which have a single transmembrane helix that is sufficient for association into a transmembrane helical bundle. These systems are being studied to determine the relative thermodynamic contributions of van der Waals interactions, conformational entropy, and polar interactions in the stabilization of membrane proteins. Although the database of thermodynamic information is not yet large, a few generalities are beginning to emerge concerning the energetic differences between membrane and water-soluble proteins: the packing of apolar side chains in the interior of helical membrane proteins plays a smaller, but nevertheless significant, role in stabilizing their structure. Polar, hydrogen-bonded interactions occur less frequently, but, nevertheless, they often provide a strong driving force for folding helix-helix pairs in membrane proteins. These studies are laying the groundwork for the design of sequence motifs that dictate the association of membrane helices.     

Turns in transmembrane helices: determination of the minimal length of a "helical hairpin" and derivation of a fine-grained turn propensity scale.

We have recently reported a first experimental turn propensity scale for transmembrane helices. This scale was derived from measurements of how efficiently a given residue placed in the middle of a 40 residue poly(Leu) stretch induces the formation of a "helical hairpin" with two rather than one transmembrane segment. We have now extended these studies, and have determined the minimum length of a poly(Leu) stretch compatible with the formation of a helical hairpin. We have also derived a more fine-grained turn propensity scale by (i) introducing each of the 20 amino acid residues into the middle of the shortest poly(Leu) stretch compatible with helical hairpin formation, and (ii) introducing pairs of residues in the middle of the 40 residue poly(Leu) stretch. The new turn propensities are consistent with the amino acid frequencies found in short hairpin loops in membrane proteins of known 3D structure.     

Interactions involved in the realignment of membrane-associated helices. An investigation using oriented solid-state NMR and attenuated total reflection Fourier transform infrared spectroscopies

A series of histidine-containing peptides (LAH4X6) was designed to investigate the membrane interactions of selected side chains. To this purpose, their pH-dependent transitions from in-plane to transmembrane orientations were investigated by attenuated total reflection Fourier transform infrared and oriented solid-state NMR spectroscopies. Peptides of the same family have previously been shown to exhibit antibiotic and DNA transfection activities. Solution NMR spectroscopy indicates that these peptides form amphipathic helical structures in membrane environments, and the technique was also used to characterize the pK values of all histidines in the presence of detergent micelles. Whereas one face of the amphipathic helix is clearly hydrophobic, the opposite side is flanked by four histidines surrounding six leucine, alanine, glycine, tryptophan, or tyrosine residues, respectively. This diversity in peptide composition causes pronounced shifts in the midpoint pH of the in-plane to transmembrane helical transition, which is completely abolished for the peptides carrying the most hydrophilic amino acid residues. These properties open up a conceptually new approach to study in a quantitative manner the hydrophobic as well as specific interactions of amino acids in membranes. Notably, the resulting scale for whole residue transitions from the bilayer interface to the hydrophobic membrane interior is obtained from extended helical sequences in lipid bilayers.     

Isolated transmembrane helices arranged across a membrane: computational studies.

A computational procedure for predicting the arrangement of an isolated helical fragment across a membrane was developed. The procedure places the transmembrane helical segment into a model triple-phase system 'water-octanol-water'; pulls the segment through the membrane, varying its 'global' position as a rigid body; optimizes the intrahelical and solvation energies in each global position by 'local' coordinates (dihedral angles of side chains); and selects the lowest energy global position for the segment. The procedure was applied to 45 transmembrane helices from the photosynthetic reaction center from Rhodopseudomonas viridis, cytochrome c oxidase from Paracoccus denitrificans and bacteriorhodopsin. In two thirds of the helical fragments considered, the procedure has predicted the vertical shifts of the fragments across the membrane with an accuracy of -0.15 +/- 3.12 residues compared with the experimental data. The accuracy for the remaining 15 fragments was 2.17 +/- 3.07 residues, which is about half of a helix turn. The procedure predicts the actual membrane boundaries of transmembrane helical fragments with greater accuracy than existing statistical methods. At the same time, the procedure overestimates the tilt values for the helical fragments.     

Dissecting membrane protein architecture: An annotation of structural complexity

α‐Helical membrane proteins exist in an anisotropic environment which strongly influences their folding, stability, and architecture, which is far more complex than a simple bundle of transmembrane helices, notably due to helix deformations, prosthetic groups and extramembrane structures. However, the role and the distribution of such heterogeneity in the supra molecular organization of membrane proteins remains poorly investigated. Using a nonredundant subset of α‐helical membrane proteins, we have annotated and analyze the statistics of several types of new elements such as incomplete helices, intramembrane loops, helical extensions of helical transmembrane domains, extracellular loops, and helices lying parallel to the membrane surface. The relevance of the annotation scheme was studied using residue composition, statistics, physical chemistry, and symmetry of their distribution in relation to the immediate membrane environment. Calculation of hydrophobicity using different scales show that different structural elements appear to have affinities coherent with their position in the membrane. Examination of the annotation scheme suggests that there is considerable information content in the amino acid compositions of the different elements suggesting that it might be useful for structural prediction. More importantly, the proposed annotation will help to decipher the complex hierarchy of interactions involved in membrane protein architecture. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 815–829, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

A simple method for predicting transmembrane alpha helices with better accuracy.

The prediction of a protein's structure from its amino acid sequence has been a long-standing goal of molecular biology. In this work, a new set of conformational parameters for membrane spanning alpha helices was developed using the information from the topology of 70 membrane proteins. Based on these conformational parameters, a simple algorithm has been formulated to predict the transmembrane alpha helices in membrane proteins. A FORTRAN program has been developed which takes the amino acid sequence as input and gives the predicted transmembrane alpha-helices as output. The present method correctly identifies 295 transmembrane helical segments in 70 membrane proteins with only two overpredictions. Furthermore, this method predicts all 45 transmembrane helices in the photosynthetic reaction center, bacteriorhodopsin and cytochrome c oxidase to an 86% level of accuracy and so is better than all other methods published to date.     

Riding the wave: structural and energetic principles of helical membrane proteins

Genome sequencing efforts have revealed that perhaps as many as 20-40% of open reading frames in complex organisms may encode proteins containing at least one helical transmembrane segment. Contrasting with this approaching tidal wave of helical membrane proteins is the fact that our understanding of the sequence-structure-function relationships for membrane proteins lags far behind that of soluble proteins. This looming reality emphasizes the tremendous biochemical and structural work that remains to be done on helical membrane proteins in order to elucidate the structural and energetic principles that specify and stabilize their folds, which define their functions. These facts are not lost on the pharmaceutical industry, where successful therapeutics and major discovery efforts are targeting membrane proteins.     

The Transmembrane Helices of Beef Heart Cytochrome Oxidase

The locations of the transmembrane helices in the 12 subunits of beef heart cytochrome oxidase were predicted with a modified form of the von Heijne-Blomberg hydrophobicity scale. Based on ~20 residues per transmembrane helix, about 480 of the estimated 660 helical residues (36.8% of 1,793 total residues) are expected to be in transmembrane helices that have their axes tilted by a small angle α from the normal to the plane of the membrane. This angle is calculated to be ~30°, based on the observed overall tilt angle θ of 39° obtained from circular dichroism (CD) measurements on multilamellar films, or about 25°, based on the observed tilt angle θ of 36° obtained from the infrared linear dichroism of films. For 21 residues per transmembrane helix, the calculated values of α become 32° and 28°, respectively, depending upon the value of θ used. Thus, a transmembrane helical tilt angle of ~30° accounts for the predicted transmembrane stretches in cytochrome oxidase if 20-21 residues are sufficient to span the membrane. Additional helical residues in the lipid head region may deviate by a larger angle from the normal to the plane of the membrane in cytochrome oxidase.     

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.     

Disposition of polar and nonpolar residues on outer surfaces of transmembrane helical segments of proteins involved in proton translocation

"Helical wheel" projections of transmembrane helical segments of membrane proteins involved in proton translocation were constructed. The particular proteins studied were the uncF protein subunit of the Escherichia coli proton-ATPase, the uncE protein subunit of the E. coli proton-ATPase, and cytochrome oxidase subunit III. Clear demarcation of polar and nonpolar regions on surfaces of transmembrane helical segments was seen in the uncF protein and in uncE protein helical segment two, but not in uncE protein helical segment one. The transmembrane segment of cytochrome oxidase subunit III which includes the dicyclohexylcarbodiimide (DCCD)-reactive residue was very similar to E. coli uncE protein helical segment two. The DCCD-reactive residue in both was clearly located on a nonpolar surface.     

Transmembrane and extramembrane contributions to membrane protein thermal stability: studies with the NaChBac sodium channel

The thermal stabilities of the extramembranous and transmembranous regions of the bacterial voltage-gated sodium channel NaChBac have been characterised using thermal-melt synchrotron radiation circular dichroism (SRCD) spectroscopy. A series of constructs, ranging from the full-length protein containing both the C-terminal cytoplasmic and the transmembranous domains, to proteins with decreasing amounts of the cytoplasmic domain, were examined in order to separately define the roles of these two types of domains in the stability and processes of unfolding of a membrane protein. The sensitivity of the SRCD measurements over a wide range of wavelengths and temperatures has meant that subtle but reproducible conformational changes could be detected with accuracy. The residues in the C-terminal extramembranous domain were highly susceptible to thermal denaturation, but for the most part the transmembrane residues were not thermally-labile and retained their helical character even at very elevated temperatures. The process of thermal unfolding involved an initial irreversible unfolding of the highly labile distal extramembranous C-terminal helical region, which was accompanied by a reversible unfolding of a small number of helical residues in the transmembrane domain. This was then followed by the irreversible unfolding of a limited number of additional transmembrane helical residues at greatly elevated temperatures. Hence this study has been able to determine the different contributions and roles of the transmembrane and extramembrane residues in the processes of thermal denaturation of this multipass integral membrane protein.     

Structural comparison and classification of alpha‐helical transmembrane domains based on helix interaction patterns

Structural classification of membrane proteins is still in its infancy due to the relative paucity of available three‐dimensional structures compared with soluble proteins. However, recent technological advances in protein structure determination have led to a significant increase in experimentally known membrane protein folds, warranting exploration of the structural universe of membrane proteins. Here, a new and completely membrane protein specific structural classification system is introduced that classifies α‐helical membrane proteins according to common helix architectures. Each membrane protein is represented by a helix interaction graph depicting transmembrane helices with their pairwise interactions resulting from individual residue contacts. Subsequently, proteins are clustered according to similarities among these helix interaction graphs using a newly developed structural similarity score called HISS. As HISS scores explicitly disregard structural properties of loop regions, they are more suitable to capture conserved transmembrane helix bundle architectures than other structural similarity scores. Importantly, we are able to show that a classification approach based on helix interaction similarity closely resembles conventional structural classification databases such as SCOP and CATH implying that helix interactions are one of the major determinants of α‐helical membrane protein folds. Furthermore, the classification of all currently available membrane protein structures into 20 recurrent helix architectures and 15 singleton proteins demonstrates not only an impressive variability of membrane helix bundles but also the conservation of common helix interaction patterns among proteins with distinctly different sequences. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Formation of cytoplasmic turns between two closely spaced transmembrane helices during membrane protein integration into the ER membrane

The helical hairpin, two closely spaced transmembrane helices separated by a short turn, is a recurring structural element in integral membrane proteins, and may serve as a compact unit that inserts into the membrane en bloc. Previously, we have determined the propensities of the 20 natural amino acids, when present in the middle of a long hydrophobic stretch, to induce the formation of a helical hairpin with a lumenally exposed turn during membrane protein assembly into the endoplasmic reticulum membrane. Here, we present results from a similar set of measurements, but with the turn placed on the cytoplasmic side of the membrane. We find that a significantly higher number of turn-promoting residues need to be present to induce a cytoplasmic turn compared to a lumenal turn, and that, in contrast to the lumenal turn, the positively charged residues Arg and Lys are the strongest turn-promoters in cytoplasmic turns. These results suggest that the process of turn formation between transmembrane helices is different for lumenal and cytoplasmic turns.     

A new coiled hollow-fiber module design for enhanced microfiltration performance in biotechnology.

The microfiltration performance of a novel membrane module design with helically wound hollow fibers is compared with that obtained with a standard commercial-type crossflow module containing linear hollow fibers. Cell suspensions (yeast, E. coli, and mammalian cell cultures) commonly clarified in the biotechnology industry are used for this comparison. The effect of variables such as transmembrane pressure, particle suspension concentration, and feed flow rate on membrane performance is evaluated. Normalized permeation fluxes versus flow rate or Dean number behave according to a heat transfer correlation obtained with centrifugal instabilities of the Taylor type. The microfiltration performance of this new module design, which uses secondary flows in helical tubes, is significantly better than an equivalent current commercial crossflow module when filtering suspensions relevant to the biotechnology industry. Flux and capacity improvements of up to 3.2-fold (constant transmembrane pressure operation) and 3.9-fold (constant flux operation), respectively, were obtained with the helical module over those for the linear module.     

Probing the local secondary structure of bacteriophage S21 pinholin membrane protein using electron spin echo envelope modulation spectroscopy

There have recently been advances in methods for detecting local secondary structures of membrane protein using electron paramagnetic resonance (EPR). A three pulsed electron spin echo envelope modulation (ESEEM) approach was used to determine the local helical secondary structure of the small hole forming membrane protein, S²¹ pinholin. This ESEEM approach uses a combination of site-directed spin labeling and ²H-labeled side chains. Pinholin S²¹ is responsible for the permeabilization of the inner cytosolic membrane of double stranded DNA bacteriophage host cells. In this study, we report on the overall global helical structure using circular dichroism (CD) spectroscopy for the active form and the negative-dominant inactive mutant form of S²¹ pinholin. The local helical secondary structure was confirmed for both transmembrane domains (TMDs) for the active and inactive S²¹ pinholin using the ESEEM spectroscopic technique. Comparison of the ESEEM normalized frequency domain intensity for each transmembrane domain gives an insight into the α-helical folding nature of these domains as opposed to a π or 3₁₀-helix which have been observed in other channel forming proteins.     

Are membrane proteins "inside-out" proteins?

One of the central paradigms of structural biology is that membrane proteins are "inside-out" proteins, in that they have a core of polar residues surrounded by apolar residues. This is the reverse of the characteristics found in water-soluble proteins. We have decided to test this paradigm, now that sufficient numbers of transmembrane alpha-helical structures are accessible to statistical analysis. We have analyzed the correlation between accessibility and hydrophobicity of both individual residues and complete helices. Our analyses reveal that hydrophobicity of residues in a transmembrane helical bundle does not correlate with any preferred location and that the hydrophilic vector of a helix is a poor indicator of the solvent exposed face of a helix. Neither polar nor hydrophobic residues show any bias for the exterior or the interior of a transmembrane domain. As a control, analysis of water-soluble helical bundles performed in a similar manner has yielded clear correlations between hydrophobicity and accessibility. We therefore conclude that, based on the data set used, membrane proteins as "inside-out" proteins is an unfounded notion, suggesting that packing of alpha-helices in membranes is better understood by maximization of van der Waal's forces, rather than by a general segregation of hydrophobicities driven by lipid exclusion.     

Understanding peptide interactions with the lipid bilayer: a guide to membrane protein engineering

In recent years, studies of the interactions of designed hydrophobic and amphipathic polypeptides with biological membranes have progressed considerably. In-plane and transmembrane helical domains have been engineered as well as sequences that exhibit dynamic distributions of different topologies. These sequences not only help our understanding of the thermodynamic interaction contributions that determine peptide orientation, but also exhibit interesting biological functions as autonomous units. In addition, helices are considered to be folding intermediates of multi-spanning membrane proteins. The specificity and thermodynamic stability of transmembrane helix-helix interactions or the assembly of beta sheets at the membrane surface have, therefore, become a focus of ongoing research.     

An implicit membrane generalized born theory for the study of structure, stability, and interactions of membrane proteins

Exploiting recent developments in generalized Born (GB) electrostatics theory, we have reformulated the calculation of the self-electrostatic solvation energy to account for the influence of biological membranes. Consistent with continuum Poisson-Boltzmann (PB) electrostatics, the membrane is approximated as an solvent-inaccessible infinite planar low-dielectric slab. The present membrane GB model closely reproduces the PB electrostatic solvation energy profile across the membrane. The nonpolar contribution to the solvation energy is taken to be proportional to the solvent-exposed surface area (SA) with a phenomenological surface tension coefficient. The proposed membrane GB/SA model requires minor modifications of the pre-existing GB model and appears to be quite efficient. By combining this implicit model for the solvent/bilayer environment with advanced computational sampling methods, like replica-exchange molecular dynamics, we are able to fold and assemble helical membrane peptides. We examine the reliability of this model and approach by applications to three membrane peptides: melittin from bee venom, the transmembrane domain of the M2 protein from Influenza A (M2-TMP), and the transmembrane domain of glycophorin A (GpA). In the context of these proteins, we explore the role of biological membranes (represented as a low-dielectric medium) in affecting the conformational changes in melittin, the tilt of transmembrane peptides with respect to the membrane normal (M2-TMP), helix-to-helix interactions in membranes (GpA), and the prediction of the configuration of transmembrane helical bundles (GpA). The present method is found to perform well in each of these cases and is anticipated to be useful in the study of folding and assembly of membrane proteins as well as in structure refinement and modeling of membrane proteins where a limited number of experimental observables are available.