Target–ligand binding affinity from single point enthalpy calculation and elemental composition

Reliable target–ligand binding thermodynamics data are essential for successful drug design and molecular engineering projects. Besides experimental methods, a number of theoretical approaches have been introduced for the generation of binding thermodynamics data. However, available approaches often neglect electronic effects or explicit water molecules influencing target–ligand interactions. To handle electronic effects within a reasonable time frame, we introduce a fast calculator QMH-L using a single target–ligand complex structure pre-optimized at the molecular mechanics level. QMH-L is composed of the semi-empirical quantum mechanics calculation of binding enthalpy with predicted explicit water molecules at the complex interface, and a simple descriptor based on the elemental composition of the ligand. QMH-L estimates the target–ligand binding free energy with a root mean square error (RMSE) of 0.94 kcal mol−1. The calculations also provide binding enthalpy values and they were compared with experimental binding thermodynamics data collected from the most reliable isothermal titration calorimetry studies of systems including various protein targets and challenging, large peptide ligands with a molecular weight of up to 2–3 thousand. The single point enthalpy calculations of QMH-L require modest computational resources and are based on short runs with open source and/or free software like Gromacs, Mopac, MobyWat, and Fragmenter. QMH-L can be applied for fast, automated scoring of drug candidates during a virtual screen, enthalpic engineering of new ligands or thermodynamic explanation of complex interactions.