Abstract:During the last decades, the application of hybrid quantum mechanical/molecular mechanical (QM/MM) methods has been extended to the field of drug design. In principle, the approximate QM/MM approach offers a more complete description of drug–receptor non-covalent interactions. This is especially true when charge or proton transfer, chelation of metal ions or strong polarization of ligand and protein or surface chemical groups are involved. The aim of this work was to assess the accuracy of calculated non-covalent ligand–protein interaction energies (∆EintQM/MM$$ \Delta {\mathrm{E}}_{\mathrm{int}}^{\mathrm{QM}/\mathrm{MM}} $$) obtained by the hybrid QM/MM approach employed in QSite/Jaguar of Schrödinger’s Small-Molecule Drug Discovery Suite on a set of small-molecule model systems when compared to rigorous QM calculations. The QM/MM approach was used at the density functional theory (DFT) level of theory with 6-31G* basis set, hybrid B3LYP functional and OPLS-2005 force field (DFT-B3LYP/6-31G*//OPLS-2005), a popular combination frequently used in studies on larger and complex biological systems such as drug–receptor complexes. In this work, we did not attempt to compute the most precise interaction energies of the model systems. We rather tried to assess the performance of the approximate QM/MM vs. full QM approach at the same computationally accessible level. For effective use of the QM/MM approach it is essential to select an appropriate QM region of the studied systems. To aid the selection of specific protein residues or functional groups to be included in the QM region, we evaluated the effect of its size, composition and symmetry on the accuracy of the QM/MM calculated ∆EintQM/MM$$ \Delta {\mathrm{E}}_{\mathrm{int}}^{\mathrm{QM}/\mathrm{MM}} $$. This was performed by means of a set of model clusters with well-defined configurations, which mimic the basic types of non-covalent interactions in proteins. Based on these systematic quantitative comparisons, recommendations for the addition of chemical groups or protein residues into the QM region are proposed for the popular DFT-B3LYP/6-31G*//OPLS-2005 QM/MM approach, leading to a more realistic description of ligand–protein interactions. These guidelines can have a significant bearing on computational drug or material research employing hybrid QM/MM methods by providing an estimate of the accuracy that can be expected from QM/MM studies.Graphical abstractAn approximate hybrid QM/MM approach at the DFT-B3LYP/6-31G*//OPLS-2005 level was systematically assessed for the accuracy of description of non-covalent interactions on a series of six small-molecule model systems using the QSite and Jaguar modules of Schrödinger. Guidelines for rational selection of receptor residues or function groups to be included in the QM region around the ligand were proposed based on the achieved accuracy of computed ligand–protein interaction energies obtained by QM/MM vs. the full QM approach.

A Posteriori Error Estimates for Adaptive QM/MM Coupling Methods.

A Posteriori Error Estimate and Adaptivity for QM/MM Models of Crystalline Defects

Adaptive Multiscale Coupling Methods of Molecular Mechanics based on a Unified Framework of a Posteriori Error Estimates

A Posteriori Error Estimation and Adaptive Algorithm for the Atomistic/Continuum Coupling in 2D

Error Analysis and Applications of A General Qm/Mm Approach

A Posteriori Error Estimate and Adaptive Mesh Refinement Algorithm for Atomistic/Continuum Coupling with Finite Range Interactions in Two Dimensions

The adaptive buffered force QM/MM method in the CP2K and AMBER software packages

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

Adaptive Finite Element Approximations for Kohn-Sham Models

A posteriori error estimates in quantities of interest for the finite element heterogeneous multiscale method

A-posteriori QMC-FEM error estimation for Bayesian inversion and optimal control with entropic risk measure

A Posteriori Error Estimator for Adaptive Local Basis Functions to Solve Kohn-Sham Density Functional Theory

Reshaping the QM Region On-the-Fly: Adaptive-Shape QM/MM Dynamic Simulations of a Hydrated Proton in Bulk Water

Enhancing Biomechanical Simulations Based on A Posteriori Error Estimates: The Potential of Dual Weighted Residual-Driven Adaptive Mesh Refinement

How accurate is the description of ligand–protein interactions by a hybrid QM/MM approach?

An improved QM/MM approach for metals

Analysis of Blended Atomistic/Continuum Hybrid Methods

Residual-Based A Posteriori Error Estimates for Symmetric Conforming Mixed Finite Elements for Linear Elasticity Problems

Recent Development In Quantum Mechanics/Molecular Mechanics Modeling For Materials

A posteriori error estimates for a dual finite element method for singularly perturbed reaction–diffusion problems

Mixed finite element methods for nonlinear equations : a priori and a posteriori error estimates