Nahd Mohamed Elmaki1, Inass A Al Sadawi2, Anton Hermann3 and Abdul M Gbaj*2
Received: September 04, 2018; Published: September 10, 2018
Corresponding author: Abdul M Gbaj, Department of Medicinal Chemistry, University of Tripoli, Libya
Molecular modeling attempts to study the function, structure and inhibition of the acetylcholinesterase enzyme due to the fact that the inhibition of this enzyme is importance to medical conditions such as Alzheimer’s disease, myasthenia gravis and Parkinson’s disease, and it is also important in eminent toxicological susceptibility to nerve agents. In this study we present an approach for forecasting the inhibitory activity of acetylcholinesterase (AChE) inhibitors by using docking studies. The docking studies were done on acetylcholinesterase to attain the conformation of the enzyme in water surroundings. The obtained conformation of the enzyme was used for docking with four inhibitors (physostigmine, neostigmine, pyridostigmine and rivastigmine). Docking analysis showed that hydrogen bonds and hydrophobic interactions play important tasks in the acetylcholinesterase -inhibitor complex. Subsequently, all inhibitors that bind at the catalytic site of acetylcholinesterase and their interactions with acetylcholinesterase were studied. In addition, conformation stability of acetylcholinesterase -inhibitors was studied using simulation docking technique. The complex showed that acetylcholinesterase conformation did not change significantly in the presence of the four inhibitors. This paper showed important studies on acetylcholinesterase and assists to illuminate the four inhibitors interdependencies using molecular modeling.
Keywords: Acetylcholinesterase; Molecular docking; Physostigmine; Neostigmine; Pyridostigmine; Rivastigmine
Acetylcholine (ACh) is a neurotransmitter broadly distributed in the central and peripheral, autonomic nervous system (CNS). In the CNS, ACh performs several functions, such as memory, learning, motor control and attention. Acetylcholinesterase (AChE), one of the most essential enzymes in the family of serine hydrolases, catalyzes the hydrolysis of neurotransmitter acetylcholine, which plays an important role in memory and cognition [1-4]. It has been suggested that acetylcholinesterase (AChE) degrades the esters of choline and has a role in neurotransmission within the autonomic and somatic motor nervous systems and it is the target of action of the drugs (inhibitors) such as physostigmine, neostigmine, pyridostigmine and rivastigmine, physostigmine [1-4]. It is well known that AChE enzyme regulated the release and entrance of ACh in cholinergic fibers . The function of AChE in neurological disorders like Myasthenia gravis (MG) , Alzheimer’s disease , Parkinson’s disease  and other ‘non-classical’ activities such as, neurite formation, network formation and cell adhesion  draw attention of many researches related to the medical field. The exact mechanisms of the interaction of physostigmine, neostigmine, pyridostigmine and rivastigmine with the acetylcholinesterase receptor complex still need further molecular modeling studies.
Acetylcholinesterase has extremely high catalytic activity and each molecule of AChE hydrolyze about twenty-five thousand molecules of acetylcholine (ACh) per second. The active site of AChE contains two subsites which are the anionic site and the esteratic subsite [10,11]. The mechanism of action and structure of AChE have been obtained from the crystal structure of the enzyme which were studies by different research groups [12- 14]. The anionic subsite housing the positive quaternary amine of the acetylcholine in addition to other cationic inhibitors and substrates. The substrates which have cationic charge do not bound by the negatively charged amino acid in the anionic site, but the binging is performed through interaction of fourteen aromatic residues that available close to the active site [15-17]. The active site of AChE is located four angstroms from the substructure of the molecule . The acetylcholine is hydrolyzed in the esteratic subsite and produces acetate and choline. The active site contains mainly a catalytic triad of three amino acids which are: serine 203, histidine 447 and glutamate 337 . The hydrolysis reaction by AChE on carboxyl ester produces acyl-enzyme and free choline. Subsequently, the acyl-enzyme undertakes nucleophilic attack by a water molecule, supported by the histidine 447 moiety and the releasing acetic acid and regenerating the free enzyme again . The main aims of the paper are: virtual screening approach, pharmacophore modeling, molecular docking, and consensus binding affinities of acetylcholine and four well know inhibitors. The obtained results will be used as a guide to identify and design novel AChE inhibitors with higher selectivity.
The starting geometry of physostigmine, neostigmine, pyridostigmine and rivastigmine was constructed using chem3D Ultra (version 8.0, Cambridgesoft Com., USA). The optimized geometry of the four inhibitors with the lowest energy was used in the molecular docking. The crystal structures of human acetylcholinesterase in a complex with a transition-state analogue were downloaded from the Protein Data Bank (https://www. rcsb.org/structure/5hfa). The molecular dockings of the four inhibitors with human acetylcholinesterase was accomplished by AutoDock 4.2 software from the Scripps Research Institute (TSRI) (http://autodock.scripps.edu/). Firstly, the polar hydrogen atoms were added into human acetylcholinesterase and the four inhibitors molecules. Then, the partial atomic charges of the human acetylcholinesterase and the four inhibitors molecule were calculated using Kollman methods . In the process of molecular docking, the grid maps of dimensions (62Å X 62Å X 62Å) with a grid-point spacing of 0.376Å and the grid boxes is centered. The number of genetic algorithms runs, and the number of evaluations was set to 100. All other parameters were default settings. Cluster analysis was performed on the results of docking by using a root mean square (RMS) tolerance of 2.0Å, and this was dependent on the binding free energy. Lastly, the dominating configuration of the binding complex of the four inhibitors and human acetylcholinesterase with minimum energy of binding can be determined.
The modeling study was performed in this paper showed great interactions between physostigmine, neostigmine, pyridostigmine and rivastigmine and human acetylcholinesterase. The binding energies of the four inhibitors and human acetylcholinesterase were shown in Table1. The geometry of docking obtained with acetylcholine and pyridostigmine with human acetylcholinesterase as shown in (Figure 1A & B), respectively. The four inhibitors were able to form hydrogen bonds (HBs) with the amino acid residues of the enzyme, pi-pi stacking and Pi-alkyl interaction (Table 1). In addition, the molecular docking results showed that other amino acids residues are involved in the interactions with the four inhibitors. The binding energies scores as shown in Table 1 were calculated by adding up a set of weighted empirical energies including van der Waals forces (VDW), electrostatic, hydrogenbonding, desolvation, entropy, and hydrophobicity. A good number of pi‒pi interactions were observed (Table 1) that show the bond distance varying between 3.40 Å to as far as 3.95 Å which is consistent with the results obtained by Avasthi et al. whom indicated that the ideal bond distance in pi‒pi interactions are within the range 3.30–4.00 Å . In addition, in order to provide another clear reason of the activity difference in terms of dynamic behavior, distances Pi-sigma monitor interactions between Ser203 and ligands were determined and compared during the simulation time (Table 1). The distances obtained were similar to the literature where the distance of Pi-sigma interaction less than 0.5nm . This indicated that the π-sigma interaction could also be one of the key interactions to elucidate the activity dissimilarity in terms of dynamic behavior. The docking results obtained in Table 1 with the acetylcholine and the various inhibitors provided insights into the essential structural elements and motifs central of the catalytic mechanism of the acetylcholinesterase. It has been reported that the X-ray analysis of the structural features of AChE has a narrow, long, hydrophobic gorge which is about 20 Å deep [11,12]. The AChE has a catalytic triad consisting of His447, Ser203, and Glu334  situated in the active site of the narrow deep gorge, the lining of which contains mostly aromatic residues that form a narrow access to the catalytic Ser203 . In addition to that a peripheral anionic site containing an additional set of aromatic residues Trp286, Tyr124, Tyr72, Tyr341, and Asp74  was reported which is located at the border of the gorge and affords a binding site for allosteric inhibitors and modulators. The binding between the acetylcholine and the four inhibitors with the enzyme is characterized by cation-π interactions between the protonated nitrogen’s and the preserved aromatic residues, phenylalanine and tryptophan. Furthermore, π-π stacking between the aromatic moieties of the inhibitors and the aromatic amino acids as shown in Table 1 has critical roles in binding and this consistent with the literature as described by Shin Hua et al . With the current obtained data, computational methods have become essential to biological investigations. Here we have used computational approach to further understand the mechanism of interactions and binding affinity between AChE with drug molecules. The present analysis permits us to draw the number of conclusions. The computational methods on the catalytic site help the researcher to design new drug molecules. The molecular docking programs are helpful in understanding the interaction between the AChE with diverse lead/drug molecules. Our analysis also demonstrates that the four inhibitors could be the prospective lead molecule for the inhibition of AChE. Hence docking results of acetylcholine and the four inhibitors could be used as the template for designing therapeutic lead molecule. We hope that the originality and success of the computational efforts in this paper offer fruitful promise for the future predictions of finding new inhibitors which could results into enormous reductions in therapeutics development time.
Table 1: Binding and intermolecular energies from molecular docking for compounds acetylcholine, physostigmine, neostigmine, pyridostigmine and rivastigmine.
Figure 1: Three-dimensional representation of acetylcholine (A) and pyridostigmine (B) interacting with active site of target macromolecule human acetylcholinesterase.
ILs are receiving more and more attention every day both in academic research and commercial applications and they seem as good replacements for volatile organic solvents. However, there is a discussion about the greenness of the ILs due to their incomplete physical, chemical and toxicological data. most of the researchers will continue to work with this new solvent, the developments of new applications utilizing ILs will increase rapidly and the number of publications will rise exponentially in the future. The present costs of the ILs are quite prohibitive in many probable commercial applications. However, there are hopes that in the near future, the cost/benefit figures of the ILs will bring economic viability to their more common use. The development of pharmaceutical applications using IL-based methodologies requires a deep understanding of ILs both in terms of their macroscopic properties and also at the molecular level, because structural aspects have been shown to play a crucial and unexpected role in a large number of situations. Despite the fact that ILs can no longer be considered a new field, the large diversity of combinations of cations and anions producing novel ILs with new specific properties is astonishing.
Bio chemistryUniversity of Texas Medical Branch, USA
Department of Criminal JusticeLiberty University, USA
Department of PsychiatryUniversity of Kentucky, USA
Department of MedicineGally International Biomedical Research & Consulting LLC, USA
Department of Urbanisation and AgriculturalMontreal university, USA
Oral & Maxillofacial PathologyNew York University, USA
Gastroenterology and HepatologyUniversity of Alabama, UK
Department of MedicineUniversities of Bradford, UK
OncologyCirculogene Theranostics, England
Radiation ChemistryNational University of Mexico, USA
Analytical ChemistryWentworth Institute of Technology, USA
Minimally Invasive SurgeryMercer University school of Medicine, USA
Pediatric DentistryUniversity of Athens , Greece
The annual scholar awards from Lupine Publishers honor a selected number Read More...