Motional dynamics of residues in a beta-hairpin peptide measured by 13C-NMR relaxation.

Structurally characterizing partially folded peptides is problematic given the nature of their transient conformational states. 13C-NMR relaxation data can provide information on the geometry of bond rotations, motional restrictions, and correlated bond rotations of the backbone and side chains and, therefore, is one approach that is useful to assess the presence of folded structure within a conformational ensemble. A peptide 12mer, R1GITVNG7KTYGR12, has been shown to partially fold in a relatively stable beta-hairpin conformation centered at NG. Here, five residues, G2, V5, G7, Y10, G11, were selectively 13C-enriched, and 13C-NMR relaxation experiments were performed to obtain auto- and cross-correlation motional order parameters, correlation times, bond rotation angular variances, and bond rotational correlation coefficients. Our results indicate that, of the three glycines, G7 within the hairpin beta-turn displays the most correlated phi(t),psi(t) rotations with its axis of rotation bisecting the angle defined by the H-C-H bonds. These positively correlated bond rotations give rise to "twisting" type motions of the HCH group. V5 and Y10 phi,psi bond rotations are also positively correlated, with their CbetaCalphaH groups undergoing similar "twisting" type motions. Motions of near-terminal residues G2 and G11 are less restricted and less correlated and are best described as wobbling-in-a-cone. V5 and Y10 side-chain motions, aside from being highly restricted, were found to be correlated with phi,psi bond rotations. At 303 K, where the hairpin is considered "unfolded," the peptide exists in a transient, collapsed state because backbone and side-chain motions of V5, G7, and Y10 remain relatively restricted, unlike their counterparts in GXG-based tripeptides. These results provide unique information toward understanding conformational variability in the unfolded state of proteins, which is necessary to solve the protein folding problem.     

Internal motional amplitudes and correlated bond rotations in an alpha-helical peptide derived from 13C and 15N NMR relaxation

Peptide GFSKAELAKARAAKRGGY folds in an alpha-helical conformation that is stabilized by formation of a hydrophobic staple motif and an N-terminal capping box (Munoz V. Blanco FJ, Serrano L, 1995, Struct Biol 2:380-385). To investigate backbone and side-chain internal motions within the helix and hydrophobic staple, residues F2, A5, L7, A8, and A10 were selectively 13C- and 15N-enriched and NMR relaxation experiments were performed in water and in water/trifluoroethanol (TFE) solution at four Larmor frequencies (62.5, 125, 150, and 200 MHz for 13C). Relaxation data were analyzed using the model free approach and an anisotropic diffusion model. In water, angular variances of motional vectors range from 10 to 20 degrees and backbone phi,psi bond rotations for helix residues A5, L7, A8, and A10 are correlated indicating the presence of Calpha-H, Calpha-Cbeta, and N-H rocking-type motions along the helix dipole axis. L7 side-chain CbetaH2 and CgammaH motions are also correlated and as motionally restricted as backbone CalphaH, suggesting considerable steric hindrance with neighboring groups. In TFE which stabilizes the fold, internal motional amplitudes are attenuated and rotational correlations are increased. For the side chain of hydrophobic staple residue F2, wobbling-in-a-cone type motions dominate in water, whereas in TFE, the Cbeta-Cgamma bond and phenyl ring fluctuate more simply about the Calpha-Cbeta bond. These data support the Daragan-Mayo model of correlated bond rotations (Daragan VA, Mayo KH, 1996, J Phys Chem 100:8378-8388) and contribute to a general understanding of internal motions in peptides and proteins.     

Backbone and side-chain dynamics of residues in a partially folded beta-sheet peptide from platelet factor-4.

Structurally characterizing partially folded states is problematic given the nature of these transient species. A peptide 20mer, T38AQLIATLKNGRKISLDLQA57 (P20), which has been shown to partially fold in a relatively stable turn/loop conformation (LKNGR) and transient beta-sheet structure, is a good model for studying backbone and side-chain mobilities in a transiently folded peptide by using 13C-NMR relaxation. Here, four residues in P20, A43, T44, G48, and 151, chosen for their positions in or near the loop conformation and for compositional variety, have been selectively 13C-enriched. Proton-coupled and decoupled 13C-NMR relaxation experiments have been performed to obtain the temperature dependencies (278 K to 343 K) of auto- and cross-correlation motional order parameters and correlation times. In order to differentiate sequence-neighbor effects from folding effects, two shorter peptides derived from P20, IATLK (P5) and NGRKIS (P6), were similarly 13C-enriched and investigated. For A43, T44, G48, and 151 residues in P20 relative to those in P5/P6, several observations are consistent with partial folding in P20: (1) C alpha H motional tendencies are all about the same, vary less with temperature, and are relatively more restricted, (2) G48 C alpha H2 phi (t) psi (t) rotations are more correlated, and (3) methyl group rotations are slower and yield lower activation energies consistent with formation of hydrophobic "pockets." In addition, T44 and 151 C beta H mobilities in P20 are more restricted at lower temperature than those of their C alpha H and display significantly greater sensitivity to temperature suggesting a larger enthalpic contribution to side-chain mobility. Moreover, at higher temperatures, side-chain methyls and methylenes in P20 are more motionally restricted than those in P5/P6, suggesting that some type of "folded" or "collapsed" structure remains in P20 for what normally would be considered an "unfolded" state.     

13C multiplet nuclear magnetic resonance relaxation-derived ring puckering and backbone dynamics in proline-containing glycine-based peptides.

13CH2-multiplet nuclear magnetic resonance relaxation studies on proline (P)-containing glycine (G)-based peptides, GP, PG, GPG, PGG, and GPGG, provided numerous dipolar auto- and cross-correlation times for various motional model analyses of backbone and proline-ring bond rotations. Molecular dynamics simulations and bond rotation energy profiles were calculated to assess which motions could contribute most to observed relaxation phenomena. Results indicate that proline restricts backbone psi 1, psi 2, and phi 2 motions by 50% relative to those found for a polyglycine control peptide. psi 1 rotations are more restricted in the trans-proline isomer state than in the cis form. A two-state jump model best approximates proline ring puckering which in water could occur either by the C gamma endo-exo or by the C2 interconversion mechanism. The temperature dependence (5 degrees to 75 degrees C) of C beta, and C gamma, and C delta angular changes is rather flat, suggesting a near zero enthalpic contribution to the ring puckering process. In lower dielectric solvents, dimethylsulfoxide and methanol, which may mimic the hydrophobic environment within a protein, the endo-exo mechanism is preferred.     

Structure and dynamics of an amphiphilic peptide in a lipid bilayer: a molecular dynamics study

A molecular dynamics simulation of a simple model membrane system composed of a single amphiphilic helical peptide (ace-K2GL16K2A-amide) in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer was performed for a total of 1060 ps. The secondary structure of the peptide and its stability were described in terms of average dihedral angles, phi and psi, and the C alpha torsion angles formed by backbone atoms; by the average translation per residue along the helix axis; and by the intramolecular peptide hydrogen bonds. The results indicated that residues 6 through 15 remain in a stable right-handed alpha-helical conformation, whereas both termini exhibit substantial fluctuations. A change in the backbone dihedral angles for residues 16 and 17 is accompanied by the loss of two intramolecular hydrogen bonds, leading to a local but long-lived disruption of the helix. The dynamics of the peptide was characterized in terms of local and global helix motions. The local motions of the N-H bond angles were described in terms of the autocorrelation functions of P2[cos thetaNH(t, t + tau)] and reflected the different degrees of local peptide order as well as a variation in time scale for local motions. The chi1 and chi2 dihedral angles of the leucine side chains underwent frequent transitions between potential minima. No connection between the side-chain positions and their mobility was observed, however. In contrast, the lysine side chains displayed little mobility during the simulation. The global peptide motions were characterized by the tilting and bending motions of the helix. Although the peptide was initially aligned parallel to the bilayer normal, during the simulation it was observed to tilt away from the normal, reaching an angle of approximately 25 degrees by the end of the simulation. In addition, a slight bend of the helix was detected. Finally, the solvation of the peptide backbone and side-chain atoms was also investigated.     

Lysine side-chain dynamics derived from 13C-multiplet NMR relaxation studies on di- and tripeptides

Summary ¹³C NMR relaxation data have been used to determine dipolar auto- and cross-correlation times for the di- and tripeptides GK, KG and GKG, primarily to analyze lysine side-chain motional dynamics. In general, correlation times are largest for backbone positions and decrease on going through the lysine side chain, consistent with the idea of increased mobility at C and C methylenes. Correlation times, however, vary with the peptide ionization state. In the zwitterionic state of GK, for example, both auto-and cross-correlation times are at their maximum values, indicating reduced internal motions probably resulting from intramolecular electrostatic interactions. Modifying the charge state increases motional fluctuations. Activation energies determined from the temperature dependence of CH rotational autocorrelation times at neutral pH are approximately equal for glycine and lysine C and lysine C and C positions (4.1±0.2 to 4.5±0.2 kcal/mol) and tend to decrease slightly for lysine C and C (3.8±0.2 to 4.3±0.2 kcal/mol). The sign of lysine side-chain cross-correlations could not be explained by using any available rotational model, including one parameterized for multiple internally restricted rotations and anisotropic overall tumbling. Molecular and stochastic dynamics calculations were performed to obtain insight into correlated internal rotations and coupled overall tumbling and internal motions. Relatively strong correlations were found for i,i+1 backbone and lysine side-chain internal bond rotations. Stochastic dynamics calculations were more successful at explaining experimentally observed correlation times. In the fully charged state, a preferred conformation was detected with an all-trans lysine side chain.Abbreviations rf radio frequency - GK dipeptide glycine-lysine - KG dipeptide lysine-glycine - GKG tripeptide glycine-lysine-glycine     

Conformational effects on tryptophan fluorescence in cyclic hexapeptides

The peptide bond quenches tryptophan fluorescence by excited-state electron transfer, which probably accounts for most of the variation in fluorescence intensity of peptides and proteins. A series of seven peptides was designed with a single tryptophan, identical amino acid composition, and peptide bond as the only known quenching group. The solution structure and side-chain chi(1) rotamer populations of the peptides were determined by one-dimensional and two-dimensional (1)H-NMR. All peptides have a single backbone conformation. The -, psi-angles and chi(1) rotamer populations of tryptophan vary with position in the sequence. The peptides have fluorescence emission maxima of 350-355 nm, quantum yields of 0.04-0.24, and triple exponential fluorescence decays with lifetimes of 4.4-6.6, 1.4-3.2, and 0.2-1.0 ns at 5 degrees C. Lifetimes were correlated with ground-state conformers in six peptides by assigning the major lifetime component to the major NMR-determined chi(1) rotamer. In five peptides the chi(1) = -60 degrees rotamer of tryptophan has lifetimes of 2.7-5.5 ns, depending on local backbone conformation. In one peptide the chi(1) = 180 degrees rotamer has a 0.5-ns lifetime. This series of small peptides vividly demonstrates the dominant role of peptide bond quenching in tryptophan fluorescence.     

Side-chain dynamics of the SAP SH2 domain correlate with a binding hot spot and a region with conformational plasticity

X-linked lymphoproliferative disease is caused by mutations in the protein SAP, which consists almost entirely of a single SH2 domain. SAP interacts with the Tyr281 site of the T<-->B cell signaling protein SLAM via its SH2 domain. Interestingly, binding is not dependent on phosphorylation but does involve interactions with residues N-terminal to the Tyr. We have used 15N and 2H NMR relaxation experiments to investigate the motional properties of the SAP SH2 domain backbone amides and side-chain methyl groups in the free protein and complexes with phosphorylated and non-phosphorylated peptides derived from the Tyr281 site of SLAM. The most mobile methyl groups are in side-chains with large RMSD values between the three crystal structures of SAP, suggesting that fast time-scale dynamics in side-chains is associated with conformational plasticity. The backbone amides of two residues which interact with the C-terminal part of the peptides experience fast time-scale motions in the free SH2 domain that are quenched upon binding of either the phosphorylated or non-phosphorylated peptide. Of most importance, the mobility of methyl groups in and around the binding site for residues in the N-terminus of the peptide is significantly restricted in the complexes, underscoring the dominance of this interaction with SAP and demonstrating a correlation between changes in rapid side-chain motion upon binding with local binding energy.     

Side-chain conformational entropy in protein unfolded states

The largest force disfavoring the folding of a protein is the loss of conformational entropy. A large contribution to this entropy loss is due to the side-chains, which are restricted, although not immobilized, in the folded protein. In order to accurately estimate the loss of side-chain conformational entropy that occurs upon folding it is necessary to have accurate estimates of the amount of entropy possessed by side-chains in the ensemble of unfolded states. A new scale of side-chain conformational entropies is presented here. This scale was derived from Monte Carlo computer simulations of small peptide models. It is demonstrated that the entropies are independent of host peptide length. This new scale has the advantage over previous scales of being more precise with low standard errors. Better estimates are obtained for long (e.g., Arg and Lys) and rare (e.g., Trp and Met) side-chains. Excellent agreement with previous side-chain entropy scales is achieved, indicating that further advancements in accuracy are likely to be small at best. Strikingly, longer side-chains are found to possess a smaller fraction of the theoretical maximum entropy available than short side-chains. This indicates that rotations about torsions after chi(2) are significantly affected by side-chain interactions with the polypeptide backbone. This finding invalidates previous assumptions about side-chain-backbone interactions. Proteins 2000;40:443-450.     

Fluorescence and 13C NMR determination of side-chain and backbone dynamics of synthetic melittin and melittin analogues in isotropic solvents

The dynamics in isotopic solvents of selectively 13C labeled synthetic melittin and three analogues have been investigated by using NMR and fluorescence techniques both separately and in combination. In conjunction with the "model-free" approach to interpretation of NMR relaxation data [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4570], the availability of steady-state fluorescence anisotropy and lifetime data augment T1, T2, and NOE data to provide quantitative information about fluorophore dynamics in these peptides. A method is presented for using combined fluorescence and NMR data to obtain technique- and model-independent values for parameters describing local motion of 13C-labeled fluorophores in peptides and proteins. The dynamics of melittin and melittin analogues are found to be consistent with structural characteristics inferred from CD, fluorescence, and NMR spectral information presented in the preceding paper (Weaver et al., 1989). In particular, the mobility of the random coil peptide monomers is shown to be quite similar, while side-chain as well as peptide backbone motion in the aggregated or oligomeric species differs markedly among the analogues. For melittin itself, experimentally determined overall rotational correlation times for the monomer and tetramer agree very well with values predicted on the basis of solvent-accessible protein surface area. The local dynamics of selectively 13C-labeled Trp-19 and Gly-12 residues of melittin are also found to be consistent with peptide structure. In random coil melittin monomer, a specific model for the motion indicates that the Trp side chain moves through an approximate angle of +/- 71 degrees about the beta-gamma bond with a correlation time of 159 +/- 24 ps. In melittin tetramer, the indole moiety is spatially more confined with a flip angle of +/- 37 degrees, yet demonstrates an increased rate of motion with a correlation time of 56 +/- 8 ps. The constrained mobility of the Trp-19 side chain is consistent with motional constraints inferred from the X-ray structure of melittin tetramer. These results show that protein side-chain motion, even of moieties as large as indole, can occur on the picosecond time scale and that these motions are reasonably similar to those inferred from molecular dynamics simulations.     

Side chain accessibility and dynamics in the molten globule state of alpha-lactalbumin: a (19)F-NMR study

The molten globule state of alpha-lactalbumin (alpha-LA) has been considered a prototype of partially folded proteins. Despite the importance of molten globules in understanding the mechanisms of protein folding and its relevance to some biological phenomena, site-specific information on the structure and dynamics of a molten globule is limited, largely because of the high conformational flexibility and heterogeneity. Here, we use selective isotope labeling and (19)F NMR to investigate the solvent accessibility and side-chain dynamics of aromatic residues in the molten globule of alpha-LA. Comparison of these properties with those of the native and unfolded protein indicates that the alpha-LA molten globule is highly heterogeneous; each residue has its unique solvent accessibility and motional environment. Many aromatic residues normally buried in the interior of native alpha-LA remain significantly buried in the molten globule and the side-chain dynamics of these residues are highly restricted. Our results suggest that hydrophobic and van der Waals interactions mediated by the inaccessible surface area could be sufficient to account for all the stability of the alpha-LA molten globule, which is approximately 50% of the value for the native protein.     

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.     

Effect of backbone cyclization on protein folding stability: chain entropies of both the unfolded and the folded states are restricted

Circular versions of a large number of proteins have been designed by connecting the N and C termini via peptide linkers. A motivation for these designs is the assumed enhancement in folding stability, because backbone cyclization reduces the chain entropy of the unfolded state. Here, it is recognized that backbone cyclization also reduces the chain entropy of a flexible peptide linker in the folded state. Specifically, the end-to-end distance of the linker is restricted to fluctuations around the average displacement between the N and C termini of the folded protein. The balance of the chain-entropy reductions in the folded and unfolded states is used to predict the change in the unfolding free energy, deltadeltaG(cycl), by backbone cyclization. Predicted values of deltadeltaG(cycl) are in quantitative agreement with results of a careful study on cyclizing the 34 residue PIN1 WW domain by linkers with two to seen residues. The experimental results of an optimal linker length l=4 and a maximum stabilization of 1.7 kcal/mol are reproduced. Calculations of deltadeltaG(cycl) for a broad selection of circular proteins suggest that the stabilizing effect of backbone cyclization is modest, reflecting entropy reductions in both the unfolded and the folded states.     

Determinants of membrane association in the SH4 domain of Fyn: Roles of N-terminus myristoylation and side-chain thioacylation

The SH4 domain of Fyn, a member of the Src family of tyrosine kinases, though rich in polar amino acid residues, anchors to the cytosolic face of membranes upon fatty acylation. In order to probe the requirement of specific fatty acylation at the N-terminus and at the side-chain of this domain for membrane-association, we have studied the interaction of peptides corresponding to the polar segment of the SH4 domain of Fyn and its mono- and dually fatty acylated analogs with model membranes. While the polar segment without covalently linked fatty acids (KDKEATKLTEW-amide) does not interact with lipid vesicles, peptides with one covalently linked fatty acid at the N-terminus or in the side-chain, associate with zwitterionic and anionic lipids to varying degrees. The interaction of dually acylated peptides (Myr-GK(ε-myr)KDKEATKLTEW-amide and Myr-GC(S-pal)KDKEATKLTEW-amide) with lipids depends on the linkage between fatty acyl side-chain and peptide backbone. The peptide chain associates with membranes only when the side-chain acylation is via an amide bond and not via a thioester bond. Our investigations indicate that acylation is essential for membrane targeting and unacylated polar stretch of the SH4 domain does not have a role in membrane-anchoring. Side-chain acylation via a thioester bond not only provides membrane anchorage but also directs the peptide chain away from the bilayer which might be important to enable the full length protein to interact with other signaling partners.     

Basis of substrate binding by the chaperonin GroEL.

The molecular chaperonins are essential proteins involved in protein folding, complex assembly, and polypeptide translocation. While there is abundant structural information about the machinery and the mechanistic details of its action are well studied, it is yet unresolved how chaperonins recognize a large number of structurally unrelated polypeptides in their unfolded or partially folded forms. To determine the nature of chaperonin-substrate recognition, we have characterized by NMR methods the interactions of GroEL with synthetic peptides that mimic segments of unfolded proteins. In previous work, we found using transferred nuclear Overhauser effect (trNOE) analysis that two 13 amino acid peptides bound GroEL in an amphipathic alpha-helical conformation. By extending the study to a variety of peptides with differing sequence motifs, we have observed that peptides can adopt conformations other than alpha-helix when bound to GroEL. Furthermore, peptides of the same composition exhibited significantly different affinities for GroEL as manifested by the magnitude of trNOEs. Binding to GroEL correlates well with the ability of the peptide to cluster hydrophobic residues on one face of the peptide, as determined by the retention time on reversed-phase (RP) HPLC. We conclude that the molecular basis of GroEL-substrate recognition is the presentation of a hydrophobic surface by an incompletely folded polypeptide and that many backbone conformations can be accommodated.     

Thermal motions of surface alpha-helices in the D-galactose chemosensory receptor. Detection by disulfide trapping.

The D-galactose chemosensory receptor of Escherichia coli is a .32 kDa globular protein possessing two distinct structural domains, each organized in an alpha/beta folding motif. Helices I and X lie at adjacent approximately parallel positions on the surface of the N-terminal domain, near the hinge region. In order to analyze the relative thermal motions of these two helices, the present study utilizes a generalizable disulfide trapping approach: first, site-directed mutagenesis is used to place a pair of cysteine residues at locations of interest on the protein surface, then disulfide bond formation is used to trap intramolecular cysteine-cysteine collisions resulting from thermal motions. Specifically, four engineered di-cysteine receptors have been constructed, each possessing one cysteine at position 26 on helix I, and a second cysteine at varying positions on helix X. A fifth control receptor possesses one cysteine at position 26, and a second on the opposite surface of the molecule. These surface cysteine substitutions have little or no effect on the measurable receptor parameters as judged by ligand binding equilibria and kinetics, protein stability, and 19F nuclear magnetic resonance, indicating that the engineered receptors are useful probes of native backbone dynamics. Spatial and kinetic features of backbone motions have been investigated by measuring intramolecular disulfide formation rates for cysteine pairs in the fully liganded receptor. The resulting rates decrease monotonically with increasing distance between cysteines in the crystal structure, while no disulfide formation is observed for the control pair unless the molecule is unfolded. The minimum translational amplitudes of the observed backbone motions range from 4.5 to 15.2 A, and the minimum rotational amplitudes are as large as 35 degrees. For each motion the rate of intramolecular sulfhydryl-sulfhydryl collision has been estimated from the measured rate of disulfide formation: the 4.5 and 15.2 A translations yield approximately 10(4) and approximately 10 collisions s-1 molecule-1, respectively. These collision rates, which are faster than ligand dissociation, likely underestimate the actual motional frequencies since only an undetermined fraction of the total motions yield collisions. The simplest plausible trajectory capable of producing such collisions is a rate-limiting translation of one or both helices along their long axes, coupled with minor helix rotations. When sugar is removed from the receptor, a substantial increase in backbone dynamics is observed, indicating the presence of new long-range backbone trajectories. Overall, the results suggest that internal motions in proteins may have larger amplitudes than previously observed.     

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I' beta-turn forming tetrapeptide and an alpha-helical tetradecapeptide.

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I'beta-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt alpha-helical conformations stabilized by 11 successive 5 --> 1 hydrogen bonds. In addition, a single 4 --> 1 hydrogen bond is also observed at the N-terminus. All five Dpg residues adopt backbone torsion angles (phi, psi) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle N-C(alpha)-C' (tau) and the observed backbone phi,psi values. For tau > 106 degrees, helices are observed, while fully extended structures are characterized by tau < 106 degrees. The mean tau values for extended and folded conformations for the Dpg residue are 103.6 degrees +/- 1.7 degrees and 109.9 degrees +/- 2.6 degrees, respectively.     

Statistical and energetic analysis of side-chain conformations in oligopeptides

The distributions of side-chain conformations in 258 crystal structures of oligopeptides have been analyzed. The sample contains 321 residues having side chains that extend beyond the C beta atom. Statistically observed preferences of side-chain dihedral angles are summarized and correlated with stereochemical and energetic constraints. The distributions are compared with observed distributions in proteins of known X-ray structures and with computed minimum-energy conformations of amino acid derivatives. The distributions are similar in all three sets of data, and they appear to be governed primarily by intraresidue interactions. In side chains with no beta-branching, the most important interactions that determine chi 1 are those between the C gamma H2 group and atoms of the neighboring peptide groups. As a result, the g- conformation (chi 1 congruent to -60 degrees) occurs most frequently for rotation around the C alpha-C beta bond in oligopeptides, followed by the t conformation (chi 1 congruent to 180 degrees), while the g+ conformation (chi 1 congruent to 60 degrees) is least favored. In residues with beta-branching, steric repulsions between the C gamma H2 or C gamma H3 groups and backbone atoms govern the distribution of chi 1. The extended (t) conformation is highly favored for rotation around the C beta-C gamma and C gamma-C delta bonds in unbranched side chains, because the t conformer has a lower energy than the g+ and g- conformers in hydrocarbon chains. This study of the observed side-chain conformations has led to a refinement of one of the energy parameters used in empirical conformational energy computations.     

Extending the accuracy limits of prediction for side-chain conformations

Current techniques for the prediction of side-chain conformations on a fixed backbone have an accuracy limit of about 1.0-1.5 A rmsd for core residues. We have carried out a detailed and systematic analysis of the factors that influence the prediction of side-chain conformation and, on this basis, have succeeded in extending the limits of side-chain prediction for core residues to about 0.7 A rmsd from native, and 94 % and 89 % of chi(1) and chi(1+2 ) dihedral angles correctly predicted to within 20 degrees of native, respectively. These results are obtained using a force-field that accounts for only van der Waals interactions and torsional potentials. Prediction accuracy is strongly dependent on the rotamer library used. That is, a complete and detailed rotamer library is essential. The greatest accuracy was obtained with an extensive rotamer library, containing over 7560 members, in which bond lengths and bond angles were taken from the database rather than simply assuming idealized values. Perhaps the most surprising finding is that the combinatorial problem normally associated with the prediction of the side-chain conformation does not appear to be important. This conclusion is based on the fact that the prediction of the conformation of a single side-chain with all others fixed in their native conformations is only slightly more accurate than the simultaneous prediction of all side-chain dihedral angles.     

A peptide model for proline isomerism in the unfolded state of staphylococcal nuclease.

Nuclear magnetic resonance spectroscopy has been used to investigate a synthetic peptide (YVYKPNNTHE) corresponding to residues 113 to 122 of staphylococcal nuclease. In the major folded state of the protein this region forms a type VIa beta-turn containing a cis Lys116-Pro117 peptide bond. There is, however, no evidence for any significant population of such a turn in the peptide in aqueous solution and the X-Pro bond is predominantly in the trans configuration. The peptide exhibits several well-resolved minor resonances due to the presence of a small fraction (4 +/- 2%) of the cis-proline isomer. The ratio of cis to trans isomer populations was found to be independent of temperature between 5 degrees C and 70 degrees C, indicating that delta H for the isomerism is close to zero. Using magnetization transfer techniques the rate of trans to cis interconversion was found to be 0.025(+/- 0.013) s-1 at 50 degrees C. The thermodynamics and kinetics of isomerism in the peptide are very similar to those estimated for the Lys116-Pro117 peptide bond in unfolded nuclease, suggesting that the cis-trans equilibrium in the unfolded protein is largely determined by the residues adjacent to Pro117 in the sequence. These results are consistent with previous suggestions that the cis-proline bond is stabilized late in the folding process and that the predominance of the cis form in folded nuclease is due to stabilizing interactions within the protein that give rise to a favorable enthalpy term.