QM/MM MD Simulations on the Origin of Donor/Acceptor Selectivity of Family GH51 α-L-Arabinofuranosidase-Catalyzed Hydrolysis and Transglycosylation Reactions
Wijitra Jitonnom, Tanchanok Wanjai, Jitrayut Jitonnom
Abstract
α- l -Arabinofuranosidase (EC 3.2.1.55) plays a crucial role in the degradation and processing of arabinoxylan and arabinoxylo-oligosaccharides (AXOS), making it a promising target for industrial biomass conversion and biotechnological applications. Despite its ability to cleave α(1→2), α(1→3), and α(1→5) glycosidic bonds, the origin of its preference for specific AXOS linkages remains unclear. Using quantum mechanics (QM) calculations and quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations, we investigated the catalytic mechanism and selectivity of hydrolysis and transglycosylation in a GH51 α- l -arabinofuranosidase from Thermobacillus xylanilyticus . We assessed various synthetic substrates and mono- and disubstituted AXOS ( p NP-Ara f, A 3, A 2 X, A 3 X, XA 2 X, XA 3 X and XA 2+3 X), predicting their binding modes and substrate specificity. Glu176 (acid/base) and Glu298 (nucleophile) facilitate AXOS hydrolysis via a two-step, double-displacement mechanism, with the second step (deglycosylation) being the rate-limiting step. Based on combined QM/MM MD and QM(ONIOM) studies, His240 is most likely protonated to maintain the positioning of the Glu176 residue within the Asp297-His240-Glu176 catalytic triad, thereby enabling efficient catalysis via an optimal electrostatic environment. Consistent with the experiment, the enzyme efficiently hydrolyzes p NP-Ara f and A 3, with free-energy barriers of 14.8–15.4 kcal mol –1 (compared to 17.9–29.3 kcal mol –1 for AXOS); it selectively cleaves Ara f -α(1→3) over Ara f -α(1→2) linkages of donor substrates and favors xylobiose as an acceptor over aryl α- d -xylopyranoside in transglycosylation due to lower reaction barriers and effective transition-state stabilization. Our study provides a theoretical foundation for the engineering of this GH51 enzyme for efficient hydrolysis and synthesis of valuable xylan-based products.