Computational Investigation of Bisphosphate Inhibitors of 3-Deoxy-d-manno-octulosonate 8-phosphate Synthase

Role of UDP-N-acetylmuramic acid in the regulation of MurA activity revealed by molecular dynamics simulations

The peptidoglycan biosynthesis pathway plays a vital role in bacterial cells, and facilitates peptidoglycan layer formation, a fundamental structural component of the bacterial cell wall. The enzymes in this pathway are candidates for antibiotic development, as most do not have mammalian homologues. The UDP-N-acetylglucosamine (UNAG) enolpyruvyl transferase enzyme (MurA) in the peptidoglycan pathway cytoplasmic step is responsible for the phosphoenolpyruvate (PEP)-UNAG catalytic reaction, forming UNAG enolpyruvate and inorganic phosphate. Reportedly, UDP-N-acetylmuramic acid (UNAM) binds tightly to MurA forming a dormant UNAM-PEP-MurA complex and acting as a MurA feedback inhibitor. MurA inhibitors are complex, owing to competitive binding interactions with PEP, UNAM, and UNAG at the MurA active site. We used computational methods to explore UNAM and UNAG binding. UNAM showed stronger hydrogen-bond interactions with the Arg120 and Arg91 residues, which help to stabilize the closed conformation of MurA, than UNAG. Binding free energy calculations using end-point computational methods showed that UNAM has a higher binding affinity than UNAG, when PEP is attached to Cys115. The unbinding process, simulated using τ-random acceleration molecular dynamics, showed that UNAM has a longer relative residence time than UNAG, which is related to several complex dissociation pathways, each with multiple intermediate metastable states. This prevents the loop from opening and exposing the Arg120 residue to accommodate UNAG and potential new ligands. Moreover, we demonstrate the importance of Cys115-linked PEP in closed-state loop stabilization. We provide a basis for evaluating novel UNAM analogues as potential MurA inhibitors. PUBLIC SIGNIFICANCE: MurA is a critical enzyme involved in bacterial cell wall biosynthesis and is involved in antibiotic resistance development. UNAM can remain in the target protein's active site for an extended time compared to its natural substrate, UNAG. The prolonged interaction of this highly stable complex known as the 'dormant complex' comprises UNAM-PEP-MurA and offers insights into antibiotic development, providing potential options against drug-resistant bacteria and advancing our understanding of microbial biology.

Computational insights for predicting the binding and selectivity of peptidomimetic plasmepsin IV inhibitors against cathepsin D

Plasmepsins (Plms) are aspartic proteases involved in the degradation of human hemoglobin by P. falciparum and are essential for the survival and growth of the parasite. Therefore, Plm enzymes are reported as an important antimalarial drug target. Herein, we have applied molecular docking, molecular dynamics (MD) simulations, and binding free energy with the Linear Interaction Energy (LIE) approach to investigate the binding of peptidomimetic PlmIV inhibitors with a particular focus on understanding their selectivity against the human Asp protease cathepsin D (CatD). The residual decomposition analysis results suggest that amino acid differences in the subsite S3 of PlmIV and CatD are responsible for the higher selectivity of the 5a inhibitor. These findings yield excellent agreement with experimental binding data and provide new details regarding van der Waals and electrostatic interactions of subsite residues as well as structural properties of the PlmIV and CatD systems.

Role of Half-of-Sites Reactivity and Inter-Subunit Communications in DAHP Synthase Catalysis and Regulation

α-Carboxyketose synthases, including 3-deoxy-d-arabinoheptulosonate 7-phosphate synthase (DAHPS), are long-standing targets for inhibition. They are challenging targets to create tight-binding inhibitors against, and inhibitors often display half-of-sites binding and partial inhibition. Half-of-sites inhibition demonstrates the existence of inter-subunit communication in DAHPS. We used X-ray crystallography and spatially resolved hydrogen-deuterium exchange (HDX) to reveal the structural and dynamic bases for inter-subunit communication in Escherichia coli DAHPS(Phe), the isozyme that is feedback-inhibited by phenylalanine. Crystal structures of this homotetrameric (dimer-of-dimers) enzyme are invariant over 91% of its sequence. Three variable loops make up 8% of the sequence and are all involved in inter-subunit contacts across the tight-dimer interface. The structures have pseudo-twofold symmetry indicative of inter-subunit communication across the loose-dimer interface, with the diagonal subunits B and C always having the same conformation as each other, while subunits A and D are variable. Spatially resolved HDX reveals contrasting responses to ligand binding, which, in turn, affect binding of the second substrate, erythrose-4-phosphate (E4P). The N-terminal peptide, M1-E12, and the active site loop that binds E4P, F95-K105, are key parts of the communication network. Inter-subunit communication appears to have a catalytic role in all α-carboxyketose synthase families and a regulatory role in some members.

Advances in UDP-N-Acetylglucosamine Enolpyruvyl Transferase (MurA) Covalent Inhibition

Peptidoglycan is a cross-linked polymer responsible for maintaining the bacterial cell wall integrity and morphology in Gram-negative and Gram-positive bacteria. The peptidoglycan pathway consists of the enzymatic reactions held in three steps: cytoplasmic, membrane-associated, and periplasmic. The Mur enzymes (MurA-MurF) are involved in a cytoplasmic stage. The UDP-N-acetylglucosamine enolpyruvyl transferase (MurA) enzyme is responsible for transferring the enolpyruvate group from phosphoenolpyruvate (PEP) to UDP-N-acetylglucosamine (UNAG) to form UDP-N-acetylglucosamine enolpyruvate (EP-UNAG). Fosfomycin is a natural product analogous to PEP that acts on the MurA target enzyme via binding covalently to the key cysteine residue in the active site. Similar to fosfomycin, other MurA covalent inhibitors have been described with a warhead in their structure that forms a covalent bond with the molecular target. In MurA, the nucleophilic thiolate of Cys115 is pointed as the main group involved in the warhead binding. Thus, in this minireview, we briefly describe the main recent advances in the design of MurA covalent inhibitors.

In silico investigation and surmounting of Lipopolysaccharide barrier in Gram-Negative Bacteria: How far has molecular dynamics Come?

Lipopolysaccharide (LPS), a main component of the outer membrane of Gram-negative bacteria, has crucial implications on both antibiotic resistance and the overstimulation of the host innate immune system. Fighting against these global concerns calls for the molecular understanding of the barrier function and immunostimulatory ability of LPS. Molecular dynamics (MD) simulations have become an invaluable tool for uncovering important findings in LPS research. While the reach of MD simulations for investigating the immunostimulatory ability of LPS has been already outlined, little attention has been paid to the role of MD simulations for exploring its barrier function and synthesis. Herein, we give an overview about the impact of MD simulations on gaining insight into the shield role and synthesis pathway of LPS, which have attracted considerable attention to discover molecules able to surmount antibiotic resistance, either circumventing LPS defenses or disrupting its synthesis. We specifically focus on the enhanced sampling and free energy calculation methods that have been combined with MD simulations to address such research. We also highlight the use of special-purpose MD supercomputers, the importance of appropriate LPS and ions parameterization to obtain reliable results, and the complementary views that MD and wet-lab experiments provide. Thereby, this work, which covers the last five years of research, apart from outlining the phenomena and strategies that are being explored, evidences the valuable insights that are gained by MD, which may be useful to advance antibiotic design, and what the prospects of this in silico method could be in LPS research.

An Inhibitor-in-Pieces Approach to DAHP Synthase Inhibition: Potent Enzyme and Bacterial Growth Inhibition

3-Deoxy-d-arabinoheptulosonate-7-phosphate (DAHP) synthase catalyzes the first step in the shikimate biosynthetic pathway and is an antimicrobial target. We used an inhibitor-in-pieces approach, based on the previously reported inhibitor DAHP oxime, to screen inhibitor fragments in the presence and absence of glycerol 3-phosphate to occupy the distal end of the active site. This led to DAHP hydrazone, the most potent inhibitor to date, K_i = 10 ± 1 nM. Three trifluoropyruvate (TFP)-based inhibitor fragments were efficient inhibitors with ligand efficiencies of up to 0.7 kcal mol^-1/atom compared with 0.2 kcal mol^-1/atom for a typical good inhibitor. The crystal structures showed the TFP-based inhibitors binding upside down in the active site relative to DAHP oxime, providing new avenues for inhibitor development. The ethyl esters of TFP oxime and TFP semicarbazone prevented E. coli growth in culture with IC₅₀ = 0.21 ± 0.01 and 0.77 ± 0.08 mg mL^-1, respectively. Overexpressing DAHP synthase relieved growth inhibition, demonstrating that DAHP synthase was the target. Growth inhibition occurred in media containing aromatic amino acids, suggesting that growth inhibition was due to depletion of some other product(s) of the shikimate pathway, possibly folate.

Targeting shikimate pathway: In silico analysis of phosphoenolpyruvate derivatives as inhibitors of EPSP synthase and DAHP synthase

The shikimate pathway consists of seven enzymatic steps involved in the conversion of erythrose-4-phosphate and phosphoenolpyruvate to chorismate and also responsible to the production of aromatic amino acids, such as phenylalanine, tyrosine, and tryptophan which are essential to the bacterial metabolism. The 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase (DAHPS) and 5-enolpyruvylshikimate 3-phosphate synthase (EPSPS) catalyze important steps in the shikimate pathway using as substrate the phosphoenolpyruvate (PEP). Due to the importance of PEP in shikimate pathway, its structure has been investigated to develop new bioinspired competitive inhibitors against DAHPS and EPSPS. In the present study, we perform a literature survey of 28 PEP derivatives, then we analyzed the selectivity and affinity of these compounds against the EPSPS and DAHPS structures using consensual molecular docking, pharmacophore prediction, molecular dynamics (MD) simulations, and binding free energy calculations. Here, we propose consistent binding modes of the selected ligands and indicate that their structures show interesting pharmacophoric properties related to multi-targets inhibitors for both enzymes. Our computational results are supported by previous experimental findings related to the interactions of PEP derivatives with DAHPS and EPSPS structures.

Evaluating the Performance of a Non-Bonded Cu2+ Model Including Jahn-Teller Effect into the Binding of Tyrosinase Inhibitors

Tyrosinase (TYR) is a metalloenzyme classified as a type-3 copper protein, which is involved in the synthesis of melanin through a catalytic process beginning with the conversion of the amino acid l-Tyrosine (l-Tyr) to l-3,4-dihydroxyphenylalanine (l-DOPA). It plays an important role in the mechanism of melanogenesis in various organisms including mammals, plants, and fungi. Herein, we used a combination of computational molecular modeling techniques including molecular dynamic (MD) simulations and the linear interaction energy (LIE) model to evaluate the binding free energy of a set of analogs of kojic acid (KA) in complex with TYR. For the MD simulations, we used a dummy model including the description of the Jahn–Teller effect for Cu²⁺ ions in the active site of this enzyme. Our results show that the LIE model predicts the TYR binding affinities of the inhibitor in close agreement to experimental results. Overall, we demonstrate that the classical model provides a suitable description of the main interactions between analogs of KA and Cu²⁺ ions in the active site of TYR.