Deuterium spin probes of backbone order in proteins: 2H NMR relaxation study of deuterated carbon alpha sites

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

Stereoselective Biocatalytic α-Deuteration of L-Amino Acids by a Pyridoxal 5'-Phosphate-Dependent Mannich Cyclase

α-Deuterated amino acids are valuable building blocks for developing deuterated drugs, and are important tools for studying biological systems. Biocatalytic deuteration represents an attractive strategy to directly access enantiopure α-deuterated amino acids. Here, we show that a PLP-dependent Mannich cyclase, LolT, involved in the biosynthesis of loline alkaloids, is capable of deuterating a diverse range of L-amino acids, including basic and acidic, nonpolar and polar, aliphatic and aromatic amino acids. Furthermore, complete deuteration of many amino acids can be achieved within minutes with exquisite control on the site- and stereoselectivity. During the course of this investigation, we also unexpectedly discovered that LolT exhibits β-elimination activity with L-cystine and O-acetyl-L-serine, confirming our previous hypothesis based on structural and phylogenetic analysis that LolT, a Cα-C bond forming enzyme, is evolved from a primordial Cβ-S lyase family. Overall, our study demonstrates that LolT is an extremely versatile biocatalyst, and can be used for not only heterocyclic quaternary amino acid biosynthesis, but also biocatalytic amino acid deuteration.

Site-specific protein backbone deuterium 2Hα quadrupolar patterns by proton-detected quadruple-resonance 3D 2HαcαNH MAS NMR spectroscopy

A novel deuterium-excited and proton-detected quadruple-resonance three-dimensional (3D) ²H_αc_αNH MAS nuclear magnetic resonance (NMR) method is presented to obtain site-specific ²H_α deuterium quadrupolar couplings from protein backbone, as an extension to the 2D version of the experiment reported earlier. Proton-detection results in high sensitivity compared to the heteronuclei detection methods. Utilizing four independent radiofrequency (RF) channels (quadruple-resonance), we managed to excite the ²H_α, then transfer deuterium polarization to its attached C_α, followed by polarization transfers to the neighboring backbone nitrogen and then to the amide proton for detection. This experiment results in an easy to interpret HSQC-like 2D ¹H-¹⁵N fingerprint NMR spectrum, which contains site-specific deuterium quadrupolar patterns in the indirect third dimension. Provided that four-channel NMR probe technology is available, the setup of the ²H_αc_αNH experiment is relatively straightforward, by using low power deuterium excitation and polarization transfer schemes we have been developing. To our knowledge, this is the first demonstration of a quadruple-resonance MAS NMR experiment to link ²H_α quadrupolar couplings to proton-detection, extending our previous triple-resonance demonstrations. Distortion-free excitation and polarization transfer of ∼160-170 kHz ²H_α quadrupolar coupling were presented by using a deuterium RF strength of ∼20 kHz. From these ²H_α patterns, an average backbone order parameter of S = 0.92 was determined on a deuterated SH3 sample, with an average η = 0.22. These indicate that SH3 backbone represents sizable dynamics in the microsecond timescale where the ²H_α lineshape is sensitive. Moreover, site-specific ²H_α T₁ relaxation times were obtained for a proof of concept. This 3D ²H_αc_αNH NMR experiment has the potential to determine structure and dynamics of perdeuterated proteins by utilizing deuterium as a novel reporter.

Capturing the Polarization Response of Solvated Proteins under Constant Electric Fields in Molecular Dynamics Simulations

We capture and compare the polarization response of a solvated globular protein ubiquitin to static electric (E-fields) using atomistic molecular dynamics simulations. We collectively follow E-field induced changes, electrical and structural, occurring across multiple trajectories using the magnitude of the protein dipole vector (P_p ). E-fields antiparallel to P_p induce faster structural changes and more facile protein unfolding relative to parallel fields of the same strength. While weak E-fields (0.1-0.5 V/nm) do not unfold ubiquitin and produce a reversible polarization, strong E-fields (1-2 V/nm) unfold the protein through a pathway wherein the helix:β-strand interactions rupture before those for the β1-β5 clamp. Independent of E-field direction, high E-field induced structural changes are also reversible if the field is switched off before P_p exceeds 2 times its equilibrium value. We critically examine the dependence of water properties, protein rotational diffusion and E-field induced protein unfolding pathways on the thermostat/barostat parameters used in our simulations.

Probing Methyl Group Dynamics in Proteins by NMR Cross-Correlated Dipolar Relaxation and Molecular Dynamics Simulations

Nuclear magnetic resonance (NMR) spin relaxation is the most informative approach to experimentally probe the internal dynamics of proteins on the picosecond to nanosecond time scale. At the same time, molecular dynamics (MD) simulations of biological macromolecules are steadily improving through better physical models, enhanced sampling methods, and increased computational power, and they provide exquisite information about flexibility and its role in protein stability and molecular interactions. Many examples have shown that MD is now adept in probing protein backbone motion, but improvements are still required toward a quantitative description of the dynamics of side chains, for example, probed by the dynamics of methyl groups. Thus far, the comparison of computation with experiment for side chain dynamics has primarily focused on the relaxation of ¹³C and ²H nuclei induced by autocorrelated variation of spin interactions. However, the cross-correlation of ¹³C-¹H dipolar interactions in methyl groups offers an attractive alternative. Here, we establish a computational framework to extract cross-correlation relaxation parameters of methyl groups in proteins from all-atom MD simulations. To demonstrate the utility of the approach, cross-correlation relaxation rates of ubiquitin are computed from MD simulations performed with the AMBER99SB*-ILDN and CHARMM36 force fields. Simulation results were found to agree well with those obtained by experiment. Moreover, the data obtained with the two force fields are highly consistent.

Stereoselective α-Deuteration of Serine, Cysteine, Selenocysteine, and 2,3-Diaminopropanoic Acid Derivatives

Efficient methodologies for synthesizing enantiopure α-deuterated derivatives of serine, cysteine, selenocysteine, and 2,3-diaminopropanoic acid have been developed. H/D exchange was achieved by deprotonation of a chiral bicyclic serine equivalent followed by selective deuteration. Additionally, diastereoselective additions of thiols, selenols, and amines to a chiral bicyclic dehydroalanine in deuterated alcohols allowed site-selective deuteration at the Cα atom of cysteine, selenocysteine, and 2,3-diaminopropanoic acid derivatives. A deuterated analogue of carbocysteine, a drug for the treatment of bronchiectasis, was synthesized.

Beyond structure: Deciphering site-specific dynamics in proteins from double histidine-based EPR measurements

Site-specific dynamics in proteins are at the heart of protein function. While electron paramagnetic resonance (EPR) has potential to measure dynamics in large protein complexes, the reliance on flexible nitroxide labels is limitating especially for the accurate measurement of site-specific β-sheet dynamics. Here, we employed EPR spectroscopy to measure site-specific dynamics across the surface of a protein, GB1. Through the use of the double Histidine (dHis) motif, which enables labeling with a Cu(II) - nitrilotriacetic acid (NTA) complex, dynamics information was obtained for both α-helical and β-sheet sites. Spectral simulations of the resulting CW-EPR report unique site-specific fluctuations across the surface of GB1. Additionally, we performed molecular dynamics (MD) simulations to complement the EPR data. The dynamics observed from MD agree with the EPR results. Furthermore, we observe small changes in gǁ values for different sites, which may be due to small differences in coordination geometry and/or local electrostatics of the site. Taken together, this work expands the utility of Cu(II)NTA-based EPR measurements to probe information beyond distance constraints.

Site-Selective Deuteration of Amino Acids through Dual-Protein Catalysis

Deuterated amino acids have been recognized for their utility in drug development, for facilitating nuclear magnetic resonance (NMR) analysis, and as probes for enzyme mechanism. Small molecule-based methods for the site-selective synthesis of deuterated amino acids typically involve de novo synthesis of the compound from deuterated precursors. In comparison, enzymatic methods for introducing deuterium offer improved efficiency, operating directly on free amino acids to achieve hydrogen-deuterium (H/D) exchange. However, site selectivity remains a significant challenge for enzyme-mediated deuteration, limiting access to desirable deuteration motifs. Here, we use enzyme-catalyzed deuteration, combined with steady-state kinetic analysis and ultraviolet (UV)-vis spectroscopy to probe the mechanism of a two-protein system responsible for the biosynthesis of l-allo-Ile. We show that an aminotransferase (DsaD) can pair with a small partner protein (DsaE) to catalyze Cα and Cβ H/D exchange of amino acids, while reactions without DsaE lead exclusively to Cα-deuteration. With conditions for improved catalysis, we evaluate the substrate scope for Cα/Cβ-deuteration and demonstrate the utility of this system for preparative-scale, selective labeling of amino acids.

A CEST NMR experiment to obtain glycine 1Hα chemical shifts in 'invisible' minor states of proteins

Chemical exchange saturation transfer (CEST) experiments are routinely used to study protein conformational exchange between a 'visible' major state and 'invisible' minor states because they can detect minor states with lifetimes varying from ~ 3 to ~ 100 ms populated to just ~ 0.5%. Consequently several ¹H, ¹⁵N and ¹³C CEST experiments have been developed to study exchange and obtain minor state chemical shifts at almost all backbone and sidechain sites in proteins. Conspicuously missing from this extensive set of CEST experiments is a ¹H CEST experiment to study exchange at glycine (Gly) ¹H^α sites as the existing ¹H CEST experiments that have been designed to study dynamics in amide ¹H-¹⁵N spin systems and methyl ¹³CH₃ groups with three equivalent protons while suppressing ¹H-¹H NOE induced dips are not suitable for studying exchange in methylene ¹³CH₂ groups with inequivalent protons. Here a Gly ¹H^α CEST experiment to obtain the minor state Gly ¹H^α chemical shifts is presented. The utility of this experiment is demonstrated on the L99A cavity mutant of T4 Lysozyme (T4L L99A) that undergoes conformational exchange between two compact conformers. The CEST derived minor state Gly ¹H^α chemical shifts of T4L L99A are in agreement with those obtained previously using CPMG techniques. The experimental strategy presented here can also be used to obtain methylene proton minor state chemical shifts from protein sidechain and nucleic acid backbone sites.

Biocatalytic, Stereoselective Deuteration of α-Amino Acids and Methyl Esters

α-²H amino acids are valuable precursors toward labeled pharmaceutical agents and tools for studying biological systems; however, these molecules are costly to purchase and challenging to synthesize in a site- and stereoselective manner. Here, we show that an α-oxo-amine synthase that evolved for saxitoxin biosynthesis, SxtA AONS, is capable of producing a range of α-²H amino acids and esters site- and stereoselectively using D₂O as the deuterium source. Additionally, we demonstrate the utility of this operationally simple reaction on preparative scale in the stereoselective chemoenzymatic synthesis of a deuterated analog of safinamide, a drug used to treat Parkinson's disease.

Structural and dynamic origins of ESR lineshapes in spin-labeled GB1 domain: the insights from spin dynamics simulations based on long MD trajectories

Site-directed spin labeling (SDSL) ESR is a valuable tool to probe protein systems that are not amenable to characterization by x-ray crystallography, NMR or EM. While general principles that govern the shape of SDSL ESR spectra are known, its precise relationship with protein structure and dynamics is still not fully understood. To address this problem, we designed seven variants of GB1 domain bearing R1 spin label and recorded the corresponding MD trajectories (combined length 180 μs). The MD data were subsequently used to calculate time evolution of the relevant spin density matrix and thus predict the ESR spectra. The simulated spectra proved to be in good agreement with the experiment. Further analysis confirmed that the spectral shape primarily reflects the degree of steric confinement of the R1 tag and, for the well-folded protein such as GB1, offers little information on local backbone dynamics. The rotameric preferences of R1 side chain are determined by the type of the secondary structure at the attachment site. The rotameric jumps involving dihedral angles χ1 and χ2 are sufficiently fast to directly influence the ESR lineshapes. However, the jumps involving multiple dihedral angles tend to occur in (anti)correlated manner, causing smaller-than-expected movements of the R1 proxyl ring. Of interest, ESR spectra of GB1 domain with solvent-exposed spin label can be accurately reproduced by means of Redfield theory. In particular, the asymmetric character of the spectra is attributable to Redfield-type cross-correlations. We envisage that the current MD-based, experimentally validated approach should lead to a more definitive, accurate picture of SDSL ESR experiments.

Aromatic Ring Dynamics, Thermal Activation, and Transient Conformations of a 468 kDa Enzyme by Specific 1H-13C Labeling and Fast Magic-Angle Spinning NMR

Aromatic residues are located at structurally important sites of many proteins. Probing their interactions and dynamics can provide important functional insight but is challenging in large proteins. Here, we introduce approaches to characterize the dynamics of phenylalanine residues using ¹H-detected fast magic-angle spinning (MAS) NMR combined with a tailored isotope-labeling scheme. Our approach yields isolated two-spin systems that are ideally suited for artifact-free dynamics measurements, and allows probing motions effectively without molecular weight limitations. The application to the TET2 enzyme assembly of ∼0.5 MDa size, the currently largest protein assigned by MAS NMR, provides insights into motions occurring on a wide range of time scales (picoseconds to milliseconds). We quantitatively probe ring-flip motions and show the temperature dependence by MAS NMR measurements down to 100 K. Interestingly, favorable line widths are observed down to 100 K, with potential implications for DNP NMR. Furthermore, we report the first ¹³C R_1ρ MAS NMR relaxation-dispersion measurements and detect structural excursions occurring on a microsecond time scale in the entry pore to the catalytic chamber and at a trimer interface that was proposed as the exit pore. We show that the labeling scheme with deuteration at ca. 50 kHz MAS provides superior resolution compared to 100 kHz MAS experiments with protonated, uniformly ¹³C-labeled samples.

Determination of Protein ps-ns Motions by High-Resolution Relaxometry

Many of the functions of biomacromolecules can be rationalized by the characterization of their conformational energy landscapes: the structures of the dominant states, transitions between states and motions within states. Nuclear magnetic resonance (NMR) spectroscopy is the technique of choice to study internal motions in proteins. The determination of motions on picosecond to nanosecond timescales requires the measurement of nuclear spin relaxation rates at multiple magnetic fields. High sensitivity and resolution are obtained only at high magnetic fields, so that, until recently, site-specific relaxation rates in biomolecules were only measured over a narrow range of high magnetic fields. This limitation was particularly striking for the quantification of motions on nanosecond timescales, close to the correlation time for overall rotational diffusion. High-resolution relaxometry is an emerging technique to investigate picosecond-nanosecond motions of proteins. This approach uses a high-field NMR spectrometer equipped with a sample shuttle device, which allows for the measurement of the relaxation rate constants at low magnetic fields, while preserving the sensitivity and resolution of a high-field NMR spectrometer. The combined analysis of high-resolution relaxometry and standard high-field relaxation data provides a more accurate description of the dynamics of proteins, in particular in the nanosecond range. The purpose of this chapter is to describe how to perform high-resolution relaxometry experiments and how to analyze the rates measured with this technique.

Enhanced spectral density mapping through combined multiple-field deuterium 13CH2D methyl spin relaxation NMR spectroscopy

Quadrupolar relaxation of ²H (D) nuclear spins is a powerful probe of conformational dynamics in biological macromolecules. Deuterium relaxation rate constants are determined by the spectral density function for reorientation of the C-D bond vector at zero, single-quantum, and double-quantum ²H frequencies. In the present work, ²H relaxation rate constants were measured for an E. coli ribonuclease H [U-²H, ¹⁵N] ILV-[¹³CH₂D] sample using 400, 500, 800, and 900 MHz NMR spectrometers and analyzed by three approaches to determine spectral density values. First, data recorded at each static magnetic field were analyzed independently. Second, data recorded at 400 and 800 MHz were analyzed jointly and data recorded at other fields were analyzed independently. Third, data recorded at 400 and 500 MHz were interpolated to 450 MHz, and the resulting two pairs of data, corresponding to 400 MHz/800 MHz and 450 MHz/900 MHz, were analyzed jointly. The second and third approaches rely on the identity between the double quantum frequency at the lower field and the single quantum frequency at the higher field. Spectral density values for 32 of the 48 resolvable ILV methyl resonances were fit by the Lipari-Szabo model-free formalism and used to validate the three methods. The three spectral density mapping methods performed equally well in cross validation with data recorded at 700 MHz. However, the third method yielded approximately 10-15% more precise estimates of model-free parameters and consequently provides a general strategy for analysis of ²H spin relaxation data in biological macromolecules.

Site-Specific Internal Motions in GB1 Protein Microcrystals Revealed by 3D ²H-¹³C-¹³C Solid-State NMR Spectroscopy

(2)H quadrupolar line shapes deliver rich information about protein dynamics. A newly designed 3D (2)H-(13)C-(13)C solid-state NMR magic angle spinning (MAS) experiment is presented and demonstrated on the microcrystalline β1 immunoglobulin binding domain of protein G (GB1). The implementation of (2)H-(13)C adiabatic rotor-echo-short-pulse-irradiation cross-polarization (RESPIRATION CP) ensures the accuracy of the extracted line shapes and provides enhanced sensitivity relative to conventional CP methods. The 3D (2)H-(13)C-(13)C spectrum reveals (2)H line shapes for 140 resolved aliphatic deuterium sites. Motional-averaged (2)H quadrupolar parameters obtained from the line-shape fitting identify side-chain motions. Restricted side-chain dynamics are observed for a number of polar residues including K13, D22, E27, K31, D36, N37, D46, D47, K50, and E56, which we attribute to the effects of salt bridges and hydrogen bonds. In contrast, we observe significantly enhanced side-chain flexibility for Q2, K4, K10, E15, E19, N35, N40, and E42, due to solvent exposure and low packing density. T11, T16, and T17 side chains exhibit motions with larger amplitudes than other Thr residues due to solvent interactions. The side chains of L5, V54, and V29 are highly rigid because they are packed in the core of the protein. High correlations were demonstrated between GB1 side-chain dynamics and its biological function. Large-amplitude side-chain motions are observed for regions contacting and interacting with immunoglobulin G (IgG). In contrast, rigid side chains are primarily found for residues in the structural core of the protein that are absent from protein binding and interactions.

NMR studies of dynamic biomolecular conformational ensembles

Multidimensional heteronuclear NMR approaches can provide nearly complete sequential signal assignments of isotopically enriched biomolecules. The availability of assignments together with measurements of spin relaxation rates, residual spin interactions, J-couplings and chemical shifts provides information at atomic resolution about internal dynamics on timescales ranging from ps to ms, both in solution and in the solid state. However, due to the complexity of biomolecules, it is not possible to extract a unique atomic-resolution description of biomolecular motions even from extensive NMR data when many conformations are sampled on multiple timescales. For this reason, powerful computational approaches are increasingly applied to large NMR data sets to elucidate conformational ensembles sampled by biomolecules. In the past decade, considerable attention has been directed at an important class of biomolecules that function by binding to a wide variety of target molecules. Questions of current interest are: "Does the free biomolecule sample a conformational ensemble that encompasses the conformations found when it binds to various targets; and if so, on what time scale is the ensemble sampled?" This article reviews recent efforts to answer these questions, with a focus on comparing ensembles obtained for the same biomolecules by different investigators. A detailed comparison of results obtained is provided for three biomolecules: ubiquitin, calmodulin and the HIV-1 trans-activation response RNA.

Indirect use of deuterium in solution NMR studies of protein structure and hydrogen bonding

A description of the utility of deuteration in protein NMR is provided with an emphasis on quantitative evaluation of the effects of deuteration on a number of NMR parameters of proteins: (1) chemical shifts, (2) scalar coupling constants, (3) relaxation properties (R1 and R2 rates) of nuclei directly attached to one or more deuterons as well as protons of methyl groups in a highly deuterated environment, (4) scalar relaxation of 15N and 13C nuclei in 15N-D and 13C-D spin systems as a measure of hydrogen bonding strength, and (5) NOE-based applications of deuteration in NMR studies of protein structure. The discussion is restricted to the 'indirect' use of deuterium in the sense that the description of NMR parameters and properties of the nuclei affected by nearby deuterons (15N, 13C, 1H) is provided rather than those of deuterium itself.

Deriving quantitative dynamics information for proteins and RNAs using ROTDIF with a graphical user interface

To facilitate rigorous analysis of molecular motions in proteins, DNA, and RNA, we present a new version of ROTDIF, a program for determining the overall rotational diffusion tensor from single- or multiple-field nuclear magnetic resonance relaxation data. We introduce four major features that expand the program's versatility and usability. The first feature is the ability to analyze, separately or together, (13)C and/or (15)N relaxation data collected at a single or multiple fields. A significant improvement in the accuracy compared to direct analysis of R2/R1 ratios, especially critical for analysis of (13)C relaxation data, is achieved by subtracting high-frequency contributions to relaxation rates. The second new feature is an improved method for computing the rotational diffusion tensor in the presence of biased errors, such as large conformational exchange contributions, that significantly enhances the accuracy of the computation. The third new feature is the integration of the domain alignment and docking module for relaxation-based structure determination of multi-domain systems. Finally, to improve accessibility to all the program features, we introduced a graphical user interface that simplifies and speeds up the analysis of the data. Written in Java, the new ROTDIF can run on virtually any computer platform. In addition, the new ROTDIF achieves an order of magnitude speedup over the previous version by implementing a more efficient deterministic minimization algorithm. We not only demonstrate the improvement in accuracy and speed of the new algorithm for synthetic and experimental (13)C and (15)N relaxation data for several proteins and nucleic acids, but also show that careful analysis required especially for characterizing RNA dynamics allowed us to uncover subtle conformational changes in RNA as a function of temperature that were opaque to previous analysis.

Four-bond deuterium isotope effects on the chemical shifts of amide nitrogens in proteins

An approach towards precision NMR measurements of four-bond deuterium isotope effects on the chemical shifts of backbone amide nitrogen nuclei in proteins is described. Three types of four-bond (15) N deuterium isotope effects are distinguished depending on the site of proton-to-deuterium substitution: (4)ΔN(N(i-1)D), (4)ΔN(N(i+1)D) and (4)ΔN(Cβ,(i-1)D). All the three types of isotope shifts are quantified in the (partially) deuterated protein ubiquitin. The (4)ΔN(N(i+1)D) and (4)ΔN(C(β,i-1)D) effects are by far the largest in magnitude and vary between 16 and 75 ppb and -18 and 46 ppb, respectively. A semi-quantitative correlation between experimental (4)ΔN(N(i+1)D) and (4)ΔN(C(β,i-1)D) values and the distances between nitrogen nuclei and the sites of (1)H-to-D substitution is noted. The largest isotope shifts in both cases correspond to the shortest inter-nuclear distances.

Accurate measurements of the effects of deuteration at backbone amide positions on the chemical shifts of ¹⁵N, ¹³Cα, ¹³Cβ, ¹³CO and ¹Hα nuclei in proteins

An approach towards accurate NMR measurements of deuterium isotope effects on the chemical shifts of all backbone nuclei in proteins ((15)N, (13)Cα, (13)CO, (1)Hα) and (13)Cβ nuclei arising from (1)H-to-D substitutions at amide nitrogen positions is described. Isolation of molecular species with a defined protonation/deuteration pattern at successive backbone nitrogen positions in the polypeptide chain allows quantifying all deuterium isotope shifts of these nuclei from the first to the fourth order. Some of the deuterium isotope shifts measured in the proteins ubiquitin and GB1 can be interpreted in terms of backbone geometry via empirical relationships describing their dependence on (φ; ψ) backbone dihedral angles. Because of their relatively large variability and notable dependence on the protein secondary structure, the two- and three-bond (13)Cα isotope shifts, (2)ΔCα(NiD) and (3)ΔCα(Ni+1D), and three-bond (13)Cβ isotope shifts, (3)ΔCβ(NiD), are useful reporters of the local geometry of the protein backbone.

Estimating side-chain order in methyl-protonated, perdeuterated proteins via multiple-quantum relaxation violated coherence transfer NMR spectroscopy

Relaxation violated coherence transfer NMR spectroscopy (Tugarinov et al. in J Am Chem Soc 129:1743-1750, 2007) is an established experimental tool for quantitative estimation of the amplitudes of side-chain motions in methyl-protonated, highly deuterated proteins. Relaxation violated coherence transfer experiments monitor the buildup of methyl proton multiple-quantum coherences that can be created in magnetically equivalent spin-systems as long as their transverse magnetization components relax with substantially different rates. The rate of this build-up is a reporter of the methyl-bearing side-chain mobility. Although the build-up of multiple-quantum (1)H coherences is monitored in these experiments, the decay of the methyl signal during relaxation delays occurs when methyl proton magnetization is in a single-quantum state. We describe a relaxation violated coherence transfer approach where the relaxation of multiple-quantum (1)H-(13)C methyl coherences during the relaxation delay period is quantified. The NMR experiment and the associated fitting procedure that models the time-dependence of the signal build-up, are applicable to the characterization of side-chain order in [(13)CH(3)]-methyl-labeled, highly deuterated protein systems up to ~100 kDa in molecular weight. The feasibility of extracting reliable measures of side-chain order is experimentally verified on methyl-protonated, perdeuterated samples of an 8.5-kDa ubiquitin at 10°C and an 82-kDa Malate Synthase G at 37°C.

Probing side-chain dynamics in proteins by the measurement of nine deuterium relaxation rates per methyl group

We demonstrate the feasibility of the measurement of up to nine deuterium spin relaxation rates in 13CHD2 and 13CH2D methyl isotopomers of small proteins. In addition to five measurable 2H relaxation rates in a 13CH2D methyl group (Millet, O.; Muhandiram, D. R.; Skrynnikov, N. R.; Kay, L. E. J. Am. Chem. Soc. 2002, 124, 6439-48), the measurement of additional four rates of (nearly) single-exponentially decaying magnetization terms in methyl groups of the 13CHD2 variety is reported. Consistency relationships between 2H spin relaxation rates measured in the two different types of methyl groups are derived and verified experimentally for a subset of methyl-containing side chains in the protein ubiquitin. A detailed comparison of methyl-bearing side-chain dynamics parameters obtained from relaxation measurements in 13CH2D and 13CHD2 methyls of ubiquitin at 10, 27, and 40 °C reveals that transverse 2H relaxation rates in 13CHD2 groups are reliable and accurate reporters of the amplitudes of methyl 3-fold axis motions (S(axis)2) for protein molecules with global molecular tumbling times τ(C) >~9 ns. For smaller molecules, simple correction of transverse 2H relaxation rates in 13CHD2 groups is sufficient for the derivation of robust measures of order. Residue-specific distributions of S(axis)2 are consistent with atomic-detail molecular dynamics (MD) results. Both 13CHD2- and 13CH2D-derived S(axis)2 values are in good overall agreement with those obtained from 1 μs MD simulations at all the three temperatures, although some differences in the site-specific temperature dependence between MD- and 2H-relaxation-derived S(axis)2 values are observed.

Toward a predictive understanding of slow methyl group dynamics in proteins

The development of the most recent generation of molecular mechanics force fields promises an increasingly predictive understanding of the protein dynamics-function relationship. Based on extensive validation against various types of experimental data, the AMBER force field ff99SB was benchmarked in recent years as a favorable force field for protein simulations. Recent improvements of the side chain and backbone potentials, made by different groups, led to the ff99SB-ILDN and ff99SBnmr1 force fields, respectively. The combination of these potentials into a unified force field, termed ff99SBnmr1-ILDN, was used in this study to perform a microsecond time scale molecular dynamics simulation of free ubiquitin in explicit solvent for validation against an extensive set of experimental NMR methyl group residual dipolar couplings. Our results show a high level of consistency between the experimental data and the values predicted from the molecular dynamics trajectory reflecting a systematically improved performance of ff99SBnmr1-ILDN over the original ff99SB force field. Moreover, the unconstrained ff99SBnmr1-ILDN MD ensemble achieves a similar level of agreement as the recently introduced EROS ensemble, which was constructed based on a large body of NMR data as constraints, including the methyl residual dipolar couplings. This suggests that ff99SBnmr1-ILDN provides a high-quality representation of the motions of methyl-bearing protein side chains, which are sensitive probes of protein-protein and protein-ligand interactions.

An optimized relaxation-based coherence transfer NMR experiment for the measurement of side-chain order in methyl-protonated, highly deuterated proteins

Relaxation violated coherence transfer NMR spectroscopy has emerged as a powerful experimental tool for the quantitative measurement of amplitudes of motion of methyl containing side-chains. Typically, the experiments, performed on proteins that are highly deuterated and methyl-protonated, monitor the build-up of methyl (1)H double-quantum magnetization. Because all three protons in a methyl group are degenerate, such coherences can only result from differential relaxation of transverse magnetization components, which in turn reflect the extent and time-scale of motion of the methyl probe [Tugarinov, V., Sprangers, R.; Kay, L.E. J. Am. Chem. Soc. 2007, 129, 1743-1750]. We show here that a 50% gain in the sensitivity of the experiment can be achieved through selection of (1)H triple-quantum coherence, thereby significantly increasing the utility of the approach. A theoretical treatment rationalizes the sensitivity gain that is subsequently verified through experiment. The utility of the methodology is demonstrated on a number of proteins, including the 360 kDa α(7)α(7) "half-proteasome".

Iterative Optimization of Molecular Mechanics Force Fields from NMR Data of Full-Length Proteins

High quality molecular mechanics force fields of proteins are key for the quantitative interpretation of experimental data and the predictive understanding of protein function based on computer simulations. A strategy is presented for the optimization of protein force fields based on full-length proteins in their native environment that is guided by experimental NMR chemical shifts and residual dipolar couplings (RDCs). An energy-based reweighting approach is applied to a long molecular dynamics trajectory, performed with a parent force field, to efficiently screen a large number of trial force fields. The force field that yields the best agreement with the experimental data is then used as the new parent force field, and the procedure is repeated until no further improvement is obtained. This method is demonstrated for the optimization of the backbone φ,ψ dihedral angle potential of the Amber ff99SB force field using six trial proteins and another 17 proteins for cross-validation using (13)C chemical shifts with and without backbone RDCs. The φ,ψ dihedral angle potential is systematically improved by the inclusion of correlation effects through the addition of up to 24 bivariate Gaussian functions of variable height, width, and tilt angle. The resulting force fields, termed ff99SB_φψ(g24;CS) and ff99SB_φψ(g8;CS,RDC), perform significantly better than their parent force field in terms of both NMR data reproduction and Cartesian coordinate root-mean-square deviations between the MD trajectories and the X-ray crystal structures. The strategy introduced here represents a powerful addition to force field optimization approaches by overcoming shortcomings of methods that are solely based on quantum-chemical calculations of small molecules and protein fragments in the gas phase.

Variation in quadrupole couplings of alpha deuterons in ubiquitin suggests the presence of C(alpha)-H(alpha)...O=C hydrogen bonds

Nuclear quadrupolar couplings are sensitive probes of hydrogen bonding. Experimental quadrupolar coupling constants of alpha deuterons (D(alpha) QCC) are reported for the residues of human ubiquitin that do not experience large-amplitude internal dynamics on the pico- to nanosecond time scale. Two different methods for D(alpha) QCC estimation are employed: (i) direct estimation of D(alpha) QCC values from R(1) and R(2) (2)H D(alpha) rates using the dynamics parameters (S(C(alpha)-H(alpha))(2)) derived from 1 micros molecular dynamics simulations as well as from (13)C(alpha) relaxation measurements and (ii) indirect measurements via scalar relaxation of the second kind that affects (13)C(alpha) relaxation rates in (13)C(alpha)-D(alpha) spin systems. A relatively large variability of D(alpha) QCC values is produced by both methods. The average value of 170.6 +/- 3 kHz is derived from the combined data set, with D(alpha) QCC values ranging from 159.2 to 177.2 kHz. The set of lowest quadrupolar couplings in all data sets corresponds to the residues that are likely to form weak C(alpha)-H(alpha)...O=C hydrogen bonds as predicted from the analysis of short H(alpha)...O distances in three-dimensional structures of ubiquitin. These D(alpha) nuclei show up to 10 kHz reduction in their QCC values, which is in agreement with earlier solid-state NMR measurements in alpha deuterons of glycine. A statistically significant correlation is observed between the QCC values of alpha-deuterons and the inverse cube of C(alpha)-H(alpha)...O=C distances in ubiquitin.

Estimating quadrupole couplings of amide deuterons in proteins from direct measurements of 2H spin relaxation rates

The measurements of longitudinal and transverse 2H spin relaxation rates of backbone amide deuterons (D(N)) in the [U-13C,15N]-labeled protein ubiquitin show that the utility of amide deuterons as probes of backbone order in proteins is compromised by substantial variability of D(N) quadrupolar coupling constants (QCC) from one amide site to another. However, using the dynamics parameters of 15N-2H bond vectors evaluated from 15N relaxation data, site-specific QCC values can be estimated directly from D(N)R1 and R2 rates providing useful information on hydrogen bonding in proteins. In agreement with previous indirect scalar relaxation-based measurements, the D(N) QCC values estimated directly from R1 and R2 2H relaxation rates correlate with the inverse cube of the X-ray structure-derived hydrogen bond distances in ubiquitin: QCC=(232+/-2.3)+(118+/-17) summation operator(i)(cosalpha)r(i)(-3) where r is the inter-nuclear hydrogen bond distance in ångströms, and alpha is the N-D....O(i) angle.