Optimal ratio between the average conformational entropy and the average energy of interaction between residues for fast protein folding

Based on available experimental data and using a theoretical model of protein folding, we demonstrate that there is an optimal ratio between the average conformational entropy and the average contact energy per residue for fast protein folding. A statistical analysis of the conformational entropy and the number of contacts per residue for 5829 protein domains from four main classes (α, β, α/β, α+β) shows that each class has its own characteristic average number of contacts per residue and average conformational entropy per residue. These class-specific characteristics determine the protein folding rates: α-proteins are the fastest to fold, β-proteins are the second fastest, α+β-proteins are the third, and α/β-proteins are the last to fold.     

Estimates of the loss of main-chain conformational entropy of different residues on protein folding.

The average contribution of conformational entropy for individual amino acid residues towards the free energy of protein folding is not well understood. We have developed empirical scales for the loss of the main-chain (torsion angles, phi and psi) conformational entropy by taking its side-chain into account. The analysis shows that the main-chain component of the total conformational entropy loss for a residue is significant and reflects intrinsic characteristics associated with individual residues. The values have direct correlation with the hydrophobicity values and this has important bearing on the folding process. Proteins 1999;36:332-339.     

The magnitude of the backbone conformational entropy change in protein folding

The magnitude of the conformational entropy change experienced by the peptide backbone upon protein folding was investigated experimentally and by computational analysis. Experimentally, two different pairs of mutants of a 33 amino acid peptide corresponding to the leucine zipper region of GCN4 were used for high-sensitivity microcalorimetric analysis. Each pair of mutants differed only by having alanine or glycine at a specific solvent-exposed position under conditions in which the differences in stability could be attributed to differences in the conformational entropy of the unfolded state. The mutants studied were characterized by different stabilities but had identical heat capacity changes of unfolding (ΔC_p), identical solvent-related entropies of unfolding (ΔS_solv), and identical enthalpies of unfolding (ΔH) at equivalent temperatures. Accordingly, the differences in stability between the different mutants could be attributed to differences in conformational entropy. The computational studies were aimed at generating the energy profile of backbone conformations as a function of the main chain dihedral angles ϕ and ϱ. The energy profiles permit a direct calculation of the probability distribution of different conformers and therefore of the conformational entropy of the backbone. The experimental results presented in this paper indicate that the presence of the methyl group in alanine reduces the conformational entropy of the peptide backbone by 2.46 ± 0.2 cal/K · mol with respect to that of glycine, consistent with a 3.4-fold reduction in the number of allowed conformations in the alanine-containing peptides. Similar results were obtained from the energy profiles. The computational analysis also indicates that the addition of further carbon atoms to the side chain had only a small effect as long as the side chains were unbranched at position β. A further reduction with respect to Ala of only 0.61 and 0.81 cal/K · mol in the backbone entropy was obtained for leucine and lysine, respectively. β-branching (Val) produces the largest decrease in conformational entropy (1.92 cal/K · mol less than Ala). Finally, the backbone entropy change associated with the unfolding of an α-helix is 6.51 cal/K · mol for glycine. These and previous results have allowed a complete estimation of the conformational entropy changes associated with protein folding. © 1996 Wiley-Liss, Inc.     

The backbone conformational entropy of protein folding: experimental measures from atomic force microscopy

The energy dissipated during the atomic force microscopy-based mechanical unfolding and extension of proteins is typically an order of magnitude greater than their folding free energy. The vast majority of the "excess" energy dissipated is thought to arise due to backbone conformational entropy losses as the solvated, random-coil unfolded state is stretched into an extended, low-entropy conformation. We have investigated this hypothesis in light of recent measurements of the energy dissipated during the mechanical unfolding of "polyproteins" comprised of multiple, homogeneous domains. Given the assumption that backbone conformational entropy losses account for the vast majority of the energy dissipated (an assumption supported by numerous lines of experimental evidence), we estimate that approximately 19(+/-2)J/(mol K residue) of entropy is lost during the extension of three mechanically stable beta-sheet polyproteins. If, as suggested by measured peak-to-peak extension distances, pulling proceeds to near completion, this estimate corresponds to the absolute backbone conformational entropy of the unfolded state. As such, it is exceedingly close to previous theoretical and semi-empirical estimates that place this value at approximately 20J/(mol K residue). The estimated backbone conformational entropy lost during the extension of two helical polyproteins, which, in contrast to the mechanically stable beta-sheet polyproteins, rupture at very low applied forces, is three- to sixfold less. Either previous estimates of the backbone conformational entropy are significantly in error, or the reduced mechanical strength of the helical proteins leads to the rupture of a subsequent domain before full extension (and thus complete entropy loss) is achieved.     

Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

Limited conformational space for early-stage protein folding simulation

MOTIVATION: The problem of early-stage protein folding is critical for protein structure prediction. The model presented introduces a common definition of protein structures which may be treated as the possible in silico early-stage form of the polypeptide chain. Limitation of the conformational space to the ellipse path on the Ramachandran map was tested as a possible sub-space to represent the early-stage structure for simulation of protein folding. The proposed conformational sub-space was developed on the basis of the backbone conformation, with side-chain interactions excluded. RESULTS: The ellipse-path-limited conformation of BPTI was created using the criterion of shortest distance between Phi, Psi angles in native form of protein and the Phi, Psi angles belonging to the ellipse. No knots were observed in the structure created according to ellipse-path conformational sub-space. The energy minimization procedure applied to ellipse-path derived conformation directed structural changes toward the native form of the protein with SS-bonds system introduced to the procedure. AVAILABILITY: Program 'Ellipse' to create the ellipse-path derived structure available on request: myroterm@cyf-kr.edu.pl     

On the entropy of protein folding.

The failure to appreciate that the hydration of polar groups is a major contribution to the entropy of protein unfolding has led to considerable underestimates for the loss of configurational freedom when a protein chain folds.     

Fractal structure and conformational entropy of protein chain

The structures of 25 proteins arbitrarily chosen are investigated by fractal geometry, and their fractal dimensions (D_f) and conformational entropies S(N₀) are calculated by Havlin—Ben Avraham and Monte Carlo method, respectively. Comparison of the D_f and S(N₀) gives the relation: D_f = 1.532 - 3.00 × 10⁻⁴ S(N₀). The entropy data obtained by Monte Carlo method for the chain of random self-avoiding walks confirm the prediction of renormalization group: S(N₀) = 1.544N₀ + 0.1667 In N₀ + 0.1570 where N₀ is the number of residues in a protein chain. Both the D_f and S(N₀) reflect the conformational properties of a protein molecular chain. The idea resulting from the present communication suggests that the thermodynamic behaviours of proteins may be related to multifractals.     

Designing potential energy functions for protein folding.

By following a consistent line of physical reasoning, some fundamental understanding about the foldability of proteins has been achieved. In recent years, this has led to the development of a number of successful algorithms for optimizing potential energy functions for folding protein models. The differences between the folding mechanisms of simple, contact-based lattice proteins and more traditional, realistic protein models, however, still call for further development of the potentials in addition to the optimization approaches.     

Inter-residue interactions in protein folding and stability

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about inter-residue interactions in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, we introduce the classification of inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes of globular and membrane proteins, and in various folds have been delineated. The development of contact potentials and the application of inter-residue contacts for predicting the structural class and secondary structures of globular proteins, solvent accessibility, fold recognition and ab initio tertiary structure prediction have been evaluated. Further, the relationship between inter-residue contacts and protein-folding rates has been highlighted. Moreover, the importance of inter-residue interactions in protein-folding kinetics and for understanding the stability of proteins has been discussed. In essence, the information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.     

Toward accurate prediction of pKa values for internal protein residues: the importance of conformational relaxation and desolvation energy

Proton uptake or release controls many important biological processes, such as energy transduction, virus replication, and catalysis. Accurate pK(a) prediction informs about proton pathways, thereby revealing detailed acid-base mechanisms. Physics-based methods in the framework of molecular dynamics simulations not only offer pK(a) predictions but also inform about the physical origins of pK(a) shifts and provide details of ionization-induced conformational relaxation and large-scale transitions. One such method is the recently developed continuous constant pH molecular dynamics (CPHMD) method, which has been shown to be an accurate and robust pK(a) prediction tool for naturally occurring titratable residues. To further examine the accuracy and limitations of CPHMD, we blindly predicted the pK(a) values for 87 titratable residues introduced in various hydrophobic regions of staphylococcal nuclease and variants. The predictions gave a root-mean-square deviation of 1.69 pK units from experiment, and there were only two pK(a)'s with errors greater than 3.5 pK units. Analysis of the conformational fluctuation of titrating side-chains in the context of the errors of calculated pK(a) values indicate that explicit treatment of conformational flexibility and the associated dielectric relaxation gives CPHMD a distinct advantage. Analysis of the sources of errors suggests that more accurate pK(a) predictions can be obtained for the most deeply buried residues by improving the accuracy in calculating desolvation energies. Furthermore, it is found that the generalized Born implicit-solvent model underlying the current CPHMD implementation slightly distorts the local conformational environment such that the inclusion of an explicit-solvent representation may offer improvement of accuracy.     

Importance of long-range interactions in protein folding

Long-range interactions play an active role in the stability of protein molecules. In this work, we have analyzed the importance of long-range interactions in different structural classes of globular proteins in terms of residue distances. We found that 85% of residues are involved in long-range contacts. The residues occurring in the range of 4-10 residues apart contribute more towards long-range contacts in all-alpha proteins while the range is 11-20 in all-beta proteins. The hydrophobic residues Cys, Ile and Val prefer the 11-20 range and all other residues prefer the 4-10 range. The residues in all-beta proteins have an average of 3-8 long-range contacts whereas the residues in other classes have 1-4 long-range contracts. Furthermore, the preference of residue pairs to the folding and stability will be discussed.     

Novel feature for catalytic protein residues reflecting interactions with other residues

Owing to their potential for systematic analysis, complex networks have been widely used in proteomics. Representing a protein structure as a topology network provides novel insight into understanding protein folding mechanisms, stability and function. Here, we develop a new feature to reveal correlations between residues using a protein structure network. In an original attempt to quantify the effects of several key residues on catalytic residues, a power function was used to model interactions between residues. The results indicate that focusing on a few residues is a feasible approach to identifying catalytic residues. The spatial environment surrounding a catalytic residue was analyzed in a layered manner. We present evidence that correlation between residues is related to their distance apart most environmental parameters of the outer layer make a smaller contribution to prediction and ii catalytic residues tend to be located near key positions in enzyme folds. Feature analysis revealed satisfactory performance for our features, which were combined with several conventional features in a prediction model for catalytic residues using a comprehensive data set from the Catalytic Site Atlas. Values of 88.6 for sensitivity and 88.4 for specificity were obtained by 10-fold cross-validation. These results suggest that these features reveal the mutual dependence of residues and are promising for further study of structure-function relationship.     

Effect of short- and long-range interactions on protein folding

The relative importance of short- and long-range interactions is examined using a Monte Carlo simulation of protein folding on bovine pancreatic trypsin inhibitor. The model of the protein and the interaction energies were parametrized using X-ray structures of 30 native proteins. A nearest neighbor Ising model is used to determine the conformational state at each stage of the Monte Carlo procedure. Long-range interactions are simulated by contact free energies which become effective as two residues, separated by four or more residues along the chain, approach each other, and by disulfide-bond energies. Short-range interactions for residues separated by one, two, or three residues along the chain are also modeled by contact free energies and by -helical hydrogen bonds. A hard-sphere model is used to represent repulsive interactions. The ratios of short- to long-range interactions studied are 1:1, 2:1, 1:2, 0:1, and 1:0; e.g., for the 2:1 ratio, short-range interactions are weighted twice as much as long-range interactions, and for the 1:0 ratio, long-range interactions are omitted. For each ratio of short- to long-range interactions, a native conformation is found by a Monte Carlo procedure, a segment of 11 residues (residue numbers 1–11) is then rotated away from the rest of the molecule [breaking the 5–55 native disulfide bond, and moving this segment so that the distance between the sulfur atoms of the 5 and 55 cystine side chains (averaged for all native conformations) increases from 3.9 to 7.3 Å], and the Monte Carlo simulation is carried out (allowing the conformation of the whole molecule to change) until equilibrium is attained. For each ratio, the refolded conformation is compared to the native one using triangular distance maps and differential geometry distance criteria. With ratios of short- to long-range interaction energies of 1:1 and 0:1, the native disulfide bond could be re-formed; with ratios of 2:1 and 1:2 it did not; and with the 1:0 ratio, even a stable native conformation was not achieved. Therefore, long-range interactions (in addition to short-range ones) are required to bring remote parts of the protein together and to stabilize its native conformation.NIH Postdoctoral Fellow, 1977–1978.     

Specific non-native hydrophobic interactions in a hidden folding intermediate: implications for protein folding

Structures of intermediates and transition states in protein folding are usually characterized by amide hydrogen exchange and protein engineering methods and interpreted on the basis of the assumption that they have native-like conformations. We were able to stabilize and determine the high-resolution structure of a partially unfolded intermediate that exists after the rate-limiting step of a four-helix bundle protein, Rd-apocyt b(562), by multidimensional NMR methods. The intermediate has partial native-like secondary structure and backbone topology, consistent with our earlier native state hydrogen exchange results. However, non-native hydrophobic interactions exist throughout the structure. These and other results in the literature suggest that non-native hydrophobic interactions may occur generally in partially folded states. This can alter the interpretation of mutational protein engineering results in terms of native-like side chain interactions. In addition, since the intermediate exists after the rate-limiting step and Rd-apocyt b(562) folds very rapidly (k(f) approximately 10(4) s(-1)), these results suggest that non-native hydrophobic interactions, in the absence of topological misfolding, are repaired too rapidly to slow folding and cause the accumulation of folding intermediates. More generally, these results illustrate an approach for determining the high-resolution structure of folding intermediates.     

Inversion of the balance between hydrophobic and hydrogen bonding interactions in protein folding and aggregation

Identifying the forces that drive proteins to misfold and aggregate, rather than to fold into their functional states, is fundamental to our understanding of living systems and to our ability to combat protein deposition disorders such as Alzheimer's disease and the spongiform encephalopathies. We report here the finding that the balance between hydrophobic and hydrogen bonding interactions is different for proteins in the processes of folding to their native states and misfolding to the alternative amyloid structures. We find that the minima of the protein free energy landscape for folding and misfolding tend to be respectively dominated by hydrophobic and by hydrogen bonding interactions. These results characterise the nature of the interactions that determine the competition between folding and misfolding of proteins by revealing that the stability of native proteins is primarily determined by hydrophobic interactions between side-chains, while the stability of amyloid fibrils depends more on backbone intermolecular hydrogen bonding interactions.