The Catalytic Mechanism of the Pyridoxal-5′-phosphate-Dependent Enzyme, Histidine Decarboxylase: A Computational Study

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.

Rapidly and accurately screening histidine decarboxylase inhibitors from Radix Paeoniae alba using ultrafiltration-high performance liquid chromatography/mass spectrometry combined with enzyme channel blocking and directional enrichment technique

Histidine Decarboxylase (HDC), an unique enzyme responsible for the synthesis of histamine, which is an important mediator in allergy. Inhibition of HDC activity to decrease histamine production is one way to alleviate allergic symptoms. Traditional Chinese medicines (TCMs) with reported anti-allergy effect is one of important source to search for natural HDC inhibitor. Ultrafiltration combined with high-performance liquid chromatography/mass spectrometry (UF-HPLC/MS) is an effective method for screening HDC inhibitor from TCMs. Nevertheless, false-positive and false-negative results caused by the non-specific binding and the neglection of the trace active compounds are major problems in this method. In this study, an integrated strategy that combined UF-HPLC/MS with enzyme channel blocking (ECB) technique and directional enrichment (DE) technique was developed to seek natural HDC inhibitors from Radix Paeoniae alba (RPA), and at the same time, to reduce false-positive and false-negative results. HDC activity was detected to determine the validity of the screened compounds by RP-HPLC-FD in vitro. Molecular docking was applied to assay the binding affinity and binding sites. As a result, three compounds were screened from low content components of RPA after the DE. Among them, two non-specific compounds were eliminated by ECB, and the specific compound was identified as catechin, which has obvious HDC inhibition activity with IC₅₀ 0.52 mM. Furthermore, gallic acid (IC₅₀ 1.8 mM) and paeoniflorin (IC₅₀>2 mM) from high content components of RPA were determined having HDC inhibitory activity. In conclusion, the integrated strategy of UF-HPLC/MS combined with ECB and DE technique is an effective mode for rapid and accurate screening and identification of natural HDC inhibitors from TCMs.

Catalytic Mechanism of Pyridoxal 5'-Phosphate-Dependent Aminodeoxychorismate Lyase: A Computational QM/MM Study

Aminodeoxychorismate lyase (ADCL) is a kind of pyridoxal-5'-phosphate (PLP)-dependent enzyme that catalyzes the conversion of 4-amino-4-deoxychorismate (ADC) to p-aminobenzoate (PABA), which is a key step for the biosynthesis of folate. To illuminate the reaction details at the atomistic level, an enzyme-substrate reactant model has been constructed, and QM/MM calculations have been performed. Our calculation results reveal that the overall catalytic cycle contains 11 elementary steps, which can be described by three stages, including the transamination reaction of PLP, the release of pyruvate and aromatization of ADC, and the recovery to the initial aldimine. During the reaction, a series of intramolecular proton transfer are involved, which are the key for the C-N bond formation and cleavage as well as the aromatization of the ADC ring. In addition to forming the Schiff base with the pocket residue Lys251 and substrate in the internal aldimine and the external aldimine, respectively, the coenzyme PLP also plays a critical role in the intramolecular proton transfer by employing its hydroxyl oxygen anion and phosphate group. These findings may provide useful information for further understanding the catalytic mechanism of other PLP-dependent enzymes.

Improvement of Fusel Alcohol Production by Engineering of the Yeast Branched-Chain Amino Acid Aminotransaminase

Branched-chain higher alcohols (BCHAs), or fusel alcohols, including isobutanol, isoamyl alcohol, and active amyl alcohol, are useful compounds in several industries. The yeast Saccharomyces cerevisiae can synthesize these compounds via the metabolic pathways of branched-chain amino acids (BCAAs). Branched-chain amino acid aminotransaminases (BCATs) are the key enzymes for BCHA production via the Ehrlich pathway of BCAAs. BCATs catalyze a bidirectional transamination reaction between branched-chain α-keto acids (BCKAs) and BCAAs. In S. cerevisiae, there are two BCAT isoforms, Bat1 and Bat2, which are encoded by the genes BAT1 and BAT2. Although many studies have shown the effects of deletion or overexpression of BAT1 and BAT2 on BCHA production, there have been no reports on the enhancement of BCHA production by functional variants of BCATs. Here, to improve BCHA productivity, we designed variants of Bat1 and Bat2 with altered enzyme activity by using in silico computational analysis: the Gly333Ser and Gly333Trp Bat1 and corresponding Gly316Ser and Gly316Trp Bat2 variants, respectively. When expressed in S. cerevisiae cells, most of these variants caused a growth defect in minimal medium. Interestingly, the Gly333Trp Bat1 and Gly316Ser Bat2 variants achieved 18.7-fold and 17.4-fold increases in isobutanol above that for the wild-type enzyme, respectively. The enzyme assay revealed that the catalytic activities of all four BCAT variants were lower than that of the wild-type enzyme. Our results indicate that the decreased BCAT activity enhanced BCHA production by reducing BCAA biosynthesis, which occurs via a pathway that directly competes with BCHA production.

IMPORTANCE Recently, several studies have attempted to increase the production of branched-chain higher alcohols (BCHAs) in the yeast Saccharomyces cerevisiae. The key enzymes for BCHA biosynthesis in S. cerevisiae are the branched-chain amino acid aminotransaminases (BCATs) Bat1 and Bat2. Deletion or overexpression of the genes encoding BCATs has an impact on the production of BCHAs; however, amino acid substitution variants of Bat1 and Bat2 that could affect enzymatic properties—and ultimately BCHA productivity—have not been fully studied. By using in silico analysis, we designed variants of Bat1 and Bat2 and expressed them in yeast cells. We found that the engineered BCATs decreased catalytic activities and increased BCHA production. Our approach provides new insight into the functions of BCATs and will be useful in the future construction of enzymes optimized for high-level production of BCHAs.

The Nitrogen Atom of Vitamin B6 Is Essential for the Catalysis of Radical Aminomutases

Radical aminomutases are pyridoxal 5'-phosphate (PLP, a B₆ vitamer)-dependent enzymes that require the generation of a 5'-deoxyadenosyl radical to initiate the catalytic cycle, to perform a 1,2 amino group shift reaction. The role of the nitrogen atom of PLP in radical aminomutases has not been investigated extensively yet. We report an alternative synthetic procedure to provide easy access to 1-deazaPLP (dAPLP), an isosteric analog of PLP which acts as a probe for studying the role of the nitrogen atom. Our results revealed that lysine 5,6-aminomutase (5,6-LAM), a radical aminomutase, reconstituted with dAPLP cannot turn over a substrate, demonstrating that the nitrogen atom is essential for radical aminomutases. In contrast, biochemical and spectroscopic studies on the S238A variant reconstituted with PLP revealed a minuscule loss of activity. This apparent anomaly can be explained by a water-mediated rescue of activity in S238A, as if mimicking the active site of lysine 2,3-aminomutase. This study leads to a better comprehension of how enzymes harness the optimum capability of PLP to realize catalysis.

Modern strategies for the diversification of the supply of natural compounds - the case of alkaloid painkillers

Plant-derived natural compounds are used for treating diseases since the beginning of humankind. The supply of many plant-derived natural compounds for medicinal purposes, such as thebaine, morphine, and codeine, is, nowadays, majorly dependent on opium poppy crop harvesting. This dependency puts an extra risk factor in ensuring the supply chain because crops are highly susceptible to environmental factors. Emerging technologies, such as biocatalysis, might help to solve this problem, by diversifying the sources of supply of these compounds. Here we review the first committed step in the production of alkaloid painkillers, the production of S-norcoclaurine, and the enzymes involved. The improvement of these enzymes can be carried out by experimental directed evolution and rational design strategies, supported by computational methods, to create variants that produce the S-norcoclaurine precursor for alkaloid painkillers in heterologous organisms, meeting the pharmaceutical industry standards and needs without depending on opium poppy crops.

The Catalytic Mechanism of Pdx2 Glutaminase Driven by a Cys-His-Glu Triad: A Computational Study

The catalytic mechanism of Pdx2 was studied with atomic detail employing the computational ONIOM hybrid QM/MM methodology. Pdx2 employs a Cys-His-Glu catalytic triad to deaminate glutamine to glutamate and ammonia - the source of the nitrogen of pyridoxal 5'-phosphate (PLP). This enzyme is, therefore, a rate-limiting step in the PLP biosynthetic pathway of Malaria and Tuberculosis pathogens that rely on this mechanism to obtain PLP. For this reason, Pdx2 is considered a novel and promising drug target to treat these diseases. The results obtained show that the catalytic mechanism of Pdx2 occurs in six steps that can be divided into four stages: (i) activation of Cys₈₇ , (ii) deamination of glutamine with the formation of the glutamyl-thioester intermediate, (iii) hydrolysis of the formed intermediate, and (iv) enzymatic turnover. The kinetic data available in the literature (19.1-19.5 kcal mol^-1 ) agree very well with the calculated free energy barrier of the hydrolytic step (18.2 kcal.mol^-11 ), which is the rate-limiting step of the catalytic process when substrate is readily available in the active site. This catalytic mechanism differs from other known amidases in three main points: i) it requires the activation of the nucleophile Cys₈₇ to a thiolate; ii) the hydrolysis occurs in a single step and therefore does not require the formation of a second tetrahedral reaction intermediate, as it is proposed, and iii) Glu₁₉₈ does not have a direct role in the catalytic process. Together, these results can be used for the synthesis of new transition state analogue inhibitors capable of inhibiting Pdx2 and impair diseases like Malaria and Tuberculosis.

New insights into the catalytic mechanism of the SARS-CoV-2 main protease: an ONIOM QM/MM approach

SARS-CoV-2 M^pro, also known as the main protease or 3C-like protease, is a key enzyme involved in the replication process of the virus that is causing the COVID-19 pandemic. It is also the most promising antiviral drug target targeting SARS-CoV-2 virus. In this work, the catalytic mechanism of M^pro was studied using the full model of the enzyme and a computational QM/MM methodology with a 69/72-atoms QM region treated at DLPNO-CCSD(T)/CBS//B3LYP/6-31G(d,p):AMBER level and including the catalytic important oxyanion-hole residues. The transition state of each step was fully characterized and described together with the related reactants and products. The rate-limiting step of the catalytic process is the hydrolysis of the thioester-enzyme adduct, and the calculated barrier closely agrees with the available kinetic data. The calculated Gibbs free energy profile, together with the full atomistic detail of the structures involved in catalysis, can now serve as valuable models for the rational drug design of transition state analogs as new inhibitors targeting the SARS-CoV-2 virus.

Structural basis for substrate specificity of l-methionine decarboxylase

l -Methionine decarboxylase (MetDC) from Streptomyces sp. 590 is a vitamin B₆ -dependent enzyme and catalyzes the non-oxidative decarboxylation of l -methionine to produce 3-methylthiopropylamine and carbon dioxide. We present here the crystal structures of the ligand-free form of MetDC and of several enzymatic reaction intermediates. Group II amino acid decarboxylases have many residues in common around the active site but the residues surrounding the side chain of the substrate differ. Based on information obtained from the crystal structure, and mutational and biochemical experiments, we propose a key role for Gln64 in determining the substrate specificity of MetDC, and for Tyr421 as the acid catalyst that participates in protonation after the decarboxylation reaction.

The complete catalytic mechanism of xanthine oxidase: a computational study

In this article, quantum mechanical/molecular mechanical (QM/MM) methods were used to study the full catalytic mechanism of xanthine oxidase (XO).

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.

Structure and substrate specificity determinants of the taurine biosynthetic enzyme cysteine sulphinic acid decarboxylase

Pyridoxal 5́-phosphate (PLP) is an important cofactor for amino acid decarboxylases with many biological functions, including the synthesis of signalling molecules, such as serotonin, dopamine, histamine, γ-aminobutyric acid, and taurine. Taurine is an abundant amino acid with multiple physiological functions, including osmoregulation, pH regulation, antioxidative protection, and neuromodulation. In mammalian tissues, taurine is mainly produced by decarboxylation of cysteine sulphinic acid to hypotaurine, catalysed by the PLP-dependent cysteine sulphinic acid decarboxylase (CSAD), followed by oxidation of the product to taurine. We determined the crystal structure of mouse CSAD and compared it to other PLP-dependent decarboxylases in order to identify determinants of substrate specificity and catalytic activity. Recognition of the substrate involves distinct side chains forming the substrate-binding cavity. In addition, the backbone conformation of a buried active-site loop appears to be a critical determinant for substrate side chain binding in PLP-dependent decarboxylases. Phe94 was predicted to affect substrate specificity, and its mutation to serine altered both the catalytic properties of CSAD and its stability. Using small-angle X-ray scattering, we further showed that CSAD presents open/close motions in solution. The structure of apo-CSAD indicates that the active site gets more ordered upon internal aldimine formation. Taken together, the results highlight details of substrate recognition in PLP-dependent decarboxylases and provide starting points for structure-based inhibitor design with the aim of affecting the biosynthesis of taurine and other abundant amino acid metabolites.

A single amino acid substitution converts a histidine decarboxylase to an imidazole acetaldehyde synthase

Histidine decarboxylase (HDC; EC 4.1.1.22), an enzyme that catalyzes histamine synthesis with high substrate specificity, is a member of the group II pyridoxal 5'-phosphate (PLP) -dependent decarboxylase family. Tyrosine is a conserved residue among group II PLP-dependent decarboxylases. Human HDC has a Y334 located on a catalytically important loop at the active site. In this study, we demonstrated that a HDC Y334F mutant is capable of catalyzing the decarboxylation-dependent oxidative deamination of histidine to yield imidazole acetaldehyde. Replacement of the active-site Tyr with Phe in group II PLP-dependent decarboxylases, including mammalian aromatic amino acid decarboxylase, plant tyrosine/DOPA decarboxylase, and plant tryptophan decarboxylase, is expected to result in the same functional change, given that a Y-to-F substitution at the corresponding residue (number 260) in the HDC of Morganella morganii, another group II PLP-dependent decarboxylase, yielded the same effect. Thus, it was suggested that the loss of the OH moiety from the active-site Tyr residue of decarboxylase uniquely converts the enzyme to an aldehyde synthase.

User-Friendly Quantum Mechanics: Applications for Drug Discovery

Quantum mechanics (QM) methods provide a fine description of receptor-ligand interactions and of chemical reactions. Their use in drug design and drug discovery is increasing, especially for complex systems including metal ions in the binding sites, for the design of highly selective inhibitors, for the optimization of bi-specific compounds, to understand enzymatic reactions, and for the study of covalent ligands and prodrugs. They are also used for generating molecular descriptors for predictive QSAR/QSPR models and for the parameterization of force fields. Thanks to the continuous increase of computational power offered by GPUs and to the development of sophisticated algorithms, QM methods are becoming part of the standard tools used in computer-aided drug design (CADD). We present the most used QM methods and software packages, and we discuss recent representative applications in drug design and drug discovery.

Expression of Histidine Decarboxylase and Its Roles in Inflammation

Histamine is a well-known mediator of inflammation that is released from mast cells and basophils. To date, many studies using histamine receptor antagonists have shown that histamine acts through four types of receptors: H1, H2, H3, and H4. Thus, histamine plays more roles in various diseases than had been predicted. However, our knowledge about histamine-producing cells and the molecular mechanisms underlying histamine production at inflammatory sites is still incomplete. The histamine producing enzyme, histidine decarboxylase (HDC), is commonly induced at inflammatory sites during the late and chronic phases of both allergic and non-allergic inflammation. Thus, histamine levels in tissues are maintained at effective concentrations for hours, enabling the regulation of various functions through the production of cytokines/chemokines/growth factors. Understanding the regulation of histamine production will allow the development of a new strategy of using histamine antagonists to treat inflammatory diseases.

Hyperconjugation promoted by hydrogen bonding between His98/His241 and a carboxyl group contributes to tyrosine decarboxylase catalysis

This article first demonstrates the utility of the σ → π* hyperconjugation by His98 and His241 in pyridoxal 5′-phosphate-dependent tyrosine decarboxylase.

PLP-dependent enzymes as important biocatalysts for the pharmaceutical, chemical and food industries: a structural and mechanistic perspective

PLP-dependent enzymes described on this review are attractive targets for enzyme engineering towards their application in an industrial biotechnology framework.

Lactobacillus brevis CGMCC 1306 glutamate decarboxylase: Crystal structure and functional analysis

Glutamate decarboxylase (GAD), which is a unique pyridoxal 5-phosphate (PLP)-dependent enzyme, can catalyze α-decarboxylation of l-glutamate (L-Glu) to γ-aminobutyrate (GABA). The crystal structure of GAD in complex with PLP from Lactobacillus brevis CGMCC 1306 was successfully solved by molecular-replacement, and refined at 2.2 Å resolution to an R_work factor of 18.76% (R_free = 23.08%). The coenzyme pyridoxal 5-phosphate (PLP) forms a Schiff base with the active-site residue Lys279 by continuous electron density map, which is critical for catalysis by PLP-dependent decarboxylase. Gel filtration showed that the active (pH 4.8) and inactive (pH 7.0) forms of GAD are all dimer. The residues (Ser126, Ser127, Cys168, Ile211, Ser276, His278 and Ser321) play important roles in anchoring PLP cofactor inside the active site and supporting its catalytic reactivity. The mutant T215A around the putative substrate pocket displayed an 1.6-fold improvement in catalytic efficiency (k_cat/K_m) compared to the wild-type enzyme (1.227 mM^-1 S^-1 versus 0.777 mM^-1 S^-1), which was the highest activity among all variants tested. The flexible loop (Tyr308-Glu312), which is positioned near the substrate-binding site, is involved in the catalytic reaction, and the conserved residue Tyr308 plays a vital role in decarboxylation of L-Glu.

molUP: A VMD plugin to handle QM and ONIOM calculations using the gaussian software

The notable advances obtained by computational (bio)chemistry provided its widespread use in many areas of science, in particular, in the study of reaction mechanisms. These studies involve a huge number of complex calculations, which are often carried out using the Gaussian suite of programs. The preparation of input files and the analysis of the output files are not easy tasks and often involve laborious and complex steps. Taking this into account, we developed molUP: a VMD plugin that offers a complete set of tools that enhance the preparation of QM and ONIOM (QM/MM, QM/QM, and QM/QM/MM) calculations. The starting structures for these calculations can be imported from different chemical formats. A set of tools is available to help the user to examine or modify any geometry parameter. This includes the definition of layers in ONIOM calculations, choosing fixed atoms during geometry optimizations, the recalculation or adjustment of the atomic charges, performing SCANs or IRC calculations, etc. molUP also extracts the geometries from the output files as well as the energies of each of them. All of these tasks are performed in an interactive GUI that is extremely helpful for the user. MolUP was developed to be easy to handle by inexperienced users, but simultaneously to be a fast and flexible graphical interface to allow the advanced users to take full advantage of this plugin. The program is available, free of charges, for macOS, Linux, and Windows at the PortoBioComp page https://www.fc.up.pt/PortoBioComp/database/doku.php?id=molup. © 2018 Wiley Periodicals, Inc.