Important components for accurate hyperfine coupling constants: electron correlation, dynamic contributions, and solvation effects

Systematic QM/MM Study for Predicting 31P NMR Chemical Shifts of Adenosine Nucleotides in Solution and Stages of ATP Hydrolysis in a Protein Environment

NMR (nuclear magnetic resonance) spectroscopy allows for important atomistic insights into the structure and dynamics of biological macromolecules; however, reliable assignments of experimental spectra are often difficult. Herein, quantum mechanical/molecular mechanical (QM/MM) calculations can provide crucial support. A major problem for the simulations is that experimental NMR signals are time-averaged over much longer time scales, and since computed chemical shifts are highly sensitive to local changes in the electronic and structural environment, sufficiently large averages over representative structural ensembles are essential. This entails high computational demands for reliable simulations. For NMR measurements in biological systems, a nucleus of major interest is 31P since it is both highly present (e.g., in nucleic acids) and easily observable. The focus of our present study is to develop a robust and computationally cost-efficient framework for simulating 31P NMR chemical shifts of nucleotides. We apply this scheme to study the different stages of the ATP hydrolysis reaction catalyzed by p97. Our methodology is based on MM molecular dynamics (MM-MD) sampling, followed by QM/MM structure optimizations and NMR calculations. Overall, our study is one of the most comprehensive QM-based 31P studies in a protein environment and the first to provide computed NMR chemical shifts for multiple nucleotide states in a protein environment. This study sheds light on a process that is challenging to probe experimentally and aims to bridge the gap between measured and calculated NMR spectroscopic properties.

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.

Tensor-Hypercontracted MP2 First Derivatives: Runtime and Memory Efficient Computation of Hyperfine Coupling Constants

We employ our recently introduced tensor-hypercontracted (THC) second-order Møller-Plesset perturbation theory (MP2) method [Bangerter, F. H., Glasbrenner, M., Ochsenfeld, C. J. Chem. Theory Comput. 2021, 17, 211-221] for the computation of hyperfine coupling constants (HFCCs). The implementation leverages the tensor structure of the THC factorized electron repulsion integrals for an efficient formation of the integral-based intermediates. The computational complexity of the most expensive and formally quintic scaling exchange-like contribution is reduced to effectively subquadratic, by making use of the intrinsic, exponentially decaying coupling between tensor indices through screening based on natural blocking. Overall, this yields an effective subquadratic scaling with a low prefactor for the presented THC-based AO-MP2 method for the computation of isotropic HFCCs on DNA fragments with up to 500 atoms and 5000 basis functions. Furthermore, the implementation achieves considerable speedups with up to a factor of roughly 600-1000 compared to previous implementations [Vogler, S., Ludwig, M., Maurer, M., Ochsenfeld, C. J. Chem. Phys. 2017, 147, 024101] for medium-sized organic radicals, while also significantly reducing storage requirements.

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.