Alpha-helix stabilization by alanine relative to glycine: roles of polar and apolar solvent exposures and of backbone entropy

The energetics of alpha-helix formation are fairly well understood and the helix content of a given amino acid sequence can be calculated with reasonable accuracy from helix-coil transition theories that assign to the different residues specific effects on helix stability. In internal helical positions, alanine is regarded as the most stabilizing residue, whereas glycine, after proline, is the more destabilizing. The difference in stabilization afforded by alanine and glycine has been explained by invoking various physical reasons, including the hydrophobic effect and the entropy of folding. Herein, the contribution of these two effects and that of hydrophilic area burial is evaluated by analyzing Ala and Gly mutants implemented in three helices of apoflavodoxin. These data, combined with available data for similar mutations in other proteins (22 Ala/Gly mutations in alpha-helices have been considered), allow estimation of the difference in backbone entropy between alanine and glycine and evaluation of its contribution and that of apolar and polar area burial to the helical stabilization typically associated to Gly-->Ala substitutions. Alanine consistently stabilizes the helical conformation relative to glycine because it buries more apolar area upon folding and because its backbone entropy is lower. However, the relative contribution of polar area burial (which is shown to be destabilizing) and of backbone entropy critically depends on the approximation used to model the structure of the denatured state. In this respect, the excised-peptide model of the unfolded state, proposed by Creamer and coworkers (1995), predicts a major contribution of polar area burial, which is in good agreement with recent quantitations of the relative enthalpic contribution of Ala and Gly residues to alpha-helix formation.     

Conformational analysis of protein structures derived from NMR data

A study is presented of the conformational characteristics of NMR-derived protein structures in the Protein Data Bank compared to X-ray structures. Both ensemble and energy-minimized average structures are analyzed. We have addressed the problem using the methods developed for crystal structures by examining the distribution of ?, Ψ, and χ angles as indicators of global conformational irregularity. All these features in NMR structures occur to varying degrees in multiple conformational states. Some measures of local geometry are very tightly constrained by the methods used to generate the structure, e.g., proline ? angles, α-helix ?, Ψ angles, ω angles, and C_α chirality. The more lightly restrained torsion angles do show increasead clustering as the number of overall experimental observations increases. ?, Ψ, and χ₁ angle conformational heterogeneity is strongly correlated with accessibility but shows additional differences which reflect the differing number of observations possible in NMR for the various side chains (e.g., many for Trp, few for Ser). In general, we find that the core is defined to a notional resolution of 2.0 to 2.3 Å. Of real interest is the behavior of surface residues and in particular the side chains where multiple rotameric states in different structures can vary from 10% to 88%. Later generation structures show a much tighter definition which correlates with increasing use of J-coupling information, stereospecific assignments, and heteronumclear techniques. A suite of programs is being developed to address the special needs of NMR-derived structures which will take into account the existence of increased mobility in solution. © 1993 Wiley-Liss, Inc.     

Proline in alpha-helix: stability and conformation studied by dynamics simulation

Free-energy simulations have been used to estimate the change in the conformational stability of short polyalanine alpha-helices when one of the alanines is replaced by a proline residue. For substituting proline in the middle of the helix the change in free energy of folding (delta delta G degrees) was calculated as 14 kJ/mol (3.4 kcal/mol), in excellent agreement with the one available experimental value. The helix containing proline was found to be strongly kinked; the free energy for reducing the angle of the kink from 40 degrees to 15 degrees was calculated, and found to be small. A tendency to alternate hydrogen bonding schemes was observed in the proline-containing helix. These observations for the oligopeptide agree well with the observation of a range of kink angles (18-35 degrees) and variety of hydrogen bonding schemes, in the rare instances where proline occurs in helices in globular proteins. For substituting proline at the N-terminus of the helix the change in free energy of folding (delta delta G degrees) was calculated as -4 kJ/mol in the first helical position (N1) and +6 kJ/mol in the second helical position (N2). The observed frequent occurrence of proline in position N1 in alpha-helices in proteins therefore has its origin in stability differences of secondary structure. The conclusion reached here that proline may be a better helix former in position N1 than (even) alanine, and thus be a helix initiator may be testable experimentally by measurements of fraction helical conformation of individual residues in oligopeptides of appropriate sequence. The relevance of these results in regards to the frequent occurrence of proline-containing helices in certain membrane proteins is discussed.     

Increasing protein stability by improving beta‐turns

Our goal was to gain a better understanding of how protein stability can be increased by improving β‐turns. We studied 22 β‐turns in nine proteins with 66–370 residues by replacing other residues with proline and glycine and measuring the stability. These two residues are statistically preferred in some β‐turn positions. We studied: Cold shock protein B (CspB), Histidine‐containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase α‐subunit, and Maltose binding protein. Of the 15 single proline mutations, 11 increased stability (Average = 0.8 ± 0.3; Range = 0.3–1.5 kcal/mol), and the stabilizing effect of double proline mutants was additive. On the basis of this and our previous work, we conclude that proteins can generally be stabilized by replacing nonproline residues with proline residues at the i + 1 position of Type I and II β‐turns and at the i position in Type II β‐turns. Other turn positions can sometimes be used if the φ angle is near ?60° for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in β‐turns in related proteins give similar results. Proline substitutions increase stability mainly by decreasing the entropy of the denatured state. In contrast, the large, diverse group of proteins considered here had almost no residues in β‐turns that could be replaced by Gly to increase protein stability. Improving β‐turns by substituting Pro residues is a generally useful way of increasing protein stability. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

β-Turn conformations in crystal structures of model peptides containing α,α-Di-n-propylglycine and α,α-Di-n-butylglycine

The crystal state conformations of three peptides containing the α,α-dialkylated residues. α,α-di-n-propylglycine (Dpg) and α,α-di-n-butylglycine (Dbg), have been established by x-ray diffraction. Boc-Ala-Dpg-Alu-OMe (I) and Boc-Ala-Dbg-Ala-OMe (III) adopt distorted type II β-turn conformations with Ala (1) and Dpg/Dbg (2) as the corner residues. In both peptides the conformational angles at the Dxg residue (I: ? = 66.2°, ψ = 19.3°; III: ? = 66.5°. ψ = 21.1°) deviate appreciably from ideal values for the i + 2 residue in a type II β-turn. In both peptides the observed (N…O) distances between the Boc CO and Ala (3) NH groups are far too long (1: 3.44 Å: III: 3.63 Å) for an intramolecular 4 → 1 hydrogen bond. Boc-Ala-Dpg-Ata-NHMe (II) crystallizes with two independent molecules in the asymmetric unit. Both molecules HA and HB adopt consecutive β-turn (type III-III in HA and type III-I in IIB) or incipient 3₁₀-helical structures, stabilized by two intramolecular 4 → 1 hydrogen bonds. In all four molecules the bond angle N-C^α-C′ (τ) at the Dxg residues are ≥ 110°. The observation of conformational angles in the helical region of ?,ψ space at these residues is consistent with theoretical predictions. © 1995 John Wiley & Sons, Inc.     

Quantitative Determination of Site-Specific Conformational Distributions in an Unfolded Protein by Solid-State Nuclear Magnetic Resonance

Solid-state nuclear magnetic resonance (NMR) techniques are used to investigate the structure of the 35-residue villin headpiece subdomain (HP35) in folded, partially denatured, and fully denatured states. Experiments are carried out in frozen glycerol/water solutions, with chemical denaturation by guanidine hydrochloride (GdnHCl). Without GdnHCl, two-dimensional solid-state ¹³C NMR spectra of samples prepared with uniform ¹³C labeling of selected residues show relatively sharp cross-peaks at chemical shifts that are consistent with the known three-helix bundle structure of folded HP35. At high GdnHCl concentrations, most cross-peaks broaden and shift, qualitatively indicating disruption of the folded structure and development of static conformational disorder in the frozen denatured state. Conformational distributions at one residue in each helical segment are probed quantitatively with three solid-state NMR techniques that provide independent constraints on backbone ? and ψ torsion angles in samples with sequential pairs of carbonyl ¹³C labels. Without GdnHCl, the combined data are well fit by α-helical conformations. At [GdnHCl] = 4.5 M, corresponding to the approximate denaturation midpoint, the combined data are well fit by a combination of α-helical and partially extended conformations at each site, but with a site-dependent population ratio. At [GdnHCl] = 7.0 M, corresponding to the fully denatured state, the combined data are well fit by a combination of partially extended and polyproline II conformations, again with a site-dependent population ratio. Two entirely different models for conformational distributions lead to nearly the same best-fit distributions, demonstrating the robustness of these conclusions. This work represents the first quantitative investigation of site-specific conformational distributions in partially folded and unfolded states of a protein by solid-state NMR.     

Sequence composition effects on denatured state loop formation in iso-1-cytochrome c variants: polyalanine versus polyglycine inserts

Protein folding is dependent on the formation and persistence of simple loops during the earliest events of the folding process. Ease of loop formation and persistence is believed to be dependent on the steric properties of the residues involved in loop formation. We have investigated this conformational factor in the denatured state of iso-1-cytchrome c using a five alanine insert in front of a unique histidine in the N-terminal region of the protein. The alanine residues have then been progressively substituted with sterically less-constrained glycine residues. Guanidine-HCl unfolding shows that all variants have a free energy of unfolding of approximately 2 kcal/mol. The low stability of these variants is well accounted for by stabilization of the denatured state by histidine-heme loop formation. The stability of the 22 residue histidine-heme loop has been measured in 3 M guanidine hydrochloride for all variants. Surprisingly, relative to alanine, glycine has only a very modest effect on equilibrium loop stability. Thus, the greater flexibility that glycine confers on the main-chain provides no advantage in terms of the persistence of simple loops early in folding. The underlying basis for the similar behavior of loops with polyalanine versus polyglycine inserts is discussed in terms of the current knowledge of the structure and loop formation kinetics of glycine versus alanine-rich peptides.     

Conformational dynamics is more important than helical propensity for the folding of the all α‐helical protein Im7

Im7 folds via an on‐pathway intermediate that contains three of the four native α‐helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94‐residue variant of the 87‐residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model‐free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea‐denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α‐helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α‐helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea‐denatured states promote formation of a three‐helix intermediate in which the residues that form helix III remain non‐helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.     

Conformational analysis of a set of peptides corresponding to the entire primary sequence of the N-terminal domain of the ribosomal protein L9: evidence for stable native-like secondary structure in the unfolded state

There is considerable interest in the structure of the denatured state and in the role local interactions play in protein stability and protein folding. Studies of peptide fragments provide one method to assess local conformational preferences which may be present in the denatured state under native-like conditions. A set of peptides corresponding to the individual elements of secondary structure derived from the N-terminal domain of the ribosomal protein L9 have been synthesized. This small 56 residue protein adopts a mixed alpha-beta topology and has been shown to fold rapidly in an apparent two-state fashion. The conformational preferences of each peptide have been analyzed by proton nuclear magnetic resonance spectroscopy and circular dichroism spectroscopy. Peptides corresponding to each of the three beta-stands and to the first alpha-helix are unstructured as judged by CD and NMR. In contrast, a peptide corresponding to the C-terminal helix is remarkably structured. This 17 residue peptide is 53 % helical at pH 5.4, 4 degrees C. Two-dimensional NMR studies demonstrate that the helical structure is distributed approximately uniformly throughout the peptide, although there is some evidence for fraying at the C terminus. Detailed analysis of the NMR spectra indicate that the helix is stabilized, in part, by a native N-capping interaction involving Thr40. A mutant peptide which lacks Thr40 is only 32 % helical. pH and ionic strength-dependent studies suggested that charge charge interactions make only a modest net contribution to the stability of the peptide. The protein contains a trans proline peptide bond located at the first position of the C-terminal helix. NMR analysis of the helical peptide and of a smaller peptide containing the proline residue indicates that only a small amount of cis proline isomer (8 %) is likely to be populated in the unfolded state.     

Effect of mutations involving charged residues on the stability of staphylococcal nuclease: a continuum electrostatics study

A continuum electrostatics model is used to calculate the relative stabilities of 117 mutants of staphylococcal nuclease (SNase) involving the mutation of a charged residue to an uncharged residue. The calculations are based on the crystallographic structure of the wild-type protein and attempt to take implicitly into account the effect of the mutations in the denatured state by assuming a linear relationship between the free energy changes caused by the mutation in the native and denatured states. A good correlation (linear correlation coefficient of approximately 0.8) is found with published experimental relative stabilities of these mutants. The results suggest that in the case of SNase (i) charged residues contribute to the stability of the native state mainly through electrostatic interactions, and (ii) native-like electrostatic interactions may persist in the denatured state. The continuum electrostatics method is only moderately sensitive to model parameters and leads to quasi-predictive results for the relative mutant stabilities (error of 2-3 kJ mol(-1) or of the order of k(B)T), except for mutants in which a charged residue is mutated to glycine.     

Correlation between changes in nuclear magnetic resonance order parameters and conformational entropy: molecular dynamics simulations of native and denatured staphylococcal nuclease

Recent work has suggested that changes in NMR order parameters may quantitatively reflect changes in the conformational entropy of a protein ensemble. The extent of the mathematical relationship between local entropy changes as seen by NMR order parameters and the full protein entropy change is a complex issue. As a step towards a fuller understanding of this problem, molecular dynamics calculations of both native and denatured staphylococcal nuclease were performed. The N-H bond vector motion, in both explicit and implicit solvent, was analyzed to estimate local and global entropy changes. The calculated N-H bond vector order parameters from simulation agreed on average with experimental values for both native and denatured structures. However, the inverted-U profile of order parameters versus residue number observed experimentally for denatured nuclease was only partially reproduced by simulation of compact denatured structures. Comparisons made across the full set of simulations revealed a correlation between the N-H order parameter-based conformational entropy change and the total quasiharmonic-based conformational entropy change between the native and denatured structures. The calculations showed that about 25% of the total entropy change was reflected by changes in simulated S2 values. This result suggests that NMR-derived order parameters may be used to provide a reasonable estimate of the total conformational entropy change on protein folding.     

The natively helical chain segment 169-188 of Escherichia coli adenylate kinase is formed in the latest phase of the refolding transition

The refolding transition of Escherichia coli adenylate kinase (AK) was investigated by monitoring the refolding kinetics of a selected 20 residue helical segment in the CORE domain of the protein. Residues 169 and 188 were labeled by 1-acetamido-methyl-pyrene, and by bimane, respectively. The experiment combines double-jump stopped-flow fast mixing initiation of refolding and time-resolved F?rster energy transfer spectroscopy for monitoring the conformational transitions (double-kinetics experiment). Two kinetic phases were found in the denaturant-induced unfolding of AK. In the first phase, the fluorescence quantum yields of both probes decreased. The distribution of the distances between them transformed from the native state's narrow distribution with the mean distance corresponding to the distance in the crystal structure, to a distribution compatible with an unordered structure. In the second, slow step of denaturation, neither the fluorescence parameters of the probes nor the distance distribution between them changed. This step appeared to be a transformation of the fast-folding species formed in the first phase, to the slow-folding species. Refolding of the fast-folding species of the denatured state of AK was also a two-phase process. During the first fast phase, within less than 5ms, the fluorescence emission of both probes increased, but the distance distribution between the labeled sites was unchanged. Only during the second slow refolding step did the intramolecular distance distribution change from the characteristic of the denatured state to the narrow distribution of the native state. This experiment shows that for the case of the CORE domain of AK, the large helical segment of residues 169-188 was not formed in the first compaction step of refolding. The helical conformation of this segment is established only in the second, much slower, refolding phase, simultaneously with the completion of the native structure.     

Insight into polyproline II helical bundle stability in an antifreeze protein denatured state

The use of polyproline II (PPII) helices in protein design is currently hindered by limitations in our understanding of their conformational stability and folding. Recent studies of the snow flea antifreeze protein (sfAFP), a useful model system composed of six PPII helices, suggested that a low denatured state entropy contributes to folding thermodynamics. Here, circular dichroism spectroscopy revealed minor populations of PPII like conformers at low temperature. To get atomic level information on the conformational ensemble and entropy of the reduced, denatured state of sfAFP, we have analyzed its chemical shifts and {¹H}-¹⁵N relaxation parameters by NMR spectroscopy at four experimental conditions. No significant populations of stable secondary structure were detected. The stiffening of certain N-terminal residues at neutral versus acidic pH and shifted pK_a values leads us to suggest that favorable charge-charge interactions could bias the conformational ensemble to favor the formation the C1-C28 disulfide bond during nascent folding, although no evidence for preferred contacts between these positions was detected by paramagnetic relaxation enhancement under denaturing conditions. Despite a high content of flexible glycine residues, the mobility of the sfAFP denatured ensemble is similar for denatured α/β proteins both on fast ps/ns as well as slower μs/ms timescales. These results are in line with a conformational entropy in the denatured ensemble resembling that of typical proteins and suggest that new structures based on PPII helical bundles should be amenable to protein design.     

Study of the compact denatured state of a protein by molecular dynamics simulation]

The native state can be considered as a unique conformation of the protein molecule with the lowest free energy of residue contacts. In this case, all other conformations correspond to the denatured state. The degree of their compactness varies significantly. Under folding conditions, the compact denatured state rather than the random coil is in equilibrium with native protein. The balance between the main forces of protein folding, the solvophobic interactions and conformational entropy, suggests that some properties of the compact denatured state are close to those of native protein, whereas other properties are close to those of the random coil. To investigate the molecular structure of the compact denatured state, the method of molecular dynamics simulation seems to be very useful.     

Influence of proline residues in transmembrane helix packing

Integral membrane proteins often contain proline residues in their alpha-helical transmembrane (TM) fragments, which may strongly influence their folding and association. Pro-scanning mutagenesis of the helical domain of glycophorin A (GpA) showed that replacement of the residues located at the center abrogates helix packing while substitution of the residues forming the ending helical turns allows dimer formation. Synthetic TM peptides revealed that a point mutation of one of the residues of the dimerization motif (L75P) located at the N-terminal helical turn of the GpA TM fragment, adopts a secondary structure and oligomeric state similar to the wild-type sequence in detergents. In addition, both glycosylation mapping in biological membranes and molecular dynamics showed that the presence of a proline residue at the lipid/water interface has as an effect the extension of the helical end. Thus, helix packing can be an important factor that determines appearance of proline in TM helices. Membrane proteins might accumulate proline residues at the two ends of their TM segments in order to modulate the exposition of key amino acid residues at the interface for molecular recognition events while allowing stable association and native folding.     

Sequence conservation provides the best prediction of the role of proline residues in p13suc1

The unique nature of the proline side-chain imposes severe constraints on the polypeptide backbone, and thus it seems likely that it plays a special structural or functional role in the architecture of proteins. We have investigated the role of proline residues in suc1, a member of the cyclin-dependent kinase (cks) family of proteins, whose known function is to bind to and regulate the activity of the major mitotic cdk. The effect on stability of mutation to alanine of all but two of the eight proline residues is correlated with their conservation within the family. The remaining two proline residues are located in the hinge loop between two beta-strands that mediates a domain-swapping process involving exchange of a beta-strand between two monomers to form a dimer pair. Mutation of these proline residues to alanine stabilises the protein. cdk binding is unaffected by these mutations, but dimerisation is altered. We propose, therefore, that the double-proline motif is conserved for the purpose of domain swapping, which suggests that this phenomenon plays a role in the function of cks proteins. Thus, the conservation of the proline residues is a good indicator of their roles in suc1, either in the stabilisation of the native state or in performing functions that are as yet unknown. In addition, the strain resulting from two of the proline residues was relieved successfully by mutation of the preceeding residue to glycine, suggesting a general method for designing more stable proteins.     

Proline for alanine substitutions in the C-peptide helix of ribonuclease A

The effect on overall alpha-helix content of substituting proline for alanine has been determined at 5 positions (1, 2, 4, 5, and 13) of a 13-residue peptide related in sequence to residues 1-13 of ribonuclease A. The helix content falls off rapidly as proline is moved inward, and the proline residue effectively truncates the helix. No helix-stabilizing effect of proline is found at positions 2 or 4 within the first turn of the helix. Proline substitution at either end position (1, 13) has little effect on overall helix content, in agreement with an earlier study of glycine for alanine substitutions. The two end residues of the helix appear to be strongly frayed.     

Action of crystalline acid carboxypeptidase from Penicillium janthinellum.

Acid carboxypeptidase (EC 3.4.12.-) crystallized from culture filtrate of Penicillium janthinellum has been investigated for its use in carboxy-terminal sequence determination of Z-Gly-Pro-Leu-Gly, Z-Gly-Pro-Leu-Gly-Pro, angiotensin I, native lysozyme, native ribonuclease T1, and reduced S-carboxy-methyl-lysozyme. The examination indicated that proline and glycine were liberated from Z-Gly-Pro-Leu-Gly-Pro. At high enzyme concentration, the enzyme catalyzed complete sequential release of amino acids from the carboxy-terminal leucine to the amino-terminal aspartic acid of angiotensin I. The enzyme released the carboxy-terminal leucine from native lysozyme, however, no release of the threonine from native ribonuclease T1 was observed after a prolonged period of incubation with the enzyme. The sequence of the first nine carboxy-terminal residues of denatured lysozyme, leucine, arginine, S-carboxymethyl-cysteine, glycine, arginine, isoleucine, tryptophane, alanine, and glutamine, could be deduced unequivocally from a time release plot of an incubation mixture with the enzyme.     

Structural and dynamic characterization of an unfolded state of poplar apo-plastocyanin formed under nondenaturing conditions

Plastocyanin is a predominantly beta-sheet protein containing a type I copper center. The conformational ensemble of a denatured state of apo-plastocyanin formed in solution under conditions of low salt and neutral pH has been investigated by multidimensional heteronuclear NMR spectroscopy. Chemical shift assignments were obtained by using three-dimensional triple-resonance NMR experiments to trace through-bond heteronuclear connectivities along the backbone and side chains. The (3)J(HN,Halpha) coupling constants, (15)N-edited proton-proton nuclear Overhauser effects (NOEs), and (15)N relaxation parameters were also measured for the purpose of structural and dynamic characterization. Most of the residues corresponding to beta-strands in the folded protein exhibit small upfield shifts of the (13)C(alpha) and (13)CO resonances relative to random coil values, suggesting a slight preference for backbone dihedral angles in the beta region of (phi,psi) space. This is further supported by the presence of strong sequential d(alphaN)(i, i + 1) NOEs throughout the sequence. The few d(NN)(i, i + 1) proton NOEs that are observed are mostly in regions that form loops in the native plastocyanin structure. No medium or long-range NOEs were observed. A short sequence, between residues 59 and 63, was found to populate a nonnative helical conformation in the unfolded state, as indicated by the shift of the (13)C(alpha), (13)CO, and (1)H(alpha) resonances relative to random coil values and by the decreased values of the (3)J(HN,Halpha) coupling constants. The (15)N relaxation parameters indicate restriction of motions on a nanosecond timescale in this region. Intriguingly, this helical conformation is present in a sequence that is close to but not in the same location as the single short helix in the native folded protein. The results are consistent with earlier NMR studies of peptide fragments of plastocyanin and confirm that the regions of the sequence that form beta-strands in the native protein spontaneously populate the beta-region of (phi,psi) space under folding conditions, even in the absence of stabilizing tertiary interactions. We conclude that the state of apo-plastocyanin present under nondenaturing conditions is a noncompact unfolded state with some evidence of nativelike and nonnative local structuring that may be initiation sites for folding of the protein.