Structural and Computational Bases for Dramatic Skeletal

Jul 4, 2018 - AndA, an Fe(II)/α-ketoglutarate (αKG)-dependent enzyme, is the key enzyme that constructs the unique and congested bridged-ring system...
0 downloads 0 Views 4MB Size
Article pubs.acs.org/JACS

Cite This: J. Am. Chem. Soc. 2018, 140, 9743−9750

Structural and Computational Bases for Dramatic Skeletal Rearrangement in Anditomin Biosynthesis Yu Nakashima,†,‡ Takaaki Mitsuhashi,†,‡ Yudai Matsuda,†,§,‡ Miki Senda,|| Hajime Sato,⊥,# Mami Yamazaki,⊥ Masanobu Uchiyama,*,†,# Toshiya Senda,*,||,□ and Ikuro Abe*,†,¶ †

Graduate School of Pharmaceutical Sciences, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan Department of Chemistry, City University of Hong Kong, 83 Tat Chee Avenue, Kowloon, Hong Kong SAR, China || Structural Biology Research Center, Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK), 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan ⊥ Graduate School of Pharmaceutical Science, Chiba University, 1-8-1, Inohana, Chuo-ku, Chiba 260-8675, Japan # Cluster of Pioneering Research (CPR), Advanced Elements Chemistry Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan □ Department of Materials Structure Science, School of High Energy Accelerator Science, The Graduate University for Advanced Studies (Soken-dai), 1−1 Oho, Tsukuba, Ibaraki 305−0801, Japan ¶ Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Yayoi 1-1-1, Bunkyo-ku, Tokyo 113-8657, Japan

Downloaded via UNIV OF WOLLONGONG on August 1, 2018 at 07:56:02 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

§

S Supporting Information *

ABSTRACT: AndA, an Fe(II)/α-ketoglutarate (αKG)-dependent enzyme, is the key enzyme that constructs the unique and congested bridged-ring system of anditomin (1), by catalyzing consecutive dehydrogenation and isomerization reactions. Although we previously characterized AndA to some extent, the means by which the enzyme facilitates this drastic structural reconstruction have remained elusive. In this study, we have solved three X-ray crystal structures of AndA, in its apo form and in the complexes with Fe(II), αKG, and two substrates. The crystal structures and mutational experiments identified several key amino acid residues important for the catalysis and provided insight into how AndA controls the reaction. Furthermore, computational calculations validated the proposed reaction mechanism for the bridged-ring formation and also revealed the requirement of a series of conformational changes during the transformation.



mechanisms for many αKG-dependent enzymes remain poorly understood, as even the mechanism for the dehydrogenation, although it is apparently simple, is still controversial.6 Therefore, it is important to utilize several different approaches to clarify the reactions catalyzed by αKG-dependent enzymes, which will facilitate future bioengineering efforts with this class of enzymes. AndA, involved in the biosynthesis of the fungal meroterpenoid anditomin (1), is an αKG-dependent enzyme that catalyzes two consecutive reactions to construct a unique bridged-ring system (Figure 1). 4c AndA first accepts preandiloid B (2) as a substrate to introduce the C−C double bond between C1 and C2 to yield the enone preandiloid C (3), which is further utilized as the substrate for the second reaction, in which it undergoes an intriguing skeletal rearrangement to produce andiconin (4). Although the first

INTRODUCTION As an excellent designer of small molecules, nature has evolved enzymes that perform diverse chemical transformations, often with the aid of a variety of cofactors and cosubstrates, thus generating the complex and intriguing frameworks of natural products.1 Non-heme iron and α-ketoglutarate (αKG)dependent enzymes are a widespread and major class of enzymes involved in both primary and secondary metabolism.2 They are mostly engaged in simple oxidative reactions, such as hydroxylation and dehydrogenation, but some atypical reactions are catalyzed by this class of enzymes, including the endoperoxidation by FtmOx1 (FtmF) and the epoxidation by AsqJ.3 Moreover, some αKG-dependent enzymes exhibit isomerase or epimerase activity, although they can also work as oxidative enzymes.4 For the catalysis by Fe(II)/αKG-dependent enzymes, there is a strong consensus that the oxidative decarboxylation of αKG first generates a highly active ferryloxo species, which then abstracts a hydrogen atom from a substrate to initiate the reaction.5 Yet, the detailed reaction © 2018 American Chemical Society

Received: June 9, 2018 Published: July 4, 2018 9743

DOI: 10.1021/jacs.8b06084 J. Am. Chem. Soc. 2018, 140, 9743−9750

Article

Journal of the American Chemical Society

meroterpenoid pathways, such as PrhA (PDB ID: 5YBM, RMSD of 0.91 Å for the 202 Cα atoms).8 Next, the binding modes of the metal cofactor, the cosubstrate αKG, and the substrate 2 or 3 were investigated by soaking them in the AndA apo crystal under an anaerobic atmosphere. As a result, we obtained two more crystal structures of AndA, complexed with Fe/αKG/2 (PDB ID: 5ZM3) and Fe/αKG/3 (PDB ID: 5ZM4) at 2.25 and 1.95 Å resolutions, respectively (Figures 2B−E, S1, and S2 and Table S2). In the active sites of both structures, unlike the apo structure (Figure S3), the iron is present at the 2-His-1-Asp facial triad site (His135, Asp137, and His213) conserved among the αKG-dependent dioxygenase family members,9 while the αKG is coordinated to the iron in a bidentate manner, together with the facial triad and a water molecule (w1), to form an octahedral iron complex (Figure 2D and E). The αKG also interacts with the His79, Gln132, Thr173, and Arg224 residues through hydrogen-bonding networks (Figure 2C). The quaternary complexes including substrates (2 or 3) showed that both substrates are bound to AndA in an almost identical fashion (Figures S4 and S5, RMSD of 0.19 Å for the 263 Cα atoms) and that they are located on the opposite side of the facial triad, across from the ferrous ion (Figure 2D and E). As often seen in this class of enzymes, the substrate-binding pocket is constructed by the engagement of both monomers of AndA: the loop between α2 and β3 (loop 58−76, loop A) from one monomer and the loop between α7′ and α8′ (loop 261′−281′, loop B′; hereafter the primes designate residues from the other chain of that without primes) from the other are adjacently located to form a binding site, together with Arg239 on another loop (Figures 2B, S6, and S7). Arg239 also interacts with the carboxyl group of Glu66. Interestingly, loop A was observed only upon substrate binding (Figures 2A,B and S6), indicating that it is highly flexible and that it serves as a lid that encapsulates the substrates in the active site of AndA. This lid-like region appears to recognize the D/E rings of the substrates, as suggested by the hydrogen bond networks observed in this area: two carbonyl oxygens (O3 and O5) both interact with the phenolic hydroxyl group of Tyr272′ as well as the backbones of Gln67 and Ile70, via one water molecule (w3 or w4), while O4 is hydrogen bonded with two water molecules (w2 and w3), which in turn interact with the carboxyl group of Glu66 and the backbones of Lys63 and Gln67 (Figures 2D,E and S7). Meanwhile, the terpenoid moiety (the A/B/C rings) of the substrates is relatively loosely bound to the enzyme, since the interaction between O1 and Asn121 is the only hydrogen bond detected in this region (Figure 2D and E). In line with the common reaction mechanism for αKGdependent dioxygenases, the αKG in the active site should initially undergo oxidative decarboxylation to generate the highly active ferryl-oxo species, in which the oxo group is located trans to His135. It was previously thought that this ferryl-oxo species directly abstracts a hydrogen atom from a substrate, as proposed for the well-studied αKG-dependent enzyme TauD.10 However, recent spectroscopic and computational studies indicated that the isomerization (or rotation) of the oxo group, which proceeds with a low energy barrier, is required before the hydrogen abstraction in the catalysis by several αKG-dependent enzymes.11 Thus, we sought to investigate the possibility that the isomerization event occurs in the AndA-catalyzed reactions. Although it is impossible to

Figure 1. Structure of anditomin (1) and AndA-catalyzed two-step reactions. The genuine “reducing agent” required for the second reaction to provide andiconin (4) has yet to be identified for the in vivo reaction.

reaction is a common dehydrogenation event, the conversion from 3 to 4 is an atypical isomerization. Our previous study revealed that the in vitro enzyme reaction of AndA requires ascorbate as an essential factor, which probably reflects the fact that AndA functions as an isomerase.4c However, the mechanisms for the synthesis of the bridged ring and how the enzyme allows for this drastic structural rearrangement have yet to be elucidated. In this study, to obtain the structural basis for the AndAcatalyzed reactions, we solved the X-ray crystal structures of AndA complexed with iron, αKG, and the individual substrates (2 and 3), which revealed the binding modes of the substrates and suggested the requirement of the ferryl-oxo species isomerization prior to the hydrogen abstraction. Further mutational studies revealed some key residues specifically important for the second reaction to produce 4. We also performed computational calculations of the rearrangement to construct the bridged ring and successfully provided a rationale for this conversion.



RESULTS AND DISCUSSION Crystallographic Studies on AndA. To obtain insight into how AndA participates in the oxidative rearrangement reaction, we sought to crystallize AndA for X-ray diffraction experiments. Since our initial attempt to crystallize the fulllength AndA was not successful, we prepared a truncated protein lacking eight residues at the N-terminus (Met1-Tyr8), which is predicted to be a disordered region. The truncated AndA was treated with EDTA prior to crystallization to completely remove any residual iron, and the metal-free enzyme was then crystallized under anaerobic conditions, leading to the X-ray crystal structure of AndA in its apo form at 2.5 Å resolution (PDB ID: 5ZM2; Figure 2A and Table S2). The crystal structure revealed that AndA exists as a symmetric homodimer, constructed via several hydrogen bond networks at the interface of the two monomeric units (Figure 2A). This observation is consistent with the apparent molecular weight of AndA in its native state predicted by a gel filtration experiment. AndA possesses a double-stranded β-helix fold and forms a funnel-like reaction chamber conserved in the jelly roll barrel, as commonly seen in αKG-dependent dioxygenases (Figure 2A).7 Accordingly, AndA exhibits a highly similar overall structure to other αKG-dependent enzymes involved in fungal 9744

DOI: 10.1021/jacs.8b06084 J. Am. Chem. Soc. 2018, 140, 9743−9750

Article

Journal of the American Chemical Society

Figure 2. (A, B) Overall structures of (A) AndA-apo and (B) the AndA-Fe/αKG/3 complex. (C−E) Close-up view of (C) the αKG-binding site in the AndA-Fe/αKG/3 complex and the substrate-binding sites in the (D) AndA-Fe/αKG/2 and (E) AndA-Fe/αKG/3 complexes. The αKG molecule and the substrates 2 and 3 were modeled into composite omit maps (mFo − DFc map) contoured to 3.5σ. Important lengths are shown as dashed lines with distances labeled (average values for the four subunits in the asymmetric unit, units; Å).

precisely determine the position of the oxo group, it would be close to the water molecule (w1) coordinated trans to His135 and to the carboxylic oxygen of αKG bound trans to His213 before and after the rotation, respectively. Therefore, we measured the distances between these oxygen atoms and the carbon atoms of the substrates from which a hydrogen atom could potentially be abstracted. In the first reaction to yield the enone 3, a hydrogen atom at C1 or C2 of 2 should be abstracted at the beginning of the reaction. In the quaternary

AndA/Fe/αKG/2 complex, the distances from C1 or C2 to w1 are 4.4 and 4.0 Å, respectively, while those to the carboxylic oxygen of αKG are 3.3 and 3.7 Å, respectively (Figure 2D), suggesting that the isomerization could occur to facilitate the hydrogen abstraction. For the second reaction, a hydrogen atom at C12 is initially abstracted,5c and the oxo rotation seems to be more important than that in the first reaction, since the carboxylic oxygen of αKG is located much closer to C12 than w1 (3.5 vs 4.9 Å) in the quaternary AndA/Fe/αKG/ 9745

DOI: 10.1021/jacs.8b06084 J. Am. Chem. Soc. 2018, 140, 9743−9750

Article

Journal of the American Chemical Society

Figure 3. Computed reaction pathways and potential energy profiles of route A and route B. Three-dimensional representations of 7a, 7b, 7c, and 7d are shown. The C8−C2′ length in 7d is represented by dashed lines with the distance labeled (units, Å).

3 complex (Figure 2E). Collectively, our crystallographic study suggests the requirement of the oxo group isomerization, but this should be confirmed by future studies. On the basis of the above discussion, the α-oriented hydrogen at C1 would initially be abstracted in the dehydrogenation reaction to afford 3. To complete the desaturation, the vicinal hydrogen at C2 could be abstracted, followed by a diradical recombination to introduce the C−C double bond (Figure S8A). However, a recent mechanistic study of the AsqJ-catalyzed desaturation indicated alternative mechanisms for this type of reaction: after the generation of the substrate radical, it could then undergo (i) the oxygen rebound and following dehydration or (ii) the oxidization by Fe(III) to a cationic species and subsequent deprotonation, to install the double bond.6b At this moment, none of the three possibilities could be ruled out, but given the structural similarity between AndA and AsqJ, AndA could also use the same strategy for the desaturation as that of AsqJ (Figure S8B). Computational Calculations of the AndA-Catalyzed Rearrangement. We then performed a density functional theory (DFT) calculation (UB3LYP/6-31+G(d,p)) to obtain deeper insight into the reaction pathway and the mechanism for the rearrangement to provide 4. First, we confirmed that the 3D structure and conformation of the substrate 3 used for the calculation are essentially the same as those of 3 in the cocrystal (Table S3 and Figure S9). After the initial hydrogen

abstraction from C12, we identified two possible multistep reaction pathways (routes A and B) for the dramatic skeletal rearrangement from 5 to 8, including (i) C8−O2 bond cleavage, (ii) C12−C5′ bond formation, and (iii) C8−C2′ bond formation, accompanied by several conformational changes (Figures 3 and S10 and Tables S4 and S5). The energy diagram suggested that both routes are thermodynamically and kinetically favorable: (i) the activation barriers are all low enough for the reactions to proceed smoothly at ambient temperature, (ii) the entire energy profile descends as the reactions proceed, and (iii) the overall exergonicity is very large (ΔG = −19.8 kcal/mol). Notably, the geometrical features of the computed 5 are in good agreement with those of the experimentally obtained 3 cocrystallized with AndA. We found that the intermediate 7a and the subsequent cascade reactions to afford 8 had identical geometry and energy in both routes. However, the key intermediate 7a can be reached from the reactant 5 via two different pathways. In route A, the C8− O2 bond cleavage occurs prior to the C12−C5′ bond formation, via the transition state (TS) TS_5−6 to give the exomethylene 6, with a large stabilization energy (ΔG = −19.4 kcal/mol), due to an intrinsically stable tertiary carbon radical. The activation energy for the C−O bond cleavage is only 6.6 kcal/mol. Conversely, route B starts with 5-endo-trig (antiBaldwin) cyclization to produce the α-oxy carbon radical intermediate 6′, with an activation barrier of 23.4 kcal/mol. 9746

DOI: 10.1021/jacs.8b06084 J. Am. Chem. Soc. 2018, 140, 9743−9750

Article

Journal of the American Chemical Society The more than 15 kcal/mol difference in the energies of the first steps between routes A and B translates into an over 1015fold difference in the reaction rates at room temperature; hence, route A should predominate over route B. From the intermediate 6 (route A), an intramolecular macrocyclization can occur smoothly along the intrinsic reaction coordinate with an activation barrier of 13.0 kcal/mol, to give the key intermediate 7a. After several conformational changes (7a → 7b → 7c → 7d) with very small activation barriers, the C8 carbon radical of the B ring in 7d approaches the potent Michael acceptor carbon center C2′ to complete the skeletal rearrangement, generating the most stable radical intermediate (8) with a large stabilization energy (ΔG = −12.3 kcal/mol). In the process of the conformational changes, the radical intermediates (7a, 7b, and 7c) can be interconverted with relatively low barriers. This is consistent with the X-ray crystal structure, as the A/B rings are relatively loosely bound to the enzyme in terms of the hydrogen bonding between the substrates and the enzyme (Figure 2D and E). To further increase the reliability of the above-proposed mechanism, we also investigated the effects of the hydrogen bonds on the radical rearrangement reaction by means of the theozyme (theoretical enzyme) calculation,12 in which some key interactions between the substrate and the enzyme are considered for the calculation. The rearrangement reaction proceeds on the B/C/D rings, and therefore, the hydrogen bonds with the A/E rings are too far from the reaction point to affect the rearrangement reaction. On the other hand, the hydrogen bonding of O-3/O-5 on the D ring with Ile70 and Tyr 272 via a molecule of water could electronically affect the radical rearrangement reaction. Thus, we carried out the DFT calculation with a molecule of water around O-3 and O-5 (Scheme S1), and it turned out that these hydrogen bonds do not greatly change the geometries of the intermediates/ transition states or the activation energies of this rearrangement reaction (Tables S6 and S7 and Figure S11). Thus, it is most likely that the bridged-ring synthesis proceeds through the route A (Figure 3 and Scheme S1). Mutational Experiment and Biochemical Characterization of AndA. To investigate how AndA is engaged in the rearrangement and to examine the importance of the amino acid residues constituting the active site, we performed mutational experiments on selected residues with side chains involved in the substrate binding, namely, Glu66, Asn121, Arg239, and Tyr272′. First, the point mutation of Asn121 (N121A) resulted in reduced activity, but did not completely abolish the conversions (Figure 4A), suggesting that the hydrogen bonding between Asn121 and O1 is not crucial for the reactions to proceed. Among the residues in the lid-like region, the mutations of Glu66 (E66A, E66D, and E66Q) only yielded insoluble proteins, which hampered the evaluation of the importance of this residue in the catalysis. Nevertheless, the mutageneses of Arg239 and Tyr272′ provided intriguing outcomes. The Y272F mutant exhibited activity comparable to that of the wild-type enzyme, but all of the other constructed mutants, R239A, R239M, R239 V, and Y272A, produced only the dehydrogenated product 3 and lost the ability to synthesize 4 when 2 was used as a substrate (Figure 4A). Consistently, 3 was not further transformed into 4 when the above-mentioned four mutants were used for the reaction (Figure 4B). In contrast to the first desaturation reaction, which proceeds relatively inside the substrate-binding pocket, the second reaction occurs near the entrance of the pocket and also

Figure 4. LC-MS profiles of the products from in vitro enzymatic reactions of the AndA mutants with (A) 2 or (B) 3. The chromatograms were extracted at m/z 413.233 (for 2; blue trace) and 411.217 (for 3 and 4; red trace). WT = wild type AndA. NC = negative control.

requires a significant structural rearrangement. Since both Arg239 and Tyr272′ interact with the lid-like region, the abolishment of the second reaction after the introduction of mutations to these residues could be attributed to the destruction of the lid function and the exposure of the D/E rings, which undergo significant structural rearrangements, to the external environment. Collectively, AndA could play a crucial role to protect the highly reactive radical intermediates by housing them in the closed active site, thus allowing the intriguing skeletal rearrangement to proceed without quenching during the transformation. Finally, we performed biochemical characterization of AndA and its mutants. Consistent with our previous in vitro enzymatic reaction with the wild-type AndA,4c the reactions by the mutants that produce 4 required ascorbate as an essential factor (Figure S12), which could be attributed to the prediction that ascorbate is utilized to finalize the reaction. Interestingly, ascorbate is absolutely essential even for the mutants that are only able to perform the first-round reaction 9747

DOI: 10.1021/jacs.8b06084 J. Am. Chem. Soc. 2018, 140, 9743−9750

Journal of the American Chemical Society



to produce 3 (Figure S12), although the dehydrogenation reaction does not apparently require a reducing agent. This observation that ascorbate is indispensable for all the chemistry by AndA indicates that ascorbate has more function(s) than just reducing the product radical. One possible explanation would be that the productions of 3 and 4 are significantly uncoupled from the αKG oxidation and that ascorbate is inevitably required to maintain the oxidation state of iron, but this needs to be further clarified in future studies. We also measured the kinetic parameters for AndA and its mutants (Table 1). The Km values for both 2 and 3 increased after the

substrate

kcat (s−1)

KM (μM)

kcat/KM (s−1 μM−1)

wild type N121A R239A R239V R239M Y272A Y272F wild type N121A R239A R239V R239M Y272A Y272F

2

6.31 ± 0.14 5.00 ± 0.10 2.13 ± 0.11 3.09 ± 0.12 3.05 ± 0.07 4.28 ± 0.09 4.69 ± 0.21 1.88 ± 0.01 3.86 ± 0.24 n.d.a n.d. n.d. n.d. 4.78 ± 0.02

52.6 ± 9.7 84.0 ± 10.7 70.1 ± 28.0 75.8 ± 22.3 85.5 ± 14.0 68.7 ± 10.7 92.3 ± 26.0 53.6 ± 9.3 78.9 ± 20.2 n.d. n.d. n.d. n.d. 116.2 ± 18.8

2.00 1.00 0.48 0.68 0.59 1.04 0.85 0.64 0.82 n.d. n.d. n.d. n.d. 0.69

3

MATERIALS AND METHODS

General Procedures. Solvent and chemicals were purchased from Hampton Research (CA, USA) and Wako Pure Chemical Industries, Ltd. (Tokyo, Japan), unless noted otherwise. Preandiloid B (2) and preandiloid C (3) were prepared according to the previous work.4c Oligonucleotide primers were purchased from Eurofins Genetics (Tokyo, Japan). Sequence analyses were performed by Eurofins Genetics. Enzyme Expression and Purification for Crystallization. The N-terminal truncated andA gene was PCR-amplified from the previously prepared full-length andA-pET28a,4c with the primer pairs andA-9-293-f and andA-9-293-r (Table S1). The DNA fragments were digested with NdeI and BamHI and ligated into pET28a using T4 DNA ligase (Wako Pure Chemical Industries, Ltd.). Constructed plasmids were introduced into Escherichia coli RosettaII(DE3)pLysS (Novagen), and the transformants were cultured with Luria−Bertani medium supplemented with 12.5 μg/mL chloramphenicol and 25 μg/mL kanamycin at 37 °C to an OD600 of 0.6. The cultures were further incubated at 16 °C for 18 h after 0.3 mM isopropyl β-D-1-thiogalactopyranoside addition. The cells were pelleted by a 20 min centrifugation at 6000g and lysed by sonication in buffer A (50 mM Tris-HCl pH 7.5, 200 mM NaCl, 5% (v/v) glycerol, and 5 mM imidazole). After the centrifugation, the supernatant was subjected to a HisTrap HP column (5 mL, GE Healthcare). The column was washed with buffer A containing 15 mM imidazole (30 CV), and then the protein was eluted with buffer A containing 300 mM imidazole (5 CV). The protein solution was loaded on Hi-Load 16/60 Superdex 200 prep grade (GE Healthcare) with gel-filtration buffer (20 mM Tris-HCl pH 7.5, 100 mM NaCl, and 2 mM dithiothreitol) after addition of 1 mM EDTA and rotation for 1 h. The protein was concentrated to 15 mg/mL using a 10 kDa MWCO spin filter (Millipore, USA). The protein concentration was measured by UV absorption at 280 nm (extinction coefficients, 0.678/ M·cm). Protein Crystallization and Structure Determination. The following procedure was undertaken under anaerobic conditions.13 PEG 3350 (30% w/v), 200 mM sodium tartrate, 20% DMF, and 15 mg/mL AndA were used for crystallization, and AndA-apo crystals were obtained with the sitting-drop vapor-diffusion method at 20 °C. Preandiloid B (2) or preandiloid C (3) was introduced into the AndA crystal by soaking in 30% (w/v) PEG3350, 8% tacsimate pH 7.0, 20% DMF, 5 mM FeSO4, 10 mM αKG, and 20 mM 2 or 3 for 6 h at 20 °C. AndA-apo crystal and substrate-soaked crystals were moved into 25% (v/v) glycerol added solution for 10 s and then frozen with liquid nitrogen. The X-ray diffraction data for AndA (BL-1A, Photon Factory, Japan) were processed and scaled using XDS14 and AIMLESS,15 respectively. The structure of PrhA8 (PDB ID: 5YBM) was used as a search model for the molecular replacement procedure (Phaser in PHENIX16). Coot17 and Phenix.Refine18 were utilized for crystallographic refinement. The structural data of 2 and 3 were prepared in Chem3D Ultra software (CambridgeSoft). Substrate occupancy was set at 1.0. The final crystallographic data are listed in Table S2. Computational Details. For all DFT calculations, the Gaussian 16 program19 and/or Reaction Plus programs20 were utilized. Geometry optimizations were conducted at the UB3LYP/6-31+G(d,p) level21 in the gas phase, without any symmetry restrictions. At the same level of theory, vibrational frequency calculations were carried out to confirm that a local minimum has no imaginary frequency and each TS possesses only a single imaginary frequency. Intrinsic reaction coordinate calculations22 for all TSs were conducted with GRRM1123 based on Gaussian 16 (Tables S3 to S7). In this study, the Gibbs free energy was adopted as the basis for discussion. Mutagenesis Analysis. Mutagenesis was performed using PrimeSTAR mutagenesis basal kit (Wako Pure Chemical Industries, Ltd.) using full-length andA in pET28a as a template. Each primer pair (andA-e66a-f and andA-e66a-r, andA-e66d-f and andA-e66d-r, andA-e66q-f and andA-e66q-r, andA-n121a-f and andA-n121a-r, andA-r239a-f and andA-r239a-r, andA-r239m-f and andA-r239m-r,

Table 1. Steady-State Kinetic Parameters for AndA and Its Mutants protein

Article

a

n.d. = not determined.

introduction of the mutations to the residues that are predicted to be involved in the substrate binding, confirming the importance of these residues in the substrate recognition. Additionally, the Km values were almost identical for both substrates 2 and 3, which further supports the crystallographic observation that the two substrates bind to the enzyme in a very similar manner.



CONCLUSION In this study, we have provided new insight into the enzymatic synthesis of the unique bicyclo[2.2.2]octane system of andiconin (4), by crystallographic and computational studies. The crystal structures of AndA revealed the presence of the uncommon lid-like region that locks the substrate in the active site, whose importance for the rearrangement reaction was further validated by the mutational experiments. Furthermore, the computational calculation based on the crystal structure of AndA delineated a plausible reaction mechanism for the bridged-ring formation by AndA. However, the mechanism for the very last step of the reaction, to reduce the product radical 8 into 4, still remains enigmatic. Similar radical quenching to complete the reaction is also required for several other αKGdependent enzymes, such as FtmOx1 (endoperoxidase), CarC (epimerase/dehydrogenase), and SnoN (epimerase), but the mechanisms for the radical reductions by these enzymes are still unclear, despite the availability of their crystal structures.3a,4a,b Thus, further studies are required to identify the reducing agent actually used in the “in vivo” reactions and the key factors that determine whether the enzyme functions as an oxidative enzyme or an isomerase. 9748

DOI: 10.1021/jacs.8b06084 J. Am. Chem. Soc. 2018, 140, 9743−9750

Journal of the American Chemical Society



andA-r239v-f and andA-r239v-r, andA-y272a-f and andA-y272a-r, and andA-y272f-f and andA-y272f-r) was used for constructing E66A, E66D, E66Q, N121A, R239A, R239M, R239 V, Y272A, and Y272F, respectively (Table S1). The plasmids were introduced into E. coli RosettaII(DE3)pLysS after DpnI (Wako Pure Chemical Industries, Ltd.) digestion. Expression and purification conditions were the same as for the wild-type AndA. For the in vitro assay, 5 μM preandiloid B (2) or 5 μM preandiloid C (3) was added to 20 mM Tris-HCl (pH 7.5) buffer containing 200 mM NaCl, 0.1 mM FeSO4, 2.5 mM αKG, 4 mM ascorbate, and 5 μM AndA (wild-type and mutants), in a final volume of 50 μL. After 1 h of enzymatic reaction at 30 °C, 50 μL of methanol was added to each reaction mixture and then analyzed by LC-MS. LC-MS analysis was conducted on a Bruker Compact qTOF mass spectrometer with a Shimadzu Prominence HPLC system, using electrospray ionization with a COSMOSIL 2.5C18-MS-II column (2.0 i.d. × 75 mm; Nacalai Tesque, Inc.), and the separation was performed using the following condition: 70% aqueous acetonitrile containing 20 mM formic acid, 0.1 mL/min. To determine the kinetics parameters of AndA and its mutants with 2 and 3 (5, 10, 50, 100, and 200 μM, duplicate), 50 μL of reaction mixture was incubated with 0.1 μM enzyme for 1 min at 30 °C. The reaction was stopped by adding 50 μL of methanol, and the intensity of consumed substrate was measured by LC-MS. GraphPad Prism 7 for Windows (GraphPad Software, Inc., La Jolla, CA, USA) was used for KM and kcat calculation.



REFERENCES

(1) (a) Tang, M.-C.; Zou, Y.; Watanabe, K.; Walsh, C. T.; Tang, Y. Chem. Rev. 2017, 117, 5226−5333. (b) Richter, M. Nat. Prod. Rep. 2013, 30, 1324−1345. (c) Cox, R. J. Nat. Prod. Rep. 2014, 31, 1405− 1424. (2) (a) Hausinger, R. P. Crit. Rev. Biochem. Mol. Biol. 2004, 39, 21− 68. (b) Loenarz, C.; Schofield, C. J. Nat. Chem. Biol. 2008, 4, 152. (c) Matsuda, Y.; Abe, I. Nat. Prod. Rep. 2016, 33, 26−53. (d) Matsuda, Y.; Awakawa, T.; Mori, T.; Abe, I. Curr. Opin. Chem. Biol. 2016, 31, 1−7. (e) Nakamura, H.; Matsuda, Y.; Abe, I. Nat. Prod. Rep. 2018, DOI: 10.1039/C7NP00055C. (3) (a) Yan, W.; Song, H.; Song, F.; Guo, Y.; Wu, C.-H.; Her, A. S.; Pu, Y.; Wang, S.; Naowarojna, N.; Weitz, A.; Hendrich, M. P.; Costello, C. E.; Zhang, L.; Liu, P.; Zhang, Y. J. Nature 2015, 527, 539−543. (b) Bräuer, A.; Beck, P.; Hintermann, L.; Groll, M. Angew. Chem., Int. Ed. 2016, 55, 422−426. (4) (a) Chang, W.-c.; Guo, Y.; Wang, C.; Butch, S. E.; Rosenzweig, A. C.; Boal, A. K.; Krebs, C.; Bollinger, J. M. Science 2014, 343, 1140− 1144. (b) Siitonen, V.; Selvaraj, B.; Niiranen, L.; Lindqvist, Y.; Schneider, G.; Metsä-Ketelä, M. Proc. Natl. Acad. Sci. U. S. A. 2016, 113, 5251−5256. (c) Matsuda, Y.; Wakimoto, T.; Mori, T.; Awakawa, T.; Abe, I. J. Am. Chem. Soc. 2014, 136, 15326−15336. (5) (a) Price, J. C.; Barr, E. W.; Glass, T. E.; Krebs, C.; Bollinger, J. M. J. Am. Chem. Soc. 2003, 125, 13008−13009. (b) Nam, W.; Lee, Y.M.; Fukuzumi, S. Acc. Chem. Res. 2014, 47, 1146−1154. (6) (a) Martinez, S.; Hausinger, R. P. J. Biol. Chem. 2015, 290, 20702−20711. (b) Liao, H.; Li, J.; Huang, J.; Davidson, M.; Kurnikov, I.; Lin, T.; Lee, J. L.; Kurnikova, M.; Guo, Y.; Chan, N.; Chang, W. Angew. Chem., Int. Ed. 2018, 57, 1831−1835. (7) (a) Clifton, I. J.; McDonough, M. A.; Ehrismann, D.; Kershaw, N. J.; Granatino, N.; Schofield, C. J. J. Inorg. Biochem. 2006, 100, 644−669. (b) Aik, W.; McDonough, M. A.; Thalhammer, A.; Chowdhury, R.; Schofield, C. J. Curr. Opin. Struct. Biol. 2012, 22, 691−700. (8) Nakashima, Y.; Mori, T.; Nakamura, H.; Awakawa, T.; Hoshino, S.; Senda, M.; Senda, T.; Abe, I. Nat. Commun. 2018, 9, 104. (9) Kovaleva, E. G.; Lipscomb, J. D. Nat. Chem. Biol. 2008, 4, 186− 193. (10) Bollinger, J. M.; Price, J. C.; Hoffart, L. M.; Barr, E. W.; Krebs, C. Eur. J. Inorg. Chem. 2005, 2005, 4245−4254. (11) (a) Quesne, M. G.; Latifi, R.; Gonzalez-Ovalle, L. E.; Kumar, D.; de Visser, S. P. Chem. - Eur. J. 2014, 20, 435−446. (b) Huang, J.; Li, C.; Wang, B.; Sharon, D. A.; Wu, W.; Shaik, S. ACS Catal. 2016, 6, 2694−2704. (c) Chang, W.-c.; Li, J.; Lee, J. L.; Cronican, A. A.; Guo, Y. J. Am. Chem. Soc. 2016, 138, 10390−10393. (d) Song, X.; Lu, J.; Lai, W. Phys. Chem. Chem. Phys. 2017, 19, 20188−20197. (12) Tantillo, D. J.; Jiangang, C.; Houk, K. N. Curr. Opin. Chem. Biol. 1998, 2, 743−750. (13) Senda, M.; Kishigami, S.; Kimura, S.; Senda, T. Acta Crystallogr., Sect. F: Struct. Biol. Cryst. Commun. 2007, 63, 311−314. (14) Kabsch, W. Acta Crystallogr., Sect. D: Biol. Crystallogr. 2010, 66, 133−144. (15) Evans, P. R.; Murshudov, G. N. Acta Crystallogr., Sect. D: Biol. Crystallogr. 2013, 69, 1204−1214. (16) McCoy, A. J.; Grosse-Kunstleve, R. W.; Adams, P. D.; Winn, M. D.; Storoni, L. C.; Read, R. J. J. Appl. Crystallogr. 2007, 40, 658−674. (17) Emsley, P.; Cowtan, K. Acta Crystallogr., Sect. D: Biol. Crystallogr. 2004, 60, 2126−2132. (18) Afonine, P. V.; Grosse-Kunstleve, R. W.; Echols, N.; Headd, J. J.; Moriarty, N. W.; Mustyakimov, M.; Terwilliger, T. C.; Urzhumtsev, A.; Zwart, P. H.; Adams, P. D. Acta Crystallogr., Sect. D: Biol. Crystallogr. 2012, 68, 352−367. (19) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A. V.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.;

ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/jacs.8b06084.



Article

Supplementary figures, scheme, and tables (PDF)

AUTHOR INFORMATION

Corresponding Authors

*[email protected] *[email protected] *[email protected] ORCID

Yudai Matsuda: 0000-0001-5650-4732 Masanobu Uchiyama: 0000-0001-6385-5944 Ikuro Abe: 0000-0002-3640-888X Author Contributions ‡

Y. Nakashima, T. Mitsuhashi, and Y. Matsuda contributed equally. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported in part by a Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology, Japan (JSPS KAKENHI Grant Numbers JP15H01836, JP16H06443, JP16H06454, and JP17H05430), JST/NSFC Strategic International Collaborative Research Program and JSPS Research Fellowships for Young Scientists (to Y.N. and T.M.), the Platform for Drug Discovery, Informatics, and Structural Life Science (PDIS), and Basis for Supporting Innovative Drug Discovery and Life Science Research (BINDS) from the Japan Agency for Medical Research and Development (AMED). The synchrotron radiation experiments were performed at the BL-1A of the Photon Factory with proposal Nos. 2014T006 and 2015G530. Allotment of computational resource (Project G18008) from HOKUSAI GreatWave (RIKEN) is gratefully acknowledged. 9749

DOI: 10.1021/jacs.8b06084 J. Am. Chem. Soc. 2018, 140, 9743−9750

Article

Journal of the American Chemical Society Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J. J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Keith, T. A.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian 16, Revision A.03; Gaussian, Inc.: Wallingford, CT, 2016. (20) (a) Reaction Plus, HPC Systems Inc., http://www.hpc.co.jp/ chem/react1.html. (b) Jónsson, H.; Mills, G.; Jacobsen, K. W. Nudged elastic band method for finding minimum energy paths of transitions. In Classical and Quantum Dynamics in Condensed Phase Simulations; World Scientific: Singapore, 1998; pp 385−404. (c) Weinan, E.; Ren, W.; Vanden-Eijnden, E. Phys. Rev. B: Condens. Matter Mater. Phys. 2002, 66, 052301. (21) Becke, A. D. J. Chem. Phys. 1993, 98, 5648−5652. (22) (a) Ishida, K.; Morokuma, K.; Komornicki, A. J. Chem. Phys. 1977, 66, 2153−2156. (b) Fukui, K. Acc. Chem. Res. 1981, 14, 363− 368. (c) Gonzalez, C.; Schlegel, H. B. J. Chem. Phys. 1989, 90, 2154− 2161. (d) Gonzalez, C.; Schlegel, H. B. J. Phys. Chem. 1990, 94, 5523−5527. (e) Page, M.; Doubleday, C.; McIver, J. W., Jr J. Chem. Phys. 1990, 93, 5634−5642. (23) (a) Ohno, K.; Maeda, S. Chem. Phys. Lett. 2004, 384, 277−282. (b) Maeda, S.; Ohno, K. J. Phys. Chem. A 2005, 109, 5742−5753. (c) Ohno, K.; Maeda, S. J. Phys. Chem. A 2006, 110, 8933−8941. (d) Maeda, S.; Osada, Y.; Morokuma, K.; Ohno, K. GRRM 11, Version 11.03; 2012. (e) Maeda, S.; Ohno, K.; Morokuma, K. Phys. Chem. Chem. Phys. 2013, 15, 3683−3701.

9750

DOI: 10.1021/jacs.8b06084 J. Am. Chem. Soc. 2018, 140, 9743−9750