Effect of proline residues on protein folding

Conformational energy calculations have been used to study the role of the proline residues in the folding of bovine pancreatic trypsin inhibitor. In the calculation, each of the four proline residues of this small protein is forced from the trans to cis peptide isomer while still part of the native folded structure. The cis proline residue can always be accommodated by small changes of the native conformation (< 1 Å root-mean-square deviation). For three of the four proline residues, Pro2, Pro9 and Pro 13, being in the cis form is calculated to destabilize the folded conformation by less than 11 kcal/mol, suggesting that rapid folding to a stable native-like conformation can occur with either isomeric form. For one of these three, Pro13, the destabilization is only 1 kcal/mol, suggesting the existence of an alternative folded native conformation with Pro13 cis. The fourth proline residue, Pro8, is calculated to destabilize the native conformation by so much (33 kcal/mol) that it will block folding in the manner proposed by Brandts et al. (1975).     

Learning about protein folding via potential functions

Over the last few years we have developed an empirical potential function that solves the protein structure recognition problem: given the sequence for an n-residue globular protein and a collection of plausible protein conformations, including the native conformation for that sequence, identify the correct, native conformation. Having determined this potential on the basis of only some 6500 native/nonnative pairs of structures for 58 proteins, we find it recognizes the native conformation for essentially all compact, soluble, globular proteins having known native conformations in comparisons with 10⁴ to 10⁶ reasonable alternative conformations apiece. In this sense, the potential encodes nearly all the essential features of globular protein conformational preference. In addition it “knows” about many additional factors in protein folding, such as the stabilization of multimeric proteins, quaternary structure, the role of disulfide bridges and ligands, proproteins vs. processed proteins, and minimal strand lengths in globular proteins. Comparisons are made with other sorts of protein folding problems, and applications in protein conformational determination and prediction are discussed. © 1994 Wiley-Liss, Inc.     

Statistical mechanics of protein folding by exhaustive enumeration

It is hard to construct theories for the folding of globular proteins because they are large and complicated molecules having enormous numbers of nonnative conformations and having native states that are complicated to describe. Statistical mechanical theories of protein folding are constructed around major simplifying assumptions about the energy as a function of conformation and/or simplifications of the representation of the polypeptide chain, such as one point per residue on a cubic lattice. It is not clear how the results of these theories are affected by their various simplifications. Here we take a very different simplification approach where the chain is accurately represented and the energy of each conformation is calculated by a not unreasonable empirical function. However, the set of amino acid sequences and allowed conformations is so restricted that it becomes computationally feasible to examine them all. Hence we are able to calculate melting curves for thermal denaturation as well as the detailed kinetic pathway of refolding. Such calculations are based on a novel representation of the conformations as points in an abstract 12-dimensional Euclidean conformation space. Fast folding sequences have relatively high melting temperatures, native structures with relatively low energies, small kinetic barriers between local minima, and relatively many conformations in the global energy minimum's watershed. In contrast to other folding theories, these models show no necessary relationship between fast folding and an overall funnel shape to the energy surface, or a large energy gap between the native and the lowest nonnative structure, or the depth of the native energy minimum compared to the roughness of the energy landscape. Proteins 32:425–437, 1998. © 1998 Wiley-Liss, Inc.     

Equilibrium folding and unfolding pathways for a model protein

The folding–unfolding process of reduced bovine pancreatic trypsin inhibitor was investigated with an idealized model employing approximate free energies. The protein is regarded to consist of only C^α and C^β atoms. The backbone dihedral angles are the only conformational variables and are permitted to take discrete values at every 10°. Intraresidue energies consist of two terms: an empirical part taken from the observed frequency distributions of (?,ψ) and an additional favorable energy assigned to the native conformation of each residue. Interresidue interactions are simplified by assuming that there is an attractive energy operative only between residue pairs in close contact in the native structure. A total of 230,000 molecular conformations, with no atomic overlaps, ranging from the native state to the denatured state, are randomly generated by changing the sampling bias. Each conformation is classified according to its conformational energy, F; a conformational entropy, S(F) is estimated for each value of F from the number of samples. The dependence of S(F) on energy reveals that the folding–unfolding transition for this idealized model is an “all-or-none” type; this is attributable to the specific long-range interactions. Interresidue contact probabilities, averaged over samples representing various stages of folding, serve to characterize folding intermediates. Most probable equilibrium pathways for the folding–unfolding transition are constructed by connecting conformationally similar intermediates. The specific details obtained for bovine pancreatic trypsin inhibitor are as follows: (1) Folding begins with the appearance of nativelike medium-range contacts at a β-turn and at the α-helix. (2) These grow to include the native pair of interacting β-strands. This state includes intact regular secondary conformations, as well as the interstrand sheet contacts, and corresponds to an activated state with the highest free energy on the pathway. (3) Additional native long-range contacts are completely formed either toward the amino terminus or toward the carboxyl terminus. (4) In a final step, the missing contacts appear. Although these folding pathways for this model are not consistent with experimental reports, it does indicate multiple folding pathways. The method is general and can be applied to any set of calculated conformational energies and furthermore permits investigation of gross folding features.     

Ubiquitin folds via a flip-twist-lock mechanism

To perform specific functional activities, the majority of proteins should fold into their distinct three-dimensional conformations. However, the biologically active conformation of a protein is generally found to be marginally stable than the other conformations that the chain can adopt. How a protein finds its native conformation from its post-synthesis unfolded structure in a complex conformational landscape is the unsolved question that still drives the protein folding community. Here, we report the folding mechanism of a globular protein, ubiquitin, from its chemically denatured state using all-atom molecular dynamics simulations. From the kinetic analysis of the simulated trajectories we show that the folding process can be described by the hydrophobic collapse mechanism, initiated by the “dewetting transition”, and subsequently assisted by the origination of an N-terminal folding nucleus, and finally supported by a native salt-bridge interaction between K11 and E34. We show that ubiquitin folds via an intermediate. Finally, we confirm the presence of “biological water” and explain its role to the folding process.     

Energy embedding of trypsin inhibitor

Energy embedding has been shown recently to be a useful extension of the distance geometry approach to conformational calculations in the case of very small molecules and simple energy functions. This paper tests the ability of energy embedding to locate low energy conformations satisfying both weak and strong geometric constraints when the molecule is the small protein, bovine pancreatic trypsin inhibitor, and the energy function is the complicated Oobatake-Crippen residue–residue potential. Using the potential function alone, the algorithm reaches a structure with energy lower than that of the native conformation, but with little resemblance to it. Aided by numerous geometric constraints, such as preformed secondary structure segments, the algorithm again finds a local minimum with energy better than that of the native, and with only 3.3 Å rms deviation from it. This is significantly closer to the native value than can be obtained using standard distance geometry and the geometric constraints alone. Thus, energy embedding using the Oobatake-Crippen potential function is a significant help in finding native conformations of proteins. However, additional trials on a hairpin bend fragment of trypsin inhibitor demonstrate the potential's shortcomings in encouraging proper secondary structure.     

Structural basis of hierarchical multiple substates of a protein. III: Side chain and main chain local conformations

An analysis is carried out of differences in the minimum energy conformations obtained in the previous paper by energy minimization starting from conformations sampled by a Monte Carlo simulation of conformational fluctuations in the native state of a globular protein, bovine pancreatic trypsin inhibitor. Main conformational differences in each pair of energy minima are found usually localized in several side chains and in a few local main chain segments. Such side chains and local main chain segments are found to take a few distinct local conformations in the minimum energy conformations. Energy minimum conformations can thus be described in terms of combinations of these multiple local conformations.     

A simple lattice model that exhibits a protein-like cooperative all-or-none folding transition

In a recent paper (D. Gront et al., Journal of Chemical Physics, Vol. 115, pp. 1569, 2001) we applied a simple combination of the Replica Exchange Monte Carlo and the Histogram methods in the computational studies of a simplified protein lattice model containing hydrophobic and polar units and sequence-dependent local stiffness. A well-defined, relatively complex Greek-key topology, ground (native) conformations was found; however, the cooperativity of the folding transition was very low. Here we describe a modified minimal model of the same Greek-key motif for which the folding transition is very cooperative and has all the features of the "all-or-none" transition typical of real globular proteins. It is demonstrated that the all-or-none transition arises from the interplay between local stiffness and properly defined tertiary interactions. The tertiary interactions are directional, mimicking the packing preferences seen in proteins. The model properties are compared with other minimal protein-like models, and we argue that the model presented here captures essential physics of protein folding (structurally well-defined protein-like native conformation and cooperative all-or-none folding transition).     

A new protein folding recognition potential function

On the study of protein inverse folding problem, one goal is to find simple and efficient potential to evaluate the compatibility between structure and a given sequence. We present here a novo empirical mean force potential to address the importance of electrostatic interactions in protein inverse folding study. It is based on protein main chain polar fraction and constructed in a way similar with Sippl's from a database of 64 known independent three-dimensional protein structures. This potential was applied to recognize the protein native conformations among a conformation pool. Calculated results show that this potential is powerful in picking out native conformations, in addition it can also find structure similarity between proteins with low sequence similarity. The success of this new potential clearly shows the importance of electrostatic factors in protein inverse folding studies. © 1995 Wiley-Liss, Inc.     

Two-dimensional 1H nuclear magnetic resonance study of the (5-55) single-disulphide folding intermediate of bovine pancreatic trypsin inhibitor

An analogue of the bovine pancreatic trypsin inhibitor (BPTI) folding intermediate that contains only the disulphide bond between Cys5 and Cys55 has been prepared in Escherichia coli by protein engineering methods, with the other four Cys residues replaced by Ser. Two-dimensional 1H nuclear magnetic resonance studies of the analogue have resulted in essentially complete resonance assignments of the folded form of the protein. The folded protein has a compact conformation that is structurally very similar to that of native BPTI, although there are subtle differences and the folded conformation is not very stable. Approximately half of the protein molecules are unfolded at 3 degrees C, and this proportion increases at higher temperatures. The folded and unfolded conformations are in slow exchange. The conformational properties of the analogue can explain many aspects of the kinetic role that the normal (5-55) intermediate plays in the folding of BPTI.     

Beta-sheet folding of fragment (16-36) of bovine pancreatic trypsin inhibitor as predicted by Monte Carlo simulated annealing.

A tertiary structure prediction is described using Monte Carlo simulated annealing for the peptide fragment corresponding to residues 16-36 of bovine pancreatic trypsin inhibitor (BPTI). The simulation starts with randomly chosen initial conformations and is performed without imposing experimental constraints using energy functions given for generic interatomic interactions. Out of 20 simulation trials, seven conformations show a sheet-like structure--two strands connected by a turn--although this sheet-like structure is not as rigid as that observed in native BPTI. It is also shown that these conformations are mostly looped and exhibit a native-like right-handed twist. Unlike the case with the C-peptide of RNase A, no conspicuous alpha-helical structure is found in any of the final conformations obtained in the simulation. However, the lowest-energy conformation does not resemble exactly the native structure. This indicates that the rigid beta-sheet conformation of native BPTI merely corresponds to a local minimum of the energy function if the fragment with residues 16-36 is isolated from the native protein. A statistical analysis of all 20 final conformations suggests that the tendency for the peptide segments to form extended beta-strands is strong for those with residues 18-24, and moderate for those with residues 30-35. The segment of residues 25-29 does not tend to form any definite structure. In native BPTI, the former segments are involved in the beta-sheet and the latter in the turn. A folding scenario is also speculated from this analysis.     

Kinetics of protein folding: Nucleation mechanism,time scales,and pathways

The kinetics and thermodynamics of protein folding is investigated using low friction Langevin simulation of minimal continuum mode of proteins. We show that the model protein has two characteristic temperatures: (a) T_θ, at which the chain undergoes a collapse transition from an extended conformation; (b) T_f(< T_θ), at which a finite size first-order transition to the folded state takes place. The kinetics of approach to the native state from initially denatured conformations is probed by several novel correlation functions. We find that the overall kinetics of approach to the native conformation occurs via a three-stage multiple pathway mechanism. The initial stage, characterized by a series of local dihedral angle transitions, eventually results in the compaction of the protein. Subsequently, the molecule acquires native-like structures during the second stage of folding. The final stage of folding involves activated transitions from one of the native-like structures to the native conformation. The first two stages are characterized by a multiplicity of pathways while relatively few paths are involved in the final stage. A detailed analysis of the dynamics of individual trajectories reveals a novel picture of protein folding. We find that afraction of the initial population reaches the native conformation without the formation of any detectable intermediates. This pathway is associated with a nucleation mechanism, i.e., once a critical number of tertiary contacts are established then the native state is reached rapidly. The remaining fraction of molecules become trapped in misfolded structures (stabilized by incorrect tertiary contacts). The slow folding involves transitions over barriers from these structures to the native conformation. The theoretical predictions are compared with recent experiments that probe protein folding kinetics by hydrogen exchange labeling technique. © 1995 John Wiley & Sons, Inc.     

Understanding the determinants of stability and folding of small globular proteins from their energetics

The results of minimal model calculations indicate that the stability and the kinetic accessibility of the native state of small globular proteins are controlled by few "hot" sites. By means of molecular dynamics simulations around the native conformation, which describe the protein and the surrounding solvent at the all-atom level, an accurate and compact energetic map of the native state of the protein is generated. This map is further simplified by means of an eigenvalue decomposition. The components of the eigenvector associated with the lowest eigenvalue indicate which hot sites are likely to be responsible for the stability and for the rapid folding of the protein. The comparison of the results of the model with the findings of mutagenesis experiments performed for four small proteins show that the eigenvalue decomposition method is able to identify between 60% and 80% of these (hot) sites.     

Structural basis of hierarchical multiple substates of a protein. I: Introduction

A computer experiment of protein dynamics is carried out, which consists of two steps: (1) A Monte Carlo simulation of thermal fluctuations in the native state of a globular protein, bovine pancreatic trypsin inhibitor; and (2) a simulation of the quick freezing of fluctuating conformations into energy minima by minimization of the energy of a number of conformations sampled in the Monte Carlo simulation. From the analysis of results of the computer experiment is obtained the following picture of protein dynamics: multiple energy minima exist in the native state, and they are distributed in clusters in the conformational space. The dynamics has a hierarchical structure which has at least two levels. In the first level, dynamics is restricted within one of the clusters of minima. In the second, transitions occur among the clusters. Local parts of a protein molecule, side chains and local main chain segments, can take multiple locally stable conformations in the native state. Many minima result from combinations of these multiple local conformations. The hierarchical structure in the dynamics comes from interactions among the local parts. Protein molecules have two types of flexibility, each associated with elastic and plastic deformations, respectively.     

Folding Trp-cage to NMR resolution native structure using a coarse-grained protein model

We develop a coarse-grained protein model with a simplified amino acid interaction potential. Using this model, we perform discrete molecular dynamics folding simulations of a small 20-residue protein--Trp-cage--from a fully extended conformation. We demonstrate the ability of the Trp-cage model to consistently reach conformations within 2-angstroms backbone root-mean-square distance from the corresponding NMR structures. The minimum root-mean-square distance of Trp-cage conformations in simulations can be <1 angstroms. Our findings suggest that, at least in the case of Trp-cage, a detailed all-atom protein model with a molecular mechanics force field is not necessary to reach the native state of a protein. Our results also suggest that the success of folding Trp-cage in our simulations and in the reported all-atom molecular mechanics simulation studies may be mainly due to the special stabilizing features specific to this miniprotein.     

Multiple interactions of a DNA-binding protein in vivo. II. Effects of host mutations on DNA replication of phage T4 gene 32 mutants.

A quantitative analysis has been made of the kinetics of disulphide bond formation, breakage, and rearrangement which occur during the folding and unfolding of the pancreatic trypsin inhibitor. The results have been used to infer the energetics of the protein conformational transitions which accompany each step.The folding transition is shown to be a co-operative process in which all intermediate states with one or two disulphide bonds are unstable relative to the unfolded, reduced protein and that in the fully folded conformation with three disulphide bonds. The approximate two-state nature of the transition at equilibrium was demonstrated experimentally. The folding transition of the inhibitor which has been determined kinetically is therefore analogous to that observed generally with other globular proteins.     

Prediction of protein folding from amino acid sequence over discrete conformation spaces

Predicting the three-dimensional structure of a protein given only its amino acid sequence is a long-standing goal in computational chemistry. In the thermodynamic approach, one needs a potential function of conformation that resembles the free energy of the real protein to the extent that the global minimum of the potential is attained by the native conformation and no other. In practice, this has never been achieved with certainty because even with greatly simplified representations of the polypeptide chain, there are an astronomical number of local minima to examine. If one chooses instead a protein representation with only a large but manageable number of discrete conformations, then the global preference of the potential for the native can be directly verified. Representing a protein as a walk on a two-dimensional square lattice makes it easy to see that simple functions of the interresidue contacts are sufficient to globally favor a given "native" conformation, as long as it is a compact, globular structure. Explicit representation of the solvent is not required. Another more realistic way to confine the conformational search to a finite set is to draw alternative conformations from fragments of larger proteins having known crystal structure. Then it is possible to construct a simple function of interresidue contacts in three dimensions such that only 8 proteins are required to determine the adjustable parameters, and the native conformations of 37 other proteins are correctly preferred over all alternative conformations. The deduced function favors short-range backbone-backbone contacts regardless of residue type and long-range hydrophobic associations. Interactions over long distances, such as electrostatics, are not required.     

Application of distance geometry to protein tertiary structure calculations

A general approach to the problem of molecular conformation is advanced. We describe a formalism that permits experimental and theoretical information to be incorporated into a set of upper and lower bounds on intramolecular distances. Structures (conformations) meeting these bounds can be readily generated and compared with each other. To illustrate the use of the method, we have employed a simple “firehose” model for protein folding to predict the long-range hydrophobic interactions in a small protein: pancreatic trypsin inhibitor. Models of this type lead to the proper hairpin turns and a reasonable set of long-range contacts for this protein. Application of the distance geometry method then yields backbone conformations with errors of 4–8 Å compared to the native structure. We discuss both the merits and shortcomings of the firehose model and the relation between distance geometry and energy minimization techniques.     

Denaturants can accelerate folding rates in a class of globular proteins.

We present a lattice Monte Carlo study to examine the effect of denaturants on the folding rates of simplified models of proteins. The two-dimensional model is made from a three-letter code mimicking the presence of hydrophobic, hydrophilic, and cysteine residues. We show that the rate of folding is maximum when the effective hydrophobic interaction epsilon H is approximately equal to the free energy gain epsilon S upon forming disulfide bonds. In the range 1 < or = epsilon H/ epsilon S < or = 3, multiple paths that connect several intermediates to the native state lead to fast folding. It is shown that at a fixed temperature and epsilon S the folding rate increases as epsilon H decreases. An approximate model is used to show that epsilon H should decrease as a function of the concentration of denaturants such as urea or guanidine hydrochloride. Our simulation results, in conjunction with this model, are used to show that increasing the concentration of denaturants can lead to an increase in folding rates. This occurs because denaturants can destabilize the intermediates without significantly altering the energy of the native conformation. Our findings are compared with experiments on the effects of denaturants on the refolding of bovine pancreatic trypsin inhibitor and ribonuclease T1. We also argue that the phenomenon of denaturant-enhanced folding of proteins should be general.