Comprehensive analysis of the numbers,lengths and amino acid compositions of transmembrane helices in prokaryotic,eukaryotic and viral integral membrane proteins of high-resolution structure

We report a comprehensive analysis of the numbers, lengths and amino acid compositions of transmembrane helices in 235 high-resolution structures of integral membrane proteins. The properties of 1551 transmembrane helices in the structures were compared with those obtained by analysis of the same amino acid sequences using topology prediction tools. Explanations for the 81 (5.2%) missing or additional transmembrane helices in the prediction results were identified. Main reasons for missing transmembrane helices were mis-identification of N-terminal signal peptides, breaks in α-helix conformation or charged residues in the middle of transmembrane helices and transmembrane helices with unusual amino acid composition. The main reason for additional transmembrane helices was mis-identification of amphipathic helices, extramembrane helices or hairpin re-entrant loops. Transmembrane helix length had an overall median of 24 residues and an average of 24.9 ± 7.0 residues and the most common length was 23 residues. The overall content of residues in transmembrane helices as a percentage of the full proteins had a median of 56.8% and an average of 55.7 ± 16.0%. Amino acid composition was analysed for the full proteins, transmembrane helices and extramembrane regions. Individual proteins or types of proteins with transmembrane helices containing extremes in contents of individual amino acids or combinations of amino acids with similar physicochemical properties were identified and linked to structure and/or function. In addition to overall median and average values, all results were analysed for proteins originating from different types of organism (prokaryotic, eukaryotic, viral) and for subgroups of receptors, channels, transporters and others.     

Amino acid screening based on structural modeling identifies critical residues for the function, ion selectivity and structure of Arabidopsis MTP1

Arabidopsis thaliana MTP1 is a vacuolar membrane Zn(2+)/H(+) antiporter of the cation diffusion facilitator family. Here we present a structure-function analysis of AtMTP1-mediated transport and its remarkable Zn(2+) selectivity by functional complementation tests of more than 50 mutant variants in metal-sensitive yeast strains. This was combined with homology modeling of AtMTP1 based on the crystal structure of the Escherichia coli broad-specificity divalent cation transporter YiiP. The Zn(2+)-binding sites of EcYiiP in the cytoplasmic C-terminus, and the pore formed by transmembrane helices TM2 and TM5, are conserved in AtMTP1. Although absent in EcYiiP, Cys31 and Cys36 in the extended N-terminal cytosolic domain of AtMTP1 are necessary for complementation of a Zn-sensitive yeast strain. On the cytosolic side of the active Zn(2+)-binding site inside the transmembrane pore, Ala substitution of either Asn258 in TM5 or Ser101 in TM2 non-selectively enhanced the metal tolerance conferred by AtMTP1. Modeling predicts that these residues obstruct the movement of cytosolic Zn(2+) into the intra-membrane Zn(2+)-binding site of AtMTP1. A conformational change in the immediately preceding His-rich cytosolic loop may displace Asn258 and permit Zn(2+) entry into the pore. This would allow dynamic coupling of Zn(2+) transport to the His-rich loop, thus acting as selectivity filter or sensor of cytoplasmic Zn(2+) levels. Individual mutations at diverse sites within AtMTP1 conferred Co and Cd tolerance in yeast, and included deletions in N-terminal and His-rich intra-molecular cytosolic domains, and mutations of single residues flanking the transmembrane pore or participating in intra- or inter-molecular domain interactions, all of which are not conserved in the non-selective EcYiiP.     

Deuterium/hydrogen exchange factors measured by solution nuclear magnetic resonance spectroscopy as indicators of the structure and topology of membrane proteins

Deuterium/hydrogen exchange factors (chi) were measured for the backbone amide sites of the membrane-bound forms of the 50-residue fd coat protein and the 23-residue magainin2 peptide in lipid micelles by solution nuclear magnetic resonance spectroscopy. By combining kinetic and thermodynamic effects, deuterium/hydrogen exchange factors overcome the principal limitations encountered in the measurements of kinetic protection factors and thermodynamic fractionation factors for membrane proteins. The magnitudes of the exchange factors can be correlated with the structure and topology of membrane-associated polypeptides. In fd coat protein, residues in the transmembrane helix have exchange factors that are substantially smaller than those in the amphipathic surface helix or the loop connecting the two helices. For the amphipathic helical peptide, magainin2, the exchange factors of residues exposed to the solvent are appreciably larger than those that face the hydrocarbon portion of membrane bilayers. These examples demonstrate that deuterium/hydrogen exchange factors can be measured by solution NMR spectroscopy and used to identify residues in transmembrane helices as well as to determine the polarity of amphipathic helices in membrane proteins.     

Conserved positioning of proline residues in membrane-spanning helices of ion-channel proteins.

Proline residues are a common feature of known and putative transmembrane helices of transport proteins. We find considerable consistency in the positioning of these residues within the structures. The proline residues are usually found on the hydrophilic (interior) faces of the pore-forming helices. This general observation adds considerable support to hypotheses concerning the structure of the ion-channels formed by alamethicin and melittin. As proline kinks helices, our observation suggests that the pores formed in ion-channel proteins tend to be funnel-shaped having a constriction near their center. Such a structure can aid in the capture of ions by the channel (an entropic effect) and should help in the gating mechanism of the channel. The observation will aid identification of putative transmembrane helices of ion-channels.     

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.     

The Escherichia coli twin-arginine translocase: conserved residues of TatA and TatB family components involved in protein transport

The Escherichia coli Tat system serves to export folded proteins harbouring an N-terminal twin-arginine signal peptide across the cytoplasmic membrane. In this report we have studied the functions of conserved residues within the structurally related TatA and TatB proteins. Our results demonstrate that there are two regions within each protein of high sequence conservation that are critical for efficient Tat translocase function. The first region is the interdomain hinge between the transmembrane and the amphipathic alpha-helices of TatA and TatB proteins. The second region is within the amphipathic helices of TatA and TatB. In particular an invariant phenylalanine residue within TatA proteins is essential for activity, whereas a string of glutamic acid residues on the same face of the amphipathic helix of TatB is important for function.     

Novel inter- and intrasubunit contacts between transport-relevant residues of the homodimeric mitochondrial phosphate transport protein

Ser158 is located near the middle of the matrix loop connecting transmembrane helices C and D of the mitochondrial phosphate transport protein (PTP). The mutant Ser158Thr PTP is transport-inactive. His32 is located near the middle of transmembrane helix A and Thr79 is located 5 residues away from transmembrane helix B and its N-terminal (matrix end). Single site mutant PTPs that have either residue replaced with Ala are transport-inactive. Based on the high resolution structure of a subunit of the bovine ADP/ATP translocase, on sequence similarities between members of the mitochondrial transport protein family, and on the PTP subunit/subunit contact site between transmembrane A helices, it is now suggested that the Ser158 site is at the PTP subunit/subunit contact site. This contact site is essential for keeping the transport cycles catalyzed by the two PTP subunits 180 degrees out of phase. The data also suggest that His32 and Thr79 of the same subunit interact and couple the phosphate and the proton transport paths.     

Single replacement constructs of all hydroxyl, basic, and acidic amino acids identify new function and structure-sensitive regions of the mitochondrial phosphate transport protein

The phosphate transport protein (PTP) catalyzes the proton cotransport of phosphate into the mitochondrial matrix. It functions as a homodimer, and thus residues of the phosphate and proton pores are somewhat scattered throughout the primary sequence. With 71 new single mutation per subunit PTPs, all its hydroxyl, basic, and acidic residues have now been replaced to identify these essential residues. We assayed the initial rate of pH gradient-dependent unidirectional phosphate transport activity and the liposome incorporation efficiency (LIE) of these mutants. Single mutations of Thr79, Tyr83, Lys90, Tyr94, and Lys98 inactivate transport. The spacings between these residues imply that they are located along the same face of transmembrane (TM) helix B, requiring an extension of its current model C-terminal domain by 10 residues. This extension superimposes very well onto the shorter bovine PTP helix B, leaving a 15-residue hydrophobic extension of the yeast helix B N-terminus. This is similar to the helix D and F regions of the yeast PTP. Only one transport-inhibiting mutation is located within loops: Ser158Thr in the matrix loop between helices C and D. All other transport-inhibiting mutations are located within the TM helices. Mutations that yield LIEs of <6% are all, except for four, within helices. The four exceptions are Tyr12Ala near the PTP N-terminus and Arg159Ala, Glu163Gln, and Glu164Gln in the loop between helices C and D. The PTP C-terminal segment beyond Thr214 at the N-terminus of helix E has 11 mutations with LIEs >20% and none with LIE <6%. Mutations with LIEs >20% are located near the ends of all the TM helices except TM helix D. Only a few mutations alter PTP structure (LIE) and also affect PTP transport activity. A novel observation is that Ser4Ala blocks the formation of PTP bacterial inclusion bodies.     

Characterization of genes encoding metal tolerance proteins isolated from Nicotiana glauca and Nicotiana tabacum

We have isolated a metal tolerance protein (MTP) gene, NgMTP1, from Nicotiana glauca (a potential phytoremediator plant) and two MTP genes, NtMTP1a and NtMTP1b, from Nicotiana tabacum. These three genes shared approximately 95% homology at the amino acid level. Heterologous expression of any of these three genes complemented Zn and Co tolerance in yeast mutants to a similar extent. In yeast, these proteins were shown to be located to vacuole membrane. These results suggest that the three MTPs operate by sequestering Zn and Co into vacuoles, thereby reducing the toxicity of these metals.     

Tertiary contacts of helix V in the lactose permease determined by site-directed chemical cross-linking in situ

The six N-terminal transmembrane helices (N6) and the six C-terminal transmembrane helices (C6) in lactose permease, each containing a single Cys residue, were coexpressed, and cross-linking was studied. The proximity of paired Cys residues in helices V and VII, VIII, or X was studied by thiol-specific chemical cross-linking. The results demonstrate that Cys residues in the periplasmic half of helix V cross-link with Cys residues in the periplasmic half of helix VII. In contrast, no cross-linking is evident with paired Cys residues in the cytoplasmic halves of helices V and VII. Moreover, Cys residues on one entire face of helix V cross-link with Cys residues on one face of helix VIII. Finally, paired Cys residues at the cytoplasmic ends of helices V and X cross-link, but no cross-linking is observed when paired Cys residues are placed at the periplasmic ends of the two helices. Taken together, the results indicate that the periplasmic halves of helices V and VII are in close proximity and that the two helices tilt away from one another toward the cytoplasmic side of the membrane. Furthermore, helices V and VIII are in close proximity throughout their lengths and do not tilt appreciably with respect to one another, and helices V and X are in close proximity at the cytoplasmic but not at the periplasmic face of the membrane.     

Mitochondrial phosphate transport protein. Reversions of inhibitory conservative mutations identify four helices and a nonhelix protein segment with transmembrane interactions and Asp39, Glu137, and Ser158 as nonessential for transport

The mitochondrial phosphate transport protein (PTP) has six (A--F) transmembrane (TM) helices per subunit of functional homodimer with all mutations referring to the subunit of the homodimer. In earlier studies, conservative replacements of several residues located either at the matrix end (Asp39/helix A, Glu137/helix C, Asp236/helix E) or at the membrane center (His32/helix A, Glu136/helix C) of TM helices yielded inactive single mutation PTPs. Some of these residues were suggested to act as phosphate ligands or as part of the proton cotransport path. We now show that the mutation Ser158Thr, not part of a TM helix but located near the center of the matrix loop (Ile141--Ser171) between TM helices C and D, inactivates PTP and is thus also functionally relevant. On the other side of the membrane, the single mutation Glu192Asp at the intermembrane space end of TM helix D yields a PTP with 33% wild-type activity. We constructed double mutants by adding this mutation to the six transport-inactivating mutations. Transport was detected only in those with Asp39Asn, Glu137Gln, or Ser158Thr. We conclude that TM helix D can interact with TM helices A and C and matrix loop Ile141--Ser171 and that Asp39, Glu137, and Ser158 are not essential for phosphate transport. Since our results are consistent with residues present in all 12 functionally identified members of the mitochondrial transport protein (MTP) family, they lead to a general rule that specifies MTP residue types at 7 separate locations. The conformations of all the double mutation PTPs (except that with the matrix loop Ser158Thr) are significantly different from those of the single mutation PTPs, as indicated by their very low liposome incorporation efficiency and their requirement for less detergent (Triton X-100) to stay in solution. These dramatic conformational differences also suggest an interaction between TM helices D and E. The results are discussed in terms of TM helix movements and changes in the PTP monomer/dimer ratio.     

Proline residues link the active site to transmembrane domain movements in human nucleoside triphosphate diphosphohydrolase 3 (NTPDase3)

The active sites of the membrane-bound nucleoside triphosphate diphosphohydrolases (NTPDases) regulate and are regulated by coordinated and spatially distant movements of their transmembrane helices, modulating enzyme activity, and substrate specificity. Using site-directed mutagenesis, the roles of the conserved proline residues (N-terminal: P52 and P53; C-terminal: P472, P476, P481, P484, and P485) of human NTPDase3, located in the “linker regions” that connect the N- and C-terminal transmembrane helices with the extracellular active site, were examined. Single cysteine substitutions were strategically placed in the transmembrane domain (N-terminal helix: V42C; C-terminal helix: G489C) to serve as cross-linking “sensors” of helical interactions. These “sensor” background mutant proteins (V42C and G489C NTPDase3) are enzymatically active and are cross-linked by copper phenanthroline less efficiently in the presence of adenosine triphosphate (ATP). Proline to alanine substitutions at P53, P481, P484, and P485 in the V42C background, as well as P53, P481, and P484 in the G489C background, exhibited decreased nucleotidase activities. More importantly, alanine substitutions at P53 and P481 in the V42C background and P481 in the G489C background no longer exhibited the ATP-induced decrease in transmembrane cross-linking efficiency. Interestingly, the P485A mutation abolished oxidative cross-linking at G489C both in the presence and absence of ATP. Taken together, these results suggest a role for proline residues 53 and 481 in the linker regions of human NTPDase3 for coupling nucleotide binding at the enzyme active site to movements and/or rearrangements of the transmembrane helices necessary for optimal nucleotide hydrolysis.     

Proline residues in transmembrane helices of channel and transport proteins: a molecular modelling study.

Proline residues are commonly found in putative transbilayer helices of many integral membrane proteins which act as transporters, channels and receptors. Intramembranous prolines are often conserved between homologous proteins. It has been suggested that such intrahelical prolines provide liganding sites for cations via exposure of the backbone carbonyl oxygen atoms of residues i-3 and i-4 (relative to the proline). Molecular modelling studies have been carried out to evaluate this proposal. Bundles of parallel proline-kinked helices are considered as simplified models of ion channels. The energetics of K+ ion-helix bundle interactions are explored. It is shown that carbonyl oxygens exposed by the proline-induced kink and at the C-terminus of the helices may provide cation-liganding sites. 'Hybrid' bundles of antiparallel helices, only some of which contain proline residues, are considered as models of transport proteins. Again, proline-exposed carbonyl oxygens are shown to be capable of liganding cations. The roles of alpha-helix dipoles and of the geometry of helix packing are considered in relation to cation-bundle interactions. Implications with respect to modelling of ion channel and transport proteins are discussed.     

Revealing the structure of the photosystem II chlorophyll binding proteins, CP43 and CP47

A review of the structural properties of the photosystem II chlorophyll binding proteins, CP47 and CP43, is given and a model of the transmembrane helical domains of CP47 has been constructed. The model is based on (i) the amino acid sequence of the spinach protein, (ii) an 8 A three-dimensional electron density map derived from electron crystallography and (iii) the structural homology which the membrane spanning region of CP47 shares with the six N-terminal transmembrane helices of the PsaA/PsaB proteins of photosystem I. Particular emphasis has been placed on the position of chlorophyll molecules assigned in the 8 A three-dimensional map of CP47 (K.-H. Rhee, E.P. Morris, J. Barber, W. Kühlbrandt, Nature 396 (1998) 283-286) relative to histidine residues located in the transmembrane regions of this protein which are likely to form axial ligands for chlorophyll binding. Of the 14 densities assigned to chlorophyll, the model predicted that five have their magnesium ions within 4 A of the imidazole nitrogens of histidine residues. For the remaining seven histidine residues the densities attributed to chlorophylls were within 4-8 A of the imidazole nitrogens and thus too far apart for direct ligation with the magnesium ion within the tetrapyrrole head group. Improved structural resolution and reconsiderations of the orientation of the porphyrin rings will allow further refinement of the model.     

Computational differentiation of N-terminal signal peptides and transmembrane helices

Signal peptides and transmembrane helices both contain a stretch of hydrophobic amino acids. This common feature makes it difficult for signal peptide and transmembrane helix predictors to correctly assign identity to stretches of hydrophobic residues near the N-terminal methionine of a protein sequence. The inability to reliably distinguish between N-terminal transmembrane helix and signal peptide is an error with serious consequences for the prediction of protein secretory status or transmembrane topology. In this study, we report a new method for differentiating protein N-terminal signal peptides and transmembrane helices. Based on the sequence features extracted from hydrophobic regions (amino acid frequency, hydrophobicity, and the start position), we set up discriminant functions and examined them on non-redundant datasets with jackknife tests. This method can incorporate other signal peptide prediction methods and achieve higher prediction accuracy. For Gram-negative bacterial proteins, 95.7% of N-terminal signal peptides and transmembrane helices can be correctly predicted (coefficient 0.90). Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 99% (coefficient 0.92). For eukaryotic proteins, 94.2% of N-terminal signal peptides and transmembrane helices can be correctly predicted with coefficient 0.83. Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 87% (coefficient 0.85). The method can be used to complement current transmembrane protein prediction and signal peptide prediction methods to improve their prediction accuracies.     

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.     

Proline-induced distortions of transmembrane helices

Proline residues in the transmembrane (TM) alpha-helices of integral membrane proteins have long been suspected to play a key role for helix packing and signal transduction by inducing regions of helix distortion and/or dynamic flexibility (hinges). In this study we try to characterise the effect of proline on the geometric properties of TM alpha-helices. We have examined 199 transmembrane alpha-helices from polytopic membrane proteins of known structure. After examining the location of proline residues within the amino acid sequences of TM helices, we estimated the helix axes either side of a hinge and hence identified a hinge residue. This enabled us to calculate helix kink and swivel angles. The results of this analysis show that proline residues occur with a significant concentration in the centre of sequences of TM alpha-helices. In this location, they may induce formation of molecular hinges, located on average about four residues N-terminal to the proline residue. A superposition of proline-containing TM helices structures shows that the distortion induced is anisotropic and favours certain relative orientations (defined by helix kink and swivel angles) of the two helix segments.     

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.     

Co-evolving residues in membrane proteins

MOTIVATION: The analysis of co-evolving residues has been exhaustively evaluated for the prediction of intramolecular amino acid contacts in soluble proteins. Although a variety of different methods for the detection of these co-evolving residues have been developed, the fraction of correctly predicted contacts remained insufficient for their reliable application in the construction of structural models. Membrane proteins, which constitute between one-fourth and one-third of all proteins in an organism, were only considered in few individual case studies. RESULTS: We present the first general study of correlated mutations in alpha-helical membrane proteins. Using seven different prediction algorithms, we extracted co-evolving residues for 14 membrane proteins having a solved 3D structure. On average, distances between correlated pairs of residues lying on different transmembrane segments were found to be significantly smaller compared to a random prediction. Covariation of residues was frequently found in direct sequence neighborhood to helix-helix contacts. Based on the results obtained from individual prediction methods, we constructed a consensus prediction for every protein in the dataset that combines obtained correlations from different prediction algorithms and simultaneously removes likely false positives. Using this consensus prediction, 53% of all predicted residue pairs were found within one helix turn of an observed helix-helix contact. Based on the combination of co-evolving residues detected with the four best prediction algorithms, interacting helices could be predicted with a specificity of 83% and sensitivity of 42%. AVAILABILITY: http://webclu.bio.wzw.tum.de/helixcorr/     

FtsH recognizes proteins with unfolded structure and hydrolyzes the carboxyl side of hydrophobic residues

FtsH of Escherichia coli is an essential membrane-integrated ATP-dependent protease. We cloned a gene for an FtsH homolog (T. FtsH) from Thermus thermophilus HB8, expressed it in E. coli, and purified the expressed protein. ATPase activity of T.FtsH was activated by proteins with unfolded structure ( alpha-casein and pepsin), and T.FtsH digested these proteins in an ATP-, Zn(2+)-dependent manner. alpha-Lactalbumin was digested by T.FtsH when it was largely unfolded, but not in its native form. Analysis of the proteolytic products revealed that, in most cases, T.FtsH cleaved the C-terminal side of hydrophobic residues and produced a characteristic set of small peptides (<30 kDa) without releasing a large intermediate. Thus, T.FtsH recognizes the unfolded structure of the proteins and progressively digests them at the expense of ATP. A soluble domain of T.FtsH, which lacked the N-terminal two transmembrane helices, was also prepared but was found to retain neither ATPase nor protease activities. Thus, the membrane segment appeared to be indispensable for these activities of T.FtsH.