First principles calculations of thermodynamics and kinetic parameters and molecular dynamics simulations of acetylcholinesterase reactivators: can mouse data provide new insights into humans?

Review about Structure and Evaluation of Reactivators of Acetylcholinesterase Inhibited with Neurotoxic Organophosphorus Compounds

BACKGROUND

Neurotoxic chemical warfare agents can be classified as some of the most dangerous chemicals for humanity. The most effective of those agents are the Organophosphates (OPs) capable of restricting the enzyme Acetylcholinesterase (AChE), which in turn, controls the nerve impulse transmission. When AChE is inhibited by OPs, its reactivation can be usually performed through cationic oximes. However, until today, it has not been developed one universal defense agent, with complete effective reactivation activity for AChE inhibited by any of the many types of existing neurotoxic OPs. For this reason, before treating people intoxicated by an OP, it is necessary to determine the neurotoxic compound that was used for contamination, in order to select the most effective oxime. Unfortunately, this task usually requires a relatively long time, raising the possibility of death. Cationic oximes also display a limited capacity of permeating the Blood-Brain Barrier (BBB). This fact compromises their capacity to reactivating AChE inside the nervous system.

METHODS

We performed a comprehensive search on the data about OPs available on the scientific literature today in order to cover all the main drawbacks still faced in the research for the development of effective antidotes against those compounds.

RESULTS

Therefore, this review about neurotoxic OPs and the reactivation of AChE, provides insights for the new agents' development. The most expected defense agent is a molecule without toxicity and effective to reactivate AChE inhibited by all neurotoxic OPs.

CONCLUSION

To develop these new agents, the application of diverse scientific areas of research, especially theoretical procedures as computational science (computer simulation, docking and dynamics), organic synthesis, spectroscopic methodologies, biology, biochemical and biophysical information, medicinal chemistry, pharmacology and toxicology, is necessary.

Molecular Modeling and In Vitro Studies of a Neutral Oxime as a Potential Reactivator for Acetylcholinesterase Inhibited by Paraoxon

The present work aimed to compare the small, neutral and monoaromatic oxime, isatin-3-oxime (isatin-O), to the commercial ones, pralidoxime (2-PAM) and obidoxime, in a search for a new potential reactivator for acetylcholinesterase (AChE) inhibited by the pesticide paraoxon (AChE/POX) as well as a novel potential scaffold for further synthetic modifications. The multicriteria decision methods (MCDM) allowed the identification of the best docking poses of those molecules inside AChE/POX for further molecular dynamic (MD) studies, while Ellman's modified method enabled in vitro inhibition and reactivation assays. In corroboration with the theoretical studies, our experimental results showed that isatin-O have a reactivation potential capable of overcoming 2-PAM at the initial moments of the assay. Despite not achieving better results than obidoxime, this molecule is promising for being an active neutral oxime with capacity of crossing the blood⁻brain barrier (BBB), to reactivate AChE/POX inside the central and peripheral nervous systems. Moreover, the fact that isatin-O can also act as anticonvulsant makes this molecule a possible multipotent reactivator. Besides, the MCDM method showed to be an accurate method for the selection of the best docking poses generated in the docking studies.

Revealing the importance of linkers in K-series oxime reactivators for tabun-inhibited AChE using quantum chemical, docking and SMD studies

Inhibition of acetylcholinesterase (AChE) with organophosphorus compounds has a detrimental effect on human life. Oxime K203 seems to be one of the promising reactivators for tabun-inhibited AChE than (K027, K127, and K628). These reactivators differ only in the linker units between the two pyridinium rings. The conformational analyses performed with quantum chemical RHF/6-31G* level for K027, K127, K203 and K628 showed that the minimum energy conformers have different orientations of the active and peripheral pyridinium rings for these reactivator molecules. K203 with (-CH₂-CH=CH-CH₂-) linker unit possesses more open conformation compared to the other reactivators. Such orientation of K203 experiences favorable interaction with the surrounding residues of catalytic anionic site (CAS) and peripheral anionic site (PAS) of tabun-inhibited AChE. From the steered molecular dynamics simulations, it has been observed that the oxygen atom of the oxime group of K203 reactivator approaches nearest to the P-atom of the SUN203 (3.75 Å) at lower time scales (less than ~1000 ps) as compared to the other reactivators. K203 experiences less number of hydrophobic interaction with the PAS residues which is suggested to be an important factor for the efficient reactivation process. In addition, K203 crates large number of H-bonding with CAS residues SUN203, Phe295, Tyr337, Phe338 and His447. K203 barely changes its conformation during the SMD simulation process and hence the energy penalty to adopt any other conformation is minimal in this case as compared to the other reactivators. The molecular mechanics and Poisson-Boltzmann surface area binding energies obtained for the interaction of K203 inside the gorge of tabun inhibited AChE is substantially higher (-290.2 kcal/mol) than the corresponding K628 reactivator (-260.4 kcal/mol), which also possess unsaturated aromatic linker unit.

Perspectives for the structure-based design of acetylcholinesterase reactivators

Rational design of active molecules through structure-based methods has been gaining adepts during the last decades due to the wider availability of protein structures, most of them conjugated with relevant ligands. Acetylcholinesterase (AChE) is a molecular target with a considerable amount of data related to its sequence and 3-dimensional structure. In addition, there are structural insights about the mechanism of action of the natural substrate and drugs used in Alzheimer's disease, organophosphorus compounds, among others. We looked for AChE structural data useful for in silico design of potential interacting molecules. In particular, we focused on information regarding the design of ligands aimed to reactivate AChE catalytic activity. The structures of 178 AChE were annotated and categorized on different subsets according to the nature of the ligand, source organisms and experimental details. We compared sequence homology among the active site from Torpedo californica, Mus musculus and Homo sapiens with the latter two species having the closest relationship (88.9% identity). In addition, the mechanism of organophosphorus binding and the design of effective reactivators are reviewed. A curated data collection obtained with information from several sources was included for researchers working on the field. Finally, a molecular dynamics simulation with human AChE indicated that the catalytic pocket volume stabilizes around 600 Å(3), providing additional clues for drug design.

Computational studies of acetylcholinesterase complexed with fullerene derivatives: a new insight for Alzheimer disease treatment

Here, we propose five fullerene (C60) derivatives as new drugs against Alzheimer's disease (AD). These compounds were designed to act as new human acetylcholinesterase (HssAChE) inhibitors by blocking its fasciculin II (FASII) binding site. Docking and molecular dynamic results show that our proposals bind to the HssAChE tunnel entrance, forming stable complex, and further binding free energy calculations suggest that three of the derivatives proposed here could be potent HssAChE inhibitors. We found a region formed by a set of residues (Tyr72, Asp74, Trp286, Gln291, Tyr341, and Pro344) which can be further exploited in the drug design of new inhibitors of HssAChE based on C60 derivatives. Results presented here report for the first time by a new class of molecules that can become effective drugs against AD.

High-affinity selective inhibitor against phospholipase A2 (PLA2): a computational study

Phospholipase A2 (PLA2) is the most abundant protein found in snake venom. PLA2 induces a variety of pharmacological effects such as neurotoxicity, myotoxicity and cardiotoxicity as well as anticoagulant, hemolytic, anti-platelet, hypertensive, hemorrhagic and edema inducing effects. In this study, the three dimensional structure of PLA2 of Naja sputatrix (Malayan spitting cobra) was modeled by I-TASSER, SWISS-MODEL, PRIME and MODELLER programs. The best model was selected based on overall stereo-chemical quality. Further, molecular dynamics simulation was performed to know the stability of the modeled protein using Gromacs software. Average structure was generated during the simulation period of 10 ns. High throughput virtual screening was employed through different databases (Asinex, Hit finder, Maybridge, TOSLab and ZINC databases) against PLA2. The top seven compounds were selected based on the docking score and free energy binding calculations. These compounds were analyzed by quantum polarized ligand docking, induced fit docking and density functional theory calculation. Furthermore, the stability of lead molecules in the active site of PLA2 was employed by MD simulation. The results show that selected lead molecules were highly stable in the active site of PLA2.

Identification of potent inhibitors against snake venom metalloproteinase (SVMP) using molecular docking and molecular dynamics studies

Snake venom metalloproteinase (SVMP) (Echis coloratus (Carpet viper) is a multifunctional enzyme that is involved in producing several symptoms that follow a snakebite, such as severe local hemorrhage, nervous system effects and tissue necrosis. Because the three-dimensional (3D) structure of SVMP is not known, models were constructed, and the best model was selected based on its stereo-chemical quality. The stability of the modeled protein was analyzed through molecular dynamics (MD) simulation studies. Structure-based virtual screening was performed, and 15 potential molecules with the highest binding energies were selected. Further analysis was carried out with induced fit docking, Prime/MM-GBSA (ΔGBind calculations), quantum-polarized ligand docking, and density functional theory calculations. Further, the stability of the lead molecules in the SVMP-active site was examined using MD simulation. The results showed that the selected lead molecules were highly stable in the active site of SVMP. Hence, these molecules could potentially be selective inhibitors of SVMP. These lead molecules can be experimentally validated, and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for snake antivenom.

Acetylcholinesterase reactivators (HI-6, obidoxime, trimedoxime, K027, K075, K127, K203, K282): structural evaluation of human serum albumin binding and absorption kinetics

Acetylcholinesterase (AChE) reactivators (oximes) are compounds predominantly targeting the active site of the enzyme. Toxic effects of organophosphates nerve agents (OPNAs) are primarily related to their covalent binding to AChE and butyrylcholinesterase (BChE), critical detoxification enzymes in the blood and in the central nervous system (CNS). After exposure to OPNAs, accumulation of acetylcholine (ACh) overstimulates receptors and blocks neuromuscular junction transmission resulting in CNS toxicity. Current efforts at treatments for OPNA exposure are focused on non-quaternary reactivators, monoisonitrosoacetone oximes (MINA), and diacylmonoxime reactivators (DAM). However, so far only quaternary oximes have been approved for use in cases of OPNA intoxication. Five acetylcholinesterase reactivator candidates (K027, K075, K127, K203, K282) are presented here, together with pharmacokinetic data (plasma concentration, human serum albumin binding potency). Pharmacokinetic curves based on intramuscular application of the tested compounds are given, with binding information and an evaluation of structural relationships. Human Serum Albumin (HSA) binding studies have not yet been performed on any acetylcholinesterase reactivators, and correlations between structure, concentration curves and binding are vital for further development. HSA bindings of the tested compounds were 1% (HI-6), 7% (obidoxime), 6% (trimedoxime), and 5%, 10%, 4%, 15%, and 12% for K027, K075, K127, K203, and K282, respectively.

Structure- and ligand-based drug design of novel p38-alpha MAPK inhibitors in the fight against the Alzheimer's disease

Alzheimer's disease (AD) is characterized microscopically by the presence of amyloid plaques, which are accumulations of beta-amyloid protein inter-neurons, and neurofibrillary tangles formed predominantly by highly phosphorylated forms of the microtubule-associated protein, tau, which form tangled masses that consume neuronal cell body, possibly leading to neuronal dysfunction and ultimately death. p38α mitogen-activated protein kinase (MAPK) has been implicated in both events associated with AD, tau phosphorylation and inflammation. p38α MAPK pathway is activated by a dual phosphorylation at Thr180 and Tyr182 residues. Drug design of p38α MAPK inhibitors is mainly focused on small molecules that compete for Adenosine triphosphate in the catalytic site. Here, we used different approaches of structure- and ligand-based drug design and medicinal chemistry strategies based on a selected p38α MAPK structure deposited in the Protein Data Bank in complex with inhibitor, as well as others reported in literature. As a result of the virtual screening experiments performed here, as well as molecular dynamics, molecular interaction fields studies, shape and electrostatic similarities, activity and toxicity predictions, and pharmacokinetic and physicochemical properties, we have selected 13 compounds that meet the criteria of low or no toxicity potential, good pharmacotherapeutic profile, predicted activities, and calculated values ​​comparable with those obtained for the reference compounds, while maintaining the main interactions observed for the most potent inhibitors.

Molecular Modeling Studies of Piperidine Derivatives as New Acetylcholinesterase Inhibitors against Neurodegenerative Diseases

Neurodegenerative disorders are related to the progressive loss of structure or function and, eventually, death of neurons. These processes are responsible for diseases like Parkinson’s, Alzheimer’s, and Huntington’s, and the main molecular target for the drug design against these illnesses today is the enzyme acetylcholinesterase (AChE). Following this line, in the present work, we applied docking techniques to study some piperidine derivative inhibitors of AChE and further propose structures of six new AChE inhibitors as potential new drugs against neurodegenerative disorders. The best inhibitor proposed was submitted to additional molecular dynamics simulations steps.

Aging mechanism of soman inhibited acetylcholinesterase

Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nervous system that hydrolyzes neurotransmitter acetylcholine (ACh) and terminates synaptic signals. The catalytic serine of AChE can be phosphonylated by soman, one of the most potent nerve agents, and subsequently undergo an aging reaction. This phosphonylation and aging process leads to irreversible AChE inhibition, results in accumulation of excess ACh at the synaptic clefts, and causes neuromuscular paralysis. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the aging mechanism of soman phosphonylated AChE and determined its free energy profile. This aging reaction starts with the scission of the O2-Cα bond, which is followed by methyl migration, and results in a tertiary carbenium intermediate. At the transition state, the scissile O2-Cα bond is already cleaved with an average O-C distance of 3.2 ± 0.3 Å and the migrating methyl group is shared between Cα and Cβ carbons with C-C distances of 1.9 ± 0.1 and 1.8 ± 0.1 Å, respectively. The negatively charged phosphonate group is stabilized by a salt bridge with the imidazole ring of the catalytic histidine. A major product of aging, 2,3-dimethyl-2-butanol can be formed swiftly by the reaction of a water molecule. Our characterized mechanism and simulation results provide new detailed insights into this important biochemical process.