Interplay of Intrinsic, Environmental, and Dynamic Effects in Tuning the EPR Parameters of Nitroxides: Further Insights from an Integrated Computational Approach

Profound Changes in Functional Structure and Dynamics of Serum Albumin in Children with Nephrotic Syndrome: An Exploratory Research Study

Patients with nephrotic syndrome (NS) suffer from urinary loss of albumin. As a cause, previous studies focused on the glomerular filter rather than analyzing the molecular properties of albumin itself. Later one was initiated by clinical observations indicating unexplained molecular alterations of human serum albumin (HSA) in an NS pediatric patient. Therefore, we examined serum from eight pediatric patients with steroid-sensitive and -resistant NS and compared it with serum from healthy subjects as well as commercial HSA. We used dynamic and electrophoretic light scattering to characterize the protein size and effective surface charge and electron paramagnetic resonance spectroscopy to measure the local environment and binding dynamics of up to seven fatty acids associated with HSA. Our findings suggest that pronounced differences in binding behavior and surface charge of HSA could enhance their filtration through the GBM, leading to direct toxicity of HSA to podocytes.

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.

Electronic and Vibrational Manifold of Tetracyanoethylene-Chloronaphthalene Charge Transfer Complex in Solution: Insights from TD-DFT and Ab Initio Molecular Dynamics

The interplay between light absorption and the molecular environment has a central role in the observed photophysics of a wide range of photoinduced chemical and biological phenomena. The understanding of the interplay between vibrational and electronic transitions is the focus of this work, since it can provide a rationale to tune the optical properties of charge transfer (CT) materials used for technological applications. A clear description of these processes poses a nontrivial challenge from both the theoretical and experimental points of view, where the main issue is how to accurately describe and probe drastic changes in the electronic structure and the ultrafast molecular relaxation and dynamics. In this work we focused on the intermolecular CT reaction that occurs upon photon absorption in a π-stacked model system in dichloromethane solution, in which the 1-chloronaphthalene (1ClN) acts as the electron donor and tetracyanoethylene (TCNE) is the electron acceptor. Density functional theory calculations have been carried out to characterize both the ground-state properties and more importantly the low-lying CT electronic transition, and excellent agreement with recently available experimental results [Mathies, R. A.; et al. J. Phys. Chem. A 2018, 122 (14), 3594] was obtained. The minima of the ground state and first singlet excited state have been accurately characterized in terms of spatial arrangements and vibrational Raman frequencies, and the CT natures of the first two low-lying electronic transitions in the absorption spectra have been addressed and clarified too. Finally, by modeling the possible coordination sites of the TCNE electron acceptor with respect to monovalent ions (Na+, K+) in an implicit solution of acetonitrile, we find that TCNE can accommodate a counterion in two different arrangements, parallel and orthogonal to the C═C axis, leading to the formation of a contact ion pair. The nature of the counterion and its relative position entail structural modifications of the TCNE radical anion, mainly the central C═C and C≡N bonds, compared to the isolated case. An important red shift of the C═C stretching frequency was observed when the counterion is orthogonal to the double bond, to a greater extent for Na+. On the contrary, in the second case, where the counterion ion lies along the internuclear C═C axis, we find that K+ polarizes the electron density of the double bond more, resulting in a greater red shift than with Na+.

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.

o-Semiquinone radical anion isolated as an amorphous porous solid

The use of metal cations is a commonly applied strategy to create S > 1/2 stable molecular systems containing semiquinone radicals. Persistent mono-semiquinonato complexes of diamagnetic metal ions (S = 1/2) have been hitherto less common and mostly limited to the complexes of heavy metal ions. In this work, a mono-semiquinonato complex of aluminum, derived from 1,2-dihydroxybenzene, is obtained using a surprisingly short and uncomplicated procedure. The isolated product is an amorphous and porous solid that exhibits very good stability under ambient conditions. To characterise its molecular and electronic structure, 9.7, 34 and 406 GHz EPR spectroscopy was used in concert with computational techniques (DFT and DLPNO-CCSD). It was revealed that the radical complex is composed of two chemically equivalent aluminum cations and two catechol-like ligands with the unpaired electron uniformly distributed between the two organic molecules. The good stability and porous structure make this complex applicable in heterogeneous aerobic reactions.

How accurate is density functional theory in predicting spin density? An insight from the prediction of hyperfine coupling constants

Electron paramagnetic resonance (EPR) spectroscopy has been proven to be an important technique for studying paramagnetic systems. Probably, the most accessible EPR parameter and the one that provides a significant amount of information about molecular structure and spin density is the hyperfine coupling constant (HFCC). Hence, accurate quantum-chemical modeling of HFCCs is frequently essential to the adequate interpretation of EPR spectra. It requires the precise spin density, which is the difference between the densities of α- and β-electrons, and thus, its quality is expected to reflect the quality of the total electron density. The question of which approximate exchange-correlation density functional yields sufficiently accurate HFCCs, and thus, the spin density remains open. To assess the performance of well-established density functionals for calculating HFCCs, we used a series of 26 small paramagnetic species and compared the obtained results to the CCSD reference values. The performance of DFT was also tested on EPR-studied o-semiquinone radical interacting with water molecules and Mg²⁺ cation. The HFCCs were additionally calculated by the DLPNO-CCSD method, and this wave function-based technique was found superior to all functionals we tested. Although some functionals were found, on average, to be fairly efficient, we found that the most accurate functional is system-dependent, and therefore, the DLPNO-CCSD method should be preferred for theoretical investigations of the HFCCs and spin density.

Towards reliable references for electron paramagnetic resonance parameters based on quantum chemistry: the case of verdazyl radicals

We present an efficient and accurate computational procedure to calculate properties measurable by EPR spectroscopy. We simulate a molecular dynamics (MD) trajectory by employing the quantum mechanically derived force field (QMDFF) [S. Grimme, J. Chem. Theory Comput., 2014, 10, 4497] and sample the trajectory at different time steps. For each snapshot EPR properties are calculated with a hybrid density functional theory (DFT) method. EPR spectra are simulated based on the averaged results. We applied the strategy to a number of previously published and novel verdazyl radicals, for which we recorded EPR spectra. The resulting simulated spectra are compatible with experiment already before employing an additional fitting step, in contrast to those from single point electronic-structure calculations. After the refinement, the experimental data are excellently reproduced, and the fitted EPR parameters do not deviate much from the calculated ones. This provides confidence in ascribing a direct physical meaning to the refined data in terms of experimental EPR parameters rather than merely considering them as mathematical fit parameters. We also find that couplings to hydrogen nuclei have a significant influence on the spectra of verdazyl radicals.

Measuring Spin⋅⋅⋅Spin Interactions between Heterospins in a Hybrid [2]Rotaxane

Use of molecular electron spins as qubits for quantum computing will depend on the ability to produce molecules with weak but measurable interactions between the qubits. Here we demonstrate use of pulsed EPR spectroscopy to measure the interaction between two inequivalent spins in a hybrid rotaxane molecule.

Spin-labelled diketopiperazines and peptide-peptoid chimera by Ugi-multi-component-reactions

For the first time, spin-labelled coumpounds have been obtained by isonitrile-based multi component reactions (IMCRs). The typical IMCR Ugi-protocols offer a simple experimental setup allowing structural variety by which labelled diketopiperazines (DKPs) and peptide-peptoid chimera have been synthesized. The reaction keeps the paramagnetic spin label intact and offers a simple and versatile route to a large variety of new and chemically diverse spin labels.

Density functional theory studies of MTSL nitroxide side chain conformations attached to an activation loop

A quantum mechanical (QM) method rooted on density functional theory (DFT) has been employed to determine conformations of the methane-thiosulfonate spin label (MTSL) attached to a fragment extracted from the activation loop of Aurora-A kinase. The features of the calculated energy surface revealed low energy barriers between isoenergetic minima, and the system could be described in a population of 76 rotamers that can be also considered for other systems since it was found that the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \chi_{3} $$\end{document}χ3, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \chi_{4} $$\end{document}χ4 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \chi_{5} $$\end{document}χ5 do not depend on the previous two dihedral angles. Conformational states obtained were seen to be comparable to those obtained in the α-helix systems studied previously, indicating that the protein backbone does not affect the torsional profiles significantly and suggesting the possibility to use determined conformations for other protein systems for further modelling studies.

Carbon and proton Overhauser DNP from MD simulations and ab initio calculations: TEMPOL in acetone

A computational analysis of the Overhauser effect is reported for the proton, methyl carbon, and carbonyl carbon nuclei of liquid acetone doped with the nitroxide radical TEMPOL. A practical methodology for calculating the dynamic nuclear polarization (DNP) coupling factors by accounting for both dipole-dipole and Fermi-contact interactions is presented. The contribution to the dipolar spectral density function of nuclear spins that are not too far from TEMPOL is computed through classical molecular dynamics (MD) simulations, whereas the contribution of distant spins is included analytically. Fermi contacts are obtained by subjecting a few molecules from every MD snapshot to ab initio quantum mechanical calculations. Scalar interaction is found to be an essential part of the (13)C Overhauser DNP. While mostly detrimental to the carbonyl carbon of acetone it is predicted to result in large enhancements of the methyl carbon signal at magnetic fields of 9 T and beyond. In contrast, scalar coupling is shown to be negligible for the protons of acetone. The additional influence of proton polarization on the carbon DNP (three-spin effect) is also analyzed computationally. Its effect, however, is concluded to be practically insignificant for liquid acetone.

Hydrogen bonding of nitroxide spin labels in membrane proteins

On the basis of experiments at 275 GHz, we reconsider the dependence of the continuous-wave EPR spectra of nitroxide spin-labeled protein sites in sensory- and bacteriorhodopsin on the micro-environment.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

First principles calculation of inhomogeneous broadening in solid-state cw-EPR spectroscopy

We present a scheme for the first-principles calculation of EPR lineshapes for continuous-wave-EPR spectroscopy (cw-EPR) of spin centers in complex chemical environments. We specifically focus on poorly characterized systems, e.g. powders and frozen glasses with variable microsolvation structures. Our approach is based on ab initio molecular dynamics simulations and ab initio calculations of the ensemble of g- and A-tensors along the trajectory. The method incorporates temperature effects as well as the full anharmonicity of the intra- and intermolecular degrees of freedom of the system. We apply this scheme to compute the lineshape of a prototypical spin probe, the nitrosodisulfonate dianionic radical (Fremy's salt), dissolved in a 50 : 50 mixture of water and methanol. We are able to determine the specific effect of variations of local solvent composition and microsolvation structure on the cw-EPR lineshape. Our molecular dynamics reveal a highly anisotropic solvation structure with distinct spatial preferences for water and methanol around Fremy's salt that can be traced back to a combination of steric and polar influences. The overall solvation structure and conformational preferences of Fremy's salt as found in our MD simulations agree very well with the results obtained from EPR and orientation-selective ENDOR spectroscopy performed on the frozen glass. The simulated EPR lineshapes show good agreement with the experimental spectra. When combined with our MD results, they characterize the lineshape dependence on local morphological fluctuations.

Extension of the AMBER Force Field for Nitroxide Radicals and Combined QM/MM/PCM Approach to the Accurate Determination of EPR Parameters of DMPO-H in Solution

A computational strategy that combines both time-dependent and time-independent approaches is exploited to accurately model molecular dynamics and solvent effects on the isotropic hyperfine coupling constants of the DMPO-H nitroxide. Our recent general force field for nitroxides derived from AMBER ff99SB is further extended to systems involving hydrogen atoms in β-positions with respect to NO-moiety. The resulting force-field has been employed in a series of classical molecular dynamics simulations, comparing the computed EPR parameters from selected molecular configurations to the corresponding experimental data in different solvents. The effect of vibrational averaging on the spectroscopic parameters is also taken into account, by second-order vibrational perturbation theory involving semidiagonal third energy derivatives together first and second property derivatives.

π-Stacking effects on the EPR parameters of a prototypical DNA spin label

The character and value of spin labels for probing environments like double-stranded DNA depend on the degree of change in the spin Hamiltonian parameters of the spin label induced by the environment. Herein we provide a systematic theoretical investigation of this issue, based on a density functional theory method applied to a spin labeled DNA model system, focusing on the dependence of the EPR properties of the spin label on the π stacking and hydrogen bonding that occur upon incorporating the spin label into the selected base pair inside DNA. It is found that the EPR spin Hamiltonian parameters of the spin label are only negligibly affected by its incorporation into DNA, when compared to its free form. This result gives a theoretical ground for the common empirical assumption regarding the behaviour of spin Hamiltonian parameters made in EPR based measurements of the distance between spin labels incorporated into DNA.

First determination of the spin relaxation properties of a nitronyl nitroxide in solution by electron spin echoes at X-band: a comparison with Tempone

We studied by electron spin echo pulse methods the spin relaxation properties of a phenyl nitronyl nitroxide radical (2-phenyl-4,4,5,5-tetramethylimidazoline-1-oxyl 3-oxide, PTIO) at X-band in fluid toluene solution in a wide temperature range, and in a water/glycerol 1:1 mixture near room temperature. The relaxation properties of PTIO have been compared with that of Tempone, as a widely used nitroxide. By a new procedure, based on experimental results on the temperature dependences of the relaxation times T(1) and T(2), and on the approximation of an isotropic brownian rotational diffusion, we separated non-secular, spin rotational and residual terms from the transverse relaxation rate to isolate secular and pseudosecular contributions. By comparing the results for the two radicals we found the differences in the magnetic properties that give rise to slower transverse (T(2)) and longitudinal (T(1)) electron spin relaxation for PTIO in the whole temperature range explored in this work.

Toward spectroscopic accuracy for open-shell systems: molecular structure and hyperfine coupling constants of H2CN, H2CP, NH2, and PH2 as test cases

In the present paper, we investigate the molecular structure and hyperfine couplings of a series of σ radicals containing first- and second-row atoms (H(2)CN,H(2)CP,NH(2),PH(2)) for which accurate gas-phase microwave results are available. The presence of α- and, especially, β-hydrogen atoms makes the evaluation of magnetic properties of these radicals particularly challenging. Geometrical parameters have been computed by the coupled-cluster ansatz in conjunction with hierarchical series of basis sets, thus accounting for extrapolation to the complete basis-set limit. Core correlation as well as higher excitations in the electronic-correlation treatment have also been taken into account. An analogous approach has been employed for evaluating hyperfine coupling constants with particular emphasis given to basis-set, correlation, and geometrical effects. The corresponding vibrational corrections, required for a meaningful comparison to experimental data, have also been investigated. The remarkable agreement with experiment confirms the reliability of the present computational approach, already validated for π radicals, thus establishing the way for setting up a benchmark database for magnetic properties.

A general approach for prediction of motional EPR spectra from Molecular Dynamics (MD) simulations: application to spin labelled protein

A general approach for the prediction of EPR spectra directly and completely from single dynamical trajectories generated from Molecular Dynamics (MD) simulations is described. The approach is applicable to an arbitrary system of electron and nuclear spins described by a general form of the spin-Hamiltonian for the entire motional range. It is shown that for a reliable simulation of motional EPR spectra only a single truncated dynamical trajectory generated until the point when correlation functions of rotational dynamics are completely relaxed is required. The simulation algorithm is based on a combination of the propagation of the spin density matrix in the Liouville space for this initial time interval and the use of well defined parameters calculated entirely from the dynamical trajectory for prediction of the evolution of the spin density matrix at longer times. A new approach is illustrated with the application to a nitroxide spin label MTSL attached to the protein sperm whale myoglobin. It is shown that simulation of the EPR spectrum, which is in excellent agreement with experiment, can be achieved from a single MD trajectory. Calculations reveal the complex nature of the dynamics of a spin label which is a superposition of the fast librational motions within dihedral states, of slow rotameric dynamics among different conformational states of the nitroxide tether and of the slow rotational diffusion of the protein itself. The significance of the slow rotameric dynamics of the nitroxide tether on the overall shape of the EPR spectrum is analysed and discussed.

DFT Calculations of Isotropic Hyperfine Coupling Constants of Nitrogen Aromatic Radicals: The Challenge of Nitroxide Radicals

The performance of DFT methodology to predict with accuracy the isotropic hyperfine coupling constants (hfccs) of aromatic radicals containing (14)N nucleus is investigated by an extensive study in which 165 hfccs, belonging to 38 radical species, are obtained from calculations with B3LYP and PBE0 functionals combined with 6-31G*, N07D, TZVP, and EPR-III basis sets, and are compared to the reported experimental data. The results indicate that the selection of the basis set is of fundamental importance in the calculation of (14)N hfccs, whereas there is not so great an influence on the accurate computation of that parameter for (1)H nuclei. The values of the calculated (14)N coupling constants of aromatic nitroxide radicals using DFT methodology are noticeably lower than the experimental ones. A very simple relation to predict these hfccs with high accuracy is proposed on the basis of the present results, as an interesting alternative to the highly computationally demanding integrated approaches so far used.

Environmental effects in computational spectroscopy: accuracy and interpretation

Spectroscopic techniques are valuable tools for understanding the structure and dynamics of complex systems, such as biomolecules or nanomaterials. Most of the current research is devoted to the development of new experimental techniques for improving the intrinsic resolution of different spectra. However, the subtle interplay of several different effects acting at different length and time scales still makes the interpretation and analysis of such spectra a very difficult task. In this respect, computational spectroscopy is becoming a needful and versatile tool for the assignment and interpretation of experimental spectra. It is in fact possible nowadays to model with relatively high accuracy the physical-chemical properties of complex molecules in different environments, and to link spectroscopic evidence directly to the structural and dynamical properties of optically or magnetically active solvated probes. In this Review, significant steps toward the simulation of entire spectra in condensed phases are presented together with some basic aspects of computational spectroscopy, which highlight how intramolecular and intermolecular degrees of freedom influence several spectroscopic parameters.

DFT study of nitroxide radicals: explicit modeling of solvent effects on the structural and electronic characteristics of 4-amino-2,2,6,6-tetramethyl-piperidine-N-oxyl

An explicit DFT modeling of water surroundings on the electron paramagnetic resonance properties of 4-amino-2,2,6,6-tetramethyl-piperidine-N-oxyl (TA) has been performed. A stepwise hydration of TA is accompanied with certain changes in geometrical parameters (bond lengths and angles) and redistribution of partial electric charges in TA. An aqueous cluster of 45 water molecules can be considered as an appropriate model for a complete aqueous shell around TA, although most of the structural and electronic characteristics of TA already converge at about 10 water molecules. Water surroundings induce an increase in electron spin density on the nitrogen atom of the nitroxide fragment due to stabilization of the polar resonance structure > N(+*)-O(-) at the expense of less polar structure > N-O*. The water-induced rise of the isotropic splitting constant a(iso), calculated from the contact term of the hyperfine interaction, comprises Deltaa(iso)(rho(N2)) = 2.2-2.5 G, which is typical of experimental value for TA. There are two contributions to the solvent effect on the a(iso)(rho(N2)) value: the redistribution of spin density in the nitroxide fragment (polarity effect) and water-induced distortions of TA geometry. Microscopic variations in a hydrogen-bonded water network cause noticeable fluctuations of the splitting constant a(iso)(rho(N2)). Calculations of the atomic spin density (sigma(N2)) allowed us to compute the splitting constant from the relationship a(iso)(sigma(N2)) = Qsigma(N2), where Q = 36.2 G. A practical advantage of using this relationship is that it gives 'smoothed' values of the splitting constant, which are sensitive to the environment polarity but remain tolerant to microscopic fluctuations of the hydrogen-bonded water network around a spin-label molecule.

Overhauser dynamic nuclear polarization and molecular dynamics simulations using pyrroline and piperidine ring nitroxide radicals

The efficiency of Overhauser dynamic nuclear polarization (DNP) depends on the local dynamics modulating the dipolar coupling between the two interacting spins. By attaching nitroxide based spin labels to molecules and by measuring the (1)H DNP response of solvent water, information about the local hydration dynamics near the spin label can be obtained. However, there are two commonly used types of nitroxide ring structures; a pyrroline based and a piperidine based molecule. It is important to know when comparing different experiments, whether changes in DNP enhancements are due to changes in local hydration dynamics or because of the different spin label structures. In this study we investigate the key parameters affecting DNP signal enhancements for 3-carbamoyl-2,2,5,5-tetramethyl-3-pyrrolin-1-oxyl, a 5-membered ring nitroxide radical, and for 4-oxo-2,2,6,6-tetramethyl-1-piperidinyloxy, a 6-membered ring nitroxide radical. Using X-Band DNP, field cycling relaxometry, and molecular dynamics simulations, we conclude that the key parameters affecting the DNP amplitude of the (1)H signal of water to be equal when using either nitroxide. Thus, experiments measuring hydration dynamics using either type of spin labels may be compared.

Dynamic nuclear polarization coupling factors calculated from molecular dynamics simulations of a nitroxide radical in water

The magnetic resonance signal obtained from nuclear spins is strongly affected by the presence of nearby electronic spins. This effect finds application in biomedical imaging and structural characterization of large biomolecules. In many of these applications nitroxide free radicals are widely used due to their non-toxicity and versatility as site-specific spin labels. We perform molecular dynamics simulations to study the electron-nucleus interaction of the nitroxide radical TEMPOL and water in atomistic detail. Correlation functions corresponding to the dipolar and scalar spin-spin couplings are computed from the simulations. The dynamic nuclear polarization coupling factors deduced from these correlation functions are in good agreement with experiment over a broad range of magnetic field strengths. The present approach can be applied to study solute-solvent interactions in general, and to characterize solvent dynamics on the surfaces of proteins or other spin-labeled biomolecules in particular.