The role of geometric complementarity in secondary structure packing: a systematic docking study

A strong similarity between the major aspects of protein folding and protein recognition is one of the emerging fundamental principles in protein science. A crucial importance of steric complementarity in protein recognition is a well-established fact. The goal of this study was to assess the importance of the steric complementarity in protein folding, namely, in the packing of the secondary structure elements. Although the tight packing of protein structures, in general, is a well-known fact, a systematic study of the role of geometric complementarity in the packing of secondary structure elements has been lacking. To assess the role of the steric complementarity, we used a docking procedure to recreate the crystallographically determined packing of secondary structure elements in known protein structures by using the geometric match only. The docking results revealed a significant percentage of correctly predicted packing configurations. Different types of pairs of secondary structure elements showed different degrees of steric complementarity (from high to low: beta-beta, loop-loop, alpha-alpha, and alpha-beta). Interestingly, the relative contribution of the steric match in different types of pairs was correlated with the number of such pairs in known protein structures. This effect may indicate an evolutionary pressure to select tightly packed elements of secondary structure to maximize the packing of the entire structure. The overall conclusion is that the steric match plays an essential role in the packing of secondary structure elements. The results are important for better understanding of principles of protein structure and may facilitate development of better methods for protein structure prediction.     

Topology is the Principal Determinant in the Folding of a Complex All-alpha Greek Key Death Domain from Human FADD

In order to elucidate the relative importance of secondary structure and topology in determining folding mechanism, we have carried out a phi-value analysis of the death domain (DD) from human FADD. FADD DD is a 100 amino acid domain consisting of six anti-parallel alpha helices arranged in a Greek key structure. We asked how does the folding of this domain compare with that of (a) other all-alpha-helical proteins and (b) other Greek key proteins? Is the folding pathway determined mainly by secondary structure or is topology the principal determinant? Our Φ-value analysis reveals a striking resemblance to the all-beta Greek key immunoglobulin-like domains. Both fold via diffuse transition states and, importantly, long-range interactions between the four central elements of secondary structure are established in the transition state. The elements of secondary structure that are less tightly associated with the central core are less well packed in both cases. Topology appears to be the dominant factor in determining the pathway of folding in all Greek key domains.     

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)     

Multi-State Proteins: Approach Allowing Experimental Determination of the Formation Order of Structure Elements in the Green Fluorescent Protein

The most complex problem in studying multi-state protein folding is the determination of the sequence of formation of protein intermediate states. A far more complex issue is to determine at what stages of protein folding its various parts (secondary structure elements) develop. The structure and properties of different intermediate states depend in particular on these parts. An experimental approach, named μ-analysis, which allows understanding the order of formation of structural elements upon folding of a multi-state protein was used in this study. In this approach the same elements of the protein secondary structure are “tested” by substitutions of single hydrophobic amino acids and by incorporation of cysteine bridges. Single substitutions of hydrophobic amino acids contribute to yielding information on the late stages of protein folding while incorporation of ss-bridges allows obtaining data on the initial stages of folding. As a result of such an μ-analysis, we have determined the order of formation of beta-hairpins upon folding of the green fluorescent protein.     

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.     

Correlation between the structural stability and aggregation propensity of proteins

Protein aggregation, being an outcome of improper protein folding, is largely dependent on the folding kinetics of a protein. Previous studies have reported a positive correlation between the stability of the secondary structural elements of a protein and their rate of folding/unfolding. In this in silico study, the secondary and tertiary structures of proteins a) that form inclusion bodies on overexpression in Escherichia coli, b) that form amyloid fibrils and c) that are soluble on overexpression in E. coli are analyzed for certain features that are known to be associated with structural stability. The study revealed that the soluble proteins seem to have a higher rate of folding (based on contact order) and a lower percentage of exposed hydrophobic residues as compared to the inclusion body forming or amyloidogenic proteins. The soluble proteins also seem to have a more favored helix and strand composition (based on the known secondary structural propensities of amino acids). The secondary structure analyses also reveal that the evolutionary pressure is directed against protein aggregation. This understanding of the positive correlation between structural stability and solubility, along with the other parameters known to influence aggregation, could be exploited in the design of mutations aimed at reducing the aggregation propensity of the proteins.     

Intermediates and the folding of proteins L and G

We use a minimalist protein model, in combination with a sequence design strategy, to determine differences in primary structure for proteins L and G, which are responsible for the two proteins folding through distinctly different folding mechanisms. We find that the folding of proteins L and G are consistent with a nucleation-condensation mechanism, each of which is described as helix-assisted beta-1 and beta-2 hairpin formation, respectively. We determine that the model for protein G exhibits an early intermediate that precedes the rate-limiting barrier of folding, and which draws together misaligned secondary structure elements that are stabilized by hydrophobic core contacts involving the third beta-strand, and presages the later transition state in which the correct strand alignment of these same secondary structure elements is restored. Finally, the validity of the targeted intermediate ensemble for protein G was analyzed by fitting the kinetic data to a two-step first-order reversible reaction, proving that protein G folding involves an on-pathway early intermediate, and should be populated and therefore observable by experiment.     

Gatekeepers in the ribosomal protein s6: thermodynamics, kinetics, and folding pathways revealed by a minimalist protein model

We investigate the effect of structural gatekeepers on the folding of the ribosomal protein S6. Folding thermodynamics and early refolding kinetics are studied for this system utilizing computer simulations of a minimalist protein model. When gatekeepers are eliminated, the thermodynamic signature of a folding intermediate emerges, and a marked decrease in folding efficiency is observed. We explain the prerequisites that determine the "strength" of a given gatekeeper. The investigated gatekeepers are found to have distinct functions, and to guide the folding and time-dependent packing of non-overlapping secondary structure elements in the protein. Gatekeepers avoid kinetic traps during folding by favoring the formation of "productive topologies" on the way to the native state. The trends in folding rates in the presence/absence of gatekeepers observed for our minimalist model of S6 are in very good agreement with experimental data on this protein.     

Toward minimalist models of larger proteins: a ubiquitin-like protein

Our recently developed off-lattice bead model capable of simulating protein structures with mixed alpha/beta content has been extended to model the folding of a ubiquitin-like protein and provides a means for examining the more complex kinetics involved in the folding of larger proteins. Using trajectories generated from constant-temperature Langevin dynamics simulations and sampling with the multiple multi-histogram method over five-order parameters, we are able to characterize the free energy landscape for folding and find evidence for folding through compact intermediates. Our model reproduces the observation that the C-terminus loop structure in ubiquitin is the last to fold in the folding process and most likely plays a spectator role in the folding kinetics. The possibility of a productive metastable intermediate along the folding pathway consisting of collapsed states with no secondary structure, and of intermediates or transition structures involving secondary structural elements occurring early in the sequence, is also supported by our model. The kinetics of folding remain multi-exponential below the folding temperature, with glass-like kinetics appearing at T/T(f) approximately 0.86. This new physicochemical model, designed to be predictive, helps validate the value of modeling protein folding at this level of detail for genomic-scale studies, and motivates further studies of other protein topologies and the impact of more complex energy functions, such as the addition of solvation forces.     

The "first in-last out" hypothesis on protein folding revisited

We calculated profiles for mean residue depth, contact order, and number of contacts in the native structure of a series of proteins for which folding has been studied extensively, the chymotrypsin inhibitor 2, the SH3 module from the src tyrosine kinase, the small ribonuclease barnase, the bacterial immunity protein Im7, and apomyoglobin. We compared these profiles with experimental data from equilibrium or pulse labeling hydrogen-deuterium exchange obtained from NMR and phi values obtained from the protein engineering approach. We find a good qualitative agreement between the hierarchy of formation of topological elements during the folding process and the ranking of secondary structure elements in terms of residue depth. Residues that are most deeply buried in the core of the native protein usually belong to stretches of secondary structure elements that are formed early in the folding pathway. Residue depth can thus provide a useful and simple tool for the design of folding experiments.     

Does trifluoroethanol affect folding pathways and can it be used as a probe of structure in transition states?

Nonaqueous co-solvents, particularly 2,2,2-trifluoroethanol (TFE), have been used as tools to study protein folding. By analyzing FKBP12, an alpha/beta-protein that folds with two-state kinetics, we have been able to address three key questions concerning the use of TFE. First, does TFE perturb the folding pathway? Second, can the observed changes in the rate of folding and unfolding in TFE be attributed to a change in free energy of a single state? Finally, can TFE be used to infer information on secondary structure formation in the transition state? Protein engineering experiments on FKBP12, coupled with folding and unfolding experiments in 0% and 9.6% TFE, conclusively show that TFE does not perturb the folding pathway of this protein. Our results also suggest that the changes in folding and unfolding rates observed in 9.6% TFE are due to a global effect of TFE on the protein, rather than the stabilization of any elements of secondary structure in the transition state. Thus, studies with TFE and other co-solvents can be accurately interpreted only when combined with other techniques.     

Role of loops in the folding and stability of yeast phosphoglycerate kinase.

Yeast phosphoglycerate kinase (yPGK) is a monomeric two domain protein used as folding model representative of large proteins. We inserted short unstructured sequences (four Gly or four Thr) into the connections between secondary structure elements and studied the consequences of these insertions on the folding process and stability of yPGK. All the mutated proteins can refold efficiently. The effect per residue on stability is larger for the first inserted residue. Insertion in two long betaalpha loops (at residue positions 71 and 129) is more destabilizing than an insertion in a short alphabeta loop (at residue position 89) located on the opposite side of the N-terminal domain. The effect on stability is mainly due to a large increase of the unfolding rate rather than a decrease of the folding rate. This suggests that these connections between secondary structure elements do not play an active role in directing the folding process. Insertion into the short alphabeta loop (position 89) has limited effects on stability and results in the detection of a kinetic phase not previously seen with the wild-type protein, suggesting that insertions in this particular loop do qualitatively affect the folding process without a large effect on folding efficiency. For the two long betaalpha loops (positions 71 and 129) located in the inner surface of the N-terminal domain, the effects on stability are possibly associated with decoupling of the two domains as observed by differential scanning calorimetry during thermal unfolding.     

Predicting <Emphasis Type="Italic">β</Emphasis>-turns and their types using predicted backbone dihedral angles and secondary structures

Background  

β-turns are secondary structure elements usually classified as coil. Their prediction is important, because of their role in protein folding and their frequent occurrence in protein chains.     

Testing the relationship between foldability and the early folding events of dihydrofolate reductase from Escherichia coli

A "folding element" is a contiguous peptide segment crucial for a protein to be foldable and is a new concept that could assist in our understanding of the protein-folding problem. It is known that the presence of the complete set of folding elements of dihydrofolate reductase (DHFR) from Escherichia coli is essential for the protein to be foldable. Since almost all of the amino acid residues known to be involved in the early folding events of DHFR are located within the folding elements, a close relationship between the folding elements and early folding events is hypothesized. In order to test this hypothesis, we have investigated whether or not the early folding events are preserved in circular permutants and topological mutants of DHFR, in which the order of the folding elements is changed but the complete set of folding elements is present. The stopped-flow circular dichroism (CD) measurements show that the CD spectra at the early stages of folding are similar among the mutants and the wild-type DHFR, indicating that the presence of the complete set of folding elements is sufficient to preserve the early folding events. We have further examined whether or not sequence perturbation on the folding elements by a single amino acid substitution affects the early folding events of DHFR. The results show that the amino acid substitutions inside of the folding elements can affect the burst-phase CD spectra, whereas the substitutions outside do not. Taken together, these results indicate that the above hypothesis is true, suggesting a close relationship between the foldability of a protein and the early folding events. We propose that the folding elements interact with each other and coalesce to form a productive intermediate(s) early in the folding, and these early folding events are important for a protein to be foldable.     

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.     

Molecular dynamics study of the folding of hydrophobin SC3 at a hydrophilic/hydrophobic interface

Hydrophobins are fungal proteins that self-assemble at hydrophilic/hydrophobic interfaces into amphipathic membranes. These assemblages are extremely stable and posses the remarkable ability to invert the polarity of the surface on which they are adsorbed. Neither the three-dimensional structure of a hydrophobin nor the mechanism by which they function is known. Nevertheless, there are experimental indications that the self-assembled form of the hydrophobins SC3 and EAS at a water/air interface is rich with beta-sheet secondary structure. In this paper we report results from molecular dynamics simulations, showing that fully extended SC3 undergoes fast (approximately 100 ns) folding at a water/hexane interface to an elongated planar structure with extensive beta-sheet secondary elements. Simulations in each of the bulk solvents result in a mainly unstructured globular protein. The dramatic enhancement in secondary structure, whether kinetic or thermodynamic in origin, highlights the role interfaces between phases with large differences in polarity can have on folding. The partitioning of the residue side-chains to one of the two phases can serve as a strong driving force to initiate secondary structure formation. The interactions of the side-chains with the environment at an interface can also stabilize configurations that otherwise would not occur in a homogenous solution.     

mRNA环结构对蛋白质折叠速率的影响

前期的相关研究发现mRNA二级结构中存在对蛋白质折叠速率的重要影响因素.而mRNA二级结构中普遍存在着各种复杂的环结构,这些环结构是否对蛋白质折叠速率也有重要的影响呢?不同的环结构对蛋白质折叠速率的影响是否相同呢?基于此想法,建立了一个包含mRNA内部环、发夹环、膨胀环和多分支环等环结构信息和相应蛋白质折叠速率的数据库.对于数据库中的每一个蛋白质,计算了mRNA二级结构中各种环结构碱基含量、配对碱基含量及单链碱基含量等参量,分析了各参量与相应蛋白质折叠速率的相关性.结果显示,各种环结构碱基含量与蛋白质折叠速率均呈极显著或显著正相关.说明mRNA环结构对蛋白质折叠速率有重要的影响.进一步,把蛋白质按照不同折叠类型或不同二级结构类型分组后,对每一组蛋白质重复上述的分析工作.结果表明,对不同类蛋白质,mRNA的各种环结构对其相应蛋白质折叠速率的影响存在着显著差异.上述研究将为进一步开展有关mRNA和蛋白质折叠速率的研究奠定理论基础.     

Residue-Specific Analysis of Frustration in the Folding Landscape of Repeat β/α Protein Apoflavodoxin

Flavodoxin adopts the common repeat β/α topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central β-strand 3 (i.e., at β-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, ?-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the β-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central β-sheet directs initial misfolding, while proper alignment on both sides of β-strand 3 is necessary for productive folding.     

A rewired green fluorescent protein: folding and function in a nonsequential, noncircular GFP permutant

The sequential order of secondary structural elements in proteins affects the folding and activity to an unknown extent. To test the dependence on sequential connectivity, we reconnected secondary structural elements by their solvent-exposed ends, permuting their sequential order, called "rewiring". This new protein design strategy changes the topology of the backbone without changing the core side chain packing arrangement. While circular and noncircular permutations have been observed in protein structures that are not related by sequence homology, to date no one has attempted to rationally design and construct a protein with a sequence that is noncircularly permuted while conserving three-dimensional structure. Herein, we show that green fluorescent protein can be rewired, still functionally fold, and exhibit wild-type fluorescence excitation and emission spectra.     

Matching simulation and experiment: a new simplified model for simulating protein folding.

Simulations of simplified protein folding models have provided much insight into solving the protein folding problem. We propose here a new off-lattice bead model, capable of simulating several different fold classes of small proteins. We present the sequence for an alpha/beta protein resembling the IgG-binding proteins L and G. The thermodynamics of the folding process for this model are characterized using the multiple multihistogram method combined with constant-temperature Langevin simulations. The folding is shown to be highly cooperative, with chain collapse nearly accompanying folding. Two parallel folding pathways are shown to exist on the folding free energy landscape. One pathway contains an intermediate--similar to experiments on protein G, and one pathway contains no intermediates-similar to experiments on protein L. The folding kinetics are characterized by tabulating mean-first passage times, and we show that the onset of glasslike kinetics occurs at much lower temperatures than the folding temperature. This model is expected to be useful in many future contexts: investigating questions of the role of local versus nonlocal interactions in various fold classes, addressing the effect of sequence mutations affecting secondary structure propensities, and providing a computationally feasible model for studying the role of solvation forces in protein folding.