Statistical analysis of ENDOR spectra

Analysis of a Cu-Doped Metal-Organic Framework, MFM-520(Zn1-x Cux ), for NO2 Adsorption

MFM-520(Zn) confines dimers of NO2 with a high adsorption of 4.52 mmol g-1 at 1 bar at 298 K. The synthesis and the incommensurate structure of Cu-doped MFM-520(Zn) are reported. The introduction of paramagnetic Cu2+ sites allows investigation of the electronic and geometric structure of metal site by in situ electron paramagnetic resonance (EPR) spectroscopy upon adsorption of NO2 . By combining continuous wave and electron-nuclear double resonance spectroscopy, an unusual reverse Berry distorted coordination geometry of the Cu2+ centers is observed. Interestingly, Cu-doped MFM-520(Zn0.95 Cu0.05 ) shows enhanced adsorption of NO2 of 5.02 mmol g-1 at 1 bar at 298 K. Whereas MFM-520(Zn) confines adsorbed NO2 as N2 O4 , the presence of monomeric NO2 at low temperature suggests that doping with Cu2+ centers into the framework plays an important role in tuning the dimerization of NO2 molecules in the pore via the formation of specific host-guest interactions.

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.

Frequency swept pulses for the enhanced resolution of ENDOR spectra detecting on higher spin transitions of Gd(III)

Half-Integer High Spin (HIHS) systems with zero-field splitting (ZFS) parameters below 1 GHz are generally dominated by the spin |─1/2>→|+1/2 > central transition (CT). Accordingly, most pulsed Electron Paramagnetic Resonance (EPR) experiments are performed at this position for maximum sensitivity. However, in certain cases it can be desirable to detect higher spin transitions away from the CT in such systems. Here, we describe the use of frequency swept Wideband, Uniform Rate, Smooth Truncation (WURST) pulses for transferring spin population from the CT, and other transitions, of Gd(III) to the neighbouring higher spin transition |─3/2>→|─1/2 > at Q- and W-band frequencies. Specifically, we demonstrate this approach to enhance the sensitivity of ¹H Mims Electron-Nuclear Double Resonance (ENDOR) measurements on two model Gd(III) aryl substituted 1,4,7,10-tetraazacyclododecane-1,4,7-triacetic acid (DO3A) complexes, focusing on transitions other than the CT. We show that an enhancement factor greater than 2 is obtained for both complexes at Q- and W-band frequencies by the application of two polarising pulses prior to the ENDOR sequence. This is in agreement with simulations of the spin dynamics of the system during WURST pulse excitation. The technique demonstrated here should allow more sensitive experiments to be measured away from the CT at higher operating temperatures, and be combined with any relevant pulse sequence.

Gauging the importance of structural parameters for hyperfine coupling constants in organic radicals

The identification of fundamental relationships between atomic configuration and electronic structure typically requires experimental empiricism or systematic theoretical studies. Here, we provide an alternative statistical approach to gauge the importance of structure parameters, i.e., bond lengths, bond angles, and dihedral angles, for hyperfine coupling constants in organic radicals. Hyperfine coupling constants describe electron-nuclear interactions defined by the electronic structure and are experimentally measurable, for example, by electron paramagnetic resonance spectroscopy. Importance quantifiers are computed with the machine learning algorithm neighborhood components analysis using molecular dynamics trajectory snapshots. Atomic-electronic structure relationships are visualized in matrices correlating structure parameters with coupling constants of all magnetic nuclei. Qualitatively, the results reproduce common hyperfine coupling models. Tools to use the presented procedure for other radicals/paramagnetic species or other atomic structure-dependent parameters are provided.

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

Ribonucleotide reductases (RNRs) catalyze the reduction of ribonucleotides to deoxyribonucleotides, thereby playing a key role in DNA replication and repair. Escherichia coli class Ia RNR is an α₂β₂ enzyme complex that uses a reversible multistep radical transfer (RT) over 32 Å across its two subunits, α and β, to initiate, using its metallo-cofactor in β₂, nucleotide reduction in α₂. Each step is proposed to involve a distinct proton-coupled electron-transfer (PCET) process. An unresolved step is the RT involving Y₃₅₆(β) and Y₇₃₁(α) across the α/β interface. Using 2,3,5-F₃Y₁₂₂-β₂ with 3,5-F₂Y₇₃₁-α₂, GDP (substrate) and TTP (allosteric effector), a Y₃₅₆^• intermediate was trapped and its identity was verified by 263 GHz electron paramagnetic resonance (EPR) and 34 GHz pulse electron–electron double resonance spectroscopies. 94 GHz ¹⁹F electron-nuclear double resonance spectroscopy allowed measuring the interspin distances between Y₃₅₆^• and the ¹⁹F nuclei of 3,5-F₂Y₇₃₁ in this RNR mutant. Similar experiments with the double mutant E₅₂Q/F₃Y₁₂₂-β₂ were carried out for comparison to the recently published cryo-EM structure of a holo RNR complex. For both mutant combinations, the distance measurements reveal two conformations of 3,5-F₂Y₇₃₁. Remarkably, one conformation is consistent with 3,5-F₂Y₇₃₁ within the H-bond distance to Y₃₅₆^•, whereas the second one is consistent with the conformation observed in the cryo-EM structure. The observations unexpectedly suggest the possibility of a colinear PCET, in which electron and proton are transferred from the same donor to the same acceptor between Y₃₅₆ and Y₇₃₁. The results highlight the important role of state-of-the-art EPR spectroscopy to decipher this mechanism.

Resolution of chemical shift anisotropy in 19F ENDOR spectroscopy at 263 GHz/9.4 T

Pulsed ¹⁹F ENDOR spectroscopy provides a selective method for measuring angstrom to nanometer distances in structural biology. Here, the performance of ¹⁹F ENDOR at fields of 3.4 T and 9.4 T is compared using model compounds containing one to three ¹⁹F atoms. CF₃ groups are included in two compounds, for which the possible occurrence of uniaxial rotation might affect the distance distribution. At 9.4 T, pronounced asymmetric features are observed in many of the presented ¹⁹F ENDOR spectra. Data analysis by spectral simulations shows that these features arise from the chemical shift anisotropy (CSA) of the ¹⁹F nuclei. This asymmetry is also observed at 3.4 T, albeit to a much smaller extent, confirming the physical origin of the effect. The CSA parameters are well consistent with DFT predicted values and can be extracted from simulation of the experimental data in favourable cases, thereby providing additional information about the geometrical and electronic structure of the spin system. The feasibility of resolving the CSA at 9.4 T provides important information for the interpretation of line broadening in ENDOR spectra also at lower fields, which is relevant for developing methods to extract distance distributions from ¹⁹F ENDOR spectra.

Statistical analysis of ENDOR spectra

Electron-nuclear double resonance (ENDOR) measures the hyperfine interaction of magnetic nuclei with paramagnetic centers and is hence a powerful tool for spectroscopic investigations extending from biophysics to material science. Progress in microwave technology and the recent availability of commercial electron paramagnetic resonance (EPR) spectrometers up to an electron Larmor frequency of 263 GHz now open the opportunity for a more quantitative spectral analysis. Using representative spectra of a prototype amino acid radical in a biologically relevant enzyme, the [Formula: see text] in Escherichia coli ribonucleotide reductase, we developed a statistical model for ENDOR data and conducted statistical inference on the spectra including uncertainty estimation and hypothesis testing. Our approach in conjunction with ¹H/²H isotopic labeling of [Formula: see text] in the protein unambiguously established new unexpected spectral contributions. Density functional theory (DFT) calculations and ENDOR spectral simulations indicated that these features result from the beta-methylene hyperfine coupling and are caused by a distribution of molecular conformations, likely important for the biological function of this essential radical. The results demonstrate that model-based statistical analysis in combination with state-of-the-art spectroscopy accesses information hitherto beyond standard approaches.