Solvent effects on hydrogen bonds in Watson–Crick, mismatched, and modified DNA base pairs

Binding energies and hydrogen bonds effects on DNA-cisplatin interactions: a DFT-xTB study

CONTEXT

A systematic study of hydrogen bonds in base pairs and the interaction of cisplatin with DNA fragments was carried out. Structure, binding energies, and electron density were analyzed. xTB has proven to be an accurate method for obtaining structures and binding energies in DNA structures. Our xTB values for DNA base binding energy were in the same order and in some cases better than CAM-B3LYP values compared to experimental values. Double-stranded DNA-cisplatin structures have been calculated and the hydrogen bonds of water molecules are a decisive factor contributing to the preference for the cisplatin-Guanine interaction. Higher values of the water hydrogen bonding energies were obtained in cisplatin-Guanine structures. Furthermore, the electrostatic potential was used to investigate and improve the analysis of DNA-cisplatin structures.

METHODS

We applied the xTB method and the CAM-B3LYP functional combined with def2-SVP basis set to perform and analyze of the bonding energies of the cisplatin interaction and the effects of the hydrogen bonds. Results were calculated employing the xTB and the ORCA software.

Medium viscosity influence on the open states genesis in a DNA molecule

The studies were carried out by the mathematical modeling of DNA mechanical deformations. Numerical calculations done for the interferon alpha 17 gene, which consists of 980 base pairs. It has been established that the genesis and dynamics of open states in the DNA molecule depends on the magnitude of the external influence (torque) and on the viscosity of the environment. In addition, it is shown that the dynamics of open states zones can have a jump-like character with a small change in the magnitude of the torque. When torque is applied to all 980 base pairs of the gene, the following effect is observed: an increase in the viscosity of the medium leads to an increase in the value of the torque necessary for the occurrence of OS and DNA unwinding, i.e. viscosity plays an important stabilizing role in DNA dynamics. Under the influence of a localized torque on different (by the content of A-T and G-C pairs and location) regions of the interferon alpha 17 gene, it was found that the magnitude of the external torque necessary for the occurrence of open states at all calculated values of viscosity depends on the nucleotide composition. The dependence of the torque magnitude required for the open states occurrence on viscosity is observed when the torque is applied to areas close to the gene boundaries. At the same time, the significance of the end effect, which weakens DNA, decreased with increasing viscosity of the medium.Communicated by Ramaswamy H. Sarma.

Complementary base pair interactions between different rare tautomers of the second-generation artificial genetic alphabets

The functionality of a semisynthetic DNA in the biological environment will depend on the base pair nature of its complementary base pairs. To understand this, base pair interactions between complementary bases of recently proposed eight second-generation artificial nucleobases are studied herein by considering their rare tautomeric conformations and a dispersion-corrected density functional theoretic method. It is found that the binding energies of two hydrogen-bonded complementary base pairs are more negative than those of the three hydrogen-bonded base pairs. However, as the former base pairs are endothermic, the semisynthetic duplex DNA would involve the latter base pairs.

Solvent effects on the NMR shieldings of stacked DNA base pairs

Stacking effects are among the most important effects in DNA. We have recently studied their influence in fragments of DNA through the analysis of NMR magnetic shieldings, firstly in vacuo. As a continuation of this line of research we show here the influence of solvent effects on the shieldings through the application of both explicit and implicit models. We found that the explicit solvent model is more appropriate for consideration due to the results matching better in general with experiments, as well as providing clear knowledge of the electronic origin of the value of the shieldings. Our study is grounded on a recently developed theoretical model of our own, by which we are able to learn about the magnetic effects of given fragments of DNA molecules on selected base pairs. We use the shieldings of the atoms of a central base pair (guanine-cytosine) of a selected fragment of DNA molecules as descriptors of physical effects, like π-stacking and solvent effects. They can be taken separately and altogether. The effect of π-stacking is introduced through the addition of some pairs above and below of the central base pair, and now, the solvent effect is considered including a network of water molecules that consist of two solvation layers, which were fixed in the calculations performed in all fragments. We show that the solvent effects enhance the stacking effects on the magnetic shieldings of atoms that belong to the external N-H bonds. The net effect is of deshielding on both atoms. There is also a deshielding effect on the carbon atoms that belong to CO bonds, for which the oxygen atom has an explicit hydrogen bond (HB) with a solvent water molecule. Solvent effects are found to be no higher than a few percent of the total value of the shieldings (between 1% and 5%) for most atoms, although there are few for which such an effect can be higher. There is one nitrogen atom, the acceptor of the HB between guanine and cytosine, that is more highly shielded (around 15 ppm or 10%) when the explicit solvent is considered. In a similar manner, the most external nitrogen atom of cytosine and the hydrogen atom that is bonded to it are highly deshielded (around 10 ppm for nitrogen and around 3 ppm for hydrogen).

Halogen Bonds Exist between Noncovalent Ligands and Natural Nucleic Acids

Because of their strong electron-rich properties, nucleic acids (NAs) can theoretically serve as halogen bond (XB) acceptors. From a PDB database survey, Kolář found that no XBs are formed between noncovalent ligands and NAs. Through statistical database analysis, quantum-mechanics/molecular-mechanics (QM/MM) optimizations, and energy calculations, we find that XBs formed between natural NAs and noncovalent ligands are primarily underestimated and that NAs can serve as XB acceptors to interact with noncovalent halogen ligands. Finally, through energy calculations, natural bond orbital analysis, and noncovalent interaction analysis, XBs are confirmed in 13 systems, among which two systems (445D and 4Q9Q) have relatively strong XBs. In addition, on the basis of energy scanning of four model systems, we explore the geometric rule for XB formation in NAs. This work will inspire researchers to utilize XBs in rational drug design targeting NAs.

B-DNA Structure and Stability: The Role of Nucleotide Composition and Order

We have quantum chemically analyzed the influence of nucleotide composition and sequence (that is, order) on the stability of double-stranded B-DNA triplets in aqueous solution. To this end, we have investigated the structure and bonding of all 32 possible DNA duplexes with Watson-Crick base pairing, using dispersion-corrected DFT at the BLYP-D3(BJ)/TZ2P level and COSMO for simulating aqueous solvation. We find enhanced stabilities for duplexes possessing a higher GC base pair content. Our activation strain analyses unexpectedly identify the loss of stacking interactions within individual strands as a destabilizing factor in the duplex formation, in addition to the better-known effects of partial desolvation. Furthermore, we show that the sequence-dependent differences in the interaction energy for duplexes of the same overall base pair composition result from the so-called "diagonal interactions" or "cross terms". Whether cross terms are stabilizing or destabilizing depends on the nature of the electrostatic interaction between polar functional groups in the pertinent nucleobases.

Matrix effects on hydrogen bonding and proton transfer in fluoropyridine - HCl complexes

We report an extensive computational and spectroscopic study of several fluoropyridine-HCl complexes, and the parent, pyridine-HCl system. Matrix-IR spectra for pentafluoropyridine-HCl, 2,6-difluororpyridine-HCl, and 3,5-difluororpyridine-HCl in solid neon exhibit shifts for the H-Cl stretching band that parallel the effects of fluorination on hydrogen-bond strength. Analogous spectral shifts observed across various host environments (solid neon, argon, and nitrogen) for pentafluoropyridine-HCl and 2,6-difluororpyridine-HCl convey a systematically varying degree of matrix stabilization on the hydrogen bonds in these complexes. An extended quantum-chemical study of pyridine-HCl and eight fluorinated analogs, including 2-, 3-, and 4-fluoropyridine-HCl, 2,6- and 3,5-difluororpyridine-HCl, 2,4,6- and 3,4,5-trifluropyridine-HCl, as well as pentafluoropyridine-HCl, was also performed. Equilibrium structures and binding energies for the gas-phase complexes illustrate two clear trends in how fluorine substitution affects hydrogen bond strength; increasing fluorination weakens these interactions, yet substitution at the 2- and 6-positions has the most pronounced effect. Bonding analyses for a select subset of these systems reveal shifts in electron density that accompany hydrogen bonding, and most notably, the values of the electron density at the N-H bond critical points among the stronger systems in this subset significantly exceed those typical for moderately strong hydrogen-bonds. We also explored the effects of dielectric media on the structural and bonding properties of these systems. For pyridine-HCl, 3-fluoropyridine-HCl, and 3,5-difluororpyridine-HCl, a transition to proton transfer-type structures is observed at ε-values of 1.2, 1.5, and 2.0, respectively. This is signaled by key structural changes, as well as an increase in the negative charge on the chorine, and dramatic shifts in topological properties of the H-Cl and N-H bonds. In the case of pentafluoropyridine-HCl, and 2,6-difluororpyridine-HCl, we do not predict proton transfer in dielectric media up to ε = 20.0. However, there are clear indications that the media enhance hydrogen-bond strength, and moreover, these observations are completely consistent with the experimental IR spectra.

An IRMPD spectroscopic and computational study of protonated guanine-containing mismatched base pairs in the gas phase

Infrared multiple photon dissociation spectroscopy has been used to probe the structures of the three protonated base-pair mismatches containing 9-ethylguanine (9eG) in the gas phase. Some of these protonated base-pairs have been identified in RNA.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

The influence of substituents and the environment on the NMR shielding constants of supramolecular complexes based on A-T and A-U base pairs

In the present study, we have theoretically analyzed supramolecular complexes based on the Watson-Crick A-T and A-U base pairs using dispersion-corrected density functional theory (DFT). Hydrogen atoms H8 and/or H6 in the natural adenine and thymine/uracil bases were replaced, respectively, by substituents X8, Y6 = NH^-, NH₂, NH₃⁺ (N series), O^-, OH, OH₂⁺ (O series), F, Cl or Br (halogen series). We examined the effect of the substituents on the hydrogen-bond lengths, strength and bonding mechanism, and the NMR shielding constants of the C2-adenine and C2-thymine/uracil atoms in the base pairs. The general belief in the literature that there is a direct connection between changes in the hydrogen-bond strength and the C2-adenine shielding constant is conclusively rejected by our computations.

DNA Polymerase λ Active Site Favors a Mutagenic Mispair between the Enol Form of Deoxyguanosine Triphosphate Substrate and the Keto Form of Thymidine Template: A Free Energy Perturbation Study

Human DNA polymerase λ is an intermediate fidelity member of the X family, which plays a role in DNA repair. Recent X-ray diffraction structures of a ternary complex of a loop-deletion mutant of polymerase λ, a deoxyguanosine triphosphate analogue, and a gapped DNA show that guanine and thymine form a mutagenic mispair with an unexpected Watson-Crick-like geometry rather than a wobble geometry. Hence, there is an intriguing possibility that either thymine in the DNA or guanine in the deoxyguanosine triphosphate analogue may spend a substantial fraction of time in a deprotonated or enol form (both are minor species in aqueous solution) in the active site of the polymerase λ mutant. The experiments do not determine particular forms of the nucleobases that contribute to this mutagenic mispair. Thus, we investigate the thermodynamics of formation of various mispairs between guanine and thymine in the ternary complex at a neutral pH using classical molecular dynamics simulations and the free energy perturbation method. Our free energy calculations, as well as a comparison of the experimental and computed structures of mispairs, indicate that the Watson-Crick-like mispair between the enol tautomer of guanine and the keto tautomer of thymine is dominant. The wobble mispair between the keto forms of guanine and thymine and the Watson-Crick-like mispair between the keto tautomer of guanine and the enol tautomer of thymine are less prevalent, and mispairs that involve deprotonated guanine or thymine are thermodynamically unlikely. These findings are consistent with the experiment and relevant for understanding mechanisms of spontaneous mutagenesis.

B-DNA model systems in non-terran bio-solvents: implications for structure, stability and replication

We have computationally analyzed a comprehensive series of Watson–Crick and mismatched B-DNA base pairs, in the gas phase and in several solvents, including toluene, chloroform, ammonia, methanol and water, using dispersion-corrected density functional theory and implicit solvation.

Ab initio study of the excited-state proton transfer mechanisms for 3-hydroxy-2-(thiophen-2-yl)chromen-4-one

The optical absorption properties and fluorescence properties of 3-hydroxy-2-(thiophen-2-yl)chromen-4-one (3-HTC) were simulated by the time-dependent density functional theory (TDDFT) method.

2-Thiouracil deprived of thiocarbonyl function preferentially base pairs with guanine rather than adenine in RNA and DNA duplexes

2-Thiouracil-containing nucleosides are essential modified units of natural and synthetic nucleic acids. In particular, the 5-substituted-2-thiouridines (S2Us) present in tRNA play an important role in tuning the translation process through codon-anticodon interactions. The enhanced thermodynamic stability of S2U-containing RNA duplexes and the preferred S2U-A versus S2U-G base pairing are appreciated characteristics of S2U-modified molecular probes. Recently, we have demonstrated that 2-thiouridine (alone or within an RNA chain) is predominantly transformed under oxidative stress conditions to 4-pyrimidinone riboside (H2U) and not to uridine. Due to the important biological functions and various biotechnological applications for sulfur-containing nucleic acids, we compared the thermodynamic stabilities of duplexes containing desulfured products with those of 2-thiouracil-modified RNA and DNA duplexes. Differential scanning calorimetry experiments and theoretical calculations demonstrate that upon 2-thiouracil desulfuration to 4-pyrimidinone, the preferred base pairing of S2U with adenosine is lost, with preferred base pairing with guanosine observed instead. Therefore, biological processes and in vitro assays in which oxidative desulfuration of 2-thiouracil-containing components occurs may be altered. Moreover, we propose that the H2U-G base pair is a suitable model for investigation of the preferred recognition of 3'-G-ending versus A-ending codons by tRNA wobble nucleosides, which may adopt a 4-pyrimidinone-type structural motif.

Theoretical modelling of epigenetically modified DNA sequences

We report herein a set of calculations designed to examine the effects of epigenetic modifications on the structure of DNA. The incorporation of methyl, hydroxymethyl, formyl and carboxy substituents at the 5-position of cytosine is shown to hardly affect the geometry of CG base pairs, but to result in rather larger changes to hydrogen-bond and stacking binding energies, as predicted by dispersion-corrected density functional theory (DFT) methods. The same modifications within double-stranded GCG and ACA trimers exhibit rather larger structural effects, when including the sugar-phosphate backbone as well as sodium counterions and implicit aqueous solvation. In particular, changes are observed in the buckle and propeller angles within base pairs and the slide and roll values of base pair steps, but these leave the overall helical shape of DNA essentially intact. The structures so obtained are useful as a benchmark of faster methods, including molecular mechanics (MM) and hybrid quantum mechanics/molecular mechanics (QM/MM) methods. We show that previously developed MM parameters satisfactorily reproduce the trimer structures, as do QM/MM calculations which treat bases with dispersion-corrected DFT and the sugar-phosphate backbone with AMBER. The latter are improved by inclusion of all six bases in the QM region, since a truncated model including only the central CG base pair in the QM region is considerably further from the DFT structure. This QM/MM method is then applied to a set of double-stranded DNA heptamers derived from a recent X-ray crystallographic study, whose size puts a DFT study beyond our current computational resources. These data show that still larger structural changes are observed than in base pairs or trimers, leading us to conclude that it is important to model epigenetic modifications within realistic molecular contexts.

DNA Stability in Ionic Liquids and Deep Eutectic Solvents

DNA molecules are known as the genetic information carriers. Recently, they are being explored as a new generation of biocatalysts or chiral scaffolds for metal catalysts. There is also a growing interest of finding alternative solvents for DNA preservation and stabilization, including two unique types of solvents: ionic liquids (ILs) and deep eutectic solvents (DES). Therefore, it is important to understand how DNA molecules interact with these novel ionic solvent systems (i.e. ILs and DES). It is well known that inorganic di- and monovalent ions preferentially bind with major and minor grooves of DNA structures. However, in the case of ILs and DES, organic cation may intrude into the DNA minor grooves; more importantly, electrostatic attraction between organic cations and the DNA phosphate backbone becomes a predominant interaction, accompanying by hydrophobic and polar interactions between ILs and DNA major and minor grooves. In addition, anions may form hydrogen-bonds with cytosine, adenine and guanine bases. Despites these strong interactions, DNA molecules maintain double helical structure in most ionic solvent systems, especially in aqueous IL solutions. Furthermore, the exciting advances of G-quadruplexe DNA structures in ILs and DES are discussed.

Investigation of hydrogen bonds properties in the terephthalic acid crystal, using molecular dynamics method

The aim of this study was to perform calculations using the method of Car-Parrinello molecular dynamics, leading to the optimized geometry of the molecules of 1,4-benzenedicarboxylic acid (terephthalic acid) in crystals, for the hydrogen form and three variants of substitution of deuterium atoms inside a carboxyl group. Based on the results, trajectories and dipole moments were calculated, what makes possible to simulate vibrations in different systems, and to make calculation of theoretical infrared spectra and atomic power spectra. Theoretical results were compared with the experimental spectra, which verifies the correctness of the method and also was compared with the results obtained by quantum-mechanical calculations using DFT for the isolated dimer. Comparison of the spectra of different forms, allowed for in-depth analysis of the effect of isotopic substitution on the frequency of vibrations and shapes of bands, and confirm the presence of possible coupling effects and intra- and intermolecular interactions. Comparison with the DFT results for the dimer show influence of the crystal structure on the spectra.

B-DNA structure and stability: the role of hydrogen bonding, π–π stacking interactions, twist-angle, and solvation

Insight into structure and stability of B-DNA is obtained through systematic quantum chemical analyses of the roles played by hydrogen bonding, π–π stacking, solvation, and twist-angle.

Tautomerisation of thymine acts against the Hückel 4N + 2 rule. The effect of metal ions and H-bond complexations on the electronic structure of thymine

Not necessarily the π-electron delocalization is responsible for the stability of thymine tautomers.