Angular variances for internal bond rotations of side chains in GXG-based tripeptides derived from (13)C-NMR relaxation measurements: Implications to protein folding

The study of backbone and side-chain internal motions in proteins and peptides is crucial to having a better understanding of protein/peptide "structure" and to characterizing unfolded and partially folded states of proteins and peptides. To achieve this, however, requires establishing a baseline for internal motions and motional restrictions for all residues in the fully, solvent-exposed "unfolded state." GXG-based tripeptides are the simpliest peptides where residue X is fully solvent exposed in the context of an actual peptide. In this study, a series of GXG-based tripeptides has been synthesized with X being varied to include all twenty common amino acid residues. Proton-coupled and -decoupled (13)C-nmr relaxation measurements have been performed on these twenty tripeptides and various motional models (Lipari-Szabo model free approach, rotational anisotropic diffusion, rotational fluctuations within a potential well, rotational jump model) have been used to analyze relaxation data for derivation of angular variances and motional correlation times for backbone and side-chain chi(1) and chi(2) bonds and methyl group rotations. At 298 K, backbone motional correlation times range from about 50 to 85 ps, whereas side-chain motional correlation times show a much broader spread from about 18 to 80 ps. Angular variances for backbone phi,psi bond rotations range from 11 degrees to 23 degrees and those for side chains vary from 5 degrees to 24 degrees for chi(1) bond rotations and from 5 degrees to 27 degrees for chi(2) bond rotations. Even in these peptide models of the "unfolded state," side-chain angular variances can be as restricted as those for backbone and beta-branched (valine, threonine, and isoleucine) and aromatic side chains display the most restricted motions probably due to steric hinderence with backbone atoms. Comparison with motional data on residues in partially folded, beta-sheet-forming peptides indicates that side-chain motions of at least hydrophobic residues are less restricted in the partially folded state, suggesting that an increase in side-chain conformational entropy may help drive early-stage protein folding. Copyright 1999 John Wiley & Sons, Inc.     

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.     

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.     

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.     

Dissecting the stability of a beta-hairpin peptide that folds in water: NMR and molecular dynamics analysis of the beta-turn and beta-strand contributions to folding.

NMR studies of the folding and conformational properties of a beta-hairpin peptide, several peptide fragments of the hairpin, and sequence-modified analogues, have enabled the various contributions to beta-hairpin stability in water to be dissected. Temperature and pH-induced unfolding studies indicate that the folding-unfolding equilibrium approximates to a two-state model. The hairpin is highly resistant to denaturation and is still significantly folded in 7 M urea at 298 K. Thermodynamic analysis shows the hairpin to fold in water with a significant change in heat capacity, however, DeltaCp degrees in 7 M urea is reduced. V/Y-->A mutations on one strand of the hairpin reduce folding to <10 %, consistent with a hydrophobic stabilisation model. We show that in a truncated peptide (residues 6-16) lacking the hydrophobic residues on one beta-strand, the type I' Asn-Gly turn in the sequence SINGKK is significantly populated in water in the absence of interstrand hydrophobic contacts. Unrestrained molecular dynamics simulations of unfolding, using an explicit solvation model, show that the conformation of the NG turn persists for longer than the AG analogue, which has a much lower propensity for type I' turn formation from a data base analysis of preferred turns. The origin of the high stability of the Asn-Gly turn is not entirely clear; data base analysis of 66 NG turns, together with molecular dynamics simulations, reveals no participation of the Asn side-chain in turn-stabilising interactions with the peptide backbone. However, hydration analysis of the molecular dynamics simulations reveals a pocket of "high density" water bridging between the Asn side-chain and peptide main-chain that suggests solvent-mediated interactions may play an important role in modulating phi,psi propensities in the NG turn region.     

A molecular dynamics investigation of the elastomeric restoring force in elastin

A repetitive polypentapeptide organized as a connected chain of beta-bends is believed to be an important structural element of elastin, the major elastomer in biological systems. Molecular dynamics simulations were carried out on hydrated polymers of (Val-Pro-Gly-Val- Gly)18 at various extensions. Analysis of the fluctuations of backbone angles in relaxed elastin showed that particularly large-amplitude torsional motions occur in phi and psi angles of residues connecting sequentially adjacent hairpin bends. Many such motions reflect peptide plane librations that result from anticorrelated crankshaft rotations of psi i and phi i+1. These effects were much reduced in stretched polymer models. The conformational entropy of relaxed and stretched elastin models was estimated using a treatment due to Meirovitch, and gave a calculated decrease in entropy of about 1 cal/mol deg when the polymer was stretched to 1.75 times its original length. There are large changes in solvent-accessible surface area during the initial stages of elastin stretching. Collectively these results suggest that hydrophobic interactions make contributions to elastin entropy at low extensions, but that librational mechanisms make larger contributions to the elastic restoring force at longer extensions.     

Polypeptide motions are dominated by peptide group oscillations resulting from dihedral angle correlations between nearest neighbors

To identify basic local backbone motions in unfolded chains, simulations are performed for a variety of peptide systems using three popular force fields and for implicit and explicit solvent models. A dominant "crankshaft-like" motion is found that involves only a localized oscillation of the plane of the peptide group. This motion results in a strong anticorrelated motion of the phi angle of the ith residue (phi(i)) and the psi angle of the residue i - 1 (psi(i-1)) on the 0.1 ps time scale. Only a slight correlation is found between the motions of the two backbone dihedral angles of the same residue. Aside from the special cases of glycine and proline, no correlations are found between backbone dihedral angles that are separated by more than one torsion angle. These short time, correlated motions are found both in equilibrium fluctuations and during the transit process between Ramachandran basins, e.g., from the beta to the alpha region. A residue's complete transit from one Ramachandran basin to another, however, occurs in a manner independent of its neighbors' conformational transitions. These properties appear to be intrinsic because they are robust across different force fields, solvent models, nonbonded interaction routines, and most amino acids.     

The utility of side-chain cyclization in determining the receptor-bound conformation of peptides: cyclic tripeptides and angiotensin II.

The effect of side-chain cyclization on accessible backbone conformations of tripeptides, X-Ala-Y (X and/or Y = Cys, Hcy (Hcy: homocysteine), cis 4-mercaptoproline (MPc), and trans 4-mercaptoproline (MPt)), was elucidated using two variants of systematic conformational search. In addition to cyclization through a disulfide bond, the thioether (-S-CH2-) and amide (-CO-NH-) side-chain analogues of Cys-Ala-Cys and Hcy-Ala-Hcy were evaluated. The number of valid backbone conformations and the allowed phi, psi space were evaluated for each compound, and the ability of the cyclic tripeptides to accommodate beta-turn conformations was examined in order to assess the value of cyclization in limiting conformational freedom. Based on the number of conformations, cyclization was highly effective in reducing the backbone degree of freedom: in order of decreasing number of conformations, Ala-Ala-Ala 1 > Hcy-Ala-Hcy 2 > Cys-Ala-Hcy 3 approximately equal to Hcy-Ala-Cys 4 > MPc-Ala-Hcy 5, 7 > Cys-Ala-Cys 6 > MPc-Ala-Cys 8 > Hcy-Ala-MPt 9 > Cys-Ala-MPt 10 approximately equal to MPc-Ala-MPt 11. Although Hcy-Ala-Hcy 2 had the greatest number of conformations of the cyclic peptides studied, it was still greatly constrained relative to its linear analogue 1. The bicyclic ring system introduced by MP was even more effective in constraining the cycle, having greater impact at position 3 than at position 1. Under the conditions of the study, cyclization of MP-containing analogues could be effected only with the cis isomer (MPc) at position 1 and/or the trans isomer (MPt) at position 3. Sterically allowed conformations of Ala2 for the cyclic tripeptides 2-4 were generally similar to those of the linear tripeptide 1, while those of Cys-Ala-Cys 6 and MPc-Ala-Hcy 7 were restricted to a smaller region of phi 2, psi 2 space: the right- and left-handed alpha-helical conformation and the beta-conformation. This trend was even more pronounced for Hcy-Ala-MPt 9, Cys-Ala-MPt 10, and MPc-Ala-MPt 11, in which Ala2 was severely restricted to a very small region of phi, psi space: the left-handed alpha-helical conformation for 9-11, plus the beta conformation for 9. This suggests that MP at the 3-position is incompatible with a right-handed alpha-helical conformation at position 2.(ABSTRACT TRUNCATED AT 400 WORDS)     

Side-chain flexibility in protein-ligand binding: the minimal rotation hypothesis

The goal of this work is to learn from nature about the magnitudes of side-chain motions that occur when proteins bind small organic molecules, and model these motions to improve the prediction of protein-ligand complexes. Following analysis of protein side-chain motions upon ligand binding in 63 complexes, we tested the ability of the docking tool SLIDE to model these motions without being restricted to rotameric transitions or deciding which side chains should be considered as flexible. The model tested is that side-chain conformational changes involving more atoms or larger rotations are likely to be more costly and less prevalent than small motions due to energy barriers between rotamers and the potential of large motions to cause new steric clashes. Accordingly, SLIDE adjusts the protein and ligand side groups as little as necessary to achieve steric complementarity. We tested the hypothesis that small motions are sufficient to achieve good dockings using 63 ligands and the apo structures of 20 different proteins and compared SLIDE side-chain rotations to those experimentally observed. None of these proteins undergoes major main-chain conformational change upon ligand binding, ensuring that side-chain flexibility modeling is not required to compensate for main-chain motions. Although more frugal in the number of side-chain rotations performed, this model substantially mimics the experimentally observed motions. Most side chains do not shift to a new rotamer, and small motions are both necessary and sufficient to predict the correct binding orientation and most protein-ligand interactions for the 20 proteins analyzed.     

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.     

A shorter peptide model from staphylococcal nuclease for the folding-unfolding equilibrium of a beta-hairpin shows that unfolded state has significant contribution from compact conformational states

It is important to understand the conformational features of the unfolded state in equilibrium with folded state under physiological conditions. In this paper, we consider a short peptide model LMYKGQPM from staphylococcal nuclease to model the conformational equilibrium between a hairpin conformation and its unfolded state using molecular dynamics simulation under NVT conditions at 300K using GROMOS96 force field. The free energy landscape has overall funnel-like shape with hairpin conformations sampling the minima. The "unfolded" state has a higher free energy of approximately 12kJ/mol with respect to native hairpin minimum and occupies a plateau region. We find that the unfolded state has significant contributions from compact conformations. Many of these conformations have hairpin-like topology. Further, these compact conformational forms are stabilized by hydrophobic interactions. Conversion between native and non-native hairpins occurs via unfolded states. Frequent conversions between folded and unfolded hairpins are observed with single exponential kinetics. We compare our results with the emerging picture of unfolded state from both experimental and theoretical studies.     

Autonomous folding of a peptide corresponding to the N-terminal beta-hairpin from ubiquitin.

The N-terminal 17 residues of ubiquitin have been shown by 1H NMR to fold autonomously into a beta-hairpin structure in aqueous solution. This structure has a specific, native-like register, though side-chain contacts differ in detail from those observed in the intact protein. An autonomously folding hairpin has previously been identified in the case of streptococcal protein G, which is structurally homologous with ubiquitin, but remarkably, the two are not in topologically equivalent positions in the fold. This suggests that the organization of folding may be quite different for proteins sharing similar tertiary structures. Two smaller peptides have also been studied, corresponding to the isolated arms of the N-terminal hairpin of ubiquitin, and significant differences from simple random coil predictions observed in the spectra of these subfragments, suggestive of significant limitation of the backbone conformational space sampled, presumably as a consequence of the strongly beta-structure favoring composition of the sequences. This illustrates the ability of local sequence elements to express a propensity for beta-structure even in the absence of actual sheet formation. Attempts were made to estimate the population of the folded state of the hairpin, in terms of a simple two-state folding model. Using published "random coil" values to model the unfolded state, and values derived from native ubiquitin for the putative unique, folded state, it was found that the apparent population varied widely for different residues and with different NMR parameters. Use of the spectra of the subfragment peptides to provide a more realistic model of the unfolded state led to better agreement in the estimates that could be obtained from chemical shift and coupling constant measurements, while making it clear that some other approaches to population estimation could not give meaningful results, because of the tendency to populate the beta-region of conformational space even in the absence of the hairpin structure.     

A survey of left-handed polyproline II helices

Left-handed polyproline II helices (PPII) are contiguous elements of protein secondary structure in which the phi and psi angles of constituent residues are restricted to around -75 degrees and 145 degrees, respectively. They are important in structural proteins, in unfolded states and as ligands for signaling proteins. Here, we present a survey of 274 nonhomologous polypeptide chains from proteins of known structure for regions that form these structures. Such regions are rare, but the majority of proteins contain at least one PPII helix. Most PPII helices are shorter than five residues, although the longest found contained 12 amino acids. Proline predominates in PPII, but Gln and positively charged residues are also favored. The basis of Gln's prevalence is its ability to form an i, i + 1 side-chain to main-chain hydrogen bond with the backbone carbonyl oxygen of the proceeding residue; this helps to fix the psi angle of the Gln and the phi and psi of the proceeding residue in PPII conformations and explains why Gln is favored at the first position in a PPII helix. PPII helices are highly solvent exposed, which explains why apolar amino acids are disfavored despite preferring this region of phi/psi space when in isolation. PPII helices have perfect threefold rotational symmetry and within these structures we find significant correlation between the hydrophobicity of residues at i and i + 3; thus, PPII helices in globular proteins can be considered to be amphipathic.     

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.     

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.     

Application of cross-correlated NMR spin relaxation to the zinc-finger protein CRP2(LIM2): evidence for collective motions in LIM domains

The solution structure of quail CRP2(LIM2) was significantly improved by using an increased number of NOE constraints obtained from a 13C,15N-labeled protein sample and by applying a recently developed triple-resonance cross-correlated relaxation experiment for the determination of the backbone dihedral angle psi. Additionally, the relative orientation of the 15N(i)-1HN(i) dipole and the 13CO(i) CSA tensor, which is related to both backbone angles phi and psi, was probed by nitrogen-carbonyl multiple-quantum relaxation and used as an additional constraint for the refinement of the local geometry of the metal-coordination sites in CRP2(LIM2). The backbone dynamics of residues located in the folded part of CRP2(LIM2) have been characterized by proton-detected 13C'(i-1)-15N(i) and 15N(i)-1HN(i) multiple-quantum relaxation, respectively. We show that regions having cross-correlated time modulation of backbone isotropic chemical shifts on the millisecond to microsecond time scale correlate with residues that are structurally altered in the mutant protein CRP2(LIM2)R122A (disruption of the CCHC zinc-finger stabilizing side-chain hydrogen bond) and that these residues are part of an extended hydrogen-bonding network connecting the two zinc-binding sites. This indicates the presence of long-range collective motions in the two zinc-binding subdomains. The conformational plasticity of the LIM domain may be of functional relevance for this important protein recognition motif.     

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     

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.     

Serine and threonine residues bend alpha-helices in the chi(1) = g(-) conformation

The relationship between the Ser, Thr, and Cys side-chain conformation (chi(1) = g(-), t, g(+)) and the main-chain conformation (phi and psi angles) has been studied in a selection of protein structures that contain alpha-helices. The statistical results show that the g(-) conformation of both Ser and Thr residues decreases their phi angles and increases their psi angles relative to Ala, used as a control. The additional hydrogen bond formed between the O(gamma) atom of Ser and Thr and the i-3 or i-4 peptide carbonyl oxygen induces or stabilizes a bending angle in the helix 3-4 degrees larger than for Ala. This is of particular significance for membrane proteins. Incorporation of this small bending angle in the transmembrane alpha-helix at one side of the cell membrane results in a significant displacement of the residues located at the other side of the membrane. We hypothesize that local alterations of the rotamer configurations of these Ser and Thr residues may result in significant conformational changes across transmembrane helices, and thus participate in the molecular mechanisms underlying transmembrane signaling. This finding has provided the structural basis to understand the experimentally observed influence of Ser residues on the conformational equilibrium between inactive and active states of the receptor, in the neurotransmitter subfamily of G protein-coupled receptors.     

Entrapping unusual folding characteristics of the beta-Ala residues in a model peptide: Boc-beta-Ala-Aib-beta-Ala-NHCH3.

Crystal structure analysis of a model peptide: Boc-beta-Ala-Aib-beta-Ala-NHCH3 (beta-Ala: 3-amino propionic acid; Aib: alpha-aminoisobutyric acid) revealed distinct conformational preferences for folded [phi approximately 136 degrees, mu approximately -62 degrees, psi approximately 100 degrees] and semifolded [phi approximately 83 degrees, mu approximately -177 degrees, psi approximately -117 degrees] structures of the N-and C-terminus beta-Ala residues, respectively. The overall folded conformation is stabilized by unusual Ni...H-Ni+1 and nonconventional C-H...O intramolecular hydrogen bonding interactions.