Abstract:The accurate prediction of protein-ligand binding free energies with tractable computational methods has the potential to revolutionize drug discovery. Modeling the protein-ligand interaction at a quantum mechanical level, instead of relying on empirical classical-mechanics methods, is an important step toward this goal. In this study, we explore the QM-PBSA method to calculate the free energies of binding of seven ligands to the T4-lysozyme L99A/M102Q mutant using linear-scaling density functional theory on the whole protein-ligand complex. By leveraging modern high-performance computing we perform over 2900 full-protein (2600 atoms) DFT calculations providing new insights into the convergence, precision and reproducibility of the QM-PBSA method. We find that even at moderate sampling over 50 snapshots, the convergence of QM-PBSA is similar to traditional MM-PBSA and that the DFT-based energy evaluations are very reproducible. We show that in the QM-PBSA framework, the physically-motivated GGA exchange-correlation functional PBE outperforms the more modern, dispersion-including non-local and meta-GGA-nonlocal functionals VV10 and B97M-rV. Different empirical dispersion corrections perform similarly well but the three-body dispersion term, as included in Grimme's D3 dispersion, is significant and improves results slightly. Inclusion of an entropy correction term sampled over less than 25 snapshots is detrimental while an entropy correction sampled over the same 50 or 100 snapshots as the enthalpies improves the accuracy of the QM-PBSA method. As full-protein DFT calculations can now be performed on modest computational resources our study demonstrates that they can be a useful addition to the toolbox of free energy calculations.

Using the Fast Fourier Transform in Binding Free Energy Calculations

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

Alchemical Transfer Approach to Absolute Binding Free Energy Estimation

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Ligand binding free-energy calculations with funnel metadynamics

Binding Energy Distribution Analysis Method: Hamiltonian Replica Exchange with Torsional Flattening for Binding Mode Prediction and Binding Free Energy Estimation

Equilibrium and Non-equilibrium Ensemble Methods for Accurate, Precise and Reproducible Absolute Binding Free Energy Calculations

Feature functional theory - binding predictor (FFT-BP) for the blind prediction of binding free energies

Protein-ligand free energies of binding from full-protein DFT calculations: convergence and choice of exchange-correlation functional

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model

Computing Ligand Binding Free Energy in a Large Flexible Pocket of a Large Protein

Nonbonded Force Field Parameters from MBIS Partitioning of the Molecular Electron Density Improve Binding Affinity Predictions of the T4-Lysozyme Double Mutant

Non-bonded force field parameters from MBIS partitioning of the molecular electron density improve binding affinity predictions of the T4-lysozyme double mutant

Free Energy Calculations Using the Movable Type Method with Molecular Dynamics Driven Protein-Ligand Sampling

Calculation of absolute free energy of binding for theophylline and its analogs to RNA aptamer using nonequilibrium work values

Efficient Calculation of Protein-Ligand Binding Free Energy Using GFN Methods: the Power of the Cluster Model.

Accurate calculation of absolute protein-ligand binding free energies

Accurate determination of protein:ligand standard binding free energies from molecular dynamics simulations

Fast, Accurate, and Reliable Protocols for Routine Calculations of Protein–Ligand Binding Affinities in Drug Design Projects Using AMBER GPU-TI with ff14SB/GAFF

Exploring the Impact of Protein Chain Selection in Binding Energy Calculations with DFT