Unexpected electron transfer mechanism upon AdoMet cleavage in radical SAM proteins (original) (raw)

Abstract

Radical _S_-adenosine-L-methionine (SAM or AdoMet) proteins are involved in chemically difficult reactions including the synthesis of cofactors, the generation of protein radicals, and the maturation of complex organometallic catalytic sites. In the first and common step of the reaction, a conserved [Fe4S4] cluster donates an electron to perform the reductive cleavage of AdoMet into methionine and a reactive radical 5′-dA· species. The latter extracts a hydrogen atom from substrate eliciting one of the about 40 reactions so far characterized for this family of proteins. It has been suggested that the radical-generating mechanism differs depending on whether AdoMet is a cofactor or a substrate. It has also been speculated that electron transfer from the [Fe4S4] cluster to AdoMet is sulfur-based. Here we have used protein crystallography and theoretical calculations to show that regardless whether AdoMet serves as a cofactor or a substrate, the 5′-dA· generating mechanism should be common to the radical SAM proteins studied so far, and that electron transfer is mediated by a unique Fe from the conserved [Fe4S4] cluster. This unusual electron transfer is determined by the sulfonium ion in AdoMet.

Keywords: density functional theory, hybrid potentials, iron sulfur cluster, redox chemistry, _S_-adenosyl-L-methionine


The catalytic cycle of all of the radical _S_-adenosyl-L-methionine (SAM or AdoMet) enzymes (13) involves the formation of a 5′-deoxyadenosyl radical species (5′-dA·), resulting from the reductive cleavage of AdoMet. Diverse reactions are catalyzed by this family of proteins, such as sulfur insertion, isomerization and protein radical formation; all initiated by the abstraction of a hydrogen atom from the substrate by the highly reactive 5′-dA·. The electron donor to form the radical species is a reduced [Fe4S4]+ cluster, coordinated by a generally conserved Cys-X3-Cys-X2-Cys motif (4). Model chemistry has shown that reductive cleavage of sulfonium-bearing compounds by reduced [Fe4S4] clusters is a facile reaction, although it is generally a two-electron process (5). In addition, electrophilic attack of sulfonium to thiolates has been also noted in these model systems but not in the enzymes. These observations underscore the role of the protein in modulating the one-electron reaction that leads to 5′-dA· formation. A 2H and 13C ENDOR study of pyruvate formate-lyase activating enzyme (PFL-AE) predicted the distances between the closest iron atom of the [Fe4S4] cluster and the methyl protons and carbon atom of AdoMet to be approximately 3.0 Å to 3.8 Å and approximately 4–5 Å, respectively (6). Using these data, a model for the interaction between AdoMet and the cubane was postulated where the carboxylate moiety of AdoMet established a bidentate interaction with an iron ion from the cluster. The ENDOR-based model was subsequently refined using additional isotopic labeling in the same enzyme to include an interaction of the carboxy and amino groups of AdoMet (7) with a unique iron site, identified by Mössbauer spectroscopy (8). An additional interaction between one of the sulfides of the cluster and the sulfonium ion (Sδ+) of AdoMet was also postulated and it was argued that this interaction could be involved in sulfur-based inner sphere electron transfer from the cluster to AdoMet (6, 7). The X-ray structures of the AdoMet-bound Radical SAM enzymes HemN (9), BioB (10), MoaA (11), and LAM (12) have confirmed the prediction that the [Fe4S4] cluster interacts with the carboxy and amino groups of AdoMet. However, as it will be discussed below, no close interaction between S (S3* in Fig. 1) and Sδ+ (6, 7) has been observed in these structures.

Fig. 1.

Fig. 1.

Depiction of AdoMet radical cleavage. Fe4* indicates the unique iron site; S3* and Sδ+/Sδ are discussed in the text ([Fe4S4] cluster numbering as in PDB code 3CIW).

Selenium K-edge X-ray absorption spectroscopy (XAS) experiments carried out with lysine 2,3-aminomutase (LAM) revealed a Se-Fe interaction at 2.7 Å after cleavage of _S_-adenosyl-L-selenomethionine (13). The absence of this interaction in PFL-AE and BioB (14), was thought to reflect the role of AdoMet during the reaction. AdoMet is a substrate for both PFL-AE that generates a glycyl radical in PFL (15) and for BioB that inserts a sulfur atom in dethiobiotin (16). Conversely, AdoMet is a cofactor for LAM that catalyzes the interconversion of α-L-lysine and β-L-lysine (17) and uses it as a reversible source of 5′-dA·. A structural basis for the XAS-based difference in selenomethionine binding after AdoMet cleavage was proposed using the ENDOR model described above: after cleavage, the methionine Sδ would interact with one iron ion from the cluster in LAM (13) and it would be either more distant or dissociated from the [Fe4S4] cluster in PFL-AE and BioB (14). In all these structures, and in the recently published structures of the [FeFe]-hydrogenase active site maturase HydE from Thermotoga maritima (18) and PFL-AE complexed with a peptide (19), AdoMet binds the conserved [Fe4S4] cluster in the same way. This observation favors methionine dissociation because it is not obvious why, upon AdoMet cleavage, methionine should bind differently in any of these proteins.

The 1.35-Å resolution HydE structure (18) provided a high resolution image of a radical AdoMet enzyme with bound _S_-adenosyl-L-homocysteine (SAH) (in the remainder the AdoMet cleavage reaction site will be called the active site). Following our report on that crystal structure, and to better define the binding of cleavage products in a Radical SAM protein, we have solved the structures of HydE bound to AdoMet and [5′-dA + Met] up to 1.62 Å and 1.25-Å resolution, respectively (Fig. 2). These high resolution structures have subsequently been used to investigate in detail the AdoMet cleavage mechanism by computational methods.

Fig. 2.

Fig. 2.

Active site X-ray structures. 2_Fo-Fc_ electron density maps on the x-ray model of (A) AdoMet-bound HydE and (B) [5′-dA + Met]-bound HydE. Map contours were drawn at the 1 σ level with a cover radius of 1.8 Å. The three cysteine residues, the iron sulfur cluster and AdoMet or [5′-dA + Met] are represented as sticks; C gray, Fe brown, N blue, O red, S yellow.

Results

Crystal Structures of the AdoMet-Bound and [5′-dA + Met]-Bound HydE.

Table 1 summarizes data collection and refinement statistics (see Materials and Methods). Fig. 2 A and B depict the active sites of AdoMet-bound and [5′-dA + Met]-bound HydE, respectively. We observe that, as in all other available AdoMet-bound X-ray structures (Table 2), the distance between the AdoMet sulfonium ion and the unique iron of the cluster is shorter than the S3*–Sδ+ distance. As shown in Fig. 2, there are some differences between the [5′-dA + Met]-HydE complex and its AdoMet counterpart. 5′-dA binds to HydE approximately as the adenosine part of AdoMet does (SI Text). However, free methionine binding differs significantly from its covalently bound counterpart in AdoMet. Whereas the Sδ from methionine is still far from the closest [Fe4S4] sulfur atom S3* at 3.50 Å, its distance to the unique iron, Fe4*, is only 2.67 Å (Fig. 1). When AdoMet-bound and [5′-dA + Met]-bound HydE structures are compared, the distortion of the cluster is more pronounced in the latter, resulting from pseudooctahedral coordination at Fe4* (SI Text). We have already discussed a similar situation in the SAH-bound HydE structure (18). A structure of the active site after AdoMet cleavage has been postulated based on a crystallographic study of MoaA in complex with its substrate 5′GTP, obtained after reduction of the [Fe4S4] clusters (20). However, our examination of the corresponding electron density strongly suggests the presence of uncleaved AdoMet in one of the two MoaA molecules of the asymmetric unit and the absence of AdoMet-related species in the other (SI Text). A very recent report by Challand et al. (21) has shown that in the case of three AdoMet-dependent proteins [BioB, lipoyl synthase (LipA) and tyrosine lyase (ThiH)] that use AdoMet as a substrate, the products of its cleavage, methionine and 5′-dA, are inhibitory in vitro. The inhibition was removed by the enzyme 5′-methylthioadenosine _S_-adenosylhomocysteine nucleosidase that hydrolyses 5′-dA. BioB, LipA, and ThiH are expected to have lower binding constant for methionine and 5′-dA than enzymes that use AdoMet as a cofactor. Indeed, although it is not explicitly stated in their paper, at equimolar concentrations (30 μM) BioB was only about 10% inhibited by products [Fig. 2C in (21)]. This is to be compared with the XAS study of LAM (13) where, at equimolar concentrations, there was a strong Fe-Se peak at d = 2.7 Å resulting from the interaction of Se-Met and the [Fe4S4] cluster. Experiments equivalent to those of Challand et al. (21) have not been reported for LAM or other reversible radical AdoMet enzymes. However, we have observed extremely tight binding of methionine and 5′-dA to HydE in our crystals. Even after a prolonged back-soak of a crystal in its crystallization solution devoid of these molecules, the electron density map was identical to the one shown in Fig. 2B. Although the role of this protein in [FeFe]-hydrogenase maturation is unknown, this observation suggests that HydE uses AdoMet as a cofactor.

Table 1.

Crystallographic data and refinement statistics

Crystal
AdoMet [5′-dA + Met]
Space group _P_212121
Data collection
Beamline ID14-eh2 ID29
Cell parameters
a, Å 50.63 51.05
b, Å 78.78 78.92
c, Å 85.94 86.19
Wavelength, Å 0.93 0.976
Resolution, Å 1.62 1.25
_R_sym 0.066 (0.369) 0.045 (0.407)
I/σ(Ι) 13.39 (2.38) 14.87 (2.10)
Completeness, % 93.5 (72.9) 95.1 (76.8)
Multiplicity 2.5 (1.4) 3.2 (1.45)
Refinement statistics
_R_cryst 0.137 0.139
_R_free 0.184 0.166
No. of reflections
Work set 38,939 83,859
Test set 2,204 4,675
rmsd from ideal geometry
Bonds, Å 0.014 0.012
Angles, ° 1.749 1.675
No. of non-hydrogen atoms 3,472 3,545
Water molecules 446 467
Iron atoms 6 4
Others 193 197

Table 2.

AdoMet binding

Protein Reference HemN (2.07) (ref. 9) BioB (3.40) (ref. 10) LAM (2.10) (ref. 12) MoaA (2.20) (ref. 11) PFL-AE (2.77) (ref. 19) HydE (1.62) (this work)
Sδ-S3* 3.64 4.77 3.79 3.84 3.97 3.66
Sδ-Fe4* 3.35 4.03 3.15 3.19 3.22 3.25
OAdoMet-Fe4* 2.13 2.51 1.98 1.97 2.17 2.25
NAdoMet-Fe4* 2.33 2.35 1.98 2.30 2.12 2.33

Calculations on Enzyme Models.

Hybrid quantum mechanical (QM)/molecular mechanical (MM) potentials have been used to investigate numerous enzymatic reaction mechanisms (22, 23). We have used the QM/MM approach as implemented in the Schrödinger Suite [(a) Maestro, version 8.5; (b) QSite, version 5.0; (c) Jaguar, version 7.5, Schrödinger, LLC, New York, NY, 2008] (see Materials and Methods), along with the crystal structures of AdoMet and [5′-dA + Met] HydE complexes, to model AdoMet cleavage (SI Text, Fig. S1). In Table 3, we report the key distances at the active site for the geometry-optimized structures of a model of the reactant (R) and the radical product (P), as depicted in Fig. 1. Table 3 also shows key distances for the transition state structure (TS). Its energy, relative to the R state, corresponds to the electron transfer activation barrier between R and P. We observed only minor deviations between the QM/MM optimized models of R and P and the starting X-ray models; thus, the N/O anchoring is maintained, the AdoMet binding does not vary and 5′-dA· and 5′-dA adopt similar conformations. The distortion of the [Fe4S4] cluster due to methionine binding to Fe4* is slightly more pronounced in the calculated models than in the X-ray structures (Table 3) while the other S-Fe iron distances in the cluster are standard (d = 2.3 Å). The calculated barrier height for AdoMet cleavage is about 54.0 kJ/mol (Fig. 3), similar to the estimated 37.7 and 54.4 kJ/mol barriers for AdoMet cleavage in LAM, in the presence and absence of substrate, respectively (24). The long Sδ-C5′ distance in the TS state (3.3 Å) indicates that the bond is already broken and the spin density at C5′ confirms that the electron has already been transferred (Fig. 4A). In addition, the TS and P structures are very similar, the only difference residing in the position of Sδ that is further away from the radical C5′ in the latter. This similarity explains the small 8.4 kJ/mol difference between the two states (Fig. 3).

Table 3.

QM/MM models

X-ray X-ray R§ TS P
NAdoMet-Fe4* 2.33 2.26 2.30 2.24 2.25
OAdoMet-Fe4* 2.25 2.27 2.15 2.13 2.14
Sδ-Fe4* 3.25 2.67 3.20 2.62 2.63
Sδ-C5′ 1.82 3.78 1.84 3.33 4.00
Sδ-S3* 3.66 3.50 3.70 3.61 3.58
S3*-Fe4* 2.28 2.35 2.42 2.48 2.47
S1-Fe4* 2.25 2.39 2.58 2.61 2.62
S4-Fe4* 2.33 2.44 2.43 2.60 2.62

Fig. 3.

Fig. 3.

Barrier height for the formation of [5′-dA· + Met] from AdoMet. Energies relative to that of R are in kJ/mol.

Fig. 4.

Fig. 4.

Transition state. (A) The spin density at the active site of the QM/MM-optimized structure of the transition state. (B) HOMO at the active site of the QM/MM-optimized structure of the transition state (surfaces in red [+] and blue [-]; atom color codes are as in Fig. 1).

The highest occupied molecular orbital (HOMO) of the TS gives an indication of how the reducing electron is delocalized when going from R to P. As can be seen in the TS HOMO of Fig. 4B, the main contributions to the orbital come from the radical-bearing C5′ (see Table S2), the intercalating Sδ from methionine and the unique Fe4*. To better understand how the reducing electron would preferentially be transferred during cleavage, we looked at the HOMO and the lowest unoccupied molecular orbital (LUMO) of R. As a first approximation, these delocalized orbitals can be conceptualized as resulting from the interaction between the putative donor (centered on Fe4*) and acceptor (centered on Sδ+ and C5′) localized orbitals (SI Text). However, technically, the perturbation induced by the protein matrix treated with MM prevents us from analyzing the LUMO_s_. Because we wanted to further investigate the role of Fe4* in the electron transfer upon AdoMet cleavage (Fig. 4B), we constructed a small model of the reactant that could be treated quantum mechanically (Fig. S2_A_). We extracted from the QM/MM optimized reactant structure the three cysteines, the iron sulfur cluster and we simplified AdoMet to _S_-methylmethionine (SI Text). Consistent with the composition of the TS HOMO (Fig. 4B), the R HOMO is mainly centered on Fe4*. As expected, both α (spin up) and β ((spin down) LUMOs are centered on the _S_-methylmethionine, with a small Fe4* contribution (Fig. 5). However, Sδ+ does not contribute significantly to the HOMO of this reactant model system. Moreover, we found little spin density on Sδ+ (Table S3), although for electron transfer to occur, non-zero orbital overlap between donor and acceptor is required.

Fig. 5.

Fig. 5.

Active site model. LUMO (A) α (spin up) and (B) β (spin down) for the system composed of [Fe4S4]+ coordinated by three methyl thiolates to model the cysteine residues and _S_-methylmethionine to model AdoMet.

This result can be rationalized by the well known fact that the B3LYP potential used for these QM calculations tends to underestimate the extent of covalency, and, as a consequence, it minimizes orbital overlap, notably in iron sulfur clusters (25). Consequently, we performed additional calculations on the same model with the more delocalized VBP potential implemented in the ADF QM code (26) (see Materials and Methods). This potential is expected to give a more accurate description of the orbitals (SI Text). Using this approach, 17% of the total spin population was found to reside on Sδ+ (Table S4). In addition, the Sδ+ pz character appears in both the HOMO and LUMO, indicating that the electron transfer to C5′ involves a direct path from Fe4* via the Sδ+ orbitals. Finally, given that the X-ray structures of SAH (18) and AdoMet superimpose well (SI Text), we built a second minimal model to evaluate the effect of changing the charge on Sδ from + 1 in the sulfonium of the R model to 0 in the thioether of SAH (Fig. S2_B_). In this model, we found no spin population on Sδ (even when using the VBP potential), indicating that when this atom is neutral it does not contribute to the HOMO/LUMO (Table S3 and Table S4). In both minimal models, the orbitals corresponding to the N/O atoms of methionine lie at a lower energy and, consequently, do not contribute to the HOMO/LUMO involving Fe4* and Sδ.

Discussion

The calculated barrier for AdoMet cleavage of 54.0 kJ/mol agrees well with an estimate of the corresponding barrier in LAM (24). This value should be compared to the reported 133.9 kJ/mol for AdoMet cleavage in solution. The reaction energy profile obtained in this study, with a TS close in energy to the P state, clearly indicates that the reaction is concerted and can easily proceed backwards, regenerating AdoMet. This is in agreement with the fact that, in all radical SAM proteins that have been studied so far, AdoMet is stable and difficult to cleave in the absence of substrate (27).

The model built from distances found between the methyl group of AdoMet and the [Fe4S4] cluster in the 2H and 13C ENDOR study of PFL-AE (6) needs to be revised. Firstly, this model postulated that only one of the iron ions (corresponding to Fe1 in Fig. 1) interacted with the C atom from the methyl group. However, in all available AdoMet-bound radical SAM enzyme crystal structures, both Fe4* and Fe1 are at about 4 Å from the methyl C atom of AdoMet and therefore, both are expected to contribute to the distant dipolar hyperfine tensor of the 13C atom. Secondly, the ENDOR-based model predicted that the redox chemistry of electron transfer from the cluster to the AdoMet was sulfur-centered (6, 28). Our calculations show that it is not the case in HydE and, given the common AdoMet binding, it should not apply to any of the radical SAM proteins of known three-dimensional structure. Indeed, there is no evidence for anything more than a van der Waals contact between S3* and Sδ (Fig. 1) in any of these proteins. The sulfur-centered mechanism (29) has been justified using data published by Noodleman and Case (30). In their study of a [Fe4S4]+ cluster coordinated by four cysteines, the reducing electron charge was found to be mainly located on the sulfide ions, suggesting they play a role in electron transfer from the cluster. Our system differs from Noodleman and Case's in that AdoMet bidentate binding to one of the Fe ions replaces cysteine ligation and places a sulfonium ion close to the cluster (Fig. S3). We found that in the bound AdoMet the Sδ+ orbital is lower in energy than the S2− orbital (SI Text) and the Sδ+ and Fe4* orbitals have matching energies, a necessary condition for electron transfer between two atoms to occur. In our [5′-dA + Met]-bound HydE structure the distance between Sδ and Fe4* is 2.67 Å. This compares well with the 2.7 Å distance between Se and the unique iron found in the LAM XAS study mentioned above (13). In addition, the reported models for AdoMet and cleavage product complexes with LAM (12, 24) are stereochemically correct, as ascertained by our HydE structure. In fact, the two HydE structures reported here satisfy all of the spectroscopic data reported so far for radical SAM proteins.

Taken together, our results indicate that the AdoMet cleavage mechanism should be common to, at least, those radical SAM proteins with known three-dimensional structures (Table 2). We also think that the most parsimonious conclusion would be that the freed methionine binds to the [Fe4S4] cluster, even if only transiently, always as observed in our [5′-dA + Met]-bound HydE complex. Consequently, we postulate that, besides the nature of the reaction, the role of AdoMet as either a cofactor or a substrate is determined by the relative affinity of the enzyme for [5′-dA + Met] and AdoMet and not by the cleavage mechanism itself. We interpret the fact that a Se-Fe interaction was not detected in the XAS study of PFL-AE and BioB using equimolar concentrations of protein and products (14), as an indication of the low affinity of methionine for the AdoMet binding site in these enzymes. The recent report by Challand et al. (21) on product inhibition of irreversible AdoMet-dependent radical proteins gives an idea of this affinity. Total inhibition of BioB with products was observed when about a 10-fold excess of the latter was used. It would be of interest to repeat the XAS studies on BioB (14) using an equivalent excess of products. Knowing the relative affinities of reversible and irreversible AdoMet enzymes for their products would help in the classification of new members of this family.

Our calculations highlight the crucial role of the unique iron site Fe4* from which the electron is transferred to AdoMet. They also underscore the influence of its Sδ+ that specifically changes the electronic distribution and favors iron-centered over the standard sulfur-centered (30) redox chemistry for electron transfer. Structural and functional data from other radical SAM proteins will be necessary to establish whether, as may be expected from our results, the cleavage mechanism postulated here is common to all of the members of this family.

Materials and Methods

For details see SI Text.

Crystal Sample Preparation.

The HydE protein from Thermotoga maritima (_Tm_HydE) protein was purified as previously described (18, 31). Before crystallization either _S_-adenosyl-L-methionine (AdoMet) (2 mM) or 5′-deoxyadenosine (5′-dA) (10 mM) was added to a 10 mg/mL _Tm_HydE protein solution. To obtain the [5′-dA + Methionine]-HydE complex, methionine was added to the pretreated 5′-dA-HydE protein solution in the crystallization drops (final concentration 30 mM). Crystals were grown and flash-cooled as previously described (18).

Data Collection and Model Refinement.

The two X-ray data sets used here were collected at the European Synchrotron Radiation Facility. Data were processed with XDS (32) (see Table 1 for statistics). Both the AdoMet and the [5′-dA + Met] complex structures were solved using our previously published 1.35-Å resolution structure (18) (PDB code 3CIW). Because all of the crystals are isomorphous, there was no need to apply an initial rigid-body refinement procedure. The _S_-adenosyl-L-homocysteine (SAH) and the water molecules filling the barrel cavity were removed from structure factor and phase calculations, before refinement. Subsequent manual building, including fitting of the AdoMet or [5′-dA + Met] into the observed electron density map, followed by refinement with REFMAC (33), led to the final models with statistics presented in Table 1.

Theoretical Calculations.

Calculations on the enzyme models were performed with the QSite code in the Schrödinger Suite [(a) Maestro, version 8.5; (b) QSite, version 5.0; (c) Jaguar, version 7.5, Schrödinger, LLC, New York, NY, 2008]. A hybrid potential combining (i) density functional theory (DFT) with the B3LYP functional to describe the atoms at the active site and (ii) molecular mechanics to include the effect of the protein matrix were used to model the reactant R and the product P of the AdoMet cleavage reaction. We used the LACVP** basis set (34) on metals and 6–31G** on all other atoms for geometry optimization. For single-point energy calculations at geometry-optimized structures, we chose the larger basis set of triple-ζ quality, LACV3P** (Schrödinger Suite) for metals and cc-pVTZ(-f) (35) for all other QM atoms. The rest of the system was treated with the OPLS2005 force field (36).

To comfort the interpretation of our results with the enzyme, we also performed the above computations on small QM models (Fig. S2_A_ and B). Along with using the Gaussian-based Jaguar code in the Schrödinger suite, we ran the Slater-based ADF code (26) (triple-ζ basis sets for all atoms) with the Vosko, Wilk, and Nusair (VWN) LDA functional (37, 38) and Becke-Perdew exchange-correlation potential (39, 40) (combination referred to as VBP). While the mixed Hartree-Fock/VWN B3LYP potential gives good energetics, the standard ‘pure’ DFT VBP potential is better in terms of orbital characterization (25, 41). Reality (so to speak) lies in between.

Supplementary Material

Supporting Information

Acknowledgments.

We thank Ms. Carole Mathevon and Pr. Marc Fontecave for providing purified HydE protein and for fruitful discussions. This work was financed by contract ANR-08-BLAN-0224-01 from the French Agence Nationale pour la Recherche and institutional funds from the Commissariat à l'Energie Atomique et le Centre National de la Recherche Scientifique.

Footnotes

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

Data deposition: The atomic coordinates and structure factors have been deposited in the Protein Data Bank, www.pdb.org (PDB ID codes 3IIZ and 3IIX).

References

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Supporting Information