A selective method for optimizing ensemble docking-based experiments on an InhA Fully-Flexible receptor model

Abstract

BackgroundIn the rational drug design process, an ensemble of conformations obtained from a molecular dynamics simulation plays a crucial role in docking experiments. Some studies have found that Fully-Flexible Receptor (FFR) models predict realistic binding energy accurately and improve scoring to enhance selectiveness. At the same time, methods have been proposed to reduce the high computational costs involved in considering the explicit flexibility of proteins in receptor-ligand docking. This study introduces a novel method to optimize ensemble docking-based experiments by reducing the size of an InhA FFR model at docking runtime and scaling docking workflow invocations on cloud virtual machines.ResultsFirst, in order to find the most affordable cost-benefit pool of virtual machines, we evaluated the performance of the docking workflow invocations in different configurations of Azure instances. Second, we validated the gains obtained by the proposed method based on the quality of the Reduced Fully-Flexible Receptor (RFFR) models produced using AutoDock4.2. The analyses show that the proposed method reduced the model size by approximately 50% while covering at least 86% of the best docking results from the 74 ligands tested. Third, we tested our novel method using AutoDock Vina, a different docking software, and showed the positive accuracy achieved in the resulting RFFR models. Finally, our results demonstrated that the method proposed optimized ensemble docking experiments and is applicable to different docking software. In addition, it detected new binding modes, which would be unreachable if employing only the rigid structure used to generate the InhA FFR model.ConclusionsOur results showed that the selective method is a valuable strategy for optimizing ensemble docking-based experiments using different docking software. The RFFR models produced by discarding non-promising snapshots from the original model are accurately shaped for a larger number of ligands, and the elapsed time spent in the ensemble docking experiments are considerably reduced.Electronic supplementary materialThe online version of this article (10.1186/s12859-018-2222-2) contains supplementary material, which is available to authorized users.

<p>ISL Induces Apoptosis and Autophagy in Hepatocellular Carcinoma via Downregulation of PI3K/AKT/mTOR Pathway in vivo and in vitro</p>

“… 26 Clusters with more conformations and the largest absolute value of binding energy were selected from the results. 27 , 28 The most stable binding configuration in the top-ranked cluster was selected for further analysis using Pymol software 29 ( https://pymol.org/2/ ) and MOE software.…”

Exploring the active mechanism of berberine against HCC by systematic pharmacology and experimental validation

“…Molecular docking was performed and analyzed using AutoDock 4.2. A Lamarckian genetic algorithm method (Runs 100) was implemented in the program suite to identify appropriate binding modes and conformation of the ligand molecules (30). A maximum of 25 million energy evaluations were applied for the experiment.…”

Key Topics in Molecular Docking for Drug Design

“…A strategy, usually referred to as Receptor Ensemble Docking (or simply Ensemble Docking) is the most frequently used to model those scenarios. It is based on the concept of Conformational Selection and consists in using multiple conformations of the receptor molecule, that can be obtained via different methods, such as MD simulations [66,67], Normal Mode Analysis [68], and even by using alternative experimentally-determined receptor conformations [69]. It is worth noting that some software, such as GOLD and Glide have implemented functionality to execute this type of analysis.…”