Abstract:The opportunity to prospectively predict ligand bound poses and free energies of binding to the Farnesoid X Receptor in the D3R Grand Challenge 2 provided a useful exercise to evaluate CHARMM based docking (CDOCKER) and \(\lambda\)-dynamics methodologies for use in “real-world” applications in computer aided drug design. In addition to measuring their current performance, several recent methodological developments have been analyzed retrospectively to highlight best procedural practices in future applications. For pose prediction with CDOCKER, when the protein structure used for rigid receptor docking was close to the crystallographic holo structure, reliable poses were obtained. Benzimidazoles, with a known holo receptor structure, were successfully docked with an average RMSD of 0.97 \(\AA\). Other non-benzimidazole ligands displayed less accuracy largely because the receptor structures we chose for docking were too different from the experimental holo structures. However, retrospective analysis has shown that when these ligands were re-docked into their holo structures, the average RMSD dropped to 1.18 \(\AA\) for all ligands. When sulfonamides and spiros were docked with the apo structure, which agrees more with their holo structure than the structures we chose, five out of six ligands were correctly docked. These docking results emphasize the need for flexible receptor docking approaches. For \(\lambda\)-dynamics techniques, including multisite \(\lambda\)-dynamics (MS\(\lambda\)D), reasonable agreement with experiment was observed for the 33 ligands investigated; root mean square errors of 2.08 and 1.67 kcal/mol were obtained for free energy sets 1 and 2, respectively. Retrospectively, soft-core potentials, adaptive landscape flattening, and biasing potential replica exchange (BP-REX) algorithms were critical to model large substituent perturbations with sufficient precision and within restrictive timeframes, such as was required with participation in Grand Challenge 2. These developments, their associated benefits, and proposed procedures for their use in future applications are discussed.

Prediction of binding poses to FXR using multi-targeted docking combined with molecular dynamics and enhanced sampling

Binding pose and affinity prediction in the 2016 D3R Grand Challenge 2 using the Wilma-SIE method

Prediction of ligand binding mode among multiple cross-docking poses by molecular dynamics simulations

CDOCKER and $$\Lambda$$-Dynamics for Prospective Prediction in D3R Grand Challenge 2

Structure-based pose prediction: Non-cognate docking extended to macrocyclic ligands

Biologically enhanced sampling geometric docking and backbone flexibility treatment with multiconformational superposition

Predicting the bioactive conformations of macrocycles: a molecular dynamics-based docking procedure with DynaDock

DOX: A New Computational Protocol for Accurate Prediction of the Protein-Ligand Binding Structures.

Synergistic Application of Molecular Docking and Machine Learning for Improved Binding Pose

Uni-Mol Docking V2: Towards Realistic and Accurate Binding Pose Prediction

The Application of the MM/GBSA Method in the Binding Pose Prediction of FGFR Inhibitors

Performance of HADDOCK and a simple contact-based protein-ligand binding affinity predictor in the D3R Grand Challenge 2

DockFormer: Efficient Multi-Modal Receptor-Ligand Interaction Prediction using Pair Transformer

A Robust Induced Fit Docking Approach with the Combination of the Hybrid All-Atom/United-Atom/Coarse-Grained Model and Simulated Annealing

Docking Challenge: Protein Sampling and Molecular Docking Performance

Quantum Mechanics-Based Fast and Reliable Prediction of Binding Pose Structures

DiffBindFR: An SE(3) Equivariant Network for Flexible Protein-Ligand Docking

Docking guidance with experimental ligand structural density improves docking pose prediction and virtual screening performance

On Evaluating Molecular-Docking Methods for Pose Prediction and Enrichment Factors

Fast docking using the CHARMM force field with EADock DSS