Deciphering the evolutionary origin of the stereoselectivity of short-chain dehydrogenases in the oxidation of the monoterpenol 1-borneol

Enzyme engineering has produced numerous methods to optimize enzymes for biotechnological processes; however, less is known about how natural evolution creates new functionalities. We investigate the evolutionary emergence of enantioselectivity in plant borneol dehydrogenases (BDHs), which feature hydrophobic active-sites and are enantioselective towards dibornane-type monoterpenols. Ancestral sequence reconstruction provided a trajectory from the oldest unselective BDH ancestor N30 (E=12) toward the youngest selective ancestor N32, involving 19 mutations: 18 mutations are peripheral, one (I111L) occurs in the active-site. The mutation L111I in the hydrophobic pocket increased the selectivity of N30, while the back-mutation I111L decreased the selectivity of N32. Additional peripheral mutations (V136L/G169A/V183I) were required for high selectivity. Crystal structures suggested that protein dynamics, rather than structural changes shape these catalytic properties. Molecular simulations with funnel-metadynamics revealed a correlation between the active-sites solvent-accessible surface area (SASA) and selectivity. This potential evolutionary pathway shapes enantioselectivity, and guides future enzyme engineering campaigns.