Hydrophobicity of transmembrane proteins: spatially profiling the distribution

A hallmark of soluble globular protein tertiary structure is a hydrophobic core and a protein exterior populated predominantly by hydrophilic residues. Recent hydrophobic moment profiling of the spatial distribution of 30 globular proteins of diverse size and structure had revealed features of this distribution that were comparable. Analogous profiling of the hydrophobicity distribution of the alpha-helical buried bundles of several transmembrane proteins, as the lipid/protein interface is approached from within the bilayer, reveals spatial hydrophobicity profiles that contrast with those obtained for the soluble proteins. The calculations, which enable relative changes of hydrophobicity to be simply identified over the entire spatial extent of the multimer within the lipid bilayer, show the accumulated zero-order moments of the bundles to be mainly inverted with respect to that found for the soluble proteins. This indicates a statistical increase in the average residue hydrophobic content as the lipid bilayer is approached. This result differs from that of a relatively recent calculation and qualitatively agrees with earlier calculations involving lipid exposed and buried residues of the alpha-helices of transmembrane proteins. Spatial profiling, over the entire spatial extent of the multimer with scaled values of residue hydrophobicity, provides information that is not available from calculations using lipid exposure alone.     

Molecular structure of an apolipoprotein determined at 2.5-A resolution

The three-dimensional structure of an apolipoprotein isolated from the African migratory locust Locusta migratoria has been determined by X-ray analysis to a resolution of 2.5 A. The overall molecular architecture of this protein consists of five long alpha-helices connected by short loops. As predicted from amino acid sequence analyses, these helices are distinctly amphiphilic with the hydrophobic residues pointing in toward the interior of the protein and the hydrophilic side chains facing outward. The molecule falls into the general category of up-and-down alpha-helical bundles as previously observed, for example, in cytochrome c'. Although the structure shows the presence of five long amphiphilic alpha-helices, the alpha-helical moment and hydrophobicity of the entire molecule fall into the range found for normal globular proteins. Thus, in order for the amphiphilic helices to play a role in the binding of the protein to a lipid surface, there must be a structural reorganization of the protein which exposes the hydrophobic interior to the lipid surface. The three-dimensional motif of this apolipoprotein is compatible with a model in which the molecule binds to the lipid surface via a relatively nonpolar end and then spreads on the surface in such a way as to cause the hydrophobic side chains of the helices to come in contact with the lipid surface, the charged and polar residues to remain in contact with water, and the overall helical motif of the protein to be maintained.     

Mammalian hemodynamics: a new similarity principle

An attempt is made to obtain information as to the spatial distribution of amino acid residues in globular proteins in terms of their chemical, physical, energetic and conformational properties. The crystallographic data on twenty-one protein molecules form the basis for the study. The properties of the residues, namely, hydrophobicity, polarity, acidity, molecular weight, bulkness, chromatographic index, refractivity, short/medium/long-range energetics, and powers to adopt a-helical, extended and bend structures are analysed by dividing the protein globule into six concentric shells containing equal numbers of residues. The results show that the decisive factor in determining the spatial position of a residue in a protein molecule is of composite nature involving a compromise between the various properties of the residue. The observed deviations from the general hydrophobic interior and hydrophilic exterior in globular proteins are nicely brought out. A valuable clue is obtained as to the directional properties of α-helical and extended structure segments in proteins. The relative buried, exposed and intermediate characters of the residues in protein globules are obtained in a very realistic approach.     

The spatial distribution of physical, chemical, energetic and conformational properties of amino acid residues in globular proteins

An attempt is made to obtain information as to the spatial distribution of amino acid residues in globular proteins in terms of their chemical, physical, energetic and conformational properties. The crystallographic data on twenty-one protein molecules form the basis for the study. The properties of the residues, namely, hydrophobicity, polarity, acidity, molecular weight, bulkness, chromatographic index, refractivity, short/medium/long-range energetics, and powers to adopt a-helical, extended and bend structures are analysed by dividing the protein globule into six concentric shells containing equal numbers of residues. The results show that the decisive factor in determining the spatial position of a residue in a protein molecule is of composite nature involving a compromise between the various properties of the residue. The observed deviations from the general hydrophobic interior and hydrophilic exterior in globular proteins are nicely brought out. A valuable clue is obtained as to the directional properties of α-helical and extended structure segments in proteins. The relative buried, exposed and intermediate characters of the residues in protein globules are obtained in a very realistic approach.     

Local hydrophobicity stabilizes secondary structures in proteins

The probability of occurrence of helix and β-sheet residues in 47 globular proteins was determined as a function of local hydrophobicity, which was defined by the sum of the Nozaki-Tanford transfer free energies at two nearest-neighbors on both sides of the amino acid sequence. In general, hydrophilic amino acids favor neither helix nor β-sheet formations when neighbor residues are also hydrophilic but favor helix formation at higher local hydrophobicity. On the other hand, some hydrophobic amino acids such as Met, Leu, and Ile favor helix formation when neighbor residues are hydrophilic. None of the hydrophobic amino acids favor β-sheet formation with hydrophilic neighbors, but most of them strongly favor β-sheet formation at high local hydrophobicity. When the average of 20 amino acids is taken, both helix and β-sheet residue probabilities are higher at higher local hydrophobicity, although the increase is steeper for β-sheets. Therefore, β-sheet formation is more influenced by local hydrophobicity than helix formation. Generally, helices are nearer the surface and tend to have hydrophilic and hydrophobic faces at opposite sides. The tendency of alternating regions of hydrophilic and hydrophobic residues in a helical sequence was revealed by calculating the correlation of the Nozaki-Tanford values. Such amphipathic helices may be important in protein–protein and protein–lipid interactions and in forming hydrophilic channels in the membrane. The choice of 30 nonhomologous proteins as the data set did not alter the above results.     

The hydrophobic moment and its use in the classification of amphiphilic structures (review)

Amphiphilic alpha-helices play a major role in membrane dependent processes and are manifested in the primary structure of a protein by the periodic appearance of hydrophobic residues. Based on these periodic sequences, the hydrophobic moment was introduced, , which essentially treats the hydrophobicity of amino acid residues as a two-dimensional vector sum and provides a measure of amphiphilicity within regular repeat structures. To identify putative amphiphilic alpha-helix forming sequences, hydrophobic moment analysis assumes an amino acid residue periodicity of 100 and scans protein primary structures to find the 11-residue window with maximal . Taken with the window's mean hydrophobicity, , hydrophobic moment plot analysis uses the coordinate pair, [, ] to classify alpha-helices as either surface active, globular or transmembrane. More recently, this latter analysis has been extended to recognize candidate oblique orientated alpha-helices. Here, the hydrophobic moment is reviewed and data to query the logic of using a fixed window length and a fixed residue angular periodicity in hydrophobic moment analysis are provided. In addition, problems associated with the use of such analysis to predict alpha-helix structure/function relationships are considered.     

Sequence analysis of phospholamban. Identification of phosphorylation sites and two major structural domains

Phospholamban is a regulatory protein in cardiac sarcoplasmic reticulum that is phosphorylated by cAMP- and Ca2+/calmodulin-dependent protein kinase activities. In this report, we present the partial amino acid sequence of canine cardiac phospholamban and the identification of the sites phosphorylated by these two protein kinases. Gas-phase protein sequencing was used to identify 20 NH2-terminal residues. Overlap peptides produced by trypsin or papain digestion extended the sequence 16 residues to give the following primary structure: Ser-Ala-Ile-Arg-Arg-Ala-Ser-Thr-Ile-Glu-Met-Pro-Gln-Gln-Ala- Arg-Gln-Asn-Leu-Gln-Asn-Leu-Phe-Ile-Asn-Phe-(Cys)-Leu-Ile-Leu-Ile-(Cys)- Leu-Leu-Leu-Ile-. Phospholamban phosphorylated by either cAMP-dependent or Ca2+/calmodulin-dependent protein kinase was cleaved with trypsin, and the major phosphorylated peptide (comprising greater than 70% of the incorporated 32P label) was purified by reverse-phase high performance liquid chromatography. The identical sequence was revealed for the radioactive peptide obtained from phospholamban phosphorylated by either kinase: Arg-Ala-Ser-Thr-Ile-Glu-Met-Pro-Gln-Gln-. The adjacent residues Ser7 and Thr8 of phospholamban were identified as the unique sites phosphorylated by cAMP- and Ca2+/calmodulin-dependent protein kinases, respectively. These results establish that phospholamban is an oligomer of small, identical polypeptide chains. A hydrophilic, cytoplasmically oriented NH2-terminal domain on each monomer contains the unique, adjacent residues phosphorylated by cAMP- and Ca2+/calmodulin-dependent protein kinase activities. Analysis by hydropathic profiling and secondary structure prediction suggests that phospholamban monomers also contain a hydrophobic domain, which could form amphipathic helices sufficiently long to traverse the sarcoplasmic reticulum membrane. A model of phospholamban as a pentamer is presented in which the amphipathic alpha-helix of each monomer is a subunit of the pentameric membrane-anchored domain, which is comprised of an exterior hydrophobic surface and an interior hydrophilic region containing polar side chains.     

Monte Carlo studies on equilibrium globular protein folding. II. Beta-barrel globular protein models

In the context of dynamic Monte Carlo simulations on a model protein confined to a tetrahedral lattice, the interplay of protein size and tertiary structure, and the requirements for an all-or-none transition to a unique native state, are investigated. Small model proteins having a primary sequence consisting of a central bend neutral region flanked by two tails having an alternating hydrophobic/hydrophilic pattern of residues are seen to undergo a continuous transition to a beta-hairpin collapsed state. On increasing the length of the tails, the beta-hairpin structural motif is found to be in equilibrium with a four-member beta-barrel. Further increase of the tail length results in the shift of the structural equilibrium to the four-member beta-barrel. The random coil to beta-barrel transition is of an all-or-none character, but while the central turn is always the desired native bend, the location of the turns involving the two external strands is variable. That is, beta-barrels having the external stands that are two residues out of register are also observed in the transition region. Introduction into the primary sequence of two additional regions that are at the very least neutral toward turn formation produces an all-or-none transition to the unique, native, four-member beta-barrel. Various factors that can augment the stability of the native conformation are explored. Overall, these folding simulations strongly indicate that the general rules of globular protein folding are rather robust--namely, one requires a general pattern of hydrophobic/hydrophilic residues that allow the protein to have a well-defined interior and exterior and the presence of regions in the amino acid sequence that at the very least are locally indifferent to turn formation. Since no site-specific interactions between hydrophobic and hydrophilic residues are required to produce a unique four-member beta-barrel, these simulations strongly suggest that site specificity is involved in structural fine-tuning.     

Mutational replacements in subtilisin 309. Val104 has a modulating effect on the P4 substrate preference.

The previous notion that the amino acid side chain at position 104 of subtilisins is involved in the binding of the side chain at position P4 of the substrate has been investigated. The amino acid residue Val104 in subtilisin 309 has been replaced by Ala, Arg, Asp, Phe, Ser, Trp and Tyr by site-directed mutagenesis. It is shown that the P4 specificity of this enzyme is not determined solely by the amino acid residue occupying position 104, as the enzyme exhibits a marked preference for aromatic groups in P4, regardless of the nature of the position-104 residue. With hydrophilic amino acid residues at this position, no involvement is seen in binding of either hydrophobic or hydrophilic amino acid residues at position P4 of the substrates. The substrate with Asp in P4 is an exception, as the preference for this substrate is increased dramatically by introduction of an arginine residue at position 104 in the enzyme, presumably due to a substrate-induced conformational change. However, when position 104 is occupied by hydrophobic residues, it is highly involved in binding of hydrophobic amino acid residues, either by increasing the hydrophobicity of S4 or by determining the size of the pocket. The results suggest that the amino acid residue at position 104 is mobile such that it is positioned in the S4 binding site only when it can interact favourably with the substrate's side chain at position P4.     

A comparative study of several membrane proteins from the nervous system

A comparative study is carried out on the amino acid composition and average hydrophobicity of several membrane proteins from the central and peripheral nervous system. For this purpose, the amino acids are grouped into different sets whose ratios are defined in terms of the proportion of hydrophilic and hydrophobic residues. The analysis results in the classification of proteolipid proteins and the cholinergic receptor protein of Electrophorus as integral proteins, and acetylcholinesterase and basic proteins as peripheral.     

On the significance of alternating patterns of polar and non-polar residues in beta-strands

A common assumption about protein sequences in beta-strands is that they have alternating patterns of polar and non-polar residues. It is thought that such patterns reflect the interior/exterior geometry of amino acid residue side-chains on a beta-sheet. Here we study the prevalence of simple hydrophobicity patterns in parallel and antiparallel beta-sheets in proteins of known structure and in the sequences of amyloidogenic proteins. The occurrence of 32 possible pentapeptide binary patterns (polar (P)/non-polar (N)) is computed in 1911 non-homologous protein structures. Despite their tendency to aggregate in experimentally designed proteins, the purely alternating hydrophobic/polar patterns (PNPNP and NPNPN) are most frequent in beta-sheets, typically occurring in antiparallel strands. The overall distribution of the pentapeptide binary patterns is significantly different in strands within parallel and antiparallel sheets. In both types of sheets, complementary patterns (where the hydrophobic and polar residues pair with one another) associate preferentially. We do not find alternating patterns to be common in amyloidogenic proteins or in short fragments involved directly in amyloid formation. However, we do note some similarities between patterns present in amyloidogenic sequences and those in parallel strands.     

Rise and fall of the delta globin gene

The complete nucleotide sequence of the gene phoE, which codes for the phosphate limitation inducible outer membrane pore protein of Escherichia coli K12 was established. The results show that PhoE protein is synthesized in a precursor form with a 21 amino acid residue amino-terminal extension. This peptide has the general characteristics of a signal sequence. The promoter region of phoE has no homlogy with the consensus sequence of E. coli promoter regions, but homologous sequences with the promoter region of phoA, the structural gene for alkaline phosphatase, were observed. The deduced amino acid sequence showed that the mature PhoE protein is composed of 330 amino acid residues with a calculated molecular weight of 36,782. A number of 81 charged amino acids was found scattered throughout the protein while no large stretches of hydrophobic amino acids were observed. Hydrophobicity and hydration profiles of PhoE protein showed five pronounced hydrophilic maxima which are all located in the region from the amino terminus to residue 212.When the deduced amino acid sequence of PhoE protein was compared with the established sequence of the OmpF pore protein, a number of 210 identical residues was found. Some aspects of the structure-function relationship of PhoE protein are discussed in view of the hydrophobicity and hydration profiles, and the homology between PhoE protein and OmpF protein.     

Relationship between structure and amino acid sequence of strongly twisted and coiled β-hairpins in globular proteins

β-Hairpins are widespread in proteins, and it is possible to find them both within β-sheets and separately. In this work, a comparative analysis of amino acid sequences of β-strands within strongly twisted β-hairpins from different structural protein subclasses has been conducted. Strongly twisted and coiled β-hairpin generates in the space a right double helix out of β-strands that are connected by a loop region (connections). The frequencies of amino acid residues on the internal (concave) and external (convex) surfaces of strongly twisted β-hairpins have been determined (220 β-hairpins from nonhomologous proteins were studied). The concave surface of these β-hairpins is mainly generated by hydrophobic residues, while the convex surface by hydrophilic residues; accordingly, the alternation of hydrophobic internal and hydrophilic external residues is observed in their amino acid sequences. Amino acid residues of glycine and alanine (especially in places of the largest twisting of the strands) were anomalously frequently found in internal positions of strongly twisted and coiled β-hairpins. It was established that internal positions never contain the proline residues, while external positions in the twisting region contain them in a relatively large amount. It was demonstrated that at least one amino acid residue in α_L- or ε-conformation is required for generation of relatively short (up to 7 amino acid residues) connection. As a rule, these positions are occupied by glycines. Thus, not only the alternation of hydrophobic and hydrophilic amino acid residues, but also the presence of one or two glycine residues in the connection region and the excess of glycines and alanines in the places of the largest strand twisting on the concave surface, as well as the presence of prolines on the convex surface, are required to generate a strongly twisted and coiled β-hairpin.     

A general model of the P2 protein of peripheral nervous system myelin based on secondary structure predictions, tertiary folding principles, and experimental observations

The amino acid sequence of the P2 protein of peripheral myelin was analyzed with regard to regions of probable alpha-helix, beta-structure, beta-turn, and unordered conformation by means of several algorithms commonly used to predict secondary structure in proteins. Because of the high beta-sheet content and virtual absence of alpha-helix shown by the circular dichroic spectra of the protein, a bias was introduced into the algorithms to favor the beta-structure over the alpha-helical conformation. In order to define those beta-sheet residues that could lie on the external hydrophilic surface of the protein and those that could lie in its hydrophobic interior, the predicted beta-strands were examined for charged and uncharged amino acids located at alternating positions in the sequence. The sequential beta-strands in the predicted secondary structure were then ordered into beta-sheets and aligned according to generally accepted tertiary folding principles and certain chemical properties peculiar to the P2 protein. The general model of the P2 protein that emerged was a "Greek key" beta-barrel, consisting of eight antiparallel beta-strands with a two-stranded ribbon of antiparallel beta-structure emerging from one end. The model has an uncharged, hydrophobic core and a highly hydrophilic surface. The two Cys residues, which form a disulfide, occur in a loop connecting two adjacent antiparallel strands. Two hydrophilic loops, each containing a cluster of acidic residues and a single Phe, protrude from one end of the molecule. The general model is consistent with many of the properties of the actual protein, including the relatively weak nature of its association with myelin lipids and the positions of amino acid substitutions. Alternative beta-strand orderings yield three specific models having different interstrand connections across the barrel ends.     

Analysis of accessible surface of residues in proteins

We analyzed the total, hydrophobic, and hydrophilic accessible surfaces (ASAs) of residues from a nonredundant bank of 587 3D structure proteins. In an extended fold, residues are classified into three families with respect to their hydrophobicity balance. As expected, residues lose part of their solvent-accessible surface with folding but the three groups remain. The decrease of accessibility is more pronounced for hydrophobic than hydrophilic residues. Amazingly, Lysine is the residue with the largest hydrophobic accessible surface in folded structures. Our analysis points out a clear difference between the mean (other studies) and median (this study) ASA values of hydrophobic residues, which should be taken into consideration for future investigations on a protein-accessible surface, in order to improve predictions requiring ASA values. The different secondary structures correspond to different accessibility of residues. Random coils, turns, and beta-structures (outside beta-sheets) are the most accessible folds, with an average of 30% accessibility. The helical residues are about 20% accessible, and the difference between the hydrophobic and the hydrophilic residues illustrates the amphipathy of many helices. Residues from beta-sheets are the most inaccessible to solvent (10% accessible). Hence, beta-sheets are the most appropriate structures to shield the hydrophobic parts of residues from water. We also show that there is an equal balance between the hydrophobic and the hydrophilic accessible surfaces of the 3D protein surfaces irrespective of the protein size. This results in a patchwork surface of hydrophobic and hydrophilic areas, which could be important for protein interactions and/or activity.     

Strong hydrophobic nature of cysteine residues in proteins.

The differences between disulfide-bonding cystine (Cys_SS) and free cysteine (Cys_SH) residues were examined by analyzing the statistical distribution of both types of residue in proteins of known structure. Surprisingly, Cys_SH residues display stronger hydrophobicity than Cys_SS residues. A detailed survey of atoms which come into contact with the sulfhydryl group (sulfur atom) of Cys_SH revealed those atoms are essentially the same in number and variety as those of the methyl group of isoleucine, but are quite different to those of the hydroxyl group of serine. Moreover, the relationships among amino acids were also determined using the 3D-profile table of known protein structures. Cys_SH was located in the hydrophobic cluster, along with residues such as Met, Trp and Tyr, and was clearly separated from Ser and Thr in the polar cluster. These results imply that free cysteines behave as strongly hydrophobic, and not hydrophilic, residues in proteins.     

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.