Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Biomedical EPR Part-B Methodology Instrumentation and Dynamics - Sandra R. Eaton

.pdf
Скачиваний:
79
Добавлен:
10.08.2013
Размер:
11.5 Mб
Скачать

322

DEREK MARSH ET AL.

Figure 4. Model axial ST-EPR powder lineshapes (heavy lines) calculated from Eq. 19, together with Eqs. 18, 19, 21 and 22, for increasing rates of rotational diffusion. From top to bottom, 0.005, 0.01, 0.02 and 0.04; in all cases and the integrated intensity is normalised to unity in the absence of rotational diffusion [ordinate in units of Inset gives the diagnostic ST-EPR lineheight ratio measured at point P', relative to the turning point P at the right-hand extremum of the spectrum (see Marsh and Horváth, 1992a).

The above analysis suggests that spectral lineheight ratios and integrated intensities, R, may have the following dependence on (cf. Eq. 19):

where is the value of R in the absence of rotational diffusion, a and b are constants to be fitted that depend only on intrinsic spectral parameters, and the ratio a/b is effectively related to the orientational degeneracy parameter

SATURATION TRANSFER SPECTROSCOPY

323

at corresponding to the diagnostic spectral position P'. Equation 23 rather well describes the dependence on rotational correlation time of the diagnostic line height ratios and intensities of the saturation transfer EPR spectra from spin-labelled haemoglobin in glycerol-water mixtures (Marsh and Horváth, 1992a). It therefore can be used to give the following simple expression for the correlation time calibrations of the experimental ST-EPR spectra:

where values of the experimental calibration constants, k, and b, for the different diagnostic spectral parameters are given in Table 1. This is a much more readily accessible form for the calibrations of rotational correlation time than hitherto was presented (e.g., Horváth and Marsh, 1988) and has the additional advantage of reflecting directly the underlying spectral diffusion process.

In principle, slow rotational diffusion also may be studied from the power saturation of the conventional EPR spectra. Squier and Thomas (1986) have done this in terms of saturation factors determined from the ratio of the integrated intensities of the conventional first-derivative EPR spectra recorded at low (subsaturating) and high (partially saturating) microwave powers. Calibrations for this method that are equivalent to those given for ST-EPR by Eq. 24 can be found in Marsh (1995).

324

DEREK MARSH ET AL.

6.APPLICATIONS: SLOW ROTATION

This section gives examples of the classical application of saturation transfer EPR to the study of slow rotational diffusion on the submillisecond timescale. Applications are chosen to illustrate specific aspects of rotational diffusion in membranes such as the effects of protein concentration, aqueous viscosity, hydrophobic matching and anisotropic rotation, rather than giving a comprehensive review.

6.1Dependence on Protein Density

A systematic study of the dependence of the ST-EPR rotational correlation times on lipid/protein ratio, LP, was undertaken by Fajer et al. (1989) with spin-labelled cytochrome c oxidase reconstituted in bilayer membranes of dimyristoyl phosphatidylcholine. Fig. 5 gives the effective rotational relaxation rate, deduced from the central ST-EPR lineheight ratio, C'/C, as a function of lipid/protein ratio in the reconstituted membranes. Protein rotational diffusion is drastically reduced in gel-phase membranes at 1°C where the lipid chains are largely frozen. Rotation is much faster in fluid-phase membranes, with effective correlation times in the tens of microsecond regime. The rate at which cytochrome oxidase rotates in the membrane decreases (i.e., the correlation time increases) progressively with increasing protein packing density.

The hindering of cytochrome oxidase rotation by protein crowding can be described with a simple collisional model based on random proteinprotein contacts. The observed diffusion coefficient, is a statistical average of that for freely rotating species that do not experience any influence from other proteins and that for highly hindered species that have other proteins immediately adjacent

where denotes the probability for the freely rotating species. This model assumes that the lifetime of protein-protein contacts is shorter than the rotation period. For a translational diffusion coefficient of this is likely to be the case (Fajer et al., 1989). The probability is obtained from a lattice model designed to calculate the frequency of lipid-protein contacts in random dispersions (Hoffmann et al., 1981). Each lipid occupies one lattice site and each protein occupies R lattice sites, where for the cytochrome c oxidase monomer (Deatherage et al., 1982) is then the probability that all N lattice sites at the protein perimeter are occupied by lipid molecules:

SATURATION TRANSFER SPECTROSCOPY

325

where LP is the lipid/protein mole ratio (see Fig. 5).

Figure 5. Dependence of the rotational relaxation rate, on lipid/protein molar ratio for cytochrome c oxidase in membranes of dimyristoyl phosphatidylcholine at 30°C (fluid phase, and 1°C (gel phase, Solid line is a non-linear least squares fit of the random collision model (Eqs. 25 and 26) to the data at 30°C with fixed R = 27 (see Fajer et al., 1989). The inset shows the lattice model used to calculate the probability, that a protein (given by the hexagons that occupy R lipid lattice sites) does not contact any other protein. The probability that a lipid (circle) does not occupy a site (e.g., the asterisked position) that is adjacent to a protein is given by R/(LP+R) where LP is the lipid/protein mole ratio.

A non-linear least squares fit with R fixed and N as the parameter to be optimised is given by the solid line in Fig. 5. This depicts the time-averaged protein-protein interactions taking place in the fluid phase. The fitted value of N = 18 corresponds to 36 first-shell lipid sites at the perimeter of the protein, in both bilayer halves of the membrane. This is somewhat smaller than the number of boundary lipids found by EPR measurements with spinlabelled lipids (Knowles et al., 1979). Presumably, the latter reflects the invaginated nature of the intramembranous surface of the protein, which is effectively smoothed when considering protein-protein contacts.

326

DEREK MARSH ET AL.

The effect of lipid/protein ratio on rotational mobility, therefore, can be reasonably described in terms of random protein collisions that occur with increasing protein density. An alternative model of a heterogeneous population of rotating species, with a varying proportion of higher oligomers is unable to explain the dependence on lipid/protein ratio (Fajer et al., 1989). Also, the increase in membrane viscosity with increasing protein concentration is predicted to be insufficient to account for the changes in rotation rate.

6.2Hydrophobic Matching

Rhodopsin has been reconstituted in bilayer membranes formed from phosphatidylcholines with different acyl chainlengths (Ryba and Marsh, 1992). Rotational diffusion of the spin-labelled protein recorded by ST-EPR was used to follow the extent of dispersal of the protein in the different chainlength lipids. A similar study with equivalent results, was undertaken earlier by Kusumi and Hyde (1982).

Figure 6. Lipid chainlength (n) dependence of the effective rotational correlation times of maleimide spin-labelled rhodopsin in membranes of different saturated diacyl phosphatidylcholines at a lipid/protein ratio of 60:1 (mol/mol). Effective correlation times are deduced from the low-field L"/L (+45° hatching) and highfield, H"/H (-45° hatching) ST-EPR diagnostic lineheight ratios (see Ryba and Marsh, 1992). Matching/mismatching of the hydrophobic lengths of lipid and protein are indicated schematically for lipid chainlengths n = 12, 15 and 17.

SATURATION TRANSFER SPECTROSCOPY

327

Figure 6 gives the effective rotational correlation times of rhodopsin in the lipids of different chainlength that are deduced from the low-field and high-field regions of the ST-EPR spectra. These measurements are all made in the fluid membrane state at equivalent temperatures above the chainmelting transition for the different lipids. Mostly, the correlation times for L"/L and H"/H are comparable, which is expected because both are reflecting the same anisotropic rotation. The longest effective rotational correlation times are obtained from recombinants with dilauroyl phosphatidylcholine, with a steep decrease on increasing the lipid chainlength, through a minimum at chainlengths of C14 to C15, and a subsequent rise on increasing the lipid chainlength to dipalmitoyl and diheptadecanoyl phosphatidylcholines, and

The pronounced increases in rotational correlation time for rhodopsin in the long and short chainlength lipids can be attributed to protein aggregation. This is driven by hydrophobic mismatch in both cases, as illustrated diagrammatically in the inset to Fig. 6. A lipid that is too short exposes part of the hydrophobic domain of the protein to a polar environment. A lipid that is too long forces contact of the hydrophobic lipid chains with polar groups on the protein. In each case, these energetically unfavourable interactions are alleviated by segregation of the protein from the lipids.

A topic of equal interest is the oligomer state of rhodopsin in the C15 chainlength lipids for which hydrophobic matching is best. The effective rotational correlation time, deduced from the ST-EPR spectra using calibrations from isotropic solutions is related to the diffusion coefficient, for uniaxial rotation by (Robinson and Dalton, 1980; Marsh and Horváth, 1989):

where is the orientation of the spin-label z-axis relative to the membrane normal, and The rotational diffusion coefficient is related to the cross-sectional dimensions, and and the intramembrane height, of the rotating species by the Stokes-Einstein equation (see e.g., Marsh and Horváth, 1989):

where is the rotational frictional coefficient, is the effective intramembrane viscosity and is a shape factor that depends weakly

328

DEREK MARSH ET AL.

on the asymmetry for (see also following section 6.3). The corresponding true rotational correlation time is: As usual, is Boltzmann’s constant and T is the absolute temperature.

The effective rotational correlation time of rhodopsin in dipentadecanoyl phosphatidylcholine, is (Fig. 6); an interpolated value for the membrane thickness is (Tardieu, 1972); and the effective membrane viscosity is in the region of (Cherry and Godfrey, 1981). This yields a value of for the intramembranous diameter of the rotating species. This is an upper estimate because it is assumed that and Cross-sectional dimensions of the dimer of frog rhodopsin are (Corless et al., 1982). Thus rhodopsin is most likely a monomer in as it is in rod outer segment disc membranes (Downer, 1985).

6.3Dependence on Extramembrane Viscosity

As is well known, the intramembrane viscosity that characterises the torque on large integral membrane proteins is much greater than the lipid microviscosity that is determined with small probe molecules by applying

the Debye equation (Cherry and Godfrey,

1981). Therefore, because the

effective viscosity in the membrane

is so much higher than that

of water, the rotational diffusion coefficients of membrane proteins are normally determined solely by the intramembranous sections of the proteins (cf., Eq. 28). Only if the extramembrane viscosity is increased considerably, e.g., by addition of sucrose or glycerol, does the rotational diffusion coefficient become dependent on the dimensions of the extramembranous sections of the protein.

The frictional torques exerted on the separate sections of the protein are additive, therefore so are also the individual contributions to the overall frictional coefficient:

where the frictional coefficients of right circular cylinders with volumes, equal to those of the different sections, i, of the protein are given by:

and the corresponding shape factors are:

SATURATION TRANSFER SPECTROSCOPY

329

where are the elliptical semi-axes of the different protein cross-sections (see Fig. 7 inset). The rotational correlation time of the protein therefore is given from Eqs. 28-30 by:

Specifically, the dependence on external viscosity is given by:

where the summation (as indicated by the prime) now extends only over the extramembrane sections of the protein, and is the rotational correlation time given by Eq. 28 above, when external viscosity can be neglected.

Figure 7 gives the dependence of the effective rotational correlation time of membranous spin-labelled Na,K-ATPase on viscosity of the external glycerol-containing medium. From the slope and intercept of the linear viscosity dependence, together with Eqs. 28 and 33, it is concluded that 5070% of the Na,K-ATPase protein is external to the membrane (Esmann et al., 1994). This conclusion obtained from hydrodynamics is consistent with the results of low-resolution structural studies on this protein (Maunsbach et al., 1989). Fig. 7 also shows that, with polyethylene glycol solutions, a pronouncedly non-linear dependence on the viscosity, is found that is much larger than the viscosity dependence obtained with glycerol solutions. This greater effect of polyethylene glycol undoubtedly corresponds to a dehydration-induced aggregation of the membrane proteins that may be related to the ability of polyethylene glycol to induce membrane fusion. The rotational correlation time reached at 50% polyethylene glycol corresponds to a degree of aggregation of the membrane proteins between two and five, depending on whether the ethylene glycol polymer is excluded from the membrane surface region (Esmann et al., 1994). Clearing of proteins from areas of apposing membrane, by aggregation, is a prerequisite for the close approach of the lipid bilayers that is needed for effective membrane fusion.

330

DEREK MARSH ET AL.

Figure 7. Dependence on extramembrane viscosity, of the effective rotational correlation time, deduced from the high-field diagnostic lineheight ratio in the ST-EPR spectra of Na,K-ATPase spin-labelled with a chloromercuri-reagent. The solid line is a linear regression for the glycerol data with intercept and gradient The dashed line for the PEG data is simply to guide the eye (see Esmann et al., 1994). The inset illustrates the different dimensions of the extramembrane and intramembrane sections of the protein that experience viscosities and respectively.

6.4Anisotropic Rotational Diffusion

Rotational correlation times are routinely deduced from experimental STEPR spectra by comparing the diagnostic lineheight ratios in the low-field, central and high-field regions of the spectrum with those obtained from isotropically rotating spin-labelled haemoglobin in solutions of known viscosity (see Table 1). The outer lineheight ratios, L"/L and H"/H, are sensitive to rotation of the nitroxide z-axis, via modulation of the hyperfine interaction, and the central lineheight ratio, C'/C, is sensitive to rotation about all three nitroxide axes, via modulation of the g-value anisotropy (see Eq. 7 and Fajer and Marsh, 1983). For anisotropic rotation, the different lineheight ratios will therefore have differential sensitivities, as illustrated in Table 2 for the rotational diffusion of a spin-labelled phospholipid in gelphase lipid bilayer membranes (Marsh, 1980). The nitroxide z-axis is oriented along the lipid long molecular axis for this particular spin probe. At low temperatures the effective correlation times, deduced from

SATURATION TRANSFER SPECTROSCOPY

331

the different lineheight ratios using the isotropic model system for calibration, are all very similar. Above the bilayer pretransition at 25°C, the effective correlation time deduced from the outer lineheight ratios is relatively unchanged, whilst that deduced from the central spectral region decreases abruptly, indicating the onset of rapid anisotropic rotational diffusion about the long axis of the lipid molecule. A differential response of the different lineheight ratios may therefore be used to diagnose anisotropic rotational diffusion (Fajer and Marsh, 1983b). Comparison of the integral of the high-field region of the saturation transfer spectrum with that of the total ST-EPR spectrum may also be similarly used (see Table 1).

Simulations of first harmonic phase-quadrature dispersion ST-EPR spectra by Robinson and Dalton (1980) give some guide to the quantitative interpretation of the effective correlation times in terms of the true rotational diffusion parameters. It was found that the effective correlation times deduced from the low-field and high-field regions of the spectrum were very similar, independent of the degree of anisotropy of the motion (cf. Table 2). If the anisotropy of the rotation is great enough and the rotational rates are slow, the dependence of the effective correlation times on the orientation, of the nitroxide z-axis with respect to the rotational diffusion axis, is given

by Eq. 27

that was introduced above. Precise measurements thus require

knowledge

of the

orientation,

However, some estimate of whether is

close to 0° or close to 90°

may be obtained from the relative sizes of

and

or by

comparing the high-field and total spectral

integrals. In addition, it may be possible to discriminate between protein monomer and oligomer formation on the basis of Eqs. 28-31 without accurate knowledge of

7. NONLINEAR EPR DISPLAYS

As will be seen from the examples given later in Section 9, the impetus for exploring the direct sensitivity of non-linear EPR to spin-lattice relaxation and cross-relaxation processes is in the study of paramagnetic