Intrahelical hydrogen bonding of serine, threonine and cysteine residues within alpha-helices and its relevance to membrane-bound proteins

A survey of known protein structures reveals that approximately 70% of serine residues and at least 85% (potentially 100%) of threonine residues in helices make hydrogen bonds to carbonyl oxygen atoms in the preceding turn of the helix. The high frequency of intrahelical hydrogen bonding is of particular significance for intrinsic membrane-bound proteins that form transmembrane helices. Hydrogen bonding within a helix provides a way for serine, threonine and cysteine residues to satisfy their hydrogen-bonding potential permitting such residues to occur in helices buried within a hydrophobic milieu.     

Motifs of two small residues can assist but are not sufficient to mediate transmembrane helix interactions

Both experimental and statistical searches for specific motifs that mediate transmembrane helix-helix interactions showed that two glycine residues separated by three intervening residues (GxxxG) provide a framework for specific interactions. Further work suggested that other motifs of small residues can mediate the interaction of transmembrane domains, so that the AxxxA-motif could also drive strong interactions of alpha-helices in soluble proteins. Thus, all these data indicate that a motif of two small residues in a distance of four might be enough to provide a framework for transmembrane helix-helix interaction. To test whether GxxxG is equivalent to (small)xxx(small), we investigated the effect of a substitution of either of the two Gly residues in the glycophorin A GxxxG-motif by Ala or Ser using the recently developed GALLEX system. The results of this mutational study demonstrate that, while a replacement of either of the two Gly by Ala strongly disrupts GpA homo-dimerization, the mutation to Ser partly stabilizes a dimeric structure. We suggest that the Ser residue can form a hydrogen bond with a backbone carbonyl group of the adjacent helix stabilizing a preformed homo-dimer. While (small)xxx(small) serves as a useful clue, the context of adjacent side-chains is essential for stable helix interaction, so each case must be tested.     

The stability of transmembrane helix interactions measured in a biological membrane

Despite some promising progress in the understanding of membrane protein folding and assembly, there is little experimental information regarding the thermodynamic stability of transmembrane helix interactions and even less on the stability of transmembrane helix-helix interactions in a biological membrane. Here we describe an approach that allows quantitative measurement of transmembrane helix interactions in a biological membrane, and calculation of changes in the interaction free energy resulting from substitution of single amino acids. Dimerization of several variants of the glycophorin A transmembrane domain are characterized and compared to the wild-type (wt) glycophorin A transmembrane helix dimerization. The calculated DeltaDeltaG(app) values are further compared with values found in the literature. In addition, we compare interactions between the wt glycophorin A transmembrane domain and helices in which critical glycine residues are replaced by alanine or serine, respectively. The data demonstrate that replacement of the glycine residues by serine is less destabilizing than replacement by alanine with a DeltaDeltaG(app) value of about 0.4 kcal/mol. Our study comprises the first measurement of a transmembrane helix interaction in a biological membrane, and we are optimistic that it can be further developed and applied.     

Re-Introduction of Transmembrane Serine Residues Reduce the Minimum Pore Diameter of Channelrhodopsin-2

Channelrhodopsin-2 (ChR2) is a microbial-type rhodopsin found in the green algae Chlamydomonas reinhardtii. Under physiological conditions, ChR2 is an inwardly rectifying cation channel that permeates a wide range of mono- and divalent cations. Although this protein shares a high sequence homology with other microbial-type rhodopsins, which are ion pumps, ChR2 is an ion channel. A sequence alignment of ChR2 with bacteriorhodopsin, a proton pump, reveals that ChR2 lacks specific motifs and residues, such as serine and threonine, known to contribute to non-covalent interactions within transmembrane domains. We hypothesized that reintroduction of the eight transmembrane serine residues present in bacteriorhodopsin, but not in ChR2, will restrict the conformational flexibility and reduce the pore diameter of ChR2. In this work, eight single serine mutations were created at homologous positions in ChR2. Additionally, an endogenous transmembrane serine was replaced with alanine. We measured kinetics, changes in reversal potential, and permeability ratios in different alkali metal solutions using two-electrode voltage clamp. Applying excluded volume theory, we calculated the minimum pore diameter of ChR2 constructs. An analysis of the results from our experiments show that reintroducing serine residues into the transmembrane domain of ChR2 can restrict the minimum pore diameter through inter- and intrahelical hydrogen bonds while the removal of a transmembrane serine results in a larger pore diameter. Therefore, multiple positions along the intracellular side of the transmembrane domains contribute to the cation permeability of ChR2.     

Structural role of glycine in amyloid fibrils formed from transmembrane alpha-helices

Amyloid fibrils associated with diseases such as Alzheimer's are often derived from the transmembrane helices of membrane proteins. It is known that the fibrils have a cross-beta-sheet structure where main chain hydrogen bonding occurs between beta-strands in the direction of the fibril axis. However, the structural basis for how the membrane-spanning helix is converted into a beta-sheet or how protofibrils associate into fibrils is not known. Here, we use a model peptide corresponding to a portion of the single transmembrane helix of glycophorin A to investigate the structural role of glycine in amyloid-like fibrils formed from transmembrane helices. Glycophorin A contains a GxxxG motif that is found in many transmembrane sequences including that of the amyloid precursor protein and prion protein. We propose that glycine, which mediates helix interactions in membrane proteins, also provides key packing motifs when it occurs in beta-sheets. We show that glycines in the glycophorin A transmembrane helix promote extended beta-strand formation when the helix partitions into aqueous environments and stabilize the packing of beta-sheets in the formation of amyloid-like fibrils. We demonstrate that fibrillization can be disrupted with a new class of inhibitors that target the molecular grooves created by glycine.     

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.     

Computer modeling of polyleucine-based coiled coil dimers in a realistic membrane environment: insight into helix-helix interactions in membrane proteins

Simulated annealing was performed to model parallel dimers of alpha-helical transmembrane peptides with the sequence L(11)XL(12), predicting left-handed coiled coil geometry in all cases. Insertion of peptides containing threonine, asparagine, alanine, phenylalanine, and leucine in position 12 into realistic model membranes showed these structures were stable for 20 ns of molecular dynamics simulation time. Threonine could participate in intermolecular hydrogen bonds, but predominantly formed hydrogen bonds to the backbone of the helix it resided on. These hydrogen bonds, although infrequent, appeared to promote closer association of polyleucine helices. Asparagine participated in multiple, rapidly fluctuating intermolecular and intramolecular hydrogen bonds, and may have slightly destabilized optimum van der Waals packing in favor of optimum hydrogen bonding. Coordinated rotations of transmembrane helices about their axes were observed, indicating helices may rotate around one another during the folding of membrane proteins or other processes. These rotations were inhibited by phenylalanine, suggesting a role for bulky residues in modulating membrane protein dynamics.     

Complex Patterns of Histidine, Hydroxylated Amino Acids and the GxxxG Motif Mediate High-affinity Transmembrane Domain Interactions

Specific interactions of transmembrane helices play a pivotal role in the folding and oligomerization of integral membrane proteins. The helix-helix interfaces frequently depend on specific amino acid patterns. In this study, a heptad repeat pattern was randomized with all naturally occurring amino acids to uncover novel sequence motifs promoting transmembrane domain interactions. Self-interacting transmembrane domains were selected from the resulting combinatorial library by means of the ToxR/POSSYCCAT system. A comparison of the amino acid composition of high-and low-affinity sequences revealed that high-affinity transmembrane domains exhibit position-specific enrichment of histidine. Further, sequences containing His preferentially display Gly, Ser, and/or Thr residues at flanking positions and frequently contain a C-terminal GxxxG motif. Mutational analysis of selected sequences confirmed the importance of these residues in homotypic interaction. Probing heterotypic interaction indicated that His interacts in trans with hydroxylated residues. Reconstruction of minimal interaction motifs within the context of an oligo-Leu sequence confirmed that His is part of a hydrogen bonded cluster that is brought into register by the GxxxG motif. Notably, a similar motif contributes to self-interaction of the BNIP3 transmembrane domain.     

The GxxxG motif: a framework for transmembrane helix-helix association

In order to identify strong transmembrane helix packing motifs, we have selected transmembrane domains exhibiting high-affinity homo-oligomerization from a randomized sequence library based on the right-handed dimerization motif of glycophorin A. Sequences were isolated using the TOXCAT system, which measures transmembrane helix-helix association in the Escherichia coli inner membrane. Strong selection was applied to a large range of sequences ( approximately 10(7) possibilities) and resulted in the identification of sequence patterns that mediate high-affinity helix-helix association. The most frequent motif isolated, GxxxG, occurs in over 80% of the isolates. Additional correlations suggest that flanking residues act in concert with the GxxxG motif, and that size complementarity is maintained at the interface, consistent with the idea that the identified sequence patterns represent packing motifs. The convergent identification of similar sequence patterns from an analysis of the transmembrane domains in the SwissProt sequence database suggests that these packing motifs are frequently utilized in naturally occurring helical membrane proteins.     

Hydroxylation-induced stabilization of the collagen triple helix. Further characterization of peptides with 4(R)-hydroxyproline in the Xaa position

4(R)-Hydroxyproline in the Yaa position of the -Gly-Xaa-Yaa-repeated sequence of collagen plays a crucial role in the stability of the triple helix. Since the peptide (4(R)-Hyp-Pro-Gly)10 does not form a triple helix, it was generally believed that polypeptides with a -Gly-4(R)-Hyp-Yaa-repeated sequence do not form a triple helix. Recently, we found that acetyl-(Gly-4(R)-Hyp-Thr)10-NH2 forms a triple helix in aqueous solutions. To further study the role of 4(R)-hydroxyproline in the Xaa position, we made a series of acetyl-(Gly-4(R)-Hyp-Yaa)10-NH2 peptides where Yaa was alanine, serine, valine, and allo-threonine. We previously hypothesized that the hydroxyl group of threonine might form a hydrogen bond to the hydroxyl group of 4(R)hydroxyproline. In water, only the threonine- and the valine-containing peptides were triple helical. The remaining peptides did not form a triple helix in water. In 1,2- and in 1,3-propanediol at 4 degrees C, all the soluble peptides were triple helical. From the transition temperature of the triple helices, it was found that among the examined residues, threonine was the most stable residue in the acetyl-(Gly-4(R)-Hyp-Yaa)10-NH2 peptide. The transition temperatures of the valine- and allo-threonine-containing peptides were 10 degrees lower than those of the threonine peptide. Surprisingly, the serine-containing peptide was the least stable. These results indicate that the stability of these peptides depends on the presence of a methyl group as well as the hydroxyl group and that the stereo configuration of the two groups is essential for the stability. In the threonine peptide, we hypothesize that the methyl group shields the interchain hydrogen bond between the glycine and the Xaa residue from water and that the hydroxyl groups of threonine and 4(R)hydroxyproline can form direct or water-mediated hydrogen bonds.     

Dissection of helix capping in T4 lysozyme by structural and thermodynamic analysis of six amino acid substitutions at Thr 59.

Threonine 59, a helix-capping residue at the amino terminus of the longest helix in T4 phage lysozyme, was substituted with valine, alanine, glycine, serine, asparagine, and aspartic acid. The valine, alanine, and glycine replacements were observed to be somewhat more destabilizing than serine, asparagine, and aspartic acid. The crystal structures of the different variants showed that changes in conformation occurred at the site of substitution, including Asp 61, which is nearby, as well as displacement of a solvent molecule that is hydrogen-bonded to the gamma-oxygen of Thr 59 in wild-type lysozyme. Neither the structures nor the stabilities of the mutant proteins support the hypothesis of Serrano and Fersht (1989) that glycine and alanine are better helix-capping residues than valine because a smaller-sized residue allows better hydration at the end of the helix. In the aspartic acid and asparagine replacements the substituted side chains form hydrogen bonds with the end of the helix, as does threonine and serine at this position. In contrast, however, the Asp and Asn side chains also make unusually close contacts with carbon atoms in Asp 61. This suggests a structural basis for the heretofore puzzling observations that asparagine is more frequently observed as a helix-capping residue than threonine [Richardson, J. S., & Richardson, D. C. (1988) Science 240, 1648-1652] yet Thr----Asn replacements at N-cap positions in barnase were found to be destabilizing [Serrano, L., & Fersht, A. R. (1989) Nature 342, 296-299].(ABSTRACT TRUNCATED AT 250 WORDS)     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Cloning and sequencing of a human pancreatic tumor mucin cDNA

A monospecific polyclonal antiserum against deglycosylated human pancreatic tumor mucin was used to select human pancreatic mucin cDNA clones from a lambda gt11 cDNA expression library developed from a human pancreatic tumor cell line. The full-length 4.4-kilobase mucin cDNA sequence included a 72-base pair 5'-untranslated region and a 307-base pair 3'-untranslated region. The predicted amino acid sequence for this cDNA revealed a protein of 122,071 daltons containing 1,255 amino acid residues of which greater than 60% were serine, threonine, proline, alanine, and glycine. Approximately two-thirds of the protein sequence consisted of identical 20-amino acid tandem repeats which were flanked by degenerate tandem repeats and nontandem repeat sequences on both the amino-terminal and carboxyl-terminal ends. The amino acid sequence also contained five putative N-linked glycosylation sites, a putative signal sequence and transmembrane domain, and numerous serine and threonine residues (potential O-linked glycosylation sites) outside and within the tandem repeat position. The cDNA and deduced amino acid sequence of the pancreatic mucin sequence was over 99% homologous with a mucin cDNA sequence derived from breast tumor mucin, even though the native forms of these molecules are quite distinct in size and degree of glycosylation.     

Accessible surfaces of beta proteins increase with increasing protein molecular mass more rapidly than those of other proteins

Here we present a systematic analysis of accessible surface areas and hydrogen bonds of 2554 globular proteins from four structural classes (all-α, all-β, α/β and α+β proteins) that is aimed to learn in which structural class the accessible surface area increases with increasing protein molecular mass more rapidly than in other classes, and what structural peculiarities are responsible for this effect. The beta structural class of proteins was found to be the leader, with the following possible explanations of this fact. First, in beta structural proteins, the fraction of residues not included in the regular secondary structure is the largest, and second, the accessible surface area of packaged elements of the beta-structure increases more rapidly with increasing molecular mass in comparison with the alpha-structure. Moreover, in the beta structure, the probability of formation of backbone hydrogen bonds is higher than that in the alpha helix for all residues of α+β proteins (the average probability is 0.73±0.01 for the beta-structure and 0.60±0.01 for the alpha-structure without proline) and α/β proteins, except for asparagine, aspartic acid, glycine, threonine, and serine (0.70±0.01 for the beta-structure and 0.60±0.01 for the alpha-structure without the proline residue). There is a linear relationship between the number of hydrogen bonds and the number of amino acid residues in the protein (Number of hydrogen bonds=0.678·number of residues-3.350).     

Sequence preference of α-helix N-terminal tetrapeptide

The α-helix is the most abundant secondary structure in proteins. Due to the specific i, i+4 hydrogen bond pattern, the two termini have unsatisfied hydrogen bonds, and are less constrained; in order to compensate for this, specific residues are preferred for the terminal positions. However, a naive combination of the statistically-preferred residues for each position may not result in a stable N-terminal helical sequence. In order to provide a set of preferable N-terminal peptides for α-helix design, we have studied the N-terminal tetrapeptide sequence motifs that are favorable for helix formation using statistical analysis and atomistic simulations. A set of tetrapeptide sequences including TEEE and TPEE were found to be favorable motifs. In addition to forming more hydrogen bonds in the helical conformation, the favorable motifs also tended to form more capping boxes. To empirically test our predictions, we obtained 10 peptides with different N-terminal motifs and measured their α-helical content by circular dichroism spectroscopy. The experimental results agreed qualitatively with the statistical and simulation results. Furthermore, some of the suggested preferable tetrapeptide sequences have been successfully applied in de novo protein design.     

Ionic Interactions Promote Transmembrane Helix-Helix Association Depending on Sequence Context

Folding and oligomerization of integral membrane proteins frequently depend on specific interactions of transmembrane helices. Interacting amino acids of helix-helix interfaces may form complex motifs and exert different types of molecular forces. Here, a set of strongly self-interacting transmembrane domains (TMDs), as isolated from a combinatorial library, was found to contain basic and acidic residues, in combination with polar nonionizable amino acids and C-terminal GxxxG motifs. Mutational analyses of selected sequences and reconstruction of high-affinity interfaces confirmed the cooperation of these residues in homotypic interactions. Probing heterotypic interaction indicated the presence of interhelical charge-charge interactions. Furthermore, simple motifs of an ionizable residue and GxxxG are significantly overrepresented in natural TMDs, and a specific combination of these motifs exhibits high-affinity heterotypic interaction. We conclude that intramembrane charge-charge interactions depend on sequence context. Moreover, they appear important for homotypic and heterotypic interactions of numerous natural TMDs.     

Phosphoserine/threonine-binding domains

Phosphorylation of proteins on serine and threonine residues has traditionally been viewed as a means to allosterically regulate catalytic activity. Research within the past five years, however, has revealed that serine/threonine phosphorylation can also directly result in the formation of multimolecular signaling complexes through specific interactions between phosphoserine/threonine (pSer/Thr)-binding modules and phosphorylated sequence motifs. pSer/Thr-binding proteins and domains currently include 14-3-3, WW domains, forkhead-associated domains, and, tentatively, WD40 repeats and leucine-rich regions. It seems likely that additional modules will be found in the future. The amino acid sequences recognized by these pSer/Thr-binding modules show partial overlap with the optimal phosphorylation motifs for different protein kinase subfamilies, allowing the formation of specific signaling complexes to be controlled through combinatorial interactions between particular upstream kinases and a particular binding module. The structural basis for pSer/Thr binding differs dramatically between 14-3-3 proteins, WW domains and forkhead-associated domains, suggesting that their pSer/Thr binding function was acquired through convergent evolution.     

Side-Chain to Main-Chain Hydrogen Bonding Controls the Intrinsic Backbone Dynamics of the Amyloid Precursor Protein Transmembrane Helix

Many transmembrane helices contain serine and/or threonine residues whose side chains form intrahelical H-bonds with upstream carbonyl oxygens. Here, we investigated the impact of threonine side-chain/main-chain backbonding on the backbone dynamics of the amyloid precursor protein transmembrane helix. This helix consists of a N-terminal dimerization region and a C-terminal cleavage region, which is processed by γ-secretase to a series of products. Threonine mutations within this transmembrane helix are known to alter the cleavage pattern, which can lead to early-onset Alzheimer’s disease. Circular dichroism spectroscopy and amide exchange experiments of synthetic transmembrane domain peptides reveal that mutating threonine enhances the flexibility of this helix. Molecular dynamics simulations show that the mutations reduce intrahelical amide H-bonding and H-bond lifetimes. In addition, the removal of side-chain/main-chain backbonding distorts the helix, which alters bending and rotation at a diglycine hinge connecting the dimerization and cleavage regions. We propose that the backbone dynamics of the substrate profoundly affects the way by which the substrate is presented to the catalytic site within the enzyme. Changing this conformational flexibility may thus change the pattern of proteolytic processing.     

Phenylalanine promotes interaction of transmembrane domains via GxxxG motifs

Interactions of transmembrane helices play a crucial role in the folding and oligomerisation of integral membrane proteins. In order to uncover novel sequence motifs mediating these interactions, we randomised one face of a transmembrane helix with a set of non-polar or moderately polar amino acids. Those sequences capable of self-interaction upon integration into bacterial inner membranes were selected by means of the ToxR/POSSYCCAT system. A comparison between low/medium-affinity and high-affinity sequences reveals that high-affinity sequences are strongly enriched in phenylalanine residues that are frequently observed at the − 3 position of GxxxG motifs, thus yielding FxxGxxxG motifs. Mutation of Phe or GxxxG in selected sequences significantly reduces self-interaction of the transmembrane domains without affecting their efficiency of membrane integration. Conversely, grafting FxxGxxxG onto unrelated transmembrane domains strongly enhances their interaction. Further, we find that FxxGxxxG is significantly over-represented in transmembrane domains of bitopic membrane proteins. The same motif contributes to self-interaction of the vesicular stomatitis virus G protein transmembrane domain. We conclude that Phe stabilises membrane-spanning GxxxG motifs. This is one example of how the role of certain side-chains in helix-helix interfaces is modulated by sequence context.