Internal proton transfer in the external pyridoxal 5'-phosphate Schiff base in dopa decarboxylase

1,3-Proton Transfer of Pyridoxal 5'-Phosphate Schiff Base in the Branched-Chain Aminotransferase: Concerted or Stepwise Mechanism?

The branched-chain aminotransferase from Mycobacterium tuberculosis (MtIlvE) is a pyridoxal 5'-phosphate (PLP) dependent enzyme, and it is essential for the synthesis of the branched-chain amino acids. Ketimine is an important intermediate in the catalytic process. We have investigated the mechanism of ketimine formation and the energy landscape using the multiple computational methods. It is found that the 1,3-proton transfer involved in ketimine formation occurs through a stepwise process rather than a one-step process. Lys204 is identified as a key residue for ligand binding and as a base that abstracts the Cα proton from the PLP-Glu Schiff base, yielding a carbanionic intermediate. The first proton transfer is the rate-limiting step with an energy barrier of 17.8 kcal mol-1. Our study disclosed the detailed pathway of the proton transfer from external aldimine to ketimine, providing novel insights into the catalytic mechanism of other PLP-dependent enzymes.

Free Energy Landscapes of the Tautomeric Interconversion of Pyridoxal 5'-Phosphate Aldimines at the Active Site of Ornithine Decarboxylase in Aqueous Media

The pyridoxal 5'-phosphate (PLP) acts as a coenzyme for a large number of biochemical reactions. It exists in mainly two bound forms at the active site of the concerned enzyme: the internal aldimine, in which the PLP is bound with the ϵ-amino group of lysine at the active site, and the external aldimine, where the PLP is bound to the substrate amino acid. Both the internal and external aldimines have Schiff base linkage (N-H-O) and can exist in two tautomeric structures of ketoenamine and enolimine forms. In this work, we have investigated the free energy landscape for the tautomeric proton transfer in the internal and external aldimines at the active site of the ornithine decarboxylase enzyme in an aqueous medium. We performed hybrid quantum-classical metadynamics and force field-based molecular dynamics simulations, which revealed that the ketoenamine tautomer is more stable than the enolimine form. The QM/MM metadynamics calculations show that the free energy difference between the ketoenamine and enolimine forms for the internal aldimine is 3.9 kcal/mol, and it is found to be 5.8 kcal/mol for the external aldimine, with the ketoenamine form being more stable in both cases. The results are further supported by calculations of the binding free energies from classical simulations and static quantum chemical calculations in different environments. We have also analyzed the configurational structure of the microenvironment at the active site in order to have better insights into the interactions of the active site residues with the PLP in its two tautomeric forms.

PLP-Dependent Enzyme Methionine γ-Lyase: Insights into the Michaelis Complex from Molecular Dynamics and Free Energy Simulations

Methionine γ-lyase (MGL) breaks down methionine, with the help of its cofactor pyridoxal-5'-phosphate (PLP), or vitamin B6. Methionine depletion is damaging for cancer cells but not normal cells, so MGL is of interest as a therapeutic protein. To increase our understanding and help engineer improved activity, we focused on the reactive, Michaelis complex M between MGL, covalently bound PLP, and substrate Met. M is not amenable to crystallography, as it proceeds to products. Experimental activity measurements helped exclude a mechanism that would bypass M . We then used molecular dynamics and alchemical free energy simulations to elucidate its structure and dynamics. We showed that the PLP phosphate has a pKa strongly downshifted by the protein, whether Met is present or not. Met binding affects the structure surrounding the reactive atoms. With Met, the Schiff base linkage between PLP and a nearby lysine shifts from a zwitterionic, keto form to a neutral, enol form that makes it easier for Met to approach its labile, target atom. The Met ligand also stabilizes the correct orientation of the Schiff base, more strongly than in simulations without Met, and in agreement with structures in the Protein Data Bank, where the Schiff base orientation correlates with the presence or absence of a co-bound anion or substrate analogue in the active site. Overall, the Met ligand helps organize the active site for the enzyme reaction by reducing fluctuations and shifting protonation states and conformational populations.

A Computational Model for the PLP-Dependent Enzyme Methionine γ-Lyase

Pyridoxal-5′-phosphate (PLP) is a cofactor in the reactions of over 160 enzymes, several of which are implicated in diseases. Methionine γ-lyase (MGL) is of interest as a therapeutic protein for cancer treatment. It binds PLP covalently through a Schiff base linkage and digests methionine, whose depletion is damaging for cancer cells but not normal cells. To improve MGL activity, it is important to understand and engineer its PLP binding. We develop a simulation model for MGL, starting with force field parameters for PLP in four main states: two phosphate protonation states and two tautomeric states, keto or enol for the Schiff base moiety. We used the force field to simulate MGL complexes with each form, and showed that those with a fully-deprotonated PLP phosphate, especially keto, led to the best agreement with MGL structures in the PDB. We then confirmed this result through alchemical free energy simulations that compared the keto and enol forms, confirming a moderate keto preference, and the fully-deprotonated and singly-protonated phosphate forms. Extensive simulations were needed to adequately sample conformational space, and care was needed to extrapolate the protonation free energy to the thermodynamic limit of a macroscopic, dilute protein solution. The computed phosphate pK_a was 5.7, confirming that the deprotonated, −2 form is predominant. The PLP force field and the simulation methods can be applied to all PLP enzymes and used, as here, to reveal fine details of structure and dynamics in the active site.

Characterization of a Novel Shewanella algae Arginine Decarboxylase Expressed in Escherichia coli

Arginine decarboxylase (ADC) catalyzes the decarboxylation of arginine to form agmatine, an important physiological and pharmacological amine, and attracts attention to the enzymatic production of agmatine. In this study, we for the first time overexpressed and characterized the marine Shewanella algae ADC (SaADC) in Escherichia coli. The recombinant SaADC showed the maximum activity at pH 7.5 and 40 °C. The SaADC displayed previously unreported substrate inhibition when the substrate concentration was higher than 50 mM, which was the upper limit of testing condition in other reports. In the range of 1-80 mM L-arginine, the SaADC showed the K_m, k_cat, K_i, and k_cat/K_m values of 72.99 ± 6.45 mM, 42.88 ± 2.63 s^-1, 20.56 ± 2.18 mM, and 0.59 s/mM, respectively, which were much higher than the K_m (14.55 ± 1.45 mM) and k_cat (12.62 ± 0.68 s^-1) value obtained by assaying at 1-50 mM L-arginine without considering substrate inhibition. Both the k_cat values of SaADC with and without substrate inhibition are the highest ones to the best of our knowledge. This provides a reference for the study of substrate inhibition of ADCs.

A Photoregulated Racemase Mimic

The racemase enzymes convert L-amino acids to their D-isomer. The reaction proceeds through a stepwise deprotonation-reprotonation mechanism that is assisted by a pyridoxal phosphate (PLP) coenzyme. This work reports a PLP-photoswitch-imidazole triad where the racemization reaction can be controlled by light by tweaking the distance between the basic residue and the reaction centre.

A QM/MM simulation study of transamination reaction at the active site of aspartate aminotransferase: Free energy landscape and proton transfer pathways

Transaminase is a key enzyme for amino acid metabolism, which reversibly catalyzes the transamination reaction with the help of PLP (pyridoxal 5^' -phosphate) as its cofactor. Here we have investigated the mechanism and free energy landscape of the transamination reaction involving the aspartate transaminase (AspTase) enzyme and aspartate-PLP (Asp-PLP) complex using QM/MM simulation and metadynamics methods. The reaction is found to follow a stepwise mechanism where the active site residue Lys258 acts as a base to shuttle a proton from α-carbon (CA) to imine carbon (C4A) of the PLP-Asp Schiff base. In the first step, the Lys258 abstracts the CA proton of the substrate leading to the formation of a carbanionic intermediate which is followed by the reprotonation of the Asp-PLP Schiff base at C4A atom by Lys258. It is found that the free energy barrier for the proton abstraction by Lys258 and that for the reprotonation are 17.85 and 3.57 kcal/mol, respectively. The carbanionic intermediate is 7.14 kcal/mol higher in energy than the reactant. Hence, the first step acts as the rate limiting step. The present calculations also show that the Lys258 residue undergoes a conformational change after the first step of transamination reaction and becomes proximal to C4A atom of the Asp-PLP Schiff base to favor the second step. The active site residues Tyr70* and Gly38 anchor the Lys258 in proper position and orientation during the first step of the reaction and stabilize the positive charge over Lys258 generated at the intermediate step.

Theoretical Evaluation of the Reaction Mechanism of Serine Hydroxymethyltransferase

Serine hydroxymethyltransferase (SHMT) is a pyridoxal phosphate (PLP)-dependent enzyme that catalyzes the reversible conversion of serine and tetrahydrofolate (THF) to glycine and 5,10-methylene THF. SHMT is a folate pathway enzyme and is therefore of considerable medical interest due to its role as an important intervention point for antimalarial, anticancer, and antibacterial treatments. Despite considerable experimental effort, the precise reaction mechanism of SHMT remains unclear. In this study, we explore the mechanism of SHMT to determine the roles of active site residues and the nature and the sequence of chemical steps. Molecular dynamics (MD) methods were employed to generate a suitable starting structure which then underwent analysis using hybrid quantum mechanical/molecular mechanical (QM/MM) simulations. The QM region consisted of 12 key residues, two substrates, and explicit solvent. Our results show that the catalytic reaction proceeds according to a retro-aldol synthetic process with His129 acting as the general base in the reaction. The rate-determining step involves the cleavage of the PLP-serine aldimine C_α-C_β bond and the formation of formaldehyde in line with experimental evidence. The pyridyl ring of the PLP-serine aldimine substrate exists in deprotonated form, being stabilized directly by Asp208 via a strong H-bond, as well as through interactions with Arg371, Lys237, and His211, and with the surrounding protein which was electrostatically embedded. This knowledge has the potential to impact the design and development of new inhibitors.

QM/MM molecular dynamics study of the galactopyranose → galactofuranose reaction catalysed by Trypanosoma cruzi UDP-galactopyranose mutase

The enzyme UDP-Galactopyranose Mutase (UGM) catalyses the conversion of galactopyranose into galactofuranose. It is known to be critical for the survival and proliferation of several pathogenic agents, both prokaryotic and eukaryotic. Among them is Trypanosoma cruzi, the parasite responsible for Chagas' disease. Since the enzyme is not present in mammals, it appears as a promising target for the design of drugs to treat this illness. A precise knowledge of the mechanism of the catalysed reaction would be crucial to assist in such design. In this article we present a detailed study of all the putative steps of the mechanism. The study is based on QM/MM free energy calculations along properly selected reaction coordinates, and on the analysis of the main structural changes and interactions taking place at every step. The results are discussed in connection with the experimental evidence and previous theoretical studies.

NMR structural study of the prototropic equilibrium in solution of Schiff bases as model compounds

An NMR titration method has been used to simultaneously measure the acid dissociation constant (pKa) and the intramolecular NHO prototropic constant ΔKNHO on a set of Schiff bases. The model compounds were synthesized from benzylamine and substituted ortho-hydroxyaldehydes, appropriately substituted with electron-donating and electron-withdrawing groups to modulate the acidity of the intramolecular NHO hydrogen bond. The structure in solution was established by 1H-, 13C- and 15N-NMR spectroscopy. The physicochemical parameters of the intramolecular NHO hydrogen bond (pKa, ΔKNHO and ΔΔG°) were obtained from 1H-NMR titration data and pH measurements. The Henderson-Hasselbalch data analysis indicated that the systems are weakly acidic, and the predominant NHO equilibrium was established using Polster-Lachmann δ-diagram analysis and Perrin model data linearization.

Free-energy computations identify the mutations required to confer trans-sialidase activity into Trypanosoma rangeli sialidase

Trypanosoma rangeli's sialidase (TrSA) and Trypanosoma cruzi's trans-sialidase (TcTS) are members of the glycoside hydrolase family 33 (GH-33). They share 70% of sequence identity and their crystallographic Cα RMSD is 0.59 Å. Despite these similarities they catalyze different reactions. TcTS transfers sialic acid between glycoconjugates while TrSA can only cleave sialic acid from sialyl-glyconjugates. Significant effort has been invested into unraveling the differences between TrSA and TcTS, and into conferring TrSA with trans-sialidase activity through appropriate point mutations. Recently, we calculated the free-energy change for the formation of the covalent intermediate (CI) in TcTS and performed an energy decomposition analysis of that process. In this article we present a similar study for the formation of the CI in TrSA, as well as in a quintuple mutant (TrSA5mut), which has faint trans-sialidase activity. The comparison of these new results with those previously obtained for TcTS allowed identifying five extra mutations to be introduced in TrSA5mut that should create a mutant (TrSA10mut ) with high trans-sialidase activity.

Aromatic L-amino acid decarboxylase (AADC) is crucial for brain development and motor functions

Aromatic L-amino acid decarboxylase (AADC) deficiency is a rare pediatric neuro-metabolic disease in children. Due to the lack of an animal model, its pathogenetic mechanism is poorly understood. To study the role of AADC in brain development, a zebrafish model of AADC deficiency was generated. We identified an aadc gene homolog, dopa decarboxylase (ddc), in the zebrafish genome. Whole-mount in situ hybridization analysis showed that the ddc gene is expressed in the epiphysis, locus caeruleus, diencephalic catecholaminergic clusters, and raphe nuclei of 36-h post-fertilization (hpf) zebrafish embryos. Inhibition of Ddc by AADC inhibitor NSD-1015 or anti-sense morpholino oligonucleotides (MO) reduced brain volume and body length. We observed increased brain cell apoptosis and loss of dipencephalic catecholaminergic cluster neurons in ddc morphants (ddc MO-injected embryos). Seizure-like activity was also detected in ddc morphants in a dose-dependent manner. ddc morphants had less sensitive touch response and impaired swimming activity that could be rescued by injection of ddc plasmids. In addition, eye movement was also significantly impaired in ddc morphants. Collectively, loss of Ddc appears to result in similar phenotypes as that of ADCC deficiency, thus zebrafish could be a good model for investigating pathogenetic mechanisms of AADC deficiency in children.

Glutamate decarboxylase-dependent acid resistance in orally acquired bacteria: function, distribution and biomedical implications of the gadBC operon

For successful colonization of the mammalian host, orally acquired bacteria must overcome the extreme acidic stress (pH < 2.5) encountered during transit through the host stomach. The glutamate-dependent acid resistance (GDAR) system is by far the most potent acid resistance system in commensal and pathogenic Escherichia coli, Shigella flexneri, Listeria monocytogenes and Lactococcus lactis. GDAR requires the activity of glutamate decarboxylase (GadB), an intracellular PLP-dependent enzyme which performs a proton-consuming decarboxylation reaction, and of the cognate antiporter (GadC), which performs the glutamatein /γ-aminobutyrateout (GABA) electrogenic antiport. Herein we review recent findings on the structural determinants responsible for pH-dependent intracellular activation of E. coli GadB and GadC. A survey of genomes of bacteria (pathogenic and non-pathogenic), having in common the ability to colonize or to transit through the host gut, shows that the gadB and gadC genes frequently lie next or near each other. This gene arrangement is likely to be important to ensure timely co-regulation of the decarboxylase and the antiporter. Besides the involvement in acid resistance, GABA production and release were found to occur at very high levels in lactic acid bacteria originally isolated from traditionally fermented foods, supporting the evidence that GABA-enriched foods possess health-promoting properties.

Triple hydrogen bonding in a circular arrangement: ab initio, DFT and first-principles MD studies of tris-hydroxyaryl enamines

First-principles Car-Parrinello molecular dynamics, ab initio (MP2) and density functional schemes have been used to explore the tautomeric equilibrium in three tris(amino(R)methylidene)cyclohexane-1,3,5-triones (R = hydrogen, methyl or phenyl group). The dynamic nature of the cyclic hydrogen bonding has been studied by the first-principles MD method. The comparison of the results obtained by aforesaid methods has been accomplished on the basis of calculations of structural and spectroscopic characteristics of the compounds. The conformational analysis of the studied compounds has been carried out at the MP2/6-31+G(d,p) and B3LYP/6-31+G(d,p) levels of theory. The influence of steric and electronic effects on the cyclic hydrogen bonding has been analysed. The extent of the proton delocalization has been modified by the substituents according to the sequence: hydrogen < phenyl < methyl. This fact is verified by the spectroscopic and structural data as well as the energy potential curve. A prevalence of the keto-enamine tautomeric form has been observed in the static ab initio and DFT models, and confirmed by the first-principles MD.

Desmotropy, polymorphism, and solid-state proton transfer: four solid forms of an aromatic o-hydroxy Schiff base

The Schiff base derived from salicylaldehyde and 2-amino-3-hydroxypyridine affords a diversity of solid forms, two polymorphic pairs of the enol-imino (D1 a and D1 b) and keto-amino (D2 a and D2 b) desmotropes. The isolated phases, identified by IR spectroscopy, X-ray crystallography, and (13)C cross-polarization/magnetic angle spinning (CP/MAS) NMR spectroscopy, display essentially planar molecular conformations characterized by strong intramolecular hydrogen bonds of the O-H⋅⋅⋅N (D1) or N-H⋅⋅⋅O (D2) type. A change in the position of the proton within this O⋅⋅⋅H⋅⋅⋅N system is accompanied by substantially different molecular conformations and, subsequently, by divergent supramolecular architectures. The appearance and interconversion conditions for each of the four phases have been established on the basis of a number of solution and solvent-free experiments, and evaluated against the results of computational studies. Solid phases readily convert into the most stable form (D1 a) upon exposure to methanol vapor, heating, or by mechanical treatment, and these transformations are accompanied by a change in the color of the sample. The course of thermally induced transformations has been monitored in detail by means of temperature-resolved powder X-ray diffraction and infrared spectroscopy. Upon dissolution, all forms equilibrate immediately, as confirmed by NMR and UV/Vis spectroscopy in several solvents, with the equilibrium shifted far towards the enol tautomer. This study reveals the significance of peripheral groups in the stabilization of metastable tautomers in the solid state.

Computation of kinetic isotope effects for enzymatic reactions

We describe a computational approach, incorporating quantum mechanics into enzyme kinetics modeling with a special emphasis on computation of kinetic isotope effects. Two aspects are highlighted: (1) the potential energy surface is represented by a combined quantum mechanical and molecular mechanical (QM/MM) potential in which the bond forming and breaking processes are modeled by electronic structure theory, and (2) a free energy perturbation method in path integral simulation is used to determine both kinetic isotope effects (KIEs). In this approach, which is called the PI-FEP/UM method, a light (heavy) isotope is mutated into a heavy (light) counterpart in centroid path integral simulations. The method is illustrated in the study of primary and secondary KIEs in two enzyme systems. In the case of nitroalkane oxidase, the enzymatic reaction exhibits enhanced quantum tunneling over that of the uncatalyzed process in water. In the dopa delarboxylase reaction, there appears to be distinguishable primary carbon-13 and secondary deuterium KIEs when the internal proton tautomerism is in the N-protonated or in the O-protonated positions. These examples show that the incorporation of quantum mechanical effects in enzyme kinetics modeling offers an opportunity to accurately and reliably model the mechanisms and free energies of enzymatic reactions.

Critical hydrogen bonds and protonation states of pyridoxal 5'-phosphate revealed by NMR

In this contribution we review recent NMR studies of protonation and hydrogen bond states of pyridoxal 5'-phosphate (PLP) and PLP model Schiff bases in different environments, starting from aqueous solution, the organic solid state to polar organic solution and finally to enzyme environments. We have established hydrogen bond correlations that allow one to estimate hydrogen bond geometries from (15)N chemical shifts. It is shown that protonation of the pyridine ring of PLP in aspartate aminotransferase (AspAT) is achieved by (i) an intermolecular OHN hydrogen bond with an aspartate residue, assisted by the imidazole group of a histidine side chain and (ii) a local polarity as found for related model systems in a polar organic solvent exhibiting a dielectric constant of about 30. Model studies indicate that protonation of the pyridine ring of PLP leads to a dominance of the ketoenamine form, where the intramolecular OHN hydrogen bond of PLP exhibits a zwitterionic state. Thus, the PLP moiety in AspAT carries a net positive charge considered as a pre-requisite to initiate the enzyme reaction. However, it is shown that the ketoenamine form dominates in the absence of ring protonation when PLP is solvated by polar groups such as water. Finally, the differences between acid-base interactions in aqueous solution and in the interior of proteins are discussed. This article is part of a special issue entitled: Pyridoxal Phosphate Enzymology.

Molecular dynamics simulations of the intramolecular proton transfer and carbanion stabilization in the pyridoxal 5'-phosphate dependent enzymes L-dopa decarboxylase and alanine racemase

Molecular dynamics simulations using a combined quantum mechanical and molecular mechanical (QM/MM) potential have been carried out to investigate the internal proton transfer equilibrium of the external aldimine species in l-dopa decarboxylase, and carbanion stabilization by the enzyme cofactor in the active site of alanine racemase. Solvent effects lower the free energy of the O-protonated PLP tautomer both in aqueous solution and in the active site, resulting a free energy difference of about -1 kcal/mol relative to the N-protonated Schiff base in the enzyme. The external aldimine provides the dominant contribution to lowering the free energy barrier for the spontaneous decarboxylation of l-dopa in water, by a remarkable 16 kcal/mol, while the enzyme l-dopa decarboxylase further lowers the barrier by 8 kcal/mol. Kinetic isotope effects were also determined using a path integral free energy perturbation theory on the primary (13)C and the secondary (2)H substitutions. In the case of alanine racemase, if the pyridine ring is unprotonated as that in the active site, there is destabilizing contribution to the formation of the α-carbanion in the gas phase, although when the pyridine ring is protonated the contribution is stabilizing. In aqueous solution and in alanine racemase, the α-carbanion is stabilized both when the pyridine ring is protonated and unprotonated. The computational studies illustrated in this article show that combined QM/MM simulations can help provide a deeper understanding of the mechanisms of PLP-dependent enzymes. This article is part of a Special Issue entitled: Pyridoxal Phosphate Enzymology.

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.

Substituent effects on electrophilic catalysis by the carbonyl group: anatomy of the rate acceleration for PLP-catalyzed deprotonation of glycine

First-order rate constants, determined by (1)H NMR, are reported for deuterium exchange between solvent D(2)O and the α-amino carbon of glycine in the presence of increasing concentrations of carbonyl compounds (acetone, benzaldehyde, and salicylaldehyde) and at different pD and buffer concentrations. These rate data were combined with (1)H NMR data that define the position of the equilibrium for formation of imines/iminium ions from addition of glycine to the respective carbonyl compounds, to give second-order rate constants k(DO) for deprotonation of α-imino carbon by DO(-). The assumption that these second-order rate constants lie on linear structure-reactivity correlations between log k(OL) and pK(a) was made in estimating the following pK(a)'s for deprotonation of α-imino carbon: pK(a) = 22, glycine-acetone iminium ion; pK(a) = 27, glycine-benzaldehyde imine; pK(a) ≈ 23, glycine-benzaldehyde iminium ion; and, pK(a) = 25, glycine-salicylaldehyde iminium ion. The much lower pK(a) of 17 [Toth, K.; Richard, J. P. J. Am. Chem. Soc. 2007, 129, 3013-3021] for carbon deprotonation of the adduct between 5'-deoxypyridoxal (DPL) and glycine shows that the strongly electron-withdrawing pyridinium ion is unique in driving the extended delocalization of negative charge from the α-iminium to the α-pyridinium carbon. This favors carbanion protonation at the α-pyridinium carbon, and catalysis of the 1,3-aza-allylic isomerization reaction that is a step in enzyme-catalyzed transamination reactions. An analysis of the effect of incremental changes in structure on the activity of benzaldehyde in catalysis of deprotonation of glycine shows the carbonyl group electrophile, the 2-O(-) ring substituent and the cation pyridinium nitrogen of DPL each make a significant contribution to the catalytic activity of this cofactor analogue. The extraordinary activity of DPL in catalysis of deprotonation of α-amino carbon results from the summation of these three smaller effects.

Proton transfer facilitated by ligand binding. An energetic analysis of the catalytic mechanism of Trypanosoma cruzi trans-sialidase

Trans-sialidase is a crucial enzyme for the infection of Trypanosoma cruzi, the protozoa responsible for Chagas' disease in humans. This enzyme catalyzes the transfer of sialic acids from mammalian host cells to parasitic cell surfaces in order to mask the infection from the host's immune system. It represents a promising target for the development of therapeutics to treat the disease and has been subject of extensive structural studies. Elaborate experiments suggested formation of a long-lived covalent intermediate in the catalytic mechanism and identified a Tyr/Glu pair as an unusual catalytic couple. This requires that the tyrosine hydroxyl proton is transferred to the carboxylate group of glutamate before the nucleophilic attack. Since the solution pK(a)s of tyrosine and glutamate are very different, this transfer can only be accomplished if the reaction environment selectively stabilizes the product state. We compute the free energy profile for the proton transfer in different environments, and our results indicate that it can take place in the active site of trans-sialidase, but only after substrate binding. By means of the energy decomposition method, we explain the influence that the active site residues exert on the reaction and how the pattern is changed when the substrate is present. This study represents an initial step that can shed light on our understanding of the catalytic mechanism of this reaction.