Asn- and Asp-mediated interactions between transmembrane helices during translocon-mediated membrane protein assembly

Inter-helix hydrogen bonding involving asparagine (Asn, N), glutamine (Gln, Q), aspartic acid (Asp, D) or glutamic acid (Glu, E) can drive efficient di- or trimerization of transmembrane helices in detergent micelles and lipid bilayers. Likewise, Asn-Asn and Asp-Asp pairs can promote the formation of helical hairpins during translocon-mediated membrane protein assembly in the endoplasmic reticulum. By in vitro translation of model integral membrane protein constructs in the presence of rough microsomes, we show that Asn- or Asp-mediated interactions with a neighbouring transmembrane helix can enhance the membrane insertion efficiency of a marginally hydrophobic transmembrane segment. Our observations suggest that inter-helix hydrogen bonds can form during Sec61 translocon-assisted insertion and thus could be important for membrane protein assembly.     

Differentiation between transmembrane helices and peripheral helices by the deconvolution of circular dichroism spectra of membrane proteins.

The interpretation of the circular dichroism (CD) spectra of proteins to date requires additional secondary structural information of the proteins to be analyzed, such as X-ray or NMR data. Therefore, these methods are inappropriate for a CD database whose secondary structures are unknown, as in the case of the membrane proteins. The convex constraint analysis algorithm (Perczel, A., Hollósi, M., Tusnády, G., & Fasman, G. D., 1991, Protein Eng. 4, 669-679), on the other hand, operates only on a collection of spectral data to extract the common spectral components with their spectral weights. The linear combinations of these derived "pure" CD curves can reconstruct the original data set with great accuracy. For a membrane protein data set, the five-component spectra so obtained from the deconvolution consisted of two different types of alpha helices (the alpha helix in the soluble domain and the alpha T helix, for the transmembrane alpha helix), a beta-pleated sheet, a class C-like spectrum related to beta turns, and a spectrum correlated with the unordered conformation. The deconvoluted CD spectrum for the alpha T helix was characterized by a positive red-shifted band in the range 195-200 nm (+95,000 deg cm2 dmol-1), with the intensity of the negative band at 208 nm being slightly less negative than that of the 222-nm band (-50,000 and -60,000 deg cm2 dmol-1, respectively) in comparison with the regular alpha helix, with a positive band at 190 nm and two negative bands at 208 and 222 nm with magnitudes of +70,000, -30,000, and -30,000 deg cm2 dmol-1, respectively.     

Sequence context strongly modulates association of polar residues in transmembrane helices

Polar residues are capable of mediating the association of membrane-embedded helices through the formation of side-chain/side-chain inter-helical hydrogen bonds. However, the extent to which native van der Waals packing of the residues surrounding the polar locus can enhance, or interfere with, the interaction of polar residues has not yet been studied. We examined the propensities of four polar residues (aspartic acid, asparagine, glutamic acid, and glutamine) to promote self-association of transmembrane (TM) domains in several biologically derived sequence environments, including (i). four naturally occurring TM domains that contain a Glu or Gln residue (Tnf5/CD40 ligand, C79a/Ig-alpha, C79b/Ig-beta, and Fut3/alpha-fucosyltransferase); and (ii). variants of bacteriophage M13 major coat protein TM segment with Asp and Asn at interfacial and non-interfacial positions. Self-association was quantified by the TOXCAT assay, which measures TM helix self-oligomerization in the Escherichia coli inner membrane. While an appropriately placed polar residue was found in several cases to significantly stabilize TM helix-helix interactions through the formation of an interhelical hydrogen bond, in other cases the strongly polar residues did not enhance the association of the two helices. Overall, these results suggest that an innate structural mechanism may operate to control non-specific association of membrane-embedded polar residues.     

A specific interface between integrin transmembrane helices and affinity for ligand

Conformational communication across the plasma membrane between the extracellular and intracellular domains of integrins is beginning to be defined by structural work on both domains. However, the role of the α and β subunit transmembrane domains and the nature of signal transmission through these domains have been elusive. Disulfide bond scanning of the exofacial portions of the integrin α_IIβ and β₃ transmembrane domains reveals a specific heterodimerization interface in the resting receptor. This interface is lost rather than rearranged upon activation of the receptor by cytoplasmic mutations of the α subunit that mimic physiologic inside-out activation, demonstrating a link between activation of the extracellular domain and lateral separation of transmembrane helices. Introduction of disulfide bridges to prevent or reverse separation abolishes the activating effect of cytoplasmic mutations, confirming transmembrane domain separation but not hinging or piston-like motions as the mechanism of transmembrane signaling by integrins.     

A periodicity analysis of transmembrane helices

Transmembrane helices and the helical bundles which they form are the major building blocks of membrane proteins. Since helices are characterized by a given periodicity, it is possible to search for patterns of traits which typify one side of the helix and not the other (e.g. amphipathic helices contain a polar and apolar sides). Using Fourier transformation we have analyzed solved membrane protein structures as well as sequences of membrane proteins from the Swiss-Prot database. The traits searched included aromaticity, volume and ionization. While a number of motifs were already recognized in the literature, many were not. One particular example involved helix VII of lactose permease which contains seven aromatic residues on six helical turns. Similarly six glycine residues in four consecutive helical turns were identified as forming a motif in the chloride channel. A tabulation of all the findings is presented as well as a possible rationalization of the function of the motif.     

Identifying interactions between transmembrane helices from the adenosine A2A receptor

The human adenosine A(2A) receptor (A(2A)R) is an integral membrane protein and a member of the G-protein-coupled receptor (GPCR) superfamily, characterized by seven transmembrane (TM) helices. Although helix-helix association in the lipid bilayer is known to be an essential step in the folding of GPCRs, the determinants of their structures, folding, and assembly in the cell membrane are poorly understood. Previous studies in our group showed that while peptides corresponding to all seven TM domains of A(2A)R form stable helical structures in detergent micelles and lipid vesicles, they display significant variability in their helical propensity. This finding suggested to us that some TM domains might need to interact with other domains to properly insert and fold in hydrophobic environments. In this study, we assessed the ability of TM peptides to interact in pairwise combinations. We analyzed peptide interactions in hydrophobic milieus using circular dichroism spectroscopy and F?rster resonance energy transfer. We find that specific interactions between TM helices occur, leading to additional helical content, especially in weakly helical TM domains, suggesting that some TM domains need a partner for proper folding in the membrane. The approach developed in this study will enable complete analysis of the TM domain interactions and the modeling of a folding pathway for A(2A)R.     

Spontaneous local membrane curvature induced by transmembrane proteins

The (local) curvature of cellular membranes acts as a driving force for the targeting of membrane-associated proteins to specific membrane domains, as well as a sorting mechanism for transmembrane proteins, e.g., by accumulation in regions of matching spontaneous curvature. The latter measure was previously experimentally employed to study the curvature induced by the potassium channel KvAP and by aquaporin AQP0. However, the direction of the reported spontaneous curvature levels as well as the molecular driving forces governing the membrane curvature induced by these integral transmembrane proteins could not be addressed experimentally.Here, using both coarse-grained and atomistic molecular dynamics (MD) simulations, we report induced spontaneous curvature values for the homologous potassium channel Kv 1.2/2.1 Chimera (KvChim) and AQP0 embedded in unrestrained lipid bicelles that are in very good agreement with experiment. Importantly, the direction of curvature could be directly assessed from our simulations: KvChim induces a strong positive membrane curvature ($\approx 0.036$ nm^?1) whereas AQP0 causes a comparably small negative curvature ($\approx ?0.019$ nm^?1).Analyses of protein-lipid interactions within the bicelle revealed that the potassium channel shapes the surrounding membrane via structural determinants. Differences in shape of the protein-lipid interface of the voltage-gating domains between the extracellular and cytosolic membrane leaflets induce membrane stress and thereby promote a protein-proximal membrane curvature. In contrast, the water pore AQP0 displayed a high structural stability and an only faint effect on the surrounding membrane environment that is connected to its wedge-like shape.     

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.     

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.     

A turn propensity scale for transmembrane helices.

Using a model protein with a 40 residue hydrophobic transmembrane segment, we have measured the ability of all the 20 naturally occurring amino acids to form a tight turn when placed in the middle of the hydrophobic segment. Turn propensities in a transmembrane helix are found to be markedly different from those of globular proteins, and in most cases correlate closely with the hydrophobicity of the residue. The turn propensity scale may be used to improve current methods for membrane protein topology prediction.     

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.     

A unifying mechanism accounts for sensing of membrane curvature by BAR domains,amphipathic helices and membrane-anchored proteins

The discovery of proteins that recognize membrane curvature created a paradigm shift by suggesting that membrane shape may act as a cue for protein localization that is independent of lipid or protein composition. Here we review recent data on membrane curvature sensing by three structurally unrelated motifs: BAR domains, amphipathic helices and membrane-anchored proteins. We discuss the conclusion that the curvature of the BAR dimer is not responsible for sensing and that the sensing properties of all three motifs can be rationalized by the physicochemical properties of the curved membrane itself. We thus anticipate that membrane curvature will promote the redistribution of proteins that are anchored in membranes through any type of hydrophobic moiety, a thesis that broadens tremendously the implications of membrane curvature for protein sorting, trafficking and signaling in cell biology.     

The measurement of transmembrane helices by the deconvolution of CD spectra of membrane proteins: A review

The interpretation of the CD spectra of proteins to date requires additional secondary structural information of the proteins to be analyzed, such as x-ray or nmr data. Therefore, these methods are inappropriate for a CD data base whose secondary structures are unknown, as in the case of the membrane proteins. The Convex Constraint Analysis algorithm [A. Perczel, M. Hollósi, G. Tusnády, and G. D. Fasman (1991) Protein Engineering, Vol. 4, 669–679], on the other hand, operates only on a collection of spectral data to extract the common spectral components with their spectral weights. The linear combinations of these derived “pure” CD curves can reconstruct the original data set with great accuracy. For a membrane protein data set, the five-component spectra so obtained from the deconvolution consisted of two different types of α-helices (the α-helix in the soluble domain and the α_T-helix, for the transmembrane α-helix), a β-pleated sheet, a class C-like spectrum related to β-turns, and a spectrum correlated with the unordered conformation. The deconvoluted CD spectrum for the α_T-helix was characterized by a positive red-shifted band in the range 195–200 nm (+95,000 deg cm² dmol^?l), with the intensity of the negative band at 208 nm being slightly less negative than that of the 222 nm band (?50,000 and ?60,000 deg cm² dmol^?1, respectively) in comparison with the regular α-helix, with a positive band at 190 nm and two negative bands at 208 and 222 nm with magnitudes of + 70,000, ?30,000, and ?30,000 deg cm² dmol^?1, respectively. © 1994 John Wiley & Sons, Inc.     

Coupling between lipid shape and membrane curvature

Using molecular dynamics simulations, we examine the behavior of lipids whose preferred curvature can be systematically varied. This curvature is imposed by controlling the headgroup size of a coarse-grained lipid model recently developed by us. To validate this approach, we examine self-assembly of each individual lipid type and observe the complete range of expected bilayer and micelle phases. We then examine binary systems consisting of lipids with positive and negative preferred curvature and find a definite sorting effect. Lipids with positive preferred curvature are found in greater proportions in outer monolayers with the opposite observed for lipids with negative preferred curvature. We also observe a similar, but slightly stronger effect for lipids in a developing spherical bud formed by adhesion to a colloid (e.g., a viral capsid). Importantly, the magnitude of this effect in both cases was large only for regions with strong mean curvature (radii of curvature <10 nm). Our results suggest that lipid shape must act in concert with other physico-chemical effects such as phase transitions or interactions with proteins to produce strong sorting in cellular pathways.     

The effect of nucleotide bias upon the composition and prediction of transmembrane helices

Transmembrane helices are the most readily predictable secondary structure components of proteins. They can be predicted to a high degree of accuracy in a variety of ways. Many of these methods compare new sequence data with the sequence characteristics of known transmembrane domains. However, the known transmembrane sequences are not necessarily representative of a particular organism. We attempt to demonstrate that parameters optimized for the known transmembrane domains are far from optimal when predicting transmembrane regions in a given genome. In particular, we have tested the effect of nucleotide bias upon the composition and hence the prediction characteristics of transmembrane helices. Our analysis shows that nucleotide bias of a genome has a strong and predictable influence upon the occurrences of several of the most important hydrophobic amino acids found within transmembrane helices. Thus, we show that nucleotide bias should be taken into account when determining putative transmembrane domains from sequence data.     

A Golgi-associated PDZ domain protein modulates cystic fibrosis transmembrane regulator plasma membrane expression.

We identified a novel cystic fibrosis transmembrane conductance regulator (CFTR)-associating, PDZ domain-containing protein, CAL (CFTR associated ligand) containing two predicted coiled-coiled domains and one PDZ domain. The PDZ domain of CAL binds to the C terminus of CFTR. Although CAL does not have any predicted transmembrane domains, CAL is associated with membranes mediated by a region containing the coiled-coil domains. CAL is located primarily at the Golgi apparatus, co-localizing with trans-Golgi markers and is sensitive to Brefeldin A treatment. Immunoprecipitation experiments suggest that CAL exists as a multimer. Overexpression of CAL reduces CFTR chloride currents in mammalian cells and decreases expression, rate of insertion and half-life of CFTR in the plasma membrane. The Na(+)/H(+) exchanger regulatory factor, NHE-RF, a subplasma membrane PDZ domain protein, restores cell surface expression of CFTR and chloride currents. In addition, NHE-RF inhibits the binding of CAL to CFTR. CAL modulates the surface expression of CFTR. CAL favors retention of CFTR within the cell, whereas NHE-RF favors surface expression by competing with CAL for the binding of CFTR. Thus, the regulation of CFTR in the plasma membrane involves the dynamic interaction between at least two PDZ domain proteins.     

GALLEX,a measurement of heterologous association of transmembrane helices in a biological membrane

Whereas a variety of two-hybrid systems are available to measure the interaction of soluble proteins, related methods are significantly less developed for the measurement of membrane protein interactions. Here we present a two-hybrid system to follow the heterodimerization of membrane proteins in the Escherichia coli inner membrane. The method is based on the repression of a reporter gene activity by two LexA DNA binding domains with different DNA binding specificities. When coupled to transmembrane domains, heterodimeric association is reported by repression of beta-galactosidase synthesis. The LexA-transmembrane chimeric proteins were found to correctly insert into the membrane, and reproducible signals were obtained measuring the homodimerization as well as heterodimerization of wild-type and mutant glycophorin A transmembrane helices. The GALLEX data were compared with data recently gained by other methods and discussed in the general context of heteroassociation of single TM helices. Additionally, the formation of heterodimers between the TM domains of the alpha(4) and the beta(7) integrin subunits were tested. The results show that both homo- and heterodimerization of membrane proteins can be measured accurately using the GALLEX system.     

Capping transmembrane helices of MscL with aromatic residues changes channel response to membrane stretch

Tyrosines and tryptophans that anchor both ends of the helices to membrane interfaces in many transmembrane proteins are not common in MscL and homologous mechanosensitive channels. This characteristic absence of two aromatic "belts" may be critical for MscL function as the opening transition is predicted to be associated with a strong helical reorientation. A single tyrosine (Y75) on the extracellular side of the M2 helix of pentameric EcoMscL is absent in TbMscL, which instead has a single tyrosine (Y87) on the cytoplasmic side of M2. Moving the tyrosine of EcoMscL to the intracellular side (Y75F/F93Y) or capping the TM2 helix on both sides (F93Y/W) slows the kinetics of gating and increases the threshold for activation, leading to a partial loss-of-function in osmotic shock survival assays. Increasing the distance between the caps (L98W, L102Y/W) partially restores channel function presumably by loosening restraints for tilting. Capping the TM2 helix with a charged residue (Y75E) causes a right shift of the activation curve ("stiff" phenotype) and abolishes function. Introducing a "cap" into the TM1 helix (I41W) decreases the activation threshold and shortens the mean open time but unexpectedly leads to a complete loss-of-function in vivo. The data are consistent with the view that restraining helical positions in MscL by introducing specific protein-lipid interactions at membrane interfaces compromises MscL function. Subtle differences in osmotic shock survival are more evident at low levels of mutant protein expression. We observed a correlation between the right shift of tension activation threshold and the loss-of-function channel phenotype, with a few exceptions that point to other parameters of gating that may define the osmotic rescuing ability in vivo.     

Effects of 'hydrophobic mismatch' on the location of transmembrane helices in the ER membrane

We have studied the effects of 'hydrophobic mismatch' between a poly-Leu transmembrane helix (TMH) and the ER membrane using a glycosylation mapping approach. The simplest interpretation of our results is that the lumenal end of the TMH is located deeper in the membrane for both short (negative mismatch) and long (positive mismatch) TMHs than for poly-Leu segments of intermediate length. We further find that the position-specific effect of Lys residues on the location of short TMHs in the membrane varies with an apparent helical periodicity when the Lys residue is moved along the poly-Leu stretch. We discuss these findings in the context of models for peptide-lipid interactions during hydrophobic mismatch.