Mechanism of glycoside hydrolysis: A comparative QM/MM molecular dynamics analysis for wild type and Y69F mutant retaining xylanases.

Computational simulations have been performed using hybrid quantum-mechanical/molecular-mechanical potentials to investigate the catalytic mechanism of the retaining endo-beta-1, 4-xylanase (BCX) from B. circulans. Two-dimensional potential-of-mean-force calculations based upon molecular dynamics with the AM1/OPLS method for wild-type BCX with a p-nitrophenyl xylobioside substrate in water clearly indicates a stepwise mechanism for glycosylation: the rate-determining step is nucleophilic substitution by Glu78 to form the covalently bonded enzyme-substrate intermediate without protonation of the leaving group by Glu172. The geometrical configuration of the transition state for the enzymic reaction is essentially the same as found for a gas-phase model involving only the substrate and a propionate/propionic acid pair to represent the catalytic glutamate/glutamic acid groups. In addition to stabilizing the (2,5)B boat conformation of the proximal xylose in the non-covalent reactant complex of the substrate with BCX, Tyr69 lowers the free-energy barrier for glycosylation by 42 kJ mol(-1) relative to that calculated for the Y69F mutant, which lacks the oxygen atom O(Y). B3LYP/6-31+G* energy corrections reduce the absolute height of the barrier to reaction. In the oxacarbenium ion-like transition state O(Y) approaches closer to the endocyclic oxygen O(ring) of the sugar ring but donates its hydrogen bond not to O(ring) but rather to the nucleophilic oxygen of Glu78. Comparison of the average atomic charge distributions for the wild-type and mutant indicates that charge separation along the bond between the anomeric carbon and O(ring) is matched in the former by a complementary separation of charge along the O(Y)-H(Y) bond, corresponding to a pair of roughly antiparallel bond dipoles, which is not present in the latter.

Computational mutagenesis reveals the role of active-site tyrosine in stabilising a boat conformation for the substrate: QM/MM molecular dynamics studies of wild-type and mutant xylanases.

Molecular dynamics simulations have been performed for non-covalent complexes of phenyl beta-xylobioside with the retaining endo-beta-1,4-xylanase from B. circulans (BCX) and its Tyr69Phe mutant using a hybrid QM/MM methodology. A trajectory initiated for the wild-type enzyme-substrate complex with the proximal xylose ring bound at the -1 subsite (adjacent to the scissile glycosidic bond) in the (4)C(1) chair conformation shows spontaneous transformation to the (2,5)B boat conformation, and potential of mean force calculations indicate that the boat is approximately 30 kJ mol(-1) lower in free energy than the chair. Analogous simulations for the mutant lacking one oxygen atom confirm the key role of Tyr69 in stabilizing the boat in preference to the (4)C(1) chair conformation, with a relative free energy difference of about 20 kJ mol(-1), by donating a hydrogen bond to the endocyclic oxygen of the proximal xylose ring. QM/MM MD simulations for phenyl beta-xyloside in water, with and without a propionate/propionic acid pair to mimic the catalytic glutamate/glutamic acid pair of the enzyme, show the (4)C(1) chair to be stable, although a hydrogen bond between the OH group at C2 of xylose and the propionate moiety seems to provide some stabilization for the (2,5)B conformation.