Importance of hydrophobic cluster formation through long-range contacts in the folding transition state of two-state proteins

Understanding the folding pathways of proteins is a challenging task. The Phi value approach provides a detailed understanding of transition-state structures of folded proteins. In this work, we have computed the hydrophobicity associated with each residue in the folded state of 16 two-state proteins and compared the Phi values of each mutant residue. We found that most of the residues with high Phi value coincide with local maximum in surrounding hydrophobicity, or have nearby residues that show such maximum in hydrophobicity, indicating the importance of hydrophobic interactions in the transition state. We have tested our approach to different structural classes of proteins, such as alpha-helical, SH3 domains of all-beta proteins, beta-sandwich, and alpha/beta proteins, and we observed a good agreement with experimental results. Further, we have proposed a hydrophobic contact network pattern to relate the Phi values with long-range contacts, which will be helpful to understand the transition-state structures of folded proteins. The present approach could be used to identify potential hydrophobic clusters that may form through long-range contacts during the transition state.     

The role of hydrophobicity patterns in prion folding as revealed by recurrence quantification analysis of primary structure

It has been suggested that the number and strength of local contacts are the major factors governing conformation accessibility of model two ground-state polypeptide chains. This phenomenology has been posed as a possible factor influencing prion folding. To test this conjecture, recurrence quantification analysis was applied to two model 36mers, and the Syrian hamster prion protein. A unique divergence of the radius function for the recurrence quantification variable %DET of hydrophobicity patterns was observed for both 36mers, and in a critical region of the hamster prion protein. This divergence suggests a partition between strong short- and long-range hydrophobicity patterns, and may be an important factor in prion phenomenology, along with other global thermodynamic factors.     

Structural analysis of residues involving cation-pi interactions in different folding types of membrane proteins

Cation-pi interactions play an important role to the stability of protein structures. In our earlier work, we have analyzed the influence and energetic contribution of cation-pi interactions in three-dimensional structures of membrane proteins. In this work, we investigate the characteristic features of residues that are involved in cation-pi interactions. We have computed several parameters, such as surrounding hydrophobicity, number of long-range contacts, conservation score and normalized B-factor for all these residues and identified their location, whether in the membrane or at surface. We found that the cation-pi interactions are mainly formed by long-range interactions. The cationic residues involved in cation-pi interactions have higher surrounding hydrophobicity than their average values in the whole dataset and an opposite trend is observed for aromatic residues. In transmembrane helical proteins, except Phe, all other residues that are responsible for cation-pi interactions are highly conserved with other related protein sequences whereas in transmembrane strand proteins, an appreciable conservation is observed only for Arg. The analysis on the flexibility of residues reveals that the cation-pi interaction forming residues are more stable than other residues. The results obtained in the present study would be helpful to understand the role of cation-pi interactions in the structure and folding of membrane proteins.     

Determination of intrinsic hydrophilicity/hydrophobicity of amino acid side chains in peptides in the absence of nearest-neighbor or conformational effects

Understanding the hydrophilicity/hydrophobicity of amino acid side chains in peptides/proteins is one the most important aspects of biology. Though many hydrophilicity/hydrophobicity scales have been generated, an "intrinsic" scale has yet to be achieved. "Intrinsic" implies the maximum possible hydrophilicity/hydrophobicity of side chains in the absence of nearest-neighbor or conformational effects that would decrease the full expression of the side-chain hydrophilicity/hydrophobicity when the side chain is in a polypeptide chain. Such a scale is the fundamental starting point for determining the parameters that affect side-chain hydrophobicity and for quantifying such effects in peptides and proteins. A 10-residue peptide sequence, Ac-X-G-A-K-G-A-G-V-G-L-amide, was designed to enable the determination of the intrinsic values, where position X was substituted by all 20 naturally occurring amino acids and norvaline, norleucine, and ornithine. The coefficients were determined by reversed-phase high-performance liquid chromatography using six different mobile phase conditions involving different pH values (2, 5, and 7), ion-pairing reagents, and the presence and absence of different salts. The results show that the intrinsic hydrophilicity/hydrophobicity of amino acid side chains in peptides (proteins) is independent of pH, buffer conditions, or whether C(8) or C(18) reversed-phase columns were used for 17 side chains (Gly, Ala, Cys, Pro, Val, nVal, Leu, nLeu, Ile, Met, Tyr, Phe, Trp, Ser, Thr, Asn, and Gln) and dependent on pH and buffer conditions, including the type of salt or ion-pairing reagent for potentially charged side chains (Orn, Lys, His, Arg, Asp, and Glu).     

Variation of amino acid properties in all-beta globular and outer membrane protein structures

Discriminating outer membrane (OM) proteins from globular proteins is an important task. The structural analysis of β-strands dominating globular (all-β) proteins and OM proteins provides useful insight to distinguish between them. In this work, we analyze the characteristic features of the 20 amino acid residues in all-β and OM proteins. We set up numerical indices for several properties of amino acid residues, such as, conformational parameters, surrounding hydrophobicity, accessible surface area and reduction in accessibility, and inter-residue contacts. We found that all the aromatic residues prefer to be in β-strands of both globular and OM proteins. The surrounding hydrophobicity of aromatic and non-polar amino acid residues in globular proteins is significantly higher than that of OM proteins. The residues Trp, Arg, Phe and Gln show a remarkable difference of reduction in accessibility between all-β globular (βG) and OM proteins. The positively charged residues, Lys and Arg in the membrane part of OM proteins have more number of contacts than globular proteins. Further, the behavior of the 20 amino acid residues in β-strand segments of globular and OM proteins have been discussed. The parameters developed in this work can be used for identifying transmembrane β-strands in OM proteins and for discriminating βG proteins from OM proteins.     

Analyses of the Sequence and Structural Properties Corresponding to Pentapeptide and Large Palindromes in Proteins

The analyses of 3967 representative proteins selected from the Protein Data Bank revealed the presence of 2803 pentapeptide and large palindrome sequences with known secondary structure conformation. These represent 2014 unique palindrome sequences. 60% palindromes are not associated with any regular secondary structure and 28% are in helix conformation, 11% in strand conformation and 1% in the coil conformation. The average solvent accessibility values are in the range between 0–155.28 Ų suggesting that the palindromes in proteins can be either buried, exposed to the solvent or share an intermittent property. The number of residue neighborhood contacts defined by interactions ≤ 3.2 Ǻ is in the range between 0–29 residues. Palindromes of the same length in helix, strand and coil conformation are associated with different amino acid residue preferences at the individual positions. Nearly, 20% palindromes interact with catalytic/active site residues, ligand or metal ions in proteins and may therefore be important for function in the corresponding protein. The average hydrophobicity values for the pentapeptide and large palindromes range between -4.3 to +4.32 and the number of palindromes is almost equally distributed between the negative and positive hydrophobicity values. The palindromes represent 107 different protein families and the hydrolases, transferases, oxidoreductases and lyases contain relatively large number of palindromes.     

Influence of long-range contacts and surrounding residues on the transition state structures of proteins

Understanding the parameters influencing the formation of transition state structures in proteins is an important problem in protein folding and kinetics. In this work, we have analyzed the structure-based parameters, surrounding hydrophobicity, secondary structure, solvent accessibility, number of medium- and long-range contacts, and surrounding residues for understanding the transition state structures of 15 proteins. The analysis of Φ-values shows that 29% of the studied 378 mutants have a Φ-value of more than 0.5. The combination of different structure-based parameters could discriminate the residues that have a Φ-value cutoff of more than 0.5 with a 5-fold cross-validation accuracy of 68%, which indicates that the surrounding residues and contacts play important roles in the formation of transition state structures. Systematic analysis on different proteins reveals that the proteins azurin, cold shock protein, and C-terminal domain of ribosomal protein L9 are influenced by the number of medium- and long-range proteins, whereas barnase, FK506 binding protein, and IM9 are influenced by surrounding residues. The discrimination accuracy lies in the ranges of 81–95% and 74–85% for these respective classes of protein. Furthermore, the combination of surrounding residues and contacts improved the accuracy up to 24% in other considered proteins. We suggest that the structure-based parameters along with noncovalent interactions and conservation of residues may aid in identifying the potential residues in the formation of transition state structures in proteins.     

Novel scales based on hydrophobicity indices for secondary protein structure

This paper is concerned with a branch of computational biology related to protein prediction and analysis of secondary structure of proteins. Although traditional methods use a simple amino acid composition to predict the secondary structure content, hydrophobicity has been recently found to improve the results in this and several related prediction tasks. To this end, we propose and analyze advantages of two new hydrophobicity index-based scales that incorporate information about long-range interactions along the protein sequence and contrast them with currently used raw hydrophobic index values. We also compare three leading hydrophobicity indices, i.e., Eisenberg's, Fauchere-Pliska's, and Cid's, using the proposed scales. The analysis is performed using fuzzy cognitive maps that quantify the strength of relation between the hydrophobicity scales/indices and the protein content values. A set of empirical tests that involve generation of fuzzy cognitive map models for a set of 200 low homology proteins have been performed. The results show that the secondary structure content along the protein sequence is characterized by about 2.5 times stronger relation with the two proposed hydrophobicity scales when compared with the currently used raw index values. The new scales exhibit stronger relation irrespective of the applied hydrobhobicity indices. Analysis of different scales shows superiority of the Eisenberg's hydrophobicity index, when used with the new scales. In contrast, the Fauchere-Pliska's index is found to perform better when compared with the two other indices when using raw hydrophobic index values that disregard the long-range interactions.     

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.     

ATP-induced conformational transition of denatured proteins.

Although the conformational change occurring in proteins upon ATP binding is important in many biological reactions, the mechanism by which ATP binding induces the conformational change is unknown. We found that ATP induces acid-unfolded (pH 2) ferricytochrome c or apomyoglobin to adopt a compact structure with a significant amount of alpha-helix and increased hydrophobicity. A very similar conformational transition was observed at neutral pH for an amphiphilic model polypeptide. The effectiveness of various adenine nucleotides in inducing the conformational transition was found to be proportional to their phosphate group contents, i.e., adenosine tetraphosphate greater than ATP greater than ADP greater than AMP. These results should be important when considering the mechanism of the ATP-induced conformational change in proteins during various biological reactions.     

Crystal structure of a collagen-like polypeptide with repeating sequence Pro-Hyp-Gly at 1.4 A resolution: implications for collagen hydration

The use of polypeptide models has proved to be a valuable tool to obtain accurate information on the collagen triple helix. Here we report the high resolution crystal structure of a collagen-like polypeptide with repeating sequence Pro-Hyp-Gly. The structure has been refined to an R(factor) of 0.137 and an R(free) of 0.163 using synchrotron diffraction data extending up to 1.4 A resolution. The polypeptide triple-helical structure binds a large number of water molecules, in contrast with a previous structure determination at lower resolution. The highly hydrated nature of this polypeptide confirms a number of previous studies conducted both in solution and in the crystal state. In addition, neighboring polypeptide triple helices are directly bound in the crystal through Hyp-Hyp hydrogen-bonding interactions. This finding supports the idea that Hyp residues may be important for the assembly of the triple helices in the collagen fibrils and may stabilize the fibrils by mediating direct contacts between neighboring molecules.     

Use of thioredoxin as a reporter to identify a subset of Escherichia coli signal sequences that promote signal recognition particle-dependent translocation

We have previously reported that the DsbA signal sequence promotes efficient, cotranslational translocation of the cytoplasmic protein thioredoxin-1 via the bacterial signal recognition particle (SRP) pathway. However, two commonly used signal sequences, those of PhoA and MalE, which promote export by a posttranslational mechanism, do not export thioredoxin. We proposed that this difference in efficiency of export was due to the rapid folding of thioredoxin in the cytoplasm; cotranslational export by the DsbA signal sequence avoids the problem of cytoplasmic folding (C. F. Schierle, M. Berkmen, D. Huber, C. Kumamoto, D. Boyd, and J. Beckwith, J. Bacteriol. 185:5706-5713, 2003). Here, we use thioredoxin as a reporter to distinguish SRP-dependent from non-SRP-dependent cleavable signal sequences. We screened signal sequences exhibiting a range of hydrophobicity values based on a method that estimates hydrophobicity. Successive iterations of screening and refining the method defined a threshold hydrophobicity required for SRP recognition. While all of the SRP-dependent signal sequences identified were above this threshold, there were also a few signal sequences above the threshold that did not utilize the SRP pathway. These results suggest that a simple measure of the hydrophobicity of a signal sequence is an important but not a sufficient indicator for SRP recognition. In addition, by fusing a number of both classes of signal sequences to DsbA, we found that DsbA utilizes an SRP-dependent signal sequence to achieve efficient export to the periplasm. Our results suggest that those proteins found to be exported by SRP-dependent signal sequences may require this mode of export because of their tendency to fold rapidly in the cytoplasm.     

Possible role of the highly conserved amino acids Trp-8 and Pro-13 in the N-terminal segment of the pigment-binding polypeptide LHI alpha of Rhodobacter capsulatus

Trp-8 and Pro-13 of the Rhodobacter capsulatus light-harvesting (LH) I alpha polypeptide are highly conserved among LHI and LHII alpha proteins of several species of the Rhodospirillaceae. Exchange of Trp-8 and Pro-13 to other amino acyl residues similar in structure and/or hydrophobicity indicates that Trp-8 is involved in the insertion of the LHI alpha polypeptide into the intracytoplasmic membrane (ICM). Pro-13, however, seems not to participate in the integration process of the LHI alpha protein but seems to be important for stable insertion of the LHI beta partner protein in the ICM.     

Correctly folded proteins make twice as many hydrophobic contacts

A novel statistical analysis of non-bonded contacts in a set of known protein structures shows that the natural residue types fall into five or six groups distinguishable by nearest neighbor preference. The observed pattern of contact specificities clearly reflects residue hydrophobicity and charge. Its most striking feature is that residues in the hydrophobic group make about twice as many contacts with one another as would be expected on a random basis. A similar increase in hydrophobic contact frequency can be observed at the level of individual proteins. Native proteins make, on average, about twice as many hydrophobic contacts as corresponding misfolded proteins, generated by computer. On the basis of these observations increased hydrophobic contact frequency is proposed as a simple model of the hydrophobic effect.     

Underlying hydrophobic sequence periodicity of protein tertiary structure

Hydropathy plots or window averages over local stretches of the sequence of residue hydrophobicity have revealed patterns related to various protein tertiary structural features. This has enabled identification of regions of the sequence that are at the surface or within the interior of globular soluble proteins, regions located within the lipid bilayer of transmembrane proteins, portions of the sequence that characterize repeating motifs, as well as motifs that usefully characterize different protein structural families. This, therefore, provides one example of the generally expressed maxim that "sequence determines structure". On the other hand, a number of previous investigations have shown the rapidly varying values of residue hydrophobicity along the sequence to be distributed approximately randomly. So one might question just how much of the sequence actually determines structure. It is, therefore, of interest to extract that part of this rapidly varying distribution of residue hydrophobicity that is responsible for the longer wavelength variations that correlate with protein tertiary structural features and to determine their prevalence within the entire distribution. This is accomplished by a finite Fourier analysis of the sequence of residue hydrophobicity and of a new measure of residue distance from the protein interior. Calculations are performed on a number of globins, immunoglobulins, cuprodoxins, and papain-like structures. The spectral power of the Fourier amplitudes of the frequencies extracted, whose inverse transforms underlie the windowed values of residue hydrophobicity is shown to be a small fraction of the total power of the hydrophobicity distribution and thereby consistent with a distribution that might appear to be predominantly random. The wide range of sequence identity between proteins having the same fold, all exhibiting similar small fractions of power amplitude that correlate with the longer wavelength inside-to-outside excursions of the amino acid residues, supports the general contention that close sequence identity is an expression of a close evolutionary relationship rather than an expression of structural similarity. Practical implications of the present analysis for protein structure prediction and engineering are also described.     

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%.     

Underlying Hydrophobic Sequence Periodicity of Protein Tertiary Structure

Abstract

Hydropathy plots or window averages over local stretches of the sequence of residue hydrophobicity have revealed patterns related to various protein tertiary structural features. This has enabled identification of regions of the sequence that are at the surface or within the interior of globular soluble proteins, regions located within the lipid bilayer of transmembrane proteins, portions of the sequence that characterize repeating motifs, as well as motifs that usefully characterize different protein structural families. This, therefore, provides one example of the generally expressed maxim that “sequence determines structure”. On the other hand, a number of previous investigations have shown the rapidly varying values of residue hydrophobicity along the sequence to be distributed approximately randomly. So one might question just how much of the sequence actually determines structure. It is, therefore, of interest to extract that part of this rapidly varying distribution of residue hydrophobicity that is responsible for the longer wavelength variations that correlate with protein tertiary structural features and to determine their prevalence within the entire distribution. This is accomplished by a finite Fourier analysis of the sequence of residue hydrophobicity and of a new measure of residue distance from the protein interior. Calculations are performed on a number of globins, immunoglobulins, cuprodoxins, and papain-like structures. The spectral power of the Fourier amplitudes of the frequencies extracted, whose inverse transforms underlie the windowed values of residue hydrophobicity is shown to be a small fraction of the total power of the hydrophobicity distribution and thereby consistent with a distribution that might appear to be predominantly random. The wide range of sequence identity between proteins having the same fold, all exhibiting similar small fractions of power amplitude that correlate with the longer wavelength inside-to- outside excursions of the amino acid residues, supports the general contention that close sequence identity is an expression of a close evolutionary relationship rather than an expression of structural similarity. Practical implications of the present analysis for protein structure prediction and engineering are also described.     

Genetic algorithm-based optimization of hydrophobicity tables

SUMMARY: The genomic abundance and pharmacological importance of membrane proteins have fueled efforts to identify them based solely on sequence information. Previous methods based on the physicochemical principle of a sliding window of hydrophobicity (hydropathy analysis) have been replaced by approaches based on hidden Markov models or neural networks which prevail due to their probabilistic orientation. In the current study, an optimization of the hydrophobicity tables used in hydropathy analysis is performed using a genetic algorithm. As such, the approach can be viewed as a synthesis between the physicochemically and statistically based methods. The resulting hydrophobicity tables lead to significant improvement in the prediction accuracy of hydropathy analysis. Furthermore, since hydropathy analysis is less dependent on the basis set of membrane proteins is used to hone the statistically based methods, as well as being faster, it may be valuable in the analysis of new genomes. Finally, the values obtained for each of the amino acids in the new hydrophobicity tables are discussed.     

Exploring the effects of sparse restraints on protein structure prediction

One of the main barriers to accurate computational protein structure prediction is searching the vast space of protein conformations. Distance restraints or inter‐residue contacts have been used to reduce this search space, easing the discovery of the correct folded state. It has been suggested that about 1 contact for every 12 residues may be sufficient to predict structure at fold level accuracy. Here, we use coarse‐grained structure‐based models in conjunction with molecular dynamics simulations to examine this empirical prediction. We generate sparse contact maps for 15 proteins of varying sequence lengths and topologies and find that given perfect secondary‐structural information, a small fraction of the native contact map (5%‐10%) suffices to fold proteins to their correct native states. We also find that different sparse maps are not equivalent and we make several observations about the type of maps that are successful at such structure prediction. Long range contacts are found to encode more information than shorter range ones, especially for α and αβ‐proteins. However, this distinction reduces for β‐proteins. Choosing contacts that are a consensus from successful maps gives predictive sparse maps as does choosing contacts that are well spread out over the protein structure. Additionally, the folding of proteins can also be used to choose predictive sparse maps. Overall, we conclude that structure‐based models can be used to understand the efficacy of structure‐prediction restraints and could, in future, be tuned to include specific force‐field interactions, secondary structure errors and noise in the sparse maps.     

Nitrite reductase from Pseudomonas aeruginosa: sequence of the gene and the protein

The gene coding for nitrite reductase of Pseudomonas aeruginosa has been cloned and its sequence determined. The coding region is 1707 bp long and contains information for a polypeptide chain of 568 amino acids. The sequence of the mature protein has been confirmed independently by extensive amino acid sequencing. The amino-terminus of the mature protein is located at Lys-26; the preceding 25 residue long extension shows the features typical of signal peptides. Therefore the enzyme is probably secreted into the periplasmic space. The mature protein is made of 543 amino acid residues and has a molecular mass of 60,204 Da. The c-heme-binding domain, which contains the only two Cys of the molecule, is located at the amino-terminal region. Analysis of the protein sequence in terms of hydrophobicity profile gives results consistent with the fact that the enzyme is fully water soluble and not membrane bound; the most hydrophilic region appears to correspond to the c-heme domain. Secondary structure predictions are in general agreement with previous analysis of circular dichroic data.