Free energy landscapes of prototropic tautomerism in pyridoxal 5′-phosphate schiff bases at the active site of an enzyme in aqueous medium

In Silico and In Vitro Potential of FDA-Approved Drugs for Antimalarial Drug Repurposing against Plasmodium Serine Hydroxymethyltransferases

Malaria has spread in many countries, with a 12% increase in deaths after the coronavirus disease 2019 pandemic. Malaria is one of the most concerning diseases in the Greater Mekong subregion, showing increased drug-resistant rates. Serine hydroxymethyltransferase (SHMT), a key enzyme in the deoxythymidylate synthesis pathway, has been identified as a promising antimalarial drug target due to its conserved folate binding pocket. This study used a molecular docking approach to screen 2509 US Food and Drug Administration (FDA)-approved drugs against seven Plasmodium SHMT structures. Eight compounds had significantly lower binding energies than the known SHMT inhibitors pyrazolopyran(+)-86, tetrahydrofolate, and antimalarial drugs, ranging from 4 to 10 kcal/mol. Inhibition assays testing the eight compounds against Plasmodium falciparum SHMT (PfSHMT) showed that amphotericin B was a competitive inhibitor of PfSHMT with a half-maximal inhibitory concentration (IC50) of 106 ± 1 μM. Therefore, a 500 ns molecular dynamics simulation of PfSHMT/PLS/amphotericin B was performed. The backbone root-mean-square deviation of the protein-ligand complex indicated the high complex stability during simulations, supported by its radius of gyration, hydrogen-bond interactions, and number of atom contacts. The appreciable binding affinity of amphotericin B for PfSHMT was indicated by their solvated interaction energy (-11.15 ± 0.09 kcal/mol) and supported by strong ligand-protein interactions (≥80% occurrences) with its essential residues (i.e., Y78, K151, N262, F266, and V365) predicted by pharmacophore modeling and per-residue decomposition free energy methods. Therefore, our findings identify a promising new PfSHMT inhibitor, albeit with less inhibitory activity, and suggest a core structure that differs from that of previous SHMT inhibitors, thus being a rational approach for novel antimalarial drug design.

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.

The yeast ALA synthase C-terminus positively controls enzyme structure and function

5-Aminolevulinic acid synthase (ALAS) is a pyridoxal 5'-phosphate (PLP)-dependent enzyme that catalyzes the first and rate-limiting step of heme biosynthesis in α-proteobacteria and several non-plant eukaryotes. All ALAS homologs contain a highly conserved catalytic core, but eukaryotes also have a unique C-terminal extension that plays a role in enzyme regulation. Several mutations in this region are implicated in multiple blood disorders in humans. In Saccharomyces cerevisiae ALAS (Hem1), the C-terminal extension wraps around the homodimer core to contact conserved ALAS motifs proximal to the opposite active site. To determine the importance of these Hem1 C-terminal interactions, we determined the crystal structure of S. cerevisiae Hem1 lacking the terminal 14 amino acids (Hem1 ΔCT). With truncation of the C-terminal extension, we show structurally and biochemically that multiple catalytic motifs become flexible, including an antiparallel β-sheet important to Fold-Type I PLP-dependent enzymes. The changes in protein conformation result in an altered cofactor microenvironment, decreased enzyme activity and catalytic efficiency, and ablation of subunit cooperativity. These findings suggest that the eukaryotic ALAS C-terminus has a homolog-specific role in mediating heme biosynthesis, indicating a mechanism for autoregulation that can be exploited to allosterically modulate heme biosynthesis in different organisms.

A Multiple Proton Transfer Mechanism for the Charging Step of the Aminoacylation Reaction at the Active Site of Aspartyl tRNA Synthetase

Aspartyl-tRNA synthetase catalyzes the attachment of aspartic acid to its cognate tRNA by the aminoacylation reaction during the initiation of the protein biosynthesis process. In the second step of the aminoacylation reaction, known as the charging step, the aspartate moiety is transferred from aspartyl-adenylate to the 3'-OH of A76 of tRNA through a proton transfer process. We have investigated different pathways for the charging step through three separate QM/MM simulations combined with the enhanced sampling method of well-sliced metadynamics and found out the most feasible pathway for the reaction at the active site of the enzyme. In the charging reaction, both the phosphate group and the ammonium group after deprotonation can potentially act as a base for proton transfer in the substrate-assisted mechanism. We have considered three possible mechanisms involving different pathways of proton transfer, and only one of them is determined to be enzymatically feasible. The free energy landscape along reaction coordinates where the phosphate group acts as the general base showed that, in the absence of water, the barrier height is 52.6 kcal/mol. The free energy barrier is reduced to 39.7 kcal/mol when the active site water molecules are also treated quantum mechanically, thus allowing a water mediated proton transfer. The charging reaction involving the ammonium group of the aspartyl adenylate is found to follow a path where first a proton from the ammonium group moves to a water in the vicinity forming a hydronium ion (H₃O⁺) and NH₂ group. The hydronium ion subsequently passes the proton to the Asp233 residue, thus minimizing the chance of back proton transfer from hydronium to the NH₂ group. The neutral NH₂ group subsequently takes the proton from the O3' of A76 with a free energy barrier of 10.7 kcal/mol. In the next step, the deprotonated O3' makes a nucleophilic attack to the carbonyl carbon forming a tetrahedral transition state with a free energy barrier of 24.8 kcal/mol. Thus, the present work shows that the charging step proceeds through a multiple proton transfer mechanism where the amino group formed after deprotonation acts as the base to capture a proton from O3' of A76 rather than the phosphate group. The current study also shows the important role played by Asp233 in the proton transfer process.

Free Energy Landscape of the Adenylation Reaction of the Aminoacylation Process at the Active Site of Aspartyl tRNA Synthetase

The process of protein biosynthesis is initiated by the aminoacylation process where a transfer ribonucleic acid (tRNA) is charged by the attachment of its cognate amino acid at the active site of the corresponding aminoacyl tRNA synthetase enzyme. The first step of the aminoacylation process, known as the adenylation reaction, involves activation of the cognate amino acid where it reacts with a molecule of adenosine triphosphate (ATP) at the active site of the enzyme to form the aminoacyl adenylate and inorganic pyrophosphate. In the current work, we have investigated the adenylation reaction between aspartic acid and ATP at the active site of the fully solvated aspartyl tRNA synthetase (AspRS) from Escherichia coli in aqueous medium at room temperature through hybrid quantum mechanical/molecular mechanical (QM/MM) simulations combined with enhanced sampling methods of well-tempered and well-sliced metadynamics. The objective of the present work is to study the associated free energy landscape and reaction barrier and also to explore the effects of active site mutation on the free energy surface of the reaction. The current calculations include finite temperature effects on free energy profiles. In particular, apart from contributions of interaction energies, the current calculations also include contributions of conformational, vibrational, and translational entropy of active site residues, substrates, and also the rest of the solvated protein and surrounding water into the free energy calculations. The present QM/MM metadynamics simulations predict a free energy barrier of 23.35 and 23.5 kcal/mol for two different metadynamics methods used to perform the reaction at the active site of the wild type enzyme. The free energy barrier increases to 30.6 kcal/mol when Arg217, which is an important conserved residue of the wild type enzyme at its active site, is mutated by alanine. These free energy results including the effect of mutation compare reasonably well with those of kinetic experiments that are available in the literature. The current work also provides molecular details of structural changes of the reactants and surroundings as the system dynamically evolves along the reaction pathway from reactant to the product state through QM/MM metadynamics simulations.

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.

Functional characterization of the HMP-P synthase of Legionella pneumophila (Lpg1565)

The production of the pyrimidine moiety in thiamine synthesis, 2-methyl-4-amino-5-hydroxymethylpyrimidine phosphate (HMP-P), has been described to proceed through the Thi5-dependent pathway in Saccharomyces cerevisiae and other yeast. Previous work found that ScThi5 functioned poorly in a heterologous context. Here we report a bacterial ortholog to the yeast HMP-P synthase (Thi5) was necessary for HMP synthesis in Legionella pneumophila. Unlike ScThi5, LpThi5 functioned in vivo in Salmonella enterica under multiple growth conditions. The protein LpThi5 is a dimer that binds pyridoxal-5'-phosphate (PLP), apparently without a solvent-exposed Schiff base. A small percentage of LpThi5 protein co-purifies with a bound molecule that can be converted to HMP. Analysis of variant proteins both in vivo and in vitro confirmed that residues in sequence motifs conserved across bacterial and eukaryotic orthologs modulate the function of LpThi5. IMPORTANCE: Thiamine is an essential vitamin for the vast majority of organisms. There are multiple strategies to synthesize and salvage this vitamin. The predominant pathway for synthesis of the pyrimidine moiety of thiamine involves the Fe-S cluster protein ThiC. An alternative pathway utilizes Thi5, a novel enzyme that uses PLP as a substrate. The Thi5-dependent pathway is poorly characterized in yeast and has not been characterized in Bacteria. Here we demonstrate that a Thi5-dependent pathway is necessary for thiamine biosynthesis in Legionella pneumophila and provide biochemical data to extend knowledge of the Thi5 enzyme, the corresponding biosynthetic pathway, and the role of metabolic network architecture in optimizing its function.

Tautomeric Equilibrium in Condensed Phases

We present an ab initio molecular dynamics (MD) investigation of the tautomeric equilibrium for the aqueous solutions of glycine and acetone under realistic experimental conditions. Metadynamics is used to accelerate proton migration among tautomeric centers. Due to the formation of complex water-ion structures involved in the proton dynamics in the aqueous environment, standard enhanced sampling approaches may face severe limitations in providing a general description of the phenomenon. Recently, we have developed a set of collective variables (CVs) designed to study protons transfer reactions in complex condensed systems [Grifoni, E. Proc. Natl. Acad. Sci. U.S.A. 2019, 116, 4054 4057]. In this work, we applied this approach to study proton dissociation dynamics leading to tautomeric interconversion of biologically and chemically relevant prototypical systems, namely, glycine and acetone in water. Although relatively simple from a chemical point of view, the results show that even for these small systems, complex reaction pathways and nontrivial conversion dynamics are observed. The generality of our method allows obtaining these results without providing any prior information on the dissociation dynamics but only the atomic species that can exchange protons in the process. Our results agree with literature estimates and demonstrate the general applicability of this method in the study of tautomeric reactions.

Direct observation of intermediates in the SufS cysteine desulfurase reaction reveals functional roles of conserved active-site residues

Iron-sulfur (Fe-S) clusters are necessary for the proper functioning of numerous metalloproteins. Fe-S cluster (Isc) and sulfur utilization factor (Suf) pathways are the key biosynthetic routes responsible for generating these Fe-S cluster prosthetic groups in Escherichia coli Although Isc dominates under normal conditions, Suf takes over during periods of iron depletion and oxidative stress. Sulfur acquisition via these systems relies on the ability to remove sulfur from free cysteine using a cysteine desulfurase mechanism. In the Suf pathway, the dimeric SufS protein uses the cofactor pyridoxal 5'-phosphate (PLP) to abstract sulfur from free cysteine, resulting in the production of alanine and persulfide. Despite much progress, the stepwise mechanism by which this PLP-dependent enzyme operates remains unclear. Here, using rapid-mixing kinetics in conjunction with X-ray crystallography, we analyzed the pre-steady-state kinetics of this process while assigning early intermediates of the mechanism. We employed H123A and C364A SufS variants to trap Cys-aldimine and Cys-ketimine intermediates of the cysteine desulfurase reaction, enabling direct observations of these intermediates and associated conformational changes of the SufS active site. Of note, we propose that Cys-364 is essential for positioning the Cys-aldimine for Cα deprotonation, His-123 acts to protonate the Ala-enamine intermediate, and Arg-56 facilitates catalysis by hydrogen bonding with the sulfhydryl of Cys-aldimine. Our results, along with previous SufS structural findings, suggest a detailed model of the SufS-catalyzed reaction from Cys binding to C-S bond cleavage and indicate that Arg-56, His-123, and Cys-364 are critical SufS residues in this C-S bond cleavage pathway.

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.