Concordance of residual dipolar couplings, backbone order parameters and crystallographic B-factors for a small alpha/beta protein: a unified picture of high probability, fast atomic motions in proteins

Using ensemble refinement of the third immunoglobulin binding domain (GB3) of streptococcal protein G (a small alpha/beta protein of 56 residues), we demonstrate that backbone (N-H, N-C', Calpha-Halpha, Calpha-C') residual dipolar coupling data in five independent alignment media, generalized order parameters from 15N relaxation data, and B-factors from a high-resolution (1.1A), room temperature crystal structure are entirely consistent with one another within experimental error. The optimal ensemble size representation is between four and eight, as assessed by complete cross-validation of the residual dipolar couplings. Thus, in the case of GB3, all three observables reflect the same low-amplitude anisotropic motions arising from fluctuations in backbone phi/psi torsion angles in the picosecond to nanosecond regime in both solution and crystalline environments, yielding a unified picture of fast, high-probability atomic motions in proteins. An understanding of these motions is crucial for understanding the impact of protein dynamics on protein function, since they provide part of the driving force for triggered conformational changes that occur, for example, upon ligand binding, signal transduction and enzyme catalysis.     

Self-consistent residual dipolar coupling based model-free analysis for the robust determination of nanosecond to microsecond protein dynamics

Residual dipolar couplings (RDCs) provide information about the dynamic average orientation of inter-nuclear vectors and amplitudes of motion up to milliseconds. They complement relaxation methods, especially on a time-scale window that we have called supra-tau(c) (tau(c) < supra-tau(c) < 50 micros). Here we present a robust approach called Self-Consistent RDC-based Model-free analysis (SCRM) that delivers RDC-based order parameters-independent of the details of the structure used for alignment tensor calculation-as well as the dynamic average orientation of the inter-nuclear vectors in the protein structure in a self-consistent manner. For ubiquitin, the SCRM analysis yields an average RDC-derived order parameter of the NH vectors 0.72 +/- 0.02 compared to = 0.778 +/- 0.003 for the Lipari-Szabo order parameters, indicating that the inclusion of the supra-tau(c) window increases the averaged amplitude of mobility observed in the sub-supra-tau(c) window by about 34%. For the beta-strand spanned by residues Lys48 to Leu50, an alternating pattern of backbone NH RDC order parameter S2(rdc)(NH) = (0.59, 0.72, 0.59) was extracted. The backbone of Lys48, whose side chain is known to be involved in the poly-ubiquitylation process that leads to protein degradation, is very mobile on the supra-tau(c) time scale (S2(rdc)(NH) = 0.59 +/- 0.03), while it is inconspicuous (S2(LS)(NH)= 0.82) on the sub-tau(c) as well as on micros-ms relaxation dispersion time scales. The results of this work differ from previous RDC dynamics studies of ubiquitin in the sense that the results are essentially independent of structural noise providing a much more robust assessment of dynamic effects that underlie the RDC data.     

Structural characterization of unfolded states of apomyoglobin using residual dipolar couplings

The conformational propensities of unfolded states of apomyoglobin have been investigated by measurement of residual dipolar couplings between (15)N and (1)H in backbone amide groups. Weak alignment of apomyoglobin in acid and urea-unfolded states was induced with both stretched and compressed polyacrylamide gels. In 8 M urea solution at pH 2.3, conditions under which apomyoglobin contains no detectable secondary or tertiary structure, significant residual dipolar couplings of uniform sign were observed for all residues. At pH 2.3 in the absence of urea, a change in the magnitude and/or sign of the residual dipolar couplings occurs in local regions of the polypeptide where there is a high propensity for helical secondary structure. These results are interpreted on the basis of the statistical properties of the unfolded polypeptide chain, viewed as a polymer of statistical segments. For a folded protein, the magnitude and sign of the residual dipolar couplings depend on the orientation of each bond vector relative to the alignment tensor of the entire molecule, which reorients as a single entity. For unfolded proteins, there is no global alignment tensor; instead, residual dipolar couplings are attributed to alignment of the statistical segments or of transient elements of secondary structure. For apomyoglobin in 8 M urea, the backbone is highly extended, with phi and psi dihedral angles favoring the beta or P(II) regions. Each statistical segment has a highly anisotropic shape, with the N-H bond vectors approximately perpendicular to the long axis, and becomes weakly aligned in the anisotropic environment of the strained acrylamide gels. Local regions of enhanced flexibility or chain compaction are characterized by a decrease in the magnitude of the residual dipolar couplings. The formation of a small population of helical structure in the acid-denatured state of apomyoglobin leads to a change in sign of the residual dipolar couplings in local regions of the polypeptide; the population of helix estimated from the residual dipolar couplings is in excellent agreement with that determined from chemical shifts. The alignment model described here for apomyoglobin can also explain the pattern of residual dipolar couplings reported previously for denatured states of staphylococcal nuclease and other proteins. In conjunction with other NMR experiments, residual dipolar couplings can provide valuable insights into the dynamic conformational propensities of unfolded and partly folded states of proteins and thereby help to chart the upper reaches of the folding landscape.     

Fast internal main-chain dynamics of human ubiquitin.

The fast internal dynamics of human ubiquitin have been studied by the analysis of 15N relaxation of backbone amide nitrogens. The amide 15N resonances have been assigned by use of heteronuclear multiple-quantum spectroscopy. Spin lattice relaxation times at 60.8 and 30.4 MHz and the steady-state nuclear Overhauser effect at 60.8 MHz have been determined for 67 amide 15N sites in the protein using two-dimensional spectroscopy. These data have been analyzed in terms of the model free treatment of Lipari and Szabo [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4559]. The global motion of the protein is shown to be isotropic and is characterized by a correlation time of 4.1 ns rad-1. The generalized order parameters (S2) of backbone amide N-H vectors in the globular region of the protein range from 0.5 to 0.95. No apparent correlation between secondary structure and generalized order parameters is observed. There is, however, a strong correlation between the magnitude of the generalized order parameters of a given N-H vector and the presence of hydrogen bonding of the amide hydrogen or its peptide bond associated carbonyl. Using a chemical shift tensor breadth of 160 ppm, the N-H vectors of peptide linkages participating in one or more hydrogen bonds to the main chain show an average generalized order parameter of 0.80 (SD 0.06), while those amide NH of peptide linkages free of hydrogen-bonding interactions with the main chain show an average order parameter of 0.69 (SD 0.06).(ABSTRACT TRUNCATED AT 250 WORDS)     

Protein backbone motions viewed by intraresidue and sequential HN–Hα residual dipolar couplings

Triple resonance E.COSY-based techniques were used to measure intra-residue and sequential H(N)-H(alpha) residual dipolar couplings (RDCs) for the third IgG-binding domain of protein G (GB3), aligned in Pf1 medium. Measurements closely correlate with values predicted on the basis of an NMR structure, previously determined on the basis of a large number of one-bond backbone RDCs measured in five alignment media. However, in particular the sequential H(N)-H(alpha) RDCs are smaller than predicted for a static structure, suggesting a degree of motion for these internuclear vectors that exceeds that of the backbone amide N-H vectors. Of all experimentally determined GB3 structures available, the best correlation between experimental (1)H-(1)H couplings is observed for a GB3 ensemble, previously derived to generate a realistic picture of the conformational space sampled by GB3 (Clore and Schwieters, J Mol Biol 355:879-886, 2006). However, for both NMR and X-ray-derived structures the (1)H-(1)H couplings are found to be systematically smaller than expected on the basis of alignment tensors derived from (15)N-(1)H amide RDCs, assuming librationally corrected N-H bond lengths of 1.041 A.     

Main chain and side chain dynamics of oxidized flavodoxin from Cyanobacterium anabaena.

Oxidized flavodoxin from Cyanobacterium anabaena PCC 7119 is used as a model system to investigate the fast internal dynamics of a flavin-bearing protein. Virtually complete backbone and side chain resonance NMR assignments of an oxidized flavodoxin point mutant (C55A) have been determined. Backbone and side chain dynamics in flavodoxin (C55A) were investigated using (15)N amide and deuterium methyl NMR relaxation methods. The squared generalized order parameters (S(NH)(2)) for backbone amide N-H bonds are found to be uniformly high ( approximately 0.923 over 109 residues in regular secondary structure), indicating considerable restriction of motion in the backbone of the protein. In contrast, methyl-bearing side chains are considerably heterogeneous in their amplitude of motion, as indicated by obtained symmetry axis squared generalized order parameters (S(axis)(2)). However, in comparison to nonprosthetic group-bearing proteins studied with these NMR relaxation methods, the side chains of oxidized flavodoxin are unusually rigid.     

Influence of the fluctuations of the alignment tensor on the analysis of the structure and dynamics of proteins using residual dipolar couplings

It has been suggested that the fluctuations of the alignment tensor can affect the results of procedures for characterizing the structure and the dynamics of proteins using residual dipolar couplings. We show here that the very significant fluctuations of the steric alignment tensor caused by the dynamics of proteins can be safely ignored when they do not correlate with those of the bond vectors. A detailed analysis of these correlations in the protein ubiquitin reveals that their effects are negligible for the analysis of backbone motions within secondary structure elements, but also that they may be significant in turns, loops and side chains, especially for bond vectors that have small residual dipolar couplings. Our results suggest that methods that explicitly consider the motions of the alignment tensor will be needed to study the large-scale structural fluctuations that take place on the millisecond timescale, which are often important for the biological function of proteins, from residual dipolar coupling measurements.     

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 dynamics and thermodynamics of Borrelia outer surface protein A

Nuclear spin relaxation experiments performed at 298K, 308K and 318K are used to characterize the intramolecular dynamics and thermodynamics of outer surface protein A (OspA), a key protein in the life-cycle of Borrelia burgdorferi, the causative agent of Lyme disease. It has recently been demonstrated that OspA specifically binds to the gut of the intermediate tick host (Ixodes scapularis), and that this interaction is mediated, at least in part, by residues in the C-terminal domain of OspA that are largely inaccessible to solvent in all X-ray structures of this protein. Our analysis of 15N relaxation parameters in OspA shows that the putative-binding region contains and is surrounded by flexible residues, which could facilitate accessibility to solvent and ligands. In addition, residues with similar activation energies are clustered in a manner that suggests locally collective motions. We have used molecular modeling to show that these collective motions are consistent with a hinge-bending mechanism that exposes residues implicated in binding. Characteristic temperatures describing the energy landscape of the OspA backbone are derived from the temperature dependence of the N-H bond vector order parameters, and a comparison is made between the N and C-terminal globular domains and the unusual single-layer beta-sheet connecting them. The average characteristic temperatures in the three regions indicate that, with an increase in temperature, a larger increase in accessible conformational states occurs for N-H bond vectors in the single-layer central beta-sheet than for bond vectors in the globular N and C-terminal domains. These conformational states are accessible without disruption of hydrogen bonds, providing a conformational entropic gain, upon increase in temperature, without a significant enthalpic penalty. This increase in heat capacity may help to explain the unexpected thermal stability of the unusual single-layer beta-sheet.     

NMR assignments, secondary structure, and global fold of calerythrin, an EF-hand calcium-binding protein from Saccharopolyspora erythraea

Calerythrin is a 20 kDa calcium-binding protein isolated from gram-positive bacterium Saccharopolyspora erythraea. Based on amino acid sequence homology, it has been suggested that calerythrin belongs to the family of invertebrate sarcoplasmic EF-hand calcium-binding proteins (SCPs), and therefore it is expected to function as a calcium buffer. NMR spectroscopy was used to obtain structural information on the protein in solution. Backbone and side chain 1H, 13C, and 15N assignments were obtained from triple resonance experiments HNCACB, HN(CO)CACB, HNCO, CC(CO)NH, and [15N]-edited TOCSY, and HCCH-TOCSY. Secondary structure was determined by using secondary chemical shifts and characteristic NOEs. In addition, backbone N-H residual dipolar couplings were measured from a spin-state selective [1H, 15N] correlation spectrum acquired from a sample dissolved in a dilute liquid crystal. Four EF-hand motifs with characteristic helix-loop-helix patterns were observed. Three of these are typical calcium-binding EF-hands, whereas site 2 is an atypical nonbinding site. The global fold of calerythrin was assessed by dipolar couplings. Measured dipolar couplings were compared with values calculated from four crystal structures of proteins with sequence homology to calerythrin. These data allowed us to recognize an overall similarity between the folds of calerythrin and sarcoplasmic calcium-binding proteins from the sandworm Nereis diversicolor and the amphioxus Branchiostoma lanceolatum.     

Calculations of NMR dipolar coupling strengths in model peptides

Ab initio MP2 and density functional quantum chemistry calculations are used to explore geometries and vibrational properties of N-methylacetamide and of the alanine dipeptide with backbone angles characteristic of helix and sheet regions in proteins. The results are used to explore one-bond direct dipolar couplings for the N–H, C–H, C–N, and C–C bonds, as well as for the two-bond C–H interaction. Vibrational averaging affects these dipolar couplings, and these effects can be expressed as effective bond lengths that are 0.5–3% larger than the true bond lengths; bending and torsion vibrations have a bigger influence on the effective coupling than do stretching vibrations. Because of zero-point motion, these effects are important even at low temperature. Hydrogen bonding interactions at the amide group also increase the N-H effective bond length. Although vibrational contributions to effective bond lengths are small, they can have a significant influence on the extraction of order parameters from relaxation data, and a knowledge of relative bond lengths is needed when several types of dipolar couplings are to be simultaneously used for refinement. The present computational results are compared to both solid- and liquid-state NMR experiments. The analysis suggests that secondary structural elements in many proteins may be more rigid than is commonly thought.     

Rotational diffusion anisotropy of proteins from simultaneous analysis of 15N and 13Cα nuclear spin relaxation

Current methods of determining the rotational diffusion tensors of proteins in solution byNMR spectroscopy exclusively utilize relaxation rate constants for backbone amide 15N spins.However, the distributions of orientations of N-H bond vectors are not isotropic in manyproteins, and correlations between bond vector orientations reduce the accuracy and precisionof rotational diffusion tensors extracted from 15N spin relaxation data. The inclusion of both13C and 15N spin relaxation rate constants increases the robustness of the diffusiontensor analysis because the orientations of the C-H bond vectors differ from theorientations of the N-H bond vectors. Theoretical and experimental results for calbindin D9k,granulocyte colony stimulating factor, and ubiquitin, three proteins with different distributionsof N-H and C-H bond vectors, are used to illustrate the advantages of thesimultaneous utilization of 13C and 15N relaxation data.     

Mapping of the spectral densities of N-H bond motions in eglin c using heteronuclear relaxation experiments.

A new strategy is used for studying the internal motions of proteins based on measurements of NMR relaxation parameters. The strategy yields values of the so-called spectral density functions J(omega) for N-H bond vectors. The spectral density functions are related to the distribution of frequencies contained in the rotational (overall and internal) motions of these NH bond vectors. No a priori model assumptions about the dynamics are required in this approach. The method involves measurements of six relaxation parameters consisting of 15N longitudinal relaxation rates, transverse relaxation rates of in-phase and antiphase coherence, the relaxation rates of heteronuclear 1H-15N two-spin order, the heteronuclear 1H-15N nuclear Overhauser effects, and longitudinal relaxation rates of the amide protons. The values of the spectral density functions at the five frequencies 0, omega N, omega H + omega N, omega H, and omega H - omega N are determined from the relaxation parameters using analytical relations derived previously [Peng & Wagner (1992) J. Magn. Reson. 98, 308-332]. Here, the method is applied to characterize the backbone dynamics of the 15N-enriched proteinase inhibitor eglin c, a protein of 70 residues. The values for J(0) and J(omega N = 50 MHz) vary significantly with the amino acid sequence, whereas the spectral densities at higher frequencies, J(450 MHz), J(500 MHz), and J(550 MHz), are typically much smaller and show no significant variation with the sequence. The collective behavior of the J(omega) values indicate greater internal motion for the proteinase binding loop residues and the first eight N-terminal residues. The additional internal motion in these regions is in the rate range below 450 MHz. The values of J(omega) are also compared with root mean square deviations (rmsds) of backbone atoms as obtained in NMR structure determinations. Low values of J(0) and J(omega N) are correlated with high rmsds. Spectral densities at higher frequencies, J(450 MHz), J(500 MHz), and J(550 MHz), are small and show no correlation with rmsds. A comparison with the spectral density functions obtained by fitting the experimental data to the functional dependence of the Lipari and Szabo formalism [Lipari & Szabo (1982a) J. Am. Chem. Soc. 104, 4546-4559] is made.     

Characterization of the overall and local dynamics of a protein with intermediate rotational anisotropy: Differentiating between conformational exchange and anisotropic diffusion in the B3 domain of protein G

Because the overall tumbling provides a major contribution to protein spectral densities measured in solution, the choice of a proper model for this motion is critical for accurate analysis of protein dynamics. Here we study the overall and backbone dynamics of the B3 domain of protein G using ¹⁵N relaxation measurements and show that the picture of local motions is markedly dependent on the model of overall tumbling. The main difference is in the interpretation of the elevated R ₂ values in the -helix: the isotropic model results in conformational exchange throughout the entire helix, whereas no exchange is predicted by anisotropic models that place the longitudinal axis of diffusion tensor almost parallel to the helix axis. Due to small size (fast tumbling) of the protein, the T ₁ values have low sensitivity to NH bond orientation. The diffusion tensor derived from orientation dependence of R ₂/R ₁ is anisotropic (D _par/D _perp=1.4), with a small rhombic component. In order to distinguish the correct picture of motion, we apply model-independent methods that are sensitive to conformational exchange and do not require knowledge of protein structure or assumptions about its dynamics. A comparison of the CSA/dipolar cross-correlation rate constants with ¹⁵N relaxation rates and the estimation of R _ex terms from relaxation data at 9.4 and 14.1 T indicate no conformational exchange in the helix, in support of the anisotropic models. The experimentally derived diffusion tensor is in excellent agreement with theoretical predictions from hydrodynamic calculations; a detailed comparison with various hydrodynamic models revealed optimal parameters for hydrodynamic calculations.     

Molecular motions and dynamics of a diunsaturated acyl chain in a lipid bilayer: implications for the role of polyunsaturation in biological membranes.

The nature and dynamics of the motions of a diunsaturated fatty acyl chain in a lipid bilayer were examined using a comprehensive simulation program for 2H NMR line shapes developed by Wittebort et al. [Wittebort, R. J., Olejniczak, E. T., & Griffin, R. G. (1987) J. Chem. Phys. 36, 5411-5420]. A motional model in which the isolinoleoyl chain (18:2 delta 6,9) adopts two conformations consistent with the low energy structures proposed for 1,4-pentadiene [Applegate, K. R., & Glomset, J. A. (1986) J. Lipid Res. 27, 658-680], but undergoes a rapid jump between these states, is sufficient to account for the experimentally observed quadrupolar couplings, the 2H-2H and 1H-2H dipolar couplings, the longitudinal relaxation times, and the changes in the average conformation of the chain that occur with a variation in temperature. The jump motion originates via rotations about the C7-C8 and the C8-C9 carbon bonds and leads to the low order parameters assigned to the C8 methylene segment (0.18) and the C9-C10 double bond (0.11). In contrast, the C6-C7 double bond, which is not involved in the two-site jump, characterized by a relatively large order parameter (0.56). Fatty acyl chains containing three or more double bonds likely cannot undergo the same jump motion and consequently will be highly ordered structures. Correlation times for diffusion of the molecular long axis of the diunsaturated acyl chain about the bilayer normal (approximately 10(-10) s) and for the local jump motion (approximately 10(-10) s) were calculated.(ABSTRACT TRUNCATED AT 250 WORDS)     

Effects of denaturants and substitutions of hydrophobic residues on backbone dynamics of denatured staphylococcal nuclease

Analysis of residual dipolar couplings (RDCs) in the Delta131Delta fragment of staphylococcal nuclease has demonstrated that its ensemble-averaged structure is resistant to perturbations such as high concentrations of urea, low pH, and substitution of hydrophobic residues, suggesting that its residual structure is encoded by local side-chain/backbone interactions. In the present study, the effects of these same perturbations on the backbone dynamics of Delta131Delta were examined through (1)H-(15)N relaxation methods. Unlike the global structure reported by RDCs, the transverse relaxation rates R(2) were quite sensitive to denaturing conditions. At pH 5.2, Delta131Delta exhibits an uneven R(2) profile with several characteristic peaks involving hydrophobic chain segments. Protonation of carboxyl side chains by lowering the pH reduces the values of R(2) along the entire chain, yet these characteristic peaks remain. In contrast, high concentrations of urea or the substitution of 10 hydrophobic residues eliminates these peaks and reduces the R(2) values by a greater amount. The combination of low pH and high urea leads to further decreases in R(2). These denaturant-induced increases in backbone mobility are also reflected in decreases in (15)N NOEs and in relaxation interference parameters, with the former reporting an increase in fast motions and the latter a decrease in slow motions. Comparison between the changes in chain dynamics and the corresponding changes in Stokes radius and the patterns of RDCs suggests that regional variations in backbone dynamics in denatured nuclease arise primarily from local contacts between hydrophobic side chains and local interactions involving charged carboxyl groups.     

Comparison of fast backbone dynamics at amide nitrogen and carbonyl sites in dematin headpiece C-terminal domain and its S74E mutant

We perform a detailed comparison of fast backbone dynamics probed at amide nitrogen versus carbonyl carbon sites for dematin headpiece C-terminal domain (DHP) and its S74E mutant (DHPS74E). Carbonyl dynamics is probed via auto-correlated longitudinal rates and transverse C′/C′-C^α CSA/dipolar and C′/C′–N CSA/dipolar cross-correlated rates, while ¹⁵N data are taken from a previous study. Resulting values of effective order parameters and internal correlation times support the conclusion that C′ relaxation reports on a different subset of fast motions compared to those probed at N–H bond vectors in the same peptide planes. ¹³C′ order parameters are on the average 0.08 lower than ¹⁵N order parameters with the exception of the flexible loop region in DHP. The reduction of mobility in the loop region upon the S74E mutation can be seen from the ¹⁵N order parameters but not from the ¹³C order parameters. Internal correlation times at ¹³C′ sites are on the average an order of magnitude longer than those at ¹⁵N sites for the well-structured C-terminal subdomains, while the more flexible N-terminal subdomains have more comparable average internal correlation times.     

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.     

Improvement of hydrogen bond geometry in protein NMR structures by residual dipolar couplings – an assessment of the interrelation of NMR restraints

We have examined how the hydrogen bond geometry in three different proteins is affected when structural restraints based on measurements of residual dipolar couplings are included in the structure calculations. The study shows, that including restraints based solely on (1)H(N)-(15)N residual dipolar couplings has pronounced impact on the backbone rmsd and Ramachandran plot but does not improve the hydrogen bond geometry. In the case of chymotrypsin inhibitor 2 the addition of (13)CO-(13)C(alpha) and (15)N-(13)CO one bond dipolar couplings as restraints in the structure calculations improved the hydrogen bond geometry to a quality comparable to that obtained in the 1.8 A resolution X-ray structure of this protein. A systematic restraint study was performed, in which four types of restraints, residual dipolar couplings, hydrogen bonds, TALOS angles and NOEs, were allowed in two states. This study revealed the importance of using several types of residual dipolar couplings to get good hydrogen bond geometry. The study also showed that using a small set of NOEs derived only from the amide protons, together with a full set of residual dipolar couplings resulted in structures of very high quality. When reducing the NOE set, it is mainly the side-chain to side-chain NOEs that are removed. Despite of this the effect on the side-chain packing is very small when a reduced NOE set is used, which implies that the over all fold of a protein structure is mainly determined by correct folding of the backbone.