Molecular recognition events in biological systems are driven by non-covalent interactions between interacting species. Here, we have studied hydrogen bonds of the CH⋅⋅⋅Y type involving electron-deficient CH donors using dispersion-corrected density functional theory (DFT) calculations applied to acetylcholinesterase-ligand complexes. The strengths of CH⋅⋅⋅Y interactions activated by a proximal cation were considerably strong; comparable to or greater than those of classical hydrogen bonds. Significant differences in the energetic components compared to classical hydrogen bonds and non-activated CH⋅⋅⋅Y interactions were observed. Comparison between DFT and molecular mechanics calculations showed that common force fields could not reproduce the interaction energy values of the studied hydrogen bonds. The presented results highlight the importance of considering CH⋅⋅⋅Y interactions when analysing protein-ligand complexes, call for a review of current force fields, and opens up possibilities for the development of improved design tools for drug discovery.
Keywords: acetylcholinesterase; density functional calculations; drug design; hydrogen bonds; quantum chemistry.
© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.