Side-Chain Dynamics in MUP
A R T I C L E S
1
3
1
2
1
2
C NMR Relaxation Measurements. Samples of [γ ,γ -
In contrast, an equivalent spectrum recorded on [γ ,γ -
1
3
1
1
2
2 2
13
1
1
2
2 2
C2,R,â,γ ,γ ,γ ,γ - H6] valine enriched MUP both alone and
C2,R,â,γ ,γ ,γ ,γ - H6] valine enriched MUP (Figure 1b) is
complexed with 2-methoxy-3-isopropylpyrazine and 2-methoxy-
essentially free from resonance overlap, and all 12 valine methyl
3
7
-isobutylpyrazine were prepared from a single sample at pH
.0 and a protein concentration of 1 mM. A further sample of
groups can be observed and assigned from the work of Abbate
3
3
et al. MUP exhibits significant promiscuity in binding small
hydrophobic ligands, and we have chosen to investigate the
thermodynamics of binding of 2-methoxy-3-isopropylpyrazine
and 2-methoxy-3-isobutylpyrazine on the basis of their global
thermodynamics of binding from isothermal titration calorimetry
experiments (data not shown). These ligands bind to MUP with
Kd’s of ∼1800 nM and ∼300 nM, respectively, at 308 K. Figure
1
3
2
methyl- C, H enriched MUP was prepared at pH 7.0 and a
protein concentration of 0.5 mM. Longitudinal (R1) and
transverse (R2) 13C relaxation rates were determined essentially
as described by Ishima et al.25 All spectra were recorded at a
proton frequency of 600 MHz and a probe temperature of 35
°
C. R1 rates were determined with relaxation delay times of
6, 64, 128, 240, 400, 720, 1120, 1600, 2240, and 2880 ms. R2
rates were determined with relaxation delay times of 16, 32,
8, 64, 96, 128, 160, 192, 240, and 288 ms and an effective
1
3
1
1
1c shows C, H HSQC spectra of complexes of these ligands
1
2 13
1
1
2
2 2
with [γ ,γ - C2,R,â,γ ,γ ,γ ,γ - H6] valine enriched MUP.
4
Significant shift perturbations are observed only for the cor-
γ1
γ2
field strength of 2 kHz. Relaxation data were fit to a single
exponential decay function I ) Io exp(-tR), where Io is the
initial resonance intensity, t is the relaxation delay time, and R
is the relaxation rate (R1 or R2). Isotropic rotational tumbling
was assumed on the basis of previous backbone N relaxation
measurements for this protein.17 Thus, relaxation data were fitted
to the Lipari-Szabo model-free spectral density3 of the form:
relations from Val 82 C and C , as anticipated since Val 82
3
7
is located within the binding pocket of MUP. Stereospecific
γ1
γ2
assignments of Val 82 C and C in these complexes were
15
confirmed by three-dimensional NOE measurements using N,
C enriched MUP and the relevant unlabeled ligands (data not
shown). In both complexes, all correlations are well-resolved,
in contrast to complexes with methyl- C, H enriched MUP
data not shown). Under these circumstances, measurement of
15
13
38
5
13
2
(
2
2
2
m
2
2 2
i
J(ω) ) S τ /(1 + ω τ ) + (1 - S )τ /(1 + ω τ ) (1)
m
i
13
C longitudinal and transverse relaxation rates (R1 and R2) can
where τi-1 ) τM-1 + τe-1, τM is the overall molecular tumbling
be undertaken in each sample at high sensitivity and with
minimal interference due to resonance overlap. Typical relax-
ation curves are shown in Figure 2. Values of Saxis were
correlation time, and τe is the effective internal correlation time.
This equation is valid for fast internal motions, τe , τM, under
which conditions the order parameter for methyl CH dipolar
2
derived from R1 and R2 data for each methyl group using the
2
axis
2
2
axis
3
5
2
model-free spectral density function given in the General
Methods section in the Supporting Information and are listed
in Table 1. From these data, it is immediately apparent that
relaxation is given by S ) S [P2(cos θH)] . S is the order
parameter of the methyl rotation axis and θH is the angle made
by this axis and the CH bond vector. The MODELFREE
program was used for all these calculations,36 kindly supplied
by Prof. Art Palmer (Columbia University) and modified in-
2
there are only minor changes in Saxis for any of the six valine
residues in MUP upon binding of either 2-methoxy-3-isopro-
pylpyrazine or 2-methoxy-3-isobutylpyrazine, most of which
13
house to include relaxation of C by the attached deuterons in
2
13
CHD2 isotopomers. A global rotational correlation time of 8.57
are within experimental error. An exception is Saxis for Val-70
γ
ns was used for these calculations, according to the previously
reported value.17 Errors in S values were estimated by Monte
Carlo simulations as implemented in the MODELFREE pack-
age. In all cases, good fits were obtained to the above spectral
density function without a contribution from exchange-broaden-
ing Rex. Although valine residues are perdeuterated at nonmethyl
positions, the protein is otherwise protonated, and dipolar
relaxation due to these “external” protons is not negligible for
C 2, which appears to increase dramatically on binding 2-meth-
2
oxy-3-isobutylpyrazine. However, this is an anomalous result:
γ
γ
Val-70 C 1 and C 2 possess similar shifts in both complexes
and differ by ∼3 Hz in the 2-methoxy-3-isobutylpyrazine
γ
γ
complex, under which conditions C 1 and C 2 will be strongly
coupled via the two bond homonuclear scalar coupling, which
is of similar magnitude. This results in a poor fit of R2 relaxation
γ
2
data for Val-70 C 2 (ø ) 49.7).
1
3
2
C relaxation. Thus, by analogy with the work of Ishima et
It is notable that measured S values for the two methyl
axis
al.,26 the contribution of these external protons to C R2 values
was estimated at 25% from the X-ray coordinates of MUP (R.
Bingham and S. E. V. Phillips, unpublished data). Consequently,
13
groups of certain valine residues are not identical within
experimental error. At first sight, this is inconsistent with the
requirement for their mobilities to be essentially the same, since
1
3
measured C R2 rates were multiplied by 0.75 to account for
they form part of the same isopropyl group. However, this is
2
axis
these external dipolar relaxation processes.
not necessarily reflected in equivalent S
previously,
values: as noted
2
2,39,40
order parameters of methyl groups from the
Results and Discussion
same isopropyl group may differ if the effective averaging axis
for this group makes different angles with the methyl 3-fold
13
1
Figure 1a shows a region from the C, H HSQC spectrum
of methyl 13C, 2H enriched mouse major urinary protein,
axes.
γ1
γ1
γ2
γ2
2
axis
containing valine C -H and C -H correlations. Significant
overlap is present in the spectrum, arising from the combined
effects of interference from resonances from 13CH3 isotopomers,
together with methyl resonances derived from residues other
than valine.
The constant values of S on complex formation for the
methyl groups of Val-82 is unexpected, since this residue is
(37) Timm, D. E.; Baker, L. J.; Mueller, H.; Zidek, L.; Novotny, M. V. Protein
Sci. 2001, 10, 997-1004.
(
38) Zwahlen, C.; Legault, P.; Vincent, S. J. F.; Greenblatt, J.; Konrat, R.; Kay,
L. E. J. Am. Chem. Soc. 1997, 119, 6711-6721.
(
34) Lee, A. L.; Urbauer, J. L.; Wand, A. J. J. Biomol. NMR 1997, 9, 437-
(39) Lemaster, D. M.; Kushlan, D. M. J. Am. Chem. Soc. 1996, 118, 9255-
4
40.
9264.
(
35) Lipari, G.; Szabo, A. J. Am. Chem. Soc. 1982, 104, 4546-4559.
(40) Mittermaier, A.; Kay, L. E.; FormanKay, J. D. J. Biomol. NMR 1999, 13,
181-185.
(
36) Mandel, A. M.; Akke, M.; Palmer, A. G. Mol. Biol. 1995, 246, 144-163.
J. AM. CHEM. SOC.
9
VOL. 125, NO. 51, 2003 15769