Vladimír Sychrovský

Nuclear magnetic resonance spin–spin coupling constants from coupled perturbed density functional theory

For the first time, a complete implementation of coupled perturbed density functional theory (CPDFT) for the calculation of NMR spin–spin coupling constants (SSCCs) with pure and hybrid DFT is presented. By applying this method to several hydrides, hydrocarbons, and molecules with multiple bonds, the performance of DFT for the calculation of SSCCs is analyzed in dependence of the XC functional used. The importance of electron correlation effects is demonstrated and it is shown that the hybrid functional B3LYP leads to the best accuracy of calculated SSCCs. Also, CPDFT is compared with sum-over-states (SOS) DFT where it turns out that the former method is superior to the latter because it explicitly considers the dependence of the Kohn–Sham operator on the perturbed orbitals in DFT when calculating SSCCs. The four different coupling mechanisms contributing to the SSCC are discussed in connection with the electronic structure of the molecule.

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 HgII-mediated T–T base pairs (thymine–HgII–thymine, T–HgII–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–HgII–T base pairs whose chemical structure was determined with the 15N 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–HgII–T base pair formation. The positive ΔS value arises owing to the HgII dehydration, which was approved with the 3D structure. The 3D structure explained extraordinary affinity of thymine towards HgII and revealed arrangement of T–HgII–T base pairs in metallo-DNA.

Raman spectroscopic detection of the T-HgII-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 HgII-mediated T-T (T-HgII-T) base pair by Raman spectroscopy, which revealed the unique physical and chemical properties of HgII. 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 (18O-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-HgII-T pair. This resonance includes the N-HgII-bonded state (HgII-N3-C4=O4) and the N-HgII-dissociated state (HgII+ N3=C4-O4−), and the latter contributor reduced the bond order of N-HgII. Consequently, the HgII nucleus in the T-HgII-T pair exhibited a cationic character. Natural bond orbital (NBO) analysis supports the interpretations of the Raman experiments.

Charge Transport in DNA Oligonucleotides with Various Base-Pairing Patterns

We combined various experimental (scanning tunneling microscopy and Raman spectroscopy) and theoretical (density functional theory and molecular dynamics) approaches to study the relationships between the base-pairing patterns and the charge transfer properties in DNA 32-mer duplexes that may be relevant for identification and repair of defects in base pairing of the genetic DNA and for DNA use in nanotechnologies. Studied were two fully Watson-Crick (W-C)-paired duplexes, one mismatched (containing three non-W-C pairs), and three with base pairs chemically removed. The results show that the charge transport varies strongly between these duplexes. The conductivity of the mismatched duplex is considerably lower than that of the W-C-paired one despite the fact that their structural integrities and thermal stabilities are comparable. Structurally and thermally much less stable abasic duplexes have still lower conductivity but not markedly different from the mismatched duplex. All duplexes are likely to conduct by the hole mechanism, and water orbitals increase the charge transport probability.

Formation of a Thymine‐HgII‐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.

Interactions of hydrated divalent metal cations with nucleic acid bases. How to relate the gas phase data to solution situation and binding selectivity in nucleic acids

Binding modes of hydrated Zn2+ and Mg2+ cations to the N7 and O6 positions of guanine have been characterized by state-of-the-art ab initio model calculations. We show how differences in the electronic structure of the cationic complexes, as revealed by gas-phase computations, translate into differences in the biological function of the two studied metals. The thermodynamic driving force of the metal binding process is estimated on the basis of interaction energies and total electronic energies. The computed results unambiguously reveal that the N7 position of guanine exhibits a greater propensity to bind Zn2+ than Mg2+ while both cations have a similar affinity to bind to O6. Contrary to the intuitive expectations, however, the computed data do not suggest any superiority of the N7 inner shell binding mode for Zn2+ compared to the O6 binding. For Mg2+ the O6 inner shell binding mode is favored over the N7 one. The gas-phase data, when properly exrapolated, provide a relevant picture of many (though not all) fundamental aspects of the diversity of cation binding to nucleic acids.

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.

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.

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.

Effect of local sugar and base geometry on 13C and 15N magnetic shielding anisotropy in DNA nucleosides

Density functional theory was employed to study the dependence of 13C and 15N magnetic shielding tensors on the glycosidic torsion angle (χ) and conformation of the sugar ring in 2′-deoxyadenosine, 2′-deoxyguanosine, 2′-deoxycytidine, and 2′-deoxythymidine. In general, the magnetic shielding of the glycosidic nitrogens and the sugar carbons was found to depend on both the conformation of the sugar ring and χ. Our calculations indicate that the magnetic shielding anisotropy of the C6 atom in pyrimidine and the C8 atom in purine bases depends strongly on χ. The remaining base carbons were found to be insensitive to both sugar pucker and χ re-orientation. These results call into question the underlying assumptions of currently established methods for interpreting residual chemical shift anisotropies and 13C and 15N auto- and cross-correlated relaxation rates and highlight possible limitations of DNA applications of these methods.