Spin-label Order Parameter Calibrations for Slow Motion.

Calibrations are given to extract orientation order parameters from pseudo-powder electron paramagnetic resonance line shapes of 14N-nitroxide spin labels undergoing slow rotational diffusion. The nitroxide z-axis is assumed parallel to the long molecular axis. Stochastic-Liouville simulations of slow-motion 9.4-GHz spectra for molecular ordering with a Maier-Saupe orientation potential reveal a linear dependence of the splittings, [Formula: see text] and [Formula: see text], of the outer and inner peaks on order parameter [Formula: see text] that depends on the diffusion coefficient [Formula: see text] which characterizes fluctuations of the long molecular axis. This results in empirical expressions for order parameter and isotropic hyperfine coupling: [Formula: see text] and [Formula: see text], respectively. Values of the calibration constants [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text] and [Formula: see text] are given for different values of [Formula: see text] in fast and slow motional regimes. The calibrations are relatively insensitive to anisotropy of rotational diffusion [Formula: see text], and corrections are less significant for the isotropic hyperfine coupling than for the order parameter.

Introduction

Spin-labelled molecules in partially ordered environments are restricted in their rotational amplitude. This is a common situation in liquid crystals or lipid membranes; also for spin labels anchored to a macromolecule, and in ligand binding sites or occlusion cavities. The degree of ordering is expressed by an order parameter that is related to the rotational amplitude and can be determined directly from the electron paramagnetic resonance (EPR) spectrum.

When the nitroxide z-axis lies close to the long molecular axis of a spin-labelled molecule, we can extract orientation order parameters from the EPR powder patterns of non-aligned samples, because spectral anisotropy is then largest [1-3]. Notable practical examples are lipid chains incorporating the DOXYL nitroxide (see e.g., [4-6]), and helical peptides containing the TOAC nitroxide residue [7-9]. Such situations also arise frequently in site-directed spin-labelling with MTSSL [10-12], and also for rigidly spin-labelled DNA bases [13].

For rotational motion that is fast relative to the spectral anisotropies, the 14N hyperfine tensor is partially averaged to give axial principal values: where and . The order parameter is given by:where is the instantaneous angle that the nitroxide z-axis makes with the director for ordering, and angular brackets indicate a time or ensemble average over the orientational distribution in .

Values for the order parameter are obtained from the splitting of the outer peaks in an experimental powder pattern, and from that of the inner peaks . These peaks in the powder spectrum correspond to turning points in the angular dependence of spectra from macroscopically aligned samples. For fast motion, is rigorously equal to [1, 14]. But to obtain we must make a correction to , because of spectral overlap in the central region of the powder pattern [1, 2, 15, 16]. When the motion is slow on the nitroxide EPR timescale, however, both outer and inner peaks shift relative to the positions specified by the motionally averaged tensor elements and . We expect this because slow motion causes shifts in the outer extrema of the powder pattern, even for unrestricted isotropic rotational diffusion [17, 18]. The main purpose here is to obtain appropriate corrections when the spectra are in the slow-motion regime for nitroxide EPR.

Results

Figure 1 shows 9.4-GHz EPR spectra for a 14N-nitroxide with fixed order parameter , and different values of the rotational diffusion coefficient that extend from the slow to fast motional regimes. Powder spectra from non-aligned samples are simulated with the stochastic-Liouville equation for Brownian diffusion [19], as implemented in EasySpin 5.1.9 [20], with molecular ordering specified by a Maier–Saupe orientational potential: [21]. In terms of the latter, the order parameter becomes:from the Boltzmann distribution, where . The nitroxide z-axis lies along the long molecular axis. The diffusion coefficient for excursions of the molecular long axis (and simultaneously the nitroxide z-axis) is , and that for rotation about the same axis is ten times faster. As decreases below , the splitting of the outer peaks in the powder pattern increases and that of the inner peaks decreases, as expected in the slow-motion region. The spectral line shapes also give us some guide as to whether a particular experimental powder pattern contains slow-motional components.

Fig. 1

Stochastic-Liouville simulations of 14N-nitroxide 9.4-GHz EPR powder spectra, for different values of the rotational diffusion coefficient , as indicated, and fixed order parameter . Vertical ticks show the positions of the outer and inner hyperfine peaks, which have splittings and respectively. Hyperfine and g-tensors are: and . Residual Gaussian broadening, FWHM = 0.18 mT

In the fast regime, the motionally averaged hyperfine constants and depend linearly on the order parameter according to Eqs. 1 and 2. So also do the experimental constants and . In the slow-motion regime, we find from simulations such as those in Fig. 1 that and depend approximately linearly on order parameter over the range (see Fig. 2):

Fig. 2

Dependence on order parameter of the shifts, (squares) and (circles), in and from their motionally averaged values and , because of slow motion and spectral overlap. Spectra are simulated as in Fig. 1, using a Maier–Saupe orientational potential, and fixed diffusion coefficient . Straight lines are linear regressions over the range (see Eqs. 5, 6)

Values of the parameters, , and , , for different values of are obtained by linear regression. Figure 3 shows the dependence of these calibration parameters in Eqs. 5 and 6 on diffusion coefficient. Note that the -axis is logarithmic. We see that the dependence for all four parameters is approximately logarithmic over the range , in the slow-motion regime. This dependence is given by the solid lines. In the fast-motion region , the linear dependence of on order parameter is slight and restricted to the range . The dotted lines in Fig. 3 are fits of the logistic equation: , with logarithmic x-axis. These provide fits over a wider range of diffusion coefficient, but with lower precision. We return to this type of fit later.

Fig. 3

Dependence on rotational diffusion coefficient of linear regression parameters for the line shifts , and , as a function of (see Fig. 2). Solid lines are linear regressions over the range for the logarithmic x-axis. Dotted lines are fits of the logistic equation (see Table 1)

Table 1

Calibration parameters in Eqs. 8 and 13 for order parameter and isotropic hyperfine coupling, for different diffusion coefficients

D_R⊥ (s−1)	s_1 (mT−1)	s_o	f_max	f_min	δa_o (mT)
1.5 × 107	0.529	1.067	1.020	1.980	−0.046
2.0 × 107	0.535	1.050	1.024	1.976	−0.058
2.4 × 107	0.522	0.973	1.019	1.981	−0.043
2.65 × 107	0.512	0.914	1.015	1.985	−0.032
3.0 × 107	0.488	0.805	1.004	1.996	0.002
4.0 × 107	0.440	0.560	0.988	2.012	0.050
5.0 × 107	0.411	0.402	0.982	2.018	0.068
6.0 × 107	0.397	0.315	0.980	2.020	0.073
7.0 × 107	0.378	0.228	0.972	2.028	0.093
8.0 × 107	0.382	0.205	0.982	2.018	0.070
1.0 × 108	0.373	0.140	0.982	2.018	0.069
1.2 × 108	0.369	0.105	0.984	2.016	0.063
1.5 × 108	0.366	0.077	0.986	2.014	0.059
	(0.353	0.044	0.974	2.026	0.087)a
2.0 × 108	0.364	0.051	0.988	2.012	0.053
	(0.357	0.035	0.982	2.018	0.067)a

Empirical logistic fits

aFits with ,

Substituting from Eqs. 5 and 6 into Eq. 3, we get the following calibration for the corrected order parameter:in terms of the measurables and . We therefore write Eq. 7 as:where the calibration parameters and are given for different values of in Table 1. Logistic fits of the dependence on diffusion coefficient (cf. dotted lines in Fig. 3), over the full range , are given in the footnote to Table 1. We get higher precision, but over a more limited range, by using linear regression with logarithmic -axis (cf. solid lines in Fig. 3): over the range . When in the fast regime, over the range , and the corresponding calibration becomes:

For comparison, previous corrections require subtraction of 0.03–0.05 for fast motion [1-3].

We also get corrected values of the isotropic hyperfine coupling, , from Eqs. 5 and 6:

This depends on , which we substitute from Eq. 7 (or 8) to give the following general form:where the calibration parameters , and are given for different values of in Table 1. Using linear regression with logarithmic -axis (cf. Fig. 3), we get the approximate dependence of the calibration parameters: over the range . Again, logistic fits over the full range of are included in the footnote to Table 1. Note that corrections to for slow motion are smaller than those to the order parameter, because the opposite shifts in and compensate. When in the fast regime, the calibration corresponding to Eq. 11 becomes:

For comparison, previous treatments make an additive correction of 0.05–0.09 mT for fast motion [1-3]

When calculating the order parameter, we use the isotropic hyperfine coupling to correct for differences between the hyperfine tensor of the experimental sample and that used to get for Eq. 3. The spin-Hamiltonian tensors used here, with  = 2.74 mT (see Fig. 1), correspond to the MTSSL spin used in site-directed spin labelling (see, e.g., [22]). They correspond reasonably closely with those for DOXYL spin labels, which also have a five-membered ring. To make the correction, we simply scale in Eq. 8 or 11 by the factor: , where comes from Eq. 13 or 17, respectively. This corrects for differences in environmental polarity, and to a first approximation for differences between nitroxide structures (e.g., TOAC spin labels with the six-membered piperidine ring).

The calibration parameters in Table 1 are obtained from simulations for fast axial diffusion: , corresponding to a prolate rigid rotator with axial ratio [23]. For segmental rotation, however, the elements of the diffusion tensor are more likely to be equal. Figure 4 compares spectra simulated for (solid lines) with those for (dashed lines), where the order parameter is now ; in each case the value of is the same. We see some differences in line shape for slow motion, but not in the fast motional regime. Table 2 gives the calibration parameters for order parameter and isotropic hyperfine coupling when . Only in the slow motion regime do these differ from those in Table 1, and mostly the differences are not large. We conclude from this that calibration parameters depend markedly on in the slow motion regime, but little on the rotational anisotropy .

Fig. 4

Stochastic-Liouville simulations of 14N-nitroxide 9.4-GHz EPR powder spectra, for different values of the rotational diffusion coefficient (solid lines) and (dashed lines), as indicated, and fixed order parameter . Simulation parameters otherwise as in Fig. 1

Table 2

Calibration parameters in Eqs. 8 and 13 for order parameter and isotropic hyperfine coupling, for different diffusion coefficients

D_R (s−1)	s_1 (mT−1)	s_o	f_max	f_min	δa_o (mT)
1.5 × 107	0.544	1.078	1.039	1.961	−0.105
2.0 × 107	0.540	1.052	1.031	1.969	−0.081
2.4 × 107	0.527	0.975	1.026	1.975	−0.064
2.65 × 107	0.516	0.916	1.021	1.980	−0.049
3.0 × 107	0.493	0.814	1.010	1.990	−0.016
4.0 × 107	0.446	0.562	0.998	2.002	0.022
5.0 × 107	0.427	0.438	0.995	2.005	0.034
6.0 × 107	0.404	0.328	0.986	2.014	0.058
7.0 × 107	0.380	0.229	0.974	2.026	0.087
8.0 × 107	0.368	0.170	0.970	2.030	0.097
1.0 × 108	0.374	0.144	0.983	2.017	0.066
1.2 × 108	0.369	0.105	0.984	2.016	0.064
1.5 × 108	0.367	0.076	0.987	2.013	0.057
	(0.353	0.044	0.975	2.025	0.086)a
2.0 × 108	0.364	0.051	0.988	2.012	0.053
	(0.357	0.035	0.982	2.018	0.067)a

aFits with ,

Discussion and Conclusions

Experience shows that spin-label spectra from systems exhibiting partial molecular order frequently contain components that are in the slow motional regime for nitroxide EPR. In such cases, we must make corrections to both and , which we obtain from the canonical turning points in pseudo-powder spectra from non-aligned samples, when calculating orientational order parameters and isotropic hyperfine couplings. Because of partial cancellation in the latter case, corrections are more important for than for . A previous treatment with an oversimplified slow-motion model gave no correction for because a near rigid-limit value of was used to emulate a powder distribution [3].

To calculate order parameters for slow motion, we must know the diffusion tensor element that characterizes angular excursions of the molecular long axis (to which the nitroxide z-axis is assumed approximately parallel). This requires simulation results for at least one spectrum in an experimental series. The spectra in Figs. 1 and 4, which correspond to two different order parameters, can be a helpful guide here. We then get order parameter and isotropic coupling constant from Eqs. 8 and 13, respectively, by using the appropriate calibration constants listed in Table 1 (or Table 2). Alternatively, the calibration constants come from Eqs. 9, 10, 14–16 for , or less precisely outside this range from the logistic fits listed at the bottom of Table 1. As we see from Fig. 4, the value of is crucial, whereas the faster rotation around the molecular axis, , is relatively unimportant.

The calibrations given here are for partially averaged powder spectra that display well defined axial anisotropy. For this, the nitroxide z-axis must lie close to the long molecular axis that is being ordered, because this is the principal axis for the N-hyperfine anisotropy. If the nitroxide z-axis is inclined at angle to the molecular axis, we know from the spherical harmonic addition theorem that the order parameter measured by EPR is given by: , where is the order parameter of the long molecular axis, and is a second rank Legendre polynomial. As already mentioned, practical examples of spin labels with the appropriately oriented nitroxide axes are DOXYL-labelled hydrocarbon chains, and TOAC in peptide helices. Complications arise with lipid chains because motion of the long molecular axis often lies in the slow regime [24], but rotational isomerism of the chain segments is much faster [25, 26]. In this case, the effective hyperfine anisotropy is reduced from to by rapid trans-gauche isomerism in the chain (cf. Eqs. 1 and 2). Order-parameter calculations for the slow long-axis motion are then made using this reduced anisotropy. Multifrequency EPR is helpful in separating fast and slow components of spin-label motion. For instance, the slower long-axis motion of lipid chains in phospholipid membranes is driven into the rigid limit at 94-GHz EPR frequency, and the high field spin-label spectrum reflects motional averaging by trans-gauche isomerism alone [27-29].

Although the slow-motion order parameter calibrations require an estimate of the diffusion coefficient, they offer significant improvements over those based purely on motional narrowing theory. Even when a full spectral simulation study is to be attempted, the calibrations given here are a useful starting point to establish initial simulation parameters.

It goes without saying that the present corrections apply only to those pseudo-powder spectra which show a well-defined inner hyperfine splitting, . Although this accounts for a major class of spectra from important spin labels, there are numerous cases when this does not apply, e.g., polymeric liquid crystals and glasses (see for instance [30]).