Experimental and Computational Investigation of Tacrine-Based Inhibitors of Acetylcholinesterase

Acetylcholinesterase (AChE) terminates cholinergic neurotransmission by catalyzing the hydrolysis of the neurotransmitter acetylcholine (ACh). Inhibition of AChE has proven an effective treatment for the memory loss exhibited by early stage Alzheimer's disease (AD) patients; four AChE inhibitors (AChEI) have been approved by the FDA for this purpose. The first AChEI approved for the palliative treatment of AD-related memory loss was 9-amino-1,2,3,4-tetrahydroacridine (tacrine).

Inhibition of AChE may present either therapeutic or toxic effects depending upon the dose administered. With the goal of discovering safe and effective pesticides to control the population of Anopheles gambiae, a malaria-transmitting mosquito indigenous to Sub-Saharan Africa, the reoptimization of the tacrine pharmacophore was undertaken. Because the optimized drug would necessarily be a poor inhibitor for human AChE (hAChE), initial ligand design focused on modification to tacrine known to negatively impact the inhibition potency for hAChE. Ultimately, an AChEI was discovered, which exhibited micromolar inhibition of Anopheles gambiae AChE (AgAChE) and essentially no potency for hAChE. Two units of this lead compound were tethered through an alkyl chain to yield a nanomolar inhibitor of AgAChE that was more than 1,100-fold selective for the mosquito enzyme over hAChE.

Dimerization of an active inhibitor is an effective strategy to increase the potency and selectivity of AChEI, and many examples of tacrine hetero- and homodimers complexed to AChE can be found in the RCSB Protein Data Bank (PDB). The bond formed between the exocyclic amine moiety and the heterocyclic ring system of tacrine is analogous to an amide bond when tacrine is protonated. Therefore, the rotational profile of protonated N-alkyltacrine should exhibit a conformational profile in which dihedral angles significantly out of the plane formed by the ring system are associated with high energies relative to those when the dihedral angles are nearly coplanar with the ring system. The barrier of rotation (ΔG^‡) produced by this phenomenon in two tacrine derivatives and two quinoline derivatives was experimentally determined using dynamic 1H NMR. These values were compared to density functional theory (DFT) derived values for the same phenomenon. Furthermore, since the ΔG^‡ proved to be impossible to experimentally determine for the optimal model compound for the active site portion of tacrine dimers, N-methyltacrine, the DFT method employed for modeling the ΔG^‡ of the tacrine and quinoline analogs was used to computationally derive the entire rotational conformation diagram of N-methyltacrine. The calculated values were then used to comment on the relative energies of adopting certain conformations found in the X-ray crystal structures of dimer/AChE complexes.