# Cross-correlated relaxation measurements under adiabatic sweeps: determination of local order in proteins

Pavel Kadeřávek, Sarina Grutsch, Nicola Salvi, Martin Tollinger, Lukáš Žídek, Geoffrey Bodenhausen, Fabien Ferrage

► To cite this version:

Pavel Kadeřávek, Sarina Grutsch, Nicola Salvi, Martin Tollinger, Lukáš Žídek, et al.. Cross-correlated relaxation measurements under adiabatic sweeps: determination of local order in proteins. Journal of Biomolecular NMR, 2015, 63 (4), pp.353-365. 10.1007/s10858-015-9994-8. hal-01223722

HAL Id: hal-01223722

https://hal.sorbonne-universite.fr/hal-01223722

Submitted on 3 Nov 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Cross-correlated relaxation measurements under adiabatic sweeps: determination of local order in proteins

Pavel Kadeřávek<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="4" />,<fnref n="5" />,<fnref n="6" /> · Sarina Grutsch<fnref n="7" /> · Nicola Salvi<fnref n="6" /> · Martin Tollinger<fnref n="7" /> · Lukáš Žídek<fnref n="4" />,<fnref n="5" /> · Geoffrey Bodenhausen<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,<fnref n="6" /> · Fabien Ferrage<fnref n="1" />,<fnref n="2" />,<fnref n="3" />

Received: 22 May 2015 / Accepted: 14 October 2015
© The Author(s) 2015. This article is published with open access at Springerlink.com

**Abstract** Adiabatically swept pulses were originally designed for the purpose of broadband spin inversion. Later, unexpected advantages of their utilization were also found in other applications, such as refocusing to excite spin echoes, studies of chemical exchange or fragment-based drug design. Here, we present new experiments to characterize fast (ps–ns) protein dynamics, which benefit from little-known properties of adiabatic pulses. We developed a strategy for measuring cross-correlated cross-relaxation (CCCR) rates during adiabatic pulses. This experiment provides a linear combination of longitudinal and transverse CCCR rates, which is offset-independent

across a typical amide ¹⁵N spectrum. The pulse sequence can be recast to provide accurate transverse CCCR rates weighted by the populations of exchanging states. Sensitivity can be improved in systems in slow exchange. Finally, the experiments can be easily modified to yield residue-specific correlation times. The average correlation time of motions can be determined with a single experiment while at least two different experiments had to be recorded until now.

**Keywords** NMR · Protein dynamics · Adiabatic sweep · Cross-correlated cross relaxation

**Electronic supplementary material** The online version of this article (doi:10.1007/s10858-015-9994-8) contains supplementary material, which is available to authorized users.

<figure><img src="image_3.png" /></figure>

Fabien Ferrage
fabien.ferrage@ens.fr

## Introduction

Protein dynamics are essential for most biological processes. The possibility of studying dynamics on various timescales at atomic resolution makes the analysis of NMR relaxation rates unique among other biophysical methods. In principle, motions of any arbitrary spin can be studied, but relaxation studies of ¹⁵N–¹H spin pairs in peptide bonds have been popular since they offer a series of benefits (i) they can be easily introduced into proteins by ¹⁵N labeling, (ii) the approximation of isolated spin pairs can be safely assumed to be fulfilled, provided all other protons are replaced by deuterons, (iii) they provide a residue-specific description of protein backbone dynamics, and finally (iv) the measurement techniques are quite robust (Korzhnev et al. 2001; Lakomek et al. 2012).

According to semi-classical NMR relaxation theory (Wangsness and Bloch 1953; Redfield 1965) all relaxation rates are determined by linear combinations of discrete values of the spectral density function J(ω) which describes the probability of finding motions at a given

<footnote n="1">¹ École Normale Supérieure - PSL Research University, Département de Chimie, 24 rue Lhomond, 75005 Paris, France</footnote>

<footnote n="2">² Sorbonne Universités UPMC Univ Paris 06, LBM, 4 place Jussieu, 75005 Paris, France</footnote>

<footnote n="3">³ UMR 7203 LBM, CNRS, 75005 Paris, France</footnote>

<footnote n="4">⁴ National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic</footnote>

<footnote n="5">5 Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic</footnote>

<footnote n="6">⁶ Institut des Sciences et Ingénierie Chimiques, École polytechnique fédérale de Lausanne, 1015 Lausanne, Switzerland</footnote>

<footnote n="7">⁷ Institute of Organic Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, Innrain 80/82, 6020 Innsbruck, Austria</footnote>
frequency. The determination of the values of the spectral
density function and their interpretation is known as
spectral density mapping (Peng and Wagner 1992a, b). A
simplified variant is called reduced spectral density map-
ping (Ishima and Nagayama 1995a, b; Farrow et al. 1995).
Later, it has been proposed to complement reduced spectral
density mapping analysis by protocols that utilize CCCR
rates (Kadeřávek et al. 2014; Kroenke et al. 1998). Robust
symmetrical reconversion methods allow accurate mea-
surement of transverse CCCR rates (Pelupessy et al. 2003),
but the precision is significantly reduced in the presence of
slow exchange when the precession frequency changes
during the exchange process.

Relaxation under adiabatic pulses has been studied in the presence of protein-ligand interactions (Auer et al. 2010), slow motions in the µs-ms regime (Mulder et al. 2001; Mangia et al. 2010; Auer et al. 2011) and in the presence of fast (ps-ns) dynamics (Konrat and Tollinger 1999). Here, we present measurements of CCCR rates under adiabatically swept pulses and show that deleterious effects of chemical exchange are significantly reduced.

In addition, we introduce variants of experiments that
can directly provide values of the spectral density func-
tion at zero frequency J(0). These J(0) values provide
interesting information about the timescales of the
dominant motion in the ps-ns range. Local disorder in a
protein can be readily identified by a decrease of J(0)
values in flexible segments where the weights of rapid
internal motions characterized by small correlation times
increase. To the best of our knowledge, this is the first
method that allows one to measure J(0) directly in a
single experiment.

**Materials and methods**

Protein samples

Uniformly ¹⁵N labeled and deuterated KID-binding domain (KIX) [residues 586–672 of human CREB-binding protein (CBP)] was expressed using *Escherichia coli* strain BL21 (DE3). Cells were grown in minimal D₂O medium containing ¹⁵N labeled ammonium chloride as sole nitrogen source. KIX was purified as described elsewhere (Tollinger et al. 2006), and the exchangeable deuterons were back-exchanged in H₂O buffer. The KIX sample contained 1 mM protein in 50 mM potassium phosphate buffer, pH 5.8, and 25 mM NaCl before adding 10 % D₂O.

Human ubiquitin was purchased from VLI. The sample
was prepared by dissolution of lyophilized ubiquitin in
50 mM acetate buffer (pH 4.6 before 10 % of D2O was
added to the solution). The final protein concentration was
0.67 mM.

NMR experiments

The experiments were performed in magnetic fields B₀ = 11.74 and 18.79 T on Bruker Avance spectrometers (¹H Larmor frequencies 500 and 800 MHz, respectively). Both spectrometers were equipped with TCI cryoprobes. The KIX domain was studied at 20 °C at 500 MHz. Experiments with ubiquitin were carried at 5 °C (500 MHz) and 30 °C (500 and 800 MHz). The temperature was calibrated by measuring the difference between ¹H chemical shifts of a 4 % solution of methanol in deuterated methanol.

Symmetrical reconversion experiments (Pelupessy et al. 2003, 2007) were used to measure the transverse ($\eta_{xy}$) and longitudinal ($\eta_z$) cross-relaxation rates due to cross-correlated fluctuations of the ¹⁵N chemical shielding anisotropy (CSA) and the ¹⁵N–¹H dipole–dipole interaction in ubiquitin. The transverse cross-relaxation rate $\eta_{xy}$ was measured with relaxation delays $T = 50$ and 70 ms, the longitudinal cross-relaxation rate $\eta_z$ with $T = 100, 175$, and 250 ms. Experiments with adiabatic sweeps employed $T_{adiab} = 60$ and 80 ms long apodized chirp pulses (Böhlen et al. 1990). The maximum amplitude of all chirp pulses was applied for 60 % of the duration of the pulses. The chirp pulses were

<figure><img src="image_4.png" /></figure>

**Fig. 1** Variants of relaxation intervals a scheme 1—with adiabatically swept pulse and b scheme 2—with single spin echo. Rectangular $\pi/2$ and $\pi$ pulses are denoted as *open and filled boxes*. All pulses are applied with phase along x axis unless indicated. The phases $\phi_1$ and $\phi_2$ follow cycles {y,−y} and {y,y,−y,−y}, respectively. The long *grey shaped* pulse represents the adiabatic inversion pulse of length $T_{adiab}$ where the amplitude is modulated by two ramps of length $T_r$ each. Since type selective pulses applied to the water resonance are represented by *semi-elliptic symbols*. The carrier frequency for all hard pulses was placed in the centre of the amide region while all flip back pulses were applied on-resonance with respect to the water signal. The amplitudes of the selective pulses preceding and following the composite pulses were calibrated separately
divided into 10,000 discrete steps to cover a sweep width of 10 kHz. The central part of the relaxation blocks ($T_{xy}$ in Fig. 1) of our experiments with a single $^{15}$N echo and non-selective $^{15}$N pulses was 50 or 70 ms long.

The transverse CCCR rates $\eta_{xy}$ of KIX were determined by a single experiment with a relaxation delay $T = 38.768$ ms and by our experiment with a chirp pulse of a length $T_{adiab} = 80$ ms (all other parameters were as for the experiments performed on ubiquitin).

The longitudinal relaxation rates $R_1(^{15}\text{N})$ in ubiquitin were measured (Korzhnev et al. 2001) with relaxation delays $T = 140$ (twice), 220, 380, 500, 540 (twice), and 660 ms at 500 MHz and with $T = 20, 60, 140$ (twice), 220, 380, 500, 540 (twice), and 660 ms at 800 MHz.

The $^{15}\text{N}-^{1}\text{H}$ steady-state nuclear Overhauser effect $\sigma\{^{1}\text{H}\}$ (ssNOE) in ubiquitin at 800 MHz and 30 °C was determined from the ratio between the intensities of spectra acquired with proton saturation achieved by a 8 s long repetition of a segment consisting of a 180 ° $^{1}\text{H}$ pulse and a 10.87 ms delay and the intensities in reference spectra acquired with a 17 s interscan delay (Ferrage et al. 2009). The steady-state and reference spectra were measured in interleaved fashion. The same experiments were performed at 500 MHz and 5 °C, but the duration of the proton saturation period was 6 s and the interscan delay in the reference experiments was 15 s.

The transverse relaxation rates $R_2(^{15}\text{N})$ of ubiquitin were determined from the relaxation rates $R_1$ and $R_{1\rho}$ (Mulder et al. 2001):

$$R_2 = \frac{R_{1\rho} - R_1 \cos^2 \theta}{\sin^2 \theta} \qquad (1)$$

where $\theta = \text{atan}(\omega_1/\Omega)$ is the tilt angle of the effective field during the spin-lock irradiation, $\Omega$ being the offset between the resonance frequency and the carrier frequency and $\omega_1$ the rf amplitude. The carrier frequency during the spin-lock was set in the middle of the amide region for durations $T = 0, 20, 40$ (twice), 60, 80, 100, 120, and 140 (twice) ms.

The rf amplitudes for spin locking and chirp pulses were calibrated using a comparison of residual and true $J(^{15}\text{N}-^{1}\text{H}^{\text{N}})$ coupling constants in spectra acquired with continuous-wave decoupling (Palmer et al. 2001). The calibrated rf amplitudes of chirp pulses were 2303.8 ± 1.1, 2343.8 ± 0.7, 2278.2 ± 1.1, and 2306.1 ± 3.0 Hz for measurements of ubiquitin at 800 MHz (30 °C), 500 MHz (30 °C), 500 MHz (5 °C), and KIX domain, respectively. The calibrated rf amplitudes used for spin locking in $R_{1\rho}$ experiments were 1617.2 ± 0.5, and 1040.6 ± 0.5 Hz for experiments at 800 MHz (30 °C), and 500 MHz (5 °C), respectively. Note that the errors were evaluated as asymptotic standard errors of the fit and are likely underestimated but this has no consequence on the analysis.

A compensation block was inserted to ensure uniform sample heating in measurements of $R_1$, $R_{1\rho}$, $\eta_{xy}$ and $\eta_z$.

All spectra were processed with NMRPIPE (Delaglio et al. 1995) and analyzed with SPARKY (Goddard and Kneller 2006). The $R_1$ and $R_{1\rho}$ relaxation rates were obtained from fits of the signal decays to mono-exponential functions using the program OCTAVE 3.2.4 (Eaton et al. 2009).

The CCCR rates $\eta_{xy}$, $\eta_z$, and their linear combinations $\eta_{ave}$ were obtained from a fit of the observable $Y(T)$ defined as:

$$Y(T) = \sqrt{\frac{I_{P,Q}(T)I_{Q,P}(T)}{I_{P,P}(T)I_{Q,Q}(T)}} \qquad (2)$$

where $I_{P,Q}(T)$ denotes a peak intensity in spectra derived from operator terms $P$ and $Q$ selected before and after the relaxation period $T$ in the manner of symmetrical reversion. $P$ and $Q$ stand for $N_z$ and $2N_zH_z$ operators, respectively. The rates $\eta_{xy}$ and $\eta_{ave}$ were determined from a linearized equation:

$$\text{atanh}(Y(T)) \approx \eta_i T, \qquad (3)$$

with $\eta_i = \eta_{xy}$ or $\eta_{ave}$. $\eta_z$ was determined from [Eq. 7 in Pelupessy et al. (2007)]

$$\text{atanh}(Y^{-1}(T)) = \eta_z T \left( 1 + \frac{\mu^2 T^2}{24} \right) \qquad (4)$$

where $\mu$ is a correction factor determined from a simultaneous fit of [Eq. 8 in Pelupessy et al. (2007)]:

$$\frac{I_{P,P}(T)}{I_{Q,Q}(T)} = \frac{\lambda \exp(-\mu T)}{6 + \eta_z^2 \mu T^3} \qquad (5)$$

where $\lambda$ is a fitted pre-exponential factor.

The errors of the peak intensities were estimated based on random noise in the spectra and the uncertainties of the relaxation rates were estimated from 2000 independent fits of the relaxation rates while the peaks intensities were varied according to their estimated errors assuming a normal error distribution.

The values of spectral density functions $J(\omega)$ were calculated using the reduced spectral density mapping approach, with the internuclear distance $r_{\text{N-H}} = 1.02$ Å, the anisotropy of the nitrogen chemical shielding tensor $\Delta\sigma = \sigma_{\parallel} - \sigma_{\perp} = -170$ ppm and the angle between the unique axis of the CSA tensor and the N-H bond to 20.6°. The choice of the parameters was justified by a comparison of the spectral density values at the $^{15}\text{N}$ frequency obtained from the CCCR rates and the reduced spectral density mapping (Kadeřávek et al. 2014).

Simulations of the evolutions of the spin system (Levitt and Bari 1992; Ghose 2000) during our pulse sequence
were performed with MATLAB R2014a [MATLAB] using SPINACH 1.5.2440 (Hogben et al. 2011). The full basis set was used with Redfield relaxation theory (Wangsness and Bloch 1953; Redfield 1965) using the secular approximation. Simulations for a two-site exchange between states on different time scales were performed by a script written in MATHEMATICA 9 (Wolfram Research Inc. 2012).

# Theory

## Cross-correlated cross-relaxation

The relaxation of an isolated ¹H–¹⁵N spin pair can be described by the semi-classical relaxation theory (Wangsness and Bloch 1953; Redfield 1965). Cross-correlation between the ¹⁵N chemical shielding anisotropy (CSA) and the dipole–dipole interaction with the HN proton is responsible for cross-relaxation pathways between the in-phase and anti-phase ¹⁵N coherences (the transverse CCCR rate ηxy) and between the polarization Nz and the two-spin order 2NzHz (the longitudinal CCCR rate ηz). The CCCR rates depend on the spectral density function:

$$ \eta_{xy} = C_{cd}(8J(0) + 6J(\omega_N)), \quad (6) $$

$$ \eta_z = 12C_{cd}J(\omega_N), \quad (7) $$

where ωN is the nitrogen Larmor frequency, Ccd = (3 cos² φ - 1)γN²γH B₀Δσμ₀ħrN⁻³/96π, γH and γN are the gyromagnetic ratios of ¹H and ¹⁵N, respectively, rN–H is the H–N internuclear distance, μ₀ is the permeability of vacuum, ħ is Planck’s constant divided by 2π, Δσ is the anisotropy of the ¹⁵N chemical shielding tensor, φ is the angle between the H–N bond and the symmetry axis of the ¹⁵N chemical shielding tensor, and B₀ is the external magnetic field. The spectral density function is described as a sum of k Lorentzian functions J(ω) = (2/5) Σk-1 i=0 aᵢτᵢ/ (1 + (ωτᵢ)²), where τᵢ are the effective correlation times and aᵢ are weighting constants. Hence, J(0) = (2/5)τ₀ for motions with a correlation function described by a mono-exponential decay characterized by a single correlation time τ₀, and J(0) = (2/5) Σk-1 i=0 aᵢτᵢ for more general cases where statistically independent motions participate. Within the frame of the popular model-free approach (Lipari and Szabo 1982a, b; Halle et al. 1981; Halle 2009) a local order parameter is defined as S² = a₀, where the index i = 0 refers to the overall tumbling of the molecule in the solution with a correlation time τ₀ and τ₁ = (τ⁻¹₀ + τ⁻¹₁)⁻¹, where τf is the timescale of the fast internal motions.

## Adiabatic sweeps

The adiabatic fast passage pulses (Abragam 1961; Böhlen and Bodenhausen 1993; Kupče and Freeman 1995) used in the relaxation period (Fig. 1) are phase-modulated pulses. The time dependence of the phase effectively sweeps the resonance frequency of the pulse in linear fashion over a defined range. The time dependence of the carrier frequency for the chirp pulses used in this study is:

$$ \Omega(t) = 2\pi \left( \frac{t}{T_{\text{adiab}}} - 0.5 \right) W \quad (8) $$

where Tadiab is the total pulse duration and W is the sweep width in Hz. The amplitude of the chirp pulse starts from zero and smoothly increases like the first quarter of the period of a sine function during a ramping interval 0<t<Tr:

$$ \omega_1(t) = 2\pi B_1 \sin \frac{\pi t}{2T_r} \quad (9) $$

where B₁ is the maximum of rf amplitude in Hz, and Tr is the length of the ramp. At time (Tadiab - Tr) the amplitude of the rf pulse is gradually decreased according to equation:

$$ \omega_1(t) = 2\pi B_1 \sin \frac{\pi(t - T_{\text{adiab}} + 2T_r)}{2T_r} \quad (10) $$

The amplitude is constant ω₁(t) = 2πB₁ for t ∈ [Tr, Tadiab - Tr]. In order to evaluate the effective relaxation of a ¹⁵N–¹H spin pair during the chirp pulse, average Liouvillian theory can be applied (Levitt and Bari 1992; Ghose 2000) as average Hamiltonian theory (Valentine et al. 2007). The relaxation can be described in a spin operator subspace {Nz, 2NzHz, Nx, 2NzHz} in the instantaneous rotating frame (assuming that the chirp pulse was applied along the x-axis of the accelerating rotating frame). The zero-order average Liouvillian over the whole length of the pulse is given by:

$$ \bar{\mathbf{L}} = \frac{1}{T_{\text{adiab}}} \int_{0}^{T_{\text{adiab}}} \mathbf{L}'(t) \mathrm{d}t. \quad (11) $$

L'(t) is the Liouvillian at a time t which can be calculated as:

$$ \mathbf{L}'(t) = \mathbf{U}^{-1}(t)\mathbf{L}\mathbf{U}(t), \quad (12) $$

where L is the Liouvillian of the initial state with the following matrix representation in the basis {Nz, 2NzHz, Nx, 2NzHz}:
$$ \mathbf{L} = \begin{pmatrix} R_1 & \eta_z & 0 & 0 \\ \eta_z & R'_1 & 0 & 0 \\ 0 & 0 & R_2 & \eta_{xy} \\ 0 & 0 & \eta_{xy} & R'_2 \end{pmatrix} \qquad (13) $$

where $R_1$ is the longitudinal $^{15}$N relaxation rate, $R'_1$ is the relaxation rate of the two spin order and $R_2$ and $R'_2$ are the transverse relaxation rates of the in-phase and anti-phase $^{15}$N coherences, respectively. $\eta_z$ and $\eta_{xy}$ are the longitudi-

interval $[0, T_{adiab}]$ vanishes. We introduce the parameter $\alpha$:

$$ \alpha = \frac{1}{T_{adiab}} \int_{0}^{T_{adiab}} \sin^2(\theta(t)) dt \qquad (18) $$

which can be evaluated numerically. The averaged Liouvillian has the form:

$$ \bar{\mathbf{L}} = \begin{pmatrix} R_1(1-\alpha) + R_2\alpha & \eta_z(1-\alpha) + \eta_{xy}\alpha & 0 & 0 \\ \eta_z(1-\alpha) + \eta_{xy}\alpha & R'_1(1-\alpha) + R'_2\alpha & 0 & 0 \\ 0 & 0 & R_2(1-\alpha) + R_1\alpha & \eta_{xy}(1-\alpha) + \eta_z\alpha \\ 0 & 0 & \eta_{xy}(1-\alpha) + \eta_z\alpha & R'_2(1-\alpha) + R'_1\alpha \end{pmatrix} \qquad (19) $$

nal and transverse CCCR rates. $U(t)$ is the transformation into the time-dependent tilted frame.

$$ U(t) = \begin{pmatrix} c & 0 & -s & 0 \\ 0 & c & 0 & -s \\ s & 0 & c & 0 \\ 0 & s & 0 & c \end{pmatrix}, \qquad (14) $$

where $c$ and $s$ stand for $\cos(\theta(t))$ and $\sin(\theta(t))$, respectively. The angle $\theta(t)$ can be calculated with the knowledge of the time-dependence of the offset and the amplitude (Eqs. 8–10)

$$ \cos \theta(t) = \frac{\Omega(t)}{\sqrt{\Omega^2(t) + \omega_1^2(t)}} \qquad (15) $$

$$ \sin \theta(t) = \frac{\omega_1(t)}{\sqrt{\Omega^2(t) + \omega_1^2(t)}} \qquad (16) $$

Eqs. 12, 13, and 14 provide a Liouvillian in the transformed base $\{N'_z(t), 2N'_z H_z(t), N_{x'}(t), 2N_{x'} H_z(t)\}$:

We provide a list of $\alpha$ parameters calculated for various sets $\{B_1, W, T_r/T_{adiab}\}$ in the supplementary material (Tab. S1). It should be noted that the calculated $\alpha$ parameters and the zero-order averaged Liouvillian presented in Eq. 19 correspond to $^{15}$N with a resonance frequency in the centre of the adiabatic sweep. For all other cases a constant offset must be added to $\Omega(t)$ in Eq. 8. Consequently, the proportion of $\eta_{xy}$ and $\eta_z$ contributions to the average cross-relaxation rates depends on the offset. In addition, the integral of the product $\cos(\theta(t)) \sin(\theta(t))$ does not vanish, so all terms in the matrix in Eq. 19 are non-zero a priori. However, errors are limited under the usual conditions (see below).

## Principles of the new experiments

The experiments that are presented in this paper combine relaxation during an adiabatic pulse with the inversion of proton magnetization and additional relaxation delays. The complete relaxation block (scheme 1) is shown in Fig 1a. It starts with a longitudinal relaxation delay $T_z$ followed by

$$ \mathbf{L}'(t) = \begin{pmatrix} R_1 c^2 + R_2 s^2 & \eta_z c^2 + \eta_{xy}s^2 & (R_2 - R_1)cs & (\eta_{xy} - \eta_z)cs \\ \eta_z c^2 + \eta_{xy}s^2 & R'_1 c^2 + R'_2 s^2 & (\eta_{xy} - \eta_z)cs & (R'_2 - R'_1)cs \\ (R_2 - R_1)cs & (\eta_{xy} - \eta_z)cs & R_2 c^2 + R_1 s^2 & \eta_{xy} c^2 + \eta_z s^2 \\ (\eta_{xy} - \eta_z)cs & (R'_2 - R'_1)cs & \eta_{xy} c^2 + \eta_z s^2 & R'_2 c^2 + R'_1 s^2 \end{pmatrix}. \qquad (17) $$

The integral of $\mathbf{L}'(t)$ over the interval $[0, T_{adiab}]$ gives the average Liouvillian. Note that the integral of the function $\cos(\theta(t)) \sin(\theta(t))$ over the integration

an adiabatic pulse of a length $T_{adiab}$ flanked by two proton inversion pulses and terminated by a second longitudinal relaxation delay $T_z$. The averaged Liouvillian can be easily
calculated as a weighted sum of three Liouvillians representing relaxation during the three delays. The Liouvillian during the first free relaxation period $T_z$ is equal to $\mathbf{L}$ in Eq. 13. The effect of the proton inversion pulses is described by a unitary transformation $\mathbf{U}_H$ in the same basis as Eq. 13:

$$ \mathbf{U}_H = \begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & -1 \end{pmatrix} \qquad (20) $$

The averaged Liouvillian over the relaxation delay $T \approx (2T_z + T_{adiab})$ where the duration of pulses can be neglected is:

$$ \bar{\mathbf{L}}_1 = \frac{1}{T} (2\mathbf{L}T_z + T_{adiab} \mathbf{U}_H^{-1} \bar{\mathbf{L}} \mathbf{U}_H) \qquad (21) $$

where $\bar{\mathbf{L}}$ is defined in Eq. 19.

A similar experiment (scheme 2) can be proposed if the adiabatic sweep is replaced by a $^{15}$N refocusing $\pi$ pulse flanked by two delays $T_{xy}/2$ and two non-selective $^{15}$N $\pi/2$ pulses (Fig 1b). The transformation corresponding to the $^{15}$N refocusing pulse is described by a simple multiplication by the factor $-1$ which has no effect on the Liouvillian. $\mathbf{U}_N$ represents the transformation of the Liouvillian due to a $^{15}$N $\pi/2$ pulse applied along y-axis:

$$ \mathbf{U}_N = \begin{pmatrix} 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & -1 \\ 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \end{pmatrix} \qquad (22) $$

The averaged Liouvillian over the relaxation interval in scheme 2 is thus given by:

$$ \bar{\mathbf{L}}_2 = \frac{1}{2T_z + T_{xy}} (2\mathbf{L}T_z + T_{xy} \mathbf{U}_N^{-1} \mathbf{U}_H^{-1} \mathbf{L} \mathbf{U}_H \mathbf{U}_N) \qquad (23) $$

which is equal to $\bar{\mathbf{L}}_1$ of Eq. 21 if $\alpha = 1$ and $T_{xy} = T_{adiab}$.

# Results and discussion

## Choice of the relaxation intervals

The cross-relaxation measurements use symmetrical reconversion (Pelupessy et al. 2003): $N_z$ or $2N_zH_z$ terms are prepared prior to the relaxation period $T$ and either $N_z$ or $2N_zH_z$ are selected at the end of $T$. All four relaxation pathways are recorded. An analysis of relaxation during our experiments shows that the observable $Y(T)$ (Eqs. 2, 3) results from a linear combination of transverse ($\eta_{xy}$) and longitudinal ($\eta_z$) CCCR rates in the zeroth order average

Liouvillian. The observed averaged CCCR rate $\eta_{ave}$ determined by our experiments is:

$$ \eta_{ave} = \frac{2T_z\eta_z - T_{centr}(\alpha\eta_{xy} + (1-\alpha)\eta_z)}{T} \qquad (24) $$

where $T$ is the total relaxation period and $T_{centr} = T_{adiab}$ or $T_{centr} = T_{xy}$ for scheme 1 or 2 as defined in Fig 1. The parameter $\alpha$ defines the proportion of time that the $^{15}$N magnetization spends in the transverse plane during the $T_{centr}$ period. $\alpha = 1$ for all N-H$^N$ groups for scheme 2 (called single echo experiment in this article) provided the pulses are sufficiently strong. $0 < \alpha < 1$ when an adiabatically swept pulse is used (scheme 1, called adiabatic experiment in this article). If the $T_z$ and $T_{centr}$ periods are chosen to fulfill $T_z = 0.5(1-\alpha)T_{centr}$, the right hand side of Eq. 24 does not depend on $\eta_z$, so that $\eta_{xy}$ can readily be determined:

$$ \eta_{ave} = \eta_{ave}^{xy} = \eta_{xy} \frac{\alpha}{\alpha - 2} \qquad (25) $$

On the other hand, if $T_z = 0.5(1 - 0.5\alpha)T_{centr}$ the right hand side of Eq. 24 does not depend on $J(\omega_N)$. Then, $J(0)$ can be calculated directly from the analysis of this single experiment. Equations 6, 7 and 24 give in this case:

$$ \eta_{ave} = \eta_{ave}^{J} = J(0) \frac{16\alpha C_{cd}}{\alpha - 4} \qquad (26) $$

With $0 < \alpha < 1$, the magnitude of both $\eta_{ave}^{xy}$ and $\eta_{ave}^{J}$ increase with $\alpha$. The $\alpha$ parameter should be maximized to increase the precision of the $\eta_{xy}$ and $J(0)$ determination.

Finally, if the flanking $T_z$ periods are omitted in the adiabatic experiment (scheme 1 in Fig. 1) and if the relaxation delay is reduced to the adiabatic pulse, a linear combination of longitudinal and transverse CCCR rates is obtained

$$ \eta_{ave} = \eta_{ave}^{xyz} = \alpha\eta_{xy} + (1-\alpha)\eta_z \qquad (27) $$

Both $\eta_z$ and $\eta_{xy}$ can be determined from at least two experiments with different $\alpha$ parameters. In principle, this is similar to the standard experiments for the separate determination of $\eta_{xy}$ and $\eta_z$. But it might be more advantageous if an exchange significantly contributes to relaxation of the $^{15}$N magnetization, as discussed below.

## Systematic errors

The $\alpha$ parameter weakly depends on the rf field strength $B_1$, so that the method is robust with respect to a possible miscalibration of the amplitude of the adiabatic pulse. The relative errors are (see Appendix):

$$ \frac{\eta_{xy} - \eta'_{xy}}{\eta_{xy}} = (1 - \frac{\alpha}{\alpha'}) F_{\eta} \qquad (28) $$
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Contributions to the relative errors of J(0) and η<sub>xy</sub> caused by miscalibration of the amplitude of the chirp pulse. **a** A dependence of the factors F<sub>η</sub> and F<sub>J</sub> on the correlation time of rotational diffusion. F<sub>η</sub> is related to the error of η<sub>xy</sub> (red lines) and F<sub>J</sub> is related to the error of J(0) (blue lines). The solid and dashed lines correspond to measurements at 500 MHz and 1 GHz, respectively. **b** A contour plot of the dependence of the error factor (1 - α/α') on the amplitude of the adiabatic pulse B'<sub>1</sub> and its deviation ΔB<sub>1</sub> = B<sub>1</sub> - B'<sub>1</sub> from the nominal value B<sub>1</sub>. The blue, magenta, black, green, and red lines represent contours of the function (1 - α/α') at levels -0.05, -0.025, 0, 0.025, and 0.05, respectively. The solid, dashed, and dotted lines show contours for chirp pulses with 10, 15, and 20 kHz sweep widths, respectively, and with a ratio of the length of the ramp and the total length of the chirp pulse T<sub>r</sub>/T<sub>adiab</sub> = 0.05. **b** The relative errors for slow motion. For faster motions the errors are scaled down by the factor F<sub>η</sub> or F<sub>J</sub> shown in (a)</figcaption></figure>

$$ \frac{J(0) - J'(0)}{J(0)} = \left(1 - \frac{\alpha}{\alpha'}\right) F_J \qquad (29) $$

where J'(0) and η'<sub>xy</sub> are obtained from the analysis of experimental data based on an α' parameter set to an erroneous value, while the correct value was α. The factors F<sub>η</sub> and F<sub>J</sub> depend only on the rotational diffusion of the molecule, its internal motions and the magnetic field of the spectrometer:

$$ F_{\eta} = \frac{4J(0) - 3J(\omega_{\text{N}})}{4J(0) + 3J(\omega_{\text{N}})} \qquad (30) $$

$$ F_J = \frac{4J(0) - 3J(\omega_N)}{4J(0)} \qquad (31) $$

For slow motions, when J(0) ≫ J(ω<sub>N</sub>), the factors converge to F<sub>η</sub> = F<sub>J</sub> = 1 as shown by two examples in Fig. 2a. The protocol employed here (Palmer et al. 2001) ensures that a miscalibration of the B<sub>1</sub> field larger than 100 Hz is unlikely. Therefore, the factor (1 - α/α') was evaluated within the limits |B<sub>1</sub> - B'<sub>1</sub>| < 100 Hz, where B'<sub>1</sub> is an erroneously determined amplitude of the chirp pulse. Fig. 2b shows that errors arising from such B<sub>1</sub> field mis-calibrations are often negligible. Errors increase with increasing sweep width, decreasing B<sub>1</sub> field, and decreasing ratio of the length of ramp and the total length of chirp pulse T<sub>r</sub>/T<sub>adiab</sub>. For clarity only the case T<sub>r</sub>/T<sub>adiab</sub> = 0.05 is shown in Fig. 2b.

As mentioned in the “Theory” section, the average Liouvillian differs from Eq. 19 for residues with <sup>15</sup>N resonance frequencies that are not in the centre of the adiabatic sweep. Therefore, a systematic error arising from off-resonance effects must be considered. The error was simulated for the case of a <sup>15</sup>N nucleus with a 15 ppm resonance offset from the centre of the adiabatic sweep. The single correlation time of the motion was varied over the range 1 < τ<sub>0</sub> < 50 ns. The dependence of the error on the correlation time between these limits was monotonic but non-linear, suggesting that the simulation for τ<sub>0</sub> = 50 ns is close to the slow motion limit. The dependence of the error on the B<sub>1</sub> amplitude of the chirp pulses is shown in Fig. 3 for τ<sub>0</sub> = 10 ns and B<sub>0</sub> = 23.5 T (1 GHz <sup>1</sup>H Larmor frequency). The results for B<sub>0</sub> = 18.79 and 11.74 T (800 and 500 MHz <sup>1</sup>H Larmor frequencies) follow the same trends but errors are scaled down so that they do not exceed 2.7

<figure><img src="image_6.png" /><figcaption>**Fig. 3** Simulated dependence of the relative errors of η<sub>xy</sub> determined by the adiabatic experiment (scheme 1) on the amplitude of the adiabatic pulse [the relative errors of J(0) were almost identical]. The simulations refer to 1 GHz, the <sup>15</sup>N resonance was shifted by 15 ppm from the centre of the frequency sweep of the adiabatic pulse. The red, green, and blue lines show the errors predicted for a 80 ms chirp pulse with sweep widths of 10, 15, and 20 kHz, respectively. The solid, dashed and dotted lines show errors for chirp pulses with the ratios of the length of the ramp and the total duration of the chirp pulse T<sub>r</sub>/T<sub>adiab</sub> = 0.05, 0.1, and 0.2, respectively</figcaption></figure>
and 1.1 %, respectively. The absolute values of the relative
errors calculated for extreme correlation times τ₀ = 1 and
50 ns do not exceed 1.8 and 4.3 %, respectively, even at 1
GHz.

The theory is valid for an isolated pair of spins. An
amide ¹⁵N-¹H spin pair in an otherwise deuterated protein
sample fulfills these assumptions to a very good approxi-
mation. However, effects of H² and other protons should be
considered for partially deuterated or fully protonated
protein samples. We performed simulations for a spin
system where dipole-dipole interactions with two addi-
tional protons were considered. The protons were placed
2.1, 2.5, and 2.9 Å from both amide nitrogen and hydrogen.
The three studied cases approximate fully protonated,
partially deuterated, and fully deuterated protein samples,
respectively. First, the simulations show (see Supplemen-
tary information) that the error rises significantly with
decreasing distance between the amide spin pair and the
additional interacting protons. Second, the accuracy
decreases with increasing correlation time of a motion of
an amide spin pair. Third, the error rises with the length of
the relaxation delay *T*, and finally, the error is larger if an
adiabatic pulse with a smaller *α* parameter is used (selected
dependences of the error on the correlation time and the
adiabatic pulse parameters are shown in Supplementary
information). Generally, the results of the simulations show
that a high level of deuteration is required for a quantitative
analysis if methods with adiabatically swept pulses
(scheme 1) are used. Although a partial deuteration can
suffice in small proteins, scheme 2 should be favoured for a
protonated or partially deuterated sample. The enhanced
performance of scheme 2 comes from the averaging of the
auto-relaxation rates of the in-phase and anti-phase terms
effectively performed by the evolution under scalar cou-
pling. The use of a conversion block between in-phase and
anti-phase terms (Kroenke et al. 1998; Ghose et al. 1999)
in scheme 1 should enhance the accuracy of the
measurements.

Finally, the effect of amide proton exchange was con-
sidered. The simulations show that the results of experi-
ments with an adiabatically swept pulse (scheme 1) are
corrupted in the presence of very fast proton exchange. The
error rises with increasing proton exchange rate and
decreasing correlation time of the studied amide spin pair.
The error is larger for longer relaxation delays *T* and/or if
an adiabatic pulse with a smaller *α* parameter is used (see
Supplementary information). For a typical experimental set
up similar to the one used in this study (simulated for an
experiment performed at 500 MHz with an adiabatic pulse
with the following parameters: *B*₁ = 2.3 kHz, *SW* = 10 kHz,
*T*<sub>adiab</sub> = 80 ms, *T*<sub>r</sub>/*T*<sub>adiab</sub> = 0.2) a maximum exchange rate
of 4 Hz was tolerable to keep the error below 3.5 % for

residues with correlation time τ = 1 ns. Because faster
proton exchange rates are typical for disordered proteins
(Croke et al. 2005) the scheme 2 should be the method of
choice in this case. The error was lower than 1.5 % even
for the simulation with a proton exchange rate equal to
20 s⁻¹, a correlation time τ = 1 ns, and a relaxation delay
T = 120 ms. Like for proton-proton dipolar relaxation, fast
averaging of in-phase and anti-phase auto-relaxation rates
makes scheme 2 more accurate in these cases.

Experimental results

The experiments were first tested on ubiquitin. Spectral
density values J(0) and experimental relaxation rates η<sub>xy</sub>,
and η<sub>ave</sub><sup>xyz</sup> were compared with the expected values derived
from transverse and longitudinal CCCR rates measured
with standard experiments. The values obtained by both
adiabatic and single echo experiments (schemes 1 and 2 in
Fig. 1) are in good agreement with the conventional
approach (Fig. 4). As expected, the J(0) values obtained
from any analysis based on CCCR rates are not contami-
nated by chemical exchange, in contrast to J(0) values
obtained from classical reduced spectral density mapping
(Ishima and Nagayama 1995b, a; Farrow et al. 1995). A
comparison of the determined J(0) values is shown in
Fig. 5.

Our method was also applied to the KIX domain of
CBP, which is known to have a significant amount of slow
conformational exchange (Tollinger et al. 2006). The
transverse CCCR rates were measured with the adiabatic
experiment (scheme 1 in Fig 1) and conventional experi-
ments. The results displayed in Fig. 6 show a very good
agreement between the two methods. However, our
experiment yields η<sub>xy</sub> values that are systematically lower
for residues in the first α helix of KIX (the average dif-
ference is 0.4 s<sup>-1</sup>), while the agreement between the data in
other parts of the sequence is satisfactory. It has been
reported (Tollinger et al. 2006) that helices 1 and 2
undergo a slow exchange between a folded and a partially
un-folded state. H/D exchange protecting factors (Schanda
et al. 2008) were shown to be higher in helix 2 than helix 1,
suggesting the persistence of residual structure for helix 2
in the partially unfolded state.

In order to identify the source of the disagreement,
simulations of both experiments were performed. Two-
state exchange was included in the Liouvillian matrix and
the exchange and relaxation parameters were chosen to
correspond to published values for KIX (Tollinger et al.
2006). Relaxation due to fast motions (ps–ns) was treated
following the semi-classical relaxation theory (Wangsness
and Bloch 1953; Redfield 1965) assuming that only an
unrestricted isotropic motion with a single correlation
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Correlations between the results of our experiments and standard experiments. Panels **a–d**, **e–h** and **i–l** correspond to the results measured for ubiquitin at 800 MHz and 30 °C, 500 MHz and 30 °C, and 500 MHz and 5 °C, respectively. The comparisons of *J*(0) measured using the single-echo variant of our pulse program with non-selective nitrogen pulses and with adiabatically swept pulse are</figcaption></figure>

shown in panels (a, e, i) and (b, f, j), respectively. The correlations of the transverse CCCR are shown in panels (c, g, k) and the linear combinations of transverse and longitudinal CCCR rates are in panels (d, h, l). Root mean square deviation (RMSD) and correlation parameter (CORR) are shown in each panel

time $\tau_0$ contributes to relaxation. Two simulations were performed in which the slow exchange parameters were identical, but the parameters defining the fast motions differed. In the first simulation, the correlation times of both exchanging states were equal ($\tau_0$ = 10 ns), while in the second case the motion of the less populated state was ten times faster ($\tau_0$ = 1 ns). The results of the simulations of both the standard experiment and of the adiabatic experiment (scheme 1 in Fig. 1) were identical in the first case but different in the second case. While the simulations of the standard experiment provide very similar results in both cases, the result of the adiabatic experiment is affected by the short correlation time of the less populated excited state. If the standard method to measure $\eta_{xy}$

is used, both in-phase and anti-phase coherences of the excited state are lost rapidly due to enhanced relaxation caused by the exchange. This effect reduces the weight of cross-relaxation in the excited state to a negligible value. So, the measured value reflects mostly the major folded state. On the contrary, the adiabatic sweep efficiently suppresses the loss of coherence due to exchange in the minor state. Therefore, the measured value represents a population-weighted $\eta_{xy}$ rate that is determined by both the ground and excited states, provided the exchange rate is fast enough to explore both states during the relaxation delay.

Theoretically, knowledge of the major state $\eta_{xy}$, of the population weighted $\eta_{xy}$, and the population of the excited
<figure><img src="image_8.png" /><figcaption>**Fig. 5** J(0) values in ubiquitin measured **a** at 800 MHz and 30°C, and **b** at 500 MHz and 5°C as a function of the residue number. The red symbols represent spectral density values obtained by reduced spectral density mapping, the green symbols correspond to data measured by our experiment using an adiabatically swept pulse</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 6** Transverse cross-correlated cross-relaxation rates for different residues in the KIX domain. The red and green symbols show data measured by the standard method and by our pulse program using adiabatically swept pulse in the relaxation period. The blue horizontal bars show the positions of three main α-helices in the KIX domain</figcaption></figure>

state [based on relaxation dispersion experiments (Tollinger et al. 2006; Schanda et al. 2008)] allows one to calculate the CCCR rate of the excited state. Unfortunately, the precision of our measurements is not sufficient to obtain reasonably precise CCCR rates of the excited state. Nevertheless, the systematic deviations of the η<sub>xy</sub> values between the adiabatic experiment (scheme 1 in Fig. 1) and

<figure><img src="image_10.png" /><figcaption>**Fig. 7** Simulated dependence of the transverse CCCR rate η<sub>xy</sub> obtained by the standard experiment (dashed lines) and our experiment (Eq. 25) using adiabatically swept pulse (solid lines) on the rate of the exchange. The simulations were performed for a static field of 11.75 T, an isolated <sup>15</sup>N–<sup>1</sup>H spin system undergoing a two-state exchange A ↔ B characterized by equilibrium populations p<sub>A</sub> = 0.9 and p<sub>B</sub> = 0.1, the <sup>15</sup>N chemical shift of the major state was in the middle of the chirp sweep width and the frequency of the minor state was shifted by 100 Hz (red), 200 Hz (green), 300 Hz (blue), 400 Hz (magenta), 500 Hz (cyan), and 600 Hz (grey). The simulations were done for a chirp pulse with sweep width W = 10 kHz, amplitude B<sub>1</sub> = 2.3 kHz, proportion of the relative length of the apodization ramps T<sub>f</sub>/T<sub>adiab</sub> = 0.2, and total length of the chirp pulse T<sub>adiab</sub> = 80 ms. The transverse CCCR rates were calculated for correlation times τ<sub>0</sub> = 10 and 1 ns for the major and minor states, respectively. The black dotted and dash and dot lines represent the transverse CCCR rate of the major state η<sub>xy</sub><sup>A</sup> and population averaged value p<sub>A</sub>η<sub>xy</sub><sup>A</sup> + p<sub>B</sub>η<sub>xy</sub><sup>B</sup>, respectively</figcaption></figure>

the standard experiment clearly demonstrate the dynamic
character of the first α helix in the excited state. It is well
known from relaxation dispersion and hydrogen/deuterium
exchange experiments (Tollinger et al. 2006; Schanda
et al. 2008) that this α helix is unfolded in the excited state,
however, to the best of our knowledge, no proof had been
obtained so far of any enhanced ps–ns motions of this α
helix in the excited state. Interestingly, we find that helix
α₂, which is expected to be in fast exchange between
unfolded and folded states (both roughly 50 % populated;
Schanda et al. 2008), does not show a significant change of
its average ps–ns dynamics in the excited state. Not sur-
prisingly, helix α₃ was found to be as rigid in the excited
state (where it remains mostly α-helical) as in the folded
state.

Several simulations were performed with various pop-
ulations, chemical shift differences, and rates of exchange.
An example is shown in Fig. 7. Generally, it was found that
our adiabatic experiment (scheme 1 in Fig. 1) provides a
good estimate of the population weighted CCCR rate
unless the exchange rate is very slow (k<sub>ex</sub>~1/T<sub>adiab</sub>), in
which case the system does not 'hop' frequently enough to
effectively average η<sub>xy</sub>. The standard experiment provides
a value close to η<sub>xy</sub> of the major state, unless the process is
fast enough to reach the fast exchange regime.
For comparison of the sensitivity, the relaxation delay $T$ of the standard experiment was set to $T = \alpha T_{\text{adiab}}$. A comparison of the intensities are shown in Supplementary information. Generally, the intensities of residues that are not affected by slow exchange should always be larger in spectra measured by the standard experiment because its total relaxation delay $T = \alpha T_{\text{adiab}}$ is shorter than the total relaxation delay $T$ of scheme 1 in Fig. 1. The difference is most pronounced for flexible residues, which typically have intense signals. Because the relaxation of the two-spin order $2N_zH_z$ is faster than the relaxation of the $N_z$ polarization, the most dramatic attenuation affects the intensities $I_{P,Q}$ where both $P$ and $Q$ represent $2N_zH_z$, which usually show the most intense signals in symmetrical reconversion. Therefore, a drop of the sensitivity compared to the standard experiment observed for almost all peaks in this spectrum is not critical. However, if slow exchange contributes significantly to relaxation the adiabatic method becomes the method of choice since the irradiation diminishes the dominant relaxation mechanism. The precision of the symmetrical reconversion method is mostly limited by the weak intensities of $I_{P,Q}$ and $I_{Q,P}$, where $P$ stands for $N_z$ and $Q$ stands for $2N_zH_z$. The intensity was improved for 46 out of 68 peaks and only 3 peaks exhibit an intensity decrease larger than 10 % in the spectra, where $N_z$ and $2N_zH_z$ terms were selected before and after the relaxation period. A similar proportion of the peaks (43 out of 68) with an increased intensity was found in spectra corresponding to the complementary cross-relaxation pathway ($2N_zH_z$ term selected before the relaxation period and $N_z$ after the relaxation period). A decrease of intensity larger than 10 % was found only for 9 peaks in this case. Finally, the intensities of 97 % peaks were increased in spectra where $N_z$ magnetization was selected both before and after the relaxation period.

Note that $\eta_{xy}$ and $\eta_z$ can be also determined from two experiments with adiabatic pulses without the flanking $T_z$ periods. In that case the relaxation during $T_z$ is avoided, while the advantage of the suppression of exchange is preserved. A measurement of $\eta_z$ by standard experiments and a linear combination of $\eta_{xy}$ and $\eta_z$ by adiabatic experiments is another possibility which keeps the advantage of an independent and direct determination of $\eta_z$.

## Conclusions

We introduced a series of experiments to quantify cross-correlated cross-relaxation CCCR rates under adiabatically swept pulses. These experiments allow one to extend applications of such pulses to the study of fast (ps–ns) motions in proteins. The adiabatic experiment for the determination of transverse CCCR rates provides results that are consistent with standard experiments. The

adiabatic sweeps ensure almost uniform conditions over the whole range of amide chemical shifts even at the highest fields accessible. The method was demonstrated for the KIX domain and human ubiquitin measured at room and low temperatures to mimic different dynamic behaviours. Two different advantages of our approach could be useful in the presence of slow chemical exchange. First, the sensitivity is improved thanks to the suppression of a loss of phase coherence due to slow exchange. Second, the measured transverse CCCR rates represent ensemble-averaged values over all states that undergo slow exchange except if the exchange rates are very small. Finally, the experiment can be modified to yield spectral density values at zero frequency $J(0)$ directly free of any bias due to exchange. In addition, a single experiment is sufficient for the determination of $J(0)$. Precise values $J(0)$ can be extracted even when the standard analysis is limited by the precision of the measurement of longitudinal CCCR rates.

The methods with an adiabatically swept pulse are best suited for deuterated protein samples because interactions with further spins may introduce significant systematic errors. In addition, the presented adiabatic methods are not well suited for rapidly moving proteins undergoing fast amide proton exchange ($>4 \text{ s}^{-1}$) which is typical of intrinsically disordered proteins close to physiological conditions. In both these cases the standard experiment for a measurement of CCCR rates (Pelupessy et al. 2003, 2007) or the variant of the presented method for a direct $J(0)$ measurement without adiabatically swept pulse should be used.

**Acknowledgments** This research has received funding from the European Research Council (ERC) under the European Communitys Seventh Framework Programme (FP7/2007–2013), ERC Grant Agreement 279519 (2F4BIODYN). This work was supported by the Scientific Exchange Programme between Switzerland and the New Member States of the EU [Sciex Foundation Project Code 12.251] (P. K., L. Ž, G. B.), SYLICA [Project Code FP7-REGPOT-2011-1] (P. K., L. Ž), and the Austrian Science Fund FWF (P26849 to M.T.). Special thanks are addressed to Shahid Khan for his preparation of a ubiquitin sample and to Vladimír Sklenář for his advice and discussions.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## Appendix

### Derivation of errors due to erroneous $\alpha$ parameter

* $\alpha$ the true parameter

* $\alpha'$ the assumed parameter
$\eta'_{xy}$ the result of the experiment performed on the assumption that $\alpha'$ is correct

$J'(0)$ the result of the experiment performed on the assumption that $\alpha'$ is correct

The experiment for the determination of the CCCR rate $\eta_{xy}$ provides $\eta'_{xy}$:

$$ \eta'_{\text{ave}} = \frac{(1 - \alpha')\eta_z - (\alpha\eta_{xy} + (1 - \alpha)\eta_z)}{1 + 1 - \alpha'} \quad (32) $$

and $\eta'_{xy}$ is calculated according to Eq. 25:

$$ \eta'_{xy} = \eta_{\text{ave}}^{xy} \frac{\alpha' - 2}{\alpha'} \quad (33) $$

Then, a relative error of $\eta'_{xy}$ with respect to the correct value is given by:

$$ \frac{\eta_{xy} - \eta'_{xy}}{\eta_{xy}} = \left(1 - \frac{\alpha}{\alpha'}\right) \left(1 - \frac{\eta_z}{\eta_{xy}}\right) \quad (34) $$

A similar approach is used for the derivation of an error $J(0)$:

$$ \eta'_{\text{ave}} = \frac{(1 - 0.5\alpha')\eta_z - (\alpha\eta_{xy} + (1 - \alpha)\eta_z)}{1 + 1 - 0.5\alpha'} \quad (35) $$

providing $J'(0)$ calculated according to Eq. 26:

$$ J'(0) = \eta'_{\text{ave}} \frac{\alpha' - 4}{16C_{cd}\alpha'} \quad (36) $$

Finally, a relative error with respect to the correct value $J(0) = (\eta_{xy} - 0.5\eta_z)/8C_{cd}$ is calculated as:

$$ \frac{J(0) - J'(0)}{J(0)} = \left(1 - \frac{\alpha}{\alpha'}\right) \left(1 - \frac{\eta_z}{2\eta_{xy} - \eta_z}\right) \quad (37) $$

The expressions in Eqs. 28 and 29 can be derived simply by substitution of $\eta_{xy}$ and $\eta_z$ according to Eqs. 6 and 7 into Eqs. 34 and 37, respectively.

## References

* Abragam A (1961) The principles of nuclear magnetism. Clarendon Press, Oxford

* Auer R, Kloiber K, Vavrinska A, Geist L, Coudevylle N, Konrat R (2010) Pharmacophore mapping via cross-relaxation during adiabatic fast passage. J Am Chem Soc 132:1480–1481. doi:10.1021/ja910098s

* Auer R, Tollinger M, Kuprov I, Konrat R, Kloiber K (2011) Mathematical treatment of adiabatic fast passage pulses for the computation of nuclear spin relaxation rates in proteins with conformational exchange. J Biomol NMR 51:35–47. doi:10.1007/s10858-011-9539-8

* Böhlen JM, Bodenhausen G (1993) Experimental aspects of chirp NMR spectroscopy. J Magn Reson Ser A 102:293–301. doi:10.1006/jmra.1993.1107

* Böhlen JM, Burghardt I, Rey M, Bodenhausen G (1990) Frequency modulated “Chirp” pulses for broadband inversion recovery in magnetic resonance. J Magn Reson 90:183–191. doi:10.1016/0022-2364(90)90377-L

* Croke RL, Sallum CO, Watson E, Watt ED, Alexandrescu AT (2005) Hydrogen exchange of monomeric $\alpha$-synuclein shows unfolded structure persists at physiological temperature and is independent of molecular crowding in *Escherichia coli*. Protein Sci. 17:1434–1445. doi:10.1110/ps.033803.107

* Delaglio F, Grzesiek S, Vuister GW, Zhu G, Pfeifer J, Bax A (1995) NMRPipe: a multidimensional spectral processing system based on UNIX pipes. J Biomol NMR 6:277–293. doi:10.1007/BF00197809

* Eaton JW, Bateman D, Hauberg S (2009) GNU Octave version 3.0.1 manual: a high-level interactive language for numerical computations. CreateSpace Independent Publishing Platform, Seattle

* Farrow NA, Zhang O, Szabo A, Torchia DA, Kay LE (1995) Spectral density function mapping using $^{15}$N relaxation data exclusively. J Biomol NMR 6:153–162

* Ferrage F, Cowburn D, Ghose R (2009) Accurate sampling of high frequency motions in proteins by steady-state $^{15}$N$^{1}$H nuclear Overhauser effect measurements in the presence of cross-correlated relaxation. J Am Chem Soc 131:6048–6049. doi:10.1021/ja809526q

* Ghose R, Eykyn TR, Bodenhausen G (1999) Average Liouvillian theory revisited: cross-correlated relaxation between chemical shift anisotropy and dipolar couplings in the rotating frame in nuclear magnetic resonance. Mol Phys 96:1281–1288

* Ghose R (2000) Average Liouvillian theory in nuclear magnetic resonance - Principles, properties, and applications. Concepts Magn Reson 12:152–172

* Goddard TD, Kneller DG (2006) Sparky. UCSF, San Francisco

* Halle B (2009) The physical basis of model-free analysis of NMR relaxation data from proteins and complex fluids. J Chem Phys 131(224):507. doi:10.1063/1.3269991

* Halle B, Andersson T, Forsén S, Lindman B (1981) Protein hydration from water oxygen-17 magnetic relaxation. J Am Chem Soc 103:500–508. doi:10.1021/ja00393a004

* Hogben HJ, Krzystyniak M, Charnock GTP, Hore PJ, Kuprov I (2011) Spinach A software library for simulation of spin dynamics in large spin systems. J Magn Reson 208:179–194. doi:10.1016/j.jmr.2010.11.008

* Ishima R, Nagayama K (1995a) Protein backbone dynamics revealed by quasi spectral density function analysis of amide N-15 nuclei. Biochemistry 34:3162–3171. doi:10.1021/bi00010a005

* Ishima R, Nagayama K (1995b) Quasi-spectral-density function analysis for nitrogen-15 nuclei in proteins. J Magn Reson Ser B 108:73–76. doi:10.1006/jmrb.1995.1104

* Kadeřávek P, Zapletal V, Rabatinová A, Krásný L, Sklenář V, Žídek L (2014) Spectral density mapping protocols for analysis of molecular motions in disordered proteins. J Biomol NMR 58:193–207. doi:10.1007/s10858-014-9816-4

* Konrat R, Tollinger M (1999) Heteronuclear relaxation in time-dependent spin systems: $^{15}$N-T$_{1}$$_{p}$ dispersion during adiabatic fast passage. J Biomol NMR 13:213–221. doi:10.1023/A:1008324721356

* Korzhnev DM, Billeter M, Arseniev AS, Orekhov VY (2001) NMR studies of Brownian tumbling and internal motions in proteins. Prog Nucl Magn Reson Spectrosc 38:197–266. doi:10.1016/S0079-6565(00)00028-5

* Kroenke CD, Loria JP, Lee LK, Rance M, Palmer AG (1998) Longitudinal and transverse $^{1}$H-$^{15}$N dipolar/$^{15}$N chemical shift anisotropy relaxation interference: Unambiguous determination of rotational diffusion tensors and chemical exchange effects in
biological macromolecules. J Am Chem Soc 120:7905-7915.
doi:10.1021/ja980832l

Kupče E, Freeman R (1995) Adiabatic pulses for wide-band inversion and broad-band decoupling. J Magn Reson Ser A 115:273-276.
doi:10.1006/jmra.1995.1179

Lakomek NA, Ying J, Bax A (2012) Measurement of ¹⁵N relaxation rates in perdeuterated proteins by TROSY-based methods. J Biomol NMR 53:209-221. doi:10.1007/s10858-012-9626-5

Levitt MH, Bari LD (1992) Steady state in magnetic-resonance pulse experiments. Phys Rev Lett 69:3124-3127. doi:10.1103/PhysRevLett.69.3124

Lipari G, Szabo A (1982a) Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules 1. Theory and range of validity. J Am Chem Soc 104:4546-4559.
doi:10.1021/ja00381a009

Lipari G, Szabo A (1982b) Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules 2. Analysis of experimental results. J Am Chem Soc 104:4559-4570. doi:10.1021/ja00381a010

Mangia S, Traaseth NJ, Veglia G, Garwood M, Michaeli S (2010) Probing slow protein dynamics by adiabatic $R_{1\rho}$ and $R_{2\rho}$ NMR experiments. J Am Chem Soc 132:9979-9981. doi:10.1021/ja1038787

MATLAB and statistics toolbox release R2014a, The MathWorks, Inc, Natick, Massachusetts, United States

Mulder FAA, Mittermaier A, Hon B, Dahlquist FW (2001) Studying excited states of proteins by NMR spectroscopy. Nat Struct Biol 8:932-935. doi:10.1038/nsb1101-932

Palmer AG III, Kroenke CD, Loria JP (2001) Nuclear magnetic resonance methods for quantifying microsecond-to-millissecond motions in biological macromolecules. Methods Enzymol 339:204-238

Pelupessy P, Espallargas GM, Bodenhausen G (2003) Symmetrical reconversion: measuring cross-correlation rates with enhanced accuracy. J Magn Res 161:258-264. doi:10.1016/S1090-7807(02)00190-8

Pelupessy P, Ferrage F, Bodenhausen G (2007) Accurate measurement of longitudinal cross-relaxation rates in nuclear magnetic resonance. J Chem Phys 126:134508. doi:10.1063/1.2715583

Peng JW, Wagner G (1992a) Mapping of spectral density-functions using heteronuclear NMR relaxation measurements. J Magn Reson 98:308-332. doi:10.1016/0022-2364(92)90135-T

Peng JW, Wagner G (1992b) Mapping of the spectral densities of N-H bond motions in Eglin c using heteronuclear relaxation experiments. Biochemistry 31:8571-8586. doi:10.1021/bi00151a027

Redfield AG (1965) The theory of relaxation processes. Adv Magn Reson 1:1-32

Schanda P, Brutscher B, Konrat R, Tollinger M (2008) Folding of the KIX domain: characterization of the equilibrium analog of a folding intermediate using ¹⁵N/¹³C relaxation dispersion and fast ¹H/²H amide exchange NMR spectroscopy. J Mol Biol 380:726-741. doi:10.1016/j.jmb.2008.05.040

Tollinger M, Kloiber K, Agoston B, Dorigoni C, Lichtenecker R, Schmid W, Konrat R (2006) An isolated helix persists in a sparsely populated form of KIX under native conditions. Biochemistry 45:8885-8893. doi:10.1021/bi0607305

Valentine ER, Ferrage F, Massi F, Cowburn D, Palmer AG III (2007) Joint composite-rotation adiabatic-sweep isotope filtration. J Biomol NMR 38:11-22. doi:10.1007/s10858-006-9131-9

Wangsness R, Bloch F (1953) The dynamical theory of nuclear induction. Phys Rev Lett 89:728-739. doi:10.1103/PhysRev.89.728

Wolfram Research Inc. (2012) Mathematica