Family 18 chitinases catalyze the hydrolysis of β-1,4-glycosidic bonds in chitin. The mechanism has been proposed to involve the formation of an oxazolinium ion intermediate via an unusual substrate-assisted mechanism, in which the substrate itself acts as an intramolecular nucleophile (instead of an enzyme residue). Here, we have modeled the first step of the chitin hydrolysis catalyzed by Serratia marcescens chitinase B for the first time using a combined quantum mechanics/molecular mechanics approach. The calculated reaction barriers based on multiple snapshots are 15.8-19.8 kcal mol(-1) [B3LYP/6-31+G(d)//AM1-CHARMM22], in good agreement with the activation free energy of 16.1 kcal mol(-1) derived from experiment. The enzyme significantly stabilizes the oxazolinium intermediate. Two stable conformations ((4)C(1)-chair and B(3,O)-boat) of the oxazolinium ion intermediate in subsite -1 were unexpectedly observed. The transition state structure has significant oxacarbenium ion-like character. The glycosyl residue in subsite -1 was found to follow a complex conformational pathway during the reaction ((1,4)B → [(4)H(5)/(4)E](++) → (4)C(1) ↔ B(3,O)), indicating complex conformational behavior in glycoside hydrolases that utilize a substrate-assisted catalytic mechanism. The D142N mutant is found to follow the same wild-type-like mechanism: the calculated barriers for reaction in this mutant (16.0-21.1 kcal mol(-1)) are higher than in the wild type, in agreement with the experiment. Asp142 is found to be important in transition state and intermediate stabilization.

Serratia marcescens Chitinase B (ChiB), belonging to the glycosidase family 18 (GH18), catalyzes the hydrolysis of β-1,4-glycosidic bond, with retention of configuration, via an unusual substrate-assisted mechanism, in which the substrate itself acts as an intramolecular nucleophile. Here, both elementary steps (glycosylation and deglycosylation) of the ChiB-catalyzed reaction are investigated by means of combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations at the SCC-DFTB/CHARMM22 level of theory. We examine the influence of the Asp142 protonation state on the reaction and the role that this residue performs in the reaction. Our simulations show that reaction with a neutral Asp142 is preferred and demonstrate that this residue provides electrostatic stabilization of the oxazolinium ion intermediate formed in the reaction. Insight into the conformational itinerary ((1,4)B↔(4)H5↔(4)C1) adopted by the substrate (bound in subsite -1) along the preferred reaction pathway is also provided by the simulations. The relative energies of the stationary points found along the reaction pathway calculated with SCC-DFTB and B3LYP were compared. The results suggest that SCC-DFTB is an accurate method for estimating the relative barriers for both steps of the reaction; however, it was found to overestimate the relative energy of an intermediate formed in the reaction when compared with the higher level of theory. Glycosylation is suggested to be a rate-determining step in the reaction with calculated overall reaction free-energy barrier of 20.5 kcal/mol, in a reasonable agreement with the 16.1 kcal/mol barrier derived from the experiment. The role of Tyr214 in catalysis was also investigated with the results, indicating that the residue plays a critical role in the deglycosylation step of the reaction. Simulations of the enzyme-product complex were also performed with an unbinding event suggested to have been observed, affording potential new mechanistic insight into the release of the product of ChiB.

No abstract

No abstract