19F Electron-Nuclear Double Resonance Reveals Interaction between Redox-Active Tyrosines across the α/β Interface of E. coli Ribonucleotide Reductase

Exploring protein structural ensembles: Integration of sparse experimental data from electron paramagnetic resonance spectroscopy with molecular modeling methods

Under physiological conditions, proteins continuously undergo structural fluctuations on different timescales. Some conformations are only sparsely populated, but still play a key role in protein function. Thus, meaningful structure-function frameworks must include structural ensembles rather than only the most populated protein conformations. To detail protein plasticity, modern structural biology combines complementary experimental and computational approaches. In this review, we survey available computational approaches that integrate sparse experimental data from electron paramagnetic resonance spectroscopy with molecular modeling techniques to derive all-atom structural models of rare protein conformations. We also propose strategies to increase the reliability and improve efficiency using deep learning approaches, thus advancing the field of integrative structural biology.

Conformational control over proton-coupled electron transfer in metalloenzymes

From the reduction of dinitrogen to the oxidation of water, the chemical transformations catalysed by metalloenzymes underlie global geochemical and biochemical cycles. These reactions represent some of the most kinetically and thermodynamically challenging processes known and require the complex choreography of the fundamental building blocks of nature, electrons and protons, to be carried out with utmost precision and accuracy. The rate-determining step of catalysis in many metalloenzymes consists of a protein structural rearrangement, suggesting that nature has evolved to leverage macroscopic changes in protein molecular structure to control subatomic changes in metallocofactor electronic structure. The proton-coupled electron transfer mechanisms operative in nitrogenase, photosystem II and ribonucleotide reductase exemplify this interplay between molecular and electronic structural control. We present the culmination of decades of study on each of these systems and clarify what is known regarding the interplay between structural changes and functional outcomes in these metalloenzyme linchpins.

Extending the Range of Distances Accessible by 19F Electron-Nuclear Double Resonance in Proteins Using High-Spin Gd(III) Labels

Fluorine electron-nuclear double resonance (19F ENDOR) has recently emerged as a valuable tool in structural biology for distance determination between F atoms and a paramagnetic center, either intrinsic or conjugated to a biomolecule via spin labeling. Such measurements allow access to distances too short to be measured by double electron-electron resonance (DEER). To further extend the accessible distance range, we exploit the high-spin properties of Gd(III) and focus on transitions other than the central transition (|-1/2⟩ ↔ |+1/2⟩), that become more populated at high magnetic fields and low temperatures. This increases the spectral resolution up to ca. 7 times, thus raising the long-distance limit of 19F ENDOR almost 2-fold. We first demonstrate this on a model fluorine-containing Gd(III) complex with a well-resolved 19F spectrum in conventional central transition measurements and show quantitative agreement between the experimental spectra and theoretical predictions. We then validate our approach on two proteins labeled with 19F and Gd(III), in which the Gd-F distance is too long to produce a well-resolved 19F ENDOR doublet when measured at the central transition. By focusing on the |-5/2⟩ ↔ |-3/2⟩ and |-7/2⟩ ↔ |-5/2⟩ EPR transitions, a resolution enhancement of 4.5- and 7-fold was obtained, respectively. We also present data analysis strategies to handle contributions of different electron spin manifolds to the ENDOR spectrum. Our new extended 19F ENDOR approach may be applicable to Gd-F distances as large as 20 Å, widening the current ENDOR distance window.

Structure and Flexibility of Copper-Modified DNA G-Quadruplexes Investigated by 19 F ENDOR Experiments at 34 GHz

DNA G-quadruplexes (GQs) are of great interest due to their involvement in crucial biological processes such as telomerase maintenance and gene expression. Furthermore, they are reported as catalytically active DNAzymes and building blocks in bio-nanotechnology. GQs exhibit remarkable structural diversity and conformational heterogeneity, necessitating precise and reliable tools to unravel their structure-function relationships. Here, we present insights into the structure and conformational flexibility of a unimolecular GQ with high spatial resolution via electron-nuclear double resonance (ENDOR) experiments combined with Cu(II) and fluorine labeling. These findings showcase the successful application of the 19 F-ENDOR methodology at 34 GHz, overcoming the limitations posed by the complexity and scarcity of higher-frequency spectrometers. Importantly, our approach retains both sensitivity and orientational resolution. This integrated study not only enhances our understanding of GQs but also expands the methodological toolbox for studying other macromolecules.

19F electron nuclear double resonance (ENDOR) spectroscopy for distance measurements using trityl spin labels in DNA duplexes

The combination of fluorine labeling and pulsed electron-nuclear double resonance (ENDOR) is emerging as a powerful technique for obtaining structural information about proteins and nucleic acids. In this work, we explored the capability of Mims 19F ENDOR experiments on reporting intermolecular distances in trityl- and 19F-labeled DNA duplexes at three electron paramagnetic resonance (EPR) frequencies (34, 94, and 263 GHz). For spin labeling, we used the hydrophobic Finland trityl radical and hydrophilic OX063 trityl radical. Fluorine labels were introduced into two positions of a DNA oligonucleotide. The results indicated that hyperfine splittings visible in the ENDOR spectra are consistent with the most populated interspin distances between 19F and the trityl radical predicted from molecular dynamic (MD) simulations. Moreover, for some cases, ENDOR spectral simulations based on MD results were able to reproduce the conformational distribution reflected in the experimental ENDOR line broadening. Additionally, MD simulations provided more detailed information about the melting of terminal base pairs of the oligonucleotides and about the configuration of the trityls relative to a DNA end.

Structural insights into the half-of-sites reactivity in homodimeric and homotetrameric metalloenzymes

Half-of-sites reactivity in many homodimeric and homotetrameric metalloenzymes has been known for half a century, yet its benefit remains poorly understood. A recently reported cryo-electron microscopy structure has given some clues on the less optimized reactivity of Escherichia coli ribonucleotide reductase with an asymmetric association of α2β2 subunits during catalysis. Moreover, nonequivalence of enzyme active sites has been reported in many other enzymes, possibly as a means of regulation. They are often induced by substrate binding or caused by a critical component introduced from a neighboring subunit in response to substrate loadings, such as in prostaglandin endoperoxide H synthase, cytidine triphosphate synthase, glyoxalase, tryptophan dioxygenase, and several decarboxylases or dehydrogenases. Overall, half-of-sites reactivity is likely not an act of wasting resources but rather a method devised in nature to accommodate catalytic or regulatory needs.

Bayesian optimization to estimate hyperfine couplings from 19F ENDOR spectra

ENDOR spectroscopy is a fundamental method to detect nuclear spins in the vicinity of paramagnetic centers and their mutual hyperfine interaction. Recently, site-selective introduction of 19F as nuclear labels has been proposed as a tool for ENDOR-based distance determination in biomolecules, complementing pulsed dipolar spectroscopy in the range of angstrom to nanometer. Nevertheless, one main challenge of ENDOR still consists of its spectral analysis, which is aggravated by a large parameter space and broad resonances from hyperfine interactions. Additionally, at high EPR frequencies and fields (⩾94 GHz/3.4 Tesla), chemical shift anisotropy might contribute to broadening and asymmetry in the spectra. Here, we use two nitroxide-fluorine model systems to examine a statistical approach to finding the best parameter fit to experimental 263 GHz 19F ENDOR spectra. We propose Bayesian optimization for a rapid, global parameter search with little prior knowledge, followed by a refinement by more standard gradient-based fitting procedures. Indeed, the latter suffer from finding local rather than global minima of a suitably defined loss function. Using a new and accelerated simulation procedure, results for the semi-rigid nitroxide-fluorine two and three spin systems lead to physically reasonable solutions, if minima of similar loss can be distinguished by DFT predictions. The approach also delivers the stochastic error of the obtained parameter estimates. Future developments and perspectives are discussed.

Radical Transport Facilitated by a Proton Transfer Network at the Subunit Interface of Ribonucleotide Reductase

Ribonucleotide reductases (RNRs) play an essential role in the conversion of nucleotides to deoxynucleotides in all organisms. The Escherichia coli class Ia RNR requires two homodimeric subunits, α and β. The active form is an asymmetric αα'ββ' complex. The α subunit houses the site for nucleotide reduction initiated by a thiyl radical (C439•), and the β subunit houses the diferric-tyrosyl radical (Y122•) that is essential for C439• formation. The reactions require a highly regulated and reversible long-range proton-coupled electron transfer pathway involving Y122•[β] ↔ W48?[β] ↔ Y356[β] ↔ Y731[α] ↔ Y730[α] ↔ C439[α]. In a recent cryo-EM structure, Y356[β] was revealed for the first time and it, along with Y731[α], spans the asymmetric α/β interface. An E52[β] residue, which is essential for Y356 oxidation, allows access to the interface and resides at the head of a polar region comprising R331[α], E326[α], and E326[α'] residues. Mutagenesis studies with canonical and unnatural amino acid substitutions now suggest that these ionizable residues are important in enzyme activity. To gain further insights into the roles of these residues, Y356• was photochemically generated using a photosensitizer covalently attached adjacent to Y356[β]. Mutagenesis studies, transient absorption spectroscopy, and photochemical assays monitoring deoxynucleotide formation collectively indicate that the E52[β], R331[α], E326[α], and E326[α'] network plays the essential role of shuttling protons associated with Y356 oxidation from the interface to bulk solvent.

Direct Proton-Coupled Electron Transfer between Interfacial Tyrosines in Ribonucleotide Reductase

Ribonucleotide reductase (RNR) regulates DNA synthesis and repair in all organisms. The mechanism of Escherichia coli RNR requires radical transfer over a proton-coupled electron transfer (PCET) pathway spanning ∼32 Å across two protein subunits. A key step along this pathway is the interfacial PCET reaction between Y356 in the β subunit and Y731 in the α subunit. Herein, this PCET reaction between two tyrosines across an aqueous interface is explored with classical molecular dynamics and quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The simulations suggest that the water-mediated mechanism involving double proton transfer through an intervening water molecule is thermodynamically and kinetically unfavorable. The direct PCET mechanism between Y356 and Y731 becomes feasible when Y731 is flipped toward the interface and is predicted to be approximately isoergic with a relatively low free energy barrier. This direct mechanism is facilitated by the hydrogen bonding of water to both Y356 and Y731. These simulations provide fundamental insights into radical transfer across aqueous interfaces.

Selective incorporation of 5-hydroxytryptophan blocks long range electron transfer in oxalate decarboxylase

Oxalate decarboxylase from Bacillus subtilis is a binuclear Mn-dependent acid stress response enzyme that converts the mono-anion of oxalic acid into formate and carbon dioxide in a redox neutral unimolecular disproportionation reaction. A π-stacked tryptophan dimer, W96 and W274, at the interface between two monomer subunits facilitates long-range electron transfer between the two Mn ions and plays an important role in the catalytic mechanism. Substitution of W96 with the unnatural amino acid 5-hydroxytryptophan leads to a persistent EPR signal which can be traced back to the neutral radical of 5-hydroxytryptophan with its hydroxyl proton removed. 5-Hydroxytryptophan acts as a hole sink preventing the formation of Mn(III) at the N-terminal active site and strongly suppresses enzymatic activity. The lower boundary of the standard reduction potential for the active site Mn(II)/Mn(III) couple can therefore be estimated as 740 mV against the normal hydrogen electrode at pH 4, the pH of maximum catalytic efficiency. Our results support the catalytic importance of long-range electron transfer in oxalate decarboxylase while at the same time highlighting the utility of unnatural amino acid incorporation and specifically the use of 5-hydroxytryptophan as an energetic sink for hole hopping to probe electron transfer in redox proteins.