Crystal structure of alpha-amino-beta-carboxymuconate-epsilon-semialdehyde decarboxylase: insight into the active site and catalytic mechanism of a novel decarboxylation reaction

Alpha-amino-beta-carboxymuconate-epsilon-semialdehyde decarboxylase (ACMSD) is a widespread enzyme found in many bacterial species and all currently sequenced eukaryotic organisms. It occupies a key position at the branching point of two metabolic pathways: the tryptophan to quinolinate pathway and the bacterial 2-nitrobenzoic acid degradation pathway. The activity of ACMSD determines whether the metabolites in both pathways are converted to quinolinic acid for NAD biosynthesis or to acetyl-CoA for the citric acid cycle. Here we report the first high-resolution crystal structure of ACMSD from Pseudomonas fluorescens which validates our previous predictions that this enzyme is a member of the metal-dependent amidohydrolase superfamily of the (beta/alpha)(8) TIM barrel fold. The structure of the enzyme in its native form, determined at 1.65 A resolution, reveals the precise spatial arrangement of the active site metal center and identifies a potential substrate-binding pocket. The identity of the native active site metal was determined to be Zn. Also determined was the structure of the enzyme complexed with cobalt at 2.50 A resolution. The hydrogen bonding network around the metal center suggests that Arg51 and His228 may play important roles in catalysis. The metal center configuration of PfACMSD is very similar to that of Zn-dependent adenosine deaminase and Fe-dependent cytosine deaminase, suggesting that ACMSD may share certain similarities in its catalytic mechanism with these enzymes. These data enable us to propose possible catalytic mechanisms for ACMSD which appear to be unprecedented among all currently characterized decarboxylases.     

Crystal structures of creatininase reveal the substrate binding site and provide an insight into the catalytic mechanism

Creatininase from Pseudomonas putida is a member of the urease-related amidohydrolase superfamily. The crystal structure of the Mn-activated enzyme has been solved by the single isomorphous replacement method at 1.8A resolution. The structures of the native creatininase and the Mn-activated creatininase-creatine complex have been determined by a difference Fourier method at 1.85 A and 1.6 A resolution, respectively. We found the disc-shaped hexamer to be roughly 100 A in diameter and 50 A in thickness and arranged as a trimer of dimers with 32 (D3) point group symmetry. The enzyme is a typical Zn2+ enzyme with a binuclear metal center (metal1 and metal2). Atomic absorption spectrometry and X-ray crystallography revealed that Zn2+ at metal1 (Zn1) was easily replaced with Mn2+ (Mn1). In the case of the Mn-activated enzyme, metal1 (Mn1) has a square-pyramidal geometry bound to three protein ligands of Glu34, Asp45, and His120 and two water molecules. Metal2 (Zn2) has a well-ordered tetrahedral geometry bound to the three protein ligands of His36, Asp45, and Glu183 and a water molecule. The crystal structure of the Mn-activated creatininase-creatine complex, which is the first structure as the enzyme-substrate/inhibitor complex of creatininase, reveals that significant conformation changes occur at the flap (between the alpha5 helix and the alpha6 helix) of the active site and the creatine is accommodated in a hydrophobic pocket consisting of Trp174, Trp154, Tyr121, Phe182, Tyr153, and Gly119. The high-resolution crystal structure of the creatininase-creatine complex enables us to identify two water molecules (Wat1 and Wat2) that are possibly essential for the catalytic mechanism of the enzyme. The structure and proposed catalytic mechanism of the creatininase are different from those of urease-related amidohydrolase superfamily enzymes. We propose a new two-step catalytic mechanism possibly common to creatininases in which the Wat1 acts as the attacking nucleophile in the water-adding step and the Wat2 acts as the catalytic acid in the ring-opening step.     

Identification of an active site ligand for a group I ribozyme catalytic metal ion

The transition state of the group I intron self-splicing reaction is stabilized by three metal ions. The functional groups within the intron substrates (guanosine and an oligoribonucleotide mimic of the 5'-exon) that coordinate these metal ions have been systematically defined through a series of metal ion specificity switch experiments. In contrast, the catalytic metal ligands within the ribozyme active site are unknown. In an effort to identify them, stereospecific (R(P) or S(P)) single-site phosphorothioate substitutions were introduced at five phosphates predicted to be in the vicinity of the catalytic center (A207, C208, A304, U305, and A306) within the Tetrahymena intron. Of the 10 ribozymes that were studied, four phosphorothioate substitutions (A207 S(P), C208 S(P), A306 R(P), and A306 S(P)) exhibited a significant reduction in the cleavage rate. Only the effect of the C208 S(P) phosphorothioate substitution could be significantly rescued by the addition of a thiophilic metal ion, either Mn(2+) or Zn(2+), when tested with an all-oxy substrate. The effect was not rescued with Cd(2+). To determine if one of the catalytic metal ions is coordinated to the C208 pro-S(P) oxygen, the phosphorothioate-substituted ribozymes were also assayed using oligonucleotide substrates with a 3'-phosphorothiolate or an S(P) phosphorothioate substitution at the scissile phosphate. This resulted in a second metal specificity switch, in that Mn(2+) or Zn(2+) no longer rescued the C208 S(P) ribozyme, but Cd(2+) provided efficient rescue in the context of either sulfur-containing substrate. The 3'-oxygen and the pro-S(P) oxygen of the scissile phosphate are both known to coordinate the same metal ion, M(A), which stabilizes the negative charge on the leaving group 3'-oxygen in the transition state. Taken together, these data suggest that metal M(A) is coordinated to the C208 pro-S(P) phosphate oxygen, which constitutes the first functional link between a specific catalytic metal ion and a particular functional group within the group I ribozyme active site.     

Evidence for a polynuclear metal ion binding site in the catalytic domain of ribonuclease P RNA

Interactions with divalent metal ions are essential for the folding and function of the catalytic RNA component of the tRNA processing enzyme ribonuclease P (RNase P RNA). However, the number and location of specific metal ion interactions in this large, highly structured RNA are poorly understood. Using atomic mutagenesis and quantitative analysis of thiophilic metal ion rescue we provide evidence for metal ion interactions at the pro-R(P) and pro-S(P) non-bridging phosphate oxygens at nucleotide A67 in the universally conserved helix P4. Moreover, second-site modifications within helix P4 and the adjacent single stranded region (J3/4) provide the first evidence for metal ion interactions with nucleotide base functional groups in RNase P RNA and reveal the presence of an additional metal ion important for catalytic function. Together, these data are consistent with a cluster of metal ion interactions in the P1-P4 multi-helix junction that defines the catalytic core of the RNase P ribozyme.     

Crystal structure of the leadzyme at 1.8 A resolution: metal ion binding and the implications for catalytic mechanism and allo site ion regulation

The leadzyme is a small ribozyme, derived from in vitro selection, which catalyzes site specific, Pb(2+)-dependent RNA cleavage. Pb(2+) is required for activity; Mg(2+) inhibits activity, while many divalent and trivalent ions enhance it. The leadzyme structure consists of an RNA duplex interrupted by a trinucleotide bulge. Here, crystal structures determined to 1.8 A resolution, both with Mg(2+) as the sole divalent counterion and with Mg(2+) and Sr(2+) (which mimics Pb(2+) with respect to binding but not catalysis), reveal the metal ion interactions with both the ground state and precatalytic conformations of the leadzyme. Mg(H(2)O)(6)(2+) ions bridge complementary strands of the duplex at multiple locations by binding tandem purines of one RNA strand in the major groove. At one site, Mg(H(2)O)(6)(2+) ligates the phosphodiester backbone of the trinucleotide bulge in the ground state conformation, but not in the precatalytic conformation, suggesting (a) Mg(2+) may inhibit leadzyme activity by stabilizing the ground state and (b) metal ions which displace Mg(2+) from this site may activate the leadzyme. Binding of Sr(2+) to the presumed catalytic Pb(2+) site in the precatalytic leadzyme induces local structural changes in a manner that would facilitate alignment of the catalytic ribose 2'-hydroxyl with the scissile bond for cleavage. These data support a model wherein binding of a catalytic ion to a precatalytic conformation of the leadzyme, in conjunction with the flexibility of the trinucleotide bulge, may facilitate structural rearrangements around the scissle phosphodiester bond favoring configurations that allow bond cleavage.     

Ca2+ binding in the active site of HincII: implications for the catalytic mechanism

The 2.8 A crystal structure of the type II restriction endonuclease HincII bound to Ca(2+) and cognate DNA containing GTCGAC is presented. The DNA is uncleaved, and one calcium ion is bound per active site, in a position previously described as site I in the related blunt cutting type II restriction endonuclease EcoRV [Horton, N. C., Newberry, K. J., and Perona, J. J. (1998) Proc. Natl. Acad. Sci. U.S.A. 95 (23), 13489-13494], as well as that found in other related enzymes. Unlike the site I metal in EcoRV, but similar to that of PvuII, NgoMIV, BamHI, BglII, and BglI, the observed calcium cation is directly ligated to the pro-S(p) oxygen of the scissile phosphate. A calcium ion-ligated water molecule is well positioned to act as the nucleophile in the phosphodiester bond cleavage reaction, and is within hydrogen bonding distance of the conserved active site lysine (Lys 129), as well as the pro-R(p) oxygen of the phosphate group 3' of the scissile phosphate, suggesting possible roles for these groups in the catalytic mechanism. Kinetic data consistent with an important role for the 3'-phosphate group in DNA cleavage by HincII are presented. The previously observed sodium ion [Horton, N. C., Dorner, L. F., and Perona, J. J. (2002) Nat. Struct. Biol. 9, 42-47] persists in the active sites of the Ca(2+)-bound structure; however, kinetic data show little effect on the single-turnover rate of DNA cleavage in the absence of Na(+) ions.     

An inhibitor binding pocket distinct from the catalytic active site on human beta-APP cleaving enzyme

Beta-APP cleaving enzyme (BACE) is responsible for the first of two proteolytic cleavages of the APP protein that together lead to the generation of the Alzheimer's disease-associated Abeta peptide. It is widely believed that halting the production of Abeta peptide, by inhibition of BACE, is an attractive therapeutic modality for the treatment of Alzheimer's disease. BACE is an aspartyl protease, and there is significant effort in the pharmaceutical community to apply traditional design methods to the development of active site-directed inhibitors of this enzyme. We report here the discovery of a ligand binding pocket within the catalytic domain of BACE that is distinct from the enzymatic active site (i.e., an exosite). Peptides, initially identified from combinatorial phage peptide libraries, contain the sequence YPYF(I/L)P(L/I) and bind specifically to this exosite, even in the presence of saturating concentrations of active site-directed inhibitors. Binding of peptides to the BACE exosite leads to a concentration-dependent inhibition of proteolysis for APP-related, protein-based substrates of BACE. The discovery of this exosite opens new opportunities for the identification and development of novel and potentially selective small molecule inhibitors of BACE that act through exosite, rather than active site, binding interactions.     

ATP conformations and ion binding modes in the active site of anthrax edema factor: A computational analysis

The Edema Factor (EF), one of the virulence factors of anthrax, is an adenylyl cyclase that promotes the overproduction of cyclic‐AMP (cAMP) from ATP, and therefore perturbs cell signaling. Crystallographic structures of EF bound to ATP analogs and reaction products, cyclic‐AMP, and Pyrophosphate (PPi), revealed different substrate conformations and catalytic‐cation binding modes, one or two cations being observed in the active site. To shed light into the biological significance of these crystallographic structures, the energetics, geometry, and dynamics of the active site are analyzed using molecular dynamics simulations. The ATP conformation observed in the one‐metal‐ion structure allows stronger interactions with the catalytic ion, and ATP is more restrained than in the structure containing two Mg²⁺ ions. Therefore, we propose that the conformation observed in the one‐ion crystal structure is a more probable starting point for the reaction. The simulations also suggest that a C3′‐endo sugar pucker facilitates nucleophilic attack. Additionally, the two‐cation binding mode restrains the mobility of the reaction products, and thus their tendency to dissociate. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural insights into the mechanism of nuclease A, a betabeta alpha metal nuclease from Anabaena

Nuclease A (NucA) is a nonspecific endonuclease from Anabaena sp. capable of degrading single- and double-stranded DNA and RNA in the presence of divalent metal ions. We have determined the structure of the delta(2-24),D121A mutant of NucA in the presence of Zn2+ and Mn2+ (PDB code 1ZM8). The mutations were introduced to remove the N-terminal signal peptide and to reduce the activity of the nonspecific nuclease, thereby reducing its toxicity to the Escherichia coli expression system. NucA contains a betabeta alpha metal finger motif and a hydrated Mn2+ ion at the active site. Unexpectedly, NucA was found to contain additional metal binding sites approximately 26 A apart from the catalytic metal binding site. A structural comparison between NucA and the closest analog for which structural data exist, the Serratia nuclease, indicates several interesting differences. First, NucA is a monomer rather than a dimer. Second, there is an unexpected structural homology between the N-terminal segments despite a poorly conserved sequence, which in Serratia includes a cysteine bridge thought to play a regulatory role. In addition, although a sequence alignment had suggested that NucA lacks a proposed catalytic residue corresponding to Arg57 in Serratia, the structure determined here indicates that Arg93 in NucA is positioned to fulfill this role. Based on comparison with DNA-bound nuclease structures of the betabeta alpha metal finger nuclease family and available mutational data on NucA, we propose that His124 acts as a catalytic base, and Arg93 participates in the catalysis possibly through stabilization of the transition state.     

Structural insight into poly(A) binding and catalytic mechanism of human PARN

Poly(A)-specific ribonuclease (PARN) is a processive, poly(A)-specific 3' exoribonuclease. The crystal structure of C-terminal truncated human PARN determined in two states (free and RNA-bound forms) reveals that PARNn is folded into two domains, an R3H domain and a nuclease domain similar to those of Pop2p and epsilon186. The high similarity of the active site structures of PARNn and epsilon186 suggests that they may have a similar catalytic mechanism. PARNn forms a tight homodimer, with the R3H domain of one subunit partially enclosing the active site of the other subunit and poly(A) bound in a deep cavity of its nuclease domain in a sequence-nonspecific manner. The R3H domain and, possibly, the cap-binding domain are involved in poly(A) binding but these domains alone do not appear to contribute to poly(A) specificity. Mutations disrupting dimerization abolish both the enzymatic and RNA-binding activities, suggesting that the PARN dimer is a structural and functional unit. The cap-binding domain may act in concert with the R3H domain to amplify the processivity of PARN.     

Engineering a catalytic metal binding site into a calcium-independent phosphatidylinositol-specific phospholipase C leads to enhanced stereoselectivity

Eukaryotic phosphatidylinositol-specific phospholipase Cs (PI-PLCs) utilize calcium as a cofactor during catalysis, whereas prokaryotic PI-PLCs use a spatially conserved guanidinium group from Arg69. In this study, we aimed to construct a metal-dependent mutant of a bacterial PI-PLC and characterize the catalytic role of the metal ion, seeking an enhanced understanding of the functional differences between these two positively charged moieties. The following results indicate that a bona fide catalytic metal binding site was created by the single arginine-to-aspartate mutation at position 69: (1) The R69D mutant was activated by Ca(2+), and the activation was specific for R69D, not for other mutants at this position. (2) Titration of R69D with Ca(2+), monitored by (15)N-(1)H HSQC (heteronuclear single quantum coherence) NMR, showed that addition of Ca(2+) to R69D restores the environment of the catalytic site analogous to that attained by the WT enzyme. (3) Upon Ca(2+) activation, the thio effect of the S(P)-isomer of the phosphorothioate analogue (k(O)/k(Sp) = 4.4 x 10(5)) approached a value similar to that of the WT enzyme, suggesting a structural and functional similarity between the two positively charged moieties (Arg69 and Asp69-Ca(2+)). The R(P)-thio effect (k(O)/k(Rp) = 9.4) is smaller than that of the WT enzyme by a factor of 5. Consequently, R69D-Ca(2+) displays higher stereoselectivity (k(Rp)/k(Sp) = 47,000) than WT (k(Rp)/k(Sp) = 7600). (4) Results from additional mutagenesis analyses suggest that the Ca(2+) binding site is comprised of side chains from Asp33, Asp67, Asp69, and Glu117. Our studies provide new insight into the mechanism of metal-dependent and metal-independent PI-PLCs.     

Two distinct binding modes of a protein cofactor with its target RNA

Like most cellular RNA enzymes, the bI5 group I intron requires binding by a protein cofactor to fold correctly. Here, we use single-molecule approaches to monitor the structural dynamics of the bI5 RNA in real time as it assembles with its CBP2 protein cofactor. These experiments show that CBP2 binds to the target RNA in two distinct modes with apparently opposite effects: a "non-specific" mode that forms rapidly and induces large conformational fluctuations in the RNA, and a "specific" mode that forms slowly and stabilizes the native RNA structure. The bI5 RNA folds though multiple pathways toward the native state, typically traversing dynamic intermediate states induced by non-specific binding of CBP2. These results suggest that the protein cofactor-assisted RNA folding involves sequential non-specific and specific protein-RNA interactions. The non-specific interaction potentially increases the local concentration of CBP2 and the number of conformational states accessible to the RNA, which may promote the formation of specific RNA-protein interactions.     

DNA binding and cleavage by the periplasmic nuclease Vvn: a novel structure with a known active site

The Vibrio vulnificus nuclease, Vvn, is a non-specific periplasmic nuclease capable of digesting DNA and RNA. The crystal structure of Vvn and that of Vvn mutant H80A in complex with DNA were resolved at 2.3 A resolution. Vvn has a novel mixed alpha/beta topology containing four disulfide bridges, suggesting that Vvn is not active under reducing conditions in the cytoplasm. The overall structure of Vvn shows no similarity to other endonucleases; however, a known 'betabetaalpha-metal' motif is identified in the central cleft region. The crystal structure of the mutant Vvn-DNA complex demonstrates that Vvn binds mainly at the minor groove of DNA, resulting in duplex bending towards the major groove by approximately 20 degrees. Only the DNA phosphate backbones make hydrogen bonds with Vvn, suggesting a structural basis for its sequence-independent recognition of DNA and RNA. Based on the enzyme-substrate and enzyme-product structures observed in the mutant Vvn-DNA crystals, a catalytic mechanism is proposed. This structural study suggests that Vvn hydrolyzes DNA by a general single-metal ion mechanism, and indicates how non-specific DNA-binding proteins may recognize DNA.     

Burst kinetics and redox transformations of the active site manganese ion in oxalate oxidase: implications for the catalytic mechanism

Oxalate oxidase (EC 1.2.3.4) catalyzes the oxidative cleavage of oxalate to carbon dioxide and hydrogen peroxide. In this study, unusual nonstoichiometric burst kinetics of the steady state reaction were observed and analyzed in detail, revealing that a reversible inactivation process occurs during turnover, associated with a slow isomerization of the substrate complex. We have investigated the underlying molecular mechanism of this kinetic behavior by preparing recombinant barley oxalate oxidase in three distinct oxidation states (Mn(II), Mn(III), and Mn(IV)) and producing a nonglycosylated variant for detailed biochemical and spectroscopic characterization. Surprisingly, the fully reduced Mn(II) form, which represents the majority of the as-isolated native enzyme, lacks oxalate oxidase activity, but the activity is restored by oxidation of the metal center to either Mn(III) or Mn(IV) forms. All three oxidation states appear to interconvert under turnover conditions, and the steady state activity of the enzyme is determined by a balance between activation and inactivation processes. In O(2)-saturated buffer, a turnover-based redox modification of the enzyme forms a novel superoxidized mononuclear Mn(IV) biological complex. An oxalate activation role for the catalytic metal ion is proposed based on these results.     

DNA binding mediates conformational changes and metal ion coordination in the active site of PcrA helicase.

Based upon the crystal structures of PcrA helicase, we have made and characterised mutations in a number of conserved helicase signature motifs around the ATPase active site. We have also determined structures of complexes of wild-type PcrA with ADPNP and of a mutant PcrA complexed with ADPNP and Mn2+. The kinetic and structural data define roles for a number of different residues in and around the ATP binding site. More importantly, our results also show that there are two functionally distinct conformations of ATP in the active site. In one conformation, ATP is hydrolysed poorly whereas in the other (activated) conformation, ATP is hydrolysed much more rapidly. We propose a mechanism to explain how the stimulation of ATPase activity afforded by binding of single-stranded DNA stabilises the activated conformation favouring Mg2+binding and a consequent repositioning of the gamma-phosphate group which promotes ATP hydrolysis. A part of the associated conformational change in the protein forces the side-chain of K37 to vacate the Mg2+binding site, allowing the cation to bind and interact with ATP.     

The first structure of dipeptidyl-peptidase III provides insight into the catalytic mechanism and mode of substrate binding

Dipeptidyl-peptidases III (DPP III) are zinc-dependent enzymes that specifically cleave the first two amino acids from the N terminus of different length peptides. In mammals, DPP III is associated with important physiological functions and is a potential biomarker for certain types of cancer. Here, we present the 1.95-A crystal structure of yeast DPP III representing the prototype for the M49 family of metallopeptidases. It shows a novel fold with two domains forming a wide cleft containing the catalytic metal ion. DPP III exhibits no overall similarity to other metallopeptidases, such as thermolysin and neprilysin, but zinc coordination and catalytically important residues are structurally conserved. Substrate recognition is accomplished by a binding site for the N terminus of the peptide at an appropriate distance from the metal center and by a series of conserved arginine residues anchoring the C termini of different length substrates.     

Model complexes of the active site of galactose oxidase. Effects of the metal ion binding sites

Model compounds of the active site of galactose oxidase have been developed by using new cofactor model ligands, L1H (2-methylthio-4-tert-butyl-6-[{bis(pyridin-2-ylmethyl)amino}methyl]phenol) and L2H (2-methylthio-4-tert-butyl-6-[{bis(6-methylpyridin-2-ylmethyl)amino}methyl]phenol). Treatment of the ligands with copper(II) and zinc(II) perchlorate in the presence of triethylamine followed by anion exchange reaction with NaPF₆ or NaBPh₄ provided the corresponding copper(II) and zinc(II) complexes, the crystal structures of which have been determined by X-ray crystallographic analysis. All the copper(II) and zinc(II) complexes have been isolated as a dimeric form in which the phenolate oxygen of each ligand acts as the bridging ligand to form a rhombic M₂(OAr)₂ core (M=Cu or Zn). The dimeric complexes can be converted into the corresponding monomer complexes by the treatment with exogenous ligand such as acetate ion. The redox potential and the spectroscopic features of the monomer complexes have also been examined. Furthermore, the copper(II)- and zinc(II)-complexes of the phenoxyl radical species of the ligands have been generated in situ by the oxidation of the phenolate complexes with (NH₄)₂[Ce^IV(NO₃)₆] (CAN) in CH₃CN, and their spectroscopic features have been explored. The structures and physicochemical properties of the phenolate and phenoxyl radical complexes of L1 and L2 have been compared to those of the previously reported copper(II) and zinc(II) complexes of L3 (2-methylthio-4-tert-butyl-6-[{bis(2-pyridin-2-ylethyl)amino}methyl]phenol) in order to get insights into the interaction between the metal ions and the organic cofactor moiety.     

Crystal structure of glycerophosphodiester phosphodiesterase (GDPD) from Thermoanaerobacter tengcongensis, a metal ion-dependent enzyme: insight into the catalytic mechanism

Glycerophosphodiester phosphodiesterase (GDPD; EC 3.1.4.46) catalyzes the hydrolysis of a glycerophosphodiester to an alcohol and glycerol 3-phosphate in glycerol metabolism. It has an important role in the synthesis of a variety of products that participate in many biochemical pathways. We report the crystal structure of the Thermoanaerobacter tengcongensis GDPD (ttGDPD) at 1.91 A resolution, with a calcium ion and glycerol as a substrate mimic coordinated at this calcium ion (PDB entry 2pz0). The ttGDPD dimer with an intermolecular disulfide bridge and two hydrogen bonds is considered as the potential functional unit. We used site-directed mutagenesis to characterize ttGDPD as a metal ion-dependent enzyme, identified a cluster of residues involved in substrate binding and the catalytic reaction, and we propose a possible general acid-base catalytic mechanism for ttGDPD. Superposing the active site with the homologous structure GDPD from Agrobacterium tumefaciens (PDB entry 1zcc), which binds a sulfate ion in the active site, the sulfate ion can represent the phosphate moiety of the substrate, simulating the binding mode of the true substrate of GDPD.     

Structural basis of fullerene derivatives as novel potent inhibitors of protein acetylcholinesterase without catalytic active site interaction: insight into the inhibitory mechanism through molecular modeling studies

Abstract

Acetylcholinesterase (AChE) is an important kind of esterase that plays a key biological role in the central and peripheral nervous systems. Recent research has demonstrated that some fullerene derivatives serve as a new nanoscale class of potent inhibitors of AChE, but the specific inhibition mechanism remains unclear. In the present article, several molecular modeling methods, such as molecular docking, molecular dynamics simulations and molecular mechanics/generalized Born surface area calculations, were used for the investigation of the binding mode and inhibition mechanism of fullerene inhibitions for AChE. Results revealed that fullerene inhibitors block the entrance of substrates by binding with the peripheral anionic site (PAS) region. Thus, fullerene derivatives might mainly serve as competitive inhibitors. The interactions of a fullerene backbone with AChE are at the same level in different single side chain systems and seem to be related to the length or aromaticity of the side chain. The inhibitor with multihydroxyl side chains shows an obviously large electrostatic interaction as it forms additional hydrogen bonds with AChE. Moreover, fullerene derivatives might exhibit noncompetitive inhibition partly by affecting some secondary structures around them. Thus, the destructions of these secondary structures can lead to conformational changes in some important regions, such as the catalytic triad and acyl pocket. The investigation is of great importance to the discovery of good fullerene inhibitors.

Communicated by Ramaswamy H. Sarma