Differential binding of bispyridinium oxime drugs with acetylcholinesterase

FDA-Approved Oximes and Their Significance in Medicinal Chemistry

Despite the scientific advancements, organophosphate (OP) poisoning continues to be a major threat to humans, accounting for nearly one million poisoning cases every year leading to at least 20,000 deaths worldwide. Oximes represent the most important class in medicinal chemistry, renowned for their widespread applications as OP antidotes, drugs and intermediates for the synthesis of several pharmacological derivatives. Common oxime based reactivators or nerve antidotes include pralidoxime, obidoxime, HI-6, trimedoxime and methoxime, among which pralidoxime is the only FDA-approved drug. Cephalosporins are β-lactam based antibiotics and serve as widely acclaimed tools in fighting bacterial infections. Oxime based cephalosporins have emerged as an important class of drugs with improved efficacy and a broad spectrum of anti-microbial activity against Gram-positive and Gram-negative pathogens. Among the several oxime based derivatives, cefuroxime, ceftizoxime, cefpodoxime and cefmenoxime are the FDA approved oxime-based antibiotics. Given the pharmacological significance of oximes, in the present paper, we put together all the FDA-approved oximes and discuss their mechanism of action, pharmacokinetics and synthesis.

Designed inhibitors with hetero linkers for gastric proton pump H+,K+-ATPase: Steered molecular dynamics and metadynamics studies

Acid suppressant SCH28080 and its derivatives reversibly reduce acid secretion activity of the H⁺,K⁺-ATPase in a K⁺ competitive manner. The results on homologation of the SCH28080 by varying the linker chain length suggested the improvement in efficacy. However, the pharmacokinetic studies reveal that the hydrophobic nature of the CH₂ linker units may not help it to function as a better acid suppressant. We have exploited the role of linker unit to enhance the efficacy of such reversible acid suppressant drug molecules using hetero linker, i.e., disulfide and peroxy linkers. The logarithm of partition coefficient defined for a drug molecule relates to the partition coefficient, which allows the optimum solubility characteristics to reach the active site. The logarithm of partition coefficient calculated for the designed inhibitors suggests that inhibitors would possibly reach the active site in sufficient concentration like in the case of SCH28080. The steered molecular dynamics studies have revealed that the Inhibitor-1 with disulfide linker unit is more stable at the active site due to greater noncovalent interactions compared to the SCH28080. Centre of mass distance analysis suggests that the Cysteine-813 amino acid residue selectively plays an important role in the inhibition of H⁺,K⁺-ATPase for Inhibitor-1. Furthermore, the quantum chemical calculations with M11L/6-31+G(d,p) level of theory have been performed to account the noncovalent interactions responsible for the stabilization of inhibitor molecules in the active site gorge of the gastric proton pump at different time scale. The hydrogen bonding and hydrophobic interaction studies corroborate the center of mass distance analysis as well. Well-tempered metadynamics free energy surface and center of mass separation analysis for the Inhibitor-1 is in good agreement with the steered molecular dynamics results. The torsional angle of the linker units seems to be crucial for better efficacy of drug molecules. The torsional angle of linker units of SCH28080 (COCH₂C) and of Inhibitor 1 (CSSC) prefers to lie within ∼60°-90° for a longer time during the simulations, whereas, the peroxy linker (COOC) of Inhibitor 2 prefers to adopt ∼120-160°. Therefore, it appears that the smaller torsion angle of linker units can achieve better interactions with the active site residues of H⁺,K⁺-ATPase to inhibit the acid secretion activity. The reversible drug molecules with disulfide linker unit would be a promising candidate as proton pump antagonist to H⁺,K⁺-ATPase.

Computational Studies on Acetylcholinesterases

Functions of biomolecules, in particular enzymes, are usually modulated by structural fluctuations. This is especially the case in a gated diffusion-controlled reaction catalyzed by an enzyme such as acetylcholinesterase. The catalytic triad of acetylcholinesterase is located at the bottom of a long and narrow gorge, but it catalyzes the extremely rapid hydrolysis of the neurotransmitter, acetylcholine, with a reaction rate close to the diffusion-controlled limit. Computational modeling and simulation have produced considerable advances in exploring the dynamical and conformational properties of biomolecules, not only aiding in interpreting the experimental data, but also providing insights into the internal motions of the biomolecule at the atomic level. Given the remarkably high catalytic efficiency and the importance of acetylcholinesterase in drug development, great efforts have been made to understand the dynamics associated with its functions by use of various computational methods. Here, we present a comprehensive overview of recent computational studies on acetylcholinesterase, expanding our views of the enzyme from a microstate of a single structure to conformational ensembles, strengthening our understanding of the integration of structure, dynamics and function associated with the enzyme, and promoting the structure-based and/or mechanism-based design of new inhibitors for it.

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.

Unbinding of fluorinated oxime drug from the AChE gorge in polarizable water: a well-tempered metadynamics study

A well-tempered metadynamics study reveals that fluorinated obidoxime is held more firmly in the AChE gorge in comparison to obidoxime.

The reactivation of tabun-inhibited mutant AChE with Ortho-7: steered molecular dynamics and quantum chemical studies

A highly toxic nerve agent, tabun, can inhibit acetylcholinesterase (AChE) at cholinergic sites, which leads to serious cardiovascular complications, respiratory compromise and death. We have examined the structural features of the tabun-conjugated AChE complex with an oxime reactivator, Ortho-7, to provide a strategy for designing new and efficient reactivators. Mutation of mAChE within the choline binding site by Y337A and F338A and its interaction with Ortho-7 has been investigated using steered molecular dynamics (SMD) and quantum chemical methods. The overall study shows that after mutagenesis (Y337A), the reactivator can approach more freely towards the phosphorylated active site of serine without any significant steric hindrance in the presence of tabun compared to the wild type and double mutant. Furthermore, the poor binding of Ortho-7 with the peripheral residues of mAChE in the case of the single mutant compared to that of the wild-type and double mutant (Y337A/F338A) can contribute to better efficacy in the former case. Ortho-7 has formed a greater number of hydrogen bonds with the active site surrounding residues His447 and Phe295 in the case of the single mutant (Y337A), and that stabilizes the drug molecule for an effective reactivation process. The DFT M05-2X/6-31+G(d) level of theory shows that the binding energy of Ortho-7 with the single mutant (Y337A) is energetically more preferred (-19.8 kcal mol(-1)) than the wild-type (-8.1 kcal mol(-1)) and double mutant (Y337A/F338A) (-16.0 kcal mol(-1)). The study reveals that both the orientation of the oxime reactivator for nucleophilic attack and the stabilization of the reactivator at the active site would be crucial for the design of an efficient reactivator.

Can hydroxylamine be a more potent nucleophile for the reactivation of tabun-inhibited AChE than prototype oxime drugs? An answer derived from quantum chemical and steered molecular dynamics studies

Organophosphorus nerve agents are highly toxic compounds which strongly inhibit acetylcholinesterase (AChE) in the blood and in the central nervous system (CNS). Tabun is one of the highly toxic organophosphorus (OP) compounds and is resistant to many oxime drugs formulated for the reactivation of AChE. The reactivation mechanism of tabun-conjugated AChE with various drugs has been examined with density functional theory and ab initio quantum chemical calculations. The presence of a lone-pair located on the amidic group resists the nucleophilic attack at the phosphorus center of the tabun-conjugated AChE. We have shown that the newly designed drug candidate N-(pyridin-2-yl)hydroxylamine, at the MP2/6-31+G*//M05-2X/6-31G* level in the aqueous phase with the polarizable continuum solvation model (PCM), is more effective in reactivating the tabun-conjugated AChE than typical oxime drugs. The rate determining activation barrier with N-(pyridin-2-yl)hydroxylamine was found to be ∼1.7 kcal mol(-1), which is 7.2 kcal mol(-1) lower than the charged oxime trimedoxime (one of the most efficient reactivators in tabun poisonings). The greater nucleophilicity index (ω(-)) and higher CHelpG charge of pyridinylhydroxylamine compared to TMB4 support this observation. Furthermore, we have also examined the reactivation process of tabun-inhibited AChE with some other bis-quaternary oxime drug candidates such as methoxime (MMB4) and obidoxime. The docking analysis suggests that charged bis-quaternary pyridinium oximes have greater binding affinity inside the active-site gorge of AChE compared to the neutral pyridinylhydroxylamine. The peripheral ligand attached to the neutral pyridinylhydroxylamine enhanced the binding with the aromatic residues in the active-site gorge of AChE through effective π-π interactions. Steered molecular dynamics (SMD) simulations have also been performed with the charged oxime (TMB4) and the neutral hydroxylamine. From protein-drug interaction parameters (rupture force profiles, hydrogen bonds, hydrophobic interactions), geometry and the orientation of the drug candidates, the hydroxylamine is suggested to orchestrate the reactivation process better than TMB4. Furthermore, the calculated log P values show the effective penetration of the neutral drug candidate through the blood-brain barrier. The toxicity measurements and the IC50 values (a measure of the intrinsic affinity toward AChE) suggest that the pyridinylhydroxylamine compound could have similar toxic behavior compared to the prototype oxime antidotes used for reactivation purposes. The newly designed pyridinylhydroxylamine drug candidate can be an effective antidote both kinetically and structurally to reactivate the tabun-inhibited enzyme.

Ortho-7 bound to the active-site gorge of free and OP-conjugated acetylcholinesterase: cation-π interactions

Despite the immense importance of cation-π interactions prevailing in bispyridinium drug acetylcholinesterase (AChE) complexes, a precise description of cation-π interactions at molecular level has remained elusive. Here, we consider a bispyridinium drug, namely, ortho-7 in three different structures of AChE, with and without complexation with organophosphorus (OP) compounds for detailed investigation using all atom molecular dynamics simulation. By quantum mechanical calculations, Y72, W86, Y124, W286, Y337, and Y341 aromatic residues of the enzyme are investigated for possible cation-π interactions with ortho-7. The cation-π interactions in each of the protein-drug complexes are studied using distance, angle, a suitable functional form of them, and electrostatic criteria. The variation of cation-π functional is remarkably consistent with that of the Columbic variation. It is clearly observed that cation-π interactions for some of the residues in the catalytic active site (CAS) and peripheral anionic site (PAS) of the enzyme are either enhanced or reduced based on the nature of OP conjugation (i.e., nerve gas, tabun or pesticide, fenamiphos) when compared with the OP-free enzyme. The strength of cation-π interaction is strongly dependent on the type OP conjugation. The effect of conjugation at CAS is also seen to influence the cation-π interaction at the PAS region. The variation of cation-π interactions on the type of conjugating OP compounds might be suggestive of a reason as to why wide spectrum drug against any OP poisoning is yet to arrive in the market.

Quantum chemical and steered molecular dynamics studies for one pot solution to reactivate aged acetylcholinesterase with alkylator oxime

Dimethyl(pyridin-2-yl)sulfonium based oxime has been designed to reverse the aging process of organophosphorus inhibited AChE and to reactivate the aged-AChE adduct. We have employed DFT M05-2X/6-31G(∗) level of theory in aqueous phase with polarizable continuum solvation model (PCM) for the methylation of phosphonate ester monoanion of the soman-aged adduct. The calculated free energy of activation for the methyl transfer process from designed dimethyl(phenyl)sulfonium (1) to aged AChE-OP adduct occurs with a barrier of 31.4kcal/mol at M05-2X/6-31G(∗) level of theory, which is 6.4kcal/mol lower compared to the aging process signifies the preferential reversal process to recover the aged AChE-OP adduct. The pyridine ring containing alkylated sulfonium species, dimethyl(pyridin-2-yl)sulfonium (2), reduced the free energy of activation by 4.4kcal/mol compared to the previously reported alkylating agent N-methyl-2-methoxypyridinium species (A) for the alkylation of aged AChE-OP adduct. The free enzyme can be liberated from the inhibited acetylcholinesterase with the sulfonium compound decorated with an oxime group to avoid the administration of oxime drugs separately. The calculated potential energy surfaces show that the oxime based sulfonium compound (3) can effectively methylate the aged phosphonate ester monoanion of soman aged-adduct. The calculated global reactivity descriptors of the oxime 3 also shed light on this observation. To gain better understanding for protein drug interaction as well as the unbinding and conformational changes of the drug candidate in the active site of cholinesterase, steered molecular dynamics (SMD) simulation with 3 has been performed. Through a protein-drug interaction parameters (rupture force profiles, hydrogen bonds, hydrophobic interactions), geometrical and the orientation of drug-like candidate, the oxime 3 suggests to orchestrate the better reactivation process. The docking studies have been performed with 3 in the aged AChE and BChE to obtain the initial geometry of the SMD studies. The docking methods adopted in this study have been verified with the available crystal geometry of 1-methyl-2-(pentafluorobenzyloxyimino)pyridinium compound in aged soman inhibited human BChE (PDB code: 4B0P). The computational study suggests that the newly designed oxime is a potential candidate to reactivate the aged-AChE adduct.

Unbinding free energy of acetylcholinesterase bound oxime drugs along the gorge pathway from metadynamics-umbrella sampling investigation

Because of the pivotal role that the nerve enzyme, acetylcholinesterase plays in terminating nerve impulses at cholinergic synapses. Its active site, located deep inside a 20 Å gorge, is a vulnerable target of the lethal organophosphorus compounds. Potent reactivators of the intoxicated enzyme are nucleophiles, such as bispyridinium oxime that binds to the peripheral anionic site and the active site of the enzyme through suitable cation-π interactions. Atomic scale molecular dynamics and free energy calculations in explicit water are used to study unbinding pathways of two oxime drugs (Ortho-7 and Obidoxime) from the gorge of the enzyme. The role of enzyme-drug cation-π interactions are explored with the metadynamics simulation. The metadynamics discovered potential of mean force (PMF) of the unbinding events is refined by the umbrella sampling (US) corrections. The bidimensional free energy landscape of the metadynamics runs are further subjected to finite temperature string analysis to obtain the transition tube connecting the minima and bottlenecks of the unbinding pathway. The PMF is also obtained from US simulations using the biasing potential constructed from the transition tube and are found to be consistent with the metadynamics-US corrected results. Although experimental structural data clearly shows analogous coordination of the two drugs inside the gorge in the bound state, the PMF of the drug trafficking along the gorge pathway point, within an equilibrium free energy context, to a multistep process that differs from one another. Routes, milestones and subtlety toward the unbinding pathway of the two oximes at finite temperature are identified.

Functional Analysis and Molecular Docking studies of Medicinal Compounds for AChE and BChE in Alzheimer's Disease and Type 2 Diabetes Mellitus

Acetylcholinesterase and Butyrylcholinesterase share unravelling link with components of metabolic syndromes that's characterised by low levels of HDL cholesterol, obesity, high fast aldohexose levels, hyper-trigliceridaemia and high blood pressure, by regulation of cholinergic transmission and therefore the enzyme activity within a living system. The phosphomotifs associated with amino acid and tyrosine binding motifs in AChE and BChE were known to be common. Phylogenetic tree was constructed to these proteins usinf UPGMA and Maximum Likelihood methods in MEGA software has shown interaction of AChE and BChE with ageing diseases like Alzheimer's disease and Diabetes. AChE has shown closely related to BChE, retinol dehydrogenase and β-polypeptide. The present studies is also accomplished that AChE, BChE, COLQ, HAND1, APP, NLGN2 and NGF proteins has interactions with diseases such as Alzheimer's and D2M using Pathwaylinker and STRING. Medicinal compounds like Ortho-7, Dibucaine and HI-6 are predicted as good targets for modeled AChE and BChE proteins based on docking studies. Hence perceptive studies of cholinesterase structure and the biological mechanisms of inhibition are necessary for effective drug development.

Energetics of Ortho-7 (oxime drug) translocation through the active-site gorge of tabun conjugated acetylcholinesterase

Oxime drugs translocate through the 20 Å active-site gorge of acetylcholinesterase in order to liberate the enzyme from organophosphorus compounds' (such as tabun) conjugation. Here we report bidirectional steered molecular dynamics simulations of oxime drug (Ortho-7) translocation through the gorge of tabun intoxicated enzyme, in which time dependent external forces accelerate the translocation event. The simulations reveal the participation of drug-enzyme hydrogen bonding, hydrophobic interactions and water bridges between them. Employing nonequilibrium theorems that recovers the free energy from irreversible work done, we reconstruct potential of mean force along the translocation pathway such that the desired quantity represents an unperturbed system. The potential locates the binding sites and barriers for the drug to translocate inside the gorge. Configurational entropic contribution of the protein-drug binding entity and the role of solvent translational mobility in the binding energetics is further assessed.

New way of direct nitrogen atom phenylation in quinoline derivatives

Comparison of ion-molecular reactions of free-phenyl cations generated by tritium β -decay with 2-methyl- and 2-phenylquinolines has been investigated. The reaction of direct nitrogen atom phenylation with the help of nucleogenic phenyl cations has been fulfilled for the first time and a new one-step synthesis of tritium-labeled N-phenyl-2-phenylquinolinium salt-lipophilic radioactive biological marker has been elaborated.