Nucleic Acids Research

Nucleic Acids Research 2004 32(14):4269-4280; doi:10.1093/nar/gkh765

This Article Abstract Print PDF (672K) Supplementary Material Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Add to My Personal Archive Download to citation manager Search for citing articles in: ISI Web of Science (67) Request Permissions Commercial Re-use Guidelines for Open Access NAR Content Google Scholar Articles by Várnai, P. Articles by Zakrzewska, K. Search for Related Content PubMed PubMed Citation Articles by Várnai, P. Articles by Zakrzewska, K. Social Bookmarking          What's this?

Published online 10 August 2004

Nucleic Acids Research, Vol. 32 No. 14 © Oxford University Press 2004; all rights reserved

DNA and its counterions: a molecular dynamics study

Péter Várnai and Krystyna Zakrzewska^*

Laboratoire de Biochimie Théorique, CNRS UPR 9080, Institut de Biologie Physico-Chimique, 13 rue Pierre et Marie Curie, Paris 75005, France

^* To whom correspondence should be addressed. Tel: +33 1 58415174; Fax: +33 1 58415026; Email: krystyna{at}ibpc.fr
Present address: Péter Várnai, University of Cambridge, Department of Chemistry, Lensfield Road, Cambridge CB2 1EW, UK

Received April 7, 2004; Revised June 15, 2004; Accepted July 26, 2004

	ABSTRACT

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS SUPPLEMENTARY MATERIAL REFERENCES

 
The behaviour of mobile counterions, Na⁺ and K⁺, was analysed around a B-DNA double helix with the sequence CCATGCGCTGAC in aqueous solution during two 50 ns long molecular dynamics trajectories. The movement of both monovalent ions remains diffusive in the presence of DNA. Ions sample the complete space available during the simulation time, although individual ions sample only about one-third of the simulation box. Ions preferentially sample electronegative sites around DNA, but direct binding to DNA bases remains a rather rare event, with highest site occupancy values of <13%. The location of direct binding sites depends greatly on the nature of the counterion. While Na⁺ binding in both grooves is strongly sequence-dependent with the preferred binding site in the minor groove, K⁺ mainly visits the major groove and binds close to the centre of the oligomer. The electrostatic potential of an average DNA structure therefore cannot account for the ability of a site to bind a given cation; other factors must also play a role. An extensive analysis of the influence of counterions on DNA conformation showed no evidence of minor groove narrowing upon ion binding. A significant difference between the conformations of the double helix in the different simulations can be attributed to extensive / transitions in the phosphate backbone during the simulation with Na⁺. These transitions, with lifetimes over tens of nanoseconds, however, appear to be correlated with ion binding to phosphates. The ion-specific conformational properties of DNA, hitherto largely overlooked, may play an important role in DNA recognition and binding.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS SUPPLEMENTARY MATERIAL REFERENCES

 
It has been known since the first experimental studies on DNA structure that both solvent and counterions play a major role in stabilizing the double helix and in determining its overall conformation (1). Solvent and counterion effects have also been evoked at a finer level in explaining the base-sequence-dependent structure of DNA oligomers, as observed in the single crystal structure of the so-called Drew–Dickerson dodecamer (2). This oligomer, containing the underlined EcoRI binding site (CGCGAATTCGCG), attracted particular attention because of the narrow minor groove formed by the central AATT sequence which constitutes a short A-tract, a feature associated both with conformational rigidity and with DNA bending. Several factors have been proposed to explain the unusual properties of A-tracts, subsequently observed in a number of different crystal structures. These include bifurcated hydrogen bonds, linking the adenine and thymine bases of adjacent base pairs (3), and the so-called ‘hydration spine’, a layer of water molecules forming a set of hydrogen bonds between the minor groove edges of the adenine–thymine bases, observed in the central part of the Drew–Dickerson dodecamer (4,5).

The structuring role of ordered water molecules in the first and higher solvation shells has been discussed at length, yet relatively little attention has been accorded to counterions. These ions are clearly necessary for DNA stability but are more difficult to study. Nevertheless, recent data obtained both by experimental and theoretical methods have described the association of monovalent cations with the double helix and their impact on DNA conformation. In crystal structures of the Drew–Dickerson dodecamer, it has been found that partially dehydrated monovalent cations (Na⁺ and K⁺) bind in the minor groove (6,7). Such sequence-specific cation binding has been interpreted as contributing to DNA conformational heterogeneity and linked notably to axial bending and to minor groove narrowing (6). A subsequent study of an atomic-resolution crystal structure of the same sequence grown under a higher Mg²⁺/Na⁺ ratio did not show such bound ions (8), and the results of the earlier works were challenged (9). However, Rb⁺ has been shown to bind at the central ApT step of the dodecamer with 50% occupancy (10), and Tl⁺ has been localized within the G-tract major groove (20–35%) and the A-tract minor groove (10%) in these crystals (11,12). Solution NMR studies of A-tract sequences show that monovalent ions penetrate into the minor groove (13), but the low occupancy binding argues against a strong structural influence on DNA (14).

Another way to approach the problem of counterion effects is by means of all-atom molecular dynamics simulations. Such simulations give us the possibility of following in detail the trajectories of individual counterions over multi-nanosecond timescales, and at the same time giving access to detailed information on DNA conformational dynamics. Although early simulations were severely limited due to both computer power and problems of instability, fast processors and the introduction of the particle–mesh Ewald (PME) approach to treat long-range electrostatic interactions have enabled simulations to be carried well beyond the nanosecond timescale. A number of specific questions concerning counterions in simulations, however, remain to be answered including the accuracy of ion–DNA interactions in solution and the sampling of the available conformational space around the DNA double helix (15).

Previous studies have addressed the question of direct, sequence-specific ion binding to DNA with contradictory results. Early molecular dynamics simulations described the movements of ions in the minor groove of the Drew–Dickerson dodecamer (16). In an extended 15 ns trajectory, it was found that the Na⁺ fractional occupancy in the AATT segment of the minor groove was 5–10%, while >10% occupancy was found at the CpG step in the major groove (17). The presence of sodium ions in the AATT region was found to exert a negligible influence on the width of the minor groove (with an apparent change of <1 Å, which is well below the thermal fluctuations of groove width), but was, however, found to be linked to major groove narrowing.

The same oligomer has been analysed in a 15 ns simulation in another study (18) with surprisingly different conclusions, although the only notable change in conditions involved a slightly different salt concentration. This analysis consisted of studying variations in groove width as a function of percentage site occupancies. The authors considered the effects of direct binding in the minor groove, especially at the central ApT step, and of binding at the entrance to the minor groove bridging the anionic oxygens of the two interstrand phosphates. It was concluded that the minor groove width was coupled to ion binding with a difference of up to 2 Å separating the cases of 0 and 100% occupancy (no error bounds were given).

Simulations of two 25mers containing A-tracts ATAGGCA₆TAGGCA₅TGG and G₅(GA₄T₄C)₂C₅ indicated a correlation between axial bending and the presence of a counterion close to the bending locus (17). Other AT-rich oligomers, C(AT)₄G, CA₄T₄G and CT₄A₄G, have also been simulated for 8 ns each (19), leading to the conclusion that Na⁺ intrusions into the first hydration shell of the minor groove of A-tracts occurred but were not coupled to axial bending.

Studies of a decamer C₅T₅ (20) described Na⁺ binding to guanines within the major groove and along the hydration spine in the minor groove of the AT base pair sequence, and reported a small (0.7 Å) reduction in groove width linked to ion binding within the AT segment. No such coupling was found for either monovalent ions or spermine in simulations of three non-A-tract B-DNA decamers in periodic conditions chosen to model a crystalline environment (21). Simulations of alternating (CG)₁₂ and (AT)₁₂ sequences with K⁺ were reported to result in direct major groove binding at purine/pyrimidine (RpY) steps, but no binding was seen in the minor groove (22,23). Finally, dependence of ion distribution on the force field used was noted for the CGATTAATCG decamer, but these simulations were too short (1.2 ns) to achieve ion convergence (24).

Although this summary is not exhaustive, it shows clearly that the influence of counterions on DNA conformation remains an open question. All the sequences studied previously had special features including either A-tracts, homopolymeric runs or alternating sequences. Therefore, it is important to gain insight into the behaviour of counterions around a sequence, which does not show any of these properties. In an attempt to answer this question, we have carried out two 50-ns simulations of a double-stranded B-DNA dodecamer, CCATGCGCTGAC, containing the underlined target sequence for HhaI cytosine methyltransferase, which we have previously used to study base pair opening (25). We note that our simulation is far longer than current standards. This is particularly important as it has been observed that ions can get trapped in particular binding sites for nanoseconds (18). The two simulations (using equivalent initial conditions and simulation protocols) were carried out with different monovalent cations, Na⁺ and K⁺, both of which are important for DNA stability and are commonly used in experimental studies. Ions can also serve as sensitive electrostatic probes to map the dynamic surface of DNA for macromolecular recognition. In this study, we analyse how well ions sample the space around DNA, their direct site binding and diffusion properties, and concurrent changes in DNA structure.

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS SUPPLEMENTARY MATERIAL REFERENCES

 
Simulation protocol
Model building and simulations were performed using the AMBER 6 suite of programs (26) and the parm99 parameter set (27). Simulations were carried out on a dodecamer with the sequence d(CCATGCGCTGAC) which was constructed in a standard B-DNA conformation. This oligomer was neutralized with 22 Na⁺ or K⁺ ions (placed following electrostatic potential values) and solvated with more than 5000 TIP3P water molecules in a truncated octahedral box, which allowed for at least a 10 Å shell of water around the solute. The minimal salt conditions were chosen because the proper description of ions in simulations represents highly controversial parts of current force fields (28–30). Molecular dynamics simulations were then carried out at constant temperature (300 K) and pressure (1 bar) using the Berendsen algorithm (31) and periodic boundary conditions. An integration time step of 2 fs was used and all bond lengths involving hydrogens were constrained using SHAKE (32). Long-range interactions were treated using the PME approach (33) with a 9 Å direct space cut-off. The non-bonded pair list was updated heuristically, whenever any atom has moved more than half of a ‘skin’ distance of 2 Å, and the centre of mass motion was removed every 10 ps during the simulation. Initially, the water molecules and ions were relaxed by energy minimization and allowed to equilibrate at 300 K around the fixed DNA oligomer for 100 ps at constant volume. The entire system was then heated from 100 to 300 K during 10 ps and equilibrated for 40 ps with harmonic restraints of 5.0 kcal/mol/Ų on the solute atoms at constant volume. Subsequently, the simulation was continued at constant pressure. The restraints were gradually removed over a period of 250 ps and the system was equilibrated for 1 ns. The following 50 ns of unrestrained simulation was used for analysis.

Counterion analysis
In order to analyse ion movements around the solute in a periodic boundary simulation, the ions that move far away need to be first translated back to the primary simulation cell. These ‘imaging’ calculations were performed with the ptraj utility of the AMBER package. We next analysed what percentage of the available space in the simulation box was effectively sampled by counterions during the course of the simulation. The volume occupied by the solute was calculated numerically using a grid with a 0.2 Å spacing and taking into account the van der Waals radii of the DNA atoms (H, 1.0 Å; C, 1.7 Å; N, 1.5 Å; O, 1.4 Å; P, 1.9 Å) incremented by the radius of the ion (Na⁺, 1.0 Å; K⁺, 1.4 Å).

The volume sampled by the counterions was calculated using the same grid spacing and ionic radius. An analysis of the ion displacement with 1 ps time intervals, corresponding to the saving frequency of coordinates, showed that the mean displacement was 1 Å, although displacements of up to 5 Å were also observed (excluding any effects caused by the imaging procedure mentioned above). Ion positions not sampled due to the use of 1 ps steps were taken into account by generating intermediate points. One such point was created by linear interpolation when the distance between two consecutively sampled points was between 2.5 and 3 Å and two equally spaced points were created for distances between 3 and 5 Å.

Diffusion coefficients
Motions of the counterions in the presence of DNA were characterized by their diffusion coefficients (D) according to Einstein's formula:

This requires calculation of the mean square displacements of each ion for a series of time intervals. The diffusion coefficient is then simply given by the slope of the plot of the average square displacements versus time. The time intervals used in this calculation were varied from 10 to 1000 ps in steps of 10 ps. For each calculation, the displacement was averaged over all time intervals along the trajectory.

Binding site analysis
An ion was considered bound to the solute if its distance from an electronegative solute atom was <3.0 and 3.5 Å for Na⁺ and K⁺, respectively. This value corresponds to a minimum in the radial distribution function of the ions around the electronegative solute atoms, which separates a sharp maximum corresponding to a direct binding and a secondary maximum corresponding to water-mediated interactions (20,22).

Site occupancy was calculated for both base- and backbone-binding sites on DNA. There are two potential binding sites in the minor groove, Y(O2) and R(N3) (where Y = pyrimidine and R = purine), and three in the major groove, R(N7), G(O6) and T(O4). For the backbone, we considered all five oxygen atoms: the phosphate anionic oxygens, O1P and O2P, the phosphate esteric oxygens, O5' and O3' and the sugar-ring oxygen, O4'. For each of these sites, we calculated the occupancy, defined as a percentage of the simulation time during which a counterion was directly bound. We also recorded the average and the maximum durations of binding. For the sites with the highest occupancies we calculated the ion residence times using a standard time correlation function (34), with residence time distribution (rtd) for a given binding site defined as

where p_i(t) is 1 if the ion i is present at time t and otherwise 0. Events when an ion fluctuates across the boundary within 1 ps are ignored. Calculations were carried out with time intervals ranging from 10 to 1000 ps. The effective residence time (R) is then calculated from the decay coefficient of an exponential function fitted to rtd:

Counterion influence on DNA conformation
In order to investigate the possible impact of ion binding on DNA conformation we calculated both helicoidal parameters and axial-bending angles throughout the simulation using the Curves program (35). Minor groove widths were measured using two different approaches. First, we calculated the interstrand distance between the phosphorous atoms P_i and P'_i+3 across the minor groove (d_PP). This distance was reduced by 5.8 Å (sum of the van der Waals radii of the two phosphate groups) in order to be comparable with the convention adopted by other authors. Second, we used the minor groove width calculated by the Curves program, which corresponds to the shortest distance between splines passing through the backbone phosphorus atoms and lying in a plane passing through the pseudodyad vector of the chosen base pair (36). DNA structures were sorted into groups depending on the presence or absence of counterions at particular locations around the DNA, and the average groove width and standard deviation were calculated for each group. The following groups were considered: (i) counterions bound in the minor groove associated with a given base pair; (ii) counterions bound in the major groove associated with a given base pair; (iii) counterions between two phosphate groups across the minor groove; and (iv) no bound counterions. Note that structures satisfying more than one criterion were rare and were subsequently excluded from our analysis.

To find structures belonging to the third group, we considered ions falling at less than (2 + d_PP/2) Å from the phosphorous atoms, P_i and P'_i+3, which defines the local minor groove width. This approach defines a convex-lens-shaped volume equidistant from the phosphate groups at the entrance to the minor groove. An alternative approach, where the volume of the lens is constant, was also tested and gave very similar results. We also calculated the relationship between the minor groove width and ion site occupancy in order to compare our results with those of Hamelberg and co-workers (18). This involved calculating the average groove width and the corresponding groove and phosphate ion occupancies for all 10 ps segments during the total simulation time.

Molecular electrostatic potential
Calculation of the electrostatic potential around the DNA oligomer was carried out with the program Delphi (37). The program was adapted to calculate an average value of the electrostatic potential over the solvent accessible area assigned to a given atom. A grid of 145 points was used and calculations were performed with radii of 1.2 or 1.4 Å. The internal and external dielectric constants were set to 2 and 80, respectively, and the ionic force was zero.

	RESULTS AND DISCUSSION

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS SUPPLEMENTARY MATERIAL REFERENCES

 
Counterion convergence
Let us consider first the problem of how efficiently ions sample the simulation box during the available simulation time. The total volume of the simulation box fluctuated around 1.79 x 10⁵ ų (±1500 ų) as a consequence of the constant pressure adjustments, and the volume excluded by the DNA was 10920 ų (±60 ų). Figure 1A presents the volume visited by the total number of counterions as a function of simulation time. The curves shown correspond to the imaging of ions with the solute molecule centred in the box. The molecule can tumble during the simulation and this movement changes the apparent volume occupied by the solute to some extent, however, this change is small in comparison with the volume of the box (see below). The figure shows that the volume sampled by K⁺ increases faster than that for Na⁺. The main increase takes place during the first 10 ns of the simulation time and a large fraction of the total volume is already sampled (practically 100% for K⁺ and 80% for Na⁺). This increase then reaches a plateau value at 1.73 x 10⁵ ų, 98% of the total box volume. The difference of 6000 ų between the total volume of the simulation box and the volume visited by counterions is slightly lower than the excluded volume of the solute because the tumbling of the DNA opens some new space for the ions.

View larger version (19K): [in this window] [in a new window]   Figure 1. Volume sampled by Na+ (dotted line) and K+ (solid line) counterions during 50 ns simulation time calculated for all the ions (A) and for the slowest and fastest ions (B).

 
Although virtually the whole simulation box is visited by at least one ion, individual ions sample only 30–40% of the available space. Figure 1B illustrates this for the ‘slowest’ and ‘fastest’ moving ions. While the fast ion moves through the box rather uniformly, the trajectory of the slower ones is characterized by long pauses. For example, the slowest Na⁺ contacts many points on the DNA molecule; starting from a location close to T₉ it moves to the opposite end of the helix, and 6 ns later, it is back at phosphate 10. It then contacts the terminal base G₁₃, before moving back to the centre of the oligomer. Owing to such pauses, the total volumes sampled by individual ions can vary by >10⁴ ų. In the simulation with K⁺, individual ions sample a higher fraction of the total volume and the curves show many short pauses with more uniform sampling of space.

An analysis of sampling of separate volume elements around DNA was also performed. This calculation shows how cylindrical shells with increasing radial distances from a linear DNA axis are visited along the simulation time (Supplementary Figure S1). A maximum in visited volume is observed at a distance of 14 Å from the DNA axis, which corresponds to a hydrated ion close to the surface of the DNA.

Ion diffusion coefficient
Movements of the ions during the simulation can be characterized by their diffusion coefficient (D), although this must be considered an approximate measure due to the existence of binding events. Nevertheless, the mean-square displacement averaged for all the ions over time intervals gives a perfect line. The slope corresponds to D = 1.7 x 10^–9 m² s^–1 for Na⁺ and 2.8 x 10^–9 m² s^–1 for K⁺. Values for individual Na⁺ ions vary between 1.3 and 2.1 x 10^–9 m² s^–1, whereas those for K⁺ are between 2.5 and 3.3 x 10^–9 m² s^–1.

Our average value for Na⁺ is identical to that reported in earlier DNA simulations under similar conditions (38), although lower values (1.3 x 10^–9 m² s^–1) have been calculated from simulations under different conditions, and, notably, using a flexible, simple point-charge water model and the CHARMM 22 force field for DNA (21). In comparison, the experimental value for Na⁺ in pure water is 1.33 x 10^–9 m² s^–1 (39), which has been well reproduced by simulation (1.28 x 10^–9 m² s^–1) (40). Diffusion properties of K⁺ in the presence of DNA have not, to our knowledge, been published. The experimental diffusion constant of K⁺ in pure water is 1.96 x 10^–9 m² s^–1 (39), and a value of 1.86 x 10^–9 m² s^–1 was calculated from simulations (40). Our calculated diffusion coefficients in the presence of DNA show a similar Na⁺/K⁺ ratio as those found in aqueous solution.

Ion mobility was calculated as a function of the mean distance of ions from the DNA axis (Supplementary Figure S2). This shows that K⁺ has a mobility of 1.2 Å/ps on average that is roughly independent of its distance from DNA. In contrast, Na⁺ mobility is low in the vicinity of DNA (0.75 Å/ps), characteristic of direct binding, and then increases to 1.1 Å/ps beyond 15 Å from the DNA axis.

Ions in the grooves
To get a clear picture of where the ions were most frequently located during the simulation, each snapshot was re-oriented by superposing the solute and the ion occupation of each volume element was counted using a 1 ų grid. Locations in the simulation box that correspond to >20% of the maximum occupancy along with the average DNA structure are shown in Figure 2.

View larger version (58K): [in this window] [in a new window]   Figure 2. Counterion distribution around an average DNA structure for Na+ (left) and K+ (right). The positions shown correspond to roughly 700 to 3500 visits per 1 Å3 over 50 ns simulation (top 80% of ion density) with red for the highest, white for the medium and blue for the lowest densities.

 
Both monovalent ions, as would be expected, preferentially visit volume elements in the vicinity of (but not necessarily directly bound to) the double helix. This does not imply that individual ions spend most of their time around the DNA and this is indeed not the case. The figure shows that, in the case of Na⁺, the minor groove has higher ion density than the major groove, even for this GC-rich sequence, although there is a favoured site in the central region of the major groove of the oligomer. The highest occupancies are located in the minor groove, mainly around T₉(O2), shown in red at the lower half of the oligomer (Figure 2). A different picture emerges from simulations with K⁺ where ions preferentially visit the central part of the oligomer. The highest ion density locations for K⁺, however, correspond to the major groove, specifically to the GpC steps in the central GCGC recognition sequence of the oligomer.

Calculated percentage ion occupancies for direct binding to bases in the grooves and to phosphates are summarized in Table 1. Note that occupancy values for various sites are not strictly additive since an ion can simultaneously interact with more than one site. The most important feature is the significantly different binding pattern for the two monovalent ions. Sodium binding to DNA is highly non-uniform along the sequence. High occupancies were found in both the minor and major grooves of the G₇–A₁₁ segment as well as around certain phosphate groups. The most occupied sites are located away from the centre of the oligomer, in contrast to what might be expected from an oligomeric electrostatic potential distribution or from sequence-dependent electrostatic properties. While direct binding to T₉(O2) dominates with a 12% occupancy, many base sites have occupancy values <1%. Interestingly, much lower occupancies were found at the ApT step at the opposite end of the oligomer, in contrast to results seen in previous studies for central ApT step in the Drew–Dickerson dodecamer (7,16). The different sequence contexts of AT pairs explain the different ion-binding patterns. Ions that are shared between different sites are not very frequent, mainly at G(O6) in the major groove and at C/T(O2) in the minor groove. All phosphate groups, and specifically the O1P site, which lies on the minor groove side of the duplex, have direct ion binding for at least 2% of the simulation time, and these values increase to 10% for phosphates close to bases with high ion occupancies. The single exception is P₂₁ with 17% direct binding, which is mainly due to the unusually high occupancy in the final 15 ns of simulation time. An analysis of the structure of bound counterions at this site revealed a rare coordination of two Na⁺ around O1P, a configuration also seen in the crystal structure of the disodium salt of adenosine triphosphate (41).

View this table: [in this window] [in a new window]   Table 1. Binding site occupancy (%) of counterions

 
In the case of K⁺, the distribution of the ion-binding sites appears to reflect the electrostatic potential distribution of an oligomeric DNA fragment. The most occupied sites are located in the centre of the sequence with a strong preference for the major groove. The highest occupancy is 13% at G₇(O6), a similar value to that observed for Na⁺ at T₉(O2). High K⁺ occupancy between interstrand G(O6) in GpC steps has been observed in other studies (22,23). Direct binding in the minor groove also takes place, although with low occupancy values. Binding to phosphates is uniform with occupancies varying between 1 and 6%.

It should be noted that overall the major groove is better than the minor groove for direct binding of either counterion. In relative terms, however, there is more Na⁺ binding in the minor groove compared with K⁺. The accumulated sodium occupancy (calculated from analysis where an ion is assigned to the closest groove atom) is 23% in the minor groove and 36% in the major groove. For potassium, these values are 14 and 45%, respectively. The preferential Na⁺ binding site at T₉(O2) in the minor groove provide an explanation for this fact. Interestingly, despite the different binding patterns of the two counterions, the average number of cations directly bound to base sites of the oligomer is roughly equal (1.2 cations per nanosecond) in the two simulations.

It is important to remark that every ion spends some time directly bound to DNA sites during the 50 ns of simulation. The reliability of our results, even for the ‘stickier’ Na⁺, is supported by the fact that many independent binding events occur at a given site along the trajectory. For example, sites on base pair T₉–A₁₆ are visited by 9 different Na⁺ ions in the minor groove and 16 different ions in the major groove. Interestingly, simultaneous occupancy of both grooves at this base pair level, as well as simultaneous binding to base sites in the same groove, were observed several times. Nevertheless, such binding events represent only a small fraction (1–3%) of the total simulation time and hence they do not influence the normal behaviour of B-DNA. On the other hand, Na⁺ can have an important structural role in other nucleic acid structures, such as RNA loops (42,43) or the G-DNA quadruplex (44,45), with correlated binding and high occupancies (>50%).

To characterize the dynamics of binding of Na⁺ and K⁺ to DNA, we calculated ion residence times, the number of binding events, the average direct-binding time and the longest duration of binding for sites with high occupancies (Table 2). As noted above, the site T₉(O2) is exceptionally attractive for Na⁺ with relatively few binding events, but with long residence times (381 ps) and a maximum binding duration of 988 ps. Other base sites show more frequent Na⁺ binding–unbinding events (between 65 and 149) associated with considerably shorter residence times (from 73 to 143 ps). Interestingly, the most occupied phosphate sites have residence times comparable with those of the groove sites, ranging from 147 to 372 ps. Overall, these values are of the same order as those obtained in other studies for Na⁺ (20,21). Importantly, however, the binding properties of K⁺ are very different from those of Na⁺. For similar total site occupancies, there were several hundred binding events with much shorter average binding times, of the order of 10 ps. This corresponds well with the picture of a fast-moving K⁺ around DNA that was also indicated by its higher diffusion coefficient.

View this table: [in this window] [in a new window]   Table 2. Ion residence time , number of binding events (nb), the average direct binding time and the longest duration of binding (ps) for selected binding sites

 
It is informative to look in detail at the evolution of the occupancy values along the 50 ns simulation time. This is especially important in judging previous, shorter molecular dynamics simulation studies of counterion behaviour around DNA. Binding occupancies are shown for consecutive 5 ns segments of the trajectories in Table 3. Cumulative occupancy values can be calculated from these data for various simulation lengths. It is remarkable that Na⁺ binding to DNA, which is restricted to a few sites, shows very irregular results based on the 5 ns occupancy sampling. For example, G₅(O6) is occupied only during the first and the last 5 ns, and C₁₂(O2) is occupied significantly only after 35 ns. The highest occupancy calculated for a single 5 ns segment is 33% at T₉(O2). In contrast, K⁺ binds to more base sites on DNA and the time distribution is more uniform. High occupancies are observed in the central part of the oligomer, specifically around G₇ and G₁₇. The highest value of 21% was calculated for G₁₇(O6). It is important to note that no bias was visible in the site occupancy values in the production phase of the simulation attributable to the initial (electrostatic potential based) placement of counterions.

View this table: [in this window] [in a new window]   Table 3. Binding site occupancy (%) of counterions for selected groove atoms averaged over 5 ns intervals

 
Binding occupancies of Na⁺ around phosphates again show large fluctuations along the trajectory (Table 4). Several phosphates, such as P₇, P₁₀, P₁₈ and P₂₁, have occupancy values well above 30% for a given 5 ns segment. The rare coordination of two Na⁺ around O1P mentioned above is responsible for an occupancy of 66% between 41 and 45 ns of simulation time. For K⁺, the phosphate occupancy values are again more uniform, rarely >10%, with the single exception of 18% for P₁₆ between 36 and 40 ns.

View this table: [in this window] [in a new window]   Table 4. Binding site occupancy (%) of counterions for selected phosphates averaged over 5 ns intervals

 
Influence of ion binding on DNA groove widths and bending
In order to analyse the effect of ions on the conformation of the double helix, we calculated the minor groove width and bending properties, both considered to be sensitive measures of DNA conformation, along the trajectories. Structures were sorted into groups depending on the presence of ions in the minor or major groove, between the phosphates at the entrance to the minor groove, or the absence of ions around the given base pair level. The minor groove width was measured either as an interstrand distance between the two phosphorous atoms, P_i and P_i+3, or by using the Curves algorithm. Both methods gave similar results and hence we present only the average P–P distances along the oligomer for the different structural groups obtained from simulations with either Na⁺ or K⁺ (Figure 3).

View larger version (22K): [in this window] [in a new window]   Figure 3. Minor groove width along the sequence calculated for structures sorted according to their counterion positions from simulation with Na+ (A) and K+ (B). In green is the average width for structures with ions located in the minor groove of the corresponding base pair, black for ions in the major grove, in blue for ions located between phosphates and average values with standard deviations in red correspond to structures without ions in the given zone. The groove width was measured as a P–P distance reduced by 5.8 Å, the sum of the van der Waals radii of two phosphates.

 
The average groove width along the sequence is rather constant with the narrowest part located at the C₈–G₁₇ base pair level for both counterions. An interesting observation is that the narrowest site of the oligomer does not correspond to the highest occupancy site for either ion. Instead, counterions bind to the neighbouring base pairs, at the T₉ and G₇ levels in the case of Na⁺ and K⁺, respectively. Binding of ions in the minor groove appears to slightly narrow the groove width, an effect varying from 0 to 1.5 Å along the sequence for both Na⁺ and K⁺. Ions in the major groove seem to have a small opening effect on the minor groove, while the presence of counterions between phosphates at the minor groove entrance does not seem to have any consistent influence on the minor groove width. A comparison of the influence of the different counterions shows that the minor groove width is somewhat narrower for Na⁺, along the entire oligomer. However, variations between the structural groups with a given ion at any base pair level or between the different ions are well within the large fluctuations of the minor groove width in the absence of bound counterions. This implies that the flexibility of the minor groove alone could explain the narrow structures seen, without evoking the presence of a bound counterion.

Hamelberg and co-workers (18) have calculated the dependence of the minor groove width on the site occupancy, and obtained a linear dependence with a difference of 2 Å between 0 and 100% occupancy. We repeated this calculation at the narrowest part of our oligomer, around the C₈–G₁₇ base pair level, where recurring Na⁺ binding resulted in long residence times and high occupancy. Our analysis, however, did not lead to any statistically significant correlation for either the base or the phosphate sites. It is important to note in this respect that even during the longest Na⁺ binding event of almost 1 ns at T₉(O2), no clear minor groove narrowing was seen. This might be attributable to the fact that during the same period, several binding events took place in the major groove, which indicates a dynamic equilibrium of monovalent ions in the DNA grooves rather than the exclusive minor or major groove binding in this sequence.

Asymmetric charge neutralization is expected to result in DNA axis bending (46), and DNA bending was found to be influenced by groove-bound counterions in long oligomers containing A-tracts (17). We also analysed the DNA curvature in our sequence using the local helical axis bending angles calculated by the Curves algorithm for each base pair level along the simulation. Bending angles were sorted into groups depending on the presence or absence of ions at the base pair. Our analysis finds only very small effects with differences of 1–2°, compared to standard deviations in bending of 3°. This reinforces the finding that the sequence studied here behaves differently from A-tract and G-tract sequences, for which the predominance of ion binding in only one of the grooves was suggested (47).

Finally, using an alternative approach to the question, we selected two groups of 1000 structures having either a very narrow or a very wide minor groove at the level of the C₈–G₁₇ base pair. The DNA structures were superposed using P₁₀ and P₁₉ atoms and the midpoint of the C1'–C1' vector of the C₈–G₁₇ base pair, and ion positions were recorded if they fell within the direct-binding distance from this plane. We then analysed the ion distributions for these two groups (Figure 4). For Na⁺, the two groups correspond to structures with d_PP < 7.9 Å and d_PP > 14.0 Å, for the narrow and wide structures, respectively. Roughly half of the narrowest structures have counterion in the minor groove, while 330 ions were located in the minor groove of the widest structures. For K⁺, the limiting d_pp for the narrowest and widest minor groove corresponds to 8.6 and 15.0 Å, respectively. There were only 160 and 120 counterions in the minor groove of the narrowest and widest structures, respectively. These results show that although fewer counterions are located in the widest minor grooves, the narrowest minor grooves appear without the presence of any counterion in more than half of the cases.

View larger version (64K): [in this window] [in a new window]   Figure 4. Counterion distribution in the minor groove of an average DNA using the 1000 narrowest (left) and widest (right) minor groove structures at site C8–G17 from simulations with Na+ (top) and K+ (bottom). Phosphates that were used to measure the minor groove width, P10 and P19, are indicated with black spheres.

 
An important factor in determining the counterion distribution is the electrostatic potential of DNA. We therefore calculated the surface electrostatic potential using average DNA structures corresponding to either the narrowest or the widest minor grooves discussed above. The results indicate that the potential of electronegative atoms in the minor groove gets more negative with the narrowing of the groove. A patch of low electrostatic potentials was seen at the level of C₆–T₉ in the minor groove of the narrow structure obtained with Na⁺. This area is significantly more negative (–23 kT/e) than the potential in the major groove (–17 kT/e). A similar patch, but with a potential value reduced to –19 kT/e, was observed for the narrow minor groove with K⁺. Finally, the potential in the minor and major grooves become comparable for the wide minor groove structures. Therefore, the preference of K⁺ to enter the major groove cannot be explained by the electrostatic potential of these structures, other factors, such as the size of the hydrated ion and the desolvation free energy must also play a role.

Influence of ion binding on DNA backbone conformations
We have shown that the different types of monovalent counterions exhibit different binding patterns around DNA and that the minor groove widths are somewhat narrower in the case of Na⁺. Nevertheless, no direct coupling of bound ions and groove narrowing could be deduced from our simulations. The narrowest minor groove was located at C₈ in simulations with either counterions and hence this appears to be a feature of this particular sequence rather than an effect of direct ion binding. The analysis of helicoidal parameters and backbone torsions, however, revealed an interesting difference between the two simulations.

Minor groove width variations at level C₈ along the simulation time showed that while the groove width is constant in the simulation with K⁺, it narrows significantly in the second half of the simulation with Na⁺ (Figure 5). Similarly, the twist between the levels C₈/T₉ is roughly constant at 34° in the simulation with K⁺; however, it drops to 20° after a 15 ns simulation with Na⁺ (Supplementary Figure S3). We compare variations in backbone torsions and for the two simulations in Figure 6. The most striking observation is that roughly half of the canonical and backbone torsions (gauche^– and gauche⁺, respectively, in canonical B-DNA) transit to alternative conformations during the simulation with Na⁺. These mainly include coupled transitions to /:g⁺/t, but /:g^–/t conformations are also seen as a result of flipping to trans. Incidentally, most backbone torsions occur in the G₇–A₁₁ part of the oligomer where increased ion binding was also found. Most transitions occurred after a simulation time of 20 ns, which would explain why this conformational heterogeneity was not discussed in earlier theoretical studies based on shorter molecular dynamics trajectories. Some of the alternative / states flipped back to the canonical form after 10–15 ns, but the rest remained stable until the end of the simulation time. Such / transitions appear to be coupled to reduced twist and minor groove width in the neighbouring base pair steps. A comparison of Table 4 and Figure 6 suggests that the counterion binding to phosphates and the presence of non-canonical backbone are correlated. Curiously, these alternative backbone conformations appeared only in very few instances during the simulation with K⁺.

View larger version (17K): [in this window] [in a new window]   Figure 5. Minor groove width at level C8 as a function of simulation time. Average values over 5 ns with standard deviation are given for simulation with Na+ (dotted line) and K+ counterions (solid line).

 

View larger version (16K): [in this window] [in a new window]   Figure 6. Non-canonical / backbone torsions (gauche+/trans) as a function of time for simulation with Na+ (A) and K+ (B); in solid line and in dotted line.

 
Non-canonical / states have been observed in crystal structures of many protein–DNA complexes, but only in some particular cases in data of B-DNA oligomers. The spontaneous appearance of alternative / conformations has been reported in early molecular dynamics simulations of B-DNA, with Nuclear Overhauser Effect distance restraints (48,49). A more recent NMR study found in trans conformation (50). The sequence studied in the present work has been the subject of an earlier study in our laboratory where the structural and energetic properties of / transitions in the central GpC step were analysed (25). It was found that the lowest free-energy barrier is >7 kcal/mol, and hence the spontaneous occurrence of alternative forms is unlikely. However, the present study clearly shows that backbone transitions are strongly counterion dependent and this raises the question of how the barrier heights are affected by the presence of direct cation binding. Another important problem is the evaluation of this backbone flexibility in different sequence contexts. Undoubtedly, more stringent evaluation of force-field parameters as well as more solution-phase experimental data are needed to clarify this point.

	CONCLUSIONS

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS SUPPLEMENTARY MATERIAL REFERENCES

 
The behaviour of two commonly studied monovalent counterions, Na⁺ and K⁺, has been analysed around a B-DNA dodecamer containing the target sequence for HhaI cytosine methyltransferase using 50 ns long molecular dynamics trajectories. Our results show that the averaged movement of counterions in the presence of DNA remains diffusive, the movement of K⁺ being much faster than that of Na⁺. Ions sample the complete space available around the solute during the simulation time, but individual ions sample only about one-third of the space.

Ions preferentially visit electronegative sites around DNA, although direct binding to DNA bases is still a rather rare event, with highest site occupancy values not >13%. All ions spend some time in the DNA grooves and the most occupied sites are visited by many different ions during the simulation. The preferred binding sites depend greatly on the nature of the ion: Na⁺ binds mostly in the minor groove around T₉(O2), while K⁺ binds in the centre of the major groove. Interestingly, the A₃pT₄ step, considered as a likely binding pocket in the minor groove, does not show strong ion binding. This indicates that the ion-binding sites may be sensitive to their sequence contexts. Phosphate site occupancy is more uniform along the sequence, but phosphates proximate to bases with high ion occupancy generally exhibit higher fractional occupancy. The Na⁺ and K⁺ exhibit very different binding dynamics. Residence times for Na⁺ range from 100 to 400 ps for highly occupied sites with the longest binding duration reaching 1 ns. In contrast, K⁺ binds and dissociates more frequently, resulting in residence times of <70 ps. These features, along with its dominance in the intracellular environment, make K⁺ particularly well adapted for simulations of DNA.

Using various analysis techniques, we have found no evidence for coupling of monovalent ion binding and changes in local helix geometry. Important structural features of the oligomer (such as the narrow minor groove at C₈–G₁₇) were present in simulations with either counterion and between different groups of structures. Therefore, the structural properties of DNA seem to be coded in its base sequence. Nevertheless, we found differences in the DNA structure depending on the type of monovalent ion used in the simulations. Specifically, non-canonical / backbone conformations appeared mainly in the simulation with Na⁺ and were responsible for reduced twist and narrower minor groove. These / transitions produce long-lived conformational substates on the nanosecond time-scale, which require particular attention in future studies. It cannot be excluded that these conformational substates arise as an artefact of the force field. However, they might have been simply overlooked in experimental studies. The final answer to this question can only come from specially designed experimental studies of nucleic acid backbones.

	SUPPLEMENTARY MATERIAL

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS SUPPLEMENTARY MATERIAL REFERENCES

 
Supplementary Material is available at NAR Online.

	ACKNOWLEDGEMENTS

 
The authors thank Richard Lavery for helpful discussions and the CINES (France) computer centre for access to their facilities. P.V. was a Wellcome Trust travelling research fellow in Paris during this work (grant 060078).

	REFERENCES

TOP ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION CONCLUSIONS SUPPLEMENTARY MATERIAL REFERENCES

 

1. Saenger,W. ( (1984) ) Principles of Nucleic Acid Structure. Springer, NY.

2. Dickerson,R.E. and Drew,H.R. ( (1981) ) Structure of a B-DNA dodecamer: influence of base sequence on helix structure. J. Mol. Biol., , 149, , 761–786.[CrossRef][Web of Science][Medline]

3. Nelson,H.C., Finch,J.T., Luisi,B.F. and Klug,A. ( (1987) ) The structure of an oligo(dA)·oligo(dT) tract and its biological implications. Nature, , 330, , 221–226.[CrossRef][Medline]

4. Drew,H.R. and Dickerson,R.E. ( (1981) ) Structure of a B-DNA dodecamer. III. Geometry of hydration. J. Mol. Biol., , 151, , 535–556.[CrossRef][Web of Science][Medline]

5. Kopka,M.L., Fratini,A.V., Drew,H.R. and Dickerson,R.E. ( (1983) ) Ordered water structure around a B-DNA dodecamer. A quantitative study. J. Mol. Biol., , 163, , 129–146.[CrossRef][Web of Science][Medline]

6. Shui,X., Sines,C.C., McFail-Isom,L., VanDerveer,D. and Williams,L.D. ( (1998) ) Structure of the potassium form of CGCGAATTCGCG: DNA deformation by electrostatic collapse around inorganic cations. Biochemistry, , 37, , 16877–16887.[CrossRef][Medline]

7. Shui,X., McFail-Isom,L., Hu,G.G. and Williams,L.D. ( (1998) ) The B-DNA dodecamer at high resolution reveals a spine of water on sodium. Biochemistry, , 37, , 8341–8355.[CrossRef][Medline]

8. Tereshko,V., Minasov,G. and Egli,M. ( (1999) ) The Dickerson–Drew B-DNA dodecamer revisited at atomic resolution. J. Am. Chem. Soc., , 121, , 470–471.[CrossRef]

9. Chiu,T.K., Kaczor-Grzeskowiak,M. and Dickerson,R.E. ( (1999) ) Absence of minor groove monovalent cations in the crosslinked dodecamer CGCGAATTCGCG. J. Mol. Biol., , 292, , 589–608.[CrossRef][Web of Science][Medline]

10. Tereshko,V., Minasov,G. and Egli,M. ( (1999) ) A ‘hydrat-ion’ spine in a B-DNA minor groove. J. Am. Chem. Soc., , 121, , 3590–3595.[CrossRef][Web of Science]

11. McFail-Isom,L., Sines,C.C. and Williams,L.D. ( (1999) ) DNA structure: cations in charge? Curr. Opin. Struct. Biol., , 9, , 298–304.[CrossRef][Web of Science][Medline]

12. Howerton,S.B., Sines,C.C., VanDerveer,D. and Williams,L.D. ( (2001) ) Locating monovalent cations in the grooves of B-DNA. Biochemistry, , 40, , 10023–10031.[CrossRef][Medline]

13. Hud,N.V., Sklenar,V. and Feigon,J. ( (1999) ) Localization of ammonium ions in the minor groove of DNA duplexes in solution and the origin of DNA A-tract bending. J. Mol. Biol., , 286, , 651–660.[CrossRef][Web of Science][Medline]

14. Denisov,V.P. and Halle,B. ( (2000) ) Sequence-specific binding of counterions to B-DNA. Proc. Natl Acad. Sci. USA, , 97, , 629–633.[Abstract/Free Full Text]

15. Cheatham,T.E.,III and Young,M.A. ( (2001) ) Molecular dynamics simulation of nucleic acids: successes, limitations, and promise. Biopolymers, , 56, , 232–256.

16. Young,M.A., Jarayam,B. and Beveridge,D.L. ( (1997) ) Intrusion of counterions into the spine of hydration in the minor groove of B-DNA: fractional occupancy of electronegative pockets. J. Am. Chem. Soc., , 119, , 59–69.

17. McConnell,K.J. and Beveridge,D.L. ( (2000) ) DNA structure: what's in charge? J. Mol. Biol., , 304, , 803–820.[CrossRef][Web of Science][Medline]

18. Hamelberg,D., McFail-Isom,L., Williams,L.D. and Wilson,W.D. ( (2000) ) Flexible structure of DNA: ion dependence of minor-groove structure and dynamics. J. Am. Chem. Soc., , 122, , 10513–10520.[CrossRef]

19. Mocci,F. and Saba,G. ( (2003) ) Molecular dynamics simulations of A·T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA. Biopolymers, , 68, , 471–485.[CrossRef][Web of Science][Medline]

20. Feig,M. and Pettitt,B.M. ( (1999) ) Sodium and chlorine ions as part of the DNA solvation shell. Biophys. J., , 77, , 1769–1781.[Web of Science][Medline]

21. Korolev,N., Lyubartsev,A.P., Laaksonen,A. and Nordenskiold,L. ( (2002) ) On the competition between water, sodium ions, and spermine in binding to DNA: a molecular dynamics computer simulation study. Biophys. J., , 82, , 2860–2875.[Web of Science][Medline]

22. Auffinger,P. and Westhof,E. ( (2000) ) Water and ion binding around RNA and DNA (C,G) oligomers. J. Mol. Biol., , 300, , 1113–1131.[CrossRef][Web of Science][Medline]

23. Auffinger,P. and Westhof,E. ( (2001) ) Water and ion binding around r(UpA)12 and d(TpA)12 oligomers. Comparison with RNA and DNA (CpG)12 duplexes. J. Mol. Biol., , 305, , 1057–1072.[CrossRef][Web of Science][Medline]

24. Reddy,S.Y., Leclerc,F. and Karplus,M. ( (2003) ) DNA polymorphism: a comparison of force fields for nucleic acids. Biophys. J., , 84, , 1421–1449.[Web of Science][Medline]

25. Várnai,P. and Lavery,R. ( (2002) ) Base flipping in DNA: pathways and energetics studied with molecular dynamics simulations. J. Am. Chem. Soc., , 124, , 7272–7273.[CrossRef][Web of Science][Medline]

26. Case,D.A., Pearlman,D.A., Caldwell,J.W., Cheatham,T.E.,III, Ross,W.S., Simmerling,C.L., Darden,T.A., Metz,K.M., Stanton,R.V., Cheng,A.L. et al. ( (1999) ) AMBER 6. University of California, San Francisco, CA.

27. Wang,J., Cieplak,P. and Kollman,P.A. ( (2000) ) How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? J. Comput. Chem., , 21, , 1049–1074.[CrossRef][Web of Science]

28. Patra,M. and Karttunen,M. ( (2004) ) Systematic comparison of force fields for microscopic simulations of NaCl in aqueous solutions: diffusion, free energy of hydration, and structural properties. J. Comput. Chem., , 25, , 678–689.[CrossRef][Web of Science][Medline]

29. Cheatham,T.E.,III ( (2004) ) Simulation and modeling of nucleic acid structure, dynamics and interactions. Curr. Opin. Struct. Biol., , 14, , 360–367.[CrossRef][Web of Science][Medline]

30. Kastenholtz,M.A. and Hünenberger,H.P. ( (2004) ) Influence of artificial periodicity and ionic strength in MD simulations of charged biomolecules employing lattice-sum methods. J. Phys. Chem. B, , 108, , 774–788.[CrossRef]

31. Berendsen,H.J.C., Postma,J.P.M., van Gunsteren,W.F., DiNola,A. and Haak,J.R. ( (1984) ) Molecular dynamics with coupling to an external bath. J. Chem. Phys., , 81, , 3684–3690.[CrossRef]

32. Ryckaert,J.P., Ciccotti,G. and Berendsen,H.J.C. ( (1977) ) Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comp. Phys., , 23, , 327–341.[CrossRef]

33. Darden,T., York,D. and Pedersen,L. ( (1993) ) Particle mesh ewald: an N·log(N) methods for Ewald sums in large systems. J. Chem. Phys., , 98, , 10089–10092.[CrossRef]

34. Brunne,R.M., Liepinsh,E., Otting,G., Wuthrich,K. and van Gunsteren,W.F. ( (1993) ) Hydration of proteins. A comparison of experimental residence times of water molecules solvating the bovine pancreatic trypsin inhibitor with theoretical model calculations. J. Mol. Biol., , 231, , 1040–1048.[CrossRef][Web of Science][Medline]

35. Lavery,R. and Sklenar,H. ( (1989) ) Defining the structure of irregular nucleic acids: convention and principles. J. Biomol. Struct. Dyn., , 6, , 655–667.[Web of Science][Medline]

36. Stofer,E. and Lavery,R. ( (1994) ) Measuring the geometry of DNA grooves. Biopolymers, , 34, , 337–346.[CrossRef][Web of Science][Medline]

37. Sharp,K.A. and Honig,B. ( (1990) ) Electrostatic interactions in macromolecules: theory and applications. Annu. Rev. Biophys. Chem., , 19, , 301–332.[CrossRef][Web of Science][Medline]

38. Makarov,V.A., Feig,M., Andrews,B.K. and Pettitt,B.M. ( (1998) ) Diffusion of solvent around biomolecular solutes: a molecular dynamics simulation study. Biophys. J., , 75, , 150–158.[Web of Science][Medline]

39. Atkins,P.W. ( (1998) ) Physical Chemistry. Oxford University Press, Oxford.

40. Koneshan,S., Rasaiah,J.C., Lynden-Bell,R.M. and Lee,S.H. ( (1998) ) Solvent structure, dynamics, and ion mobility in aqueous solutions at 25°C. J. Phys. Chem. B, , 102, , 4193–4204.[CrossRef]

41. Kennard,O., Isaacs,N.W., Motherwell,W.D., Coppola,J.C., Wampler,D.L., Larson,A.C. and Watson,D.G. ( (1971) ) The crystal and molecular structure of adenosine triphosphate. Proc. R. Soc. Lond. A, , 325, , 401–436.[Abstract/Free Full Text]

42. Reblova,K., Spackova,N., Sponer,J.E., Koca,J. and Sponer,J. ( (2003) ) Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets. Nucleic Acids Res., , 31, , 6942–6952.[Abstract/Free Full Text]

43. Reblova,K., Spackova,N., Stefl,R., Csaszar,K., Koca,J., Leontis,N.B. and Sponer,J. ( (2003) ) Non-Watson–Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E. Biophys. J., , 84, , 3564–3582.[Web of Science][Medline]

44. Stefl,R., Cheatham,T.E.,III, Spackova,N., Fadrna,E., Berger,I., Koca,J. and Sponer,J. ( (2003) ) Formation pathways of a guanine-quadruplex DNA revealed by molecular dynamics and thermodynamic analysis of the substates. Biophys. J., , 85, , 1787–1804.[Web of Science][Medline]

45. Spackova,N., Berger,I. and Sponer,J. ( (2001) ) Structural dynamics and cation interactions of DNA quadruplex molecules containing mixed guanine/cytosine quartets revealed by large-scale MD simulations. J. Am. Chem. Soc., , 123, , 3295–3307.[CrossRef][Web of Science][Medline]

46. Kosikov,K.M., Gorin,A.A., Lu,X.-J., Olson,W.K. and Manning,G.S. ( (2002) ) Bending of DNA by asymmetric charge neutralization: all-atom energy simulations. J. Am. Chem. Soc., , 124, , 4838–4847.[CrossRef][Web of Science][Medline]

47. Hud,N.V. and Plavec,J. ( (2003) ) A unified model for the origin of DNA sequence-directed curvature. Biopolymers, , 69, , 144–159.[CrossRef][Web of Science][Medline]

48. Pearlman,D.A. and Kollman,P.A. ( (1991) ) Are time-averaged restraints necessary for nuclear magnetic resonance refinement? A model study for DNA. J. Mol. Biol., , 220, , 457–479.[CrossRef][Web of Science][Medline]

49. Schmitz,U., Ulyanov,N.B., Kumar,A. and James,T.L. ( (1993) ) Molecular dynamics with weighted time-averaged restraints for a DNA octamer. Dynamic interpretation of nuclear magnetic resonance data. J. Mol. Biol., , 234, , 373–389.[CrossRef][Web of Science][Medline]

50. Ojha,R.P., Dhingra,M.M., Sarma,M.H., Shibata,M., Farrar,M., Turner,C.J. and Sarma,R.H. ( (1999) ) DNA bending and sequence-dependent backbone conformation NMR and computer experiments. Eur. J. Biochem., , 265, , 35–53.[Web of Science][Medline]

CiteULike    Connotea    Del.icio.us    What's this?

This article has been cited by other articles:

				R. Lavery, K. Zakrzewska, D. Beveridge, T. C. Bishop, D. A. Case, T. Cheatham III, S. Dixit, B. Jayaram, F. Lankas, C. Laughton, et al. A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA Nucleic Acids Res., November 6, 2009; (2009) gkp834v2. [Abstract] [Full Text] [PDF]

				S. Yamasaki, T. Terada, K. Shimizu, H. Kono, and A. Sarai A generalized conformational energy function of DNA derived from molecular dynamics simulations Nucleic Acids Res., November 1, 2009; 37(20): e135 - e135. [Abstract] [Full Text] [PDF]

				J. Curuksu, M. Zacharias, R. Lavery, and K. Zakrzewska Local and global effects of strong DNA bending induced during molecular dynamics simulations Nucleic Acids Res., June 1, 2009; 37(11): 3766 - 3773. [Abstract] [Full Text] [PDF]

				S. Jiranusornkul and C. A Laughton Destabilization of DNA duplexes by oxidative damage at guanine: implications for lesion recognition and repair J R Soc Interface, December 6, 2008; 5(Suppl_3): 191 - 198. [Abstract] [Full Text] [PDF]

				S. Khalid, P. J Bond, J. Holyoake, R. W Hawtin, and M. S.P Sansom DNA and lipid bilayers: self-assembly and insertion J R Soc Interface, December 6, 2008; 5(Suppl_3): 241 - 250. [Abstract] [Full Text] [PDF]

				C. Mura and J. A. McCammon Molecular dynamics of a {kappa}B DNA element: base flipping via cross-strand intercalative stacking in a microsecond-scale simulation Nucleic Acids Res., September 1, 2008; 36(15): 4941 - 4955. [Abstract] [Full Text] [PDF]

				A. Perez, F. Lankas, F. J. Luque, and M. Orozco Towards a molecular dynamics consensus view of B-DNA flexibility Nucleic Acids Res., April 1, 2008; 36(7): 2379 - 2394. [Abstract] [Full Text] [PDF]

				Y. Timsit and S. Bombard The 1.3 A resolution structure of the RNA tridecamer r(GCGUUUGAAACGC): Metal ion binding correlates with base unstacking and groove contraction RNA, December 1, 2007; 13(12): 2098 - 2107. [Abstract] [Full Text] [PDF]

				S. Fujii, H. Kono, S. Takenaka, N. Go, and A. Sarai Sequence-dependent DNA deformability studied using molecular dynamics simulations Nucleic Acids Res., September 25, 2007; 35(18): 6063 - 6074. [Abstract] [Full Text] [PDF]

				J. Sefcikova, M. V. Krasovska, J. Sponer, and N. G. Walter The genomic HDV ribozyme utilizes a previously unnoticed U-turn motif to accomplish fast site-specific catalysis Nucleic Acids Res., March 19, 2007; 35(6): 1933 - 1946. [Abstract] [Full Text] [PDF]

				N. Spackova and J. Sponer Molecular dynamics simulations of sarcin-ricin rRNA motif Nucleic Acids Res., February 2, 2006; 34(2): 697 - 708. [Abstract] [Full Text] [PDF]

				Y. Cheng, N. Korolev, and L. Nordenskiold Similarities and differences in interaction of K+ and Na+ with condensed ordered DNA. A molecular dynamics computer simulation study Nucleic Acids Res., January 30, 2006; 34(2): 686 - 696. [Abstract] [Full Text] [PDF]

				J. Zhou, S. Krueger, and S. K. Gregurick A coarse graining approach to determine nucleic acid structures from small angle neutron scattering profiles in solution Nucleic Acids Res., November 10, 2005; 33(19): 6361 - 6371. [Abstract] [Full Text] [PDF]

				K. Range, E. Mayaan, L. J. Maher III, and D. M. York The contribution of phosphate-phosphate repulsions to the free energy of DNA bending Nucleic Acids Res., March 1, 2005; 33(4): 1257 - 1268. [Abstract] [Full Text] [PDF]

				S. Y. Ponomarev, K. M. Thayer, and D. L. Beveridge Ion motions in molecular dynamics simulations on DNA PNAS, October 12, 2004; 101(41): 14771 - 14775. [Abstract] [Full Text] [PDF]

This Article Abstract Print PDF (672K) Supplementary Material Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Add to My Personal Archive Download to citation manager Search for citing articles in: ISI Web of Science (67) Request Permissions Commercial Re-use Guidelines for Open Access NAR Content Google Scholar Articles by Várnai, P. Articles by Zakrzewska, K. Search for Related Content PubMed PubMed Citation Articles by Várnai, P. Articles by Zakrzewska, K. Social Bookmarking          What's this?

Online ISSN 1362-4962 - Print ISSN 0305-1048

Copyright © 2009 Oxford University Press

Oxford Journals Oxford University Press

• Site Map
• Privacy Policy
• Frequently Asked Questions

Other Oxford University Press sites: Oxford University Press Oxford Journals China Oxford Journals Japan American National Biography Booksellers' Information Service Children's Fiction and Poetry Children's Reference Corporate & Special Sales Dictionaries Dictionary of National Biography Digital Reference English Language Teaching Higher Education Textbooks Humanities International Education Unit Law Medicine Music Online Products Oxford English Dictionary Reference Rights and Permissions Science School Books Social Sciences Very Short Introductions World's Classics