Chemical analysis of processing of spiralin, the major lipoprotein of Spiroplasma melliferum

The plasma membrane of Spiroplasma melliferum contains a major membrane-associated lipoprotein called spiralin. In this study, the processing pathway of spiralin was investigated by chemical analysis of the purified protein and by using [35S]cysteine, [35S]methionine, [14C]myristic acid (14C-14:0), [14C]palmitic acid (14C-16:0), and globomycin. SDS-PAGE analysis of membrane proteins showed the leader peptide cleavage of prospiralin and provided evidence for an apparent selectivity in the acylation: the unprocessed protein was labelled with 14C-16:0 only (O-ester-linked acyl chains), and the mature form with both 14C-labelled fatty acids (O-ester-linked + amide-linked chains). Chemical analysis of the purified protein revealed that spiralin contains S-glycerylcysteine and is covalently modified with two O-ester-linked acyl chains and one amide-linked fatty acid chain. However, a specific selectivity in the O- and the N-acylations was not confirmed; palmitate and stearate were the major components. The amounts of O-ester- and amide-linked acyl chains, the resistance to Edman degradation and the presence of S-glycerylcysteine together indicate that spiralin is a "classical" lipoprotein (i.e. is triacylated) and is probably processed by a mechanism similar to that described for gram-negative eubacteria. On the basis of these findings, a biogenesis pathway for spiralin is proposed.     

Selective acylation of membrane proteins in Acholeplasma laidlawii

In membranes of the cell-wall-less prokaryote Acholeplasma laidlawii most proteins are of the integral type. A substantial fraction of these proteins are enriched in hydrophilic amino acid residues. Approximately 20 different major as well as minor proteins were found to be covalently modified with acyl chains. The same set of proteins are acylated when cells are grown in different fatty-acid-supplemented media. In individual proteins the ratio of palmitoyl/oleoyl acyl chains was 12-14 times larger than the acyl chain ratio in polar membrane lipids. The transmembrane protein D12 has close to two acyl chains per molecule. Proteins T2 and T4a, localized in the outer and inner leaflet of the membrane, respectively, occur each as pairs with a difference in relative molecular mass within each pair of approximately 2000. Each of these proteins as well as the other acyl proteins, except the light form of T4a, has close to one acyl chain per molecule. The extent of acylation was increased for certain proteins and decreased for others by treatment with globomycin or phenethylalcohol. The relative amounts of the T2 and T4a pairs were affected by these drugs. It is concluded that the mechanism of acylation is different from that in Escherichia coli lipoprotein and Bacillus penicillinase. The mean hydrophobicity [Kyte & Doolittle (1982) J. Mol. Biol. 157, 105-132] of the A. laidlawii acyl proteins are similar to those of other bacterial acyl proteins but significantly lower than for non-acylated integral membrane proteins, supporting an anchoring function of the acyl chains. The number of membrane acyl proteins in A. laidlawii and two other mycoplasmas are at least twice that in other bacteria.     

Secondary structure of spiralin in solution, at the air/water interface, and in interaction with lipid monolayers

The surface of spiroplasmas, helically shaped pathogenic bacteria related to the mycoplasmas, is crowded with the membrane-anchored lipoprotein spiralin whose structure and function are unknown. In this work, the secondary structure of spiralin under the form of detergent-free micelles (average Stokes radius, 87.5 A) in water and at the air/water interface, alone or in interaction with lipid monolayers was analyzed. FT-IR and circular dichroism (CD) spectroscopic data indicate that spiralin in solution contains about 25+/-3% of helices and 38+/-2% of beta sheets. These measurements are consistent with a consensus predictive analysis of the protein sequence suggesting about 28% of helices, 32% of beta sheets and 40% of irregular structure. Brewster angle microscopy (BAM) revealed that, in water, the micelles slowly disaggregate to form a stable and homogeneous layer at the air/water interface, exhibiting a surface pressure up to 10 mN/m. Polarization modulation infrared reflection absorption spectroscopy (PMIRRAS) spectra of interfacial spiralin display a complex amide I band characteristic of a mixture of beta sheets and alpha helices, and an intense amide II band. Spectral simulations indicate a flat orientation for the beta sheets and a vertical orientation for the alpha helices with respect to the interface. The combination of tensiometric and PMIRRAS measurements show that, when spiroplasma lipids are used to form a monolayer at the air/water interface, spiralin is adsorbed under this monolayer and its antiparallel beta sheets are mainly parallel to the polar-head layer of the lipids without deep perturbation of the fatty acid chains organization. Based upon these results, we propose a 'carpet model' for spiralin organization at the spiroplasma cell surface. In this model, spiralin molecules anchored into the outer leaflet of the lipid bilayer by their N-terminal lipid moiety are composed of two colinear domains (instead of a single globular domain) situated at the lipid/water interface. Owing to the very high amount of spiralin in the membrane, such carpets would cover most if not all the lipids present in the outer leaflet of the bilayer.     

Anomalous orientation-dependent effective pair interaction among histidine and other amino acid residues in metalloproteins: breakdown of the hydropathy scale index

Hydropathy scale is widely used to obtain a measure of the effective interaction between any two amino acid residues in proteins and is based on the assumption that attraction between two hydrophobic groups and repulsion between hydrophilic groups (in water) can be translated straightforwardly to protein environment. Here we employ a recently developed statistical mechanical approach combined with the Protein Data Bank to obtain both distance- and orientation-dependent potential of mean force (ODPMF). This allowed us to explore effective pair potential among many amino acid residues and to examine the validity of the hydropathy scale in modeling the interaction among amino acid residues. We find that in some cases, like Phe-Phe and Lys-Lys, the hydropathy scale approach is largely obeyed. However, we also observe many unexpected pair interactions which defy the trend given by published hydropathy scales. An example of the former is the arginine-arginine (Arg-Arg) pair interaction which is found to be strongly and surprisingly attractive at short separation, even though it is the most hydrophilic residue. Here the head-to-head (see text) interaction is also stabilized. Tryptophan residues also exhibit strong attractive interaction. Equally important, we find strong influence of metal in determining effective interaction among the amino acid residues. It is the behavior of the histidine (His) which is found to be the most unusual. It exhibits a strong attractive interaction with itself which gets significantly enhanced in metalloproteins. These results highlight the important (sometime hidden) role of metals in protein structure and folding.     

Transmembrane topology of pmt1p, a member of an evolutionarily conserved family of protein O-mannosyltransferases

The identification of the evolutionarily conserved family of dolichyl-phosphate-D-mannose:protein O-mannosyltransferases (Pmts) revealed that protein O-mannosylation plays an essential role in a number of physiologically important processes. Strikingly, all members of the Pmt protein family share almost identical hydropathy profiles; a central hydrophilic domain is flanked by amino- and carboxyl-terminal sequences containing several putative transmembrane helices. This pattern is of particular interest because it diverges from structural models of all glycosyltransferases characterized so far. Here, we examine the transmembrane topology of Pmt1p, an integral membrane protein of the endoplasmic reticulum, from Saccharomyces cerevisiae. Structural predictions were directly tested by site-directed mutagenesis of endogenous N-glycosylation sites, by fusing a topology-sensitive monitor protein domain to carboxyl-terminal truncated versions of the Pmt1 protein and, in addition, by N-glycosylation scanning. Based on our results we propose a seven-transmembrane helical model for the yeast Pmt1p mannosyltransferase. The Pmt1p amino terminus faces the cytoplasm, whereas the carboxyl terminus faces the lumen of the endoplasmic reticulum. A large hydrophilic segment that is oriented toward the lumen of the endoplasmic reticulum is flanked by five amino-terminal and two carboxyl-terminal membrane spanning domains. We could demonstrate that this central loop is essential for the function of Pmt1p.     

Spiralins of Spiroplasma citri and Spiroplasma melliferum: amino acid sequences and putative organization in the cell membrane

Spiralin is the major membrane protein of the helical mollicute Spiroplasma citri. A similar protein occurs in the membrane of Spiroplasma melliferum, an organism related to S. citri. The gene encoding spiralin has been sequenced. A restriction fragment of the spiralin gene has been used as a probe to detect the gene encoding S. melliferum spiralin. A 4.6-kilobase-pair ClaI DNA fragment from S. melliferum strongly hybridized with the probe. This fragment was inserted in pBR322 and cloned in Escherichia coli. It was further subcloned in the replicative forms of M13mp18 and M13mp19, and its nucleotide sequence was determined (GenBank accession number M33991). An open reading frame showing 88.6% base sequence homology with the S. citri spiralin gene could be identified and was assumed to be the gene encoding S. melliferum spiralin. The deduced amino acid sequence of the protein had 75% homology with the spiralin sequence. In particular, the two proteins possess a stretch of 20 amino acids which can form an alpha-helix, in which all polar amino acids occupy approximately one-third of the axial projection down the helix. On the basis of these data and published data, we propose a topological model for the structural organization of the spiralin in the cell membrane of spiroplasmas.     

The role of leucine and isoleucine in tuning the hydropathy of class A GPCRs

Leucine and Isoleucine are two amino acids that differ only by the positioning of one methyl group. This small difference can have important consequences in α-helices, as the β-branching of Ile results in helix destabilization. We set out to investigate whether there are general trends for the occurrences of Leu and Ile residues in the structures and sequences of class A GPCRs (G protein-coupled receptors). GPCRs are integral membrane proteins in which α-helices span the plasma membrane seven times and which play a crucial role in signal transmission. We found that Leu side chains are generally more exposed at the protein surface than Ile side chains. We explored whether this difference might be attributed to different functions of the two amino acids and tested if Leu tunes the hydrophobicity of the transmembrane domain based on the Wimley-White whole-residue hydrophobicity scales. Leu content decreases the variation in hydropathy between receptors and correlates with the non-Leu receptor hydropathy. Both measures indicate that hydropathy is tuned by Leu. To test this idea further, we generated protein sequences with random amino acid compositions using a simple numerical model, in which hydropathy was tuned by adjusting the number of Leu residues. The model was able to replicate the observations made with class A GPCR sequences. We speculate that the hydropathy of transmembrane domains of class A GPCRs is tuned by Leu (and to some lesser degree by Lys and Val) to facilitate correct insertion into membranes and/or to stably anchor the receptors within membranes.     

Organization and nucleotide sequences of the Spiroplasma citri genes for ribosomal protein S2, elongation factor Ts, spiralin, phosphofructokinase, pyruvate kinase, and an unidentified protein

The gene for spiralin, the major membrane protein of the helical mollicute Spiroplasma citri, was cloned in Escherichia coli as a 5-kilobase-pair (kbp) DNA fragment. The complete nucleotide sequence of the 5.0-kbp spiroplasmal DNA fragment was determined (GenBank accession no. M31161). The spiralin gene was identified by the size and amino acid composition of its translational product. Besides the spiralin gene, the spiroplasmal DNA fragment was found to contain five additional open reading frames (ORFs). The translational products of four of these ORFs were identified by their amino acid sequence homologies with known proteins: ribosomal protein S2, elongation factor Ts, phosphofructokinase, and pyruvate kinase, respectively encoded by the genes rpsB, tsf, pfk, and pyk. The product of the fifth ORF remains to be identified and was named protein X (X gene). The order of the above genes was tsf--X--spiralin gene--pfk--pyk. These genes were transcribed in one direction, while the gene for ribosomal protein S2 (rpsB) was transcribed in the opposite direction.     

Spiralin polymorphism in strains of Spiroplasma citri is not due to differences in posttranslational palmitoylation.

Spiralin is defined as the major membrane protein of the helical mollicute Spiroplasma citri. According to the S. citri strain used, spiralin shows polymorphism in its electrophoretic mobility. The spiralin gene sequences of eight S. citri strains were determined by direct sequencing of the PCR-amplified genes. All spiralins were found to be 241 amino acids long, except for the spiralin of strain Palmyre, which is 242 amino acids long. The molecular masses calculated from these sequences did not explain the differences observed in the electrophoretic mobilities. In all of the spiralins examined, the first 24 N-terminal amino acids were conserved, including a cysteine at position 24, and had the features of typical signal peptides of procaryotic lipoproteins. When S. citri strains were grown in the presence of [3H]palmitic acid, at least 10 proteins, including spiralin, became labeled. In the presence of globomycin, a lipoprotein signal peptidase inhibitor in eubacteria, apparently unprocessed spiralin could be detected. Formic acid hydrolysis of the [3H]palmitic acid-labeled spiralins of four representative S. citri strains yielded two peptide fragments for each spiralin, as expected from the gene sequence. On fragment was [3H]palmitic acid labeled, and it had almost the same electrophoretic mobility irrespective of the spiralins used. Samples of the unlabeled peptide fragments from the four representative strains had slightly different electrophoretic mobilities (delta Da approximately equal to 800 Da); however, these were much smaller than those of the whole spiralins before formic acid hydrolysis (delta Da approximately equal to 8,000 Da). These results suggest that spiralin polymorphism in S. citri is not due to differences in posttranslational modification by palmitic acid and is certainly a structural property of the whole protein or could result from an unidentified posttranslational modification of spiralin.     

Hydrophobic membrane thickness and lipid-protein interactions of the leucine transport system of Lactococcus lactis

The effect of the phospholipid acyl chain carbon number on the activity of the branched-chain amino acid transport system of Lactococcus lactis has been investigated. Major fatty acids identified in a total lipid extract of L. lactis membranes are palmitic acid (16:0), oleic acid (18:1) and the cyclopropane-ring containing lactobacillic acid (19 delta). L. lactis membrane vesicles were fused with liposomes prepared from equimolar mixtures of synthetic phosphatidylethanolamine (PE) and phosphatidylcholine (PC) with cis mono-unsaturated acyl chains. The activity of the branched-chain amino acid carrier is determined by the bulk properties of the membrane (Driessen, A.J.M., Zheng, T., In 't Veld, G., Op den Kamp, J.A.F. and Konings, W.N. (1988) Biochemistry 27, 865-872). PE acts as an activator and PC is ineffective. Counterflow and protonmotive-force driven transport of leucine is sensitive to changes in the acyl chain carbon number of both phospholipids and maximal with dioleoyl-PE/dioleoyl-PC. Above the gel to liquid-crystalline phase transition temperature of the lipid species, membrane fluidity decreased with increasing acyl chain carbon number. Our data suggest that the carbon number of the acyl chains of PE and PC determine to a large extent the activity of the transport system. This might be relevant for the interaction of PE with the transport protein. Variations in the acyl chain composition of PC exert a more general effect on transport activity. The acyl chain composition of phospholipids determines the membrane thickness (Lewis, B.A. and Engelman, D.M. (1983) J. Mol. Biol. 166, 211-217). We therefore propose that the degree of matching between the lipid-bilayer and the hydrophobic thickness of the branched-chain amino acid carrier is an important parameter in lipid-protein interactions.     

Role of mutational bias and natural selection on genome-wide nucleotide bias in prokaryotic organisms

Correlations between genomic GC contents and amino acid frequencies were studied in the homologous sequences of 12 eubacterial genomes. Results show that amino acids encoded by GC-rich codons increases significantly with genomic GC contents, whereas opposite trend was observed in case of amino acids encoded by GC-poor codons. Further studies show all the amino acids do not change in the predicted direction according to their genomic GC pressure, suggesting that protein evolution is not entirely dictated by their nucleotide frequencies. Amino acid substitution matrix calculated among hydrophobic, amphipathic and hydrophilic amino acid groups' shows that amphipathic and hydrophilic amino acids are more frequently substituted by hydrophobic amino acids than from hydrophobic to hydrophilic or amphipathic amino acids. This indicates that nucleotide bias induces a directional changes in proteome composition in such a way that underwent strong changes in hydropathy values. In fact, significant increases in hydrophobicity values have also been observed with the increase of genomic GC contents. Correlations between GC contents and amino acid compositions in three different predicted protein secondary structures show that hydropathy values increases significantly with GC contents in aperiodic and helix structures whereas strand structure remains insensitive with the genomic GC levels. The relative importance of mutation and selection on the evolution of proteins have been discussed on the basis of these results.     

Purification, amino acid composition and N-terminal sequence of the major protein (protein H) of the outer membrane of Pasteurella multocida]

The major protein (protein H) of the outer membrane of Pasteurella multocida was purified by size-exclusion chromatography after selective extraction with detergents. The protein forms homotrimers which are stable in the presence of SDS at room temperature. Upon treatment at 100 degrees C, the protein is fully dissociated by the detergent into monomers exhibiting an apparent molecular mass of 37 kDa as estimated by electrophoresis. The amino acid composition of protein H is characterized by a low hydropathy index (HI = -0.40) and is strongly related to the compositions of bacterial porins, notably porins P2 (Haemophilus influenzae), PIA (Neisseria gonorrhoeae) and Cl.2 ("class 2 porin" of N. meningitidis). The N-terminal amino acid sequence of protein H shares a strong homology with those of porins OmpC (Escherichia coli) and P2. These data indicate that protein H of P. multocida is a porin belonging to the superfamily of the non-specific porins of Gram-negative eubacteria outer membrane.     

Amphilphilic nature of spiralin, the major protein of the Spiroplasma citri cell membrane.

Spiralin could not be solubilized in the absence of detergents, and it was shown by charge-shift crossed immunoelectrophoresis that this protein was capable of binding detergents under nondenaturing conditions. These properties indicate the amphiphilic nature of spiralin, which therefore should be regarded as an intrinsic membrane protein. The efficiency of mild (ionic and neutral) detergents to solubilize spiralin was as follows: deoxycholate greater than lauroyl sarcosinate, cholate, taurocholate, taurodeoxycholate greater than Triton X-100 greater than Brij 58 greater than Tween 20, indicating that mild ionic detergents were more effective than neutral ones. Solubilization of spiralin was quantitative with sodium deoxycholate. It was also shown that although a membrane protein is not extractable by a given detergent from the membrane, this does not necessarily mean that the protein is not soluble in this detergent.     

The bacterial porin superfamily: sequence alignment and structure prediction

The porins of Gram-negative bacteria are responsible for the 'molecular sieve' properties of the outer membrane. They form large water-filled channels which allow the diffusion of hydrophilic molecules into the periplasmic space. Owing to the strong hydrophilicity of their amino acid sequence and the nature of their secondary structure (beta strands), conventional hydropathy methods for predicting membrane topology are useless for this class of protein. The large number of available porin amino acid sequences was exploited to improve the accuracy of the prediction in combination with tools detecting amphipathicity of secondary structure. Using the constraints of beta-sheet structure these porins are predicted to contain 16 membrane-spanning strands, 14 of which are common to the two (enteric and the neisserial) porin subfamilies.     

Structural analysis of a peptide fragment of transmembrane transporter protein bilitranslocase

Using a combination of genomic and post-genomic approaches is rapidly altering the number of identified human influx carriers. A transmembrane protein bilitranslocase (TCDB 2.A.65) has long attracted attention because of its function as an organic anion carrier. It has also been identified as a potential membrane transporter for cellular uptake of several drugs and due to its implication in drug uptake, it is extremely important to advance the knowledge about its structure. However, at present, only the primary structure of bilitranslocase is known. In our work, transmembrane subunits of bilitranslocase were predicted by a previously developed chemometrics model and the stability of these polypeptide chains were studied by molecular dynamics (MD) simulation. Furthermore, sodium dodecyl sulfate (SDS) micelles were used as a model of cell membrane and herein we present a high-resolution 3D structure of an 18 amino acid residues long peptide corresponding to the third transmembrane part of bilitranslocase obtained by use of multidimensional NMR spectroscopy. It has been experimentally confirmed that one of the transmembrane segments of bilitranslocase has alpha helical structure with hydrophilic amino acid residues oriented towards one side, thus capable of forming a channel in the membrane.     

A fast method for the quantitative estimation of the distribution of hydrophobic and hydrophilic segments in alpha-helices of membrane proteins

The work presents a fast quantitative approach for estimating the orientations of hydrophilic and hydrophobic regions in the helical wheels of membrane-spanning alpha-helices of transmembrane proteins. The common hydropathy analysis provides an estimate of the integral hydrophobicity in a moving window which scans an amino acid sequence. The new parameter, orientation hydrophobicity, is based on the estimate of hydrophobicity of the angular segment that scans the helical wheel of a given amino acid sequence. The corresponding procedure involves the treatment of transmembrane helices as cylinders with equal surface elements for each amino acid residue. The orientation hydrophobicity, P(phi), phi = 0-360 degrees, of a helical cylinder is given as a sum of hydrophobicities of individual amino acids which are taken as the S-shaped functions of the angle between the centre of amino acid surface element and the centre of the segment. Non-zero contribution to P(phi) comes only from the amino acids belonging to the angular segment for a given angle phi. The size of the angular segment is related to the size of the channel pore. The amplitudes of amino acid S-functions are calibrated in the way that their maximum values (reached when the amino acid is completely exposed into the pore) are equal to the corresponding hydropathy index in the selected scale (here taken as Goldman-Engelman-Steitz hydropathy scale). The given procedure is applied in the studies of three ionic channels with well characterized three-dimensional structures where the channel pore is formed by a bundle of alpha-helices: cholera toxin B, nicotinic acetylcholine homopentameric alpha7 receptor, and phospholamban. The estimated maximum of hydrophilic properties at the helical wheels are in a good agreement with the spatial orientations of alpha-helices in the corresponding channel pores.     

Nucleotide sequence of the pldB gene and characteristics of deduced amino acid sequence of lysophospholipase L2 in Escherichia coli

The nucleotide sequence of the pldB gene of Escherichia coli K-12, which codes for lysophospholipase L2 located in the inner membrane, was determined. The deduced amino acid sequence of lysophospholipase L2 contains 340 amino acid residues, resulting in a protein with a molecular weight of 38,934. It is characterized by a high content of arginine residues (36 out of 340 residues). The amino acid sequence near the NH2-terminus of the protein is composed of a large number of polar or charged amino acid residues, suggesting that this region cannot be a signal peptide. The hydropathy profile of the deduced amino acid sequence of lysophospholipase L2 was studied. Most of the region was rather hydrophilic, and there was no stretch of hydrophobic amino acid region, such as might be predicted to traverse the lipid bilayer. These results are consistent with the experimental observation that lysophospholipase L2 is extracted by salt solution from the membrane fraction, and it may be classified as a peripheral membrane protein. Computer analysis showed that there is no homology in amino acid sequences between lysophospholipase L2 and other extracellular phospholipases, as well as detergent-resistant phospholipase A, which is another membrane-bound phospholipase in E. coli and whose DNA sequence was determined (Homma, H., Kobayashi, T., Chiba, N., Karasawa, K., Mizushima, H., Kudo, I., Inoue, K., Ideka, H., Sekiguchi, M., & Nojima, S. (1984) J. Biochem. 96, 1655-1664). This is the first report of the primary structure of a lysophospholipase.     

Molecular cloning and sequencing of cDNA for yeast porin, an outer mitochondrial membrane protein: a search for targeting signal in the primary structure.

We have cloned a full-length cDNA for yeast porin, the major outer mitochondrial membrane protein from Saccharomyces cerevisiae, and determined its nucleotide sequence. The primary structure of the protein, deduced from the nucleotide sequence, consisted of 283 amino acid residues and its NH2-terminal sequence, Met-Ser-Pro-Pro-Val-Tyr-Ser, coincided with that determined by Edman degradation for yeast porin, except that the initiator methionine was missing in the mature protein. The deduced sequence had an overall polarity index of 46.3%, a value which falls in the normal range for soluble proteins. An evaluation of hydropathy of the protein indicated that the NH2-terminal one third was relatively hydrophilic and the rest of the molecule was rather hydrophobic. An interesting finding was that the NH2-terminal region of yeast porin (consisting of some 50 amino acid residues) shows structural features that resemble those of the corresponding portion of 70-kd protein, which is also a yeast outer mitochondrial membrane protein. We postulate that this NH2-terminal sequence, like that of 70-kd protein, is required for targeting the porin to the outer mitochondrial membrane.     

Protein acylation in Tetrahymena

Examination of exhaustively delipidated Tetrahymena mimbres cells by sodium dodecyl sulfate-polyacrylamide gel electrophoresis revealed the presence of several protein bands containing covalently linked fatty acids. Palmitic (16:0) and stearic (18:0) acids together accounted for approximately 90% of the protein-linked acyl chains, with myristic acid (14:0) comprising most of the remainder. Each of these three fatty acids was present mainly in alkali-stable linkage, indicating that unlike most other systems examined, fatty acids are attached to proteins of Tetrahymena principally by amide bonds. Smaller proportions of the acyl chains were susceptible to release by hydroxylaminolysis or by alkaline hydrolysis as would be expected from an ester linkage. The protein-bound acyl chains accounted for 0.3% of the cells' total fatty acids. They closely resembled in composition the highly saturated free fatty acid pool but not the vast pool of glycerolipid-associated fatty acids, which were mainly unsaturated. Cells subjected to thermal stress by rapid chilling from 39 to 15 degrees C responded by sharply increasing the ratio of palmitate to stearate in covalent association with proteins.     

Energetics, stability, and prediction of transmembrane helices

We show that the peptide backbone of an alpha-helix places a severe thermodynamic constraint on transmembrane (TM) stability. Neglect of this constraint by commonly used hydrophobicity scales underlies the notorious uncertainty of TM helix prediction by sliding-window hydropathy plots of membrane protein (MP) amino acid sequences. We find that an experiment-based whole-residue hydropathy scale (WW scale), which includes the backbone constraint, identifies TM helices of membrane proteins with an accuracy greater than 99 %. Furthermore, it correctly predicts the minimum hydrophobicity required for stable single-helix TM insertion observed in Escherichia coli. In order to improve membrane protein topology prediction further, we introduce the augmented WW (aWW) scale, which accounts for the energetics of salt-bridge formation. An important issue for genomic analysis is the ability of the hydropathy plot method to distinguish membrane from soluble proteins. We find that the method falsely predicts 17 to 43 % of a set of soluble proteins to be MPs, depending upon the hydropathy scale used.