An explanation of the induction of mutations by 2-aminopurine from an ab initio molecular orbital study

The molecular mechanism of induction of mutations by 2-aminopurine (AP) was studied by an ab initio molecular orbital method. Cytosine (C) is converted to its disfavored imino tautomer more easily than AP, judging from the calculated total energies of the bases and the base analogue. This suggests that a stable AP:C base mispair via two hydrogen bonds can be formed with the imino tautomer of C. These results stress the importance of the imino form of C in AP-induced mutagenesis and support the 'trigger mechanism', in which formation of one hydrogen bond between AP and C is considered to stimulate the tautomeric shift of AP or C. The calculated relative stabilities of various base pairs and mispairs were in good agreement with experimental findings.     

A double-deletion method to quantifying incremental binding energies in proteins from experiment: example of a destabilizing hydrogen bonding pair

The contribution of a specific hydrogen bond in apoflavodoxin to protein stability is investigated by combining theory, experiment and simulation. Although hydrogen bonds are major determinants of protein structure and function, their contribution to protein stability is still unclear and widely debated. The best method so far devised to estimate the contribution of side-chain interactions to protein stability is double mutant cycle analysis, but the interaction energies so derived are not identical to incremental binding energies (the energies quantifying net contributions of two interacting groups to protein stability). Here we introduce double-deletion analysis of 'isolated' residue pairs as a means to precisely quantify incremental binding. The method is exemplified by studying a surface-exposed hydrogen bond in a model protein (Asp96/Asn128 in apoflavodoxin). Combined substitution of these residues by alanines slightly destabilizes the protein due to a decrease in hydrophobic surface burial. Subtraction of this effect, however, clearly indicates that the hydrogen-bonded groups in fact destabilize the native conformation. In addition, molecular dynamics simulations and classic double mutant cycle analysis explain quantitatively that, due to frustration, the hydrogen bond must form in the native structure because when the two groups get approximated upon folding their binding becomes favorable. We would like to remark that 1), this is the first time the contribution of a specific hydrogen bond to protein stability has been measured by experiment; and 2), more hydrogen bonds need to be analyzed to draw general conclusions on protein hydrogen bond energetics. To that end, the double-deletion method should be of help.     

Fidelity of mispair formation and mispair extension is dependent on the interaction between the minor groove of the primer terminus and Arg668 of DNA polymerase I of Escherichia coli

The hydrogen bonding interactions between the Klenow fragment of Escherichia coli DNA polymerase I with the proofreading exonuclease inactivated (KF(-)) and the minor groove of DNA were examined with modified oligodeoxynucleotides in which 3-deazaguanine (3DG) replaced guanine. This substitution would prevent a hydrogen bond from forming between the polymerase and that one site on the DNA. If the hydrogen bonding interaction were important, then we should observe a decrease in the rate of reaction. The steady-state and pre-steady-state kinetics of DNA replication were measured with 10 different oligodeoxynucleotide duplexes in which 3DG was placed at different positions. The largest decrease in the rate of replication was observed when 3DG replaced guanine at the 3'-terminus of the primer. The effect of this substitution on mispair extension and formation was then probed. The G to 3DG substitution at the primer terminus decreased the k(pol) for the extension past G/C, G/A, and G/G base pairs but not the G/T base pair. The G to 3DG substitution at the primer terminus also decreased the formation of correct base pairs as well as incorrect base pairs. However, in all but two mispairs, the effect on correct base pairs was much greater than that of mispairs. These results indicate that the hydrogen bond between Arg668 and the minor groove of the primer terminus is important in the fidelity of both formation and extension of mispairs. These experiments support a mechanism in which Arg668 forms a hydrogen bonding fork between the minor groove of the primer terminus and the ring oxygen of the deoxyribose moiety of the incoming dNTP to align the 3'-hydroxyl group with the alpha-phosphate of the dNTP. This is one mechanism by which the polymerase can use the geometry of the base pairs to modulate the rate of formation and extension of mispairs.     

Structure, energetics, vibrational frequencies and charge transfer of base pairs, nucleoside pairs, nucleotide pairs and B-DNA pairs of trinucleotides: ab initio HF/MINI-1 and empirical force field study

Geometries, interaction energies and vibrational frequencies of base pairs, nucleoside pairs and nucleotide pairs were studied by ab initio Hartree-Fock (HF) method using MINI-1 basis set and empirical Cornell et al. force field (AMBER 4.1). A good agreement was found between HF/MINI-1 and AMBER results. In addition, both methods provide reasonable agreement with available high-level ab initio data. Finally, AMBER potential was used to determine the structure, energetics and vibrational frequencies of B-DNA pairs of trinucleotides. Stabilization energies of clusters are lowered when passing from base pairs to nucleoside pairs, nucleotide pairs and to pairs of trinucleotides. The lowest vibrations of base pairs and nucleoside pairs correspond to intermolecular motions of bases, specifically to buckle and propeller motions. In the case of pairs of larger subunits the lowest vibrations are of intramolecular nature (rotation around glycosidic bond, sugar and phosphate vibration). The spectra of these clusters became more complicated and quasi-degenerate. Intermolecular charge transfer between bases in H-bonded and stacked pairs is negligible, while a significant intramolecular charge transfer was observed.     

The influence of Li+, Na+, Mg2+, Ca2+, and Zn2+ ions on the hydrogen bonds of the Watson-Crick base pairs

The interaction of mono- and divalent metal ions with the nucleic acid base pairs A:T and G:C has been studied using ab initio self-consistent field Hartree-Fock computations with minimal basis sets. Energy-optimized structures of the two base pairs with a final base-base distance of L = 10.35 A have been determined and were further used in calculations on ternary complexes Mn+ - A:B together with previously computed coordination geometries of the cations at adenine (Ade), thymine (Thy), and guanine (Gua). Besides the binding energy of the various metal ions to the base pairs, changes in the stability of the H bonds between Ade and Thy or Gua and Cyt have been determined. Polarization effects of the metal ion on the ligand turned out to increase the binding between complementary bases. Regardless of the metal species, cation binding to Gua N(3) and Thy O(2) leads to a special increase in H-bond stability, whereas binding to Ade N(3) changes the H-bond stability least. Situated in between are the stabilizing effects caused by Gua and Ade N(7) coordination. A remarkable relation between the stability of the H bond and the distance from metal binding site to H bonds was found. This relationship has been rationalized in terms of partial charges of the atoms participating in H bonding, which can reveal the trend in the electrostatic part of total H bond energy. It can be shown that a short distance between coordination site and acceptor hydrogen increases the H-bond strength substantially, while a long distance shows minor effects as supposed. On the other hand, the opposite effect is observed for the influence of the distance between binding site and donor atom. A comparison of our findings with a new model of transition metal ion facilitated rewinding of denatured DNA proposed by S. Miller, D. VanDerveer, and L. Marzilli is given [(1985) J. Am. Chem. Soc. 107, 1048-1055].     

Structure and the energy of base pairing in non-natural bases of nucleic acids: the azaguanine-cytosine and azaadenine-thymine base pairs

Watson-Crick optimized geometries and the energies of base pairing for the natural pairs of nucleic bases: adenine-thymine (AT) and guanine-cytosine (GC) have been recalculated by ab initio methods in order to compare results to those found for the non-natural azaadenine-thymine (AAT) and azaguanine-cytosine (AGC) pairs. Geometry optimizations carried out at the HF/6-31G** level and energies obtained at MP2/6-31G**, show that AAT and AGC have hydrogen bonding patterns similar to the natural AT and GC and that the interaction energies (DeltaH0int) for the former are ca. 7 kcal/mol more stable than the latter. Accordingly, the pairs based on azapurines would be favored with respect to the natural pairs. Some possible explanations why nature does not use extensively the azabases in base pairing are given.     

Stability of XGCGCp, GCGCYp, and XGCGCYp helixes: an empirical estimate of the energetics of hydrogen bonds in nucleic acids

The stabilizing effects of dangling ends and terminal base pairs on the core helix GCGC are reported. Enthalpy and entropy changes of helix formation were measured spectrophotometrically for AGCGCU, UGCGCA, GGCGCCp, CGCGCGp, and the corresponding pentamers XGCGCp and GCGCYp containing the GCGC core plus a dangling end. Each 5' dangling end increases helix stability at 37 degrees C roughly 0.2 kcal/mol and each 3' end from 0.8 to 1.7 kcal/mol. The free energy increments for dangling ends on GCGC are similar to the corresponding increments reported for the GGCC core [Freier, S. M., Alkema, D., Sinclair, A., Neilson, T., & Turner, D. H. (1985) Biochemistry 24, 4533-4539], indicating a nearest-neighbor model is adequate for prediction of stabilization due to dangling ends. Nearest-neighbor parameters for prediction of the free energy effects of adding dangling ends and terminal base pairs next to G.C pairs are presented. Comparison of these free energy changes is used to partition the free energy of base pair formation into contributions of "stacking" and "pairing". If pairing contributions are due to hydrogen bonding, the results suggest stacking and hydrogen bonding make roughly comparable favorable contributions to the stability of a terminal base pair. The free energy increment associated with forming a hydrogen bond is estimated to be -1 kcal/mol of hydrogen bond.     

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

NMR studies of the stable mismatch purine-thymine in the self-complementary d(CGPuAATTTCG) duplex in solution

One- and two-dimensional nuclear Overhauser effect experiments demonstrate that a single hydrogen bond between a T imino proton and purine N3 is sufficient to hold the base pair dPu.dT in d(CGPuAATTTCG) by a Watson-Crick fashion rather than a Hoogsteen type. In addition, the dPu.dT base pair is well stacked with neighboring base pairs. The spin-lattice relaxation measurements at 30 and 35 degrees C of two decamers, d(CGPuAATTTCG) and d(CGAAATTTCG), reveal that the elimination of two single hydrogen bonds of dA.dT base pairs (due to the substitution of adenine for purine) in the sequence results in an increase in the overall imino proton exchange rate from 7 to 36 s-1 at the site of mismatch.     

Sequence dependence for the energetics of dangling ends and terminal base pairs in ribonucleic acid

Stability increments of terminal unpaired nucleotides (dangling ends) and terminal base pairs on the core helixes AUGCAU and UGCGCA are reported. Enthalpy, entropy, and free energy changes of helix formation were measured spectrophotometrically for 18 oligoribonucleotides containing the core sequences. The results indicate 3' dangling purines add more stability than 3' dangling pyrimidines. In most cases, the additional stability from a 3' dangling end on an AU base pair is less than that on a GC base pair [Freier, S.M., Burger, B.J., Alkema, D., Neilson, T., & Turner, D.H. (1985) Biochemistry 22, 6198-6206]. The sequence dependence provides a test for the importance of dangling ends for various RNA interactions. Correlations are suggested with codon context effects and with the three-dimensional structure of yeast phenylalanine transfer RNA. In the latter case, all terminal unpaired nucleotides having stability increments more favorable than -1 kcal/mol are stacked on the adjacent base pair. All terminal unpaired nucleotides having stability increments less favorable than -0.3 kcal/mol are not stacked on the adjacent base pair. In several cases, this lack of stacking is associated with a turn in the sugar-phosphate backbone. This suggests stability increments measured on oligoribonucleotides may be useful for predicting tertiary structure in large RNA molecules. Comparison of the stability increments for terminal dangling ends and base pairs, and of terminal GC and AU base pairs, indicates the free energy increment associated with forming a hydrogen bond can be about -1 kcal/mol of hydrogen bond.     

Probing the role of hydrogen bonds in the stability of base pairs in double-helical DNA

Aromatic stacking and hydrogen bonding between nucleobases are two of the key interactions responsible for stabilization of DNA double-helical structures. The present work aims at defining the specific contributions of these interactions to the stability of individual base pairs in DNA. The two DNA double helices investigated are formed, respectively, by the palindromic base sequences 5'-dCCAACGTTGG-3' and 5'-dCGCAGATCTGCG-3'. The strength of the N==H...N inter-base hydrogen bond in each base pair is characterized from the measurement of the protium-deuterium fractionation factor of the corresponding imino proton using NMR spectroscopy. The structural stability of each base pair is evaluated from the exchange rate of the imino proton, measured by NMR. The results reveal that the fractionation factors of the imino protons in the two DNA double helices investigated fall within a narrow range of values, between 0.92 and 1.0. In contrast, the free energies of structural stabilization for individual base pairs span 3.5 kcal/mol, from 5.2 to 8.7 kcal/mol (at 15 degrees C). These findings indicate that, in the two DNA double helices investigated, the strength of N==H...N inter-base hydrogen bonds does not change significantly depending on the nature or the sequence context of the base pair. Hence, the variations in structural stability detected by proton exchange do not involve changes in the strength of inter-base hydrogen bonds. Instead, the results suggest that the energetic identity of a base pair is determined by the number of inter-base hydrogen bonds, and by the stacking interactions with neighboring base pairs.     

Effects of base substituents on the hydration of B- and Z-DNA: correlations to the B- to Z-DNA transition.

We present a study of how substituent groups of naturally occurring and modified nucleotide bases affect the degree of hydration of right-handed B-DNA and left-handed Z-DNA. A comparison of poly(dG-dC) and poly(dG-dm5C) titrations with the lipotropic salts of the Hofmeister series infers that the methyl stabilization of cytosines as Z-DNA is primarily a hydrophobic effect. The hydration free energies of various alternating pyrimidine-purine sequences in the two DNA conformations were calculated as solvent free energies from solvent accessible surfaces. Our analysis focused on the N2 amino group of purine bases that sits in the minor groove of the double helix. Removing this amino group from guanine to form inosine (I) destabilizes Z-DNA, while adding this group to adenines to form 2-aminoadenine (A') stabilizes Z-DNA. These predictions were tested by comparing the salt concentrations required to crystallize hexanucleotide sequences that incorporate d(CG), d(CI), d(TA) and d(TA') base pairs as Z-DNA. Combining the current results with our previous analysis of major groove substituents, we derived a thermodynamic cycle that relates the systematic addition, deletion, or substitution of each base substituent to the B- to Z-DNA transition free energy.     

基于密度泛函理论研究Watson–Crick碱基对中的氢键特征(英文)

基于密度泛函理论(DFT)研究腺嘌呤、胸腺嘧啶、鸟嘌呤、胞嘧啶以及腺嘌呤胸腺嘧啶碱基对、鸟嘌呤胞嘧啶碱基对。在DFT-B3LYP/6-31G**水平上利用自然键轨道理论分析研究结果显示,互补碱基对的结构和电子特征有利于氢键的形成。本文中讨论几何结构、电子结构、分子轨道和能量对于氢键形成的影响。此研究结果将有助于更好的理解AT和GC碱基对中氢键与它们的结构特性之间的关系。     

Vibrational fluctuations of hydrogen bonds in a DNA double helix with alternating TA type inserts.

The Green's function technique is applied to a study of breathing modes in a DNA double helix which contains a region of different base pairs from the rest of the double helix. The calculation is performed on an alternating poly(dC-dG).poly(dC-dG) helix in the B conformation with four consecutive base pairs replaced by a model of a biological promoter region with four alternating T-A,A-T base pairs, henceforth referred to as (TATA)2. The average stretch of interbase hydrogen bonds is found to be amplified around the insert. This is likely related to the (TATA)2 insert having a lower stability against hydrogen bond melting than the two semi-infinite poly(dC-dG).poly(dC-dG) helices. The insert region may be considered to be a site of enhanced tendency to melt in such a helix. The results show that an alternating AT insert of four base pairs has a larger average hydrogen bond stretch inside and outside the insert region than the average hydrogen bond stretch inside and outside an insert of four consecutive A-T base pairs, henceforth referred to as (AAAA).(TTTT). Calculations are performed which show that the enhancement of the average hydrogen bond stretch around an alternating TA type insert is greatly dependent upon the local modes and not the inband modes. The amount of local mode enhanced average stretch is explored as a function of insert size.     

The thermal stability of RNA duplexes containing modified base pairs placed at internal and terminal positions of the oligoribonucleotides

The presence of various modifications within oligomers changes their thermodynamic stability. To get more systematic data, we measured effects of 5- and 6-substituted uridine on thermal stability of (AUCU(Mod.)AGAU)2 and (AUCUAGAU(Mod.))2. Collected results lead to the following conclusions: (i) 5-halogenated and 5-alkylated substituents of the uridine affect thermal stability of the RNA duplexes differently. Moreover, the 5-fluorouridine changes stability of the RNA duplexes opposite to remaining 5-halogenouridines; (ii) for oligomers containing 5-chloro, 5-bromo or 5-iodouridine stronger hydrogen bond formed between oxygen-4 of the 5-halogenated uracil and 6-amino group of the adenine is presumably responsible for stabilizing effect; (iii) placing of A-U(5R) base pairs closer to the end of the duplex enhance thermal stability relatively to oligomer with central position of this base pair; (iv) the effects of 5-substituents are additive, particularly for substituents which stabilize RNA duplexes; (v) 6-methyluridines (N1 and N3 isomers) as well as 3N-methyluridine present at internal position of A-U(Mod.) inhibit duplexes formation; (vi) 6-methyluridines (N1 and N3 isomers) as well as 3N-methyluridine placed as terminal base pairs stabilize the duplexes mostly via 3'-dangling end effect.     

Synthesis and characterization of 8-methoxy-2'- deoxyadenosine-containing oligonucleotides to probe the syn glycosidic conformation of 2'-deoxyadenosine within DNA.

The synthesis of 8-methoxy-2'-deoxyadenosine (moA) protected at N6 as an N,N-dimethylformamidine derivative and incorporation of the modified nucleoside into oligodeoxynucleotides via the phosphoramidite method are described. UV thermal denaturation studies were conducted on duplexes containing moA:G, moA:C and moA:T base pairs to determine the thermodynamic stability of duplexes containing moA relative to their adenosine (A)-containing counterparts. In the case of moA:G base pairs the effect of moA substitution is sequence dependent. In A:G mismatch-containing sequences, which have been shown by structural characterization to have a syn conformational preference at the glycosidic bond of A, moA substitution results in stabilization of the duplex. In contrast, in sequences where the A in the A:G mismatch has been shown to prefer the anti conformation moA substitution is destabilizing to the duplex. Thus moA may be a useful probe for investigating the conformational preferences of the N-glycosidic bond of adenosine within DNA. In addition, moA nucleoside is more resistant to acid-catalyzed depurination than previously described 8-bromo-2'-deoxyadenosine, allowing for facile incorporation into oligonucleotides via automated solid phase DNA synthesis.     

Substitution, cooperative, and solvent effects on π pnicogen bonds in the FH2P and FH2As complexes

Ab initio calculations have been carried out to study the substitution effect on the π pnicogen bond in ZH(2)P-C(2)HM (Z?=?H, H(3)C, NC, F; M?=?H, CH(3), Li) dimer, cooperative effect of the π pnicogen bond and hydrogen bond in XH-FH(2)Y-C(2)H(4) (X?=?HO, NC, F; Y?=?P and As) trimer, and solvent effect on the π pnicogen bond in FH(2)P-C(2)H(2), FH(2)P-C(2)H(4), FH(2)As-C(2)H(2), and FH(2)As-C(2)H(4) dimers. The interaction energy of π pnicogen bond increases in magnitude from -1.51?kcal?mol(-1) in H(3)P-C(2)H(2) dimer to -7.53?kcal?mol(-1) in FH(2)P-C(2)HLi dimer at the MP2/aug-cc-pVTZ level. The π pnicogen bond is enhanced by 12-30?% due to the presence of hydrogen bond in the trimer. The π pnicogen bond is also enhanced in solvents. The natural bond orbital analysis and symmetry adapted perturbation theory (SAPT) were used to unveil the source of substitution, cooperative, and solvent effects.