Relaxation-induced dipolar exchange with recoupling (RIDER) distortions in CODEX experiments

Abstract Chemical shift anisotropy (CSA) and dipolar CODEX (Cenralband Only Detection of EXchange) experiments enable abundant quantitative information on the reorientation of the CSA and dipolar tensors to be obtained on millisecond–second timescales. At the same time, proper performance of the experiments and data analysis can often be a challenge since CODEX is prone to some interfering effects that may lead to incorrect interpretation of the experimental results. One of the most important such effects is RIDER (relaxation-induced dipolar exchange with recoupling). It appears due to the dipolar interaction of the observed X nuclei with some other nuclei, which causes an apparent decay in the mixing time dependence of the signal intensity reflecting not molecular motion, but spin flips of the adjacent nuclei. This may hamper obtaining correct values of the parameters of molecular mobility. In this contribution we consider in detail the reasons why the RIDER distortions remain even under decoupling conditions and propose measures to eliminate them. That is, we suggest (1) using an additional Z filter between the cross-polarization (CP) section and the CODEX recoupling blocks that suppresses the interfering anti-phase coherence responsible for the X -H RIDER and (2) recording only the cosine component of the CODEX signal since it is less prone to the RIDER distortions in comparison to the sine component. The experiments were conducted on rigid model substances as well as microcrystalline 2 H  /   15 N-enriched proteins (GB1 and SH3) with a partial back-exchange of labile protons. Standard CSA and dipolar CODEX experiments reveal a fast-decaying component in the mixing time dependence of 15 N nuclei in proteins, which can be misinterpreted as a slow overall protein rocking motion. However, the RIDER-free experimental setup provides flat mixing time dependences, meaning that the studied proteins do not undergo global motions on the millisecond timescale.


Introduction
CODEX (Cenralband Only Detection of EXchange) (deAzevedo et al., 1999(deAzevedo et al., , 2000Luz et al., 2002;Reichert and Krushelnitsky, 2018) is a powerful nuclear magnetic resonance (NMR) tool for studying molecular dynamics in millisecond to second timescales under magic angle spinning (MAS). It is based on the stimulated echo principle; the simplified pulse sequence is shown in Fig. 1. Depending on the phases of the rf pulses and receiver, one may record a signal, which is proportional to sin( 1 ) · sin( 2 ) (SIN component) or cos( 1 ) · cos( 2 ) (COS component), where 1 and 2 are the phases accumulated by the magnetization vector during the precession under recoupling conditions in the dephasing and rephasing periods, respectively. The sum of the two signals (COS and SIN components) is proportional to cos( 1 − 2 ). The classical CODEX experiment was designed for observing the reorientation of the chemical shift anisotropy (CSA) tensor: the REDOR-like (Gullion and Schaefer, 1989) train of rotor-synchronized recoupling π pulses applied with a half-rotor period spacing on the X nuclei reintroduces the CSA interaction and, thus, the phases 1 and 2 are determined by the precession under the influence of the CSA interaction during the de(re)phasing periods. Potentially interfering dipolar interactions with protons are supposed to be averaged out by proton decoupling during the de(re)phasing periods. However, CODEX can be easily modified for observing motionally modulated Black vertical bars denote π/2 pulses; ϕ CP , ϕ 1 , ϕ 2 and ϕ Rec are the phases of the X-channel cross-polarization (CP) pulse, two π/2 pulses and the receiver, respectively. The COS component is recorded when the phase differences are ϕ CP − ϕ 1 = π/2 and ϕ 2 −ϕ Rec = π/2; the SIN component corresponds to ϕ CP = ϕ 1 and ϕ 2 = ϕ Rec . dipolar interaction or isotropic chemical shift (i.e. chemical exchange). This can be achieved by a corresponding modification of the recoupling pulses (Krushelnitsky et al., 2013;Reichert and Krushelnitsky, 2018).
In the CODEX experiment, one can measure the signal intensity as a function of both mixing time and the length of the de(re)phasing periods N T R (T R is the MAS period and N is the number of rotor cycles in the de(re)phasing periods), which provides the information on both timescale and geometry of a molecular motion (Luz et al., 2002). Thus, the CODEX experiment enables more abundant quantitative information on molecular dynamics to be obtained in comparison to standard NMR relaxation studies. At the same time, CODEX is prone to some interfering effects that may distort the information on molecular dynamics and that should be taken into account in the data analysis. Two most important effects are the proton-driven spin diffusion between X nuclei and RIDER (relaxation-induced dipolar exchange with recoupling) (Saalwächter and Schmidt-Rohr, 2000). Spin diffusion reveals itself as a signal decay in the mixing time dependence, which can be in some cases erroneously attributed to a molecular motion process. Suppressing the spin diffusion by proton decoupling during the mixing time is in principle possible but is rather difficult and not always reliable and effective (Reichert and Krushelnitsky, 2018). The most robust way of removing the undesirable spin-diffusion effect is a spin dilution, e.g. using natural abundance 13 C or perdeuterated samples.
RIDER also leads to an additional decay in the mixing time dependence. Dipolar interaction of X nuclei (S) with either protons or some other magnetic nuclei present in a sample (I ) adds two terms of the precessing X-nuclei magnetization -in-phase S x cos(ωt) and anti-phase 2S y I z sin(ωt). The last term is the origin of RIDER, which can be simplistically explained as follows: flips of I z during the mix-ing time change the sign of the inter-nuclear dipolar interaction (for 1/2 nuclei) and thus change the sign of the dipolar contribution to the precession frequency. This in turn leads to incomplete rephasing of the S magnetization at the end of the rephasing period and thus to decrease in the signal. Therefore, the characteristic time of the decay in the mixing time dependence due to RIDER is determined by the timescale of I z flips, that is, T 1 relaxation of nuclei I . In addition, if the homonuclear dipolar interaction between I spins is significant, spin diffusion (flip-flops) also contributes to the timescale of RIDER, which can be much shorter than T 1 of I spins. The standard way of suppressing RIDER in the CODEX experiments is heteronuclear I -S decoupling during the de(re)phasing periods. For some I nuclei with a large quadrupolar moment, e.g. 14 N, decoupling is not effective, and in this case, the only way of removing the undesirable RIDER influence is isotopic editing of a sample.
Our interest in the methodological problems of the CODEX experiments was stimulated by the study of slow motions in solid proteins. Recently, it was shown by means of R 1ρ relaxometry that proteins in a solid state undergo slow overall rocking motion (Ma et al., 2015;Lamley et al., 2015;Kurauskas et al., 2017;Krushelnitsky et al., 2018). The timescale of this motion is tens of microseconds, which is the limit of the time window accessible with R 1ρ relaxation experiments. What happens on the (sub)millisecond timescale up to now remained unclear, and the CODEX experiments could answer the question of whether the rocking motion extends to longer correlation times or not.
We have thus conducted CSA and dipolar CODEX experiments on 15 N nuclei in 15 N, 2 H-enriched microcrystalline proteins (SH3 and GB1) with a partial back-exchange of labile protons. These experiments were conducted with a sitespecific resolution in a 2D 1 H-15 N correlation spectrum using indirect proton detection of a signal (Chevelkov et al., 2006;Krushelnitsky et al., 2009). Surprisingly, all peaks in 2D spectra without exception reveal decays in the mixing time dependences as shown in Fig. 2. The amplitude of the decay and the apparent correlation time of the fast component (around 20 ms) are the same for all residues. This component cannot be due to the proton-driven spin diffusion since the timescale of the spin diffusion between 15 N nuclei even in fully protonated proteins is much longer (Krushelnitsky et al., 2006). In the CSA CODEX, this could be the RIDER effect arising due to the dipolar interaction between 15 N and abundant 2 H nuclei. On the other hand, in the dipolar CODEX experiment, we observe very similar shapes of the mixing time dependences with very similar parameters of the fast component. This was observed for both SH3 and GB1 microcrystalline proteins. In the dipolar CODEX experiment, the recoupling π pulses are applied on the proton channel and, thus, the 15 N-2 H dipolar interaction should be effectively averaged out by MAS. From this one could conclude that the observed fast component of the mixing time dependences is not an artefact and does report on a real over- all protein motion on the millisecond timescale. This would mean that the rocking motion of a protein in a crystal has a very wide correlation time distribution, from microseconds to milliseconds.
However, it turned out that the fast-decaying component in the mixing time dependences is actually a highly non-trivial artefact. Its nature proved to be more complicated than the simple 15 N-2 H RIDER effect. Below we explain the details of the effects responsible for the appearance of this component and suggest some measures for correctly conducting CODEX experiments and avoiding misinterpretations of CODEX data in proteins as well as other samples with complex isotopic composition.

Theory
Here we consider the time evolution of spin coherences in the CSA CODEX experiments using product operator formalism. It is well known that after I → S cross-polarization (CP), both in-phase S x and anti-phase −2S y I z terms appear; see e.g. Schmidt-Rohr and Spiess (1994). The anti-phase term is usually neglected since in standard CP/MAS experiments it is suppressed by the heteronuclear proton decoupling during the FID acquisition. In the CSA CODEX, it is supposed to be suppressed by the proton decoupling during the de(re)phasing periods as well. However, in the CODEX pulse sequence this suppression is much less effective. The train of the X-channel recoupling π pulses applied during the de(re)phasing periods restores not only CSA, but also dipolar X-1 H interaction. Hence, the proton decoupling affects not the residual (after MAS) dipolar interaction, but the restored (recoupled) value of this interaction. For this reason, the small but appreciable dipolar X-1 H interaction survives during the de(re)phasing periods, which will be demonstrated experimentally below, and we have to take it into account in our analysis.
Let us consider the time evolution of the in-phase and antiphase terms in the CSA CODEX experiment under the simultaneous influence of the CSA and (not completely suppressed) dipolar interactions during the de(re)phasing periods. The phases acquired during the dephasing period under the influence of the CSA and dipolar interactions we denote as CSA and D , respectively. We assume for simplicity that CSA remains the same for both the dephasing and rephasing periods, but D can change due to RIDER. Thus, for the rephasing period, the acquired phase will be denoted as D + D . In-phase term, dephasing period: − 2S x I z sin ( CSA ) sin ( D ) + S y sin ( CSA ) cos ( D ) + 2S y I z cos ( CSA ) sin ( D ) . (1) The first two terms are picked up in the COS component and two second terms in the SIN component of the CODEX signal. At the end of the rephasing period, we have the COS component, for spin I = 1/2, since D can be either 0 or −2 D . This means that for the in-phase term, the effect of the incomplete suppression of the dipolar X-1 H interaction is almost negligible: it leads only to a small decrease in the signal, proportional to cos 2 ( D ).
Let us now consider the time evolution of the anti-phase term. At the end of the dephasing period we have Analogously to Eq. (1), the first two terms in Eq. (5) correspond to the COS component and the second two terms to the SIN component. After the rephasing period, the COS component reads as and the SIN component is Again, in Eqs. (7) and (8) only the observable terms are left that correspond to the COS (Eq. 7) and SIN (Eq. 8) components. It is seen from these equations that for the anti-phase term, the RIDER effect is not negligible, and the inequality analogous to Eq. (4) cannot be written if D is small but appreciable. But how can the RIDER effect arising from the anti-phase term be recognized in the analysis of experimental data? This is relatively simple: one may compare the shapes of the mixing time dependence of the COS and SIN components. If these curves, namely the ratio S ∞ /S 0 (S 0 and S ∞ are the signal amplitudes at very short and very long mixing times, respectively), are not similar, then RIDER is relevant. In general, the ratio S ∞ /S 0 for the COS and SIN components should be exactly the same if only molecular motions and/or spin diffusion are present in a sample. This can be proved as follows. Let us denote the phases acquired during the dephasing and rephasing periods as and + , respectively. At short mixing times, = 0; then the ratio S ∞ /S 0 for different cases would be as follows.
Classical CODEX (COS+SIN components): COS component: SIN component: Next, we have to recall that ij = − j i (i and j are the numbers of the exchanging sites), and since it is always assumed that we are dealing with dynamic equilibrium (i.e. the populations of the exchanging sites are constant in time), then obviously sin( ) = 0. Thus, in all cases S ∞ S 0 = cos( ) , that is, the shapes of the COS and SIN components should be the same, although the absolute amplitudes in the general case are of course different. Now, let us estimate the ratio S ∞ /S 0 for the COS and SIN components described in Eqs. (7) and (8) taking into account Eq. (5) and the equation sin( D ) = sin( D + D ) = 0 (note that this is valid only for I = 1/2). COS component: SIN component: Hence, it is clear that the RIDER effect leads to different shapes of the mixing time dependence of the COS and SIN components. Note that if D is not small, the ratio S ∞ /S 0 would be different for the COS and SIN components, also for the in-phase term; see Eqs. (2) and (3). From Eqs.
(2), (3), (7) and (8) it can also be deduced that the SIN component is about twice as prone to the RIDER distortions as the COS component. This follows from the comparison of the amplitudes of different coherences. The amplitudes of the COS component of the in-phase and anti-phase terms (see Eqs. 2 and 7) are proportional to cos 2 ( CSA ) · cos 2 ( D ) and cos 2 ( CSA )·cos( D )·sin( D ), respectively (here we assume CSA to be not large). Hence, the ratio of the anti-phase-term amplitude to the in-phase-term amplitude is proportional to tan( D ) or even smaller if the second terms in the parentheses in Eqs. (2) and (7) are taken into account. The same ratio for the SIN component (see Eqs. 3 and 8) is proportional to 2 · tan( D ). Thus, the contribution of the anti-phase coherence to the total signal is larger for the SIN component.
The analysis presented above is valid only for an isolated I -S spin pair. For multinuclear spin systems, the description would be much more complicated since many types of multiple coherences with a complex network of homo-and hetero-nuclear dipolar interactions should be taken into account. Quantifying this is outside the scope of our work; still, we believe that on a qualitative level, two most important points remain valid: first, the anti-phase term appearing after the CP pulses may cause RIDER distortions of the mixing time dependences, and second, the RIDER effect can be recognized from the comparison of the shapes of the COS and SIN components. This will be proven experimentally below.

Samples
In our work we used four different samples. Model substances: 15 N-enriched BOC (N-(tert-Butoxycarbonyl)) glycine and 15 N-enriched glycine, which were purchased from Sigma-Aldrich. Proteins: small GB1 and SH3 proteins in a form of microcrystals, 15 N and 2 H enriched with a partial back-exchange of labile protons. The GB1 sample was purchased from Giotto Biotech (Florence, Italy); the SH3 sample was prepared in Bernd Reif's lab (FMP, Berlin). These are the same samples that were used in our recent work on R 1ρ relaxometry (Krushelnitsky et al., 2018). Both protein samples were prepared according to the protocol ensuring 20 % of the back-exchange of labile protons. However, we believe that in reality this percentage is somewhat different: in GB1 it is higher, which is indicated by a stronger signal and faster proton-driven spin diffusion between 15 N nuclei (see Figs. 12 and 13 below). The quantitative estimation of this difference is as yet rather difficult and uncertain. Since the GB1 sample provides a better signal-to-noise ratio, most of the experiments were conducted with this sample.

NMR experiments
The experiments were performed on a Bruker AVANCE II NMR spectrometer (600 MHz) with a 3.2 mm MAS probe. In the CODEX experiments with the protein samples, the integral intensity of the entire signal was determined without site-selective specification (except for the data shown in Fig. 2). One-dimensional double CP ( 1 H→ 15 N→ 1 H) proton-detected spectra for SH3 and GB1 proteins were shown in Krushelnitsky et al. (2018). For the BOC glycine and glycine samples direct 15 N or 13 C signal detection was employed, and for the protein samples we used indirect 1 H signal detection of the 15 N CODEX mixing time dependences. This was implemented by using a back CP section ( 15 N→ 1 H) at the end of the pulse sequence, according to the approach described earlier (Chevelkov et al., 2006;Krushelnitsky et al., 2009). We have checked -the direct 15 N and indirect 1 H signal detections in the protein samples provide the same shape of the CODEX mixing time dependences; in the latter case the signal-to-noise ratio was however better.
To exclude the effect of spin-lattice relaxation during the mixing time, each CODEX mixing time dependence was T 1normalized. For that, for each mixing time dependence two experiments were performed: measuring the mixing time dependence itself and measuring a T 1 -relaxation curve within the same time range. After that, the mixing time dependence was divided by a T 1 -relaxation curve. This is a routine procedure described earlier (deAzevedo et al., 1999(deAzevedo et al., , 2000Reichert et al., 2001;Reichert and Krushelnitsky, 2018). Below are shown only the T 1 -normalized mixing time dependences for all CODEX experiments in protein samples. For BOC glycine, the T 1 normalization was not performed since 15 N T 1 in this sample was extremely long (800-900 s).
The pulse sequences of the CSA and dipolar CODEX are shown in Figs. 3 and 4. To measure the mixing time dependence, τ m was variable and τ r was fixed at 1 ms; to measure the T 1 -relaxation curve, τ m was fixed at 1 ms and τ r was variable. The phase cycle for both the COS and SIN components consists of 64 steps: 2× spin-temperature inversion (ensuring that the signal decays to zero; Torchia, 1978) for T 1 relaxation, 2× spin-temperature inversion for mixing time, 4× CYCLOPS for the π/2 pulse after mixing time, and 4× CYCLOPS for the π/2 pulse after τ r delay (Reichert et al., 2001). Typical values for the π/2 pulse for 1 H and 15 N channels were 1.4-1.8 and 6.0-6.5 µs, respectively.
The experimental error in estimation of the signal amplitude was less than 1 % for BOC glycine, 1 %-2 % for GB1, 2 %-4 % for SH3 and 5 %-10 % for natural abundance 13 C in glycine. On top of the signal noise, a certain contribution to the experimental error in the mixing time dependences comes from the instability of the MAS controller; this was however significant only for BOC glycine. The final error of the mixing time dependences for this sample was around 1 %-2 %. For better visual distinguishing between the mixing time dependences shown in Figs. 5, 7 and 8, the adjacent Solid and open bars denote π/2 and π pulses, respectively. The mixing time τ m is an integer multiple of the MAS period, which is achieved by MAS rotor triggering before and at the end of the mixing time (see details in Reichert and Krushelnitsky, 2018). Rotor synchronization during the τ r delay is not necessary. Waltz decoupling in the indirect detection sequence aims to suppress only J coupling between X and 1 H nuclei; therefore, it has a low amplitude (a few hundred Hz). An additional initial Z filter and 2 H decoupling (see below) are not shown. Phase cycle: ϕ 1 = x; ϕ 2 = y; ϕ 3 = x; ϕ 4 = x; ϕ 5 = (y, −y) (COS component); averaging over a five-point filter was applied to the experimental curves in these figures, which significantly decreased the noise spread of the points without a change in the overall shape of the curves.   Fig. 3. π pulses applied on the X channel are set in the middle of the de(re)phasing periods; therefore, the duration of these periods should be an even multiple of the MAS period. The phase cycle is identical to that shown in Fig. 3. The phases of the π pulses applied during the de(re)phasing periods on the 1 H channel have no critical significance.
diffusion between 15 N nuclei in BOC glycine is very slow (Krushelnitsky et al., 2006). First, we demonstrate that the anti-phase term discussed above does really cause RIDER distortions in the CSA CODEX mixing time dependence. The anti-phase term appears in the course of CP; thus, its contribution to the total CODEX signal should depend on the CP contact time. The CSA CODEX mixing time dependences at various CP times are shown in Fig. 5. These data fully confirm the qualitative theoretical analysis presented above. One may see that the amplitude of the RIDER decay depends on the CP time, that the COS and SIN components of the mixing time dependences are different and that the SIN component is more prone to the RIDER distortions than the COS component. It is interesting to mention that the mixing time dependences shown in Fig. 5 reveal the decays on two different timescales: a few milliseconds and a few hundred milliseconds. Such a two-component shape of the decays reflects two different mechanisms that cause proton spin flips mentioned in the introduction above: spin diffusion (flip-flops) and spin-lattice relaxation.
If the heteronuclear proton decoupling during the de(re)phasing periods was effective enough, than the RIDER effect caused by the anti-phase coherence could have been of course avoided. However, this is not always possible for practical reasons because of the hardware limitations for the power of the long proton pulses. We tried to optimize the proton decoupling by the maximum signal at short mixing times. Different decoupling schemes were checked (TPPI, WALZ, SPINAL) at maximum proton power around 100-130 kHz, but none of them provided much better efficiency than simple CW decoupling (which is not the case for 1 H decoupling during FID, where CW is not the best choice). Therefore, in the experiments shown here we used CW 1 H decoupling during the de(re)phasing periods in the CSA CODEX measurements. We do not claim that CW decoupling is the best option for this purpose. It is quite possible that some other decoupling schemes specifically designed for the case of the recoupling X pulses can perform better. However, even having such a decoupling scheme at hand, one should carefully optimize it for different MAS rates and 1 H field strengths. We suggest here another, more simple and robust way of suppressing the undesired RIDER effect.
The anti-phase term can simply be suppressed by an additional Z filter between the CP pulses and the dephasing period, as illustrated in Fig. 6. The delay of this Z filter should be short compared to 15 N T 1 and long compared to T 2 . Thus, after such a Z filter one would have only an in-phase component. Figure 7 shows the mixing time dependences of the COS and SIN components at different delays of the Z filter. It is clearly seen that the Z filter fully removes the contribution of the anti-phase coherence.
Still, it is seen that even at long delays of the Z filter, the mixing time dependences are not completely flat, as they should be. The observed distortions are obviously the RIDER effect of the in-phase coherence. The Z filter eliminates the anti-phase coherence (Eqs. 7 and 8), but it does not improve the efficiency of the proton decoupling during the de(re)phasing periods, and thus the phase D remains non-zero. If the second terms in the parentheses in Eqs. (2) and (3) are not negligibly small in comparison to the first terms, then the RIDER is present also in the in-phase coherence and the Z filter obviously cannot remove it. Figure 8 presents the COS and SIN components of the mixing time dependences at different durations of the de(re)phasing periods measured with the additional Z filter. This is clearly seen: Figure 6. Initial part of the CODEX pulse sequence (Figs. 1 and 3) with the additional Z filter between the CP section and the dephasing period. ϕ CP − ϕ Z1 = ±π/2, ϕ Z2 = −ϕ Z1 . All the dependences were normalized to the same initial amplitude. MAS 20 kHz, 108 kHz CW 1 H decoupling during the de(re)phasing periods, NT R 2 ms, CP contact time 3 ms. the longer N T R is, the larger the RIDER distortions. This is reasonable since D is proportional to N T R . Note that the COS component is less prone to distortions, not only in the case of the "anti-phase", but also in the case of the "in-phase" RIDER. We are not able at the moment to explain the unusual bell-shaped form of the mixing time dependences. It is likely that the network of multi-nuclear dipolar interactions should be taken into account, and thus the explanation will not be simple. We also cannot exclude the possibility that transient NOE effects may play a certain role.
However, in any case this is the unwanted distortion, and regardless of the exact nature of this distortion, it should be maximally suppressed in real experiments. For this, the efficiency of the 1 H decoupling during the de(re)phasing periods must be optimized as far as possible. As mentioned above, the standard heteronuclear decoupling schemes used for FID detection do not help much for the case of de(re)phasing periods. This is illustrated by the example of the SPINAL sequence; see Fig. 8. Still, one may minimize the "in-phase" RIDER effect by keeping N T R as short as possible and by recording only the COS component of the mixing time dependence. Anyway, the "in-phase" RIDER is much smaller than the "anti-phase" one, and in most real experiments it can be safely neglected, as we will see below by the example of the protein samples.
At the end of this section, we demonstrate the application of the additional Z filter to the natural abundance 13 C CSA CODEX experiment performed on carbonyl carbons in 15 Nenriched glycine ( 15 N enrichment is necessary to avoid the 13 C-14 N RIDER effect). We see the same effect as in the case of 15 N CSA CODEX (Fig. 9). The dependences measured with the Z filter (red points in Fig. 9) are not completely flat; however, this is hardly due to the "in-phase" RIDER since the shapes of the SIN and COS components are very similar (in the case of RIDER they should be different) and the time constant of the decay (about 50-60 s) is obviously too long compared to the proton T 1 (a few seconds). We suspect that this decay is a manifestation of the proton-driven spin diffusion between natural abundance 13 C nuclei. Its time constant is roughly of the same order of magnitude as spin diffusion times between natural abundance 13 C nuclei measured in other organic solids; see e.g. Reichert et al. (1998). Spin diffusion however has no direct relevance to the topic of this work, and we did not analyse this in detail.
In summary, the theoretical and experimental results presented above show that the proton decoupling under the influence of the recoupling π pulses in the CSA CODEX is not fully efficient. This leads to the evolution of both in-phase and anti-phase coherences during the de(re)phasing periods under the influence of the residual 15 N( 13 C)-1 H dipolar coupling, that is, to the RIDER effect. The dominant contribution to the RIDER distortions of a mixing time dependence arises from the anti-phase term. This contribution can be suppressed by the additional Z filter between CP and the dephasing period. The RIDER distortion of the in-phase term  (Fig. 6). All decays were normalized to the same initial amplitude, the real ratio between the amplitudes of SIN and COS components is 0.7 for both experiments; 80 kHz CW 1 H decoupling and 35 kHz CW 15 N decoupling during the de(re)phasing periods were applied. MAS 22 kHz, NT R 2 ms, Zfilter delay 20 ms, CP contact time 3 ms. is smaller but still appreciable at long de(re)phasing periods. This interfering effect cannot be suppressed completely, but it can be significantly minimized if only the COS component of the mixing time dependence is measured and analysed since this component is less prone to RIDER in comparison to the SIN component.

Protein samples
In the protein samples we have three types of nuclei that we need to take into account -15 N, 1 H and 2 H. The direct and inverse 1 H-15 N CP sections employed in the CODEX pulse sequence ensure that in the experiment we observe only those nitrogens that have protons attached, and all 15 N's coupled to 2 H in the protein backbone remain invisible. Still, the interactions between protonated 15 N's and many remote 2 H's can be sufficient to induce RIDER-type distortions in the CODEX experiment. To demonstrate the hierarchy of the inter-nuclear interactions in our samples, we measured 15 N Hahn-echo decays (Fig. 10) and the initial signal S 0 (the signal at short mixing time) in the CSA and dipolar CODEX experiments as a function of N T R (Fig. 11) for various combinations of 1 H and 2 H decoupling schemes. Note that the S 0 vs. N T R dependence is in fact the analogue of the Hahn-echo experiment, the only difference is that either CSA or dipolar interaction is reintroduced by means of recoupling pulses during the transverse relaxation.
The conclusions that can be deduced from these data are as follows. First, despite the proton dilution, the 15 N-1 H dipolar line broadening at the MAS frequency 20 kHz remains quite  appreciable and strong proton decoupling is necessary to suppress the 15 N-1 H dipolar interaction. The comparison of the S 0 vs. NT R dependences of dipolar CODEX in fully protonated BOC glycine and the deuterated protein shows that the proton dilution reduces of course the inter-proton interaction (flip-flops) and thus, the rate of the 15 N decay: slower 1 H flips ensure slower 15 N-1 H coupling modulation and, hence, better refocus the signal after the end of the rephasing period. Still, the rate of the proton flip-flops in the protein sample remains in the millisecond range. This is an important point which will be discussed below. This result corresponds well to the proton line width estimations made by Reif and coworkers (Chevelkov et al., 2006). Second, it is clearly seen that the 130 kHz CW decoupling performs much worse in comparison to the SPINAL scheme (Fig. 10). As mentioned above, under the influence of the 15 N recoupling π pulses during the de(re)phasing periods, SPINAL does not provide significant advantage in comparison to CW. This confirms our previous statement that the proton decoupling efficiency under the influence of the Xchannel recoupling pulses is much different in comparison to FID detection.
Third, the 15 N-2 H interaction is non-negligible. 2 H decoupling does not lead to longer the Hahn-echo decays since it is effectively (but not completely; see below) reduced by MAS even without decoupling. However, the reintroduction of the 15 N-2 H dipolar coupling by the REDOR pulse train applied on deuterons appreciably shortens the decays, see Fig. 10. In the CSA CODEX experiment, the 15 N-2 H interaction is initially reintroduced by means of REDOR pulse train applied on 15 N's. In this case, the 2 H decoupling has the effect and makes the decay slower (Fig. 11). Now the recipe for a methodologically correct CSA CODEX experiment is evident. In deuterated proteins, there are two simultaneous RIDER effects arising from the 15 N-1 H and 15 N-2 H dipolar interactions, and one has to take care of both of them. Coincidentally, both RIDERs have similar, although not exactly the same, time constants. The time constant for the proton flip-flops can be estimated directly from the Hahn decay, which gives the value of about 10 ms (Fig. 10). As for the 2 H T 1 , it has a value of 25 ms for aliphatic deuterons in the SH3 protein sample, which was measured by a simple inversion-recovery method. Both these values are quite close to the time constant of the short component of the CODEX mixing time dependences observed in proteins (Fig. 2).
The 15 N-1 H and 15 N-2 H RIDER effects can be suppressed by the additional Z filter between CP and the dephasing period (see above) and the rf decoupling, respectively. We remind the reader that the Z filter suppresses only the "antiphase" 15 N-1 H RIDER, but not the "in-phase" one. However, the "in-phase" RIDER distortion of the COS component at reasonably short N T R is practically negligible, as our data demonstrate. Figures 12 and 13 present the mixing time dependences of the CSA CODEX at various combinations of the 15 N-1 H and 15 N-2 H suppression tools for GB1 and SH3 protein samples. It is seen that the dominant contribution to the short component in the mixing time dependence (Fig. 2) comes from the 15 N-2 H RIDER. Still 2 H decoupling alone cannot ensure the artefact-free experiment, and only combination of Z filter and 2 H decoupling provides the flat mixing time dependence on the millisecond timescale for both proteins. This demonstrates that both SH3 and GB1 proteins in microcrystalline form do not undergo global motions on the millisecond timescale, and the overall rocking motion is limited to the microsecond range only.
The mixing time dependences in Figs. 12 and 13 also reveal a rather slow decay with a time constant in the second  range. This is spin diffusion between 15 N nuclei, which is easy to prove. The spin-diffusion rate should not depend on temperature and should depend on the MAS rate (Reichert et al., 2001;Krushelnitsky et al., 2006). We measured the mixing time dependence for the GB1 sample at two temperatures and two MAS rates; see Fig. 14. The results shown in this figure leave no doubts that this is the ordinary proton-driven spin diffusion. The rate of these decays is approximately 3-4 times slower than the spin-diffusion rate in a fully protonated protein (Krushelnitsky et al., 2006); still, it is quite appreciable. Spin diffusion rate in SH3 protein is noticeably slower; we believe this is due to the lower proton density in this sample, which is confirmed by a somewhat weaker signal from SH3 compared to GB1. Figure 14. 15 N CSA CODEX mixing time dependences measured in GB1 at different MAS rates and temperatures. In all cases Z filter and 2 H CW decoupling during the de(re)phasing periods were applied (the parameters are the same as mentioned in the caption to Fig. 12). NT R 2 ms.

Dipolar CODEX
The principal problem of the dipolar CODEX is that the 15 N-1 H interaction cannot be decoupled for obvious reasons and thus the anti-phase term responsible for the RIDER effect emerges explicitly during the de(re)phasing periods even without CP. To solve this problem, in our first paper on dipolar CODEX (Krushelnitsky et al., 2009) we suggested to measure only the COS component of the mixing time dependence. The COS component must be RIDER-free, which directly follows from Eq. (2). In the dipolar CODEX CSA = 0, and since cos( D ) = cos( D + D ) (we repeat again that this is valid only for I = 1/2), the COS component of the dipolar CODEX mixing time dependence should not be affected by the 15 N-1 H RIDER. However, this is only true under the condition that we did not pay a proper attention to at that time. This condition is that the dipolar interaction must be constant during the de(re)phasing periods, i.e. the timescale of the I -spin flips should be much longer than N T R . If this is not so, then cos( D ) = cos( D + D ) since D and D + D are randomly modulated by I -spin flips within the de(re)phasing periods. Thus, the COS component at this condition is not RIDER-free anymore.
The comparability of N T R and the timescale of proton spins flips is exactly our case. We have estimated above the characteristic time of the protons flip-flops, which is about 10 ms (Fig. 10). The duration of the de(re)phasing periods in the CODEX experiments is usually from few hundred microseconds to several milliseconds. This is shorter than 10 ms but still comparable, which is enough for the RIDER effect. From this we pessimistically conclude that the X-H dipolar CODEX experiment even in proton-diluted samples like deuterated proteins with a partial back-exchange of labile protons is not suitable for studying slow molecular dynamics -there will always be RIDER distortions. This means that the decay in the dipolar CODEX mixing time dependences of backbone nitrogens in SH3 protein that we observed earlier (Krushelnitsky et al., 2009) is not due to molecular motions but due to the RIDER effect and that these data were misinterpreted. The dipolar CODEX experiment, however, can be implemented using other nuclei pairs, e.g. 13 C-15 N (McDermott and Li, 2009), ensuring that the flip-flop time is much longer than the duration of the de(re)phasing periods.
The last point that deserves to be discussed here is the influence of the 15 N-2 H interaction on the dipolar CODEX results. At first sight, there should be no influence, since this interaction is not reintroduced in the dipolar CODEX and it should be simply suppressed by MAS. However, this is not the case. Figure 15 presents the mixing time dependences in GB1 measured at different powers of the CW-2 H decoupling during the de(re)phasing periods. The data demonstrate that in spite of MAS, the 15 N-2 H interaction has a small but very visible contribution to the short component, i.e. RIDER effect, of the mixing time dependence. The residual 15 N-2 H interaction is rather small since only a few kHz of CW decoupling is enough to suppress it completely. So, why does MAS not do its job alone, without the rf decoupling? The reason is the protein mobility on the microsecond timescale. If the 15 N-2 H interaction is modulated by a molecular motion on a timescale of the MAS period (for 20 kHz it is 50 µs), then MAS cannot suppress this interaction completely, which leads to the increased line widths of the MAS centerband (Suwelack et al., 1980). As we know, the correlation time of the protein rocking motion is few tens microseconds (Krushelnitsky et al., 2018). On top of that, there can be interaction of protein nitrogens with deuterons of solvent molecules, and these molecules can also reveal a mobility on the microsecond timescale. Thus, the appearance of the residual 15 N-2 H interaction after MAS can be reasonably explained.
In summary, we have shown that both 15 N-2 H and 15 N-1 H RIDER effects contribute to the short component of the mixing time dependences of both CSA and dipolar CODEX experiments in the protein samples. However, the dominant contributions in these two experiments are different: in the CSA CODEX the dominant source of the short component is the 15 N-2 H interaction, and in the dipolar CODEX it is the 15 N-1 H interaction. As estimated above, the time constants of the two RIDER effects are similar but not the same: 2 H spin-lattice relaxation is somewhat slower than the proton flip-flop rate. Therefore, the apparent decay rate of the short component in the CSA and dipolar CODEX experiments should also be somewhat different. This is illustrated in Fig. 16, which presents the fast RIDER components of the CSA and dipolar CODEX experiments after subtraction of the spin-diffusion component and normalization of the decay amplitudes to the same value. The direct comparison of these decays is in a perfect agreement with the findings described above. It is interesting to note that in SH3, the difference of the apparent correlation times of the short component for the   Fig. 12 (CSA CODEX, without Z filter and without 2 H decoupling) and Fig. 15 (dipolar CODEX, no 2 H decoupling). Red solid lines -the exponential fits of the initial parts of the spin-diffusion components. (b) Fast initial components of the decays after subtraction the spin-diffusion components and normalization to the same initial amplitude.
CSA and dipolar CODEX is much smaller; see Fig. 2 (τ c as a function of the residue number). This can also be reasonably explained by the different proton density in the GB1 and SH3 samples: the less the proton density, the slower the flip-flop rate and thus the smaller the difference between the rates of proton spin diffusion and deuteron spin-lattice relaxation.