An in silico hierarchal approach for drug candidate mining and validation of natural product inhibitors against pyrimidine biosynthesis enzyme in the antibiotic-resistant Shigella flexneri

Shigella flexneri is the main causative agent of the communicable diarrheal disease, shigellosis. It is estimated that about 80-165 million cases and > 1 million deaths occur every year due to this disease. S. flexneri causes dysentery mostly in young children, elderly and immunocompromised patients, all over the globe. Recently, due to the emergence of S. flexneri antibiotic resistance strains, it is a dire need to predict novel therapeutic drug targets in the bacterium and screen natural products against it, which could eliminate the curse of antibiotic resistance. Therefore, in current study, available antibiotic-resistant genomes (n = 179) of S. flexneri were downloaded from PATRIC database and a pan-genome and resistome analysis was conducted. Around 5059 genes made up the accessory, 2469 genes made up the core, and 1558 genes made up the unique genome fraction, with 44, 34, and 13 antibiotic-resistant genes in each fraction, respectively. Core genome fraction (27% of the pan-genome), which was common to all strains, was used for subtractive genomics and resulted in 384 non-homologous, and 85 druggable targets. Dihydroorotase was chosen for further analysis and docked with natural product libraries (Ayurvedic and Streptomycin compounds), while the control was orotic acid or vitamin B13 (which is a natural binder of this protein). Dynamics simulation of 50 ns was carried out to validate findings for top-scored inhibitors. The current study proposed dihydroorotase as a significant drug target in S. flexneri and 4-tritriacontanone & patupilone compounds as potent drugs against shigellosis. Further experiments are required to ascertain validity of our findings.