Loop dynamics redesign is an important strategy to manipulate protein function. Cellobiose 2-epimerase (CE) and other members of its superfamily are widely used for diverse industrial applications. The structural feature of the loops connecting barrel helices contributes greatly to the differences in their functional characteristics. Inspired by the in-silico mutation with molecular dynamics (MD) simulation analysis, we propose a strategy for identifying disulfide bond mutation candidates based on the prediction of protein flexibility and residue-residue interaction. The most beneficial mutant with the newly introduced disulfide bond would simultaneously improve both its thermostability and its reaction propensity to the targeting isomerization product. The ratio of the isomerization/epimerization catalytic rate was improved from 4:103 to 9:22. MD simulation and binding free energy calculations were applied to provide insights into molecular recognition upon mutations. The comparative analysis of enzyme/substrate binding modes indicates that the altered catalytic reaction pathway is due to less efficient binding of the native product. The key residue responsible for the observed phenotype was identified by energy decomposition and was further confirmed by the mutation experiment. The rational design of the key loop region might be a promising strategy to alter the catalytic behavior of all (α/α)6-barrel-like proteins.
Read full abstract