|Home | About | Journals | Submit | Contact Us | Français|
Benzoylformate decarboxylase (BFDC) is a thiamin diphosphate (ThDP)-dependent enzyme acting on aromatic substrates. In addition to its metabolic role in the mandelate pathway, BFDC shows broad substrate specificity coupled with tight stereocontrol in the carbon-carbon bond-forming reverse reaction, making it a useful biocatalyst for the production of chiral α-hydroxyketones. The reaction of methyl benzoylphosphonate (MBP), an analogue of the natural substrate benzoylformate, with BFDC results in the formation of a stable analogue (C2α-phosphonomandelylThDP) of the covalent ThDP-substrate adduct C2α-mandelylThDP. Formation of the stable adduct is confirmed both by formation of a circular dichroism band characteristic of the 1’,4’-iminopyrimidine tautomeric form of ThDP (commonly observed when ThDP forms tetrahedral complexes with its substrates), and by high resolution mass spectrometry of the reaction mixture. In addition, the structure of BFDC with the MBP inhibitor was solved by X-ray crystallography, to a spatial resolution of 1.37 Å. The electron density clearly shows formation of a tetrahedral adduct between the C2 atom of ThDP and the carbonyl carbon atom of the MBP. This adduct resembles the intermediate from the penultimate step of the carboligation reaction between benzaldehyde and acetaldehyde. The combination of real-time kinetic information via stopped-flow circular dichroism with steady-state data from equilibrium circular dichroism measurements and X-ray crystallography reveals details of the first step of the reaction catalyzed by BFDC. The MBP-ThDP adduct on BFDC is compared to the recently solved structure of the same adduct on benzaldehyde lyase, another ThDP-dependent enzyme capable of catalyzing aldehyde condensation with high stereospecificity.
Characterizing the discrete steps of a chemical reaction remains an exceptionally difficult task. One strategy for enzyme-catalyzed reactions has been to obtain the structure of the protein catalyst in the presence of a reaction intermediate. In most cases, though, structural determination is relatively slow and the lifetime of relevant intermediates fleeting. The use of molecules that resemble the intermediate but that are longer lasting is an established method for obtaining detailed information on the molecular species formed during an enzyme-catalyzed reaction pathway(1–6).
Benzoylformate decarboxylase (BFDC, EC 126.96.36.199) catalyzes the decarboxylation of benzoylformate to benzaldehyde (Scheme 1a). A thiamin diphosphate (ThDP)-dependent enzyme, BFDC plays a fundamental role in mandelate metabolism(7, 8). One of several ThDP-dependent enzymes that catalyze the decarboxylation of 2-oxo acids, it has received considerable interest in studies aimed at understanding the substrate specificity in enzymes of this class(9–12), because of the enzyme’s ability to catalyze stereospecific carboligation reactions(13–16) (Scheme 1b–c).
BFDC was incubated with methyl benzoylphosphonate (MBP, Scheme 2), an analogue of benzoylformate, forming an MBP-ThDP adduct that cannot be decarboxylated(17–21). The formation of a covalent MBP-ThDP adduct on BFDC was initially followed by time-resolved CD spectroscopy, and its formation was confirmed by mass spectrometry. Finally, the X-ray structure of the MBP-ThDP adduct on BFDC was determined and compared with that of the same adduct on benzaldehyde lyase (BAL), a related ThDP-dependent enzyme also acting on aromatic substrates(22). BAL has also been investigated for its utility in producing stereospecific 2-hydroxy ketones(10, 16, 23–25). Although both enzymes form (R)-benzoin in > 99% ee, BFDC catalyzes formation of 2-(S)-hydroxypropiophenone (2-HPP) from benzaldehyde and acetaldehyde in 90% ee. Conversely, BAL produces 2-(R)-HPP in 99% ee. The causes of this difference in chiral induction have been the subject of some speculation(9, 10). Here we examine, for the first time, structural evidence for the source of the observed stereospecificity on BFDC and BAL.
Methyl benzoylphosphonate was synthesized as reported elsewhere(22). Benzoylformate, NADH, and horse liver alcohol dehydrogenase (HLADH) were purchased from Sigma-Aldrich (St. Louis, MO). All other chemicals purchased were of the highest purity grade available and were used without further purification.
The enzyme was expressed, purified and stored as described elsewhere(26). The activity of BFDC was measured using a NADH/horse liver alcohol dehydrogenase (HLADH) coupled assay(26). In a typical assay, 1 mL of the reaction mixture contained in 50 mM KH2PO4 (pH 6.0): benzoylformate (3.5 mM), MgSO4 (2.5 mM), ThDP (0.5 mM), NADH (0.28 mM) and 0.25 mg (10 units) of HLADH. The reaction was initiated by the addition of BFDC. The disappearance of NADH was monitored at 340 nm at 30 °C.
The inhibition constant (Ki) for methyl benzoylphosphonate was determined by varying the concentrations of benzoylformate (0.3–3×Km) and methyl benzoylphosphonate (0.3–4×Ki). Initial velocities were then fit to eq 1 using the Enzyme Kinetics package of SigmaPlot 9.0 (Systat Software Inc).
To the BFDC (1.7 mg/mL, concentration of active centers = 30.2 µM) in 20 mM KH2PO4 (pH 6.0), containing MgCl2 (2.5 mM) and ThDP (0.50 mM) aliquots of MBP (0.001– 0.50 mM) or MAP (1–30 mM) were added and CD spectra were recorded in the near-UV region using an Aviv Circular Dichroism Spectrometer (Model 202). Difference spectra were obtained on subtraction of the spectrum of BFDC in the absence of MBP using SigmaPlot 2001 (v.7.101). Values of Kd,MBP were calculated by fitting the data to a Hill equation: CD = (CDmax * [MBP]n) / (Kd n +[MBP]n).
Experiments were carried out on a π*-180 CDF Spectrometer from Applied Photophysics. In a typical experiment, BFDC (83 µM active-site concentration) in 50 mM KH2PO4 (pH 6.0) containing 2.5 mM MgSO4 and 0.5 mM ThDP (buffer A) was mixed with an equal volume of substrate of intended concentration in the same buffer. The reaction was monitored for varied time (0.01– 200 s) and data points were collected at varied time intervals (2.5–10 ms). A slit width of 2 mm and a path length of 10 mm were used. The temperature was maintained at 30 °C. The data were analyzed using Pro-K Global Analysis software. SigmaPlot 2001 version 7.101 was used for data analysis.
BFDC was concentrated to 30–40 mg/mL in its storage buffer (100 mM potassium phosphate pH 6.0, 1 mM MgSO4, 0.5 mM ThDP, 10% glycerol) and could be stored at −80 °C for up to several months prior to crystallization. Under the conditions of Hasson for hanging drop vapor diffusion, the protein crystallized from heavy precipitate to give very large, well-formed prisms in a matter of days(27). A 1 M solution of MBP was utilized at a 1:50 dilution for co-crystallization studies. The MBP was diluted 1:5 into protein storage buffer, and then incubated 1:10 with protein solution on ice for 10–15 min. 1 µL of this BFDC-MBP solution was mixed with 1 µL of well solution [100 mM TrisCl (pH 8.5), 150 mM CaCl2, 0.5% (v/v) MPD, 22% (v/v) PEG 400] and left to stand over the well solution in a hanging drop vapor diffusion tray. As with protein alone, large crystals formed in 3–5 days. After transfer to fresh well solution containing 2 mM MBP and incubation for approximately 60 s, the crystals were transferred again to a cryoprotectant well solution made up with 2 mM MBP and 30% v/v glycerol and flash-frozen in liquid nitrogen.
Crystals were irradiated at the Advanced Photon Source at Argonne National Laboratory (Argonne, IL), using the 23-IDD beamline administered by GM/CA-CAT. The data were scaled in I222, the same space group previously reported for wild-type BFDC in the absence of inhibitor. The HKL 2000 software package and the CCP4 suite of programs were used for data reduction and processing(28, 29). A search model for molecular replacement was obtained from the crystal structure of Hasson (PDB ID: 1BFD)(27). As was the case for BFDC itself, the asymmetric unit of the BFDC-MBP complex was determined to contain a single monomer. Subsequent refinement of the model was carried out as described below, and crystallization and refinement statistics are given in Table 1.
The initial model resulting from molecular replacement was refined using Refmac5(30). Placement of water molecules into difference density peaks with an intensity greater than 3σ was carried out by the water-picking algorithm in Coot 0.4.1(31). Individual water molecules were checked for spherical electron density and the existence of at least one appropriately positioned hydrogen-bonding partner. Water molecules in the vicinity of the active site were removed. After refinement with a model including all non-active site water molecules, placement of the ThDP cofactor into the positive difference density proceeded, again using Coot. The model containing non-active site water molecules, the polypeptide chain and the ThDP cofactor was refined to generate an omit map for the MBP. Coordinate files and appropriate restraints for the unreacted MBP, as well as both enantiomers of the likely adduct were generated using the Sketcher program of the CCP4 suite and the PRODRG server(29, 33). Given the fairly high spatial resolution of the structure, the best-fitting species was unequivocally the MBP-ThDP adduct (see Discussion). In order to independently cross-validate the model, the structure was independently solved using the automated molecular replacement routine of Phenix1.3b(32), starting with the model from 1BFD. The simulated-annealing omit map generated in Phenix was in excellent agreement with that produced from CCP4. The final structural solution was subjected to validation using SFCheck and ProCheck 3.4.4, deposited in the Protein Data Bank, and assigned a PDB ID of xxxx(34, 35).
The steady-state kinetic data in Figure 1 show that MBP, a close structural and electrostatic analog of the true substrate benzoylformate, is a competitive inhibitor of BFDC. The data provided a Ki value for MBP of 0.38±0.04 mM. By comparison, the Km value for benzoylformate was found to be 0.76±0.09 mM.
Incremental addition of MBP to BFDC produced a positive CD band with λmax at 299 nm in the difference CD spectra of BFDC, corresponding to the 1’,4’-iminopyrimidine tautomer of C2α-phosphonomandelyl-ThDP (PMThDP) (Figure 2A). The amplitude of this band increased with increasing concentration of MBP and displayed saturation with a Kd,MBP = 0.113 mM (average from three independent experiments; see Figure 2B for an example) as compared with Ki,MBP = 0.38±0.04 mM obtained from kinetics. The CD band at 299 nm was not detected on BFDC after the overnight dialysis, indicating that the formation of 1’,4’-imino PMThDP is reversible. Next, stopped-flow CD was used to determine the rate of formation of 1’,4’-imino PMThDP on BFDC (Scheme 1a and Figure 2C). A positive CD band at 305 nm developed with time, which followed an exponential growth to maximum with a rate constant of 2.08 s−1.
Methyl acetylphosphonate (MAP) is similar to MBP characterized above, but with a methyl group replacing the phenyl ring in MBP. A comparison of the behavior of MAP and MBP would throw light on the substrate specificity of BFDC. A CD titration of BFDC with MAP revealed the development of the same positive CD band at 299 nm (Figure 3A), a characteristic of 1’,4’-iminophosphonolactylThDP (PLThDP), the covalent adduct between MAP and ThDP, already reported on yeast pyruvate decarboxylase (YPDC) and the E1 component of E. coli pyruvate dehydrogenase(36–40). The Kd for MAP was 5.49 mM (Figure 3A inset), 55-times larger than that for MBP. Overnight dialysis again eliminated the CD band at 299 nm (Figure 3B), indicating that the formation of the PLThDP is reversible.
The formation of the adduct between ThDP and MBP on BFDC was confirmed by mass spectral analysis. The BFDC was titrated by MBP and formation of the positive CD band at 299 nm was monitored. Once the titration was completed, 12.5% TCA was added to precipitate the protein. The supernatant was analyzed for products by FTICR mass spectrometry. The species with mass of 625.06276 Da (M+) and mass of 626.06613 (M+1) were identified as PMThDP, confirming its formation on BFDC (Figure 4).
The structural consequence of addition of the ThDP ylide to the MBP is the formation of a tetrahedral center at the 2’ carbon (Scheme 2). The electron density after refinement of the structural model, including the cofactor, clearly shows tetrahedral positive difference density above the plane of the cofactor, along with a ring of density at an angle to the plane of the thiazolium ring (Figure 5). This electron density is relatively well fit by a tetrahedral adduct of ThDP and MBP, but not by the planar sp2-hybridized unbound MBP. In the initial refinement stages, the adduct C2α-phosphonomandelylThDP (PMThDP) was modeled independently of the cofactor. Given the covalent adduct PMThDP that is formed, subsequent rounds of refinement involved fitting PMThDP. The best fit is obtained when the C–C distance between C2 and C2α is somewhat longer than a typical C–C bond. Thus, it is likely that the observed electron density represents a mixture of species in the active site. Some positive difference density is observable near the phenyl ring and phosphonate group of bound MBP, again suggesting multiple conformations. This C–C distance between C2 and C2α is in contrast to the observed electron density for BAL with MBP, where the same distance was observed to be 1.5 Å(22). This difference may result from differences in experimental conditions. In any case, it is clear that the observed electron density is best fit by the covalent adduct between ThDP and MBP.
The compound MBP behaved as a mechanism-based, competitive, reversible inhibitor of BFDC (Ki = 0.38±0.04 µM) (Figure 1). The spectroscopic signal for the covalent adduct formed was characteristic of C2α-phosphonomandelylThDP, a stable MThDP analogue. Earlier, the use of MAP with ThDP-dependent enzymes YPDC and E coli PDHc E1 enabled the observation of a positive CD band in the 300–305 nm region(41, 42). Model studies showed(37, 43) that formation of 1’,4’-iminopyrimidine tautomer accounts for this electronic transition(37–40, 42, 44–46). Accordingly, the band with an absorption maximum at 299 nm (Figure 2A) was attributed to the 1’,4’-iminoPLThDP. Addition of both MBP and MAP produced the same positive band (λmax 299 nm; Figure 2 and Figure 3). A rate constant of 2.08 s−1 was obtained by stopped-flow CD of BFDC with MBP for formation of the band and, by extension, for covalent bond formation (Figure 2C).
The role of the 4’-aminopyrimidine ring in the catalytic cycle of ThDP-dependent enzymes has been a central point of interest for a long time(37–40, 42–45, 47–51). The MBP-ThDP adduct observed here (PMThDP) is a stable analogue of MThDP, so the observations above likely apply to reaction with the native substrate. The formation of the 1’,4’-imino tautomer of MThDP could be explained by the mechanism shown in Scheme 2. The 4’-aminopyrimidine ring undergoes protonation at N1’ by a conserved glutamate, forming 4’-aminopyrimidinium-ThDP, which loses a proton to form 1’,4’-imino-ThDP (IP form of the cofactor, λmax = 299 nm). The imino nitrogen then abstracts the proton from C2-H of thiazolium producing the ylide and the 4’-aminopyrimidiniumThDP (APH+). The ylide carries out a nucleophilic attack on the incoming α-keto acid at its carbonyl carbon, producing the C2α-oxyanion which picks up the proton to form the 1’,4’-iminoMThDP. We have suggested that the tetrahedral ThDP-bound intermediates on all ThDP dependent enzymes exist in their 1’,4’-imino tautomeric forms(38).
Several hypotheses have been put forward to explain how BFDC accelerates the rate of decomposition over the uncatalyzed, spontaneous decarboxylation of the C2α-mandelylthiamin diphosphate intermediate(52). The current thinking postulates that the primary source of this acceleration is not maintenance of the adduct in a hydrophobic environment, but rather positioning of acidic residues to protonate the carbanion that forms upon decarboxylation(53). These residues were proposed to be His70 and His281, on the grounds of proximity and that replacement of the imidazole ring of His70 and His281 can lead to significant reduction in kcat(54). The present crystal structure of the MBP adduct shows that the closest heavy-atom approach between the C2α anion and atoms of either histidine is > 5 Å (Figure 6). In addition, the catalyzed rate is a million-fold greater than the uncatalyzed rate, where replacement of His281 by alanine is accompanied by a roughly 200-fold rate reduction(54). Recent saturation mutagenesis experiments have thoroughly examined the sidechain tolerance of His70 and His281, as well as Ser26(26). Surprisingly, these protic active site residues can be replaced by hydrophobic amino acids with relatively modest reduction in catalytic efficiency (on the order of ten-fold), as measured by kcat/Km. Together with this mutational data, the current structure leads us to propose a potential role for promotion of decarboxylation by stabilization of a particular substrate geometry rather than proton transfer. Interaction of the substrate carboxylate with His281, for example, could position the carboxylate over the thiazolium ring, whereas interaction with Ser26 would place the plane of the carboxylate perpendicular to the thiazolium-C2α bond, and out of range of interaction with the thiazolium π system. These conformations may promote decarboxylation to a different extent or via a slightly different reaction coordinate. Future examination of the relative energetics of these conformations may provide insight into the role of polar active site residues in positioning benzoylformate for decarboxylation in the BFDC active site. The replacement of these residues with large hydrophobic sidechains, such as leucine or phenylalanine, may influence the localization of negative charge onto the two carboxylate oxygens, also an issue that lends itself to theoretical study.
Benzaldehyde lyase is a ThDP-dependent enzyme that carries out the reversible conversion of (R)-benzoin to benzaldehyde(21, 55–57) (Scheme 1). Although BFDC and BAL are homologous (27% sequence identity) few active site residues are conserved, and BFDC is essentially inactive toward benzoin(21, 56). Likewise, BAL has very little decarboxylase activity, although its A28S variant does exhibit significant BFDC activity(56). Nonetheless, the two enzymes carry out very similar carboligation reactions in that both condense two equivalents of benzaldehyde to give (R)-benzoin with very high enantioselectivity (>98% ee). However, one intriguing distinction between BFDC and BAL is that they give products with the opposite stereochemistry in a mixed acyloin condensation reaction. For example, BFDC condenses benzaldehyde and acetaldehyde to produce 2-(S)-hydroxypropiophenone (HPP) with >90% ee, whereas BAL generates the (R)-enantiomer with >99% ee(58).
Alignment of the two structures shows the differences in active site residues directly contacting the MBP (Figure 7). A notable distinction is the orientation of the phosphonate hydrogen bonding network in the BFDC active site relative to that of BAL(22). Where BFDC has the possibility for orienting the substrate toward Ser26 or toward His281, BAL supports only orientation toward Ala28 (the structural homologue of Ser26) and the ordered water positioned there.
The complex of BFDC and BAL with MBP mimics the result of addition of the donor benzaldehyde to acceptor acetaldehyde, the penultimate and stereogenic step of the carboligation reaction (Scheme 3). Taking the three phosphonate oxygens to represent three possible oxygen binding sites of the single oxygen of the reaction intermediate, considering both configurations for each conformation leads to six possibilities. Disregarding those that would give rise to 2-(R)-HPP (which is not observed) leaves three possible conformations to consider (Figure 8a). In addition to the likely hydrogen bonds indicated above, Figure 8b shows potential stabilizing hydrophobic interactions for a substrate methyl group (yellow lines showing heavy atom distances to nearest atom of the relevant hydrophobic residue)(9, 10). Stabilization of the carbinol oxygen by two sidechains, a proper fit into the small hydrophobic pocket of the methyl group, and the presence of His281 as a potential proton donor/acceptor support the likelihood of a conformation like that of panel B. Note that the MBP-ThDP adduct lacks the methyl group to nestle into the methyl-binding pocket. We would expect the best fit to result from a ~15° clockwise rotation around the C2α-P bond (viewed from above). This rotation brings the hydroxyl group into a position that preserves the interaction with Ser26 but loses that with His281 (Figure 8c). His281 is not a necessary residue for catalysis of 2-HPP formation. Indeed, mutation to alanine results in an approximately 125-fold increase in catalytic rate(15), perhaps because the removal of the imidazole sidechain favors a more productive binding mode wherein the methyl binding pocket and Ser26 dominate substrate positioning.
A similar analysis may be undertaken for condensation of benzaldehyde with acetaldehyde catalyzed by BAL, which gives the opposite enantiomer of HPP to that given by BFDC(58) (Figure 9). Of the three pro-R conformations available to BAL, D would seem to be preferred, as it is the only one of the three in which a hydrogen bonding interaction is available to the 2-hydroxyl group of the ThDP-bound HPP intermediate. This conformation also has the effect of positioning the methyl group of acetaldehyde in the more open portion of the active site cleft. The methyl binding pocket observed in BFDC is occluded in BAL by the γ-methyl group of Thr481, which perhaps sets up a steric clash that leads to preferential positioning of the aldehydic proton in this small pocket(9, 10). This explanation leads naturally to the prediction that the BAL T481S mutant may well produce 2-(S)-HPP, a prediction which is currently under investigation.
This analysis also serves to show why both enzymes produce (R)-benzoin. In the case of BAL, both the methyl group of acetaldehyde and the phenyl of benzaldehyde are too bulky to fit in the active site. So, BAL produces (R)-benzoin for the same reason it produces 2-(R)-HPP. In the case of BFDC, we argued that binding the methyl group of acetaldehyde in the methyl-binding pocket gave rise to the (S)-product. This pocket, however, fails to accommodate the phenyl ring of benzaldehyde. Now, when the acceptor benzaldehyde binds, it positions its aldehydic proton to enter the methyl binding pocket, the carbonyl oxygen to interact with Ser26, and the phenyl group to occupy the mouth of the active site cleft (Figure 8). This is a similar binding mode to that proposed for BAL and correspondingly generates the R-product.
The same explanation holds for recent studies showing that BFDC-catalyzed condensation of benzaldehyde with propanal overwhelmingly produces the R-product, while a point mutant of the enzyme gives the S-product(9). Increasing the space available in the active site by converting Leu461 to Ala permits the substrate to orient itself in such a way that the propane sidechain can be tucked into the created space, generating the S-product.
The work presented here seeks to characterize formation of the pre-decarboxyation covalent adduct between the inhibitor MBP and the ThDP cofactor of the enzyme BFDC. MBP was shown to be a mechanism-based inhibitor of BFDC, by determining the parameters of its inhibition of the enzyme and also by mass spectrometric analysis of the cofactor adduct isolated from the reaction mixture. In addition, the rate of formation the covalent adduct was followed in real time by stopped-flow CD, and the results were interpreted in light of the increasing understanding of ThDP activation in a variety of ThDP-dependent enzymes. Finally, the crystal structure of the covalent adduct was solved. The BFDC-MBP adduct is of particular interest, since an analysis of this protein-inhibitor complex may shed light on the detailed mechanism of BFDC. This enzyme has been employed as a biocatalyst for the asymmetric production of 2-hydroxyketones, as has the related enzyme BAL. The two enzymes catalyze the production of opposite enantiomers for the condensation of benzaldehyde with acetaldehyde. The source of this stereospecificity is examined in light of structures available for the covalent complex of both enzymes with the same MBP inhibitor. Altogether, this work provides a detailed analysis of the first step in the mechanism of benzoylformate decarboxylation and its functional consequences.
Supported by NIH NRSF F32GM069057 (GSB, Brandeis), NIH GM050380 (Rutgers) and NSF EF-0425719 (Brandeis, Michigan).