On the calculation of G tensors of organic radicals

Manipulating D-A interaction to achieve stable photoinduced organic radicals in triphenylphosphine crystals

New strategies for the design and synthesis of stable organic radicals without additives are highly desirable. Herein, we design a series of donor-acceptor structured triarylphosphines and disclose the fast color change triggered by UV-irradiation in the crystalline state. Photoinduced organic radicals are undoubtedly verified and proved to be the reason for the color change by time-dependent and quantitative electron paramagnetic resonance analysis, X-ray crystallographic analysis, and theoretical calculations. It is revealed that the intrinsic symmetry breaking of peripheral architecture helps to form continuous molecular chains by donor-acceptor counterpart pairing. Intermolecular electron-transfer occurs among molecular chains and results in radical ion pairs upon photoirradiation.

Structure of the biliverdin radical intermediate in phycocyanobilin:ferredoxin oxidoreductase identified by high-field EPR and DFT

The cyanobacterial enzyme phycocyanobilin:ferredoxin oxidoreductase (PcyA) catalyzes the two-step four-electron reduction of biliverdin IXalpha to phycocyanobilin, the precursor of biliprotein chromophores found in phycobilisomes. It is known that catalysis proceeds via paramagnetic radical intermediates, but the structure of these intermediates and the transfer pathways for the four protons involved are not known. In this study, high-field electron paramagnetic resonance (EPR) spectroscopy of frozen solutions and single crystals of the one-electron reduced protein-substrate complex of two PcyA mutants D105N from the cyanobacteria Synechocystis sp. PCC6803 and Nostoc sp. PCC7120 are examined. Detailed analysis of Synechocystis D105N mutant spectra at 130 and 406 GHz reveals a biliverdin radical with a very narrow g tensor with principal values 2.00359(5), 2.00341(5), and 2.00218(5). Using density-functional theory (DFT) computations to explore the possible protonation states of the biliverdin radical, it is shown that this g tensor is consistent with a biliverdin radical where the carbonyl oxygen atoms on both the A and the D pyrrole rings are protonated. This experimentally confirms the reaction mechanism recently proposed (Tu, et al. Biochemistry 2007, 46, 1484).

EPR and ENDOR study of radiation-induced radical formation in purines: hypoxanthine hydrochloride monohydrate crystals X-irradiated at 10 K

Electron paramagnetic resonance (EPR) and electron-nuclear double resonance (ENDOR) study of hypoxanthine.HCl.H(2)O crystals irradiated at low temperatures (10 K) identified three radical species. In these crystals, the parent molecules exist in a cationic form with a proton at N7. R1 was the product of net hydrogen addition to N3 and exhibited alpha-proton hyperfine couplings to HC2, HN1, HC8, and HN3. The coupling to HC2 has an isotropic component smaller than usual, evidently an indication that the bonds to C2 are nonplanar. R2 was the product of net hydrogen loss from N7, equivalent to the one-electron oxidation product of neutral hypoxanthine, and exhibited alpha-proton hyperfine couplings to HC2 and HC8. Both couplings are characteristic of planar bonding arrangements at the centers of spin. R3 was provisionally identified as the product of net hydrogen addition to O6 and exhibited hyperfine alpha-proton couplings to HC8 and NH1. To identify the set of radicals, the experiments employed four crystal types: normal, deuterated only at NH positions, deuterated at HC8 and NH positions, and deuterated at HC8 only. The low-temperature data also showed clear evidence for H/D isotope effects in formation and/or stabilization of all radicals. To aid and support the identifications, the experimental results were compared to DFT calculations performed on a variety of radical structures plausible for the parent molecule and molecular packing within the crystal.

Substituent effect on g-tensor: multifrequency ESR study and DFT calculation of polycrystalline phenoxyl radicals in diamagnetic crystals

The substituent effect on the g-tensor of polycrystalline 2,6-di-tert-butyl phenoxyl radical derivatives diluted in diamagnetic crystals was investigated using multifrequency ESR spectroscopy and DFT calculations. It was revealed that the g-tensors of the series of phenoxyl radical derivatives essentially have an orthorhombic symmetry. For some radicals, the hyperfine-splitting tensors from the para groups were resolved. The interpretations and the assignments of the spin-Hamiltonian parameters were confirmed with computer simulations in all bands. The DFT-calculated g-tensors were consistent with the experimental g-tensors. Furthermore, the shifts Delta(g) from the free electron ge were analyzed in details as the sum of three contributions. The spin-orbit interactions were found to be the dominant factor with regard to the Delta(g). With a focus on the s-o term, thus, the relationship of the g-values and the electronic excited states was explained by visualizing the molecular orbitals of the phenoxyl radical derivatives. This study thus showed the very significant potential of the combination of a multi-frequency ESR approach and a DFT calculation to advanced ESR analysis, particularly, g-tensor analysis, even for a powder-sample radical.

Structure of radicals from X-irradiated guanine derivatives: an experimental and computational study of sodium guanosine dihydrate single crystals

In sodium guanosine dihydrate single crystals, the guanine moiety is deprotonated at N1 due to growth from high-pH (>12) solutions. Electron paramagnetic resonance (EPR) and electron-nuclear double resonance (ENDOR) studies of crystals X-irradiated at 10 K detected evidence for three radical forms. Radical R1, characterized by two proton and two nitrogen hyperfine interactions, was identified as the product of net hydrogenation at N7 of the N1-deprotonated guanine unit. R1 exhibited an unusually distorted structure leading to net positive isotropic components of the hydrogen alpha-couplings. Radical R2, characterized by one proton and one nitrogen hyperfine coupling, was identified as the primary electron-loss product. This product is equivalent to that of deprotonation at N1 by the guanine cation and represents the first ENDOR characterization of that product. Radical R3, characterized by a single hydrogen hyperfine coupling, was identified as the product of net dehydrogenation at C1' of the ribose moiety. The identification of radicals R1-R3 was supported by density functional theory (DFT) calculations on several possible structures using the B3LYP/6-311G(2df,p)//6-31G(d,p) approach. Radical R4, detected after warming the crystals to room temperature, was identified as the well-known product of net hydrogenation of C8 of the (N1-deprotonated) guanine component. Radical R1, evidently formed by protonation of the primary electron addition product, was present as roughly 60% of the total radicals detected at 10 K. Radical R2 was present as roughly 27% of the total yield, and the concentration of R3 contributed the remaining 13%. R3 is evidently the product of one-electron oxidation followed by deprotonation; thus, the balance of oxidation and reduction products is approximately equal within experimental uncertainty.

Protein–cofactor interactions in bacterial reaction centers from Rhodobacter sphaeroides R-26: Effect of hydrogen bonding on the electronic and geometric structure of the primary quinone. A density functional theory study

The effect of hydrogen bonding to the primary quinone (Q(A) and Q(*)(-)(A)) in bacterial reaction centers was studied using density functional theory (DFT) calculations. The charge neutral state Q(A) was investigated by optimizing the hydrogen atom positions of model systems extracted from 15 different X-ray structures. From this analysis, mean values of the H-bond lengths and directions were derived. It was found that the N(delta)-H of His M219 forms a shorter H-bond to Q(A) than the N-H of Ala M260. The H-bond of His M219 is linear and more twisted out of the quinone plane. The radical anion Q(*)(-)(A) in the protein environment was investigated by using a mixed quantum mechanics/molecular mechanics (QM/MM) approach. Two geometry optimizations with a different number of flexible atoms were performed. H-bond lengths were obtained and spectroscopic parameters calculated, i.e. the hyperfine and nuclear quadrupole couplings of magnetic nuclei coupled to the radical. Good agreement was found with the results provided by EPR/ENDOR spectroscopy. This implies that the calculated lengths and directions of the H-bonds to Q(*)(-)(A) are reliable values. From a comparison of the neutral and reduced state of Q(A) it was concluded that the H-bond distances are shortened by approximately 0.17 Angstroms (His M219) and approximately 0.13 Angstroms (Ala M260) upon single reduction of the quinone. It is shown that the point-dipole approximation can not be used for an estimation of H-bond lengths from measured hyperfine couplings in a system with out-of-plane H-bonding. In contrast, the evaluation of the nuclear quadrupole couplings of (2)H nuclei substituted in the hydrogen bonds yields H-bond lengths close to the values that were deduced from DFT geometry optimizations. The significance of hydrogen bonding to the quinone cofactors in biological systems is discussed.

Combining high-field EPR with site-directed spin labeling reveals unique information on proteins in action

In the last decade, joint efforts of biologists, chemists and physicists have helped in understanding the dominant factors determining specificity and directionality of transmembrane transfer processes in proteins. In this endeavor, electron paramagnetic resonance (EPR) spectroscopy has played an important role. Characteristic examples of such determining factors are hydrogen-bonding patterns and polarity effects of the microenvironment of protein sites involved in the transfer process. These factors may undergo characteristic changes during the reaction and, thereby, control the efficiency of biological processes, e.g. light-induced electron and proton transfer across photosynthetic membranes or ion-channel formation of bacterial toxins. In case the transfer process does not involve stable or transient paramagnetic species or states, site-directed spin labeling with suitable nitroxide radicals still allows EPR techniques to be used for studying structure and conformational dynamics of the proteins in action. By combining site-directed spin labeling with high-field/high-frequency EPR, unique information on the proteins is revealed, which is complementary to that of X-ray crystallography, solid-state NMR, FRET, fast infrared and optical spectroscopic techniques. The main object of this publication is twofold: (i) to review our recent spin-label high-field EPR work on the bacteriorhodopsin light-driven proton pump from Halobacterium salinarium and the Colicin A ion-channel forming bacterial toxin produced in Escherichia coli, (ii) to report on novel high-field EPR experiments for probing site-specific pK(a) values in protein systems by means of pH-sensitive nitroxide spin labels. Taking advantage of the improved spectral and temporal resolution of high-field EPR at 95 GHz/3.4 T and 360 GHz/12.9 T, as compared to conventional X-band EPR (9.5 GHz/0.34 T), detailed information on the transient intermediates of the proteins in biological action is obtained. These intermediates can be observed and characterized while staying in their working states on biologically relevant timescales. The paper concludes with an outlook of ongoing high-field EPR experiments on site-specific protein mutants in our laboratories at FU Berlin and Osnabrück.

Gauge invariance of the spin-other-orbit contribution to the g-tensors of electron paramagnetic resonance

The spin-other-orbit (SOO) contribution to the g-tensor (DeltagSOO) of electron paramagnetic resonance arises due to the interaction of electron-spin magnetic moment with the magnetic field produced by the orbital motion of other electrons. A similar mechanism is responsible for the leading term in nuclear magnetic-shielding tensors sigma. We demonstrate that analogous to sigma, paramagnetic DeltagSOO contribution exhibits a pronounced dependence on the choice of the magnetic-field gauge. The gauge corrections to DeltagSOO are similar in magnitude, and opposite in sign, to the paramagnetic SOO term. We calculate gauge-invariant DeltagSOO values using gauge-including atomic orbitals and density-functional theory. For organic radicals, complete gauge-invariant DeltagSOO values typically amount to less than 500 parts per million (ppm), and are small compared to other g-tensor contributions. For the first-row transition-metal compounds, DeltagSOO may contribute several thousand ppm to the g-tensor, but are negligible compared to the remaining deviations from experiment. With popular choices for the magnetic-field gauge, the individual gauge-variant contributions may be an order of magnitude higher, and do not provide a reliable estimation of DeltagSOO.

Ab Initio Molecular Dynamics Simulations and g-Tensor Calculations of Aqueous Benzosemiquinone Radical Anion: Effects of Regular and “T-Stacked” Hydrogen Bonds

Car-Parrinello molecular dynamics (CP-MD) simulations of the benzosemiquinone radical anion in aqueous solution have been performed at ambient conditions. Analysis of the trajectory shows not only extensive hydrogen bonding to the carbonyl oxygen atoms (ca. 4-5.6 water molecules depending on distance criteria), but also relatively long-lived "T-stacked" hydrogen bonds to the semiquinone pi-system. These results are discussed in the context of recent findings on semiquinone-protein interactions in photosynthetic reaction centers, and of EPR and vibration spectroscopical data for the aqueous system. Snapshots from the CP-MD trajectory are used for the first quantum chemical analyses of dynamical effects on electronic g-tensors, using cluster models and a recently developed density functional method. In particular, the effects of intermolecular hydrogen-bond dynamics on the g-tensor components are examined, in comparison with recent EPR and ENDOR studies.

Hydrogen Bond Geometries from Electron Paramagnetic Resonance and Electron−Nuclear Double Resonance Parameters: Density Functional Study of Quinone Radical Anion−Solvent Interactions

Density functional theory was used to study the impact of hydrogen bonding on the p-benzosemiquinone radical anion BQ(*-) in coordination with water or alcohol molecules. After complete geometry optimizations, (1)H, (13)C, and (17)O hyperfine as well as (2)H nuclear quadrupole coupling constants and the g-tensor were computed. The suitability of different model systems with one, two, four, and 20 water molecules was tested; best agreement between theory and experiment could be obtained for the largest model system. Q-band pulse (2)H electron-nuclear double resonance (ENDOR) experiments were performed on BQ(*-) in D(2)O. They compare very well with the spectra simulated by use of the theoretical values from density functional theory. For BQ(*-) in coordination with four water or alcohol molecules, rather similar hydrogen-bond lengths between 1.75 and 1.78 A were calculated. Thus, the computed electron paramagnetic resonance (EPR) parameters are hardly distinguishable for the different solvents, in agreement with experimental findings. Furthermore, the distance dependence of the EPR parameters on the hydrogen-bond length was studied. The nuclear quadrupole and the dipolar hyperfine coupling constants of the bridging hydrogens show the expected dependencies on the H-bond length R(O.H). A correlation was obtained for the g-tensor. It is shown that the point-dipole model is suitable for the estimation of hydrogen-bond lengths from anisotropic hyperfine coupling constants of the bridging (1)H nuclei for H-bond lengths larger than approximately 1.7 A. Furthermore, the estimation of H-bond lengths from (2)H nuclear quadrupole coupling constants of bridging deuterium nuclei by empirical relations is discussed.

Quantum chemical calculations of spectroscopic properties of metalloproteins and model compounds: EPR and Mössbauer properties

Recently developed theoretical methods to predict EPR and Mössbauer parameters open the way for close interactions between theorists and experimentalists to elucidate the geometric and electronic structures of metalloenzymes and model complexes and to obtain insight into their reactive properties. Spectral calculations (g-values, hyperfine couplings, zero-field splittings, isomer shifts and quadrupole splittings) are also a means to validate theoretical approaches and therefore complement the prediction of geometries, reaction energies and transition states.

Density functional calculations of electronic g-tensors for semiquinone radical anions. The role of hydrogen bonding and substituent effects

A recently developed density functional approach has been used to carry out a systematic computational study of electronic g-tensors for a series of 1,4-semiquinone radical anions. Good agreement with high-field EPR data in frozen 2-propanol is achieved only after taking into account the significant reduction of g-tensor anisotropy caused by hydrogen bonding to solvent molecules. The comparison of various model systems for the first solvation shell suggests two hydrogen bonds from 2-propanol molecules to each of the carbonyl groups of the radical anions, and one additional hydrogen bond to each of the methoxy groups in ubiquinone systems. 2-Propanol makes stronger hydrogen bonds than water and thus influences g-tensor anisotropy more strongly. Substituent effects at the semiquinone are reproduced quantitatively by the calculations. The g-tensor anisotropy is influenced significantly by the conformations of methyl and methoxy substituents, with opposite contributions. Analyses and interpretations of the interrelations between structure, bonding, and spectroscopic data are provided. The relevance of the computational results for the EPR spectroscopy of semiquinone radical anions in photosynthetic reaction centers is discussed.