Physicochemical factors for discriminating between soluble and membrane proteins: hydrophobicity of helical segments and protein length

The average hydrophobicity of a polypeptide segment is considered to be the most important factor in the formation of transmembrane helices, and the partitioning of the most hydrophobic (MH) segment into the alternative nonpolar environment, a membrane or hydrophobic core of a globular protein may determine the type of protein produced. In order to elucidate the importance of the MH segment in determining which of the two types of protein results from a given amino acid sequence, we statistically studied the characteristics of MH helices, longer than 19 residues in length, in 97 membrane proteins whose three-dimensional structure or topology is known, as well as 397 soluble proteins selected from the Protein Data Bank. The average hydrophobicity of MH helices in membrane proteins had a characteristic relationship with the length of the protein. All MH helices in membrane proteins that were longer than 500 residues had a hydrophobicity greater than 1.75 (Kyte and Doolittle scale), while the MH helices in membrane proteins smaller than 100 residues could be as hydrophilic as 0.1. The possibility of developing a method to discriminate membrane proteins from soluble ones, based on the effect of size on the type of protein produced, is discussed.     

Prediction of natively unfolded regions in protein chain

We have shown that the ability of a protein to be in globular or in natively unfolded state (under native conditions) may be determined (besides low overall hydrophobicity and a large net charge) by such a property as the average environment density, the average number of residues enclosed at the given distance. A statistical scale of the average number of residues enclosed at the given distance for 20 types of amino acid residues in globular state has been created on the basis of 6626 protein structures. Using this scale for separation of 80 globular and 90 natively unfolded proteins we fail only in 11% of proteins (compared with 17% of errors which are observed if to use hydrophobicity scale). The present scale may be used both for prediction of form (folded or unfolded) of the native state of protein and for prediction of natively unfolded regions in protein chains. The results of comparison of our method of predicting natively unfolded regions with the other known methods show that our method has the highest fraction of correctly predicted natively unfolded regions (that is 87% and 77% if to make averaging over residues and over proteins correspondingly).     

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.     

TM Finder: a prediction program for transmembrane protein segments using a combination of hydrophobicity and nonpolar phase helicity scales

Based on the principle of dual prediction by segment hydrophobicity and nonpolar phase helicity, in concert with imposed threshold values of these two parameters, we developed the automated prediction program TM Finder that can successfully locate most transmembrane (TM) segments in proteins. The program uses the results of experiments on a series of host-guest TM segment mimic peptides of prototypic sequence KK AAAXAAAAAXAAWAAXAAAKKKK-amide (where X = each of the 20 commonly occurring amino acids) through which an HPLC-derived hydropathy scale, a hydrophobicity threshold for spontaneous membrane insertion, and a nonpolar phase helical propensity scale were determined. Using these scales, the optimized prediction algorithm of TM Finder defines TM segments by first searching for competent core segments using the combination of hydrophobicity and helicity scales, and then performs a gap-joining operation, which minimizes prediction bias caused by local hydrophilic residues and/or the choice of window size. In addition, the hydrophobicity threshold requirement enables TM Finder to distinguish reliably between membrane proteins and globular proteins, thereby adding an important dimension to the program. A full web version of the TM Finder program can be accessed at http://www.bioinformatics-canada.org/TM/.     

Location of ''continuous'' antigenic determinants in the protruding regions of proteins.

A simple method is described to locate 'antigenic' peptides from the alpha-carbon co-ordinates of a protein, based on protrusion from the protein's globular surface. A good correlation is found between those parts of a protein which protrude and the experimentally determined antigenic peptides in myoglobin, lysozyme and myohemerythrin. A comparison is made between the use of protrusion index, mobility, solvent accessibility and hydrophilicity for predicting the most likely antigenic peptides.     

Predicting alternate structure attainment and amyloidogenesis: a nonlinear signal analysis approach

Chain hydrophobicity values have been used in prediction of alternate structure attainment by a polypeptide. Nonlinear signal analysis on the hydrophobicity values gives important clues about the propensities of particular stretches of a protein to form local or nonlocal contacts. These contacts determine the folding behavior of a polypeptide and helps in predicting the final structure that can be attained. A nonlinear signal analysis called the recurrent quantification analysis has been carried out using the hydrophobicity values on a wide range of proteins obtained from human, plant, and fungal sources. Here, we show that such an analysis gives us an easy handle in determining sequences within the proteins that may be important in beta-sheet formation leading to amyloidosis.     

The globular domains of type VI collagen are related to the collagen-binding domains of cartilage matrix protein and von Willebrand factor.

Type VI collagen is a transformation-sensitive glycoprotein of the extracellular matrix of fibroblasts. We have isolated and sequenced several overlapping cDNA clones (4153 bp) which encode the entire alpha 2 subunit of chicken type VI collagen. The deduced amino acid sequence predicts that the alpha 2(VI) polypeptide consists of 1015 amino acid residues that are arranged in four domains: a hydrophobic signal peptide of 20 residues, an amino-terminal globular domain of 228 residues, a collagenous segment of 335 residues and a carboxy-terminal globular domain of 432 residues. The collagenous domain contains seven Arg-Gly-Asp tripeptide units, some of which are likely to be used as cell-binding sites. The globular domains contain three homologous repeats with an average length of 180 amino acid residues. These repeats show a striking similarity to the collagen-binding motifs found in von Willebrand factor and cartilage matrix protein. We therefore speculate that the globular domains of the alpha 2(VI) polypeptide may interact with collagenous structures.     

An algorithm for determining the conformation of polypeptide segments in proteins by systematic search

The feasibility of determining the conformation of segments of a polypeptide chain up to six residues in length in globular proteins by means of a systematic search through the possible conformations has been investigated. Trial conformations are generated by using representative sets of phi, psi, and chi angles that have been derived from an examination of the distributions of these angles in refined protein structures. A set of filters based on simple rules that protein structures obey is used to reduce the number of conformations to a manageable total. The most important filters are the maintenance of chain integrity and the avoidance of too-short van der Waals contacts with the rest of the protein and with other portions of the segment under construction. The procedure is intended to be used with approximate models so that allowance is made throughout for errors in the rest of the structure. All possible main chains are first constructed and then all possible side-chain conformations are built onto each of these. The electrostatic energy, including a solvent screening term, and the exposed hydrophobic area are evaluated for each accepted conformation. The method has been tested on two segments of chain in the trypsin like enzyme from Streptomyces griseus. It is found that there is a wide spread of energies among the accepted conformations, and the lowest energy ones have satisfactorily small root mean square deviations from the X-ray structure.     

A turn propensity scale for transmembrane helices.

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

Origins of globular structure in proteins

Thermodynamic incompatibility of polymers in a common solvent is possibly a driving force for formation and evolution of globular protein structures. Folding of polypeptide chains leads to a decrease in both excluded volume of molecules and chemical differences between surfaces of globular molecules with chemical information hidden in the hydrophobic interior. Folding of polypeptide chains results in 'molecular or thermodynamic mimicry' of globular proteins and in at least more than 10-fold higher phase separation threshold values of mixed protein solutions compared to those of classical polymers. Unusually high co-solubility might be necessary for efficient biological functioning of proteins, e.g. enzymes, enzyme inhibitors, etc.     

Determination of the dioxygen quenching constant for protein and model indole triplets

The decay of the indole triplet of single tryptophan-containing proteins and model compounds can be readily determined at room temperature in solution by monitoring the triplet absorption or emission following an exciting laser pulse. The dioxygen triplet quenching constants, can be measured for all these molecules and compared to the analogous singlet values determined by fluorescence methods. The dioxygen triplet quenching constant (tkq) ranged from a high of 5.1.10(9) M-1.s-1 for the exposed indole of corticotropin to a low of 0.1.10(9) M-1.s-1 for the buried indole of asparaginase. The ratio of these values with their respective dioxygen singlet quenching constants (skq), tkq/skq, ranged from 0.3 to 0.6 for aqueous exposed polypeptide indoles. For globular proteins the tkq/skq value is observed to be 0.2 +/- 0.1. This lower value for protein indoles is not attributable to 'bulk' environmental or hydrogen bonding effects, since the magnitude of tkq/skq (= 0.5 +/- 0.1) for model indoles was independent of solvent dielectric constant, polarity, and proticity. Temperature-dependence studies were done to test whether tkq could be used to characterize the nature of the protein matrix. The activation energy (Ea) for tkq was found to be 11 +/- 2 kcal/mol for most proteins. This Ea was independent of whether the indole side-chain was solvent exposed or buried in the non-aqueous protein interior. Large Ea values were also obtained for model indoles, naphthalene and nalidixic acid, dissolved in water, whereas the same compounds dissolved in 95% ethanol exhibited much smaller Ea values. These data, in combination with the observation that the tkq of model indoles is insensitive to changes in solvent viscosity, indicate that dioxygen quenching at the triplet level can not be easily used to characterize the dynamics of proteins.     

Prediction of natively unfolded regions in protein chains

Analysis showed that the globular or natively unfolded state of a protein can be inferred not only from a lower hydrophobicity or a higher charge, but also from the average environment density (average number of close residues located within a certain distance of a given one) of its residues. A database of 6626 protein structures was used to construct a statistical scale of the average number of close residues in globular structures for the 20 amino acids. The portion of false predictions in distinguishing between 80 globular and 90 natively unfolded proteins was 11% with the new scale and 17% with a hydrophobicity scale. The new scale proved suitable for predicting the folded or unfolded state for native proteins or the natively unfolded regions for protein chains. In comparisons with the available algorithms, the new method yielded the highest portion of true predictions (87 and 77% with averaging over residues and over proteins, respectively).     

A new approach to predict the helix/strand content of globular proteins

An improved multiple linear regression method has been proposed to predict the content of alpha-helix and beta-strand of a globular protein based on its primary sequence and structural class. The amino acid composition and the auto-correlation functions derived from the hydrophobicity profile of the primary sequence have been taken into account. However, only the compositions of a part of the amino acids and a part of the auto-correlation functions are selected as the regression terms, which lead to the least prediction error. The resubstitution test shows that the average absolute errors are 0.052 and 0.047 with the standard deviations 0.050 and 0.047 for the prediction of helix/strand content, respectively. A rigorous cross-validation test, the jackknife test shows that the average absolute errors are 0.058 and 0.053 with the standard deviations 0.057 and 0.053 for the prediction of helix/strand content, respectively. Both tests indicate the self-consistency and the extrapolating effectiveness of the new method. The high prediction accuracy means that the method is suitable for practical applications.     

Structural features of transmembrane helices

A total of 160 transmembrane helices of 15 non-homologous high-resolution X-ray protein structures have been analyzed in respect of their structural features. The dihedral angles and hydrogen bonds of the helical sections that span the hydrophobic interior of the lipid bilayer have been investigated. The Ramachandran plot of protein channels and solute transporters exhibit a significant shift Delta (phi- and psi-angles) of Delta mean (+4.5 degrees and -5.4 degrees ), compared to a reference group of 151 alpha-helices of the same average length derived from water-soluble globular proteins. At the C-termini of transmembrane helices structural motifs equivalent to the Gly-caps of helices in globular proteins have been found, with two third of the transmembrane Gly-caps taking up a primary structure that is typically not found at helix termini exposed to a polar solvent. The structural particularities reported here are relevant for the three-dimensional modelling of membrane protein structures.     

Local sequence patterns of hydrophobicity and solvent accessibility in soluble globular proteins

We examined the variation in the solvent accessibility and hydrophobicity of the amino acids along the sequences of 58 soluble globular proteins with known tertiary structure. We found that there is a significant tendency for the accessibilities to run in clusters along the sequence but that the hydrophobicities are distributed without such nonrandom clusters. Theseresults suggest severe limitations on the power of sequence analysis tools that use average hydrophobicity scores of overlapping subsequences to predict accessibility.     

To be folded or to be unfolded?

The lack of ordered structure in “natively unfolded” proteins raises a general question: Are there intrinsic properties of amino acid residues that are responsible for the absence of fixed structure at physiological conditions? In this article, we demonstrate that the competence of a protein to be folded or to be unfolded may be determined by the property of amino acid residues to form a sufficient number of contacts in a globular state. The expected average number of contacts per residue calculated from the amino acid sequence alone (using the average number of contacts for 20 amino acid residues in globular proteins) can be used as one of the simple indicators of natively unfolded proteins. The prediction accuracy for the sets of 80 folded and 90 natively unfolded proteins reaches 89% if the expected average number of contacts is used as a parameter and 83% in the case of hydrophobicity. An optimal set of artificial parameters for 20 amino acid residues obtained by Monte Carlo algorithm to maximally separate the sets of 90 natively unfolded and 80 folded proteins demonstrates the upper limit for prediction accuracy, which is 95%.     

PresCont: predicting protein-protein interfaces utilizing four residue properties

An important task of computational biology is to identify those parts of a polypeptide chain, which are involved in interactions with other proteins. For this purpose, we have developed the program PresCont, which predicts in a robust manner amino acids that constitute protein-protein interfaces (PPIs). PresCont reaches state-of-the-art classification quality on the basis of only four residue properties that can be readily deduced from the 3D structure of an individual protein and a multiple sequence alignment (MSA) composed of homologs. The core of PresCont is a support vector machine, which assesses solvent-accessible surface area, hydrophobicity, conservation, and the local environment of each amino acid on the protein surface. For training and performance testing, we compiled three nonoverlapping datasets consisting of permanently formed or transient complexes, respectively. A comparison with SPPIDER, ProMate, and meta-PPISP showed that PresCont compares favorably with these highly sophisticated programs, and that its prediction quality is less dependent on the type of protein complex being considered. This balance is due to a mutual compensation of classification weaknesses observed for individual properties: For PPIs of permanent complexes, solvent-accessible surface and hydrophobicity contribute most to classification quality, for PPIs of transient complexes, the assessment of the local environment is most significant. Moreover, we show that for permanent complexes a segmentation of PPIs into core and rim residues has only a moderate influence on prediction quality. PresCont is available as a web service at http://www-bioinf.uni-regensburg.de/.     

A simple method for displaying the hydropathic character of a protein

A computer program that progressively evaluates the hydrophilicity and hydrophobicity of a protein along its amino acid sequence has been devised. For this purpose, a hydropathy scale has been composed wherein the hydrophilic and hydrophobic properties of each of the 20 amino acid side-chains is taken into consideration. The scale is based on an amalgam of experimental observations derived from the literature. The program uses a moving-segment approach that continuously determines the average hydropathy within a segment of predetermined length as it advances through the sequence. The consecutive scores are plotted from the amino to the carboxy terminus. At the same time, a midpoint line is printed that corresponds to the grand average of the hydropathy of the amino acid compositions found in most of the sequenced proteins. In the case of soluble, globular proteins there is a remarkable correspondence between the interior portions of their sequence and the regions appearing on the hydrophobic side of the midpoint line, as well as the exterior portions and the regions on the hydrophilic side. The correlation was demonstrated by comparisons between the plotted values and known structures determined by crystallography. In the case of membrane-bound proteins, the portions of their sequences that are located within the lipid bilayer are also clearly delineated by large uninterrupted areas on the hydrophobic side of the midpoint line. As such, the membrane-spanning segments of these proteins can be identified by this procedure. Although the method is not unique and embodies principles that have long been appreciated, its simplicity and its graphic nature make it a very useful tool for the evaluation of protein structures.     

A two-component nucleation model of protein hydrophobicity

A fundamental characteristic of soluble globular protein structure is a hydrophobic core and protein exterior comprised predominantly of hydrophilic residues. This distribution of amino acid residue hydrophobicity, from protein interior to exterior, has recently been profiled with the use of hydrophobic moments. The calculations enable comparison of the radial hydrophobicity distribution of different proteins and had revealed two features common to 30 proteins of diverse size and structure. One, a global feature, is the overall shape of the second-order ellipsoidal hydrophobic moment. The second, a specific feature, is a quasi-invariant hydrophobic-ratio of distances. Both features are dependent upon the rates of increase, from protein interior to exterior, of the accumulated numbers of hydrophobic and hydrophilic amino acid residues. These rates can be simulated simply with a two-component nucleation model of protein hydrophobicity. The model provides insight into the origin of the shape of the observed hydrophobic moment profiles and of the observed range of hydrophobic ratios. Consistent with observation, it is shown that a relatively wide range of hydrophobic and hydrophilic rates of increase yield a relatively narrow range of hydrophobic ratios. Furthermore, the model identifies one factor, the decrease in residue density with increasing distance from the protein interior, that is critical in providing the range of values that is comparable with the observed range.     

What does it mean to be natively unfolded?

Natively unfolded or intrinsically unstructured proteins constitute a unique group of the protein kingdom. The evolutionary persistence of such proteins represents strong evidence in the favor of their importance and raises intriguing questions about the role of protein disorders in biological processes. Additionally, natively unfolded proteins, with their lack of ordered structure, represent attractive targets for the biophysical studies of the unfolded polypeptide chain under physiological conditions in vitro. The goal of this study was to summarize the structural information on natively unfolded proteins in order to evaluate their major conformational characteristics. It appeared that natively unfolded proteins are characterized by low overall hydrophobicity and large net charge. They possess hydrodynamic properties typical of random coils in poor solvent, or premolten globule conformation. These proteins show a low level of ordered secondary structure and no tightly packed core. They are very flexible, but may adopt relatively rigid conformations in the presence of natural ligands. Finally, in comparison with the globular proteins, natively unfolded polypeptides possess 'turn out' responses to changes in the environment, as their structural complexities increase at high temperature or at extreme pH.