Searching for bioactive conformations of drug-like ligands with current force fields: how good are we?

Drug-like ligands obtained from protein–ligand complexes deposited in the Protein Databank were subjected to conformational searching using various force fields and solvation settings. For each ligand, the resulting conformer pool was examined for the presence of the bioactive (crystal pose-like) conformation. Similarity of conformers toward the crystal-pose was quantified as the best achievable root mean squared deviation (RMSD, heavy atoms only). Analyzing the conformer pools generated by various force fields revealed only small differences in the likelihood of finding a crystal pose-like conformation. However, employing different solvents in the conformational search was found to be very important for achieving RMSDs below 1.0 Å. The best statistical values of likelihood were observed with a recently released force field covering a large portion of dihedral angles occurring in drug-like compounds in combination with the water as solvent. In order to enable computational chemists and modelers to efficiently use available software tools, we have additionally performed several focused analyses on ligands, grouped according to descriptors most relevant for the rational drug design.