Active-site gorge and buried water molecules in crystal structures of acetylcholinesterase from Torpedo californica

Derivatives of 9-phosphorylated acridine as butyrylcholinesterase inhibitors with antioxidant activity and the ability to inhibit β-amyloid self-aggregation: potential therapeutic agents for Alzheimer's disease

We investigated the inhibitory activities of novel 9-phosphoryl-9,10-dihydroacridines and 9-phosphorylacridines against acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and carboxylesterase (CES). We also studied the abilities of the new compounds to interfere with the self-aggregation of β-amyloid (Aβ42) in the thioflavin test as well as their antioxidant activities in the ABTS and FRAP assays. We used molecular docking, molecular dynamics simulations, and quantum-chemical calculations to explain experimental results. All new compounds weakly inhibited AChE and off-target CES. Dihydroacridines with aryl substituents in the phosphoryl moiety inhibited BChE; the most active were the dibenzyloxy derivative 1d and its diphenethyl bioisostere 1e (IC50 = 2.90 ± 0.23 µM and 3.22 ± 0.25 µM, respectively). Only one acridine, 2d, an analog of dihydroacridine, 1d, was an effective BChE inhibitor (IC50 = 6.90 ± 0.55 μM), consistent with docking results. Dihydroacridines inhibited Aβ42 self-aggregation; 1d and 1e were the most active (58.9% ± 4.7% and 46.9% ± 4.2%, respectively). All dihydroacridines 1 demonstrated high ABTS•+-scavenging and iron-reducing activities comparable to Trolox, but acridines 2 were almost inactive. Observed features were well explained by quantum-chemical calculations. ADMET parameters calculated for all compounds predicted favorable intestinal absorption, good blood-brain barrier permeability, and low cardiac toxicity. Overall, the best results were obtained for two dihydroacridine derivatives 1d and 1e with dibenzyloxy and diphenethyl substituents in the phosphoryl moiety. These compounds displayed high inhibition of BChE activity and Aβ42 self-aggregation, high antioxidant activity, and favorable predicted ADMET profiles. Therefore, we consider 1d and 1e as lead compounds for further in-depth studies as potential anti-AD preparations.

Resveratrol–Maltol and Resveratrol–Thiophene Hybrids as Cholinesterase Inhibitors and Antioxidants: Synthesis, Biometal Chelating Capability and Crystal Structure

New resveratrol–thiophene and resveratrol–maltol hybrids were synthesized as cholinesterase inhibitors and antioxidants. As with photostability experiments, biological tests also found remarkable differences in the properties and behavior of thiophene and maltol hybrids. While resveratrol–thiophene hybrids have excellent inhibitory and antioxidant properties (similar to the activity of reference drug galantamine), maltols have been proven to be weaker inhibitors and antioxidants. The molecular docking of selected active ligands gave insight into the structures of docked enzymes. It enabled the identification of interactions between the ligand and the active site of both cholinesterases. The maltols that proved to be active cholinesterase inhibitors were able to coordinate Fe³⁺ ion, forming complexes of 1:1 composition. Their formation constants, determined by spectrophotometry, are very similar, lgK = 11.6–12.6, suggesting that Fe³⁺ binds to the common hydroxy-pyranone moiety and is hardly affected by the other aromatic part of the ligand. Accordingly, the characteristic bands in their individual absorption spectra are uniformly red-shifted relative to those of the free ligands. The crystal structures of two new resveratrol–maltol hybrids were recorded, giving additional information on the molecules’ intermolecular hydrogen bonds and packing. In this way, several functionalities of these new resveratrol hybrids were examined as a necessary approach to finding more effective drugs for complicated neurodegenerative diseases.

Novel 3-aminobenzofuran derivatives as multifunctional agents for the treatment of Alzheimer’s disease

A novel multifunctional series of 3-aminobenzofuran derivatives 5a-p were designed and synthesized as potent inhibitors of acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE). The target compounds 5a-p were prepared via a three-step reaction, starting from 2-hydroxy benzonitrile. In vitro anti-cholinesterase activity exhibited that most of the compounds had potent acetyl- and butyrylcholinesterase inhibitory activity. In particular, compound 5f containing 2-fluorobenzyl moiety showed the best inhibitory activity. Furthermore, this compound showed activity on self- and AChE-induced Aβ-aggregation and MTT assay against PC12 cells. The kinetic study revealed that compound 5f showed mixed-type inhibition on AChE. Based on these results, compound 5f can be considered as a novel multifunctional structural unit against Alzheimer’s disease.

Multiadducts of C60 Modulate Amyloid-β Fibrillation with Dual Acetylcholinesterase Inhibition and Antioxidant Properties: In Vitro and In Silico Studies

Background: Amyloid-β (Aβ) fibrils induce cognitive impairment and neuronal loss, leading to onset of Alzheimer’s disease (AD). The inhibition of Aβ aggregation has been proposed as a therapeutic strategy for AD. Pristine C60 has shown the capacity to interact with the Aβ peptide and interfere with fibril formation but induces significant toxic effects in vitro and in vivo. Objective: To evaluate the potential of a series of C60 multiadducts to inhibit the Aβ fibrillization. Methods: A series of C60 multiadducts with four to six diethyl malonyl and their corresponding disodium-malonyl substituents were synthesized as individual isomers. Their potential on Aβ fibrillization inhibition was evaluated in vitro, in cellulo, and silico. Antioxidant activity, acetylcholinesterase inhibition capacity, and toxicity were assessed in vitro. Results: The multiadducts modulate Aβ fibrils formation without inducing cell toxicity, and that the number and polarity of the substituents play a significant role in the adducts efficacy to modulate Aβ aggregation. The molecular mechanism of fullerene-Aβ interaction and modulation was identified. Furthermore, the fullerene derivatives exhibited antioxidant capacity and reduction of acetylcholinesterase activity. Conclusion: Multiadducts of C60 are novel multi-target-directed ligand molecules that could hold considerable promise as the starting point for the development of AD therapies.

Computational Studies on Acetylcholinesterase Inhibitors: From Biochemistry to Chemistry

Acetylcholinesterase inhibitors are the most promising therapeutics for Alzheimer's disease treatment as these prevent the loss of acetylcholine and slows the progression of the disease. The drugs approved for the management of Alzheimer's disease by the FDA are acetylcholinesterase inhibitors but are associated with side effects. Consistent and stringent efforts by the researchers with the help of computational methods opened new ways of developing novel molecules with good acetylcholinesterase inhibitory activity. In this manuscript, we reviewed the studies that identified the essential structural features of acetylcholinesterase inhibitors at the molecular level as well as the techniques like molecular docking, molecular dynamics, quantitative structure-activity relationship, virtual screening, and pharmacophore modelling that were used in designing these inhibitors.

A Comprehensive Review of Cholinesterase Modeling and Simulation

The present article reviews published efforts to study acetylcholinesterase and butyrylcholinesterase structure and function using computer-based modeling and simulation techniques. Structures and models of both enzymes from various organisms, including rays, mice, and humans, are discussed to highlight key structural similarities in the active site gorges of the two enzymes, such as flexibility, binding site location, and function, as well as differences, such as gorge volume and binding site residue composition. Catalytic studies are also described, with an emphasis on the mechanism of acetylcholine hydrolysis by each enzyme and novel mutants that increase catalytic efficiency. The inhibitory activities of myriad compounds have been computationally assessed, primarily through Monte Carlo-based docking calculations and molecular dynamics simulations. Pharmaceutical compounds examined herein include FDA-approved therapeutics and their derivatives, as well as several other prescription drug derivatives. Cholinesterase interactions with both narcotics and organophosphate compounds are discussed, with the latter focusing primarily on molecular recognition studies of potential therapeutic value and on improving our understanding of the reactivation of cholinesterases that are bound to toxins. This review also explores the inhibitory properties of several other organic and biological moieties, as well as advancements in virtual screening methodologies with respect to these enzymes.

Torpedo californica acetylcholinesterase is stabilized by binding of a divalent metal ion to a novel and versatile 4D motif

Stabilization of Torpedo californica acetylcholinesterase by the divalent cations Ca⁺², Mg⁺², and Mn⁺² was investigated. All three substantially protect the enzyme from thermal inactivation. Electron paramagnetic resonance revealed one high‐affinity binding site for Mn⁺² and several much weaker sites. Differential scanning calorimetry showed a single irreversible thermal transition. All three cations raise both the temperature of the transition and the activation energy, with the transition becoming more cooperative. The crystal structures of the Ca⁺² and Mg⁺² complexes with Torpedo acetylcholinesterase were solved. A principal binding site was identified. In both cases, it consists of four aspartates (a 4D motif), within which the divalent ion is embedded, together with several water molecules. It makes direct contact with two of the aspartates, and indirect contact, via waters, with the other two. The 4D motif has been identified in 31 acetylcholinesterase sequences and 28 butyrylcholinesterase sequences. Zebrafish acetylcholinesterase also contains the 4D motif; it, too, is stabilized by divalent metal ions. The ASSAM server retrieved 200 other proteins that display the 4D motif, in many of which it is occupied by a divalent cation. It is a very versatile motif, since, even though tightly conserved in terms of RMSD values, it can contain from one to as many as three divalent metal ions, together with a variable number of waters. This novel motif, which binds primarily divalent metal ions, is shared by a broad repertoire of proteins. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:Protein_Science:3.

PDB‐ID(s): 7B38, 7B8E and 7B2W;

Molecular modeling-guided optimization of acetylcholinesterase reactivators: A proof for reactivation of covalently inhibited targets

Covalent drugs have been intensively studied in some very important fields such as anti-tumor and anti-virus, including the currently global-spread SARS-CoV-2. However, these drugs may interact with a variety of biological macromolecules and cause serious toxicology, so how to reactivate the inhibited targets seems to be imperative in the near future. Organophosphate was an extreme example, which could form a covalent bound easily with acetylcholinesterase and irreversibly inhibited the enzyme, causing high toxicology. Some nucleophilic oxime reactivators for organophosphate poisoned acetylcholinesterase had been developed, but the reactivation process was still less understanding. Herein, we proposed there should be a pre-reactivated pose during the reactivating process and compounds whose binding pose was easy to transfer to the pre-reactivated pose might be efficient reactivators. Then we refined the previous reactivators based on the molecular dynamic simulation results, the resulting compounds L7R3 and L7R5 were proven as much more efficient reactivators for organophosphate inhibited acetylcholinesterase than currently used oximes. This work might provide some insights for constructing reactivators of covalently inhibited targets by using computational methods.

Chiral Separation, X-ray Structure, and Biological Evaluation of a Potent and Reversible Dual Binding Site AChE Inhibitor

Acetylcholinesterase (AChE) inhibitors (AChEIs) still remain the leading therapeutic options for the symptomatic treatment of cognitive deficits associated with mild-to-moderate Alzheimer's disease. The search for new AChEIs benefits from well-established knowledge of the molecular interactions of selective AChEIs, such as donepezil and related dual binding site inhibitors. Starting from a previously disclosed coumarin-based inhibitor (±)-cis-1, active as racemate in the nanomolar range toward AChE, we proceeded on a double track by (i) achieving chiral resolution of the enantiomers of 1 by HPLC and (ii) preparing two close achiral analogues of 1, i.e., compounds 4 and 6. An eudismic ratio as high as 20 was observed for the (-) enantiomer of cis-1. The X-ray crystal structure of the complex between the (-)-cis-1 eutomer (coded as MC1420) and T. californica AChE was determined at 2.8 Å, and docking calculation results suggested that the eutomer in (1R,3S) absolute configuration should be energetically more favored in binding the enzyme than the eutomer in (1S,3R) configuration. The achiral analogues 4 and 6 were less effective in inhibiting AChE compared to (±)-cis-1, but interestingly butylamide 4 emerged as a potent inhibitor of butyrylcholinesterase (BChE).

Molecular docking and in vitro evaluation of a new hybrid molecule (JM-20) on cholinesterase activity from different sources

The main function of AChE is the hydrolysis of the neurotransmitter acetylcholine (ACh) at the neuromuscular and in cholinergic brain synapses. In some pathologies, loss of cholinergic neurons may be associated with a deficiency of ACh in specific brain areas. Consequently, the study of new safe drugs that inhibit AChE is important, because they can increase ACh levels in the synaptic cleft without adverse effects. Here, we evaluated the effects of JM-20 (a benzodiazepine-dihydropyridine hybrid molecule) on cholinesterase (ChE) activities from distinct sources (AChE from Electrophorus electricus (EeAChE), human erythrocyte membranes (HsAChE (ghost)), total erythrocyte (HsAChE (erythrocyte)) and BChE from plasma (HsBChE) and purified enzyme from the horse (EcBChE)). Kinetic parameters were determined in the presence of 0.05-1.6 mM of substrate concentration. The interactions ChEs with JM-20 were performed using molecular docking simulations. JM-20 inhibited all tested AChE but not BChE. The IC₅₀ values were 123 nM ± 0.2 (EeAChE), 158 nM ± 0.1 (ghost HsAChE), and 172 nM ± 0.2 (erythrocytic HsAChE). JM-20 caused a mixed type of inhibition (it altered K_m and V_max of AChE). The molecular docking indicated the binding poses and the most plausible active isomer of JM-20. Besides giving important data for future drug design, our results help us understand the mode of action of JM-20 as a specific inhibitor of AChE enzymes.

Multidomain Convergence of Argonaute during RISC Assembly Correlates with the Formation of Internal Water Clusters

Despite the relevance of Argonaute proteins in RNA silencing, little is known about the structural steps of small RNA loading to form RNA-induced silencing complexes (RISCs). We report the 1.9 Å crystal structure of human Argonaute4 with guide RNA. Comparison with the previously determined apo structure of Neurospora crassa QDE2 revealed that the PIWI domain has two subdomains. Binding of guide RNA fastens the subdomains, thereby rearranging the active-site residues and increasing the affinity for TNRC6 proteins. We also identified two water pockets beneath the nucleic acid-binding channel that appeared to stabilize the mature RISC. Indeed, mutating the water-pocket residues of Argonaute2 and Argonaute4 compromised RISC assembly. Simulations predict that internal water molecules are exchangeable with the bulk solvent but always occupy specific positions at the domain interfaces. These results suggest that after guide RNA-driven conformational changes, water-mediated hydrogen-bonding networks tie together the converged domains to complete the functional RISC structure.

Water structure changes in oxime-mediated reactivation process of phosphorylated human acetylcholinesterase

The role of water in oxime-mediated reactivation of phosphylated cholinesterases (ChEs) has been asked with recurrence. To investigate oximate water structure changes in this reaction, reactivation of paraoxon-inhibited human acetylcholinesterase (AChE) was performed by the oxime asoxime (HI-6) at different pH in the presence and absence of lyotropic salts: a neutral salt (NaCl), a strong chaotropic salt (LiSCN) and strong kosmotropic salts (ammonium sulphate and phosphate HPO42-). At the same time, molecular dynamic (MD) simulations of enzyme reactivation under the same conditions were performed over 100 ns. Reactivation kinetics showed that the low concentration of chaotropic salt up to 75 mM increased the percentage of reactivation of diethylphosphorylated AChE whereas kosmotropic salts lead only to a small decrease in reactivation. This indicates that water-breaker salt induces destructuration of water molecules that are electrostricted around oximate ions. Desolvation of oximate favors nucleophilic attack on the phosphorus atom. Effects observed at high salt concentrations (>100 mM) result either from salting-out of the enzyme by kosmotropic salts (phosphate and ammonium sulphate) or denaturing action of chaotropic LiSCN. MDs simulations of diethylphosphorylated hAChE complex with HI-6 over 100 ns were performed in the presence of 100 mM (NH4)2SO4 and 50 mM LiSCN. In the presence of LiSCN, it was found that protein and water have a higher mobility, i.e. water is less organized, compared with the ammonium sulphate system. LiSCN favors protein solvation (hydrophobic hydration) and breakage of elelectrostricted water molecules around of oximate ion. As a result, more free water molecules participated to reaction steps accompanying oxime-mediated dephosphorylation.

Kinetic and structural studies on the interactions of Torpedo californica acetylcholinesterase with two donepezil-like rigid analogues

Acetylcholinesterase inhibitors were introduced for the symptomatic treatment of Alzheimer’s disease (AD). Among the currently approved inhibitors, donepezil (DNP) is one of the most preferred choices in AD therapy. The X-ray crystal structures of Torpedo californica AChE in complex with two novel rigid DNP-like analogs, compounds 1 and 2, have been determined. Kinetic studies indicated that compounds 1 and 2 show a mixed-type inhibition against TcAChE, with K_i values of 11.12 ± 2.88 and 29.86 ± 1.12 nM, respectively. The DNP rigidification results in a likely entropy-enthalpy compensation with solvation effects contributing primarily to AChE binding affinity. Molecular docking evidenced the molecular basis for the binding of compounds 1 and 2 to the active site of β-secretase-1. Overall, these simplified DNP derivatives may represent new structural templates for the design of lead compounds for a more effective therapeutic strategy against AD by foreseeing a dual AChE and BACE-1 inhibitory activity.

Comparison of the Binding of Reversible Inhibitors to Human Butyrylcholinesterase and Acetylcholinesterase: A Crystallographic, Kinetic and Calorimetric Study

Acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) hydrolyze the neurotransmitter acetylcholine and, thereby, function as coregulators of cholinergic neurotransmission. Although closely related, these enzymes display very different substrate specificities that only partially overlap. This disparity is largely due to differences in the number of aromatic residues lining the active site gorge, which leads to large differences in the shape of the gorge and potentially to distinct interactions with an individual ligand. Considerable structural information is available for the binding of a wide diversity of ligands to AChE. In contrast, structural data on the binding of reversible ligands to BChE are lacking. In a recent effort, an inhibitor competition approach was used to probe the overlap of ligand binding sites in BChE. Here, we extend this study by solving the crystal structures of human BChE in complex with five reversible ligands, namely, decamethonium, thioflavin T, propidium, huprine, and ethopropazine. We compare these structures to equivalent AChE complexes when available in the protein data bank and supplement this comparison with kinetic data and observations from isothermal titration calorimetry. This new information now allows us to define the binding mode of various ligand families and will be of importance in designing specific reversible ligands of BChE that behave as inhibitors or reactivators.

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.

Assessing the binding of cholinesterase inhibitors by docking and molecular dynamics studies

In this report we assessed by docking and molecular dynamics the binding mechanisms of three FDA-approved Alzheimer drugs, inhibitors of the enzyme acetylcholinesterase (AChE): donepezil, galantamine and rivastigmine. Dockings by the softwares Autodock-Vina, PatchDock and Plant reproduced the docked conformations of the inhibitor-enzyme complexes within 2Å of RMSD of the X-ray structure. Free-energy scores show strong affinity of the inhibitors for the enzyme binding pocket. Three independent Molecular Dynamics simulation runs indicated general stability of donepezil, galantamine and rivastigmine in their respective enzyme binding pocket (also referred to as gorge) as well as the tendency to form hydrogen bonds with the water molecules. The binding of rivastigmine in the Torpedo California AChE binding pocket is interesting as it eventually undergoes carbamylation and breaks apart according to the X-ray structure of the complex. Similarity search in the ZINC database and targeted docking on the gorge region of the AChE enzyme gave new putative inhibitor molecules with high predicted binding affinity, suitable for potential biophysical and biological assessments.

Recent developments in structural studies on acetylcholinesterase

This review focuses on several recent developments concerning structure-function relationships in vertebrate acetylcholinesterase. These include studies on high-resolution structures of human acetylcholinesterase and its complexes; the first crystal structure of a snake venom acetylcholinesterase, in which open and closed states of the 'back door' are visualized; a powerful algorithm for redesigning proteins for enhanced expression in prokaryotic systems, as applied to human acetylcholinesterase, which has hitherto been an intractable target; in situ implementation of 'click chemistry' in crystalline acetylcholinesterase, which yields novel insights into the steric and dynamic changes involved in the reaction within the active-site gorge; and a study that demonstrates the effect of crystallization conditions on ligand alignment within a protein complex, in this case the methylene blue-Torpedo californica acetylcholinesterase complex, which highlights the relevance of the precipitant employed to structure-based drug design. This is an article for the special issue XVth International Symposium on Cholinergic Mechanisms.

On the Influence of the Protonation States of Active Site Residues on AChE Reactivation: A QM/MM Approach

Acetylcholinesterase (AChE), an enzyme of the serine hydrolase superfamily, is a mediator of signal transmission at cholinergic synapses by catalyzing acetylcholine cleavage into acetate and choline. This enzyme is vulnerable to covalent inhibition by organophosphate compounds (like VX). Covalent inhibition of AChE does not revert spontaneously. Known reactivator compounds have limited action in restoring catalytic activity. QM/MM simulations of VX-inhibited AChE reactivation by pralidoxime (2-PAM), a classical reactivator, were performed. These afforded a broad view of the effect of protonation states of active-site residues, and provide evidence for the role of Glu202, which needs to be protonated for reactivation to occur. In situ deprotonation of 2-PAM for both protonation states of Glu202 showed that His447 is able to deprotonate 2-PAM with the assistance of Glu202. Because the active site of serine hydrolases is highly conserved, this work provides new insights on the interplay between the catalytic triad residues and this glutamate, newly identified as protonatable.

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.

Automated Structure- and Sequence-Based Design of Proteins for High Bacterial Expression and Stability

Upon heterologous overexpression, many proteins misfold or aggregate, thus resulting in low functional yields. Human acetylcholinesterase (hAChE), an enzyme mediating synaptic transmission, is a typical case of a human protein that necessitates mammalian systems to obtain functional expression. We developed a computational strategy and designed an AChE variant bearing 51 mutations that improved core packing, surface polarity, and backbone rigidity. This variant expressed at ∼2,000-fold higher levels in E. coli compared to wild-type hAChE and exhibited 20°C higher thermostability with no change in enzymatic properties or in the active-site configuration as determined by crystallography. To demonstrate broad utility, we similarly designed four other human and bacterial proteins. Testing at most three designs per protein, we obtained enhanced stability and/or higher yields of soluble and active protein in E. coli. Our algorithm requires only a 3D structure and several dozen sequences of naturally occurring homologs, and is available at http://pross.weizmann.ac.il.

•

A new computational method is used to stabilize five recalcitrant proteins

•

Designed variants show higher expression and stability with unmodified function

•

A designed human acetylcholinesterase variant expresses solubly in bacteria

•

The method is fully automated and implemented on a webserver

The P450 CYP6Z1 confers carbamate/pyrethroid cross-resistance in a major African malaria vector beside a novel carbamate-insensitive N485I acetylcholinesterase-1 mutation

Carbamates are increasingly used for vector control notably in areas with pyrethroid resistance. However, a cross‐resistance between these insecticides in major malaria vectors such as Anopheles funestus could severely limit available resistance management options. Unfortunately, the molecular basis of such cross‐resistance remains uncharacterized in An. funestus, preventing effective resistance management. Here, using a genomewide transcription profiling, we revealed that metabolic resistance through upregulation of cytochrome P450 genes is driving carbamate resistance. The P450s CYP6P9a,CYP6P9b and CYP6Z1 were the most upregulated detoxification genes in the multiple resistant mosquitoes. However, in silico docking simulations predicted CYP6Z1 to metabolize both pyrethroids and carbamates, whereas CYP6P9a and CYP6P9b were predicted to metabolize only the pyrethroids. Using recombinant enzyme metabolism and inhibition assays, we demonstrated that CYP6Z1 metabolizes bendiocarb and pyrethroids, whereas CYP6P9a and CYP6P9b metabolize only the pyrethroids. Other upregulated gene families in resistant mosquitoes included several cuticular protein genes suggesting a possible reduced penetration resistance mechanism. Investigation of the target‐site resistance in acetylcholinesterase 1 (ace‐1) gene detected and established the association between the new N485I mutation and bendiocarb resistance (odds ratio 7.3; P < 0.0001). The detection of multiple haplotypes in single mosquitoes after cloning suggested the duplication of ace‐1. A TaqMan genotyping of the N485I in nine countries revealed that the mutation is located only in southern Africa with frequency of 10–15% suggesting its recent occurrence. These findings will help in monitoring the spread and evolution of carbamate resistance and improve the design of effective resistance management strategies to control this malaria vector.

Structure of a prereaction complex between the nerve agent sarin, its biological target acetylcholinesterase, and the antidote HI-6

Organophosphorus nerve agents interfere with cholinergic signaling by covalently binding to the active site of the enzyme acetylcholinesterase (AChE). This inhibition causes an accumulation of the neurotransmitter acetylcholine, potentially leading to overstimulation of the nervous system and death. Current treatments include the use of antidotes that promote the release of functional AChE by an unknown reactivation mechanism. We have used diffusion trap cryocrystallography and density functional theory (DFT) calculations to determine and analyze prereaction conformers of the nerve agent antidote HI-6 in complex with Mus musculus AChE covalently inhibited by the nerve agent sarin. These analyses reveal previously unknown conformations of the system and suggest that the cleavage of the covalent enzyme-sarin bond is preceded by a conformational change in the sarin adduct itself. Together with data from the reactivation kinetics, this alternate conformation suggests a key interaction between Glu202 and the O-isopropyl moiety of sarin. Moreover, solvent kinetic isotope effect experiments using deuterium oxide reveal that the reactivation mechanism features an isotope-sensitive step. These findings provide insights into the reactivation mechanism and provide a starting point for the development of improved antidotes. The work also illustrates how DFT calculations can guide the interpretation, analysis, and validation of crystallographic data for challenging reactive systems with complex conformational dynamics.

The impact of crystallization conditions on structure-based drug design: A case study on the methylene blue/acetylcholinesterase complex

Structure-based drug design utilizes apoprotein or complex structures retrieved from the PDB. >57% of crystallographic PDB entries were obtained with polyethylene glycols (PEGs) as precipitant and/or as cryoprotectant, but <6% of these report presence of individual ethyleneglycol oligomers. We report a case in which ethyleneglycol oligomers' presence in a crystal structure markedly affected the bound ligand's position. Specifically, we compared the positions of methylene blue and decamethonium in acetylcholinesterase complexes obtained using isomorphous crystals precipitated with PEG200 or ammonium sulfate. The ligands' positions within the active-site gorge in complexes obtained using PEG200 are influenced by presence of ethyleneglycol oligomers in both cases bound to W84 at the gorge's bottom, preventing interaction of the ligand's proximal quaternary group with its indole. Consequently, both ligands are ∼3.0Å further up the gorge than in complexes obtained using crystals precipitated with ammonium sulfate, in which the quaternary groups make direct π-cation interactions with the indole. These findings have implications for structure-based drug design, since data for ligand-protein complexes with polyethylene glycol as precipitant may not reflect the ligand's position in its absence, and could result in selecting incorrect drug discovery leads. Docking methylene blue into the structure obtained with PEG200, but omitting the ethyleneglycols, yields results agreeing poorly with the crystal structure; excellent agreement is obtained if they are included. Many proteins display features in which precipitants might lodge. It will be important to investigate presence of precipitants in published crystal structures, and whether it has resulted in misinterpreting electron density maps, adversely affecting drug design.

Synthesis and structure-activity relationship study of benzofuran-based chalconoids bearing benzylpyridinium moiety as potent acetylcholinesterase inhibitors

A series of benzofuran-based chalconoids 6a-v were designed and synthesized as new potential AChE inhibitors. The in vitro assay of synthesized compounds 6a-v showed that most compounds had significant anti-AChE activity at micromolar or sub-micromolar levels. Among the tested compounds, 3-pyridinium derivative 6m bearing N-(2-bromobenzyl) moiety and 7-methoxy substituent on the benzofuran ring exhibited superior activity. This compound with IC₅₀ value of 0.027 μM was as potent as standard drug donepezil.

Statistical survey of the buried waters in the Protein Data Bank

The structures of buried water molecules were studied in an ensemble of high-quality and non-redundant protein crystal structures. Buried water molecules were clustered and classified in lake-like clusters, which are completely isolated from the bulk solvent, and bay-like clusters, which are in contact with the bulk solvent through a surface water molecule. Buried water molecules are extremely common: lake-like clusters are found in 89 % of the protein crystal structures and bay-like clusters in 93 %. Clusters with only one water molecule are much more common than larger clusters. Both cluster types incline to be surrounded by loop residues, and to a minor extent by residues in extended secondary structure. Helical residues on the contrary do not tend to surround clusters of buried water molecules. One buried water molecule is found every 30-50 amino acid residues, depending on the secondary structures that are more abundant in the protein. Both main- and side-chain atoms are in contact with buried waters; they form four hydrogen bonds with the first water and 1-1.5 additional hydrogen bond for each additional water in the cluster. Consequently, buried water molecules appear to be firmly packed and rigid like the protein atoms. In this regard, it is remarkable to observe that prolines often surround water molecules buried in the protein interior. Interestingly, clusters of buried water molecules tend to be just beneath the protein surface. Moreover, water molecules tend to form a one-dimensional wire rather than more compact arrangements. This agrees with recent evidence of the mechanisms of solvent exchange between internal cavities and bulk solvent.

Synthesis and Evaluation of Chroman-4-One Linked to N-Benzyl Pyridinium Derivatives as New Acetylcholinesterase Inhibitors

A novel series of chroman-4-one derivatives containing the N-benzyl pyridinium moiety were designed, synthesized, and evaluated for their acetylcholinesterase (AChE) inhibitory activities. Among the various synthesized compounds, (E)-1-(2,3-dibromobenzyl)-4-((7-ethoxy-4-oxochroman-3-ylidene)methyl)pyridinium bromide (8l) depicted the most potent anti-AChE activity (IC50 = 0.048 μM). In addition, the molecular modeling study allowed us to detect possible binding modes that are in full compliance with the observed results through in vitro experiments.

Understanding the molecular mechanism of aryl acylamidase activity of acetylcholinesterase - An in silico study

Acetylcholinesterase (AChE) exhibits two different activities, namely esterase and aryl acylamidase (AAA). Unlike esterase, AAA activity of AChE is inhibited by the active site inhibitors while remaining unaffected by the peripheral anionic site inhibitors. This differential inhibitory pattern of active and peripheral anionic site inhibitors on the AAA activity remains unanswered. To answer this, we investigated the mechanism of binding and trafficking of AAA substrates using in silico tools. Molecular docking of serotonin and AAA substrates (o-nitroacetanilide, and o-nitrotrifluoroacetanilide,) onto AChE shows that these compounds bind at the side door of AChE. Thus, we conceived that the AAA substrates prefer the side door to reach the active site for their catalysis. Further, steered molecular dynamics simulations show that the force required for binding and trafficking of the AAA substrate through the side door is comparatively lesser than their dissociation (900kJ/mol/nm). Among the two substrates, o-nitrotrifluoroacetanilide required lesser force (380kJ/mol/nm) than o-nitroacetanilide the (550kJ/mol/nm) for its binding, thus validating o-nitrotrifluoroacetanilide as a better substrate. With these observations, we resolve that the AAA activity of AChE is mediated through its side door. Therefore, binding of PAS inhibitors at the main door of AChE remain ineffective against AAA activity.

Analysis of the activation of acetylcholinesterase by carbon nanoparticles using a monolithic immobilized enzyme microreactor: role of the water molecules in the active site gorge

A biochromatographic system was used to study the direct effect of carbon nanoparticles (CNPs) on the acetylcholinesterase (AChE) activity. The AChE enzyme was covalently immobilized on a monolithic CIM-disk via its NH2 residues. Our results showed an increase in the AChE activity in presence of CNPs. The catalytic constant (k(cat)) was increased while the Michaelis constant (K(m)) was slightly decreased. This indicated an increase in the enzyme efficiency with increase of the substrate affinity to the active site. The thermodynamic data of the activation mechanism of the enzyme, i.e. ΔH* and ΔS*, showed no change in the substrate interaction mechanism with the anionic binding site. The increase of the enthalpy (ΔH*) and the entropy (ΔS*) with decrease in the free energy of activation (Ea) was related to structural conformation change in the active site gorge. This affected the stability of water molecules in the active site gorge and facilitated water displacement by substrate for entering to the active site of the enzyme.

Indolinone-based acetylcholinesterase inhibitors: synthesis, biological activity and molecular modeling

A series of indolinone-based compounds bearing benzylpyridinium moiety was designed as dual-binding inhibitors of acetylcholinesterase (AChE). The target compounds 3a-u were synthesized by condensation of oxindole and pyridin-4-carbalehyde, and subsequent N-benzylation. The anti-cholinesterase activity evaluation of synthesized compounds revealed that most of them had very potent inhibitory activity against AChE, superior to standard drug donepezil. Particularly, 2-chlorobenzyl derivative 3c was the most potent compound against AChE with IC50 value of 0.44 nM, being 32-fold more potent than donepezil. Also, most of compounds were more potent than standard drug donepezil against butyrylcholinesterase (BuChE). Docking study revealed that the hydrophobic aromatic part (indoline) of representative compound 3c binds to the PAS and the N-benzylpyridinium residue binds to the CAS of AChE.

Crystal structures of human cholinesterases in complex with huprine W and tacrine: elements of specificity for anti-Alzheimer's drugs targeting acetyl- and butyryl-cholinesterase

The multifunctional nature of Alzheimer's disease calls for MTDLs (multitarget-directed ligands) to act on different components of the pathology, like the cholinergic dysfunction and amyloid aggregation. Such MTDLs are usually on the basis of cholinesterase inhibitors (e.g. tacrine or huprine) coupled with another active molecule aimed at a different target. To aid in the design of these MTDLs, we report the crystal structures of hAChE (human acetylcholinesterase) in complex with FAS-2 (fasciculin 2) and a hydroxylated derivative of huprine (huprine W), and of hBChE (human butyrylcholinesterase) in complex with tacrine. Huprine W in hAChE and tacrine in hBChE reside in strikingly similar positions highlighting the conservation of key interactions, namely, π-π/cation-π interactions with Trp86 (Trp82), and hydrogen bonding with the main chain carbonyl of the catalytic histidine residue. Huprine W forms additional interactions with hAChE, which explains its superior affinity: the isoquinoline moiety is associated with a group of aromatic residues (Tyr337, Phe338 and Phe295 not present in hBChE) in addition to Trp86; the hydroxyl group is hydrogen bonded to both the catalytic serine residue and residues in the oxyanion hole; and the chlorine substituent is nested in a hydrophobic pocket interacting strongly with Trp439. There is no pocket in hBChE that is able to accommodate the chlorine substituent.

Effects of viscosity and osmotic stress on the reaction of human butyrylcholinesterase with cresyl saligenin phosphate, a toxicant related to aerotoxic syndrome: kinetic and molecular dynamics studies

CSP (cresyl saligenin phosphate) is an irreversible inhibitor of human BChE (butyrylcholinesterase) that has been involved in the aerotoxic syndrome. Inhibition under pseudo-first-order conditions is biphasic, reflecting a slow equilibrium between two enzyme states E and E′. The elementary constants for CSP inhibition of wild-type BChE and D70G mutant were determined by studying the dependence of inhibition kinetics on viscosity and osmotic pressure. Glycerol and sucrose were used as viscosogens. Phosphorylation by CSP is sensitive to viscosity and is thus strongly diffusion-controlled (kon≈108 M−1·min−1). Bimolecular rate constants (ki) are about equal to kon values, making CSP one of the fastest inhibitors of BChE. Sucrose caused osmotic stress because it is excluded from the active-site gorge. This depleted the active-site gorge of water. Osmotic activation volumes, determined from the dependence of ki on osmotic pressure, showed that water in the gorge of the D70G mutant is more easily depleted than that in wild-type BChE. This demonstrates the importance of the peripheral site residue Asp70 in controlling the active-site gorge hydration. MD simulations provided new evidence for differences in the motion of water within the gorge of wild-type and D70G enzymes. The effect of viscosogens/osmolytes provided information on the slow equilibrium E⇌E′, indicating that alteration in hydration of a key catalytic residue shifts the equilibrium towards E′. MD simulations showed that glycerol molecules that substitute for water molecules in the enzyme active-site gorge induce a conformational change in the catalytic triad residue His438, leading to the less reactive form E′.

Structure-based search for new inhibitors of cholinesterases

Cholinesterases are important biological targets responsible for regulation of cholinergic transmission, and their inhibitors are used for the treatment of Alzheimer’s disease. To design new cholinesterase inhibitors, of different structure-based design strategies was followed, including the modification of compounds from a previously developed library and a fragment-based design approach. This led to the selection of heterodimeric structures as potential inhibitors. Synthesis and biological evaluation of selected candidates confirmed that the designed compounds were acetylcholinesterase inhibitors with IC₅₀ values in the mid-nanomolar to low micromolar range, and some of them were also butyrylcholinesterase inhibitors.

Design, synthesis and evaluation of novel heterodimers of donepezil and huperzine fragments as acetylcholinesterase inhibitors

Four series of novel heterodimers comprised of donepezil and huperzine A (HupA) fragments were designed, synthesized, and evaluated in search of potent acetylcholinesterase (AChE) inhibitors as potential therapeutic treatment for Alzheimer's disease. Heterodimers comprised of dimethoxyindanone (from donepezil), hupyridone (from HupA), and connected with a multimethylene linker, were identified as potent and selective inhibitors of AChE. Diastereomeric heterodimers (RS,S)-17b (with a tetramethylene linker) exhibited the highest potency of inhibition towards AChE with an IC(50) value of 9 nM and no detectable inhibitory effect on butyrylcholinesterase at 1mM.

Back-scattering interferometry: an ultrasensitive method for the unperturbed detection of acetylcholinesterase-inhibitor interactions

A series of inhibitors of acetylcholinesterase (AChE) have been screened by back-scattering interferometry (BSI). Enzyme levels as low as 100 pM (22,000 molecules of AChE) can be detected. This method can be used to screen for mixed AChE inhibitors, agents that have shown high efficacy against Alzheimer's disease, by detecting dual-binding interactions. E = enzyme, I = inhibitor, S = substrate.

Mapping of the interaction sites of galanthamine: a quantitative analysis through pairwise potentials and quantum chemistry

A quantitative analysis of the interaction sites of the anti-Alzheimer drug galanthamine with molecular probes (water and benzene molecules) representative of its surroundings in the binding site of acetylcholinesterase (AChE) has been realized through pairwise potentials calculations and quantum chemistry. This strategy allows a full and accurate exploration of the galanthamine potential energy surface of interaction. Significantly different results are obtained according to the distances of approaches between the various molecular fragments and the conformation of the galanthamine N-methyl substituent. The geometry of the most relevant complexes has then been fully optimized through MPWB1K/6-31 + G(d,p) calculations, final energies being recomputed at the LMP2/aug-cc-pVTZ(-f) level of theory. Unexpectedly, galanthamine is found to interact mainly from its hydrogen-bond donor groups. Among those, CH groups in the vicinity of the ammonium group are prominent. The trends obtained provide rationales to the predilection of the equatorial orientation of the galanthamine N-methyl substituent for binding to AChE. The analysis of the interaction energies pointed out the independence between the various interaction sites and the rigid character of galanthamine. The comparison between the cluster calculations and the crystallographic observations in galanthamine-AChE co-crystals allows the validation of the theoretical methodology. In particular, the positions of several water molecules appearing as strongly conserved in galanthamine-AChE co-crystals are predicted by the calculations. Moreover, the experimental position and orientation of lateral chains of functionally important aminoacid residues are in close agreement with the ones predicted theoretically. Our study provides relevant information for a rational drug design of galanthamine based AChE inhibitors.