Detection of Water Molecules on the Radical Transfer Pathway of Ribonucleotide Reductase by 17O Electron-Nuclear Double Resonance Spectroscopy.

The role of water in biological proton-coupled electron transfer (PCET) is emerging as a key for understanding mechanistic details at atomic resolution. Here we demonstrate 17O high-frequency electron-nuclear double resonance (ENDOR) in conjunction with H217O-labeled protein buffer to establish the presence of ordered water molecules at three radical intermediates in an active enzyme complex, the α2β2 E. coli ribonucleotide reductase. Our data give unambiguous evidence that all three, individually trapped, intermediates are hyperfine coupled to one water molecule with Tyr-O···17O distances in the range 2.8-3.1 Å. The availability of this structural information will allow for quantitative models of PCET in this prototype enzyme. The results also provide a spectroscopic signature for water H-bonded to a tyrosyl radical.

Water is no longer known as just the solvent in which biochemical reactions take place but has been recognized as an essential player in these reactions.[1] Of particular interest is water involvement in electron transfer processes,[2−5] its action as a proton wire[6−8] or its role in proton-coupled electron transfer (PCET).[9−12] The identification of internal water in proteins can be achieved by X-ray diffraction.[13−15] However, the crystallization of transient protein complexes is difficult. One key approach for detection of water in biological systems has been the use of 17O-enriched water in conjunction with magnetic resonance spectroscopy.[16−21] Among these methods, electron paramagnetic resonance (EPR) can take advantage of high selectivity, as it detects nuclei only in the ligand sphere (r ≲ 1.5 nm)[22] of paramagnetic centers.

EPR-based 17O hyperfine (hf) spectroscopy has been established for the detection of water binding to transition-metal ions, where the oxygen usually coordinates to the ion and large hyperfine couplings (several MHz) can be observed.[23−25] However, the most common water coordination motif to biological radicals occurs via H-bond interactions. The hf coupling to 17O is diminished in comparison to the metal ion coordination, due to a longer interspin distance. In addition, the small 17O gyromagnetic ratio (γH/γ ≈ 7.4)26 and high nuclear spin (I = 5/2) have rendered the 17O hf coupling difficult to resolve.

Here we illustrate that high-frequency (94 and 263 GHz) electron–nuclear double resonance (ENDOR) spectroscopy can detect the 17O signal of ordered water molecules at an H-bond distance to radical intermediates in E. coli ribonucleotide reductase (RNR). The enzyme uses a long-range (32 Å) radical transfer (RT) to initiate nucleotide reduction (Scheme ).[27] Three tyrosines (Y356, Y731, and Y730) are essential pathway residues, which form transient intermediates in the active complex α2β2, consisting of the two homodimeric subunits α2 and β2.[13,27] Water has been observed only crystallographically in inactive α2s without β2.[13,14,28,29] Using site-selectively inserted tyrosine analogues to trap Y intermediates,[30] our previous 1H/2H ENDOR and DFT studies[10,11,31] revealed H-bonds attributed to water molecules and proposed a role of water in RT. However, all active sites of proteins have exchangeable protons, and thus alternative interpretations to our water proposal were possible. Recently, a cryo-EM structure of α2β2 was reported but the resolution was insufficient for water observation.[32] Since our original proposals, studies using photo-RNRs and MD simulations implying waters in α2β2 have appeared.[33,34] However, water has never been directly detected.

Scheme 1

Current PCET Model of the 32 Å (Y122 to C439) RT in E. coli RNR[27,32],

Redox-active tyrosines 356, 731, and 730 are shown in cyan, electron transfer steps as red arrows, and proton transfer steps as blue arrows. Water molecules revealed in this study in respective site-selective mutants are shown in boldface.

Therefore, we explored the capability of H217O ENDOR spectroscopy by exchanging the RNR buffer with H217O. α2β2-Y356• was generated by a 2,3,5-F3Y122• mutation in β2,[35] whereas radicals at Y731 and Y730 were trapped by replacing the respective residue with 3-aminotyrosine (NH2Y),[36] leading to α2β2-NH2Y731• and α2β2-NH2Y730•. The individual variants were mixed with the complementary α2 or β2 protein, CDP as a substrate, and ATP as an effector. The reaction was then quenched after a few seconds inside EPR tubes. Details on the sample preparation are given in sections SI1 and SI2.

Figure displays representative 94 GHz 17O Mims[37] ENDOR spectra of the radical intermediates.

Figure 1

94 GHz 17O Mims[37] ENDOR spectra of (A) the intermediate Y356• and (B) NH2Y731• and NH2Y730• at B0∥g in the EPR line (T = 50 K, τMims = 390 ns). Acquisition time: 46 h (Y356•), 40 h (NH2Y731•), and 18 h (NH2Y730•). Y356• is from β2-F3Y122•/α2-Y730F, which gives the highest radical yield (section SI2). Experimental spectra are shown in gray, with a Savitzky–Golay filter (fourth-order polynomial, 20-point window) shown in color. Simulations used Easyspin[38] (section SI1.6) with parameters given in Table and section S3. Solid lines (teal) represent transitions among m > 0 manifolds and dashed lines (orange) those among m < 0 manifolds. The simulation does not distinguish between dihedral θ = 0°or θ = 180°.

Table 1

Simulation and DFT Parameters for 17O and 1H hf Couplings of Water in RNR Intermediatesa

	Y356• sim/DFTsmall	NH2Y731• sim	NH2Y730• sim/DFTlargeb
Ax (17O)	0.43/0.19	0.70	0.65/0.24
Ay (17O)	0.66/0.59	0.84	0.80/0.6
Az (17O)	0.70/0.65	0.89	0.89/0.6
A(H1)	6.2[31]/7.4	≲2.5b	2.7b/4.2
ρ(17O)c (%)	0.05		0.03
rO···17O (Å)	2.9 ± 0.1	∼3.0	∼3.0

Except as noted, values are in MHz. Simulated quadrupole values for 17O were {P;P;P} = {−0.02;–0.32;0.34} MHz with eqQ/h = 6.8 MHz and η = 0.93.[41]

Values from 2H couplings in refs (11 and 10) using γ/γ ≈ 6.5.[26]

Loewdin spin density[42] from DFT. Uncertainties in coupling constants are less than 10% for simulations and up to 20% for DFT.

Each spectrum shows a sharp doublet centered on the 17O Larmor frequency (19.3 MHz at 3.4 T), which can be assigned to the central spin transition (m(17O) = +1/2 → −1/2) of one coupled 17O nucleus. As 17O is contained only in the water of the protein buffer, these sharp signals must arise from water molecules coupled to the radicals. Control experiments with only β2 protein confirmed that the signal is associated to the radicals generated in α2β2 (section SI3). The broad resonances at ±2.5 MHz are attributed to other nuclear transitions of the I = 5/2 spin system, broadened by nuclear quadrupole coupling (Figure A). Additionally, we note asymmetry of the doublet, which arises from second-order effects of the quadrupole coupling (section SI4). A comparison of the ENDOR spectra at the low (B0∥g) and high-field (B0∥g) edges of the EPR line (section SI5) indicates an almost isotropic hf coupling, with the dipolar contribution dominating the line width of the central doublet. The 17O ENDOR spectra could be simulated with one 17O nucleus, from which the asymmetry of the central peaks resulted using full diagonalization of the spin Hamiltonian (Figure and section SI8). Parameters are given in Table and section SI3. The spectra of Y356• and NH2Y731• additionally contain signals close to the Larmor frequency not reproduced in the simulations, which likely originate from second-sphere water molecules at the subunit interface. Additional broadening is also observed, particularly at NH2Y731•. It might be caused by conformational distribution of this residue, which was found to have flexibility.[39,40,33]

To rationalize the coupling, we began with a DFT-optimized small model (25 atoms, details in section SI1) of Y356•, as previous ENDOR spectra revealed 1H couplings consistent with one water at the H-bond distance rO–H ≈ 1.8 Å.[31] The 17O coupling from this model was Amax(17O) ≈ 1 MHz, slightly exceeding the present experimental value of 0.6 ± 0.05 MHz. To optimize the model, we computed dihedral θ (C3–C4–O···H) and distance scans for 17O couplings, including the quadrupole tensor and the relative energies (section SI6). The DFT equilibrium distance always resulted in rO–H ≈ 1.8 Å. We found that hf couplings and energies vary significantly with θ, while the quadrupole coupling is less affected (Figure S9A–C). Axyz values of ≲1 MHz are found for θ in the range ≲±30° (or equivalently 150° ≲ θ ≲ 240°): i.e., close to the ring plane. Water coordination in the ring plane also results in minimal relative energies (Figure S9B). Importantly, predicted spin densities on 17O are <0.1% but are sufficient for producing a marked 17O isotropic splitting. The spin density transfer or spin polarization is likely related to the H-bond nature. A distance scan for the optimized dihedral of +20° predicts Amax(17O) in the range 0.75–0.56 MHz (Figure S10A) for rO–H ≈ 1.8–2.0 Å. Consideration of the DFT-predicted 1H couplings (Figure S10B) and comparison with the experimental values[31] of ∼6.2 (H1) and ∼1.6 MHz (H2) indicates that the water is located at rTyr-O··· = 2.9 ± 0.1 Å, corresponding to an rO–H value of 1.9 ± 0.1 Å. Notably, the DFT-predicted dipolar coupling (T|| ≈ 0.3 MHz, section SI6) is consistent with the point-dipole model and the aforementioned broadening of the sharp peaks.

Analogous DFT calculations were performed on the isolated amino tyrosyl NH2Y•.[10,11,36] We observed a trend for the 17O hf coupling in the dihedral and distance scans (section SI7) very similar to the Y• model. The calculation predicts that Aiso(17O) of NH2Y• is slightly larger (10–15%) than that of Y• at similar Tyr-O···17O distances and orientations, which could explain the experimental observation. The amino group introduces an asymmetry in the radical, and the energetically most favored water orientation is found at the opposite side of the amino group (Figure S11B). Nevertheless, this small model could not account simultaneously for the 17O and 1H couplings observed for these two intermediates (Figure S12). As noted in a previous g calculation,[10] the coordination of the water molecule to NH2Y•s is influenced by the surrounding second-sphere residues, as these two intermediates are buried in α2β2 (Scheme ).

Having established that at least one water molecule is hf-coupled to each of the three intermediates, we examine their current molecular models in light of this finding. First, we consider the radical site Y356• (Scheme ). To explain the unprecedented g value of Y356• (g = 2.0062), we previously proposed that two almost equivalent waters might be simultaneously bonded to Y356•.[31] While the present results are most consistent with the distance and orientation proposed for one water, the 94 GHz 17O ENDOR spectra (Figure A) cannot resolve a second water. We note that the spectral line shape and 17O hf coupling in Figure A are conserved in other RNR constructs that generate Y356• (section SI8), including the F3Y122•/E52Q-β2 double mutant used to solve a recent cryo-EM structure.[43]

To gain spectral resolution, we recorded 17O ENDOR spectra of Y356• at 263 GHz/9.4 T (Figure ).[44−46]

Figure 2

Comparison of 94 and 263 GHz Mims ENDOR of Y356• at the three canonical positions in the EPR line. Total acquisition time for 263 GHz (T = 20 K): 18 h (B0∥g), 10 h (B0∥g), and 11 h (B0∥g). Experimental spectra are shown in gray, with a Savitzky–Golay filter (fourth-order polynomial, 10 point window) in color. Simulations of 263 GHz spectra are in black with parameters as for 94 GHz (see Table and Table S4).

The results illustrate that the line width of the central doublet substantially narrows, particularly at B0∥g (Figure ). Despite the narrowing, a factor of approximately 2 from 94 to 263 GHz, we cannot discern two distinct 17O contributions. Simulations of the 263 GHz spectra with the same parameters used at 94 GHz reproduce the line narrowing and support the analysis at 94 GHz. The lack of evidence for a second, almost equivalent water H-bonded to Y356• strongly suggests that the two-water model has become very unlikely and alternative explanations for the shifted g value of Y356• will have to be examined. The precise location of second-sphere residues might play a role,[12] which will require further experimental and computational investigation.

For the radical intermediates in the subunit α, a previous combined ENDOR/DFT model of NH2Y730• proposed a water molecule coordinated in plane at a distance rNH ≈ 3.0 Å.[10] The present results are consistent with this model and provide direct evidence for this postulated water in the enzyme complex α2β2-NH2Y730•. The DFT-predicted hf parameters (DFTlarge) for this large model (140 atoms) are reported in Table , and the model is displayed in section SI9.

Finally, for α2β2-NH2Y731•, large-scale (215 atoms) DFT calculations previously proposed three models of the trapped intermediate (section SI10). Among these models, only one (model 3, Figure S15) contained a water molecule at an H-bond distance. The DFT-predicted 17O hf couplings of model 3 (∼2.5 MHz), however, largely exceed the present experimental values (Table S5). However, this DFT model did not include residues from the β subunit, which we now know are close to this residue in the active complex.[43] Therefore, the model will require further refinement. Nevertheless, the present results give evidence for a water molecule coordinated almost in the plane of NH2Y731•.

In conclusion, we have reported the capability of 17O high-frequency ENDOR to detect water H-bonded to tyrosyl radicals. The spectroscopic approach led to the first detection of ordered water molecules at three trapped radicals proposed to be representative of Y• intermediates in the PCET of E. coli RNR. These results verify previous hypotheses on the presence and role of water in the RNR mechanism and provide a new starting point for computational studies. Knowledge of this 17O signature will also be generally useful for many other biological systems, in which tyrosyl radicals are involved.