A combined quantum mechanical and molecular mechanical study of the reaction mechanism and alpha-amino acidity in alanine racemase

Combined quantum mechanical/molecular mechanical simulations have been carried out to investigate the origin of the carbon acidity enhancement in the alanine racemization reaction catalyzed by alanine racemase (AlaR). The present study shows that the enhancement of carbon acidity of alpha-amino acids by the cofactor pyridoxal 5'-phosphate (PLP) with an unusual, unprotonated pyridine is mainly due to solvation effects, in contrast to the intrinsic electron-withdrawing stabilization by the pyridinium ion to form a quinonoid intermediate. Alanine racemase further lowers the alpha-proton acidity and provides an overall 14-17 kcal/mol transition-state stabilization. The second key finding of this study is that the mechanism of racemization of an alanine zwitterion in water is altered from an essentially concerted process to a stepwise reaction by formation of an external aldimine adduct with the PLP cofactor. Finally, we have used a centroid path integral method to determine the intrinsic kinetic isotope effects for the two proton abstraction reactions, which are somewhat greater than the experimental estimates.     

A structure-based design approach for the identification of novel inhibitors: application to an alanine racemase

We report a new structure-based strategy for the identification of novel inhibitors. This approach has been applied to Bacillus stearothermophilus alanine racemase (AlaR), an enzyme implicated in the biosynthesis of the bacterial cell wall. The enzyme catalyzes the racemization of l- and d-alanine using pyridoxal 5-phosphate (PLP) as a cofactor. The restriction of AlaR to bacteria and some fungi and the absolute requirement for d-alanine in peptidoglycan biosynthesis make alanine racemase a suitable target for drug design. Unfortunately, known inhibitors of alanine racemase are not specific and inhibit the activity of other PLP-dependent enzymes, leading to neurological and other side effects.This article describes the development of a receptor-based pharmacophore model for AlaR, taking into account receptor flexibility (i.e. a `dynamic' pharmacophore model). In order to accomplish this, molecular dynamics (MD) simulations were performed on the full AlaR dimer from Bacillus stearothermophilus (PDB entry, 1sft) with a d-alanine molecule in one active site and the non-covalent inhibitor, propionate, in the second active site of this homodimer. The basic strategy followed in this study was to utilize conformations of the protein obtained during MD simulations to generate a dynamic pharmacophore model using the property mapping capability of the LigBuilder program. Compounds from the Available Chemicals Directory that fit the pharmacophore model were identified and have been submitted for experimental testing.The approach described here can be used as a valuable tool for the design of novel inhibitors of other biomolecular targets.     

Role of the pyridine nitrogen in pyridoxal 5'-phosphate catalysis: activity of three classes of PLP enzymes reconstituted with deazapyridoxal 5'-phosphate

Pyridoxal 5'-phosphate (PLP; vitamin B(6))-catalyzed reactions have been well studied, both on enzymes and in solution, due to the variety of important reactions this cofactor catalyzes in nitrogen metabolism. Three functional groups are central to PLP catalysis: the C4' aldehyde, the O3' phenol, and the N1 pyridine nitrogen. In the literature, the pyridine nitrogen has traditionally been assumed to be protonated in enzyme active sites, with the protonated pyridine ring providing resonance stabilization of carbanionic intermediates. This assumption is certainly correct for some PLP enzymes, but the structures of other active sites are incompatible with protonation of N1, and, consequently, these enzymes are expected to use PLP in the N1-unprotonated form. For example, aspartate aminotransferase protonates the pyridine nitrogen for catalysis of transamination, while both alanine racemase and O-acetylserine sulfhydrylase are expected to maintain N1 in the unprotonated, formally neutral state for catalysis of racemization and β-elimination. Herein, kinetic results for these three enzymes reconstituted with 1-deazapyridoxal 5'-phosphate, an isosteric analogue of PLP lacking the pyridine nitrogen, are compared to those for the PLP enzyme forms. They demonstrate that the pyridine nitrogen is vital to the 1,3-prototropic shift central to transamination, but not to reactions catalyzed by alanine racemase or O-acetylserine sulfhydrylase. Not all PLP enzymes require the electrophilicity of a protonated pyridine ring to enable formation of carbanionic intermediates. It is proposed that modulation of cofactor electrophilicity plays a central role in controlling reaction specificity in PLP enzymes.     

Alanine racemase free energy profiles from global analyses of progress curves

Free energy profiles for alanine racemase from Bacillus stearothermophilus have been determined at pH 6.9 and 8.9 from global analysis of racemization progress curves. This required a careful statistical design due to the problems in finding the global minimum in mean square for a system with eight adjustable parameters (i.e., the eight rate constants that describe the stepwise chemical mechanism). The free energy profiles obtained through these procedures are supported by independent experimental evidence: (1). steady-state kinetic constants, (2). solvent viscosity dependence, (3). spectral analysis of reaction intermediates, (4). equilibrium overshoots for progress curves measured in D(2)O, and (5). the magnitudes of calculated intrinsic kinetic isotope effects. The free energy profiles for the enzyme are compared to those of the uncatalyzed and the PLP catalyzed reactions. At pH 6.9, PLP lowers the free energy of activation for deprotonation by 8.4 kcal/mol, while the inclusion of apoenzyme along with PLP additionally lowers it by 11 kcal/mol.     

Kinetic isotope effects of L-Dopa decarboxylase

A mixed centroid path integral and free energy perturbation method (PI-FEP/UM) has been used to investigate the primary carbon and secondary hydrogen kinetic isotope effects (KIEs) in the amino acid decarboxylation of L-Dopa catalyzed by the enzyme L-Dopa decarboxylase (DDC) along with the corresponding uncatalyzed reaction in water. DDC is a pyridoxal 5'-phosphate (PLP) dependent enzyme. The cofactor undergoes an internal proton transfer between the zwitterionic protonated Schiff base configuration and the neutral hydroxyimine tautomer. It was found that the cofactor PLP makes significant contributions to lowering the decarboxylation barrier, while the enzyme active site provides further stabilization of the transition state. Interestingly, the O-protonated configuration is preferred both in the Michaelis complex and at the decarboxylation transition state. The computed kinetic isotope effects (KIE) on the carboxylate C-13 are consistent with that observed on decarboxylation reactions of other PLP-dependent enzymes, whereas the KIEs on the α carbon and secondary proton, which can easily be validated experimentally, may be used as a possible identification for the active form of the PLP tautomer in the active site of DDC.     

Covalent catalysis by pyridoxal: evaluation of the effect of the cofactor on the carbon acidity of glycine

First-order rate constants for deprotonation of the alpha-imino carbon of the adduct between 5'-deoxypyridoxal (1) and glycine were determined as the rate constants for Claisen-type addition of glycine to 1 where deprotonation is rate determining for product formation. There is no significant deprotonation at pH 7.1 of the form of the 1-glycine iminium ion with the pyridine nitrogen in the basic form. The value of kHO for hydroxide ion-catalyzed deprotonation of the alpha-imino carbon increases from 7.5 x 10(2) to 3.8 x 10(5) to 3.0 x 10(7) M(-1) s(-1), respectively, with protonation of the pyridine nitrogen, the phenoxide oxyanion, and the carboxylate anion of the 1-glycine iminium ion. There is a corresponding decrease in the pKas for deprotonation of the alpha-imino carbon from 17 to 11 to 6. It is proposed that enzymes selectively bind and catalyze the reaction of the iminium ion with pKa = 17. A comparison of kB = 1.7 x 10(-3) s(-1) for deprotonation of the alpha-imino carbon of this cofactor-glycine adduct (pKa = 17 by HPO4(2-) with k(cat)/K(m) = 4 x 10(5) M(-1) s(-1) for catalysis of amino-acid racemization by alanine racemase shows that the enzyme causes a ca 2 x 10(8)-fold acceleration of the rate of deprotonation the alpha-imino carbon. This corresponds to about one-half of the burden borne by alanine racemase in catalysis of deprotonation of alanine.     

Quantum Mechanics/Molecular Mechanics Studies on the Mechanism of Action of Cofactor Pyridoxal 5′‐Phosphate in Ornithine 4,5‐Aminomutase

A computational study was performed on the experimentally elusive cyclisation step in the cofactor pyridoxal 5′‐phosphate (PLP)‐dependent D ‐ornithine 4,5‐aminomutase (OAM)‐catalysed reaction. Calculations using both model systems and a combined quantum mechanics/molecular mechanics approach suggest that regulation of the cyclic radical intermediate is achieved through the synergy of the intrinsic catalytic power of cofactor PLP and the active site of the enzyme. The captodative effect of PLP is balanced by an enzyme active site that controls the deprotonation of both the pyridine nitrogen atom (N1) and the Schiff‐base nitrogen atom (N2). Furthermore, electrostatic interactions between the terminal carboxylate and amino groups of the substrate and Arg297 and Glu81 impose substantial “strain” energy on the orientation of the cyclic intermediate to control its trajectory. In addition the “strain” energy, which appears to be sensitive to both the number of carbon atoms in the substrate/analogue and the position of the radical intermediates, may play a key role in controlling the transition of the enzyme from the closed to the open state. Our results provide new insights into several aspects of the radical mechanism in aminomutase catalysis and broaden our understanding of cofactor PLP‐dependent reactions.     

A model for enzyme-substrate interaction in alanine racemase.

We report on a theoretical model for the complex of the enzyme alanine racemase with its natural substrate (L-alanine) and cofactor (pyridoxal 5'-phosphate). Electrostatic potentials were calculated and ionization states were predicted for all of the ionizable groups in alanine racemase. Some rather unusual charge states were predicted for certain residues. Tyr265' has an unusually low predicted pK(a) of 7.9 and at pH 7.0 has a predicted average charge of -0.37, meaning that 37% of the Tyr265' residues in an ensemble of enzyme molecules are in the phenolate form. At pH 8-9, the majority of Tyr265' side groups will be in the phenolate form. This lends support to the experimental evidence that Tyr265' is the catalytic base involved in the conversion of L-alanine to D-alanine. Residues Lys39 and Lys129 have predicted average charges of +0.91 and +0.14, respectively, at pH 7.0. Lys39 is believed to be the catalytic base for the conversion of D-alanine to L-alanine, and the present results show that, at least some of the time, it is in the unprotonated amine form and thus able to act as a base. Cys311', which is located very close to the active site, has an unusually low predicted pK(a) of 5.8 and at pH 7.0 has a predicted average charge of -0.72. The very low predicted charge for Lys129 is consistent with experimental evidence that it is carbamylated, since an unprotonated amine group is available to act as a Lewis base and form the carbamate with CO(2). Repeating the pK(a) calculations on the enzyme with Lys129 in carbamylated form predicts trends similar to those of the uncarbamylated enzyme. It appears that the enzyme has the ability to stabilize negative charge in the region of the active site. Implications for selective inhibitor design are discussed.     

Rapid photodynamics of vitamin B6 coenzyme pyridoxal 5'-phosphate and its Schiff bases in solution

The active form of vitamin B6, pyridoxal 5'-phosphate (PLP), is an important cofactor for numerous enzymes in amine and amino acid metabolism. Presented here is the first femtosecond transient absorption study of free PLP and two Schiff bases, PLP-valine and PLP-alpha-aminoisobutyric acid (AIB), in solution. Photoexcitation of free PLP leads to efficient triplet formation with an internal conversion rate that increases with increasing pH. The measured excited-state kinetics of the PLP-valine Schiff base exhibits a dramatic deuterium dependence as a result of excited-state proton transfer (ESPT) of the Calpha hydrogen in the amino acid substrate. This is consistent with formation of the key reaction carbanionic intermediate (quinonoid), which is resonance stabilized by the electron-deficient, conjugated pi system of the Schiff base/pyridine ring. The transient absorption signals of the PLP-Schiff base with alpha-methylalanine (2-aminoisobutyric acid), which does not have a Calpha proton, lack an observable deuterium effect, verifying ESPT formation of the quinonoid intermediate. In contrast to previous studies, no dependence on the excitation wavelength of the femtosecond kinetics is observed with PLP or PLP-valine, which suggests that a rapid (<250 fs) tautomerization occurs between the enolimine (absorbing at 330 nm) and ketoenamine (absorbing at 410 nm) tautomers in solution.     

A QM/MM study of the racemization of vinylglycolate catalyzed by mandelate racemase enzyme.

The experimentally postulated mechanism for the interconversion between (S)-vinylglycolate and (R)-vinylglycolate catalyzed by mandelate racemase enzyme consists of a two-step quite symmetric process through a dianionic enolic intermediate that is formed after the abstraction of the alpha-proton of vinylglycolate by a basic enzymatic residue and is then reprotonated by another residue. The challenging problem behind this reaction is how the enzyme manages to stabilize such an intermediate, that is, how it lowers enough the high pK(a) of the alpha-proton for the reaction to take place. The QM/MM simulations performed in this paper indicate that catalysis is based on the stabilization of the negative charge developed on the substrate along the reaction. We have identified three different reaction mechanisms starting from different quasi-degenerate structures of the substrate-enzyme complex. In two of them the stabilizing role is done by means of a catalytic proton transfer that avoids the formation of a dianionic intermediate, and they involve six steps instead of the two experimentally proposed. On the contrary, the third mechanism passes through a dianionic species stabilized by the concerted approach of a protonated enzymatic residue during the proton abstraction. The potential energy barriers theoretically found along these mechanisms are qualitatively in good agreement with the experimental free energy barriers determined for racemization of vinylglycolate and mandelate. The theoretical study of the effect of the mutation of Glu317 by Gln317 in the kinetics of the reaction reveals the important role in the catalysis of the hydrogen bond formed by Glu317 in the native enzyme, as only one of the mechanisms, the slower one, is able to produce the racemization in the active site of the mutant. However, we have found that this hydrogen bond is not an LBHB within our model.     

Functional Models of the Nickel Pincer Nucleotide Cofactor of Lactate Racemase

A novel nickel pincer cofactor was recently discovered in lactate racemase. Reported here are three synthetic nickel pincer complexes that are both structural and functional models of the pincer cofactor in lactate racemase. DFT computations suggest the ipso‐carbon atom of the pyridinium pincer ligands act as a hydride acceptor for lactate isomerization, whereas an organometallic pathway involving nickel‐mediated β‐hydride elimination is less favored.     

Theoretical studies on pyridoxal 5'-phosphate-dependent transamination of alpha-amino acids

Density functional methods have been applied to investigate the irreversible transamination between glyoxylic acid and pyridoxamine analog and the catalytic mechanism for the critical [1,3] proton transfer step in aspartate aminotransferase (AATase). The results indicate that the catalytic effect of pyridoxal 5'-phosphate (PLP) may be attributed to its ability to stabilize related transition states through structural resonance. Additionally, the PLP hydroxyl group and the carboxylic group of the amino acid can shuttle proton, thereby lowering the barrier. The rate-limiting step is the tautomeric conversion of the aldimine to ketimine by [1,3] proton transfer, with a barrier of 36.3 kcal/mol in water solvent. A quantum chemical model consisting 142 atoms was constructed based on the crystal structure of the native AATase complex with the product L-glutamate. The electron-withdrawing stabilization by various residues, involving Arg386, Tyr225, Asp222, Asn194, and peptide backbone, enhances the carbon acidity of 4'-C of PLP and Calpha of amino acid. The calculations support the proposed proton transfer mechanism in which Lys258 acts as a base to shuttle a proton from the 4'-C of PLP to Calpha of amino acid. The first step (proton transfer from 4'-C to lysine) is shown to be the rate-limiting step. Furthermore, we provided an explanation for the reversibility and specificity of the transamination in AATase.     

Metal ion inhibition of nonenzymatic pyridoxal phosphate catalyzed decarboxylation and transamination.

Nonenzymatic pyridoxal phosphate (PLP) catalyzed decarboxylations and transaminations have been revisited experimentally. Metal ions are known to catalyze a variety of PLP-dependent reactions in solution, including transamination. It is demonstrated here that the rate accelerations previously observed are due solely to enhancement of Schiff base formation under subsaturating conditions. A variety of metal ions were tested for their effects on the reactivity of the 2-methyl-2-aminomalonate Schiff bases. All were found to have either no effect or a small inhibitory one. The effects of Al(3+) were studied in detail with the Schiff bases of 2-methyl-2-aminomalonate, 2-aminoisobutyrate, alanine, and ethylamine. The decarboxylation of 2-methyl-2-aminomalonate is unaffected by metalation with Al(3+), while the decarboxylation of 2-aminoisobutyrate is inhibited 125-fold. The transamination reaction of ethylamine is 75-fold slower than that of alanine. Ethylamine transamination is inhibited 4-fold by Al(3+) metalation, while alanine transamination is inhibited only 1.3-fold. Metal ion inhibition of Schiff base reactivity suggests a simple explanation for the lack of known PLP dependent enzymes that make direct mechanistic use of metal ions. A comparison of enzyme catalyzed, PLP catalyzed, and uncatalyzed reactions shows that PLP dependent decarboxylases are among the best known biological rate enhancers: decarboxylation occurs 10(18)-fold faster on the enzyme surface than it does free in solution. PLP itself provides the lion's share of the catalytic efficiency of the holoenzyme: at pH 8, free PLP catalyzes 2-aminoisobutyrate decarboxylation by approximately 10(10)-fold, with the enzyme contributing an additional approximately 10(8)-fold.     

Identification of iron(III) peroxo species in the active site of the superoxide reductase SOR from Desulfoarculus baarsii

The active site of superoxide reductase SOR consists of an Fe2+ center in an unusual [His4 Cys1] square-pyramidal geometry. It specifically reduces superoxide to produce H2O2. Here, we have reacted the SOR from Desulfoarculus baarsii directly with H2O2. We have found that its active site can transiently stabilize an Fe3+-peroxo species that we have spectroscopically characterized by resonance Raman. The mutation of the strictly conserved Glu47 into alanine results in a stabilization of this Fe3+-peroxo species, when compared to the wild-type form. These data support the hypothesis that the reaction of SOR proceeds through such a Fe3+-peroxo intermediate. This also suggests that Glu47 might serve to help H2O2 release during the reaction with superoxide.     

A Combined Effect of Thermal and Enzymatic Racemization for d‐Aspartic Acid to l‐Aspartic Acid by High‐Performance Liquid Chromatography Assay

The racemization of d ‐aspartic acid to l ‐aspartic acid has been successfully performed with a coupled enzyme system at 90 °C and a pH of about 4.0 by the assay of high‐performance liquid chromatography. This coupled enzymatic racemization is a successive two‐step reaction first induced by d ‐amino acid oxidase and a subsequent coupled reaction by an aminotransferase clonezyme with the help of coenzyme pyridoxal 5′‐phosphate and cosubstrate l ‐glutamate. Due to the very high temperature, part of the l ‐aspartic acid is produced by the thermal effect. In fact the thermal racemization for aspartic acid can proceed from either d ‐ or l ‐aspartic acid via an intermediate fumaric acid and leads to the formation of d ,l ‐malic acid. The formation of α‐oxalacetic acid formed irreversibly from d ‐aspartic acid with d ‐amino acid oxidase can induce a side reaction to l ‐alanine. The thermal effect may also be responsible for the production of d ‐, and l ‐alanine.     

Reaction mechanism of the mandelate anion racemization catalyzed by mandelate racemase enzyme: a QM/MM molecular dynamics free energy study

The present work studies the reaction mechanism of the racemization of mandelate substrate by mandelate racemase enzyme. The reaction has some intriguing aspects such as the deprotonation of a nonacid hydrogen and the achievement of the pseudosymmetry necessary to obtain the racemic mixture. We will make use of a QM/MM potential energy surface to compute the free energy profiles associated with the reaction. The most favorable reaction mechanism consists of two proton transfers and the configuration inversion of the stereogenic carbon taking place in a concerted manner. We have also designed a suitable reaction coordinate to compute the free energy profiles for this rather complicated reaction. In addition, analysis of the electrostatic effects and bond distances along the reaction will explain how the enzyme accomplishes the catalysis. Finally, the enzymatic reaction will be compared to a model of the uncatalyzed reaction and the catalytic effect of mandelate racemase will be evaluated.     

Intrinsic primary and secondary hydrogen kinetic isotope effects for alanine racemase from global analysis of progress curves

The pyridoxal phosphate dependent alanine racemase catalyzes the interconversion of L- and D-alanine. The latter is an essential component of peptidoglycan in cell walls of Gram-negative and -positive bacteria, making alanine racemase an attractive target for antibacterials. Global analysis of protiated and deuterated progress curves simultaneously enables determination of intrinsic kinetic and equilibrium isotope effects for alanine racemase. The intrinsic primary kinetic isotope effects for Calpha hydron abstraction are 1.57 +/- 0.05 in the D --> L direction and 1.66 +/- 0.09 in the L --> D direction. Secondary kinetic isotope effects were found for the external aldimine formation steps in both the L --> D (1.13 +/- 0.05, forward; 0.90 +/- 0.03, reverse) and D --> L (1.13 +/- 0.06, forward; 0.89 +/- 0.03, reverse) directions. The secondary equilibrium isotope effects calculated from these are 1.26 +/- 0.07 and 1.27 +/- 0.07 for the L --> D and D --> L directions, respectively. These equilibrium isotope effects imply substantial ground-state destabilization of the C-H bond via hyperconjugation with the conjugated Schiff base/pyridine ring pi system. The magnitudes of the intrinsic primary kinetic isotope effects, the lower boundary on the energy of the quinonoid intermediate, and the protonation states of the active site catalytic acids/bases (K39-epsilonNH2 and Y265-OH) suggest that the pKa of the substrate Calpha-H bond in the external aldimine lies between those of the two catalytic bases, such that the proton abstraction transition state is early in the D --> L direction and late in the L --> D direction.     

On the ionization state of the substrate in the active site of glutamate racemase. A QM/MM study about the importance of being zwitterionic

Computer simulations on a QM/MM potential energy surface have been carried out to gain insights into the catalytic mechanism of glutamate racemase (MurI). Understanding such a mechanism is a challenging task from the chemical point of view because it involves the deprotonation of a low acidic proton by a relatively weak base to give a carbanionic intermediate. First, we have examined the dependency of the kinetics and thermodynamics of the racemization process catalyzed by MurI on the ionization state of the substrate (glutamate) main chain. Second, we have employed an energy decomposition procedure to study the medium effect on the enzyme-substrate electrostatic and polarization interactions along the reaction. Importantly, the present theoretical results quantitatively support the mechanistic proposal by Rios et al. [J. Am. Chem. Soc. 2000, 122, 9373-9385] for the PLP-independent amino acid racemases.     

Large pKa shifts of alpha-carbon acids induced by copper(II) complexes

A series of synthetic receptors (4-6) incorporating metal ions, specifically copper(II), were examined for their ability to enhance the acidity of active methylene compounds. The copper(II) complexes were observed to reduce the pKa of 1,3-diketone carbon acids in acetonitrile by as much as 12 pKa units. The relatively large pKa reduction achieved by the complex is attributed to the electrostatic interaction between the anionic pi system of the enolate and the copper(II) ions. The cage structure and hydrogen bonding sites in receptors 4 and 5 lead to a very modest further enhancement of the acidity relative to that with 6. This study provides insight into the way in which metalloenzymes stabilize an enolate intermediate.     

In search of radical intermediates in the reactions catalyzed by lysine 2,3-aminomutase and lysine 5,6-aminomutase

High-level ab initio calculations have been used to study radical intermediates in the reactions catalyzed by lysine 2,3-aminomutase (2,3-LAM) and lysine 5,6-aminomutase (5,6-LAM). The reactions of these enzymes with the substrate analogues 4-oxalysine (4-OL), 4-thialysine (4-TL), or trans-4,5-dehydrolysine (t-4,5-DL) are rationalized in terms of stabilization provided by the substituent to the adjacent radical center. Large changes in the exothermicity accompanying the initial H-abstraction are observed relative to the lysine reference values that follow the series 4-OL < 4-TL < t-4,5-DL. These changes have the primary effect of increasing the endothermicity for subsequent ring-closure to form the putative aziridinylcarbinyl radical intermediate. Such stabilization is consistent with experimental observations of the substrate-derived radical (S*) in the reaction of 2,3-LAM with 4-TL as well as the ability of t-4,5-DL to act as an irreversible inhibitor of 2,3-LAM. Our calculations suggest that 4-TL and trans-3,4-dehydrolysine may also permit experimental characterization of S* radicals in the reactions catalyzed by 5,6-LAM. Strategies for modifying PLP are presented that might lead to the first observation of the aziridinylcarbinyl radical intermediate (I*) in the aminomutase-catalyzed reactions.