Understanding heme proteins with hyperfine spectroscopy

Characterization of the active site in the thiocyanate-forming protein from Thlaspi arvense (TaTFP) using EPR spectroscopy

Glucosinolates are plant thioglucosides, which act as chemical defenses. Upon tissue damage, their myrosinase-catalyzed hydrolysis yields aglucones that rearrange to toxic isothiocyanates. Specifier proteins such as thiocyanate-forming protein from Thlaspi arvense (TaTFP) are non-heme iron proteins, which capture the aglucone to form alternative products, e.g. nitriles or thiocyanates. To resolve the electronic state of the bound iron cofactor in TaTFP, we applied continuous wave electron paramagnetic resonance (CW EPR) spectroscopy at X-and Q-band frequencies (∼9.4 and ∼34 GHz). We found characteristic features of high spin and low spin states of a d 5 electronic configuration and local rhombic symmetry during catalysis. We monitored the oxidation states of bound iron during conversion of allylglucosinolate by myrosinase and TaTFP in presence and absence of supplemented Fe2+. Without added Fe2+, most high spin features of bound Fe3+ were preserved, while different g'-values of the low spin part indicated slight rearrangements in the coordination sphere and/or structural geometry. We also examined involvement of the redox pair Fe3+/Fe2 in samples with supplemented Fe2+. The absence of any EPR signal related to Fe3+ or Fe2+ using an iron-binding deficient TaTFP variant allowed us to conclude that recorded EPR signals originated from the bound iron cofactor.

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.

The effect of pH and nitrite on the haem pocket of GLB-33, a globin-coupled neuronal transmembrane receptor of Caenorhabditis elegans

Out of the 34 globins in Caenorhabditis elegans, GLB-33 is a putative globin-coupled transmembrane receptor with a yet unknown function. The globin domain (GD) contains a particularly hydrophobic haem pocket, that rapidly oxidizes to a low-spin hydroxide-ligated haem state at physiological pH. Moreover, the GD has one of the fastest nitrite reductase activity ever reported for globins. Here, we use a combination of electronic circular dichroism, resonance Raman and electron paramagnetic resonance (EPR) spectroscopy with mass spectrometry to study the pH dependence of the ferric form of the recombinantly over-expressed GD in the presence and absence of nitrite. The competitive binding of nitrite and hydroxide is examined as well as nitrite-induced haem modifications at acidic pH. Comparison of the spectroscopic results with data from other haem proteins allows to deduce the important effect of Arg at position E10 in stabilization of exogenous ligands. Furthermore, continuous-wave and pulsed EPR indicate that ligation of nitrite occurs in a nitrito mode at pH 5.0 and above. At pH 4.0, an additional formation of a nitro-bound haem form is observed along with fast formation of a nitri-globin.

A unified description for polarization-transfer mechanisms in magnetic resonance in static solids: Cross polarization and DNP

Polarization transfers are crucial building blocks in magnetic resonance experiments, i.e., they can be used to polarize insensitive nuclei and correlate nuclear spins in multidimensional nuclear magnetic resonance (NMR) spectroscopy. The polarization can be transferred either across different nuclear spin species or from electron spins to the relatively low-polarized nuclear spins. The former route occurring in solid-state NMR can be performed via cross polarization (CP), while the latter route is known as dynamic nuclear polarization (DNP). Despite having different operating conditions, we opinionate that both mechanisms are theoretically similar processes in ideal conditions, i.e., the electron is merely another spin-1/2 particle with a much higher gyromagnetic ratio. Here, we show that the CP and DNP processes can be described using a unified theory based on average Hamiltonian theory combined with fictitious operators. The intuitive and unified approach has allowed new insights into the cross-effect DNP mechanism, leading to better design of DNP polarizing agents and extending the applications beyond just hyperpolarization. We explore the possibility of exploiting theoretically predicted DNP transients for electron–nucleus distance measurements—such as routine dipolar-recoupling experiments in solid-state NMR.

A dimeric chlorophyll electron acceptor differentiates type I from type II photosynthetic reaction centers

This research addresses one of the most compelling issues in the field of photosynthesis, namely, the role of the accessory chlorophyll molecules in primary charge separation. Using a combination of empirical and computational methods, we demonstrate that the primary acceptor of photosystem (PS) I is a dimer of accessory and secondary chlorophyll molecules, Chl_2A and Chl_3A, with an asymmetric electron charge density distribution. The incorporation of highly coupled donors and acceptors in PS I allows for extensive delocalization that prolongs the lifetime of the charge-separated state, providing for high quantum efficiency. The discovery of this motif has widespread implications ranging from the evolution of naturally occurring reaction centers to the development of a new generation of highly efficient artificial photosynthetic systems.

Video abstract

Two-dimensional HYSCORE spectroscopy reveals a histidine imidazole as the axial ligand to Chl3A in the M688HPsaA genetic variant of Photosystem I

Recent studies on Photosystem I (PS I) have shown that the six core chlorophyll a molecules are highly coupled, allowing for efficient creation and stabilization of the charge-separated state. One area of particular interest is the identity and function of the primary acceptor, A₀, as the factors that influence its ultrafast processes and redox properties are not yet fully elucidated. It was recently shown that A₀ exists as a dimer of the closely-spaced Chl₂/Chl₃ molecules wherein the reduced A₀^- state has an asymmetric distribution of electron spin density that favors Chl₃. Previous experimental work in which this ligand was changed to a hard base (histidine, M688H_PsaA) revealed severely impacted electron transfer processes at both the A₀ and A₁ acceptors; molecular dynamics simulations further suggested two distinct conformations of PS I in which the His residue coordinates and forms a hydrogen bond to the A₀ and A₁ cofactors, respectively. In this study, we have applied 2D HYSCORE spectroscopy in conjunction with molecular dynamics simulations and density functional theory calculations to the study of the M688H_PsaA variant. Analysis of the hyperfine parameters demonstrates that the His imidazole serves as the axial ligand to the central Mg²⁺ ion in Chl_3A in the M688H_PsaA variant. Although the change in ligand identity does not alter delocalization of electron density over the Chl₂/Chl₃ dimer, a small shift in the asymmetry of delocalization, coupled with the electron withdrawing properties of the ligand, most likely accounts for the inhibition of forward electron transfer in the His-ligated conformation.

Two-dimensional 67Zn HYSCORE spectroscopy reveals that a Zn-bacteriochlorophyll aP′ dimer is the primary donor (P840) in the type-1 reaction centers of Chloracidobacterium thermophilum

Using pulsed EPR spectroscopy and isotopic labeling we demonstrate that reaction centers of Chloracidobacterium thermophilum have an unusual primary donor that is a dimer of Zn-bacteriochlorophyll a_P′ molecules.

AMOEBA Polarizable Force Field Parameters of the Heme Cofactor in Its Ferrous and Ferric Forms

We report the first parameters of the heme redox cofactors for the polarizable AMOEBA force field in both the ferric and ferrous forms. We consider two types of complexes, one with two histidine side chains as axial ligands and one with a histidine and a methionine side chain as ligands. We have derived permanent multipoles from second-order Møller-Plesset perturbation theory (MP2). The sets of parameters have been validated in a first step by comparison of AMOEBA interaction energies of heme and a collection of biologically relevant molecules with MP2 and Density Functional Theory (DFT) calculations. In a second validation step, we consider interaction energies with large aggregates comprising around 80 H₂O molecules. These calculations are repeated for 30 structures extracted from semiempirical PM7 DM simulations. Very encouraging agreement is found between DFT and the AMOEBA force field, which results from an accurate treatment of electrostatic interactions. We finally report long (10 ns) MD simulations of cytochromes in two redox states with AMOEBA testing both the 2003 and 2014 AMOEBA water models. These simulations have been carried out with the TINKER-HP (High Performance) program. In conclusion, owing to their ubiquity in biology, we think the present work opens a wide array of applications of the polarizable AMOEBA force field on hemeproteins.