Structural Insights from Molecular Dynamics Simulations of Tryptophan 7-Halogenase and Tryptophan 5-Halogenase

Identifying and Engineering Flavin Dependent Halogenases for Selective Biocatalysis

ConspectusOrganohalogen compounds are extensively used as building blocks, intermediates, pharmaceuticals, and agrochemicals due to their unique chemical and biological properties. Installing halogen substituents, however, frequently requires functionalized starting materials and multistep functional group interconversion. Several classes of halogenases evolved in nature to enable halogenation of a different classes of substrates; for example, site-selective halogenation of electron rich aromatic compounds is catalyzed by flavin-dependent halogenases (FDHs). Mechanistic studies have shown that these enzymes use FADH2 to reduce O2 to water with concomitant oxidation of X- to HOX (X = Cl, Br, I). This species travels through a tunnel within the enzyme to access the FDH active site. Here, it is believed to interact with an active site lysine proximal to bound substrate, enabling electrophilic halogenation with selectivity imparted via molecular recognition, rather than directing groups or strong electronic activation.The unique selectivity of FDHs led to several early biocatalysis efforts, preparative halogenation was rare, and the hallmark catalyst-controlled selectivity of FDHs did not translate to non-native substrates. FDH engineering was limited to site-directed mutagenesis, which resulted in modest changes in site-selectivity or substrate preference. To address these limitations, we optimized expression conditions for the FDH RebH and its cognate flavin reductase (FRed), RebF. We then showed that RebH could be used for preparative halogenation of non-native substrates with catalyst-controlled selectivity. We reported the first examples in which the stability, substrate scope, and site selectivity of a FDH were improved to synthetically useful levels via directed evolution. X-ray crystal structures of evolved FDHs and reversion mutations showed that random mutations throughout the RebH structure were critical to achieving high levels of activity and selectivity on diverse aromatic substrates, and these data were used in combination with molecular dynamics simulations to develop predictive model for FDH selectivity. Finally, we used family wide genome mining to identify a diverse set of FDHs with novel substrate scope and complementary regioselectivity on large, three-dimensionally complex compounds.The diversity of our evolved and mined FDHs allowed us to pursue synthetic applications beyond simple aromatic halogenation. For example, we established that FDHs catalyze enantioselective reactions involving desymmetrization, atroposelective halogenation, and halocyclization. These results highlight the ability of FDH active sites to tolerate different substrate topologies. This utility was further expanded by our recent studies on the single component FDH/FRed, AetF. While we were initially drawn to AetF because it does not require a separate FRed, we found that it halogenates substrates that are not halogenated efficiently or at all by other FDHs and provides high enantioselectivity for reactions that could only be achieved using RebH variants after extensive mutagenesis. Perhaps most notably, AetF catalyzes site-selective aromatic iodination and enantioselective iodoetherification. Together, these studies highlight the origins of FDH engineering, the utility and limitations of the enzymes developed to date, and the promise of FDHs for an ever-expanding range of biocatalytic halogenation reactions.

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.

Asymmetric catalysis by flavin-dependent halogenases

In nature, flavin-dependent halogenases (FDHs) catalyze site-selective chlorination and bromination of aromatic natural products. This ability has led to extensive efforts to engineer FDHs for selective chlorination, bromination, and iodination of electron rich aromatic compounds. On the other hand, FDHs are unique among halogenases and haloperoxidases that exhibit catalyst-controlled site selectivity in that no examples of enantioselective FDH catalysis in natural product biosynthesis have been characterized. Over the past several years, our group has established that FDHs can catalyze enantioselective reactions involving desymmetrization, atroposelective halogenation, and halocyclization. Achieving high activity and selectivity for these reactions has required extensive mutagenesis and mitigation of problems resulting from hypohalous acid generated during FDH catalysis. The single-component flavin reductase/FDH AetF is unique among the wild type enzyme we have studied in that it provides high activity and selectivity toward several asymmetric transformations. These results highlight the ability of FDH active sites to tolerate different substrate topologies and suggest that they could be useful for a broad range of oxidative halogenations.

The Single-Component Flavin Reductase/Flavin-Dependent Halogenase AetF is a Versatile Catalyst for Selective Bromination and Iodination of Arenes and Olefins

Flavin-dependent halogenases (FDHs) natively catalyze selective halogenation of electron rich aromatic and enolate groups. Nearly all FDHs reported to date require a separate flavin reductase to supply them with FADH2 , which complicates biocatalysis applications. In this study, we establish that the single component flavin reductase/flavin dependent halogenase AetF catalyzes halogenation of a diverse set of substrates using a commercially available glucose dehydrogenase to drive its halogenase activity. High site selectivity, activity on relatively unactivated substrates, and high enantioselectivity for atroposelective bromination and bromolactonization was demonstrated. Site-selective iodination and enantioselective cycloiodoetherification was also possible using AetF. The substrate and reaction scope of AetF suggest that it has the potential to greatly improve the utility of biocatalytic halogenation.

Mechanism of Action of Flavin-Dependent Halogenases

To rationally engineer the substrate scope and selectivity of flavin-dependent halogenases (FDHs), it is essential to first understand the reaction mechanism and substrate interactions in the active site. FDHs have long been known to achieve regioselectivity through an electrophilic aromatic substitution at C7 of the natural substrate Trp, but the precise role of a key active-site Lys residue remains ambiguous. Formation of hypochlorous acid (HOCl) at the cofactor-binding site is achieved by the direct reaction of molecular oxygen and a single chloride ion with reduced FAD and flavin hydroxide, respectively. HOCl is then guided 10 Å into the halogenation active site. Lys79, located in this site, has been proposed to direct HOCl toward Trp C7 through hydrogen bonding or a direct reaction with HOCl to form an -NH₂Cl⁺ intermediate. Here, we present the most likely mechanism for halogenation based on molecular dynamics (MD) simulations and active-site density functional theory "cluster" models of FDH PrnA in complex with its native substrate l-tryptophan, hypochlorous acid, and the FAD cofactor. MD simulations with different protonation states for key active-site residues suggest that Lys79 directs HOCl through hydrogen bonding, which is confirmed by calculations of the reaction profiles for both proposed mechanisms.

Analysis of Laboratory-Evolved Flavin-Dependent Halogenases Affords a Computational Model for Predicting Halogenase Site Selectivity

Flavin-dependent halogenases (FDHs) catalyze selective halogenation of electron-rich aromatic compounds without the need for harsh oxidants required by conventional oxidative halogenation reactions. Predictive models for halogenase site selectivity could greatly improve their utility for chemical synthesis. Toward this end, we analyzed the structures and selectivity of three halogenase variants evolved to halogenate tryptamine with orthogonal selectivity. Crystal structures and reversion mutations revealed key residues involved in altering halogenase selectivity. Density functional theory calculations and molecular dynamics simulations are both consistent with hypohalous acid as the active halogenating species in FDH catalysis. This model was used to accurately predict the site selectivity of halogenase variants toward different synthetic substrates, providing a valuable tool for implementing halogenases in biocatalysis efforts.

Directed Evolution of Flavin-Dependent Halogenases for Site- and Atroposelective Halogenation of 3-Aryl-4(3H)-Quinazolinones via Kinetic or Dynamic Kinetic Resolution

In this study, we engineer a variant of the flavin-dependent halogenase RebH that catalyzes site- and atroposelective halogenation of 3-aryl-4(3H)-quinazolinones via kinetic or dynamic kinetic resolution. The required directed evolution uses a combination of random and site-saturation mutagenesis, substrate walking using two probe substrates, and a two-tiered screening approach involving the analysis of variant conversion and then enantioselectivity of improved variants. The resulting variant, 3-T, provides >99:1 e.r. for the (M)-atropisomer of the major brominated product, 25-fold improved conversion, and 91-fold improved site selectivity relative to the parent enzyme on the probe substrate used in the final rounds of evolution. This high activity and selectivity translate well to several additional substrates with varied steric and electronic properties. Computational modeling and docking simulations are used to rationalize the effects of key mutations on substrate binding. Given the range of substrates that have been used for atroposelective synthesis via electrophilic halogenation in the literature, these results suggest that flavin-dependent halogenases (FDHs) could find many additional applications for atroposelective catalysis. More broadly, this study highlights how RebH can be engineered to accept structurally diverse substrates that enable its use for enantioselective catalysis.

Engineered RebH Halogenase Variants Demonstrating a Specificity Switch from Tryptophan towards Novel Indole Compounds

Activating industrially important aromatic hydrocarbons by installing halogen atoms is extremely important in organic synthesis and often improves the pharmacological properties of drug molecules. To this end, tryptophan halogenase enzymes are potentially valuable tools for regioselective halogenation of arenes, including various industrially important indole derivatives and similar scaffolds. Although endogenous enzymes show reasonable substrate scope towards indole compounds, their efficacy can often be improved by engineering. Using a structure-guided semi-rational mutagenesis approach, we have developed two RebH variants with expanded biocatalytic repertoires that can efficiently halogenate several novel indole substrates and produce important pharmaceutical intermediates. Interestingly, the engineered enzymes are completely inactive towards their natural substrate tryptophan in spite of their high tolerance to various functional groups in the indole ring. Computational modelling and molecular dynamics simulations provide mechanistic insights into the role of gatekeeper residues in the substrate binding site and the dramatic switch in substrate specificity when these are mutated.

Comparative molecular dynamics simulations identify a salt-sensitive loop responsible for the halotolerant activity of GH5 cellulases

Halotolerant glycoside hydrolases (GH) have broad application potentials in biorefinery industries. Elucidating the structure-activity relationship underlying the halotolerant catalysis is essential to design superior biocatalysts. Here, we performed molecular dynamics simulations to investigate the structural dynamics of two GH5 cellulases, namely the halotolerant Cel5R and non-halotolerant TfCel5A. Through characterizing the physical properties at different salt concentrations, the results revealed that the overall structures of Cel5R and TfCel5A were marginally affected by the increase in salt concentrations. However, a salt-sensitive loop was identified from both Cel5R and TfCel5A based on its significantly increased flexibility at high salt concentrations. Importantly, compared to TfCel5A the salt-sensitive loop of Cel5R engaged more sodium ions and water molecules around the active site of the enzyme. Besides, the unique residue motif of the salt-sensitive loop in Cel5R formed more intramolecular hydrogen bonds, stabilizing the active architecture of Cel5R at high salt concentrations. Collectively, the structural and dynamic differences may contribute to the various catalytic halotolerance of Cel5R and TfCel5A. These findings provide mechanistic insight into the halotolerant catalysis and will guide the ration design of GH5 cellulases with improved catalytic properties.Communicated by Ramaswamy H. Samy.

Genomic Determinants Encode the Reactivity and Regioselectivity of Flavin-Dependent Halogenases in Bacterial Genomes and Metagenomes

Halogenases create diverse natural products by utilizing halide ions and are of great interest in the synthesis of potential pharmaceuticals and agrochemicals. An increasing number of halogenases discovered in microorganisms are annotated as flavin-dependent halogenases (FDHs), but their chemical reactivities are markedly different and the genomic contents associated with such functional distinction have not been revealed yet. Even though the reactivity and regioselectivity of FDHs are essential in the halogenation activity, these FDHs are annotated inaccurately in the protein sequence repositories without characterizing their functional activities. We carried out a comprehensive sequence analysis and biochemical characterization of FDHs. Using a probabilistic model that we built in this study, FDHs were discovered from 2,787 bacterial genomes and 17 sediment metagenomes. We analyzed the essential genomic determinants that are responsible for substrate binding and subsequent reactions: four flavin adenine dinucleotide-binding, one halide-binding, and four tryptophan-binding sites. Compared with previous studies, our study utilizes large-scale genomic information to propose a comprehensive set of sequence motifs that are related to the active sites and regioselectivity. We reveal that the genomic patterns and phylogenetic locations of the FDHs determine the enzymatic reactivities, which was experimentally validated in terms of the substrate scope and regioselectivity. A large portion of publicly available FDHs needs to be reevaluated to designate their correct functions. Our genomic models establish comprehensive links among genotypic information, reactivity, and regioselectivity of FDHs, thereby laying an important foundation for future discovery and classification of novel FDHs.

IMPORTANCE Halogenases are playing an important role as tailoring enzymes in biosynthetic pathways. Flavin-dependent tryptophan halogenases (Trp-FDHs) are among the enzymes that have broad substrate scope and high selectivity. From bacterial genomes and metagenomes, we found highly diverse halogenase sequences by using a well-trained profile hidden Markov model built from the experimentally validated halogenases. The characterization of genotype, steady-state activity, substrate scope, and regioselectivity has established comprehensive links between the information encoded in the genomic sequence and reactivity of FDHs reported here. By constructing models for accurate and detailed sequence markers, our work should guide future discovery and classification of novel FDHs.

AoiQ Catalyzes Geminal Dichlorination of 1,3-Diketone Natural Products

Enzymes that can perform halogenation of aliphatic carbons are of significant interest to the synthetic and biocatalysis communities. Here we describe the characterization of AoiQ, a single-component flavin-dependent halogenase (FDH) that catalyzes gem-dichlorination of 1,3-diketone substrates in the biosynthesis of dichlorodiaporthin. AoiQ represents the first biochemically reconstituted FDH that can halogenate an enolizable sp3-hybridized carbon atom.

Two Novel, Flavin-Dependent Halogenases from the Bacterial Consortia of Botryococcus braunii Catalyze Mono- and Dibromination

Halogen substituents often lead to a profound effect on the biological activity of organic compounds. Flavin-dependent halogenases offer the possibility of regioselective halogenation at non-activated carbon atoms, while employing only halide salts and molecular oxygen. However, low enzyme activity, instability, and narrow substrate scope compromise the use of enzymatic halogenation as an economical and environmentally friendly process. To overcome these drawbacks, it is of tremendous interest to identify novel halogenases with high enzymatic activity and novel substrate scopes. Previously, Neubauer et al. developed a new hidden Markov model (pHMM) based on the PFAM tryptophan halogenase model, and identified 254 complete and partial putative flavin-dependent halogenase genes in eleven metagenomic data sets. In the present study, the pHMM was used to screen the bacterial associates of the Botryococcus braunii consortia (PRJEB21978), leading to the identification of several putative, flavin-dependent halogenase genes. Two of these new halogenase genes were found in one gene cluster of the Botryococcus braunii symbiont Sphingomonas sp. In vitro activity tests revealed that both heterologously expressed enzymes are active flavin-dependent halogenases able to halogenate indole and indole derivatives, as well as phenol derivatives, while preferring bromination over chlorination. Interestingly, SpH1 catalyses only monohalogenation, while SpH2 can catalyse both mono- and dihalogenation for some substrates.

Structure of apo flavin-dependent halogenase Xcc4156 hints at a reason for cofactor-soaking difficulties

Flavin-dependent halogenases regioselectively introduce halide substituents into electron-rich substrates under mild reaction conditions. For the enzyme Xcc4156 from Xanthomonas campestris, the structure of a complex with the cofactor flavin adenine dinucleotide (FAD) and a bromide ion would be of particular interest as this enzyme exclusively brominates model substrates in vitro. Apo Xcc4156 crystals diffracted to 1.6 Å resolution. The structure revealed an open substrate-binding site lacking the loop regions that close off the active site and contribute to substrate binding in tryptophan halogenases. Therefore, Xcc4156 might accept larger substrates, possibly even peptides. Soaking of apo Xcc4156 crystals with FAD led to crumbling of the intergrown crystals. Around half of the crystals soaked with FAD did not diffract, while in the others there was no electron density for FAD. The FAD-binding loop, which changes its conformation between the apo and the FAD-bound form in related enzymes, is involved in a crystal contact in the apo Xcc4156 crystals. The conformational change that is predicted to occur upon FAD binding would disrupt this crystal contact, providing a likely explanation for the destruction of the apo crystals in the presence of FAD. Soaking with only bromide did not result in bromide bound to the catalytic halide-binding site. Simultaneous soaking with FAD and bromide damaged the crystals more severely than soaking with only FAD. Together, these latter two observations suggest that FAD and bromide bind to Xcc4156 with positive cooperativity. Thus, apo Xcc4156 crystals provide functional insight into FAD and bromide binding, even though neither the cofactor nor the halide is visible in the structure.

Conformational flexibility influences structure-function relationships in nucleic acid N-methyl demethylases

N-Methylation of DNA/RNA bases can be regulatory or damaging and is linked to diseases including cancer and genetic disorders. Bacterial AlkB and human FTO are DNA/RNA demethylases belonging to the Fe(ii) and 2-oxoglutarate oxygenase superfamily. Modelling studies reveal conformational dynamics influence structure-function relationships of AlkB and FTO, e.g. why 1-methyladenine is a better substrate for AlkB than 6-methyladenine. Simulations show that the flexibility of the double stranded DNA substrate in AlkB influences correlated motions, including between the core jelly-roll fold and an active site loop involved in substrate binding. The FTO N- and C-terminal domains move in respect to one another in a manner likely important for substrate binding. Substitutions, including clinically observed ones, influencing catalysis contribute to the network of correlated motions in AlkB and FTO. Overall, the calculations highlight the importance of the overall protein environment and its flexibility to the geometry of the reactant complexes.

Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms

The combined quantum mechanics/molecular mechanics (QM/MM) methods have become a valuable tool in computational biochemistry and received versatile applications for studying the reaction mechanisms of enzymes. The approach combines the calculations of the electronic structure of the active site by QM, with modeling of the protein environment using MM force field, which allows the long-range electrostatics and steric effects on the enzyme reactivity to be accounted for. In this review, we review some key theoretical and computational aspects of the method and we also present some applications to particular enzymatic reactions such as tryptophan-7-halogenase, cyclooxygenase-1, and the epidermal growth factor receptor.