Spectroscopically Guided Simulations Reveal Distinct Strategies for Positioning Substrates to Achieve Selectivity in Nonheme Fe(II)/α-Ketoglutarate-Dependent Halogenases

Discovery and substrate specificity engineering of nucleotide halogenases

C2'-halogenation has been recognized as an essential modification to enhance the drug-like properties of nucleotide analogs. The direct C2'-halogenation of the nucleotide 2'-deoxyadenosine-5'-monophosphate (dAMP) has recently been achieved using the Fe(II)/α-ketoglutarate-dependent nucleotide halogenase AdaV. However, the limited substrate scope of this enzyme hampers its broader applications. In this study, we report two halogenases capable of halogenating 2'-deoxyguanosine monophosphate (dGMP), thereby expanding the family of nucleotide halogenases. Computational studies reveal that nucleotide specificity is regulated by the binding pose of the phosphate group. Based on these findings, we successfully engineered the substrate specificity of these halogenases by mutating second-sphere residues. This work expands the toolbox of nucleotide halogenases and provides insights into the regulation mechanism of nucleotide specificity.

Proton-triggered chemoselective halogenation of aliphatic C-H bonds with nonheme FeIV-oxo complexes

Halogenation of aliphatic C-H bonds is a chemical transformation performed in nature by mononuclear nonheme iron dependent halogenases. The mechanism involves the formation of an iron(IV)-oxo-chloride species that abstracts the hydrogen atom from the reactive C-H bond to form a carbon-centered radical that selectively reacts with the bound chloride ligand, a process commonly referred to as halide rebound. The factors that determine the halide rebound, as opposed to the reaction with the incipient hydroxide ligand, are not clearly understood and examples of well-defined iron(IV)-oxo-halide compounds competent in C-H halogenation are scarce. In this work we have studied the reactivity of three well-defined iron(IV)-oxo complexes containing variants of the tetradentate 1-(2-pyridylmethyl)-1,4,7-triazacyclononane ligand (Pytacn). Interestingly, these compounds exhibit a change in their chemoselectivity towards the functionalization of C-H bonds under certain conditions: their reaction towards C-H bonds in the presence of a halide anionleads to exclusive oxygenation, while the addition of a superacid results in halogenation. Almost quantitative halogenation of ethylbenzene is observed when using the two systems with more sterically congested ligands and even the chlorination of strong C-H bonds such as those of cyclohexane is performed when a methyl group is present in the sixth position of the pyridine ring of the ligand. Mechanistic studies suggest that both reactions, oxygenation and halogenation, proceed through a common rate determining hydrogen atom transfer step and the presence of the acid dictates the fate of the resulting alkyl radical towards preferential halogenation over oxygenation.

Controllable multi-halogenation of a non-native substrate by SyrB2 iron halogenase

Geminal, multi-halogenated functional groups are widespread in natural products and pharmaceuticals, yet no synthetic methodologies exist that enable selective multi-halogenation of unactivated C-H bonds. Biocatalysts are powerful tools for late-stage C-H functionalization, as they operate with high degrees of regio-, chemo-, and stereoselectivity. 2-oxoglutarate (2OG)-dependent non-heme iron halogenases chlorinate and brominate aliphatic C-H bonds offering a solution for achieving these challenging transformations. Here, we describe the ability of a non-heme iron halogenase, SyrB2, to controllably halogenate non-native substrate alpha-aminobutyric acid (Aba) to yield mono-chlorinated, di-chlorinated, and tri-chlorinated products. These chemoselective outcomes are achieved by controlling the loading of 2OG cofactor and SyrB2 biocatalyst. By using a ferredoxin-based biological reductant for electron transfer to the catalytic center of SyrB2, we demonstrate order-of-magnitude enhancement in the yield of tri-chlorinated product that were previously inaccessible using any single halogenase enzyme. We also apply these strategies to broaden SyrB2's reactivity scope to include multi-bromination and demonstrate chemoenzymatic conversion of the ethyl side chain in Aba to an ethylyne functional group. We show how steric hindrance induced by the successive addition of halogen atoms on Aba's C4 carbon dictates the degree of multi-halogenation by hampering C3-C4 bond rotation within SyrB2's catalytic pocket. Overall, our work showcases the synthetic potential of iron halogenases to facilitate multi-C-H functionalization chemistry.

How Do Differences in Electronic Structure Affect the Use of Vanadium Intermediates as Mimics in Nonheme Iron Hydroxylases?

We study active-site models of nonheme iron hydroxylases and their vanadium-based mimics using density functional theory to determine if vanadyl is a faithful structural mimic. We identify crucial structural and energetic differences between ferryl and vanadyl isomers owing to the differences in their ground electronic states, i.e., high spin (HS) for Fe and low spin (LS) for V. For the succinate cofactor bound to the ferryl intermediate, we predict facile interconversion between monodentate and bidentate coordination isomers for ferryl species but difficult rearrangement for vanadyl mimics. We study isomerization of the oxo intermediate between axial and equatorial positions and find the ferryl potential energy surface to be characterized by a large barrier of ca. 10 kcal/mol that is completely absent for the vanadyl mimic. This analysis reveals even starker contrasts between Fe and V in hydroxylases than those observed for this metal substitution in nonheme halogenases. Analysis of the relative bond strengths of coordinating carboxylate ligands for Fe and V reveals that all of the ligands show stronger binding to V than Fe owing to the LS ground state of V in contrast to the HS ground state of Fe, highlighting the limitations of vanadyl mimics of native nonheme iron hydroxylases.

Substrate Positioning Dynamics Involves a Non-Electrostatic Component to Mediate Catalysis

Substrate positioning dynamics (SPD) orients the substrate in the active site, thereby influencing catalytic efficiency. However, it remains unknown whether SPD effects originate primarily from electrostatic perturbation inside the enzyme or can independently mediate catalysis with a significant non-electrostatic component. In this work, we investigated how the non-electrostatic component of SPD affects transition state (TS) stabilization. Using high-throughput enzyme modeling, we selected Kemp eliminase variants with similar electrostatics inside the enzyme but significantly different SPD. The kinetic parameters of these mutants were experimentally characterized. We observed a valley-shaped, two-segment linear correlation between the TS stabilization free energy (converted from kinetic parameters) and substrate positioning index (a metric to quantify SPD). The energy varies by approximately 2 kcal/mol. Favorable SPD was observed for the distal mutant R154W, increasing the proportion of reactive conformations and leading to the lowest activation free energy. These results indicate the substantial contribution of the non-electrostatic component of SPD to enzyme catalytic efficiency.

Why Nonheme Iron Halogenases Do Not Fluorinate C-H Bonds: A Computational Investigation

Selective halogenation is necessary for a range of fine chemical applications, including the development of therapeutic drugs. While synthetic processes to achieve C-H halogenation require harsh conditions, enzymes such as nonheme iron halogenases carry out some types of C-H halogenation, i.e., chlorination or bromination, with ease, while others, i.e., fluorination, have never been observed in natural or engineered nonheme iron enzymes. Using density functional theory and correlated wave function theory, we investigate the differences in structural and energetic preferences of the smaller fluoride and the larger chloride or bromide intermediates throughout the catalytic cycle. Although we find that the energetics of rate-limiting hydrogen atom transfer are not strongly impacted by fluoride substitution, the higher barriers observed during the radical rebound reaction for fluoride relative to chloride and bromide contribute to the difficulty of C-H fluorination. We also investigate the possibility of isomerization playing a role in differences in reaction selectivity, and our calculations reveal crucial differences in terms of isomer energetics of the key ferryl intermediate between fluoride and chloride/bromide intermediates. While formation of monodentate isomers believed to be involved in selective catalysis is shown for chloride and bromide intermediates, we find that formation of the fluoride monodentate intermediate is not possible in our calculations, which lack additional stabilizing interactions with the greater protein environment. Furthermore, the shorter Fe-F bonds are found to increase isomerization reaction barriers, suggesting that incorporation of residues that form a halogen bond with F and elongate Fe-F bonds could make selective C-H fluorination possible in nonheme iron halogenases. Our work highlights the differences between the fluoride and chloride/bromide intermediates and suggests potential steps toward engineering nonheme iron halogenases to enable selective C-H fluorination.

Mutexa: A Computational Ecosystem for Intelligent Protein Engineering

Protein engineering holds immense promise in shaping the future of biomedicine and biotechnology. This Review focuses on our ongoing development of Mutexa, a computational ecosystem designed to enable "intelligent protein engineering". In this vision, researchers will seamlessly acquire sequences of protein variants with desired functions as biocatalysts, therapeutic peptides, and diagnostic proteins through a finely-tuned computational machine, akin to Amazon Alexa's role as a versatile virtual assistant. The technical foundation of Mutexa has been established through the development of a database that combines and relates enzyme structures and their respective functions (e.g., IntEnzyDB), workflow software packages that enable high-throughput protein modeling (e.g., EnzyHTP and LassoHTP), and scoring functions that map the sequence-structure-function relationship of proteins (e.g., EnzyKR and DeepLasso). We will showcase the applications of these tools in benchmarking the convergence conditions of enzyme functional descriptors across mutants, investigating protein electrostatics and cavity distributions in SAM-dependent methyltransferases, and understanding the role of nonelectrostatic dynamic effects in enzyme catalysis. Finally, we will conclude by addressing the future steps and fundamental challenges in our endeavor to develop new Mutexa applications that assist the identification of beneficial mutants in protein engineering.

Bioinspired complexes confined in well-defined capsules: getting closer to metalloenzyme functionalities

Reproducing the key features offered by metalloprotein binding cavities is an attractive approach to overcome the main bottlenecks of current open artificial models (in terms of stability, efficiency and selectivity). In this context, this featured article brings together selected examples of recent developments in the field of confined bioinspired complexes with an emphasis on the emerging hemicryptophane caged ligands. In particular, we focused on (1) the strategies allowing the insulation and protection of complexes sharing similarities with metalloprotein active sites, (2) the confinement-induced improvement of catalytic efficiencies and selectivities and (3) very recent efforts that have been made toward the development of bioinspired complexes equipped with weakly binding artificial cavities.

Data-driven enzyme engineering to identify function-enhancing enzymes

Identifying function-enhancing enzyme variants is a 'holy grail' challenge in protein science because it will allow researchers to expand the biocatalytic toolbox for late-stage functionalization of drug-like molecules, environmental degradation of plastics and other pollutants, and medical treatment of food allergies. Data-driven strategies, including statistical modeling, machine learning, and deep learning, have largely advanced the understanding of the sequence-structure-function relationships for enzymes. They have also enhanced the capability of predicting and designing new enzymes and enzyme variants for catalyzing the transformation of new-to-nature reactions. Here, we reviewed the recent progresses of data-driven models that were applied in identifying efficiency-enhancing mutants for catalytic reactions. We also discussed existing challenges and obstacles faced by the community. Although the review is by no means comprehensive, we hope that the discussion can inform the readers about the state-of-the-art in data-driven enzyme engineering, inspiring more joint experimental-computational efforts to develop and apply data-driven modeling to innovate biocatalysts for synthetic and pharmaceutical applications.

Convergence in determining enzyme functional descriptors across Kemp eliminase variants

Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSD_{active_site} and SASA_ratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EF_C-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.

What Drives Radical Halogenation versus Hydroxylation in Mononuclear Nonheme Iron Complexes? A Combined Experimental and Computational Study

Nonheme iron halogenases are unique enzymes in nature that selectively activate an aliphatic C–H bond of a substrate to convert it into C–X (X = Cl/Br, but not F/I). It is proposed that they generate an Fe^III(OH)(X) intermediate in their catalytic cycle. The analogous Fe^III(OH) intermediate in nonheme iron hydroxylases transfers OH^• to give alcohol product, whereas the halogenases transfer X^• to the carbon radical substrate. There remains significant debate regarding what factors control their remarkable selectivity of the halogenases. The reactivity of the complexes Fe^III(BNPA^Ph2O)(OH)(X) (X = Cl, Br) with a secondary carbon radical (R^•) is described. It is found that X^• transfer occurs with a secondary carbon radical, as opposed to OH^• transfer with tertiary radicals. Comprehensive computational studies involving density functional theory were carried out to examine the possible origins of this selectivity. The calculations reproduce the experimental findings, which indicate that halogen transfer is not observed for the tertiary radicals because of a nonproductive equilibrium that results from the endergonic nature of these reactions, despite a potentially lower reaction barrier for the halogenation pathway. In contrast, halogen transfer is favored for secondary carbon radicals, for which the halogenated product complex is thermodynamically more stable than the reactant complex. These results are rationalized by considering the relative strengths of the C–X bonds that are formed for tertiary versus secondary carbon centers. The computational analysis also shows that the reaction barrier for halogen transfer is significantly dependent on secondary coordination sphere effects, including steric and H-bonding interactions.

Structural and Functional Insights into a Nonheme Iron- and α-Ketoglutarate-Dependent Halogenase That Catalyzes Chlorination of Nucleotide Substrates

Nonheme iron- and α-ketoglutarate (αKG)-dependent halogenases (NHFeHals), which catalyze the regio- and stereoselective halogenation of the unactivated C(sp³)-H bonds, exhibit tremendous potential in the challenging asymmetric halogenation. AdeV from Actinomadura sp. ATCC 39365 is the first identified carrier protein-free NHFeHal that catalyzes the chlorination of nucleotide 2'-deoxyadenosine-5'-monophosphate (2'-dAMP) to afford 2'-chloro-2'-deoxyadenosine-5'-monophosphate. Here, we determined the complex crystal structures of AdeV/Fe^II/Cl and AdeV/Fe^II/Cl/αKG at resolutions of 1.76 and 1.74 Å, respectively. AdeV possesses a typical β-sandwich topology with H194, H252, αKG, chloride, and one water molecule coordinating Fe^II in the active site. Molecular docking, mutagenesis, and biochemical analyses reveal that the hydrophobic interactions and hydrogen bond network between the substrate-binding pocket and the adenine, deoxyribose, and phosphate moieties of 2'-dAMP are essential for substrate recognition. Residues H111, R177, and H192 might play important roles in the second-sphere interactions that control reaction partitioning. This study provides valuable insights into the catalytic selectivity of AdeV and will facilitate the rational engineering of AdeV and other NHFeHals for synthesis of halogenated nucleotides. IMPORTANCE Halogenated nucleotides are a group of important antibiotics and are clinically used as antiviral and anticancer drugs. AdeV is the first carrier protein-independent nonheme iron- and α-ketoglutarate (αKG)-dependent halogenase (NHFeHal) that can selectively halogenate nucleotides and exhibits restricted substrate specificity toward several 2'-dAMP analogues. Here, we determined the complex crystal structures of AdeV/Fe^II/Cl and AdeV/Fe^II/Cl/αKG. Molecular docking, mutagenesis, and biochemical analyses provide important insights into the catalytic selectivity of AdeV. This study will facilitate the rational engineering of AdeV and other carrier protein-independent NHFeHals for synthesis of halogenated nucleotides.

Molecular Dynamics-Derived Descriptor Informs the Impact of Mutation on the Catalytic Turnover Number in Lactonase Across Substrates

Molecular dynamics simulations have been extensively employed to reveal the roles of protein dynamics in mediating enzyme catalysis. However, simulation-derived predictive descriptors that inform the impacts of mutations on catalytic turnover numbers remain largely unexplored. In this work, we report the identification of molecular modeling-derived descriptors to predict mutation effect on the turnover number of lactonase SsoPox with both native and non-native substrates. The study consists of 10 enzyme-substrate complexes resulting from a combination of five enzyme variants with two substrates. For each complex, we derived 15 descriptors from molecular dynamics simulations and applied principal component analysis to rank the predictive capability of the descriptors. A top-ranked descriptor was identified, which is the solvent-accessible surface area (SASA) ratio of the substrate to the active site pocket. A uniform volcano-shaped plot was observed in the distribution of experimental activation free energy against the SASA ratio. To achieve efficient lactonase hydrolysis, a non-native substrate-bound enzyme variant needs to involve a similar range of the SASA ratio to the native substrate-bound wild-type enzyme. The descriptor reflects how well the enzyme active site pocket accommodates a substrate for reaction, which has the potential of guiding optimization of enzyme reaction turnover for non-native chemical transformations.

Mechanistic analysis of carbon-carbon bond formation by deoxypodophyllotoxin synthase

The completion of the tetracyclic core of etoposide, classified by the World Health Organization as an essential medicine, by the Fe/2OG oxygenase deoxypodophyllotoxin synthase follows a hybrid radical-polar pathway not previously seen in other members of this enzyme class. The implication of a substrate-based benzylic carbocation in this mechanism will inform ongoing efforts to create analogs of this important drug with improved or emergent properties and represents a new route for resolution of the initial substrate radical that is common to members of the class. This study adds to our understanding on a growing number of biochemical transformations in which carbocation intermediates are likely to be crucial.

Deoxypodophyllotoxin contains a core of four fused rings (A to D) with three consecutive chiral centers, the last being created by the attachment of a peripheral trimethoxyphenyl ring (E) to ring C. Previous studies have suggested that the iron(II)- and 2-oxoglutarate–dependent (Fe/2OG) oxygenase, deoxypodophyllotoxin synthase (DPS), catalyzes the oxidative coupling of ring B and ring E to form ring C and complete the tetracyclic core. Despite recent efforts to deploy DPS in the preparation of deoxypodophyllotoxin analogs, the mechanism underlying the regio- and stereoselectivity of this cyclization event has not been elucidated. Herein, we report 1) two structures of DPS in complex with 2OG and (±)-yatein, 2) in vitro analysis of enzymatic reactivity with substrate analogs, and 3) model reactions addressing DPS’s catalytic mechanism. The results disfavor a prior proposal of on-pathway benzylic hydroxylation. Rather, the DPS-catalyzed cyclization likely proceeds by hydrogen atom abstraction from C7', oxidation of the benzylic radical to a carbocation, Friedel–Crafts-like ring closure, and rearomatization of ring B by C6 deprotonation. This mechanism adds to the known pathways for transformation of the carbon-centered radical in Fe/2OG enzymes and suggests what types of substrate modification are likely tolerable in DPS-catalyzed production of deoxypodophyllotoxin analogs.

Reaction pathway engineering converts a radical hydroxylase into a halogenase

Fe^II/α-ketoglutarate (Fe^II/αKG)-dependent enzymes offer a promising biocatalytic platform for halogenation chemistry owing to their ability to functionalize unactivated C-H bonds. However, relatively few radical halogenases have been identified to date, limiting their synthetic utility. Here, we report a strategy to expand the palette of enzymatic halogenation by engineering a reaction pathway rather than substrate selectivity. This approach could allow us to tap the broader class of Fe^II/αKG-dependent hydroxylases as catalysts by their conversion to halogenases. Toward this goal, we discovered active halogenases from a DNA shuffle library generated from a halogenase-hydroxylase pair using a high-throughput in vivo fluorescent screen coupled to an alkyne-producing biosynthetic pathway. Insights from sequencing halogenation-active variants along with the crystal structure of the hydroxylase enabled engineering of a hydroxylase to perform halogenation with comparable activity and higher selectivity than the wild-type halogenase, showcasing the potential of harnessing hydroxylases for biocatalytic halogenation.