Ribonucleic acids (RNAs), particularly the noncoding RNAs, play key roles in cancer, making them attractive drug targets. While conventional methods such as high throughput screening are resource-intensive, computational methods such as RNA-ligand docking can be used as an alternative. However, currently available docking methods are fine-tuned to perform protein-ligand and protein-protein docking. In this work, we evaluated three commonly used docking methods─AutoDock Vina, HADDOCK, and HDOCK─alongside RLDOCK, which is specifically designed for RNA-ligand docking. Our evaluation was based on several criteria including cognate docking, blind docking, scoring potential, and ranking potential. In cognate docking, only RLDOCK showed a success rate of 70% for the top-scoring docked pose. Despite this, all four docking methods did not achieve an overall success rate exceeding 50% amidst our attempt to refine the top-scoring docked poses using molecular dynamics simulations. Meanwhile, all four docking methods showed poor performance in scoring potential evaluation. Although AutoDock Vina achieved an area under the receiver operating characteristic curve of 0.70, it showed poor performance in terms of Matthews' correlation coefficient, precision, enrichment factors, and normalized enrichment factors at 1, 2, and 5%. These results highlight the growing need for further optimization of docking methods to assess RNA-ligand interactions.
The rampant spread of multidrug-resistant Pseudomonas aeruginosa strains severely threatens global health. This severity is compounded against the backdrop of a stagnating antibiotics development pipeline. Moreover, with many promising therapeutics falling short of expectations in clinical trials, targeting the las quorum sensing (QS) system remains an attractive therapeutic strategy to combat P. aeruginosa infection. Thus, our primary goal was to develop a drug prediction algorithm using machine learning to identify potent LasR inhibitors. In this work, we demonstrated using a Multilayer Perceptron (MLP) algorithm boosted with AdaBoostM1 to discriminate between active and inactive LasR inhibitors. The optimal model performance was evaluated using 5-fold cross-validation and test sets. Our best model achieved a 90.7% accuracy in distinguishing active from inactive LasR inhibitors, an area under the Receiver Operating Characteristic Curve value of 0.95, and a Matthews correlation coefficient value of 0.81 when evaluated using test sets. Subsequently, we deployed the model against the Enamine database. The top-ranked compounds were further evaluated for their target engagement activity using molecular docking studies, Molecular Dynamics simulations, MM-GBSA analysis, and Free Energy Landscape analysis. Our data indicate that several of our chosen top hits showed better ligand-binding affinities than naringenin, a competitive LasR inhibitor. Among the six top hits, five of these compounds were predicted to be LasR inhibitors that could be used to treat P. aeruginosa-associated infections. To our knowledge, this study provides the first assessment of using an MLP-based QSAR model for discovering potent LasR inhibitors to attenuate P. aeruginosa infections.