Helix propensities of the amino acids measured in alanine-based peptides without helix-stabilizing side-chain interactions.

Helix propensities of the amino acids have been measured in alanine-based peptides in the absence of helix-stabilizing side-chain interactions. Fifty-eight peptides have been studied. A modified form of the Lifson-Roig theory for the helix-coil transition, which includes helix capping (Doig AJ, Chakrabartty A, Klingler TM, Baldwin RL, 1994, Biochemistry 33:3396-3403), was used to analyze the results. Substitutions were made at various positions of homologous helical peptides. Helix-capping interactions were found to contribute to helix stability, even when the substitution site was not at the end of the peptide. Analysis of our data with the original Lifson-Roig theory, which neglects capping effects, does not produce as good a fit to the experimental data as does analysis with the modified Lifson-Roig theory. At 0 degrees C, Ala is a strong helix former, Leu and Arg are helix-indifferent, and all other amino acids are helix breakers of varying severity. Because Ala has a small side chain that cannot interact significantly with other side chains, helix formation by Ala is stabilized predominantly by the backbone ("peptide H-bonds"). The implication for protein folding is that formation of peptide H-bonds can largely offset the unfavorable entropy change caused by fixing the peptide backbone. The helix propensities of most amino acids oppose folding; consequently, the majority of isolated helices derived from proteins are unstable, unless specific side-chain interactions stabilize them.     

Helix propagation and N-cap propensities of the amino acids measured in alanine-based peptides in 40 volume percent trifluoroethanol.

The helix propagation and N-cap propensities of the amino acids have been measured in alanine-based peptides in 40 volume percent trifluoroethanol (40% TFE) to determine if this helix-stabilizing solvent uniformly affects all amino acids. The propensities in 40% TFE are compared with revised values of the helix parameters of alanine-based peptides in water. Revision of the propensities in water is the result of redefining the capping statistical weights and evaluating the helix nucleation constant with N-capping explicitly included in the helix-coil model. The propagation propensities of all amino acids increase in 40% TFE relative to water, but the increases are highly variable. In water, all beta-branched and beta-substituted amino acids are helix breakers. In 40% TFE, the propagation propensities of the nonpolar amino acids increase greatly, leaving charged and neutral polar, beta-substituted amino acids as helix breakers. Glycine and proline are strong helix breakers in both solvents. Free energy differences for helix propagation (delta delta G) between alanine and other nonpolar amino acids are twice as large in water as predicted from side-chain conformational entropies, but delta delta G values in 40% TFE are close to those predicted from side-chain entropies. This dependence of delta delta G on the solvent points to a specific role of water in determining the relative helix propensities of the nonpolar amino acids. The N-cap propensities converge toward a common value in 40% TFE, suggesting that differential solvation by water contributes to the diversity of N-cap values shown by the amino acids.     

N- and C-capping preferences for all 20 amino acids in alpha-helical peptides.

We have determined the N- and C-capping preferences of all 20 amino acids by substituting residue X in the peptides NH2-XAKAAAAKAAAAKAAGY-CONH2 and in Ac-YGAAKAAAAKAAAAKAX-CO2H. Helix contents were measured by CD spectroscopy to obtain rank orders of capping preferences. The data were further analyzed by our modified Lifson-Roig helix-coil theory, which includes capping parameters (n and c), to find free energies of capping (-RT ln n and -RT ln c), relative to Ala. Results were obtained for charged and uncharged termini and for different charged states of titratable side chains. N-cap preferences varied from Asn (best) to Gln (worst). We find, as expected, that amino acids that can accept hydrogen bonds from otherwise free backbone NH groups, such as Asn, Asp, Ser, Thr, and Cys generally have the highest N-cap preference. Gly and acetyl group are favored, as are negative charges in side chains and at the N-terminus. Our N-cap preference scale agrees well with preferences in proteins. In contrast, we find little variation when changing the identity of the C-cap residue. We find no preference for Gly at the C-cap in contrast to the situation in proteins. Both N-cap and C-cap results for Tyr and Trp are inaccurate because their aromatic groups affect the CD spectrum. The data presented here are of value in rationalizing mutations at capping sites in proteins and in predicting the helix contents of peptides.     

α-Helix-forming propensities in peptides and proteins

Much effort has been invested in seeking to understand the thermodynamic basis of helix stability in both peptides and proteins. Recently, several groups have measured the helix-forming propensities of individual residues (Lyu, P. C., Liff, M. I., Marky, L. A., Kallenbach, N. R. Science 250:669–673, 1990; O'Neil, K. T., DeGrado, W. F. Science 250:646–651, 1990; Padmanabhan, S., Marqusee, S., Ridgeway, T., Laue, T. M., Baldwin, R. L. Nature (London) 344:268–270, 1990). Using Monte Carlo computer simulations, we tested the hypothesis that these differences in measured helix-forming propensity are due primarily to loss of side chain conformational entropy upon helix formation (Creamer, T. P., Rose, G. D. Proc. Natl. Acad. Sci. U.S.A. 89:5937–5941, 1992). Our previous study employed a rigid helix backbone, which is here generalized to a completely flexible helix model in order to ensure that earlier results were not a methodological artifact. Using this flexible model, side chain rotamer distributions and entropy losses are calculated and shown to agree with those obtained earlier. We note that the side chain conformational entropy calculated for Trp in our previous study was in error; a corrected value is presented. Extending earlier work, calculated entropy losses are found to correlate strongly with recent helix propensity scales derived from substitutions made within protein helices (Horovitz, A., Matthews, J. M., Fersht, A. R. J. Mol. Biol. 227:560–568, 1992; Blaber, M., Zhang, X.-J., Matthews, B. M. Science 260:1637–1640, 1993). In contrast, little correlation is found between these helix propensity scales and the accessible surface area buried upon formation of a model polyalanyl α-helix. Taken in sum, our results indicate that loss of side chain entropy is a major determinant of the helix-forming tendency of residues in both peptide and protein helices. © 1994 Wiley-Liss, Inc.     

Intrinsic secondary structure propensities of the amino acids,using statistical ϕ–ψ matrices: Comparison with experimental scales

Today there are several different experimental scales for the intrinsic α-helix as well as β-strand, propensities of the 20 amino acids obtained from the thermodynamic analysis of various model systems. These scales do not compare well with those extracted from statistical analysis of three-dimensional structure databases. Possible explanations for this could be the limited size of the databases used, the definitions of intrinsic propensities, or the theoretical approach. Here we report a statistical determination of α-helix and β-strand propensities derived from the analysis of a database of 279 three-dimensional structures. Contrary to what has been generally done, we have considered a particular residue as in α-helix or β-strand conformation by looking only at its dihedral angles (?–ψ matrices). Neither the identity nor the conformation of the surrounding residues in the amino acid sequence has been taken into consideration. Pseudoenergy empirical scales have been calculated from the statistical propensities. These scales agree very well with the experimental ones in relative and absolute terms. Moreover, its correlation with the average of the experimental scales for α-helix or β-strand is as good as the correlations of the individual experimental scales with the average. These results show that by using a large enough database and a proper definition for the secondary structure propensities, it is possible to obtain a scale as good as any of experimental origin. Interestingly the ?–ψ analysis of the Ramachandran plot suggests that the amino acids could have different β-strand propensities in different subregions of the β-strand area. © 1994 Wiley-Liss, Inc.     

Effect of the N1 residue on the stability of the alpha-helix for all 20 amino acids

N1 is the first residue in an alpha-helix. We have measured the contribution of all 20 amino acids to the stability of a small helical peptide CH(3)CO-XAAAAQAAAAQAAGY-NH(2) at the N1 position. By substituting every residue into the N1 position, we were able to investigate the stabilizing role of each amino acid in an isolated context. The helix content of each of the 20 peptides was measured by circular dichroism (CD) spectroscopy. The data were analyzed by our modified Lifson-Roig helix-coil theory, which includes the n1 parameter, to find free energies for placing a residue into the N1 position. The rank order for free energies is Asp(-), Ala > Glu(-) > Glu(0) > Trp, Leu, Ser > Asp(0), Thr, Gln, Met, Ile > Val, Pro > Lys(+), Arg, His(0) > Cys, Gly > Phe > Asn, Tyr, His(+). N1 preferences are clearly distinct from preferences for the preceding N-cap and alpha-helix interior. pK(a) values were measured for Asp, Glu, and His, and protonation-free energies were calculated for Asp and Glu. The dissociation of the Asp proton is less favorable than that of Glu, and this reflects its involvement in a stronger stabilizing interaction at the N terminus. Proline is not energetically favored at the alpha-helix N terminus despite having a high propensity for this position in crystal structures. The data presented are of value both in rationalizing mutations at N1 alpha-helix sites in proteins and in predicting the helix contents of peptides.     

Effect of the N2 residue on the stability of the α-helix for all 20 amino acids

N2 is the second position in the alpha-helix. All 20 amino acids were placed in the N2 position of a synthetic helical peptide (CH(3)CO-[AXAAAAKAAAAKAAGY]-NH(2)) and the helix content was measured by circular dichroism spectroscopy at 273K. The dependence of peptide helicity on N2 residue identity has been used to determine a free-energy scale by analysis with a modified Lifson-Roig helix-coil theory that includes a parameter for the N2 energy (n2). The rank order of DeltaDeltaG((relative to Ala)) is Glu(-), Asp(-) > Ala > Glu(0), Leu, Val, Gln, Thr, Ile, Ser, Met, Asp(0), His(0), Arg, Cys, Lys, Phe > Asn, > Gly, His(+), Pro, Tyr. The results correlate very well with N2 propensities in proteins, moderately well with N1 and helix interior preferences, and not at all with N-cap preferences. The strongest energetic effects result from interactions with the helix dipole, which favors negative charges at the helix N terminus. Hydrogen bonds to side chains at N2, such as Gln, Ser, and Thr, are weak, despite occurring frequently in protein crystal structures, in contrast to the N-cap position. This is because N-cap hydrogen bonds are close to linear, whereas N2 hydrogen bonds have poor geometry. These results can be used to modify protein stability rationally, help design helices, and improve prediction of helix location and stability.     

Microfolding: conformational probability map for the alanine dipeptide in water from molecular dynamics simulations

A direct attack on the protein-folding problem has been initiated with the free energy perturbation methods of molecular dynamics. The complete conformational probability map for the alanine dipeptide is presented. This work uses the SPC model for the explicit hydration of the dipeptide. Free energy differences for the four observed minima (beta, alpha R, alpha L, C7ax) are given, and the free energy barriers between minima are outlined.     

Position-resolved free energy of solvation for amino acids in lipid membranes from molecular dynamics simulations

Studies of insertion and interactions of amino acids in lipid membranes are pivotal to our understanding of membrane protein structure and function. Calculating the insertion cost as a function of transmembrane helix sequence is thus an important step towards improved membrane protein prediction and eventually drug design. Here, we present position-dependent free energies of solvation for all amino acid analogs along the membrane normal. The profiles cover the entire region from bulk water to hydrophobic core, and were produced from all-atom molecular dynamics simulations. Experimental differences corresponding to mutations and costs for entire segments match experimental data well, and in addition the profiles provide the spatial resolution currently not available from experiments. Polar side-chains largely maintain their hydration and assume quite ordered conformations, which indicates the solvation cost is mainly entropic. The cost of solvating charged side-chains is not only significantly lower than for implicit solvation models, but also close to experiments, meaning these could well maintain their protonation states inside the membrane. The single notable exception to the experimental agreement is proline, which is quite expensive to introduce in vivo despite its hydrophobicity--a difference possibly explained by kinks making it harder to insert helices in the translocon.     

Similarities in the thermodynamics and kinetics of aggregation of disease-related Abeta(1-40) peptides

Increasing evidence indicates that polypeptide aggregation often involves a nucleation and a growth phase, although the relationship between the factors that determine these two phases has not yet been fully clarified. We present here an analysis of several mutations at different sites of the Abeta(1-40) peptide, including those associated with early onset forms of the Alzheimer's disease, which reveals that the effects of specific amino acid substitutions in the sequence of this peptide are strongly modulated by their structural context. Nevertheless, mutations at different positions perturb in a correlated manner the free energies of aggregation as well as the lag times and growth rates. We show that these observations can be rationalized in terms of the intrinsic propensities for aggregation of the Abeta(1-40) sequence, thus suggesting that, in the case of this peptide, the determinants of the thermodynamics and of the nucleation and growth of the aggregates have a similar physicochemical basis.     

Removal of kinetic traps and enhanced protein folding by strategic substitution of amino acids in a model alpha-helical hairpin peptide

The presence of non-native kinetic traps in the free energy landscape of a protein may significantly lengthen the overall folding time so that the folding process becomes unreliable. We use a computational model alpha-helical hairpin peptide to calculate structural free energy landscapes and relate them to the kinetics of folding. We show how protein engineering through strategic changes in only a few amino acid residues along the primary sequence can greatly increase the speed and reliability of the folding process, as seen experimentally. These strategic substitutions also prevent the formation of long-lived misfolded configurations that can cause unwanted aggregations of peptides. These results support arguments that removal of kinetic traps, obligatory or nonobligatory, is crucial for fast folding.     

Comparison of binding energies of SrcSH2-phosphotyrosyl peptides with structure-based prediction using surface area based empirical parameterization

The prediction of binding energies from the three-dimensional (3D) structure of a protein-ligand complex is an important goal of biophysics and structural biology. Here, we critically assess the use of empirical, solvent-accessible surface area-based calculations for the prediction of the binding of Src-SH2 domain with a series of tyrosyl phosphopeptides based on the high-affinity ligand from the hamster middle T antigen (hmT), where the residue in the pY+ 3 position has been changed. Two other peptides based on the C-terminal regulatory site of the Src protein and the platelet-derived growth factor receptor (PDGFR) are also investigated. Here, we take into account the effects of proton linkage on binding, and test five different surface area-based models that include different treatments for the contributions to conformational change and protein solvation. These differences relate to the treatment of conformational flexibility in the peptide ligand and the inclusion of proximal ordered solvent molecules in the surface area calculations. This allowed the calculation of a range of thermodynamic state functions (deltaCp, deltaS, deltaH, and deltaG) directly from structure. Comparison with the experimentally derived data shows little agreement for the interaction of SrcSH2 domain and the range of tyrosyl phosphopeptides. Furthermore, the adoption of the different models to treat conformational change and solvation has a dramatic effect on the calculated thermodynamic functions, making the predicted binding energies highly model dependent. While empirical, solvent-accessible surface area based calculations are becoming widely adopted to interpret thermodynamic data, this study highlights potential problems with application and interpretation of this type of approach. There is undoubtedly some agreement between predicted and experimentally determined thermodynamic parameters: however, the tolerance of this approach is not sufficient to make it ubiquitously applicable.     

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.     

Helix-coil transition of alanine peptides in water: force field dependence on the folded and unfolded structures

The force fields used in classical modeling studies are semiempirical in nature and rely on their validation by comparison of simulations with experimental data. The all-atom replica-exchange molecular dynamics (REMD) methodology allows us to calculate the thermodynamics of folding/unfolding of peptides and small proteins, and provides a way of evaluating the reliability of force fields. We apply the REMD to obtain equilibrium folding/unfolding thermodynamics of a 21-residue peptide containing only alanine residues in explicit aqueous solution. The thermodynamics of this peptide is modeled with both the OPLS/AA/L and the A94/MOD force fields. We find that the helical content and the values for the helix propagation and nucleation parameters for this alanine peptide are consistent with measurements on similar peptides and with calculations using the modified AMBER force field (A94/MOD). The nature of conformations, both folded and unfolded, that contributes to the helix-coil transition profile, however, is quite different between these two force fields.     

Fuzzy cluster analysis of simple physicochemical properties of amino acids for recognizing secondary structure in proteins.

Fuzzy cluster analysis has been applied to the 20 amino acids by using 65 physicochemical properties as a basis for classification. The clustering products, the fuzzy sets (i.e., classical sets with associated membership functions), have provided a new measure of amino acid similarities for use in protein folding studies. This work demonstrates that fuzzy sets of simple molecular attributes, when assigned to amino acid residues in a protein''s sequence, can predict the secondary structure of the sequence with reasonable accuracy. An approach is presented for discriminating standard folding states, using near-optimum information splitting in half-overlapping segments of the sequence of assigned membership functions. The method is applied to a nonredundant set of 252 proteins and yields approximately 73% matching for correctly predicted and correctly rejected residues with approximately 60% overall success rate for the correctly recognized ones in three folding states: alpha-helix, beta-strand, and coil. The most useful attributes for discriminating these states appear to be related to size, polarity, and thermodynamic factors. Van der Waals volume, apparent average thickness of surrounding molecular free volume, and a measure of dimensionless surface electron density can explain approximately 95% of prediction results. hydrogen bonding and hydrophobicity induces do not yet enable clear clustering and prediction.     

The barrier for proton transport in aquaporins as a challenge for electrostatic models: the role of protein relaxation in mutational calculations

The origin of the barrier for proton transport through the aquaporin channel is a problem of general interest. It is becoming increasingly clear that this barrier is not attributable to the orientation of the water molecules across the channel but rather to the electrostatic penalty for moving the proton charge to the center of the channel. However, the reason for the high electrostatic barrier is still rather controversial. It has been argued by some workers that the barrier is due to the so-called NPA motif and/or to the helix macrodipole or to other specific elements. However, our works indicated that the main reason for the high barrier is the loss of the generalized solvation upon moving the proton charge from the bulk to the center of the channel and that this does not reflect a specific repulsive electrostatic interaction but the absence of sufficient electrostatic stabilization. At this stage it seems that the elucidation and clarification of the origin of the electrostatic barrier can serve as an instructive test case for electrostatic models. Thus, we reexamine the free-energy surface for proton transport in aquaporins using the microscopic free-energy perturbation/umbrella sampling (FEP/US) and the empirical valence bond/umbrella sampling (EVB/US) methods as well as the semimacroscopic protein dipole Langevin dipole model in its linear response approximation version (the PDLD/S-LRA). These extensive studies help to clarify the nature of the barrier and to establish the "reduced solvation effect" as the primary source of this barrier. That is, it is found that the barrier is associated with the loss of the generalized solvation energy (which includes of course all electrostatic effects) upon moving the proton charge from the bulk solvent to the center of the channel. It is also demonstrated that the residues in the NPA region and the helix dipole cannot be considered as the main reasons for the electrostatic barrier. Furthermore, our microscopic and semimacroscopic studies clarify the problems with incomplete alternative calculations, illustrating that the effects of various electrostatic elements are drastically overestimated by macroscopic calculations that use a low dielectric constant and do not consider the protein reorganization. Similarly, it is pointed out that microscopic potential of mean force calculations that do not evaluate the electrostatic barrier relative to the bulk water cannot be used to establish the origin of the electrostatic barrier. The relationship between the present study and calculations of pK(a)s in protein interiors is clarified, pointing out that approaches that are applied to study the aquaporin barrier should be validated by pK(a)s calculations. Such calculations also help to clarify the crucial role of solvation energies in establishing the barrier in aquaporins.