Computational design of proteins stereochemically optimized in size, stability, and folding speed

Artificial proteins potentially barrier-free in the folding kinetics are approached computationally under the guidance of protein-folding theories. The smallest and fastest folding globular protein triple-helix-bundle (THB) is so modified as to minimize or eliminate its presumed barriers in folding speed. As the barriers may reside in the ordering of either secondary or tertiary structure, the elements of both secondary and tertiary structure in the protein are targeted for prenucleation with suitable stereochemically constrained amino acid residues. The required elements of topology and sequence for the THB are optimized independently; first the topology is optimized with simulated annealing in polypeptides of highly simplified alphabet; next, the sequence in side chains is optimized using the standard inverse design methods. The resultant three best-adapted THBs, variable in topology and distinctive in sequences, are assessed by comparing them with a few benchmark proteins. The results of mainly molecular dynamics (MD) comparisons, undertaken in explicit water at different temperatures, show that the designed sequences are favorably placed against the chosen benchmarks as THB proteins potentially thermostable in the native folds. Folding simulation experiments with MD establish that the designed sequences are rapid in the folding of individual helices, but not in the evolution of tertiary structure; energetic cum topological frustrations remain but could be the artifacts of the starting conformations that were chosen in the THBs in the folding simulations. Overall, a practical high-throughput approach for de novo protein design has been developed that may have fruitful application for any type of tertiary structure.     

Mechanism of protein folding

The high structural resolution of the main transition states for the formation of native structure for the six small proteins of which Phi-values for a large set of mutants have become available, barstar, barnase, chymotrypsin inhibitor 2, Arc repressor, the src SH3 domain, and a tetrameric p53 domain reveals that for the first 5 of these proteins: (1) Residues that belong to regular secondary structure have a significantly larger average fraction of native structural consolidation than residues in loops; (2) on the other hand, secondary and tertiary structures have built up to the same degree, or at least a high degree, but nonuniformly distributed over the molecule; (3) the most consolidated parts of each protein molecule in the transition state cluster together, and these clusters contain a significantly higher percentage of residues that belong to regular secondary structure than the rest of the molecule. These observations further reconcile the framework model with the nucleation-condensation mechanism for folding: The amazing speed of protein folding can be understood as caused by the catalytic effect of the formation of clusters of residues which have particularly high preferences for the early formation of regular secondary structure in the presence of significant amounts of tertiary structure interactions.     

Systematic circular permutation of an entire protein reveals essential folding elements

The importance of chain connectivity in determining protein function and stability can be examined by breaking the peptide backbone using a technique such as circular permutation. Cleavage at certain positions results in a complete loss of the ability of the protein to fold. When such cleavage sites occur sequentially in the primary structure, we call the region a 'folding element', a new concept that could assist in our understanding of the protein folding problem. The folding elements of dihydrofolate reductase have been assigned by conducting a systematic circular permutation analysis in which the peptide backbone was sequentially broken between every pair of residues in the protein. The positions of folding elements do not appear to correspond to secondary structure motifs, substrate or coenzyme binding sites, or accessible surface area. However, almost all of the amino acid residues known to be involved in early folding events are located within the folding elements.     

The folding of an enzyme. VI. The folding pathway of barnase: comparison with theoretical models.

The sequence of events in the refolding pathway of barnase has been analysed to search for general principles in protein folding. There appears to be a correlation between burying hydrophobic surface area and early folding events. All the regions that fold early interact extensively with the beta-sheet. These interactions involve predominantly hydrophobic interactions and the burial of very extensive hydrophobic areas in which multiple, close, hydrophobic-hydrophobic contacts are established around a central group of aliphatic residues. There is no burial of hydrophilic residues in these regions; those that are partly screened from the solvent make hydrogen bonds. All the regions or interactions that are made late in the folding pathway do not make extensive contacts with the beta-sheet. Their buried hydrophobic regions lack a central hydrophobic residue or residues around which other hydrophobic residues pack. Further, in some of these regions there is an extensive burial of hydrophilic residues. The results are consistent with one of the earlier events in protein folding being the local formation of native-like secondary structure elements driven by local hydrophobic surface burial. A possible candidate for an initiation site is a beta-hairpin between beta-strands 3 and 4 that is conserved in the microbial ribonuclease family. A comparison of structures in this family shows that those regions that can be superimposed, or have sequence homology, correspond to elements of structure that are formed and interact with each other early in the folding pathway, suggesting that some of these residues could be involved in directing the folding process. The data on barnase combined with results from other laboratories suggest the following tentative conclusions for the refolding of small monomeric proteins. (1) The refolding pathway is, at least in part, sequential and of compulsory order. (2) Secondary structure formation is driven by local hydrophobic surface burial and precedes the formation of most tertiary interactions. These elements are then stabilized and sometimes elongated by tertiary interactions. It is plausible that there are stop signals encoded in the linear sequence that prevent the elongation of isolated secondary structure elements in solution to a larger extent than is found in the folded protein. (3) Many tertiary interactions are not very constrained in the intermediate but become more and more defined as the hydrophobic cores consolidate, loop structures form and the configuration of surface residues takes place. The interactions between different elements of secondary structure are the last ones to be consolidated while the interactions within the secondary structure elements are consolidated earlier.(ABSTRACT TRUNCATED AT 400 WORDS)     

Using motion planning to study protein folding pathways.

We present a framework for studying protein folding pathways and potential landscapes which is based on techniques recently developed in the robotics motion planning community. Our focus in this work is to study the protein folding mechanism assuming we know the native fold. That is, instead of performing fold prediction, we aim to study issues related to the folding process, such as the formation of secondary and tertiary structure, and the dependence of the folding pathway on the initial denatured conformation. Our work uses probabilistic roadmap (PRM) motion planning techniques which have proven successful for problems involving high-dimensional configuration spaces. A strength of these methods is their efficiency in rapidly covering the planning space without becoming trapped in local minima. We have applied our PRM technique to several small proteins (~60 residues) and validated the pathways computed by comparing the secondary structure formation order on our paths to known hydrogen exchange experimental results. An advantage of the PRM framework over other simulation methods is that it enables one to easily and efficiently compute folding pathways from any denatured starting state to the (known) native fold. This aspect makes our approach ideal for studying global properties of the protein's potential landscape, most of which are difficult to simulate and study with other methods. For example, in the proteins we study, the folding pathways starting from different denatured states sometimes share common portions when they are close to the native fold, and moreover, the formation order of the secondary structure appears largely independent of the starting denatured conformation. Another feature of our technique is that the distribution of the sampled conformations is correlated with the formation of secondary structure and, in particular, appears to differentiate situations in which secondary structure clearly forms first and those in which the tertiary structure is obtained more directly. Overall, our results applying PRM techniques are very encouraging and indicate the promise of our approach for studying proteins for which experimental results are not available.     

Identification and ab initio simulations of early folding units in proteins

The location of protein subunits that form early during folding, constituted of consecutive secondary structure elements with some intrinsic stability and favorable tertiary interactions, is predicted using a combination of threading algorithms and local structure prediction methods. Two folding units are selected among the candidates identified in a database of known protein structures: the fragment 15-55 of 434 cro, an all-alpha protein, and the fragment 1-35 of ubiquitin, an alpha/beta protein. These units are further analyzed by means of Monte Carlo simulated annealing using several database-derived potentials describing different types of interactions. Our results suggest that the local interactions along the chain dominate in the first folding steps of both fragments, and that the formation of some of the secondary structures necessarily occurs before structure compaction. These findings led us to define a prediction protocol, which is efficient to improve the accuracy of the predicted structures. It involves a first simulation with a local interaction potential only, whose final conformation is used as a starting structure of a second simulation that uses a combination of local interaction and distance potentials. The root mean square deviations between the coordinates of predicted and native structures are as low as 2-4 A in most trials. The possibility of extending this protocol to the prediction of full proteins is discussed. Proteins 2001;42:164-176.     

Investigation of the folding pathway of the TEM-1 β-lactamase

The TEM-1 β-lactamase is a globular protein containing 12 proline residues. The folding mechanism of this enzyme was investigated by kinetic and equilibrium experiments with the help of fluorescence spectroscopy and circular dichroism. The equilibrium denaturation of the protein induced by guanidine hydrochloride occurs in two discrete steps, indicating the existence of a thermodynamically stable intermediate state. Thisstate is 5.2 ± 0.4 kcal/mol less stable than the native conformation and 5.7 ± 0.2 kcal/mol more stable than the fully denaturedprotein. This intermediate state exhibits a high content of native secondary structure elements but is devoid of specific tertiary organization; its relation to the “molten globule” is discussed. Refolding kinetic experimentsrevealed the existence of a transient intermediate conformation between thethermodynamically stable intermediate and the native protein. This transient intermediate appears rapidly during the folding reaction. It exhibits a secondary structure content very similar to that of the native protein and has also recovered a significant amount of tertiary organisation. The final refolding step of the TEM-1 β-lactamase, leading to the native enzyme, is dominated by two major slow kinetic phases which probablyreflect a very complex process kinetically limited by proline cis/transisomerization. © 1995 Wiley-Liss, Inc.     

Secondary structure prediction and folding of globular protein: refolding of ferredoxin

The physicochemical mechanism of protein folding has been elucidated by the island model, describing a growth type of folding. The folding pathway is closely related with nucleation on the polypeptide chain and thus the formation of small local structures or secondary structures at the earliest stage of folding is essential to all following steps. The island model is applicable to any protein, but a high precision of secondary structure prediction is indispensable to folding simulation. The secondary structures formed at the earliest stage of folding are supposed to be of standard form, but they are usually deformed during the folding process, especially at the last stage, although the degree of deformation is different for each protein. Ferredoxin is an example of a protein having this property. According to X-ray investigation (1FDX), ferredoxin is not supposed to have secondary structures. However, if we assumed that in ferredoxin all the residues are in a coil state, we could not attain the correct structure similar to the native one. Further, we found that some parts of the chain are not flexible, suggesting the presence of secondary structures, in agreement with the recent PDB data (1DUR). Assuming standard secondary structures (-helices and -strands) at the nonflexible parts at the early stage of folding, and deforming these at the final stage, a structure similar to the native one was obtained. Another peculiarity of ferredoxin is the absence of disulfide bonds, in spite of its having eight cysteines. The reason cysteines do not form disulfide bonds became clear by applying the lampshade criterion, but more importantly, the two groups of cysteines are ready to make iron complexes, respectively, at a rather later stage of folding. The reason for poor prediction accuracy of secondary structure with conventional methods is discussed.     

The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H

Elucidation of the high-resolution structures of folding intermediates is a necessary but difficult step toward the ultimate understanding of the mechanism of protein folding. Here, using hydrogen-exchange-directed protein engineering, we populated the folding intermediate of the Thermus thermophilus ribonuclease H, which forms before the rate-limiting transition state, by removing the unfolded regions of the intermediate, including an α-helix and two β-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 Å). It has a native-like backbone topology and shows some local deviations from the native structure, revealing that the structure of the folded region of an early folding intermediate can be as well defined as the native structure. The topological parameters calculated from the structures of the intermediate mimic and the native state predict that the intermediate should fold on a millisecond time scale or less and form much faster than the native state. Other factors that may lead to the slow folding of the native state and the accumulation of the intermediate before the rate-limiting transition state are also discussed.     

Probing minimal independent folding units in dihydrofolate reductase by molecular dissection.

Molecular dissection was employed to identify minimal independent folding units in dihydrofolate reductase (DHFR) from Escherichia coli. Eight overlapping fragments of DHFR, spanning the entire sequence and ranging in size from 36 to 123 amino acids, were constructed by chemical cleavage. These fragments were designed to examine the effect of tethering multiple elements of secondary structure on folding and to test if the secondary structural domains represent autonomous folding units. CD and fluorescence spectroscopy demonstrated that six fragments containing up to a total of seven alpha-helices or beta-strands and, in three cases, the adenine binding domain (residues 37-86), are largely disordered. A stoichiometric mixture of the two fragments comprising the large discontinuous domain, 1-36 and 87-159, also showed no evidence for folding beyond that observed for the isolated fragments. A fragment containing residues 1-107 appears to have secondary and tertiary structure; however, spontaneous self-association made it impossible to determine if this structure solely reflects the behavior of the monomeric form. In contrast, a monomeric fragment spanning residues 37-159 possesses significant secondary and tertiary structure. The urea-induced unfolding of fragment 37-159 in the presence of 0.5 M ammonium sulfate was found to be a well-defined, two-state process. The observation that fragment 37-159 can adopt a stable native fold with unique, aromatic side-chain packing is quite striking because residues 1-36 form an integral part of the structural core of the full-length protein.     

Cooperative folding of the isolated alpha-helical domain of hen egg-white lysozyme.

Proteins in the alpha-lactalbumin and c-type lysozyme family have been studied extensively as model systems in protein folding. Early formation of the alpha-helical domain is observed in both alpha-lactalbumin and c-type lysozyme; however, the details of the kinetic folding pathways are significantly different. The major folding intermediate of hen egg-white lysozyme has a cooperatively formed tertiary structure, whereas the intermediate of alpha-lactalbumin exhibits the characteristics of a molten globule. In this study, we have designed and constructed an isolated alpha-helical domain of hen egg-white lysozyme, called Lyso-alpha, as a model of the lysozyme folding intermediate that is stable at equilibrium. Disulfide-exchange studies show that under native conditions, the cysteine residues in Lyso-alpha prefer to form the same set of disulfide bonds as in the alpha-helical domain of full-length lysozyme. Under denaturing conditions, formation of the nearest-neighbor disulfide bonds is strongly preferred. In contrast to the isolated alpha-helical domain of alpha-lactalbumin, Lyso-alpha with two native disulfide bonds exhibits a well-defined tertiary structure, as indicated by cooperative thermal unfolding and a well-dispersed NMR spectrum. Thus, the determinants for formation of the cooperative side-chain interactions are located mainly in the alpha-helical domain. Our studies suggest that the difference in kinetic folding pathways between alpha-lactalbumin and lysozyme can be explained by the difference in packing density between secondary structural elements and support the hypothesis that the structured regions in a protein folding intermediate may correspond to regions that can fold independently.     

Molecular basis of co-operativity in protein folding. III. Structural identification of cooperative folding units and folding intermediates.

The hierarchical partition function formalism for protein folding developed earlier has been extended through the use of three-dimensional polar and apolar contact plots. For each amino acid residue in the protein, these plots indicate the apolar and polar surfaces that are buried from the solvent, the identity of all amino acid residues that contribute to this shielding, and the magnitude of their contributions. These contact plots are then used to examine the distribution of the free energy of stabilization throughout the protein molecule. Analysis of these data allows identification of co-operative folding units and their hierarchical levels, and the identification of partially folded intermediates with a significant probability of being populated. The overall folding/unfolding thermodynamics of 12 globular proteins, for which crystallographic and experimental thermodynamics are available, is predicted within error. An energetic classification of partially folded intermediates is presented and the results compared to those cases for which structural and thermodynamic experimental information is available. Four different types of partially folded states and their structural energies are considered. (1) Local intermediates, in which only a local region of the protein loses secondary and tertiary interactions, while the rest of the protein remains intact. (2) Global intermediates, corresponding to the standard molten globule definition, in which significant secondary structure is maintained but native-like tertiary structure contacts are disrupted. (3) Extended intermediates characterized by the existence of secondary structure elements (e.g. alpha-helices) exposed to solvent. (4) Folding intermediates in proteins with two structural domains. The structure and energetics of folding intermediates of apo-myoglobin, alpha-lactalbumin, phosphoglycerate kinase and arabinose-binding protein are considered in detail.     

Contribution of disulfide bonds to stability of the folding intermediate of alpha-lactalbumin

The secondary structure formed in disulfide reduced alpha-lactalbumin is investigated by CD spectrum and is compared with that of the folding intermediate of the disulfide intact protein. The peptide backbone structure of the reduced protein depends strongly on salt concentration in contrast to that of the intermediate. It is close to a random coil in the absence of salt, but it is almost the same as that of the intermediate at a high concentration of salt. The secondary structures of both the proteins undergo broad unfolding transitions when temperature is raised or when urea is added. The secondary structure of the reduced protein is less stable against both heat and urea. These results show that the disulfide bonds are not a determinant of the secondary structure formed at an early stage of folding, and they stabilize the secondary structure of the folding intermediate.     

Rapid formation of secondary structure framework in protein folding studied by stopped-flow circular dichroism

Kinetic refolding reactions of ferricytochrome c and beta-lactoglobulin have been studied by stopped-flow circular dichroism by monitoring rapid ellipticity changes of peptide backbone and side-chain chromophores. In both proteins, a transient intermediate accumulates within the dead time of stopped-flow mixing (18 ms), and the intermediate has an appreciable amount of secondary structure but possesses an unfolded tertiary structure. It is suggested that the rapid formation of a secondary structure framework in protein folding is a common property observed in a variety of globular proteins.     

Unravelling the folding of bacteriorhodopsin

The folding mechanism of integral membrane proteins has eluded detailed study, largely as a result of the inherent difficulties in folding these proteins in vitro. The seven-transmembrane helical protein bacteriorhodopsin has, however, allowed major advances to be made, not just on the folding of this particular protein, but also on the factors governing folding of transmembrane alpha-helical proteins in general. This review focusses on kinetic and equilibrium studies of bacteriorhodopsin folding in vitro. It covers what is currently known about secondary and tertiary structure formation as well as the events accompanying retinal binding, for protein in detergent and lipid systems, including native membrane samples.     

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.     

Understanding the mechanism of beta-sheet folding from a chemical and biological perspective

Perturbing the structure of the Pin1 WW domain, a 34-residue protein comprised of three beta-strands and two intervening loops has provided significant insight into the structural and energetic basis of beta-sheet folding. We will review our current perspective on how structure acquisition is influenced by the sequence, which determines local conformational propensities and mediates the hydrophobic effect, hydrogen bonding, and analogous intramolecular interactions. We have utilized both traditional site-directed mutagenesis and backbone mutagenesis approaches to alter the primary structure of this beta-sheet protein. Traditional site-directed mutagenesis experiments are excellent for altering side-chain structure, whereas amide-to-ester backbone mutagenesis experiments modify backbone-backbone hydrogen bonding capacity. The transition state structure associated with the folding of the Pin1 WW domain features a partially H-bonded, near-native reverse turn secondary structure in loop 1 that has little influence on thermodynamic stability. The thermodynamic stability of the Pin1 WW domain is largely determined by the formation of a small hydrophobic core and by the formation of desolvated backbone-backbone H-bonds enveloped by this hydrophobic core. Loop 1 engineering to the consensus five-residue beta-bulge-turn found in most WW domains or a four-residue beta-turn found in most beta-hairpins accelerates folding substantially relative to the six-residue turn found in the wild type Pin1 WW domain. Furthermore, the more efficient five- and four-residue reverse turns now contribute to the stability of the three-stranded beta-sheet. These insights have allowed the design of Pin1 WW domains that fold at rates that approach the theoretical speed limit of folding.     

Prelude and Fugue, predicting local protein structure, early folding regions and structural weaknesses

Prelude&Fugue are bioinformatics tools aiming at predicting the local 3D structure of a protein from its amino acid sequence in terms of seven backbone torsion angle domains, using database-derived potentials. Prelude(&Fugue) computes all lowest free energy conformations of a protein or protein region, ranked by increasing energy, and possibly satisfying some interresidue distance constraints specified by the user. (Prelude&)Fugue detects sequence regions whose predicted structure is significantly preferred relative to other conformations in the absence of tertiary interactions. These programs can be used for predicting secondary structure, tertiary structure of short peptides, flickering early folding sequences and peptides that adopt a preferred conformation in solution. They can also be used for detecting structural weaknesses, i.e. sequence regions that are not optimal with respect to the tertiary fold. AVAILABILITY: http://babylone.ulb.ac.be/Prelude_and_Fugue.     

Engineering a compact non-native state of intestinal fatty acid-binding protein

The last three C-terminal residues (129-131) of intestinal fatty acid-binding protein (IFABP) participate in four main-chain hydrogen bonds and two electrostatic interactions to sequentially distant backbone and side-chain atoms. To assess if these interactions are involved in the final adjustment of the tertiary structure during folding, we engineered an IFABP variant truncated at residue 128. An additional mutation, Trp-6-->Phe, was introduced to simplify the conformational analysis by optical methods. Although the changes were limited to a small region of the protein surface, they resulted in an IFABP with altered secondary and tertiary structure. Truncated IFABP retains some cooperativity, is monomeric, highly compact, and has the molecular dimensions and shape of the native protein. Our results indicated that residues 129-131 are part of a crucial conformational determinant in which several long-range interactions, essential for the acquisition of the native state, are established. This work suggests that carefully controlled truncation can populate equilibrium non-native states under physiological conditions. These non-native states hold a great promise as experimental models for protein folding.