Conformational diversity and enantioconvergence in potato epoxide hydrolase 1

Evidence of a Distinctive Enantioselective Binding Mode for the Photoinduced Radical Cyclization of α-Chloroamides in Ene-Reductases

We demonstrate here through molecular simulations and mutational studies the origin of the enantioselectivity in the photoinduced radical cyclization of α-chloroacetamides catalyzed by ene-reductases, in particular the Gluconobacter oxidans ene-reductase and the Old Yellow Enzyme 1, which show opposite enantioselectivity. Our results reveal that neither the π-facial selectivity model nor a protein-induced selective stabilization of the transition states is able to explain the enantioselectivity of the radical cyclization in the studied flavoenzymes. We propose a new enantioinduction scenario according to which enantioselectivity is indeed controlled by transition-state stability; however, the relative stability of the prochiral transition states is not determined by direct interaction with the protein but is rather dependent on an inherent degree of freedom within the substrate itself. This intrinsic degree of freedom, distinct from the traditional π-facial exposure mode, can be controlled by the substrate conformational selection upon binding to the enzyme.

Engineered P450 Atom-Transfer Radical Cyclases are Bifunctional Biocatalysts: Reaction Mechanism and Origin of Enantioselectivity

New-to-nature radical biocatalysis has recently emerged as a powerful strategy to tame fleeting open-shell intermediates for stereoselective transformations. In 2021, we introduced a novel metalloredox biocatalysis strategy that leverages the innate redox properties of the heme cofactor of P450 enzymes, furnishing new-to-nature atom-transfer radical cyclases (ATRCases) with excellent activity and stereoselectivity. Herein, we report a combined computational and experimental study to shed light on the mechanism and origins of enantioselectivity for this system. Molecular dynamics and quantum mechanics/molecular mechanics (QM/MM) calculations revealed an unexpected role of the key beneficial mutation I263Q. The glutamine residue serves as an essential hydrogen bond donor that engages with the carbonyl moiety of the substrate to promote bromine atom abstraction and enhance the enantioselectivity of radical cyclization. Therefore, the evolved ATRCase is a bifunctional biocatalyst, wherein the heme cofactor enables atom-transfer radical biocatalysis, while the hydrogen bond donor residue further enhances the activity and enantioselectivity. Unlike many enzymatic stereocontrol rationales based on a rigid substrate binding model, our computations demonstrate a high degree of rotational flexibility of the allyl moiety in an enzyme-substrate complex and succeeding intermediates. Therefore, the enantioselectivity is controlled by the radical cyclization transition states rather than the substrate orientation in ground-state complexes in the preceding steps. During radical cyclization, anchoring effects of the Q263 residue and steric interactions with the heme cofactor concurrently control the π-facial selectivity, allowing for highly enantioselective C-C bond formation. Our computational findings are corroborated by experiments with ATRCase mutants generated from site-directed mutagenesis.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)₈-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)₈-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

Structure-function relationship between soluble epoxide hydrolases structure and their tunnel network

Enzymes with buried active sites maintain their catalytic function via a single tunnel or tunnel network. In this study we analyzed the functionality of soluble epoxide hydrolases (sEHs) tunnel network, by comparing the overall enzyme structure with the tunnel's shape and size. sEHs were divided into three groups based on their structure and the tunnel usage. The obtained results were compared with known substrate preferences of the studied enzymes, as well as reported in our other work evolutionary analyses data. The tunnel network architecture corresponded well with the evolutionary lineage of the source organism and large differences between enzymes were observed from long fragments insertions. This strategy can be used during protein re-engineering process for large changes introduction, whereas tunnel modification can be applied for fine-tuning of enzyme.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Ground-State Destabilization by Active-Site Hydrophobicity Controls the Selectivity of a Cofactor-Free Decarboxylase

Bacterial arylmalonate decarboxylase (AMDase) and evolved variants have become a valuable tool with which to access both enantiomers of a broad range of chiral arylaliphatic acids with high optical purity. Yet, the molecular principles responsible for the substrate scope, activity, and selectivity of this enzyme are only poorly understood to date, greatly hampering the predictability and design of improved enzyme variants for specific applications. In this work, empirical valence bond and metadynamics simulations were performed on wild-type AMDase and variants thereof to obtain a better understanding of the underlying molecular processes determining reaction outcome. Our results clearly reproduce the experimentally observed substrate scope and support a mechanism driven by ground-state destabilization of the carboxylate group being cleaved by the enzyme. In addition, our results indicate that, in the case of the nonconverted or poorly converted substrates studied in this work, increased solvent exposure of the active site upon binding of these substrates can disturb the vulnerable network of interactions responsible for facilitating the AMDase-catalyzed cleavage of CO2. Finally, our results indicate a switch from preferential cleavage of the pro-(R) to the pro-(S) carboxylate group in the CLG-IPL variant of AMDase for all substrates studied. This appears to be due to the emergence of a new hydrophobic pocket generated by the insertion of the six amino acid substitutions, into which the pro-(S) carboxylate binds. Our results allow insight into the tight interaction network determining AMDase selectivity, which in turn provides guidance for the identification of target residues for future enzyme engineering.

Manipulating regioselectivity of an epoxide hydrolase for single enzymatic synthesis of (R)-1,2-diols from racemic epoxides

Both the activity and regioselectivity of Phaseolus vulgaris epoxide hydrolase were remarkably improved via reshaping two substrate tunnels based on rational design.

Exploring Solanum tuberosum Epoxide Hydrolase Internal Architecture by Water Molecules Tracking

Several different approaches are used to describe the role of protein compartments and residues in catalysis and to identify key residues suitable for the modification of the activity or selectivity of the desired enzyme. In our research, we applied a combination of molecular dynamics simulations and a water tracking approach to describe the water accessible volume of Solanum tuberosum epoxide hydrolase. Using water as a molecular probe, we were able to identify small cavities linked with the active site: (i) one made up of conserved amino acids and indispensable for the proper positioning of catalytic water and (ii) two others in which modification can potentially contribute to enzyme selectivity and activity. Additionally, we identified regions suitable for de novo tunnel design that could also modify the catalytic properties of the enzyme. The identified hot-spots extend the list of the previously targeted residues used for modification of the regioselectivity of the enzyme. Finally, we have provided an example of a simple and elegant process for the detailed description of the network of cavities and tunnels, which can be used in the planning of enzyme modifications and can be easily adapted to the study of any other protein.

Epoxide hydrolysis as a model system for understanding flux through a branched reaction scheme

The epoxide hydrolase StEH1 catalyzes the hydrolysis of trans-methylstyrene oxide to 1-phenyl-propane-1,2-diol. The (S,S)-epoxide is exclusively transformed into the (1R,2S)-diol, while hydrolysis of the (R,R)-epoxide results in a mixture of product enantiomers. In order to understand the differences in the stereoconfigurations of the products, the reactions were studied kinetically during both the pre-steady-state and steady-state phases. A number of closely related StEH1 variants were analyzed in parallel, and the results were rationalized by structure-activity analysis using the available crystal structures of all tested enzyme variants. Finally, empirical valence-bond simulations were performed in order to provide additional insight into the observed kinetic behaviour and ratios of the diol product enantiomers. These combined data allow us to present a model for the flux through the catalyzed reactions. With the (R,R)-epoxide, ring opening may occur at either C atom and with similar energy barriers for hydrolysis, resulting in a mixture of diol enantiomer products. However, with the (S,S)-epoxide, although either epoxide C atom may react to form the covalent enzyme intermediate, only the pro-(R,S) alkylenzyme is amenable to subsequent hydrolysis. Previously contradictory observations from kinetics experiments as well as product ratios can therefore now be explained for this biocatalytically relevant enzyme.

Exploring the origins of selectivity in soluble epoxide hydrolase from Bacillus megaterium

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. In this work, the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative are computationally investigated using density functional theory (DFT).

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. These enzymes present different enantio- and regioselectivity depending upon either the substrate structure or the substitution pattern of the epoxide ring. In this study, we computationally investigate the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative using density functional theory (DFT) and an active site cluster model consisting of 195 and 197 atoms, respectively. Full reaction mechanisms for epoxide ring opening were evaluated considering the attack at both oxirane carbons and considering two possible orientations of the substrate at the BmEH active site. Our results indicate that for both SO and p-NSO substrates the BmEH enantio- and regioselectivity is opposite to the inherent (R)-BmEH selectivity, the attack at the benzylic position (C1) of the (S)-enantiomer being the most favoured chemical outcome.

Active-Site Flexibility and Substrate Specificity in a Bacterial Virulence Factor: Crystallographic Snapshots of an Epoxide Hydrolase

Pseudomonas aeruginosa secretes an epoxide hydrolase with catalytic activity that triggers degradation of the cystic fibrosis transmembrane conductance regulator (CFTR) and perturbs other host defense networks. Targets of this CFTR inhibitory factor (Cif) are largely unknown, but include an epoxy-fatty acid. In this class of signaling molecules, chirality can be an important determinant of physiological output and potency. Here we explore the active-site chemistry of this two-step α/β-hydrolase and its implications for an emerging class of virulence enzymes. In combination with hydrolysis data, crystal structures of 15 trapped hydroxyalkyl-enzyme intermediates reveal the stereochemical basis of Cif's substrate specificity, as well as its regioisomeric and enantiomeric preferences. The structures also reveal distinct sets of conformational changes that enable the active site to expand dramatically in two directions, accommodating a surprising array of potential physiological epoxide targets. These new substrates may contribute to Cif's diverse effects in vivo, and thus to the success of P. aeruginosa and other pathogens during infection.

CADEE: Computer-Aided Directed Evolution of Enzymes

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

Laboratory-Evolved Enzymes Provide Snapshots of the Development of Enantioconvergence in Enzyme-Catalyzed Epoxide Hydrolysis

Engineered enzyme variants of potato epoxide hydrolase (StEH1) display varying degrees of enrichment of (2R)‐3‐phenylpropane‐1,2‐diol from racemic benzyloxirane. Curiously, the observed increase in the enantiomeric excess of the (R)‐diol is not only a consequence of changes in enantioselectivity for the preferred epoxide enantiomer, but also to changes in the regioselectivity of the epoxide ring opening of (S)‐benzyloxirane. In order to probe the structural origin of these differences in substrate selectivity and catalytic regiopreference, we solved the crystal structures for the evolved StEH1 variants. We used these structures as a starting point for molecular docking studies of the epoxide enantiomers into the respective active sites. Interestingly, despite the simplicity of our docking analysis, the apparent preferred binding modes appear to rationalize the experimentally determined regioselectivities. The analysis also identifies an active site residue (F33) as a potentially important interaction partner, a role that could explain the high conservation of this residue during evolution. Overall, our experimental, structural, and computational studies provide snapshots into the evolution of enantioconvergence in StEH1‐catalyzed epoxide hydrolysis.