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].     

Backdoor opening mechanism in acetylcholinesterase based on X-ray crystallography and molecular dynamics simulations

The transient opening of a backdoor in the active‐site wall of acetylcholinesterase, one of nature's most rapid enzymes, has been suggested to contribute to the efficient traffic of substrates and products. A crystal structure of Torpedo californica acetylcholinesterase in complex with the peripheral‐site inhibitor aflatoxin is now presented, in which a tyrosine at the bottom of the active‐site gorge rotates to create a 3.4‐Å wide exit channel. Molecular dynamics simulations show that the opening can be further enlarged by movement of Trp84. The crystallographic and molecular dynamics simulation data thus point to the interface between Tyr442 and Trp84 as the key element of a backdoor, whose opening permits rapid clearance of catalysis products from the active site. Furthermore, the crystal structure presented provides a novel template for rational design of inhibitors and reactivators, including anti‐Alzheimer drugs and antidotes against organophosphate poisoning.     

Structure-Function Studies of a Melanocarpus albomyces Laccase Suggest a Pathway for Oxidation of Phenolic Compounds

Melanocarpus albomyces laccase crystals were soaked with 2,6-dimethoxyphenol, a common laccase substrate. Three complex structures from different soaking times were solved. Crystal structures revealed the binding of the original substrate and adducts formed by enzymatic oxidation of the substrate. The dimeric oxidation products were identified by mass spectrometry. In the crystals, a 2,6-dimethoxy-p-benzoquinone and a C-O dimer were observed, whereas a C-C dimer was the main product identified by mass spectrometry. Crystal structures demonstrated that the substrate and/or its oxidation products were bound in the pocket formed by residues Ala191, Pro192, Glu235, Leu363, Phe371, Trp373, Phe427, Leu429, Trp507 and His508. Substrate and adducts were hydrogen-bonded to His508, one of the ligands of type 1 copper. Therefore, this surface-exposed histidine most likely has a role in electron transfer by laccases. Based on our mutagenesis studies, the carboxylic acid residue Glu235 at the bottom of the binding site pocket is also crucial in the oxidation of phenolics. Glu235 may be responsible for the abstraction of a proton from the OH group of the substrate and His508 may extract an electron. In addition, crystal structures revealed a secondary binding site formed through weak dimerization in M. albomyces laccase molecules. This binding site most likely exists only in crystals, when the Phe427 residues are packed against each other.     

Substitutions at the opening of the Rubisco central solvent channel affect holoenzyme stability and CO2/O2 specificity but not activation by Rubisco activase

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) catalyzes the initial step of carbon metabolism in photosynthesis. The holoenzyme comprises eight large subunits, arranged as a tetramer of dimers around a central solvent channel that defines a fourfold axis of symmetry, and eight small subunits, arranged as two tetramers at the poles of the axis. The phylogenetically divergent small-subunit loops between β-strands A and B form the entrance to the solvent channel. In the green alga Chlamydomonas reinhardtii, Ile-58 from each of the four small-subunit βA–βB loops defines the minimal diameter of the channel opening. To understand the role of the central solvent channel in Rubisco function, directed mutagenesis and transformation of Chlamydomonas were employed to replace Ile-58 with Ala, Lys, Glu, Trp, or three Trp residues (I58W3) to close the entrance to the channel. The I58E, I58K, and I58W substitutions caused only small decreases in photosynthetic growth at 25 and 35 °C, whereas I58W3 had a substantial effect at both temperatures. The mutant enzymes had decreased carboxylation rates, but the I58W3 enzyme had decreases in both carboxylation and CO₂/O₂ specificity. The I58E, I58W, and I58W3 enzymes were inactivated at lower temperatures than wild-type Rubisco, and were degraded at slower rates under oxidative stress. However, these mutant enzymes were activated by Rubisco activase at normal rates, indicating that the structural transition required for carboxylation is not affected by altering the solvent channel opening. Structural dynamics alone may not be responsible for these distant effects on the Rubisco active site.     

Influence of size and polarity of residue 31 in porcine pancreatic phospholipase A2 on catalytic properties

Residue 31 of porcine pancreatic phospholipase A2 (PLA2) is located at the entrance to the active site. To study the role of residue 31 in PLA2, six mutant enzymes were produced by site-directed mutagenesis, replacing Leu by either Trp, Arg, Ala, Thr, Ser or Gly. Direct binding studies indicated a three to six times greater affinity of the Trp31 PLA2 for both monomeric and micellar substrate analogs, relative to the wild-type enzyme. The other five mutants possess an unchanged affinity for monomers of the product analog n-decylphosphocholine and for micelles of the diacyl substrate analog rac-1,2-dioctanoylamino-dideoxy-glycero-3-phosphocholine. The affinities for micelles of the monoacyl product analog n-hexadecylphosphocholine were decreased 9-20 times for these five mutants. Kinetic studies with monomeric substrates showed that the mutants have Vmax values which range between 15 and 70% relative to the wild-type enzyme. The Vmax values for micelles of the zwitterionic substrate 1,2-dioctanoyl-sn-glycero-3-phosphocholine were lowered 3-50 times. The Km values for the monomeric substrate and the Km values for the micellar substrate were hardly affected in the case of five of the six mutants, but were considerably decreased when Trp was present at position 31. The results of these investigations point to a versatile role for the residue at position 31: involvement in the binding and orientating of monomeric substrate (analogs), involvement in the binding of the enzyme to micellar substrate analogs and possibly involvement in shielding the active site from excess water.     

Unveiling Prolyl Oligopeptidase Ligand Migration by Comprehensive Computational Techniques

Prolyl oligopeptidase (POP) is a large 80 kDa protease, which cleaves oligopeptides at the C-terminal side of proline residues and constitutes an important pharmaceutical target. Despite the existence of several crystallographic structures, there is an open debate about migration (entrance and exit) pathways for ligands, and their coupling with protein dynamics. Recent studies have shown the capabilities of molecular dynamics and classical force fields in describing spontaneous binding events and nonbiased ligand migration pathways. Due to POP’s size and to the buried nature of its active site, an exhaustive sampling by means of conventional long enough molecular dynamics trajectories is still a nearly impossible task. Such a level of sampling, however, is possible with the breakthrough protein energy landscape exploration technique. Here, we present an exhaustive sampling of POP with a known inhibitor, Z-pro-prolinal. In >3000 trajectories Z-pro-prolinal explores all the accessible surface area, showing multiple entrance events into the large internal cavity through the pore in the β-propeller domain. Moreover, we modeled a natural substrate binding and product release by predicting the entrance of an undecapeptide substrate, followed by manual active site cleavage and nonbiased exit of one of the products (a dipeptide). The product exit shows preference from a flexible 18-amino acid residues loop, pointing to an overall mechanism where entrance and exit occur in different sites.     

Structural insights into substrate traffic and inhibition in acetylcholinesterase

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

Importance of van der Waals contact between Glu 35 and Trp 109 to the catalytic action of human lysozyme.

The importance of van der Waals contact between Glu 35 and Trp 109 to the active-site structure and the catalytic properties of human lysozyme (HL) has been investigated by site-directed mutagenesis. The X-ray analysis of mutant HLs revealed that both the replacement of Glu 35 by Asp or Ala, and the replacement of Trp 109 by Phe or Ala resulted in a significant but localized change in the active-site cleft geometry. A prominent movement of the backbone structure was detected in the region of residues 110 to 120 and in the region of residues 100 to 115 for the mutations concerning Glu 35 and Trp 109, respectively. Accompanied by the displacement of the main-chain atoms with a maximal deviation of C alpha atom position ranging from 0.7 A to 1.0 A, the mutant HLs showed a remarkable change in the catalytic properties against Micrococcus luteus cell substrate as compared with native HL. Although the replacement of Glu 35 by Ala completely abolished the lytic activity, HL-Asp 35 mutant retained a weak but a certain lytic activity, showing the possible involvement of the side-chain carboxylate group of Asp 35 in the catalytic action. The kinetic consequence derived from the replacement of Trp 109 by Phe or Ala together with the result of the structural change suggested that the structural detail of the cleft lobe composed of the residues 100 to 115 centered at Ala 108 was responsible for the turnover in the reaction of HL against the bacterial cell wall substrate. The results revealed that the van der Waals contact between Glu 35 and Trp 109 was an essential determinant in the catalytic action of HL.     

Mechanistic studies on N-acetylmuramic acid 6-phosphate hydrolase (MurQ): an etherase involved in peptidoglycan recycling

Peptidoglycan recycling is a process in which bacteria import cell wall degradation products and incorporate them back into either peptidoglycan biosynthesis or basic metabolic pathways. The enzyme MurQ is an N-acetylmuramic acid 6-phosphate (MurNAc 6-phosphate) hydrolase (or etherase) that hydrolyzes the lactyl side chain from MurNAc 6-phosphate and generates GlcNAc 6-phosphate. This study supports a mechanism involving the syn elimination of lactate to give an alpha,beta-unsaturated aldehyde with (E)-stereochemistry, followed by the syn addition of water to give product. The observation of both a kinetic isotope effect slowing the reaction of [2-(2)H]MurNAc 6-phosphate and the incorporation of solvent-derived deuterium into C2 of the product indicates that the C2-H bond is cleaved during catalysis. The observation that the solvent-derived (18)O isotope is incorporated into the C3 position of the product, but not the C1 position, provides evidence of the cleavage of the C3-O bond and argues against imine formation. The finding that 3-chloro-3-deoxy-GlcNAc 6-phosphate serves as an alternate substrate is also consistent with an elimination-addition mechanism. Upon extended incubations of MurQ with GlcNAc 6-phosphate, the alpha,beta-unsaturated aldehydic intermediate accumulates in solution, and (1)H NMR analysis indicates it exists as the ring-closed form of the (E)-alkene. A structural model is developed for the Escherichia coli MurQ and is compared to that of the structural homologue glucosamine-6-phosphate synthase. Putative active site acid/base residues are probed by mutagenesis, and Glu83 and Glu114 are found to be crucial for catalysis. The Glu83Ala mutant is essentially inactive as an etherase yet is capable of exchanging the C2 proton of substrate with solvent-derived deuterium. This suggests that Glu83 may function as the acidic residue that protonates the departing lactate.     

Energetic and mechanistic studies of glucoamylase using molecular recognition of maltose OH groups coupled with site-directed mutagenesis

Molecular recognition using a series of deoxygenated maltose analogues was used to determine the substrate transition-state binding energy profiles of 10 single-residue mutants at the active site of glucoamylase from Aspergillus niger. The individual contribution of each substrate hydroxyl group to transition-state stabilization with the wild type and each mutant GA was determined from the relation Delta(DeltaG()) = -RT ln[(k(cat)/K(M))(x)/(k(cat)/K(M))(y)], where x represents either a mutant enzyme or substrate analogue and y the wild-type enzyme or parent substrate. The resulting binding energy profiles indicate that disrupting an active site hydrogen bond between enzyme and substrate, as identified in crystal structures, not only sharply reduces or eliminates the energy contributed from that particular hydrogen bond but also perturbs binding contributions from other substrate hydroxyl groups. Replacing the active site acidic groups, Asp55, Glu180, or Asp309, with the corresponding amides, and the neutral Trp178 with the basic Arg, all substantially reduced the binding energy contribution of the 4'- and 6'-OH groups of maltose at subsite -1, even though both Glu180 and Asp309 are localized at subsite 1. In contrast, the substitution, Asp176 --> Asn, located near subsites -1 and 1, did not substantially perturb any of the individual hydroxyl group binding energies. Similarly, the substitutions Tyr116 --> Ala, Ser119 --> Tyr, or Trp120 --> Phe also did not substantially alter the energy profiles even though Trp120 has a critical role in directing conformational changes necessary for activity. Since the mutations at Trp120 and Asp176 reduced k(cat) values by 50- and 12-fold, respectively, a large effect on k(cat) is not necessarily accompanied by changes in hydroxyl group binding energy contributions. Two substitutions, Asn182 --> Ala and Tyr306 --> Phe, had significant though small effects on interactions with 3- and 4'-OH, respectively. Binding interactions between the enzyme and the glucosyl group in subsite -1, particularly with the 4'- and 6'-OH groups, play an important role in substrate binding, while subsite 1 interactions may play a more important role in product release.     

Engineering of the pH-dependence of thermolysin activity as examined by site-directed mutagenesis of Asn112 located at the active site of thermolysin

Asn112 is located at the active site of thermolysin, 5-8 A from the catalytic Zn2+ and catalytic residues Glu143 and His231. When Asn112 was replaced with Ala, Asp, Glu, Lys, His, and Arg by site-directed mutagenesis, the mutant enzymes N112D and N112E, in which Asn112 is replaced with Asp and Glu, respectively, were secreted as an active form into Escherichia coli culture medium, while the other four were not. In the hydrolysis of a neutral substrate N-[3-(2-furyl)acryloyl]-Gly-L-Leu amide, the kcat/Km values of N112D and N112E exhibited bell-shaped pH-dependence, as did the wild-type thermolysin (WT). The acidic pKa of N112D was 5.7 +/- 0.1, higher by 0.4 +/- 0.2 units than that of WT, suggesting that the introduced negative charge suppressed the protonation of Glu143 or Zn2+-OH. In the hydrolysis of a negatively charged substrate, N-carbobenzoxy-l-Asp-l-Phe methyl ester (ZDFM), the pH-dependence of kcat/Km of the mutants decreased with increase in pH from 5.5 to 8.5, while that of WT was bell-shaped. This difference might be explained by the electrostatic repulsion between the introduced Asp/Glu and ZDFM, suggesting that introducing ionizing residues into the active site of thermolysin might be an effective means of modifying its pH-activity profile.     

A multinuclear NMR study of the active site of an endoglucanase from a strain of Bacillus. Use of Trp residues as structural probes.

In the hydrolytic reaction catalyzed by an endoglucanase from a Bacillus strain (endoglucanase K), 2 of 12 Trp residues, Trp174 and Trp243, are responsible for binding of the substrate and/or for the catalysis (Kawaminami, S., Ozaki, K., Sumitomo, N., Hayashi, Y., Ito, S., Shimada, I., and Arata, Y. (1994) J. Biol. Chem. 269, 28752-28756). Here we report results of a stable isotope-aided NMR analysis of the active site of endoglucanase K, using Trp174 and Trp243 as structural probes. Hydrogen-deuterium exchange experiments performed for the NH protons of main and side chains of Trp residues revealed that Trp174 and Trp243 are located in the hydrophilic and hydrophobic microenvironments in the active site, respectively. We also carried out pH titration experiments for indole C2 proton resonances of Trp residues and measured the pH dependence of specific activities for wild-type endoglucanase K and its mutants in which Glu or Asp residues are replaced with their respective amide forms. On the basis of the results obtained from the present study, we conclude that (a) Glu130 and Asp191, which are in spatial proximity to Trp174 and Trp243 in the active site, play a crucial role in the enzymatic activity; (b) Glu130 and Asp191 interact with each other in the active site, leading to an increase in the pKa values to 5.5 for both amino acid residues; and (c) the pKa values of Glu130 and Asp191 would lead to an unusually narrow pH-activity profile of the endoglucanase K.     

Aromatic stacking in the sugar binding site of the lactose permease

Major determinants for substrate recognition by the lactose permease of Escherichia coli are at the interface between helices IV (Glu126, Ala122), V (Arg144, Cys148), and VIII (Glu269). We demonstrate here that Trp151, one turn of helix V removed from Cys148, also plays an important role in substrate binding probably by aromatic stacking with the galactopyranosyl ring. Mutants with Phe or Tyr in place of Trp151 catalyze active lactose transport with time courses nearly the same as wild type. In addition, apparent K(m) values for lactose transport in the Phe or Tyr mutants are only 6- or 3-fold higher than wild type, respectively, with a comparable V(max). Surprisingly, however, binding of high-affinity galactoside analogues is severely compromised in the mutants; the affinity of mutant Trp151-->Phe or Trp151-->Tyr is diminished by factors of at least 50 or 20, respectively. The results demonstrate that Trp151 is an important component of the binding site, probably orienting the galactopyranosyl ring so that important H-bond interactions with side chains in helices IV, V, and VIII can be realized. The results are discussed in the context of a current model for the binding site.     

The apolar channel in Cerebratulus lacteus hemoglobin is the route for O2 entry and exit

The major pathway for O2 binding to mammalian myoglobins (Mb) and hemoglobins (Hb) involves transient upward movement of the distal histidine (His-64(E7)), allowing ligand capture in the distal pocket. The mini-globin from Cerebratulus lacteus (CerHb) appears to have an alternative pathway between the E and H helices that is made accessible by loss of the N-terminal A helix. To test this pathway, we examined the effects of changing the size of the E7 gate and closing the end of the apolar channel in CerHb by site-directed mutagenesis. Increasing the size of Gln-44(E7) from Ala to Trp causes variation of association (k'O2) and dissociation (kO2) rate coefficients, but the changes are not systematic. More significantly, the fractions (Fgem approximately 0.05-0.19) and rates (kgem approximately 50-100 micros(-1)) of geminate CO recombination in the Gln-44(E7) mutants are all similar. In contrast, blocking the entrance to the apolar channel by increasing the size of Ala-55(E18) to Phe and Trp causes the following: 1) both k'O2 and kO2 to decrease roughly 4-fold; 2) Fgem for CO to increase from approximately 0.05 to 0.45; and 3) kgem to decrease from approximately 80 to approximately 9 micros(-1), as ligands become trapped in the channel. Crystal structures and low temperature Fourier-transform infrared spectra of Phe-55 and Trp-55 CerHb confirm that the aromatic side chains block the channel entrance, with little effect on the distal pocket. These results provide unambiguous experimental proof that diatomic ligands can enter and exit a globin through an interior channel in preference to the more direct E7 pathway.     

Molecular mechanism of ADP-ribose hydrolysis by human NUDT5 from structural and kinetic studies

Human NUDT5 (hNUDT5) is an ADP-ribose (ADPR) pyrophosphatase (ADPRase) that plays important roles in controlling the intracellular levels of ADPR and preventing non-enzymatic ADP-ribosylation of proteins by hydrolyzing ADPR to AMP and ribose 5′-phosphate. We report the crystal structure of hNUDT5 in complex with a non-hydrolyzable ADPR analogue, α,β-methyleneadenosine diphosphoribose, and three Mg^2 + ions representing the transition state of the enzyme during catalysis. Analysis of this structure and comparison with previously reported hNUDT5 structures identify key residues involved in substrate binding and catalysis. In the transition-state structure, three metal ions are bound at the active site and are coordinated by surrounding residues and water molecules. A conserved water molecule is at an ideal position for nucleophilic attack on the α-phosphate of ADPR. The side chain of Glu166 on loop L9 changes its conformation to interact with the conserved water molecule compared with that in the substrate-bound structure and appears to function as a catalytic base. Mutagenesis and kinetic studies show that Trp28 and Trp46 are important for the substrate binding; Arg51 is involved in both the substrate binding and the catalysis; and Glu112 and Glu116 of the Nudix motif, Glu166 on loop L9, and Arg111 are critical for the catalysis. The structural and biochemical data together reveal the molecular basis of the catalytic mechanism of ADPR hydrolysis by hNUDT5. Specifically, Glu166 functions as a catalytic base to deprotonate a conserved water molecule that acts as a nucleophile to attack the α-phosphate of ADPR, and three Mg^2 + ions are involved in the activation of the nucleophile and the binding of the substrate. Structural comparison of different ADPRases also suggests that most dimeric ADPRases may share a similar catalytic mechanism of ADPR hydrolysis.     

Molecular dynamics simulation of the processive endocellulase Cel48F from Clostridium cellulolyticum: a novel “water‐control mechanism” in enzymatic hydrolysis of cellulose

Glycoside hydrolase of Cel48F from Clostridium cellulolyticum is an important processive cellulose, which can hydrolyze cellulose into cellobiose. Molecular dynamics simulations were used to investigate the hydrolysis mechanism of cellulose. The two conformations of the Cel48F‐cellotetrose complex in which the cellotetroses are bound at different sites (known as the sliding conformation and the hydrolyzing conformation) were simulated. By comparing these two conformations, a water‐control mechanism is proposed, in which the hydrolysis proceeds by providing a water molecule for every other glucosidic linkage. The roles of certain key residues are determined: Glu55 and Asp230 are the most probable candidates for acid and base, respectively, in the mechanism of inverting anomeric carbon. Met414 and Trp417 constitute the water‐control system. Glu44 might keep the substrate at a certain location within the active site or help the substrate chain to move from the sliding conformation to the hydrolyzing conformation. The other hydrophobic residues around the substrate can decrease the sliding energy barrier or provide a hydrophobic environment to resist entry of the surrounding water molecules into the active site, except for those coming from a specific water channel. Copyright © 2014 John Wiley & Sons, Ltd.     

Complexes of mutants of Escherichia coli aminopeptidase P and the tripeptide substrate ValProLeu

Aminopeptidase P (APPro) is a manganese-containing enzyme that catalyses the hydrolysis of the N-terminal residue of a polypeptide if the second residue is proline. Structures of APPro mutants with reduced or negligible activity have been determined in complex with the tripeptide substrate ValProLeu. In the complex of Glu383Ala APPro with ValProLeu one of the two metal sites is only partly occupied, indicating an essential role for Glu383 in metal binding in the presence of substrate. His361Ala APPro clearly possesses residual activity as the ValProLeu substrate has been cleaved in the crystals; difference electron density consistent with bound ProLeu dipeptide and a disordered Val amino acid is present at the active site. Contrary to previous suggestions, the His243Ala mutant is capable of binding substrate. The structure of the His243Ala APPro complex with ValProLeu shows that the peptide interacts with one of the active-site metal atoms via its terminal amino group. The implications of these complexes for the roles of the respective residues in APPro catalysis are discussed.     

Characterization of the channel constriction allowing the access of the substrate to the active site of yeast oxidosqualene cyclase

In oxidosqualene cyclases (OSCs), an enzyme which has been extensively studied as a target for hypocholesterolemic or antifungal drugs, a lipophilic channel connects the surface of the protein with the active site cavity. Active site and channel are separated by a narrow constriction operating as a mobile gate for the substrate passage. In Saccharomyces cerevisiae OSC, two aminoacidic residues of the channel/constriction apparatus, Ala525 and Glu526, were previously showed as critical for maintaining the enzyme functionality. In this work sixteen novel mutants, each bearing a substitution at or around the channel constrictions, were tested for their enzymatic activity. Modelling studies showed that the most functionality-lowering substitutions deeply alter the H-bond network involving the channel/constriction apparatus. A rotation of Tyr239 is proposed as part of the mechanism permitting the access of the substrate to the active site. The inhibition of OSC by squalene was used as a tool for understanding whether the residues under study are involved in a pre-catalytic selection and docking of the substrate oxidosqualene.