- 1Department of Physics, University of Warwick, Coventry, United Kingdom
- 2Department of Chemistry, University of Warwick, Coventry, United Kingdom
Order parameters are a useful tool for quantifying amplitudes of molecular motions. Here we measure dipolar order parameters by recoupling heteronuclear dipole-dipole couplings under fast spinning. We apply symmetry based recoupling methods to samples spinning under magic angle at 60 kHz by employing a variable flip angle compound inversion pulse. We validate the methods by measuring site-specific 15N-1H order parameters of a microcrystalline protein over a small temperature range and the same protein in a large, precipitated complex with antibody. The measurements of the order parameters in the complex are consistent with the observed protein undergoing overall motion within the assembly.
Introduction
The astounding number of structures found in the protein databank speaks to the usefulness of structural data to provide insights into the structure-function relationship in biology and biochemistry (Berman et al., 2000; Burley et al., 2021). With advent of powerful computational structure prediction approaches such as AlphaFold (Jumper et al., 2021) and RoseTTAFold (Baek et al., 2021) there is an almost exponential increase of systems for which a reasonable quality structure or models become available. However, a structure is a snapshot that does not necessarily capture the choreography of the protein it needs to execute in order to perform its function (Koshland, 1958; Frauenfelder et al., 1991). The motion of a protein is often intrinsic to its activity. Understanding the dynamics and how structure changes in time is sometimes nearly as important as knowing a single, even high-resolution, snapshot. An ultimate example of this idea are intrinsically disordered proteins (IDP) and intrinsically disordered regions (IDR), which are involved in controlling countless processes in eukaryotic but also prokaryotic organisms, e.g. biosynthetic steps in production of bioactive natural products (Jenner et al., 2018; Kosol et al., 2019; Fage et al., 2021).
Nuclear magnetic resonance (NMR), can be used to find molecular motions under near physiological conditions at atomic resolution over a several orders of magnitude of the time scale, from as fast as picoseconds, to as slow as months, but will only report on local conditions and typically over small distance scales (Palmer, 2004; Kovermann et al., 2016; Sekhar and Kay, 2019; Alderson and Kay, 2021). However, to access such a vast range of dynamics a battery of different tools reporting on different parameters of motion in different regimes is required. For example, NMR relaxation is sensitive to both amplitudes and time scales of motions typically in the picoseconds-nanoseconds range in solution and picoseconds-milliseconds range in the solid state (Lewandowski, 2013), which provides some unique opportunities for characterizing protein motions (Castellani et al., 2002; Chevelkov et al., 2003; Chevelkov et al., 2006; Lewandowski et al., 2011; Asami and Reif, 2012; Lamley et al., 2014; Lamley et al., 2015a; Lamley et al., 2015b; Sternberg et al., 2018; Öster et al., 2019). However, the extended range of time scales of motions, which influence relaxation in the solid state, comes also at a price: reliable quantification of motional amplitudes with relaxation rates alone is challenging and sometimes impossible. Often reliable quantification of dynamics using relaxation rates requires them to be combined with measurements of order parameters (typically dipolar order parameters), which constrain the overall amplitude of motions (Schanda and Ernst, 2016).
Order parameters can be obtained by recoupling specific terms of the NMR Hamiltonian in a separated local field (SLF) experiment (Hester et al., 1976). SLF techniques use pulses to create a Hamiltonian where one term of the full Hamiltonian is recoupled, and the other terms are averaged to zero. Examples of such methods include C7 (Hohwy et al., 1998), RFDR (Bennett et al., 1992), REDOR (Gullion and Schaefer, 1989a; Gullion and Schaefer, 1989b), TMREV (Hohwy et al., 2000; Franks et al., 2005; Franks et al., 2006a), and many others. The recoupling portion of SLF experiments have been summarized into a uniform theory using symmetry principles (Levitt, 2002). The development has mainly been focused on slowly spinning samples, as fast rotation was not available at the time. However, fast magic angle spinning (MAS) NMR experiments, introduced since the main formulation of this theory, have improved the process of assignment and structure calculation of large proteins and complexes that were very difficult to solve using solid-state NMR otherwise (Zhou et al., 2007a; Barbet-Massin et al., 2010; Knight et al., 2012; Barbet-Massin et al., 2014). For example, the membrane protein OMPG had been actively studied for almost 15 years using carbon detection but the assignment and structure were finally solved with the use of fast spinning and proton detection (Hiller et al., 2005; Retel et al., 2017). Thus, it is desirable to extend symmetry methods to this attractive new regime. Unfortunately, symmetry sequences require applied fields that scale linearly with the spinning rate, and thus the applied field requirements cannot be usually achieved under fast spinning conditions. The application of symmetry principles to heteronuclear dipole-dipole recoupling was previously demonstrated under 40 kHz spinning (Hou et al., 2011), but did not engender optimism for application at higher spinning rates. However, spinning rates of ∼60 kHz are routine at the time of writing of this manuscript, with 100 kHz spinning becoming more common, and current cutting-edge probes reach rates on the order of 150 kHz (Penzel et al., 2019; Schledorn et al., 2020) and even 200 kHz. Consequently, symmetry methods have not been applied extensively to fast spinning samples (Brinkmann and Levitt, 2001; Levitt, 2002). SLF experiments undertaken with spinning frequencies of 60 kHz or greater have been cross-polarization-based (CP) (Chevelkov et al., 2009; Chevelkov et al., 2010; Paluch et al., 2013; Paluch et al., 2018) or use phase modulated rotary resonance pulses (Liang et al., 2021). Symmetry-based recoupling is comparable to CP based SLF experiments and has many of the same disadvantages, but symmetry be advantageous in a few ways. First, the symmetry sequences can be constructed to be very selective of the terms allowed, where only a few terms in the NMR Hamiltonian are still active. Second, there is only one channel that has high-power pulses applied which limits the power deposition, where the CP methods apply high fields on both channels simultaneously.
In this work, we introduce an approach to generate pulse sequences with optimized recoupling at fast spinning given probe performance requirements. Candidate symmetry sequences are generated, the scaling factor (κ) is optimized in silico using variable flip angle pulse sequence elements, and then the highest performing sequences are selected. The most promising sequences are tested against B1 inhomogeneity and match condition mis-set. The experimental performance of several candidate sequences is evaluated on a favorable model sample, the micro-crystalline protein GB1 (β1 immunoglobulin binding domain of protein G) which is uniformly labelled with 2H, 13C, 15N, and then back-exchanged with 1H at all exchangeable sites, and on a more challenging >300 kDa precipitated complex of GB1 with immunoglobulin G (IgG). We have previously investigated differences in protein dynamics for protein GB1 in these two environments. The analysis of various relaxation and relaxation dispersion experiments has indicated that while ps-ns motions and some µs motions appear to be largely similar for GB1 in the two environments, there appears to be an additional overall motional mode present only in GB1 in the complex with IgG (Lamley et al., 2015a; Öster et al., 2019). Thus we decided to investigate whether we see the presence of this additional dynamic mode reflected in the measured order parameters. In addition, we have previously performed variable temperature dynamics measurements on crystalline GB1 and in the analysis assumed the observed trends are dominated by changes to the time scales of the motions rather than changes in amplitude and assumed a constant order parameter (Lewandowski et al., 2015; Busi et al., 2018). The initial variable temperature molecular dynamics simulations suggested that indeed the 15N-1H order parameter changes only very slightly with temperature in the explored 30 °C range but we thought this study to be a good opportunity to begin to explore validity of such approximation experimentally. Incidentally, crystalline GB1 and GB1 in the complex with IgG cover the range of a favorable model sample and a challenging “real” sample.
Symmetry Based Pulse Sequences
Symmetry-based sequences allow for the selection of portions of the full NMR Hamiltonian (Levitt, 2002). The performance of the pulse sequence with regards to the extent of the reintroduction for an interaction is indicated by a scaling factor (κ). The scaling factor is the magnitude of any coupling when compared to the static limit, which can vary between 0 and 1, where a larger κ indicates a more efficient recoupling/reintroduction. The recoupling performance can be altered in two ways: by using a different symmetry, or by using a different rotation element. This work demonstrates a strategy to find high performance heteronuclear recoupling pulse sequences by exploring the possible variations of symmetry derived pulse sequences.
The primary limitation for the application of symmetry at high spinning frequencies is the electronic performance of the probe. The nutation of the spins which are induced by the applied radio frequency field in the probe must match the conditions specified by the symmetry sequence, which is dependent on the spinning rate. The same symmetry will require higher radio frequency fields with higher spinning rates. For example, in the well-known R1817 sequence (Levitt, 2002; Wylie and Rienstra, 2008; Wylie et al., 2011) the match condition is 9 times the spinning rate, which means that ν1 = 90 kHz at 10 kHz magic angle spinning (MAS) rate, but ν1 = 540 kHz at 60 kHz spinning frequency. The electronics in the probe will likely break down due to the high voltages needed to generate such strong pulses, or alternatively, the protein sample will denature when the temperature gets too high from radiative heating. Most modern probes are specified to work with an applied 1H field of ∼100 kHz for long pulses (i.e., decoupling during acquisition), with fields on the order of 250 kHz available for short periods (and nonconductive samples). Since the prior art does not produce many appropriate choices, a method to identify and optimize symmetry-based recoupling to measure the heteronuclear dipole-dipole coupling was devised.
Initially, we intended to apply a series of symmetries identified for use at 40 kHz MAS (Hou et al., 2011). We simulated these schemes at 60 kHz MAS with a standard π pulse as the R-element (see Figure 1; Supplementary Figure S1, yellow circles). The curve of scaling factor against match condition matches the literature well but the scaling factor is smaller at 60 kHz spinning than at 40 kHz (Hou et al., 2011). That is, the scaling increases with increasing field, approaching some upper limit (here, κ ∼ 0.3). We chose the best candidate symmetries to test experimentally where the criterion for selection was an applied field ν1 < 130 kHz. After careful calibration of the applied field using long nutation experiments, the performance of the candidate sequences was found to be poor, with little to no recoupling apparent. The disappointing performance was attributed to poor 1H channel B1 homogeneity of the probe. The ratio of the NMR signal intensity at the 90°, 450° and 810° pulse is used as a proxy of an actual B1 homogeneity measurement. The probes used in this study typically showed I450/I90 ∼70%, and I810/I90 ∼55% for the 1H channel.
FIGURE 1. Heteronuclear dipole-dipole recoupling scaling factor for symmetries found in Hou et al. (2011). The full range of symmetry match conditions for standard composite pulses. Yellow circles are for R-symmetries with a standard rectangular π pulse. Red triangles report the scaling with a compound R-element of a π/2(0)−π/2(90)−π/2(0) with the equivalent match of a 3π/2 pulse. Orange squares report the scaling with a compound R-element of a π/2(0)−3π/2(180) with the equivalent match of a 2π pulse. See Supplementary Figure S1 for full B1 range.
Composite pulses were implemented to compensate for the probe deficiencies since the standard symmetry recoupling experiments were unsuccessful. Two composite-π pulses are prominent in the literature, the [90(0)90(90)90(0)] and [90(0)270(180)] where the flip angle of a portion of the pulse is denoted by the large number and the phase of that portion is denoted by the subscripted number in parentheses (Levitt, 2002). The comparison of the scaling factors determined in silico by simulations with the SIMPSON program (method described below) shows a dramatic change in the scaling factor when the symmetry element is changed. The performance of the [90(0)90(90)90(0)] composite pulse (Figure 1; Supplementary Figure S1, orange squares) follows the same general trend as the standard π pulse (yellow circles), but with worse efficiency and higher applied field. However, the [90(0)270(180)] composite pulse (Figure 1; Supplementary Figure S1, red triangles) does not follow the same trend. The shape of the curve produced by the [90(0)270(180)] composite pulse has a maximum in the curve, whereas the other pulses asymptotically rise. The maximum scaling factor is found at a relatively low field, with performance similar to the π pulse variant.
While it is made clear in the literature that the specifics of the R-elements contribute to the efficiency of the recoupling by altering the scaling factor (κ), the magnitude of this contribution was underappreciated. These preliminary simulations demonstrate that more symmetries than those found in the literature should be tried. Those with lower match conditions can be a viable option with composite pulse rotation elements. Also, a variety of symmetry elements will be useful to identify the best candidate sequences given the desired experimental conditions. To these ends, we present tools to generate appropriate symmetry lists, tools to test these symmetries, and experimentally test the best candidates.
Generating Candidate Symmetry Sequences
An
where kp is determined by the specifics of the R or C-element.
A custom program was written in Python 3 (see Supplementary Material: SymmetrySelector/SymSelect.py) to generate new symmetry sequences that are more relevant to faster spinning rates. This program reproduces the full list of sequences from Levitt (2002) (excepting 3 minor book-keeping errors, see symmetry selection rules section of the Supplementary Information; Supplementary Tables S2–S4). An arbitrarily large number of candidate symmetry sequences for application at 60 kHz spinning frequency is then generated. Symmetries for most terms in the Hamiltonian can be generated, with the possibility to limit the output based on experimental considerations. 203 symmetries fit the selection criteria (Supplementary Table S5); the sequence must allow heteronuclear dipole-dipole couplings, disallow homonuclear dipole-dipole coupling, and the base applied field is between 0.1 and 150 kHz. It should be noted that this program can be used to generate symmetry sequences for most spin ½ Hamiltonians, not just the heteronuclear dipole-dipole coupling.
Composite Rotation Element Pulses
Complex inversion pulses have been used for R-sequences, such as adiabatic inversions (Herbst et al., 2011; Herbst et al., 2015), numerically optimized optimal control pulses (Nielsen et al., 2009), and composite inversion pulses (Levitt, 1982). However, the R-element seems to have previously been chosen for some desirable property of the composite pulse and the expected shortcomings of the symmetry sequence, not explicitly for the performance of the sequence.
Since the timing and phase behavior of the multiple-pulse sequence must still fulfill the symmetry requirements the match field grows in proportion to the total arc swept out by the composite pulse. This term is the pulse contribution to the match condition “kp,” where:
For example, a simple π inversion pulse sweeps an arc of 180°, so kp = 1. For the composite pulses [90(0)180(90)90(0)] and [90(0)270(180)] the RF sweeps out 360°, but the result is only an inversion of the polarization. These specific examples result in kp = 2. The match field B1 is, likewise, twice that of the symmetry alone.
It is unclear which composite pulses produce the highest scaling factor for the smallest applied field. The sweep angles for compound pulses have previously been constrained to only use the principal directions, i.e., they only include integer multiples of 90°, but such a constraint is not enforced in this study. For example, when we modify the [90(0)270(180)] composite pulse we trade the initial 90(0) portion for a θ(0) pulse, and the 270(180) becomes (180 + θ) (180). The portions of the composite pulse are consecutively numbered τ1,2…n.The field, timing dependencies, and a diagram of this element, named “1a,” and two others are shown in Table 1. The “2a” element is a slight variation on 1a, it has an extra “wiggle” before finishing. The “4a” element is a variation of the 90(0)180(90)90(0) element where the middle, out of plane, pulse is allowed to vary. It is possible to change θ continuously, and smoothly for these compound pulses and the end point (inversion) will not be changed. In total, 22 R-elements were constructed where the sweep angle (and thus the applied field) can be varied continuously, but which always produces a traceable inversion pulse. A simple inversion occurs when θ is zero for many, but not all the R-elements. All 22 of the R-elements with their pulse timings, field match dependence, and a visualization is found in Supplementary Table S6.
TABLE 1. Selected variable flip angle R-elements, with field and timing dependencies. Red arrows represent the first portion of the composite pulse τ1, blue is the second, τ2, and yellow is the third portion, τ3.
Simpson Simulations
Numerical simulations were conducted to determine the scaling factor with the SIMPSON NMR calculation software (Bak et al., 2000). The numerical scaling factor is comparable to the analytical solution presented in Levitt (2002), but differs in magnitude by up to 0.05, where the numerical method always over-estimates the scaling in comparison to the analytical solution. The spin system for testing the scaling factor is an isolated two spin (1H-15N) system with the dipole-dipole coupling defined to be 10 kHz which corresponds to an internuclear distance of 1.07 Å. The single crystallite crystal file “alpha0beta90” was used so that the maximum dipole-dipole coupling will be produced using the method of Hou et al. (2011) (via personal correspondence).
A series of simulations were performed where the applied field was varied from the lowest match condition up to 250 kHz for the composite R-elements to produce dipole recoupling efficiency curves. These curves are well-defined, as demonstrated by the simulations of the symmetry sequence R22219 in Figure 2. It is convenient to vary the applied field and back calculate the variable flip angles and pulse durations to maintain a valid range for the applied field. The choice of the recoupling element results in decoupling in a sequence which is supposed to recouple, which may point to further selection rules for symmetry elements since there are several instances which result in a coupling of zero. These zero-coupling points might be useful for other applications such as decoupling. These zero points are also present with more crystallites, and thus do not appear to be an artifact of the single crystal simulations. The shape of the curve is an indication of the sensitivity to, for example, B1 inhomogeneity or a mis-set match condition, where a flatter curve should mean less sensitive condition. There is a significant zero-frequency in the simulated spectra indicating that some fraction of the polarization is not recoupled, although this component is lost when more crystallites are included in the calculation when the calculations are repeated with crystal averaging, as seen in Supplementary Figure S2.
FIGURE 2. Dipole-dipole recoupling spectra of the R22219 symmetry as a function of ν1 match condition for different variable flip angle pulse sequence elements.
Dipole recoupling efficiency curves are generated for all 22 R-elements with 101 applied fields for each of the 203 symmetries. The highest scaling factor for each symmetry and element pair are plotted in Figure 3. Here, the recoupling for the “m” and “μ” (space and spin component) quantum numbers are either correlated (orange X), or anticorrelated (blue +), where this correlation indicates a second order dependency on frequency offset from the carrier (correlated) or on a field mis-set as in B1 inhomogeneity (anticorrelated), where the frequency offset dependence is the more favorable deficiency.
FIGURE 3. Maximum scaling factor for heteronuclear dipole-dipole scaling factor determined in silico of all test symmetries and pulse sequence elements as a function of (a) applied (B1) field and of additional tilt angle (θ). Blue“+”’s indicate sequences where “m” and “μ” are correlated, and orange “x”s are anti-correlated. Black circles indicate the position of sequences chosen for experimental verification.
In Figure 3A, there is a clear limit to the performance as a function of the applied field, where the maximum scaling rises quickly up to a maximum of about κ = 0.27 and ω1 = 175 kHz and then slowly reduces as the applied field increases. There are many candidates with κ > 0.225 and ω1 < 150 kHz, which are suitable for further testing. When the scaling is plotted against the flip angle, we find that there are certain flip angles that are favored. The first local maximum is at zero, indicating that at least some sequences do not improve using composite pulses. There are local maxima at about 15° and 35°, but the global maximum occurs when there is a flip angle of about 125°, and then approximately every 90° after.
The SIMPSON input files used to evaluate the pulse sequences and Python scripts to process the simulated datasets (i.e., find the maximum in a simple spectrum and report the scaling factor of a 2D simulation) are available online (see Materials and Methods).
The 1a element was found to produce the largest scaling factor (or is tied for the largest) for any given symmetry sequence. Additionally, the 1a element has the lowest requirements of any element that performs similarly. That is, the 1a element consistently produces the highest scaling factor for a given symmetry for the least applied field amongst all other competitive R-elements. Therefore, we only considered the 1a element further. The plot of the scaling factors for only the 1a elements as a function applied field and flip angle can be found in Supplementary Figure S3. The intensity of the applied field during the recoupling period is the largest concern, however the maximum duration is short so we felt that we could push the limits of the probe and chose an upper limit of about 150 kHz. Since compensating for poor B1 homogeneity was the motivation to use compound pulses, we devised a fast test to show the B1, or mis-set, dependence. The 2-spin simulation is run with the optimum flip angle, but the applied field is multiplied by 0.9 or 1.1. About half of the sequences respond strongly and move by hundreds of Hz, while the other half move less than ∼50 Hz. The 0.9 mis-set spectra generally shift more than the 1.1. We found that the mis-set dependence depends on the specific components that are recoupled, specifically if the space component “m” and the spin component “μ” are correlated (both positive or both negative within the same allowed set), the mis-set dependence is typically small. This is due to a “pulse” term in the second order Hamiltonian which is allowed in the anticorrelated set of symmetries and disallowed in the correlated set. The “pulse” term is replaced by a second order frequency offset term in the correlated set of symmetries, which can be demonstrated in the simulations by observing the dependence of the scaling when introducing a frequency offset. This is not an absolute rule, though, since some anti-correlated symmetries are not greatly affected by mis-set, such as R1485.
The criteria to choose a symmetry sequence are to find a sequence that works and will not damage either the sample or the hardware. Initially, the criteria were that the maximum applied field should be less than ∼150 kHz, we should use the “1a” element for homonuclear recoupling, the scaling factor should be as large as possible, the spectral width should be large enough that the spectrum does not fold, and the echo time should be short. Ideally, the allowed “m” and “μ” components should be the same sign and the R-element should not stray too far from those previously devised. On close inspection, one notices a small gap at around 150 kHz, after which the scaling factors no longer greatly improve. The sequence immediately after this gap is R2098(124(0)304(180)). Reducing the flip angle to 115° [making R2098(115(0)295(180))] reduces the applied field to 152 kHz, but does not affect the scaling factor [the scaling is the same down to a 105° flip angle (B1 = 145 kHz)]. R2098(115(0)295(180)) thus fulfills all our desired traits. There are many candidate sequences with good scaling, and less demanding match conditions that meet the criteria, but were not tested experimentally. However, the sequences R22219 (300(0),480(180)) and R14225 (460(0),640(180)) were chosen to test the robustness of the simulations approach, since the unusual flip-angles in the composite pulses change the scaling factors from almost 0 to above 0.2.
This same in silico analysis can, of course, be made under 100 kHz spinning (Supplementary Figure S4). 310 gamma encoded candidate symmetry sequences were generated where N = 10 through 42, n = 1 through 37, and ν ≤ N/2 and the base match field is limited so that ω1 < 200 kHz. The recoupling is less efficient at 100 kHz spinning, as the curve equivalent to the one shown in Figure 3A is shifted to higher match conditions at 100 kHz spinning (Supplementary Figure S4A). The dependence on the flip angle is the same as at 60 kHz spinning, where there are local maxima at ∼15°, 35°, 125°, and then approximately every 90° afterwards (Supplementary Figure S4B). Still, there are several candidate sequences with scaling between κ = 0.15 and κ = 0.20 with relevant match conditions. If we limit the search to the “1a” element, with a match condition of less than ω1 = 150 kHz we find 8 candidates, R22161(θ = 104), R26191(θ = 105), R30221(θ = 106), R32232(θ = 103), R34251(θ = 104), R38281(θ = 102), R40292(θ = 104), and R42311 (θ = 97).
Experimental Verification
A cross-polarization-based 1H-detected 1H-15N correlation experiment was converted into a 3D experiment, where the third dimension consists of a constant time echo period on the low-gamma frequency (Figure 4). The recoupling time increases to create the third, separated local field dimension. The echo period (τecho) is calculated according to the chosen symmetry and data sampling by the following equation:
FIGURE 4. The pulse sequence used for the experiments. Thin lines represent 90° pulses, and thick lines represent 180° pulses, boxes indicate cross polarization spinlock pulses. The sequence is modified from the usual solid-state NMR CP-HSQC by the addition of an echo period, and the application of the recoupling period. ϕ0 = 1; ϕ1 = 0 0 2 2; ϕ2 = 0*4 2*4; ϕ3 = 1*8 2*8; ϕ4 = 3 1; ϕ5 = 1; ϕ6 = 0; ϕ7 = 1; ϕrec = 1 3 3 1 3 1 1 3 3 1 1 3 1 3 3 1 where X*4 indicates that phase X is repeated 4 times. The phase shifts for the recoupling experiment are determined by the symmetry sequence and must be adjusted by hand.
Where “n” is the space winding number from the symmetry sequence, “τr” is the rotor period, “k” is the number of repeated symmetry sequences (usually 1) to better fit the spectral width, and “points-1” is the zero-time-point inclusive number of points. We have chosen to apply the recoupling sequence only during the first half of the echo, and a standard decoupling pulse train for the rest of the echo and chemical shift dimensions. That is, the recoupling sequence increments up to the echo pulse, and the echo time is constant throughout the experiment. This is not the ideal case since the echo period should ideally be as short as possible so as little signal is lost as possible. However, the logic needed to loop the symmetry elements before, after and during the π pulse, to properly invert the phases of the symmetry pulses after the echo, and to maintain the proper timings in the Bruker scripting language was too cumbersome, so a simple echo was settled on since the Nitrogen T2* is quite long for these samples. It might be possible to use “compound pulse decoupling” pulses to simplify the logic, but we were unsuccessful in our attempts. Alternatively, the π pulse could be used for chemical shift evolution in a constant time evolution experiment, but this will reduce the sensitivity further.
Three recoupling sequences with good theoretical scaling factors and appropriate match conditions were chosen to validate our approach, and to determine which candidate scheme is the most promising. These sequences: R2098 (115(0),295(180)); R22219 (300(0),480(180)); and R14225 (460(0),640(180)) were tested with both their π-pulse version and the numerically optimized sequence elements (Supplementary Figure S5). The sequences were chosen partially because of the diversity of the flip angle, applied field, and the difference in scaling between standard and optimized sequences. The in-silico performance of these three sequences, along with R1485 (115(0),295(180)) are summarized in Table 2. The experimental performance closely follows the in-silico performance as demonstrated by strong recoupling with the variable flip angle pulses which lends credence to the approach, and the robustness of symmetry theory in general.
While the R2098 (115(0),295(180)) seems to be a great option, its match condition is considerable at ν1 = 152 kHz. During the power calibration with long nutation experiments, the signal disappeared in a few ms at the fields needed for R2098 (115(0),295(180)). The probe detuned indicating that such high match conditions would likely damage the sample or the probe, however fields up to about 130 kHz were long lived. The R14225 (460(0),640(180)) was also tested, but did not produce any recoupling, so the R22219 (300(0),480(180)) sequence was also thrown out since both have quite long echo times, and quite large flip angles. Other candidate sequences were identified from the scaling curve that have more conventional R-elements and lower match conditions. Amongst a handful of candidate sequences, the R1485 (115(0),295(180)) sequence was the first one that worked (it was also the first one tested). The R1485 (115(0),295(180)) has a lower match condition than the R2098(115(0),295(180)), the R-element is not far from a canonical R-element, and it seems to have a high tolerance for field mis-set in numerical simulations. The match condition (ν1 = 119.5 kHz) for R1485 (115(0),295(180)) is very near the ν1/2p = 2ωr rotary resonance condition, although this does not seem detrimental to the quality of the data in the microcrystalline sample. If the match is a greater concern, the R-element could be adjusted for a lower match condition as it was for the R2098 (115(0),295(180)). The scaling factor is the same for all “1a” R-elements between θ = 100° and θ = 115°, where ν1 ranges from 111.8 to 119.5 kHz. The ability to turn down the power was not appreciated at the time the experiment was conducted.
Temperature Dependent Order Parameters
To validate the designed experiments, we have first applied them to a favorable model sample of crystalline 100% H2O [U-2H,13C,15N]GB1, which means that the protein is uniformly 13C and 15N labelled and perdeuterated with only exchangeable protons being reintroduced at 100%. We used R1485 (115(0),295(180)) sequence to measure 15N-1H order parameters at three different temperatures, nominally 263.2, 273.2, and 283.2 K, which correspond to sample temperatures of ∼302 K, ∼309 and ∼315 K (larger temperature differences are difficult to obtain on our 1.3 mm probes at 60 kHz spinning; see Materials and Methods). A representative 1H-15N correlation spectrum taken from the first plane of the 3D can found in the Supplementary Figure S9.
The quality of the recoupled line-shapes for all temperatures is excellent as evidenced in Figure 5A. There is surprisingly little intensity at the zero-frequency, and the dipole line-shape is clear and strong. As expected, there is little difference in the experimentally observed (fitted) coupling over the explored temperature range and the obtained 15N-1H order parameters are generally very similar for all three measurements (see Supplementary Table S7). We only observe a noticeable change of 15N-1H dipole-dipole coupling (from 11.1 to 10.2 kHz) as a function of temperature for L5, D36, and T49 which are low intensity peaks in the 2D fingerprint spectrum. In the ∼309 K data, the order parameter for D36 is spuriously low for reasons which are not clear (See Supplementary Figure S8). The precision of the experiment is very good: the 1σ standard deviation in the dipole couplings is around ±220 Hz, with some down to ±100 Hz, and outliers ranging from ±500 to 1500 Hz. These errors are comparable to previous symmetry-based methods on the same protein (Franks et al., 2005; Franks et al., 2006b). However, the precision may be overstated due to poor noise estimation in the Monte Carlo error analysis. Monte Carlo error analysis first finds the best fit for the experimental data to a simulation, here the time-domain dipole-dipole recoupling curve. Noise is then added to each point of the experimental data and the best fit is found again, and the values saved. Noise is added to the original experimental data several times (here, 5,000 times) to estimate the amount of spread in the simulation values found in the experimental data. The method we have used to add the noise could be improved. First, most peaks have very similar initial intensities, so a constant noise value was used for all peaks in the Monte Carlo analysis, which was ±7.5% of the total. All trajectories were normalized to 1 during the integrations, which results in an undesirable loss of information regarding the intensity. The result is that a resonance that is 100 intensity units high will have noise ranging from −7.5 to +7.5 added, while a peak that is only 10 Intensity units high will only have noise added that ranges from −0.75 to +0.75 units. The noise estimate works well for peaks with a typical intensity, and may even be larger than necessary, but the error analysis fails for peaks with poor intensity. A typical Monte Carlo fit is shown for residue K28 in Supplementary Figure S6. Similar figures for all fits can be found in their corresponding datasets in the online materials. Those residues with worse sensitivity will have spuriously good fits, such as found for the K13 peak in the GB1+IgG complex (Supplementary Figure S8).
FIGURE 5. 15N-1H dipolar couplings as a function of temperature in crystalline GB1. (A) Representative dipole-dipole recoupling spectra using the R1485 (115(0),295(180)) pulse sequence at ∼302 K, ∼309 and ∼315 K at 60 kHz spinning. (B) S2NH obtained from the measured dipolar couplings, error bars are drawn at 1σ. (C) S2NH obtained from variable temperature molecular dynamics simulations. MD derived order parameters were extracted from 400 ns molecular dynamics simulations of a 3 × 3 × 3 supercell of GB1 containing 108 monomers.
The 15N-1H dipolar order parameters shown in Figure 5B fit well with those reported previously (Franks et al., 2005; Franks et al., 2006b). The 15N-1H dipolar order parameters are generally around 0.9 (order parameters, S2, take values between 0 and 1, which mean unrestricted motion and no motion respectively) with an occasional dip to around 0.8 near the loops. The largest amplitudes of motions are observed for residue G41, which is in a loop between the alpha helix and beta strand 3.
We have also compared the experimentally determined 15N-1H order parameters to those obtained from molecular dynamics simulations performed in the 280–310 K temperature range (Figure 5C). The simulated rates are in general good agreement with the experimental ones (see Supplementary Figure S7). As before (Busi et al., 2018), the MD simulations predict that there should be little change in the order parameters of the expected temperature range.
Order Parameters of the GB1-IgG Complex
The complex of 100% H2O [U-2H,13C,15N]GB1 with IgG is a much more challenging sample compared to crystalline GB1 to apply the described methods. As a precipitate it is more heterogenous, >90% of the sample volume is taken by the antibody resulting in lower sensitivity and GB1 in the complex exhibits much more pronounced slow motions in the microsecond range. (Lamley et al., 2014; Lamley et al., 2015a; Öster et al., 2019). A paramagnetic doping agent, 2 mM Gd (DTPA-BMA), was added to speed up the measurements by reducing the required relaxation delays (Linser et al., 2007), although in this particular application sample heating is a large concern so the recovery delay remains quite long at 1s. The 1H-15N correlation spectrum taken from the first plane of the 3D can found in the Supplementary Figure S9.
Despite the more challenging nature of the sample, the GB1 in the complex still produces high-quality dipole-recoupling spectra, as seen in Figure 6A. There is a significant zero-frequency component in most of these spectra, and the sensitivity is generally worse (especially for T16). The origin of the zero-frequency component is possibly due to the larger amount of 1H atoms in the sample, less efficient 1H-1H homo-decoupling, increased dynamics, sample heating, and/or probe detuning. The determined 15N-1H order parameters are generally lower in the GB1 complex (average S2NH ∼0.7) than in crystalline GB1 (average S2NH ∼0.9). If one uses expression for order parameter in diffusion in a cone model this difference corresponds on average to ∼23° additional motional amplitude for most residues. This provides further support for presence of a microsecond range overall motion of GB1 in the complex with IgG as proposed previously (Lamley et al., 2015a).
FIGURE 6. 15N-1H dipolar order parameter data for GB1 in a crystal (black line) and in a precipitated complex with IgG (red line). (A) Overlay of representative dipole recoupling spectra. (B) 15N-1H order parameters as a function of residue, error bars are drawn at 1σ. Both data sets were recorded at the nominal temperature of 263.2 K corresponding to sample temperature of ∼302 K at 60 kHz spinning frequency.
Conclusion
We have presented a method to apply symmetry-based recoupling theory to fast MAS experiments using variable flip angle compound pulses. The method generates many candidate sequences that have a reasonably high scaling factor and applicable match conditions. Being able to apply symmetry principles under fast-MAS makes SLF methods applicable to faster spinning and the other benefits that usually comes with it. The heteronuclear dipole-dipole coupling can be measured site-specifically in microcrystalline GB1 and in the GB1-IgG complex. The order parameters found using the variable flip angle method are consistent with previous datasets and molecular dynamics in the microcrystalline sample. The GB1 in the complex shows both localized differences in dynamics and global increase of cumulative motional amplitudes.
Materials and Methods
Sample Preparation
Uniformly (2H, 13C, 15N) GB1 was produced as described previously (Franks et al., 2005; Zhou et al., 2007b). After production in 2H buffer with 13C glucose and 15N NH4Cl, the protein is placed in 1H containing buffer, and heated so that the exchangeable 1H sites are 1H labelled. The protein is then either crystallized or incubated with natural abundance IgG in an equimolar ratio (Lamley et al., 2014). All back-bone amide sites on the GB1 molecules are thus labeled, but the strong 1H dipole coupling network is disrupted locally by the 2H labelling of the sidechain carbons. The buffer of the IgG-GB1 complex contained 2 mM of paramagnetic Gd (DTPA-BMA).
NMR Spectroscopy
All experiments were performed on either a Bruker Avance III spectrometer at 700.13 MHz 1H Larmor frequency or a Bruker Avance II spectrometer or at a 599.4 MHz 1H Larmor frequency. A Bruker 1.3 mm HCN Probe operating in HCN triple resonance mode with a sample spinning rate of 60 kHz ± 3 Hz was used with both instruments. 1,200 L/h of cooling gas was used at the nominal temperatures of 263.2, 273.2, and 283.2 K. The nutation frequencies for the 90° pulses were calibrated so that the hard pulses for 1H were 2.1 μs (ν1 = 120 kHz); 13C, 2.5 μs (ν1 = 100 kHz); and 15N, 3.25 μs (ν1 = 77 kHz). The 1H carrier radiofrequency (RF) was centred on the H2O signal (∼4.7 ppm), 15N at 120 ppm, and 13C at 100 ppm. Heteronuclear 1H decoupling (∼10 kHz SPINAL-16) (Fung et al., 2000) was used during the indirect chemical shift dimension, and during the echo period when the recoupling sequence was not being applied, and approximately 10 kHz WALTZ-64 (Zhou et al., 2007c) 15N decoupling was used during 1H acquisition. The States-TPPI method was employed for quadrature detection in the indirect chemical shift dimension (Marion et al., 1989) and only the real portion of the dipole-coupling dimension is acquired. The MISSISSIPPI (Zhou and Rienstra, 2008) solvent suppression scheme was applied with a spinlock field of ∼30 kHz for four 15 ms intervals before detection. The symmetry match condition was calibrated to the theoretical value by varying the applied power in a 2D, nitrogen edited, 1H nutation experiment until the nutation experiment is within 1 Hz.
The total experiment time for each temperature point of the crystalline GB1 was 16 h each 1H free induction decay was acquired for 40 ms with a spectral width of 30 ppm with 32 coadded transients. The 15N dimension for the microcrystalline protein were acquired with 80 rows with a dwell of 333.33 µs for a total of 13.3 ms in the indirect dimensions. The R1485 (115(0),295(180)) dimension was acquired for 15 real points with an increment of 8*τr = 133.33 µs for a total of 1.87 ms (3.73 ms total echo time). The recovery delay was 1.5 s.
The spectrum of the GB1-IgG complex was collected in four blocks of 34.1 h each which were later summed together, for a total of 5 days 16.5 h. There was 1,200 L/h of variable temperature gas flow at the nominal temperature of 263.2 K. Each 1H free induction decay was acquired for 30 ms with a spectral width of 30 ppm in four blocks of 128 coadded transients (512 total), the 15N dimension for the microcrystalline protein were acquired with 64 rows with a dwell of 333.33 µs, with a spectral width of 42 ppm (15N) for a total of 10.7 ms in the indirect dimensions, the R1485 (115(0),295(180)) dimension was acquired for 15 real points with an increment of 8*τr = 133.33 µs for a total of 1.87 ms (3.73 ms total echo time) with a relaxation delay of 1 s All 3D data was processed using NMRPipe and the four blocks were added using NMRPipe (Delaglio et al., 1995). The dipole-recoupling dimension was not Fourier-transformed in NMRPipe so that peak volumes could be extracted. The Fourier-transform of the dipole recoupling dimension was performed on the peak volumes extracted by NMRPipe by a the fast fourier transform routine found in python’s numpy package. The imaginary components of the trajectories were filled with zeroes, and the trajectory was zero filled to 512 points. The 2D datasets with a dipole recoupling dimension (Supplementary Figure S5) were processed in Topspin with by performing a Hilbert-transform to fill the imaginary portion of the dipole trajectory before Fourier transforming the dipole recoupled dimension. (See Supplementary Material: BrukerMacros/2DHN).
External KBr (Thurber and Tycko, 2009) and neat methanol (Ammann et al., 1982) were used as external standards to calibrate the temperature. The samples did not have adequate resolution to unambiguously identify the bulk water signal from the isopropanol and methyl-pentane-diol OH signals, which precluded temperature calibration by the chemical shift difference between water and DSS (Hoogen et al., 1988; Wishart et al., 1995). Temperatures derived using either the chemical shift or the T1 of KBr (Thurber and Tycko, 2009) were not self-consistent on the 700 MHz instrument at either 10 or 60 kHz spinning. There is approximately 10°C difference between the two methods. The T1 method indicates a 20°C difference across the nominal temperatures, where the chemical shift method indicates a difference of 13.6°C. Calibrating the temperature by the chemical shift difference in the 1H spectrum of methanol (Ammann et al., 1982) indicates temperatures of 301.9 (263.2), 309.0 (273.2), and 314.9 K (283.2) or 28.7, 35.8, and 41.7°C under the experimental conditions, for a range of 12.9°C.
Data Processing, Fitting and Simulations
The GB1 resonances were identified in the 2D finger-print spectra from the literature and previous work (Franks et al., 2005; Zhou et al., 2007a). Peaks were integrated using NMRPipe (Delaglio et al., 1995), where the peak volumes were converted into dipole recoupling trajectories. The volume of the peak on the first plane is used to normalize the rest of the curve, so all curves range between ±1, and start with an intensity of 1.00. Overlapping resonances were fit, but these resonances were not included in the figures or analysis since the resonance could not be unambiguously identified.
A library of numerical simulations was created in SIMPSON 4.2.1 on an Apple MacBook Pro for use in the Monte-Carlo fitting routine. The simulation library was created using 251 dipole-dipole couplings ranging from 7,500 Hz to 12,500 Hz in steps of 20 Hz, the spin rate was 60 kHz, the calculation method was “direct”, the crystal file was “zcw376,” and 16 gamma angles were used. The time-domain trajectories and frequency domain spectra were saved as a 2D SIMPSON file.
5000 Monte Carlo steps were used for error analysis for all datasets only using the time domain SIMPSON library. The library was expanded using simple operations on the time-domain. The DC parameter is used to add a constant is to all data points (varied between −0.2 and +0.5). The scaling factor multiplies each point by a constant (varied from 0.9 to 1.1). Finally, relaxation is simulated by applying line-broadening for each trajectory, that is, the simulation is multiplied by a time dependent exponential function (from 0 Hz to 2,500 Hz). To Fourier transform the dipole trajectories, the imaginary time portion was filled with zeroes, and the trajectory was zero-filled to 128 points. For the crystalline GB1, most line shapes fit with a small, negative DC offset (−0.05), a scaling multiplier of 1.00, and less than 300 Hz of line-broadening. For the IgG-GB1 complex, there is a small positive DC offset (+0.14), a scale of 1.00, and approximately 600 Hz of line-broadening on average. The rigid limit for the N-H dipolar coupling is taken as 11,477.3 Hz to determine the Order Parameters, which corresponds to an N-H bond length of 1.02 Å.
Molecular Dynamics Simulations
A molecular dynamics trajectory for a 3 × 3 × 3 supercell of GB1 containing 108 monomers was computed using AMBER MD (Case et al., 2005; Doshi and Hamelberg, 2009; Maier et al., 2015; Tian et al., 2019). The coordinates of the X-ray structure of GB1 (PDB:2gi9 Franks et al., 2006b) were taken as a starting conformation. To the supercell, 108 PO43− counter ions were added. 12,852 explicit water molecules were added, followed by charge balancing with sodium ions giving an overall box size of 75.591 Å × 107.152 Å × 150.822 Å. The ff19SB (Tian et al., 2019) forcefield was used for the GB1 proteins, with OPC water (Izadi et al., 2014) and GAFF cocrystallites (Wang et al., 2004). After minimization, the system was replicated and heated to the temperatures indicated in the figures (280, 290, 300, 310 K). The systems were then simulated for full 400 ns runs. For each, a 2 fs timestep was used with a cut-off of 11 Å for non-bonded interactions. Temperatures were maintained using a Langevin thermostat, and the SHAKE algorithm (Ryckaert et al., 1977) was applied to all bond lengths involving a hydrogen atom. Anisotropic pressure scaling was used with periodic boundary conditions.
Prior to processing, the Cα carbons between timesteps were aligned using cpptraj (Roe and Cheatham, 2013). Then, correlation functions for each N-H vector were calculated according to the iRED framework using cpptraj (Prompers and Brüschweiler, 2002). The median was calculated for each residue over all GB1s in the supercell for which order parameters could be extracted, and the error taken as twice the median absolute difference.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
WF designed, optimized, and performed the experiment, analyzed the simulated and experimental NMR data, and prepared the manuscript. BT performed and analyzed the molecular dynamics simulations. JT performed and analyzed NMR simulations. JL designed the experiment, analyzed the data, and prepared the manuscript.
Funding
The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement 639907, which supported JL and WF. JL also acknowledges funding from EPSRC EP/L025906/1 and BBSRC Grant BB/R010218/1.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
BT thanks Bruker and the University of Warwick for PhD funding as part of the Warwick Centre for Doctoral Training in Analytical Science.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.791026/full#supplementary-material
References
Alderson, T. R., and Kay, L. E. (2021). NMR Spectroscopy Captures the Essential Role of Dynamics in Regulating Biomolecular Function. Cell 184, 577–595. doi:10.1016/j.cell.2020.12.034
Ammann, C., Meier, P., and Merbach, A. (1982). A Simple Multinuclear NMR Thermometer. J. Magn. Reson. (1969) 46, 319. doi:10.1016/0022-2364(82)90147-0
Asami, S., and Reif, B. (2012). Assignment Strategies for Aliphatic Protons in the Solid-State in Randomly Protonated Proteins. J. Biomol. NMR 52, 31–39. doi:10.1007/s10858-011-9591-4
Baek, M., DiMaio, F., Anishchenko, I., Dauparas, J., Ovchinnikov, S., Lee, G. R., et al. (2021). Accurate Prediction of Protein Structures and Interactions Using a Three-Track Neural Network. Science 373, 871–876. doi:10.1126/science.abj8754
Bak, M., Rasmussen, J. T., and Nielsen, N. C. (2000). SIMPSON: A General Simulation Program for Solid-State NMR Spectroscopy. J. Magn. Reson. 147, 296–330. doi:10.1006/jmre.2000.2179
Barbet-Massin, E., Pell, A. J., Retel, J. S., Andreas, L. B., Jaudzems, K., Franks, W. T., et al. (2014). Rapid Proton-Detected NMR Assignment for Proteins with Fast Magic Angle Spinning. J. Am. Chem. Soc. 136, 12489–12497. doi:10.1021/ja507382j
Barbet-Massin, E., Ricagno, S., Lewandowski, J. R., Giorgetti, S., Bellotti, V., Bolognesi, M., et al. (2010). Fibrillar vs Crystalline Full-Length β-2-Microglobulin Studied by High-Resolution Solid-State NMR Spectroscopy. J. Am. Chem. Soc. 132, 5556–5557. doi:10.1021/ja1002839
Bennett, A. E., Griffin, R. G., Ok, J. H., and Vega, S. (1992). Chemical Shift Correlation Spectroscopy in Rotating Solids: Radio Frequency‐driven Dipolar Recoupling and Longitudinal Exchange. J. Chem. Phys. 96, 8624–8627. doi:10.1063/1.462267
Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The Protein Data Bank. Nucleic Acids Res. 28, 235–242. doi:10.1093/nar/28.1.235
Brinkmann, A., and Levitt, M. H. (2001). Symmetry Principles in the Nuclear Magnetic Resonance of Spinning Solids: Heteronuclear Recoupling by Generalized Hartmann-Hahn Sequences. J. Chem. Phys. 115, 357–384. doi:10.1063/1.1377031
Burley, S. K., Bhikadiya, C., Bi, C., Bittrich, S., Chen, L., Crichlow, G. V., et al. (2021). RCSB Protein Data Bank: Powerful New Tools for Exploring 3D Structures of Biological Macromolecules for Basic and Applied Research and Education in Fundamental Biology, Biomedicine, Biotechnology, Bioengineering and Energy Sciences. Nucleic Acids Res. 49, D437. doi:10.1093/nar/gkaa1038
Busi, B., Yarava, J. R., Hofstetter, A., Salvi, N., Cala-De Paepe, D., Lewandowski, J. R., et al. (2018). Probing Protein Dynamics Using Multifield Variable Temperature NMR Relaxation and Molecular Dynamics Simulation. J. Phys. Chem. B 122, 9697–9702. doi:10.1021/acs.jpcb.8b08578
Case, D. A., Cheatham, T. E., Darden, T., Gohlke, H., Luo, R., Merz, K. M., et al. (2005). The Amber Biomolecular Simulation Programs. J. Comput. Chem. 26, 1668–1688. doi:10.1002/jcc.20290
Castellani, F., van Rossum, B., Diehl, A., Schubert, M., Rehbein, K., and Oschkinat, H. (2002). Structure of a Protein Determined by Solid-State Magic-Angle-Spinning NMR Spectroscopy. Nature 420, 98–102. doi:10.1038/nature01070
Chevelkov, V., Fink, U., and Reif, B. (2009). Accurate Determination of Order Parameters from 1H,15N Dipolar Couplings in MAS Solid-State NMR Experiments. J. Am. Chem. Soc. 131, 14018–14022. doi:10.1021/ja902649u
Chevelkov, V., Rehbein, K., Diehl, A., and Reif, B. (2006). Ultrahigh Resolution in Proton Solid-State NMR Spectroscopy at High Levels of Deuteration. Angew. Chem. Int. Ed. 45, 3878–3881. doi:10.1002/anie.200600328
Chevelkov, V., van Rossum, B. J., Castellani, F., Rehbein, K., Diehl, A., Hohwy, M., et al. (2003). H Detection in MAS Solid-State NMR Spectroscopy of Biomacromolecules Employing Pulsed Field Gradients for Residual Solvent Suppression. J. Am. Chem. Soc. 125, 7788–7789. doi:10.1021/ja029354b
Chevelkov, V., Xue, Y., Linser, R., Skrynnikov, N. R., and Reif, B. (2010). Comparison of Solid-State Dipolar Couplings and Solution Relaxation Data Provides Insight into Protein Backbone Dynamics. J. Am. Chem. Soc. 132, 5015–5017. doi:10.1021/ja100645k
Delaglio, F., Grzesiek, S., Vuister, G. W., Zhu, G., Pfeifer, J., and Bax, A. (1995). NMRPipe: A Multidimensional Spectral Processing System Based on UNIX Pipes. J. Biomol. NMR 6, 277–293. doi:10.1007/BF00197809
Doshi, U., and Hamelberg, D. (2009). Reoptimization of the AMBER Force Field Parameters for Peptide Bond (Omega) Torsions Using Accelerated Molecular Dynamics. J. Phys. Chem. B 113, 16590–16595. doi:10.1021/jp907388m
Fage, C. D., Kosol, S., Jenner, M., Öster, C., Gallo, A., Kaniusaite, M., et al. (2021). Communication Breakdown: Dissecting the COM Interfaces between the Subunits of Nonribosomal Peptide Synthetases. ACS Catal. 11, 10802–10813. doi:10.1021/acscatal.1c02113
Franks, W. T., Wylie, B. J., Stellfox, S. A., and Rienstra, C. M. (2006). Backbone Conformational Constraints in a Microcrystalline U-15N-Labeled Protein by 3D Dipolar-Shift Solid-State NMR Spectroscopy. J. Am. Chem. Soc. 128, 3154–3155. doi:10.1021/ja058292x
Franks, W. T., Wylie, B. J., Stellfox, S. A., and Rienstra, C. M. (2006). Backbone Conformational Constraints in a Microcrystalline U-15N-Labeled Protein by 3D Dipolar-Shift Solid-State NMR Spectroscopy. J. Am. Chem. Soc. 128, 3154–3155. doi:10.1021/ja058292x
Franks, W. T., Zhou, D. H., Wylie, B. J., Money, B. G., Graesser, D. T., Frericks, H. L., et al. (2005). Magic-Angle Spinning Solid-State NMR Spectroscopy of the β1 Immunoglobulin Binding Domain of Protein G (GB1): 15N and 13C Chemical Shift Assignments and Conformational Analysis. J. Am. Chem. Soc. 127, 12291–12305. doi:10.1021/ja044497e
Frauenfelder, H., Sligar, S. G., and Wolynes, P. G. (1991). The Energy Landscapes and Motions of Proteins. Science 254, 1598–1603. doi:10.1126/science.1749933
Fung, B. M., Khitrin, A. K., and Ermolaev, K. (2000). An Improved Broadband Decoupling Sequence for Liquid Crystals and Solids. J. Magn. Reson. 142, 97–101. doi:10.1006/jmre.1999.1896
Gullion, T., and Schaefer, J. (1989). Detection of Weak Heteronuclear Dipolar Coupling by Rotational-Echo Double-Resonance Nuclear Magnetic Resonance. Adv. Magn. Opt. Reson. 13, 57–83. doi:10.1016/B978-0-12-025513-9.50009-4
Gullion, T., and Schaefer, J. (1989). Rotational-echo Double-Resonance NMR. J. Magn. Reson. (1969) 81, 196–200. doi:10.1016/0022-2364(89)90280-1
Herbst, C., Bellstedt, P., Görlach, M., and Ramachandran, R. (2015). MAS Solid State NMR of Proteins: Simultaneous 15N-13CA and 15N-13CO Dipolar Recoupling via Low-Power Symmetry-Based RF Pulse Schemes. J. Biomol. NMR 62, 7–15. doi:10.1007/s10858-015-9910-2
Herbst, C., Herbst, J., Leppert, J., Ohlenschläger, O., Görlach, M., and Ramachandran, R. (2011). Chemical Shift Correlation at High MAS Frequencies Employing Low-Power Symmetry-Based Mixing Schemes. J. Biomol. NMR 50, 277–284. doi:10.1007/s10858-011-9516-2
Hester, R. K., Ackerman, J. L., Neff, B. L., and Waugh, J. S. (1976). Separated Local Field Spectra in NMR: Determination of Structure of Solids. Phys. Rev. Lett. 36, 1081–1083. doi:10.1103/physrevlett.36.1081
Hiller, M., Krabben, L., Vinothkumar, K. R., Castellani, F., van Rossum, B.-J., Kühlbrandt, W., et al. (2005). Solid-state Magic-Angle Spinning NMR of Outer-Membrane Protein G from Escherichia coli. Chembiochem 6, 1679–1684. doi:10.1002/cbic.200500132
Hohwy, M., Jakobsen, H. J., Edén, M., Levitt, M. H., and Nielsen, N. C. (1998). Broadband Dipolar Recoupling in the Nuclear Magnetic Resonance of Rotating Solids: A Compensated C7 Pulse Sequence. J. Chem. Phys. 108, 2686–2694. doi:10.1063/1.475661
Hohwy, M., Jaroniec, C. P., Reif, B., Rienstra, C. M., and Griffin, R. G. (2000). Local Structure and Relaxation in Solid-State NMR: Accurate Measurement of Amide N−H Bond Lengths and H−N−H Bond Angles. J. Am. Chem. Soc. 122, 3218–3219. doi:10.1021/ja9913737
Hoogen, Y. T., Treurniet, S. J., Roelen, H. C. P. F., Vroom, E., Marel, G. A., Boom, J. H., et al. (1988). Conformational Analysis of the Tetranucleotides m62A-m62A-U-m62A (m62A = N6-Dimethyladenosine) and U-m62A-U-m62A and of the Hybrid dA-r(U-A). A One- and Two-Dimensional NMR Study. Eur. J. Biochem. 171, 155–162. doi:10.1111/j.1432-1033.1988.tb13771.x
Hou, G., Byeon, I.-J. L., Ahn, J., Gronenborn, A. M., and Polenova, T. (2011). 1H-13C/1H-15N Heteronuclear Dipolar Recoupling by R-Symmetry Sequences under Fast Magic Angle Spinning for Dynamics Analysis of Biological and Organic Solids. J. Am. Chem. Soc. 133, 18646–18655. doi:10.1021/ja203771a
Izadi, S., Anandakrishnan, R., and Onufriev, A. V. (2014). Building Water Models: A Different Approach. J. Phys. Chem. Lett. 5, 3863–3871. doi:10.1021/jz501780a
Jenner, M., Kosol, S., Griffiths, D., Prasongpholchai, P., Manzi, L., Barrow, A. S., et al. (2018). Mechanism of Intersubunit Ketosynthase-Dehydratase Interaction in Polyketide Synthases. Nat. Chem. Biol. 14, 270–275. doi:10.1038/nchembio.2549
Jumper, J., Evans, R., Pritzel, A., Green, T., Figurnov, M., Ronneberger, O., et al. (2021). Highly Accurate Protein Structure Prediction with AlphaFold. Nature 596, 583–589. doi:10.1038/s41586-021-03819-2
Knight, M. J., Pell, A. J., Bertini, I., Felli, I. C., Gonnelli, L., Pierattelli, R., et al. (2012). Structure and Backbone Dynamics of a Microcrystalline Metalloprotein by Solid-State NMR. Proc. Natl. Acad. Sci. 109, 11095–11100. doi:10.1073/pnas.1204515109
Koshland, D. E. (1958). Application of a Theory of Enzyme Specificity to Protein Synthesis. Proc. Natl. Acad. Sci. 44, 98–104. doi:10.1073/pnas.44.2.98
Kosol, S., Gallo, A., Griffiths, D., Valentic, T. R., Masschelein, J., Jenner, M., et al. (2019). Structural Basis for Chain Release from the Enacyloxin Polyketide Synthase. Nat. Chem. 11, 913–923. doi:10.1038/s41557-019-0335-5
Kovermann, M., Rogne, P., and Wolf-Watz, M. (2016). Protein Dynamics and Function from Solution State NMR Spectroscopy. Quart. Rev. Biophys. 49, e6. doi:10.1017/s0033583516000019
Lamley, J. M., Iuga, D., Öster, C., Sass, H.-J., Rogowski, M., Oss, A., et al. (2014). Solid-state NMR of a Protein in a Precipitated Complex with a Full-Length Antibody. J. Am. Chem. Soc. 136, 16800–16806. doi:10.1021/ja5069992
Lamley, J. M., Lougher, M. J., Sass, H. J., Rogowski, M., Grzesiek, S., and Lewandowski, J. R. (2015). Unraveling the Complexity of Protein Backbone Dynamics with Combined 13C and 15N Solid-State NMR Relaxation Measurements. Phys. Chem. Chem. Phys. 17, 21997–22008. doi:10.1039/c5cp03484a
Lamley, J. M., Öster, C., Stevens, R. A., and Lewandowski, J. R. (2015). Intermolecular Interactions and Protein Dynamics by Solid‐State NMR Spectroscopy. Angew. Chem. Int. Ed. 54, 15374–15378. doi:10.1002/anie.201509168
Levitt, M. H. (1982). Symmetrical Composite Pulse Sequences for NMR Population Inversion. II. Compensation of Resonance Offset. J. Magn. Reson. (1969) 50, 95–110. doi:10.1016/0022-2364(82)90035-x
Levitt, M. H. (2002). Symmetry-Based Pulse Sequences in Magic-Angle Spinning Solid-State NMR. Encyclopedia Magn. Reson. 9, 165–196.
Lewandowski, J. R. (2013). Advances in Solid-State Relaxation Methodology for Probing Site-specific Protein Dynamics. Acc. Chem. Res. 46, 2018–2027. doi:10.1021/ar300334g
Lewandowski, J. R., Halse, M. E., Blackledge, M., and Emsley, L. (2015). Direct Observation of Hierarchical Protein Dynamics. Science 348, 578–581. doi:10.1126/science.aaa6111
Lewandowski, J. R., Sass, H. J., Grzesiek, S., Blackledge, M., and Emsley, L. (2011). Site-specific Measurement of Slow Motions in Proteins. J. Am. Chem. Soc. 133, 16762–16765. doi:10.1021/ja206815h
Liang, L., Ji, Y., Zhao, Z., Quinn, C. M., Han, X., Bao, X., et al. (2021). Accurate Heteronuclear Distance Measurements at All Magic-Angle Spinning Frequencies in Solid-State NMR Spectroscopy. Chem. Sci. 12, 11554–11564. doi:10.1039/d1sc03194e
Linser, R., Chevelkov, V., Diehl, A., and Reif, B. (2007). Sensitivity Enhancement Using Paramagnetic Relaxation in MAS Solid-State NMR of Perdeuterated Proteins. J. Magn. Reson. 189, 209–216. doi:10.1016/j.jmr.2007.09.007
Maier, J. A., Martinez, C., Kasavajhala, K., Wickstrom, L., Hauser, K. E., and Simmerling, C. (2015). ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theor. Comput. 11, 3696–3713. doi:10.1021/acs.jctc.5b00255
Marion, D., Ikura, M., Tschudin, R., and Bax, A. (1989). Rapid Recording of 2D NMR Spectra without Phase Cycling. Application to the Study of Hydrogen Exchange in Proteins. J. Magn. Reson. (1969) 85, 393–399. doi:10.1016/0022-2364(89)90152-2
Nielsen, A. B., Bjerring, M., Nielsen, J. T., and Nielsen, N. C. (2009). Symmetry-based Dipolar Recoupling by Optimal Control: Band-Selective Experiments for Assignment of Solid-State NMR Spectra of Proteins. J. Chem. Phys. 131, 025101. doi:10.1063/1.3157737
Öster, C., Kosol, S., and Lewandowski, J. R. (2019). Quantifying Microsecond Exchange in Large Protein Complexes with Accelerated Relaxation Dispersion Experiments in the Solid State. Sci. Rep. 9, 11082. doi:10.1038/s41598-019-47507-8
Palmer, A. G. (2004). NMR Characterization of the Dynamics of Biomacromolecules. Chem. Rev. 104, 3623–3640. doi:10.1021/cr030413t
Paluch, P., Pawlak, T., Amoureux, J.-P., and Potrzebowski, M. J. (2013). Simple and Accurate Determination of X-H Distances under Ultra-fast MAS NMR. J. Magn. Reson. 233, 56–63. doi:10.1016/j.jmr.2013.05.005
Paluch, P., Pawlak, T., Ławniczak, K., Trébosc, J., Lafon, O., Amoureux, J.-P., et al. (2018). Simple and Robust Study of Backbone Dynamics of Crystalline Proteins Employing 1H-15N Dipolar Coupling Dispersion. J. Phys. Chem. B 122, 8146–8156. doi:10.1021/acs.jpcb.8b04557
Penzel, S., Oss, A., Org, M.-L., Samoson, A., Böckmann, A., Ernst, M., et al. (2019). Spinning Faster: Protein NMR at MAS Frequencies up to 126 kHz. J. Biomol. NMR 73, 19–29. doi:10.1007/s10858-018-0219-9
Prompers, J. J., and Brüschweiler, R. (2002). General Framework for Studying the Dynamics of Folded and Nonfolded Proteins by NMR Relaxation Spectroscopy and MD Simulation. J. Am. Chem. Soc. 124, 4522–4534. doi:10.1021/ja012750u
Retel, J. S., Nieuwkoop, A. J., Hiller, M., Higman, V. A., Barbet-Massin, E., Stanek, J., et al. (2017). Structure of Outer Membrane Protein G in Lipid Bilayers. Nat. Commun. 8, 2073. doi:10.1038/s41467-017-02228-2
Roe, D. R., and Cheatham, T. E. (2013). PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theor. Comput. 9, 3084–3095. doi:10.1021/ct400341p
Ryckaert, J.-P., Ciccotti, G., and Berendsen, H. J. C. (1977). Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of N-Alkanes. J. Comput. Phys. 23, 327–341. doi:10.1016/0021-9991(77)90098-5
Schanda, P., and Ernst, M. (2016). Studying Dynamics by Magic-Angle Spinning Solid-State NMR Spectroscopy: Principles and Applications to Biomolecules. Prog. Nucl. Magn. Reson. Spectrosc. 96, 1–46. doi:10.1016/j.pnmrs.2016.02.001
Schledorn, M., Malär, A. A., Torosyan, A., Penzel, S., Klose, D., Oss, A., et al. (2020). Protein NMR Spectroscopy at 150 kHz Magic‐Angle Spinning Continues to Improve Resolution and Mass Sensitivity. ChemBioChem 21, 2540–2548. doi:10.1002/cbic.202000341
Sekhar, A., and Kay, L. E. (2019). An NMR View of Protein Dynamics in Health and Disease. Annu. Rev. Biophys. 48, 297–319. doi:10.1146/annurev-biophys-052118-115647
Sternberg, U., Witter, R., Kuprov, I., Lamley, J. M., Oss, A., Lewandowski, J. R., et al. (2018). 1H Line Width Dependence on MAS Speed in Solid State NMR - Comparison of experiment and Simulation. J. Magn. Reson. 291, 32–39. doi:10.1016/j.jmr.2018.04.003
Thurber, K. R., and Tycko, R. (2009). Measurement of Sample Temperatures under Magic-Angle Spinning from the Chemical Shift and Spin-Lattice Relaxation Rate of 79Br in KBr Powder. J. Magn. Reson. 196, 84–87. doi:10.1016/j.jmr.2008.09.019
Tian, C., Kasavajhala, K., Belfon, K. A. A., Raguette, L., Huang, H., Migues, A. N., et al. (2019). ff19SB: Amino-acid-specific Protein Backbone Parameters Trained against Quantum Mechanics Energy Surfaces in Solution. J. Chem. Theor. Comput. 16, 528–552. doi:10.1021/acs.jctc.9b00591
Wang, J., Wolf, R. M., Caldwell, J. W., Kollman, P. A., and Case, D. A. (2004). Development and Testing of a General Amber Force Field. J. Comput. Chem. 25, 1157–1174. doi:10.1002/jcc.20035
Wishart, D. S., Bigam, C. G., Yao, J., Abildgaard, F., Dyson, H. J., Oldfield, E., et al. (1995). 1H, 13C and 15N Chemical Shift Referencing in Biomolecular NMR. J. Biomol. NMR 6, 135–140. doi:10.1007/bf00211777
Wylie, B. J., and Rienstra, C. M. (2008). Multidimensional Solid State NMR of Anisotropic Interactions in Peptides and Proteins. J. Chem. Phys. 128, 052207. doi:10.1063/1.2834735
Wylie, B. J., Sperling, L. J., Nieuwkoop, A. J., Franks, W. T., Oldfield, E., and Rienstra, C. M. (2011). Ultrahigh Resolution Protein Structures Using NMR Chemical Shift Tensors. Proc. Natl. Acad. Sci. U. S. A. 108, 16974–16979. doi:10.1073/pnas.1103728108
Zhou, D. H., Shea, J. J., Nieuwkoop, A. J., Franks, W. T., Wylie, B. J., Mullen, C., et al. (2007). Solid-state Protein-Structure Determination with Proton-Detected Triple-Resonance 3D Magic-Angle-Spinning NMR Spectroscopy. Angew. Chem. Int. Ed. Engl. 46, 8380–8383. doi:10.1002/anie.200702905
Zhou, D. H., and Rienstra, C. M. (2008). High-performance Solvent Suppression for Proton Detected Solid-State NMR. J. Magn. Reson. 192, 167–172. doi:10.1016/j.jmr.2008.01.012
Zhou, D. H., Shea, J. J., Nieuwkoop, A. J., Franks, W. T., Wylie, B. J., Mullen, C., et al. (2007). Solid-state Protein-Structure Determination with Proton-Detected Triple-Resonance 3D Magic-Angle-Spinning NMR Spectroscopy. Angew. Chem. Int. Ed. 46, 8380–8383. doi:10.1002/anie.200702905
Keywords: symmetry, dynamics, fast MAS NMR, recoupling, order parameter, proton detection
Citation: Franks WT, Tatman BP, Trenouth J and Lewandowski JR (2021) Dipolar Order Parameters in Large Systems With Fast Spinning. Front. Mol. Biosci. 8:791026. doi: 10.3389/fmolb.2021.791026
Received: 07 October 2021; Accepted: 05 November 2021;
Published: 09 December 2021.
Edited by:
Mark Pfuhl, King’s College London, United KingdomReviewed by:
Johanna Becker-Baldus, Goethe University Frankfurt, GermanySergey Dvinskikh, Royal Institute of Technology, Sweden
Copyright © 2021 Franks, Tatman, Trenouth and Lewandowski. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: W. Trent Franks, t.franks@warwick.ac.uk; Józef R. Lewandowski, j.r.lewandowski@warwick.ac.uk