The global antibiotic resistance problem necessitates fast and effective approaches to finding novel inhibitors to treat bacterial infections. In this study, we propose a computational workflow to identify plausible high-affinity compounds from FDA-approved, investigational, and experimental libraries for the decoding center on the small subunit 30S of the E. coli ribosome. The workflow basically consists of two molecular docking calculations on the intact 30S, followed by molecular dynamics (MD) simulations coupled with MM-GBSA calculations on a truncated ribosome structure. The parameters used in the molecular docking suits, Glide and AutoDock Vina, as well as in the MD simulations with Desmond were carefully adjusted to obtain expected interactions for the ligand-rRNA complexes. A filtering procedure was followed, considering a fingerprint based on aminoglycoside's binding site on the 30S to obtain seven hit compounds either with different clinical usages or aminoglycoside derivatives under investigation, suggested for in vitro studies. The detailed workflow developed in this study promises an effective and fast approach for the estimation of binding free energies of large protein-RNA and ligand complexes.
Keywords: 16S rRNA; Antibiotic resistance; Binding free energy; Drug repurposing; Molecular docking; Molecular dynamics simulations.
Copyright © 2024 Elsevier Inc. All rights reserved.