Previous Article |
Table of Contents
| Next Article
BIOLOGICAL SCIENCES / BIOPHYSICS
Intramolecular domaindomain association/dissociation and phosphoryl transfer in the mannitol transporter of Escherichia coli are not coupled
Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892-0520
Edited by Alan R. Fersht, University of Cambridge, Cambridge, United Kingdom, and approved January 11, 2007 (received for review October 17, 2006)
| Abstract |
|---|
|
|
|---|
500 s1 for A
B and B
A phosphoryl transfer, derived from qualitative line-shape analysis of 1H-15N correlation spectra taken during the course of active catalysis. Thus, on average,
80 association/dissociation events are required to effect a single phosphoryl-transfer reaction. We conclude that intramolecular phosphoryl transfer between the A and B domains of IIMtl is rate-limited by chemistry and not by the rate of formation or dissociation of a stereospecific complex in which the active sites are optimally apposed.
domain motion | NMR | protein dynamics | relaxation dispersion | phosphotransferase system
Relaxation measurements, including relaxation dispersion, have been used to study dynamics of enzyme function and protein folding on milli- to microsecond time scales (210). Recent work has suggested that the dynamics of atomic motions observed by relaxation dispersion are directly linked to catalysis (5, 7, 8). For adenylate kinase, the rate of lid opening, which involves large-amplitude correlated motions, seems to represent the rate-limiting step in catalysis (5). This observation led to the suggestion that many enzymes have evolved such that the catalytic reaction itself is so fast that catalytic power is limited not by chemistry but by the rate at which conformational rearrangements can take place to optimally align the reactive atoms (5). However, this generalization has been called into question, because the correlation between conformational changes and kinetics observed for adenylate kinase does not necessarily imply a direct link between dynamics, correlated motions, and catalysis (11).
Previously, we determined the solution structure of an analog of the product complex of the phosphoryl-transfer reaction involving the isolated IIAMtl domain and a stably phosphorylated IIBMtl domain from Escherichia coli (12). IIBMtl was stably phosphorylated by mutating the active-site cysteine (Cys-384) to serine (13), whereas the active-site histidine (His-554) of IIAMtl was mutated to glutamine to prevent transfer of the phosphoryl group back to IIAMtl in the complex. The interaction between the isolated IIAMtl and phosphoIIBMtl domains is weak, with an equilibrium dissociation constant (Kd) of
3.7 mM. In the native enzyme, however, the two domains are connected by an
21-residue flexible linker (1). The linker renders the equilibrium between associated (closed) and dissociated (open) forms of intact IIBAMtl concentration independent and is expected to stabilize the associated state, because the reduction in configurational space that accompanies tethering is equivalent to raising the effective concentration of the domains (12, 14). In this article we investigate the dynamics of intramolecular domaindomain association and dissociation between the A and B domains of single-chain IIBAMtl of E. coli by using 1H relaxation dispersion spectroscopy in the rotating frame and compare these rates with estimates of the reversible phosphoryl-transfer rates between the A and B domains derived from qualitative line-shape analysis of cross-peaks in 1H-15N correlation spectra during the course of active catalysis. We show that the intramolecular domaindomain association and dissociation rates are not rate-limiting and that phosphoryl transfer occurs, on average, once for every
80 association and dissociation events.
| Results |
|---|
|
|
|---|
The profile of 1HN/15N chemical-shift perturbations relative to the free phosphoIIBMtl and free IIAMtl is essentially the same for the mixture of isolated IIAMtl and phosphoIIBMl (comprising
30% complex) and the intact phosphoIIBAMtl, indicating that the interaction of the two domains is the same in both systems (Fig. 1). The magnitudes of the perturbations, however, are larger for phosphoIIBAMtl than for the mixture, indicating that the fraction of closed (associated) form in the single-chain phosphoIIBAMtl is >30%. The 1HN/15N chemical shifts of phosphoIIBAMtl are not affected by concentration (over a range of 50500 µM), indicating that the chemical-shift perturbation relative to the isolated IIAMtl and phosphoIIBMtl domains is solely caused by intramolecular association of the two domains.
|
|
![]() |
where IIAMtl · P-IIBMtl is the closed state in which the A and B domains are associated to form a stereospecific complex, IIAMtlP-IIBMtl is the open, dissociated, state (which, in reality, is an ensemble of states) in which there is no interaction between the A and B domains, and kOC and kCO are the unimolecular rate constants for the open-to-closed and closed-to-open transitions, respectively. For this two-state system, the overall relaxation rate in the rotating frame, R1
, is given by the sum of the intrinsic relaxation rate, R
, and the exchange of chemical-shifts contribution, Rex (9):
|
|
Because exchange is fast on the relaxation time scale, (|R
R
| << kex), R
is given by
|
|
where R
and R
are the intrinsic relaxation rates for the open and closed forms, respectively, and popen and pclosed are the populations of the open and closed states, respectively, with popen = (1 pclosed). Note that the values of R1
for the open and closed states in the absence of chemical exchange are not accessible experimentally because this is a unimolecular system in equilibrium.
In the fast-exchange limit, (
/kex)2 << 1, on the chemical-shift time scale, the Rex contribution to R1
is given by (10, 19, 20)
|
|
where 
(=
closed
open) is the frequency difference between the chemical shifts of the two states; kex is the sum of the two rate constants, kOC and kCO, for exchange between the two states; and
SL is the spin-lock frequency. In general, pclosed and (
)2 cannot be determined independently of one another. In this instance, however, the chemical shifts of the open state,
open, are known, because they are the same as those for the two isolated domains (
free). The fraction of the closed state, therefore, is given by pclosed = (
obs
free)/(
closed
free), where
obs is the observed shift measured for phosphoIIBAMtl. Thus, Rex can be expressed as
|
|
thereby enabling pclosed to be determined independently, because both
obs and
free are directly accessible experimentally.
Mapping of the residues used for analysis of the relaxation dispersion data onto the structure of the previously determined IIAMtlphosphoIIBMtl complex (12) is shown in Fig. 3. The titration spectra obtained with the isolated IIAMtl and IIBMtl domains were used to estimate the lower limit of kex and the range of 
. Because the values of 
were larger in the 1H than 15N dimension and stronger spin-lock field strengths can be applied for 1H than 15N to suppress the exchange contribution to R1
, we made use of backbone amide proton (1HN) relaxation dispersion measurements in the rotating frame (20) by using perdeuterated phosphoIIBAMtl. Among the titrating residues of IIBAMtl, 10 residues with both well resolved cross-peaks in the 2D 1H-15N correlation spectrum and large 1HN chemical-shift differences (
H/2
> 60 Hz at a 1H frequency of 750 MHz) relative to the free protein were selected for further analysis. Note that simulations of Eq. 4, varying the 
H, kex (for values larger than the estimated lower limit of kex), and
SL, indicate that a 60-Hz cutoff is required to ensure that Rex values larger than the 34% experimental errors in R1
can be detected.
|
, by using the error function
|
|
where R
(i, j) and R
(i, j) are the calculated and observed R1
values for the ith spin-lock field and jth residue, respectively, and
ij is the error obtained from Monte Carlo simulation for the single exponential fitting of R1
.
Fig. 4a shows typical R1
data measured for the backbone amide proton of Asp-454 at three different field strengths, illustrating that the data follow single exponential decays. Selected relaxation dispersion profiles for residues with (Asp-454 and Gly-394) and without (Ala-513 and Gly-592) a chemical-exchange contribution to R1
are displayed in Fig. 4b. Global fitting of the data for the 10 selected residues yielded values of kex = 3.5(±0.5) x 104 s1 and pclosed = 0.48 ± 0.06 with
2/N = 0.8 (where N is the number of degrees of freedom). The value of pclosed obtained from the relaxation dispersion data are in agreement with the approximate estimate derived from the observed chemical shifts of phospho-IIBAMtl, the chemical shifts of the isolated IIAMtl and phospho-IIBMtl domains free in solution, and the extrapolated chemical shifts of the fully saturated IIAMtlphosphoIIBMtl complex derived from titration data with the isolated domains (note >9095% saturation is difficult to achieve given the Kd of 3.7 mM). Because kex = kOC + kCO and pclosedkCO = popenkOC at equilibrium, the unimolecular association (kOC) and dissociation (kCO) rate constants in Eq. 1 are given by pclosedkex and popenkex, respectively, yielding values of kOC = 1.7(±0.3) x 104 and kOC = 1.8(±0.4) x 104 s1.
|
Kinetics of Reversible Intramolecular Phosphoryl Transfer Between the A and B Domains of IIABMtl.
Under conditions of catalytic amounts of enzyme I and HPr and excess PEP, both the A and B domains of wild-type IIBAMtl are fully phosphorylated. Preliminary experiments indicated that the isolated phosphorylated IIAMtl domain is relatively unstable at neutral pH, and irreversible dephosphorylation by hydrolysis occurs with a half-life of
20 min. In contrast, the isolated phosphorylated IIBMtl domain is rather stable, and the phosphorylated form has a half-life of
17 h. Thus, PEP is continually consumed and inorganic phosphorus generated during the time course of the reaction. The 31P-NMR spectrum of fully phosphorylated IIBAMtl displays two resonances of equal intensity: one at 12.2 ppm corresponding to phospho-Cys-384 (B domain) and the other at 6.7 ppm corresponding to phospho-His-554 (A domain) (referenced relative to inorganic phosphate at 0 ppm), in perfect agreement with previous reports on full-length IIMtl [ref. 22; see supporting information (SI) Fig. 6]. Once PEP is fully consumed, hydrolysis of phospho-His-554 occurs. The intensity of the resonances of phospho-His-554 and phospho-Cys-384 remain equal and decrease in concert to finally disappear after
4 h (SI Fig. 6). The two phosphorus resonances remain in slow exchange on the chemical-shift time scale throughout the reaction. Therefore, one can conclude that the overall phosphoryl-transfer reaction rate, k
(given by the sum of the A
B and B
A phosphoryl-transfer reaction rates), is at least 10-fold slower than 
P = 2.9 x 104 s1 (calculated from the chemical-shift difference between the two signals, 18.9 ppm, and the 31P frequency of 242.94 MHz).
To further probe the rate of reversible phosphoryl transfer between the A and B domains we recorded a series of 1H-15N heteronuclear single quantum coherence (HSQC) spectra (at a 1H frequency of 800 MHz) during active phosphoryl transfer between the two domains subsequent to the depletion of PEP. The results are shown in Fig. 5. At time 0 h (the point at which PEP is fully depleted), the 1H-15N HSQC spectrum reflects fully phosphorylated IIBAMtl, that is, the species in which both the B and A domains are phosphorylated. At 4 h subsequent to PEP depletion, the 1H-15N HSQC spectrum corresponds to the fully unphosphorylated IIBAMtl species. The intermediate time points reflect a mixture of biphosphorylated, monophosphorylated, and unphosphorylated species, with the monophosphorylated species reaching a maximum between 1 and 2 h subsequent to PEP depletion. In some instances, the chemical shifts of the cross-peaks progressively shift during the course of the reaction from the fully phosphorylated to the unphosphorylated positions without any significant line broadening. This is the fast-exchange regime on the chemical-shift time scale and is exemplified by the cross-peak for Ala-383 (Fig. 5a). 
HN for Ala-383 is 151 s1, indicating that the overall phosphoryl-transfer reaction rate k
is in excess of
750 s1 (kex > 5
HN). In the case of Asp-432 (Fig. 5b), a continual shift is also observed, but extensive line broadening is seen in the spectra taken at 1 and 2 h, characteristic of an exchange rate on the fast side of intermediate exchange: 
HN for Asp-432 is 437 s1, suggesting that k
1,000 s1 (2
HN < kex < 3
HN). This estimate is fully confirmed by the behavior of the cross-peaks of Gly-555 (
HN = 905 s1; Fig. 5c) and Thr-556 (
HN = 1,116 s1; Fig. 5d) where the cross-peaks completely disappear in the 1H-15N HSQC spectra taken at 1 and 2 h, which is indicative of the intermediate exchange regime where k

HN. Because the intensities of the 31P resonances of phospho-Cys-384 and phospho-His-554 decrease in concert and remain approximately equal to each other throughout the reaction, one can conclude that the forward and backward rate constants for phosphoryl transfer are approximately equal with a value of
500 s1.
|
| Discussion |
|---|
|
|
|---|
3.7 mM (12). Because the occupancy of the closed (associated) state of phosphoIIBAMtl determined from the relaxation dispersion measurements is
48%, one can conclude that the 21-residue linker connecting the A and B domains in IIBAMtl results in an effective local concentration of
4 mM for each domain, a value that is in good accord with that predicted from polymer chain theory (12, 23) and configurational entropy calculations (14). Moreover, if the unimolecular association rate constant, kOC (1.7 x 104 s1), is divided by the effective local concentration, the value of
5 x 106 M1·s1 for the apparent bimolecular rate constant is reasonably close to what one might expect for a diffusion-limited reaction between two proteins, particularly when factors such as the chain-diffusion coefficient of the linker are taken into account. Thus, intradomain association in the IIBAMtl system is tuned to function, which requires reasonably high occupancy coupled with rapid association and dissociation (
2 x 104 s1 in both directions) to efficiently carry out three sequential phosphoryl-transfer reactions, namely from HPr to IIAMtl, from IIAMtl to IIBMtl, and finally from IIBMtl onto the incoming sugar bound to the cytosolic side of IICMtl.
Relationship Between Phosphoryl Transfer and Domain Dynamics.
The NMR data presented in this article clearly indicate that the rate constants for the forward (A
B) and backward (B
A) phosphoryl-transfer reactions are approximately equal with values of
500 s1, which is
40-fold lower than the rate constants for intramolecular domaindomain association and dissociation. Thus, on average,
80 association/dissociation events take place for every phosphoryl-transfer reaction. One can conclude, therefore, that the rate-limiting step for intramolecular phosphoryl transfer in IIBAMtl is governed by the chemistry of the phosphoryl-transfer reaction rather than the time it takes for the two domains to form a stereospecific complex in which the active sites of the A and B domains are optimally positioned for phosphoryl transfer. In this regard, it is worth noting that in the case of all of the bimolecular complexes of the PTS solved to date, including the N-terminal domain of enzyme I (EIN)HPr (24), IIAGlcHPr (25), IIAMtlHPr (26), IIAManHPr (27), IIAGlcIIBGlc (28), and IIAMtlIIBMtl (12) complexes, either minimal or no backbone changes are required to form an optimal pentacoordinate phosphoryl-transition state.
Phosphoryl-Transfer Rates in the PTS.
There have been a number of kinetic studies on the glucose-branch PTS using rapid-quench techniques (29, 30). (Note that the kinetics of the intramolecular phosphoryl-transfer reaction studied here is not accessible to classical biochemical methodology.) The reported apparent second-order rate constants for reversible phosphoryl transfer along the enzyme I, HPr, IIAGlc, and IICBGlc pathway range from
4 x 106 to 2 x 108 M1·s1 (29). By using estimated intracellular protein concentrations of
5, 20100, 2060, and 10 µM for enzyme I (monomer), HPr, IIAGlc, and IIBCGlc, respectively (29), one can deduce that the upper limits of the pseudo-first-order rate constants for the forward and backward phosphoryl-transfer reactions in vivo are
1,000 and 160800 s1, respectively, between enzyme I and HPr; 1,2006,000 and 1,0003,000 s1, respectively, between HPr and IIAGlc, and 200600 and 40 s1 between IIAGlc and IICBGlc, respectively. Thus, the intramolecular phosphoryl-transfer rates of
500 s1 observed here for IIBAMtl are quite comparable to those in the glucose branch of the PTS.
Comparison with Adenylate Kinase.
Adenylate kinase catalyzes the reversible transfer of phosphorus from two molecules of ADP to form ATP and AMP. Its structure comprises a core and two mobile domains (the lid and the AMP-binding domains) that close over the nucleotide-binding sites (31). The two mobile domains exist in either open or closed forms, and the transition between these two states involves large-scale correlated motions within the hinge regions connecting the two domains to the core (31). The opening rates for both a thermophilic and mesophilic adenylate kinase (
40 and 290 s1, respectively) are essentially identical to the respective catalytic rates, whereas the closing rate (1,4001,600 s1) is much faster (5). These data suggest that the conformational transition from closed to open states required for product release is rate-limiting in catalysis (5). In contrast, in the IIBAMtl system, the intramolecular domaindomain opening (dissociation/product release) and closing (association/substrate binding) rates are comparable (
2 x 104 s1), and it is the phosphoryl-transfer reaction that seems to be rate-limiting (
500 s1).
The different behaviors observed for IIBAMtl and adenylate kinase may be rationalized in structural terms. In the case of IIBAMtl, the two domains are connected by a long 21-residue linker that behaves essentially as a random-coil polymer. Thus, concerted backbone motions within the linker are not required to bring the two domains together or move them apart; rather the function of the linker is to restrict diffusion of one domain relative to another to a volume limited by the average rms end-to-end distance of the linker. For adenylate kinase, however, motion of the lid and AMP-binding domains relative to the core domain is determined by concerted backbone conformational changes of a few residues at the hinge points between the domains. The concerted nature of these changes is presumably more complex and, therefore, energetically more costly and slower relative to the essentially random, largely uncorrelated nature of the motions within the IIBAMtl linker.
Biological Implications for the PTS. In the absence of external sugar, the PTS is in the resting state with all signal transducers phosphorylated (32). When the bacterium encounters external sugar in the medium, the first event involves the transfer of the phosphoryl group from IIB to the incoming sugar located on IIC, thereby switching on the PTS signaling pathway. Dephosphorylated IIB and subsequently dephosphorylated upstream mediators (e.g., IIA and enzyme I) interact with several transcription factors to turn on or off a variety of proteins involved in sugar uptake. The fast intramolecular domaindomain association between the A and B domains of IIBAMtl and subsequent phosphoryl transfer represents an efficient means for rapid initiation and amplification of sugar-uptake signaling in bacteria.
| Materials and Methods |
|---|
|
|
|---|
-D-thiogalactopyranoside at an A600 of
0.8 and harvested by centrifugation after 4 h of induction. Wild-type IIBAMtl and IIBAMtl(C384S/H554Q/C571A) were purified as described (12). Phosphorylation of the latter was carried out as described (12, 13). The isolated IIBMtl(C384S) and IIAMtl(H554QA/C571A) domains were expressed and purified as described (12).
NMR Spectroscopy.
NMR samples contained 0.5 mM 2H,15N- or 13C,15N-labeled phosphoIIBAMtl in 20 mM Tris-d11 (pH 7.4), 0.01% (wt/vol) sodium azide, and 10% 2H2O (vol/vol). 13C,15N-labeled phosphoIIBAMtl protein was used to assign the linker region and also to confirm assignments of each domain by using 3D HNCACB and CBCA(CO)NH experiments (33). The 1HN/15N assignments of the A and B domains of wild-type IIBAMtl in unphosphorylated and phosphorylated states were obtained directly by reference to those of the isolated IIAMtl (26) and IIBMtl (13, 34) domains in their respective unphosphorylated and phosphorylated states. 1DNH RDCs were obtained by taking the difference in the corresponding 1JNH couplings measured in aligned and isotropic (water) medium (15). The alignment medium comprised a neutral 5% polyacrylamide gel, prepared as described (13), and 1JNH couplings were measured from a 2D in-phase/antiphase 1H-15N HSQC spectrum (15). Calculated 1DNH values were obtained by fitting the experimental RDC data to the coordinates of each domain (12) individually by using singular value decomposition (15). Agreement between observed and calculated RDCs is expressed as an RDC R factor that scales between 0% and 100% and is defined as the ratio of the rms deviation between observed and calculated values to the expected rms deviation if all of the N-H vectors were randomly distributed; the latter is given by [2Da2(4 +
2)/5]1/2, where Da is the magnitude of the axial component of the alignment tensor, and
is the rhombicity (35).
R1
-unlike spin relaxation rates in the rotating frame for backbone amide protons were measured on 2H,15N-labeled phospho-IIBAMt1 by using the pulse sequence described in ref. 20 with minor modifications. A 3-9-19 pulse train was used for water suppression in the WATERGATE scheme, and heat compensation was used during the relaxation delay (36). Proton R1
measurements were carried out by using seven spin-lock times (4, 8, 13, 17, 25, 33, and 42 ms) on a Bruker DMX750 spectrometer at 30°C. The 1H spin-lock field strengths (
H/2
) used for the relaxation dispersion measurements were 2, 4, 6, 8, 10, and 12 kHz, at three 1H spin-lock frequencies of interest (7.03, 8.56, and 9.38 ppm) to ensure the on-resonance R1
condition. With the offsets of amide proton frequencies from the spin-lock field frequency <300 Hz, the contribution of longitudinal relaxation to R1
is <3% even for the lowest spin-lock field strength of 2 kHz, so that the on-resonance condition holds throughout the experiments (37). The relaxation dispersion curves were fit simultaneously as described in Results. Errors in optimized parameters were estimated by the jackknife method, in which one data set was removed in turn.
Time-course experiments during phosphoryl transfer between the A and B domains of wild-type IIBAMtl were carried out by using 2D 1H-15N HSQC (1H frequency = 800.13 MHz) and 1D 31P-NMR (31P frequency = 242.94 MHz) spectroscopy. 15N-labeled IIBAMtl (0.5 mM) was prepared in 20 mM Tris-d11 buffer (pH 7.4) and phosphorylated by addition of 5 µM enzyme I/5 µM HPr/5 mM MgCl2/1.5 mM PEP. The PEP stock solution was prepared in Trizma base, and the Trizma base buffer maintained the pH of the solution at 7.4 throughout the course of the reaction.
| Acknowledgements |
|---|
|
|
|---|
| Footnotes |
|---|
Abbreviations: PEP, phosphoenolpyruvate; PTS, PEP:sugar phosphotransferase system; HPr, histidine phosphocarrier protein; RDC, residual dipolar coupling; HSQC, heteronuclear single quantum coherence.
*To whom correspondence should be addressed. E-mail: mariusc{at}intra.niddk.nih.gov
Author contributions: J.-Y.S. and G.M.C. designed research; J.-Y.S. performed research; J.-Y.S., J.I., and G.M.C. analyzed data; and J.-Y.S. and G.M.C. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS direct submission.
This article contains supporting information online at www.pnas.org/cgi/content/full/0609103104/DC1.
© 2007 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||