A computational method for the systematic screening of reaction barriers in enzymes: Searching for Bacillus circulans xylanase mutants with greater activity towards a synthetic substrate

TL;DR

This paper introduces a semi-empirical computational method to systematically screen enzyme mutants for improved activity, demonstrated on Bacillus circulans xylanase, identifying promising mutants for experimental validation.

Contribution

The paper presents a novel semi-empirical approach for estimating enzyme mutation effects on reaction barriers, enabling efficient mutant screening without extensive quantum calculations.

Findings

PM6 estimates align well with experimental data.

Identified mutants with significantly increased reaction barriers.

Developed a comprehensive mutant library for Bacillus circulans xylanase.

Abstract

We present a semi-empirical (PM6-based) computational method for systematically estimating the effect of all possible single mutants, within a certain radius of the active site, on the barrier height of an enzymatic reaction. The intent of this method is not a quantitative prediction of the barrier heights, but rather to identify promising mutants for further computational or experimental study. The method is applied to identify promising single and double mutants of Bacillus circulans xylanase (BCX) with increased hydrolytic activity for the artificial substrate ortho-nitrophenyl \beta-xylobioside (ONPX$_2$). The estimated reaction barrier for wild-type (WT) BCX is 18.5 kcal/mol, which is in good agreement with the experimental activation free energy value of 17.0 kcal/mol extracted from the observed k$_\text{cat}$ using transition state theory (Joshi et al., Biochemistry 2001, 40, 10115). The PM6 reaction profiles for eight single point mutations are recomputed using FMO-MP2/PCM/6-31G(d) single points. PM6 predicts an increase in barrier height for all eight mutants while FMO predicts an increase for six of the eight mutants. Both methods predict that the largest change in barrier occurs for N35F, where PM6 and FMO predict a 9.0 and 15.8 kcal/mol increase, respectively. We thus conclude that PM6 is sufficiently accurate to identify promising mutants for further study. We prepared a set of all theoretically possible (342) single mutants in which every amino acid of the active site (except for the catalytically active residues E78 and E172) was mutated to every other amino acid. Based on results from the single mutants we construct a set of 111 double mutants consisting of all possible pairs of single mutants with the lowest barrier for a particular position and compute their reaction profile. None of the mutants have, to our knowledge, been prepared experimentally[...].