Carboxymuconolactone decarboxylase is a prospective molecular target for multi-drug resistant Acinetobacter baumannii-computational modeling, molecular docking and dynamic simulation studies

Abstract

Multidrug-resistant Acinetobacter baumannii (MDRAb), a priority-I pathogen declared by the World Health Organization, became a potential healthcare concern worldwide with a high mortality rate. Thus, the identification of putative molecular targets and potential lead molecules is an important concern in healthcare. The present study aimed to screen a prospective molecular target and effectual binders for the drug discovery of MDRAb by computational virtual screening approach. Based on the functional role, γ-carboxymuconolactone decarboxylase (CMD) was prioritized as the target and its three-dimensional (3D) structure was computationally modeled. Based on the availability of the 3D structure, twenty-five herbal molecules were selected by database search, and their drug-likeliness, pharmacokinetic, and toxicity features were predicted. The effectual binding of the selected molecules towards CMD was predicted by molecular docking. The stability of the best-docked complexes was predicted by molecular dynamics (MD) simulation for 100 ns and binding energy calculations were carried out by molecular mechanics generalized Born and surface area solvation (MM/GBSA) method. Out of twenty-five molecules screened, hirsutine (an indole alkaloid of Uncaria rhynchophylla) and thymoquinone (a phytochemical of Nigella sativa) were qualified for drug likeliness, pharmacokinetic, and toxicity features and demonstrated significant effectual binding to CMD when compared with the binding of co-crystallized inhibitor and CMD (control). The docked complexes of hirsutine and thymoquinone, and CMD were stabilized by the binding energies of −8. 30 and −8. 46 kcal/mol respectively. These molecules were qualified in terms of ideal drug likeliness, ADME, and toxicity properties. MD simulation studies showed that the ligand-protein complexes were stable throughout the simulation. The binding free energies of the complexes by MMGBSA were estimated to be −42.08157745 kcal/mol and −36.58618242 kcal/mol for hirsutine and thymoquinone respectively when compared with the calculated binding free energy of the control (−28.75032666 kcal/mol). This study concluded that hirsutine and thymoquinone can act as potential lead molecules against CMD and the present hypothesis can be scaled up to develop potential inhibitors against MDRAb.