Crystal packing mediates enantioselective ligand recognition at the peripheral site of acetylcholinesterase

Recently, alkylene-linked heterodimers of tacrine (1) and 5-amino-5,6,7,8-tetrahydroquinolinone (2, hupyridone) were shown to exhibit higher acetylcholinesterase (AChE) inhibition than either monomeric 1 or 2. Such inhibitors are potential drug candidates for ameliorating the cognitive decrements in early Alzheimer patients. In an attempt to understand the inhibition mechanism of one such dimer, (RS)-(+/-)-N-9-(1,2,3,4-tetrahydroacridinyl)-N'-5-[5,6,7,8-tetrahydro-2'(1'H)-quinolinonyl]-1,10-diaminodecane [(RS)-(+/-)-3] bisoxalate, the racemate was soaked in trigonal Torpedo californica AChE (TcAChE) crystals, and the X-ray structure of the resulting complex was solved to 2.30 A resolution. Its structure revealed the 1 unit bound to the "anionic" subsite of the active site, near the bottom of the active-site gorge, as seen for the 1/TcAChE complex. Interestingly, only the (R)-enantiomer of the 2 unit was seen in the peripheral "anionic" site (PAS) at the top of the gorge, and was hydrogen-bonded to the side chains of residues belonging to an adjacent, symmetry-related AChE molecule covering the gorge entrance. When the same racemate was soaked in orthorhombic crystals of TcAChE, in which the entrance to the gorge is more exposed, the crystal structure of the corresponding complex revealed no substantial enantiomeric selectivity. This observation suggests that the apparent enantiomeric selectivity of trigonal crystals of TcAChE for (R)-3 is mainly due to crystal packing, resulting in preferential binding of one enantiomeric inhibitor both to its "host" enzyme and to its neighbor in the asymmetric unit, rather than to steric constraints imposed by the geometry of the active-site gorge.     

Crystal structure of thioflavin T bound to the peripheral site of Torpedo californica acetylcholinesterase reveals how thioflavin T acts as a sensitive fluorescent reporter of ligand binding to the acylation site

Acetylcholinesterase plays a key role in cholinergic synaptic transmission by hydrolyzing the neurotransmitter acetylcholine with one of the highest known catalytic rate constants. Hydrolysis occurs in a narrow and deep gorge that contains two sites of ligand binding: A peripheral site, or P-site, near the gorge entrance that contributes to catalytic efficiency both by transiently trapping substrate molecules as they enter the gorge and by allosterically accelerating the transfer of the substrate acyl group to a serine hydroxyl in an acylation site or A-site at the base of the gorge. Thioflavin T is a useful reporter of ligand interactions with the A-site. It binds specifically to the P-site with fluorescence that is enhanced approximately 1000-fold over that of unbound thioflavin T, and the enhanced fluorescence is quenched 1.5- to 4-fold when another ligand binds to the A-site in a ternary complex. To clarify the structural basis of this advantageous signal change, we here report the X-ray structure of the complex of thioflavin T with Torpedo californica acetylcholinesterase. The two aromatic rings in thioflavin T are coplanar and are packed snugly parallel to the aromatic side chains of Trp279, Tyr334, and Phe330. Overlays of this structure with the crystal structures of Torpedo californica acetylcholinesterase complexes with either edrophonium or m-( N, N, N-trimethylammonio)-2,2,2-trifluoroacetophenone, two small aromatic ligands that bind specifically to the A-site, indicate that the phenyl side chain of Phe330 must rotate to sterically accommodate both thioflavin T and the A-site ligand in the ternary complex. This rotation may allow some relaxation of the strict coplanarity of the aromatic rings in the bound thioflavin T and result in partial quenching of its fluorescence.     

Localization of a novel adhesion-promoting site on acetylcholinesterase using catalytic antiacetylcholinesterase antibodies displaying cholinesterase-like activity

A monoclonal antibody (MAb) raised against human acetylcholinesterase was found to have catalytic activity. A similar phenomenon was observed in a polyclonal antibody raised against the same antigen. The antibodies were demonstrated to be pure, and no contamination with either acetylcholinesterase or butyrylcholinesterase was found. Both antibodies hydrolyzed acetylthiocholine, anacetylcholinesterasesubstrate, and the MAb followed Michaelis-Menten kinetics. Six other MAbs and one other polyclonal antibody showed no evidence of catalytic activity. Acetylcholinesterase is a key component in the transmission of the nerve impulse, and is also expressed nonsynaptically during embryonic development, and abnormalities in expression are seen in neural tumors and degenerative disorders. This unusual expression is believed to be associated with a novel function of the enzyme related to differentiation and cell adhesion. Autoantibodies to acetylcholinesterase have been observed in a variety of neurologic, muscular, and autoimmune disorders. In an investigation of the possible role of acetylcholinesterase in cell adhesion, we showed that the enzyme promoted neurite outgrowth, in neuroblastoma cell lines, and conversely, that certain antiacetylcolinesterase antibodies abrogated cell-substrate adhesion. Interestingly, the antibodies most effective in this regard were catalytic. Preliminary epitopeanalysis indicated a conformational epitope in the N-terminal domain. This domain contains the active site, within a deep gorge, and the peripheral anionic site at the rim of the gorge. Peripheral-site inhibitors, but not active-site inhibitors, also interfered with adhesion, and competed with the catalytic monoclonal binding to acetylcholinesterase, indicating that the epitope recognized is associated with the peripheral anionic site. The inhibitor data also support the supposition that catalysis in these antibodies may have arisen from stable complexation of acetylcholinesterase with an inhibitor. We conclude that the catalytic antiacetylcholinesterase antibody interacts with structures associated with the peripheral anionic site, thus defining a novel site on the molecule involved, in cell adhesion. This finding has implications for our understanding of the potential importance of this peripheral site in a variety of congenital, neoplastic, and degenerative conditions.     

Crystal structure of human carboxylesterase 1 complexed with the Alzheimer's drug tacrine: from binding promiscuity to selective inhibition

Human carboxylesterase 1 (hCE1) is a broad-spectrum bioscavenger that plays important roles in narcotic metabolism, clinical prodrug activation, and the processing of fatty acid and cholesterol derivatives. We determined the 2.4 A crystal structure of hCE1 in complex with tacrine, the first drug approved for treating Alzheimer's disease, and compare this structure to the Torpedo californica acetylcholinesterase (AcChE)-tacrine complex. Tacrine binds in multiple orientations within the catalytic gorge of hCE1, while it stacks in the smaller AcChE active site between aromatic side chains. Our results show that hCE1's promiscuous action on distinct substrates is enhanced by its ability to interact with ligands in multiple orientations at once. Further, we use our structure to identify tacrine derivatives that act as low-micromolar inhibitors of hCE1 and may provide new avenues for treating narcotic abuse and cholesterol-related diseases.     

Flexibility of active-site gorge aromatic residues and non-gorge aromatic residues in acetylcholinesterase

The presence of an unusually large number of aromatic residues in the active site gorge of acetylcholinesterase is a subject of great interest. Flexibility of these residues has been suspected to be a key player in controlling the ligand traversal in the gorge. This raises the question of whether the over-representation of aromatic residues in the gorge implies higher-than-normal flexibility of these residues. The current study suggests that it does not. Large changes in the hydrophobic cross-sectional area due to dihedral oscillations are probably the reason of their presence in the gorge.     

Acetylcholinesterase: enhanced fluctuations and alternative routes to the active site in the complex with fasciculin-2

A 15 ns molecular dynamics simulation is reported for the complex of mouse acetylcholinesterase (mAChE) and the protein neurotoxin fasciculin-2. As compared to a 15 ns simulation of apo-mAChE, the structural fluctuations of the enzyme are substantially increased in magnitude for the enzyme in the complex. Fluctuations of part of the long omega loop (residues 69-96) are particularly enhanced. This loop forms one wall of the active site, and the enhanced fluctuations lead to additional routes of access to the active site.     

How does huperzine A enter and leave the binding gorge of acetylcholinesterase? Steered molecular dynamics simulations

The entering and leaving processes of Huperzine A (HupA) binding with the long active-site gorge of Torpedo californica acetylcholinesterase (TcAChE) have been investigated by using steered molecular dynamics simulations. The analysis of the force required along the pathway shows that it is easier for HupA to bind to the active site of AChE than to disassociate from it, which for the first time interprets at the atomic level the previous experimental result that unbinding process of HupA is much slower than its binding process to AChE. The direct hydrogen bonds, water bridges, and hydrophobic interactions were analyzed during two steered molecular dynamics (SMD) simulations. Break of the direct hydrogen bond needs a great pulling force. The steric hindrance of bottleneck might be the most important factor to produce the maximal rupture force for HupA to leave the binding site but it has a little effect on the binding process of HupA with AChE. Residue Asp72 forms a lot of water bridges with HupA leaving and entering the AChE binding gorge, acting as a clamp to take out HupA from or put HupA into the active site. The flip of the peptide bond between Gly117 and Gly118 has been detected during both the conventional MD and SMD simulations. The simulation results indicate that this flip phenomenon could be an intrinsic property of AChE and the Gly117-Gly118 peptide bond in both HupA bound and unbound AChE structures tends to adopt the native enzyme structure. At last, in a vacuum the rupture force is increased up to 1500 pN while in water solution the greatest rupture force is about 800 pN, which means water molecules in the binding gorge act as lubricant to facilitate HupA entering or leaving the binding gorge.     

The complex of a bivalent derivative of galanthamine with torpedo acetylcholinesterase displays drastic deformation of the active-site gorge: implications for structure-based drug design

Bifunctional derivatives of the alkaloid galanthamine, designed to interact with both the active site of the enzyme acetylcholinesterase (AChE) and its peripheral cation binding site, have been assayed with Torpedo californica AChE (TcAChE), and the three-dimensional structures of their complexes with the enzyme have been solved by X-ray crystallography. Differences were noted between the IC(50) values obtained for TcAChE and those for Electrophorus electricus AChE. These differences are ascribed to sequence differences in one or two residues lining the active-site gorge of the enzyme. The binding of one of the inhibitors disrupts the native conformation of one wall of the gorge, formed by the loop Trp279-Phe290. It is proposed that flexibility of this loop may permit the binding of inhibitors such as galanthamine, which are too bulky to penetrate the narrow neck of the gorge formed by Tyr121 and Phe330 as seen in the crystal structure.     

Binding free energy calculations to rationalize the interactions of huprines with acetylcholinesterase

In the present study, the binding free energy of a family of huprines with acetylcholinesterase (AChE) is calculated by means of the free energy perturbation method, based on hybrid quantum mechanics and molecular mechanics potentials. Binding free energy calculations and the analysis of the geometrical parameters highlight the importance of the stereochemistry of huprines in AChE inhibition. Binding isotope effects are calculated to unravel the interactions between ligands and the gorge of AChE. New chemical insights are provided to explain and rationalize the experimental results. A good correlation with the experimental data is found for a family of inhibitors with moderate differences in the enzyme affinity. The analysis of the geometrical parameters and interaction energy per residue reveals that Asp72, Glu199, and His440 contribute significantly to the network of interactions between active site residues, which stabilize the inhibitors in the gorge. It seems that a cooperative effect of the residues of the gorge determines the affinity of the enzyme for these inhibitors, where Asp72, Glu199, and His440 make a prominent contribution.     

Dynamic mechanism of E2020 binding to acetylcholinesterase: a steered molecular dynamics simulation

The unbinding process of E2020 ((R,S)-1-benzyl-4-[(5,6-dimethoxy-1-indanon)-2-yl]-methylpiperidine) leaving from the long active site gorge of Torpedo californica acetylcholinesterase (TcAChE) was studied by using steered molecular dynamics (SMD) simulations on a nanosecond scale with different velocities, and unbinding force profiles were obtained. Different from the unbinding of other AChE inhibitors, such as Huperzine A that undergoes the greatest barrier located at the bottleneck of the gorge, the major resistance preventing E2020 from leaving the gorge is from the peripheral anionic site where E2020 interacts intensively with several aromatic residues (e.g., Tyr70, Tyr121, and Trp279) through its benzene ring and forms a strong direct hydrogen bond and a water bridge with Ser286 via its O24. These interactions cause the largest rupture force, approximately 550 pN. It was found that the rotatable bonds of the piperidine ring to the benzene ring and dimethoxyindanone facilitate E2020 to pass the bottleneck through continuous conformation change by rotating those bonds to avoid serious conflict with Tyr121 and Phe330. The aromatic residues lining the gorge wall are the major components contributing to hydrophobic interactions between E2020 and TcAChE. Remarkably, these aromatic residues, acting in three groups as "sender" and "receiver", compose a "conveyer belt" for E2020 entering and leaving the TcAChE gorge.     

Influence of ligand binding on the conformation of Torpedo californica acetylcholinesterase

With the aim of identifying structural changes in acetylcholinesterase, induced by ligand binding, we use a completely automatic procedure to analyse the differences between the backbone conformation of the free enzyme and those in eight complexes of Torpedo californica acetylcholinesterase, with various quaternary ammonium ligands, and with the protein inhibitor fasciculin. In order to discriminate between structural changes due to ligand binding and those arising from model imprecision, we also examine protein–ligand and protein–water contacts. Except for the peptide flip in the complex with huperzine A, the backbones of other complexes with quaternary ammonium ligands display negligible changes relative to the free enzyme. Another exception is the complex with the bisquaternary ammonium ligand decamethonium, where several loops display above average deformations, but only two, those spanning residues 334–348 and residues 277–304, seem to move as a result of ligand binding. Movement of the ω loop (residues 61–95) is detected only in the complex with the protein fasciculin. Received: 21 July 2000 / Accepted: 18 September 2000 / Published online: 28 February 2001     

Conformational remodeling of femtomolar inhibitor-acetylcholinesterase complexes in the crystalline state

The active center of acetylcholinesterase (AChE), a target site for competitive inhibitors, resides centrosymmetric to the subunit at the base of a deep, narrow gorge lined by aromatic residues. At the gorge entry, a peripheral site encompasses overlapping binding loci for noncompetitive inhibitors, which alter substrate access to the gorge. The click-chemistry inhibitor TZ2PA6 links the active center ligand, tacrine, to the peripheral site ligand, propidium, through a biorthogonal reaction of an acetylene and an azide that forms either a syn1 or an anti1 triazole. Compared with wild-type mouse AChE, a Tyr337Ala mutant displays full catalytic activity, albeit with 2-3 orders of magnitude higher affinities for the TZ2PA6 syn1 and anti1 regioisomers, reflected in low femtomolar K(d) values, diffusion-limited association, and dissociation half-times greater than 1 month and 1 week, respectively. Three structures of each of the co-crystallized syn1 and anti1 complexes of the Tyr337Ala mutant were solved at three distinct times of crystal maturation, consistent with or exceeding the half-lives of the complexes in solution, while crystalline complexes obtained from soaked Tyr337Ala crystals led to picturing "freshly formed" complexes. The structures, at 2.55-2.75 ? resolution, reveal a range of unprecedented conformations of the bound regioisomers, not observed in the wild-type AChE complexes, associated with concerted positional rearrangements of side chains in the enzyme gorge. Moreover, time-dependent conformational remodeling of the crystalline complexes appears to correlate with the dissociation half-times of the solution complexes. Hence, for the tight-binding TZ2PA6 inhibitors, the initial complexes kinetically driven in solution slowly form more stable complexes governed by thermodynamic equilibrium and observable in mature crystals.     

Activation mechanisms of butyrylcholinesterase by 2,4,6-trinitrotoluene, 3,3-dimethylbutyl-N-n-butylcarbamate, and 2-trimethylsilyl-ethyl-N-n-butylcarbamate

The goal of this work was to propose a possible mechanism for the butyrylcholinesterase activation by 2,4,6-trinitrotoluene (TNT), 3,3-dimethylbutyl-N-n-butylcarbamate (1), and 2-trimethylsilyl-ethyl-N-n-butylcarbamate (2). Kinetically, TNT, and compounds 1 and 2 were characterized as the nonessential activators of butyrylcholinesterase. TNT, and compounds 1 and 2 were hydrophobic compounds and were proposed to bind to the hydrophobic activator binding site, which was located outside the active site gorge of the enzyme. The conformational change from a normal active site gorge to a more accessible active site gorge of the enzyme was proposed after binding of TNT, and compounds 1 and 2 to the activator binding site of the enzyme. Therefore, TNT, and compounds 1 and 2 may act as the excess of butyrylcholine in the substrate activator for the butyrylcholinesterase catalyzed reactions.     

Kinetics and Mechanisms of Acetylcholinesterase and Butyrylcholinesterase Inhibition by Cardiovascular Drugs and Benzodiazepines

The goal of this work is to determine enzyme kinetics and mechanisms of acetylcholinesterase and butyrylcholinesterase inhibition by five cardiovascular drugs, lovastatin, simvastatin, amlodipine besylate, nifedipine, and hydralazine hydrochloride, and two benzodiazepines, diazepam and chlordiazepoxide hydrochloride. All drugs in this study are reversible mixed‐type inhibitors of acetylcholinesterase and butyrylcholinesterase. The pK_i values for acetylcholinesterase and butyrylcholinesterase inhibition by the cardiovascular drugs are linearly correlated with the molecular weights of the drugs with the slopes of 0.005 and 0.0021, respectively. Therefore, van der Waals' interactions between acetylcholinesterase and the cardiovascular drugs are stronger than those between butyrylcholinesterase and the drugs. This is probably due to a smaller active site gorge and a more significant peripheral anionic substrate binding site of acetylcholinesterase than those of butyrylcholinesterase. The fact that the pK_i values for both butyrylcholinesterase and acetylcholinesterase inhibition by the cardiovascular drugs are linearly correlated with each other suggests that both enzyme inhibition reactions proceed via a common mechanism. Furthermore, amlodipine besylate may be useful in Alzheimer's disease treatment similar to donepezil.     

Acetylcholinesterase complexed with bivalent ligands related to huperzine a: experimental evidence for species-dependent protein-ligand complementarity

Acetylcholinesterase (AChE) inhibitors improve the cognitive abilities of Alzheimer patients. (-)-Huperzine A [(-)-HupA], an alkaloid isolated from the club moss, Huperzia serrata, is one such inhibitor, but the search for more potent and selective drugs continues. Recently, alkylene-linked dimers of 5-amino-5,6,7,8-tetrahydroquinolinone (hupyridone, 1a), a fragment of HupA, were shown to serve as more potent inhibitors of AChE than (-)-HupA and monomeric 1a. We soaked two such dimers, (S,S)-(-)-bis(10)-hupyridone [(S,S)-(-)-2a] and (S,S)-(-)-bis(12)-hupyridone [(S,S)-(-)-2b] containing, respectively, 10 and 12 methylenes in the spacer, into trigonal TcAChE crystals, and solved the X-ray structures of the resulting complexes using the difference Fourier technique, both to 2.15 A resolution. The structures revealed one HupA-like 1a unit bound to the "anionic" subsite of the active-site, near the bottom of the active-site gorge, adjacent to Trp84, as seen for the TcAChE/(-)-HupA complex, and the second 1a unit near Trp279 in the "peripheral" anionic site at the top of the gorge, both bivalent molecules thus spanning the active-site gorge. The results confirm that the increased affinity of the dimeric HupA analogues for AChE is conferred by binding to the two "anionic" sites of the enzyme. Inhibition data show that (-)-2a binds to TcAChE approximately 6-7- and > 170-fold more tightly than (-)-2b and (-)-HupA, respectively. In contrast, previous data for rat AChE show that (-)-2b binds approximately 3- and approximately 2-fold more tightly than (-)-2a and (-)-HupA, respectively. Structural comparison of TcAChE with rat AChE, as represented by the closely related mouse AChE structure (1maa.pdb), reveals a narrower gorge for rat AChE, a perpendicular alignment of the Tyr337 ring to the gorge axis, and its conformational rigidity, as a result of hydrogen bonding between its hydroxyl group and that of Tyr341, relative to TcAChE Phe330. These structural differences in the active-site gorge explain the switch in inhibitory potency of (-)-2a and 2b and the larger dimer/(-)-HupA potency ratios observed for TcAChE relative to rat AChE. The results offer new insights into factors affecting protein-ligand complementarity within the gorge and should assist the further development of improved AChE inhibitors.     

Central cholinergic agents. I. Potent acetylcholinesterase inhibitors, 2-[omega-[N-alkyl-N-(omega-phenyl-alkyl)amino]alkyl]-1H-isoindole- 1,3(2H)-diones, based on a new hypothesis of the enzyme's active site.

It has been suggested that the active site of acetylcholinesterase contains a hydrophobic binding site (HBS-1), which is closely adjacent to both the anionic and the esteratic sites. In this paper, we assumed that there exists another hydrophobic binding site (HBS-2), some distance removed from the anionic site. On this assumption, a new working hypothesis was proposed for the design of acetylcholinesterase inhibitors. A series of 2-[omega-[N-alkyl-N-(omega-phenyl-alkyl)amino]alkyl]-1H- isoindole-1,3(2H)-diones was designed based on this hypothesis and tested for its inhibitory activities on acetylcholinesterase. Some in this series were revealed to be more potent than physostigmine. Optimum activity was found to be associated with a five carbon chain length separating the benzylamino group from the 1H-isoindole-1,3(2H)-dione (phthalimide) moiety. Quantitative study of substitution effect on the phthalimide moiety revealed that hydrophilic and electron-withdrawing groups enhance the activity.     

The role of the peripheral anionic site and cation-pi interactions in the ligand penetration of the human AChE gorge

We study the ligand (tetramethylammonium) recognition by the peripheral anionic site and its penetration of the human AChE gorge by using atomistic molecular dynamics simulations and our recently developed metadynamics method. The role of both the peripheral anionic site and the formation of cation-pi interactions in the ligand entrance are clearly shown. In particular, a simulation with the W286A mutant shows the fundamental role of this residue in anchoring the ligand at the peripheral anionic site of the enzyme and in positioning it prior to the gorge entrance. Once the ligand is properly oriented, the formation of specific and synchronized cation-pi interactions with W86, F295, and Y341 enables the gorge penetration. Eventually, the ligand is stabilized in a free energy basin by means of cation-pi interactions with W86.     

Hydrogen-bonding interactions in the binding of loop 1 of fasciculin 2 to Torpedo californica acetylcholinesterase: a density functional theory study

In the present study, the interactions of model complexes at the interface between loop1 of fasciculin 2 (Fas2) and acetylcholinesterase (AChE) are theoretically explored. Three interaction models based upon the crystal structure of the complex of Fas2 with AChE from Torpedo californica (PDB code ) were fully optimized at the B3LYP/6-311G(d,p) level of theory. The atoms-in-molecules (AIM) approach was employed to characterize the corresponding noncovalent hydrogen bonds through the densities and the Laplacians of electron densities at the bond critical points. The total interaction energy of loop 1 (Fas2) with AChE is predicted to be -99.4 kcal/mol. It is concluded that the Fas2 residue Thr8, which contributes more than half of the total binding energy, plays the most important role among the three binding sites. The energy decomposition results through the Kitaura-Morokuma scheme suggest that the electrostatic term is the major component of the entire interaction energy. The positive cooperativity effect revealed in the Thr8(F)-related models was confirmed through the geometry characteristics, AIM results, and the energy decomposition analysis.     

Prediction of the binding sites of huperzine A in acetylcholinesterase by docking studies

Summary We have performed docking studies with the SYSDOC program on acetylcholinesterase (AChE) to predict the binding sites in AChE of huperzine A (HA), which is a potent and selective, reversible inhibitor of AChE. The unique aspects of our docking studies include the following: (i) Molecular flexibility of the guest and the host is taken into account, which permits both to change their conformations upon binding. (ii) The binding energy is evaluated by a sum of energies of steric, electrostatic and hydrogen bonding interactions. In the energy calculation no grid approximation is used, and all hydrogen atoms of the system are treated explicitly. (iii) The energy of cation- interactions between the guest and the host, which is important in the binding of AChE, is included in the calculated binding energy. (iv) Docking is performed in all regions of the host's binding cavity. Based on our docking studies and the pharmacological results reported for HA and its analogs, we predict that HA binds to the bottom of the binding cavity of AChE (the gorge) with its ammonium group interacting with Trp⁸⁴, Phe³³⁰, Glu¹⁹⁹ and Asp⁷² (catalytic site). At the the opening of the gorge with its ammonium group partially interacting with Trp²⁷⁹ (peripheral site). At the catalytic site, three partially overlapping subsites of HA were identified which might provide a dynamic view of binding of HA to the catalytic site.