Car-Parrinello molecular dynamics simulations and EPR property calculations on aqueous ubisemiquinone radical anion

Car-Parrinello molecular dynamics (CPMD) simulations have been performed on ubisemiquinone radical anion in aqueous solution. The different types of hydrogen bonding formed between the semiquinone and the solvent were studied in terms of frequency and directionality, in comparison with the parent benzosemiquinone radical anion. The EPR parameters (g-tensors and hyperfine coupling constants) were obtained from quantum chemical property calculations performed on snapshots along the MD trajectory, and the contributions of different solvation effects to the EPR parameters have been evaluated. The influence of the anion’s conformational behaviour was examined, including the orientation-dependent effects of hyperconjugation on side-chain hyperfine coupling.     

Phylloquinone and related radical anions studied by pulse electron nuclear double resonance spectroscopy at 34 GHz and density functional theory

1H hyperfine (hf) coupling constants of semiquinone radical anions of 1,4-naphthoquinone, 2-methyl-1,4-naphthoquinone, and 2-methyl-3-phytyl-1,4-naphthoquinone in frozen alcoholic solutions were measured using pulse Q-band electron nuclear double resonance spectroscopy. The resolved signals of the quinone protons as well as from hydrogen bond and solvent shell protons were analyzed and assigned. Both in-plane and out-of-plane hydrogen bonding with respect to the pi-plane of the radical is observed. Interactions with nonexchangeable protons from the surrounding matrix are detected and assigned to solvent protons above and below the quinone plane. Density functional theory was used to calculate spin Hamiltonian parameters of the radical anions. Solvent molecules of the first solvent shell that provide hydrogen bonds to the quinones were included in the geometry optimization. The conductor-like screening model was employed to introduce additional effects of the solvent cage. From a comparison of the experimental and calculated hf tensors it is concluded that four solvent molecules are coordinated via hydrogen bonds to the quinone oxygens. For all radicals very good agreement between experimental and calculated data is observed. The influence of different substituents on the spin density distribution and hydrogen bond geometries is discussed.     

Protein-cofactor interactions and EPR parameters for the Q(H) quinone binding site of quinol oxidase. A density functional study

Recent multifrequency EPR studies of the "high-affinity" quinone binding site of quinol oxidase (Q(H) site) have suggested a very asymmetric hydrogen-bonding environment for the semiquinone radical anion state. Single-sided hydrogen bonding to the O1 carbonyl position was one of the proposals, which contrasts with some previous experimental indications. Here density functional calculations of the EPR parameters (g-tensors, 13C, 1H, and 17O hyperfine tensors) for a wide variety of supermolecular model complexes have been used to provide insight into the detailed relations among structure, environment, and EPR parameters of ubisemiquinone radical anions. A single-sided binding model is not able to account for the experimentally observed low g(x) component of the g-tensor or for the observed magnitude of the asymmetry of the 13C carbonyl HFC tensors. Based on the detailed comparison between computation and experiment, a model with two hydrogen bonds to O1 and one hydrogen bond to O4 is suggested for the Q(H) site, but a model with one more hydrogen bond on each side cannot be excluded. Several general conclusions on the interrelations between EPR parameters and hydrogen bond patterns of ubisemiquinones in proteins are provided.     

Study of rhamnose radicals in the solid state adopting a density functional theory cluster approach

A theoretical study is performed on the radiation-induced radicals in crystalline alpha-l-rhamnose, using density functional theory (DFT) calculations. Irrespective of earlier structural assignments, a host of possible radical models is examined in search for a structure that accurately reproduces experimental electron paramagnetic resonance (EPR) properties. A cluster approach is followed, incorporating all hydrogen bond interactions between radical and crystalline environment. Hyperfine coupling tensors as well as g tensors are determined and a comparison is made with available experimental data. Three carbon-centered hydroxyalkyl radicals are validated, in accordance with experimental suggestions for their structure. The occurrence of a carbon-centered oxygen anion radical for one of the radical species is rejected on theoretical grounds, and instead an altered hydroxyalkyl structure is suggested. Our cluster calculations are able to determine g and hyperfine tensors for the oxygen-centered alkoxy radical in rhamnose, in accordance with one of the two measurements for this species. For all radical models, quantitative agreement with experimental hyperfine tensors is obtained by performing full cluster DFT calculations. The inclusion of the molecular environment for the determination of this EPR property proved to be essential.     

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.     

Electron magnetic resonance and density functional theory study of room temperature X-irradiated β-D-fructose single crystals

Stable free radical formation in fructose single crystals X-irradiated at room temperature was investigated using Q-band electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR), and ENDOR induced EPR (EIE) techniques. ENDOR angular variations in the three main crystallographic planes allowed an unambiguous determination of 12 proton HFC tensors. From the EIE studies, these hyperfine interactions were assigned to six different radical species, labeled F1-F6. Two of the radicals (F1 and F2) were studied previously by Vanhaelewyn et al. [Vanhaelewyn, G. C. A. M.; Pauwels, E.; Callens, F. J.; Waroquier, M.; Sagstuen, E.; Matthys, P. J. Phys. Chem. A 2006, 110, 2147.] and Tarpan et al. [Tarpan, M. A.; Vrielinck, H.; De Cooman, H.; Callens, F. J. J. Phys. Chem. A 2009, 113, 7994.]. The other four radicals are reported here for the first time and periodic density functional theory (DFT) calculations were used to aid their structural identification. For the radical F3 a C3 carbon centered radical with a carbonyl group at the C4 position is proposed. The close similarity in HFC tensors suggests that F4 and F5 originate from the same type of radical stabilized in two slightly different conformations. For these radicals a C2 carbon centered radical model with a carbonyl group situated at the C3 position is proposed. A rather exotic C2 centered radical model is proposed for F6.     

Accurate Backbone 13C and 15N Chemical Shift Tensors in Galectin-3 Determined by MAS NMR and QM/MM: Details of Structure and Environment Matter

Chemical shift tensors obtained from solid-state NMR spectroscopy are very sensitive reporters of structure and dynamics in proteins. While accurate ¹³C and ¹⁵N chemical shift tensors are accessible by magic angle spinning (MAS) NMR, their quantum mechanical calculations remain challenging, particularly for ¹⁵N atoms. Here we compare experimentally determined backbone ¹³C^α and ¹⁵N^H chemical shift tensors by MAS NMR with hybrid quantum mechanics/molecular mechanics/molecular dynamics (MD-QM/MM) calculations for the carbohydrate-binding domain of galectin-3. Excellent agreement between experimental and computed ¹⁵N^H chemical shift anisotropy values was obtained using the Amber ff15ipq force field when solvent dynamics was taken into account in the calculation. Our results establish important benchmark conditions for improving the accuracy of chemical shift calculations in proteins and may aid in the validation of protein structure models derived by MAS NMR.     

Thermal and solvent effects on NMR indirect spin-spin coupling constants of a prototypical Chagas disease drug

NMR J-couplings across hydrogen bonds reflect the static and dynamic character of hydrogen bonding. They are affected by thermal and solvent effects and can therefore be used to probe such effects. We have applied density functional theory (DFT) to compute the NMR (n)J(N,H) scalar couplings of a prototypical Chagas disease drug (metronidazole). The calculations were done for the molecule in vacuo, in microsolvated cluster models with one or few water molecules, in snapshots obtained from molecular dynamics simulations with explicit water solvent, and in a polarizable dielectric continuum. Hyperconjugative and electrostatic effects on spin-spin coupling constants were assessed through DFT calculations using natural bond orbital (NBO) analysis and atoms in molecules (AIM) theory. In the calculations with explicit solvent molecules, special attention was given to the nature of the hydrogen bonds formed with the solvent molecules. The results highlight the importance of properly incorporating thermal and solvent effects into NMR calculations in the condensed phase.     

Molecular mechanical devices based on quinone-pyrrole and quinone-indole dyads: a computational study

A set of intramolecularly connected dyads consisting of a quinone unit and a pyrrole or indole moiety have been designed and evaluated in quantum-chemical calculations. It is shown computationally for several systems, depending on the length and attachment points of the interconnecting chains, that a reduction of the quinone to the semiquinone radical anion or quinolate dianion state leads to a reversible intramolecular reorientation from a pi-stacked to a T-stacked arrangement. In the rearranged structures, a hydrogen bond from the pyrrole or indole N-H function to the semiquinone or quinolate pi-system is created upon reduction. In some systems, hydrogen bonds to the semiquinone or quinolate oxygen atoms are partly feasible and will be preferred over T-stacking. The choice of systems has been based on recent computational observations related to photosystem I. Systems with pyrrole or indole units should provide a better basis for the envisioned molecular motor than recently proposed quinone-benzene dyads. The intramolecular interactions modify the quinone redox potentials. Electronic g-tensors have been computed for the semiquinone states. These reflect characteristically the presence and nature of hydrogen bonds to the semiquinone and represent suitable electron paramagnetic resonance spectroscopic probes for the preferred structures. Intramolecular proton transfer is possible in the dianionic state.     

Theoretical determination of standard oxidation and reduction potentials of chlorophyll-a in acetonitrile

QM/MM calculations were performed on ethyl chlorophyllide-a and its radical cation and anion, by using the density functional (DF) B3LYP method to determine the molecular characteristics, and a molecular mechanics (MM) method to simulate the solvating medium. The presence of the solvent was accounted for during the optimization of the geometry of the 85-atom chlorophyll-a system by using an ONIOM methodology. A total of 24 solvent molecules were explicitly considered during the optimization process, and these were treated by the universal force field (UFF) method. Initially, the split-valence 3-21G basis set was used for optimizing the geometry of the 85-atom species, neutral, cation and anion. Electronic energies were then determined for the optimized species by making use of the polarized 6-31G(d) basis set. The ionization energy calculated (6.0 eV) is in very good agreement with the observed one (6.1 eV). The MM+ force field was used to investigate the dynamics of the acetonitrile molecules around the neutral species as well as the radical ions of chlorophyll. The required atomic charges on all the atoms were obtained from calculations on all involved molecules at the DFT/6-31G(d) level. Randomly sampled configurations were used to determine the first solvation layer contribution to the free energy of solvation of various species. A truncated 46-atom model of ethyl chlorophyllide-a was used to evaluate the thermal energies of neutral chlorophyll molecule relative to its two radical ions in the gas phase. Born energy, Onsager energy, and the Debye-Huckel energy of the chlorophyll-solvent aggregate were added as perturbative corrections to the free energy of solvation that was initially obtained through molecular dynamics method for the same complex. These calculations yield the oxidation potential as 0.75 +/- 0.32 V and the reduction potential -1.18 +/- 0.31 V at 298.15 K. The calculated values are in good agreement with the experimental midpoint potentials of +0.76 and -1.04 V, respectively.     

Simulations of structure and vibrational spectra of deoxyoctanucleotides

Combined molecular dynamics and ab initio computations were applied for analysis of infrared absorption and vibrational circular dichroism spectra of deoxyoctanucleotides. Unlike for previous idealized models, molecular geometries of these shorter DNA fragments in solutions were obtained as dynamic averages from simulations in a periodic water box. Vibrational spectra for the whole octamers including hydrogen-bonded solvent molecules were simulated on the basis of density-functional calculations on small fragments and subsequent transfer of molecular property tensors. Explicit and continuum solvent models were compared. Apparently, the DNA segments retain an approximate B-conformation in the aqueous solutions, but the terminal base pairs significantly deviate from the planar arrangement and the vibrational circular dichroism spectrum for (CG)(2) nucleotide indicates a larger average helical twist. Sodium counterions moved freely around the molecule during the simulation and do not influence spectral intensities. Simulated absorption spectra faithfully reproduced the experimental signal of principal functional groups, while only qualitative agreement was obtained for the dependence on the basis sequence.     

Radiation products at 77 K in trehalose single crystals: EMR and DFT analysis

The radicals obtained in trehalose dihydrate single crystals after 77 K X-irradiation have been investigated at the same temperature using X-band electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR), and ENDOR-induced EPR (EIE) techniques. Five proton hyperfine coupling tensors were unambiguously determined from the ENDOR measurements and assigned to three carbon-centered radical species (T1, T1*, and T2) based on the EIE spectra. EPR angular variations revealed the presence of four additional alkoxy radical species (T3 to T6) and allowed determination of their g tensors. Using periodic density functional theory (DFT) calculations, T1/T1*, T2, and T3 were identified as H-loss species centered at C4, C1', and O2', respectively. The T4 radical is proposed to have the unpaired electron at O4, but considerable discrepancies between experimental and calculated HFC values indicate it is not simply the (net) H-loss species. No suitable models were found for T5 and T6. These exhibit a markedly larger g anisotropy than T3 and T4, which were not reproduced by any of our DFT calculations.     

Hydration of hydrogentungstate anions at different ph conditions: a car-parrinello molecular dynamics study

Standard density functional theory calculations with a continuous model of solvation as well as Car-Parrinello molecular dynamics simulations with explicit solvent molecules are carried out to analyze the effect of the pH of the solution on the coordination sphere of the W (VI) ion. Both methodologies agree in predicting an expansion of the coordination sphere of the W (VI) ion upon a decrease in the pH. Continuous solvation models, however, are unable to predict as stable some structural isomers of a hydrated hydrogentungstate anion and tungstic acid.     

Reaction of superoxide radical with quinone molecules

When the superoxide radical O(2)(?-) is generated on reaction of KO(2) with water in dimethyl sulfoxide, the decay of the radical is dramatically accelerated by inclusion of quinones in the reaction mix. For quinones with no or short hydrophobic tails, the radical product is a semiquinone at much lower yield, likely indicating reduction of quinone by superoxide and loss of most of the semiquinone product by disproportionation. In the presence of ubiquinone-10, a different species (I) is generated, which has the EPR spectrum of superoxide radical. However, pulsed EPR shows spin interaction with protons in fully deuterated solvent, indicating close proximity to the ubinquinone-10. We discuss the nature of species I, and possible roles in the physiological reactions through which ubisemiquinone generates superoxide by reduction of O(2) through bypass reactions in electron transfer chains.     

Solvation of N3- at the water surface: the polarizable continuum model approach

We present a new quantum mechanical model to introduce Pauli repulsion interaction between a molecular solute and the surrounding solvent in the framework of the Polarizable Continuum Model. The new expression is derived in a way to allow naturally for a position-dependent solvent density. This development makes it possible to employ the derived expression for the calculation of molecular properties at the interface between two different dielectrics. The new formulation has been tested on the azide anion (N3-) for which we have calculated the solvation energy, the dipole moment, and the static polarizability at the interface as a function of the ion position. The calculations have been carried out for different ion-surface orientations, and the results have also been compared with the parallel electrostatic-only solvation model.     

An ONIOM study of the QA site semiquinone in the Rhodobacter sphaeroides photosynthetic reaction centre

ONIOM (QM/MM) calculations are performed to investigate the spin density distribution for the ubisemiquinone anion radical in the Q_A binding site of the photosynthetic bacterium Rhodobacter sphaeroides. The calculated spin density in the Q_A site model suggests that differential hydrogen bonding strength to the O1 and O4 oxygen atoms of the radical results in an asymmetric spin density distribution in the semiquinone anion free radical form. The origin of the spin density asymmetry is attributed to the presence of the divalent iron or zinc ion situated between the Q_A and Q_B sites.     

Investigation of the correlation between the energy of the highest occupied molecular orbital (HOMO) and the logarithm of the OH rate constant of hydrofluorocarbons and hydrofluoroethers

A regression based model was developed to determine whether highest occupied molecular orbital (HOMO) energies, calculated using Kohn-Sham orbital density functional theory (DFT), could be used to estimate the OH rate constants of hydrofluorocarbons (HFCs) and hydrofluoroethers (HFEs), proposed substitutes for stratospheric O₃ depleting chlorofluorocarbons. The goodness of fit of the DFT model was compared with a second regression model, derived using recently reported HOMO energies obtained from Hartree Fock theory (HFT). Both models were employed to predict OH rate constants for a number of HFCs and HFEs whose OH rate constants have not been measured, thus providing data on the types of chemical structures that may increase the OH reactivity of the substitute and hence decrease its contribution to global warming. The estimated percent standard errors in the OH rate constant HFT and DFT regression models were 72% and 78%, respectively. The goodness of fits were such that the models can differentiate between two rate constants only when their ratio exceeds about a factor of four. Model results suggest that (1) only a limited number of HFEs will have OH rate constants that are more than an order of magnitude greater than the value for their corresponding HFCs and (2) the strategy of introducing an ether linkage into an HFC to dramatically enhance its reactivity will be most effective for the least fluorinated HFCs. © 1994 John Wiley & Sons, Inc.

1 This article is a US Government work and, as such, is in the public domain in the United States of America.

     

DNA oligonucleotide-cis-platin binding: Ab initio interpretation of the vibrational spectra

The cis-platin binding to the d(CCTGGTCC)*d(GGACCAGG) model DNA octamer was monitored with infrared absorption (IR) and vibrational circular dichroism (VCD) spectroscopies. The spectra were modeled with the aid of density functional computations and a Cartesian coordinate-based transfer of molecular property tensors from smaller DNA fragments. Because of the fragmentation, the tensors could be calculated with a higher precision. Environmental effects, such as the presence of the solvent or the cis-platin ligand, could be included in the modeling. The solvent was modeled by an explicit inclusion of hydrogen-bound water molecules, positions of which were estimated from a molecular dynamics simulation, or by the polarized continuum COSMO model. The B3LYP and BPW91 functionals used for the calculations of the spectral parameters were combined with the relativistic LANL2DZ platinum pseudo-potentials. The simulations reproduced the main IR and VCD DNA spectral features and explained most of the changes observed experimentally upon metal binding. The results confirmed that the influence of the ligand on DNA vibrational properties is quite complex; it originates in the geometry deformation and normal mode coupling pattern changes of the platinated octamer, as well as in local perturbations of the electronic structure and force field of the GC base pairs to which the platinum is bound. Many of the local effects could be accounted for by a point charge used in place of the metal in the GC complex.     

Nuclear magnetic shielding constants of liquid water: insights from hybrid quantum mechanics/molecular mechanics models

We present a gauge-origin independent method for the calculation of nuclear magnetic shielding tensors of molecules in a structured and polarizable environment. The method is based on a combination of density functional theory (DFT) or Hartree-Fock wave functions with molecular mechanics. The method is unique in the sense that it includes three important properties that need to be fulfilled in accurate calculations of nuclear magnetic shielding constants: (i) the model includes electron correlation effects, (ii) the model uses gauge-including atomic orbitals to give gauge-origin independent results, and (iii) the effect of the environment is treated self-consistently using a discrete reaction-field methodology. The authors present sample calculations of the isotropic nuclear magnetic shielding constants of liquid water based on a large number of solute-solvent configurations derived from molecular dynamics simulations employing potentials which treat solvent polarization either explicitly or implicitly. For both the (17)O and (1)H isotropic shielding constants the best predicted results compare fairly well with the experimental data, i.e., they reproduce the experimental solvent shifts to within 4 ppm for the (17)O shielding and 1 ppm for the (1)H shielding.     

Low-field optically detected EPR spectroscopy of transient photoinduced radical pairs

The effects of simultaneously applied weak static and weak radio frequency magnetic fields on the recombination of transient (<100 ns) radical pairs have been investigated using a low-field optically detected electron paramagnetic resonance technique. Measurements on the photoinduced electron-transfer reaction of perdeuterated pyrene with 1,3-dicyanobenzene using a approximately 0.3 mT radio frequency field at three separate frequencies (5, 20, and 65 MHz) in the presence of 0-4 mT static fields yield spectra that are strikingly sensitive to the frequency of the time-dependent field, to the strength of the static field, and to the relative orientation of the two fields. The spectra are simulated using a modified form of the gamma-COMPUTE algorithm originally devised for calculating magic angle spinning NMR spectra of polycrystalline samples. The essential features of the spectra are consistent with the radical pair mechanism and were satisfactorily simulated using parameters whose values are either known independently or for which estimates are readily available. The calculations included hyperfine couplings to four deuterons in the pyrene cation radical and three protons in the 1,3-dicyanobenzene anion radical. Spin-selective recombination was modeled using an exponential distribution of radical encounter times. The results are discussed in the context of the proposal that radical pair chemistry forms the basis of the magnetoreceptor that allows birds to sense the Earth's magnetic field as a source of compass information during migration.