Theoretical study of the interaction pattern and the binding affinity between procaine and DNA bases

The interacting patterns and mechanism of the binding affinity between the local anaesthetic procaine and four DNA bases (adenine, cytosine, guanine and thymine) in neutral form have been investigated in gas phase using the Austin Model l and density functional methods. The results show that the complexes are mainly stabilized by the H-bonding interactions. The bond critical point properties of the optimized complexes were analyzed by using the atoms in molecules theory with DFT method and the results show that the presence of the C?H···O or C?H···N hydrogen bonding. The natural bond orbital analysis was performed to quantitatively evaluate the hydrogen bonding interaction. The interacting energy shows that the binding of procaine with guanine is the most strong, whereas its binding to cytosine exhibits relatively weaker stability. The strength order of the relevant transferred charge between procaine and DNA base with natural population analysis are consist with the HOMO–LUMO gap results for each complex. And the order is accord with the relevant electrochemical experimental results.     

Effect of hydrogen bonding on the photo-oxidation of DNA

The one-electron oxidation of DNA has been extensively studied as it leads to the formation of oxidative lesions that cause carcinogenesis and aging. In this paper, experimental results specifically addressing the effect of hydrogen bonding on the one-electron oxidation rate of nucleosides are presented. To separate the hydrogen bonding from the π-stacking effect, experiments were performed in dichloromethane, in which base-pair formation is possible at the monomer level. The effect of base pairing of guanine with cytosine on the rate constant of the electron transfer from guanine to electron acceptor molecules in the triplet excited state was investigated, and a selective enhancement of the electron transfer was observed for the guanine:cytosine base pair. By introducing a methyl or bromo group to the C5 position of cytosine, acceleration or suppression, respectively, of the one-electron oxidation relative to the guanine:cytosine base pair was observed. The results demonstrate that the one-electron oxidation rate of guanine in DNA can be regulated by introducing a substituent on the base-pairing cytosine.     

Insights into DNA binding of ruthenium arene complexes: role of hydrogen bonding and pi stacking

Density functional theory (DFT) methods are used to investigate the binding of ruthenium arene complexes, proposed as promising anticancer drugs, to isolated nucleobases. This shows a clear preference for binding at guanine over any other base and an approximately 100 kJ mol (-1) difference in binding between guanine and adenine in the gas phase, while binding to cytosine and inosine are intermediate in energy between these extremes. Solvation reduces binding energies and the discrimination between bases but maintains the overall pattern of binding. DFT and ab initio data on arene-base interactions in the absence of ruthenium show that stacking and hydrogen-bonding interactions play a significant role but cannot account for all of the energy difference between bases observed. Atoms-in-molecules analysis allows further decomposition of binding energies into contributions from covalent-binding, hydrogen-bonding, and pi-stacking interactions. Larger arenes undergo stabilizing stacking interactions, whereas N-H...X hydrogen bonding is independent of arene. Pairing of guanine to cytosine is affected by ruthenium complexation, with individual hydrogen-bonding energies being altered but the overall pairing energy remaining almost constant.     

Strikingly Different Effects of Hydrogen Bonding on the Photodynamics of Individual Nucleobases in DNA: Comparison of Guanine and Cytosine

Ab initio surface hopping dynamics calculations were performed to study the photophysical behavior of cytosine and guanine embedded in DNA using a quantum mechanical/molecular mechanics (QM/MM) approach. It was found that the decay rates of photo excited cytosine and guanine were affected in a completely different way by the hydrogen bonding to the DNA environment. In case of cytosine, the geometrical restrictions exerted by the hydrogen bonds did not influence the relaxation time of cytosine significantly due to the generally small cytosine ring puckering required to access the crossing region between excited and ground state. On the contrary, the presence of hydrogen bonds significantly altered the photodynamics of guanine. The analysis of the dynamics indicates that the major contribution to the lifetime changes comes from the interstrand hydrogen bonds. These bonds considerably restricted the out-of-plane motions of the NH(2) group of guanine which are necessary for the ultrafast decay to the ground state. As a result, only a negligible amount of trajectories decayed into the ground state for guanine embedded in DNA within the simulation time of 0.5 ps, while for comparison, the isolated guanine relaxed to the ground state with a lifetime of about 0.22 ps. These examples show that, in addition to phenomena related to electronic interactions between nucleobases, there also exist relatively simple mechanisms in DNA by which the lifetime of a nucleobase is significantly enhanced as compared to the gas phase.     

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.     

Charge transfer interaction and terahertz studies of a nonlinear optical material L-glutamine picrate: A DFT study

Charge transfer interaction, vibrational spectra, and DFT computation of l-glutamine picrate has been analyzed. The equilibrium geometry, bonding features, and harmonic vibrational wavenumbers have been investigated with the help of B3LYP density functional theory method. The natural bond orbital analysis confirms the occurrence of strong intramolecular hydrogen bonding in the molecule. Terahertz time-domain spectroscopy was used to detect the absorption spectra in the frequency range from 0.025 to 2.8 THz. The vibrational modes found in molecular crystalline materials should be described as phonon modes with strong coupling to the intramolecular vibrations.     

Homopairing possibilities of the DNA bases cytosine and guanine: an ab initio DFT study

All the planar homopairings of cytosine and guanine are reported for the first time in this study. The idea of binding sites suggested for the simple case of adenine homopairs (J. Phys. Chem. B 2005, 109, 11933) is shown to be applicable to more complicated molecules binding to each other via multiple hydrogen bonds and can be considered as a general method for constructing hydrogen bonding structures. As an example we consider homopairs formed by DNA bases cytosine and guanine, suggesting that there may be 13 cytosine and 17 guanine homopairs. However, only 11 cytosine and 15 guanine homopairs remain after atomic relaxation performed using ab initio density functional theory. Most of the homopairs obtained have not been studied before. The homopairs have significant binding energies, varying from -0.19 to -1.12 eV, that are explained by multiple hydrogen bonds formed between monomers in the pairs, up to four hydrogen bonds in most energetically favorable cases. The detailed information on all guanine and cytosine planar homopairs contained in this work can be used to construct various cytosine and guanine superstructures observed on different surfaces.     

Regulation of one-electron oxidation rate of guanine by base pairing with cytosine derivatives

Effects of base pairing on the one-electron oxidation rate of guanine derivatives, guanine, 8-bromoguanine, and 8-oxo-7,8-dihydroguanine have been studied. The one-electron oxidation rate of guanine derivatives was determined by triplet-quenching experiments, using N,N'-dibutylnaphthaldiimide (NDI) in the triplet excited state (3NDI*) and fullerene (C(60)) in the triplet excited state ((3)C(60*)) as oxidants. In all three guanine derivatives studied here, acceleration of the one-electron oxidation was observed upon hydrogen bonding with cytosine, which demonstrates lowering of the oxidation potential of guanine derivatives by base pairing with cytosine. When a methyl or bromo group was introduced to the C5 position of cytosine, acceleration or suppression of the one-electron oxidation relative to the guanine:cytosine base pair was observed, respectively. The results demonstrate that the one-electron oxidation rate of guanine in DNA can be regulated by introducing a substituent on base pairing cytosine.     

A Comparative Study of Hydrogen Bonding Using Density Functional Theory

The hydrogen‐bond energies of water dimer and water‐formaldehyde complexes have been studied using density functional theory (DFT). Basis sets up to aug‐cc‐pVXZ (X=D, T, Q) were used. It was found that counterpoise corrected binding energies using the aug‐cc‐pVDZ basis set are very close to those predicted with the aug‐cc‐pVQZ set. Comparative studies using various DFT functionals on these two systems show that results from B3LYP, mPW1PW91 and PW91PW91 functionals are in better agreements with those predicted using high‐level ab initio methods. These functionals were applied to the study of hydrogen bonding between guanine (G) and cytosine (C), and between adenine (A) and thy mine (T) base pairs. With the aug‐cc‐pVDZ basis set, the predicted binding energies of base pairs are in good agreement with the most elaborate ab initio results.     

Gas-phase DNA oligonucleotide structures. A QM/MM and atoms in molecules study

QM/MM calculations have been employed to investigate the role of hydrogen bonding and pi-stacking in single- and double-stranded DNA oligonucleotides. DFT calculations and Atoms in Molecules analysis on QM/MM-optimized structures allow characterization and estimation of the energies of pi-stacking and hydrogen-bond interactions. This shows that pi-stacking interactions depend on the number and the nature of the DNA bases for single-stranded nucleotides; for instance, guanines are found to be involved in strong hydrogen bonds, whereas adenines interact mainly via stacking interactions. The role of interbase hydrogen bonding was explored: the -NH2 groups of guanine, adenine, and cytosine participate in N-H...O and N-H...N interactions. These are much stronger in single-strand oligonucleotides, where the -NH2 groups are highly nonplanar. In double-stranded DNA, the strong base-pairing hydrogen bonds of complementary bases lead to more planar -NH2 groups, which tend to be involved in pi-stacking interactions rather than H-bonds. The use of AIM also allows us to evaluate the interplay of pi-stacking and H-bonding, suggesting that cooperativity does occur, but is generally limited to about 1-2 kcal/mol.     

How well can new-generation density functional methods describe stacking interactions in biological systems?

We compare the performance of four recently developed DFT methods (MPW1B95, MPWB1K, PW6B95, and PWB6K) and two previous, generally successful DFT methods (B3LYP and B97-1) for the calculation of stacking interactions in six nucleic acid bases complexes and five amino acid pairs and for the calculation of hydrogen bonding interactions in two Watson-Crick type base pairs. We found that the four newly developed DFT methods give reasonable results for the stacking interactions in the DNA base pairs and amino acid pairs, whereas the previous DFT methods fail to describe interactions in these stacked complexes. We conclude that the new generation of DFT methods have greatly improved performance for stacking interaction as compared to previously available methods. We recommend the PWB6K method for investigating large DNA or protein systems where stacking plays an important role.     

Nuclear quantum effect and temperature dependency on the hydrogen‐bonded structure of base pairs

The structure of Watson–Crick‐type adenine‐thymine and guanine‐cytosine pairs has been studied by hybrid Monte Carlo (HMC) and path integral hybrid Monte Carlo (PIHMC) simulations with the use of semiempirical PM6‐DH+ method in the gas phase. We elucidated the nuclear quantum effect and temperature dependency on the hydrogen‐bonded moiety of base pairs. It was shown that the contribution of nuclear quantum effect on the hydrogen‐bonded structure is significant not only at low temperature 150 K but also at temperature as high as 450 K. The relative position of hydrogen‐bonded proton between two heavy atoms and the nuclear quantum nature of the proton are also shown. Furthermore, we have applied principal component analysis to HMC and PIHMC simulations to analyze the nuclear quantum effect on intermolecular motions. We found that the ratio of Buckle mode (lowest vibrational mode from normal mode analysis) decreases due to the nuclear quantum effect, whereas that of Propeller mode (second lowest vibrational mode) increases. In addition, nonplanar structures of base pairs were found to become stable due to the nuclear quantum effect from two‐dimensional free energy landscape along Buckle and Propeller modes. © 2013 Wiley Periodicals, Inc.     

Cis Watson-Crick/Watson-Crick Base Pairs in DNA: A Density Functional Theory(DFT) Study

The four nucleic acid DNA bases(adenine, thymine, guanine, cytosine) and ten cis Watson-Crick/Watson-Crick(cis WC/WC) DNA base pairs were investigated by density functional theory(DFT) quantum chemical calculations. Geometry optimizations were carried out on the four bases and ten base pairs at the B3LYP level with 6-31G~(**) basis set. All the optimizations were performed within Cs symmetry. The optimum structures for the four bases and seven cis WC/WC base pairs were obtained, and Natural Bond Orbital analysis(NBO) was based on these structures. The possibilities of matches between any two of the four bases through their Watson-Crick(WC) edges were discussed. The structures of seven cis WC/WC base pairs change to a certain extent relative to these of the four bases due to the formation of hydrogen bonds. These base pairs existing in DNA have an important influence on the structural stability of the double helix. The analysis of the electronic structures and molecular orbitals for seven cis WC/WC base pairs can provide significant information about the relationship between charge transfer along the hydrogen bond and the Frontier orbitals of these base pairs.     

Halogen bonding in DNA base pairs

Halogen bonding (R-X···Y) is a qualitative analogue of hydrogen bonding that may prove useful in the rational design of artificial proteins and nucleotides. We explore halogen-bonded DNA base pairs containing modified guanine, cytosine, adenine and thymine nucleosides. The structures and stabilities of the halogenated systems are compared to the normal hydrogen bonded base pairs. In most cases, energetically stable, coplanar structures are identified. In the most favorable cases, halogenated base pair stabilities are within 2 kcal mol(-1) of the hydrogen bonded analogues. Among the halogens X = Cl, Br, and I, bromine is best suited for inclusion in these biological systems because it possesses the best combination of polarizability and steric suitability. We find that the most stable structures result from a single substitution of a hydrogen bond for a halogen bond in dA:dT and dG:dC base pairs, which allows 1 or 2 hydrogen bonds, respectively, to complement the halogen bond.     

MP2 study on the hydrogen bonding interaction between 5-hydroxy-5-methylhydantoin and DNA bases: A,C, G,T

The 5-hydroxy-5-methylhydantoin (5-OH-5-Me-dHyd) is a nucleobase lesion induced by the action of ionizing radiation on thymine residue in DNA. In this study, we present the hydrogen bonding base pairs involving 5-OH-5-Me-dHyd bound to the four bases in DNA: adenine (A), cytosine (C), guanine (G), and thymine (T). Full geometry optimizations have been performed for the studied complexes by MP2 method. The interaction energies were corrected for the basis-set superposition error (BSSE), using the full Boys–Bernardi counterpoise correction scheme. Hydrogen bonding patterns of these base pairs were characterized using NBO analysis and AIM analysis. According to the calculated binding energies and structural parameters, the stability of the base pairs decrease in the following order: 5-OH-5-Me-dHyd:G>5-OH-5-Me-dHyd:A>5-OH-5-Me -dHyd:C~5-OH-5-Me-dHyd:T.     

DFT and MP2 investigations on the hydrogen bonding interaction between 5,6-dihydrothymine and dna bases: A, C, G, T

5,6-Dihydrothymine (DHT) is a nucleobase lesion induced by the action of ionizing radiation on thymine residue in DNA. In this work, we present the hydrogen bonding base pairs involving 5,6-dihydrothymine bound to four bases in DNA: adenine (A), cytosine (C), guanine (G), and thymine (T). Full geometry optimizations are performed for the studied complexes by the B3LYP method. Interaction energies are corrected for the basis set superposition error, using the full Boys-Bernardi counterpoise correction scheme. Hydrogen bonding patterns of these base pairs are characterized using NBO and AIM analysis. According to the calculated binding energies and structural parameters, the stability of the base pairs decrease in the following order: DHT:G ～ DHT:A > DHT:C > DHT:T.     

Assignment of the lowest-lying THz absorption signatures in biotin and lactose monohydrate by solid-state density functional theory

The narrow terahertz (THz) features in crystalline biotin and lactose monohydrate observed in recent experimental studies are considered by solid-state density functional theory (DFT) calculations. The lowest-frequency THz features in both solid-state biotin and lactose monohydrate are assigned to external hindered rotational modes and not to the lowest-frequency internal modes predicted from isolated-molecule calculations. The motions of the molecules associated with these narrow THz features and the interactions between molecules in the hydrogen-bonded networks of these molecular crystals are discussed, and comparisons are made to similar studies on molecular crystals not exhibiting strong intermolecular interactions.     

Oligonucleotide cross-linking with copper ions. Spectral and quantum chemical study

Copper(II) complexes with synthetic oligonucleotides consisting of repeating adenine–thymine and guanine–cytosine complementary base pairs have been studied by UV spectroscopy and simulated by DFT quantum chemical calculations at the B3LYP/6-311G++(d,p) level of theory with inclusion of solvation (hydration) effects. The obtained data suggest selective interaction of copper(II) ions with guanine–cytosine complementary pairs, followed by DNA cross-linking at those sites.     

Theoretical analysis of the terahertz spectrum of the high explosive PETN.

The experimental solid-state terahertz (THz) spectrum (3 to 120 cm(-1)) of the high explosive pentaerythritol tetranitrate (PETN, C(5)H(6)N(4)O(12)) has been modeled using solid-state density functional theory (DFT) calculations. Solid-state DFT, employing the BP density functional, is in best qualitative agreement with the features in the previously reported THz spectrum. The crystal environment of PETN includes numerous intermolecular hydrogen-bonding interactions that contribute to large (up to 80 cm(-1)) calculated shifts in molecular normal-mode positions in the solid state. Comparison of the isolated-molecule and solid-state normal-mode calculations for a series of density functionals reveals the extent to which the inclusion of crystal-packing interactions and the relative motions between molecules are required for correctly reproducing the vibrational structure of solid-state THz spectra. The THz structure below 120 cm(-1) is a combination of both intermolecular (relative rotations and translations) and intramolecular (torsions, large amplitude motions) vibrational motions. Vibrational-mode analyses indicate that the first major feature (67.2 cm(-1)) in the PETN THz spectrum contains all of the optical rotational and translational cell modes and no internal (molecular) vibrational modes.     

MP2 study on the hydrogen-bonding interaction between 5-hydroxymethyl-uracil and DNA bases: A,C, G,T

The 5-hydroxymethyl-uracil (HmU) is a product of oxidative attack on the methyl group of thymine in DNA. In this work, we present the hydrogen bonding complexes formation involving HmU bound to the four bases in DNA: adenine (A), cytosine (C), guanine (G), and thymine (T). Full geometry optimizations have been performed for the studied complexes by MP2 method. The interaction energies were corrected for the basis-set superposition error (BSSE), using the full Boys-Bernardi counterpoise correction scheme. Hydrogen bonding patterns of these base pairs were characterized using NBO analysis and AIM analysis. According to the calculated binding energies and structural parameters, the stability of the base pairs decrease in the following order: HmU:A > HmU:G > HmU:C > HmU:T.