The Ad-MD method to calculate NMR shift including effects due to conformational dynamics: The 31 P NMR shift in DNA.

A method for averaging of NMR parameters by molecular dynamics (MD) has been derived from the method of statistical averaging in MD snapshots, benchmarked and applied to structurally dynamic interpretation of the 31 P NMR shift (δ31P ) in DNA phosphates. The method employs adiabatic dependence of an NMR parameter on selected geometric parameter(s) that is weighted by MD-calculated probability distribution(s) for the geometric parameter(s) (Ad-MD method). The usage of Ad-MD for polymers is computationally convenient when one pre-calculated structural dependence of an NMR parameter is employed for all chemically equivalent units differing only in dynamic behavior. The Ad-MD method is benchmarked against the statistical averaging method for δ31P in the model phosphates featuring distinctively different structures and dynamic behavior. The applicability of Ad-MD is illustrated by calculating 31 P NMR spectra in the Dickerson-Drew DNA dodecamer. δ31P was calculated with the B3LYP/IGLO-III/PCM(water) and the probability distributions for the torsion angles adjacent to the phosphorus atoms in the DNA phosphates were calculated using the OL15 force field.

Raman spectroscopy and DFT calculations of PEDOT:PSS in a dipolar field.

The conductive polymer-electrolyte interface plays an important role in many electrochemical devices. An unusual situation arises when a solvent-free ionic liquid (SF-IL) is used as the electrolyte because it behaves as a molten salt rather than an electrolyte solution. On the basis of Raman spectra, it was found that the presence of ion pairs of SF-IL in the vicinity of poly(3,4-ethylenedioxythiophene):poly(styrenesulfonate) (PEDOT:PSS) results in a decrease in the oxidation level of the polymer and an increase in the HOMO-LUMO gap. The process of polymer "dedoping" and the modification of the electronic structure of the polymer are illustrated by quantum chemical calculations.

QM and QM/MM umbrella sampling MD study of the formation of Hg(II)-thymine bond: Model for evaluation of the reaction energy profiles in solutions with constant pH.

The formation of the Hg-N3(T) bond between the 1-methylthymine (T) molecule and the hydrated Hg2+ cation was explored with the combined quantum mechanics/molecular mechanics (QM/MM) method including Free Energy Perturbation corrections. The thermodynamic properties were determined in the whole pH range, when these systems were explicitly investigated and considered as the QM part: (1) T + [Hg(H2 O)6 ]2+ , (2) T + [Hg(H2 O)5 (OH)]+ , (3) T + Hg(H2 O)4 (OH)2 , and (4) N3-deprotonated T + Hg(H2 O)4 (OH)2 . The MM part contained only solvent molecules and counterions. As a result, the dependence of Gibbs-Alberty reaction free energy on pH was obtained along the reaction coordinate. We found that an endoergic reaction in acidic condition up to pH < 4-5 becomes exoergic for a higher pH corresponding to neutral and basic solutions. The migration of the Hg2+ cation between N3 and O4/2 positions in dependence on pH is discussed as well. For the verification, DFT calculations of stationary points were performed confirming the qualitative trends of QM/MM MD simulations and NMR parameters were determined for them.

The mechanism of the glycosylase reaction with hOGG1 base-excision repair enzyme: concerted effect of Lys249 and Asp268 during excision of 8-oxoguanine.

The excision of 8-oxoguanine (oxoG) by the human 8-oxoguanine DNA glycosylase 1 (hOGG1) base-excision repair enzyme was studied by using the QM/MM (M06-2X/6-31G(d,p):OPLS2005) calculation method and nuclear magnetic resonance (NMR) spectroscopy. The calculated glycosylase reaction included excision of the oxoG base, formation of Lys249-ribose enzyme-substrate covalent adduct and formation of a Schiff base. The formation of a Schiff base with ΔG# = 17.7 kcal/mol was the rate-limiting step of the reaction. The excision of the oxoG base with ΔG# = 16.1 kcal/mol proceeded via substitution of the C1΄-N9 N-glycosidic bond with an H-N9 bond where the negative charge on the oxoG base and the positive charge on the ribose were compensated in a concerted manner by NH3+(Lys249) and CO2-(Asp268), respectively. The effect of Asp268 on the oxoG excision was demonstrated with 1H NMR for WT hOGG1 and the hOGG1(D268N) mutant: the excision of oxoG was notably suppressed when Asp268 was mutated to Asn. The loss of the base-excision function was rationalized with QM/MM calculations and Asp268 was confirmed as the electrostatic stabilizer of ribose oxocarbenium through the initial base-excision step of DNA repair. The NMR experiments and QM/MM calculations consistently illustrated the base-excision reaction operated by hOGG1.

Ascorbigen A-NMR identification.

The connectivities of all atoms in ascorbigen A, an important metabolite, were determined unambiguously for the first time. The connectivity between carbon atoms was established by 2D INADEQUATE, and one-bond 13 C-13 C coupling constants were determined for all pairs of directly connected carbon atoms except for two strongly coupled carbon pairs. The 13 C-13 C coupling in one of the pairs was proved by a modification of standard INADEQUATE; however, the signals from the other pair were too weak to be observed. The connectivity within the two strongly coupled C-C pairs was confirmed by a combination of COSY and gHSQC; the latter experiment also identified all C-H bonds. The proton nuclear magnetic resonance (1 H NMR) spectra in dry dimethyl sulfoxide allowed identification and assignment of the signals due to NH and OH protons. The derived structure, 3-((1H-indol-3-yl)methyl)-3,3a,6-trihydroxytetrahydrofuro[3,2-b]furan-2(5H)-one, agrees with the structure suggested for ascorbigen A in 1966. The density functional theory (DFT) calculations showed that among 16 possible stereoisomers, only two complied with the almost zero value of the measured 3 J(H6-H6a). Of the two stereoisomers, 3S,3aS,6S,6aR and 3R,3aR,6R,6aS, the latter was excluded on synthetic grounds. The nuclear Overhauser effect measurements unveiled close proximity between H2' proton of the indole and the H6a proton of the tetrahydrofuro[3,2-b]furan part. Detailed structural interpretation of the measured NMR parameters by means of DFT NMR was hampered by rotational flexibility of the indole and tetrahydrofuro[3,2-b]furan parts and inadequacy of Polarizable Continuum Model (PCM) solvent model.

The benchmark of 31P NMR parameters in phosphate: a case study on structurally constrained and flexible phosphate.

A benchmark for structural interpretation of the 31P NMR shift and the 2JP,C NMR spin-spin coupling in the phosphate group was obtained by means of theoretical calculations and NMR measurements in diethylphosphate (DEP) and 5,5-dimethyl-2-hydroxy-1,3,2-dioxaphosphinane 2-oxide (cDEP). The NMR parameters were calculated employing the B3LYP, BP86, BPW91, M06-2X, PBE0, KT2, KT3, MP2, and HF methods, and the 6-31+G(d), Iglo-n (n = II, III), cc-pVnZ (n = D, T, Q, 5), aug-cc-pVnZ (n = D, T and Q), and pcS-n and pcJ-n (n = 1, 2, 3, 4) bases, including the solvent effects described with explicit water molecules and/or the implicit Polarizable Continuum Model (PCM). The effect of molecular dynamics (MD) on NMR parameters was MD-calculated using the GAFF force field inclusive of explicit hydration with TIP3P water molecules. Both the optimal geometries and the dynamic behaviors of the DEP and cDEP phosphates differed notably, which allowed a reliable theoretical benchmark of the 31P NMR parameters for highly flexible and structurally constrained phosphate in a one-to-one relationship with the corresponding experiment. The calculated 31P NMR shifts were referenced employing three different NMR reference schemes to highlight the effect of the 31P NMR reference on the accuracy of the calculated 31P NMR shift. The relative Δδ(31P) NMR shift calculated employing the MD/B3LYP/Iglo-III/PCM method differed from the experiment by 0.16 ppm while the NMR shifts referenced to H3PO4 and/or PH3 deviated from the experiment notably more, which illustrated the superior applicability of the relative NMR reference scheme. The 2JP,C coupling in DEP and cDEP calculated employing the MD/B3LYP/Iglo-III(DSO,PSO,SD)/cc-PV5Z(FC)/PCM method inclusive of correction due to explicit hydration differed from the experiment by 0.32 Hz and 0.15 Hz, respectively. The NMR calculations demonstrated that reliable structural interpretation of the 31P NMR parameters in phosphate must involve both the structural and the dynamical components.

Impact of nucleic acid self-alignment in a strong magnetic field on the interpretation of indirect spin-spin interactions.

Heteronuclear and homonuclear direct (D) and indirect (J) spin-spin interactions are important sources of structural information about nucleic acids (NAs). The Hamiltonians for the D and J interactions have the same functional form; thus, the experimentally measured apparent spin-spin coupling constant corresponds to a sum of J and D. In biomolecular NMR studies, it is commonly presumed that the dipolar contributions to Js are effectively canceled due to random molecular tumbling. However, in strong magnetic fields, such as those employed for NMR analysis, the tumbling of NA fragments is anisotropic because the inherent magnetic susceptibility of NAs causes an interaction with the external magnetic field. This motional anisotropy is responsible for non-zero D contributions to Js. Here, we calculated the field-induced D contributions to 33 structurally relevant scalar coupling constants as a function of magnetic field strength, temperature and NA fragment size. We identified two classes of Js, namely (1)JCH and (3)JHH couplings, whose quantitative interpretation is notably biased by NA motional anisotropy. For these couplings, the magnetic field-induced dipolar contributions were found to exceed the typical experimental error in J-coupling determinations by a factor of two or more and to produce considerable over- or under-estimations of the J coupling-related torsion angles, especially at magnetic field strengths >12 T and for NA fragments longer than 12 bp. We show that if the non-zero D contributions to J are not properly accounted for, they might cause structural artifacts/bias in NA studies that use solution NMR spectroscopy.

Structure Determination of an Ag(I) -Mediated Cytosine-Cytosine Base Pair within DNA Duplex in Solution with (1) H/(15) N/(109) Ag NMR Spectroscopy.

The structure of an Ag(I) -mediated cytosine-cytosine base pair, C-Ag(I) -C, was determined with NMR spectroscopy in solution. The observation of 1-bond (15) N-(109) Ag J-coupling ((1) J((15) N,(109) Ag): 83 and 84 Hz) recorded within the C-Ag(I) -C base pair evidenced the N3-Ag(I) -N3 linkage in C-Ag(I) -C. The triplet resonances of the N4 atoms in C-Ag(I) -C demonstrated that each exocyclic N4 atom exists as an amino group (-NH2 ), and any isomerization and/or N4-Ag(I) bonding can be excluded. The 3D structure of Ag(I) -DNA complex determined with NOEs was classified as a B-form conformation with a notable propeller twist of C-Ag(I) -C (-18.3±3.0°). The (109) Ag NMR chemical shift of C-Ag(I) -C was recorded for cytidine/Ag(I) complex (δ((109) Ag): 442 ppm) to completed full NMR characterization of the metal linkage. The structural interpretation of NMR data with quantum mechanical calculations corroborated the structure of the C-Ag(I) -C base pair.

Limits in Proton Nuclear Singlet-State Lifetimes Measured with para-Hydrogen-Induced Polarization.

The synthesis of a hyperpolarized molecule was developed, where the polarization and the singlet state were preserved over two controlled chemical steps. Nuclear singlet-state lifetimes close to 6 min for protons are reported in dimethyl fumarate. Owing to the high symmetry (AA'X3 X3 ' and A2 systems), the singlet-state readout requires either a chemical desymmetrization or a long and repeated spin lock. Using DFT calculations and relaxation models, we further determine nuclear spin singlet lifetime limiting factors, which include the intramolecular dipolar coupling mechanism (proton-proton and proton-deuterium), the chemical shift anisotropy mechanism (symmetric and antisymmetric), and the intermolecular dipolar coupling mechanism (to oxygen and deuterium). If the limit of paramagnetic relaxation caused by residual oxygen could be lifted, the intramolecular dipolar coupling to deuterium would become the limiting relaxation mechanism and proton lifetimes upwards of 26 min could become available in the molecules considered here (dimethyl maleate and dimethyl fumarate).

The structure of metallo-DNA with consecutive thymine-HgII-thymine base pairs explains positive entropy for the metallo base pair formation.

We have determined the three-dimensional (3D) structure of DNA duplex that includes tandem Hg(II)-mediated T-T base pairs (thymine-Hg(II)-thymine, T-Hg(II)-T) with NMR spectroscopy in solution. This is the first 3D structure of metallo-DNA (covalently metallated DNA) composed exclusively of 'NATURAL' bases. The T-Hg(II)-T base pairs whose chemical structure was determined with the (15)N NMR spectroscopy were well accommodated in a B-form double helix, mimicking normal Watson-Crick base pairs. The Hg atoms aligned along DNA helical axis were shielded from the bulk water. The complete dehydration of Hg atoms inside DNA explained the positive reaction entropy (ΔS) for the T-Hg(II)-T base pair formation. The positive ΔS value arises owing to the Hg(II) dehydration, which was approved with the 3D structure. The 3D structure explained extraordinary affinity of thymine towards Hg(II) and revealed arrangement of T-Hg(II)-T base pairs in metallo-DNA.

Synthesis of locked cyclohexene and cyclohexane nucleic acids (LCeNA and LCNA) with modified adenosine units.

We describe here the preparation of conformationally locked cyclohexane nucleic acids designed as hybrids between locked nucleic acids (LNAs) and cyclohexene nucleic acids (CeNAs), both of which excel in hybridization with complementary RNAs. We have accomplished the synthesis of these adenine derivatives starting from a simple ketoester and installed all four chiral centres by means of total synthesis. The acquired monomers were incorporated into nonamer oligonucleotides.

Raman spectroscopic detection of the T-Hg II-T base pair and the ionic characteristics of mercury.

Developing applications for metal-mediated base pairs (metallo-base-pair) has recently become a high-priority area in nucleic acid research, and physicochemical analyses are important for designing and fine-tuning molecular devices using metallo-base-pairs. In this study, we characterized the Hg(II)-mediated T-T (T-Hg(II)-T) base pair by Raman spectroscopy, which revealed the unique physical and chemical properties of Hg(II). A characteristic Raman marker band at 1586 cm(-1) was observed and assigned to the C4=O4 stretching mode. We confirmed the assignment by the isotopic shift ((18)O-labeling at O4) and density functional theory (DFT) calculations. The unusually low wavenumber of the C4=O4 stretching suggested that the bond order of the C4=O4 bond reduced from its canonical value. This reduction of the bond order can be explained if the enolate-like structure (N3=C4-O4(-)) is involved as a resonance contributor in the thymine ring of the T-Hg(II)-T pair. This resonance includes the N-Hg(II)-bonded state (Hg(II)-N3-C4=O4) and the N-Hg(II)-dissociated state (Hg(II+) N3=C4-O4(-)), and the latter contributor reduced the bond order of N-Hg(II). Consequently, the Hg(II) nucleus in the T-Hg(II)-T pair exhibited a cationic character. Natural bond orbital (NBO) analysis supports the interpretations of the Raman experiments.

Influence of the O-phosphorylation of serine, threonine and tyrosine in proteins on the amidic ¹5N chemical shielding anisotropy tensors.

Density functional theory was employed to study the influence of O-phosphorylation of serine, threonine, and tyrosine on the amidic (15)N chemical shielding anisotropy (CSA) tensor in the context of the complex chemical environments of protein structures. Our results indicate that the amidic (15)N CSA tensor has sensitive responses to the introduction of the phosphate group and the phosphorylation-promoted rearrangement of solvent molecules and hydrogen bonding networks in the vicinity of the phosphorylated site. Yet, the calculated (15)N CSA tensors in phosphorylated model peptides were in range of values experimentally observed for non-phosphorylated proteins. The extent of the phosphorylation induced changes suggests that the amidic (15)N CSA tensor in phosphorylated proteins could be reasonably well approximated with averaged CSA tensor values experimentally determined for non-phosphorylated amino acids in practical NMR applications, where chemical surrounding of the phosphorylated site is not known a priori in majority of cases. Our calculations provide estimates of relative errors to be associated with the averaged CSA tensor values in interpretations of NMR data from phosphorylated proteins.

Formation of a thymine-Hg(II)-thymine metal-mediated DNA base pair: proposal and theoretical calculation of the reaction pathway.

A reaction mechanism that describes the substitution of two imino protons in a thymine:thymine (T:T) mismatched DNA base pair with a Hg(II) ion, which results in the formation of a (T)N3-Hg(II)-N3(T) metal-mediated base pair was proposed and calculated. The mechanism assumes two key steps: The formation of the first Hg(II)-N3(T) bond is triggered by deprotonation of the imino N3 atom in thymine with a hydroxo ligand on the Hg(II) ion. The formation of the second Hg(II)-N3(T) bond proceeds through water-assisted tautomerization of the remaining, metal-nonbonded thymine base or through thymine deprotonation with a hydroxo ligand of the Hg(II) ion already coordinated to the thymine base. The thermodynamic parameters ΔGR =-9.5 kcal mol(-1), ΔHR =-4.7 kcal mol(-1), and ΔSR =16.0 cal mol(-1) K(-1) calculated with the ONIOM (B3LYP:BP86) method for the reaction agreed well with the isothermal titration calorimetric (ITC) measurements by Torigoe et al. [H. Torigoe, A. Ono, T. Kozasa, Chem. Eur. J. 2010, 16, 13218-13225]. The peculiar positive reaction entropy measured previously was due to both dehydration of the metal and the change in chemical bonding. The mercury reactant in the theoretical model contained one hydroxo ligand in accord with the experimental pKa  value of 3.6 known for an aqua ligand of a Hg(II) center. The chemical modification of T:T mismatched to the T-Hg(II)-T metal-mediated base pair was modeled for the middle base pair within a trinucleotide B-DNA duplex, which ensured complete dehydration of the Hg(II) ion during the reaction.

Chemical shifts in nucleic acids studied by density functional theory calculations and comparison with experiment.

NMR chemical shifts are highly sensitive probes of local molecular conformation and environment and form an important source of structural information. In this study, the relationship between the NMR chemical shifts of nucleic acids and the glycosidic torsion angle, χ, has been investigated for the two commonly occurring sugar conformations. We have calculated by means of DFT the chemical shifts of all atoms in the eight DNA and RNA mono-nucleosides as a function of these two variables. From the DFT calculations, structures and potential energy surfaces were determined by using constrained geometry optimizations at the BP86/TZ2P level of theory. The NMR parameters were subsequently calculated by single-point calculations at the SAOP/TZ2P level of theory. Comparison of the (1)H and (13)C NMR shifts calculated for the mono-nucleosides with the shifts determined by NMR spectroscopy for nucleic acids demonstrates that the theoretical shifts are valuable for the characterization of nucleic acid conformation. For example, a clear distinction can be made between χ angles in the anti and syn domains. Furthermore, a quantitative determination of the χ angle in the syn domain is possible, in particular when (13)C and (1)H chemical shift data are combined. The approximate linear dependence of the C1' shift on the χ angle in the anti domain provides a good estimate of the angle in this region. It is also possible to derive the sugar conformation from the chemical shift information. The DFT calculations reported herein were performed on mono-nucleosides, but examples are also provided to estimate intramolecularly induced shifts as a result of hydrogen bonding, polarization effects, or ring-current effects.

Guanine bases in DNA G-quadruplex adopt nonplanar geometries owing to solvation and base pairing.

The effect of base pairing and solvation on pyramidalization of the glycosidic nitrogen found in the residues of parallel G-quadruplex with NDB ID UDF062 is analyzed and explained with theoretical calculations. The extent of the pyramidalization depends on the local geometry of the 2'-deoxyguanosine residues, namely on their glycosidic torsion and sugar pucker, which are predetermined by the 3D-architecture of G-quadruplex. Pyramidal inversion of the glycosidic nitrogen found in 2'-deoxyguanosines of G-quadruplex is induced owing to site-specifically coordinated solvent. Different adiabatic structural constraints used for fixing the base-to-sugar orientation of 2'-deoxyguanosine in geometry optimizations result in different extents of pyramidalization and induce pyramidal inversion of the glycosidic nitrogen. These model geometry constraints helped us analyze the effect of real constraints represented by explicit molecular environment of selected residues of the G-quadruplex. The maximal extent of the glycosidic nitrogen pyramidalization found in the high-resolution crystal structure corresponds to the calculation to deformation energy of only 1 kcal mol(-1). The out-of-plane deformations of nucleobases thus provide a way for compensating the site-specific external environmental stress on the G-quadruplex.

Detection of mercury-TpT dinucleotide binding by Raman spectra: a computational study.

The Hg(2+) ion stabilizes the thymine-thymine mismatched base pair and provides new ways of creating various DNA structures. Recently, such T-Hg-T binding was detected by the Raman spectroscopy. In this work, detailed differences in vibrational frequencies and Raman intensity patterns in the free TpT dinucleotide and its metal-mediated complex (TpT·Hg)(2) are interpreted on the basis of quantum chemical modeling. The computations verified specific marker Raman bands indicating the effect of mercury binding to DNA. Although the B3LYP functional well-describes the Raman frequencies, a dispersion correction had to be added for all atoms including mercury to obtain realistic geometry of the (TpT·Hg)(2) dimer. Only then, the DFT complex structure agreed with those obtained with the wave function-based MP2 method. The aqueous solvent modeled as a polarizable continuum had a minor effect on the dispersion interaction, but it stabilized conformations of the sugar and phosphate parts. A generalized definition of internal coordinate force field was introduced to monitor covalent bond mechanical strengthening and weakening upon the Hg(2+) binding. Induced vibrational frequency shifts were rationalized in terms of changes in electronic structure. The simulations thus also provided reliable insight into the complex structure and stability.

Probabilistic Interpretation of NMR J-Couplings Determines BI-BII State Equilibria in DNA.

Interpretation of 3JP,H3' NMR scalar spin-spin coupling constants in DNA becomes more reliable by including distinct structural states such as BI and BII, using the weighted-static or, better still, the recently implemented adiabatic-MD (Ad-MD) method. The calculation method employs an adiabatic ("Ad") dependence of 3JP,H3' coupling on NMR-assigned torsion angle, ε, weighted by P(ε) probability distribution calculated by molecular dynamics (MD). Ad-MD calculations enable cross-validation of the bsc1, OL15, and OL21 force fields and various parametrizations of the Karplus equation describing the dependence of 3JP,H3' coupling on ε torsion (KE). The mean absolute deviation of Ad-MD 3JP,H3' couplings from the experimental values in Dickerson-Drew DNA is comparable to the scatter of 3JP,H3' couplings among four separate NMR experiments. A commonly accepted assumption of homogeneity of one kind of structure-dynamic state within DNA (BI or BII) is questionable because the principal characteristics of relevant P(ε) probabilities (shapes and positioning) vary with DNA sequence. The theory outlined in the present work sets limits to future reparameterization of MD force fields, as relevant to NMR data.

On the role of mercury in the non-covalent stabilisation of consecutive U-Hg(II)-U metal-mediated nucleic acid base pairs: metallophilic attraction enters the world of nucleic acids.

Metal atoms with a closed-shell electronic structure and positive charge as for example the Au(I), Pt(II), Ag(I), Tl(I) or Hg(II) atoms do not in some compounds repel each other due to the so-called metallophilic attraction (P. Pyykkö, Chem. Rev., 1997, 97, 597-636). Here we highlight the role of the Hg(II)Hg(II) metallophilic attraction between the consecutive metal-mediated mismatched base pairs of nucleic acids. Usually, the base stacking dominates the non-covalent interactions between steps of native nucleic acids. In the presence of metal-mediated base pairs these non-covalent interactions are enriched by the metal-base interactions and the metallophilic attraction. The two interactions arising due to the metal linkage of the mismatches were found in this study to have a stabilizing effect on nucleic acid structure. The calculated data are consistent with recent experimental observations. The stabilization due to the metallophilic attraction seems to be a generally important concept for the nucleic acids containing heavy metals with short contacts.