A residue in helical conformation in the native state adopts a β-strand conformation in the folding transition state despite its high and canonical Φ-value

Delineating structures of the transition states in protein folding reactions has provided great insight into the mechanisms by which proteins fold. The most common method for obtaining this information is Φ-value analysis, which is carried out by measuring the changes in the folding and unfolding rates caused by single amino acid substitutions at various positions within a given protein. Canonical Φ-values range between 0 and 1, and residues displaying high values within this range are interpreted to be important in stabilizing the transition state structure, and to elicit this stabilization through native-like interactions. Although very successful in defining the general features of transition state structures, Φ-value analysis can be confounded when non-native interactions stabilize this state. In addition, direct information on backbone conformation within the transition state is not provided. In the work described here, we have investigated structure formation at a conserved β-bulge (with helical conformation) in the Fyn SH3 domain by characterizing the effects of substituting all natural amino acids at one position within this structural motif. By comparing the effects on folding rates of these substitutions with database-derived local structure propensity values, we have determined that this position adopts a non-native backbone conformation in the folding transition state. This result is surprising because this position displays a high and canonical Φ-value of 0.7. This work emphasizes the potential role of non-native conformations in folding pathways and demonstrates that even positions displaying high and canonical Φ-values may, nevertheless, adopt a non-native conformation in the transition state.     

Local interactions drive the formation of nonnative structure in the denatured state of human alpha-lactalbumin: a high resolution structural characterization of a peptide model in aqueous solution.

There are a small number of peptides derived from proteins that have a propensity to adopt structure in aqueous solution which is similar to the structure they possess in the parent protein. There are far fewer examples of protein fragments which adopt stable nonnative structures in isolation. Understanding how nonnative interactions are involved in protein folding is crucial to our understanding of the topic. Here we show that a small, 11 amino acid peptide corresponding to residues 101-111 of the protein alpha-lactalbumin is remarkably structured in isolation in aqueous solution. The peptide has been characterized by 1H NMR, and 170 ROE-derived constraints were used to calculate a structure. The calculations yielded a single, high-resolution structure for residues 101-107 that is nonnative in both the backbone and side-chain conformations. In the pH 6.5 crystal structure, residues 101-105 are in an irregular turn-like conformation and residues 106-111 form an alpha-helix. In the pH 4.2 crystal structure, residues 101-105 form an alpha-helix, and residues 106-111 form a loopike structure. Both of these structures are significantly different from the conformation adopted by our peptide. The structure in the peptide model is primarily the result of local side-chain interactions that force the backbone to adopt a nonnative 310/turn-like structure in residues 103-106. The structure in aqueous solution was compared to the structure in 30% trifluoroethanol (TFE), and clear differences were observed. In particular, one of the side-chain interactions, a hydrophobic cluster involving residues 101-105, is different in the two solvents and residues 107-111 are considerably more ordered in 30% TFE. The implications of the nonnative structure for the folding of alpha-lactalbumin is discussed.     

Helix, sheet, and polyproline II frequencies and strong nearest neighbor effects in a restricted coil library

A central issue in protein folding is the degree to which each residue's backbone conformational preferences stabilize the native state. We have studied the conformational preferences of each amino acid when the amino acid is not constrained to be in a regular secondary structure. In this large but highly restricted coil library, the backbone preferentially adopts dihedral angles consistent with the polyproline II conformation rather than alpha or beta conformations. The preference for the polyproline II conformation is independent of the degree of solvation. In conjunction with a new masking procedure, the frequencies in our coil library accurately recapitulate both helix and sheet frequencies for the amino acids in structured regions, as well as polyproline II propensities. Therefore, structural propensities for alpha-helices and beta-sheets and for polyproline II conformations in unfolded peptides can be rationalized solely by local effects. In addition, these propensities are often strongly affected by both the chemical nature and the conformation of neighboring residues, contrary to the Flory isolated residue hypothesis.     

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.     

Conformational mapping of the N-terminal segment of surfactant protein B in lipid using 13C-enhanced Fourier transform infrared spectroscopy.

Synthetic peptides based on the N-terminal domain of human surfactant protein B (SP-B1-25; 25 amino acid residues; NH2-FPIPLPYCWLCRALIKRIQAMIPKG) retain important lung activities of the full-length, 79-residue protein. Here, we used physical techniques to examine the secondary conformation of SP-B1-25 in aqueous, lipid and structure-promoting environments. Circular dichroism and conventional, 12C-Fourier transform infrared (FTIR) spectroscopy each indicated a predominate alpha-helical conformation for SP-B1-25 in phosphate-buffered saline, liposomes of 1-palmitoyl-2-oleoyl phosphatidylglycerol and the structure-promoting solvent hexafluoroisopropanol; FTIR spectra also showed significant beta- and random conformations for peptide in these three environments. In further experiments designed to map secondary structure to specific residues, isotope-enhanced FTIR spectroscopy was performed with 1-palmitoyl-2-oleoyl phosphatidylglycerol liposomes and a suite of SP-B1-25 peptides labeled with 13C-carbonyl groups at either single or multiple sites. Combining these 13C-enhanced FTIR results with energy minimizations and molecular simulations indicated the following model for SP-B1-25 in 1-palmitoyl-2-oleoyl phosphatidylglycerol: beta-sheet (residues 1-6), alpha-helix (residues 8-22) and random (residues 23-25) conformations. Analogous structural motifs are observed in the corresponding homologous N-terminal regions of several proteins that also share the 'saposin-like' (i.e. 5-helix bundle) folding pattern of full-length, human SP-B. In future studies, 13C-enhanced FTIR spectroscopy and energy minimizations may be of general use in defining backbone conformations at amino acid resolution, particularly for peptides or proteins in membrane environments.     

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.     

163 Enhanced sampling of peptides and proteins with a new biasing replica exchange method

Classical MD simulations (cMD) are limited by the sampling of relevant states of the peptides. Replica exchange (REMD) methods aim to search the conformational space of proteins more efficiently (reviewed in Ostermeir & Zacharias, 2013). We have developed a Hamiltonian REMD method that takes advantage of an intrinsic property of proteins, the specific Φ ? dihedral angle combinations along the polymer backbone. By employing a coupled two-dimensional biasing potential the energy barriers along the polymer backbone are reduced more effectively than by a previous approach based on a one-D biasing potential (Kannan & Zacharias, 2007). Thus, adjacent amino acids along the polymers backbone can easily switch between favourable regions in the Ramachandran plot. Additionally, energy barriers of rotameric states of amino acid side chains of proteins are also biased in the replica runs. The method improves the sampling of conformational substates of proteins at a modest number of replicas (nine replicas in the standard set-up with one replica running without biasing potential) compared to much larger numbers necessary in the case of standard temperature (T)-REMD simulations. A further improvement is achieved by a dynamical adjustment of the penalty potential levels in the replicas such that high exchange rates and improved mixing of conformations between different replicas are guaranteed. The biasing potential (BP)-REMD method turns out to be suitable to speed up both the folding of spaghetti-like test peptides and the refinement of loop decoy structures. Starting from extended structures, an α-helical oligo-alanine and β-hairpin chignolin and the Trp-cage protein fold more rapidly in near-native structures than in cMD simulations. The BP-REMD simulations not only accelerate the folding process of test proteins but also enlarge the variety of sampled configurations in conformational space. Since flexible parts of the protein can be penalized selectively, this method provides a precise tool to investigate regions of interest of the protein.     

A new approach to the rapid determination of protein side chain conformations

Two efficient algorithms have been developed which allow amino acid side chain conformations to be optimized rapidly for a given peptide backbone conformation. Both these approaches are based on the assumption that each side chain can be represented by a small number of rotameric states. These states have been obtained by a dynamic cluster analysis of a large data base of known crystallographic structures. Successful applications of these algorithms to the prediction of known protein conformations are presented.     

Conformational Entropy of the RNA Phosphate Backbone and Its Contribution to the Folding Free Energy

While major contributors to the free energy of RNA tertiary structures such as basepairing, base-stacking, and charge and counterion interactions have been studied extensively, little is known about the intrinsic free energy of the backbone. To assess the magnitude of the entropic strains along the phosphate backbone and their impact on the folding free energy, we have formulated a mathematical treatment for computing the volume of the main-chain torsion-angle conformation space between every pair of nucleobases along any sequence to compute the corresponding backbone entropy. To validate this method, we have compared the computed conformational entropies against a statistical free energy analysis of structures in the crystallographic database from several-thousand backbone conformations between nearest-neighbor nucleobases as well as against extensive computer simulations. Using this calculation, we analyzed the backbone entropy of several ribozymes and riboswitches and found that their entropic strains are highly localized along their sequences. The total entropic penalty due to steric congestions in the backbone for the native fold can be as high as 2.4 cal/K/mol per nucleotide for these medium and large RNAs, producing a contribution to the overall free energy of up to 0.72 kcal/mol per nucleotide. For these RNAs, we found that low-entropy high-strain residues are predominantly located at loops with tight turns and at tertiary interaction platforms with unusual structural motifs.     

Folding and stability of helical proteins: Carp myogen

In this work we use our very simple general representation of protein structures to study the mainly helical protein carp myogen. The representation, which treats the amino acid side-chains as simple spheres, is further simplified by rigidly fixing residues in α-helices. With this model we are able to reproduce the geometry and energetic stability of the native myogen conformation. Studies of the formation of α-helical sub-assemblies showed that the simulated folding of two and four-helix systems worked well, reaching compact native-like conformations with a good rate of success. Greater problems were encountered with the whole molecule (six helices), possibly due to the omission of entropic effects or to simulating the folding too rapidly. Finally, studies of the conformation of a pair of helices when isolated and when part of the whole molecule native conformation showed that long-range interactions have an unexpectedly strong influence on the conformation of the pair of helices.     

Molecular thermodynamic model to predict the alpha-helical secondary structure of polypeptide chains in solution.

The native state of a protein molecule in aqueous solutions represents one of the lowest states of Gibbs energy [Anfinsen, C.B. (1973) Science 181, 223-230]. Much progress has been made about the rules of protein folding [King, J. (1989) Chem. Eng. News 67, 32-54] and the dominant forces in protein folding [Dill, K.A. (1990) Biochemistry 29, 7133-7155]. However, the quantitative contributions of different Gibbs energy terms to protein stability remains a controversial issue [Moult, J., & Unger, R. (1991) Biochemistry 30, 3816-3824]. A molecular thermodynamic model has been proposed for the Gibbs energy of folding a residue in aqueous homopolypeptides from a random-coiled state to either the alpha-helix state or the beta-sheet state [Chen, C.-C., Zhu, Y., King, J.A., & Evans, L.B. (1992) Biopolymers 32, 1375-1392]. In this work, we present a generalization of the molecular thermodynamic model for the Gibbs energy of folding natural and synthetic heteropolypeptides from random-coiled conformations into alpha-helical conformations. The generalized model incorporates the intrinsic folding potential due to residue-solvent interactions, the cooperative folding effect due to residue-residue interactions, and the location and length of alpha-helices. The utility of the model was demonstrated by examining the stability of alpha-helical conformations of a number of natural polypeptides including C-peptide (residues 1-13) and S-peptide (residues 1-20) of RNase A (bovine pancreatic ribonuclease A), the P alpha fragment in BPTI (bovine pancreatic trypsin inhibitor), and synthetic polypeptides (the copolymers of different amino acid residues) including alanine-based peptides (16 or 17 residues long) in water. The computed Gibbs energies correspond well with the experimental data on helicity. The results also accounted for the effects of amino acid substitution and temperature on the stability of alpha-helical conformations of the test polypeptides.     

Two peptide fragments G55-I72 and K97-A109 from staphylococcal nuclease exhibit different behaviors in conformational preferences for helix formation

Two synthetic peptides, SNasealpha1 and SNasealpha2, corresponding to residues G55-I72 and K97-A109, respectively, of staphylococcal nuclease (SNase), are adopted for detecting the role of helix alpha1 (E57-A69) and helix alpha2 (M98-Q106) in the initiation of folding of SNase. The helix-forming tendencies of the two SNase peptide fragments are investigated using circular dichroism (CD) and two-dimensional (2D) nuclear magnetic resonance (NMR) methods in water and 40% trifluoroethanol (TFE) solutions. The coil-helix conformational transitions of the two peptides in the TFE-H2O mixture are different from each other. SNasealpha1 adopts a low population of localized helical conformation in water, and shows a gradual transition to helical conformation with increasing concentrations of TFE. SNasealpha2 is essentially unstructured in water, but undergoes a cooperative transition to a predominantly helical conformation at high TFE concentrations. Using the NMR data obtained in the presence of 40% TFE, an ensemble of alpha-helical structures has been calculated for both peptides in the absence of tertiary interactions. Analysis of all the experimental data available indicates that formation of ordered alpha-helical structures in the segments E57-A69 and M98-Q106 of SNase may require nonlocal interactions through transient contact with hydrophobic residues in other parts of the protein to stabilize the helical conformations in the folding. The folding of helix alpha1 is supposed to be effective in initiating protein folding. The formation of helix alpha2 depends strongly on the hydrophobic environment created in the protein folding, and is more important in the stabilization of the tertiary conformation of SNase.     

Preferred conformations of a linear RGD tripeptide.

Conformational study of RGD tripeptides in the nonhydrated and hydrated states was carried out using an empirical potential function ECEPP/3 and the hydration shell model in order to investigate preferred conformations and factors responsible for their stability. RGD tripeptides in the nonhydrated and hydrated states can be interpreted as existing as an ensemble of feasible conformations rather than as a single dominant conformation from the analysis of distributions of backbone conformations, hydrogen bonds and beta-turns. The different distributions of conformations for the neutral and zwitterionic RGD tripeptides in both states may indicate that the conformation of the RGD tripeptide is liable to depend on solvent polarity and pH values. beta-Turn populations for the neutral tripeptide in both states are reasonably consistent with NMR measurements on linear RGD-containing peptides. The degradation of RGD tripeptide seems to be attributed mainly to the hydrogen bonds between the Asp side-chain and the backbone of Asp residue or C-terminal NHMe group, rather than to the flexible backbones of Gly and Asp residues.     

Extracting information on folding from the amino acid sequence: accurate predictions for protein regions with preferred conformation in the absence of tertiary interactions.

A recently developed procedure to predict backbone structure from the amino acid sequence [Rooman, M., Kocher, J. P., & Wodak, S. (1991) J. Mol. Biol, 221, 961-979] is fine tuned to identify protein segments, of length 5-15 residues, that adopt well-defined conformations in the absence of tertiary interactions. These segments are obtained by requiring that their predicted lowest energy structures have a sizable energy gap relative to other computed conformations. Applying this procedure to 69 proteins of known structure, we find that regions with largest energy gaps--those having highly preferred conformations--are also the most accurately predicted ones. On the basis of previous findings that such regions correlate well with sites that become structured early during folding, our approach provides the means of identifying such sites in proteins without prior knowledge of the tertiary structure. Furthermore, when predictions are performed so as to ignore the influence of residues flanking each segment along the sequence, a situation akin to excising the considered peptide from the rest of the chain, they offer the possibility of identifying protein segments liable to adopt well-defined conformations on their own. The described approach should have useful applications in experimental and theoretical investigations of protein folding and stability, and aid in designing peptide drugs and vaccines.     

Conformational Entropy of the RNA Phosphate Backbone and Its Contribution to the Folding Free Energy

While major contributors to the free energy of RNA tertiary structures such as basepairing, base-stacking, and charge and counterion interactions have been studied extensively, little is known about the intrinsic free energy of the backbone. To assess the magnitude of the entropic strains along the phosphate backbone and their impact on the folding free energy, we have formulated a mathematical treatment for computing the volume of the main-chain torsion-angle conformation space between every pair of nucleobases along any sequence to compute the corresponding backbone entropy. To validate this method, we have compared the computed conformational entropies against a statistical free energy analysis of structures in the crystallographic database from several-thousand backbone conformations between nearest-neighbor nucleobases as well as against extensive computer simulations. Using this calculation, we analyzed the backbone entropy of several ribozymes and riboswitches and found that their entropic strains are highly localized along their sequences. The total entropic penalty due to steric congestions in the backbone for the native fold can be as high as 2.4 cal/K/mol per nucleotide for these medium and large RNAs, producing a contribution to the overall free energy of up to 0.72 kcal/mol per nucleotide. For these RNAs, we found that low-entropy high-strain residues are predominantly located at loops with tight turns and at tertiary interaction platforms with unusual structural motifs.     

PII structure in the model peptides for unfolded proteins: studies on ubiquitin fragments and several alanine-rich peptides containing QQQ, SSS, FFF, and VVV

A great deal of attention has been paid lately to the structures in unfolded proteins due to the recent discovery of many biologically functional but natively unfolded proteins and the far-reaching implications of order in unfolded states for protein folding. Recently, studies on oligo-Ala, oligo-Lys, oligo-Asp, and oligo-Glu, as well as oligo-Pro, have indicated that the left-handed polyproline II (PII) is the major local structure in these short peptides. Here, we show by NMR and CD studies that ubiquitin fragments, model unfolded peptides composed of nonrepeating amino acids, and four alanine-rich peptides containing QQQ, SSS, FFF, and VVV sequences are all present in aqueous solution predominantly in the extended PII or beta conformation. The results from this and related studies indicate that PII might be a major backbone conformation in unfolded proteins. The presence of defined local backbone structure in unfolded proteins is inconsistent with predictions from random coil models.     

Coupling between conformation and proton binding in proteins

Interest centers here on whether the use of a fixed charge distribution of a protein solute, or a treatment that considers proton-binding equilibria by solving the Poisson equation, is a better approach to discriminate native from non-native conformations of proteins. In this analysis of the charge distribution of 7 proteins, we estimate the solvation free energy contribution to the total free energy by exploring the 2(zeta) possible ionization states of the whole molecule, with zeta being the number of ionizable groups in the amino acid sequence, for every conformation in the ensembles of 7 proteins. As an additional consideration of the role of electrostatic interactions in determining the charge distribution of native folds, we carried out a comparison of alternative charge assignment models for the ionizable residues in a set of 21 native-like proteins. The results of this work indicate that (1) for 6 out of 7 proteins, estimation of solvent polarization based on the Generalized Born model with a fixed charge distribution provides the optimal trade-off between accuracy, with respect to the Poisson equation, and speed when compared to the accessible surface area model; for the seventh protein, consideration of all possible ionization states of the whole molecule appears to be crucial to discriminate the native from non-native conformations; (2) significant differences in the degree of ionization and hence the charge distribution for native folds are found between the different charge models examined; (3) the stability of the native state is determined by a delicate balance of all the energy components, and (4) conformational entropy, and hence the dynamics of folding, may play a crucial role for a successful ab initio protein folding prediction.     

Tryptophan rich peptides: influence of indole rings on backbone conformation

Synthetic peptides with defined secondary structure scaffolds, namely hairpins and helices, containing tryptophan residues, have been investigated in this study to probe the influence of a large number of aromatic amino acids on backbone conformations. Solution NMR investigations of Boc-W-L-W-(D)P-G-W-L-W-OMe (peptide 1), designed to form a well-folded hairpin, clearly indicates the influence of flanking aromatic residues at the (D)Pro-Gly region on both turn nucleation and strand propagation. Indole-pyrrolidine interactions in this peptide lead to the formation of the less-frequent type I' turn at the (D)Pro-Gly segment and frayed strand regions, with the strand residues adopting local helical conformations. An analog of peptide 1 with an Aib-Gly turn-nucleated hairpin (Boc-W-L-W-U-G-W-L-W-OMe (peptide 2)) shows a preference for helical structures in solution, in both chloroform and methanol. Peptides with either one (Boc-W-L-W-U-W-L-W-OMe (peptide 3)) or two (Boc-U-W-L-W-U-W-L-W-OMe (peptide 4)) helix-nucleating Aib residues give rise to the well-folded helical conformations in the chloroform solution. The results are indicative of a preference for helical folding in peptides containing a large number of Trp residues. Investigation of a tetrapeptide analog of peptide 2, Boc-W-U-G-W-OMe (peptide 5), in solution and in the crystal state (by X-ray diffraction), also indicates a preference for a helical fold. Additionally, peptide 5 is stabilized in crystals by both aromatic interactions and an array of weak interactions. Examination of Trp-rich sequences in protein structures, however, reveals no secondary structure preference, suggesting that other stabilizing interactions in a well-folded protein may offset the influence of indole rings on backbone conformations.     

A small tripeptide AFA undergoes two state cooperative conformational transitions: implications for conformational biases in unfolded states

It is important to understand the conformational biases that are present in unfolded states to understand protein folding. In this context, it is surprising that even a short tripeptide like AFA samples folded/ordered conformation as demonstrated recently by NMR experiments of the peptide in aqueous solution at 280 K. In this paper, we present molecular dynamics simulation of the peptide in explicit water using OPLS-AA/L all-atom force field. The results are in overall agreement with NMR results and provide some further insights. The peptide samples turn and extended conformational forms corresponding to minima in free energy landscape. Frequent transitions between the minima are observed due to modest free energy barriers. The turn conformation seems to be stabilized by hydrophobic interactions and possibly by bridging water molecules between backbone donors and acceptors. Thus the peptide does not sample conformations randomly, but samples well defined conformations. The peptide served as a model for folding-unfolding equilibrium in the context of peptide folding. Further, implications for drug design are also discussed.