Virtual Screening, Molecular Docking and Dynamic Simulation of Shikimate Kinase from Mycobacterium Tuberculosis Using in Silico Approach
Preprints are manuscripts made publicly available before they have been submitted for formal peer review and publication. They might contain new research findings or data. Preprints can be a draft or final version of an author's research but must not have been accepted for publication at the time of submission.
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.