CIFDock: A novel CHARMM‐based flexible receptor–flexible ligand docking protocol

Docking studies play a critical role in the current workflow of drug discovery. However, limitations may often arise through factors including inadequate ligand sampling, a lack of protein flexibility, scoring function inadequacies (e.g., due to metals, co‐factors, etc.), and difficulty in retaining explicit water molecules. Herein, we present a novel CHARMM‐based induced fit docking (CIFDock) workflow that can circumvent these limitations by employing all‐atom force fields coupled to enhanced sampling molecular dynamics procedures. Self‐guided Langevin dynamics simulations are used to effectively sample relevant ligand conformations, side chain orientations, crystal water positions, and active site residue motion. Protein flexibility is further enhanced by dynamic sampling of side chain orientations using an expandable rotamer library. Steps in the procedure consisting of fixing individual components (e.g., the ligand) while sampling the other components (e.g., the residues in the active site of the protein) allow for the complex to adapt to conformational changes. Ultimately, all components of the complex—the protein, ligand, and waters—are sampled simultaneously and unrestrained with SGLD to capture any induced fit effects. This modular flexible docking procedure is automated using CHARMM scripting, interfaced with SLURM array processing, and parallelized to use the desired number of processors. We validated the CIFDock procedure by performing cross‐docking studies using a data set comprised of 21 pharmaceutically relevant proteins. Five variants of the CHARMM‐based SWISSDOCK scoring functions were created to quantify the results of the final generated poses. Results obtained were comparable to, or in some cases improved upon, commercial docking program data. This work presents a novel CHARMM‐based flexible receptor–flexible ligand docking (CIFDock) approach that has both an accuracy and performance on par with the current commercial induced fit docking programs. CIFDock can be used as an effective tool for drug design through its ability to customize the current protocol with added functionalities present in CHARMM and to predict final ligand poses with the latest CHARMM‐based consensus scoring. Its availability as a free (including source) program to all academic institutions has the potential to transform how computational drug discovery is carried out in academic settings and introduced to students.