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

Atomic-level structure of the amorphous drug atuliflapon via NMR crystallography

We determine the complete atomic-level structure of the amorphous form of the drug atuliflapon, a 5-lipooxygenase activating protein (FLAP) inhibitor, via chemical-shift-driven NMR crystallography. The ensemble of preferred structures allows us to identify a number of specific conformations and interactions that stabilize the amorphous structure. These include preferred hydrogen-bonding motifs with water and with other drug molecules, as well as conformations of the cyclohexane and pyrazole rings that stabilize structure by indirectly allowing for optimization of hydrogen bonding.

Charge Relay Without Proton Transfer: Coupling of Two Short Hydrogen Bonds via Imidazole in Models of Catalytic Triad of Serine Protease Active Site

A homologous series of 20 substituted alcohol-imidazole-acetate model complexes imitating the charge relay system in Ser-His-Asp catalytic triad of serine proteases is considered quantum-chemically. We show qualitatively that the geometries of alcohol-imidazole and imidazole-acetate short hydrogen bonds are strongly coupled via the central imidazole and such complexes are capable of effectively relaying the charge from acetate to alcohol moiety upon relatively small concerted proton displacements. We hypothesize an alternative catalytic mechanism of serine proteases that does not require two complete proton transfers or hydrogen bond breakage between Ser and His residues.

Complexes of Gold(III) with Hydrazones Derived from Pyridoxal: Stability, Structure, and Nature of UV-Vis Spectra

Pyridoxal and pyridoxal 5'-phosphate are aldehyde forms of B6 vitamin that can easily be transformed into each other in the living organism. The presence of a phosphate group, however, provides the related compounds (e.g., hydrazones) with better solubility in water. In addition, the phosphate group may sometimes act as a binding center for metal ions. In particular, a phosphate group can be a strong ligand for a gold(III) ion, which is of interest for researchers for the anti-tumor and antimicrobial potential of gold(III). This paper aims to answer whether the phosphate group is involved in the complex formation between gold(III) and hydrazones derived from pyridoxal 5'-phosphate. The answer is negative, since the comparison of the stability constants determined for the gold(III) complexes with pyridoxal- and pyridoxal 5'-phosphate-derived hydrazones showed a negligible difference. In addition, quantum chemical calculations confirmed that the preferential coordination of two series of phosphorylated and non-phosphorylated hydrazones to gold(III) ion is similar. The preferential protonation modes for the gold(III) complexes were also determined using experimental and calculated data.

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.

Mechanism-based and computational modeling of hydrogen sulfide biogenesis inhibition: interfacial inhibition

Hydrogen sulfide (H₂S) is a gaseous signaling molecule that participates in various signaling functions in health and diseases. The tetrameric cystathionine γ-lyase (CSE) contributes to H₂S biogenesis and several investigations provide evidence on the pharmacological modulation of CSE as a potential target for the treatment of a multitude of conditions. D-penicillamine (D-pen) has recently been reported to selectively impede CSE-catalyzed H₂S production but the molecular bases for such inhibitory effect have not been investigated. In this study, we report that D-pen follows a mixed-inhibition mechanism to inhibit both cystathionine (CST) cleavage and H₂S biogenesis by human CSE. To decipher the molecular mechanisms underlying such a mixed inhibition, we performed docking and molecular dynamics (MD) simulations. Interestingly, MD analysis of CST binding reveals a likely active site configuration prior to gem-diamine intermediate formation, particularly H-bond formation between the amino group of the substrate and the O3' of PLP. Similar analyses realized with both CST and D-pen identified three potent interfacial ligand-binding sites for D-pen and offered a rational for D-pen effect. Thus, inhibitor binding not only induces the creation of an entirely new interacting network at the vicinity of the interface between enzyme subunits, but it also exerts long range effects by propagating to the active site. Overall, our study paves the way for the design of new allosteric interfacial inhibitory compounds that will specifically modulate H₂S biogenesis by cystathionine γ-lyase.

Vitamin B6-Based Biomimetic Asymmetric Catalysis

ConspectusOne of the fundamental goals of chemists is to develop highly efficient methods for producing optically active compounds, given their wide range of applications in chemistry, pharmaceutical industry, chemical biology, and material science. Biomimetic asymmetric catalysis, which imitates the structures and functions of enzymes, has emerged as an extremely attractive strategy for producing chiral compounds. This field has drawn tremendous research interest and has led to various protocols for constructing complex molecular scaffolds. The Vitamin B₆ family, including pyridoxal, pyridoxamine, pyridoxine, and the corresponding phosphorylated derivatives, serves as the cofactors to catalyze more than 200 enzymatic functions, accounting for ∼4% of all enzyme activities. Although significant progress has been made in simulating the biological roles of vitamin B₆ during the past several decades, its extraordinary catalytic power has not yet been successfully applied into asymmetric synthesis. In recent years, our group has been devoted to developing vitamin B₆-based biomimetic asymmetric catalysis using chiral pyridoxals/pyridoxamines as catalysts. We are particularly interested in mimicking the processes of enzymatic transamination and biological aldol reaction of glycine, respectively, developing asymmetric biomimetic transamination and carbonyl catalysis enabled α-C-H transformation of primary amines. Using a chiral α,α-diarylprolinol-derived pyridoxal as the catalyst, we reported the first chiral pyridoxal catalyzed asymmetric transamination of α-keto acids in 2015. A significant breakthrough in biomimetic transamination was achieved by using an axially chiral biaryl pyridoxamine catalyst that bears a lateral amine side arm. The amine side arm acts as an intramolecular base, accelerating the transamination and proving highly effective for transamination of α-keto acids and α-keto amides. In addition, we discovered the catalytic power of chiral pyridoxals as carbonyl catalysts for asymmetric biomimetic Mannich/aldol reactions of glycinates. These chiral pyridoxals also enabled more α-C-H conversions of glycinates, such as asymmetric 1,4-addition toward α,β-unsaturated esters and asymmetric α-allylation with Morita-Baylis-Hillman acetates. Moreover, carbonyl catalysis can be further applied to highly challenging primary amines with inert α-C-H bonds, such as propargylamines and benzylamines, which represents a powerful strategy for direct asymmetric α-C-H functionalization of various primary amines without protection of the NH₂ group. These biomimetic/bioinspired transformations provide efficient new protocols for the synthesis of chiral amines. Herein, we summarize our recent efforts on the development of the vitamin B₆-based biomimetic asymmetric catalysis.

Enzyme-independent catabolism of cysteine with pyridoxal-5'-phosphate

Pyridoxal-5'-phosphate (PLP) is a versatile cofactor that assists in different types of enzymatic reactions. PLP has also been reported to react with substrates and catalyze some of these reactions independent of enzymes. One such catalytic reaction is the breakdown of cysteine to produce hydrogen sulfide (H₂S) in the presence of multivalent metal ions. However, the enzyme-independent catalytic activity of PLP in catabolizing cysteine in the absence of multivalent ions is unknown. In this study, we show that PLP reacts with cysteine to form a thiazolidine product, which is supported by quantum chemical calculations of the absorption spectrum. The reaction of PLP with cysteine is dependent on ionic strength and pH. The thiazolidine product slowly decomposes to produce H₂S and the PLP regenerates to its active form with longer reaction times (> 24 h), suggesting that PLP can act as a catalyst. We propose an enzyme-independent plausible reaction mechanism for PLP catalyzed cysteine breakdown to produce H₂S, which proceeds through the formation of thiazolidine ring intermediates that later hydrolyzes slowly to regenerate the PLP. This work demonstrates that PLP catalyzes cysteine breakdown in the absence of enzymes, base, and multivalent metal ions to produce H₂S.

Shared Hydrogen Atom Location and Chemical Composition in Picolinic Acid and Pyridoxal Hydrochloride Derivatives Determined by X-ray Crystallography

Three new single-crystal structures were isolated for picolinic acid (2), the trifluoroacetate salt of picolinic acid (1), and pyridoxal hydrochloride (3). These compounds displayed unconventional crystallographic features that must be considered when structural refinements are carried out. Thus, the generated Fourier differences map obtained with the diffraction data collected at 100 K was crucial to visualize electron densities, which were balanced by either one hydrogen atom or a hydrogen atom with an occupancy factor of 1/2 located between either two carboxylate moieties, two phenolic oxygen atoms, or two pyridinic nitrogen atoms. Moreover, NMR studies were conducted to analyze the bulk chemical composition of single crystals of 2-pyridinecarboxylic acid obtained from the gem-diol/hemiacetal forms and the polymerization products after the treatment of 2-pyridinecarboxaldehyde with TFA:H₂O (1) or a diluted Cu(NO₃)₂ solution (2). The quantitative yield of the pyridoxal hydrochloride crystalline material (3) obtained from a diluted CuCl₂ solution was exhaustively characterized by solid-state NMR methods. These methods allowed the resolution of the signals corresponding to the protons of the hydroxyl moiety of the intramolecular hemiacetal group and the phenolic hydrogen. Theoretical calculations using DFT methods were done to complement the atomic location of the hydrogen atoms obtained from the X-ray analysis.

An N⋯H⋯N low-barrier hydrogen bond preorganizes the catalytic site of aspartate aminotransferase to facilitate the second half-reaction

Pyridoxal 5'-phosphate (PLP)-dependent enzymes have been extensively studied for their ability to fine-tune PLP cofactor electronics to promote a wide array of chemistries. Neutron crystallography offers a straightforward approach to studying the electronic states of PLP and the electrostatics of enzyme active sites, responsible for the reaction specificities, by enabling direct visualization of hydrogen atom positions. Here we report a room-temperature joint X-ray/neutron structure of aspartate aminotransferase (AAT) with pyridoxamine 5'-phosphate (PMP), the cofactor product of the first half reaction catalyzed by the enzyme. Between PMP NSB and catalytic Lys258 Nζ amino groups an equally shared deuterium is observed in an apparent low-barrier hydrogen bond (LBHB). Density functional theory calculations were performed to provide further evidence of this LBHB interaction. The structural arrangement and the juxtaposition of PMP and Lys258, facilitated by the LBHB, suggests active site preorganization for the incoming ketoacid substrate that initiates the second half-reaction.

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.

Free Energy Landscape and Proton Transfer Pathways of the Transimination Reaction at the Active site of the Serine Hydroxymethyltransferase Enzyme in Aqueous Medium

Serine hydroxymethyltransferase (SHMT) is a ubiquitous enzyme belonging to the fold type I or aspartate aminotransferase (AspAT) family of the pyridoxal 5'-phosphate (PLP)-dependent enzymes. Like other PLP-dependent enzymes, SHMT also undergoes the so-called transimination reaction before exhibiting its enzymatic activity. The transimination process constitutes an important pre-step for all PLP-dependent enzymes, where an internal aldimine of the PLP-enzyme complex gets converted to an external aldimine of the substrate-PLP complex at the active site of the enzyme. In case of the transimination reaction involving SHMT, the PLP molecule bound to the active site lysine residue of SHMT (internal aldimine) gets detached from the enzyme by a serine substrate to produce an external aldimine complex, where the PLP is now bound to the serine substrate. In the current study, the free energy surfaces and reaction pathways of different steps of the transimination reaction at the active site of SHMT are investigated by employing hybrid quantum mechanical/molecular mechanical (QM/MM) simulations combined with metadynamics methods of rare event sampling. It is found that the process of transimination involving serine and PLP at the active site of the SHMT enzyme takes place through different elementary steps such as the formation of the first geminal diamine intermediate (GDI1), transfer of a proton from the substrate serine to the phenolic oxygen of PLP, followed by another proton transfer from PLP to the amine nitrogen of lysine with the formation of the second geminal diamine intermediate (GDI2), and finally, detachment of the active site lysine residue from PLP to produce the external aldimine.

Asymmetric biomimetic transamination of α-keto amides to peptides

Peptides are important compounds with broad applications in many areas. Asymmetric transamination of α-keto amides can provide an efficient strategy to synthesize peptides, however, the process has not been well developed yet and still remains a great challenge in both enzymatic and catalytic chemistry. For biological transamination, the high activity is attributed to manifold structural and electronic factors of transaminases. Based on the concept of multiple imitation of transaminases, here we report N-quaternized axially chiral pyridoxamines 1 for enantioselective transamination of α-keto amides, to produce various peptides in good yields with excellent enantio- and diastereoselectivities. The reaction is especially attractive for the synthesis of peptides made of unnatural amino acids since it doesn’t need great efforts to make chiral unnatural amino acids before amide bond formation.

Asymmetric transamination of α-keto amides could provide an efficient strategy to synthesise peptides, but has not been well developed yet. Here, the authors design chiral pyridoxamine catalyst and realize the asymmetric biomimetic transamination of α-keto amides, providing access to various peptides with excellent enantiopurities.

Actual Symmetry of Symmetric Molecular Adducts in the Gas Phase, Solution and in the Solid State

This review discusses molecular adducts, whose composition allows a symmetric structure. Such adducts are popular model systems, as they are useful for analyzing the effect of structure on the property selected for study since they allow one to reduce the number of parameters. The main objectives of this discussion are to evaluate the influence of the surroundings on the symmetry of these adducts, steric hindrances within the adducts, competition between different noncovalent interactions responsible for stabilizing the adducts, and experimental methods that can be used to study the symmetry at different time scales. This review considers the following central binding units: hydrogen (proton), halogen (anion), metal (cation), water (hydrogen peroxide).

1,3,5-Triaza-7-Phosphaadamantane (PTA) as a 31P NMR Probe for Organometallic Transition Metal Complexes in Solution

Due to the rigid structure of 1,3,5-triaza-7-phosphaadamantane (PTA), its ³¹P chemical shift solely depends on non-covalent interactions in which the molecule is involved. The maximum range of change caused by the most common of these, hydrogen bonding, is only 6 ppm, because the active site is one of the PTA nitrogen atoms. In contrast, when the PTA phosphorus atom is coordinated to a metal, the range of change exceeds 100 ppm. This feature can be used to support or reject specific structural models of organometallic transition metal complexes in solution by comparing the experimental and Density Functional Theory (DFT) calculated values of this ³¹P chemical shift. This approach has been tested on a variety of the metals of groups 8-12 and molecular structures. General recommendations for appropriate basis sets are reported.

Modeling of Solute-Solvent Interactions Using an External Electric Field-From Tautomeric Equilibrium in Nonpolar Solvents to the Dissociation of Alkali Metal Halides

An implicit account of the solvent effect can be carried out using traditional static quantum chemistry calculations by applying an external electric field to the studied molecular system. This approach allows one to distinguish between the effects of the macroscopic reaction field of the solvent and specific solute-solvent interactions. In this study, we report on the dependence of the simulation results on the use of the polarizable continuum approximation and on the importance of the solvent effect in nonpolar solvents. The latter was demonstrated using experimental data on tautomeric equilibria between the pyridone and hydroxypyridine forms of 2,6-di-tert-butyl-4-hydroxy-pyridine in cyclohexane and chloroform.

Accurate location of hydrogen atoms in hydrogen bonds of tizoxanide from the combination of experimental and theoretical models

To obtain detailed information about the position of hydrogen atoms in hydrogen bonds, HBs, of crystalline organic molecular compounds is not an easy task. In this work we propose a combination of ssNMR experimental data with theoretical procedures to get such information. Furthermore, the combination of experimental and theoretical models provides us with well-defined grounds to analyse the strength of π-stacking interactions between layers of hydrogen bonded molecules. Two different theoretical models were considered, both approaches being quite different. The first one is a solid-state model, so that the periodicity of a crystalline system underlies calculations of the electronic energy, the electronic density and NMR parameters. The other one is a molecular model in which molecules are taken as isolated monomers, dimers and tetramers. These two models were applied to the tizoxanide, TIZ, molecular crystal though it can widely be applied to any other molecular crystal. By the application of the quantum molecular model it was possible to learn about the way the intermolecular HBs affect the position of hydrogen atoms that belong to HBs in TIZ. This molecule has two intermolecular HBs that stabilize the structure of a basic dimer, but it also has an intramolecular HB in each monomer whose position should be optimized together with the other ones. We found that by doing this it is possible to obtain reliable results of calculations of NMR spectroscopic parameters. Working with the solid-state model we found that any local variation of the TIZ crystalline structure is correlated with the variation of the values of the NMR parameters of each nucleus. The excellent agreement between experimental and calculated chemical shifts leads to the conclusion that the N10-H10 bond distance should be (1.00 ± 0.02) Å.

Crystallographic Snapshots of the Dunathan and Quinonoid Intermediates provide Insights into the Reaction Mechanism of Group II Decarboxylases

PLP-dependent enzymes catalyze a plethora of chemical reactions affecting diverse physiological functions. Here we report the structural determinants of the reaction mechanism in a Group II PLP-dependent decarboxylase by assigning two early intermediates. The in-crystallo complexes of the PLP bound form, and the Dunathan and quinonoid intermediates, allowed direct observation of the active site interactions. The structures reveal that a subtle rearrangement of a conserved Arg residue in concert with a water-mediated interaction with the carboxylate of the Dunathan intermediate, appears to directly stabilize the alignment and facilitate the release of CO₂ to yield the quinonoid. Modeling indicates that the conformational change of a dynamic catalytic loop to a closed form controls a conserved network of hydrogen bond interactions between catalytic residues to protonate the quinonoid. Our results provide a structural framework to elucidate mechanistic roles of residues that govern reaction specificity and catalysis in PLP-dependent decarboxylation.

Catalytic Roles of Coenzyme Pyridoxal-5'-phosphate (PLP) in PLP-dependent Enzymes: Reaction Pathway for Methionine-γ-lyase-catalyzed L-methionine Depletion

Pyridoxal-5'-phosphate (PLP), the active form of vitamin B₆, is an important and versatile coenzyme involved in a variety of enzymatic reactions, accounting for about 4% of all classified activities. However, the detailed catalytic reaction pathways for PLP-dependent enzymes remain to be explored. Methionine-γ-lyase (MGL), a promising alternative anti-tumor agent to conventional chemotherapies whose catalytic mechanism is highly desired for guiding further development of re-engineered enzymes, was used as a representative PLP-dependent enzyme, and the catalytic mechanism for L-Met elimination by MGL was explored at the first-principles quantum mechanical/molecular mechanical (QM/MM) level with umbrella sampling. The QM/MM calculations revealed that the enzymatic reaction pathway consists of 4 stages for a total of 19 reaction steps with five intermediates captured in available crystal structures. Furthermore, the more comprehensive role of PLP was revealed. Besides the commonly known role of "electron sink", coenzyme PLP can also assist proton transfer and temporarily store the excess proton generated in some intermediate states by using its hydroxyl group and phosphate group. Thus, PLP is participated in most of the 19 steps. This study not only provided a theoretical basis for further development and re-engineering MGL as a potential anti-tumor agent, but also revealed the comprehensive role of PLP which could be used to explore the mechanisms of other PLP-dependent enzymes.

Synthesis of 15 N-labelled 3,5-dimethylpyridine

¹⁵N‐labelled pyridines are liquid‐ and solid‐state nuclear magnetic resonance (NMR) probes for chemical and biological environments because their ¹⁵N chemical shifts are sensitive to hydrogen‐bond and protonation states. By variation of the type and number of substituents, different target pyridines can be synthesized exhibiting different pK_a values and molecular volumes. Various synthetic routes have been described in the literature, starting from different precursors or modification of other ¹⁵N‐labelled pyridines. In this work, we have explored the synthesis of ¹⁵N ¹⁵N‐labelled pyridines using a two‐step process via the synthesis of alkoxy‐3,4‐dihydro‐2H‐pyran as precursor exhibiting already the desired pyridine substitution pattern. As an example, we have synthesized 3,5‐dimethylpyridine‐¹⁵N (lutidine‐¹⁵N) as demonstrated by ¹⁵N‐NMR spectroscopy. That synthesis starts from methacrolein, propenyl ether, and ¹⁵N‐labelled NH₄Cl as nitrogen source.

Structural insights into the mechanism of internal aldimine formation and catalytic loop dynamics in an archaeal Group II decarboxylase

Formation of the internal aldimine (LLP) is the first regulatory step that activates pyridoxal 5'-phosphate (PLP) dependent enzymes. The process involves a nucleophilic attack on PLP by an active site Lys residue, followed by proton transfers resulting in a carbinolamine (CBA) intermediate that undergoes dehydration to form the aldimine. Despite a general understanding of the pathway, the structural basis of the mechanistic roles of specific residues in each of these steps is unclear. Here we determined the crystal structure of the LLP form (holo-form) of a Group II PLP-dependent decarboxylase from Methanocaldococcus jannaschii (MjDC) at 1.7 Å resolution. By comparing the crystal structure of MjDC in the LLP form with that of the pyridoxal-P (non-covalently bound aldehyde) form, we demonstrate structural evidence for a water-mediated mechanism of LLP formation. A conserved extended hydrogen-bonding network around PLP coupled to the pyridinyl nitrogen influences activation and catalysis by affecting the electronic configuration of PLP. Furthermore, the two cofactor bound forms revealed open and closed conformations of the catalytic loop (CL) in the absence of a ligand, supporting a hypothesis for a regulatory link between LLP formation and CL dynamics. The evidence suggests that activation of Group II decarboxylases involves a complex interplay of interactions between the electronic states of PLP, the active site micro-environment and CL dynamics.

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.

Enamine/Dienamine and Brønsted Acid Catalysis: Elusive Intermediates, Reaction Mechanisms, and Stereoinduction Modes Based on in Situ NMR Spectroscopy and Computational Studies

Over the years, the field of enantioselective organocatalysis has seen unparalleled growth in the development of novel synthetic applications with respect to mechanistic investigations. Reaction optimization appeared to be rather empirical than rational. This offset between synthetic development and mechanistic understanding was and is generally due to the difficulties in detecting reactive intermediates and the inability to experimentally evaluate transition states. Thus, the first key point for mechanistic studies is detecting elusive intermediates and characterizing them in terms of their structure, stability, formation pathways, and kinetic properties. The second key point is evaluating the importance of these intermediates and their properties in the transition state.

In the past 7 years, our group has addressed the problems with detecting elusive intermediates in organocatalysis by means of NMR spectroscopy and eventually theoretical calculations. Two main activation modes were extensively investigated: secondary amine catalysis and, very recently, Brønsted acid catalysis. Using these examples, we discuss potential methods to stabilize intermediates via intermolecular interactions; to elucidate their structures, formation pathways and kinetics; to change the kinetics of the reactions; and to address their relevance in transition states. The elusive enamine in proline-catalyzed aldol reactions is used as an example of the stabilization of intermediates via inter- and intramolecular interactions; the determination of kinetics on its formation pathway is discussed. Classical structural characterization of intermediates is described using prolinol and prolinol ether enamines and dienamines. The Z/E dilemma for the second double bond of the dienamines shows how the kinetics of a reaction can be changed to allow for the detection of reaction intermediates. We recently started to investigate substrate–catalyst complexes in the field of Brønsted acid catalysis. These studies on imine/chiral phosphoric acid complexes show that an appropriate combination of highly developed NMR and theoretical methods can provide detailed insights into the complicated structures, exchange kinetics, and H-bonding properties of chiral ion pairs. Furthermore, the merging of these structural investigations and photoisomerization even allowed the active transition state combinations to be determined for the first time on the basis of experimental data only, which is the gold standard in mechanistic investigations and was previously thought to be exclusively the domain of theoretical calculations.

Thus, this Account summarizes our recent mechanistic work in the field of organocatalysis and explains the potential methods for addressing the central questions in mechanistic studies: stabilization of intermediates, elucidation of structures and formation pathways, and addressing transition state combinations experimentally.

NMR Study of Solvation Effect on the Geometry of Proton-Bound Homodimers of Increasing Size

Hydrogen bond geometries in the proton-bound homodimers of quinoline and acridine derivatives in an aprotic polar solution have been experimentally studied using ¹H NMR at 120 K. The reported results show that an increase of the dielectric permittivity of the medium results in contraction of the N···N distance. The degree of contraction depends on the homodimer's size and its substituent-specific solvation features. Neither of these effects can be reproduced using conventional implicit solvent models employed in computational studies. In general, the N···N distance in the homodimers of pyridine, quinoline, and acridine derivatives decreases in the sequence gas phase > solid state > polar solvent.

Direct visualization of critical hydrogen atoms in a pyridoxal 5'-phosphate enzyme

Enzymes dependent on pyridoxal 5′-phosphate (PLP, the active form of vitamin B₆) perform a myriad of diverse chemical transformations. They promote various reactions by modulating the electronic states of PLP through weak interactions in the active site. Neutron crystallography has the unique ability of visualizing the nuclear positions of hydrogen atoms in macromolecules. Here we present a room-temperature neutron structure of a homodimeric PLP-dependent enzyme, aspartate aminotransferase, which was reacted in situ with α-methylaspartate. In one monomer, the PLP remained as an internal aldimine with a deprotonated Schiff base. In the second monomer, the external aldimine formed with the substrate analog. We observe a deuterium equidistant between the Schiff base and the C-terminal carboxylate of the substrate, a position indicative of a low-barrier hydrogen bond. Quantum chemical calculations and a low-pH room-temperature X-ray structure provide insight into the physical phenomena that control the electronic modulation in aspartate aminotransferase.

Pyridoxal 5’-phosphate (PLP) is a ubiquitous co factor for diverse enzymes, among them aspartate aminotransferase. Here the authors use neutron crystallography, which allows the visualization of the positions of hydrogen atoms, and computation to characterize the catalytic mechanism of the enzyme.

Enormous Hydrogen Bond Strength Enhancement through π-Conjugation Gain: Implications for Enzyme Catalysis

Surprisingly large resonance-assistance effects may explain how some enzymes form extremely short, strong hydrogen bonds to stabilize reactive oxyanion intermediates and facilitate catalysis. Computational models for several enzymic residue-substrate interactions reveal that when a π-conjugated, hydrogen bond donor (XH) forms a hydrogen bond to a charged substrate (Y^-), XH can become significantly more π-electron delocalized, and this "extra" stabilization may boost the [XH···Y^-] hydrogen bond strength by ≥15 kcal/mol. This reciprocal relationship departs from the widespread pK_a concept (i.e., the idea that short, strong hydrogen bonds form when the interacting moieties have matching pK_a values), which has been the rationale for enzymic acid-base reactions. The findings presented here provide new insight into how short, strong hydrogen bonds could form in enzymes.

Direct evidence that an extended hydrogen-bonding network influences activation of pyridoxal 5'-phosphate in aspartate aminotransferase

Pyridoxal 5'-phosphate (PLP) is a fundamental, multifunctional enzyme cofactor used to catalyze a wide variety of chemical reactions involved in amino acid metabolism. PLP-dependent enzymes optimize specific chemical reactions by modulating the electronic states of PLP through distinct active site environments. In aspartate aminotransferase (AAT), an extended hydrogen bond network is coupled to the pyridinyl nitrogen of the PLP, influencing the electrophilicity of the cofactor. This network, which involves residues Asp-222, His-143, Thr-139, His-189, and structural waters, is located at the edge of PLP opposite the reactive Schiff base. We demonstrate that this hydrogen bond network directly influences the protonation state of the pyridine nitrogen of PLP, which affects the rates of catalysis. We analyzed perturbations caused by single- and double-mutant variants using steady-state kinetics, high resolution X-ray crystallography, and quantum chemical calculations. Protonation of the pyridinyl nitrogen to form a pyridinium cation induces electronic delocalization in the PLP, which correlates with the enhancement in catalytic rate in AAT. Thus, PLP activation is controlled by the proximity of the pyridinyl nitrogen to the hydrogen bond microenvironment. Quantum chemical calculations indicate that Asp-222, which is directly coupled to the pyridinyl nitrogen, increases the pK_a of the pyridine nitrogen and stabilizes the pyridinium cation. His-143 and His-189 also increase the pK_a of the pyridine nitrogen but, more significantly, influence the position of the proton that resides between Asp-222 and the pyridinyl nitrogen. These findings indicate that the second shell residues directly enhance the rate of catalysis in AAT.

Some Brief Notes on Theoretical and Experimental Investigations of Intramolecular Hydrogen Bonding

A review of selected literature data related to intramolecular hydrogen bonding in ortho-hydroxyaryl Schiff bases, ortho-hydroxyaryl ketones, ortho-hydroxyaryl amides, proton sponges and ortho-hydroxyaryl Mannich bases is presented. The paper reports on the application of experimental spectroscopic measurements (IR and NMR) and quantum-mechanical calculations for investigations of the proton transfer processes, the potential energy curves, tautomeric equilibrium, aromaticity etc. Finally, the equilibrium between the intra- and inter-molecular hydrogen bonds in amides is discussed.

NMR Crystallography of a Carbanionic Intermediate in Tryptophan Synthase: Chemical Structure, Tautomerization, and Reaction Specificity

Carbanionic intermediates play a central role in the catalytic transformations of amino acids performed by pyridoxal-5'-phosphate (PLP)-dependent enzymes. Here, we make use of NMR crystallography-the synergistic combination of solid-state nuclear magnetic resonance, X-ray crystallography, and computational chemistry-to interrogate a carbanionic/quinonoid intermediate analogue in the β-subunit active site of the PLP-requiring enzyme tryptophan synthase. The solid-state NMR chemical shifts of the PLP pyridine ring nitrogen and additional sites, coupled with first-principles computational models, allow a detailed model of protonation states for ionizable groups on the cofactor, substrates, and nearby catalytic residues to be established. Most significantly, we find that a deprotonated pyridine nitrogen on PLP precludes formation of a true quinonoid species and that there is an equilibrium between the phenolic and protonated Schiff base tautomeric forms of this intermediate. Natural bond orbital analysis indicates that the latter builds up negative charge at the substrate Cα and positive charge at C4' of the cofactor, consistent with its role as the catalytic tautomer. These findings support the hypothesis that the specificity for β-elimination/replacement versus transamination is dictated in part by the protonation states of ionizable groups on PLP and the reacting substrates and underscore the essential role that NMR crystallography can play in characterizing both chemical structure and dynamics within functioning enzyme active sites.

Molecular dynamics simulations of apo, holo, and inactivator bound GABA-at reveal the role of active site residues in PLP dependent enzymes

The pyridoxal 5-phosphate (PLP) cofactor is a significant organic molecule in medicinal chemistry. It is often found covalently bound to lysine residues in proteins to form PLP dependent enzymes. An example of this family of PLP dependent enzymes is γ-aminobutyric acid aminotransferase (GABA-AT) which is responsible for the degradation of the neurotransmitter GABA. Its inhibition or inactivation can be used to prevent the reduction of GABA concentration in brain which is the source of several neurological disorders. As a test case for PLP dependent enzymes, we have performed molecular dynamics simulations of GABA-AT to reveal the roles of the protein residues and its cofactor. Three different states have been considered: the apoenzyme, the holoenzyme, and the inactive state obtained after the suicide inhibition by vigabatrin. Different protonation states have also been considered for PLP and two key active site residues: Asp298 and His190. Together, 24 independent molecular dynamics trajectories have been simulated for a cumulative total of 2.88 µs. Our results indicate that, unlike in aqueous solution, the PLP pyridine moiety is protonated in GABA-AT. This is a consequence of a pKa shift triggered by a strong charge-charge interaction with an ionic "diad" formed by Asp298 and His190 that would help the activation of the first half-reaction of the catalytic mechanism in GABA-AT: the conversion of PLP to free pyridoxamine phosphate (PMP). In addition, our MD simulations exhibit additional strong hydrogen bond networks between the protein and PLP: the phosphate group is held in place by the donation of at least three hydrogen bonds while the carbonyl oxygen of the pyridine ring interacts with Gln301; Phe181 forms a π-π stacking interaction with the pyridine ring and works as a gate keeper with the assistance of Val300. All these interactions are hypothesized to help maintain free PMP in place inside the protein active site to facilitate the second half-reaction in GABA-AT: the regeneration of PLP-bound GABA-AT (i.e., the holoenzyme). Proteins 2016; 84:875-891. © 2016 Wiley Periodicals, Inc.

Catalytic roles of βLys87 in tryptophan synthase: (15)N solid state NMR studies

The proposed mechanism for tryptophan synthase shows βLys87 playing multiple catalytic roles: it bonds to the PLP cofactor, activates C4' for nucleophilic attack via a protonated Schiff base nitrogen, and abstracts and returns protons to PLP-bound substrates (i.e. acid-base catalysis). ε-¹⁵N-lysine TS was prepared to access the protonation state of βLys87 using ¹⁵N solid-state nuclear magnetic resonance (SSNMR) spectroscopy for three quasi-stable intermediates along the reaction pathway. These experiments establish that the protonation state of the ε-amino group switches between protonated and neutral states as the β-site undergoes conversion from one intermediate to the next during catalysis, corresponding to mechanistic steps where this lysine residue has been anticipated to play alternating acid and base catalytic roles that help steer reaction specificity in tryptophan synthase catalysis. This article is part of a Special Issue entitled: Cofactor-dependent proteins: evolution, chemical diversity and bio-applications. Guest Editors: Andrea Mozzarelli and Loredano Pollegioni.

Long-range proton relay shows an inverse linear free energy relationship depending on the pKa of the hydrogen-bonded wire

The long-range proton transfer dependence on the pK_a of hydroxyl molecules in hydrogen (H)-bonded wires was investigated using quantum mechanical calculations.

NMR studies of active-site properties of human carbonic anhydrase II by using (15) N-labeled 4-methylimidazole as a local probe and histidine hydrogen-bond correlations

By using a combination of liquid and solid-state NMR spectroscopy, (15) N-labeled 4-methylimidazole (4-MI) as a local probe of the environment has been studied: 1) in the polar, wet Freon CDF3 /CDF2 Cl down to 130 K, 2) in water at pH 12, and 3) in solid samples of the mutant H64A of human carbonic anhydrase II (HCA II). In the latter, the active-site His64 residue is replaced by alanine; the catalytic activity is, however, rescued by the presence of 4-MI. For the Freon solution, it is demonstrated that addition of water molecules not only catalyzes proton tautomerism but also lifts its quasidegeneracy. The possible hydrogen-bond clusters formed and the mechanism of the tautomerism are discussed. Information about the imidazole hydrogen-bond geometries is obtained by establishing a correlation between published (1) H and (15) N chemical shifts of the imidazole rings of histidines in proteins. This correlation is useful to distinguish histidines embedded in the interior of proteins and those at the surface, embedded in water. Moreover, evidence is obtained that the hydrogen-bond geometries of His64 in the active site of HCA II and of 4-MI in H64A HCA II are similar. Finally, the degeneracy of the rapid tautomerism of the neutral imidazole ring His64 reported by Shimahara et al. (J. Biol. Chem.- 2007, 282, 9646) can be explained with a wet, polar, nonaqueous active-site conformation in the inward conformation, similar to the properties of 4-MI in the Freon solution. The biological implications for the enzyme mechanism are discussed.

Protonation states of the tryptophan synthase internal aldimine active site from solid-state NMR spectroscopy: direct observation of the protonated Schiff base linkage to pyridoxal-5'-phosphate

The acid-base chemistry that drives catalysis in pyridoxal-5'-phosphate (PLP)-dependent enzymes has been the subject of intense interest and investigation since the initial identification of PLP's role as a coenzyme in this extensive class of enzymes. It was first proposed over 50 years ago that the initial step in the catalytic cycle is facilitated by a protonated Schiff base form of the holoenzyme in which the linking lysine ε-imine nitrogen, which covalently binds the coenzyme, is protonated. Here we provide the first (15)N NMR chemical shift measurements of such a Schiff base linkage in the resting holoenzyme form, the internal aldimine state of tryptophan synthase. Double-resonance experiments confirm the assignment of the Schiff base nitrogen, and additional (13)C, (15)N, and (31)P chemical shift measurements of sites on the PLP coenzyme allow a detailed model of coenzyme protonation states to be established.

Aspartate aminotransferase: an old dog teaches new tricks

Aspartate aminotransferase (AAT) is a prototypical pyridoxal 5'-phosphate (PLP) dependent enzyme that catalyzes the reversible interconversion of l-aspartate and α-ketoglutarate with oxalacetate and l-glutamate via a ping-pong catalytic cycle in which the pyridoxamine 5'-phosphate enzyme form is an intermediate. There is a bountiful literature on AAT that spans approximately 60years, and much fundamental mechanistic information on PLP dependent reactions has been gained from its study. Here, we review our recent work on AAT, where we again used it as a test bed for fundamental concepts in PLP chemistry. First, we discuss the role that coenzyme protonation state plays in controlling reaction specificity, then ground state destabilization via hyperconjugation in the external aldimine intermediate is examined. The third topic is light enhancement of catalysis of Cα-H deprotonation by PLP in solution and in AAT, which occurs through a triplet state of the external aldimine intermediate. Lastly, we consider recent advances in our analyses of enzyme multiple sequence alignments for the purpose of predicting mutations that are required to interconvert structurally similar but catalytically distinct enzymes, and the application of our program JANUS to the conversion of AAT into tyrosine aminotransferase.