| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
* Laboratorium für Physikalische Chemie, ETH Hönggerberg, HCI, Zürich, Switzerland; and Department of Biology, Massachusetts Institute of Technology, Cambridge, Massachusetts
Correspondence: Address reprint requests to P. H. Hünenberger, Tel.: 41-1-632-5503; E-mail: phil{at}igc.phys.chem.ethz.ch.
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONTHEORYRESULTSCONCLUSIONSREFERENCES |
|---|
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONTHEORYRESULTSCONCLUSIONSREFERENCES |
|---|
), numerous experimental studies have investigated the relative stability of this DNA form, the influence of environmental factors on this stability, and the possible role of Z-DNA conformations in biological processes (2). As a result of this research, the occurrence of Z-DNA in vivo is now widely recognized, and the major thermodynamic factors influencing its stability are well understood. Yet, the conformational pathway(s) followed during the interconversion between B- and Z-DNA remain experimentally unknown at the atomic level.
Under physiological (low-salt) conditions, DNA preferentially adopts B or A conformations (along with 
20 sequence-dependent variations of these two forms (3–5)), characterized by Watson-Crick basepairing, anti orientation of the nucleotide bases, and a right-handed helicity (3,6). The A-DNA conformation is the one favored under low hydration or moderately high ionic-strength conditions, while the B-DNA conformation is preferred otherwise.
However, already in 1972, Pohl and Jovin (7) were able to observe an almost complete inversion of the ultraviolet circular dichroism (CD) spectrum of poly(GC) double-helical oligonucleotides upon addition of large amounts of salt (e.g., 4 m NaCl). The explanation for this effect was provided by the subsequent determination of the corresponding three-dimensional structure through x-ray diffraction (1,8), resulting in a left-handed helicity with an alternating anti (pyrimidine) and syn (purine) orientation of the nucleotide bases together with Watson-Crick basepairing (1). Due to the particular zigzag pattern adopted by the sugar-phosphate backbone, this new DNA form was termed Z-DNA. Apart from the inverted helicity, the Z-DNA structure is more compact along the double-helical axis (compared to the B form), with an average of one additional basepair per turn (3). As a consequence, the successive negatively-charged phosphate groups along the backbone reside on average closer to each other. At low ionic strengths, the resulting increase in electrostatic repulsion significantly destabilizes the Z form relative to the B form, even if the stacking interactions are expected to be more favorable in the Z-DNA conformation (9). Upon increasing the ionic strength, however, counterions provide more efficient electrostatic screening, and the relative stabilities of the two forms may be inverted (10). Indeed, it is now widely accepted that the dominant contribution to the free energy difference between the two forms is due to electrostatic interactions (11–15). The structural difference between the B- and Z-DNA conformations also leads to a different stacking pattern of the nucleotide bases. While the B-DNA conformation permits a continuous 
-stacking of successive bases, the stacking in Z-DNA involves discontinuous series of four-base clusters along each of the two strands (16).
Besides the ionic strength, the relative stabilities of the B- and Z-DNA forms are influenced by the nature of the (local) nucleotide sequence and the temperature. The Z-DNA conformation appears to be favored by sequences of repeating AT or CG doublets along one strand or in supercoiled DNA (17). Furthermore, due to its lower conformational entropy (compared to B-DNA), the relative stability of the Z-DNA conformation generally decreases upon raising the temperature (18). This property, together with the intrinsic differences in -stacking patterns between the two DNA forms (9), has recently been exploited in the design of a nanothermometer (19). The design of a nanomechanical motor device based on the B-to-Z transition has also been proposed (20,21) (under the basic assumption that the DNA unwinds during the transition), as well as that of a temperature-dependent reversible biomolecular switch (22).
Since the discovery of Z-DNA, much attention has been focused on the detection of Z-DNA in vivo and on the possible biological function of this particular DNA morphology (2,23–26). Because Z-DNA is intrinsically unstable under physiological conditions, this research is inherently difficult. However, it is now recognized that Z-DNA is not just a chemical oddity, but does occur in vivo (27). It is, for instance, formed transiently behind a moving RNA polymerase, stabilized by the negative supercoiling in the wake of the enzyme (28–31). In addition, there is growing evidence that Z-DNA segments or junctions may play a role in gene expression (16), recombination (32–34) and regulation (35–37), or in destructive effects of terminal differentiation (38), and that Z-DNA stability may be enhanced by etiological factors related with Alzheimer's disease (39,40).
Despite all these advances, one of the main unanswered questions remains how the very large conformational change associated with the B-to-Z transition actually takes place at the atomic level (41). Apart from the inversion of the helicity, the most striking conformational change occurring during the transition is the 180° rotation of the nucleotide bases with respect to their initial positions. The rotational mechanism involved during the transition appears to differ for the purine and pyrimidine bases. In the B-DNA conformation, the nucleotide base 
-angle of both nucleotide types, defined as the torsional angle O4'-C1'-N1-C2 (pyrimidine) or O4'-C1'-N9-C4 (purine), is in anti orientation, i.e., in the range 180 ± 90°. This orientation is required for the formation of canonical Watson-Crick hydrogen bonds, and bases with a syn orientation (0 ± 90°) in B-DNA serve as an indicator for a distortion of the double helix due to basepair opening or the presence of mismatched basepairs (3,42). However, in Z-DNA, purine bases generally adopt a syn conformation, i.e., they must rotate with respect to their associated deoxyribose unit along the B-to-Z transition pathway. In contrast, pyrimidine bases remain in anti orientation within Z-DNA due to steric constraints (3,42,43). Therefore, to keep the Watson-Crick basepairing intact, the whole pyrimidine nucleotide (including the sugar) needs to rotate, which places the deoxyribose 5'-carbon atom before the 3'-oxygen in terms of projection along the helix axis (with the 3'-end of the strand on the left and the 5'-end on the right), introducing a local half-twist in the backbone (44). This feature is commonly referred to as the chain-sense paradox, because the chain sense (3' 
5' or 5' 3') becomes locally inverted (44).
Four main hypotheses have been formulated concerning the mechanism of the B-to-Z transition.
The first hypothesis (known as the Wang model) suggests that the transition involves basepair opening before the rotation of the nucleotide bases and of the backbone dihedral angles (1,7,10,45,46).
The second hypothesis (known as the Harvey model) postulates that the transition consists of a series of nondisruptive transfiguration steps (correlated rotations of the base and backbone dihedral angles) permitting a transition without disruption of the Watson-Crick hydrogen bonds (44,47–49).
A third hypothesis assumes a transition through an intermediate (A-type) conformation without disruption of Watson-Crick hydrogen bonds and without sterical interference (50,51).
The fourth hypothesis was formulated recently (52), based on the first crystallographic structure of a B-Z DNA junction (cocrystallized with an enzyme). In this structure, two basepairs are flipped out of the duplex between the two DNA forms (analogous to B-DNA structures cocrystallized with DNA editing enzymes such as methyltransferases (53)). Based on this observation, the authors suggested that the B-to-Z transition may be accomplished by sequentially extruding basepairs from the duplex, permitting their free rotation by 180°, before reentrance at their original positions (in contrast to the Wang model, where base rotation is assumed to occur within the double-helix core).
Since the B-to-Z transition is known to be highly cooperative, it is not unreasonable to assume the involvement of domain formation during the interconversion (54), especially when long sequences are considered (e.g., negative supercoiling effects in plasmids). The transition pathway(s) may therefore also depend on the length and nature of the DNA sequence, as well as on the detailed experimental conditions. This view is in principle contained in the zipper model (55), where after an initial high-barrier nucleation step resulting in a short Z-DNA segment enclosed between two B-Z junctions, the junctions migrate in opposing directions as the intervening basepairs flip and the backbone winds in the left-handed direction to form Z-DNA. The zipper model is in principle compatible with any of the four above-mentioned hypotheses for the actual base-flipping mechanism.
To determine or rule out possible mechanisms, it is necessary to acquire time-resolved information at the atomic level, either by experiment or through computer simulations, which is in both cases highly nontrivial. On the one hand, experimental techniques typically used to probe the structure of (bio-)molecular systems at atomic resolution, such as x-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy, offer limited information in the context of dynamical structural changes. On the other hand, computer simulation methods, such as explicit-solvent molecular dynamics (MD) simulations based on empirical force fields, are currently limited to timescales of the order of 10–100 ns for oligonucleotides (56). This is clearly insufficient to investigate the B-to-Z transition by straightforward MD, since this transition occurs on a timescale of seconds to minutes (either induced by the addition of salt (7,10) or of a Z-DNA stabilizing protein (57)), with energy barriers up to hundreds of kJ x mol–1 (7,46,57,58). Furthermore, computational methods commonly used to improve the sampling along a transition pathway (e.g., umbrella sampling (59,60)) are not applicable here, because the reaction coordinate for the transition is unknown. Traditional free energy perturbation methods (61–64) are also out of question due to the large structural dissimilarity between the two conformations, so that only the effect of base modifications on the relative stability of Z-DNA (compared to B-DNA) can be investigated (65,66). Given these difficulties together with the inherent structural complexity of the transition, theoretical studies have so far mostly relied on simplified models, which commonly use a coarse-grained representation of the oligomer and model the solvent via a mean-field approach. Examples of such models are the Soumpasis potential-of-mean-force approach (67) and others (12,14,68–74), which are able to reproduce reasonably well the experimental free energy difference and its dependency on the charge, size, and concentration of the counterions (73). However, simplified representations of this kind cannot give detailed insight on the transition pathway(s) at the atomic level.
Recently, simulations at atomistic resolution were applied using the approach of the stochastic difference equation with length of Elber et al. (75) to investigate the B-to-Z transition of the 3'-(CGCGCG)-5' duplex. The results suggested yet another model for the B-to-Z transition pathway (76,77). In these simulations, the transition was observed to occur through an unwinding and overstretching of the oligomer, where the stretched conformation (increased terminal separation) could be viewed as an explicit intermediate structure for the transition (stretched-intermediate model). The main strength of the proposed model is that it provides a simple solution for the need to overcome very large sterical constraints to achieve the rotation of the nucleotide bases (purines) or sugar units (pyrimidines) during the transition. Conceptually, this suggested mechanism is related to the previously postulated mode of longitudinal breathing (44), which assumes that local backbone (longitudinal) stretching facilitates the flipping mechanism. The authors also proposed an experimental verification of the suggested mechanism by means of a fluorescence resonance energy transfer experiment (78) through the addition of two dyes at the oligomer termini. A decrease in the fluorescence signal during the transition would indicate an increase in the average dye separation and thus, provide evidence for the proposed stretching mechanism. However, three questions remain open with respect to the applied computational methodology and simulation setup. First, despite the conceptual elegance of variational methods based on minimizing the classical action integral, it is not clear whether the observed transition pathway has any physical significance at room temperature (since it corresponds in principle to a zero-temperature pathway). Second, the approach does not treat solvent molecules explicitly, but through a generalized-Born implicit (mean-field) solvation model (79). While the performance of the generalized Born approach has been excellent for some systems, in particular for studies on peptide folding (80) (with some notable exceptions (81,82)), its validity remains questionable in the context of DNA simulations, due to the key role of specific solute-solvent interactions (83–90). Third, the above simulations were undertaken in the absence of explicit counterions, and did not permit the investigation of ionic-strength effects on the transition.
Any computational approach to investigate the B-to-Z transition pathway currently requires the application of special methods for bridging of the timescale gap between simulation and experiment. However, this should be done by imposing as little constraint as possible on the dynamics of the system. A simple approach to this purpose is targeted molecular dynamics (TMD) (91–94). This algorithm allows for the description of a progressive conformational change between a starting and a target structure, by linearly decreasing the (mass-weighted) root-mean-square deviation (RMSD) in atomic position between the current structure and the target one (RMSD constraint) along the simulation. By performing the transformation reasonably slowly, the system should ideally remain at equilibrium with respect to all degrees of freedom orthogonal to the RMSD constraint. Thus, the sampled conformational pathway should provide a reasonable first guess to the true transition pathway, under the assumption that this true pathway involves a monotonic decrease of the RMSD to the target structure. In addition, the free energy difference between the initial and final states may be estimated by integrating the constraint force along the forced transition pathway (59), in a way analogous to the slow-growth protocol (59,61,95). Finally, this approach is easily combined with an explicit representation of the solvent molecules and counterions, allowing for a realistic modeling of solvation and ionic-strength effects.
The main goal of the present work is to:
,77), suggesting a stretched-intermediate model for this transition.
To this purpose, MD and TMD simulations are carried out for the double-helical hexanucleotide 3'-(CGCGCG)-5' in water. The self-complementary (CGCGCG)2 duplex (melting temperature: Tm 
320 K (96)) has been investigated previously in a number of experimental and theoretical studies (1,97–102), and is known to show a high propensity to form Z-DNA due to the presence of CG repeats (103). Six independent MD simulations at room temperature involve the hexamer duplex starting from either the B or the Z conformation, and either without counterions, with a neutralizing amount of sodium counterions, or with an excess amount (4 m) of sodium and chloride counterions. In addition, three independent TMD simulations of the B-to-Z transition (starting from the B conformation) are performed under the same temperature and ionic-strength conditions.
|
THEORY |
|---|
|
TOPABSTRACTINTRODUCTIONTHEORYRESULTSCONCLUSIONSREFERENCES |
|---|
,105) together with the GROMOS96 45A4 force field for DNA (106) and the SPC water model (107). For the six MD simulations, a modeled B-DNA or Z-DNA hexamer duplex with sequence 3'-(CGCGCG)-5' and net charge –10e was solvated in a cubic box of dimension 5 · 5 · 5 nm3 with explicit water molecules at a density of 1 g x cm–3, resulting in a total of 4111 (B-DNA) or 4102 (Z-DNA) solvent molecules. For each conformer, one system was simulated in the absence of counterions, one system with a neutralizing amount of 10 sodium ions and one system with an amount of sodium and chloride counterions corresponding to a 4-molal solution (for an overall neutral system). The initial ion locations were selected by successively replacing water molecules with sodium counterions at the positions of lowest electrostatic potential or with chloride ions at the positions of highest electrostatic potential (sodium and chloride ions were inserted in alternation), while obeying an ion-ion exclusion distance of 0.35 nm. Since the modeled B- and Z-DNA structures do not exactly correspond to stable energy minima in the GROMOS96 force field (106), the starting and target structures for the three TMD simulations were chosen as the equilibrated structures of the corresponding MD simulations with a neutralizing amount of counterions after 1 ns. The simulation conditions and codes employed in this article are summarized in Table 1 for the nine simulations.
|
,109) for the omitted interactions beyond the long-range cutoff distance (the relative permittivity of the medium outside the cutoff sphere was set to 66.6 (110)). Although lattice-sum methods (e.g., Ewald, P3M, or PME) are often quoted as the proper remedy to the numerous problems associated with straight-cutoff truncation of electrostatic interactions (111–118), the reaction-field scheme also appears to provide excellent (and often qualitatively similar) results (106,115,119). Furthermore, the application of lattice-sum methods in explicit-solvent simulations of solutions may lead (at least in some cases) to periodicity-induced artifacts (120–126).
All simulations were carried out at constant volume and at a constant temperature of 280.15 K (40 K below the melting temperature of the hexamer (96) to avoid any possible denaturation), maintained by weak-coupling to a heat bath (127) using a coupling time of 0.1 ps. For equilibration, the MD simulations were slowly heated up over a period of 0.1 ns while applying a position-restraining potential with a force constant of 2.5 · 104 kJ x mol–1 · nm–2 on all solute atoms, so as to relax the surrounding solvent molecules. Bond lengths were constrained by application of the SHAKE procedure (128) using a relative geometric tolerance of 10–4. A timestep of 2 or 1 fs was used to integrate the equations of motion for the MD and TMD simulations, respectively. The system center-of-mass motion was removed every 5000 time steps to prevent the buildup of center of mass velocity. All MD and TMD simulations were carried out for 2 ns.
For the calculation of the RMSD constraint forces (implemented analogously to the SHAKE algorithm (128) via a Lagrange multiplier approach), the same relative geometric tolerance of 10–4 was used. Since the bond-length (SHAKE) and RMSD (TMD) constraints may interfere with each other during a given integration timestep, the two protocols were iterated sequentially until convergence was reached for both. The atoms considered for the RMSD constraint were the NT = 84 heavy atoms of the 12 nucleotide bases, which should capture the largest conformational change during the B-to-Z transition. At every integration timestep, the target structure was superimposed onto the current structure by means of a (mass-weighted) translational and rotational least-squares fit (129) involving the NT atoms, followed by a calculation of the (mass-weighted) RMSD value (R), as
 | (1) |
the corresponding coordinate vector in the target structure after fitting onto r, and 
, with mi being the mass of atom i and 
the average atomic mass of the NT atoms. The use of mass-weighted coordinates for the least-squares fitting and RMSD calculation results in a constraint force exempt of any overall translational and rotational component (94). The constraint to be satisfied by the coordinates at any timestep is written as
 | (2) |
, is added to the physical (unconstrained) potential energy Vuc, i.e.,
 | (3) |
on each atom reads
 | (4) |
denotes the common Lagrange multiplier determining the strength of the constraint force at the given timestep. In the case of the leap-frog integrator the coordinates are propagated in time as
 | (5) |
 | (6) |
 | (7) |
t is the integration timestep, 
, vi, and 
denote the free-flight coordinate, velocity, and unconstrained force associated with atom i and 
corresponds to the coordinate shift required to enforce the constraint at time t + t. An expression for the Lagrange multiplier (t) is obtained by inserting Eqs, 1, 4, 5, and 7 into Eq. 2 for r(t + t). This leads to the quadratic equation
 | (8) |
 | (9) |
 | (10) |
 | (11) |
 | (12) |
). The corresponding coordinate shift 
to be applied to each atom i may then be calculated at each timestep as
 | (13) |
Over the course of the TMD simulation, the RMSD distance Rref(t) to be satisfied is linearly reduced from an initial RMSD difference RBZ = 0.55 nm between the B- and Z-form down to zero, via a constant decrease Rref = 2.75 · 10–7 nm at each timestep t = 1 fs, which corresponds to a total simulation time of 2 ns.
The (Helmholtz) free-energy profile A(R) along the transition pathway defined by the RMSD reaction coordinate may be determined by application of the thermodynamic integration formula (130)
 | (14) |
...
R denotes ensemble averaging at a specified value R of the RMSD and the second equality follows from Eqs. 2 and 3. Using a slow-growth approximation (59), this profile may be evaluated numerically as (94)
 | (15) |
is chosen here to enforce A(RBZ) = 0, i.e., the zero of the profile is set to the B-DNA conformation.
Due to the use of a mass-weighted RMSD as reaction coordinate, the configurational probability distribution is not affected by metric-tensor effects (94,131,132), and Eq. 15 becomes exact in the limit of an infinitesimally slow variation of Rref. However, due to the use of a polar coordinate system, the free energy profile A(R) may appear biased toward high R-values because the Cartesian coordinate volume associated with a given value of R decreases with R (94,133). In particular, A(R) diverges to infinity when R tends toward zero (i.e., by the end of the simulation). This effect can be viewed as arising from an entropic contribution related to the restriction of the accessible Cartesian volume upon decreasing R.
The analysis of the simulations was performed in terms of the following basic properties:
In addition, the rotation of the nucleotide bases during the TMD simulations was monitored through the time-series of the -dihedral angle around the sugar-base bond, defined by the atoms O4-C1'-N1-C2 (cytidine) and O4'-C1'-N9-C4 (guanosine). Because the cytidine bases do not rotate around the -dihedral angle during the transition (see Introduction), it is useful to monitor an alternate internal coordinate that is also sensitive to the rotation of the whole nucleotide. This was done by following the time-series of the angle 
between the plane of the nucleotide base in the current configuration and the corresponding plane in the starting structure (after an all-atom non-mass-weighted translational and rotational least-squares fit of the two structures).
The ultraviolet (UV) absorbance (and in particular CD spectra in the UV range) are known to serve as a sensitive indicator of DNA conformation (134). Nucleic acids exhibit a strong * transition in the far UV, due to the aromaticity of the bases, as well as an n * transition at longer wavelengths. The CD spectra of B- and Z-DNA show distinct differences, the hallmark for Z-DNA being an intense negative band at 195 nm (7). It is possible to estimate the rotational-strength R10 of a 0 1 transition based on a specific system configuration from ab initio calculations, as the imaginary part of the scalar product of the corresponding electronic µ and magnetic m transition dipole moments R10 = Im{µ01 · m10}. However, this calculation requires knowledge of the wavefunction for the ground- and excited states, which is computationally intractable for the present system (292 atoms). An alternative approach is to describe the ground and excited states at various wavelengths by monopole sets, and to construct a Hamiltonian matrix from which eigenvalues, eigenfunctions, and ultimately µ01 and m10 can be estimated (135,136). Such monopole sets for the two types of nucleotide bases considered here were derived by fitting to the direction and magnitudes of experimentally obtained transition dipole moments as described earlier (134). This approach was applied to calculate far-UV CD spectra based on selected configurations sampled during the present TMD simulations. The molar CD signal at a given wavelength , E(), represents the difference in absorption between left- and right-circularly polarized light, reported in units of M–1xcm–1. It may be calculated from the rotational-strength as 
, where 10 is the half-bandwidth for the given 0 1 transition and R10 is expressed in the same units as E(). In any monopole approximation, charge transfer effects, which can be of importance in the longer wavelength regions of DNA spectra, cannot be accounted for. In the far-UV region below 200 nm, however, these effects should be negligible.
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONTHEORYRESULTSCONCLUSIONSREFERENCES |
|---|
0.5 ns. For the simulations starting from Z-DNA, the system without counterions is clearly unstable and loses its double-helical structure almost immediately, reaching an RMSD value of 0.50 nm after 2 ns. Due to the large repulsive electrostatic interactions within the Z-DNA backbone, this behavior is not surprising. In contrast, the two Z-DNA systems with explicit counterions are stable, and show maximum deviations of 0.25–0.30 nm.
|
,137,138).
The analysis of the interstrand hydrogen-bonding patterns during these six MD simulations is reported in Fig. 2, a–f, and Table 2. For the sequence considered, a total of 18 canonical Watson-Crick (WC) hydrogen-bond pairs are possible (Table 2). The simulations starting from B-DNA exhibit a high occurrence of canonical hydrogen bonds under all conditions, i.e., these simulations conserve the initial WC basepairing pattern to a high extent. For the central basepairs, these hydrogen bonds are typically present over 85–100% of the simulation time. Their occurrence may be reduced in the strand termini due to increased mobility (see RMSF in Fig. 1, a–c) and more effective competition with hydrogen bonds to water molecules. Note also that the overall occurrence of canonical hydrogen bonds slightly increases upon increasing the ionic strength of the solution. For the simulations starting from Z-DNA, the system without counterions is characterized by a progressive disruption of the canonical hydrogen-bonding pattern, confirming the result of the RMSD analysis (Fig. 1 d). The hydrogen-bonding pattern is better preserved for the system with a neutralizing set of counterions, and essentially stable for the simulation with an excess of counterions. The increase in ionic strength leads to a decrease in the repulsive electrostatic interactions between backbone phosphate groups, which in turn results in a stabilization of the Z-DNA conformation, in agreement with independent experimental (7,10) and theoretical (73) results.
|
|
|
0.7 ns, i.e., the simulations are progressing more rapidly away from the B-DNA form than they are moving toward the Z-DNA form. Due to the large conformational change involved in these simulations, the corresponding RMSF values averaged over the course of the transition are very high (because they account not only for the structural fluctuations but also for the overall structural change).
The analysis of the interstrand hydrogen-bonding patterns during these three TMD simulations is reported in Fig. 2, g–i, and Table 2. A major disruption of WC hydrogen bonds occurs after 1 ns upon leaving the initial B-DNA conformation. However, after an additional time period of 0.5 ns, most hydrogen bonds are recovered and maintained up to reaching the target Z-DNA conformation. Note that this transient disruption of WC hydrogen bonds is not accompanied by major distortions of the basepairing geometries (in particular, the flipping-out of bases into solution reported in the crystallographic structure of a B-Z junction (52), is never observed in the present simulations). There is no obvious difference with regard to the time evolution of the hydrogen-bonding pattern under different ionic-strength conditions.
The time-series of the total number of oligomer to solvent hydrogen bonds is shown in Fig. 3, g–i, for the three TMD simulations. These time-series show a significant increase between 1.0 and 1.5 ns, coinciding with the transient disruption of interstrand WC hydrogen bonds (Fig. 2, g–i). This suggests that the energetic cost involved in the transient breakdown of intrasolute hydrogen bonds is partly compensated by a concomitant increase in solute-solvent hydrogen-bonding. The total number of hydrogen bonds are characterized by mean values and standard deviations (excluding a total of 0.2 ns equilibration time) of 114.7 ± 6.6 (B0), 117.2 ± 7.5 (BZN), and 120.91 ± 6.1 (BZE). These average values are systematically higher than values in the corresponding MD simulations starting from either B- or Z-DNA. However, in contrast to the observations made in the context of these MD simulations, the number of solute-solvent hydrogen bonds is larger at higher ionic strength.
The transition pathway emerging from the TMD simulation with a neutralizing amount of counterions is illustrated in Fig. 4. Along the transition, the oligomer starts to partially unwind (first 0.25 ns), then undergoes a slight bending (0.25–0.5 ns) until reaching an extended or stretched conformation after 1 ns. The end-to-end distance of this stretched conformation at 1 ns (defined here as the average of the 1C-6G and 7C-12G phosphate atom distances of both strands) is 2.39 nm, to be compared with the corresponding values of 2.09 nm and 2.12 nm for the equilibrated B- and Z-DNA structures, respectively, corresponding to a stretching by 13.5%. At this time point, the nucleotide bases start to lose their canonical interstrand hydrogen bonds (Fig. 2, g–i) at the benefit of solute-solvent hydrogen bonding (Fig. 3, g–i), and begin to flip. After 1.5 ns, the WC hydrogen-bond pairing is entirely restored and the final Z-DNA conformation is progressively reached. The two other simulations (without and with an excess of counterions) follow a very similar pathway (not shown). The passage through a stretched configuration seems natural, because it facilitates the base-flipping process by allowing the nucleotide bases to occupy more space and enables the formation of the (local) chain-sense paradox. The transition pathway observed in the present explicit-solvent TMD simulations is in good qualitative agreement with the recent findings of Lim et al. (76,77), based on the simulation of a least-action path using an implicit solvent model.
|
). Nevertheless, the calculated single-structure spectra can be interpreted in a qualitative way. The individual structures along the pathway followed by the TMD simulation give rise to a spectral evolution, which is in good agreement with the experimentally observed signal inversion at 195 nm (7), i.e., the emergence of the spectral trademark of the Z-DNA configuration. The continuous progression of the band inversion over the course of the transition is also compatible with the experimental evolution observed in time-resolved measurements (7).
|
-angle by 180°, changing orientation from anti to syn, the cytosine bases remain in anti orientation and the flipping is accomplished through a rotation of the whole nucleoside involving the angle . The corresponding - and -angle values of the starting (equilibrated) B-DNA and target (equilibrated) Z-DNA conformations of the TMD simulations are reported in Table 3.
|
-angles are displayed in Fig. 6 for the three TMD simulations. For all simulations, the -angle values of the target Z-DNA structure (Table 3) are reached after 2 ns, i.e., all bases adopt the expected syn (guanine) or anti (cytosine) orientations with respect to the deoxyribose sugar moiety. The evolution of the -angle appears to be very similar in the three TMD simulations. The guanosine bases are distributed into two groups, i.e., those for which the rotation occurs gradually through an increase in (so-called high-anti path; 2G, 8G) and those for which rotation occurs in a stepwise fashion through a decrease in (4G, 6G, 10G, 12G). Although the cytosine bases paired with a given guanosine do not flip, all of them show significant variations in their -angle during the time in which their partner guanosine flips.
|
between the nucleotide-base plane in the current and starting (B-DNA) structure (after least-squares-fit superimposition) are displayed in Fig. 7 for the three TMD simulations. Due to possible ambiguities related to the fitting procedure, the results may only be interpreted in a qualitative fashion. The flipping of the cytosine bases, mediated by the rotation of the corresponding deoxyribose moiety, is evident here. Here also, the bases are roughly distributed into two groups, i.e., those showing gradual rotation (1C, 2G, 6G, 12G) and those flipping more suddenly at a later stage of the transition (3C, 4G, 5C, 7C, 8G, 9C, 10G, 11C).
|
1C in strand 1 and 7C 12G in strand 2, whereby bases 3C, 2G and 11C rotate the latest. In other words, the flipping process progresses in direction 5' 3' along strand 1, and 3' 5' along strand 2. This kind of propagation mechanism may be viewed as a (cooperative) zipperlike transition (55).
The free energy profiles A(R) calculated from the integral of the constraint force over the sampled pathways, are shown in Fig. 8 for the three TMD simulations under different ionic-strength conditions. The free energy profiles for the three simulations exhibit the same sigmoidal shape, with a transition point at R 0.25 nm (i.e., after 1 ns simulation time), which corresponds to the stretched intermediate conformation in Fig. 4. Only the simulation with an excess of counterions shows a clear minimum for low values of the reaction coordinate, suggesting again that the Z-DNA conformation is intrinsically unstable under lower ionic-strength conditions. The observation that the A(R) value for Z-DNA is higher compared to B-DNA, even at high salt concentration, can have two origins:
).
|
125–150 kJ x mol–1, which is of the same order of magnitude as reported literature values for the activation energy of this process (159 kJ x mol–1 (57)). |
CONCLUSIONS |
|---|
|
TOPABSTRACTINTRODUCTIONTHEORYRESULTSCONCLUSIONSREFERENCES |
|---|
,77) using a different computational methodology and an implicit-solvation model. The main features are the formation of a stretched (13.5%) intermediate configuration followed by a transient disruption of canonical hydrogen bonds, partly compensated energetically by a concomitant increase in the number of solute-solvent hydrogen bonds, and concluded by reformation of Watson-Crick hydrogen bonds and relaxation to the final Z-DNA structure. Two of the six guanine bases (2G, 8G) were found to undergo rotation about the deoxyribose -angle (to the syn orientation) via the high-anti pathway. The flipping process was found to progress along the 5' 3' direction in strand 1, and along the 3' 5' direction in strand 2, in accordance with a zipperlike progression mechanism (55). The stretched intermediate configuration was found to approximately coincide with the maximum of the corresponding free energy curve. The fundamental assumption underlying the realism of the mechanism suggested by the present TMD simulations is that the relaxation of the degrees of freedom orthogonal to the RMSD reaction coordinate (for which the TMD speed-up is of the order of 1010 compared to experiment) occurs on the subnanosecond timescale. Although there is currently no way to prove the validity of this assumption, the suggested pathway appears physically plausible, because it is compatible with the sterical constraints involved in fully flipping the bases (guanosine) and the deoxyribose-moieties (cytosine) during the transition, and provides a simple solution to the chain-sense paradox. This transition pathway also appears to be qualitatively insensitive to the ionic-strength conditions. However, the Z-DNA conformation only corresponds to a (shallow) free-energy minimum (along the profile) for the simulation under high ionic-strength conditions. Submitted on February 22, 2006; accepted for publication July 10, 2006.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONTHEORYRESULTSCONCLUSIONSREFERENCES |
|---|
2. Rich, A., and Z. Shuguang. 2003. Z-DNA: the long road to biological function. Nature Rev. Gen. 4:566–573.[Medline]
3. Saenger, W. 1984. Principles of Nucleic Acid Structure. Springer Verlag, New York.
4. Wells, R., and S. Harvey, editors. 1989. Unusual DNA Structures. Springer, New York.
5. Ghosh, A., and M. Bansal. 2003. A glossary of DNA structures from A to Z. Acta Crystallogr. D. 59:620–626.[CrossRef][Medline]
6. Watson, J. D., and F. H. C. Crick. 1953. Molecular structure of nucleic acids—a structure for deoxyribose nucleic acid. Nature. 171:737–738.[CrossRef][Medline]
7. Pohl, F. M., and T. M. Jovin. 1972. Salt-induced co-operative conformational change of a synthetic DNA: equilibrium and kinetic studies with poly(dG-dC). J. Mol. Biol. 67:375–396.[CrossRef][Medline]
8. Ho, P. S., and B. H. M. Mooers. 1997. Z-DNA crystallography. Biopolymers. 44:65–90.[CrossRef]
9. Sponer, J., H. A. Gabb, J. Leszczynski, and P. Hobza. 1997. Base-base and deoxyribose-base stacking interactions in B-DNA and Z-DNA: a quantum-chemical study. Biophys. J. 73:76–87.
10. Pohl, F. M. 1983. Salt-induced transition between two double-helical forms of oligo(DC-DG). Cold Spring Harb. Symp. Quant. Biol. 47:113–117.[Medline]
11. Kollman, P. A., P. K. Weiner, and A. Dearing. 1981. Theoretical studies of the structure and energies of basepaired nucleotides and the dissociation kinetics of a proflavine-dinucleotide complex. Ann. N. Y. Acad. Sci. 367:250–268.[Abstract]
12. Gueron, M., and J. P. Demaret. 1992. A simple explanation of the B-to-Z transition of DNA. Proc. Natl. Acad. Sci. USA. 89:5740–5743.
13. Misra, V. K., and B. Honig. 1996. The electrostatic contribution to the B-to-Z transition of DNA. Biochemistry. 35:1115–1124.[CrossRef][Medline]
14. Gueron, M., J. P. Demaret, and M. Filoche. 2000. A unified theory of the B-Z transition of DNA in high and low concentrations of multivalent ions. Biophys. J. 78:1070–1083.
15. Demaret, J. P., and M. Gueron. 2000. The sensitivity of the B-Z transition of DNA to multivalent cations in high and low concentrations: a unified electrostatic interpretation. J. Biomol. Struct. Dyn. S2:377–382.
16. Rich, A., A. Nordheim, and A. H. Wang. 1984. The chemistry and biology of left-handed Z-DNA. Annu. Rev. Biochem. 53:791–846.[CrossRef][Medline]
17. Haniford, D. B., and D. E. Pulleyblank. 1983. The in vivo occurrence of Z DNA. J. Biomol. Struct. Dyn. 1:593–609.[Medline]
18. Sugiyama, H., K. Kawai, A. Matsunaga, K. Fujimoto, I. Saito, H. Robinson, and A. H. Wang. 1996. Synthesis, structure and thermodynamic properties of 8-methylguanine-containing oligonucleotides: Z-DNA under physiological salt conditions. Nucleic Acids Res. 24:1272–1278.
19. Tashiro, R., and H. Sugiyama. 2003. A nanothermometer based on the different -stackings of B- and Z-DNA. Ang. Chem. Intl. Ed. 42:6018–6020.[CrossRef]
20. Mao, C., W. Sun, Z. Shen, and N. C. Seeman. 1999. A nanomechanical device based on the B-Z transition of DNA. Nature. 397:144–146.[CrossRef][Medline]
21. Seeman, N. C., and A. M. Belcher. 2002. Emulating biology: building nanostructures from the bottom up. Proc. Natl. Acad. Sci. USA. 99:6451–6455.
22. Tashiro, R., and H. Sugiyama. 2005. Biomolecule-based switching devices that respond inversely to thermal stimuli. J. Am. Chem. Soc. 127:2094–2097.[CrossRef][Medline]
23. Rich, A. 1996. The biology of left-handed Z-DNA. J. Biol. Chem. 271:11595–11598.
24. Stollar, B. D. 1997. Why the difference between B-DNA and Z-DNA? Lupus. 6:327–328.
25. Belmont, P., J. F. Constant, and M. Demeunynck. 2001. Nucleic acid conformation diversity: from structure to function and regulation. Chem. Soc. Rev. 30:70–81.[CrossRef]
26. Gagna, C. E., and W. C. Lambert. 2003. The halting arrival of left-handed Z-DNA. Med. Hypotheses. 60:418–423.[CrossRef][Medline]
27. Kim, Y. G., K. Lowenhaupt, D. Y. Oh, K. K. Kim, and A. Rich. 2004. Evidence that Vaccinia virulence factor ER binds to Z-DNA in vivo: implications for development of a therapy for poxvirus infection. Proc. Natl. Acad. Sci. USA. 101:1514–1518.
28. Liu, L. F., and J. C. Wang. 1987. Supercoiling of the DNA template during transcription. Proc. Natl. Acad. Sci. USA. 84:7024–7027.
29. Wittig, B., T. Dorbic, and A. Rich. 1991. Transcription is associated with Z-DNA formation in metabolically active permeabilized mammalian cell nuclei. Proc. Natl. Acad. Sci. USA. 88:2259–2263.
30. Wittig, B., S. Wolfl, T. Dorbic, W. Vahrson, and A. Rich. 1992. Transcription of human c-Myc in permeabilized nuclei is associated with formation of Z-DNA in three discrete regions of the gene. EMBO J. 11:4653–4663.[Medline]
(
