Identification of Potential Inhibitors of MurD Enzyme of Staphylococcus aureus from a Marine Natural Product Library

Staphylococcus aureus is an opportunistic pathogen that can cause fatal bacterial infections. MurD catalyzes the formation of peptide bond between UDP-N-acetylehyl-l-alanine and d-glutamic acid, which plays an important role in the synthesis of peptidoglycan and the formation of cell wall by S. aureus. Because S. aureus is resistant to most existing antibiotics, it is necessary to develop new inhibitors. In this study, Schrodinger 11.5 Prime homology modeling was selected to prepare the protein model of MurD enzyme, and its structure was optimized. We used a virtual screening program and similarity screening to screen 47163 compounds from three marine natural product libraries to explore new inhibitors of S. aureus. ADME provides analysis of the physicochemical properties of the best performing compounds during the screening process. To determine the stability of the docking effect, a 100 ns molecular dynamics was performed to verify how tightly the compound was bound to the protein. By docking analysis and molecular dynamics analysis, both 46604 and 46608 have strong interaction with the docking pocket, have good pharmacological properties, and maintain stable conformation with the target protein, so they have a chance to become drugs for S. aureus. Through virtual screening, similarity screening, ADME study and molecular dynamics simulation, 46604 and 46608 were selected as potential drug candidates for S. aureus.     

Combining molecular docking and molecular dynamics studies for modelling Staphylococcus aureus MurD inhibitory activity

The ATP-dependent bacterial MurD enzyme catalyses the formation of the peptide bond between cytoplasmic intermediate UDP-N-acetylmuramoyl-L-alanine and D-glutamic acid. This is essential for bacterial cell wall peptidoglycan synthesis in both Gram-positive and Gram-negative bacteria. MurD is recognized as an important target for the development of new antibacterial agents. In the present study we prepared the 3D-stucture of the catalytic pocket of the Staphylococcus aureus MurD enzyme by homology modelling. Extra-precision docking, binding free energy calculation by the MM–GBSA approach and a 40 ns molecular dynamics (MD) simulation of 2-thioxothiazolidin-4-one based inhibitor $1 was carried out to elucidate its inhibition potential for the S. aureus MurD enzyme. Molecular docking results showed that Lys19, Gly147, Tyr148, Lys328, Thr330 and Phe431 residues are responsible for the inhibitor–protein complex stabilization. Binding free energy calculation revealed electrostatic solvation and van der Waals energy components as major contributors for the inhibitor binding. The inhibitor-modelled S. aureus protein complex had a stable conformation in response to the atomic flexibility and interaction, when subjected to MD simulation at 40 ns in aqueous solution. We designed some molecules as potent inhibitors of S. aureus MurD, and to validate the stability of the designed molecule D1-modelled protein complex we performed a 20 ns MD simulation. Results obtained from this study can be utilized for the design of potent S. aureus MurD inhibitors.     

Extra precision docking,free energy calculation and molecular dynamics studies on glutamic acid derivatives as MurD inhibitors

The binding modes of well known MurD inhibitors have been studied using molecular docking and molecular dynamics (MD) simulations. The docking results of inhibitors 1-30 revealed similar mode of interaction with Escherichia coli-MurD. Further, residues Thr36, Arg37, His183, Lys319, Lys348, Thr321, Ser415 and Phe422 are found to be important for inhibitors and E. coli-MurD interactions. Our docking procedure precisely predicted crystallographic bound inhibitor 7 as evident from root mean square deviation (0.96 Å). In addition inhibitors 2 and 3 have been successfully cross-docked within the MurD active site, which was pre-organized for the inhibitor 7. Induced fit best docked poses of 2, 3, 7 and 15/2Y1O complexes were subjected to 10 ns MD simulations to determine the stability of the predicted binding conformations. Induce fit derived docked complexes were found to be in a state of near equilibrium as evident by the low root mean square deviations between the starting complex structure and the energy minimized final average MD complex structures. The results of molecular docking and MD simulations described in this study will be useful for the development of new MurD inhibitors with high potency.     

Potential drug targets in Mycobacterium tuberculosis through metabolic pathway analysis

The emergence of multidrug resistant varieties of Mycobacterium tuberculosis has led to a search for novel drug targets. We have performed an insilico comparative analysis of metabolic pathways of the host Homo sapiens and the pathogen M. tuberculosis. Enzymes from the biochemical pathways of M. tuberculosis from the KEGG metabolic pathway database were compared with proteins from the host H. sapiens, by performing a BLASTp search against the non-redundant database restricted to the H. sapiens subset. The e-value threshold cutoff was set to 0.005. Enzymes, which do not show similarity to any of the host proteins, below this threshold, were filtered out as potential drug targets. We have identified six pathways unique to the pathogen M. tuberculosis when compared to the host H. sapiens. Potential drug targets from these pathways could be useful for the discovery of broad spectrum drugs. Potential drug targets were also identified from pathways related to lipid metabolism, carbohydrate metabolism, amino acid metabolism, energy metabolism, vitamin and cofactor biosynthetic pathways and nucleotide metabolism. Of the 185 distinct targets identified from these pathways, many are in various stages of progress at the TB Structural Genomics Consortium. However, 67 of our targets are new and can be considered for rational drug design. As a case study, we have built a homology model of one of the potential drug targets MurD ligase using WHAT IF software. The model could be further explored for insilico docking studies with suitable inhibitors. The study was successful in listing out potential drug targets from the M. tuberculosis proteome involved in vital aspects of the pathogen's metabolism, persistence, virulence and cell wall biosynthesis. This systematic evaluation of metabolic pathways of host and pathogen through reliable and conventional bioinformatic methods can be extended to other pathogens of clinical interest.