Abstract

Molecular simulations have made great progresses in predicting koff values--the kinetic constant of drug unbinding, a key parameter for modern pharmacology--yet computed values under- or over-estimate experimental data in a system- and/or technique-dependent way. In an effort at gaining insights on this issue, here we used an established method to calculate koff values--frequency-adaptive metadynamics with force field-- and a subsequent QM/MM descriptions of the interactions. First, using force field-based metadynamics, we calculate koff of the Positron Emission Tomography (PET) ligand iperoxo targeting the human muscarinic acetylcholine receptor M2. In line with previously performed in silico studies, the prediction (3.7 {+/-} 0.7 * 10-4 s-1) turned out to differ significantly from the experimentally measured value (1.0 {+/-} 0.2 * 10-2 s-1). Next, we use DFT-based QM/MM simulations to show that this discrepancy arises from erroneous force field energetics at the transition state. It turns out that this discrepancy is partly caused by lack of electronic polarization and/or charge transfer in commonly employed force field. We expect these issues to arise also in other systems where charged portions of the system play a pivotal role, such as protein- or DNA-protein complexes. Graphical TOC Entry O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=200 SRC="FIGDIR/small/015396v1_ufig1.gif" ALT="Figure 1"> View larger version (32K): org.highwire.dtl.DTLVardef@1ee4e1aorg.highwire.dtl.DTLVardef@49ba0forg.highwire.dtl.DTLVardef@565f25org.highwire.dtl.DTLVardef@80940e_HPS_FORMAT_FIGEXP M_FIG C_FIG