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.

Molecular docking with conformer-dependent charges

Combining Classical Molecular Docking with Self-Consistent Charge Density-Functional Tight-Binding Computations for the Efficient and Quality Prediction of Ligand Binding Structure

Assessing small molecule conformational sampling methods in molecular docking

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

Conformation search across multiple-level potential-energy- surfaces (CSAMP): A strategy for accurate prediction of protein-ligand binding structures.

Improved Modeling of Peptide-Protein Binding Through Global Docking and Accelerated Molecular Dynamics Simulations

ApoDock: Ligand-Conditioned Sidechain Packing for Flexible Molecular Docking

CovalentDock: Automated Covalent Docking with Parameterized Covalent Linkage Energy Estimation and Molecular Geometry Constraints.

PackDock: a Diffusion Based Side Chain Packing Model for Flexible Protein-Ligand Docking

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

DeltaDock: A Unified Framework for Accurate, Efficient, and Physically Reliable Molecular Docking

Assessing the Performance of the Molecular Mechanics/poisson Boltzmann Surface Area and Molecular Mechanics/generalized Born Surface Area Methods. II. the Accuracy of Ranking Poses Generated from Docking

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

Dockformer: A transformer-based molecular docking paradigm for large-scale virtual screening

Toward On-The-Fly Quantum Mechanical/Molecular Mechanical (QM/MM) Docking: Development and Benchmark of a Scoring Function

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

Progress in molecular docking

Molecular Dynamics Scoring of Protein–Peptide Models Derived from Coarse-Grained Docking

How to improve docking accuracy of AutoDock4.2: a case study using different electrostatic potentials.

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

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