QM-Polarized Ligand Docking


Accurate treatment of electrostatic charges is crucial to the success of any docking algorithm. Although contemporary force fields are capable of modeling partial atomic charges on ligands with reasonable accuracy, they are generally incapable of considering charge polarization induced by the protein environment. The greater the role charge polarization plays in determining a ligand's bound conformation, the more difficult it will be for MM docking algorithms to perceive the correct binding mode. For research applications that demand the highest level of docking accuracy, Schrödinger introduces QM-Polarized Ligand Docking (QPLD), which uses ab inito charge calculations to overcome this limitation.

QPLD combines the docking power of Glide with the accuracy of QSite, Schrödinger's respected QM/MM software. The QPLD algorithm begins with a Glide docking job that generates several geometrically unique protein-ligand complexes. QSite then performs a single-point energy calculation on each complex, treating the ligand with ab initio methods and deriving partial atomic charges using electrostatic potential fitting. Glide then re-docks the ligand using each of the ligand charge sets calculated by QSite, and the QPLD algorithm returns the most energetically favorable pose. The fully automated algorithm is calibrated to provide useful default settings that can be modified at the user's discretion.

In keeping with Schrödinger's tradition of pairing innovation and practicality, QPLD calculations are effortlessly set up and launched using a single panel within the Maestro interface. Calculations are easily parallelized across multiple processors, and results are automatically incorporated into Maestro for visualization and analysis.


Rigorously tested throughout its development, QPLD has been shown to provide substantial improvements in docking accuracy over pure MM methods. Across a test set of 271 protein-ligand complexes taken from the PDB, QPLD succeeded in docking the ligand with an average RMSD of less than one angstrom (below).

In validation experiments comparing QPLD against pure force field methods, ligands were docked into their native crystal structure using both QPLD and Glide's standard precision (SP) mode. Across the 271 complex test set, QPLD offered superior docking accuracy: the average RMSD of QPLD-docked ligands was nearly half that obtained by Glide SP. The above histogram illustrates the distribution of ligand RMSDs obtained by QPLD and force field docking. QPLD returns extremely high accuracy poses (RMSD below 0.5 Å) for over 120 ligands in the test set, while reducing the overall number of low accuracy poses (RMSD above 3.0 Å) to less than 20.

"Quenching-Resistant Solid-State Photoluminescence of Graphene Quantum Dots: Reduction of π−π Stacking by Surface Functionalization with POSS, PEG, and HDA"

Park, M.; Jeong, Y.; Kim, H.S.; Lee, W.; Nam, S. H.; Lee, S.; Yoon, H.; Kim, J.; Yoo, S.; Jeon, S., , 2021, ,

"Selective targeting of dendritic cell-specific intercellular adhesion molecule-3-grabbing nonintegrin (DC-SIGN) with mannose-based glycomimetics: Synthesis and interaction studies of bis(benzylamide) derivatives of a pseudomannobioside "

Varga, N.; Sutkeviciute, I.; Guzzi, C.; McGeagh, J.; Petit-Haertlein, I.; Gugliotta, S.; Weiser, J.; Angulo, J.; Fieschi, F.; Bernardi, A., Chem. Eur. J., 2013, 19, 4786-4797

"Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and Prime/MM-GBSA calculation"

Du, J.; Sun, H.; Xi, L.; Li, J.; Yang, Y.; Liu, H.; Yao, X., J. Comp. Chem., 2011, 32(13), 2800-2809

"Correlation between Performance of QM/MM Docking and Simple Classification of Binding Sites"

Chung, J. Y.; Hah, J.; Cho, A. E., J. Chem. Inf. Model., 2009, 49, 2382–2387

ö "Extension of QM/MM docking and its applications to metalloproteins"

Cho, A. E.; Rinaldo, D., J. Comput. Chem., 2009, 30, 2609–2616

"Glucose-based spiro-isoxazolines: A new family of potent glycogen phosphorylase inhibitors"

Benltifa, M.; Hayes, J. M.; Vidal, S.; Gueyrard, D.; Goekjian, P. G.; Praly, J.; Kizilis, G.; Tiraidis, C.; Alexacou, K.; Chrysina, E. D.; Zographos, S. E.; Leonidas, D. D.; Archontis, G.; Oikonomakos, N. G., Bioorg. Med. Chem., 2009, 17, 7368-7380

"Enantioselective Proteins: Selection, Binding Studies and Molecular Modeling of Antibodies with Affinity towards Hydrophobic BINOL Derivatives"

Rasmussen, B. S.; Pedersen, J. M.; Sørensen, J.; Egebjerg, J.; Schiøtt, B.; Mortensen, K. K.; Skrydstrup, T., ChemBioChem, 2007, 8, 1974-1980

ö "Importance of Accurate Charges in Molecular Docking: Quantum Mechanical/Molecular Mechanical (QM/MM) Approach"

Cho, A. E.; Guallar, V.; Berne, B.; Friesner, R. A., J. Comput. Chem., 2005, 26, 915–931

Citations and Acknowledgements

Schrödinger Release 2021-4: QM-Polarized Ligand Docking protocol; Glide, Schrödinger, LLC, New York, NY, 2021; Jaguar, Schrödinger, LLC, New York, NY, 2021; QSite, Schrödinger, LLC, New York, NY, 2021.

Back To Top