Attack Conformation Research Articles

Computing the role of near attack conformations in an enzyme-catalyzed nucleophilic bimolecular reaction.

Near attack conformations (NACs) are conformations extending from the ground state (GS) that lie on the transition path of a chemical reaction. Here, we develop a method for computing the thermodynamic contribution to catalysis due to NAC formation in bimolecular reactions, within the limit of a classical molecular dynamics force field. We make use of the Bürgi–Dunitz theory applied to large-scale unbiased all-atom ensemble molecular dynamics simulations. We apply this to HIV-1 protease peptide hydrolysis, known to achieve a rate enhancement of ∼1011 (ΔGcat⧧ ∼ 15 kcal/mol) over the uncatalyzed bimolecular reaction (ΔGnon⧧ ∼ 30 kcal/mol). The ground state consists of a nucleophilic water molecule bound to an octapeptide substrate in the active site. We first observe multiple and reversible binding of a nucleophilic water molecule into the active site giving a free energy of binding of ΔG = −1 kcal/mol to form the GS. The free energy barriers for catalyzed and uncatalyzed NAC formation are both equivalent: ΔGNAC⧧ = 4.6 kcal/mol, constituting ∼30% and ∼15% of the overall barriers, respectively. Therefore, not only does adoption of NACs only account for minor progress along the transition path in both catalyzed and uncatalyzed reactions, but there is no preferential formation of them in the catalyzed reaction. Analysis of the catalytic hydrogen bond network reveals interactions that stabilize the GS; however, subsequent NAC formation does not preferentially favor any of the possible hydrogen bond configurations. This supports the view that the catalytic power of HIV-1 protease is not due to NAC formation.

Impact of Magnesium(II) on Beryllium Fluorides in Solutions Studied by 19F NMR Spectroscopy

AbstractBeryllium fluorides are widely used in protein phosphorylation studies to get stable transition state analogs or near attack conformers, which has attracted much attention. BeF3− is one of the optimal phosphoryl (PO3−) analogs for its identical geometry and charge, and Mg2+ naturally participates in the phosphoryl binding in biological systems. In solutions, BeF3− coexists with other beryllium fluorides (BeF42−, BeF2 and BeF+) and magnesium fluorides, and there are equilibriums between these species. In this article, 19F NMR spectroscopy was applied to the investigation of the impact of magnesium(II) on beryllium fluorides. It has been found that when Mg2+ was introduced into the solutions, the chemical shifts, the intensities and the line widths of 19F signals of various beryllium fluoride complexes were changed. After ionic strength correction, these effects were remarkable only for BeF42− and BeF3−, especially BeF42−, when the concentration of the fluoride ion is relatively low. Mechanism of the effects is proposed which involves ion pair formation between Mg2+ and beryllium fluorides.

Computationally efficient and accurate enantioselectivity modeling by clusters of molecular dynamics simulations.

Computational approaches could decrease the need for the laborious high-throughput experimental screening that is often required to improve enzymes by mutagenesis. Here, we report that using multiple short molecular dynamics (MD) simulations makes it possible to accurately model enantioselectivity for large numbers of enzyme-substrate combinations at low computational costs. We chose four different haloalkane dehalogenases as model systems because of the availability of a large set of experimental data on the enantioselective conversion of 45 different substrates. To model the enantioselectivity, we quantified the frequency of occurrence of catalytically productive conformations (near attack conformations) for pairs of enantiomers during MD simulations. We found that the angle of nucleophilic attack that leads to carbon-halogen bond cleavage was a critical variable that limited the occurrence of productive conformations; enantiomers for which this angle reached values close to 180° were preferentially converted. A cluster of 20-40 very short (10 ps) MD simulations allowed adequate conformational sampling and resulted in much better agreement to experimental enantioselectivities than single long MD simulations (22 ns), while the computational costs were 50-100 fold lower. With single long MD simulations, the dynamics of enzyme-substrate complexes remained confined to a conformational subspace that rarely changed significantly, whereas with multiple short MD simulations a larger diversity of conformations of enzyme-substrate complexes was observed.

Evidence for the role of active site residues in the hairpin ribozyme from molecular simulations along the reaction path.

The hairpin ribozyme accelerates a phosphoryl transfer reaction without catalytic participation of divalent metal ions. Residues A38 and G8 have been implicated as playing roles in general acid and base catalysis, respectively. Here we explore the structure and dynamics of key active site residues using more than 1 μs of molecular dynamics simulations of the hairpin ribozyme at different stages along the catalytic pathway. Analysis of results indicates hydrogen bond interactions between the nucleophile and proR nonbridging oxygen are correlated with active inline attack conformations. Further, the simulation results suggest a possible alternative role for G8 to promote inline fitness and facilitate activation of the nucleophile by hydrogen bonding, although this does not necessarily exclude an additional role as a general base. Finally, we suggest that substitution of G8 with N7- or N3-deazaguanosine which have elevated pKa values, both with and without thio modifications at the 5′ leaving group position, would provide valuable insight into the specific role of G8 in catalysis.

Ab initiodynamics trajectory study of the heterolytic cleavage of H2 by a Lewis acid [B(C6F5)3] and a Lewis base [P(tBu)3

Activation of H2 by a "frustrated Lewis pair" (FLP) composed of B(C6F5)3 and P(tBu)3 species has been explored with high level direct ab initio molecular dynamics (AIMD) simulations at finite temperature (T = 300 K) in gas phase. The initial geometrical conditions for the AIMD trajectory calculations, i.e., the near attack conformations of FLP + H2, were devised using the host-guest model in which suitable FLP conformations were obtained from the dynamics of the B(C6F5)3∕P(tBu)3 pair in gas phase. AIMD trajectory calculations yielded microscopic insight into effects which originate from nuclear motion in the reacting complex, e.g., the alternating compression∕elongation of the boron-phosphorous distance and the change of the pyramidality of boron in B(C6F5)3. The ensemble averaged trajectory analysis has been compared with the minimum energy path (MEP) description of the reaction. Similar to MEP, AIMD shows that an attack of the acid∕base pair on the H-H bond gives rise to the polarization of the H2 molecule and as a consequence generates a large dipole moment of the reacting complex. The MEP and AIMD portrayals of the reaction are fundamentally different in terms of the magnitude of the motion of nuclei in B(C6F5)3 and P(tBu)3 during the H2 cleavage. In the AIMD trajectory simulations, geometries of B(C6F5)3 and P(tBu)3 appear as nearly "frozen" on the short time scale of the H2 cleavage. This is contrary to the MEP picture. Several of the concepts which arise from this work, e.g., separation of time scales of nuclear motion and the time-dependence of the donor-acceptor interactions in the reacting complex, are important for the understanding of chemical reactivity and catalysis.

Water assisted reaction mechanism of OH− with CCl4 in aqueous solution – Hybrid quantum mechanical and molecular mechanics investigation

The OH−(H2O)+CCl4 reaction in aqueous solution was investigated using the combined quantum mechanical and molecular mechanics approach. The reaction mechanism of OH−(H2O)+CCl4 consists of two concerted steps – formation of OH− in the favorable attack conformation via the proton transfer process, and the nucleophilic substitution process in which the newly formed OH− attacks the CCl4. The free energy activation barrier is 38.2kcal/mol at CCSD(T)/MM level of theory for this reaction, which is about 10.3kcal/mol higher than that of the direct nucleophilic substitution mechanism of the OH−+CCl4 reaction in aqueous solution.

Reflections on biocatalysis involving phosphorus

Early studies on chemical synthesis of biological molecules can be seen to progress to preparation and biological evaluation of phosphonates as analogues of biological phosphates, with emphasis on their isosteric and isopolar character. Work with such mimics progressed into structural studies with a range of nucleotide-utilising enzymes. The arrival of metal fluorides as analogues of the phosphoryl group, PO(3)(-), for transition state (TS) analysis of enzyme reactions stimulated the symbiotic deployment of (19)F NMR and protein crystallography. Characteristics of enzyme transition state analogues are reviewed for a range of reactions. From the available MF(x) species, trifluoroberyllate gives tetrahedral mimics of ground states (GS) in which phosphate is linked to carboxylate and phosphate oxyanions. Tetrafluoroaluminate is widely employed as a TS mimic, but it necessarily imposes octahedral geometry on the assembled complexes, whereas phosphoryl transfer involves trigonal bipyramidal (tbp) geometry. Trifluoromagnesate (MgF(3)(-)) provides the near-ideal solution, delivering tbp geometry and correct anionic charge. Some of the forty reported tbp structures assigned as having AlF(3)(0) cores have been redefined as trifluoromagnesate complexes. Transition state analogues for a range of kinases, mutases, and phosphatases provide a detailed description of mechanism for phosphoryl group transfer, supporting the concept of charge balance in their TS and of concerted-associative pathways for biocatalysis. Above all, superposition of GS and TS structures reveals that in associative phosphoryl transfer, the phosphorus atom migrates through a triangle of three, near-stationary, equatorial oxygens. The extension of these studies to near attack conformers further illuminates enzyme catalysis of phosphoryl transfer.

Initial Events in the Degradation of Hyaluronan Catalyzed by Hyaluronate Lyase from Spectrococcus pneumoniae: QM/MM Simulation

Hyaluronate lyase from Spectrococcus pneumonia can degrade hyaluronic acid, which is one of the major components in the extracellular matrix. The major functions of hyaluronan are to regulate water balance and osmotic pressure and act as an ion-exchange resin. In this work, we focus on the prerequisite issue of the enzymatic reaction, i.e., the initial reactive conformer. Based on the quantum mechanical and molecular mechanical molecular dynamic simulations and free energy profiles, a near attack conformer was obtained for the degradation of hyaluronan catalyzed by the hyaluronate lyase. Along with the substrate binding, the phenylhydroxyl hydrogen atom of Tyr408 will transfer to nearby His399 via a near barrierless transition state, which results in a negatively charged Tyr408 and positively charged His399. The Tyr408, rather than the previously proposed His399, was suggested to act as the general base for the subsequent β-elimination reaction. The His399 was suggested to have the function of neutralizing the C5-carboxyl group.

Substrate Specificity of Fluoroacetate Dehalogenase: An Insight from Crystallographic Analysis, Fluorescence Spectroscopy, and Theoretical Computations

The high substrate specificity of fluoroacetate dehalogenase was explored by using crystallographic analysis, fluorescence spectroscopy, and theoretical computations. A crystal structure for the Asp104Ala mutant of the enzyme from Burkholderia sp. FA1 complexed with fluoroacetate was determined at 1.2 Å resolution. The orientation and conformation of bound fluoroacetate is different from those in the crystal structure of the corresponding Asp110Asn mutant of the enzyme from Rhodopseudomonas palustris CGA009 reported recently (J. Am. Chem. Soc. 2011, 133, 7461). The fluorescence of the tryptophan residues of the wild-type and Trp150Phe mutant enzymes from Burkholderia sp. FA1 incubated with fluoroacetate and chloroacetate was measured to gain information on the environment of the tryptophan residues. The environments of the tryptophan residues were found to be different between the fluoroacetate- and chloroacetate-bound enzymes; this would come from different binding modes of these two substrates in the active site. Docking simulations and QM/MM optimizations were performed to predict favorable conformations and orientations of the substrates. The F atom of the substrate is oriented toward Arg108 in the most stable enzyme-fluoroacetate complex. This is a stable but unreactive conformation, in which the small O-C-F angle is not suitable for the S(N)2 displacement of the F(-) ion. The cleavage of the C-F bond is initiated by the conformational change of the substrate to a near attack conformation (NAC) in the active site. The second lowest energy conformation is appropriate for NAC; the C-O distance and the O-C-F angle are reasonable for the S(N) 2 reaction. The activation energy is greatly reduced in this conformation because of three hydrogen bonds between the leaving F atom and surrounding amino acid residues. Chloroacetate cannot reach the reactive conformation, due to the longer C-Cl bond; this results in an increase of the activation energy despite the weaker C-Cl bond.

Near attack conformers dominate β-phosphoglucomutase complexes where geometry and charge distribution reflect those of substrate

Experimental observations of fluoromagnesate and fluoroaluminate complexes of β-phosphoglucomutase (β-PGM) have demonstrated the importance of charge balance in transition-state stabilization for phosphoryl transfer enzymes. Here, direct observations of ground-state analog complexes of β-PGM involving trifluoroberyllate establish that when the geometry and charge distribution closely match those of the substrate, the distribution of conformers in solution and in the crystal predominantly places the reacting centers in van der Waals proximity. Importantly, two variants are found, both of which satisfy the criteria for near attack conformers. In one variant, the aspartate general base for the reaction is remote from the nucleophile. The nucleophile remains protonated and forms a nonproductive hydrogen bond to the phosphate surrogate. In the other variant, the general base forms a hydrogen bond to the nucleophile that is now correctly orientated for the chemical transfer step. By contrast, in the absence of substrate, the solvent surrounding the phosphate surrogate is arranged to disfavor nucleophilic attack by water. Taken together, the trifluoroberyllate complexes of β-PGM provide a picture of how the enzyme is able to organize itself for the chemical step in catalysis through the population of intermediates that respond to increasing proximity of the nucleophile. These experimental observations show how the enzyme is capable of stabilizing the reaction pathway toward the transition state and also of minimizing unproductive catalysis of aspartyl phosphate hydrolysis.

Theoretical investigation of experimentally determined α-εffects — The role of electronics

Recent experimental reports involving both α-nucleophiles and normal nucleophiles have reported both the presence and absence of an α-effect. In ester systems, such as dimethylmethylphosphonate (DMMP), a small α-effect is reported, but the reference point is a stationary point of the potential energy surface that must rearrange to acquire the near attack conformation (NAC) necessary for the Sn2 pathway to proceed. The second type of study involves use of highly fluorinated alkoxides as normal nucleophiles and reports no α-effect. This paper employs linear free energy plots in an investigation of electronic effects in methyl formate SN2 reactions, using high-level computations of transition states for determination of energy barriers.

Protonation states of the key active site residues and structural dynamics of the glmS riboswitch as revealed by molecular dynamics.

The glmS catalytic riboswitch is part of the 5'-untranslated region of mRNAs encoding glucosamine-6-phosphate (GlcN6P) synthetase (glmS) in numerous gram-positive bacteria. Binding of the cofactor GlcN6P induces site-specific self-cleavage of the RNA. However, the detailed reaction mechanism as well as the protonation state of the glmS reactive form still remains elusive. To probe the dominant protonation states of key active site residues, we carried out explicit solvent molecular dynamic simulations involving various protonation states of three crucial active site moieties observed in the available crystal structures: (i) guanine G40 (following the Thermoanaerobacter tengcongensis numbering), (ii) the GlcN6P amino/ammonium group, and (iii) the GlcN6P phosphate moiety. We found that a deprotonated G40(-) seems incompatible with the observed glmS active site architecture. Our data suggest that the canonical form of G40 plays a structural role by stabilizing an in-line attack conformation of the cleavage site A-1(2'-OH) nucleophile, rather than a more direct chemical role. In addition, we observe weakened cofactor binding upon protonation of the GlcN6P phosphate moiety, which explains the experimentally observed increase in K(m) with decreasing pH. Finally, we discuss a possible role of cofactor binding and its interaction with the G65 and G1 purines in structural stabilization of the A-1(2'-OH) in-line attack conformation. On the basis of the identified dominant protonation state of the reaction precursor, we propose a hypothesis of the self-cleavage mechanism in which A-1(2'-OH) is activated as a nucleophile by the G1(pro-R(p)) nonbridging oxygen of the scissile phosphate, whereas the ammonium group of GlcN6P acts as the general acid protonating the G1(O5') leaving group.

The gem-disubstituent effect—a computational study that exposes the relevance of existing theoretical models

The DFT B3LYP/6-31G (d, p) and HF/6-31G calculation results for the cyclization reactions of six different di-substituted 4-bromobutylamines (Brown’s system) reveal that the driving force for accelerations in the rate for intramolecular cyclization compared to the linear unsubstituted compound is due to strain effects (difference in strain energy in the ground state and in the transition state) and not to proximity in the ground state stemming from the distribution of rotational conformations ‘near attack conformation’ (time of residence at a reactive distance). Further, the results indicate that the reaction rate (log k rel) for the intramolecular cyclization process is strongly correlated with the distance between the two reacting centers ( r) and the attack angle ( α). The latter result is the first to correlate strain and reaction rate with geometrical parameters ( r and α).

Threshold Occupancy and Specific Cation Binding Modes in the Hammerhead Ribozyme Active Site are Required for Active Conformation

The relationship between formation of active in-line attack conformations and monovalent (Na +) and divalent (Mg 2+) metal ion binding in hammerhead ribozyme (HHR) has been explored with molecular dynamics simulations. To stabilize repulsions between negatively charged groups, different requirements of the threshold occupancy of metal ions were observed in the reactant and activated precursor states both in the presence and in the absence of a Mg 2+ in the active site. Specific bridging coordination patterns of the ions are correlated with the formation of active in-line attack conformations and can be accommodated in both cases. Furthermore, simulation results suggest that the HHR folds to form an electronegative recruiting pocket that attracts high local concentrations of positive charge. The present simulations help to reconcile experiments that probe the metal ion sensitivity of HHR catalysis and support the supposition that Mg 2+, in addition to stabilizing active conformations, plays a specific chemical role in catalysis.

Role of Mg2+ in Hammerhead Ribozyme Catalysis from Molecular Simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.

Dramatic Acceleration of the Menschutkin Reaction and Distortion of Halide Leaving-Group Order

Compared to structurally related linear trialkylamines, a simple macrocyclic amine with an anion-binding cavity exhibits very large rate enhancements (>10(5)) for stoichiometric N-alkylation with primary alkyl, allyl, and benzyl halides in the weakly polar solvent CDCl3. There is also a major distortion of the halide leaving-group order. For example, with benzyl halides the relative leaving-group order with a control amine is Cl (1) < Br (71) < I (160), whereas the leaving-group order with the macrocyclic amine is I (0.4) < Cl (1) < Br (8.5). Reaction with the macrocyclic amine is inhibited by the addition of DMSO, which is unusual because the Menschutkin reaction is normally enhanced by the presence of a polar aprotic solvent. Competitive inhibition studies indicate that the reaction proceeds through a prereaction complex. Effective molarities for the subsequent unimolecular N-alkylation step with 4-t-butylbenzyl halides are 4-t-BuBnCl (62,000 M) > 4-t-BuBnBr (2200 M) > 4-t-BuBnI (35 M); thus, the free energy of activation is selectively decreased for organohalides having smaller and more charge dense leaving groups. Likely reasons for this selective enhancement effect are: (a) increased transition-state stabilization due to hydrogen bonding in the macrocyclic pocket and (b) reduced entropic penalty in the transition state due to an increased fraction of prereaction complexes that are oriented in a near attack conformation. The study suggests that it should be possible to develop highly reactive macrocyclic amines that selectively sense or scavenge carcinogenic haloalkanes from the atmosphere.

Phosphorylation Reaction in cAPK Protein Kinase-Free Energy Quantum Mechanical/Molecular Mechanics Simulations

We present results of a theoretical analysis of the phosphorylation reaction in cAMP-dependent protein kinase using a combined quantum mechanical and molecular mechanics (QM/MM) approach. Detailed analysis of the reaction pathway is provided using a novel QM/MM implementation of the nudged elastic band method, finite temperature fluctuations of the protein environment are taken into account using free energy calculations, and an analysis of hydrogen bond interactions is performed on the basis of calculated frequency shifts. The late transfer of the substrate proton to the conserved aspartate (D166), the activation free energy of 15 kcal/mol, and the slight exothermic (-3 kcal/mol) character of the reaction are all consistent with the experimental data. The near attack conformation of D166 in the reactant state is maintained by interactions with threonine-201, asparagine-177, and most notably by a conserved water molecule serving as a strong structural link between the primary metal ion and the D166. The secondary Mg ion acts as a Lewis acid, attacking the beta-gamma bridging oxygen of ATP. This interaction, along with a strong hydrogen bond between the D166 and the substrate, contributes to the stabilization of the transition state. Lys-168 maintains a hydrogen bond to a transferring phosphoryl group throughout a reaction process. This interaction increases in the product state and contributes to its stabilization.

Strain and Near Attack Conformers in Enzymic Thiamin Catalysis: X-ray Crystallographic Snapshots of Bacterial Transketolase in Covalent Complex with Donor Ketoses Xylulose 5-phosphate and Fructose 6-phosphate, and in Noncovalent Complex with Acceptor Aldose Ribose 5-phosphate

Transketolase is a prominent thiamin diphosphate-dependent enzyme in sugar metabolism that catalyzes the reversible transfer of a 2-carbon dihydroxyethyl fragment between a donor ketose and an acceptor aldose. The X-ray structures of transketolase from E. coli in a covalent complex with donor ketoses d-xylulose 5-phosphate (X5P) and d-fructose 6-phosphate (F6P) at 1.47 A and 1.65 A resolution reveal significant strain in the tetrahedral cofactor-sugar adducts with a 25-30 degrees out-of-plane distortion of the C2-Calpha bond connecting the substrates' carbonyl with the C2 of the cofactor's thiazolium part. Both intermediates adopt very similar extended conformations in the active site with a perpendicular orientation of the scissile C2-C3 sugar bond relative to the thiazolium ring. The sugar-derived hydroxyl groups of the intermediates form conserved hydrogen bonds with one Asp side chain, with a cluster of His residues and with the N4' of the aminopyrimidine ring of the cofactor. The phosphate moiety is held in place by electrostatic and hydrogen-bonding interactions with Arg, His, and Ser side chains. With the exception of the thiazolium part of the cofactor, no structural changes are observable during intermediate formation indicating that the active site is poised for catalysis. DFT calculations on both X5P-thiamin and X5P-thiazolium models demonstrate that an out-of-plane distortion of the C2-Calpha bond is energetically more favorable than a coplanar bond. The X-ray structure with the acceptor aldose d-ribose 5-phosphate (R5P) noncovalently bound in the active site suggests that the sugar is present in multiple forms: in a strained ring-closed beta-d-furanose form in C2-exo conformation as well as in an extended acyclic aldehyde form, with the reactive C1 aldo function held close to Calpha of the presumably planar carbanion/enamine intermediate. The latter form of R5P may be viewed as a near attack conformation. The R5P binding site overlaps with those of the leaving group moieties of the covalent donor-cofactor adducts, demonstrating that R5P directly competes with the donor-derived products glyceraldehyde 3-phosphate and erythrose 4-phosphate, which are substrates of the reverse reaction, for the same docking site at the active site and reaction with the DHEThDP enamine.

Stereoselective iodocyclisation of 3-acylamino-2-methylene alkanoates: a computational insight

In order to explain the high stereocontrol occurring in the iodocyclisation of 3-acylamino-2-methylenealkanoates, either the conformational space of the starting products or the cyclisation reaction potential energy surface (PES) was explored at DFT level of theory, the polarised continuum formalism (PCM) for chloroform being used in order to consider the solvent effect. The observed stereoselection was ascribed to both the near attack conformations (NACs) distribution and the energy differences between the two possible and competitive cyclisation pathways leading to cis- and trans-diastereomers.

Low-frequency normal modes in horse liver alcohol dehydrogenase and motions of residues involved in the enzymatic reaction

Normal mode analysis using the elastic network model has provided characteristics and directions of the low-frequency large domain motions of horse liver alcohol dehydrogenase. Three normal modes (mode 1, mode 7, and mode 8) were identified as representative domain motions that may promote the onset of Near Attack Conformers or facilitate the product to be released. The pattern of the atomic displacement for some key residues (such as Val292 and Val203) revealed in this study is in line with experimental structural and kinetic studies and theoretical simulations.