In the present study, the binding free energy of some classical inhibitors (DMT, DNP, GNT, HUP, THA) with acetylcholinesterase (AChE) is calculated by means of the free energy perturbation (FEP) method based on hybrid quantum mechanics and molecular mechanics (QM/MM) potentials. The results highlight the key role of the van der Waals interaction for the inhibition process, since the contribution of this term to the binding free energy is almost as decisive as the electrostatic one. The analysis of the geometrical parameters and the interaction energy per residue along the QM/MM molecular dynamics (MD) simulations highlights the most relevant interactions in the different AChE-ligand systems, showing that the charged residues with a more prominent contribution to the interaction energy are Asp72 and Glu199, although the relative importance depends on the molecular size of the ligand. A correlation between the binding free energy and the number of cation-π interactions present in the systems has been established, DMT being the most potent inhibitor, capable of forming four cation-π interactions. A layer of water molecules surrounding the inhibitors has been observed, which act as bridges along a network formed by the ligands and the residues of the gorge and also between different residues. Although several hydrogen bonds between ligands and AChE do appear, no significant values of BIEs have been recorded. This behavior can be accounted for by the special features of AChE, such as the presence of several subsites of different natures in the gorge or the existence of several water molecules that act as bridges in the electrostatic interactions.

Hartree-Fock and density functional methods were used to analyze electronic and structural properties of known drugs to evaluate the influence of these data on acetylcholinesterase inhibition. The energies of the frontier orbitals and the distances between the more acidic hydrogen species were investigated to determine their contributions to the activity of a group of acetylcholinesterase inhibitors. Electrostatic potential maps indicated suitable sites for drugs-enzyme interactions. In this study, the structural, electronic and spatial properties of nine drugs with known inhibitory effects on acetylcholinesterase were examined. The data were obtained based on calculations at the B3LYP/6-31 + G(d,p) level. Multivariate principal components analysis was applied to 18 parameters to determine the pharmacophoric profile of acetylcholinesterase inhibitors. Desirable features for acetylcholinesterase inhibitor molecules include aromatic systems or groups that simulate the surface electrostatic potential of aromatic systems and the presence of a sufficient number of hydrogen acceptors and few hydrogen donors. PCA showed that electronic properties, including the HOMO-1 orbital energy, logP and aromatic system quantity, as well as structural data, such as volume, size and H-H distance, are the most significant properties.

No abstract

No abstract