Abstract:One approach to incorporate protein flexibility in molecular docking is the use of an ensemble consisting of multiple protein structures. Sequentially docking each ligand into a large number of protein structures is computationally too expensive to allow large-scale database screening. It is challenging to achieve a good balance between docking accuracy and computational efficiency. In this work, we have developed a fast, novel docking algorithm utilizing multiple protein structures, referred to as ensemble docking, to account for protein structural variations. The algorithm can simultaneously dock a ligand into an ensemble of protein structures and automatically select an optimal protein structure that best fits the ligand by optimizing both ligand coordinates and the conformational variable m, where m represents the m-th structure in the protein ensemble. The docking algorithm was validated on 10 protein ensembles containing 105 crystal structures and 87 ligands in terms of binding mode and energy score predictions. A success rate of 93% was obtained with the criterion of root-mean-square deviation < 2.5 angstrom if the top five orientations for each ligand were considered, comparable to that of sequential docking in which scores for individual docking are merged into one list by re-ranking, and significantly better than that of single rigid-receptor docking (75% on average). Similar trends were also observed in binding score predictions and enrichment tests of virtual database screening. The ensemble docking algorithm is computationally efficient, with a computational time comparable to that for docking a ligand into a single protein structure. In contrast, the computational time for the sequential docking method increases linearly with the number of protein structures in the ensemble. The algorithm was further evaluated using a more realistic ensemble in which the corresponding bound protein structures of inhibitors were excluded. The results show that ensemble docking successfully predicts the binding modes of the inhibitors, and discriminates the inhibitors from a set of noninhibitors with similar chemical properties. Although multiple experimental structures were used in the present work, our algorithm can be easily applied to multiple protein conformations generated by computational methods, and helps improve the efficiency of other existing multiple protein structure(MPS)-based methods to accommodate protein flexibility.

Protein-Ligand conformation parallel search method based on pose sharing

PaFlexPepDock: Parallel Ab-Initio Docking of Peptides Onto Their Receptors with Full Flexibility Based on Rosetta.

Enhancing Ligand Pose Sampling for Molecular Docking

Multiple ligand simultaneous docking: Orchestrated dancing of ligands in binding sites of protein

S4MPLE—Sampler for Multiple Protein-Ligand Entities: Methodology and Rigid-Site Docking Benchmarking

A filter enhanced sampling and combinatorial scoring study for protein docking in CAPRI

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

Assessing small molecule conformational sampling methods in molecular docking

Ensemble Docking of Multiple Protein Structures: Considering Protein Structural Variations in Molecular Docking

Parallel multi-swarm cooperative particle swarm optimization for protein–ligand docking and virtual screening

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

Advancing membrane-associated protein docking with improved sampling and scoring in Rosetta

Rosettadock In Capri Rounds 6-12

Docking Challenge: Protein Sampling and Molecular Docking Performance

EDock: blind protein-ligand docking by replica-exchange monte carlo simulation

XDock: A General Docking Method for Modeling Protein-Ligand and Nucleic Acid-Ligand Interactions.

FitDock: protein–ligand docking by template fitting

Prl-Dock: Protein-Ligand Docking Based on Hydrogen Bond Matching and Probabilistic Relaxation Labeling

SPOT-Ligand: Fast and effective structure-based virtual screening by binding homology search according to ligand and receptor similarity.