Calculating Distortions of Short DNA Duplexes with Base Pairing Between an Oxidatively Damaged Guanine and a Guanine

DNA is constantly being oxidized, and oxidized DNA is prone to mutation; moreover, guanine is highly sensitive to several oxidative stressors. Several oxidatively damaged forms of guanine—including 2,2,4-triamino-5(2H)-oxazolone (Oz), iminoallantoin (Ia), and spiroiminodihydantoin (Sp)—can be paired with guanine, and cause G:C-C:G transversions. Previous findings indicate that guanine is incorporated more efficiently opposite Oz than opposite Ia or Sp, and that these differences in efficiency cannot be explained by differences in the stabilities of G:Oz, G:Ia, and G:Sp base pairs calculated ab initio. Here, to explain previous experimental result, we used a 3-base-pair model DNA duplex to calculate the difference in the stability and the distortion of DNA containing a G:Oz, G:Ia, or G:Sp base pair. We found that the stability of the structure containing 5' and 3' base pairs adjacent to G:Oz was more stable than that containing the respective base pairs adjacent to G:Ia or G:Sp. Moreover, the distortion of the structure in the DNA model duplex that contained a G:Oz was smaller than that containing a G:Ia or G:Sp. Therefore, our discussion can explain the previous results involving translesion synthesis past an oxidatively damaged guanine.


Introduction
Mutations of genomic information are commonly caused by oxidative stress, and oxidative stressors have been implicated in aging, carcinogenesis, and many other diseases. Guanine is particularly sensitive to several oxidative stressors because it has the lowest oxidation potential among the four bases. Several different oxidative stressors can cause G:C-T:A and G:C-C:G transversions [1]. In fact, G:C-T:A and G:C-C:G transversions caused by passive smoking are detected at high frequencies in codon 12 of the K-ras gene [2].

Scheme 1. Oxidation products of guanine and 8-oxoG.
Guanidinohydantoin (Gh) and spiroiminodihydantoin (Sp) can be formed via oxidation of 8-oxoG (Scheme 1). Gh is a principal product of 8-oxoG oxidation under acidic conditions, and Sp is a major product under basic conditions [8][9][10][11]. Gh isomerizes to iminoallantoin (Ia) [12], but whether the predominant isomer is Gh or Ia during DNA polymerization is unknown [12]. Based on our previously published calculations of the stability of G:Ia and G:Gh base pairs, Gh is predicted to tautamerize to Ia when guanine is inserted opposite Gh/Ia [13]. Hence, in this study, we consider only Ia. Incorporation of guanine opposite Oz is more efficient than that opposite Gh/Ia or Sp, and translesion synthesis past Oz is also more efficient than that past Gh/Ia or Sp [5,7,14]. However, depending on the number of hydrogen bonds, a G:Oz base pair was less stable than either a G:Ia or the G:Sp base pair in our calculated results (Figure 1) [13]. Generally, since replicative DNA polymerases are highly susceptible to geometric distortions in DNA, they have difficulty synthesizing past distorting DNA damage [15]. The previous study showed that the geometry of DNA containing nonplanar DNA damage readily destabilizes and distorts relative to DNA that includes natural Watson-Crick base pairs [16]. Therefore, to resolve the discrepancy between experimental results and calculated results, we focused this study on differences among the distortion of DNA containing these lesions. As depicted in Figure 1, the G:Oz base pair has no sp 3 carbon and is planar, while G:Ia and G:Sp base pair have a sp 3 carbon and are nonplanar. Ia and Sp break the planarity of the adjacent base pairs due to a sp 3 carbon, and as a result, DNA containing Ia or Sp is easily distorted compared to that containing Oz.
In this study, we tried to explain the experimental results by calculating the differences in the stabilities and distortions of DNA containing a G:Oz, G:Ia, or G:Sp base pair relative to DNA containing a G:C. We use 3-base-pair model DNA duplexes that each contained a central G:Oz, G:Ia, or G:Sp base pair as a way to calculate the differences in stability and distortion.

The Minimized Structures of DNA Duplex Containing an X:G Base Pair
We previously calculated stabilization energies of oxidatively guanine damage (Oz, Ia or Sp) paired with guanine [13]. Since Ia and Sp have a sp 3 carbon, they have S and R stereoisomers. Our calculated results showed that the stabilization energies of the base pairs, "G:X (where X = Oz, S-Ia, R-Ia, S-Sp, or R-Sp)", are ordered as follows: G:S-Ia > G:R-Ia >> G:R-Sp > G:S-Sp >> G:Oz [13]. Here, we investigated differences in stabilities of duplexes containing one G:X base pair in the active site of DNA polymerase; each Pol β-DNA complex containing a G:X base pair was built by modification of the Pol β ternary complex (PDB entry 1BPY) [17]. We used the structure of G:X base pair formed the most stable base pair with guanine based on our previous report [13].
Quantum mechanical (QM) calculations are more quantitative than molecular mechanics (MM) calculations. In this study, we used QM calculations only to optimize the most important part of each G:X base pair. However, the Pol β-DNA complex model had too many atoms to use QM calculations for minimizing the model in an experimentally acceptable amount of time. Therefore, we used MM calculations, which are quantitatively inferior to QM calculations, to minimize the Pol β-DNA model. Only the most important part of each G:X base pair was optimized by QM calculation, and the geometries of Pol β-DNA complex model were minimized via MM calculations. We focused on the G:X base pair and each of two adjacent base pairs, one 5ꞌ and one 3ꞌ to the central G:X base pair, in each minimized structure (Figures 2 and 3). In Figure 2, the central G:X base pair was designated "G 2 X 2 "; the 5ꞌ A:T base pair was designated "A 1 T 1 ", and the 3ꞌ G:C base pair was designated "G 3 C 3 ". Minimized geometries of A 1 T 1 , G 2 X 2 and G 3 C 3 base pairs are depicted in Figure 3. We wanted to compare the impact of each X 2 on the calculated energies of the duplexes. However, the X 2 variants differ from one another in the number of atoms; therefore, the energies of the structures containing different X 2 could not be compared directly. Therefore, we focused on the energies of the elements common to each duplex (i.e., A 1 T 1 and G 3 C 3 ) to compare among the energies of the duplexes. By calculating the energies of the common elements and excluding each "G 2 X 2 " base pair within each structure, we assessed the destabilization caused by the oxidatively damaged guanine.

The Destabilization Energies of the Base Pairs on 5ꞌ-and 3ꞌ-Side of X
In this section, in order to evaluate the destabilization of a DNA duplex caused by particular forms of oxidatively damaged guanine, we calculated the destabilization energies of elements common to each structure and excluded each "G 2 X 2 " base pair.  Figure 3. A:T base pair on the 5ꞌ-side of X was designated "A 1 T 1 ", G:X base pair was designated "G 2 X 2 ", and G:C base pair on the 3ꞌ-side of X was designated "G 3 C 3 ". The destabilization energies of "A 1 T 1 " (∆E 1 ), "G 3 C 3 " (∆E 3 ), and "A 1 T 1 + G 3 C 3 " (∆E 1 + 3 ) were calculated ab initio as the parts common to each model duplex; each G 2 X 2 base pair was excluded from the calculations. We calculated the destabilization energies of each "A 1 T 1 " (∆E 1 ) on the 5ꞌ-side of each X 2 . The calculated ∆E 1 value associated with duplexes that contained Oz, S-Ia, R-Ia, S-Sp, or R-Sp were defined by Equation (1) Table 1. Destabilization energies (kcal/mol) of "A 1 T 1 " (∆E 1 ), "G 3 C 3 " (∆E 3 ), and "A 1 T 1 + G 3 C 3 " (∆E 1 + 3 ), each value was calculated with minimized geometries. Using kinetic parameters, Kornyushyna et al., demonstrated that incorporation of guanine opposite Gh/Ia is more efficient than that opposite Sp [14]. Our calculated data indicated that "A 1 T 1 " on the 5ꞌ-side of Ia was more stable than that of Sp regardless of the stereoisomer, and our data corresponded with these published results. However, our data indicated that "A 1 T 1 " on 5ꞌ-side of Oz was less stable than "A 1 T 1 " 5ꞌ to R-Ia, and this finding could not account for the previous experimental finding that translesion synthesis past Oz is more efficient than that past Ia [5]. Therefore, we calculated the destabilization energies for each G 3 C 3 base pair on the 3ꞌ-side of each X 2 .

The Destabilization Energies of the Base Pairs on the 3ꞌ-Side of X
The destabilization energies (in vacuo; ∆E 3 DFT and in water; ∆E 3 SCRF ) of "G 3 C 3 " were defined by Equation (2) (Experimental Section), and were calculated individually ( Table 1). The values of ∆E 3 DFT and ∆E 3 SCRF were ordered as follows with regard to each X 2 : Surprisingly, the calculated values for ∆E 3 DFT and ∆E 3 SCRF of Oz were the lowest among the values for any X 2 , and "G 3 C 3 " on 3ꞌ-side of Oz was the most stable. The stabilities of the "G 3 C 3 " base pairs on the 3ꞌ-side of X 2 were consistent with previous experimental results indicating that translesion synthesis past Oz is more efficient than that past Ia or Sp [5,14]. The "G 3 C 3 " base pairs on 3ꞌ-side of S-Sp was more stable than that of R-Ia; consequently, it was not possible to fully explain the previous experimental results based on the stabilities of G 3 C 3 base pairs on the 3ꞌ-side of X 2 . Consequently, we then considered the stability of the larger structure by examining the energies of the structure containing both base pairs (5ꞌ and 3ꞌ) adjacent to each X 2 .

The Destabilization Energy of Each Model Duplex was Calculated by Including both Adjacent
Base Pairs and Excluding the Central "G 2 X 2 " Base Pair The destabilization energy of A 1 T 1 + G 3 C 3 (in vacuo; ∆E 1 + 3 DFT and in water; ∆E 1 + 3 SCRF ) was calculated for each structure, including both adjacent base pairs and excluding the central "G 2 X 2 " base pair. These values were calculated with Equation (3) "A 1 T 1 + G 3 C 3 " with Oz was more stable than that with any other X 2 ; this finding indicated that a DNA structure containing G 2 Oz 2 was more similar to DNA containing G 2 C 2 than was that with any other G 2 :X 2 base pair. As mentioned in the introduction; the geometry of DNA that includes nonplanar DNA damage is easy to destabilize and distort relative to that containing natural Watson-Crick base pairs [16]; and DNA polymerases have difficulty synthesizing through distorted DNA [15]. Therefore; our calculated values could account for the previous experimental findings that translesion synthesis across Oz is more efficient than that across Ia or Sp [5,14].
Notably, the calculated ∆E 1 + 3 DFT values for R-Ia and S-Sp were not consistent with previous results [14]. Actual experiments involving DNA polymerization have been performed in aqueous solution [5,14]; therefore, the stability of DNA duplex during DNA polymerization is attributed to the energy in water (∆E 1 + 3 SCRF ). ∆E 1 + 3 SCRF for R-Ia was lower than ∆E 1 + 3 SCRF for S-Sp, which can account for the previous result [14]. Furthermore, the calculated stabilities of "A 1 T 1 + G 3 C 3 " with S-Sp or R-Sp could explain why bypass of S-Sp is more efficient than bypass of R-Sp in vivo [18].
In conclusion, although calculating either energies of "A 1 T 1 " on the 5ꞌ-side of X 2 or energies of "G 3 C 3 " on the 3ꞌ-side of X 2 was not sufficient, calculating energies of the structure while including both base pairs, "A 1 T 1 + G 3 C 3 ", adjacent to "G 2 X 2 " was sufficient to evaluate the stability of the DNA duplex containing an G:Oz, G:Ia, or G:Sp base pair. In the next section, we used parameterized distortion of DNA to examine the effects of the non-planarity of Oz, S-Ia, R-Ia S-Sp, or R-Sp on A 1 T 1 base pairs and G 3 C 3 base pairs in DNA duplexes that each contain one central G:X base pair.

The Degree of Distortion from DNA Duplex Containing a G:C Base Pair
Based on the images in Figure 3, we received the visual impression that there were differences in the degree of distortion among the DNA duplexes that contained C, Oz, S-Ia, R-Ia S-Sp, or R-Sp. Because abstract visual impressions were insufficient for detailed comparisons among the structures, it was necessary to parameterize actual values for those visual impressions. In this section, we examined differences in the degree of distortion caused by each G:X base pair to parameterize the distortion of each DNA duplex.

The Degree of Distortion at the 5ꞌ-Side of X
In order to determine the degree of distortion at the 5ꞌ-side of each G 2 X 2 base pair, we calculated the dihedral angles (θ) between G 2 and A 1 (θ (G 2 -A 1 )) or G 2 and T 1 (θ (G 2 -T 1 )) or A 1 and T 1 (θ (A 1 -T 1 )) (red arrows in Figure 4). Values for θ (G 2 -A 1 ), θ (G 2 -T 1 ), and θ (A 1 -T 1 ) were each calculated individually via Equations (5)-(7), respectively (Experimental Section). The degree of distortion (δ 1 ) on the 5ꞌ-side of G 2 X 2 was defined as the sum of θ (G 2 -A 1 ), θ (G 2 -T 1 ), and θ (A 1 -T 1 ) (Equation (11) in Experimental Section). In Table 2, the values of δ 1 were ordered as follows with regard to X 2 variants:  A 1 ), θ (G 2 -T 1 ) and θ (A 1 -T 1 ) showed red arrows, and the calculated dihedral angle θ (G 2 -G 3 ), θ (G 2 -C 3 ) and θ (G 3 -C 3 ) showed blue arrows. Table 2. Dihedral angles θ (G 2 -A 1 ), θ (G 2 -T 1 ), and θ (A 1 -T 1 ) (red arrows in Figure 4), and the degree of distortion δ 1 . Based on the values of δ 1 , the distortion of "A 1 T 1 " base pair adjacent to R-Ia is the smallest, and this structure was similar to the structure of the "A 1 T 1 " derived from DNA duplex containing a G 2 C 2 ( Ia and for Sp. The C1ꞌ methyl of S-Ia or S-Sp was oriented in the direction of the "A 1 T 1 " base pair (Figure 3c,e); therefore, the effects of the S configurations on 5ꞌ-side of X 2 were larger than those of the R configurations. In summary, δ 1 values could not explain previous experimental results, nor could they completely account for ∆E 1 values.

The Degree of Distortion at the 3ꞌ-Side of X
As with 5ꞌ-side of X 2 (in Section 2.3.1), we calculated the dihedral angle θ between G 2 and G 3 (θ (G 2 -G 3 )), or G 2 and C 3 (θ (G 2 -C 3 )), or G 3 and C 3 (θ (G 3 -C 3 )) (blue arrows in Figure 4). Values for θ (G 2 -G 3 ), θ (G 2 -C 3 ), and θ (G 3 -C 3 ) were each calculated individually via Equations (8)-(10), respectively (Experimental Section). The degree of distortion (δ 3 ) at the 3ꞌ-side of X 2 was defined as the sum of θ (G 2 -G 3 ), θ (G 2 -C 3 ), and θ (G 3 -C 3 ) (Equation (12) in Experimental Section). All numeric values are listed in Table 3. For the X 2 residues, the order of values for δ 3 was as follows: R-Ia > R-Sp > S-Sp > S-Ia > C > Oz. Table 3. Dihedral angles θ (G 2 -G 3 ), θ (G 2 -C 3 ), and θ (G 3 -C 3 ) (blue arrows in Figure 4), and the degree of distortion δ 3 . When comparing among the X 2 variants, the "G 3 C 3 " adjacent to Oz was the least distorted based on the values of δ 3 , and the most similar in structure to the "G 3 C 3 in a DNA duplex containing a G 2 C 2 (Figure 3a,b). Additionally, the calculated ∆E 3 with Oz was lowest, and the order of X 2 variants with regard to calculated ∆E 3 values-S-Sp > S-Ia > Oz (Table 1)-was consistent with the order of δ 3 values. As was the case with δ 1 (showed in section 2.3.1), the destabilization of the G 3 C 3 base pair could be accounted for by the distortion caused by the G 2 X 2 base pair. In contrast to δ 1 values, δ 3 values associated with R configurations were larger than those associated with S configurations for both Ia and Sp; moreover, the C1ꞌ methyl of R-Ia or R-Sp was oriented in the direction of the "G 3 C 3 " base pair (Figure 3d,f). Consequently, the C 3 base of R-Ia or R-Sp was particularly distorted by the C1ꞌ methyl of R-Ia or R-Sp, respectively. Taken together, δ 3 values could not account for ∆E 3 values or previous experimental results, as was the case for δ 1 values.

The Total Degree of Distortion at the 5ꞌ-Side and 3ꞌ-Side of X
The total degree of distortion ("δ 1 + δ 3 ") at the 5ꞌ-and 3ꞌ-side of X 2 was defined as the sum of δ 1 and δ 3 ( Table 4). The "δ 1 + δ 3 " values were ordered as follows with regard to X 2 variants: R-Sp > S-Sp > R-Ia > S-Ia > C > Oz. This order was consistent with the order of the calculated energies (∆E 1 + 3 ).
Although the value of "δ 1 + δ 3 " for Oz was smaller than that for C, the structure of "A 1 T 1 + G 3 C 3 " adjacent to "G 2 Oz 2 " was less stable than that to"G 2 C 2 " according to ∆E 1 + 3 : The stabilities of the structure of "A 1 T 1 + G 3 C 3 " were only partially explained by the degree of distortion.

The Difference Between the Degree of Distortion for Oxidatively Damaged Guanine and that for C
We tried to eliminate the discordance between the distortion and the stability in "G 2 Oz 2 " and "G 2 C 2 ". Instead of a simple comparison among the values of the degree of distortion for any X 2 , we calculated the difference between the value of δ 1 + δ 3 for C and that for Oz, S-Ia, R-Ia, S-Sp or R-Sp as well as the destabilization energy. We found that the values of δ 1 + δ 3 for S-Ia or R-Ia differed from that for C by 30.7° or 35.4°, respectively; in contrast, there was only slight difference of 4.0° between that for Oz and that for C. The distortion of the structure of "A 1 T 1 + G 3 C 3 " adjacent to "G 2 Oz 2 " was smaller than that adjacent to "G 2 S-Ia 2 " or "G 2 R-Ia 2 ", which seemed consistent with visual impressions.
Additionally, the values of δ 1 + δ 3 for S-Sp or R-Sp differed from that of C by 81.1° or 103.3°, respectively. These values indicated that the structures of "A 1 T 1 + G 3 C 3 " adjacent to "G 2 S-Sp 2 " or "G 2 R-Sp 2 " were the most distorted from that adjacent to "G 2 C 2 "; these findings seemed consistent with the visual impression of the images in Figure 3. These discussions accounted for the differences among the calculated energies (∆E 1 + 3 ) of the structure of "A 1 T 1 + G 3 C 3 " adjacent to "G 2 X 2 ".
In addition to the values of δ 1 + δ 3 , the values of δ 1 or δ 3 for Oz, S-Ia, R-Ia, S-Sp, or R-Sp were compared with that for C. The distortions of each X 2 -containing structure compared to the value of δ 1 for C was ordered as follows with regard to X 2 variants: S-Sp > R-Ia > R-Sp > S-Ia > Oz: That of δ 3 was ordered as follows: R-Ia > R-Sp > S-Sp > Oz > S-Ia. Therefore, the distortion for each X 2 compared the value of δ 1 or δ 3 for C was not consistent with the visual impression, unlike the values of δ 1 + δ 3 , which were consistent with the visual impressions.
In summary, the previous experimental results were only partially explained by either δ 1 or δ 3 (Sections 2.3.1 and 2.3.2); however, "δ 1 + δ 3 " compensated for the weaknesses of each single measure by summing δ 1 and δ 3 . Consequently, "δ 1 + δ 3 " was able to account for the previous experimental results and ∆E 1 + 3 . According to our data, the structure of "A 1 T 1 + G 3 C 3 " adjacent to "G 2 X 2 " was distorted by the effects of the non-planarity of each lesion; in other words, the instability of a DNA duplex containing oxidatively damaged guanine was determined by the non-planarity of each oxidatively damaged guanine. Thus, it seemed reasonable that the projected distortion of a DNA duplex containing a G:Oz base pair would be the smallest.

Molecular Modeling
In this study, we constructed models of Pol β-DNA complexes containing a G:X (where X = C, Oz, S-Ia, R-Ia, S-Sp, or R-Sp) base pair using the structure (PDBID:1BPY) [17]. The details of method used to construct these models are presented below. In each model, the incoming nucleotide (dCTP) was replaced with dATP, and the template G for the incoming nucleotide was replaced with T. Each complex containing Pol β and a G:X base pair was built by replacing the G:C base pair at the 3ꞌ terminus of the template-primer with G:X base pair optimized previously [13]. The amino acid residues and other DNA base sequences were unchanged. The geometry of each atom was minimized at the OPLS2005/water level using Macromodel 9.0 (Schrödinger LLC, New York) with the fixed G:X base pair.

ab Initio Calculations
For each minimized structure, all atoms were removed from the Pol β-DNA model except for (1) constituents of the bases of each G:X base pair, (2) constituents of the bases of each base pair adjacent to the G:X base pair, (3) the 2-deoxyribose C1ꞌ carbon, and (4) the C1ꞌ H. Two H atoms were then attached to the C1ꞌ methine to complete the N-methylated nucleobases (Figure 2) [5]. As mentioned in Section 2.1, the designation "A 1 T 1 ", "G 2 X 2 ", and "G 3 C 3 " were used for the consecutive base pairs (Figures 2 and 3).
Gaussian 03 (Gaussian Inc., Wallingford, CT, USA) [19] was used to calculate the destabilization energy of "A 1 T 1 " (∆E 1 ) of each G:X complex in vacuo at the B3LYP/6-31G** level. Moreover, to estimate the energy of each "A 1 T 1 " base pair in water, the Onsanger reaction field model and a dielectric constant of 78.39 were used to calculate the SCRF value of each of these base pairs. The destabilization energies of each "G 3 C 3 " (∆E 3 ) were calculated in vacuo and in water, just as ∆E 1 values were calculated. The destabilization energies of "A 1 T 1 + G 3 C 3 " (∆E 1 + 3 ) were calculated in vacuo and in water. The calculated destabilization energies (∆E 1 , ∆E 3 , and ∆E 1 + 3 ) are defined in Equations (1)-(3).

Conclusions
Previous experimental results show that incorporation of guanine is more efficient when the template base is Oz than when it is Ia or Sp; however, the difference between Oz and Ia or Sp is not explained by the stabilities of G:Oz, G:Ia, and G:Sp as calculated ab initio. Thus, we focused on differences among G:Oz, G:Ia, and G:Sp with regard to planarity of base pair, and tried to explain the previous experimental results based on differences in distortion of DNA containing these DNA lesion.
Using models of 3-base-pair DNA duplexes that each contained a variable central "G 2 X 2 " base pair, we calculated the destabilization energies of the parts common to each model relative to energy of the structure derived from DNA containing a central G 2 C 2 base pair. For the five X 2 variants, we found that the stabilities of the "A 1 T 1 " base pairs on the 5ꞌ-side of X 2 were ordered as follows: R-Ia > Oz > S-Ia > S-Sp > R-Sp, and that the stabilities of the "G 3 C 3 " base pairs on the 3ꞌ-side of X 2 were ordered as follows: Oz > S-Ia > S-Sp > R-Ia > R-Sp. These data could not explain the previous experimental results. We also calculated the stability of the structure including both adjacent base pairs and excluding the central "G 2 X 2 " ("A 1 T 1 + G 3 C 3 "); we found that, for the five X 2 variants, the stability of "A 1 T 1 + G 3 C 3 " was ordered as follows: Oz > S-Ia > R-Ia > S-Sp > R-Sp. The previous experimental results could be explained by considering the stability of the structure when accounting for both the 5ꞌ and 3ꞌ sides of X 2 .
Additionally, we parameterized the distortion of DNA duplex in order to examine the effects of Oz, S-Ia, R-Ia S-Sp, or R-Sp. The degree of distortion at the 5ꞌ-side of X 2 (δ 1 ) or at the 3ꞌ-side of X 2 (δ 3 ) could not explain the stability of "A 1 T 1 " or "G 3 C 3 ", respectively, or the previous experimental results. However, compared with the structure of DNA containing a central G 2 C 2 base pair, the total degree of distortion ("δ 1 + δ 3 ") at the 5ꞌ and 3ꞌ sides of X 2 were ordered as follows for the five X 2 variants: R-Sp >> S-Sp >> R-Ia > S-Ia >> Oz. Values "δ 1 + δ 3 " coincided with the stabilities of "A 1 T 1 + G 3 C 3 ", and either calculation could account for the previous experimental results. Taken together, our discussions indicated that DNA containing a nonplanar lesion was susceptible to distortion and to destabilization, and consequently the bypass efficiency of DNA polymerase at the more nonplanar lesion was reduced relative to that at the more planar lesion. We showed that calculating energies of the structure accounting for both base pairs adjacent to G:Oz, G:Ia, or G:Sp was sufficient to evaluate the stability of the DNA duplex containing any of these base pairs. The distortion of the structure of DNA a duplex containing a central G:Oz was smaller than that of DNA containing a G:Ia or G:Sp; therefore, a G:Oz base pair was more easily bypassed by DNA polymerase than was a G:Ia or G:Sp base pair.
In this study, to minimize by MM calculation, we used the models contained not only the residues of protein and base pairs but also the sugar and phosphate components. There are not, however, the sugar and phosphate components in the DNA models used QM calculation of the stability, because models containing these components cannot calculate in a realistically acceptable amount of time due to many atoms. Therefore, to explain the more accurate effect of distortion caused these DNA damages, it will need to be considered how distortion within the bases can lead to changes in sugar conformation in the future.