BI BII Substate Transitions Induce Changes in the Hydration of B-DNA

The BII substate distribution within one strand is shown to give an alternating pattern ... Furthermore, BI h BII conformational transitions are shown...
0 downloads 0 Views 306KB Size
J. Phys. Chem. B 2001, 105, 10379-10387

10379

BI h BII Substate Transitions Induce Changes in the Hydration of B-DNA, Potentially Mediating Signal Transduction from the Minor to Major Groove Wolfgang Flader, Bernd Wellenzohn, Rudolf H. Winger, Andreas Hallbrucker, Erwin Mayer, and Klaus R. Liedl* Institute of General, Inorganic and Theoretical Chemistry, UniVersity of Innsbruck, Innrain 52a, A-6020 Innsbruck, Austria ReceiVed: October 31, 2000; In Final Form: July 26, 2001

The EcoRI DNA dodecamer d(CGCGAATTCGCG)2 was investigated by a 10 ns molecular dynamics simulation. The BII substate distribution within one strand is shown to give an alternating pattern and thus avoids two successive phosphates remaining simultaneously in the BII substate. Increased BII population (BI/ BII ratio of 3.7 compared to 9.2 in a previous simulation) was observed by employing improved force field parameters derived recently by Cheatham et al. [Cheatham, T. E., III; Cieplak, P.; Kollman, P. A. J. Biomol. Struct. Dyn. 1999, 16, 845-862]. Furthermore, BI h BII conformational transitions are shown to correlate with the hydration pattern of the phosphate group, the sugar oxygen, and the minor and the major groove. Therefore, BI/BII substates can act as an additional conversation tool to transfer binding information from a drug ligand in the minor groove to a protein binding site in the major groove.

1. Introduction Hydration and dynamical behavior of DNA play a decisive role for its biological function.1-12 Especially the interconversion between the conformational substates of B-DNA, BI and BII, may participate significantly in ligand binding within the backbone phosphate groups, which has been shown to be important in numerous investigations.2,9,10,12-23 In the BI state the corresponding  and ζ angles derived from X-ray structures are between 120-210° (trans) and 235-295° (gauche-), respectively; for BII  lies between 210° and 300° (gauche-) and ζ between 150° and 210° (trans).24 So a BI h BII interconversion represents a small but significant reorientation of the charged phosphate group of the DNA backbone. Figure 1 shows a part of the DNA dodecamer with one phosphate group in both the BI and BII conformation. The upper picture is a snapshot at 900 ps with the highlighted phosphate in the BII substate; the lower picture is a snapshot at 2100 ps, therein the same phosphate is in the BI state. The EcoRI DNA dodecamer d(CGCGAATTCGCG)2 was subject of numerous experimental and theoretical studies.25-49 The vast amount of data available for this dodecamer and its complexes with different ligands provide an ideal field for comparative investigations regarding structure and dynamics of free and complexed B-DNA. A detailed knowledge of molecular recognition processes between DNA and proteins or small ligands can give insight into the regulation of the expression of genes23,50-67 and therefore explains function and misfunction of regulatory processes in the cell. This knowledge is of crucial importance for the development of strategies and highly selective drugs for the treatment of cancer and viral or microbial diseases.31,35,38,51,54,58,60,68-82 2. Methods Molecular dynamics (MD) simulations provide an excellent tool for the investigation of structure and dynamics of biomolecules in their biologically active form.83-86 Improved param-

etrization of the force fields and description of the long-range interactions using the Ewald summation (implemented in the AMBER program87 in the form of the particle mesh Ewald method88,89) enable MD simulations to provide stable B-DNA trajectories39,90,91 and to supplement experimental data by information of unique space and time resolution. The fast interconversion between the BI/BII substates of DNA, which are assumed to interconvert on the subnanosecond time scale,92-98 allow an exhaustive description of this conformational aspect by computer simulations. This gives an excellent completion to the current experimental methods like X-ray crystallography, and NMR and IR spectroscopy, which provide valuable experimental data as starting information. These experimental methods suffer from crystallization effects or insufficient resolution in time and space, which make a detailed description of the BI h BII substate interconversion behavior very difficult. Nevertheless, they are used to confirm the accuracy of the computer simulations. 3. Computational Details As a starting structure the canonical B-form99,100 of the DNA oligonucleotide d(CGCGAATTCGCG)2 was used. To achieve electroneutrality of the system containing 11 PO4- anions in each strand, 22 Na+ counterions were added using the CION program of the AMBER package. TIP3P101 Monte Carlo water boxes were added to form a solvation box at least 12 Å thick in all directions, resulting in a system of 68 × 45 × 45 Å3 containing 3874 water molecules. The corresponding Γ value (i.e., water molecules per nucleotide) is about 161. The simulation was carried out using the AMBER5102 suite of programs employing the all atom force field103 and the modified parameters from Cheatham et al.104 The simulation protocol was the same as in the work of Winger et al.105 4. Results In the course of the simulation the number of phosphate groups being simultaneously in the BII substate ranged from 1

10.1021/jp004046q CCC: $20.00 © 2001 American Chemical Society Published on Web 10/03/2001

10380 J. Phys. Chem. B, Vol. 105, No. 42, 2001

Flader et al.

Figure 2. Conformational substates as a function of time: The time the respective -angle of the 22 phosphates is in substate BII is marked by a black line/dot. X-axis: Enumeration of the interconverted phosphate groups according to Kopka et al.120 Y-axis: Time (ps).

Figure 1. Part of the DNA dodecamer with one phosphate group in both the BI and BII conformations. In the upper picture (snapshot at 900 ps) the highlighted phosphate is in the BII substate; in the lower picture (snapshot at 2100 ps) the same phosphate is in BI.

to 9. The distribution of the BII substates over the whole dodecamer shows a clear preference of the first and last third of every strand (see Figure 2). The ATTC regions contain very few BII states, especially both of the TpT phosphates stay in the BI state over the whole 10 ns. The time-averaged number of phosphate groups in the BII state of 4.7 corresponds to a BI/BII ratio of 3.7 (about 21% BII). This increased BII population (previous simulations105 showed a BI/BII ratio of 9.2) results from the use of improved force field parameters derived recently from Cheatham et al.104 and gives a completely new picture of the BI h BII interconversion behavior. Although the increased BII population leads to regions with a relatively high concentration of this substate, the distribution of BII states within one strand shows an alternating pattern avoiding two successive phosphates remaining simultaneously in the BII substate. For phosphates preferentially in the BII substate (P5, P11, P17, and P23) this anticorrelation with the substate of adjacent phosphates (P4, P6, P10, P16, P18, and P22) results in a mean probability of 0.5% to find a BII-BII sequence. (Phosphates of endstanding base pairs were not considered in this summary because their artificial flexibility may influence the results.)

The complete set of all -angles and their conformation as a function of time is depicted in Figure 2 for the 22 phosphates of the dodecamer’s 22 base steps. As observed in previous investigations,105 analysis of counterions did not show direct correlations between BI h BII substate transitions and dynamics of the sodium ions. For selected phosphates that show a representative BI h BII dynamic (P4, P5, P6, P10, and P11), the occupation with counterions was rather low. The area within 5 Å arround their phosphate oxygens, O1 and O2, was occupied by a sodium ion only 0.8% of the time. Within 5 Å arround the minor groove atoms N3 of G4, A5, A6, G10, and O2 of C11, the mean occupation with counterions was less than 0.7%. Within an 8 Å area of these phosphates and minor groove atoms the mean presence of a sodium ion was 3.9 and 3.3%, respectively. Anyway, also with this increased occupation no correlations with BI h BII substate transitions were observed. This is not in contradiction with the findings of Wilson et al. that minor groove structure and dynamics show direct ion dependence.106 They observed a correlation of minor groove width (which is surely related to BI/BII substates) with ion positions especially in the AATT region of the EcoRI DNA dodecamer. In this part the BII activity is rather low. Therefore, this sequence may be unsuitable for a detailed analysis of possible correlations between BI h BII substate transitions and the mentioned counterion-influenced minor groove width. The correlation between the BI/BII interconversions and sugar puckering and destacking of adjacent base pairs is shown in Figure 3 for the G4-A5 (left column) and G22-C23 base steps (right column). These two representative base steps were chosen because of the high BII amount with very distinct BI h BII interconversions. The -angles are shown as a function of time in pictures a and b, in pictures c and d the destacking of adjacent base pairs is represented by the projected C1′ (sugar carbon attached to the base) distances107 as criterion (Dxy), the sugar puckering of the preceding base is shown in pictures e/i and f/j, respectively, and the one of the following base is in pictures g/k and h/l, respectively. As expected, the BI h BII interconversions show a strong correlation with both the sugar pucker and the destacking of both adjacent base pairs. Thereby, the correlations are stronger for the sugar previous to a phosphate than for the following one. Correlation factors between the BI/ BII substates of the phosphate P5 and sugar puckering of G4 are +0.56 and -0.24 for the amplitude and the angle, respectively. The sugar puckering of A5 shows correlation

BI h BII Substate Transitions in B-DNA

J. Phys. Chem. B, Vol. 105, No. 42, 2001 10381

Figure 3. Substate interconversion of the G4-A5 (left column) and the G22-C23 base step (right column). (a) and (b): -angles as a function of time. (c) and (d): destacking (Dxy) of adjacent residues by using the projected C1′ (sugar carbon attached to the base) distances107 as criterion. Sugar puckering angle of the pseudorotation cycle of the furanose ring in G4 (e), G22 (f), A5 (g), and C23 (h). Pucker amplitude of the furanose ring in G4 (i), G22 (j), A5 (k), and C23 (l).

factors of -0.28 and +0.40 to the substates of the same phosphate P5. In the case of the phosphate P23 the difference is even stronger. The puckering of the previous sugar (G22) correlates with factors of 0.68 for the amplitude and -0.59 for the angle, the respective values for the sugar of C23 are -0.20 and +0.32 (correlation factors were calculated using the corrcoef function implemented in the matlab108 program).

The puckering angle of the pseudorotation cycle of the furanose ring fluctuates around a mean value of about 180° in the BI state and around a mean value of 150° in the BII state. Very seldom and without assignable correlation it changes to a value of about 350° for a few tens of picoseconds by decreasing from the 180° region to about -10° ()350°). The pucker amplitude of the furanose ring fluctuates around a mean value

10382 J. Phys. Chem. B, Vol. 105, No. 42, 2001

Flader et al.

Figure 4. Substate interconversion of the G4-A5 (left column) and the G22-C23 base step (right column). (a) and (b): -angles as a function of time. A time dependent hydration pattern is shown by the number of water oxygens closer than the specified distance to the following atoms: within 3.3 Å of N2 of G4 (c) and G22 (d), within 3.1 Å of N3 of G4 (e) and N3 of G22 (f), within 3.3 Å of the phosphate oxygen O1 of P5 (g) and P23 (h), within 3.0 Å of O3′ of G4 (i) and G22 (j), within 3.2 Å of sugar ogygen O4′ of A5 (k) and C23 (l). The different distances were chosen in order to give an optimal representation of the correlations.

of about 0.33 Å in the BI state, and around a mean value of 0.45 Å in the BII state. The destacking reference C1′-C1′ distance between the two adjacent residues, Dxy, fluctuates around a mean value of about 4.8 Å in the BII state, and a broad area with mean values from 2 to 3.5 Å in the BI state. The changes in the hydration of adjacent base pairs in correlation with the BI h BII interconversions are shown in Figure 4 (minor groove, ester oxygen O3′, phosphate oxygen

O1, and sugar oxygen O4′) and Figure 5 (major groove and phosphate oxygen O2). The hydration of some edges of the bases in the minor groove shows a moderate correlation to the BI h BII interconversions: N2 (G4/G22) and N3 (G4/G22) are slightly more hydrated in the BI than in the BII state (pictures c/d and e/f in Figure 4). The hydration of N3 (A5) and O2 (C23) (not shown in Figure 4) is not correlated to the substate interconversions (correlation

BI h BII Substate Transitions in B-DNA

J. Phys. Chem. B, Vol. 105, No. 42, 2001 10383

Figure 5. Substate interconversion of the G4-A5 (left column) and the G22-C23 base step (right column). (a) and (b): -angles as a function of time. A time dependent hydration pattern is shown by the number of water oxygens closer than the specified distance to the following atoms: within 3.5 Å of O6 of G4 (c) and G22 (d), within 3.5 Å of N7 of G4 (e) and G22 (f), within 3.5 Å of N6 of A5 (g) and N4 of C23 (h), within 3.5 Å of N7 of A5 (i), within 3.1 Å of the phosphate oxygen O2 of P5 (j) and P23 (k).

factor +0.03 and -0.07). The phosphate oxygens O1 and O3′, and the sugar oxygen O4′ show a strong correlation to the BI h BII interconversions: the charged phosphate oxygen pointing toward the minor groove (O1) and the phosphate ester oxygen O3′ are most hydrated in the BI state (pictures g/h and i/j in Figure 4), whereas the sugar oxygen O4′ in this state is less hydrated (pictures k/l in Figure 4). In the major groove the correlations between the hydration of the edges of the bases and the BI/BII substates are as

follows: O6 (G4/G22) and N7 (G4/G22) are more hydrated in the BII substate (pictures c/d and e/f in Figure 5), N6 and N7 (A5) are more hydrated in the BI substate (pictures g and i in Figure 5). The hydration of N4 (C23) (picture h in Figure 5) shows only a poor correlation to the substate interconversions, it is slightly better in the BI state. The base edge of cytosine has only one hydrophilic atom (N4) pointing toward the major groove (compared to N6 and N7 of adenine), so there is one picture missing in the right column of Figure 5. There is a

10384 J. Phys. Chem. B, Vol. 105, No. 42, 2001

Flader et al.

TABLE 1: Correlation between BI h BII Substate Transitions and the Hydration Shown in Figures 4 and 5 Represented by the Correlation Factor (Corr), Calculated with a Smoothed Hydration Description Using the Mean Value of 10 ps for Every Stepa minor groove

G4-A5 atom corr

N2 (G4) -0.33

N3 (G4) -0.14

a

O1 (P5) -0.75

O3′(G4) -0.77

O4′(A5) +0.57

O6 (G4) +0.37

N7 (G4) +0.25

minor groove

G22-C23 atom corr

N3 (A5) +0.03

major groove

N2 (G22) -0.62

N3 (G22) -0.18

O2 (C23) -0.07

O1 (P23) -0.77

N6 (A5) -0.48

N7 (A5) -0.57

O2 (P5) +0.40

major groove O3′(G22) -0.82

O4′(C23) +0.48

O6 (G22) +0.58

N7 (G22) +0.63

N4 (C23) -0.18

NA

O2 (P23) +0.45

A positive sign means that the atom is more hydrated in the BII substate; with a negative sign the hydration is better in the BI substate.

moderate correlation with a better hydration in the BII substate for the charged phosphate oxygen pointing toward the major groove (O2) (pictures j/k in Figure 5), although the matches of the changes in the hydration with the BI h BII substate transitions are not as evident as for the phosphate oxygen O1. The correlations of all these changes in the hydration with the BI h BII substate transitions are summarized in Table 1 by a correlation factor calculated between the time dependent -angle and a smoothed hydration description represented by the mean value of 10 ps of the water population within a given distance for every step (correlation factors were calculated using the corrcoef function implemented in the matlab108 program). To give a better insight into the chronology of the observed effects in correlation with the BI h BII substate transitions, Figure 6 shows an enlarged part of the previous pictures for the G4-A5 step. It represents the time period from 1500 to 2500 ps, containing a transition from BII to BI at about 1720 ps and a transition from BI to BII at about 2100 ps (a few trials to change from the BII to the BI substate at about 1580, 2190, 2330, and 2460 ps are too short and do not reach an equilibrated BI substate). The pictures are arranged in the same manner as in Figures 3-5. So the -angle is shown as a function of time in the first row, the left column represents destacking and sugar puckering, and the middle and the right column show the hydration of the minor and the major groove, respectively. For better visibility the hydration herein is represented by the mean value of 10 ps for every step (as used for the calculation of the correlation factor), instead of the single values shown in Figures 4 and 5. The short changes from BII to the BI substate that turn back after a few picoseconds do not have a sufficient effect on the observed parameters to allow a significant conclusion. The transition from BII to BI at 1720 ps initiates the loss of destacking, whereas the changes of sugar puckering began a few tens of picoseconds earlier and are just complete at this time. The transition from BI to BII at 2100 ps is preceded by a slower convergence to the typical values for the BII substate of both destacking and sugar pucker (left column of Figure 6). Changes in the hydration of the minor and major grooves are difficult to correlate exactly to the substate transitions at this resolution (in a detailed view there are strong local fluctuations and in a more global view the graphs are relatively smooth, which makes an unambiguous assignment impossible). For the oxygens of the DNA backbone (phosphate group and sugar) the correlation is more evident and can be summarized in the conclusion that the change in the hydration starts about contemporary with the substate transition and in some cases needs a few tens of picoseconds to reach the new equilibrated state (middle and right column of Figure 6). 5. Discussion The use of improved force field parameters derived recently from Cheatham et al. leads to an increased BII population104

(BI/BII ratio of 3.7 compared to 9.2 in a previous simulation). Despite this relatively high BII population concentrated in the first and the last third of the dodecamer, the BII states within one strand are alternating. This avoids two successive phosphates remaining simultaneously in the BII state, and it confirms the observation of Grzeskowiak et al. that “BII conformations are seldom found at adjacent steps along a chain”.107 The hierarchy of BII forming tendencies predicted in this same work107 was later extended by Winger et al.,105 leading to the conclusion that the BII substate is highly unfavorable only for pyrimidinepyrimidine steps (which are in the BII substate in only 10 of 353 analyzed steps of this type). This is consistent with our simulation where the ATTC regions show only a few transitions to BII for very short periods (see Figure 2). Another observation made on analysis of single-crystal structures of DNA decamers at high resolution107,109 and confirmed later by molecular dynamics simulations of the dodecamer in aqueous solution105 regards the correlation between the BI/BII substates and the destacking of the bases on either side of the phosphate. Our simulation confirms again the rule that destacking is a necessary but not a sufficient condition for an adoption of the BII phosphate conformation. A detailed observation of the destacking (Figure 6: second picture of the left column) shows the fast decrease of this parameter after the transition from BII to BI, whereas the following increase is slower and the typical value for the BII substate is reached before the change from BI to BII occurs. The hydration of the whole B-DNA is influenced by the BI/ BII substate pattern, and the observed changes in the hydration shell are in agreement with spectroscopic data showing BI h BII substate interconversions by IR experiments.44,46 ,110 A detailed analysis of the accessible surface for the hydrogenbonding atoms of the major groove was presented recently by Tisne´ et al.111 They describe an increased accessible surface for most of the O6 and N7 atoms of the guanines in the BII substate, which agrees with increasing hydration of these atoms in the BII substate observed in our simulation (see Table 1). But, their results show also that sometimes the correlation between BI h BII substate interconversions and the accessible surface for a certain atom type can be ambiguous even when the sequence is the same over four base pairs. So one should be careful to derive simple predictive rules from a few data. Instead, these dynamics of the BI/BII hydration may be seen in general as a transfer of this substate information to protein or drug binding sites in the major or minor groove and therefore may give an additional meaning to the sequence dependent BI/ BII behavior.24,112-114 Besides this, such specific movements in the first hydration shell may be able to reduce the thermodynamic penalty for the dehydration of hydrophilic groups and hence support the formation of hydrogen bond contacts in DNA-ligand complexes.

BI h BII Substate Transitions in B-DNA

J. Phys. Chem. B, Vol. 105, No. 42, 2001 10385

Figure 6. Substate interconversion of the G4-A5 base step in the period between 1500 and 2500 ps. The left column shows details of the pictures in Figure 3 (-angle, destacking and sugar puckering), the middle and the right columns present details corresponding to the left columns of Figures 4 and 5, respectively (time dependent hydration represented herein by the mean value of 10 ps for every step).

While the correlation of the hydration in the major groove with the BI h BII substate interconversions is surprisingly high, in the minor groove only a few of the edges of the bases show a hydration pattern with significant dependence on these substate transitions. Nevertheless, the evident correlation of the hydration

of the sugar oxygen O4′ and the phosphate ester oxygen O3′ (both pointing toward the minor groove) with the BI/BII interconversions allows us to widen the influence of the substate pattern also to the ligand binding sites in the minor groove. Especially changes in the hydration of the phosphate groups

10386 J. Phys. Chem. B, Vol. 105, No. 42, 2001 show a strong correlation with BI h BII substate transitions. These differences in the shielding of the charges of the phosphates may lead to a dynamical behavior of the electrostatic field in the environment of B-DNA, which is responsible for the orientation of ligands in the first steps of complex formation.115 Certainly the phosphate groups of the DNA backbone contribute the major part of DNA-protein contacts not only in the so-called nonspecific DNA-ligand complexes but also in complexes with proteins of high specificity,20,116 and therefore the sequence dependent BI/BII pattern and dynamics may represent an additional tool to increase their sequence specificity. Finally, this dynamical behavior of the B-DNA backbone and the correlated hydration may play a crucial role in the mechanism by which minor groove binding drugs act on the binding site of proteins in the major groove. In a complex with a small drug molecule in the minor groove, the BI/BII substate pattern of B-DNA can be frozen-in in a form completely different from that of the free DNA in solution.117-119 Usually this is a direct consequence of ligand contacts with the DNA backbone phosphates, but partially it may also result from the substitution of water molecules by the ligand and the involved loss of flexibility in the minor groove. This fixed unusual BI/BII pattern leads to the discussed changes in the hydration, and so the electrostatic behavior and the hydration of potential binding sites in the major groove may be influenced also by a minor groove binding drug in relatively rigid DNA tracts where larger structural deviations like bending or changes in the groove width are negligible. 6. Summary and Conclusions A 10 ns state-of-the-art molecular dynamics simulation of the EcoRI DNA dodecamer was performed and analyzed with special attention on BI/BII substates and hydration. The DNA backbone exhibited a highly dynamical behavior with respect to BI h BII interconversions, and only the two TpT steps remained in the BI substate over the whole 10 ns. Despite the high BII population (about 21%) the substate distribution showed an alternating pattern. This avoids BII conformations being found simultaneously at adjacent steps along a chain. The hydration of the whole B-DNA performed a dynamic that is correlated to the BI h BII substate transitions. This proposes a possible additional function of the backbone substates for sequence specific recognition and for signal transductions from minor to major groove. Acknowledgment. This work was supported by a grant of the Austrian Science Fund (grant number P13845-TPH). References and Notes (1) Berman, H. M. Curr. Opin. Struct. Biol. 1994, 4, 345-350. (2) Beamer, L. J.; Pabo, C. O. J. Mol. Biol. 1992, 227, 177-196. (3) Hartmann, B.; Lavery, R. Q. ReV. Biophys. 1996, 29, 309-368. (4) Prive´, G. G.; Heinemann, U.; Chandrasegaran, S.; Kan, L.-S.; Kopka, M. L.; Dickerson, R. E. Science 1987, 238, 498-504. (5) Robinson, C. R.; Sligar, S. G. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 2186-2191. (6) Egli, M.; Tereshko, V.; Teplova, M.; Minasov, G.; Joachimiak, A.; Sanishvili, R.; Weeks, C. M.; Miller, R.; Maier, M. A.; An, H.; Cook, P. D.; Manoharan, M. Biopolymers 2000, 48, 234-252. (7) Kunkel, T. A.; Wilson, S. H. Nature Struct. Biol. 1998, 5, 95-99. (8) Dickerson, R. E. Nucl. Acids Res. 1998, 26, 1906-1926. (9) Somers, W. S.; Phillips, S. E. Nature 1992, 359, 387-393. (10) Carrondo, M. A.; Coll, M.; Aymami, J.; Wang, A. H.; van der Marel, G. A.; van Boom, J. H.; Rich, A. Biochemistry 1989, 28, 78497859. (11) Gupta, G.; Bansal, M.; Sasiskharan, V. Proc. Natl. Acad. Sci. U.S.A. 1980, 77, 6486-6490.

Flader et al. (12) Luger, K.; Ma¨der, A. W.; Richmond, R. K.; Sargent, D. F.; Richmond, T. J. Nature 1997, 389, 251-260. (13) Suck, D.; Lahm, A.; Oefner, C. Nature 1988, 332, 464-468. (14) Otwinowski, Z.; Schevitz, R. W.; Zhang, R. G.; Lawson, C. L.; Joachimiak, A.; Marmorstein, R. Q.; Luisi, B. F.; Sigler, P. B. Nature 1988, 335, 321-329. (15) Wolberger, C.; Dong, Y. C.; Ptashne, M.; Harrison, S. C. Nature 1988, 335, 789-795. (16) Jordan, S. R.; Pabo, C. O. Science 1988, 242, 893-899. (17) Luisi, B. F.; Xu, W. X.; Otwinowski, Z.; Freedman, L. P.; Yamamoto, K. R.; Sigler, P. B. Nature 1991, 352, 497-505. (18) Beese, L. S.; Derbyshire, V.; Steitz, T. A. Science 1993, 260, 352355. (19) Shakked, Z.; Guzikevich-Guerstein, G.; Frolow, F.; Rabinovich, D.; Joachimiak, A.; Sigler, P. B. Nature 1994, 368, 469-473. (20) Mandel-Gutfreund, Y.; Schueler, O.; Margalit, H. J. Mol. Biol. 1995, 253, 370-382. (21) Rhodes, D. Nature 1997, 389, 231-233. (22) Agback, P.; Baumann, H.; Knapp, S.; Ladenstein, R.; Haerd, T. Nature Struct. Biol. 1998, 5, 579-584. (23) Chen, Y.-Q.; Gosh, S.; Gosh, G. Nature Struct. Biol. 1998, 5, 6773. (24) Schneider, B.; Neidle, S.; Berman, H. M. Biopolymers 1997, 42, 113-124. (25) Wing, R.; Drew, H.; Takano, T.; Broka, C.; Tanaka, S.; Itakura, K.; Dickerson, R. E. Nature 1980, 287, 755-758. (26) Drew, H. R.; Wing, R. M.; Takanao, T.; Broka, C.; Tanaka, S.; Itakura, K.; Dickerson, R. E. Proc. Natl. Acad. Sci. U.S.A. 1981, 78, 21792183. (27) Drew, H. R.; Dickerson, R. E. J. Mol. Biol. 1981, 151, 535-556. (28) Dickerson, R. E.; Drew, H. R. J. Mol. Biol. 1981, 149, 761-786. (29) Drew, H. R.; Samson, S.; Dickerson, R. E. Proc. Natl. Acad. Sci. 1982, 79, 4040-4044. (30) Fratini, A. V.; Kopka, M. L.; Drew, H. R.; Dickerson, R. E. J. Biol. Chem. 1982, 257, 14686-14707. (31) Patel, D. J. Proc. Natl. Acad. Sci. U.S.A. 1982, 21, 6424-6428. (32) Ravishanker, G.; Swaminathan, S.; Beveridge, D. L.; Lavery, R.; Sklenar, H. J. Biomol. Struct. Dyn. 1989, 6, 669-699. (33) Lane, A. N.; Jenkins, T. C.; Brown, T.; Neidle, S. Biochemistry 1991, 30, 1372-1385. (34) Whitka, J. M.; Swaminathan, S.; Srinivasan, J.; Beveridge, D. L.; Bolton, P. H. Science 1992, 255, 597-599. (35) Perree-Fauvet, M.; Gresh, N. J. Biomol. Struct. Dyn. 1994, 11, 1203-1224. (36) Goodsell, D. S.; Kopka, M. L.; Dickerson, R. E. Biochemistry 1995, 34, 4983-4992. (37) Trent, J. O.; Clark, G. R.; Kumar, A.; Wilson, W. D.; Boykin, D. W.; Hall, J. E.; Tidwell, R. R.; Blagburn, B. L.; Neidle, S. J. Med. Chem. 1996, 39, 4554-4562. (38) Monaco, R. R.; Polkosnik, W.; Dwarakanath, S. J. Biomol. Struct. Dyn. 1997, 15, 63-67. (39) Young, M. A.; Ravishanker, G.; Beveridge, D. L. Biophys. J. 1997, 73, 2313-2336. (40) Srinivasan, J.; Cheatham, T. E., III.; Cieplak, P.; Kollman, P. A.; Case, D. A. J. Am. Chem. Soc. 1998, 120, 9401-9409. (41) Shui, X.; McFail-Isom, L.; Hu, G. G.; Williams, L. D. Biochemistry 1998, 37, 8341-8355. (42) Tereshko, V.; Minasov, G.; Egli, M. J. Am. Chem. Soc. 1999, 121, 470-471. (43) Tereshko, V.; Minasov, G.; Egli, M. J. Am. Chem. Soc. 1999, 121, 3590-3595. (44) Pichler, A.; Ru¨disser, S.; Mitterbo¨ck, M.; Huber, C. G.; Winger, R. H.; Liedl, K. R.; Hallbrucker, A.; Mayer, E. Biophys. J. 1999, 77, 398409. (45) Pichler, A.; Ru¨disser, S.; Winger, R. H.; Liedl, K. R.; Hallbrucker, A.; Mayer, E. J. Am. Chem. Soc. 2000, 122, 716-717. (46) Pichler, A.; Ru¨disser, S.; Winger, R. H.; Liedl, K. R.; Hallbrucker, A.; Mayer, E. Chem. Phys. 2000, 258, 391-404. (47) Woods, K. K.; McFail-Isom, L.; Sines, C. C.; Howerton, S. B.; Stephens, R. K.; Williams, L. D. J. Am. Chem. Soc. 2000, 122, 15461547. (48) Halle, B.; Denisov, V. P. Biopolymers 2000, 48, 210-233. (49) Clark, G. R.; Squire, C. J.; Baker, L. J.; Martin, R. F.; White, J. Nucl. Acids Res. 2000, 28, 1259-1265. (50) Pabo, C. O.; Sauer, R. T. Annu. ReV. Biochem. 1992, 61, 10531095. (51) Gottesfeld, J. M.; Neely, L.; Trauger, J. W.; Baird, E. E.; Dervan, P. B. Nature 1997, 387, 202-205. (52) Beerli, R. R.; Dreier, B.; III, C. F. B. Proc. Natl. Acad. Sci. U.S.A. 2000, 97 (4), 1495-1500. (53) Dickinson, L. A.; Trauger, J. W.; Baird, E. E.; an, P. B. D.; Graves, B. J.; Gottesfeld, J. M. J. Biol. Chem. 1999, 30 (18), 12765-12773. (54) Choo, Y.; Sanchez, I.; Klug, A. Nature 1994, 372, 642-645.

BI h BII Substate Transitions in B-DNA (55) Gehring, W. J.; Qian, Y. Q.; Furukubo-Tokunaga, M. K.; Schlier, F.; Resendez-Perez, D.; Affolter, M.; Otting, G.; Wu¨thrich, K. Cell 1994, 78, 211-223. (56) Tan, S.; Richmond, T. J. Nature 1998, 391, 660-666. (57) Welch, J. J.; Rauscher, F. J., III; Beerman, T. A. J. Biol. Chem. 1994, 269, 31051-31058. (58) Ho, S. N.; Boyer, S. H.; Schreiber, S. L.; Danishefsky, S. J.; Crabtree, G. R. Proc. Natl. Acad. Sci. U.S.A. 1994, 91, 9203-9207. (59) Chen, L. Curr. Opin. Struct. Biol. 1999, 9, 48-55. (60) Steinmetzer, K.; Reinert, K. E. J. Biomol. Struct. Dyn. 1998, 15, 779-791. (61) Chiang, S. Y.; Welch, J.; Rauscher, F. J.; Beerman, T. A. Biochemistry 1994, 33, 7033-7040. (62) Chiang, S. Y.; Azizkhan, J. C.; Beerman, T. A. Biochemistry 1998, 37, 3109-3115. (63) Henderson, D.; Hurley, L. H. Nat. Med. 1995, 1, 525-527. (64) Raumann, B. E.; Brown, B. M.; Sauer, R. T. Curr. Opin. Struct. Biol. 1994, 4, 36-43. (65) Dorn, A.; Affolter, M.; Mu¨ller, M.; Gehring, W. J.; Leupin, W. EMBO J. 1992, 11, 279-286. (66) Harrison, S. C.; Aggarwal, A. K. Annu. ReV. Biochem. 1990, 59, 933-969. (67) Becker, S.; Groner, B.; Mu¨ller, C. W. Nature 1998, 394, 145151. (68) Neidle, S. Molecular Aspects of Anticancer Drug-DNA Interactions; Macmillan Press: London, 1993. (69) Krugh, T. R. Curr. Opin. Struct. Biol. 1994, 4, 351-364. (70) Zakrzewska, K.; Lavery, R.; Pullman, B. Nucl. Acids Res. 1983, 11, 8825-8839. (71) Broggini, M.; Erba, E.; Ponti, M.; Ballinari, D.; Geroni, C.; Spreafico, F.; De´Incalci, M. Cancer Res. 1991, 51, 199-204. (72) Xu, Y.; Zhen, Y.; Goldberg, I. H. J. Am. Chem. Soc. 1997, 119, 1133-1134. (73) Xu, Y.; Zhen, Y.; Goldberg, I. H. Biochemistry 1997, 36, 1497514984. (74) Bellorini, M.; Moncollin, V.; De´Incalci, M.; Mongelli, N.; Mantovani, R. Nucl. Acids Res. 1995, 23, 1657-1663. (75) Walker, W. L.; Kopka, M. L.; Filipowski, M. E.; Dickerson, R. E.; Goodsell, D. S. Biopolymers 1995, 35, 543-553. (76) Chen, A. Y.; Yu, C.; Gatto, B.; Liu, L. F. Proc. Natl. Acad. Sci. U.S.A. 1993, 90, 8131-8135. (77) Wemmer, D. E.; Dervan, P. B. Curr. Opin. Struct. Biol. 1997, 7, 355-361. (78) Geierstanger, B.; Wemmer, D. E. Annu. ReV. Biophys. Biomol. Struct. 1995, 24, 463-493. (79) Clark, G. R.; Gray, E. J.; Neidle, S. Biochemistry 1996, 35, 1374513752. (80) Dervan, P. B. Science 1986, 238, 464-471. (81) Chaires, J. B. Curr. Opin. Struct. Biol. 1998, 8, 314-320. (82) Sen, S.; Nilsson, L. J. Am. Chem. Soc. 1998, 120, 619-631. (83) von Kitzing, E. Methods Enzymol. 1992, 211, 449-467. (84) Beveridge, D. L.; Swaminathan, S.; Ravishanker, G.; Withka, J. M.; Srinivasan, J.; Prevost, C.; Louise-May, S.; Langley, D. R.; DiCapua, F. M.; Bolton, P. H. Molecular Dynamics Simulations on the Hydration, Structure and Motions of DNA Oligomers. In Water and Biological Molecules; The Macmillan Press Ltd.: London, 1993; pp 165-225. (85) Beveridge, D. L.; Ravishanker, G. Curr. Opin. Struct. Biol. 1994, 4, 246-255. (86) Louise-May, S.; Auffinger, P.; Westhof, E. Curr. Opin. Struct. Biol. 1996, 6, 289-298. (87) Pearlman, D. A.; Case, D. A.; Caldwell, J. W.; Ross, W. S.; Cheatham, T. E., III; Ferguson, D. M.; Seibel, G. L.; Singh, U. C.; Weiner, P. K.; Kollman, P. A. AMBER 4.1; University of California, San Francisco, 1995. (88) York, D. M.; Wlodawer, A.; Pedersen, L. G.; Darden, T. A. Proc. Natl. Acad. Sci. U.S.A. 1994, 91, 8715-8718.

J. Phys. Chem. B, Vol. 105, No. 42, 2001 10387 (89) Cheatham, T.; Miller, J.; Fox, T.; Darden, T.; Kollman, P. J. Am. Chem. Soc. 1995, 117, 4193-4194. (90) Young, M. A.; Jayaram, B.; Beveridge, D. L. J. Am. Chem. Soc. 1997, 119, 59-69. (91) Duan, Y.; Wilkosz, P.; Crowley, M.; Rosenberg, J. M. J. Mol. Biol. 1997, 272, 553-572. (92) Shindo, H.; Fujiwara, T.; Akutsu, H.; Masumoto, U.; Shimidzu, M. J. Mol. Biol. 1984, 174, 221-229. (93) Sklenar, V.; Bax, A. J. Am. Chem. Soc. 1987, 109, 7525-7526. (94) Gorenstein, D. G. Methods Enzymol. 1992, 211, 254-286. (95) Gorenstein, D. G. Chem. ReV. 1994, 94, 1315-1338. (96) Szyperski, T.; Ono, A.; Fernandez, C.; Iwai, H.; Tate, S.; Wu¨thrich, K.; Kainosho, M. J. Am. Chem. Soc. 1997, 119, 9901-9902. (97) Hogan, M. E.; Jardetzky, O. Proc. Natl. Acad. Sci. U.S.A. 1979, 76, 6341-6345. (98) Roongta, V. A.; Jones, C. R.; Gorenstein, D. G. Biochemistry 1990, 29, 5245-5258. (99) Arnott, S.; Campbell-Smith, P. J.; Chandrasekaran, R. Atomic coordinates and molecular conformation for DNA-DNA, RNA-RNA, and DNA-RNA helices. In CRC Handbook of Biochemistry and Molecular Biology; CRC Press: Boca Raton, FL, 1976; pp 411-422. (100) Chandrasekaran, R.; Arnott, S. J. Biomol. Struct. Dyn. 1996, 13, 1015-1027. (101) Jorgensen, W. L. J. Am. Chem. Soc. 1981, 103, 335-340. (102) Case, D. A.; Pearlman, D. A.; Caldwell, J. W.; Cheatham, T., III; Ross, W. S.; Simmerling, C. L.; Darden, T. A.; Merz, K. M.; Stanton, R. V.; Cheng, A. L.; Vincent, J. J.; Crowley, M.; Ferguson, D. M.; Radmer, R. J.; Seibel, G. L.; Singh, U. C.; Weiner, P. K.; Kollman, P. A. AMBER 5; University of California, San Francisco, 1997. (103) Cornell, W. D.; Cieplak, P.; Bayly, C. I.; Gould, I. R.; Merz, K. M., Jr.; Ferguson, D. M.; Spellmeyer, D. C.; Fox, T.; Caldwell, J. W.; Kollman, P. A. J. Am. Chem. Soc. 1995, 117, 5179-5197. (104) Cheatham, T. E., III; Cieplak, P.; Kollman, P. A. J. Biomol. Struct. Dyn. 1999, 16, 845-862. (105) Winger, R. H.; Liedl, K. R.; Ru¨disser, S.; Pichler, A.; Hallbrucker, A.; Mayer, E. J. Phys. Chem. B 1998, 102, 8934-8940. (106) Hamelberg, D.; McFail-Isom, L.; Williams, L. D.; Wilson, W. D. J. Am. Chem. Soc. 2000, 122, 10513-10520. (107) Grzeskowiak, K.; Yanagi, K.; Prive´, G. G.; Dickerson, R. E. J. Biol. Chem. 1991, 266, 8861-8883. (108) MATLAB 5.3.1. The MathWorks, Inc., 24 Prime Park Way, Natrick, MA 01760, 1999. (109) Prive´, G. G.; Yanagi, K.; Dickerson, R. E. J. Mol. Biol. 1991, 217, 177-199. (110) Ru¨disser, S.; Hallbrucker, A.; Mayer, E. J. Am. Chem. Soc. 1997, 119, 12251-12256. (111) Tisne´, C.; Delepierre, M.; Hartmann, B. J. Mol. Biol. 1999, 293, 139-150. (112) Hartmann, B.; Piazzola, D.; Lavery, R. Nucl. Acids Res. 1993, 21, 561-568. (113) Lefebvre, A.; Mauffret, O.; Lescot, E.; Hartmann, B.; Fermandjian, S. Biochemistry 1996, 35, 12560-12569. (114) Bertrand, H. O.; Ha-Duong, T.; Fermandjian, S.; Hartmann, B. Nucl. Acids Res. 1998, 26, 1261-1267. (115) Wlodek, S. T.; Shen, T. S.; McCammon, J. A. Biopolymers 2000, 53, 265-271. (116) Schildbach, J. F.; Karzai, A. W.; Rauman, B. E.; Sauer, R. T. Proc. Natl. Acad. Sci. U.S.A. 1999, 96, 811-817. (117) Wellenzohn, B.; Flader, W.; Winger, R. H.; Hallbrucker, A.; Mayer, E.; Liedl, K. R. J. Am. Chem. Soc. 2001, 123, 5044-5049. (118) Kielkopf, C. L.; Ding, S.; Kuhn, P.; Rees, D. C. J. Mol. Biol. 2000, 296, 787-801. (119) Kielkopf, C. L.; Baird, E. E.; Dervan, P. B.; Rees, D. C. Nature Struct. Biol. 1998, 5, 104-109. (120) Kopka, M. L.; Fratini, A. V.; Drew, H. R.; Dickerson, R. E. J. Mol. Biol. 1983, 163, 129-146.