Molecular Dynamics Study of Oligonucleotides Containing Difluorotoluene
Molecular Dynamics Study of Oligonucleotides Containing Difluorotoluene
Molecular Dynamics Study of Oligonucleotides Containing Difluorotoluene
6891
Abstract: Extended molecular dynamics (MD) and thermodynamic integration (MD-TI) calculations have been
used to determine the structural and energetic changes in DNA that accompany the replacement of thymine
(T) by the nonnatural isostere difluorotoluene (F). In a duplex DNA oligonucleotide, it is found that the TfF
mutation leads to only small changes in the average structure, but to important alterations in flexibility, hydration,
and recognition properties. The TfF mutation in the Watson-Crick or Hoogsteen position of a pyrimidine
purinepyrimidine type DNA triplex does not lead to dramatic changes in the general structure of the triplex,
but again, detailed analysis shows some alterations in flexibility, hydration, and recognition properties. MDTI calculations on the TfF mutation in duplex DNA reproduce the experimentally determined free energy
differences with good accuracy, and detailed analyses of the trajectories have enabled us to rationalize these.
Finally, MD-TI simulations have been used to predict the changes in stability of a triplex due to a TfF
mutation in either the Watson-Crick or Hoogsteen-binding pyrimidine strands. We predict that in either case
the mutation will reduce stability, being most unfavorable in the Watson-Crick strand.
Introduction
DNA replication by DNA polymerase is a high-fidelity
process with an error rate of the order of one mismatch per 104
to 105 bases.1,2 It might be thought that this fidelity has its origins
in the specificity of Watson-Crick hydrogen bonding, but it
has been argued that this is not the case. The free energy
difference between matched and mismatched terminal base pairs
in solution is estimated3 to range between 0.2 and 0.4 kcal
mol-1. This difference is too low to account for the observed
error rate, suggesting that shape-complementarity must also play
an important role in ensuring the fidelity of replication.4
However, other workers5 have argued that within the environment of the polymerase, the free energy difference between
matched and mismatched base pairs may be much greater due
to reduced solvation and so the hydrogen-bonding argument is
sufficient.
These competing hypotheses have been tested using the
thymine mimic difluorotoluene (F, Figure 1). F was designed
as a nonpolar homologue of thymine,6 which lacks H-bonding
capability, even in apolar solvents such as chloroform.7 Despite
its nonpolar nature, it was found that DNA polymerase I would
Cubero et al.
Figure 1. Chemical structure of 1,3-difluorotoluene (F) and thymine bound to adenine in the Watson-Crick and Hoogsteen orientations.
system
SS (3 steeps)
SS (3 steeps)
SS (3 steeps)
SS (3 steeps)
SS (5 steeps)
duplex d(AT)
duplex d(AT)
duplex d(AF)
triplex d(AT-T)
triplex d(AT-T)
triplex d(AF-T)
triplex d(AF-T)
triplex d(AT-F)
triplex d(AT-F)
structure
length of MD
psb
420 ps
420 ps
840 ps
420 ps
420 ps
420 ps
840 ps
420 ps
420 ps
420 ps
420 ps
840 ps
420 ps
840 ps
Figure 3. RMS deviations in the trajectories of the d(AAGAAAGAAAAG) duplex containing T or F bound to the adenine in position
6. The continuous line shows the deviation with respect to the average
structure obtained from the same trajectory. The dashed line shows
the deviation with respect to the average structure from the other
trajectory.
change associated with a particular mutation allows us to obtain a good
estimate of the statistical errors in the free energy estimates.
Trajectories for MD-TI calculations were performed considering
long-range contributions as introduced by the PME method, and
considering all inter- and nonbonded intramolecular contributions. All
the technical details of simulations used for TI calculations are identical
with those used in standard MD simulations (see above).
-1.4(0.6)
-8.5((4.9)
3.5(0.2)
6.9(2.3)
30.0(1.5)
-11.1(4.5)
116 (35)
40.7(6.2)
22.1(1.7)
12.8(0.9)
a
Averages are done using the last 5 ns of the trajectory of the duplex
with the d(AF) step and the last 1 ns of the trajectory for the control
duplex.
Cubero et al.
contain a TfF mutation: d(AF-T) in 4 and d(AT-F) in 5.
In all cases MD simulations produce stable trajectories over
the nanosecond time scale at room temperature and nearphysiological conditions, suggesting that triplexes containing a
single TfF mutation can be stable in physiological conditions.
The triplex conformations remain essentially B-type. Thus,
the RMS deviations from a canonical B-type model average
1.3 while the RMS deviations from a canonical A-type model
lie in the range 2.3-2.4 . In a manner similar to the duplexes,
convergence of the three trajectories is evident from comparison
of RMS deviations. This is shown in Table 3, where the cross
RMS deviations between triplexes containing F and T are 1.01.1 . These values are similar in magnitude to the RMS
deviations found between individual structures from a trajectory
and the MD-averaged structure derived from it.
The TfF mutation does not have any major impact on the
helical characteristics of the triplexes, which are those expected
for a B-type triplex20,24,2524,25 (see Table 4). In all the cases the
twist is around 30 with a rise around 3.3-3.4 . No change
in groove dimensions or sugar puckering (average phase angle
around 110-112) is apparent as a consequence of the TfF
mutation either in the Watson-Crick or Hoogsteen positions.
The general recognition and hydration characteristics of the
triplex remain unaltered upon the TfF mutation in the WatsonCrick position. However, such a mutation weakens the region
of negative potential in the minor groove, and a partial disruption
of the spine of hydration (see Figure 7). There is also a slight
decrease in the apparent density of water in the minor part of
the major groove (for nomenclature see ref 20). As for the
duplexes (see above), these observations can be accounted for
by considering two factors: first the loss of the lone pair at
position 2 of the pyrimidine, resulting from the TfF mutation
(see Figure 5), and second the breathing motions in the AF
pair (see below).
The TfF mutation in the Hoogsteen position does not lead
to any important change in the electronegativity in the grooves
or in their apparent water density (see Figure 7). The lack of
difference in the minor part of the major groove is particularly
interesting, since it should experience a partial loss in H-bonding
capabilities as a result of the TfF mutation. Comparing with
the results for triplex 4, this suggests that the integrity of the
spines of hydration is more sensitive to the existence of
breathing and opening movements (see ref 13 for nomenclature) than to the partial loss of H-bonding capability.
Breathing movements are detected in triplexes containing F
(4 and 5), while they are not found in the reference triplex 3.
The triplex with the TfF mutation in the Watson-Crick
position (4) shows reversible partial opening and breathing
movements (see ref 13 for nomenclature), which lead to partial
or total loss of WC hydrogen bonds due to the displacement of
F to the major groove. Analysis of the trajectory (see Supporting
Information for details) suggests free energy changes of 0.91.3 kcal/mol for partial opening, and 2.0-2.4 kcal/mol for
breathing. These values are only around 0.5 kcal/mol higher
than those found for a duplex DNA with a single TfF
mutation,13 which suggests that the breathing of WatsonCrick base pairs is not dramatically modified by the presence
of a third strand.
The triplex with the TfF mutation in the Hoogsteen position
(5) shows reversible breathing movements, where the F is
displaced to the major groove, leading to the partial or total
(24) Ragunnathan, G.; Miles, H. T.; Sasisekharan, V. Biochemistry 1993,
32, 455.
(25) Howard, F. N.; Miles, H. T.; Liu, K.; Frazier, J.; Raghunathan, G.;
Sasisekharan, V. Biochemistry 1993, 31, 10671.
Figure 4. Molecular interaction potential (MIP) and solvation maps obtained for the duplexes containing T or F. The MIP map is contoured at
-5.0 kcal/mol, and the solvation map is contoured at a water density of around 2.5 g/cm3.
Cubero et al.
Figure 5. Ab initio HF/6-31G(d) molecular electrostatic potential map computed at 3.4 above the aromatic plane. Colored areas correspond to
MEP values from -10 (deep red) to +10 (deep blue) kcal/mol.
d(AT-T)
d(AF-T)
d(AT-F)
0.97(0.16)
1.11(0.21)
1.03(0.18)
1.05(0.18)
1.10(0.21)
1.07(0.17)
1.08(0.22)
1.11(0.19)
1.01(0.18)
triplex
d(AT-T)
triplex
d(AF-T)
triplex
d(AT-F)
X-displacement ()
inclination (deg)
rise ()
roll (deg)
twist (deg)
propeller twist (deg)
phase angle (deg)
puckering amplitude (deg)
MM groove width ()
mM groove width ()
m groove width ()
-3.0(0.5)
-3.2(4.7)
3.4(0.1)
3.2(1.8)
29.6(0.7)
-12.2(4.7)
111(29)
41.5(5,8)
15.6(1.2)
9.0(0.3)
11.5(0.5)
-3.0(0.6)
-1.6(5.0)
3.4(0.1)
3.6(1.8)
29.5(0.6)
-12.8(4.7)
110(38)
41.5(5.8)
15.1(1.0)
9.0(0.3)
12.1(0.5)
-3.1(0.5)
-1.4(4.3)
3.4(0.1)
3.1(1.8)
29.7(0.7)
-13.5(4.8)
113(26)
41.7(5.7)
15.3(0.8)
9.0(0.3)
11.5(0.5)
a Averages are done using the last 1.5 (reference triplex) and 2.0
(triplexes with mutation) ns of the trajectories. MM means major-major
groove, mM means minor-major groove, and m refers to the minor
groove. See ref 20 for nomenclature.
Figures 2 and 8). This value agrees well with SCRF estimates
(-8.6 kcal/mol) obtained at the AM1/MST level,27 which gives
us further confidence in the quality of the MD/TI calculations.
Simple algebra with free energy differences in Figure 2 gives
us all the thermodynamic data associated with the effect of TfF
mutations on duplex and triplex stability. Results in Table 5
show that the TfF mutation destabilizes the duplex DNA by
about 5 kcal/mol, in good agreement with experimental values
found by Kool and co-workers (values around 4 kcal/mol from
refs 6, 7, and 9). Analysis of the central d(GpApA) sequence
of the MD-averaged structures (see Table 6) shows that the
largest difference between T and F in terms of duplex stability
lies in the H-bonding term, which is near 10 kcal/mol more
stable for duplexes containing the d(AT) pair.28
In summary, MD and MD/TI simulations strongly suggest
that the presence of a single d(AF) step destabilizes the duplex,
but the resulting DNA is still stable at room temperature, in
good agreement with all known experimental evidence.6,7,9 In
practice, the largest differences between a DNA duplex containing a d(AF) pair and the natural one should be related to (i)
the change in the intrinsic reactivity in the minor groove, which
might modify the ability of DNA to interact with small
molecules or some DNA minor-groove binding proteins, and
(ii) the existence of breathing events which might also affect
the reactivity of DNA, and contribute to the activation of DNArepairing systems.
The formation of most parallel-motif triplexes is a two-step
process. First a DNA duplex is formed, and second the triplexforming oligonucleotide is bound29 to give the triplex. The
transition from a duplex with a d(AT) pair to a triplex with a
d(AT-T) trio is around 3.2 kcal/mol more favorable than the
transition from the same duplex to a triplex with a d(AT-F)
trio (see Table 5). Interestingly, if the duplex contains a d(A
F) pair, the formation of a d(AF-T)-containing triplex is not
disfavored with respect to the formation of a triplex with a d(A
T-T) trio from a duplex with the d(AT) pair (differential free
(27) (a) Luque, F. J.; Negre, M.; Orozco, M. J. Phys. Chem. 1993, 97,
4386. (b) Luque, F. J.; Bachs, M.; Orozco, M. J. Comput. Chem. 1994, 15,
847. (c) Orozco, M.; Bachs, M.; Luque, F. J. J. Comput. Chem. 1995, 16,
564-575. (d) Miertus, S.; Scrocco, E.; Tomasi, J. Chem. Phys. 1981, 55,
117. (e) Miertus, J.; Tomasi, J. Chem. Phys. 1982, 65, 239.
(28) AM1/MST calculations show that hydration disfavors AT and A
F pairings by around 8.1 and 3.6 kcal/mol. This implies that solvation favors
relative binding of F to A, with respect to the binding of T to A, by around
4.5 kcal/mol, thus reducing the apparent large difference in stability
suggested from H-bonding considerations alone.
(29) Marky, L. A.; Breslauer, K. J. Biopolymers 1987, 26, 1601.
Figure 7. Molecular interaction potential (MIP) and solvation maps obtained for the triplexes containing T or F. The MIP maps are contoured at
-6.0 kcal/mol, and the solvation maps are contoured at a water density of around 3.5 g/cm3
energy -0.6 kcal/mol). This finding can be explained considering that the characteristics of the major groove of the duplex
containing the d(AF) pair are very similar to those of the
reference duplex (see above and Table 2). The similarity of
triplex-forming characteristics of the parent and TfF mutated
duplex strongly suggests that a duplex DNA containing a d(A
F) pair would have, a priori, a similar ability to interact with
most major groove-binding DNA proteins as the parent
duplex.
Perhaps more interesting is the study of the effect of the TfF
mutation on the thermodynamics of the whole assembly process,
from single-stranded species to triplex. As seen in Figure 8 and
Table 5, the presence of one F destabilizes the resultant triplex
with respect to the reference triplex (d(ATT)) by around 3.2
kcal/mol when F is in the Hoogsteen position, and 4.7 kcal/
Cubero et al.
Figure 8. Free energy profiles obtained from different trajectories. Average values (in kcal/mol) are shown with their standard errors. See also
Figure 2.
Table 5. Changes in the Free Energy of the DuplexfTriplex,
SinglefDuplex, and SinglefTriplex Process Associated with the
TfF Mutationa
G
SE
(kcal/mol) (kcal/mol)
folding
process
sequence
duplexf triplex
duplexf triplex
duplexf triplex
singlef duplex
singlef duplex
singlef triplex
singlef triplex
singlef triplex
d(T)+d(AT)f d(AT-T)
d(T)+d(AF)f d(AF-T)
d(F)+d(AT)f d(AT-F)
d(T)+d(A)fd(AT)
d(F)+d(A)fd(AF)
d(T)+d(A)+d(T)f d(AT-T)
d(F)+d(A)+d(T)fd(AF-T)
d(T)+d(A)+d(F)fd(AT-F)
0.0
-0.6
3.2
0.0
5.3
0.0
4.7
3.2
0.0
0.3
0.2
0.0
0.1
0.0
0.3
0.2
Standard errors for the different estimates are displayed. See Figure
8 and text for details.
Table 6. Interaction Terms (in kcal/mol) for the Central Three
Steeps (d(GpApA)) of the Duplexes Containing the d(AT) and
d(AF) Paira
interaction
d(AT)
d(AF)
stacking (intrastrand)
stacking (interstrand)
stacking (total)
hydrogen bonding
total interaction
-27.6
-4.6
-32.3
-54.5
-86.8
-24.8
-6.2
-31.0
-44.2
-75.2
MD-TI calculations strongly suggest that the TfF replacement in Watson-Crick or Hoogsteen positions destabilizes the
triplex DNAs by 3-5 kcal/mol with respect to the parent triplex
3. With use of the empirical approach of Roberts and Crothers,31
(30) B3LYP/6-31G(d, p) geometry optimization suggests stabilization
energies (after BSSE correction) of -3.1, and -2.9 kcal/mol for Hoogsteen
and Watson-Crick AF pairs. These values are remarkably similar to those
obtained using AMBER parameters: -3.3 (A-F, H) and -3.2 (A-F, WC)
kcal/mol.
(31) Roberts, R. W.; Crothers, D. M. Proc. Natl. Acad. Sci. U.S.A. 1996,
93, 4320.
d(AT-T)
d(AF-T)
d(AT-F)
stacking (intrastrand)
stacking (interstrand)
stacking (total)
hydrogen bonding
total interaction
-40.1
-20.6
-60.7
-115.0
-175.7
-35.9
-20.1
-56.0
-96.5
-152.5
-38.7
-17.5
-56.2
-104.9
-161.1