Potato epoxide hydrolase 1 (StEH1) is a biocatalytically important enzyme that exhibits rich enantio- and regioselectivity in the hydrolysis of chiral epoxide substrates. In particular, StEH1 has been demonstrated to enantioconvergently hydrolyze racemic mixes of styrene oxide (SO) to yield (R)-1-phenylethanediol. This work combines computational, crystallographic and biochemical analyses to understand both the origins of the enantioconvergent behavior of the wild-type enzyme, as well as shifts in activities and substrate binding preferences in an engineered StEH1 variant, R-C1B1, which contains four active site substitutions (W106L, L109Y, V141K and I155V). Our calculations are able to reproduce both the enantio- and regioselectivities of StEH1, and demonstrate a clear link between different substrate binding modes and the corresponding selectivity, with the preferred binding modes being shifted between the wild-type enzyme and the R-C1B1 variant. Additionally, we demonstrate that the observed changes in selectivity and the corresponding enantioconvergent behavior are due to a combination of steric and electrostatic effects that modulate both the accessibility of the different carbon atoms to the nucleophilic side chain of D105, as well as the interactions between the substrate and protein amino acid side chains and active site water molecules. Being able to computationally predict such subtle effects for different substrate enantiomers, as well as to understand their origin and how they are affected by mutations, is an important advance towards the computational design of improved biocatalysts for enantioselective synthesis.
Potato epoxide hydrolase 1 exhibits rich enantio- and regioselectivity in the hydrolysis of a broad range of substrates. The enzyme can be engineered to increase the yield of optically pure products as a result of changes in both enantio- and regioselectivity. It is thus highly attractive in biocatalysis, particularly for the generation of enantiopure fine chemicals and pharmaceuticals. The present work aims to establish the principles underlying the activity and selectivity of the enzyme through a combined computational, structural, and kinetic study using the substrate trans-stilbene oxide as a model system. Extensive empirical valence bond simulations have been performed on the wild-type enzyme together with several experimentally characterized mutants. We are able to computationally reproduce the differences between the activities of different stereoisomers of the substrate and the effects of mutations of several active-site residues. In addition, our results indicate the involvement of a previously neglected residue, H104, which is electrostatically linked to the general base H300. We find that this residue, which is highly conserved in epoxide hydrolases and related hydrolytic enzymes, needs to be in its protonated form in order to provide charge balance in an otherwise negatively charged active site. Our data show that unless the active-site charge balance is correctly treated in simulations, it is not possible to generate a physically meaningful model for the enzyme that can accurately reproduce activity and selectivity trends. We also expand our understanding of other catalytic residues, demonstrating in particular the role of a noncanonical residue, E35, as a 'backup base' in the absence of H300. Our results provide a detailed view of the main factors driving catalysis and regioselectivity in this enzyme and identify targets for subsequent enzyme design efforts.
