Binding interaction of protoberberine alkaloids against acetylcholinesterase (AChE) using molecular dynamics simulations and QM/MM calculations

Acetylcholinesterase (AChE) plays a vital role in Alzheimer's disease (AD), which is one of the most common causes of dementia. Discovering new effective inhibitors against AChE activity is seen to be one of the effective approaches to reduce the suffering from AD. Protoberberine alkaloids isolated from natural resources have previously been reported as potent AChE inhibitors. In order to gain insights into how these alkaloids could inhibit AChE, berberine, palmatine, and cyclanoline were selected to investigate in terms of binding orientation and their key interactions with AChE using molecular docking and molecular dynamics simulations and quantum chemical calculations. The results revealed that the molecular dynamics structures of palmatine and berberine indicated that their equilibrated structures did not occupy the gorge but they slightly moved away from the catalytic site (CAS). For cyclanoline, the binding mode was quite different from those of donepezil and the other protoberberine alkaloids: it preferred to stay deeper in the CAS site. Interaction energies and residual interaction energies confirmed that the key interactions for palmatine and berberine were π–π interactions with Trp286 and Tyr341 and H-bond interactions with Tyr124. Cyclanoline formed π–π interactions with Trp86 and H-bonds to the amino acids in the CAS site. The results suggested the importance of aromaticity in the core structure and the flexibility of the core structure or the substituents in order to fit into the narrow gorge. The HOMO, LUMO, bioavailability, drug-likeness and pharmacokinetics were also predicted. The results obtained will be useful for further AD drug development.