Three-dimensional structures of Drosophila melanogaster acetylcholinesterase and of its complexes with two potent inhibitors

We have crystallized Drosophila melanogaster acetylcholinesterase and solved the structure of the native enzyme and of its complexes with two potent reversible inhibitors, 1,2,3,4-tetrahydro-N-(phenylmethyl)-9-acridinamine and 1,2,3,4-tetrahydro-N-(3-iodophenyl-methyl)-9-acridinamine--all three at 2.7 A resolution. The refined structure of D. melanogaster acetylcholinesterase is similar to that of vertebrate acetylcholinesterases, for example, human, mouse, and fish, in its overall fold, charge distribution, and deep active-site gorge, but some of the surface loops deviate by up to 8 A from their position in the vertebrate structures, and the C-terminal helix is shifted substantially. The active-site gorge of the insect enzyme is significantly narrower than that of Torpedo californica AChE, and its trajectory is shifted several angstroms. The volume of the lower part of the gorge of the insect enzyme is approximately 50% of that of the vertebrate enzyme. Upon binding of either of the two inhibitors, nine aromatic side chains within the active-site gorge change their conformation so as to interact with the inhibitors. Some differences in activity and specificity between the insect and vertebrate enzymes can be explained by comparison of their three-dimensional structures.     

Flexibility of aromatic residues in the active-site gorge of acetylcholinesterase: X-ray versus molecular dynamics

The high aromatic content of the deep and narrow active-site gorge of acetylcholinesterase (AChE) is a remarkable feature of this enzyme. Here, we analyze conformational flexibility of the side chains of the 14 conserved aromatic residues in the active-site gorge of Torpedo californica AChE based on the 47 three-dimensional crystal structures available for the native enzyme, and for its complexes and conjugates, and on a 20-ns molecular dynamics (MD) trajectory of the native enzyme. The degree of flexibility of these 14 aromatic side chains is diverse. Although the side-chain conformations of F330 and W279 are both very flexible, the side-chain conformations of F120, W233, W432, Y70, Y121, F288, F290 and F331 appear to be fixed. Residues located on, or adjacent to, the Ω-loop (C67-C94), namely W84, Y130, Y442, and Y334, display different flexibilities in the MD simulations and in the crystal structures. An important outcome of our study is that the majority of the side-chain conformations observed in the 47 Torpedo californica AChE crystal structures are faithfully reproduced by the MD simulation on the native enzyme. Thus, the protein can assume these conformations even in the absence of the ligand that permitted their experimental detection. These observations are pertinent to structure-based drug design.     

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.     

Monitoring the reaction of carbachol with acetylcholinesterase by thioflavin T fluorescence and acetylthiocholine hydrolysis

Acetylcholinesterase (AChE) contains a narrow and deep active site gorge with two sites of ligand binding, an acylation site (or A-site) at the base of the gorge and a peripheral site (or P-site) near the gorge entrance. The P-site contributes to catalytic efficiency by transiently binding substrates on their way to the acylation site, where a short-lived acyl enzyme intermediate is produced. Carbamates are very poor substrates that, like other AChE substrates, form an initial enzyme-substrate complex and proceed to an acylated enzyme intermediate which is then hydrolyzed. However, the hydrolysis of the carbamoylated enzyme is slow enough to resolve the acylation and deacylation steps on the catalytic pathway. Here we show that the reaction of carbachol (carbamoylcholine) with AChE can be monitored both with acetylthiocholine as a reporter substrate and with thioflavin T as a fluorescent reporter group. The fluorescence of thioflavin T is strongly enhanced when it binds to the P-site of AChE, and this fluorescence is partially quenched when a second ligand binds to the A-site to form a ternary complex. These fluorescence changes allow not only the monitoring of the course of the carbamoylation reaction but also the determination of carbachol affinities for the A- and P-sites.     

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.     

A new proposal for urease mechanism based on the crystal structures of the native and inhibited enzyme from Bacillus pasteurii: why urea hydrolysis costs two nickels.

BACKGROUND: Urease catalyzes the hydrolysis of urea, the final step of organic nitrogen mineralization, using a bimetallic nickel centre. The role of the active site metal ions and amino acid residues has not been elucidated to date. Many pathologies are associated with the activity of ureolytic bacteria, and the efficiency of soil nitrogen fertilization with urea is severely decreased by urease activity. Therefore, the development of urease inhibitors would lead to a reduction of environmental pollution, to enhanced efficiency of nitrogen uptake by plants, and to improved therapeutic strategies for treatment of infections due to ureolytic bacteria. Structure-based design of urease inhibitors would require knowledge of the enzyme mechanism at the molecular level. RESULTS: The structures of native and inhibited urease from Bacillus pasteurii have been determined at a resolution of 2.0 A by synchrotron X-ray cryogenic crystallography. In the native enzyme, the coordination sphere of each of the two nickel ions is completed by a water molecule and a bridging hydroxide. A fourth water molecule completes a tetrahedral cluster of solvent molecules. The enzyme crystallized in the presence of phenylphosphorodiamidate contains the tetrahedral transition-state analogue diamidophosphoric acid, bound to the two nickel ions in an unprecedented mode. Comparison of the native and inhibited structures reveals two distinct conformations of the flap lining the active-site cavity. CONCLUSIONS: The mode of binding of the inhibitor, and a comparison between the native and inhibited urease structures, indicate a novel mechanism for enzymatic urea hydrolysis which reconciles the available structural and biochemical data.     

Crystallographic comparison of manganese- and iron-dependent homoprotocatechuate 2,3-dioxygenases

The X-ray crystal structures of homoprotocatechuate 2,3-dioxygenases isolated from Arthrobacter globiformis and Brevibacterium fuscum have been determined to high resolution. These enzymes exhibit 83% sequence identity, yet their activities depend on different transition metals, Mn2+ and Fe2+, respectively. The structures allow the origins of metal ion selectivity and aspects of the molecular mechanism to be examined in detail. The homotetrameric enzymes belong to the type I family of extradiol dioxygenases (vicinal oxygen chelate superfamily); each monomer has four betaalphabetabetabeta modules forming two structurally homologous N-terminal and C-terminal barrel-shaped domains. The active-site metal is located in the C-terminal barrel and is ligated by two equatorial ligands, H214NE1 and E267OE1; one axial ligand, H155NE1; and two to three water molecules. The first and second coordination spheres of these enzymes are virtually identical (root mean square difference over all atoms, 0.19 A), suggesting that the metal selectivity must be due to changes at a significant distance from the metal and/or changes that occur during folding. The substrate (2,3-dihydroxyphenylacetate [HPCA]) chelates the metal asymmetrically at sites trans to the two imidazole ligands and interacts with a unique, mobile C-terminal loop. The loop closes over the bound substrate, presumably to seal the active site as the oxygen activation process commences. An "open" coordination site trans to E267 is the likely binding site for O2. The geometry of the enzyme-substrate complexes suggests that if a transiently formed metal-superoxide complex attacks the substrate without dissociation from the metal, it must do so at the C-3 position. Second-sphere active-site residues that are positioned to interact with the HPCA and/or bound O2 during catalysis are identified and discussed in the context of current mechanistic hypotheses.     

Inhibition of Drosophila melanogaster acetylcholinesterase by high concentrations of substrate.

Acetylcholine hydrolysis by acetylcholinesterase is inhibited at high substrate concentrations. To determine the residues involved in this phenomenon, we have mutated most of the residues lining the active-site gorge but mutating these did not completely eliminate hydrolysis. Thus, we analyzed the effect of a nonhydrolysable substrate analogue on substrate hydrolysis and on reactivation of an analogue of the acetylenzyme. Analyses of various models led us to propose the following sequence of events: the substrate initially binds at the rim of the active-site gorge and then slides down to the bottom of the gorge where it is hydrolyzed. Another substrate molecule can bind to the peripheral site: (a) when the choline is still inside the gorge - it will thereby hinder its exit; (b) after choline has dissociated but before deacetylation occurs - binding at the peripheral site increases deacetylation rate but (c) if a substrate molecule bound to the peripheral site slides down to the bottom of the active-site before the catalytic serine is deacetylated, its new position will prevent the approach of water, thus blocking deacetylation.     

Structural Insights into the Substrate Specificity of (S)-Ureidoglycolate Amidohydrolase and Its Comparison with Allantoate Amidohydrolase

In plants, the ureide pathway is a metabolic route that converts the ring nitrogen atoms of purine into ammonia via sequential enzymatic reactions, playing an important role in nitrogen recovery. In the final step of the pathway, (S)-ureidoglycolate amidohydrolase (UAH) catalyzes the conversion of (S)-ureidoglycolate into glyoxylate and releases two molecules of ammonia as by-products. UAH is homologous in structure and sequence with allantoate amidohydrolase (AAH), an upstream enzyme in the pathway with a similar function as that of an amidase but with a different substrate. Both enzymes exhibit strict substrate specificity and catalyze reactions in a concerted manner, resulting in purine degradation. Here, we report three crystal structures of Arabidopsis thaliana UAH (bound with substrate, reaction intermediate, and product) and a structure of Escherichia coli AAH complexed with allantoate. Structural analyses of UAH revealed a distinct binding mode for each ligand in a bimetal reaction center with the active site in a closed conformation. The ligand directly participates in the coordination shell of two metal ions and is stabilized by the surrounding residues. In contrast, AAH, which exhibits a substrate-binding site similar to that of UAH, requires a larger active site due to the additional ureido group in allantoate. Structural analyses and mutagenesis revealed that both enzymes undergo an open-to-closed conformational transition in response to ligand binding and that the active-site size and the interaction environment in UAH and AAH are determinants of the substrate specificities of these two structurally homologous enzymes.     

Hydration energy landscape of the active site cavity in cytochrome P450cam

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein–ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol⁻¹) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable. Proteins 32:381–396, 1998. © 1998 Wiley-Liss, Inc.     

4-Aryl-4-oxo-N-phenyl-2-aminylbutyramides as acetyl- and butyrylcholinesterase inhibitors. Preparation,anticholinesterase activity,docking study,and 3D structure–activity relationship based on molecular interaction fields

Synthesis and anticholinesterase activity of 4-aryl-4-oxo-N-phenyl-2-aminylbutyramides, novel class of reversible, moderately potent cholinesterase inhibitors, are reported. Simple substituent variation on aroyl moiety changes anti-AChE activity for two orders of magnitude; also substitution and type of hetero(ali)cycle in position 2 of butanoic moiety govern AChE/BChE selectivity. The most potent compounds showed mixed-type inhibition, indicating their binding to free enzyme and enzyme–substrate complex. Alignment-independent 3D QSAR study on reported compounds, and compounds having similar potencies obtained from the literature, confirmed that alkyl substitution on aroyl moiety of molecules is requisite for inhibition activity. The presence of hydrophobic moiety at close distance from hydrogen bond acceptor has favorable influence on inhibition potency. Docking studies show that compounds probably bind in the middle of the AChE active site gorge, but are buried deeper inside BChE active site gorge, as a consequence of larger BChE gorge void.     

X-ray structures of small ligand-FKBP complexes provide an estimate for hydrophobic interaction energies

A new crystal form of native FK506 binding protein (FKBP) has been obtained which has proved useful in ligand binding studies. Three different small molecule ligand complexes and the native enzyme have been determined at higher resolution than 2.0 A. Dissociation constants of the related small molecule ligands vary from 20 mM for dimethylsulphoxide to 200 microM for tetrahydrothiophene 1-oxide. Comparison of the four available crystal structures shows that the protein structures are identical to within experimental error, but there are differences in the water structure in the active site. Analysis of the calculated buried surface areas of these related ligands provides an estimated van der Waals contribution to the binding energy of -0.5 kJ/A(2) for non-polar interactions between ligand and protein.     

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     

Structural insights into substrate traffic and inhibition in acetylcholinesterase

Acetylcholinesterase (AChE) terminates nerve-impulse transmission at cholinergic synapses by rapid hydrolysis of the neurotransmitter, acetylcholine. Substrate traffic in AChE involves at least two binding sites, the catalytic and peripheral anionic sites, which have been suggested to be allosterically related and involved in substrate inhibition. Here, we present the crystal structures of Torpedo californica AChE complexed with the substrate acetylthiocholine, the product thiocholine and a nonhydrolysable substrate analogue. These structures provide a series of static snapshots of the substrate en route to the active site and identify, for the first time, binding of substrate and product at both the peripheral and active sites. Furthermore, they provide structural insight into substrate inhibition in AChE at two different substrate concentrations. Our structural data indicate that substrate inhibition at moderate substrate concentration is due to choline exit being hindered by a substrate molecule bound at the peripheral site. At the higher concentration, substrate inhibition arises from prevention of exit of acetate due to binding of two substrate molecules within the active-site gorge.     

pK values for active site residues of carboxypeptidase A

The phenolic group of active site residue Tyr-248 in carboxypeptidase A has a pKa value of 10.06, as determined from the pH dependence of its rate of nitration by tetranitromethane. The decrease in enzyme activity (kcat/Km) in alkaline solution, characterized by a pKa value of approximately 9.0 (for cobalt carboxypeptidase A), is associated with the protonation state of an imidazole ligand of the active-site metal ion, as indicated by a selective pH dependence of the 1H NMR spectrum of the enzyme. Inhibition of the cobalt-substituted enzyme by 2-(1-carboxy-2-phenylethyl)phenol and its 4,6-dichloro- and 4-phenylazo-derivatives confirms that the decrease in enzyme activity (kcat/Km) in acidic solution, characterized by a pKa value of 5.8, is due to the protonation state of a water molecule bound to the active-site metal ion in the absence of substrate. Changes in the coordination number of the active-site metal ion are seen in its visible absorption spectrum as a consequence of binding of the phenolic inhibitors. Conventional concepts regarding the mechanisms of the enzyme are brought into question.     

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 NH₂ 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.     

Magnetic resonance investigation of ionizable residues at the active site of thermolysin.

The details of the pH dependence of the thermodynamic and magnetic interactions of the active-site region of thermolysin in which manganese has replaced the active-site zinc atom and the inhibitor N-trifluoroacetyl-D-phenylalanine have been examined. These show a number of ionizable groups in the active-site region. A cooperative displacement of manganese at the catalytic site is observed as pH is lowered. This appears to be the result of the protonation of histidine-142 and -146 which act as metal ligands. The metal is 50% displaced at pH 6.0. At higher pH values, the environment of the bound manganese changes as a result of the ionization of at least two groups of approximate pKa = 8.5 and 9.5. These values are assigned to tyrosine-157 and to the water molecule which acts as a metal ligand at the active site. The binding behavior of the inhibitor strongly suggests that two molecules of inhibitor bind to the enzyme. The weaker site is competitive with the synthetic substrate FAGLA (furylacryloylglycyl-leucinamide), while the strong site has no effect on FAGLA hydrolysis. This second site is in the vicinity of the active site with a distance of 8 A or less between the trifluoromethyl group and manganese bound at the active site.     

Structural analyses of covalent enzyme-substrate analog complexes reveal strengths and limitations of de novo enzyme design

We report the cocrystal structures of a computationally designed and experimentally optimized retro-aldol enzyme with covalently bound substrate analogs. The structure with a covalently bound mechanism-based inhibitor is similar to, but not identical with, the design model, with an RMSD of 1.4 Å over active-site residues and equivalent substrate atoms. As in the design model, the binding pocket orients the substrate through hydrophobic interactions with the naphthyl moiety such that the oxygen atoms analogous to the carbinolamine and β-hydroxyl oxygens are positioned near a network of bound waters. However, there are differences between the design model and the structure: the orientation of the naphthyl group and the conformation of the catalytic lysine are slightly different; the bound water network appears to be more extensive; and the bound substrate analog exhibits more conformational heterogeneity than typical native enzyme–inhibitor complexes. Alanine scanning of the active-site residues shows that both the catalytic lysine and the residues around the binding pocket for the substrate naphthyl group make critical contributions to catalysis. Mutating the set of water-coordinating residues also significantly reduces catalytic activity. The crystal structure of the enzyme with a smaller substrate analog that lacks naphthyl ring shows the catalytic lysine to be more flexible than in the naphthyl–substrate complex; increased preorganization of the active site would likely improve catalysis. The covalently bound complex structures and mutagenesis data highlight the strengths and weaknesses of the de novo enzyme design strategy.     

Human glutathione-dependent formaldehyde dehydrogenase. Structures of apo,binary, and inhibitory ternary complexes

The human glutathione-dependent formaldehyde dehydrogenase is unique among the structurally studied members of the alcohol dehydrogenase family in that it follows a random bi bi kinetic mechanism. The structures of an apo form of the enzyme, a binary complex with substrate 12-hydroxydodecanoic acid, and a ternary complex with NAD+ and the inhibitor dodecanoic acid were determined at 2.0, 2.3, and 2.3 A resolution by X-ray crystallography using the anomalous diffraction signal of zinc. The structures of the enzyme and its binary complex with the primary alcohol substrate, 12-hydroxydodecanoic acid, and the previously reported binary complex with the coenzyme show that the binding of the first substrate (alcohol or coenzyme) causes only minor changes to the overall structure of the enzyme. This is consistent with the random mechanism of the enzyme where either of the substrates binds to the free enzyme. The catalytic-domain position in these structures is intermediate to the "closed" and "open" conformations observed in class I alcohol dehydrogenases. More importantly, two different tetrahedral coordination environments of the active site zinc are observed in these structures. In the apoenzyme, the active site zinc is coordinated to Cys44, His66 and Cys173, and a water molecule. In the inhibitor complex, the coordination environment involves Glu67 instead of the solvent water molecule. The coordination environment involving Glu67 as the fourth ligand likely represents an intermediate step during ligand exchange at the active site zinc. These observations provide new insight into metal-assisted catalysis and substrate binding in glutathione-dependent formaldehyde dehydrogenase.