The effect of methylation on the hydrogen-bonding and stacking interaction of nucleic acid bases

Methylated nucleosides play an important role in DNA/RNA function, and may affect republication by interrupting the base-pairing and base-stacking. In order to investigate the effect of methylation on the interaction between nucleic acid bases, this work presents the hydrogen-bonding and stacking interactions between 5-methylcytosine and guanine (G), cytosine (C) and G, 1-methyladenine and thymine (T), as well as adenine and T. Geometry optimization and potential energy surface scan have been performed for the involved complexes by MP2 calculations. The interaction energies, which were corrected for the basis-set superposition error by the full Boys–Bernardi counterpoise correction scheme, were used to evaluate the interaction intensity of these nucleic acid bases. The atoms in molecules theory and natural bond orbital analysis have been performed to study the hydrogen bonds in these complexes. The result shows that the methyl substitute contributes the stability to these complexes because it enhances either the hydrogen bonding or the staking interaction between nucleic acid bases studied.     

Density Functional Theory Study of the Interaction between Guanine and Catechin

The interacting patterns and mechanism of the catechin and guanine have been investigated with the density functional theory B3LYP method by 6‐31G* basis set. Fourteen stable structures for the catechin‐guanine complexes have been found which form two hydrogen bonds at least. The results indicate that the complexes are mainly stabilized by the hydrogen bonding interactions. At the same time, the number and strength of hydrogen bond play a co‐determinant parts in the stability of the complexes which can form two or more hydrogen bonds. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) have been adopted to investigate the hydrogen bonds involved in all systems. The interaction energies of all complexes have been corrected for basis set superposition error (BSSE), ranging from ?38.86 to ?14.56 kJ/mol. The results showed that the hydrogen bonding contributes to the interaction energies dominantly. The corresponding bonds stretching motions in all complexes are red‐shifted relative to that of the monomer, which is in agreement with experimental results.     

Role of hydrogen bonds in protein-DNA recognition: effect of nonplanar amino groups

Amino groups are one of the various types of hydrogen bond donors, abundantly found in protein main chains, protein side chains, and DNA bases. The polar hydrogen atoms of these groups exhibit short ranged, specific, and directional hydrogen bonds, which play a decisive role in the specificity and stability of protein-DNA complexes. To date, planar amino groups are only considered for the analysis of protein-DNA interfacial hydrogen bonds. This assumption regarding hydrogen atom positions possibly failed to establish the expected role of hydrogen bonds in protein-DNA recognition. We have performed ab initio quantum chemical studies on amino acid side chains and DNA bases containing amino groups as well as on specific hydrogen bonded residue pairs selected from high-resolution protein-DNA complex crystal structures. Our results suggest that occurrences of pyramidal amino groups are more probable in comparison with the usually adopted planar geometry. This increases the quality of the existing hydrogen bonds in almost all cases. Further, detailed analysis of protein-DNA interfacial hydrogen bonds in 107 crystal structures using the in-house program "pyrHBfind" indicates that consideration of energetically more preferred nonplanar amino groups improves the geometry of hydrogen bonds and also gives rise to new contacts amounting to nearly 14.5% of the existing interactions. Large improvements have been observed specifically for the amino groups of guanine, which faces the DNA minor groove and thus helps to resolve the problem of insufficient directional contacts observed in many minor groove binding complexes. Apart from guanine, improvement observed for asparagine, glutamine, adenine, or cytosine also indicates that the consideration of nonplanar amino groups leads to a more realistic scenario of hydrogen bonds occurring between protein and DNA residues.     

DNA-conducting polymer complexes: a computational study of the hydrogen bond between building blocks

Ab initio quantum mechanical calculations at the MP2 level were used for an extensive study concerning the stability of hydrogen-bonded complexes formed by pyrrole and thiophene, which are the most common building blocks of conducting polymers, and DNA bases. Results indicated that very stable complexes were formed with pyrrole, which shows a clear tendency to form specific hydrogen-bonding interactions with nucleic acid bases. Furthermore, the strength of such interactions depends significantly on the base, growing in the following order: thymine < adenine approximately equal to cytosine < guanine. On the contrary, thiophene formed complexes stabilized by nonspecific interactions between the pi-cloud of the ring and the N-H groups of the nucleic acid bases rather than specific hydrogen bonds. Overall, these results are fully consistent with experimental observations: polypyrrole is able not only to stabilize adducts with DNA but also to interact specifically, while the interactions of the latter with polythiophene and their derivatives are weaker and nonspecific.     

Hydrogen bonding and the conformations of poly(alkyl acrylamides)

The conformations of poly(alkyl acrylamide) oligomers in nonpolar solvents were studied using molecular dynamics techniques. Poly(methyl acrylamide) was found to collapse to a globule-like conformation at low temperatures; however, excluded volume effects inhibited the collapse of poly(octadecyl acrylamide). A high density of structured units, characterized by a trans-gauche-trans-trans-gauche-trans torsional sequence along the backbone, was noted in all simulations. Such units were found to create a particularly stable set of intramolecular hydrogen bonds. An oligomer constructed with these stable units was found to have significantly lower minimized energy than both the all-trans and the helical backbone conformations. The constructed conformation had lower Coulomb energy (more hydrogen bonds) than the all-trans conformation and lower dihedral energy (less backbone distortion) than the helical conformation. The propensity for poly(octadecyl acrylamide) to form hydrogen bonds introduced significant disorder into the orientation of the alkyl side chains. This disorder would inhibit crystallization and restrict the ability of such polymers to form epitaxial seeds for nucleating paraffin crystals.     

Hydrogen bonding and covalent effects in binding of cisplatin to purine bases: ab initio and atoms in molecules studies

Ab initio and density functional calculations are employed to investigate the role of hydrogen bonding in the binding of cisplatin to the purine bases guanine and adenine. Through the use of the theory of atoms in molecules (AIM), it is shown that hydrogen bonds are ubiquitous in such systems, with N-H...N and N-H...Cl interactions present in addition to the expected N-H...O. This in turn means that the known stability of cisplatin-guanine complexes cannot be ascribed solely to hydrogen bonding and allows decomposition of total binding energy into contributions from covalent and hydrogen bonds. To do so, a new method for predicting hydrogen bond energies from bond critical point properties is proposed, employing partial least-squares analysis to remove the family dependence of simple models. Still more hydrogen bond motifs are found in bifunctional complexes of the general type purine-[Pt(NH(3))(2)](2+)-purine, including purine...purine contacts, though again the energetics of these are insufficient to explain the observed trends in stability. Finally, the effect of platination on the pairing of guanine with cytosine is studied in a similar manner, revealing large redistributions of hydrogen bonding but surprisingly small overall changes in pairing energy.     

Calculation of the stabilization energies of oxidatively damaged guanine base pairs with guanine

DNA is constantly exposed to endogenous and exogenous oxidative stresses. Damaged DNA can cause mutations, which may increase the risk of developing cancer and other diseases. G:C-C:G transversions are caused by various oxidative stresses. 2,2,4-Triamino-5(2H)-oxazolone (Oz), guanidinohydantoin (Gh)/iminoallantoin (Ia) and spiro-imino-dihydantoin (Sp) are known products of oxidative guanine damage. These damaged bases can base pair with guanine and cause G:C-C:G transversions. In this study, the stabilization energies of these bases paired with guanine were calculated in vacuo and in water. The calculated stabilization energies of the Ia:G base pairs were similar to that of the native C:G base pair, and both bases pairs have three hydrogen bonds. By contrast, the calculated stabilization energies of Gh:G, which form two hydrogen bonds, were lower than the Ia:G base pairs, suggesting that the stabilization energy depends on the number of hydrogen bonds. In addition, the Sp:G base pairs were less stable than the Ia:G base pairs. Furthermore, calculations showed that the Oz:G base pairs were less stable than the Ia:G, Gh:G and Sp:G base pairs, even though experimental results showed that incorporation of guanine opposite Oz is more efficient than that opposite Gh/Ia and Sp.     

Theoretical Investigation on the Geometries and Properties of Guanine-BX3 （X = F, Cl） Complex

Geometries and binding energies were predicted at the B3LYP/6-311+G* level for the guanine-BX3 (X = F, Cl) systems and four isomers with no imaginary frequencies have been obtained for both guanine-BF3 and guanine-BCl3, respectively. Single energy calculations using much larger basis sets (6-311+G(2df,p) and aug-cc-pVDZ were carried out as well. It was found that the most stable isomer of guanine-BF3 is BF3 connected to N3 of guanine with the stabilization energy of –19.93 kcal/mol (BSSE corrected), while that of guanine-BCl3 is BCl3 connected to O10 of guanine having stabilization energy of –15.02 kcal/mol at the same level. The analyses for the combining interaction between BX3 and guanine with the atom-in-molecules theory (AIM) and natural bond orbital (NBO) methods have been performed. The results indicated that all the isomers are formed with σ-p type interactions between guanine and BX3, in which pyridine-type nitrogen or carbonyl oxygen or nitrogen atom of amino group offers its lone pair electrons to the empty p orbital of boron atom and the concomitance of charge transfer from guanine to BX3 has occurred. Still, one or two hydrogen bonds exist in some isomers of guanine-BX3 system and contribute to the stability of complex systems. Frequency analysis suggested that the stretching vibration of BX3 undergoes a red shift in complexes. Guanine-BF3 complex is more stable than guanine-BCl3 although the B–Y (Y=N, O) bond distance in the latter is shorter.     

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.     

Analysis of the structures,energetics, and vibrational frequencies for the hydrogen-bonded interaction of nucleic acid bases with Carmustine pharmaceutical agent: a detailed computational approach

In the present study, it is attempted to scrutinize the hydrogen bonding interaction between Carmustine drug and DNA pyrimidine bases by means of density functional theory calculations regarding their geometries, binding energies, vibrational frequencies, and topological features of the electron density in the gas phase and the water solution. Based on the density functional theory results, it is found that the process of intermolecular interaction between Carmustine drug and nucleobases is exothermic and all of the optimized configurations are stable. Furthermore, the negative stability energy represented by a polarizable continuum model shows the significant increase in the solubility of the nucleobase after hydrogen bonding intermolecular interaction in the presence of water solvent. It is also found that the intermolecular hydrogen bonds between drug and the nucleobases play the significant role in the stability of the physisorption configurations. Hydrogen bond energies for hydrogen-bonded complexes are obtained from Espinosa method and the atoms-in-molecules theory are also applied to get a more precise insight into the nature of the intermolecular hydrogen bond interactions.     

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.     

Hydrogen‐bonding Structures and Energetics of Acrylamide Isomers,Tautomers, and Dimers: An ab initio Study and Spectral Analysis

Hydrogen‐bonding patterns and energetic profiles of acrylamide isomers (syn‐ and skew‐), tautomers (amide and imidic acid forms) and 13 stable dimers have been studied using the second‐order Møller–Plesset perturbation theory with basis sets up to aug‐cc‐pVTZ. Syn‐acrylamide is the most stable monomer with a reaction barrier of 4.15 kcal/mol for the syn–skew isomerization reaction. The direct amide–imidic acid tautomerization reaction is separated by too high a barrier to surpass. The most stable dimer corresponds to the planar double‐hydrogen‐bonded configuration, indicating its crucial role in determining the stability of the formed complex. Moreover, hydrogen bonds have significant effects on the infrared spectral features, which can be consistently explained solely based on the acrylamide dimeric structures and energetics without monomeric and dimeric tautomer forms. The results are useful for studying the stability of the acrylamide clusters in condensed‐phase samples such as those in food chemistry studies.     

Minor Groove Binding between Norfloxacin and DNA Duplexes in Solution： A Molecular Dynamics Study

Quinolones,a series of nalidixic acid analogues,represent a class of extremely potent antibacterial drugs.Although the functional target of these agents such as norfloxacin has been identified as the enzyme DNA gyrase,the direct binding site of the drug is the DNA itself1.Interactions between quinolones and diversified DNA have been studied by several biochemical,spectroscopic and computational methods2-8.Shen et al.proposed a cooperative drug-DNA binding model.The essential feature of th…     

带电组氨酸侧链与DNA碱基间非键作用强度的理论研究

采用MP2方法和6-31+G(d,p)基组优化得到了带有一个正电荷的组氨酸侧链与4个DNA碱基间形成的18个氢键复合物的气相稳定结构, 从文献中获取了组氨酸侧链与DNA碱基间形成的12个堆积和T型复合物的气相稳定结构, 使用包含基组重叠误差(BSSE)校正的MP2方法和aug-cc-pVTZ基组及密度泛函理论M06-2X-D3方法和aug-cc-pVDZ基组计算了这些复合物的结合能. 研究结果表明, 包含BSSE校正的M06-2X-D3方法和aug-cc-pVDZ基组能够给出较准确的结合能; 气相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用明显强于堆积作用和T型作用, 组氨酸侧链最易通过离子氢键与胞嘧啶C和鸟嘌呤G作用形成氢键复合物, 组氨酸与胞嘧啶C和鸟嘌呤G间的T型作用强于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用; 水相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用仍明显强于堆积作用和T型作用, 组氨酸侧链更易与胞嘧啶C和鸟嘌呤G相互作用形成氢键复合物, 但是最强的组氨酸侧链与胞嘧啶C间的T型作用明显弱于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用, 说明水相条件下组氨酸侧链与DNA碱基间主要通过离子氢键作用形成氢键复合物.     

鸟嘌呤与氨基酸残基体系的极化力场

发展了应用于鸟嘌呤G和氨基酸残基体系的浮动电荷力场, 该力场明确定义了孤对电子和键的电荷和位置, 通过电荷随着环境的浮动来体现极化效应; 通过氢键拟合函数k_HB描绘了氢键键能. 应用量子化学方法, 对G与氨基酸残基体系从氢键、 几何结构及电荷分布3个方面展开计算及分析, 并以其为基准, 确定参数发展了适用于G与氨基酸残基氢键体系的ABEEMσπ PFF. 采用3种不同力场模拟目标分子的结构和性质. 模拟结果表明, 发展的ABEEMσπ PFF与量子化学方法具有最好的一致性, 可用于模拟生物大分子体系.     

Branch migration of Holliday junction in RuvA tetramer complex studied by umbrella sampling simulation using a path‐search algorithm

Branch migration of the Holliday junction takes place at the center of the RuvA tetramer. To elucidate how branch migration occurs, umbrella sampling simulations were performed for complexes of the RuvA tetramer and Holliday junction DNA. Although conventional umbrella sampling simulations set sampling points a priori, the umbrella sampling simulation in this study set the sampling points one by one in order to search for a realistic path of the branch migration during the simulations. Starting from the X‐ray structure of the complex, in which the hydrogen bonds between two base‐pairs were unformed, the hydrogen bonds between the next base‐pairs of the shrinking stems were observed to start to disconnect. At the intermediate stage, three or four of the eight unpaired bases interacted closely with the acidic pins from RuvA. During the final stage, these bases moved away from the pins and formed the hydrogen bonds of the new base‐pairs of the growing stems. The free‐energy profile along this reaction path showed that the intermediate stage was a meta‐stable state between two free‐energy barriers of about 10 to 15 kcal/mol. These results imply that the pins play an important role in stabilizing the interactions between the pins and the unpaired base‐pairs. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Structure and hydrogen bonding in polyhydrated complexes of guanine

Abstract Molecular structure of complexes of guanine with 12, 13, 16, and 17 water molecules were calculated using B3LYP/6-311G(d,p) level of theory. Interaction with water results in some deformation of geometrical parameters of guanine, which can be described as contribution of zwitter-ionic resonant form into the structure of DNA base. Saturation of water binding sites within guanine creates possibilities for the formation of the N···H–O hydrogen bond where the nitrogen atom of amino group acts as proton acceptor. The NBO analysis of guanine–water interactions reveals that hydrogen bonds involving the N(3) and N(7) atoms of guanine represent a case of mixed N···H–O/π···H–O hydrogen bonds where contribution of π-system into total energy of interaction varies from 3% to 41%. This contribution significantly depends on orientation of the hydrogen atom of water molecule with respect to plane of purine bicycle and influence of neighboring water molecules. Graphical Abstract      

Hydrogen-bonding interaction of urea with DNA bases: A density functional theory study

This work deals with the interaction between urea and DNA bases (adenine, thymine, guanine, and cytosine). The optimized geometries, binding energies, and harmonic vibrational frequencies are calculated using the DFT/B3LYP functional combined with the 6–31+G(d,p) basis set. Their interactions are studied aiming to understand more about the nature of the intercalation binding forces between urea and DNA. Fourteen stable complexes are found on the potential energy surface. The structures are cyclic; they are stabilized by NH...O/N and CH...O interactions. The binding energies range from −19.9 kJ·mol⁻¹ to −74.0 kJ·mol⁻¹. The obtained formation energies indicate that Urea:G and Urea:C are more favorable than Urea:T and Urea:A. In addition, the Atoms in Molecules theory is performed to study the hydrogen bonds in the complexes.     

Association Complexes of Calix[6]arenes with Amino Acids Explained by Energy-Partitioning Methods

Intermolecular complexes with calixarenes are intriguing because of multiple possibilities of noncovalent binding for both polar and nonpolar molecules, including docking in the calixarene cavity. In this contribution calix[6]arenes interacting with amino acids are studied with an additional aim to show that tools such as symmetry-adapted perturbation theory (SAPT), functional-group SAPT (F-SAPT), and systematic molecular fragmentation (SMF) methods may provide explanations for different numbers of noncovalent bonds and of their varying strength for various calixarene conformers and guest molecules. The partitioning of the interaction energy provides an easy way to identify hydrogen bonds, including those with unconventional hydrogen acceptors, as well as other noncovalent bonds, and to find repulsive destabilizing interactions between functional groups. Various other features can be explained by energy partitioning, such as the red shift of an IR stretching frequency for some hydroxy groups, which arises from their attraction to the phenyl ring of calixarene. Pairs of hydrogen bonds and other noncovalent bonds of similar magnitude found by F-SAPT explain an increase in the stability of both inclusion and outer complexes.     

The hydrogen bonding properties of cytosine: a computational study of cytosine complexed with hydrogen fluoride, water, and ammonia

Density functional theory is used to study the hydrogen bonding pattern in cytosine, which does not contain alternating proton donor and acceptor sites and therefore is unique compared with the other pyrimidines. Complexes between various small molecules (HF, H(2)O, and NH(3)) and four main binding sites in (neutral and (N1) anionic) cytosine are considered. Two complexes (O2(N1) and N3(N4)) involve neighboring cytosine proton acceptor and donor sites, which leads to cooperative interactions and bidendate hydrogen bonds. The third (less stable) complex (N4) involves a single cytosine donor. The final (O2-N3) complex involves two cytosine proton acceptors, which leads to an anticooperative hydrogen bonding pattern for H(2)O and NH(3). On the neutral surface, the anticooperative O2-N3 complex is less stable than those involving bidentate hydrogen bonds, and the H(2)O complex cannot be characterized when diffuse functions are included in the (6-31G(d,p)) basis set. On the contrary, the anionic O2-N3 structure is the most stable complex, while the HF and H(2)O N3(N4) complexes cannot be characterized with diffuse functions. B3LYP and MP2 potential energy surface scans are used to consider the relationship between the water N3(N4) and O2-N3 complexes. These calculations reveal that diffuse functions reduce the conversion barrier between the two complexes on both the neutral and anionic surfaces, where the reduction leads to a (O2-N3) energy plateau on the neutral surface and complete (N3(N4)) complex destabilization on the anionic surface. From these complexes, the effects of hydrogen bonds on the (N1) acidity of cytosine are determined, and it is found that the trends in the effects of hydrogen bonds on the (N1) acidity are similar for all pyrimidines.