Toxins TcdA and TcdB from Clostridioides difficile glucosylate human colon Rho GTPases. TcdA and TcdB glucosylation of RhoGTPases results in cytoskeletal changes, causing cell rounding and loss of intestinal integrity. Clostridial toxins TcdA and TcdB are proposed to catalyze glucosylation of Rho GTPases with retention of stereochemistry from UDP-glucose. We used kinetic isotope effects to analyze the mechanisms and transition-state structures of the glucohydrolase and glucosyltransferase activities of TcdB. TcdB catalyzes Rho GTPase glucosylation with retention of stereochemistry, while hydrolysis of UDP-glucose by TcdB causes inversion of stereochemistry. Kinetic analysis revealed TcdB glucosylation via the formation of a ternary complex with no intermediate, supporting an SNi mechanism with nucleophilic attack and leaving group departure occurring on the same face of the glucose ring. Kinetic isotope effects combined with quantum mechanical calculations revealed that the transition states of both glucohydrolase and glucosyltransferase activities of TcdB are highly dissociative. Specifically, the TcdB glucosyltransferase reaction proceeds via an SNi mechanism with the formation of a distinct oxocarbenium phosphate ion pair transition state where the glycosidic bond to the UDP leaving group breaks prior to attack of the threonine nucleophile from Rho GTPase.
Toxins TcdA and TcdB from Clostridioides difficile glucosylate human colon Rho GTPases. TcdA and TcdB glucosylation of RhoGTPases results in cytoskeletal changes, causing cell rounding and loss of intestinal integrity. Clostridial toxins TcdA and TcdB are proposed to catalyze glucosylation of Rho GTPases with retention of stereochemistry from UDP-glucose. We used kinetic isotope effects to analyze the mechanisms and transition-state structures of the glucohydrolase and glucosyltransferase activities of TcdB. TcdB catalyzes Rho GTPase glucosylation with retention of stereochemistry, while hydrolysis of UDP-glucose by TcdB causes inversion of stereochemistry. Kinetic analysis revealed TcdB glucosylation via the formation of a ternary complex with no intermediate, supporting an SNi mechanism with nucleophilic attack and leaving group departure occurring on the same face of the glucose ring. Kinetic isotope effects combined with quantum mechanical calculations revealed that the transition states of both glucohydrolase and glucosyltransferase activities of TcdB are highly dissociative. Specifically, the TcdB glucosyltransferase reaction proceeds via an SNi mechanism with the formation of a distinct oxocarbenium phosphate ion pair transition state where the glycosidic bond to the UDP leaving group breaks prior to attack of the threonine nucleophile from Rho GTPase.
Glycosyltransferases commonly catalyze
glycosidic bond formation
via the transfer of a sugar moiety from the activated nucleoside diphosphate
sugar donor to acceptor substrates, including monosaccharides, polysaccharides,
lipids, proteins, nucleic acids, and small organic molecules.[1] Glycosyltransferases can be classified as inverting
or retaining depending on stereochemical outcome at the anomeric carbon.[1,2] Inverting glycosyltransferases use a single-displacement mechanism
(SN2) with a single oxocarbenium ion-like transition state.
The inverting mechanism has been supported by mechanistic studies
and kinetic isotope effect measurements.[1,3,4]The reaction mechanism of retaining glycosyltransferases
can occur
by several proposed mechanisms (Figure ).[1,5] Initially, retaining glycosyltransferases
were proposed to proceed through double-displacement mechanisms, where
an active site amino acid acts as a nucleophile, attacking the anomeric
carbon to form a covalently linked glycosyl-enzyme intermediate before
the deprotonated acceptor nucleophile can attack, resulting in retained
stereochemistry of the newly formed glycosidic bond (Figure a).[1,5] However,
retaining glycosyltransferases often lack a suitably positioned amino
acid to act as the catalytic nucleophile.[1,5−11] Retention of stereochemistry without a covalently linked glycosyl-enzyme
intermediate led to the proposal for an “internal return”
SNi-like mechanism, where nucleophilic
attack occurs at the same face as leaving group departure.[7,9,12−14] The SNi-like mechanism can proceed through a short-lived
oxocarbenium ion intermediate (stepwise-mechanism) or a single oxocarbenium
ion-like transition state (concerted mechanism) (Figure b, c).[5] The SNi mechanism has been supported
by kinetic isotope effect (KIE) measurements on the solvolysis of
glycosyl fluorides in hexafluoro-2-propanol which proposed a short-lived
oxocarbenium ion pair intermediate.[15] KIE
measurements have also been reported for trehalose-6-phosphate synthetase
(OtsA), also supporting an SNi mechanism.
However, the magnitude of the anomeric carbon KIE could not distinguish
between stepwise and concerted SNi mechanisms.[16]
Figure 1
Proposed mechanisms for retaining glycosyltransferases:
(a) Double
displacement mechanism with the formation of covalently linked glycosyl
intermediate. (b) Concerted SNi mechanism.
(c) Stepwise SNi mechanism with the formation
of the oxocarbenium-phosphate ion pair.
Proposed mechanisms for retaining glycosyltransferases:
(a) Double
displacement mechanism with the formation of covalently linked glycosyl
intermediate. (b) Concerted SNi mechanism.
(c) Stepwise SNi mechanism with the formation
of the oxocarbenium-phosphate ion pair.Recent quantum mechanics/molecular mechanics (QM/MM)
studies have
analyzed the retaining glycosyltransferases OtsA,[17] α-1,2-mannosyltransferase,[18] GalNAc-transferase 2,[19] mannoglycerate
synthetase,[20] lipopolysaccharyl α-galactosyl
transferase,[21] and glucosyl-3-phosphoglycerate
synthetase.[11] With the exception of lipopolysaccharyl
α-galactosyl transferase, the enzymes were proposed to catalyze
glycosylation via stepwise SNi mechanisms
with short-lived oxocarbenium ion pair intermediates. These transition
states were highly dissociative with the bond to the UDP leaving group
being fully dissociated prior to nucleophilic attack. These studies
suggest that the stepwise SNi mechanism
is dominant for retaining glycosyltransferases.Despite the
wealth of knowledge for the glycosyltransferases, protein
O-glucosyltransferases are less common. Here, we extend studies of
glycosyltransferases to TcdB, a bacterial glucosyltransferase by the
analysis of combined experimental KIEs and QM calculations.The TcdA and TcdB toxins from Clostridioides difficile contain a glucosyltransferase (GT) domain that when processed and
released from endosomes into the cytosol of mammalian cells glucosylate
and inactivate Rho GTPases.[22−24] The TcdA and TcdB toxins are
the primary determinants of C. difficile pathogenesis and as such are attractive therapeutic targets.[25−28] Detailed knowledge of the reaction mechanism and transition states
of TcdA and TcdB defines reaction chemistry and is useful in the design
of transition state analogues. TcdA and TcdB catalyze transfer of
a glucosyl residue to a specific threonine residue in the switch I
effector regions of Rho GTPases, including Rac1.[23] The glucosyltransferase domain of TcdA and TcdB belong
to the GT44 family of glycosyltransferases, which also contains the
bacterial toxin glucosyltransferases.[29] The TcdA and TcdB glucosyltransferase domains belong to the GT-A
fold family of glycosyltransferases. The TcdA and TcdB glucosyltransferase
domains show the highest catalytic activity in the presence of a combination
of potassium, magnesium, and manganese ions.[30−32] Similar to
other retaining glycosyltransferases, in the absence of an appropriately
positioned acceptor to provide a catalytic nucleophile, the TcdA and
TcdB toxins catalyze the slow hydrolysis of UDP-glucose.[32,33]Previous KIE measurements on the glucohydrolase reactions
of TcdA
and TcdB demonstrated a highly dissociative glucocationic transition
state. This information was used to identify the iminosugars isofagomine
and noeuromycin as transition state analogue inhibitors.[33] Here, we characterize the mechanisms of both
glucohydrolase (GH) and glucosyltransferase (GT) reactions of the
TcdB GT domain. Retention of stereochemistry is confirmed for the
TcdB GT reaction. In contrast, the TcdB GH reaction is demonstrated
to proceed with the inversion of stereochemistry. Kinetic analysis
of the TcdB GT reaction demonstrated that the reaction proceeds with
the formation of a ternary complex. While positional isotope exchange
analysis demonstrated that TcdB-catalyzed glucosyltransfer occurs
without reversible rotation of the leaving phosphoryl group.KIE analysis combined with QM calculations support an SNi mechanism with an oxocarbenium ion pair transition
state for the TcdB GT reaction. Knowledge of the TcdB GH and GT transition
states is useful for the design of transition analogues with improved
potency to act as antitoxins.
Materials and Methods
Gene Expression and Protein Purification
Expression
plasmid pET28a containing C-terminal His-6 tagged TcdB-GTD (amino
acids 1–546) was kindly provided by Tam et al.[34] Conditions for the heterologous expression of TcdB-GTD
in Escherichia coli (One Shot BL21
Star (DE3) cell line were adapted from Tam et al[34] and followed the method of Paparella et al.,[33] with purification by Ni-NTA chromatography and
protein preparation in 20 mM Tris pH 7.5, 150 mM NaCl, 15% v/v glycerol,
freezing, and storage at −80 °C.Expression of Rac1
GTPase (amino acids 1–192) also followed the method of Paparella
et al.[33] Expression from an open reading
frame with an N-terminal His6 tag was cloned into pD454-SR:375415
(Atum Bio). Rac1 GTPase was expressed in E. coli (One shot BL21 Star, DE3 cell line), and recombinant protein expression
was induced with 1 mM IPTG for 4 h at 30 °C. Purification of
Rac1 was performed by Ni-NTA chromatography. The protein preparation
was equilibrated in 20 mM Tris pH 7.5, 150 mM NaCl, 15% v/v glycerol
for freezing and storage at −80 °C.
Synthesis of Isotopically Labeled Substrates
Isotopically
labeled UDP-glucose was synthesized by coupled enzymatic reactions,
as described previously.[33] In summary,
individual 1-3H, 2-3H, 6-3H, 1-14C, 1-13C, and 6-14C labeled glucoses
were used as starting materials to prepare 1″-3H,
2-3H, 6″-3H, 1″-14C,
1″-13C, and 6″-14C UDP-glucose,
respectively (see Table ). Synthesis from labeled glucose was accomplished by the sequential
actions of hexokinase, pyruvate kinase, phosphoglucomutase, UDP-glucose
pyrophosphorylase, and inorganic pyrophosphatase. Isotopically labeled
UDP-glucoses were purified by Mono Q anion exchange HPLC in ammonium
formate. Purity was established by coelution with commercially available
UDP-glucose. 1″-18O UDP-glucose and 1″-18O, 1″-13C UDP-glucose were prepared by
oxygen exchange at the anomeric carbon in H218O (Cambridge isotope labs), detailed in the study of Paparella et
al.[33,35] Isotopic enrichment of 1″-18O, 1″-13C UDP-glucose was assessed by mass spectrometry.
Isotopic enrichment of 1″-18O, 6″-14C double labeled UDP-glucose, the isotope enrichment was assessed
via mass spectrometry on a reaction without 14C that was
run in parallel. Isotope enrichment was determined to be >95%.
Mass
spectra were acquired on a Shimazu LCMS-2010EV spectrometer.
Table 1
V/K KIEs, Intrinsic KIEs, Calculated
KIEs, and Calculated EIEs for the TcdB-GTD Glucosyltransferase Reaction
Positions of isotopic labels in
UDP-glucose designed for multiple KIE measurement and the resulting
experimental and theoretical KIEs.
Positions of isotopic labels in
UDP-glucose designed for multiple KIE measurement and the resulting
experimental and theoretical KIEs.
Stereochemistry of TcdB-GTD Glucohydrolase Reaction
The stereochemistry of the TcdB-GTD glucohydrolase reaction was monitored
by nuclear magnetic resonance (NMR) spectroscopy. Reaction conditions
were as follows: 50 mM HEPES pH 7.5, 100 mM KCl, 4 mM MgCl2, 1 mM MnCl2, 112.5 μM TcdB-GTD, and 10% D2O. UDP-glucose was added to a final concentration of 20 mM, and the
mixture was placed in a 3 mm NMR tube. The first NMR spectra was acquired
15 min after the addition of UDP-glucose. All NMR data were acquired
at 25 °C using a Bruker AVIII 600 MHz spectrometer running TopSpin
3.6 and equipped with a 5 mm H/F-TCI cryogenic probe. The anomeric
(C1′) cross peak volumes of α-d-glucose and
β-d-glucose were determined in a 2D 1H-13C HSQC to monitor the rate of glucose formation coupled with
the subsequent mutarotation that occurs with glucose in solution.
Each HSQC timepoint was run for 10 min with 32 scans acquired for
each increment using a spectral width of 14 and 80 ppm for 1H and 13C, respectively, and a recycle delay of 1 s. The
data sets were collected using 2048 and 8 complex points for 1H and 13C, respectively, and the time-domain NMR
data were linear predicted and zero filled in the indirect dimension
and multiplied with a shifted sine-bell function (SSB = 2) in each
dimension prior to Fourier transformation.
Stereochemistry of TcdB-GTD Glucosyltransferase Reaction
Stereochemistry of the TcdB-GTD glucosyltransferase reaction was
assessed by NMR spectroscopy. Rac1 was purified as described above
and glucosylated by TcdB-GTD using 1-13C UDP-glucose. A
glucosylation reaction was carried out for 90 min at room temperature
and contained 50 mM HEPES pH 7.5, 100 mM KCl, 4 mM MgCl2, 1 mM MnCl2, 1 mM DTT, 200 μM 1-13C
UDP-glucose, 160 μM Rac1, and 0.5 μM TcdB-GTD. After 90
min, the solution was concentrated to 500 μL using a 10 kDa
Amicon spin column. D2O was added to a final concentration
of 10% and transferred to a 5 mm NMR tube. The final concentration
of Rac1 for NMR analysis was 280 μM. All NMR data were acquired
at 25 °C using a Bruker AVIII 600 MHz spectrometer running TopSpin
3.6 and equipped with a 5 mm H/F-TCI cryogenic probe. A 2D 1H-13C HSQC was acquired to characterize the products.
The HSQC was run for 16 h with 512 scans acquired for each increment
using a spectral width of 12 and 80 ppm for 1H and 13C, respectively, and a recycle delay of 1.5 s. The data sets
were collected using 2048 and 32 complex points for 1H
and 13C, respectively, and the time-domain NMR data were
linear-predicted and zero-filled in the indirect dimension and multiplied
with a shifted sine-bell function (SSB = 2) in each dimension prior
to Fourier transformation.
Kinetic Analysis of TcdB-GTD Glucosyltransferase Reaction
Initial reaction rate experiments for TcdB-GTD glucosytransferase
activity were measured using 6-3H UDP-glucose as a substrate
and by capturing radiolabeled glucosylated Rac1 protein by the precipitation
of the protein product using the method described in the study of
Bensadoun and Weinstein[36] and Paparella
et al.[33] Initial reaction rate studies
used fixed UDP-glucose/6-3H UDP-glucose or Rac1 at different
concentrations while varying the concentrations of the other substrate.
Terminated reaction samples were added to a precipitation mix of 120
μg/mL sodium deoxycholate, 6% trichloroacetic acid, and 10 μg/mL
BSA. Samples were incubated at room temperature for 15 min before
being centrifuged for 20 min at 13,000 rpm to separate the glucosylated
protein product. The supernatant was removed, and the protein pellet
was resuspended in 500 μL of 200 mM Tris–HCl pH 7.5,
5% SDS, and 20 mM NaOH. Samples were vortexed briefly and added to
20 mL scintillation vials to which 10 mL of Ultima Gold scintillation
fluid was added (Perkin Elmer). The amount of glucosylated Rac1 was
determined by scintillation counting of each sample for 1 min using
a Tri-carb 2910TR scintillation counter (Perkin Elmer). The initial
rates of the reaction were fit to the Michaelis Menten equation (eq ), and a double-reciprocal
Lineweaver–Burk plot was created. The data were best fit with
a convergence of lines with a shared common x intercept.
Positional Isotope Exchange Experiment
Reaction mixtures
contained 50 mM HEPES pH 7.5, 100 mM KCl, 4 mM MgCl2, 1
mM MnCl2, 1 mM DTT, 500 μM [1-13C, 1-18O] UDP-glucose, 500 μM Rac1, and 0.1 μM TcdB-GTD.
At 50, 100, 200, and 400 s, 150 μL aliquots of the reaction
were added to 20 μL of 0.4 M EDTA (final concentration = 50
mM) to quench the reaction. This corresponded to 5, 20, and 40% conversion.
Before NMR analysis, the reaction was diluted by 37% with D2O, TSP was added as an internal chemical shift reference, and samples
placed into 3 mm NMR tubes. NMR data was acquired at 25 °C using
a Bruker AVIII 600 MHz spectrometer running TopSpin 3.6 and equipped
with a 5 mm H/F-TCI cryogenic probe. The anomeric cross peak (C1′)
of UDP-glucose was determined in a high-resolution 2D 1H-13C HSQC to resolve the cross-peak belonging to the
C1′ bonded to 18O from that bonded to 16O expected to appear if positional isotope exchange (PIX) occurred
during the enzymatic reaction. Each HSQC timepoint was run for 30
min with 8 scans acquired for each increment using a spectral width
of 14 and 2 ppm for 1H and 13C, respectively,
and a recycle delay of 1.1 s. The data sets were collected using 4096
and 128 complex points for 1H and 13C, respectively,
and the time domain NMR data were linear-predicted in the indirect
dimension and multiplied with a shifted sine-bell function (SSB =
3) in each dimension prior to Fourier transformation. The extent of
the enzyme reaction was quantitated by 2D 1H-13C HSQC analysis of samples without enzymes and with enzymes to cause
40% reaction to quantitate the cross peak for glucosylated Rac1. Each
HSQC was run for 9 h with 320 scans acquired for each increment using
a spectral width of 12 and 80 ppm for 1H and 13C, respectively, and a recycle delay of 1.3 s. The data sets were
collected using 4096 and 64 complex points for 1H and 13C, respectively, and the time-domain NMR data were linear-predicted
in the indirect dimension and multiplied with a shifted sine-bell
function (SSB = 2) in each dimension prior to Fourier transformation.
Attempted Glucocation Capture in the TcdB-GT Reaction
Reaction mixtures contained 50 mM HEPES pH 7.5, 100 mM KCl, 4 mM
MgCl2, 1 mM MnCl2, 1 mM DTT, 500 μM [1-13C] UDP-glucose, 500 μM Rac1, 5 μM TcdB-GTD, and
10% methanol. After reaction completion (2 h), proteins were removed
by ultrafiltration (10 kDa Amicon spin column). The sample was made
to 500 μL with 10% D2O and transferred to a 5 mm
NMR tube. NMR data were acquired as described for analysis of the
TcdB-GTD glucosyltransferase reaction.
V/K KIEs for TcdB-GTD Glucosyltransferase Reaction
Measurements of V/K KIEs for the TcdB-GTD glucosyltransferase reaction
were performed at room temperature using the competitive radiolabeled
approach. Reaction conditions were as follows: 50 mM postassium phosphate,
pH 6.0, 100 mM KCl, 4 mM MgCl2, 1 mM MnCl2,
1 mM DTT, 30 μM Rac1, 70 μM UDP-glucose ([heavy
label] UDP-glucose + [remote label] UDP-glucose
+ unlabeled carrier), and 1 nM TcdB-GTD. For a typical KIE measurement,
a master mix containing all reaction components (except enzyme) was
prepared, and equal volume samples were aliquoted into tubes and designated
as “control reactions,” which include the no enzyme
control and 3 replicates of the starting heavy/remote label ratio.
TcdB-GTD was added to the remaining master mix, and the reaction was
allowed to proceed to 10–20% completion. Equal volume aliquots
of the TcdB-GTD reaction were added to the same volume of 100 mM EDTA
pH 8.0 to quench the reaction. Labeled glucosylated Rac1 protein was
precipitated as described above. Scintillation counting was performed
on each sample over 10 cycles at 10 min/cycle using a Tri-carb 2910TR
scintillation counter (Perkin Elmer). The dual channel instrument
registers 3H in channel A and 14C in channels
A and B. The raw data were analyzed to determine the total counts
for 3H-labeled glucosylated Rac1 and 14C-labeled
glucosylated Rac1. A control sample of 14C-labeled glucose
was used as a standard to establish the signal overlap between channel
A and channel B for 14C, as defined by eq :Spectral deconvolution
of the KIE data was achieved for 3H and 14C
using eqs and 4, respectively:V/K KIE values were
calculated from eq ,
where R0 and Rf are the ratios of [heavy label] UDP-glucose
to [remote label] UDP-glucose prior to the reaction and at partial
conversion, respectively, and f is the fraction of
substrate conversion for KIE measurements involving 1″-14C, 1″-18O, and 6″-14C
UDP-glucose:The observed KIE values
were corrected for the 6″-3H remote isotope effect
using eq :
TcdB-GTD Glucosyltransferase Forward Commitment
Forward
commitment for UDP-glucose in the TcdB-GTD glucosyltransferase reaction
was measured using the isotope trapping method.[37] A 50 μL sample of the TcdB-GTD:UDP-glucose equilibrium
complex was formed by incubating 5 μM TcdB-GTD with 70 μM
6″-14C UDP-glucose for 5 s in 50 mM phosphate pH
6.0, 100 mM KCl, 4 mM MgCl2, 1 mM MnCl2, and
1 mM DTT at room temperature. A chase solution (450 μL of 5
mM unlabeled UDP-glucose, 30 μM Rac1, 50 mM phosphate buffer
pH 6.0, 100 mM KCl, 4 mM MgCl2, 1 mM MnCl2,
and 1 mM DTT) was rapidly mixed with the ES complex solution, and
4 × 50 μL aliquots were removed from the mixture and quenched
with 50 μL of 100 mM EDTA pH 8.0 at 5, 10, 15, and 20 s after
addition of chase solution. Control reactions containing no enzyme
were processed in parallel to correct for background levels. Labeled
glucosylated Rac1 protein was precipitated as described above. Scintillation
counting was performed as described for KIE measurements. The amount
of glucosylated Rac1 produced after the addition of chase solution
was plotted as a function of time and extrapolated to time = 0. The
concentration of the enzyme substrate (ES) complex was calculated
using eq , where E represents the enzyme concentration, S represents the concentration of 6″-14C UDP-glucose,
and KM is the Michaelis constant for UDP-glucose
in the TcdB-GTD glucosyltransferase reaction.Cf was calculated using eq , where Y is the ratio of moles of glucose
product P (y-intercept at x = 0) per mole of the TcdB-GTD:UDP-glucose ES complex (eq ).
Solvent KIE
A solvent KIE for the TcdB-GTD glucosyltransferase
activity was measured using the TcdB-GTD glucosyltransferase assay
described above for V/K KIE measurements. Reactions were carried out
in 0, 16, 32, 49, 65, and 82% D2O. Glucosylated Rac1 was
isolated using the protein precipitation method described above. Initial
rates were used to calculate the solvent KIE using eq , where k is the
reaction rate, n is the atom fraction of deuterium,
and C is commitment to catalysis. The intrinsic solvent
KIE was calculated from the term 1/ϕ3.
Computational Methods
Glucohydrolase and glucosyltransferase
transition state structures for TcdB-GTD were determined from the
experimentally measured intrinsic KIE values as constraints for theoretical
transition state structures determined via density functional theory
in Gaussian 09 using the B3LYP level of theory and the 6-31 g(d) basis
set.[38] Calculated KIE values for theoretical
transition states were calculated at 25 °C (TKELV = 298.15) using
ISOEFF98[39] from scaled vibrational frequencies
(SCFACT = 0.977) of optimized structures of UDP-glucose in the ground
state and each transition state.For TcdB-GTD glucosyltransferase
transition state calculations, a modified threonine acceptor nucleophile
was included in optimized structures of the transition state. The
input coordinates for the optimization of UDP-glucose in the ground
state were taken from the x-ray crystal structure of the TcdA-GTD
in complex with UDP-glucose (PDB ID: 3SRZ). The fully elaborated UDP-glucose
structure was truncated to glucose-1-PO4. Truncation of
reactants for Gaussian calculations (glucose-1-PO4 and
a cut-off threonine nucleophile) differs for the reality of the UDP-Mg2+ and UDP-Mn2+ leaving groups, resulting from the
combined use of Mg2+ and Mn2+ in KIE experiments.
Thus, the leaving group potential of the calculated and actual leaving
groups should be noted. All structures were assigned a formal charge
of 0 and a multiplicity of 1. The optimized ground-state structure
was located as the global minimum and displayed no imaginary frequencies.For the TcdB-GTD glucohydrolase reaction, initial transition state
calculations were performed by fixing the 1–C–1–O(Pi) and 1–C–O(H2O) bonds and varying
their lengths from 1.8 to 3.0 Å in 0.2 Å increments. Transition
state structures were refined from the results of the first calculation
by varying the 1–C–1–O(Pi) bond length
from 2.2 to 2.4 Å in 0.01 Å increments while keeping the
1–C–O(H2O) bond length fixed at 2.2 Å.
The final transition state structure had a single imaginary frequency.
For the TcdB-GTD glucohydrolase reaction, a stepwise SNi mechanism was investigated, where the 1–C–1–O(Pi) bond was fixed at 3.0 Å and the 1–C–O(Thr)
and 1–O(Pi)–H(Thr) bonds were fixed and varied
from 1.8 to 3.0 Å in 0.2 Å increments. The final transition
state structure was consistent with a distinct oxocarbenium ion intermediate.
Equilibrium isotope effects (EIEs) for the intermediate were calculated
with ISOEFF98 using the optimized structure of the UDP-glucose ground
state and the oxocarbenium ion. Electrostatic potential maps for the
TcdB-GTD glucohydrolase and glucosyltransferase ground state and transition
state structures were visualized in GaussView 6.0 (isovalue = 0.04)
from the electron density and potential cubes acquired from the checkpoint
files of the geometry optimization and frequency calculations.
Results and Discussion
Stereochemistry of TcdB Glucohydrolase and Glucosyltransferase
Reactions
Transition state analysis of the TcdB glucosyltransferase
domain (TcdB-GTD) was aided by the stereochemical outcome of the TcdB-GTD
GH and GT reactions. We used 2D 1H-13C HSQC
NMR spectroscopy to monitor the formation of α-d-glucose
or β-d-glucose from TcdB catalyzed hydrolysis of UDP-glucose.
However, d-glucose mutarotates between α- and β-anomeric
forms with an equilibrium mixture of 36% α-d-glucose
and 64% of β-d-glucose at physiological pH values.[40] The mutarotation of glucose in aqueous solution
reaches equilibrium in approximately 2 hours.[41−43] For α-d-glucose and β-d-glucose, the 13C1-1H anomeric peaks are found at approximately 5.3 and 4.7 ppm,
respectively. TcdB-GTD (112.5 μM) and UDP-glucose (20 mM) were
incubated at 25 °C, and NMR spectra were acquired 15 min after
enzyme addition (approximately 90 turnovers). At this early time point,
the anomeric carbon peak ratio corresponding to β-d-glucose/α-d-glucose = 3.9. Following equilibration
of the glucose product at 180 min, the ratio stabilized at β-d-glucose/α-d-glucose = 1.4 (Figures a and S1). The TcdB-GTD GH reaction results in the formation of
β-d-glucose, establishing inversion of stereochemistry
for TcdB-GTD GH activity.
Figure 2
NMR analysis of stereochemistry of TcdB-GTD
glucohydrolase and
glucosyltransferase reactions: (a) 2D 1H-13C
HSQC spectrum was enriched in β-d-glucose early in
the TcdB-GTD glucohydrolase reaction (blue) relative to mutarotation
equilibrium (red). (b) 2D 1H-13C HSQC spectrum
for the C1’ anomeric carbon of glucosylated Rac1 formed in
the TcdB-GTD glucosyltransferase reaction. The C1 and C1’ anomeric
carbons of α-glucose, β-glucose, α-methyl-glucoside,
β-methyl-glucoside, and α-glucosylated threonine are indicated.
NMR analysis of stereochemistry of TcdB-GTD
glucohydrolase and
glucosyltransferase reactions: (a) 2D 1H-13C
HSQC spectrum was enriched in β-d-glucose early in
the TcdB-GTD glucohydrolase reaction (blue) relative to mutarotation
equilibrium (red). (b) 2D 1H-13C HSQC spectrum
for the C1’ anomeric carbon of glucosylated Rac1 formed in
the TcdB-GTD glucosyltransferase reaction. The C1 and C1’ anomeric
carbons of α-glucose, β-glucose, α-methyl-glucoside,
β-methyl-glucoside, and α-glucosylated threonine are indicated.Stereochemical outcome of the TcdB-GTD GT reaction
was monitored
by 2D 1H-13C HSQC NMR for formation of the C1-H
anomeric peak of the glucosylated Rho GTPase acceptor protein (Rac1).
Based on earlier reports of the related C. sordellii lethal toxin, also a member of the large clostridial toxins, we
anticipated TcdB-GTD to catalyze glucosylation of Rho GTPases with
retention of stereochemistry.[43,45] TcdB-GTD was incubated
with Rac1 and [1-13C glucose]UDP-glucose, and the C1-H
anomeric peak corresponding to glucosylated threonine of Rac1 was
monitored (Figure b). The peak corresponding to the C1-H anomeric peak of glucosylated
Rac1 appeared at approximately 5.0 ppm, consistent with previous studies.[44] As controls, we acquired spectra of α-methyl-glucoside,
β-methyl-glucoside, and α-glucosylated threonine. The
overlay of the NMR spectra indicates α-glucosylated threonine,
thus confirming retention of stereochemistry for the TcdB-GTD GT reaction
(Figure b).
Kinetics of the TcdB Glucosyltransferase Reaction
Retention
of stereochemistry for the TcdB-GTD GT reaction implies a double-displacement
mechanism or an SNi mechanism. Steady
state kinetic mechanisms can indicate whether TcdB GT activity proceeds
via a double displacement mechanism with a covalent intermediate or
through a sequential mechanism with ternary complex formation.[46] Steady state kinetic measurements on the TcdB-GTD
GT reaction were used to distinguish between these mechanisms. Double
displacement mechanisms predict double reciprocal plots with varied
substrates to be parallel lines while for a sequential mechanism (SNi) with ternary complex formation, the lines
on the double reciprocal plot intersect on the left of the x-axis.[47] A sequential mechanism
would support TcdB glucosyltransfer via an SNi mechanism. Double reciprocal plots indicated a sequential mechanism,
indicating that both UDP-glucose and Rac1 need to bind to TcdB-GTD,
forming a ternary complex before reacting through an SNi mechanism (Figure ). It has been proposed that TcdB-GTD binds to UDP-glucose
prior to Rac1 binding, as indicated by X-ray crystal structures of
TcdA-GTD in complex with UDP and a recently published crystal structure
of the TcdB-GTD:Cdc42:UDP complex.[47,48]
Figure 3
Steady state
kinetic measurements on the TcdB-GTD glucosyltransferase
reaction. Initial rates were fit to the Michaelis–Menten equation,
and double reciprocal plots show lines intersecting on the x-axis. (a) Varying UDP-glucose concentrations at fixed
Rac1 concentrations. (b) Varying Rac1 concentrations at fixed UDP-glucose
concentrations.
Steady state
kinetic measurements on the TcdB-GTD glucosyltransferase
reaction. Initial rates were fit to the Michaelis–Menten equation,
and double reciprocal plots show lines intersecting on the x-axis. (a) Varying UDP-glucose concentrations at fixed
Rac1 concentrations. (b) Varying Rac1 concentrations at fixed UDP-glucose
concentrations.
Positional Isotope Exchange
The chemical steps of the
TcdB-GTD GT reaction were analyzed by 18O-UDP-[13C]glucose PIX experiments. PIX experiments detect reversible enzyme
reaction steps and can indicate the lifetimes of potential intermediates.[3] The SNi mechanism
suggests loss of the bond to UDP prior to attack by the incoming threonine
nucleophile. If there is a stabilized glucocation intermediate, the
β-phosphoryl group of UDP may be free to rotate, recombine with
the glucocation, and equilibrate with the substrate UDP-glucose, to
exchange the bridging phosphoryl 18O for 16O
(Figure a). Chemical
reversibility for an oxocarbenium ion pair intermediate in the GT
reaction of TcdB-GTD is tested by this PIX experiment.
Figure 4
Positional isotope exchange
analysis of TcdB glucosyltransferase:
(a) Positional isotope exchange reaction scheme. A glucocation with
sufficient lifetime permits the weakly-bonded β-phosphoryl group
of UDP-glucose to rotate and return to the reactant state with reposition
of the bridging 18O to nonbridging 18O with
a consequent NMR shift of the 13C-glucose anomeric carbon.
(b) 2D 1H-13C HSQC spectrum demonstrating the
transfer of glucose from UDP-glucose to form glucosyltated Rac1. Red
is TcdB-GTD glucosyltransferase products at 40% completion and blue
is the no enzyme control. c) 2D 1H-13C HSQC
spectrum showing the anomeric 13C1’ carbon of [13C1’-bridge-18O]UDP-glucose from the TcdB-GTD
glucosyltransferase reaction quenched at 5% (blue), 20% (red), and
40% (green) completion. A TcdB-GTD glucosyltransferase reaction quenched
at 40% completion (observed). As a PIX control for observation to
nonbridge-18O, the 2D 1H-13C HSQC
spectrum of unlabeled UDP-glucose is compared (purple).
Positional isotope exchange
analysis of TcdB glucosyltransferase:
(a) Positional isotope exchange reaction scheme. A glucocation with
sufficient lifetime permits the weakly-bonded β-phosphoryl group
of UDP-glucose to rotate and return to the reactant state with reposition
of the bridging 18O to nonbridging 18O with
a consequent NMR shift of the 13C-glucose anomeric carbon.
(b) 2D 1H-13C HSQC spectrum demonstrating the
transfer of glucose from UDP-glucose to form glucosyltated Rac1. Red
is TcdB-GTD glucosyltransferase products at 40% completion and blue
is the no enzyme control. c) 2D 1H-13C HSQC
spectrum showing the anomeric 13C1’ carbon of [13C1’-bridge-18O]UDP-glucose from the TcdB-GTD
glucosyltransferase reaction quenched at 5% (blue), 20% (red), and
40% (green) completion. A TcdB-GTD glucosyltransferase reaction quenched
at 40% completion (observed). As a PIX control for observation to
nonbridge-18O, the 2D 1H-13C HSQC
spectrum of unlabeled UDP-glucose is compared (purple).Double-labeled [1″-13C, bridge-β-18O]UDP-glucose was synthesized to test for PIX in the TcdB
reaction.
The presence of the 18O label at the bridging position
causes an upfield 13C NMR shift in the doublet from 98.33/98.38
to 98.30/98.35 ppm relative to 16O label at the bridging
position (Figure c).
The anomeric carbon is split into a doublet because it is coupled
to the 31P via 2J (13C, 31P) scalar coupling – 31P is spin 1/2, causing the 13C signal to split.TcdB-GTD was incubated with excess
[1″-13C, 1″-18O]UDP-glucose and
Rac1, and samples were quenched at 5, 20,
and 40% conversion to UDP and glucosylated Rac1. For each sample,
no shift in the UDP-glucose anomeric carbon peaks were detected (Figure ), indicating that
no PIX occurs for the TcdB-GTD GT reaction. Thus, the 18O label does not scramble from the bridging to the nonbridging position
during the glucosyltransferase reaction.A lack of PIX exchange
establishes that the glucocation formed
in the reaction chemistry is not free to permit β-phosphoryl
rotation and reversible reaction to reform UDP-glucose. Chemical possibilities
are that glucosyl loss is not complete prior to attack of the nucleophilic
threonine of Rac1 and/or the β-phosphate of the UDP leaving
group is restricted in the catalytic site and not free to rotate and/or
that the glucocation transition state has low bond orders to both
the leaving group and attacking nucleophile, but its lifetime is too
short to permit phosphoryl rotational equilibrium. This result is
consistent with other retaining glycosyltransferases including galactosyltransferase,[3] glycogen synthase,[49] and sucrose synthetase,[50] where no PIX
effect has been reported.
Attempted Glucocation Capture
PIX analysis suggested
a glucocationic transition state too sequestered to too short-lived
to permit rotational equilibration of the β-phosphoryl group
of UDP-glucose. An independent analysis of the presence of a stabilized
glucocation is its reactivity (capture) by small nucleophiles. Capture
of an oxocarbenium ion intermediate was attempted in the TcdB-GTD
GT reaction in the presence of 10% methanol. An oxocarbenium ion intermediate
with lifetime and geometry to permit diffusional approach will react
with methanol as a nucleophilic acceptor to form α-methyl or
β-methyl glucose. No methyl-glucose was detected (Figures and S2). Consistent with PIX, the TcdB-GTD GT reaction does not contain
an oxocarbenium ion pair with a chemically reactive environment or
reversible exchange.
KIE on the TcdB Glucosyltransferase Reaction
A glucosyltransfer
SNi mechanism can proceed through an oxocarbenium
ion intermediate or a concerted mechanism with bond making to the
Rac1 threonine simultaneous with loss of the UDP leaving group. We
measured KIEs on the TcdB-GTD GT reaction to resolve these differences.
KIEs were measured using the competitive radiolabeled approach with
isotopically labeled UDP-glucose substrates. This method yields (V/K)
KIEs, which report on all steps from substrate binding up to and including
the first irreversible chemical step.[51,52] For the TcdB-GTD
GT reaction, departure of the UDP leaving group is considered to be
irreversible, as reversal of enzyme-bound UDP would result in PIX
exchange. Experimental KIEs can be altered or obscured by forward
commitment factors (Cf). However, intrinsic
KIEs that reflect only the chemical step are required to provide transition
state information. Correction of experimental KIEs to intrinsic KIEs
involves the measurement of commitment factors, the probability that
molecules in the Michaelis complex will partition to product or be
released to the unreacted substrate pool. Experimental KIEs on the
TcdB-GTD GT reaction with excess Rac1 and UDP-glucose concentrations
at physiological pH resulted in KIEs close to unity, indicating significant
forward commitment (data not shown). Therefore, we measured the forward
commitment factor associated with the TcdB-GTD GT reaction using the
substrate trapping method developed by Rose.[37] The forward commitment factor was measured to be 0.285 ± 0.012
and was used to correct the measured KIEs for the TcdB-GTD GT reaction
(Figure S3).UDP-glucose was synthesized
enzymatically or chemoenzymatically from appropriately labeled glucose
as the starting material to yield [1″-3H], [2″-3H], [1″-14C], [1″-18O
6″-14C], [6″-3H], and [6″-14C]UDP-glucose substrates (Figure S4). [6″-14C]Glucose locates the heavy atom 3 bonds
away from the anomeric carbon, a position where 14C KIE
is known to be unity and therefore served as the reference. Therefore,
[1″-3H], [2″-3H], and [6″-3H] KIEs were measured in pairs with [6″-14C]UDP-glucose serving as the control (Table ). KIEs for [1″-14C], [1″-18O 6″-14C], and [6″-14C] were measured in pairs with [6″-3H]UDP-glucose
serving as the remote label and were corrected for [6″-3H] remote effects by the latter (Table and Figure S4). Previously reported intrinsic KIEs for the TcdB-GTD GH reaction
are also shown for comparison (Table ).[33]
Table 2
Summary of V/K KIEs, Intrinsic KIEs
and Calculated KIEs from QM Calculations for the TcdB-GTD Glucohydrolase
Reaction
Positions of isotopic labels in
the substrate UDP-glucose designed for multiple KIE measurement and
the resulting experimental and theoretical KIEs. Colored labels correspond
to those of Table .
Experimental and intrinsic
KIEs
were previously measured and taken from ref (28). Six independent measurements
were performed for each KIE, and the data represent the mean values
± SEM.
Experimental
KIEs were corrected
for forward commitment (eqs –9) to give intrinsic KIE values.
The remote 6″-14C value of unity provides the standard
unity reference for all other isotope effects.
Positions of isotopic labels in
the substrate UDP-glucose designed for multiple KIE measurement and
the resulting experimental and theoretical KIEs. Colored labels correspond
to those of Table .Experimental and intrinsic
KIEs
were previously measured and taken from ref (28). Six independent measurements
were performed for each KIE, and the data represent the mean values
± SEM.Experimental
KIEs were corrected
for forward commitment (eqs –9) to give intrinsic KIE values.
The remote 6″-14C value of unity provides the standard
unity reference for all other isotope effects.The [1″-14C] primary KIE reports
on bonding and
reaction coordinate motion of the anomeric carbon at the transition
state. This isotope effect is typically 1.00–1.02 (klight/kheavy) for
an SN1 mechanism with a discrete oxocarbenium ion transition
state, 1.025–1.06 for an SN1-like concerted mechanism
with partial bonds remaining to the anomeric carbon at the transition
state, and 1.06–1.16 for an SN2 associative concerted
mechanism.[53] The [1″-14C] KIE was 1.010 ± 0.001, indicating the presence of an oxocarbenium
ion pair between the anionic UDP leaving group oxygen and the cationic
anomeric carbon, thereby supporting an SNi mechanism with an oxocarbenium ion transition state.The [1″-18O] primary KIE reports on the extent
of cleavage of the UDP-glucose glycosidic bond at the transition state.
The [1″-18O] KIE for the TcdB-GTD GT reaction was
measured to be 1.037 ± 0.002. A fully cleaved glycosidic bond
to generate an anionic leaving group is expected to generate a maximum
KIE of 1.047.[54] For TcdB-GTD, the KIE indicates
substantial cleavage of the UDP-glucose glycosidic bond at the transition
state. The magnitude of the TcdB-GTD [1″-18O] KIE
could also reflect full departure of the UDP leaving group with partial
protonation of the departing oxygen, a transition state structure
precedented by trehalose-6-phosphate synthetase (OtsA).[15] Therefore, the [1″-18O] KIE
for the TcdB-GTD GT reaction is interpreted to reflect substantial
leaving group bond cleavage with partial protonation of the leaving
group oxygen.The α-secondary [1″-3H]
KIE reports on
the changes in bond hybridization at the anomeric carbon (sp3 vs sp2) reflecting the degree of oxocarbenium ion character
at the transition state.[53,55] An sp2 geometry
at the transition state creates out-of-plane bending freedom of the
hydrogen atom, generating a normal isotope effect. The [1″-3H] KIE was measured to be 1.064 ± 0.003. Oxocarbenium
ionlike transition states typically result in α-secondary KIEs
of magnitude >1.15.[53,55] In comparison, the [1″-3H] KIE for the TcdB-GTD GH reaction was measured to be 1.216,
and for the retaining glycosyltransferase OtsA, the α-secondary
KIE is reported to be 1.284.[16] The relatively
small α-secondary KIE for the TcdB-GTD GT reaction could reflect
a crowded environment at the transition state, where the out-of-plane
bending mode that would normally occur for an sp2 hydrogen
atom is restricted by the encroaching Rac1 threonine acceptor nucleophile.The β-secondary [2″-3H] KIE reports on
the degree of hyperconjugation between the σ(C–H) orbital
at C2 to the σ*(C–O) orbital from the anomeric carbon
to the UDP leaving group.[53,56] The [2″-3H] KIE for the TcdB-GTD GT reaction of 1.032 ± 0.002
suggests that the C2–H2 bond is near-perpendicular to the C1-UDP
bond at the transition state.[56] β-secondary
KIEs of magnitude >1.07 have been reported for other N-glycohydrolases
and glycosyltransferases with oxocarbenium ion-like transition states.[16,53] As the primary [1″-14C] and [1″-18O] KIEs support an oxocarbenium ionlike transition state, the [2″-3H] KIE reflects an unusual transition state geometry with
little hyperconjugation from the σ(C–H) orbital at C2″.
Oxocarbenium ion transition states are expected to generate C5-O5-C1-C2
atoms of the glucose ring in a near coplanar geometry. Four possible
conformations of the hexopyranose ring can accommodate the co-planar
arrangement, 2,5B, B2,5, 3H4, and 4H3 (Figure S5). Of the four possible structures, only the B2,5 hexopyranose
structure is consistent with the small β-secondary KIEs.[12,57] The B2,5 hexopyranose structure is likely to be a feature
of the transition state geometry based on the small β-secondary
[2″-3H] KIE of (1.014) for the TcdB-GTD GH reaction.A remote [6″-3H] KIE was measured to be 1.071
± 0.002 for the TcdB-GTD GT reaction. Remote tritium binding
isotope effects contribute to V/K isotope effects and typically arise
from enzyme-induced distortion of the sp3 geometry at C6″.
This is a common feature among GH and GT enzymes. A similar remote
tritium KIE was also observed for the TcdB-GTD GH reaction (Table ).[52,53]The solvent deuterium KIE for the TcdB-GTD GT reaction was
measured
in the presence of varying amounts of D2O. The solvent
KIE was calculated to be 2.24 ± 0.15 and a linear function of
deuterium concentration. A single proton transfer occurring at the
TcdB-GTD GT transition state is indicated (Figure S6).Kinetic analysis and the KIEs support an SNi mechanism that occurs via an oxocarbenium ion
pair transition state.
Unusually small α-secondary and β-secondary KIEs are attributed
to a B2,5 hexopyranose geometry and restricted environment
for the ternary complex at the transition state.
Computational Analysis of the TcdB Glucohydrolase and Glucosyltransferase
Transition States
Transition states for the TcdB-GTD GH and
GT reactions were modeled using the experimental intrinsic KIEs as
constraints for computational chemistry. Geometry for the transition
state calculations used the initial coordinates for UDP-glucose (ground
state) taken from a crystal structure of the TcdA GT domain in complex
with UDP-glucose (PDB ID: 3SRZ). A simplified atomic cutoff model
containing 29 atoms was used in the calculations. For optimization
of the GH transition state structure, a water molecule was included
on the β-face of glucose (32 atoms) (Figure b). Currently, there is no available structure
of the ternary complex of TcdB-GTD in complex with UDP-glucose and
an acceptor Rho GTPase protein. Although a structure of a TcdB-GTD:Cdc42:UDP
complex has been recently published, glucose is absent as UDP-glucose
was hydrolyzed during crystallization.[48] Analysis of the GT transition state included a truncated nucleophilic
threonine side chain placed on the α-face of the glucose ring
(38 atoms) (Figure d). Calculations were performed in Gaussian 09 using the B3LYP level
of theory and the 6-31g(d) basis set.[38]
Figure 5
Optimized
ground state and transition structures for TcdB glucohydrolase
and glucosyltransferase reactions. (a) UDP-glucose structure is taken
from PDB 3SRZ. The glycosidic bond is indicated by an asterisk. (b)
Cut-off model for the ground state used to model the TcdB-GTD glucohydrolase
transition state structure. Calculations were made by varying the
lengths of the C–O bonds indicated by d1 and d2. (c) Calculated TcdB-GTD glucohydrolase transition state
structure where calculated KIEs provided the best fit to experimental
KIEs. (d) Cut-off model for the optimized ground state structure used
to model the TcdB-GTD glucosyltransferase transition state. Calculations
fixed the C–O bond (d1) at 3.0 Å and
varied the lengths of the C–O bond and O–H bonds (d2 and d3, respectively). (e) Calculated
TcdB-GTD glucosyltransferase transition state where the theoretical
and intrinsic KIEs were in the best agreement.
Optimized
ground state and transition structures for TcdB glucohydrolase
and glucosyltransferase reactions. (a) UDP-glucose structure is taken
from PDB 3SRZ. The glycosidic bond is indicated by an asterisk. (b)
Cut-off model for the ground state used to model the TcdB-GTD glucohydrolase
transition state structure. Calculations were made by varying the
lengths of the C–O bonds indicated by d1 and d2. (c) Calculated TcdB-GTD glucohydrolase transition state
structure where calculated KIEs provided the best fit to experimental
KIEs. (d) Cut-off model for the optimized ground state structure used
to model the TcdB-GTD glucosyltransferase transition state. Calculations
fixed the C–O bond (d1) at 3.0 Å and
varied the lengths of the C–O bond and O–H bonds (d2 and d3, respectively). (e) Calculated
TcdB-GTD glucosyltransferase transition state where the theoretical
and intrinsic KIEs were in the best agreement.Theoretical transition states for the TcdB-GTD
GH reaction were
generated by varying the lengths of the C–O bond between the
anomeric carbon and phosphate leaving group and the C–O bond
between the anomeric carbon and incoming water nucleophile, marked
as d1 and d2, respectively (Figure b). Theoretical KIEs
for each potential transition state were calculated from scaled vibrational
frequencies in ISOEFF98.[39] Bond lengths
from the anomeric carbon to the phosphate leaving group or to the
incoming nucleophilic water were varied from 1.8–3.0 Å
to find a theoretical transition state that would generate calculated
KIEs to best match the experiment determinations (Table ).For the GH reaction,
a theoretical transition state structure with
C–O bond lengths of 2.33 and 2.20 Å to the phosphate leaving
group and incoming water nucleophile, respectively, gave the best
match between experimental and measured intrinsic KIEs (Figure c and Table ). The [1″-14C] calculated
KIE was outside the experimental limits of the intrinsic KIE. However,
experimental and calculated [1″-14C] KIEs remains
in agreement with a dissociative transition state with oxocarbenium
ion character. The C5–O5–C1–C2 glucosyl atoms
of the TcdB-GTD GH transition state are in a coplanar arrangement
(Figure S7). The best-fit GH transition
state from the family of KIEs predict the hexopyranose ring to be
in the 4H3 half-chair conformation (Figures S5 and S7).[58] The C4 atom is 0.1 Å above the plane on the β-face while
the C3 atom is 0.6 Å below the plane on the α-face. Only
the small 2-3H β-secondary KIE is at variance with
this transition state, which predicts a B2,5 transition
state. It is possible that the geometry of the β-secondary C2-3H atom could be influenced by a strong transition state interaction
with TcdB-GTD catalytic site amino acids.Calculation of the
TcdB-GTD GT transition state structure fixed
the C–O bond between the anomeric carbon and phosphate leaving
group at 3.0 Å for the test transition states (d1 in Figure d). This
reflects full loss of the C1-UDP bond to permit initiation of nucleophilic
attack from the same face by the Rac1 threonine. The C–O bond
distances to the nucleophilic oxygen atom and for proton addition
to the leaving group O atom were varied from 1.8–3.0 Å
and 2.4–1.0 Å, respectively (d2 and d3 in Figure d). The computed GT transition state that best matched the [1″-14C] and [1″-18O] KIEs had C–O bond
lengths of 3.0 Å to the anomeric carbon and an O–H bond
length of 2.4 Å for protonation of the phosphate leaving group.
Transition state geometry for the TcdB-GTD GT reaction supports an
oxocarbenium ion pair at the transition state, causing the glucosyl
C5–O5–C1–C2 atoms to be coplanar (Figure S7). The best-fit transition state from
the family of KIEs places the hexopyranose ring in an envelope conformation
(4E), where the C4 atom is 0.7 Å above the plane on
the β-face and the C3 atom is within the plane.[58] Similar to the TcdB-GTD GH transition state, the small
experimental α and β secondary KIEs may result from a
restricted active site when the UDP leaving group and Rac1 both crowd
the same face of the anomeric carbon. These transition state structures
were calculated without contributions from the active site amino acid
interactions that may influence bond polarization.Forward motion
through the reaction coordinate for the GT reaction
is supported by the PIX experiment, revealing an irreversible and/or
restricted catalytic environment for the GT reaction, eliminating
chemical reversal with torsional β-phosphoryl motion. The theoretical
structure most closely related to the KIE constraints (Figure d) had imaginary frequencies
along the O–H bond corresponding to proton transfer (d3) and the C–O bond between the anomeric carbon
and the nucleophilic oxygen (d2) and is thus not
a true transition state structure by the definition of a single imaginary
frequency. Accordingly, equilibrium isotope effect (EIE) values were
calculated. The intrinsic [1″-14C] KIE matched a
calculated EIE of 1.011 between UDP-glucose and the UDP + oxocarbenium
ion, where the oxocarbenium ion has a lifetime that is adequate to
permit bond vibrational modes to equilibrate. The intrinsic [1″-3H] and [2″-3H] KIEs also matched the calculated
EIEs, consistent with a short-lived oxocarbenium ion intermediate.
The [1″-18O] calculated EIE was smaller than the
experimental KIE, reflecting uncertainty in the degree of protonation
for the leaving group oxygen atom.Electrostatic potential surface
(EPS) maps for the optimized TcdB-GTD
ground state and GH and GT transition state structures indicate the
distribution of partial positive (blue) and partial negative (red)
charge (Figure ).
Key features of the TcdB-GTD GH and GT transition states center on
the charge distribution around the cationic anomeric carbon of the
glucose ring (Figure e–i). The positive charge on the anomeric carbon is most pronounced
for the TcdB-GTD GT transition state (Figure i). Isofagomine, a 2,5 dideoxy-1-imino mimic
of glucose, is a nanomolar inhibitor of the TcdA and TcdB toxins and
binds only in the presence of UDP33. The cationic isofagomine
imino group forms a 2.6 Å ion pair with the β-phosphoryl
group of UDP to mimic the transition state geometry. An EPS map from
the crystal structure of isofagomine-TcdB was also generated (Figure j). The positive
charge observed at the protonated endocyclic nitrogen atom of isofagomine
is isosteric to the positively charged anomeric carbon in the TcdB-GTD
GT transition state (Figure j).
Figure 6
Structures and electrostatic potential surfaces for the TcdB-GTD
ground state, GH and GT transition states and isofagomine, a TcdB-GTD
transition state analogue. Chemical structures are shown on top and
are the (a) TcdB-GTD ground state. (b) TcdB-GTD GH transition state.
(c) TcdB-GTD GT transition state. (d) Protonated isofagomine. Geometric
and electrostatic potential surfaces (EPS) maps for the TcdB-GTD ground
state and transition state structures and transition state analogue,
isofagomine are shown in panels e to j. Electron-rich charge is red,
and electron-deficient charge is blue. Transparent EPS maps with imbedded
structures are as follows: (e) TcdB-GTD ground state, (f) view of
the anomeric carbon from the β-face of glucose.
(g) TcdB-GTD glucohydrolase transition state structure, (h) TcdB-GTD
glucosyltransferase transition state structure, (i) TcdB-GTD glucosyltransferase
transition state structure viewed from the anomeric carbon on the
β-face of glucose, and (j) structure of the
transition state analogue isofagomine with view of the protonated
nitrogen atom.
Structures and electrostatic potential surfaces for the TcdB-GTD
ground state, GH and GT transition states and isofagomine, a TcdB-GTD
transition state analogue. Chemical structures are shown on top and
are the (a) TcdB-GTD ground state. (b) TcdB-GTD GH transition state.
(c) TcdB-GTD GT transition state. (d) Protonated isofagomine. Geometric
and electrostatic potential surfaces (EPS) maps for the TcdB-GTD ground
state and transition state structures and transition state analogue,
isofagomine are shown in panels e to j. Electron-rich charge is red,
and electron-deficient charge is blue. Transparent EPS maps with imbedded
structures are as follows: (e) TcdB-GTD ground state, (f) view of
the anomeric carbon from the β-face of glucose.
(g) TcdB-GTD glucohydrolase transition state structure, (h) TcdB-GTD
glucosyltransferase transition state structure, (i) TcdB-GTD glucosyltransferase
transition state structure viewed from the anomeric carbon on the
β-face of glucose, and (j) structure of the
transition state analogue isofagomine with view of the protonated
nitrogen atom.
Conclusions
Organisms express several hundred glycosyltransferases,
a handful
of which have been characterized for transition state analysis.[1,2] In contrast, humans express only three UDP-glucose based protein
O-glucosyltransferases, and they are also rare in bacterial functions.[59] Mechanistic and KIE analysis of the TcdB-GTD
GT transition state expands mechanistic knowledge of the glycosyltransferases
to the clinically important bacterial toxins. Knowledge of the transition
states for TcdB-GTD glucosyltransferases has also provided the first
generation of small molecule inhibitors.Glucosylation of Rho
GTPases by the TcdB toxin from C. difficile occurs with overall retention of stereochemistry
and KIE measurements support an SNi mechanism
with a distinct oxocarbenium phosphate ion pair transition state.
The glucocation complex has a lifetime long enough to conform to EIE
and thus shares characteristics of an intermediate. However, the complex
is too transient or too constrained to permit torsional equilibration
of the UDP-phosphoryl group. Combined KIE and QM calculations demonstrate
that the glycosidic bond to the leaving group is cleaved before the
approach of the nucleophilic threonine and that the distances from
the anomeric carbon to the leaving group and the nucleophile are large
(3.0 Å), indicating the presence of an oxocarbenium ion pair.
The TcdB glucosyltransferase is therefore similar to some other retaining
glycosyltransferases reported to react by an SNi mechanism and ion pair transition states. Enzymes forming
glucocation-like transition states are susceptible to inhibition by
iminosugar analogues that mimic the cationic nature of the transition
state. Thus, isofagomine and noeuromycin mimic the distinct oxocarbenium
ion present in the TcdB-GTD GT transition state and are nanomolar
inhibitors of the enzyme. The transition state structure of the TcdB
glycosyltransferase domain will be useful in designing second-generation
transition state analogues toward potential therapeutics for C. difficile infections.
Authors: Jeffrey J Letourneau; Ilana L Stroke; David W Hilbert; Andrew G Cole; Laurie J Sturzenbecker; Brett A Marinelli; Jorge G Quintero; Joan Sabalski; Yanfang Li; Linh Ma; Igor Pechik; Philip D Stein; Maria L Webb Journal: Bioorg Med Chem Lett Date: 2018-10-29 Impact factor: 2.823
Authors: Kristina Oresic Bender; Megan Garland; Jessica A Ferreyra; Andrew J Hryckowian; Matthew A Child; Aaron W Puri; David E Solow-Cordero; Steven K Higginbottom; Ehud Segal; Niaz Banaei; Aimee Shen; Justin L Sonnenburg; Matthew Bogyo Journal: Sci Transl Med Date: 2015-09-23 Impact factor: 17.956