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.

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.

HERMES - A Software Tool for the Prediction and Analysis of Magnetic-Field-Induced Residual Dipolar Couplings in Nucleic Acids.

Field-Induced Residual Dipolar Couplings (fiRDC) are a valuable source of long-range information on structure of nucleic acids (NA) in solution. A web application (HERMES) was developed for structure-based prediction and analysis of the (fiRDCs) in NA. fiRDC prediction is based on input 3D model structure(s) of NA and a built-in library of nucleobase-specific magnetic susceptibility tensors and reference geometries. HERMES allows three basic applications: (i) the prediction of fiRDCs for a given structural model of NAs, (ii) the validation of experimental or modeled NA structures using experimentally derived fiRDCs, and (iii) assessment of the oligomeric state of the NA fragment and/or the identification of a molecular NA model that is consistent with experimentally derived fiRDC data. Additionally, the program's built-in routine for rigid body modeling allows the evaluation of relative orientation of domains within NA that is in agreement with experimental fiRDCs.

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.

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 structural model of Zika virus RNA-dependent RNA polymerase in complex with RNA for rational design of novel nucleotide inhibitors.

Zika virus is a global health threat due to significantly elevated risk of fetus malformations in infected pregnant women. Currently, neither an effective therapy nor a prophylactic vaccination is available for clinical use, desperately necessitating novel therapeutics and approaches to obtain them. Here, we present a structural model of the Zika virus RNA-dependent RNA polymerase (ZIKV RdRp) in complex with template and nascent RNAs, Mg2+ ions and accessing nucleoside triphosphate. The model allowed for docking studies aimed at effective pre-screening of potential inhibitors of ZIKV RdRp. Applicability of the structural model for docking studies was illustrated with the NITD008 artificial nucleotide that is known to effectively inhibit the function of the ZIKV RdRp. The ZIKV RdRp - RNA structural model is provided for all possible variations of the nascent RNA bases pairs to enhance its general utility in docking and modelling experiments. The developed model makes the rational design of novel nucleosides and nucleotide analogues feasible and thus provides a solid platform for the development of advanced antiviral therapy.

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.

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.

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).

Benchmark Theoretical and Experimental Study on (15)N NMR Shifts of Oxidatively Damaged Guanine.

The (15)N NMR shifts of 9-ethyl-8-oxoguanine (OG) were calculated and measured in liquid DMSO and in crystal. The OG molecule is a model for oxidatively damaged 2'-deoxyguanosine that occurs owing to oxidative stress in cell. The DNA lesion is repaired with human 8-oxoguanine glycosylase 1 (hOGG1) base-excision repair enzyme, however, the exact mechanism of excision of damaged nucleobase with hOGG1 is currently unknown. This benchmark study on (15)N NMR shifts of OG aims their accurate structural interpretation and calibration of the calculation protocol utilizable in future studies on mechanism of hOGG1 enzyme. The effects of NMR reference, DFT functional, basis set, solvent, structure, and dynamics on calculated (15)N NMR shifts were first evaluated for OG in crystal to calibrate the best performing calculation method. The effect of large-amplitude motions on (15)N NMR shifts of OG in liquid was calculated employing molecular dynamics. The B3LYP method with Iglo-III basis used for B3LYP optimized geometry with 6-311++G(d,p) basis and including effects of solvent and molecular dynamic was the calculation protocol used for calculation of (15)N NMR shifts of OG. The NMR shift of N9 nitrogen of OG was particularly studied because the atom is involved in an N-glycosidic bond that is cleaved with hOGG1. The change of N9 NMR shift owing to oxidation of 9-ethylguanine (G) measured in liquid was -27.1 ppm. The calculated N9 NMR shift of OG deviated from experiment in crystal and in liquid by 0.45 and 0.65 ppm, respectively. The calculated change of N9 NMR shift owing to notable N9-pyramidalization of OG in one previously found polymorph was 20.53 ppm. We therefore assume that the pyramidal geometry of N9 nitrogen that could occur for damaged DNA within hOGG1 catalytic site might be detectable with (15)N NMR spectroscopy. The calculation protocol can be used for accurate structural interpretation of (15)N NMR shifts of oxidatively damaged guanine DNA residue.

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.

Structures, physicochemical properties, and applications of T-Hg(II)-T, C-Ag(I)-C, and other metallo-base-pairs.

Recently, metal-mediated base-pairs (metallo-base-pairs) have been studied extensively with the aim of exploring novel base-pairs; their structures, physicochemical properties, and applications have been studied. This trend has become more evident after the discovery of Hg(II)-mediated thymine-thymine (T-Hg(II)-T) and Ag(I)-mediated cytosine-cytosine (C-Ag(I)-C) base-pairs. In this article, we focus on the basic science and applications of these metallo-base-pairs, which are composed of natural bases.

Direct detection of the mercury-nitrogen bond in the thymine-Hg(II)-thymine base-pair with (199)Hg NMR spectroscopy.

We have observed the 1-bond (199)Hg-(15)N J-coupling ((1)J((199)Hg,(15)N) = 1050 Hz) within the Hg(II)-mediated thymine-thymine base pair (T-Hg(II)-T). This strikingly large (1)J((199)Hg,(15)N) is the first one for canonical sp(2)-nitrogen atoms, which can be a sensitive structure-probe of N-mercurated compounds and a direct evidence for N-mercuration.

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.

Theoretical and experimental study of charge transfer through DNA: impact of mercury mediated T-Hg-T base pair.

DNA-Hg complexes may play an important role in sensing DNA defects or in detecting the presence of Hg in the environment. A fundamental way of characterizing DNA-Hg complexes is to study the way the electric charge is transferred through the molecular chain. The main goal of this contribution was to investigate the impact of a mercury metal cation that links two thymine bases in a DNA T-T mismatched base pair (T-Hg-T) on charge transfer through the DNA molecule. We compared the charge transfer efficiencies in standard DNA, DNA with mismatched T-T base pairs, and DNA with a T-Hg(II)-T base pair. For this purpose, we measured the temperature dependence of steady-state fluorescence and UV-vis of the DNA molecules. The experimental results were confronted with the results obtained employing theoretical DFT methods. Generally, the efficiency of charge transfer was driven by mercury changing the spatial overlap of bases.

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.

Charge transfer through DNA/DNA duplexes and DNA/RNA hybrids: complex theoretical and experimental studies.

Oligonucleotides conduct electric charge via various mechanisms and their characterization and understanding is a very important and complicated task. In this work, experimental (temperature dependent steady state fluorescence spectroscopy, time-resolved fluorescence spectroscopy) and theoretical (Density Functional Theory) approaches were combined to study charge transfer processes in short DNA/DNA and RNA/DNA duplexes with virtually equivalent sequences. The experimental results were consistent with the theoretical model - the delocalized nature of HOMO orbitals and holes, base stacking, electronic coupling and conformational flexibility formed the conditions for more effective short distance charge transfer processes in RNA/DNA hybrids. RNA/DNA and DNA/DNA charge transfer properties were strongly connected with temperature affected structural changes of molecular systems - charge transfer could be used as a probe of even tiny changes of molecular structures and settings.