These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.

Virtual Screening, Molecular Docking and Dynamic Simulation of Shikimate Kinase from Mycobacterium Tuberculosis Using in Silico Approach

submitted on 15.06.2020 and posted on 17.06.2020 by Mustafa Alhaji Isa

Shikimate kinase (SK) is an enzyme that catalyzes the fifth steps in the shikimate pathway. The enzyme facilitate the transfer of phosphoryl from ATP to shikimate, to produce ADP and shikimate-3-phosphate from Mycobacterium tuberculosis (MTB). The 3D structure of SK bound ligands (4-(2-Hydroxyethyl)-1-Piperazine Ethanesulfonic Acid (EPE)), ADP and metals (Mg2+, Cl- and Pt+) obtained from PDB (PDB ID: 1L4U and resolution 1.8Å). The structural analysis of the SK revealed that it has a substrate or shikimate binding site (Asp34, Arg58, and Lys136) and substrate binding via amide nitrogen (Gly80). It also possessed nucleotide binding region (Gly12─Thr17), the ATP binding site (Arg117 and Arg153) and metallic ion (Mg2+) binding site (Ser16 and Asp32). All these residues mentioned above play an essential role in the catalytic activity of the SK. Therefore inhibition any of these residues serve as a stumbling block for the normal function of the enzyme. A total of eleven thousand three hundred and twenty-three (11323) compounds obtained from two public databases (Zinc Database and PubChem) capable of binding to SK with good binding affinities. These compounds further filtered for Lipinski’s rule of five, drug-likeness, molecular docking analysis, and ADME and toxicity analysis. Three compounds with minimum binding energies─ PubChem15478 (─11.75 kcal/mol), ZINC02838601 (─11.52 kcal/mol), and ZINC11790367 (─9.88 kcal/mol) ─were selected and used for the MD simulation analysis. Also, MD simulation of the SK bound to EPE, ADP, and Mg2+ were carried out to compare their stabilities with the selected protein-ligand complexes. The result showed that the two compounds (ZINC11790367 and PubChem15478) formed stable and rigid complexes comparable to the bound ligand and the cofactors during the 50ns MD simulation. Therefore, it concluded that the above mentioned two compounds capable of inhibiting SK considered as prospective drugs for MTB after successful experimental validation.


Email Address of Submitting Author


University of Maiduguri



ORCID For Submitting Author

Declaration of Conflict of Interest

I declare I have no conflict of interest

Version Notes

I have submitted