Rapid binding of a cationic active site inhibitor to wild type and mutant mouse acetylcholinesterase: Brownian dynamics simulation including diffusion in the active site gorge

It is known that anionic surface residues play a role in the long-range electrostatic attraction between acetylcholinesterase and cationic ligands. In our current investigation, we show that anionic residues also play an important role in the behavior of the ligand within the active site gorge of acetylcholinesterase. Negatively charged residues near the gorge opening not only attract positively charged ligands from solution to the enzyme, but can also restrict the motion of the ligand once it is inside of the gorge. We use Brownian dynamics techniques to calculate the rate constant k_on for wild type and mutant acetylcholinesterase with a positively charged ligand. These calculations are performed by allowing the ligand to diffuse within the active site gorge. This is an extension of previously reported work in which a ligand was allowed to diffuse only to the enzyme surface. By setting the reaction criteria for the ligand closer to the active site, better agreement with experimental data is obtained. Although a number of residues influence the movement of the ligand within the gorge, Asp74 is shown to play a particularly important role in this function. Asp74 traps the ligand within the gorge, and in this way helps to ensure a reaction. © 1998 John Wiley & Sons, Inc. Biopoly 46: 465–474, 1998     

Properties of water molecules in the active site gorge of acetylcholinesterase from computer simulation

A 10-ns trajectory from a molecular dynamics simulation is used to examine the structure and dynamics of water in the active site gorge of acetylcholinesterase to determine what influence water may have on its function. While the confining nature of the deep active site gorge slows down and structures water significantly compared to bulk water, water in the gorge is found to display a number of properties that may aid ligand entry and binding. These properties include fluctuations in the population of gorge waters, moderate disorder and mobility of water in the middle and entrance to the gorge, reduced water hydrogen-bonding ability, and transient cavities in the gorge.     

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.     

Pathways of ligand clearance in acetylcholinesterase by multiple copy sampling

The clearance of seven different ligands from the deeply buried active-site of Torpedo californica acetylcholinesterase is investigated by combining multiple copy sampling molecular dynamics simulations, with the analysis of protein-ligand interactions, protein motion and the electrostatic potential sampled by the ligand copies along their journey outwards. The considered ligands are the cations ammonium, methylammonium, and tetramethylammonium, the hydrophobic methane and neopentane, and the anionic product acetate and its neutral form, acetic acid. We find that the pathways explored by the different ligands vary with ligand size and chemical properties. Very small ligands, such as ammonium and methane, exit through several routes. One involves the main exit through the mouth of the enzyme gorge, another is through the so-called back door near Trp84, and a third uses a side door at a direction of approximately 45 degrees to the main exit. The larger polar ligands, methylammonium and acetic acid, leave through the main exit, but the bulkiest, tetramethylammonium and neopentane, as well as the smaller acetate ion, remain trapped in the enzyme gorge during the time of the simulations. The pattern of protein-ligand contacts during the diffusion process is highly non-random and differs for different ligands. A majority is made with aromatic side-chains, but classical H-bonds are also formed. In the case of acetate, but not acetic acid, the anionic and neutral form, respectively, of one of the reaction products, specific electrostatic interactions with protein groups, seem to slow ligand motion and interfere with protein flexibility; protonation of the acetate ion is therefore suggested to facilitate clearance. The Poisson-Boltzmann formalism is used to compute the electrostatic potential of the thermally fluctuating acetylcholinesterase protein at positions actually visited by the diffusing ligand copies. Ligands of different charge and size are shown to sample somewhat different electrostatic potentials during their migration, because they explore different microscopic routes. The potential along the clearance route of a cation such as methylammonium displays two clear minima at the active and peripheral anionic site. We find moreover that the electrostatic energy barrier that the cation needs to overcome when moving between these two sites is small in both directions, being of the order of the ligand kinetic energy. The peripheral site thus appears to play a role in trapping inbound cationic ligands as well as in cation clearance, and hence in product release.     

Structure and dynamics of the active site gorge of acetylcholinesterase: synergistic use of molecular dynamics simulation and X-ray crystallography.

The active site of acetylcholinesterase (AChE) from Torpedo californica is located 20 A from the enzyme surface at the bottom of a narrow gorge. To understand the role of this gorge in the function of AChE, we have studied simulations of its molecular dynamics. When simulations were conducted with pure water filling the gorge, residues in the vicinity of the active site deviated quickly and markedly from the crystal structure. Further study of the original crystallographic data suggests that a bis-quaternary decamethonium (DECA) ion, acquired during enzyme purification, residues in the gorge. There is additional electron density within the gorge that may represent small bound cations. When DECA and 2 cations are placed within the gorge, the simulation and the crystal structure are dramatically reconciled. The small cations, more so than DECA, appear to stabilize part of the gorge wall through electrostatic interactions. This part of the gorge wall is relatively thin and may regulate substrate, product, and water movement through the active site.     

The role of Li+, Na+, and K+ in the ligand binding inside the human acetylcholinesterase gorge

Alkali cations can affect the catalytic efficiency of enzymes. This is particularly true when dealing with enzymes whose substrate bears a formal positive charge. Computational and biochemical approaches have been combined to shed light on the atomic aspects of the role of Li(+), Na(+), and K(+) on human acetylcholinesterase (hAChE) ligand binding. In this respect, molecular dynamics simulations and our recently developed metadynamics method were applied to study the entrance of the three cations in the gorge of hAChE, and their effect on the dynamical motion of a ligand (tetramethylammonium) from the bulk of the solvent into the deep narrow enzyme gorge. Furthermore, in order to support the theoretical results, K(M) and k(cat) for the acetylcholine hydrolysis in the presence of the three cations were evaluated by using an approach based on the Ellman's method. The combination of computational and biochemical experiments clearly showed that Li(+), Na(+), and K(+) may influence the ligand binding at the hAChE gorge.     

Long Route or Shortcut? A Molecular Dynamics Study of Traffic of Thiocholine within the Active-Site Gorge of Acetylcholinesterase

The principal role of acetylcholinesterase is termination of nerve impulse transmission at cholinergic synapses, by rapid hydrolysis of the neurotransmitter acetylcholine to acetate and choline. Its active site is buried at the bottom of a deep and narrow gorge, at the rim of which is found a second anionic site, the peripheral anionic site. The fact that the active site is so deeply buried has raised cogent questions as to how rapid traffic of substrate and products occurs in such a confined environment. Various theoretical and experimental approaches have been used to solve this problem. Here, multiple conventional molecular dynamics simulations have been performed to investigate the clearance of the product, thiocholine, from the active-site gorge of acetylcholinesterase. Our results indicate that thiocholine is released from the peripheral anionic site via random pathways, while three exit routes appear to be favored for its release from the active site, namely, along the axis of the active-site gorge, and through putative back- and side-doors. The back-door pathway is that via which thiocholine exits most frequently. Our results are in good agreement with kinetic and kinetic-crystallography studies. We propose the use of multiple molecular dynamics simulations as a fast yet accurate complementary tool in structural studies of enzymatic trafficking.     

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

Buried water molecules and the water molecules in the active-site gorge are analyzed for five crystal structures of acetylcholinesterase from Torpedo californica in the resolution range 2.2-2.5 A (native enzyme, and four inhibitor complexes). A total of 45 buried hydration sites are identified, which are populated with between 36 and 41 water molecules. About half of the buried water is located in a distinct region neighboring the active-site gorge. Most of the buried water molecules are very well conserved among the five structures, and have low displacement parameters, B, of magnitudes similar to those of the main-chain atoms of the central beta-sheet structure. The active-site gorge of the native enzyme is filled with over 20 water molecules, which have poor hydrogen-bond coordination with an average of 2.9 polar contacts per water molecule. Upon ligand binding, distinct groups of these water molecules are displaced, whereas the others remain in positions similar to those that they occupy in the native enzyme. Possible roles of the buried water molecules are discussed, including their possible action as a lubricant to allow large-amplitude fluctuations of the loop structures forming the gorge wall. Such fluctuations are required to facilitate traffic of substrate, products and water molecules to and from the active-site. Because of their poor coordination, the gorge water molecules can be considered as "activated" as compared to bulk water. This should allow their easy displacement by incoming substrate. The relatively loose packing of the gorge water molecules leaves numerous small voids, and more efficient space-filling by substrates and inhibitors may be a major driving force of ligand binding.     

Probing the active center gorge of acetylcholinesterase by fluorophores linked to substituted cysteines

To examine the influence of individual side chains in governing rates of ligand entry into the active center gorge of acetylcholinesterase and to characterize the dynamics and immediate environment of these residues, we have conjugated reactive groups with selected charge and fluorescence characteristics to cysteines substituted by mutagenesis at specific positions on the enzyme. Insertion of side chains larger than in the native tyrosine at position 124 near the constriction point of the active site gorge confers steric hindrance to affect maximum catalytic throughput (k(cat)/K(m)) and rates of diffusional entry of trifluoroketones to the active center. Smaller groups appear not to present steric constraints to entry; however, cationic side chains selectively and markedly reduce cation ligand entry through electrostatic repulsion in the gorge. The influence of side chain modification on ligand kinetics has been correlated with spectroscopic characteristics of fluorescent side chains and their capacity to influence the binding of a peptide, fasciculin, which inhibits catalysis peripherally by sealing the mouth of the gorge. Acrylodan conjugated to cysteine was substituted for tyrosine at position 124 within the gorge, for histidine 287 on the surface adjacent to the gorge and for alanine 262 on a mobile loop distal to the gorge. The 124 position reveals the most hydrophobic environment and the largest hypsochromic shift of the emission maximum with fasciculin binding. This finding likely reflects a sandwiching of the acrylodan in the complex with the tip of fasciculin loop II. An intermediate spectral shift is found for the 287 position, consistent with partial occlusion by loops II and III of fasciculin in the complex. Spectroscopic properties of the acrylodan at the 262 position are unaltered by fasciculin addition. Hence, combined spectroscopic and kinetic analyses reveal distinguishing characteristics in various regions of acetylcholinesterase that influence ligand association.     

Mouse acetylcholinesterase unliganded and in complex with huperzine A: A comparison of molecular dynamics simulations

A 1 ns molecular dynamics simulation of unliganded mouse acetylcholinesterase (AChE) is compared to a previous simulation of mouse AChE complexed with huperzine A (HupA). Several common features are observed. In both simulations, the active site gorge fluctuates in size during the 1 ns trajectory and is completely pinched off several times. Many of the residues in the gorge that formed hydrogen bonds with HupA in the simulation of the complex now form hydrogen bonds with other protein residues and water molecules in the gorge. The opening of a “backdoor” entrance to the active site that was found in the simulation of the complex is also observed in the unliganded simulation. Differences between the two simulations include overall lower structural rms deviations for residues in the gorge in the unliganded simulation, a smaller diameter of the gorge in the absence of HupA, and the disappearance of a side channel that was frequently present in the liganded simulation. The differences between the two simulations can be attributed, in part, to the interaction of AChE with HupA. © 1999 John Wiley & Sons, Inc. Biopoly 50: 35–43, 1999     

Acetylcholinesterase: Role of the enzyme's charge distribution in steering charged ligands toward the active site

The electrostatic steering of charged ligands toward the active site of Torpedo californica acetylcholinesterase is investigated by Brownian dynamics simulations of wild type enzyme and several mutated forms, in which some normally charged residues are neutralized. The simulations reveal that the total ligand influx through a surface of 42 Å radius centered in the enzyme monomer and separated from the protein surface by I-14 Å is not significantly influenced by electrostatic interactions. Electrostatic effects are visible for encounters with a surface of 32 Å radius, which is partially hidden inside the protein, but mostly within the solvent. A clear accumulation of encounter events for that sphere is observed in the area directly above the entrance to the active site gorge. In this area, the encounter events are increased by 40% compared to the case of a neutral ligand. However, the differences among the encounter rates for the various mutants considered here are not pronounced, all rate constants being within ±10% of the average value. The enzyme charge distribution becomes more important as the charged ligand moves toward the bottom of the gorge, where the active site is located. We show that neither the enzyme's total charge, nor its dipole moment, fully account for the electrostatic steering of ligand to the active site. Higher moments of the enzyme's charge distribution are also important. However, for a series of mutations for which the direction of the enzyme dipole moment is constant within a few degrees, one observes a gradual decrease in the diffusional encounter rate constant with the number of neutralized residues. On the other hand, for other mutants that change the direction of the dipole moment from that of the wild type, the calculated encounter rate constants can be very close to that of the wild type. The present work yields two new insights to the kinetics of acetylcholinesterase. First, evolution appears to have built a redundant electrostatic steering capability into this important enzyme through the overall distribution of its thousands of partially charged atoms. And second, roughly half of the rate enhancement due to electrostatics arises from steering of the substrate outside the enzyme; the other half of the rate enhancement arises from improved trapping of the substrate after it has entered the gorge. The computational results reproduce qualitatively, and help to rationalize, many surprising experimental results obtained recently for human acetylcholinesterase. © 1996 John Wiley & Sons, Inc.     

Nanosecond dynamics of acetylcholinesterase near the active center gorge

To delineate the role of peptide backbone flexibility and rapid molecular motion in acetylcholinesterase catalysis and inhibitor association, we investigated the decay of fluorescence anisotropy at three sites of fluorescein conjugation to cysteine-substitution mutants of the enzyme. One cysteine was placed in a loop at the peripheral site near the rim of the active center gorge (H287C); a second was in a helical region outside of the active center gorge (T249C); a third was at the tip of a small, flexible omega loop well separated from the gorge (A262C). Mutation and fluorophore conjugation did not appreciably alter catalytic or inhibitor binding parameters of the enzyme. The results show that each site examined was associated with a high degree of segmental motion; however, the A262C and H287C sites were significantly more flexible than the T249C site. Association of the active center inhibitor, tacrine, and the peripheral site peptide inhibitor, fasciculin, had no effect on the anisotropy decay of fluorophores at positions 249 and 262. Fasciculin, but not tacrine, on the other hand, dramatically altered the decay profile of the fluorophore at the 287 position, in a manner consistent with fasciculin reducing the segmental motion of the peptide chain in this local region. The results suggest that the motions of residues near the active center gorge and across from the Cys(69)-Cys(96) omega loop are uncoupled and that ligand binding at the active center or the peripheral site does not influence acetylcholinesterase conformational dynamics globally, but induces primarily domain localized decreases in flexibility proximal to the bound ligand.     

Analysis of a 10-ns molecular dynamics simulation of mouse acetylcholinesterase.

A 10-ns molecular dynamics simulation of mouse acetylcholinesterase was analyzed, with special attention paid to the fluctuation in the width of the gorge and opening events of the back door. The trajectory was first verified to ensure its stability. We defined the gorge proper radius as the measure for the extent of gorge opening. We developed an expression of an inter-atom distance representative of the gorge proper radius in terms of projections on the principal components. This revealed the fact that collective motions of many scales contribute to the opening behavior of the gorge. Covariance and correlation results identified the motions of the protein backbone as the gorge opens. In the back-door region, side-chain dihedral angles that define the opening were identified.     

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. Proteins 2014; 82:1799–1818. © 2014 Wiley Periodicals, Inc.     

Molecular dynamics simulations of human butyrylcholinesterase

Herein, we present results from molecular dynamics (MD) simulations of the human butyrylcholinesterase (BuChE) enzyme in aqueous solution. Two configurations of the unbound form of BuChE differing in the presence or absence of a sodium ion inside the protein gorge were simulated for 10 and 5 ns, respectively. Besides complementing the structural information provided by X-ray data, the MD simulations give insight into the structure of the native BuChE enzyme. For example, it is shown that: the nucleophilic Ser(198) residue and the various binding subsites in the BuChE catalytic cavity are readily accessible from the exterior of the protein; the presence of the sodium ion dynamically explores two different binding sites in the gorge leading to the active site and stabilizes the productive conformation of the Glu(325)/His(438)/Ser(198) catalytic triad; several long-lived water bridges are fully integrated into the architecture of the active site; the positions of the residues at the rim of the gorge region display large deviations with respect to the crystal structure; and two side doors, constituted by residues situated at the tip of the acyl- and Omega-loops, respectively, open wide enough to allow the passage of water molecules. In conclusion, we compare our theoretical results with those from previous work on mouse acetylcholinesterase and discuss their implications for substrate binding and catalysis in BuChE.     

Induced‐fit or preexisting equilibrium dynamics? Lessons from protein crystallography and MD simulations on acetylcholinesterase and implications for structure‐based drug design

Crystal structures of acetylcholinesterase complexed with ligands are compared with side-chain conformations accessed by native acetylcholinesterase in molecular dynamics (MD) simulations. Several crystallographic conformations of a key residue in a specific binding site are accessed in a simulation of native acetylcholinesterase, although not seen in rotomer plots. Conformational changes upon ligand binding thus involve preexisting equilibrium dynamics. Consequently, rational drug design could benefit significantly from conformations monitored by MD simulations of native targets.     

Molecular dynamics simulation of Axillaridine–A: A potent natural cholinesterase inhibitor

Molecular Dynamics (MD) simulations were carried out for human acetylcholinesterase (hAChE) and its complex with Axillaridine–A, in order to dynamically explore the active site of the protein and the behaviour of the ligand at the peripheral binding site. Simulation of the enzyme alone showed that the active site of AChE is located at the bottom of a deep and narrow cavity whose surface is lined with rings of aromatic residues while Tyr72 is almost perpendicular to the Trp286, which is responsible for stable π -π interactions. The complexation of AChE with Axillaridine-A, results in the reduction of gorge size due to interaction between the ligand and the active site residues. The gorge size was determined by the distance between the center of mass of Glu81 and Trp286. As far as the geometry of the active site is concerned, the presence of ligand in the active site alters its specific conformation, as revealed by stable hydrogen bondings established between amino acids. With the increasing interaction between ligand and the active amino acids, size of the active site of the complex decreases with respect to time. Axillaridine-A, forms stable π -π interactions with the aromatic ring of Tyr124 that results in inhibition of catalytic activity of the enzyme. This π -π interaction keeps the substrate stable at the edge of the catalytic gorge by inhibiting its catalytic activity. The MD results clearly provide an explanation for the binding pattern of bulky steroidal alkaloids at the active site of AChE.     

Molecular dynamics of mouse acetylcholinesterase complexed with huperzine A

Two molecular dynamics simulations were performed for a modeled complex of mouse acetylcholinesterase liganded with huperzine A (HupA). Analysis of these simulations shows that HupA shifts in the active site toward Tyr 337 and Phe 338, and that several residues in the active site area reach out to make hydrogen bonds with the inhibitor. Rapid fluctuations of the gorge width are observed, ranging from widths that allow substrate access to the active site, to pinched structures that do not allow access of molecules as small as water. Additional openings or channels to the active site are found. One opening is formed in the side wall of the active site gorge by residues Val 73, Asp 74, Thr 83, Glu 84, and Asn 87. Another opening is formed at the base of the gorge by residues Trp 86, Val 132, Glu 202, Gly 448, and Ile 451. Both of these openings have been observed separately in the Torpedo californica form of the enzyme. These channels could allow transport of waters and ions to and from the bulk solution. © 1999 John Wiley & Sons, Inc. Biopoly 50: 347–359, 1999     

Dynamic structure based pharmacophore modeling of the Acetylcholinesterase reveals several potential inhibitors

Acetylcholinesterase is a critical enzyme that regulates neurotransmission by catalyzing the breakdown of neurotransmitter acetylcholine in synapses of the nervous system. It is an important target for therapeutic drugs that treat Alzheimer’s disease. Since, the degree of flexibility of the side chains of the residues in the active-site gorge of Acetylcholinesterase is diverse it results in different bound ligand conformations. The side-chain conformations of Ser293, Tyr341, Leu76, and Val73 are flexible, while the side-chain conformations of Tyr72, Tyr 124, Ser125, Phe295, and Arg296 appear to be fixed. In this study, multi-conformation dynamic pharmacophore models from the donepezyl-binding pocket based on highly populated structures chosen from molecular dynamics simulations were used for screening compounds that can properly bind acetylcholinesterase. Based on these structures, three pharmacophore models were generated. Consequently, 14 hits were retrieved as final candidates by utilizing virtual screening of ZINC database and molecular docking.     

Backdoor opening mechanism in acetylcholinesterase based on X-ray crystallography and molecular dynamics simulations

The transient opening of a backdoor in the active‐site wall of acetylcholinesterase, one of nature's most rapid enzymes, has been suggested to contribute to the efficient traffic of substrates and products. A crystal structure of Torpedo californica acetylcholinesterase in complex with the peripheral‐site inhibitor aflatoxin is now presented, in which a tyrosine at the bottom of the active‐site gorge rotates to create a 3.4‐Å wide exit channel. Molecular dynamics simulations show that the opening can be further enlarged by movement of Trp84. The crystallographic and molecular dynamics simulation data thus point to the interface between Tyr442 and Trp84 as the key element of a backdoor, whose opening permits rapid clearance of catalysis products from the active site. Furthermore, the crystal structure presented provides a novel template for rational design of inhibitors and reactivators, including anti‐Alzheimer drugs and antidotes against organophosphate poisoning.