精氨酸侧链和核酸碱基间离子氢键作用强度分析

采用MP2/6-31+G(d,p)方法优化得到了22个由精氨酸侧链与碱基尿嘧啶、 胸腺嘧啶、 胞嘧啶、 鸟嘌呤及腺嘌呤形成的氢键复合物的气相稳定结构, 使用包含BSSE校正的MP2/aug-cc-pVTZ方法计算得到了复合物的气相结合能, 通过MP2/6-31+G(d,p)方法和PCM模型优化得到了复合物的水相稳定结构, 采用MP2/aug-cc-pVTZ方法和PCM模型计算得到了复合物的水相结合能. 研究发现, 精氨酸侧链与碱基间的离子氢键作用强度与单体间电荷转移量、 氢键临界点电子密度及二阶作用稳定化能密切相关. 与中性氢键相比, 离子氢键作用具有更显著的共价作用成分. 研究还发现, 精氨酸侧链和碱基间形成的氢键复合物的稳定性次序可以通过氢键受体碱基分子上氧原子和氮原子的质子化反应焓变进行预测, 质子化反应焓变越负, 形成的氢键复合物越稳定.     

Modelling protein-RNA interactions: an electron density study of the guanidinium and formate complexes with RNA bases

The complexes formed by the double interaction established between RNA bases and guanidinium and formate ions, as a model for the interacting groups of arginine and glutamic or aspartic amino acid side chains, have been theoretically studied. A density functional theory method (B3LYP/6-31 + G**) has been used for this study. The range of interaction energies obtained allowed for a distinction between bidentate and bifurcate hydrogen bond interactions. The analysis of the electron density and the natural bond orbital analysis shows that these complexes are bound by double hydrogen bonds established between the donor and acceptor groups of guanidinium and formate respectively and those of the RNA bases. Comparisons are made with the results obtained in some previous theoretical and experimental studies.     

Energetic,Topological and Electric Field Analyses of Cation-Cation Nucleic Acid Interactions in Watson-Crick Disposition

A theoretical study of the effect of the diprotonation on the nucleic acid bases (A : U, A : T and G : C) in Watson-Crick conformation has been carried out by means of DFT computational methods in vacuum. In addition, the corresponding neutral and monoprotonated binary complexes have been considered. Most of the diprotonated species studied are stable, even though the binding energy is positive due to the overall repulsive electrostatic term. Local electrostatic attractive forces in the regions of hydrogen bonds (HBs) are responsible for equilibrium geometries, as shown by the electric field lines connecting the electrophilic and nucleophilic sites involved in the HB interactions. Secondary electrostatic effects also affect the assembling of the nucleic acid complexes in either neutral or cationic form. In particular, the electric field lines flowing from electrophilic sites in one base to nucleophilic sites in the other reinforce the linking between them. Hence, when the nucleophilic site concerns the free lone pair of the heteroatom involved in the HB interaction as acceptor, the HB distance shortens. However, if the free lone pair of the HB acceptor interacts with an electrophilic site in the same molecule, the HB distance elongates, weakening the HB interaction. The topological analysis of the electron density distribution in HB regions indicates that neutral, monoprotonated and diprotonated complexes show no differences in the nature of their HB's.     

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.     

Computational study of hydrogen-bonded complexes of HOCO with acids: HOCO⋯HCOOH, HOCO⋯H(2)SO(4), and HOCO⋯H(2)CO(3)

Quantum chemistry calculations at the density functional theory (DFT) (B3LYP), MP2, QCISD, QCISD(T), and CCSD(T) levels in conjunction with 6-311++G(2d,2p) and 6-311++G(2df,2p) basis sets have been performed to explore the binding energies of open-shell hydrogen bonded complexes formed between the HOCO radical (both cis-HOCO and trans-HOCO) and trans-HCOOH (formic acid), H(2)SO(4) (sulfuric acid), and cis-cis-H(2)CO(3) (carbonic acid). Calculations at the CCSD(T)∕6-311++G(2df,2p) level predict that these open-shell complexes have relatively large binding energies ranging between 9.4 to 13.5 kcal∕mol and that cis-HOCO (cH) binds more strongly compared to trans-HOCO in these complexes. The zero-point-energy-corrected binding strengths of the cH?Acid complexes are comparable to that of the formic acid homodimer complex (～13-14 kcal∕mol). Infrared fundamental frequencies and intensities of the complexes are computed within the harmonic approximation. Infrared spectroscopy is suggested as a potential useful tool for detection of these HOCO?Acid complexes in the laboratory as well as in various planetary atmospheres since complex formation is found to induce large frequency shifts and intensity enhancement of the H-bonded OH stretching fundamental relative to that of the corresponding parent monomers. Finally, the ability of an acid molecule such as formic acid to catalyze the inter-conversion between the cis- and trans-HOCO isomers in the gas phase is also discussed.     

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.     

水溶液中氨基酸侧链对G∶C碱基对间氢键影响的理论研究

利用量子化学方法研究了气相和水溶液下,氨基酸侧链与鸟嘌呤和胞嘧啶间的氢键作用.应用B3LYP/6-31+G(d,p)方法优化复合物几何结构,使用MP2/aug-cc-p VDZ方法进行复合物能量、自然键轨道(NBO)电荷和二阶稳定化能的计算.结果表明,水溶液可使氨基酸侧链与碱基或碱基对之间氢键键能显著减小;带电复合物气相和水溶液氢键键能之差范围为50.63~146.48 k J/mol,中性为0.17~24.94 k J/mol;电荷的转移量与氢键键能成正比,电荷转移量越多,复合物越稳定;二阶稳定化能与氢键键长成反比,与电荷转移量成正比,且气相与水溶液氢键二阶稳定化能之比约为两相的电荷转移量之比.水溶液对该类体系中氢键作用具有明显影响.     

复合物Ca+-RNA碱基的结构与稳定性

Gas-phase metal ion affinities and optimized structures of RNA nucleic acid bases for the Ca+ were determined at a density functional level employing the hybrid B3LYP exchange correlation potential in connection with the 6-311+G(2df,2p) basis set. All the molecular complexes, obtained by the interaction between several low-lying tautomers of RNA nucleic acid and Ca+ on the different binding sites, were considered. For Cytosine, the most stable complex was obtained starting from the most stable tautomer of the free nucleic acid base tautomers. As to thymine, the bond energy of the ion with the most stable tautomer of the free nucleic acid base is the weakest among the three tautomer’s complexes, and that of the ion with least stable tautomer of the free nucleic acid base is the strongest . Uracil is similar to thymine. The two kinds of relation, bond energy and total energy for the complex, are in disagreement, as the metal affinities of RNA bases for the Ca+ depend on binding sites, and total energy of complex (Ca+-RNA base) relies on all atoms and their relative positions in the complex.     

碱基与氮甲基乙酰胺相互作用的从头算和ABEEMσπ研究

选取氮甲基乙酰胺(N-Methylacetamide,NMA)和DNA碱基为研究对象,应用ABEEMσπ/MM方法研究了NMA与碱基之间形成的二聚体的几何构型、电荷分布及结合能等性质,并进行了相应的从头算(abinitio)MP2水平研究,在B3LYP/6-311 ++G(d,p)水平下对结构进行了优化.将2种方法的研...     

A polarizable dipole–dipole interaction model for evaluation of the interaction energies for NH···OC and CH···OC hydrogen‐bonded complexes

In this article, a polarizable dipole–dipole interaction model is established to estimate the equilibrium hydrogen bond distances and the interaction energies for hydrogen‐bonded complexes containing peptide amides and nucleic acid bases. We regard the chemical bonds N? H, C?O, and C? H as bond dipoles. The magnitude of the bond dipole moment varies according to its environment. We apply this polarizable dipole–dipole interaction model to a series of hydrogen‐bonded complexes containing the N? H···O?C and C? H···O?C hydrogen bonds, such as simple amide‐amide dimers, base‐base dimers, peptide‐base dimers, and β‐sheet models. We find that a simple two‐term function, only containing the permanent dipole–dipole interactions and the van der Waals interactions, can produce the equilibrium hydrogen bond distances compared favorably with those produced by the MP2/6‐31G(d) method, whereas the high‐quality counterpoise‐corrected (CP‐corrected) MP2/aug‐cc‐pVTZ interaction energies for the hydrogen‐bonded complexes can be well‐reproduced by a four‐term function which involves the permanent dipole–dipole interactions, the van der Waals interactions, the polarization contributions, and a corrected term. Based on the calculation results obtained from this polarizable dipole–dipole interaction model, the natures of the hydrogen bonding interactions in these hydrogen‐bonded complexes are further discussed. © 2013 Wiley Periodicals, Inc.     

槲皮素与尿嘧啶和胸腺嘧啶氢键的结合位点

深入理解药物分子和核酸碱基间的相互作用机制对合理设计研发新型高效药物有重要意义. 本文运用密度泛函理论B3LYP方法对核酸碱基尿嘧啶和胸腺嘧啶与药物分子槲皮素间的氢键相互作用位点进行了研究. 使用B3LYP/6-31G(d)方法优化得到了30个稳定的氢键复合物结构, 使用B3LYP/6-311++G(3df,2p)方法计算了这些复合物的结合能. 研究结果表明, 槲皮素可以使用5个不同的结合位点与尿嘧啶或胸腺嘧啶形成氢键复合物, 尿嘧啶或胸腺嘧啶可以使用3个不同的结合位点与槲皮素形成氢键复合物. 当槲皮素的结合位点固定时, 槲皮素与尿嘧啶的位点u1或胸腺嘧啶的位点t1形成的氢键作用最强, 与位点u2或位点t2形成的氢键强度最弱; 当尿嘧啶或胸腺嘧啶的作用位点固定时, 二者与槲皮素的位点qu1 形成的氢键作用最强, 与位点qu5 作用强度次之, 与位点qu3的作用强度最弱. 分子中原子(AIM)和自然键轨道(NBO)分析计算结果表明, 轨道作用在氢键中起重要作用.     

Bond energies and attachments sites of sodium and potassium cations to DNA and RNA nucleic acid bases in the gas phase

Gas-phase metal affinities of DNA and RNA bases for the Na(+) and K(+) ions were determined at density functional level employing the hybrid B3LYP exchange correlation potential in connection with the 6-311+G(2df,2p) basis set. All the molecular complexes, obtained by the interaction between several low-lying tautomers of nucleic acid bases and the alkali ions on the different binding sites, were considered. Structural features of the sodium and potassium complexes were found to be similar except in some uracil and thymine compounds in which the tendency of potassium ion toward monocoordination appeared evident. B3LYP bond energies for both metal ions were in agreement with the available experimental results in the cases of uracil and thymine for which the most stable complex was obtained starting from the most stable tautomer of the free nucleic acid base. For adenine, although the interaction of the ions with the most stable free tautomer generated the least stable molecular complex, the best agreement with experiment was found in just this case. For the remaining cytosine and guanine bases, our calculations indicated that the metal ion affinity value closest to experiment should be determined taking into account the role played by the different tautomers of the free bases with similar energy and all the possible complexes obtained by them.     

Structural and thermodynamic features of complexes formed by DNA and synthetic polynucleotides with dodecylamine and dodecyltrimethylammonium bromide

Complex formation of native and denatured DNA, single-stranded polyribonucleotides poly(A) and poly(U), as well as double-stranded poly(A).poly(U) with dodecylamine (DDA) and dodecyltrimethylammonium bromide (DTAB) has been studied by UV-, CD-, IR-spectroscopy and fluorescence analysis of hydrophobic probe pyrene. DDA and DTAB were shown to bind cooperatively with DNA and polyribonucleotides, resulting in the formation of complexes containing hydrophobic micelle-like clusters. Critical aggregation concentration (CAC) of DDA and DTAB shifts sharply to lower values (30-50 times) in the presence of DNA and polynucleotides as compared to critical micelle concentration (CMC) of free DDA and DTAB in solution. The analysis of binding isotherms within the frame of the model of cooperative binding of low-molecular ligands to linear polymers allowed us to determine the thermodynamic parameters of complex formation and estimate the contribution of electrostatic interaction of positively charged heads of amphiphiles with negatively charged phosphate groups of DNA and polyribonucleotides, and hydrophobic interaction of aliphatic chains to complex stability. Electrostatic interaction was shown to make the main contribution to the stability of DNA complexes with DDA, while preferential contribution of hydrophobic interactions is characteristic of DTAB complexes with DNA. The opposite effect of DDA and DTAB on the thermal stability of DNA double helix was demonstrated from UV-melting of DNA-while DTAB stabilizes the DNA helix, DDA, to the contrary, destabilizes it. The destabilizing effect of DDA seems to originate from the displacement of intramolecular hydrogen bonds in complementary Watson-Crick A.T and G.C base pairs with intermolecular H-bonds between unsubstituted DDA amino groups and proton-accepting sites of nucleic bases.     

Interactions of Safranin T with Nucleic Acids

The interactions of Safranin T (ST) with several nucleic acids have been investigated by electrochemical, UV‐visible and CD spectroscopic techniques. The form of the nucleic acid‐ST complexes is sensitive to the ratio of the two species. Two electrochemically inactive complexes such as, nucleic acid‐ST and nucleic acid‐2ST, were formed while ST interacts with nucleic acids. Two processes were obtained from spectral experiments: (1) at the high value of R (R is defined as the ratio of the total concentration of ST to that of nucleic acid), ST is groove‐binding with stacking, (2) at the low value of R, ST is groove‐binding without stacking. Intrinsic binding constants were obtained by spectral methods. The experiments also show that electrostatic binding plays an important role in the interaction of ST with nucleic acids.     

Gas-Phase Hydration Thermochemistry of Sodiated and Potassiated Nucleic Acid Bases

Hydration reactions of sodiated and potassiated nucleic acid bases (uracil, thymine, cytosine, and adenine) produced by electrospray have been studied in a gas phase using the pulsed ion-beam high-pressure mass spectrometer. The thermochemical properties, ΔH ( o ) ( n ), ΔS ( o ) ( n ), and ΔG ( o ) ( n ), for the hydrated systems were obtained from hydration equilibrium measurement. The structural aspects of the hydrated complexes are discussed in conjunction with available literature data. The correlation between water binding energies in the hydrated complexes and the corresponding metal ion affinities of nucleobases suggests that a significant (if not dominant) amount of the canonical structure of cytosine undergoes tautomerization during electrospray ionization, and the thermochemical values for cationized cytosine probably correspond to a mixture of tautomeric complexes.     

Physical nature of ethidium and proflavine interactions with nucleic acid bases in the intercalation plane

On the basis of the crystallographic structures of three nucleic acid intercalation complexes involving ethidium and proflavine, we have analyzed the interaction energies between intercalator chromophores and their four nearest bases, using a hybrid variation-perturbation method at the second-order M?ller-Plesset theory level (MP2) with a 6-31G(d,p) basis set. A total MP2 interaction energy minimum precisely reproduces the crystallographic position of the ethidium chromophore in the intercalation plane between UA/AU bases. The electrostatic component constitutes the same fraction of the total energy for all three studied structures. The multipole electrostatic interaction energy, calculated from cumulative atomic multipole moments (CAMMs), was found to converge only after including components above the fifth order. CAMM interaction surfaces, calculated on grids in the intercalation planes of these structures, reasonably reproduce the alignment of intercalators in crystal structures; they exhibit additional minima in the direction of the DNA grooves, however, which also need to be examined at higher theory levels if no crystallographic data are given.     

Measurement and theory of hydrogen bonding contribution to isosteric DNA base pairs

We address the recent debate surrounding the ability of 2,4-difluorotoluene (F), a low-polarity mimic of thymine (T), to form a hydrogen-bonded complex with adenine in DNA. The hydrogen bonding ability of F has been characterized as small to zero in various experimental studies, and moderate to small in computational studies. However, recent X-ray crystallographic studies of difluorotoluene in DNA/RNA have indicated, based on interatomic distances, possible hydrogen bonding interactions between F and natural bases in nucleic acid duplexes and in a DNA polymerase active site. Since F is widely used to measure electrostatic contributions to pairing and replication, it is important to quantify the impact of this isostere on DNA stability. Here, we studied the pairing stability and selectivity of this compound and a closely related variant, dichlorotoluene deoxyriboside (L), in DNA, using both experimental and computational approaches. We measured the thermodynamics of duplex formation in three sequence contexts and with all possible pairing partners by thermal melting studies using the van't Hoff approach, and for selected cases by isothermal titration calorimetry (ITC). Experimental results showed that internal F-A pairing in DNA is destabilizing by 3.8 kcal/mol (van't Hoff, 37 °C) as compared with T-A pairing. At the end of a duplex, base-base interactions are considerably smaller; however, the net F-A interaction remains repulsive while T-A pairing is attractive. As for selectivity, F is found to be slightly selective for adenine over C, G, T by 0.5 kcal mol, as compared with thymine's selectivity of 2.4 kcal/mol. Interestingly, dichlorotoluene in DNA is slightly less destabilizing and slightly more selective than F, despite the lack of strongly electronegative fluorine atoms. Experimental data were complemented by computational results, evaluated at the M06-2X/6-31+G(d) and MP2/cc-pVTZ levels of theory. These computations suggest that the pairing energy of F to A is ~28% of that of T-A, and most of this interaction does not arise from the F···HN interaction, but rather from the CH···N interaction. The nucleobase analogue shows no inherent selectivity for adenine over other bases, and L-A pairing energies are slightly weaker than for F-A. Overall, the results are consistent with a small favorable noncovalent interaction of F with A offset by a large desolvation cost for the polar partner. We discuss the findings in light of recent structural studies and of DNA replication experiments involving these analogues.     

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.     

四甲基吡啶卟啉与核酸相互作用的稳态和时间分辨光谱

采用稳态吸收和荧光光谱、圆二色谱和皮秒时间分辨荧光光谱手段, 研究了5,10,15,20-四[4-(N-甲基吡啶)]卟啉(TMPyP4)与腺嘌呤(A)、鸟嘌呤(G)、胸腺嘧啶(T)和胞嘧啶(C)等4种碱基, 以及相应的核苷、核苷酸和单链DNA的结合能力和光谱学性质. 研究结果发现, 嘌呤与TMPyP4的结合能力比嘧啶的强. 对于某一碱基系列, 结合能力强弱顺序依次为: 碱基~核苷<核苷酸<单链DNA. 时间分辨荧光谱研究发现, 除鸟嘌呤外, 核酸和TMPyP4复合物的荧光动力学均含有快(1~2 ns)和慢(约10 ns)两个衰减过程, 它们分别是由激基复合体和环境极性对激发态TMPyP4分子的影响所致. 单链DNA能诱导TMPyP4产生诱导圆二色信号, 而单分子(碱基、核苷、核苷酸)则无此功能.