PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Biochem J. Author manuscript; available in PMC 2010 October 4.
Published in final edited form as:
PMCID: PMC2949283
NIHMSID: NIHMS239142

Accommodation of physostigmine and its analogs by acetylcholinesterase is dominated by hydrophobic interactions

Abstract

The role of the functional architecture of the human acetylcholinesterase (HuAChE) in reactivity toward the carbamates pyridostigmine, rivastigmine and several analogs of physostigmine, that are currently used or considered for use as drugs for Alzheimer’s disease, was analysed using over 20 mutants of residues that constitute the interaction subsites in the active center. Both steps of the HuAChE carbamylation reaction, formation of the Michaelis complex as well as the nucleophilic process, are sensitive to accommodation of the ligand by the enzyme. For certain carbamate - HuAChE combinations, the mode of inhibition shifted from a covalent to a noncovalent type, according to the balance between dissociation and covalent reaction rates. Whereas the charged moieties of pyridostigmine and rivastigmine contribute significantly to the stability of the corresponding HuAChE complexes, no such effect was observed for physostigmine and its analogs, phenserine and cymserine. Moreover, physostigmine-like ligands carrying oxygen instead of nitrogen at position – 1 of the tricyclic moiety (physovenine and tetrahydrofurobenzofuran analogs) displayed comparable structure – function characteristics toward the various HuAChE enzymes. The essential role of the HuAChE hydrophobic pocket, comprising mostly of residues Trp86 and Tyr337, in accommodating (−)-physostigmine and in conferring ~300-fold stereoselectivity toward physostigmines, was elucidated through examination of the reactivity of selected HuAChE mutations toward enantiomeric pairs of different physostigmine analogs. The present study demonstrates that certain charged and uncharged ligands, like analogs of physostigmine and physovenine, seem to be accommodated by the enzyme mostly through hydrophobic interactions.

1. INTRODUCTION

The enzyme acetylcholinesterase (AChE) is presently the most important molecular target for therapeutic intervention in symptomatic treatment of senile dementia of the Alzheimer’s type (AD) [1]. The ongoing effort to develop more therapeutically efficacious AChE inhibitors is currently driven by the remarkable progress, made during the last 15 years, in elucidating the structural and functional properties of the enzyme through x-ray crystallography [2,3] and site directed mutagenesis [4-8]. Combination of these two powerful techniques allowed for the detailed mapping of the human acetylcholinesterase (HuAChE) active center, delineating the functional subsites involved in reactivity toward substrates and other covalent modifiers as well as toward noncovalent ligands specific for the active center. These include the catalytic triad (Ser203, His447, Glu334), the “oxyanion hole” consisting of residues Gly120, Gly121 and Ala204 as well as different combinations of the 14 aromatic amino acids which line about 40% of the HuAChE active-center gorge surface, e.g. the acyl pocket (Phe295 and Phe297); the ‘hydrophobic subsite’ (Trp86, Tyr133, Tyr337 and Phe338); the cation - π interaction locus for charged moieties of substrates and other ligands at the active center (Trp86) and the peripheral anionic site (PAS, Tyr72; Tyr124 and Trp286).

Further examination of the functional architecture of the HuAChE active center revealed that reactivity of the enzyme toward substrates and other ligands can be also affected through perturbation of functional domains which may include multiple subsites in the active center. Thus, enhanced conformational mobility of the catalytic histidine was recently implicated in the activity differences between human butyrylcholinesterase (HuBChE) and the hexamutant HuAChE carrying aliphatic replacements of all the active site gorge aromatic residues (Tyr72, Tyr124, Trp286, Phe295, Phe297, Tyr337) distinguishing between the two enzymes [9,10]. Modulation of ligand interactions with the enzyme can be also effected through disruption of polar networks in the active center. One of these may include residues Ser229 and the catalytic triad residue Glu334 [11].

Most of the AChE inhibitors approved for clinical use as AD drugs (Cognex, Aricept, Nivalin) are noncovalent inhibitors and hence their AChE complexes are amenable for crystallographic analysis [12-14]. For the recently approved covalent AChE modifier, the carbamate rivastigmine (Exelon), such analysis can be carried out only on the carbamylated enzyme and is therefore relevant primarily to the enzyme regeneration step [15]. Yet, the overall inhibition process by carbamates is determined by properties of both the carbamylated enzyme and the transient Michaelis complex. These two species determine the rates of decarbamylation and carbamylation, respectively, and hence both contribute to the efficacy of the carbamate as a drug. Therefore, dissection of the affinity characteristics toward carbamates, through functional analysis of the carbamate-HuAChE Michaelis complexes, should provide information relevant to the design of more efficacious carbamate AD therapeutics [16]. In the past, we have shown that functional analysis of such Michaelis complexes can be carried out much in the same manner as for the noncovalent ligands [8].

In the present study, we examined the reactivity of HuAChE enzymes modified at relevant binding subsites, toward the carbamates rivastigmine and analogs of physostigmine, that are currently used [17,18] or considered for use as AD drugs [19-22]. Elements of the binding environment that determine the fate (carbamylation or dissociation; see scheme 1) of the particular Michaelis complexes have been identified. We conclude that interactions of HuAChE enzymes with analogs of physostigmine are dominated by hydrophobic interactions of the tricyclic eseroline moiety, and therefore the properties of the corresponding Michaelis complexes are quite different from those of rivastigmine and pyridostigmine. Thus, functional analysis appears to be a tool of choice for analysis of molecular complexes, which are too unstable for structural studies and yet are important as templates for drug design.

Scheme 1
Under conditions where k−1>>k2: kd= k−1/k1 (dissociation constant of the Michaelis complex); ki= k2/kd (bimolecular rate constant of corbamylation)

MATERIALS AND METHODS

Materials and enjymes

Acetylthiocholine iodide (ATC), 5,5′-dithiobis-(2-nitrobenzoic acid), pyridostigmine bromide and (−)physostigmine salicylate were purchased from Sigma, and rivastigmine was obtained from Teva Ltd. (+)-Physostigimine, enantiomers of physovenine, phenserine as well as the enantiomeric pairs of cymserine, cymyl carbamate of physovenol and cymyl carbamate of tetrahydrofurobenzofurol were synthesized according to published procedures [23-26]. The chemical structures of the AChE inhibitors are shown in Figure 1, and full chemical characterization was performed to insure chemical and chiral purity.

Figure 1
Chemical structures of the carbamates used in this study. Only the (−)-3aS-enantiomers of the various physostigmine analogues are shown. The numbering shown for eseroline is applicable to all the other physostigmine analogues.

Expression of recombinant enzymes, as well as the construction of all the HuAChE mutants, was described previously [4,6,7,9,27-31]. Construction of the double mutant W86A/Y337A was carried out by replacement of the appropriate DNA fragments of the AChE-w7 variant [4] with the respective fragments of the W86A and Y337A variants. Stable recombinant cell clones expressing high levels of each of the mutants were established according to the procedure described previously [27]. Enzymes were purified (over 90% purity) as described previously [27,32].

Determination of HuAChE activity and analyisis of kinetic data

HuAChE activity was assayed according to the method of Ellman et al [33] in the presence of 0.1 mg/ml BSA, 0.3 mM DTNB, 50 mM sodium phosphate buffer (pH 8.0) and various concentrations of ATC at 27 °C and monitored with a Thermomax microplate reader (Molecular Devices).

The rate constants of progression of the carbamylation reactions (see Scheme 1) were estimated for at least four different concentrations (and at least a 10-fold in ligand concentration, around the estimated value of Kd) of carbamate (CR), by adding substrate at various time intervals and measuring the enzyme residual activity (E) (enzyme concentration was ca. 1.0nM). To avoid interference from regeneration of enzyme activity due to dissociation of enzyme carbamate conjugates, the initial velocity was used to determine kobs (V=kobs [E]) at each carbamate concentration. Thus, values of kobs were calculated from the slope of the straight lines obtained from the plots of 1/ln(E) vs. time of incubation prior to addition of substrate (cf. Figure 2 middle panels A and B). Double reciprocal plot of 1/kobs vs. 1/[CR] were used to compute k2 from the intercept, ki from the slope and Kd from the ratio between the slop and the intercept according to scheme 1 and equation 1 [8] (cf. Figure 2 right panels A and B )

1kobs=1k2+(1ki)×(1[CR])=1k2+(Kdk2)×(1[CR])
Equation 1

In cases where steady state with respect to E was formed rapidly (within few minutes, cf. panels C,D in Figure 2), and immediate recovery of full enzymatic activity was observed upon dilution (cf. Figure 3), the inhibition was treated kinetically as reversible. Thus, Lineweaver –Burk plots in the absence and in the presence of different carbamate concentrations (cf. Figure 2 middle panels C and D) yielded values of relative slopes Rs (note that the Rs = (1+ (1/ Kd)[CR]) according to ref 6 and 31). The Rs values were plotted against the carbamate concentration and the reciprocal of the slope provided the Kd values (cf. Fig. 2 right panels C and D).

Fig. 2
Derivation of inhibition rate and equilibrium constants of wild type HuAChE and its mutant F295A/F338A by the carbamates physostigmine, phenserine and cymserine. Left panels: Inhibition time curve at indicated concentrations of the various carbamates ...
Fig. 3
Regeneration of enzymatic activity, following 50-fold dilution, from HuAChE phenserine and cymserine conjugates.

Molecular modeling

Building and optimization of three-dimensional models of the HuAChE adducts with the various carbamates were performed on a Silicon Graphics workstation Octane2, using the SYBYL modeling software (Tripos Inc.). The initial models were constructed by manual docking of the ligands into the HuAChE active centre guided by interactions with residue Trp86, the active-site nucleophile (Ser203) and residues of the oxyanion hole. The initial models were optimized by molecular mechanics using the MAXMIN force field (and AMBER charge parameters for the enzyme) and zone-refined, including 122 amino acids [1.5 nm (15 Å) substructure sphere around γO-Ser203]. Initial optimization included restriction of the distances between the carbonyl oxygen and the amide nitrogen atoms of residues Gly121 and Gly122, that between the carbonyl carbon and γO-Ser203 as well as positions of residues Cys69 and Cys96, the ends of the omega loop. Those constraints were relieved in the subsequent refinement [8].

RESULTS AND DISCUSSION

Reactivity of rivastigmine toward HuAChE enzymes modified at various binding subsites in the active center

In a previous study, comparison of the reactivity characteristics of rivastigmine and pyridostigmine toward HuAChE enzymes suggested that accommodation of these carbamates in the active center is analogous to that of noncovalent inhibitors like edrophonium [16]. Results presented here indicate that although pyridostigmine and rivastigmine share the same binding subsites in the HuAChE active center, their distinct orientations with respect to the active site seem to influence the outcome of the carbamylation process. These results are summarized in Table 1 that includes HuAChE enzymes that carry replacements at the hydrophobic pocket, H-bond network, oxyanion hole, acyl pocket and the peripheral anionic site [6-8,10,11,30,31]. In addition, we report on reactivities of both carbamates toward HuAChEs, which were engineered to resemble HuBChE active center [9,11].

Table 1
Effects of mutations at the various subsites of the HuAChE active center on reactivity toward pyridostigmine and rivastigminea

Replacements of aromatic residues comprising the HuAChE active center hydrophobic pocket had a similar effect on the rates of carbamylation by rivastigmine and by pyridostigmine, implying that in both cases the positively charged moiety interacts with the cation binding subsite - Trp86. The pronounced increase of the respective dissociation constants, due to replacement of Trp86 (4400- and 6150-fold for pyridostigmine and rivastigmine respectively), resembles that for all the charged active center inhibitors [8,34,35]. Like for pyridostigmine [8], replacement of Tyr133 by alanine but not by phenylalanine had a pronounced effect on the affinity of Y133A HuAChE toward rivastigmine. As noted previously, replacement of Tyr337 by alanine had little effect on interactions of cationic ligands, although corresponding crystal structures of AChE complexes [12-14] and molecular model of HuAChE-pyridostigmine Michaelis complex [8] indicate close proximity of this residue to the ligand charged moiety.

Perturbations of the H-bond network, through replacements of residues Tyr133, Glu202 and Glu450 [10], had a relatively uniform effect on the corresponding rates of carbamylation by pyridostigmine (see Table 1). Yet, for two of those enzymes, E202Q and E450A, interaction with rivastigmine did not lead to carbamylation but rather to a regular, albeit low affinity (for corresponding values of the dissociation constants Kd see Table 1), noncovalent inhibition. This observation suggests that the balance between the rates of carbamylation and of dissociation of the corresponding Michaelis complexes can be easily tipped away from the covalent reaction. This facet of carbamate reactivity toward HuAChEs will become even more evident for certain analogs of physostigmine.

Structural modification of the oxyanion hole through replacement of residue Gly121 by alanine [31] alters the reactivity characteristics of both carbamates, converting them into noncovalent inhibitors. We have already shown that interactions of the acyl oxygen (acetyl, carbamyl or phosphoryl) with the oxyanion hole are important for both stabilization of the Michelis complex and activation of the acyl moiety for nucleophilic attack by the catalytic Ser203 [31]. Thus, dissociation of the G121A HuAChE- rivastigmine complex is probably much faster than that of the corresponding complex of the wild type enzyme (the pronounced increase of the value of Kd is mostly due to increase of the dissociation rate constant k−1), while its conversion to carbamylated enzyme is slower.

Replacements at the peripheral anionic site had only a limited effect on the carbamylation rate constants by rivastigmine. In particular, the corresponding value of ki for carbamylation of D74N HuAChe was 3-fold lower than that of the wild type enzyme. On the other hand, carbamylation of this enzyme by pyridostigmine was nearly 150-fold slower, with the corresponding value of Kd being 200-fold higher than that of the wild type HuAChE. It is already reported that this replacement resulted in a 50-fold increase of the dissociation constant for tacrine, while having only a small effect on the corresponding constant for edrophonium (5-fold increase) and no effect on huperzine A [8]. The reason for these uneven effects on affinities of the D74N enzyme toward the various charged (at the experimental pH) active center ligands still remains elusive.

While replacement of acyl pocket residue Phe295 by alanine, had little effect on the carbamylation rate by pyridostigmine, the corresponding rate for rivastigmine was 30-fold higher. This observation seems consistent with the size of the substituents on the carbamyl nitrogen. Namely, while interaction of rivastigmine with residue Phe295 of the WT HuAChE may perturb the “aromatic trap” and affect the carbamylation step, such perturbation is avoided in binding to the F295A enzyme. Accordingly, analogous substitution of the second acyl pocket residue Phe297 had a similar effect on the carbamylation rates by both carbamates. Reactivity of rivastigmine toward HuBChE, in which the corresponding acyl pocket is lined by aliphatic residues, has been found to be 15-fold higher than toward HuAChE. Most of this difference was due to the 9-fold higher value of the carbamylated enzyme formation step rate constant k2 (see Table 1). Such reactivity enhancement was not observed for HuAChE enzymes in which the active center was engineered to resemble that of HuBChE. For the “butyryl – like” enzyme carrying replacements of aromatic residues vicinal to the active site (F295L/F297V/Y337A HuAChE) as well as for that substituted at both the active center and the peripheral anionic site (Y72N/Y124Q/W286A /F295L/F297V/Y337A HuAChE) [9], carbamylation rates by rivastigmine resembled those of the wild type enzyme. On the other hand, the corresponding rates of carbamylation by pyridostigmine were 15-fold and 73-fold lower, respectively.

Analogues of physostigmine display distinct inhibition profiles toward HuAChE enzymes

Unlike the pronounced effects (over 4000-fold) of replacing Trp86 by alanine on the inhibitory activities of pyridostigmine and rivastigmine, the inhibition rate constant of physostigmine toward the W86A HuAChE was only less than 50-fold lower than toward the wild type enzyme (see Table 2). Replacement of Tyr133 by alanine had larger effect (700-fold), implying that steric congestion rather than cation-π interaction is the dominant factor in the accommodation of physostigmine in the hydrophobic pocket [28]. Perturbations of the H-bond network affect reactivities of the corresponding enzymes toward physostigmine to a similar extent as for pyridostigmine and rivastigmine. Modifications of the acyl pocket had lower effect on reactivities toward physostigmine than toward rivastigmine (see Tables Tables11,,22 and reference 8)

Table 2
Effects of mutations at the various subsites of the HuAChE active center on reactivity toward physostigmine phenserine and cymserinea

To further explore the reactivity characteristics of HuAChE toward physostigmine, we now examine physostigmine analogues differing in substitution at the carbamyl nitrogen as well as analogs with a modified tricyclic moiety. Thus, physostigmine, phenserine and cymserine (see Figure 1) were expected to display similar accommodation in the HuAChE hydrophobic pocket while differing with respect to the acyl pocket and the peripheral anionic site. From the results exemplified in Figure 2 and the regeneration of HuAChE activity from the respective enzyme-inhibitor conjugate (Figure 3), it appears that while the inhibition characteristics of phenserine toward the HuAChE enzymes resembles that of physostigmine, the reactivity of cymserine is that of a noncovalent inhibitor

Notwithstanding the difference in the inhibitory activity of physostigmine and phenserine as compared to cymserine, it seems that the three compounds are similarly accommodated in the hydrophobic pocket. In all cases affinities were affected by replacements at positions 86 and 133 by alanine, while not being sensitive to substitutions of Tyr337. Substitution of position 338 by alanine had some effect on the values of Kd for phenserine and cymserine (30-fold and 6-fold, respectively) but not on the corresponding value for physostigmine (see Table 2). Thus it appears that the different reactivity characteristics of cymserine, as compared to physostigmine and phenserine is not due to interactions of the respective eseroline moieties with the HuAChE hydrophobic pocket. Therefore, it seems reasonable to assume that in the HuAChE- cymserine Michaelis complexes the ligand is sub-optimally oriented with respect to the catalytic machinery of the enzyme.

Substitutions at the H-bond network had similar effects on affinities toward phenserine and cymserine, as was the case for most of the replacements at the peripheral anionic site. The outstanding case was the failure of phenserine to carbamylate W286A HuAChE. Similar noncovalent inhibition was observed for interactions of phenserine with F295A/F338A and with the "butyryl like" HuAChEs (see Table 2).

Perturbations of the acyl pocket had comparable effects on the affinities toward phenserine and cymserine, as well as minor effects on the carbamylation rate constant (k2) for phenserine, Thus, it does not appear that the failure of cymserine to carbamylate HuAChE enzymes results from delocalization of the catalytic His447 due to perturbations of the acyl pocket [10]. Moreover, while the prototypical perturbation of the His447 positioning (F295A/F338A HuAChE see [29]), abolished carbamylation by phenserine, the reactivity was restored by additional replacement at position 407. Such a reactivity profile [29] could not be observed in the case of cymserine (Table 2).

In order to gain further insight into the unique inhibitory characteristics of cymserine toward all the HuAChE enzymes, molecular modelling experiments of the Michaelis complexes of wild type HuAChE with physostigmine, phenserine and cymserine have been constructed (Figure 4). While in the model of phenserine, interactions of the N-aryl moiety with residues of the acyl pocket and with Phe338 could be observed, the disposition of the carbamyl moiety with respect to the active site residues resembled that of physostigmine. On the other hand, interactions of the 4-isopropyl aromatic substituent of cymserine with aromatic residues lining the mouth of the active center gorge forced an alternative conformation of the N-aryl moiety and consequently the removal of the carbamyl oxygen from the oxyanion hole. This finding seems consistent with the reversible noncovalent nature of cymserine interaction with wild type HuAChE and any of its mutant derivatives.

Figure 4
Relative orientations of (−)-physostigmine and (−)-cymserine in models of their respective Michaelis complexes with HuAChE. Superposition of the carbamyl moieties is emphasized by the shaded area. Please notice that the carbamyl oxygen ...

The notion that cymserine fails to carbamylate HuAChEs (see Figures Figures22 and and3)3) due to impaired polarization of the carbamyl oxygen, is consistent with the carbamylation profile of the G121A enzyme with physostigmine [31] as well as with pyridostigmine and rivastigmine.

Accommodation of the tricyclic moiety in the hydrophobic pocket

Another manifestation of the tight accommodation of the tricyclic moiety in the hydrophobic pocket is the marked stereoselectivity toward the (−)-physostigmine (3aS-diastereomer) enantiomeric form [20]. Since the main difference between the diastereomers is in the eseroline moiety, stereoselectivity should originate from the asymmetric interactions in the hydrophobic pocket. It was therefore reasonable to assume that through identifying the specific interactions leading to stereoselectivity toward (−)-physostigmine, a better understanding of its accommodation in the hydrophobic pocket would be achieved.

Replacement of the tricyclic eseroline group in physostigmine by physovenyl moiety (see Figure 1) resulted in an analogous physovenine [24], with similar chirality due to asymmetric carbon at position - 3a. Yet, the overall inhibitory activities, and in particular the affinities of both diastereomers of physovenine toward HuAChE were similar to that of (−)-physostigmine (see Table 3). It has been already proposed that the low affinity of AChE toward (+)-physostigmine is due to its N1-methyl group interfering with Trp84 (86 in HuAChE) [36]. However this residue is usually thought to interact with N-methyl groups, as in the case for its endogenous substrate – acetylcholine [4,6].

Table 3
Stereoselectivity of HuAChE enzymes mutated at the hydrophobic pocket and at the acyl pocket toward (+)- and (−)-enantiomers of physostigmine and physoveninea

As already reported [8,11], replacement of Trp86 by alanine had a moderate effect (~25-fold) on the affinity toward the natural [(−)-physostigmine]. Here we find that this modification has practically no effect on the affinity toward the (+)-diastereomer (see Table 3). Thus, the diminished stereoselectivity exhibited by W86A HuAChE toward physostigmine diastereomers, as compared to the wild type enzyme (20-fold, see Table 3), was solely due to a decrease in affinity toward (−)-physostigmine. This suggests that residue Trp86 does not participate at all in the interactions of HuAChE hydrophobic pocket with (+)-physostigmine (see Figure 5). On the other hand, W86A HuAChE displayed a similar decrease in affinity toward both diastereomers of physovenine and therefore, both the wild type and the W86A HuAChEs show nearly no stereoselectivity toward the physovenines.

Figure 5
Gradual displacement of the cation binding residue Trp86 from the hydrophobic pocket. Models of HuAChE Michaelis complexes with pyridostigmine and physostigmine enantiomers are shown in panels A-C. A, The positioning of residue Trp86 in the Michaelis ...

Substitution of Tyr337 by alanine maintained the stereoselectivity toward physostigmines while inducing limited stereoselectivity toward the physovenines (18-fold, see Table 3). For both cases, stereoselectivity is completely abolished in the double mutant W86A/Y337A HuAChE. We note that while replacement of the aromatic residues Trp86 and Tyr337 led to a moderate decrease in affinity toward (−)-physostigmine (28- and 3-fold respectively), the corresponding effect for the double mutant is quite dramatic (4100-fold). In contrast, only a 10-fold decline in affinity of the W86/Y337 HuAChE toward the (+)-enantiomer has been observed, demonstrating that interactions with the hydrophobic pocket determine stereoselectivity toward physostigmines (see Table 3). For both physovenine enantiomers, affinities of the double mutant are comparable to those of the W86A enzyme. Thus, physostigmines and physovenines seem to be somewhat differently oriented with respect to the hydrophobic patch in the active center [8]. While residue Trp86 is essential in accommodation of (−)-physovenine, residues Trp86 and Tyr337 seem to compensate for one another in the case of (−)-physostigmine. Such compensation seems to account for the intriguing observation that removal of the aromatic moiety from position 86 had a larger effect on affinity toward the uncharged physovenines than toward the charged (−)-physostigmine.

The results described above seem consistent with the idea that steric congestion of Trp86 and the N1-methyl of (+)-physostigmine, indeed, interferes with accommodation of this diastereomer in the hydrophobic pocket. To further examine this hypothesis, affinities of HuAChE enzymes modified at the hydrophobic pocket, toward cymserine and cymyl carbamates of physovenol and of tetrahydrofurobenzofuran [26], have been compared. As for cymserine, all the cymyl analogs are noncovalent inhibitors. While HuAChE displayed ~50-fold stereoselectivity toward (−)-cymserine, practically no stereoselectivity was observed toward the diastereomers of analogs bearing the physovenol and tetrahydrofurobenzofuran moieties. Replacement of Trp86 by alanine practically abolished stereoselectivity toward diastereomers of cymserine yet had only a limited effect on their binding affinities. Other cymyl carbamates were also little affected by residue replacements at the hydrophobic pocket (see Table 4).

Table 4
Stereoselectivity of HuAChE enzymes mutated at the hydrophobic pocket toward (+)- and (−)-enantiomers of cymserine and its analoguesa

The notion that AChE stereoselectivity toward physostigmine is mainly due to interactions of the alkyl substituent at position -1 is supported also by previous studies on physostigmine analogues [36-38]. In particular, it was interesting to examine the low AChE stereoselectivity toward analogs of 8-carbaphysostigmine, since these structures do contain the N(1)- alkyl substituent [39]. Examination of molecular models of the corresponding Michaels complexes indicates that due to bending of the tricyclic moiety at the sp3-C8, both enantiomers could be accommodated in the hydrophobic pocket without steric occlusion of Trp86. Thus the structural features of the eseroline moiety that contribute to AChE stereoselectivity toward physostigmine, are the alkyl substituent at position -1 combined with the planar disposition of the tricyclic ring system.

Accommodation of carbamates in the HuAChE active center

Carbamates are unique HuAChE inhibitors, binding both as covalent and noncovalent ligands to the different HuAChE enzymes. It appears that this property of carbamates originates from the particular dependence of the carbamyl moiety reactivity on its juxtaposition with the elements of the enzyme catalytic machinery. Carbamylation of AChEs involves nucleophilic attack on a relatively nonreactive carbonyl group, and therefore its rate depends critically upon the stability of the corresponding Michaelis complex, which manifests itself predominantly by variation in the values of the dissociation rate constant k−1 (under equilibrium conditions Kd = k−1/k1). Thus, the efficiency of the carbamylation process depends mainly on the relative values of the rate constants k−1 and k2 (see scheme 1), with the latter displaying rather limited variability [40]. Accommodation of the carbamylating agent in the AChE active center is hence the most significant molecular event in the carbamylation process. The finding that the affinity of HuAChE toward the charged physostigmine is remarkably similar to that toward the structurally similar yet uncharged physovenine, indicates that both inhibitors are accommodated mainly through hydrophobic interactions.

ACKNOWLEDGEMENT

This works was supported by the US Army Medical Research and Material Command under contract DAMD17-00-C-0021 (to A.S.), grant from Life Science Research Israel Ltd and Intramural Research Program of the National Institute on Aging, NIH.

Abbreviations

ATC
acetylthiocholine
AChE
acetylcholinesterase
AD
Alzheimer’s disease
HuAChE
human acetylcholinesterase
HuBChE
human butyrylcholinesterase

REFERENCES

1. Giacobini E. Cholinesterase inhibitors: new roles and therapeutic alternatives. Pharmacol. Res. 2004;50:433–440. [PubMed]
2. Sussman JL, Harel M, Frolow F, Oefner C, Goldman A, Toker L, Silman I. Atomic structure of acetylcholinesterase from Torpedo californica: a prototypic acetylcholine-binding protein. Science. 1991;253:872–9. [PubMed]
3. Kryger G, Harel M, Giles K, Toker L, Velan B, Lazar A, Kronman C, Barak D, Ariel N, Shafferman A, Silman I, Sussman JL. Structures of recombinant native and E202Q mutant human acetylcholinesterase complexed with the snake-venom toxin fasciculin-II. Acta Crystallogr. 2000;56:1385–1394. [PubMed]
4. Shafferman A, Velan B, Ordentlich A, Kronman C, Grosfeld H, Leitner M, Flashner Y, Cohen S, Barak D, Ariel N. Substrate inhibition of acetylcholinesterase residues affecting signal transduction from the surface to the catalytic center. EMBO J. 1992;11:3561–3568. [PubMed]
5. Taylor P, Radic Z. The cholinesterases. Annu. Rev. Pharmacol. Toxicol. 1994;34:281–320. [PubMed]
6. Ordentlich A, Barak D, Kronman C, Flashner Y, Leitner M, Segall Y, Ariel N, Cohen S, Velan B, Shafferman A. Dissection of the human residues constituting the anionic site, the hydrophobic site, and the acyl pocket. J. Biol. Chem. 1993;268:17083–17095. [PubMed]
7. Barak D, Kronman C, Ordentlich A, Ariel N, Bromberg A, Marcus D, Lazar A, Velan B, Shafferman A. Acetylcholinesterase peripheral anionic site degeneracy conferred by amino acid arrays sharing a common core. J. Biol. Chem. 1994;269:6296–630. [PubMed]
8. Ariel N, Ordentlich A, Barak D, Bino T, Velan B, Shafferman A. The ‘aromatic patch’ of three proximal residues in the human acetylcholinesterase active centre allows for versatile interaction modes with inhibitors. Biochem. J. 1998;335:95–102. [PubMed]
9. Kaplan D, Orentlich A, Barak D, Ariel N, Kronman C, Velan B, Shafferman A. Does “butyrylization” of acetylcholinesterase through substitution of the six divergent aromatic amino acids in the active center gorge generate an enzyme mimic of butyrylcholinesterase. Biochemistry. 2002;41:8245–8252. [PubMed]
10. Ordentlich A, Barak D, Kronman C, Benschop HP, De Jong LP, Ariel N, Barak R, Segall Y, Velan B, Shafferman A. Exploring the active center of human Acetylcholinesterase with stereoisomers of an organophosphorus inhibitor with two chiral centers. Biochemistry. 1999;38:3055–3066. [PubMed]
11. Shafferman A, Barak D, Kaplan D, Ordentlich A, Kronman C, Velan B. Functional requirements for the optimal catalytic configuration of the AChE active center. Chem. Biol. Interact. 2005;157-158:123–31. [PubMed]
12. Harel M, Schalk I, Ehert-Sabatier L, Bouet F, Goeldner M, Hirth C, Axelson PH, Silman I, Sussman JL. Quaternary ligand binding to aromatic residues in the active-site gorge of acetylcholinesterase. Proc.Natl. Acad. Sci. USA. 1993;90:9031–35. [PubMed]
13. Kryger G, Silman I, Sussman JL. Structure of acetylcholinesterase complexed with E2020 (Aricept): implications for the design of new anti-Alzheimer drugs. Structure. 1999;7:297–307. [PubMed]
14. Greenblatt HM, Dvir H, Silman I, Sussman JL. Acetyl-cholinesterase. J. Molec. Neurosci. 2003;20:369–383. [PubMed]
15. Bar-On P, Millard CB, Harel M, Dvir H, Enz A, Sussman JL, Silman I. Kinetic and structural studies on the interaction of cholinesterases with the anti-Alzheimer drug rivastigmine. Biochemistry. 2002;41:3555–64. [PubMed]
16. Barak D, Ordentlich A, Kaplan D, Kronman C, Velan B, Shafferman A. Lessons from functional analysis of AChE covalent and noncovalent inhibitors for design of AD therapeutic agents. Chem. Biol. Interact. 2005;157-158:219–26. [PubMed]
17. Jann MW. Rivastigmine, a new generation cholinesterase inhibitor for the treatment of Alzheimer’s disease. Pharmacotherapy. 2000;20:1–12. [PubMed]
18. Williams BR, Nazarians A, Gill MA. A review of rivastigmine: a reversible cholinesterase inhibitor. Clin. Ther. 2003;25:1634–1652. [PubMed]
19. Groner E, Ashani Y, Schorer-Apelbaum D, Sterling J, Herzig Y, Weinstock M. The kinetics of inhibition of human acetylcholinesterase and butyrylcholinesterase by two series of novel carbamates. Mol. Pharmacol. 2007;71:1610–1617. [PubMed]
20. Atack JR, Yu Q-S, Soncrant TT, Brossi A, Rapoport RI. Comparative inhibitory effects of various physostigmine analogs against acetyl- and butyrylcholinesterase. J. Pharmacol. Exp. Ther. 1989;249:194–202. [PubMed]
21. Greig NH, Pei XF, Soncrant TT, Ingram DK, Brossi A. Phenserine and ring C hetero-analogues drug candidates for trearment of Alzheimer’s disease. Med. Res. Rev. 1995;15:3–31. [PubMed]
22. Bartolucci C, Siotto M, Ghidini E, Amari G, Bolzoni PT, Racchi M, Villetti G, Delcanale M, Lamba D. Structural determinants of Torpedo californica acetylcholinesterase inhibition by the novel and orally active carbamate based anti-Alzheimer drug ganstigmine (CHF-2819) J. Med. Chem. 2006;49:5051–5058. [PubMed]
23. Dale FJ, Robinson B. The synthesis and anti-cholinesterase activities of (+)-physostigmine and (+)-physovenine. J. Pharm. Pharmacol. 1979;22:889–896. [PubMed]
24. Yu Q, Liu C, Brzostowska M, Crisley L, Brossi A, Greig NH, Attack JR, Soncrant TT, Rapoport SI, Radunz HE. Physovenines: Efficient synthesis of (−)- and (+)-physovenine and synthesis of carbamate analogues of (−)- physovenine. Anticholinesterase activity and analgesic properties of optically active physovenines. Helv. Chim. Acta. 1991;74:761–766.
25. Yu Q-S, Holloway HW, Utsuki T, Brossi A, Greig NH. Synthesis of novel phenserine-based-selective inhibitors of butyrylcholinesterase for Alzheimer’s disease. J. Med. Chem. 1999;42:1855–1861. [PubMed]
26. Luo W, Yu Q-S, Zhan M, Parrish D, Deschamps JR, Kulkarani SS, Holloway HW, Allay GM, Lahiri DK, Brossi A, Greig NH. Novel anticholinesterases based on the molecular skeletons of furobenzofuran and mehtanobenzodioxepine. J. Med. Chem. 2005;48:986–994. [PubMed]
27. Kronman C, Velan B, Gozes Y, Leitner M, Flashner Y, Lazar A, Marcus D, Sery T, Papier A, Grosfeld H, Cohen S, Shafferman A. Production and secretion of high levels of recombinant human acetylcholinesterase in cultured cell lines: microheterogeneity of the catalytic subunit. Gene. 1992;121:295–304. [PubMed]
28. Ordentlich A, Barak D, Kronman C, Ariel N, Segall Y, Velan B, Shafferman A. Contribution of aromatic moieties of Tyr-133 and of the “anionic subsite” Trp-86 to catalytic efficiency and allosteric modulation of acetylcholinesterase. J. Biol. Chem. 1995;270:2082–2091. [PubMed]
29. Barak D, Kaplan D, Ordentlich A, Ariel N, Velan B, Shafferman A. The “caging” of the catalytic histidine by aromatic residues is essential for efficient catalysis of acetylcholinesterase. Biochemistry. 2002;41:8245–8252. [PubMed]
30. Ordentlich A, Kronman C, Barak D, Stein D, Ariel N, Marcus D, Velan B, Shafferman A. Engineering resistance to ‘aging’ in phosphylated human acetylcholinesterase - role of hydrogen bond network in the active center. FEBS Lett. 1993;334:215–220. [PubMed]
31. Ordentlich A, Barak D, Kronman C, Ariel N, Segall Y, Velan B, Shafferman A. Functional characteristics of the oxyanion hole in human acetylcholinesterase. J. Biol. Chem. 1998;273:19509–19517. [PubMed]
32. Barak D, Ordentlich A, Bromberg A, Kronman C, Marcus D, Lazar A, Ariel N, Velan B, Shafferman A. Allosteric modulation of acetylcholinesterase activity by peripheral ligands involves a conformational transition of the anionic subsite. Biochemistry. 1995;34:15444–15452. [PubMed]
33. Ellman G,L, Courtney KD, Andres V, Featherstone RM. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem. Pharmacol. 1961;7:88–95. [PubMed]
34. Ordentlich A, Barak D, Sod-Moria G, Kaplan D, Mizrahi D, Segal Y, Kronman C, Karton Y, Lazar A, Marcus D, Velan B, Shafferman A. Stereoselectivity toward VX is determined by interactions with residues of the acyl pocket as well as of the peripheral anionic site of AChE. Biochemistry. 2004;43:11255–11265. [PubMed]
35. Silman I, Sussman JL. Acetylcholinesterase: ‘classical’ and ‘non-classical’ functions and pharmacology. Curr. Opin. Pharmacol. 2005;5:293–302. [PubMed]
36. Yu Q-S, Atack JR, Rapoport SI, Brossi A. Synthesis and acetylcholinesterase activity of (−)-N1-Norphysostigmine (−)-Eseramine and other (N1)-substituted analogues of (−)-Physostigmine. J. Med. Chem. 1988;31:2297–2300. [PubMed]
37. Yu Q-S, Pei X-F, Holloway HW, Greig NH. Total synthesis and anticholinesterase activities of (3aS)-N(8)-Norphysostigmine (3aS)-N(8)-Norphenserine, their antipodal isomers and other N(8)-substituted analogues. J. Med. Chem. 1997;40:2895–2901. [PubMed]
38. Yu Q-S, Luo W, Holloway HW, Utsuki T, Perry T, Lahiri DK, Greig NH, Brossi A. Racemic N1-Phenserine and its enantiomers: unpredicted inhibition of human acetyl- and butyrylcholinesterase and β-amyloid precursor protein in vitro. Heterocycles. 2003;61:529–539.
39. Chen YL, Nielsen J, Hedberg K, Dunaiskis A, Jones S, Russo L, Johnson J, Ives J, Liston D. Syntheses, resolution, and structure-activity relationships of potent acetylcholinesterase inhibitors: 8-carba-physostigmine analogs. J. Med.Chem. 1992;35:1429–1434. [PubMed]
40. Hetnarski B, O’Brien R. Electron-donor and affinity constants and their application to the inhibition of acetylcholinesterase by carbamates. J. Agric. Food Chem. 1975;23:709–713. [PubMed]