Nanosecond Dynamics of Acetylcholinesterase Near the Active Center Gorge

Connectivity between surface and interior in catalytic subunits of acetylcholinesterases inferred from their X-ray structures

Catalytic activity and function of acetylcholinesterase (AChE; EC 3.1.1.7) have been recognized and studied for over a century and its quaternary and primary structures for about half a century, and its tertiary structure has been known for about 33 years. Clear understanding of relationships between the structure and the function is still pending for this enzyme. Hundreds of crystallographic, static snapshots of AChEs from different sources reveal largely one general backbone conformation with narrow entry into the active center gorge, tightly fit to accept one acetylcholine (ACh) molecule, in contrast to its high catalytic turnover. This short review of available X-ray structures of AChEs from electric ray Torpedo californica, mouse and human, finds some limited, yet consistent deviations in conformations of selected secondary structure elements of AChE relevant for its function. The observed conformational diversity of the acyl pocket loop of AChE, unlike the large Ω-loop, appears consistent with structurally dynamic INS data and solution-based SAXS experiments to explain its dominant role in controlling the size of the active center gorge opening, as well as connectivity between the immediate surroundings of the buried active Ser, and catalytically relevant sites on the AChE surface.

Computational Studies on Acetylcholinesterase Inhibitors: From Biochemistry to Chemistry

Acetylcholinesterase inhibitors are the most promising therapeutics for Alzheimer's disease treatment as these prevent the loss of acetylcholine and slows the progression of the disease. The drugs approved for the management of Alzheimer's disease by the FDA are acetylcholinesterase inhibitors but are associated with side effects. Consistent and stringent efforts by the researchers with the help of computational methods opened new ways of developing novel molecules with good acetylcholinesterase inhibitory activity. In this manuscript, we reviewed the studies that identified the essential structural features of acetylcholinesterase inhibitors at the molecular level as well as the techniques like molecular docking, molecular dynamics, quantitative structure-activity relationship, virtual screening, and pharmacophore modelling that were used in designing these inhibitors.

Elucidating the potential neurotoxicity of chiral phenthoate: Molecular insight from experimental and computational studies

Chiral organophosphorus pollutants are existed ubiquitously in the ecological environment, but the enantioselective toxicities of these nerve agents to humans and their molecular bases have not been fully elucidated. Using experimental and computational approaches, this story was to explore the neurotoxic response process of the target acetylcholinesterase (AChE) to chiral phenthoate and further decipher the microscopic mechanism of such toxicological effect at the enantiomeric level. The results showed that the toxic reaction of AChE with chiral phenthoate exhibited significant enantioselectivity, and (R)-phenthoate (K＝1.486 × 10⁵ M^-1) has a bioaffinity for the nerve enzyme nearly three times that of (S)-phenthoate (K＝4.503 × 10⁴ M^-1). Dynamic research outcomes interpreted the wet experiments, and the inherent conformational flexibility of the target enzyme has a great influence on the enantioselective neurotoxicological action processes, especially reflected in the conformational changes of the three key loop regions (i.e. residues His-447, Gly-448, and Tyr-449; residues Gly-122, Phe-123, and Tyr-124; and residues Thr-75, Leu-76, and Tyr-77) around the reaction patch. This was supported by the quantitative results of conformational studies derived from circular dichroism spectroscopy (α-helix: 34.7%→30.2%/31.6%; β-sheet: 23.6%→19.5%/20.7%; turn: 19.2%→22.4%/21.9%; and random coil: 22.5%→27.9%/25.8%). Meanwhile, via analyzing the modes of toxic action and free energies, we can find that (R)-phenthoate has a strong inhibitory effect on the enzymatic activity of AChE, as compared with (S)-phenthoate, and electrostatic energy (-23.79/－17.77 kJ mol^-1) played a critical role in toxicological reactions. These points were the underlying causes of chiral phenthoate displaying different degrees of enantioselective neurotoxicity.

Early Detection of Growth Hormone Secretagogue Receptor Antagonists Exploiting Their Atypical Behavior in Competitive Assays

We report the proof-of-concept of a bioaffinity format designed for the early detection of growth hormone secretagogue receptor (GHS-R1a) antagonists in urine samples. We exploit here their atypical behavior in competitive experiments with labeled ghrelin (GHR), namely, the strong promoting effect on the GHR/GHS-R1a interaction at low molar ratios GHR/antagonist. The antagonists potentiate the GHR/GHS-R1a interaction, and they display the same effect on the interaction of GHS-R1a with other agonists listed as doping agents. The developed assay allows the estimation of affinity constants of ligand/receptor and antagonist/receptor binding and is amenable to optical, electrochemical, and mass-sensitive detection. The estimated affinity constants for GHR/GHS-R1a and antagonist/GHS-R1a in the absence of G proteins are in good agreement with recently reported data.

Significance Testing and Multivariate Analysis of Datasets from Surface Plasmon Resonance and Surface Acoustic Wave Biosensors: Prediction and Assay Validation for Surface Binding of Large Analytes

In this study, we performed uni- and multivariate data analysis on the extended binding curves of several affinity pairs: immobilized acetylcholinesterase (AChE)/bioconjugates of aflatoxin B₁(AFB₁) and immobilized anti-AFB₁ monoclonal antibody/AFB₁-protein carriers. The binding curves were recorded on three mass sensitive cells operating in batch configurations: one commercial surface plasmon resonance (SPR) sensor and two custom-made Love wave surface-acoustic wave (LW-SAW) sensors. We obtained 3D plots depicting the time-evolution of the sensor response as a function of analyte concentration using real-time SPR binding sensograms. These "calibration" surfaces exploited the transient periods of the extended kinetic curves, prior to equilibrium, creating a "fingerprint" for each analyte, in considerably shortened time frames compared to the conventional 2D calibration plots. The custom-made SAW sensors operating in different experimental conditions allowed the detection of AFB₁-protein carrier in the nanomolar range. Subsequent statistical significance tests were performed on unpaired data sets to validate the custom-made LW-SAW sensors.

Limitations in current acetylcholinesterase structure-based design of oxime antidotes for organophosphate poisoning

Acetylcholinesterase (AChE; EC 3.1.1.7), an essential enzyme of cholinergic neurotransmission in vertebrates, is a primary target in acute nerve agent and organophosphate (OP) pesticide intoxication. Catalytically inactive OP-AChE conjugates formed between the active-center serine and phosphorus of OPs can, in principle, be reactivated by nucleophilic oxime antidotes. Antidote efficacy is limited by the structural diversity of OP-AChE conjugates resulting from differences in the structure of the conjugated OP, the different active-center volumes they occupy when conjugated to the active-center serine of AChE, and the distinct chemical characteristics of both OPs and oximes documented in numerous X-ray structures of OP-conjugated AChEs. Efforts to improve oxime reactivation efficacy by AChE structure-based enhancement of oxime structure have yielded only limited success. We outline here the potential limitations of available AChE X-ray structures that preclude an accurate prediction of oxime structures, which are necessary for association in the OP-AChE gorge and nucleophilic attack of the OP-conjugated phosphorus.

A wrench in the works of human acetylcholinesterase: soman induced conformational changes revealed by molecular dynamics simulations

Irreversible inactivation of human acetylcholinesterase (hAChE) by organophosphorous pesticides (OPs) and chemical weapon agents (CWA) has severe morbidity and mortality consequences. We present data from quantum mechanics/molecular mechanics (QM/MM) and 80 classical molecular dynamics (MD) simulations of the apo and soman-adducted forms of hAChE to investigate the effects on the dynamics and protein structure when the catalytic Serine 203 is phosphonylated. We find that the soman phosphonylation of the active site Ser203 follows a water assisted addition-elimination mechanism with the elimination of the fluoride ion being the highest energy barrier at 6.5 kcal/mole. We observe soman-dependent changes in backbone and sidechain motions compared to the apo form of the protein. These alterations restrict the soman-adducted hAChE to a structural state that is primed for the soman adduct to be cleaved and removed from the active site. The altered motions and resulting structures provide alternative pathways into and out of the hAChE active site. In the soman-adducted protein both side and back door pathways are viable for soman adduct access. Correlation analysis of the apo and soman adducted MD trajectories shows that the correlation of gorge entrance and back door motion is disrupted when hAChE is adducted. This supports the hypothesis that substrate and product can use two different pathways as entry and exit sites in the apo form of the protein. These alternative pathways have important implications for the rational design of medical countermeasures.

Kinetic mechanism of indole-3-glycerol phosphate synthase

The (βα)(8)-barrel enzyme indole-3-glycerol phosphate synthase (IGPS) catalyzes the multistep transformation of 1-(o-carboxyphenylamino)-1-deoxyribulose 5-phosphate (CdRP) into indole-3-glycerol phosphate (IGP) in tryptophan biosynthesis. Mutagenesis data and crystal structure analysis of IGPS from Sulfolobus solfataricus (sIGPS) allowed for the formulation of a plausible chemical mechanism of the reaction, and molecular dynamics simulations suggested that flexibility of active site loops might be important for catalysis. Here we developed a method that uses extrinsic fluorophores attached to active site loops to connect the kinetic mechanism of sIGPS to structure and conformational motions. Specifically, we elucidated the kinetic mechanism of sIGPS and correlated individual steps in the mechanism to conformational motions of flexible loops. Pre-steady-state kinetic measurements of CdRP to IGP conversion monitoring changes in intrinsic tryptophan and IGP fluorescence provided a minimal three-step kinetic model in which fast substrate binding and chemical transformation are followed by slow product release. The role of sIGPS loop conformational motion during substrate binding and catalysis was examined via variants that were covalently labeled with fluorescent dyes at the N-terminal extension of the enzyme and mobile active site loop β1α1. Analysis of kinetic data monitoring dye fluorescence revealed a conformational change that follows substrate binding, suggesting an induced-fit-type binding mechanism for the substrate CdRP. Global fitting of all kinetic results obtained with wild-type sIGPS and the labeled variants was best accommodated by a four-step kinetic model. In this model, both the binding of CdRP and its on-enzyme conversion to IGP are accompanied by conformational transitions. The liberation of the product from the active site is the rate-limiting step of the overall reaction. Our results confirm the importance of flexible active loops for substrate binding and catalysis by sIGPS.

Limitations of time-resolved fluorescence suggested by molecular simulations: assessing the dynamics of T cell receptor binding loops

Time-resolved fluorescence anisotropy (TRFA) has a rich history in evaluating protein dynamics. Yet as often employed, TRFA assumes that the motional properties of a covalently tethered fluorescent probe accurately portray the motional properties of the protein backbone at the probe attachment site. In an extensive survey using TRFA to study the dynamics of the binding loops of a αβ T cell receptor, we observed multiple discrepancies between the TRFA data and previously published results that led us to question this assumption. We thus simulated several of the experimentally probed systems using a protocol that permitted accurate determination of probe and protein time correlation functions. We found excellent agreement in the decays of the experimental and simulated correlation functions. However, the motional properties of the probe were poorly correlated with those of the backbone of both the labeled and unlabeled protein. Our results warrant caution in the interpretation of TRFA data and suggest further studies to ascertain the extent to which probe dynamics reflect those of the protein backbone. Meanwhile, the agreement between experiment and computation validates the use of molecular dynamics simulations as an accurate tool for exploring the molecular motion of T cell receptors and their binding loops.

Disparate degrees of hypervariable loop flexibility control T-cell receptor cross-reactivity, specificity, and binding mechanism

αβ T-cell receptors (TCRs) recognize multiple antigenic peptides bound and presented by major histocompatibility complex molecules. TCR cross-reactivity has been attributed in part to the flexibility of TCR complementarity-determining region (CDR) loops, yet there have been limited direct studies of loop dynamics to determine the extent of its role. Here we studied the flexibility of the binding loops of the αβ TCR A6 using crystallographic, spectroscopic, and computational methods. A significant role for flexibility in binding and cross-reactivity was indicated only for the CDR3α and CDR3β hypervariable loops. Examination of the energy landscapes of these two loops indicated that CDR3β possesses a broad, smooth energy landscape, leading to rapid sampling in the free TCR of a range of conformations compatible with different ligands. The landscape for CDR3α is more rugged, resulting in more limited conformational sampling that leads to specificity for a reduced set of peptides as well as the major histocompatibility complex protein. In addition to informing on the mechanisms of cross-reactivity and specificity, the energy landscapes of the two loops indicate a complex mechanism for TCR binding, incorporating elements of both conformational selection and induced fit in a manner that blends features of popular models for TCR recognition.

Hydrolysis of acetylthiocoline, o-nitroacetanilide and o-nitrotrifluoroacetanilide by fetal bovine serum acetylcholinesterase

Besides esterase activity, acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) hydrolyze o-nitroacetanilides through aryl acylamidase activity. We have reported that BuChE tetramers and monomers of human blood plasma differ in o-nitroacetanilide (ONA) hydrolysis. The homology in quaternary structure and folding of subunits in the prevalent BuChE species (G4(H)) of human plasma and AChE forms of fetal bovine serum prompted us to study the esterase and amidase activities of fetal bovine serum AChE. The k(cat)/K(m) values for acetylthiocholine (ATCh), ONA and its trifluoro derivative N-(2-nitrophenyl)-trifluoroacetamide (F-ONA) were 398 x 10(6) M(-1) min(-1), 0.8 x 10(6) M(-1) min(-1), and 17.5 x 10(6) M(-1) min(-1), respectively. The lack of inhibition of amidase activity at high F-ONA concentrations makes it unlikely that there is a role for the peripheral anionic site (PAS) in F-ONA degradation, but the inhibition of ATCh, ONA and F-ONA hydrolysis by the PAS ligand fasciculin-2 points to the transit of o-nitroacetalinides near the PAS on their way to the active site. Sedimentation analysis confirmed substrate hydrolysis by tetrameric 10.9S AChE. As compared with esterase activity, amidase activity was less sensitive to guanidine hydrochloride. This reagent led to the formation of 9.3S tetramers with partially unfolded subunits. Their capacity to hydrolyze ATCh and F-ONA revealed that, despite the conformational change, the active site architecture and functionality of AChE were partially retained.

Parameters for Carbamate Pesticide QSAR and PBPK/PD Models for Human Risk Assessment

Our interest in providing parameters for the development of quantitative structure physiologically based pharmacokinetic/pharmacodynamic (QSPBPK/PD) models for assessing health risks to carbamates (USEPA 2005) comes from earlier work with organophosphorus (OP) insecticides (Knaak et al. 2004). Parameters specific to each carbamate are needed in the construction of PBPK/PD models along with their metabolic pathways. Parameters may be obtained by (1) development of QSAR models, (2) collecting pharmacokinetic data, and (3) determining pharmacokinetic parameters by fitting to experimental data. The biological parameters are given in Table 1 (Blancato et al. 2000). Table 1 Biological Parameters Required for Carbamate Pesticide Physiologically Based Pharmacokinetic/Pharmacodynamic (PBPK/PD) Models.(a).

Dynamics of the acetylcholinesterase tetramer

Acetylcholinesterase rapidly hydrolyzes the neurotransmitter acetylcholine in cholinergic synapses, including the neuromuscular junction. The tetramer is the most important functional form of the enzyme. Two low-resolution crystal structures have been solved. One is compact with two of its four peripheral anionic sites (PAS) sterically blocked by complementary subunits. The other is a loose tetramer with all four subunits accessible to solvent. These structures lacked the C-terminal amphipathic t-peptide (WAT domain) that interacts with the proline-rich attachment domain (PRAD). A complete tetramer model (AChEt) was built based on the structure of the PRAD/WAT complex and the compact tetramer. Normal mode analysis suggested that AChEt could exist in several conformations with subunits fluctuating relative to one another. Here, a multiscale simulation involving all-atom molecular dynamics and C alpha-based coarse-grained Brownian dynamics simulations was carried out to investigate the large-scale intersubunit dynamics in AChEt. We sampled the ns-mus timescale motions and found that the tetramer indeed constitutes a dynamic assembly of monomers. The intersubunit fluctuation is correlated with the occlusion of the PAS. Such motions of the subunits "gate" ligand-protein association. The gates are open more than 80% of the time on average, which suggests a small reduction in ligand-protein binding. Despite the limitations in the starting model and approximations inherent in coarse graining, these results are consistent with experiments which suggest that binding of a substrate to the PAS is only somewhat hindered by the association of the subunits.

Conformational dynamics in loop swap mutants of homologous fibronectin type III domains

Fibronectin type III (FN-III) domains are autonomously folded modules found in a variety of multidomain proteins. The 10th FN-III domain from fibronectin (fnFN10) and the 3rd FN-III domain from tenascin-C (tnFN3) have 27% sequence identity and the same overall fold; however, the CC' loop has a different pattern of backbone hydrogen bonds and the FG loop is longer in fnFN10 compared to tnFN3. To examine the influence of length, sequence, and context in determining dynamical properties of loops, CC' and FG loops were swapped between fnFN10 and tnFN3 to generate four mutant proteins and backbone conformational dynamics on ps-ns and mus-ms timescales were characterized by solution (15)N-NMR spin relaxation spectroscopy. The grafted loops do not strongly perturb the properties of the protein scaffold; however, specific effects of the mutations are observed for amino acids that are proximal in space to the sites of mutation. The amino acid sequence primarily dictates conformational dynamics when the wild-type and grafted loop have the same length, but both sequence and context contribute to conformational dynamics when the loop lengths differ. The results suggest that changes in conformational dynamics of mutant proteins must be considered in both theoretical studies and protein design efforts.

A model study of interactions between TcAChE peripheral site segment Tyr70Val71 and loop 1 of Fasciculin 2

The continuous chain of residues (Thr7 to Ala12) of Loop1 of Fas2 (F1) and its interaction with the peripheral binding sites (Tyr70-Val71) of AChE (P1) has been studied. Our results suggest that the flexibility of Loop1 might be caused by either the partially protonated guanidine group of Arg11 under experimental conditions or by the interaction with the negatively charged center of substrates. The binding energy of F1-P1 is predicted to be -16.6 kcal/mol at the B3LYP/6-311G(d,p) level, which is assumed to originate from one isolated O7...HN10 H-bond, one possible O10...HC71 unconventional O...HC type H-bonding, and the improved pi-bonding cooperativity around the peptide group of the AChE segment Tyr70-Val71. The classical Kitaura-Morokuma energy decomposition analysis, the NPA charge analysis, and the AIM analysis consistently reveal that the peptide group in segment P1 is more polarizable, which might play the key role in the interactions between F1 and P1. The PCM solvent effect corrected results reveal decrease of the interaction energy of the considered model. The importance of Thr8 of Fas2 in the P-site binding of AChE is also concluded. Site-directed mutations on either the Fas2 residue of Thr8 or the AChE residue of Tyr70 are expected to alter the binding behavior of the Loop1 of Fas2 with AChE.

Influence of agonists and antagonists on the segmental motion of residues near the agonist binding pocket of the acetylcholine-binding protein

Using the Lymnaea acetylcholine-binding protein as a surrogate of the extracellular domain of the nicotinic receptor, we combined site-directed labeling with fluorescence spectroscopy to assess possible linkages between ligand binding and conformational dynamics. Specifically, 2-[(5-fluoresceinyl)aminocarbonyl]ethyl methanethiosulfonate was conjugated to a free cysteine on loop C and to five substituted cysteines at strategic locations in the subunit sequence, and the backbone flexibility around each site of conjugation was measured with time-resolved fluorescence anisotropy. The sites examined were in loop C (Cys-188 using a C187S mutant), in the beta9 strand (T177C), in the beta10 strand (D194C), in the beta8-beta9 loop (N158C and Y164C), and in the beta7 strand (K139C). Conjugated fluorophores at these locations show distinctive anisotropy decay patterns indicating different degrees of segmental fluctuations near the agonist binding pocket. Ligand occupation and decay of anisotropy were assessed for one agonist (epibatidine) and two antagonists (alpha-bungarotoxin and d-tubocurarine). The Y164C and Cys-188 conjugates were also investigated with additional agonists (nicotine and carbamylcholine), partial agonists (lobeline and 4-hydroxy,2-methoxy-benzylidene anabaseine), and an antagonist (methyllycaconitine). With the exception of the T177C conjugate, both agonists and antagonists perturbed the backbone flexibility of each site; however, agonist-selective changes were only observed at Y164C in loop F where the agonists and partial agonists increased the range and/or rate of the fast anisotropy decay processes. The results reveal that agonists and antagonists produced distinctive changes in the flexibility of a portion of loop F.

Substrate and Product Trafficking through the Active Center Gorge of Acetylcholinesterase Analyzed by Crystallography and Equilibrium Binding

Hydrolysis of acetylcholine catalyzed by acetylcholinesterase (AChE), one of the most efficient enzymes in nature, occurs at the base of a deep and narrow active center gorge. At the entrance of the gorge, the peripheral anionic site provides a binding locus for allosteric ligands, including substrates. To date, no structural information on substrate entry to the active center from the peripheral site of AChE or its subsequent egress has been reported. Complementary crystal structures of mouse AChE and an inactive mouse AChE mutant with a substituted catalytic serine (S203A), in various complexes with four substrates (acetylcholine, acetylthiocholine, succinyldicholine, and butyrylthiocholine), two non-hydrolyzable substrate analogues (m-(N,N,N-trimethylammonio)-trifluoroacetophenone and 4-ketoamyltrimethylammonium), and one reaction product (choline) were solved in the 2.05-2.65-A resolution range. These structures, supported by binding and inhibition data obtained on the same complexes, reveal the successive positions and orientations of the substrates bound to the peripheral site and proceeding within the gorge toward the active site, the conformations of the presumed transition state for acylation and the acyl-enzyme intermediate, and the positions and orientations of the dissociating and egressing products. Moreover, the structures of the AChE mutant in complexes with acetylthiocholine and succinyldicholine reveal additional substrate binding sites on the enzyme surface, distal to the gorge entry. Hence, we provide a comprehensive set of structural snapshots of the steps leading to the intermediates of catalysis and the potential regulation by substrate binding to various allosteric sites at the enzyme surface.

Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus

NADH oxidase (NOX) from Thermus thermophilus is a member of a structurally homologous flavoprotein family of nitroreductases and flavin reductases. The importance of local conformational dynamics in the active site of NOX has been recently demonstrated. The enzyme activity was increased by 250% in the presence of 1 M urea with no apparent perturbation of the native structure of the protein. The present in silico results correlate with the in vitro data and suggest the possible explanation about the effect of urea on NOX activity at the molecular level. Both, X-ray structure and molecular dynamics (MD) simulations, show open conformation of the active site represented by approximately 0.9 nm distance between the indole ring of Trp47 and the isoalloxazine ring of FMN412. In this conformation, the substrate molecule can bind in the active site without sterical restraints. MD simulations also indicate more stable conformation of the active site called "closed" conformation. In this conformation, Trp47 and the isoalloxazine ring of FMN412 are so close to each other (approximately 0.5 nm) that the substrate molecule is unable to bind between them without perturbing this conformation. The open/close transition of the active site between Trp47 and the flavin ring is accompanied by release of the "tightly" bound water molecule from the active site--cofactor assisted gating mechanism. The presence of urea in aqueous solutions of NOX prohibits closing of the active site and even unlocks the closed active site because of the concomitant binding of a urea molecule in the active site cavity. The binding of urea in the active site is stabilized by formation of one/two persistent hydrogen bonds involving the carbonyl group of the urea molecule. Our report represents the first MD study of an enzyme from the novel flavoprotein family of nitroreductases and flavin reductases. The common occurrence of aromatic residues covering the active sites in homologous enzymes suggests the possibility of a general gating mechanism and the importance of local dynamics within this flavoprotein family.

Effect of dendrimers on pure acetylcholinesterase activity and structure

The effect of polyamidoamine (PAMAM) dendrimers on activity and fluorescence of pure acetylcholinesterase (EC 3.1.1.7.) was studied. It has been shown that all dendrimers studied decreased the enzymatic activity of acetylcholinesterase. This effect depended on the type of dendrimers. The data on the intrinsic fluorescence have shown that the dendrimers changed acetylcholinesterase conformation and the strongest effect was induced by PAMAM G3.5 dendrimer.

Molecular docking study on the “back door” hypothesis for product clearance in acetylcholinesterase

Acetylcholinesterase (AChE) is one of the fastest enzymes known, even though the active site is buried inside the protein at the end of a 20-A deep narrow gorge. Among the great variety of crystal structures of this enzyme, both in the absence and presence of various ligands and proteins, the structure of a complex of AChE with the pseudo-irreversible inhibitor Mf268 is of particular interest, as it assists in the proposal of a back door for product clearance from the active site. Binding of Mf268 to AChE results in the carbamoylation of Ser200 and liberation of an eseroline-fragment as the leaving group. The crystal structure of the AChE-Mf268 complex, however, proves that eseroline has escaped from the enzyme, despite the fact that the Ser-bound inhibitor fragment blocks the gorge entrance. The existence of alternative routes other than through the gorge for product clearance has been postulated but is still controversially discussed in the literature, as an experimental proof for such a back door is still missing. We have used Monte Carlo-based molecular docking methods in order to examine possible alternative pathways that could allow eseroline to be released from the protein after being cleaved from the substrate by Ser200. Based on our results, a short channel at the bottom of the gorge seems to be the most probable back-door site, which begins at amino acid Trp84 and ends at the enzyme surface in a cavity close to amino acid Glu445. [Figure: see text].