Original Research ARTICLE
Does Substrate Positioning Affect the Selectivity and Reactivity in the Hectochlorin Biosynthesis Halogenase?
- 1The Manchester Institute of Biotechnology and School of Chemical Engineering and Analytical Science, University of Manchester, Manchester, United Kingdom
- 2The Manchester Institute of Biotechnology and School of Chemistry, University of Manchester, Manchester, United Kingdom
- 3Institute of Biochemistry, Graz University of Technology, Graz, Austria
- 4Institute of Molecular Biosciences, Graz University, Graz, Austria
In this work we present the first computational study on the hectochlorin biosynthesis enzyme HctB, which is a unique three-domain halogenase that activates non-amino acid moieties tethered to an acyl-carrier, and as such may have biotechnological relevance beyond other halogenases. We use a combination of small cluster models and full enzyme structures calculated with quantum mechanics/molecular mechanics methods. Our work reveals that the reaction is initiated with a rate-determining hydrogen atom abstraction from substrate by an iron (IV)-oxo species, which creates an iron (III)-hydroxo intermediate. In a subsequent step the reaction can bifurcate to either halogenation or hydroxylation of substrate, but substrate binding and positioning drives the reaction to optimal substrate halogenation. Furthermore, several key residues in the protein have been identified for their involvement in charge-dipole interactions and induced electric field effects. In particular, two charged second coordination sphere amino acid residues (Glu223 and Arg245) appear to influence the charge density on the Cl ligand and push the mechanism toward halogenation. Our studies, therefore, conclude that nonheme iron halogenases have a chemical structure that induces an electric field on the active site that affects the halide and iron charge distributions and enable efficient halogenation. As such, HctB is intricately designed for a substrate halogenation and operates distinctly different from other nonheme iron halogenases.
Enzymatic C–Cl bond formation is a rare process in Nature, yet over the past few decades a range of haloperoxidases and halogenases have been discovered (Gribble, 2003; Vaillancourt et al., 2006; van Pée et al., 2006; Butler and Sandy, 2009; Wagner et al., 2009; Weichold et al., 2016; Agarwal et al., 2017; Schnepel and Sewald, 2017; Timmins and de Visser, 2018). Their catalytic mechanism, however, is still subject to controversies and understanding the fundamental details of these processes may have an impact on biotechnological advances as well as drug development. Three different classes of halogenation enzymes appear in Nature, namely the heme haloperoxidases, the vanadium-dependent nonheme haloperoxidases and the α-ketoglutarate dependent nonheme iron halogenases (Gribble, 2003; Vaillancourt et al., 2006; van Pée et al., 2006; Butler and Sandy, 2009; Wagner et al., 2009; Weichold et al., 2016; Agarwal et al., 2017; Schnepel and Sewald, 2017; Timmins and de Visser, 2018).
Heme haloperoxidases, (Sundaramoorthy et al., 1995; Wagenknecht and Woggon, 1997; Green et al., 2004; Kim et al., 2006) such as chloroperoxidase bind H2O2 on an iron(III)-heme center, which is then converted into an iron(IV)-oxo heme cation radical species called Compound I, (Meunier et al., 2004; Denisov et al., 2005; Shaik et al., 2005; Rittle and Green, 2010) with the help of a proton shuttle machinery. Compound I subsequently reacts with chloride to form OCl− products. The product drifts out of the active site and reacts with substrates through halogenation. A second class of haloperoxidases has a nonheme vanadium co-factor that also utilizes hydrogen peroxide and halide in a catalytic cycle to form OCl− products (Messerschmidt and Wever, 1996; Chen and van Pée, 2008). The vanadium haloperoxidases have been characterized in marine algae and are believed to have functions related to natural product synthesis associated with defense mechanisms (Martinez et al., 2001).
The final class of halogenases are the α-ketoglutarate (αKG) dependent nonheme iron halogenases, (Blasiak et al., 2006; Buongiorno and Straganz, 2013; Huang and Groves, 2017) which show structural and functional similarities to the corresponding nonheme iron—α-ketoglutarate-dependent dioxygenases (Solomon et al., 2000; Bugg, 2001; Ryle and Hausinger, 2002; Costas et al., 2004; Abu-Omar et al., 2005; Bruijnincx et al., 2008). These nonheme iron halogenases contain an iron(II) center in the resting state that is bound to the side chains of two histidine residues. During the catalytic cycle, halide and αKG bind to the iron(II) ion prior to molecular oxygen binding. It is believed dioxygen attacks αKG to form succinate and a high-valent iron(IV)-oxo species similarly to the αKG-dependent hydroxylases (Schofield and Zhang, 1999; Bollinger et al., 2005). The iron(IV)-oxo species of several nonheme iron halogenases have been spectroscopically trapped and characterized and shown to react with substrates with a rate-determining hydrogen atom abstraction barrier (Galonić Fujimori et al., 2007; Neidig et al., 2007; Matthews et al., 2009; Wong et al., 2013; Srnec et al., 2016; Srnec and Solomon, 2017). The subsequent pathway leading to the halogenated product; however, is controversial as a thermodynamically much more favorable hydroxylation pathway is prevented. How the enzyme manages to perform this unfavorable thermodynamic reaction pathway is under much debate. One proposal suggested links to substrate binding and orientation (Matthews et al., 2006). Thus, in the antibiotic biosynthesis protein SyrB2 the substrate L-Thr is linked to an acyl-carrier protein (SyrB1) via a phosphopantetheinyl (PPT) bridge. Replacement of L-Thr by L-norvaline changed the chemoselectivity of substrate halogenation to hydroxylation due to substrate positioning in the active site (Matthews et al., 2006). On the other hand, computational modeling of the Borowski group proposed a rotation step in the iron(IV)-oxo(halide) group, whereby after hydrogen atom abstraction, the positions of the hydroxo and halide ions switched position leading to easier halide rebound (Borowski et al., 2010). Clearly, the mechanism of aliphatic halogenation remains controversial. Because of this, synthetic (biomimetic) model complexes of these nonheme iron halogenases and haloperoxidases have been developed and studied with experimental (Podgoršek et al., 2009; Comba and Wunderlich, 2010; Liu and Groves, 2010; Chatterjee and Paine, 2016; Puri et al., 2016; Wang et al., 2016) and computational approaches (Noack and Siegbahn, 2007; de Visser and Latifi, 2009; Kulik et al., 2009; Pandian et al., 2009; Quesne and de Visser, 2012; Senn, 2014; Huang et al., 2016; Timmins et al., 2018). Thus, if the halogenation pathway is dependent on substrate positioning then biomimetic models that lack the protein and do not bind substrates in fixed orientations may not be able to react via substrate halogenation efficiently.
A recently discovered halogenase (HctB) from Lyngbya majuscula is involved in the biosynthesis of hectochlorin, whereby a fatty acyl substrate is dihalogenated on a nonheme iron center (Ramaswamy et al., 2007; Pratter et al., 2014a). In contrast to the nonheme iron halogenase SyrB2, HctB activates a non-amino acid group as a substrate. This would give HctB biotechnological applicability that is beyond that of halogenases like SyrB2. As such, HctB appears to have certain flexibility in the activation of an alkyl chain; therefore, we decided to investigate its mechanism using computational methods and compare the mechanism and active site features with those described for SyrB2 previously. The mechanism of O2 activation for this mononuclear nonheme iron halogenase was characterized with stopped-flow and spectroscopic (circular dichroism and magnetic circular dichroism) studies (Pratter et al., 2014b).
Functional characterization of the enzyme together with analysis of its primary structure reveals that it is a unique three-domain halogenase containing an acyl carrier protein (ACP), which binds the substrate covalently via a phosphopantetheinyl bridge (Figure 1). The other terminus of the ACP group is connected to the C-terminus of an acyl-Coenzyme-A subunit, while the N-terminus is linked to the halogenase domain. As such the protein has an intricate set-up and arrangement that enables dihalogenation of the tethered substrate. Interestingly, no evidence of substrate hydroxylation is available, but the activation of hexanoic acid by HctB apart from dihalogenation also gave products from oxygenation leading to vinyl-chloride and ketone products (Pratter et al., 2014a). HctB shows a certain degree of sequence similarity with SyrB2 (see Supporting Information for an overlay and sequence alignment), (Pratter et al., 2014a) and binds iron in a nonheme ligand configuration and utilizes α-ketoglutarate (αKG). The metal in HctB is linked to the protein via two histidine groups (His111 and His227) in a pentacoordinate environment. However, there are major differences in the overall structure of SyrB2 vs. HctB, which are currently understood and warrant a detailed computational study.
Figure 1. (A) Structure of the resting state of HctB as taken from Pratter et al. (2014b). The different domains of the secondary structure are identified and a close-up of the active site structure is given. (B) General reaction mechanism of the first halogenation step in HctB investigated here.
After αKG and halide binding, the metal binds molecular oxygen, and it has been hypothesized to react with αKG to form an iron(IV)-oxo species and succinate upon release of CO2. Unfortunately, the iron(IV)-oxo species in HctB has never been trapped and characterized and details of the halogenation mechanism are unknown. In the next stage of the catalytic cycle it is expected that a hydrogen atom abstraction occurs to form an iron(III)-hydroxo species, which is elusive as well. Technically, however, the halide should rebound to the substrate radical to form halogenated products, but thermodynamically it has been shown that hydroxo rebound is energetically favorable (Timmins and de Visser, 2015). How the enzyme avoids the low-energy hydroxylation pathway in favor of the higher-energy halogenation remains a matter of discussion. In order to understand the effect of substrate binding on the chemoselectivity of substrate halogenation vs. hydroxylation in HctB we employed a detailed molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) study. We located two substrate entrance channels and have calculated the substrate halogenation and hydroxylation pathways with a substrate in these positions. The work shows that dramatic differences in halogenation vs. hydroxylation product ratios should be expected based upon substrate positioning. Moreover, our calculations predict that the protein induces an electric field effect that withdraws electron density from the halide toward the metal during the reaction in order to make the halogenation process favorable.
Model Selection and Reactant Benchmarking
So far, no computational modeling studies have been reported on HctB and little is known how different (if at all) it is to the well-studied SyrB2. Obviously, both nonheme iron halogenases utilize a different substrate although both are tethered to a carrier protein. Conversely, SyrB2 catalyzes a monohalogenation of the substrate, whereas HctB performs a dihalogenation instead. A comparison of the amino acid sequences of HctB vs. SyrB2 shows major deviations which must result in differences in secondary structure. Therefore, their mechanisms may be quite different and we decided to do a computational study on the mechanism of halogenation vs. hydroxylation of HctB and how it compares with previous experimental and computational studies of analogous nonheme iron halogenases. In particular, we focused our computational study on the first halogenation step of the fatty acyl tethered substrate of a HctB model and took the structure reported by Pratter et al. (2014a). Our set-up of QM/MM models was reviewed thoroughly recently; therefore, we will summarize the main issues only briefly (Hernández-Ortega et al., 2015; Quesne et al., 2016a; Hofer and de Visser, 2018). The model was altered from an iron(II)-water α-ketoglutarate bound structure into an iron(IV)-oxo with ligated succinate. Prior to the full set-up of the complete QM/MM chemical system, however, we investigated possible substrate binding positions of the tethered hexanoyl-PPT moiety. In particular, we searched for alternative substrate entrance channels into the active site.
Figure 2 shows highlights of the two substrate entrance channels we identified, which are narrow channels that should fit the linear terminal chain of the substrate-carrier protein. The acyl-carrier group of the substrate protein then latches onto the surface of the protein and inserts the tethered hexanoyl-PPT group into the active site. Thus, substrate entrance channel II represents the model from Pratter et al. (2014a) with substrate position 1. In model 1 with the substrate entering through channel II, the substrate is located parallel to the iron(IV)-oxo group in a position very similar to the methylated DNA strand in the AlkB repair enzyme (Quesne et al., 2014). The substrate approaches the iron(IV)-oxo group in between the Val113 and Glu223 side chains in the corner where the halide group is also located. The entrance channel is located well below the iron(IV)-oxo group.
Figure 2. Substrate entrance channels I (in purple) and II (in cyan) into the active site of the halogenase domain of HctB. The yellow circle highlights the active site of the halogenase in the protein. The right-hand-side of the figure gives an overlay of the four substrate binding positions, where we only show the hexanoyl group and highlight the transferring hydrogen atom in white. Substrate binding position 1 is located in entrance channel II and substrate binding position 2, 3, and 4 are bound in entrance channel I.
During analysis of the structure, we identified another entrance channel and attempted to latch the acyl carrier protein onto this and manually inserted the tethered hexanoyl-PPT chain. In the analogous halogenase SyrB2 a similar substrate entrance channel is seen in the same position, (Matthews et al., 2006) but due to the lack of a substrate bound crystal structure of HctB we will consider all possible substrate orientations and entrance channels. Entrance channel I is located above the iron(IV)-oxo group and inserts the tethered hexanoyl group into a large open space (possibly filled with water molecules) and hence we created three starting orientations for the substrate, namely models 2, 3, and 4. In model 2 the terminus of the tethered hexanoyl group points down and hangs in between the Val113 and Arg245 residues, whereas in model 3 it is found in between Lys126 and Pro127. The final substrate position (4) through channel I is located on the side of the iron(IV)-oxo group and also brings the terminus of the tethered hexanoyl group into the active site nearby the Val113 and Glu223 residues. The substrate binding orientation through channel I model 4 brings it in close proximity to both the oxo and halide groups. These three substrate binding positions are distinctly different and guided by hydrogen bonding interactions of the thioester moiety of the hexanoyl-PPT moiety with amino acid residues aligning the channel walls. As such we do not expect easy interconversion between the three substrate binding models.
Subsequently, we investigated the reaction mechanism from each of these different substrate starting orientations and set up QM/MM models with the substrate bound in positions 1, 2, 3 and 4. The set-up of the QM/MM models follows previously reported and benchmarked methods (Quesne et al., 2014, 2016a; Hernández-Ortega et al., 2015; Hofer and de Visser, 2018). Firstly, we added hydrogen atoms using known pKa values and visually inspected polar residues for correct protonation states. Thereafter, we applied an iterative solvation procedure and followed this with an equilibration and heating run to a temperature of 298 K. During these set-up steps the protein backbone was fixed, but in the final molecular dynamics simulation all atoms were allowed to move. We initially ran MD simulations for a period of up to 200 ns; however, the run stabilized in all cases after about 5 ns. Clearly, HctB is a very rigid protein with substrate and active site in tight binding orientation with little flexibility as seen from the MD runs. Therefore, for all subsequent systems described here only a 10 ns MD simulation was taken, see Supporting Information Figure S1. From the MD simulations we selected a low energy snapshot after 5 ns as starting structures for the actual QM/MM calculations.
Finally, we bisected the full chemical structure of protein, substrate and water layer into a QM and MM region and included key residues in the QM region that form covalent or hydrogen bonding interactions with substrate and oxidant. Initial, exploratory calculations were done with a small QM region containing only the iron(IV)-oxo(chloro) group and the first-coordination sphere of ligands to the metal, namely the imidazole groups of His111 and His227, the acetate terminus of succinate (Succ) and the thiohexanoic acid arm of the substrate as our minimal QM region A, see Figure 3. This model contains 48 QM atoms and is overall charge neutral. To test the effect of the second coordination sphere we also calculated the full mechanism with a larger QM region that in addition to QM region A included the amino acid side chains of residues within 6 Å of the iron(IV)-oxo-chloro structure, i.e. QM region AB. For model 1, the large QM region AB contained the amino acid side chains of Ile108, Val113, Asp200, Asp202, Glu223, Val225, Met226, Arg245, and six water molecules and has a total of 160 atoms in the QM region.
Figure 3. Description of QM region A and AB for Model 1. Ball-and-stick model of QM region AB is shown on the right-hand-side.
Thereafter, QM/MM geometry optimizations of the iron(IV)-oxo(chloro) species (Re) were performed in the singlet, triplet and quintet spin states. Note that in the label of the structure we give the spin multiplicity in superscript before the label and the substrate binding position and the QM region (A or AB) in subscript after the label. The full set of results that were obtained is given in the Supporting Information, while we focus in the main text on the low-energy pathways only.
We started the work with extensive validation and benchmarking of the methods. Unfortunately, there are no experimental rate constants and spectroscopic data for the iron(IV)-oxo species. Previously, we calibrated thioanisole sulfoxidation free energies of activation of a biomimetic nonheme iron(IV)-oxo complex against experimental data and tested 50 different computational methods and techniques (Cantú Reinhard et al., 2016a). The best agreement with experiment was found for the PBE0/BS2 and B3LYP/BS2 methods with a solvent model included. In particular, free energies of activation were reproduced within 4 kcal mol−1 from experiment (Vardhaman et al., 2011, 2013; de Visser et al., 2014; Sainna et al., 2015; Barman et al., 2016a; Cantú Reinhard et al., 2016b). Furthermore, the methods reproduced experimental product distributions of bifurcation processes well (Ji et al., 2015; Kaczmarek et al., 2018). Finally, for the nonheme iron dioxygenase prolyl-4-hydroxylase six hydrogen atom abstraction barriers from substrate were investigated and the QM/MM predicted the correct regioselectivity and therefore the methods are expected to predict regio- and chemoselectivities well (Karamzadeh et al., 2010; Pratter et al., 2013; Timmins and de Visser, 2017; Timmins et al., 2017).
In agreement with the experimental studies on the iron(IV)-oxo species of the halogenase SyrB2, (Galonić Fujimori et al., 2007) we find the quintet spin state as the ground state for all chemical systems. Interestingly, the quintet spin state is below the triplet spin state by a considerable margin of well over 13 kcal mol−1; hence the triplet and singlet spin states of HctB will not play a major role during the reaction and are high in energy. Therefore, we will focus on the quintet spin state structures and mechanism only here. The Supporting Information gives all results of the alternative spin states investigated. Nevertheless, the reaction is expected to proceed through single-state-reactivity on the dominant quintet spin state surface in agreement with previous studies on pentacoordinated iron(IV)-oxo complexes reported before for nonheme iron enzymatic and model complexes (de Visser, 2006a,b, 2010; Hirao et al., 2006; Bernasconi and Baerends, 2008; Latifi et al., 2009; Ye and Neese, 2011; Ansari et al., 2013; Tang et al., 2013; Saouma and Mayer, 2014; Cantú Reinhard and de Visser, 2017a).
The quintet spin state of the iron(IV)-oxo(chloro) species in HctB has four unpaired electrons located in the , , and molecular orbitals and a virtual orbital. Thus, the and orbitals represent the antibonding interactions of the metal 3dxz/3dyz atomic orbitals with a 2px/2py orbital on the oxo group. These two orbitals are orthogonal and close in degeneracy in nonheme iron(IV)-oxo and located in the xz and yz molecular planes, whereby the z-axis is defined by the Fe–O axis. The and orbitals are both in the plane perpendicular to the Fe–O axis (xy-plane) and represent antibonding orbitals with ligands in the equatorial plane, namely the nitrogen of His111, the carboxylate of succinate and the halide atom. The final metal-type 3d molecular orbital is the orbital for the σ-type antibonding interaction along the Fe–O bond, which is virtual.
During the equilibration and MD simulation as well as the subsequent QM/MM geometry optimization the positions of the hexanoyl substrate chains have moved slightly. In particular this was the case for model 3 where the terminus moved away from the direction of the iron(IV)-oxo(chloro) group. As a result, the distance between substrate and iron(IV)-oxo(chloro) is quite large (>5 Å) in model 3. Figure 4 shows the eight different optimized geometries of the iron(IV)-oxo(chloro) reactant species as calculated with QM/MM for models 1, 2, 3 and 4 with either QM region A or AB. As can be seen, all optimized structures, regardless of the substrate binding position, give very similar optimized geometries with Fe–O distances ranging from 1.621 to 1.641 Å and Fe–Cl distances between 2.278 and 2.326 Å. These values match previous calculations on nonheme iron(IV)-oxo(chloro) and nonheme iron(IV)-oxo complexes excellently (de Visser, 2006a,b, 2010; Hirao et al., 2006, 2011; Matthews et al., 2006; Noack and Siegbahn, 2007; Bernasconi and Baerends, 2008, 2013; de Visser and Latifi, 2009; Kulik et al., 2009; Latifi et al., 2009; Pandian et al., 2009; Dey, 2010; Ye and Neese, 2011; Quesne and de Visser, 2012; Ansari et al., 2013; Liu et al., 2013; Tang et al., 2013; Usharani et al., 2013; Kumar et al., 2014; Saouma and Mayer, 2014; Senn, 2014; Huang et al., 2016; Zhao et al., 2016; Cantú Reinhard and de Visser, 2017a; Timmins et al., 2018).
Figure 4. QM/MM optimized geometries of the quintet spin iron(IV)-oxo reactant complexes as calculated at UB3LYP/BS1 level of theory in Turbomole:Charmm. Structures were optimized with either QM region A or QM region AB. Reactant structures for models 1–4 were calculated and are grouped by model: top row model 1 and model 2, bottom row model 3 and 4. Bond lengths are given in angstroms. Water molecules and aliphatic residues are shaded.
The QM/MM results from Figure 4 show that the secondary environment and the substrate positioning have little or no effect on the Fe–O and Fe–Cl distances in the iron(IV)-oxo(chloro) species although it may affect the kinetics as shown below. In addition, no major electronic differences are seen for the eight optimized structures and all converge to a quintet spin ground state. To be specific, group spin densities of all structures are close with values of around ρFe = 3.12 and ρO = 0.63. The orbital analysis and group spin densities implicate that most spin density is located along the Fe–O bond and point to an orbital occupation of for all optimized geometries. Our calculations are in agreement with experimental EPR and Mössbauer spectroscopy studies on the analogous structure in SyrB2 that also reported a high-spin ground state (Galonić Fujimori et al., 2007).
The only noticeable difference between the four substrate binding positions relates to the orientation of the active site Arg245 residue, which is locked in hydrogen bonding interactions to the dangling terminal carboxylate group of succinate via a bridging water molecule in model 1 and 2. By contrast, in model 3 the Arg245 group forms a direct salt bridge with the carboxylate of the succinate moiety and also is in close proximity to the oxo group. Finally, in structure 4, the Arg245 forms a salt bridge with Glu223, whereby it appears to close entrance channel II. It may very well be, therefore, that the Arg245 residue is involved in α-ketoglutarate binding and/or succinate release from the active site. Indeed, Arg245 is a conserved residue in most reported nonheme iron halogenases and therefore, is expected to play a key role in catalysis and/or substrate positioning. On the other hand, Glu223 is not conserved in the majority of reported nonheme iron halogenases and only found in HctB. It would be interesting to see how mutation of either Glu223 or Arg245 affects the enzyme function, chemical catalysis and product distributions. But that will have to await a future experimental study.
Next, the chemoselectivity of substrate halogenation vs. hydroxylation was investigated for the pathways described and defined as in Scheme 1. The reaction starts from the reactant complex of iron(IV)-oxo(chloro) with substrate (Re) and proceeds via a stepwise mechanism with an initial hydrogen atom abstraction transition state (TSHA) leading to an iron(III)-hydroxo complex and a substrate radical (IHA). In the next step the pathways diverge and either OH rebound (via transition state TSOH) or Cl rebound (via transition state TSCl) occurs. These mechanisms then lead to the hydroxylated products (POH) and halogenated products (PCl), respectively. Technically, the reaction apart from being chemoselective is also enantioselective where only one of the two isomers is expected. Thus, due to tight substrate binding and positioning, the hydrogen atom abstraction will be selective from only one of the two hydrogen atoms from the ω-1 position of the substrate. Hence the hydrogen atom abstraction will guide the enantioselectivity as seen before on related enzymes (Karamzadeh et al., 2010; Pratter et al., 2013; Timmins and de Visser, 2017; Timmins et al., 2017). Nevertheless, in HctB the first halogenation of substrate is followed by a second halogen transfer through the binding of another molecule of O2 and Cl− to repeat the cycle and the formation of the dihalogenated product.
We tested the mechanisms for the first halogenation step in HctB on all low lying spin states (singlet, triplet, quintet), for each of the substrate binding positions (1, 2, 3, 4) and with different QM regions (A or AB). Previously, we found the septet spin state of nonheme iron(IV)-oxo to be at least 7 kcal mol−1 higher in energy than the quintet spin state and, therefore, we did not investigate this state further (Latifi et al., 2009). However, we tested the models and methods by using different optimization techniques (basis set, DFT method, snapshot).
Firstly, changing the basis set from SV(P) to TV(P) on all atoms had little effect on the optimized geometries and lowered the hydrogen atom abstraction barrier by only a small amount. In particular, geometry optimizations at UB3LYP/BS1 and UB3LYP/BS2 of 5Re1 and 5TSHA, 1 using model 1 and the large QM region AB gave hydrogen atom abstraction barriers of 23.5 kcal mol−1 with UB3LYP/BS2 and 21.5 kcal mol−1 with UB3LYP/BS2//UB3LYP/BS1. In addition, the optimized geometries obtained with UB3LYP/BS1 vs. UB3LYP/BS2 are very similar (vide infra). Therefore, we continued with geometry optimizations at UB3LYP/BS1 level of theory for the remainder of the project.
Secondly, choosing alternative snapshots from the MD simulation also had little effect on the optimized geometries of the reactant complexes of nonheme iron complexes and did not change spin state orderings and relative energies dramatically. This probably is as nonheme iron dioxygenases often have very rigid structures where the substrate is tightly bound (Aluri and de Visser, 2007; Kumar et al., 2011; Tchesnokov et al., 2016; Faponle et al., 2017). Previous QM/MM studies on analogous nonheme iron dioxygenases indeed reproduced experimental chemoselectivities excellently using the same methods and procedures described here (Aluri and de Visser, 2007; Karamzadeh et al., 2010; Kumar et al., 2011; Pratter et al., 2013; Quesne et al., 2014, 2016a; Hernández-Ortega et al., 2015; Tchesnokov et al., 2016; Faponle et al., 2017; Timmins and de Visser, 2017; Timmins et al., 2017; Hofer and de Visser, 2018). Moreover, changing the QM region, basis set or density functional method did not change the ordering of the transition states.
Let us first focus on the hydrogen atom abstraction step, which is rate-determining in the overall reaction mechanism and optimized geometries for all four substrate binding positions are shown in Figure 5. Hydrogen atom abstraction from a quintet spin iron(IV)-oxo species with configuration leads to a sextet spin iron(III)-hydroxo with a nearby radical, or alternatively a quartet spin iron(III)-hydroxo with a nearby radical. Thus, the abstracted electron from the hydrogen atom either pairs up with the lowest electron to form a quartet spin iron(III)-hydroxo species or enters the virtual orbital to give a sextet spin iron(III)-hydroxo species. Previously, (de Visser, 2006b; Bernasconi and Baerends, 2008; Ye and Neese, 2011; Ansari et al., 2013; Tang et al., 2013; Cantú Reinhard and de Visser, 2017a) it was shown that hydrogen atom abstraction by a quintet spin pentacoordinated nonheme iron(IV)-oxo species leads to electron transfer to the metal that preferentially fills the orbital with one electron, designated the 5σ-pathway. Alternatively, the electron can move into one of the lower-lying and singly occupied π* orbitals: the so-called 5π-pathway. A spin density and molecular orbital analysis for all hydrogen atom abstraction transition states revealed that all are of 5σ-type leading to a radical intermediate with five unpaired electrons in metal-type orbitals ( ) and a down-spin electron on the substrate. Thus, in all transition states there is negative spin density on the substrate (ranging from −0.27 to −0.47) while the spin density on iron increases to 4.01–4.14. These values are in agreement with hydrogen atom abstraction reactions calculated before for nonheme iron dioxygenases and synthetic model complexes (Karamzadeh et al., 2010; Hirao et al., 2011; Pratter et al., 2013; Kumar et al., 2014; Zhao et al., 2016; Timmins and de Visser, 2017; Timmins et al., 2017). We attempted to find a transition state for the 5π-pathway for several models, but found them considerably higher in energy and we did not manage to optimize those structures with QM/MM. Filling of the orbital with one electron, as happens in the 5σ-pathway, usually means that the hydrogen atom abstraction is performed along the molecular z-axis, i.e., a linear Fe–O–C angle is seen in the transition state. However, due to constraints in substrate position and orientation and the location of amino acid residues in that region, in the protein the substrates are unable to approach the metal-oxo under this “ideal” angle. As a matter of fact, only in model 3 the substrate can approach the iron(IV)-oxo(chloro) from the top and an almost linear Fe–O–C5 angle is found, while the angle is considerably less for the other structures. As a result, the barrier is the lowest for model 3. Nevertheless, even with the transition states in the bent orientation they are still favored in the 5σ-pathway and no 5π-pathway is found. Thus, even though the transition states are not in the geometric ideal position, the lowest energy pathway is the 5σ-pathway. In a previous study on a biomimetic model complex, we investigated hydrogen atom abstraction from dihydroanthracene by a pentacoordinated iron(IV)-oxo species, where the z-axis was blocked by ligands, so that the substrate had to approach from the side. For that system, analogously to what is seen here, the 5σ-pathway was still the lowest in energy (Latifi et al., 2013). Therefore, geometric constraints cannot overcome the orbital energy differences between the 5σ− and 5π pathways for the HctB active site.
Figure 5. (A) Hydrogen atom abstraction barriers (structures and energies, in kcal mol−1) for hydrogen atom abstraction from the C5 position of the hexanoyl substrate terminus by the iron(IV)-oxo(chloro) species of HctB as calculated with QM/MM at UB3LYP/BS2/ /UB3LYP/BS1:Charmm with Turbomole using QM region AB. (B) Optimized geometries of the hydrogen atom abstraction barriers with bond lengths in angstroms, the Fe–O–C5 angle in degrees and the imaginary frequency in wave numbers.
As can be seen from Figure 5, the approach of the substrate on the iron(IV)-oxo center has a profound effect on the barrier height for hydrogen atom abstraction. In particular, pathway 3 appears to be catalytically most effective with a low hydrogen atom abstraction barrier of only 4.9 kcal mol−1 (UB3LYP/BS2//UB3LYP/BS1) for QM region AB. This is thanks to its geometric orientation, whereby substrate approach to the catalytic center is from the top and the donating hydrogen atom is aligned with the Fe–O axis. With substrates approaching through channels 1, 2, and 4, however, the Fe–O–H angle is strongly bent and the hydrogen atom abstraction barriers are dramatically raised in energy, although they still should be accessible at room temperature.
All transition state geometries display short C–H and long O–H bonds as well as relatively long O–C distances. It may very well be that the enzyme has evolved to do this on purpose in order to guide the chemoselectivity of the reaction and prevent unwanted (more favorable) side reactions. Thus, alternative nonheme iron dioxygenases, such as taurine/αKG dioxygenase have the substrate located in a pocket nearby the iron(IV)-oxo species and a QM/MM study provided a substantially lower hydrogen atom abstraction barrier than found for HctB (Borowski et al., 2004; de Visser, 2006c; Nemukhin et al., 2006; Cicero et al., 2007; Sinnecker et al., 2007; Godfrey et al., 2008; Chen et al., 2011; Bushnell et al., 2012; Mai and Kim, 2016; Wójcik et al., 2016; Álvarez-Barcia and Kästner, 2017). All hydrogen atom abstraction transition states display a large imaginary frequency for the C–H–O stretch vibration; however, distinct differences in the magnitude are seen. In particular, the barriers 5 TSHA, 2 and 5 TSHA, 3 have imaginary frequencies well over i1,000 cm−1 and hence considerable amount of tunneling can be expected. Indeed, previous computational studies on hydrogen atom abstraction reactions by iron(IV)-oxo complexes predicted significant tunneling but also a huge change in barrier upon replacement of the transferring hydrogen atom by deuterium (Kumar et al., 2003, 2004; de Visser, 2006d; Quesne et al., 2016b; Cantú Reinhard et al., 2017b). In particular, the kinetic isotope effect was shown to linearly increase with the value of the imaginary frequency in the transition state.
Table 1 summarizes all hydrogen atom abstraction barriers for models 1, 2, and 4. In general, for model 2 and 4 the UB3LYP/BS1 and UB3LYP/BS2//UB3LYP/BS1 energies are very close, whereas a change from the small QM region to the larger one leads to a drop of about 4 kcal mol−1 in the transition state energy. More dramatic changes are seen for model 1. Thus, with the small QM region large variations in the transition state energy are found dependent on the calculation method and basis set. Moreover, the energies obtained for the small QM region are quite different from those for the larger QM region, where little fluctuation upon changing the basis set is observed. Therefore, the large QM region is the most suitable chemical system for our calculations and hence we used that for the rest of the work.
Table 1. Transition state energies of 5TSHA as calculated with different QM sizes and basis sets for model 1, 2, and 4(a).
We even did a full geometry optimization at UB3LYP/BS2:Charmm with QM/MM for the same reaction mechanism using both QM region A and AB for model 1. The UB3LYP/BS2 QM/MM optimized geometry of 5TSHA, 1 is located at a value of 23.5 kcal mol−1 above the reactant complex, which is within 2 kcal mol−1 of the UB3LYP/BS1 and UB3LYP/BS2//UB3LYP/BS1 results. Moreover, the effect of the basis set on the structures of the optimized geometries is small as well and geometrically they are very close (Supporting Information). Thus, the C–H and H–O distances change from 1.22/1.45 Å calculated with UB3LYP/BS1:Charmm to 1.22/1.43 Å as obtained with UB3LYP/BS2:Charmm for model 1.
There is quite a bit of fluctuation in the hydrogen atom abstraction barrier depending on the position of the substrate (Figure 5). The lowest barrier is found for model 3 and is 4.9 kcal mol−1 in energy (UB3LYP/BS2//UB3LYP/BS1:Charmm with QM region AB), while the substrate approach via models 1, 2, and 4 are higher in energy by 16.6, 11.5, and 15.2 kcal mol−1, respectively. The latter three, although somewhat higher in energy than model 3, have barrier heights that are still within range of what would be an accessible hydrogen atom abstraction barrier at room temperature. It is very interesting to see that models 1, 2, and 4 give hydrogen atom abstraction barriers to within a couple of kcal mol−1 of each other and hence, the substrate entrance and positioning appears to have a small effect on the ability of the iron(IV)-oxo(chloro) species to react with substrate. As the substrate is connected to an acyl-carrier protein that is latched on the surface of the protein, the hydrogen atom abstraction barrier is dependent on the distance the substrate terminus can be inserted into the protein. We show here that in all four orientations (models 1, 2, 3, and 4) the substrate entrance leads to a viable hydrogen atom abstraction channel; however, the ease of hydrogen atom abstraction is strongly influenced by its positioning. Nevertheless, a low hydrogen atom abstraction barrier not necessarily correlates with the correct chemoselectivity of reaction, which is determined in the subsequent reaction step. Therefore, even though the hydrogen atom abstraction is rate determining, it does not decide the product distributions.
The hydrogen atom abstraction reaction leads to an iron(III)-hydroxo(chloro) radical intermediate (IHA) and proceeds to products through either OH or Cl rebound to the substrate radical to form the products POH and PCl. For all four substrate position orientations (1, 2, 3, and 4), we then investigated the mechanisms for either Cl rebound or OH rebound. As can be seen from Figure 6, dramatic changes in rebound as well as chemoselectivity of the reaction are obtained depending on the position of the substrate. The ordering and relative energies are not dependent on the choice of the basis set or the size of the QM region: all results point in the same direction and give the same conclusions (Supporting Information). With the substrate entering from a channel above the iron(IV)-oxo species, like through models 1 and 2, the chemoselectivity is in favor of halogenation over hydroxylation by more than 10 kcal mol−1. In particular, the halide rebound step was found to have a negligible barrier in both cases; therefore, halide transfer will be fast in both cases.
Figure 6. Energy diagram (in kcal mol−1) for the pathways leading to halogenation or hydroxylation products from the iron(III)-hydroxo(chloro) intermediate (5IHA) as calculated with QM/MM in Turbomole:Charmm. Energies (in kcal mol−1) obtained with QM/MM at UB3LYP/BS2//UB3LYP/BS1:Charmm and are the result of a full geometry optimization of the local minima and transition states. Energies are taken relative to the value of the radical intermediate state.
A complete reversal of the chemoselectivity is seen in model 3, but now the OH rebound is barrierless. By contrast, almost identical OH and Cl rebound barriers are found in model 4. Clearly, the substrate position and orientation has a major effect in the bifurcation pathways and the chemoselectivity of the reaction. Moreover, the calculations presented in Figure 6 implicate that HctB has two viable substrate entrance channels (I and II in Figure 2) and models with the substrate in either of these entrance channels give low-energy hydrogen atom abstraction barriers and preferential halogenation over hydroxylation. However, substrate entrance through channel I in HctB can lead to the location of the substrate in various positions, whereby the one in model 4 probably will give a mixture of products and position 3 gives substrate hydroxylation. As such, substrate positioning and orientation is vital for efficient halogenation in HctB and we predict the enzyme to latch the acyl carrier into substrate entrance channel II and position the substrate in the orientation as shown in model 1. The structure and orientation of the substrate in model 1 leads to a reaction mechanism with thermochemically accessible hydrogen atom abstraction barrier at room temperature and chemoselectivity preference for substrate halogenation as expected from a halogenase. Most likely, therefore, the other channel located from our structure analysis (substrate entrance channel I) could be involved in guiding halide, dioxygen or α-ketoglutarate reactants into the active site and release of succinate and CO2.
The hydrogen atom abstraction barriers reported above in Figure 5 already give a clue on the potential chemoselectivity of the reaction. Thus, in structure 5TSHA, 1 and 5TSHA, 2 the Fe–Cl bond has elongated from about 2.28 Å in the reactant complexes to 2.42 and 2.33 Å, respectively. In both cases, the substrate radical is positioned in the quadrant in between the hydroxo and chloride ligands. On the other hand, in 5TSHA, 3 the radical is located above the hydroxo group and far away from chloride, thereby making chloride rebound unlikely. To understand the large chemoselectivity preference of halogenation over hydroxylation by model 2, we show in Figure 7 the optimized geometries of the radical intermediates 5IHA as obtained with QM region A and AB. The structures are very similar regarding of the size of the QM region chosen. The Fe–O bond has elongated from 1.62 Å in 5Re2 to 1.90 (1.88) Å in 5IHA using QM region A (AB), respectively. This is due to occupation of the orbital with one electron in the radical intermediate, which is antibonding along the Fe–O bond and hence leads to its elongation. Although the radical has dissipated from the metal center slightly, it is nearer to the halide than to the hydroxo group by almost 1 Å (Figure 7). Moreover, the hydrogen atom of the hydroxo group points in the direction of the radical. However, for hydroxo rebound it needs to move out of the way to make the C–O bond formation possible. In the gas-phase and in apolar environments this O–H rotation costs little energy and hydroxo rebound tends to be low in energy as seen in many aliphatic hydroxylation studies reported previously (Ogliaro et al., 2000; Dey, 2010; Hirao et al., 2011; Bernasconi and Baerends, 2013; Liu et al., 2013; Usharani et al., 2013; Kumar et al., 2014; Zhao et al., 2016). By contrast, in the HctB models 1 and 2 the OH rebound barrier is high in energy due to the largely polar active site with lots of water molecules that prevent the OH group from rotation and keep it in a fixed orientation.
Figure 7. QM/MM optimized geometries of the radical intermediates (5IHA) of model 2 as calculated with UB3LYP/BS1:Charmm in Turbomole:Charmm. Data are given for QM region A (AB) with bond lengths in angstroms.
Recent QM/MM studies on the nonheme iron halogenase SyrB2 proposed a strong hydrogen bonding interaction between the hydroxo group in the radical intermediate with an active site Arg residue (Huang et al., 2016). In our optimized geometries of the radical intermediates (5IHA) no direct interaction was found for any substrate binding position, although in several cases a bridging hydrogen bonded water molecule was located in between, see Figure 7. The structure displayed in Figure 7 shows that the water channel that enters the active site pocket serves as a drive to prevent OH rebound and trigger an alternative reaction pathway, namely halogenation. A recent QM/MM study on a cytochrome P450 peroxygenase also had an optimized geometry for an iron(IV)-hydroxo(heme) with an active site whereby the hydroxo group underwent several hydrogen bonding interactions with crystal water molecules in the pocket (Faponle et al., 2016). Studies on the bifurcation pathways for substrate decarboxylation vs. hydroxylation also gave a preference for the energetically unfavorable decarboxylation mechanism similar to what is seen here.
Of course, it should be mentioned that the iron(II)-hydroxo complex with halogenated substrate PCl is not the final step in the catalytic cycle. One could envisage a proton transfer to generate iron(II)-water; however, our structure analysis did not characterize a proton-transfer channel, so that this pathway may not be feasible. Instead, we predict that after the formation of PCl, another chloride ion enters the substrate binding pocket and binds to the iron(II) center. A subsequent hydrogen atom abstraction by the iron(II)-hydroxo(chloro) and rebound of halide then gives the iron(II)-water resting state and a dihalogenated substrate as a product.
To gain further insight into the factors that determine the bifurcation patterns of hydroxylation vs. halogenation, we analyzed the chemoselectivity-determining transition states and investigated the role of the protein environment on their ordering and relative energies. In particular, we compare our QM/MM results with small model DFT calculations from previous work, which should give an indication how the protein constraints to the active site structure and how the long-range electrostatic interactions through the protein affect the chemoselectivity of the reaction. Specifically, small gas-phase model complexes of the nonheme iron halogenase reaction mechanism failed to find preferential halogenation over hydroxylation, (de Visser and Latifi, 2009; Kulik et al., 2009; Pandian et al., 2009; Huang et al., 2016) and only the inclusion of the full enzymatic structure and the QM/MM approach led to the correct chemoselectivity for SyrB2 (Borowski et al., 2010).
Therefore, we took the DFT optimized geometries from de Visser and Latifi (2009) calculated at UB3LYP/BS1 with a polarized continuum model with dielectric constant of ε = 5.7 included and perturb the small model complexes with external charges, electric fields and dipole moments and investigate the effects on the product distributions. As will become clear in the discussion, these external perturbations have a profound effect on the chemoselectivity of the reactions. Thus, we started with a detailed comparison of the QM/MM optimized geometries with the gas-phase DFT models from de Visser and Latifi (2009), see Figure 8. Of course, the gas-phase DFT model ignores the secondary coordination sphere of atoms and has large flexibility due to the absence of structural constraints. Indeed, we see from the optimized geometries that the substrate radical in the DFT model approaches the oxidant under an ideal angle (with little or no stereochemical repulsions of active site residues) to the iron(III)-hydroxide(halide) complex 5IHA.
Figure 8. Structural differences of QM/MM optimized geometries and gas-phase DFT models. Overlay of OH and Cl rebound transition state structures (5TSOH and 5TSCl) for model 1 (left) and model 3 (right). DFT structures (in blue) taken from de Visser and Latifi (2009) and QM/MM structures are in green. (A) Overlay of model 1 structures. (B) Overlay of model 3 structures.
Figure 8 shows overlays of 5TSCl/5TSOH structures of models 1 and 3 as compared to those obtained with the DFT model from (de Visser and Latifi, 2009). Overlays of models 2 and 4 with these structures are given in the Supporting Information and show a similar pattern. As can be seen the optimized 5TSCl, 1 structure overlaps the gas-phase DFT geometry almost perfectly with the substrate in the same location and under almost the same angle. However, the substrate location is far from the gas-phase orientation for 5TSOH, 1. Therefore, the substrate binding pocket appears to be designed to accommodate the substrate for efficient halogen transfer and not OH rebound as would be expected of a halogenase enzyme. Despite the good geometric agreement between 5TSCl, 1 and the analogous DFT model, the barrier heights relative to their precursor complexes are quite different. Thus, an almost negligible Cl transfer barrier is found for model 1 (Figure 6), whereas the DFT model gave a much larger barrier of about 10 kcal mol−1 (de Visser and Latifi, 2009). Clearly, even though the 5TSCl, 1 structures are almost the same, the QM/MM barrier is well lower in energy. Accordingly, the protein must stabilize the transition state for halogen transfer dramatically through long-range electrostatic interactions.
We then investigated the overlays of 5TSCl, 3/5TSOH, 3 with those obtained with the DFT model (Figure 8B). As can be seen, the OH rebound structure shows a reasonable match between the two models and hence rebound should stay low in energy. The overlay between the halogen transfer barriers of model 3 and the DFT model are much less good than seen for model 1 and explain why the chemoselectivity is reversed. The comparison of the gas-phase and QM/MM optimized structures, therefore, shows that the substrate should bind in a favorable position for halide rebound. In addition, as shown in Figure 7, water molecules prevent a low-energy OH-rotation pathway and lock the OH group in a non-rebound position. Consequently, substrate positioning in combination with a strong hydrogen bonding network around the hydroxo group destabilize OH rebound and drive the reaction toward halide transfer. A similar conclusion was reached by Mitchell et al. (2017) who reported experimental studies on protein engineering of the nonheme iron halogenase WelO5.
Subsequently, we explored the contribution of the protein environment to the chemoselectivity patterns. We initially investigated perturbations affecting the energetics of the 5TSCl and 5TSOH barriers for the small model complex and particularly focused on electric field effects using previously described procedures (Shaik et al., 2004). We took the gas-phase DFT model and calculated the electronic energy under the influence of an applied electric field with magnitude 0.0050, 0.0100, and 0.0150 au along each individual coordinate axis in both positive and negative directions. Thus, the x-axis is along the Fe–Cl bond, the y-axis along the Fe–OH bond and the z-axis along the Fe–N(His) bond. Relative energies of the two bifurcation transition states 5TSCl over 5TSOH were investigated with electric field effects for the small model and Figure 9 shows the ordering and relative energies of 5TSCl/5TSOH through perturbations along the x-axis, i.e., along the Fe–Cl bond. Similar calculations were also performed along the y- and z-axis, whereby the y-axis is defined along the Fe–OH bond and the z-axis is located along the Fe–N(His) bond. The results for the field along the y- and z-axis is given in the Supporting Information; however, in those cases the energies did not change by more than 3 kcal mol−1 and in all cases 5TSOH stayed well below 5TSCl. The results in Figure 9 show that without an electric field or a field along the negative x-axis, the hydroxylation pathway will be favored over the halogenation pathway by more than 7 kcal mol−1. With an electric field pointing along the positive x-axis, by contrast, a dramatic change in relative energies between the two barriers is seen and a chemoselective halogenation over hydroxylation is observed.
Figure 9. Relative energies (UB3LYP/BS2) between halide and OH rebound barriers (5TSCl vs. 5TSOH) as calculated with gas-phase DFT calculations in Gaussian with an applied electric field (with magnitude E0 in au) included. Relative energies (in italics) are in kcal mol−1 and the Table gives group charges (in atomic units) of transition states under the influence of an applied electric field. The electric field perturbations are along the x-axis (along the Fe—Cl bond).
To understand the chemoselectivity reversal through the addition of an electric field, we analyzed the unpaired group spin densities and charges of the individual complexes. The group charges under the addition of an electric field to the transition states are given at the bottom of Figure 9. A clear trend emerges from these calculations. Thus, in both TS's charge density (Q) is removed from Cl upon increasing the electric field along the positive x-axis. In particular, the negative charge of the chloride atom decreases from QCl = −0.61 with a field of E0, X = −0.015 au to a value of QCl = −0.31 with a field of +0.015 au, while an increase of the negative charge with 0.18 units is seen for the same field strengths for 5TSCl. At the same time, the charge on the OH group stays virtually the same in 5TSOH, whereas it increases from −0.32 to −0.46 in 5TSCl. Therefore, an applied electric field reduces the charge on chloride and puts more radical character on Cl and enables a Cl∙ transfer. At the same time, charge is transferred to iron that is reduced from iron (III) to iron (II).
As a result the 5TSCl barrier drops below the 5TSOH barrier when a positive electric field is applied. It also appears that the 5TSOH barrier is more strongly affected than the 5TSCl barrier to an external electric field and hence the chemoselectivity switch can occur due to destabilization of the hydroxyl rebound barrier. The calculation with an electric field of 0.015 au along the positive x-axis shows that the substrate group loses most of its spin density with respect to the unperturbed system, namely the spin density on the substrate moiety drops from ρSub = −0.68 for the system without an electric field to ρSub = −0.12 with an electric field of E0, X = +0.015 au. At the same time also the Cl and OH groups lose radical character but to a much lesser extent than the substrate. These results are in line with our previous work on the chemoselectivity of aliphatic hydroxylation vs. epoxidation by a cytochrome P450 model Compound I complex, where we also showed a chemoselectivity switch upon addition of an electric field along a specific axis (Shaik et al., 2004). Further, DFT calculations on a Compound I model of cytochrome c peroxidase with an applied electric field identified radical character on the porphyrin in one direction, whereas with the electric field in the opposite direction a tryptophan radical was obtained (de Visser, 2005).
Our studies, therefore, highlight the importance of the shape and size of the protein of HctB that create an induced electric field on the active site and thereby favorably lead to substrate halogenation to substrate halogenation of the substrate while minimizing substrate hydroxylation. Evidently, an applied electric field from the dipole moment of the protein affects the charge distributions and spin densities of intermediates and transition states in the reaction mechanism and thereby affects the product distributions. In the HctB halogenase enzyme the protein perturbs the active site structure with an induced electric field in such a way that it stabilizes the halogenation pathway and lowers it below the hydroxylation pathway. This effect is the result of considerable destabilization of the OH rebound barrier and a smaller stabilization of the Cl rebound as based upon the changes of the charge distributions between no field and a positive electric field.
To find further evidence of environmental perturbations influencing the halogenase vs. hydroxylase activity of the protein, we searched for charged residues in the vicinity of the iron(IV)-oxo(chloride) complex in the pdb structures of 5Re1, 5Re2, 5Re3, and 5Re4. Nearby the halide atom we locate the Glu223 and Arg245 side chains, whereby the latter may be involved in α-ketoglutarate binding and succinate release. The top part of Figure 10 gives the relative orientations of the Glu223 and Arg245 amino acids with respect to the iron(IV)-oxo(chloride) in the reactant complexes. As can be seen, the three enzyme models show considerable differences in the position of these amino acids that may incur Coulombic interactions with the active site atoms resulting in chemoselectivity changes. Thus, in model 1 the carboxylate group of Glu223 is aligned to the axis that goes through the Fe–Cl bond and appears to be set up to push electron density away from Cl and onto iron. By contrast, in 4 the Glu223 and Arg245 side chains have formed a salt bridge and now the positively charged Arg is on the axis of the Fe–Cl bond. Consequently, in model 4 charge density will be withdrawn from iron and move to chloride instead making it more difficult to transfer the Cl to substrate. In model 3 the Glu223 is moved to the outside part of the protein and hence is located at a relatively large distance from the active center. In addition, the Arg245 is positioned much closer to the halide but on the side rather than along the axis. Indeed, model 3 shows a regioselective hydroxylation and no halogenation, whereas in model 1 a preference for halogenation over hydroxylation is seen.
Figure 10. Top: Relative positions of Glu223 and Arg245 with respect to the Fe–Cl axis in 5Re optimized geometries with QM/MM. Optimized geometries calculated at UB3LYP/BS1: Charmm with Turbomole:Charmm. Structure displayed along the Fe–Cl axis by looking through the Fe–O bond. Bottom: Group charges (in atomic units) of the halogenation and hydroxylation transition states of the model system upon addition of a point charge (Q) in the position of the oxygen atoms of Glu223 in model 1.
To find out whether the Glu223 amino acid residue is near enough to the iron(IV)-oxo(chloride) to affect the charge distributions of the atoms, we took the DFT models of 5TSOH and 5TSCl and included a point charge to the system in the position of the carboxylate oxygen atom of Glu223 from model 1. Single point calculations with a point charge of Q = +1 and −1 were performed and the results are given at the bottom of Figure 10. Surprisingly, only small electronic changes are observed to 5TSCl upon addition of a point charge. By contrast, the hydroxyl rebound barrier sees major changes in the charge of the halogen atom under the influence of a point charge. Specifically, the halogen charge increases from −0.48 to −0.41 with a point charge with magnitude −1 included. In addition, the barrier 5TSCl drops below 5TSOH by 4.6 kcal mol−1, whereas the energy difference is −20.0 kcal mol−1 with a point charge of magnitude +1 included. Therefore, a point charge located at a distance of over 5 Å from the reaction center can incur an electrostatic perturbation that affects the charge distributions in the catalytic reaction center and consequently the chemoselectivity of the reaction. This is an important caveat, given that contrary to the active site and core region of the halogenase domain, the structure of HctB is less well defined by the model, with the positioning of the two other domains and interactions of the halogenase domain not clear/not well defined.
To finally test the effect of the salt bridge between Arg245 and Glu223 we performed a series of electrostatic potential calculations on model 4 whereby in one case the salt bridge was rearranged so that either Arg245 or Glu223 is closer to the chloride. Calculations of electrostatic potential fields were made with the Protein-sol server software package (Hebditch et al., 2017) and the results are given in Figure 11, displayed using PyMOL (Schrödinger, LLC). Thus, the salt-bridge between Arg245 and Glu223 in model 4 (lower panel of Figure 11) induces a change in negative charge on the halide that is pushed toward the iron. Swapping the salt-bridge and changing the position of the Arg and Glu residues has a dramatic effect on the charge donation to the active site and now electron density is withdrawn from iron and moved toward the chloride group. Therefore, the salt-bridge is essential for creating the right charge distribution between the iron and chloride group and induces an electric field effect on the active site that enables efficient substrate halogenation as predicted by the DFT studies reported above.
Figure 11. Induced electric field of the salt-bridge between Arg245 and Glu223 on the iron(IV)-oxo(chloride) active site. Electrostatic potential was calculated for model 5Re4 (A) and model 5Re4,mod (B), which is model 4 with rotated Arg245 and Glu223 residues with respect to the Cl–Fe axis. Color coding for the amino stick representations changes from red (negative potential) to blue (positive potential). The direction of the induced electric field is given by the arrow.
In summary, the QM/MM and DFT calculations presented here highlight that for effective substrate halogenation, radical character on the halogen atom is needed that can pair up with the substrate radical. Perturbations from the protein, i.e. from nearby located anionic amino acid residues, e.g., Glu223, can lead to a push-effect of charge density from the halogen to the metal to achieve this. At the same time, the active site Arg245 should not be aligned with the Fe–Cl bond. Our observations are in excellent agreement with recent computational studies on SyrB2 using model complexes (Wong et al., 2013; Srnec et al., 2016; Srnec and Solomon, 2017). In particular, it was proposed that ideal orbital overlap between the substrate radical and halogen atom would lower the halogen transfer transition state below the OH rebound barrier and thereby overcome the large hydroxylation driving force.
In this work we describe a detailed MM, MD, QM/MM, and DFT study on the nonheme iron halogenase HctB for the first time. Our initial analysis of a structure from the literature revealed several substrate entrance channels. We, therefore, created four models with the substrate in different positions and ran QM/MM calculations on snapshots from these MD simulations. We confirm a reaction mechanism starting from an iron(IV)-oxo species that reacts via hydrogen atom abstraction followed by halogen rebound. Our calculated hydrogen atom abstraction barriers vary dramatically with substrate binding position and a low barrier is found for model 3, while much higher barriers are seen for models 1 and 2. Interestingly, two models give preferential halogenation, one model preferential hydroxylation and the fourth one is expected to give a mixture of products. The two models that give preferential halogenation in our analyses use distinct substrate access channels. We analyzed the individual structures and identify an active site Glu residue (Glu223), which is unique for HctB as other halogenases typically have a Ser or Lys in that position. Finally, we compared the QM/MM calculations with small gas-phase model complexes and find that even though the gas-phase structure would imply preferential hydroxylation over halogenation, a complete chemoselectivity reversal can be achieved through external perturbations. Calculations on model systems using either an applied electric field along the positive x-axis or with a point charge with magnitude Q = −1 give dominant halogenation, whereas in all other cases, hydroxylation is predicted. It follows, therefore, that the HctB enzyme structure is designed in such a way, so as to destabilize the hydroxylation pathway and give favorable halogenation products. Overall, our calculations show that substrate binding position is essential for an optimal halogenation reaction. This may come at a cost through higher hydrogen atom abstraction barriers, but the enlarged hydrogen atom abstraction barriers may hamper optimal hydroxo rebound and the formation of alcohol products.
QM/MM calculations were performed on the mechanism of substrate hydroxylation vs. halogenation in the HctB halogenase using methods and procedures utilized previously on analogous enzymatic systems (Aluri and de Visser, 2007; Kumar et al., 2011; Hernández-Ortega et al., 2015; Quesne et al., 2016a; Tchesnokov et al., 2016; Faponle et al., 2017; Hofer and de Visser, 2018).
The work started from the halogenase domain reported by Pratter et al. (2014a,b) which contains an in silico docked hexanoyl phosphopentatheinyl (PPT) moiety (channel II, Figure 2). Upon our analysis of the structure; however, another substrate entrance channel was discovered similar to that identified in SyrB2 (channel I, Figure 2) (Blasiak et al., 2006). We manually docked the substrate into each substrate entrance channel and replaced the iron(II)-water group by iron(IV)-oxo and shortened the Fe–O distance to 1.63 Å. In addition, α-ketoglutarate was replaced by succinate (Succ) and the acyl protein carrier was removed. The substrate in channel I was positioned in three different orientations relative to the iron(IV)-oxo. The substrates were positioned in such a way that stereochemical clashes with protein residues along the substrate entrance channels were minimized and after insertion of the substrate the system was solvated. During solvation and molecular dynamics simulations no additional water molecules entered the active site pocket. The structures were charge neutralized with 18 chloride and 13 magnesium ions on the surface of the protein. For example, model 1 contained a total of 32,458 atoms including 9,157 water molecules.
Hydrogen atoms were added to models 1–4 using the PDB2PQR software package at pH = 7, (Dolinsky et al., 2007) whereby all Glu/Asp side chains were in their basic form and all Arg/Lys residues in the acidic forms. The protonation states of histidine residues was decided upon visual inspection of their local environment and His67, His111, His157, His185, His227, and His256 were singly protonated at the Nδ position, whereas His64, His163, His268, and His273 were protonated at the Nε position. Histidine residue His8 was chosen as being doubly protonated at both the Nδ- and Nε-positions.
A constraint geometry optimization with all protein backbone atoms fixed was performed in Charmm, (Brooks et al., 1983) and subsequently the system was solvated in a sphere with a radius of 40 Å. An iterative solvation procedure (Supporting Information Figure S1) with fixed protein backbone was followed by an equilibration and heating procedure to 298 K. For all chemical systems described here, a full molecular dynamics simulation with the Charmm forcefield was performed without constraints for a minimum of 10 ns (Supporting Information Figure S2). All systems relaxed to stable structures and a random low-energy snapshot was chosen as the starting geometry for the QM/MM calculations.
We selected a QM region containing the iron(IV)-oxo(chloro) group and included the imidazole groups of His111 and His227, the acetate terminus of succinate (Succ) and the thiohexanoic acid arm of the substrate as our minimal QM region A. In addition, a larger QM region was tested that included the full or partial amino acid R-groups within 6 Å of the iron(IV)-oxo(chloro) structure: QM region AB. For model 1, the large QM region AB contained the amino acid side chains of Ile108, Val113, Asp200, Asp202, Glu223, Val225, Met226, Arg245, and six water molecules. QM region AB for Model 2 is described with the same amino acids except Asp202, Val225, and Met245 and in addition also includes Ser141 and Phe221 and four water molecules. Model 3 in its most elaborate form covers Ile108, Val113, Ser141, Asp200, Phe221, Asn243, Arg245, and four water molecules. Finally, the large QM region for Model 4 includes Ile108, Val113, Phe221, Arg245, and five water molecules.
Subsequent QM/MM calculations of the different substrate bound iron(IV)-oxo(chloro) enzyme intermediates were performed in Turbomole:Charmm and linked via ChemShell (Ahlrichs et al., 1989; Smith and Forester, 1996; Sherwood et al., 2003). Density functional theory methods at the unrestricted B3LYP level of theory (Lee et al., 1988; Becke, 1993) were applied to the QM region in Turbomole with an SV(P) basis set on all atoms: basis set BS1 (Schafer et al., 1992). Single point calculations on the optimized structures in QM/MM were performed with an all-electron Wachters-type basis set on iron and def2-TZV(P) on the rest of the atoms: basis set BS2 (Wachters, 1970; Hay, 1977; Krishnan et al., 1980; Bauschlicher et al., 1989). For several systems we tested a full geometry optimization of using basis set BS2 but found similar structures and relative energies to those obtained with BS1; hence BS1 was used for all systems. The MM region was described with the Charmm forcefield (Brooks et al., 1983). A detailed benchmark study comparing computational free energies of activation with experimental data from Eyring plots gave good agreement between experiment and theory for the methods described here (Cantú Reinhard et al., 2016a). Furthermore, these methods were also seen to excellently reproduce reduction potentials of copper proteins, (Fowler et al., 2017) and regio- and chemoselectivities of reaction mechanisms (Jastrzebski et al., 2014; Barman et al., 2016b; Yang et al., 2016).
Geometries were optimized at UB3LYP/BS1 in Turbomole:Charmm and a frequency calculation on the QM region only was done to test whether it was a local minimum or first order saddle point. Geometry scans with one degree of freedom fixed were performed along a specified reaction coordinate. The maximum of these scans were subjected to a full geometry optimization of a first-order saddle point and the starting and end points of the scans established the proposed reaction paths.
AT and SdV devised the project. AT ran the QM/MM calculations. NF and JW performed the electrostatic calculations. AT, NF, GS, and SdV wrote the paper.
The BBSRC is acknowledged for a studentship for AT and NF under grant code BB/J014478/1. The EU-COST Network Explicit Control Over Spin-states in Technology and Biochemistry (ECOSTBio, CM1305) is acknowledged for support.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00513/full#supplementary-material
Agarwal, V., Miles, Z. D., Winter, J. M., Eustáquio, A. S., El Gamal, A. A., and Moore, B. S. (2017). Enzymatic halogenation and dehalogenation reactions: pervasive and mechanistically diverse. Chem. Rev. 117, 5619–5674. doi: 10.1021/acs.chemrev.6b00571
Ahlrichs, R., Bär, M., Häser, M., Horn, H., and Kölmel, C. (1989). Electronic structure calculations on workstation computers: the program system turbomole. Chem. Phys. Lett. 162, 165–169. doi: 10.1016/0009-2614(89)85118-8
Álvarez-Barcia, S., and Kästner, J. (2017). Atom tunneling in the hydroxylation process of taurine/α-ketoglutarate dioxygenase identified by quantum mechanics/molecular mechanics simulations. J. Phys. Chem. B. 121, 5347–5354. doi: 10.1021/acs.jpcb.7b03477
Ansari, A., Kaushik, A., and Rajaraman, G. (2013). Mechanistic insights on the ortho-hydroxylation of aromatic compounds by non-heme iron complex: a computational case study on the comparative oxidative ability of ferric-hydroperoxo and high-valent FeIV-O and FeV-O intermediates. J. Am. Chem. Soc. 135, 4235–4249. doi: 10.1021/ja307077f
Barman, P., Faponle, A. S., Vardhaman, A. K., Angelone, D., A.-,Löhr, M., Browne, W. R., et al. (2016b). Influence of ligand architecture in tuning reaction bifurcation pathways for chlorite oxidation by non-heme iron complexes. Inorg. Chem. 55, 10170–10181. doi: 10.1021/acs.inorgchem.6b01384
Barman, P., Upadhyay, P., Faponle, A. S., Kumar, J., Nag, S. S., Kumar, D., et al. (2016a). Deformylation reaction by a nonheme manganese(III)–peroxo complex via initial hydrogen-atom abstraction. Angew. Chem. Int. Ed. 55, 11091–11095. doi: 10.1002/anie.201604412
Bauschlicher, C. W. Jr., Langhoff, S. R., Partridge, H., and Barnes, L. A. (1989). Theoretical studies of the first- and second-row transition-metal methyls and their positive ions. J. Chem. Phys. 91, 2399–2411. doi: 10.1063/1.456998
Bernasconi, L., and Baerends, E. J. (2008). The EDTA complex of oxidoiron(IV) as realisation of an optimal ligand environment for high activity of FeO2+. Eur. J. Inorg. Chem. 2008, 1672–1681. doi: 10.1002/ejic.200701135
Bernasconi, L., and Baerends, J. (2013). A frontier orbital study with ab initio molecular dynamics of the effects of solvation on chemical reactivity: solvent-induced orbital control in FeO-activated hydroxylation reactions. J. Am. Chem. Soc. 135, 8857–8867. doi: 10.1021/ja311144d
Blasiak, L. C., Vaillancourt, F. H., Walsh, C. T., and Drennan, C. L. (2006). Crystal structure of the non-haem iron halogenase SyrB2 in syringomycin biosynthesis. Nature 440, 368–371. doi: 10.1038/nature04544
Bollinger, J M. Jr., Price, J. C, Hoffart, L. M., Barr, E. W., and Krebs, C. (2005). Mechanism of taurine: α-ketoglutarate dioxygenase (TauD) from Escherichia coli. Eur. J. Inorg. Chem. 2005, 4245–4254. doi: 10.1002/ejic.200500476
Borowski, T., Bassan, A., and Siegbahn, P. E. M. (2004). Mechanism of dioxygen activation in 2-oxoglutarate-dependent enzymes: a hybrid DFT study. Chem. Eur. J. 10, 1031–1041. doi: 10.1002/chem.200305306
Borowski, T., Noack, H., Radon, M., Zych, K., and Siegbahn, P. E. M. (2010). Mechanism of selective halogenation by SyrB2: a computational study. J. Am. Chem. Soc. 132, 12887–12898. doi: 10.1021/ja101877a
Brooks, B. R., Bruccoleri, R. E., Olafson, B. D., States, D. J., Swaminathan, S., and Karplus, M. (1983). CHARMM: a program for macromolecular energy, minimization, and dynamics calculation. J. Comput. Chem. 4, 187–217. doi: 10.1002/jcc.540040211
Bruijnincx, P. C. A., van Koten, G., and Klein Gebbink, R. J. M. (2008). Mononuclear non-heme iron enzymes with the 2-His-1-carboxylate facial triad: recent developments in enzymology and modeling studies. Chem. Soc. Rev. 37, 2716–2744. doi: 10.1039/b707179p
Buongiorno, D., and Straganz, G. D. (2013). Structure and function of atypically coordinated enzymatic mononuclear non-heme-Fe(II) centers. Coord. Chem. Rev. 257, 541–563. doi: 10.1016/j.ccr.2012.04.028
Bushnell, E. A. C., Fortowsky, G. B., and Gauld, J. W. (2012). Model iron–oxo species and the oxidation of imidazole: insights into the mechanism of OvoA and EgtB? Inorg. Chem. 51, 13351–13356. doi: 10.1021/ic3021172
Cantú Reinhard, F. G., Barman, P., Mukherjee, G., Kumar, J., Kumar, D., Kumar, D., et al. (2017b). Keto–enol tautomerization triggers an electrophilic aldehyde deformylation reaction by a nonheme manganese(III)-peroxo complex. J. Am. Chem. Soc. 139, 18328–18338. doi: 10.1021/jacs.7b10033
Cantú Reinhard, F. G., and de Visser, S. P. (2017a). Oxygen atom transfer using an iron(IV)-oxo embedded in a tetracyclic N-heterocyclic carbene system: how does the reactivity compare to cytochrome P450 compound I? Chem. Eur. J. 23, 2935–2944. doi: 10.1002/chem.201605505
Cantú Reinhard, F. G., Faponle, A. S., and de Visser, S. P. (2016a). Substrate Sulfoxidation by an Iron(IV)-Oxo complex: benchmarking computationally calculated barrier heights to experiment. J. Phys. Chem. A. 120, 9805–9814. doi: 10.1021/acs.jpca.6b09765
Cantú Reinhard, F. G., Sainna, M. A., Upadhyay, P., Balan, G. A., Kumar, D., Fornarini, S., et al. (2016b). A systematic account on aromatic hydroxylation by a cytochrome P450 model compound i: a low-pressure mass spectrometry and computational study. Chem. Eur. J. 22, 18608–18619. doi: 10.1002/chem.201604361
Chatterjee, S., and Paine, T. K. (2016). Hydroxylation versus halogenation of aliphatic C-H bonds by a dioxygen-derived iron-oxygen oxidant: functional mimicking of iron halogenases. Angew. Chem. Int. Ed. 55, 7717–7722. doi: 10.1002/anie.201509914
Chen, X., and van Pée, H. (2008). Catalytic mechanisms, basic roles, and biotechnological and environmental significance of halogenating enzymes. Acta. Biochim. Biophys. Sin. 40, 183–193. doi: 10.1111/j.1745-7270.2008.00390.x
Cicero, G., Carbonera, C., Valegård, K., Hajdu, J., Andersson, I., and Ranghino, G. (2007). Study of the oxidative half-reaction catalyzed by a non-heme ferrous catalytic center by means of structural and computational methodologies. Int. J. Quant. Chem. 107, 1514–1522. doi: 10.1002/qua.21275
Comba, P., and Wunderlich, S. (2010). Iron-catalyzed halogenation of alkanes: modeling of nonheme halogenases by experiment and DFT calculations. Chem. Eur. J. 16, 7293–7299. doi: 10.1002/chem.201000092
Costas, M., Mehn, M. P., Jensen, M. P., and Que Jr, L. (2004). Dioxygen activation at mononuclear nonheme iron active sites: enzymes, models, and intermediates. Chem. Rev. 104, 939–986. doi: 10.1021/cr020628n
de Visser, S. P. (2006a). Differences in and comparison of the catalytic properties of heme and non-heme enzymes with a central Oxo–iron group. Angew. Chem. Int. Ed. 45, 1790–1793. doi: 10.1002/anie.200503841
de Visser, S. P. (2006b). What factors influence the ratio of CH hydroxylation versus CC epoxidation by a nonheme cytochrome P450 biomimetic? J. Am. Chem. Soc. 128, 15809–15818. doi: 10.1021/ja065365j
de Visser, S. P. (2006c). Propene activation by the Oxo-iron active species of taurine/α-ketoglutarate dioxygenase (TauD) enzyme. How Does the Catalysis Compare to Heme-Enzymes? J. Am. Chem. Soc. 128, 9813–9824. doi: 10.1021/ja061581g
de Visser, S. P. (2006d). Substitution of Hydrogen by deuterium changes the regioselectivity of ethylbenzene hydroxylation by an oxo–iron–porphyrin catalyst. Chem. Eur. J. 12, 8168–8177. doi: 10.1002/chem.200600376
de Visser, S. P. (2010). Trends in substrate hydroxylation reactions by Heme and nonheme iron(IV)-oxo oxidants give correlations between intrinsic properties of the oxidant with barrier height. J. Am. Chem. Soc. 132, 1087–1097. doi: 10.1021/ja908340j
de Visser, S. P., and Latifi, R. (2009). Carbon dioxide: a waste Product in the catalytic cycle of α-ketoglutarate dependent halogenases prevents the formation of hydroxylated by-products. J. Phys. Chem. B. 113, 12–14. doi: 10.1021/jp8097632
de Visser, S. P., Quesne, M. G., Martin, B., Comba, P., and Ryde, U. (2014). Computational modelling of oxygenation processes in enzymes and biomimetic model complexes. Chem. Commun. 50, 262–282. doi: 10.1039/C3CC47148A
Dey, A. (2010). Density functional theory calculations on the mononuclear non-heme iron active site of hmd hydrogenase: role of the internal ligands in tuning external ligand binding and driving H2 heterolysis. J. Am. Chem. Soc. 132, 13892–13901. doi: 10.1021/ja1041918
Dolinsky, T. J., Czodrowski, P., Li, H., Nielsen, J. E., Jensen, J. H., Klebe, G., et al. (2007). PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations. Nucl. Acids Res. 35, 522–525. doi: 10.1093/nar/gkm276
Faponle, A. S., Quesne, M. G., and de Visser, S. P. (2016). Origin of the regioselective fatty-acid hydroxylation versus decarboxylation by a cytochrome p450 peroxygenase: what drives the reaction to biofuel production? Chem. Eur. J. 22, 5478–5483. doi: 10.1002/chem.201600739
Faponle, A. S., Seebeck, F. P., and de Visser, S. P. (2017). Sulfoxide synthase versus cysteine dioxygenase reactivity in a nonheme iron enzyme. J. Am. Chem. Soc. 139, 9259–9270. doi: 10.1021/jacs.7b04251
Fowler, N. J., Blanford, C. F., Warwicker, J., and de Visser, S. P. (2017). Prediction of reduction potentials of copper proteins with continuum electrostatics and density functional theory. Chem. Eur. J. 23, 15436–15445. doi: 10.1002/chem.201702901
Galonić Fujimori, D., Barr, E. W., Matthews, M. L., Koch, G. M., Yonce, J. R., Walsh, C. T., et al. (2007). Spectroscopic evidence for a high-spin Br-Fe(IV)-Oxo intermediate in the α-ketoglutarate-dependent halogenase CytC3 from Streptomyces. J. Am. Chem. Soc. 129, 13408–13409. doi: 10.1021/ja076454e
Godfrey, E., Porro, C. S., and de Visser, S. P. (2008). Comparative quantum mechanics/molecular mechanics (QM/MM) and density functional theory calculations on the Oxo–iron species of taurine/α-ketoglutarate dioxygenase. J. Phys. Chem. A. 112, 2464–2468. doi: 10.1021/jp710999v
Hebditch, M., Carballo-Amador, M. A., Charonis, S., Curtis, R., and Warwicker, J. (2017). Protein–Sol: a web tool for predicting protein solubility from sequence. Bioinformatics 33, 3098–3100. doi: 10.1093/bioinformatics/btx345
Hernández-Ortega, A., Quesne, M. G., Bui, S., Heyes, D. J., Steiner, R. A., Scrutton, N. S., et al. (2015). Catalytic mechanism of cofactor-free dioxygenases and how they circumvent spin-forbidden oxygenation of their substrates. J. Am. Chem. Soc. 137, 7474–7487. doi: 10.1021/jacs.5b03836
Hirao, H., Li, F., Que Jr, L., and Morokuma, K. (2011). Theoretical study of the mechanism of oxoiron(IV) formation from H2O2 and a nonheme iron(II) complex: O–O cleavage involving proton-coupled electron transfer. Inorg. Chem. 50, 6637–6648. doi: 10.1021/ic200522r
Hofer, T. S., and de Visser, S. P. (2018). Editorial: quantum mechanical/molecular mechanical approaches for the investigation of chemical systems – recent developments and advanced applications. Front. Chem. 6:357. doi: 10.3389/fchem.2018.00357
Huang, J., Li, C., Wang, B., Sharon, D. A., Wu, W., and Shaik, S. (2016). Selective chlorination of substrates by the halogenase SyrB2 is controlled by the protein according to a combined quantum mechanics/molecular mechanics and molecular dynamics study. ACS Catal. 6, 2694–2704. doi: 10.1021/acscatal.5b02825
Jastrzebski, R., Quesne, M. G., Weckhuysen, B. M., de Visser, S. P., and Bruijnincx, P. C. A. (2014). Experimental and computational evidence for the mechanism of intradiol catechol dioxygenation by non-heme iron(III) complexes. Chem. Eur. J. 20, 15686–15691. doi: 10.1002/chem.201404988
Ji, L., Faponle, A. S., Quesne, M. G., Sainna, M. A., Zhang, J., Franke, A., et al. (2015). Drug metabolism by cytochrome P450 enzymes: what distinguishes the pathways leading to substrate hydroxylation over desaturation? Chem. Eur. J. 21, 9083–9092. doi: 10.1002/chem.201500329
Kaczmarek, M. A., Malhotra, A., Balan, G. A., Timmins, A., and de Visser, S. P. (2018). Nitrogen reduction to ammonia on a biomimetic mononuclear iron centre: insights into the nitrogenase enzyme. Chem. Eur. J. 24, 5293–5302. doi: 10.1002/chem.201704688
Karamzadeh, B., Kumar, D., Sastry, G. N., and de Visser, S. P. (2010). Steric factors override thermodynamic driving force in regioselectivity of proline Hydroxylation by Prolyl-4-hydroxylase enzymes. J. Phys. Chem. A. 114, 13234–13243. doi: 10.1021/jp1089855
Kim, S. H., Perera, R., Hager, L. P., Dawson, J. H., and Hoffman, B. M. (2006). Rapid freeze-quench ENDOR study of chloroperoxidase compound I: the site of the radical. J. Am. Chem. Soc. 128, 5598–5599. doi: 10.1021/ja060776l
Krishnan, R., Binkley, J. S., Seeger, R., and Pople, J. A. (1980). Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions. J. Chem. Phys. 72, 650–654. doi: 10.1063/1.438955
Kumar, D., de Visser, S. P., and Shaik, S. (2003). How does product isotope effect prove the operation of a two-state “rebound” mechanism in C–H hydroxylation by cytochrome P450? J. Am. Chem. Soc. 125, 13024–13025. doi: 10.1021/ja036906x
Kumar, D., de Visser, S. P., Sharma, P. K., Cohen, S., and Shaik, S. (2004). Radical clock substrates, their C–H hydroxylation mechanism by cytochrome P450, and other reactivity patterns: what does theory reveal about the clocks' behavior? J. Am. Chem. Soc. 126, 1907–1920. doi: 10.1021/ja039439s
Kumar, D., Thiel, W., and de Visser, S. P. (2011). Theoretical study on the mechanism of the oxygen activation process in cysteine dioxygenase enzymes. J. Am. Chem. Soc. 133, 3869–3882. doi: 10.1021/ja107514f
Kumar, S., Faponle, A. S., Barman, P., Vardhaman, A. K., Sastri, C. V., Kumar, D., et al. (2014). Long-range electron transfer triggers mechanistic differences between Iron(IV)-Oxo and Iron(IV)-imido oxidants. J. Am. Chem. Soc. 136, 17102–17115. doi: 10.1021/ja508403w
Latifi, R., Bagherzadeh, M., and de Visser, S. P. (2009). Origin of the correlation of the rate constant of substrate hydroxylation by nonheme iron(IV)–oxo complexes with the bond-dissociation energy of the C-H bond of the substrate. Chem. Eur. J. 15, 6651–6662. doi: 10.1002/chem.200900211
Latifi, R., Sainna, M. A., Rybak-Akimova, E. V., and de Visser, S. P. (2013). Does hydrogen-bonding donation to manganese(IV)–Oxo and iron(IV)–Oxo oxidants affect the oxygen-atom transfer ability? A computational study. chem. Eur. J. 19, 4058–4068. doi: 10.1002/chem.201202811
Liu, Y., Guan, X. E. L., Wong, M., Liu, P. J., Huang, S. C., and Che, M. (2013). Nonheme Iron-mediated amination of C(sp3)–H bonds. Quinquepyridine-supported iron-imide/nitrene intermediates by experimental studies and DFT calculations. J. Am. Chem. Soc. 135, 7194–7204. doi: 10.1021/ja3122526
Mai, B. K., and Kim, Y. (2016). Is it Fe(III)-oxyl radical that abstracts hydrogen in the C–H activation of TauD? A theoretical study based on the DFT potential energy surfaces. Inorg. Chem. 55, 3844–3852. doi: 10.1021/acs.inorgchem.5b02939
Martinez, J. S., Carroll, G. L., Tschirret-Guth, R. A., Altenhoff, G., Little, D. R., and Butler, A. (2001). On the Regiospecificity of Vanadium Bromoperoxidase. J. Am. Chem. Soc. 123, 3289–3294. doi: 10.1021/ja004176c
Matthews, M. L., Krest, C. M., Barr, E. W., Vaillancourt, F. H., Walsh, C. T., Green, M. T., et al. (2009). Substrate-triggered formation and remarkable stability of the C–H bond-cleaving chloroferryl intermediate in the aliphatic halogenase, SyrB2. Biochemistry 48, 4331–4343. doi: 10.1021/bi900109z
Matthews, M. L., Neumann, C. S., Miles, L. A., Grove, T. L., Booker, S. J., Krebs, C., et al. (2006). Substrate positioning controls the partition between halogenation and hydroxylation in the aliphatic halogenase, SyrB2. Proc. Natl. Acad. Sci. USA. 106, 17723–17728. doi: 10.1073/pnas.0909649106
Messerschmidt, A., and Wever, R. (1996). X-ray structure of a vanadium-containing enzyme: chloroperoxidase from the fungus Curvularia inaequalis. Proc. Natl. Acad. Sci. U.S.A. 93, 392–396. doi: 10.1073/pnas.93.1.392
Mitchell, A. J., Dunham, N. P., Bergman, J. A., Wang, B., Zhu, Q., Chen, W., et al. (2017). Structure-guided reprogramming of a hydroxylase to halogenate its small molecule substrate. Biochemistry 56, 441–444. doi: 10.1021/acs.biochem.6b01173
Neidig, M. L., Brown, C. D., Light, K. M., Galonić Fujimori, D., Nolan, E. M., Price, J. C., et al. (2007). CD and MCD of CytC3 and taurine dioxygenase: role of the facial triad in α-KG-dependent oxygenases. J. Am. Chem. Soc. 129, 14224–14231. doi: 10.1021/ja074557r
Nemukhin, A. V., Topol, I. A., Cachau, R. E., and Burt, S. K. (2006). On the nature of oxoiron (IV) intermediate in dioxygen activation by non-heme enzymes. Theor. Chem. Acc. 115, 348–353. doi: 10.1007/s00214-005-0047-3
Noack, H., and Siegbahn, P. E. M. (2007). Theoretical investigation on the oxidative chlorination performed by a biomimetic non-heme iron catalyst. J. Biol. Inorg. Chem. 12, 1151–1162. doi: 10.1007/s00775-007-0284-0
Ogliaro, F., Harris, N., Cohen, S., Filatov, M., de Visser, S. P., and Shaik, S. (2000). A model “rebound” mechanism of hydroxylation by cytochrome P450: stepwise and effectively concerted pathways, and their reactivity patterns. J. Am. Chem. Soc. 122, 8977–8989. doi: 10.1021/ja991878x
Pandian, S., Vincent, M. A., Hillier, I. H., and Burton, N. A. (2009). Why does the enzyme SyrB2 chlorinate, but does not hydroxylate, saturated hydrocarbons? A density functional theory (DFT) study. Dalton Trans. 2009, 6201–6207. doi: 10.1039/b906866j
Pratter, S. M., Ivkovic, J., Birner-Gruenberger, R., Breinbauer, R., Zangger, K., and Straganz, G. D. (2014a). More than just a halogenase: modification of fatty acyl moieties by a trifunctional metal enzyme. ChemBioChem. 15, 567–574. doi: 10.1002/cbic.201300345
Pratter, S. M., Konstantinovics, C., DiGiuro, C. L. M, Leitner, E., Kumar, D., de Visser, S. P., et al. (2013). Inversion of enantioselectivity of a mononuclear non-heme iron(II)-dependent hydroxylase by tuning the interplay of metal-center geometry and protein structure. Angew. Chem. Int. Ed. 52, 9677–9681. doi: 10.1002/anie.201304633
Pratter, S. M., Light, K. M., Solomon, E. I., and Straganz, G. D. (2014b). The role of chloride in the mechanism of O2 activation at the mononuclear nonheme Fe(II) center of the halogenase HctB. J. Am. Chem. Soc. 136, 9385–9395. doi: 10.1021/ja503179m
Puri, M., Biswas, A. N., Fan, R., Guo, Y., and Que Jr, L. (2016). Modeling non-heme iron halogenases: high-spin oxoiron(IV)–halide complexes that halogenate C–H Bonds. J. Am. Chem. Soc. 138, 2484–2487. doi: 10.1021/jacs.5b11511
Quesne, M. G., Borowski, T., and de Visser, S. P. (2016a). Quantum mechanics/molecular mechanics modeling of enzymatic processes: caveats and breakthroughs. Chem. Eur. J. 22, 2562–2581. doi: 10.1002/chem.201503802
Quesne, M. G., and de Visser, S. P. (2012). Regioselectivity of substrate hydroxylation versus halogenation by a nonheme iron(IV)-oxo complex: possibility of rearrangement pathways. J. Biol. Inorg. Chem. 17, 841–852. doi: 10.1007/s00775-012-0901-4
Quesne, M. G., Latifi, R., Gonzalez-Ovalle, L. E., Kumar, D., and de Visser, S. P. (2014). Quantum mechanics/molecular mechanics study on the oxygen binding and substrate hydroxylation step in AlkB repair enzymes. Chem. Eur. J. 20, 435–446. doi: 10.1002/chem.201303282
Quesne, M. G., Senthilnathan, D., Singh, D., Kumar, D., Maldivi, P., Sorokin, A. B., et al. (2016b). Origin of the enhanced reactivity of μ-nitrido-bridged diiron(IV)-oxo porphyrinoid complexes over cytochrome P450 compound, I. ACS Catal. 6, 2230–2243. doi: 10.1021/acscatal.5b02720
Ramaswamy, A. V., Sorrels, C. M., and Gerwick, W. H. (2007). Cloning and biochemical characterization of the hectochlorin biosynthetic gene cluster from the marine cyanobacterium lyngbya majuscule. J. Nat. Prod. 70, 1977–1986. doi: 10.1021/np0704250
Sainna, M. A., Kumar, S., Kumar, D., Fornarini, S., Crestoni, M. E., and de Visser, S. P. (2015). A comprehensive test set of epoxidation rate constants for iron(IV)–oxo porphyrin cation radical complexes. Chem. Sci. 1516–1529. doi: 10.1039/C4SC02717E
Schofield, C. J., and Zhang, Z. (1999). Structural and mechanistic studies on 2-oxoglutarate-dependent oxygenases and related enzymes. Curr. Opin. Struc. Biol. 9, 722–731. doi: 10.1016/S0959-440X(99)00036-6
Shaik, S., Kumar, D., de Visser, S. P., Altun, A., and Thiel, W. (2005). Theoretical perspective on the structure and mechanism of cytochrome P450 enzymes. Chem. Rev. 105, 2279–2328. doi: 10.1021/cr030722j
Sherwood, P., de Vries, A. H., Guest, M. F., Schreckenbach, G., Catlow, C. R. A., French, S. A., et al. (2003). QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis. J. Mol. Struct. 632, 1–28. doi: 10.1016/S0166-1280(03)00285-9
Sinnecker, S., Svensen, N., Barr, E. W., Ye, S., Bollinger Jr, J. M., Neese, F., et al. (2007). Spectroscopic and Computational evaluation of the structure of the high-Spin Fe(IV)-Oxo intermediates in taurine: α-ketoglutarate dioxygenase from Escherichia coli and its His99Ala ligand variant. J. Am. Chem. Soc. 129, 6168–6179. doi: 10.1021/ja067899q
Solomon, E. I., Brunold, T. C., Davis, M. I., Kemsley, J. N., Lee, K., Lehnert, K., et al. (2000). Geometric and electronic structure/function correlations in non-heme iron enzymes. Chem. Rev. 100, 235–349. doi: 10.1021/cr9900275
Srnec, M., and Solomon, E. I. (2017). Frontier molecular orbital contributions to chlorination versus hydroxylation selectivity in the non-heme iron halogenase SyrB2. J. Am. Chem. Soc. 139, 2396–2407. doi: 10.1021/jacs.6b11995
Srnec, M., Wong, S. D., Matthews, M. L., Krebs, C., Bollinger, J. M. Jr., and Solomon, E. I. (2016). Electronic structure of the ferryl intermediate in the α-ketoglutarate dependent non-heme iron halogenase SyrB2: contributions to H atom abstraction reactivity. J. Am. Chem. Soc. 138, 5110–5122. doi: 10.1021/jacs.6b01151
Sundaramoorthy, M., Terner, J., and Poulos, T. L. (1995). The crystal structure of chloroperoxidase: a heme peroxidase–cytochrome P450 functional hybrid. Structure. 3, 1367–1377. doi: 10.1016/S.0969-2126(01)00274-X.
Tang, H., Guan, J., Liu, H., and Huang, X. (2013). Analysis of an alternative to the H-atom abstraction mechanism in methane C–H bond activation by nonheme iron(IV)-oxo oxidants. Dalton Trans. 42, 10260–10270. doi: 10.1039/c3dt50866h
Tchesnokov, E. P., Faponle, A. S., Davies, C. G., Quesne, M. G., Turner, R., Fellner, M., et al. (2016). An iron–oxygen intermediate formed during the catalytic cycle of cysteine dioxygenase. Chem. Commun. 52, 8814–8817. doi: 10.1039/C6CC03904A
Timmins, A., and de Visser, S. P. (2015). Enzymatic halogenases and haloperoxidases: computational studies on mechanism and function. Adv. Protein Chem. Struct. Biol. 100, 113–151. doi: 10.1016/bs.apcsb.2015.06.001
Timmins, A., and de Visser, S. P. (2017). How are substrate binding and catalysis affected by mutating Glu127 and Arg161 in Prolyl-4-hydroxylase? A QM/MM and MD Study. Front.Chem. 5:94. doi: 10.3389/fchem.2017.00094
Timmins, A., Quesne, M. G., Borowski, T., and de Visser, S. P. (2018). Group transfer to an aliphatic bond: a biomimetic study inspired by nonheme iron halogenases. ACS Catal. 8, 8685–8698. doi: 10.1021/acscatal.8b01673
Timmins, A., Saint-André, M., and de Visser, S. P. (2017). Understanding how prolyl-4-hydroxylase structure steers a ferryl oxidant toward scission of a Strong C–H bond. J. Am. Chem. Soc. 139, 9855–9866. doi: 10.1021/jacs.7b02839
Usharani, D., Lacy, D. C., Borovik, A. S., and Shaik, S. (2013). Dichotomous hydrogen atom transfer vs. proton-coupled electron transfer during activation of X–H bonds (X = C, N, O) by nonheme iron–oxo complexes of variable basicity. J. Am. Chem. Soc. 135, 17090–17104. doi: 10.1021/ja408073m
Vaillancourt, F. H., Yeh, E., Vosburg, D. A., Garneaur-Tsodikova, S., and Walsh, C. T. (2006). nature's inventory of halogenation catalysts: oxidative strategies predominate. Chem. Rev. 106, 3364–3378. doi: 10.1021/cr050313i,
van Pée, H., Dong, C., Flecks, S., Naismith, J., Patallo, E. P., and Wage, T. (2006). Biological halogenation has moved far beyond haloperoxidases. Adv. Appl. Microbiol. 59, 127–157. doi: 10.1016/S0065-2164(06)59005-7
Vardhaman, A. K., Barman, P., Kumar, S., Sastri, C. V., Kumar, D., and de Visser, S. P. (2013). Comparison of the reactivity of nonheme iron(IV)-oxo versus iron(IV)-imido complexes: which is the better oxidant? Angew. Chem. Int. Ed. 52, 12288–12292. doi: 10.1002/anie.201305370
Vardhaman, A. K., Sastri, C. V., Kumar, D., and de Visser, S. P. (2011). Nonheme ferric hydroperoxo intermediates are efficient oxidants of bromide oxidation. Chem. Commun. 47, 11044–11046. doi: 10.1039/c1cc13775a
Wang, Y., Li, X., Yang, G., He, G., and Chen, G. (2016). A visible-light-promoted radical reaction system for azidation and halogenation of tertiary aliphatic C–H bonds. Chem. Sci. 7, 2679–2683. doi: 10.1039/C5SC04169D
Weichold, V., Milbredt, D., and van Pée, K. H. (2016). Specific enzymatic halogenation—from the discovery of halogenated enzymes to their applications in vitro and in vivo. Angew. Chem. Int. Ed. 55, 6374–6389. doi: 10.1002/anie.201509573
Wójcik, A., Radon, M., and Borowski, T. (2016). Mechanism of O2 activation by α-ketoglutarate dependent oxygenases revisited. A Quantum Chemical Study. J. Phys. Chem. A 120, 1261–1274. doi: 10.1021/acs.jpca.5b12311
Wong, S. D., Srnec, M., Matthews, M. L., Liu, L. V., Kwak, Y., Park, K., Solomon, E. I., et al. (2013). Elucidation of the iron(IV)–oxo intermediate in the non-haem iron halogenase SyrB2. Nature 499, 320–323. doi: 10.1038/nature12304
Yang, T., Quesne, M. G., Neu, H. M., Cantú Reinhard, F. G., Goldberg, D. P., and de Visser, S. P. (2016). Singlet versus triplet reactivity in an mn(v)–oxo species: testing theoretical predictions against experimental evidence. J. Am. Chem. Soc. 138, 12375–12386. doi: 10.1021/jacs.6b05027
Ye, S., and Neese, F. (2011). Nonheme oxo-iron(IV) intermediates form an oxyl radical upon approaching the C–H bond activation transition state. Proc. Natl. Acad. Sci. U.S.A. 108, 1228–1233. doi: 10.1073/pnas.1008411108
Keywords: nonheme iron, enzyme catalysis, reaction mechanism, QM/MM, density functional theory, halogenation, hydroxylation
Citation: Timmins A, Fowler NJ, Warwicker J, Straganz GD and de Visser SP (2018) Does Substrate Positioning Affect the Selectivity and Reactivity in the Hectochlorin Biosynthesis Halogenase? Front. Chem. 6:513. doi: 10.3389/fchem.2018.00513
Received: 01 September 2018; Accepted: 04 October 2018;
Published: 30 October 2018.
Edited by:Fahmi Himo, Stockholm University, Sweden
Reviewed by:Rongzhen Liao, Huazhong University of Science and Technology, China
Robert S. Paton, Colorado State University, United States
Copyright © 2018 Timmins, Fowler, Warwicker, Straganz and de Visser. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sam P. de Visser, email@example.com