A case study of density functional theory and domain-based local pair natural orbital coupled cluster for vibrational effects on EPR hyperfine coupling constants: vibrational perturbation theory versus ab initio molecular dynamics

Theoretical terahertz spectroscopy of free radical solutes in solution: an EPR spin probe in water

Free radical species are used as spin labels in electron paramagnetic resonance (EPR) spectroscopy of biomolecular systems in water, for instance in the frame of Overhauser dynamic nuclear polarization (ODNP) relaxometry to probe the local hydration water dynamics close to protein surfaces in aqueous environments. Widely used in this context are nitroxide spin probes such as TEMPO, PROXYL or MTSL derivatives. Here, we study the THz spectroscopy of HMI (2,2,3,4,5,5-HexaMethylImidazolidin-1-oxyl) in water at ambient conditions which has been recently investigated as to how its EPR properties depend on its solvation pattern in water. To enable theoretical THz spectroscopy of molecular radicals in solution, we have generalized well-established methodologies for THz spectral decomposition of closed-shell systems, namely the supermolecular solvation complex (SSC) and cross-correlation analysis (CCA) techniques, to open-shell polyatomic solute species in water. Based on this methodological advance, we have decomposed and assigned the THz response of HMI including its solvation shell by employing the generalized SSC and CCA methods to cope with the open-shell character of this free radical solute, in particular its unpaired electron localized at the nitroxy group. We reveal that the main modulations of the far-IR spectrum of HMI are dominated by the low-frequency intramolecular modes of the spin probe molecule itself while the solvation of its two hydrogen bonding sites contribute much less intensely in this spectral window. Finally, we have computed THz spectra of HMI with its local solvation water with an aim to provide a theoretical analogue of experimental THz difference spectroscopy. Beyond the specific case, our decomposition methodology that now is able to include open shells can be applied in future work to analyze the low-frequency vibrational response of the solvation shell of other free radicals in solution.

Bidentate Substrate Binding Mode in Oxalate Decarboxylase

Oxalate decarboxylase is an Mn- and O2-dependent enzyme in the bicupin superfamily that catalyzes the redox-neutral disproportionation of the oxalate monoanion to form carbon dioxide and formate. Its best-studied isozyme is from Bacillus subtilis where it is stress-induced under low pH conditions. Current mechanistic schemes assume a monodentate binding mode of the substrate to the N-terminal active site Mn ion to make space for a presumed O2 molecule, despite the fact that oxalate generally prefers to bind bidentate to Mn. We report on X-band 13C-electron nuclear double resonance (ENDOR) experiments on 13C-labeled oxalate bound to the active-site Mn(II) in wild-type oxalate decarboxylase at high pH, the catalytically impaired W96F mutant enzyme at low pH, and Mn(II) in aqueous solution. The ENDOR spectra of these samples are practically identical, which shows that the substrate binds bidentate (κO, κO') to the active site Mn(II) ion. Domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) calculations of the expected 13C hyperfine coupling constants for bidentate bound oxalate predict ENDOR spectra in good agreement with the experiment, supporting bidentate bound substrate. Geometry optimization of a substrate-bound minimal active site model by density functional theory shows two possible substrate coordination geometries, bidentate and monodentate. The bidentate structure is energetically preferred by ~4.7 kcal/mol. Our results revise a long-standing hypothesis regarding substrate binding in the enzyme and suggest that dioxygen does not bind to the active site Mn ion after substrate binds. The results are in agreement with our recent mechanistic hypothesis of substrate activation via a long-range electron transfer process involving the C-terminal Mn ion.

Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory

Nitroxides are common EPR sensors of microenvironmental properties such as polarity, numbers of H-bonds, pH, and so forth. Their solvation in an aqueous environment is facilitated by their high propensity to form H-bonds with the surrounding water molecules. Their g- and A-tensor elements are key parameters to extracting the properties of their microenvironment. In particular, the gxx value of nitroxides is rich in information. It is known to be characterized by discrete values representing nitroxide populations previously assigned to have different H-bonds with the surrounding waters. Additionally, there is a large g-strain, that is, a broadening of g-values associated with it, which is generally correlated with environmental and structural micro-heterogeneities. The g-strain is responsible for the frequency dependence of the apparent line width of the EPR spectra, which becomes evident at high field/frequency. Here, we address the molecular origin of the gxx heterogeneity and of the g-strain of a nitroxide moiety (HMI: 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water. To treat the solvation effect on the g-strain, we combined a multi-frequency experimental approach with ab initio molecular dynamics simulations for structural sampling and quantum chemical EPR property calculations at the highest realistically affordable level, including an explicitly micro-solvated HMI ensemble and the embedded cluster reference interaction site model. We could clearly identify the distinct populations of the H-bonded nitroxides responsible for the gxx heterogeneity experimentally observed, and we dissected the role of the solvation shell, H-bond formation, and structural deformation of the nitroxide in the creation of the g-strain associated with each nitroxide subensemble. Two contributions to the g-strain were identified in this study. The first contribution depends on the number of hydrogen bonds formed between the nitroxide and the solvent because this has a large and well-understood effect on the gxx-shift. This contribution can only be resolved at high resonance frequencies, where it leads to distinct peaks in the gxx region. The second contribution arises from configurational fluctuations of the nitroxide that necessarily lead to g-shift heterogeneity. These contributions cannot be resolved experimentally as distinct resonances but add to the line broadening. They can be quantitatively analyzed by studying the apparent line width as a function of microwave frequency. Interestingly, both theory and experiment confirm that this contribution is independent of the number of H-bonds. Perhaps even more surprisingly, the theoretical analysis suggests that the configurational fluctuation broadening is not induced by the solvent but is inherently present even in the gas phase. Moreover, the calculations predict that this broadening decreases upon solvation of the nitroxide.

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.

Electron correlation and vibrational effects in predictions of paramagnetic NMR shifts

Electronic structure calculations are fundamentally important for the interpretation of nuclear magnetic resonance (NMR) spectra from paramagnetic systems that include organometallic and inorganic compounds, catalysts, or metal-binding sites in proteins. Prediction of induced paramagnetic NMR shifts requires knowledge of electron paramagnetic resonance (EPR) parameters: the electronic g tensor, zero-field splitting D tensor, and hyperfine A tensor. The isotropic part of A, called the hyperfine coupling constant (HFCC), is one of the most troublesome properties for quantum chemistry calculations. Yet, even relatively small errors in calculations of HFCC tend to propagate into large errors in the predicted NMR shifts. The poor quality of A tensors that are currently calculated using density functional theory (DFT) constitutes a bottleneck in improving the reliability of interpretation of the NMR spectra from paramagnetic systems. In this work, electron correlation effects in calculations of HFCCs with a hierarchy of ab initio methods were assessed, and the applicability of different levels of DFT approximations and the coupled cluster singles and doubles (CCSD) method was tested. These assessments were performed for the set of selected test systems comprising an organic radical, and complexes with transition metal and rare-earth ions, for which experimental data are available. Severe deficiencies of DFT were revealed but the CCSD method was able to deliver good agreement with experimental data for all systems considered, however, at substantial computational costs. We proposed a more computationally tractable alternative, where the A was computed with the coupled cluster theory exploiting locality of electron correlation. This alternative is based on the domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) method. In this way the robustness and reliability of the coupled cluster theory were incorporated into the modern formalism for the prediction of induced paramagnetic NMR shifts, and became applicable to systems of chemical interest. This approach was verified for the bis(cyclopentadienyl)vanadium(II) complex (Cp₂V; vanadocene), and the metal-binding site of the Zn²⁺ → Co²⁺ substituted superoxide dismutase (SOD) metalloprotein. Excellent agreement with experimental NMR shifts was achieved, which represented a substantial improvement over previous theoretical attempts. The effects of vibrational corrections to orbital shielding and hyperfine tensor were evaluated and discussed within the second-order vibrational perturbation theory (VPT2) framework.

High Level Electronic Structure Calculation of Molecular Solid-State NMR Shielding Constants

In this work, we present a quantum mechanics/molecular mechanics (QM/MM) approach for the computation of solid-state nuclear magnetic resonance (SS-NMR) shielding constants (SCs) for molecular crystals. Besides applying standard-DFT functionals like GGAs (PBE), meta-GGAs (TPSS), and hybrids (B3LYP), we apply a double-hybrid (DSD-PBEP86) functional as well as MP2, using the domain-based local pair natural orbital (DLPNO) formalism, to calculate the NMR SCs of six amino acid crystals. All the electronic structure methods used exhibit good correlation of the NMR shieldings with respect to experimental chemical shifts for both 1H and 13C. We also find that local electronic structure is much more important than the long-range electrostatic effects for these systems, implying that cluster approaches using all-electron/Gaussian basis set methods might offer great potential for predictive computations of solid-state NMR parameters for organic solids.

A Joint Venture of Ab Initio Molecular Dynamics, Coupled Cluster Electronic Structure Methods, and Liquid-State Theory to Compute Accurate Isotropic Hyperfine Constants of Nitroxide Probes in Water

The isotropic hyperfine coupling constant (HFCC, Aiso) of a pH-sensitive spin probe in a solution, HMI (2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water, is computed using an ensemble of state-of-the-art computational techniques and is gauged against X-band continuous wave electron paramagnetic resonance (EPR) measurement spectra at room temperature. Fundamentally, the investigation aims to delineate the cutting edge of current first-principles-based calculations of EPR parameters in aqueous solutions based on using rigorous statistical mechanics combined with correlated electronic structure techniques. In particular, the impact of solvation is described by exploiting fully atomistic, RISM integral equation, and implicit solvation approaches as offered by ab initio molecular dynamics (AIMD) of the periodic bulk solution (using the spin-polarized revPBE0-D3 hybrid functional), embedded cluster reference interaction site model integral equation theory (EC-RISM), and polarizable continuum embedding (using CPCM) of microsolvated complexes, respectively. HFCCs are obtained from efficient coupled cluster calculations (using open-shell DLPNO-CCSD theory) as well as from hybrid density functional theory (using revPBE0-D3). Re-solvation of "vertically desolvated" spin probe configuration snapshots by EC-RISM embedding is shown to provide significantly improved results compared to CPCM since only the former captures the inherent structural heterogeneity of the solvent close to the spin probe. The average values of the Aiso parameter obtained based on configurational statistics using explicit water within AIMD and from EC-RISM solvation are found to be satisfactorily close. Using either such explicit or RISM solvation in conjunction with DLPNO-CCSD calculations of the HFCCs provides an average Aiso parameter for HMI in aqueous solution at 300 K and 1 bar that is in good agreement with the experimentally determined one. The developed computational strategy is general in the sense that it can be readily applied to other spin probes of similar molecular complexity, to aqueous solutions beyond ambient conditions, as well as to other solvents in the longer run.

The distance between g-tensors of nitroxide biradicals governs MAS-DNP performance: The case of the bTurea family

Bis-nitroxide radicals are common polarizing agents (PA), used to enhance the sensitivity of solid-state NMR experiments via Magic Angle Spinning Dynamic Nuclear Polarization (MAS-DNP). These biradicals can increase the proton spin polarization through the Cross-Effect (CE) mechanism, which requires PAs with at least two unpaired electrons. The relative orientation of the bis-nitroxide moieties is critical to ensure efficient polarization transfer. Recently, we have defined a new quantity, the distance between g-tensors, that correlates the relative orientation of the nitroxides with the ability to polarize the surrounding nuclei. Here we analyse experimentally and theoretically a series of biradicals belonging to the bTurea family, namely bcTol, AMUPol and bcTol-M. They differ by the degree of substitution on the urea bridge that connects the two nitroxides. Using quantitative simulations developed for moderate MAS frequencies, we show that these modifications mostly affect the relative orientations of the nitroxide, i.e. the length and distribution of the distance between the g-tensors, that in turn impacts both the steady state nuclear polarization/depolarization as well as the build-up times. The doubly substituted urea bridge favours a large distance between the g-tensors, which enables bcTol-M to provide ∊_on/off>200 at 14.1 T/600 MHz/395 GHz with build-up times of 3.8 s using a standard homogenous solution. The methodology described herein was used to show how the conformation of the spirocyclic rings flanking the nitroxide function in the recently described c- and o-HydrOPol affects the distance between the g-tensors and thereby polarization performance.

Performance of the DLPNO-CCSD and recent DFT methods in the calculation of isotropic and dipolar contributions to 14N hyperfine coupling constants of nitroxide radicals

In the present study, the performance of a set of density functionals: BP86, PBE, OLYP, BEEF, PBEpow, TPSS, SCAN, PBEGXPBE, M06L, MN15L, B3LYP, PBE0, mPW1PW, B97, BHandHLYP, mPW1PW, B98, TPSS0, PBE1KCIS, SCAN0, M06, M06-2X, MN15, CAM-B3LYP, ωB97x, B2PLYP, and the B3LYP/N07D and PBE/N07D schemes in the calculation of the ¹⁴N anisotropic hyperfine coupling (HFC) constants of a set of 23 nitroxide radicals is evaluated. The results are compared with those obtained with the DLPNO-CCSD method and experimental HFC values. Harmonic contribution to the ¹⁴N HFC vibrational correction was calculated at the revPBE0/def2-TZVPP level and included in the evaluation. With the vibrational correction, the DLPNO-CCSD method yielded HFC values in good agreement with the experiment (mean absolute deviation (MAD) = 0.3 G for the dipole-dipole contribution and MAD = 0.8 G for the contact coupling contribution). The best DFT results are obtained using the M06 functional with MAD = 0.2 G for the dipole-dipole contribution and MAD = 0.7 G for the contact coupling contribution. In general, vibrational correction significantly improved most DFT functionals' performance but did not change its overall ranking.

A Quantum Chemistry View on Two Archetypical Paramagnetic Pentacoordinate Nickel(II) Complexes Offers a Fresh Look on Their NMR Spectra

Quantum chemical methods for calculating paramagnetic NMR observables are becoming increasingly accessible and are being included in the inorganic chemistry practice. Here, we test the performance of these methods in the prediction of proton hyperfine shifts of two archetypical high-spin pentacoordinate nickel(II) complexes (NiSAL-MeDPT and NiSAL-HDPT), which, for a variety of reasons, turned out to be perfectly suited to challenge the predictions to the finest level of detail. For NiSAL-MeDPT, new NMR experiments yield an assignment that perfectly matches the calculations. The slightly different hyperfine shifts from the two “halves” of the molecules related by a pseudo-C₂ axis, which are experimentally divided into two well-defined spin systems, are also straightforwardly distinguished by the calculations. In the case of NiSAL-HDPT, for which no X-ray structure is available, the quality of the calculations allowed us to refine its structure using as a starting template the structure of NiSAL-MeDPT.

State-of-the-art quantum chemical methods and paramagnetism-tailored NMR experiments provide a deep insight on the relation between the spectra and the electronic structure for two paramagnetic pentacoordinate nickel(II) complexes.