help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Originally published as Biophys J. BioFAST on November 11, 2005.
doi:10.1529/biophysj.105.059535
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.105.059535v1
90/3/1000    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Cubero, E.
Right arrow Articles by Orozco, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Cubero, E.
Right arrow Articles by Orozco, M.
Biophysical Journal 90:1000-1008 (2006)
© 2006 The Biophysical Society

Theoretical Study of the Hoogsteen–Watson-Crick Junctions in DNA

Elena Cubero *, F. Javier Luque {dagger} and Modesto Orozco * {ddagger}

* Molecular Modelling & Bioinformatic Unit, Institut de Recerca Biomèdica-Parc Científic de Barcelona, Barcelona 08028, Spain; and {dagger} Departament de Fisicoquímica, Facultat de Farmacia, and {ddagger} Departament de Bioquímica i Biologia Molecular, Facultat de Química, Universitat de Barcelona, Barcelona 08028, Spain

Correspondence: Address reprint requests to M. Orozco, E-mail: modesto{at}mmb.pcb.ub.es or F. J. Luque, E-mail: javier{at}far1.far.ub.es.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
A series of d (AT)n oligonucleotides containing mixtures of normal B-type Watson-Crick and antiparallel Hoogsteen helices have been studied using molecular dynamics simulation techniques to analyze the structural and thermodynamic impact of the junction between Watson-Crick and antiparallel Hoogsteen structures. Analysis of molecular dynamics simulations strongly suggests that for all oligonucleotides studied the antiparallel Hoogsteen appears as a reasonable conformation, only slightly less stable than the canonical B-type Watson-Crick one. The junctions between the Watson-Crick and Hoogsteen structures introduces a priori a sharp discontinuity in the helix, because the properties of each type of conformation are very well preserved in the corresponding fragments. However, and quite counterintuitively, junctions do not largely distort the duplex in structural, dynamics or energetic terms. Our results strongly support the possibility that small fragments of antiparallel Hoogsteen duplex might be embedded into large fragments of B-type Watson-Crick helices, making possible protein-DNA interactions that are specific of the antiparallel Hoogsteen conformation.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
DNA is very polymorphic and can adapt its structure depending on the sequence, the presence of ligands, or changes in the environment (1Go–3Go). Repetitive sequences have an especially rich conformational space (1Go–4Go) showing a marked tendency to appear in left-handed duplexes, triplexes, and tetraplexes (4Go–8Go), or even duplexes with parallel arrangements of the complementary strands (2Go,9Go,10Go). Particularly, d (AT) sequences have been experimentally detected in: i), partly or largely distorted B-type conformation (11Go–13Go), ii), in the right-handed C (14Go) and D (15Go) forms, and iii), even left-handed models have been suggested (16Go). Recently, Subirana and co-workers (17Go,18Go) solved crystals of the d (AT) polymer where the duplex appeared as an antiparallel right-handed (apH) double helix with Hoogsteen d (A-T) pairings and all the adenosines in the syn conformation. The structure has not been yet directly verified in solution by NMR experiments, but recent molecular dynamics (MD) simulations (19Go) suggested that the apH conformation is not much with respect to the canonical Watson-Crick B-form (B helix) in aqueous solution. In fact, MD and molecular mechanics-Poisson Boltzmann/surface area (MM-PB/SA) calculations suggested that both B and apH duplexes are almost isoenergetic and only the greater rigidity of the latter explain the larger population of B helix in physiological environments (19Go).

Hoogsteen d (A-T) pairings are known to stabilize triads in parallel triplexes (3Go,4Go,20Go,21Go), as well as in parallel duplexes (9Go). Furthermore, Hoogsteen A-T pairings are commonly found inserted in the middle of B-type Watson-Crick duplexes when the DNA interacts with some drugs or proteins (22Go–26Go). Because d (AT)n regions have a key role in the control of gene expression (24Go,27Go,28Go), the formation of Hoogsteen islands in the middle of tracks of d (AT)n B-DNA might be involved in the existence of a subtle regulatory mechanism of gene function (19Go,22Go,24Go).

In a previous article (19Go) we analyzed the properties of pure apH duplexes made with tracts from 2 to 8 d (AT), i.e., duplexes from 4 to 16 mer. These calculations can provide evidence for the existence of conformation derived from crystal structures reported by Subirana's group (17Go,18Go) in solution. Our calculations also suggested that segments of apH helix might exist in conditions where the flexibility of the DNA is externally restricted. Preliminary calculations indicated that, at least for poly-D (A) sequences, the junctions between B and apH helices are not sterically hindered. In this article, we will analyze systematically the structure of chimeric duplexes formed by portions of B and apH helices in d (AT) tracks. For this purpose we have performed molecular dynamics simulations of 12-, 14-, and 16-mer duplexes containing each 2, 4, or 6 apH steps in the center of the helix and surrounded by segments of normal B helix. Results are compared with those previously obtained for pure B and apH duplexes of the same global length. The analysis of >40 ns of equilibrated MD simulation on the same structural motif allowed us to describe with detail the characteristics of the B/apH chimeras.


   METHODS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
To analyze the structure and stability of chimeras containing fragments of apH and B helices for d (AT)n sequences, we built starting structures for duplexes d (AT)n (n = 12, 14, and 16) using standard Arnott's parameters (29Go) for B-DNA (see a summary of simulations in Table 1). In the center of each structure, apH helices with length 2, 4, or 6 were inserted using geometrical data from crystal structure (17Go). The junction steps were optimized by restricted MD to avoid unreliable starting geometries. The nine different starting structures were neutralized by adding a suitable number of Na+ ions and immersed in rectangular boxes containing between 3000 and 3600 water molecules. All of the systems were optimized, thermalized (298 K), and equilibrated for 200 ps using our standard multistage protocol (30Go–32Go). The equilibrated structures were then subjected to 2 and 5 ns (after the referee's suggestions some simulations were extended to 10 ns for control purposes; see Table 1) of MD simulation at constant temperature (298 K; after the referee's suggestions we performed three additional simulations at T = 310; see Table 1) and pressure (1 atm). Periodic boundary conditions were used to simulate a diluted environment, and the particle mesh Ewald method (33Go). SHAKE (34Go) was used to constrain all of the bonds at their equilibrium positions, which allowed us to use a 2-fs time step for integration of Newton equations. AMBER-98 (35Go,36Go) and TIP3P (37Go) force fields were used to describe DNA and water. All MD simulations were performed using the AMBER6.0 suite of programs (38Go). This simulation protocol is identical to that used in our previous simulations (19Go), which allowed us to include in our comparison the B and apH duplexes (n = 12, 14, and 16) described there.


View this table:
[in this window]
[in a new window]
  		TABLE 1  Scheme of MD simulations of pure helices (B and apH), and chimeras containing both types of helices

 
The free energy of every structure was estimated from Eq. 1 (9Go,19Go,39Go,40Go), where averages were performed for snapshots collected after the first nanosecond of the corresponding trajectory. The intramolecular term (Eintra) was computed using the standard AMBER-98 force field (35Go,36Go). The solvation contribution (Gsol) was determined using: i), finite difference Poisson-Boltzmann (PB) calculations as implemented in the MEAD program (41Go,42Go) for initial and final grid spacing of 1 and 0.4 Å, with external and internal dielectric constants of 80 and 1 (GB/SA solvation calculations (43Go) provide the same trends and accordingly results are not shown here ((19Go,39Go,44Go) for discussion). Intramolecular DNA entropy was determined using Schlitter's and Andricioaei-Karplus methods (45Go–47Go). As done in previous studies, end basepairs were removed from the analysis (9Go,19Go,39Go) to reduce noise in the calculations.

(1)

The free energy computed from Eq. 1 does not provide a direct measure of the stability of a structure, but for two oligonucleotides of the same sequence the difference in free energies gives a direct measure of the relative stability (see Eq. 2). Equivalent equations can be formulated for the entropic and effective energy contributions.

(2)

For a regular helix there is a linear relationship between the free energy (as computed in Eq. 1) and the total length of the oligonucleotide (see Eq. 3; (9Go,19Go,39Go)). This allowed us to obtain accurate estimates of the relative stability of two helical conformations by averaging data obtained from series of different MD simulations, thus reducing the statistical noise intrinsic to the use of Eq. 1.

(3)
where Gelon and Gnuc represent free-energy contributions associated for the elongation and nucleation of the nucleotide, respectively, and n is the length of the oligonucleotide in conformation X.

Equation 3 cannot be directly used in these systems, because in principle the free energy of the chimera depends not only on the length of the oligonucleotide, but also on the relative content of each type of structure (B, apH) and the number of junctions. A reasonable extension of Eq. 3 can be derived assuming that the thermodynamics of each constituent structure of the chimera is independent of the other. This yields to Eq. 4, where the total oligonucleotide (containing p steps) is divided in n steps of X and m steps of Y conformations, and z junctions linking the two types of structures. In a first approximation, the effective nucleation free energy of the chimera ({alpha}x/y) can be interpolated from those of the pure oligonucleotides using Eq. 5. We should note that the nucleation free energies of B and apH helices are very similar (19Go), reducing the uncertainties derived from the use of Eq. 5. We should also notice that Eqs. 4 and 5 can be easily rewritten in terms of effective energies.

(4)

(5)

Because the elongation and nucleation free energies of the B and apH helices are known from previous work (19Go) the free energy of junction can be easily derived using Eqs. 4 and 5 from the analysis of the free energies of the nine oligonucleotides of different lengths and compositions studied here.

Interaction properties of the structures were obtained from classical molecular interaction potential (cMIP) calculations for an O+ probe molecule (30Go–32Go,48Go). Solvation was represented by integrating the water population along the trajectory (30Go–32Go). The molecular flexibility was analyzed using principal component analysis method (49Go). The eigenvalues {{lambda}i} were manipulated (50Go) to derive information on the strength of the harmonic constant associated to the essential movements (Eq. 6). The eigenvectors ({nu}i) were used to describe the nature of the essential deformation modes and to compare the pattern of flexibility of two trajectories using our standard algorithms ((39Go,50Go–52Go); Eqs. 7 and 8), or the recently developed metrics (53Go), which introduces a Boltzmann-weighting scheme in the absolute (Eq. 9) and relative (Eq. 10) comparison indexes.

(6)
where k is Boltzmann's constant, T is the absolute temperature, and {lambda}i is the eigenvalue describing the essential movement i in Å2.

(7)
where A and B stand for two different trajectories of equal length, {nu} stands for the eigenvectors, and n is the number of motions that account for a given value of structural variance in the trajectories (typically only 10 motions account for ~80% of the variance). Only backbone atoms (up to C1') were considered in the comparison.

(8)
where self-similarity indexes for A and B trajectories ( and ) are obtained by comparing the first and last halves of the same trajectory. Both {gamma} and {kappa} are 1 for two identical trajectories and 0 when they sample orthogonal movements.

(9)

(10)
where {lambda}i is the eigenvalue (in Å2) associated to eigenvector {nu}i and {Delta}x is a displacement (in Å) common to all the modes, that is selected (J. R. Blas, A. Pérez, F. J. Luque, and M. Orozco, unpublished data) as the minimum one that makes negligible to the calculation of {xi} (Eq. 9) of eigenvectors (i, j) associated to high frequencies (those that are not needed to explain 99% of the variance). For coherence with the other indexes the sums are extended to the same active space as in Eqs. 9 and 10. Note that when {Delta}x = 0, Eqs. 9 and 10 converge to Eqs. 7 and 8.

Geometrical analysis was performed using the ptraj module in AMBER6.0 (38Go) and in-house programs, and helical analysis was carried out using Olson's X3DNA program (54Go) and Curves5.3 (55Go). Basepairs at the ends of the duplex were removed in all the cases.


   RESULTS AND DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Stability of trajectories of B/apH chimeras
The MD trajectories for the nine chimeras studied here (Table 1) are stable in the simulation time (2–5 ns) considered and the same is found for selected simulations extended to 10 ns or performed at 310 K. In all cases the root mean-square deviations (RMSDs) determined with regard to the MD-averaged structures are small, with values intermediate between those of B and apH helices (see Table 2), thus indicating that even the 2-ns trajectories are well converged. Interestingly, the content of B and apH duplex in the chimeras does not introduce major alterations in the RMSD with respect to the MD-averaged conformation (Table 2). The RMSD with respect to the starting structures (built from experimental data of the constituent fragments; see Table 2) are small (~2 Å in all the cases), and intermediate between those found for B and apH helices. In summary, the RMSD data suggest that the structures of chimeras seem to be very stable under these simulation conditions. Moreover, such structural stability is not affected by the length of the oligonucleotide nor the ratio between B and apH helices. The presence of junction, therefore, does not seem to introduce any dramatic effect on the structure of the duplex.


View this table:
[in this window]
[in a new window]
  		TABLE 2  All-atom root mean square deviations (Å) for duplexes containing B or apH helices (data from Cubero et al. (19Go)) or B/apH chimeras with respect to the MD-averaged and starting structures

 
The hydrogen-bond (H-bond) pattern for Watson-Crick and Hoogsteen pairings is well preserved for all the chimeras. The percentage of total hydrogen bonding is >93% in all cases (Table 3), as also found in B and apH helices (19Go), thus indicating that mixing of the two helices in the same structure does not alter the hydrogen bonding pattern. The hydrogen bonding at the junctions is also well preserved, the percentage of hydrogen bonds being similar to that found for the rest of the duplex. As found in our previous study (19Go), partial breathing movements were detected in the subnanosecond timescale for both Watson-Crick and Hoogsteen pairings. In general these movements last for ~0.5 ns and imply interchanges of hydrogen-bond donors and acceptors between adenine and thymine. For one of the trajectories (apH6B10; see Table 3), the "partial opening" movement lasts >2 ns and implies that at one Hoogsteen basepair of the junction the N6 (A) -> O4 (T) and N3 (T) -> N7 (A) interactions are lost while one N6 (A) -> O2 (T) hydrogen bond is formed leading to an overall reduction in the percentage of H-bonding detected during the entire trajectory (see Table 3). These "partial opening" events, nevertheless, do not introduce any major alteration in the general geometry of the helix.


View this table:
[in this window]
[in a new window]
  		TABLE 3  Occurrence (in % with respect to the maximum number of hydrogen-bond interactions) of canonical Watson-Crick and Hoogsteen hydrogen bonds for the different trajectories

 
In summary, all the chimeras seem to be as stable as the parent B and apH helices. No major alterations in the global structure nor in the specific pattern of basepair interactions are detected in the trajectories related to the existence of chimeras of two different helices in a common duplex.

Structural description of the chimeras
The general shape of the chimeras is very similar to that of B and apH helices (see Fig. 1). The different ratio of B and apH helices or the presence of B/apH junctions does not alter the general shape nor the helical parameters of the duplex. The periodicity and rise of the apH helix is similar to that of B-DNA (17Go–19Go) and do not change much in the different chimeras (see Table 4). The fact that phase angles for sugars in apH helices are typically smaller (5–7°) than those found in B-DNA is also visible in chimeras, because the average phase angle generally decreases as the content of apH helix increases. However, all the sugars are found in the South-Southeast region. All the structures (see Fig. 1 and Table 4) have a groove topology similar to that of B-DNA, but the average size of the minor groove decreases as the percentage of apH helix increases, as expected from the fact that apH helices show in general smaller widths of the minor groove than B-DNAs.



View larger version (27K):
[in this window]
[in a new window]
  		FIGURE 1  MD-averaged structure of 14-mer DNA duplexes containing (from left to right) none, two, four, six, and 14 antiparallel Hoogsteen (apH) steps. The B fragment always appears in black and the apH one in gray, and framed by the dotted box.

 

View this table:
[in this window]
[in a new window]
  		TABLE 4  Average values for selected helical parameters obtained in MD simulations of B and apH duplexes as well as for different chimeras of the 14-mer oligonucleotides

 
Very interestingly, the helical parameters describing B and apH helices are well maintained in the respective fragments (see Table 4). For example, the RMSD between the central 6-mer apH fragment in apH12, apH6B6, apH6B8, and apH6B10 varies between 1.0 and 1.1 Å. The maintenance of the structural characteristics of B and apH helices in the chimeric duplexes is especially clear in the shortest intrastrand C1'–C1' distance, which is the best single parameter to discriminate between B and apH helices. Thus, this parameter adopts the values expected for the canonical B and apH helices in the corresponding chimeric fragments. The junctions are characterized by a reduced twist and a slightly increased rise with respect to B and apH helices (Table 4). Apart from these changes, they do not appear to be structurally stressed nor produce any noticeable change in the global helical parameters of the duplex.

The three-dimensional distribution of classical molecular interaction potential (CMIP (48Go)) in chimeras is intermediate between those found in both B and apH helices (see Fig. 2). For B-DNA there is a strong region of favorable interaction with small cations along the minor groove. The size of such a region becomes progressively diminished as the content of apH helix in the duplex increases. All the helices are well solvated by 25–27 waters per step (water molecules at <3.5 Å from any heteroatom of DNA), with no major difference for B, apH, or chimeric duplexes. As described elsewhere (19Go), despite the different orientation of the adenine, the characteristic spine of hydration in the minor groove of B-DNA (56Go) is also clear in the apH helix and in all the chimeras (see Fig. 2). Furthermore, the specific hydration pattern of apH helices in the major groove (17Go–19Go) is found also here even for the smallest apH helical fragments (see Fig. 3 and compare with Fig. 5 in Cubero et al. (19Go)). In summary, there is a strong memory of the intrinsic interaction properties of B and apH helices in the chimeras. Accordingly, we suggest that even small apH fragments inserted in large B-DNA duplexes will display recognition properties very close to those of an apH helix, thus opening the possibility to expand (for a common sequence) the possibilities for specific DNA-protein recognition in the cell.



View larger version (46K):
[in this window]
[in a new window]
  		FIGURE 2  (Top) Classical molecular interaction potential (isocontour at –5 kcal/mol) using an O+ probe particle (see Methods). (Bottom) Solvation maps (isocontour for an apparent water density of 2.5 g/cm3) for the different 14-mer duplexes studied here.

 


View larger version (21K):
[in this window]
[in a new window]
  		FIGURE 3  Detail of the solvation along the grooves in the apH fragment of the apH2B12 duplex. The isocontour corresponds to an apparent water density of 2.5 g/cm3.

 


View larger version (8K):
[in this window]
[in a new window]
  		FIGURE 5  Dependence of the intramolecular entropy with the length of the apH segment in different duplexes. Entropy values (in cal/Kmol) were obtained using the Andricioaei-Karplus method, but Schlitter's approach yields equivalent results.

 
Dynamics properties of chimeras
The essential dynamics of B and apH duplexes and all the chimeras is dominated by bending and twisting deformation of the helices. Analysis of the 5-ns trajectories collected for 12-mer duplexes shows a large degree of similarity between pure duplexes and chimeras, as noted in relative similarity indexes ({kappa} and {delta}; see Eqs. 8 and 10) around 90% (see Table 5). Comparison of the nature of the essential movements for the 14- and 16-mer oligonucleotides requires some caution, because deformation modes might not be completely converged after 2-ns trajectories. With this caution in mind, we must notice that the results support the occurrence of similar similarity essential movements in B, apH, and chimeric duplexes (see data in Table 5 for the 16-mer duplexes). At this point, it is worth noting the similarity between normal ({kappa}) and Boltzmann's-weighted ({delta}) relative similarity indexes, which indicates that not only the nature of the essential movements, but also their relative contribution to the structural variance, is retained in the different duplexes. Overall, the presence of junctions and apH pairs does not seem to induce major changes in the natural deformation modes of the duplex.


View this table:
[in this window]
[in a new window]
  		TABLE 5  Comparison indexes for 2-ns trajectories for 16-mer duplexes and 5-ns trajectories for 12-mer duplexes

 
The essential movements of B and apH duplexes and chimeras are associated to very weak (<10 cal/mol Å2) harmonic force constants (see Eq. 6 and Fig. 4), which allow large deformations in the helices under room temperature. Interestingly, the increase in the percentage of apH steps leads to a parallel increase in all the harmonic constants (see Fig. 5). This finding suggests that the nature of the easiest deformation modes is similar for apH and B helices, but that deformation along these modes is more difficult for apH than for B fragments. Thus, B/apH chimeras will display an intermediate deformation pattern, but they will become more rigid as the apH content increases.



View larger version (14K):
[in this window]
[in a new window]
  		FIGURE 4  Harmonic force constants associated to deformations along the 10 first essential movements of B and apH duplexes and different chimeras (for 16-mer duplexes).

 
Stability of B/apH chimeras
The effective energy (Eq. 1) of the chimeras is mostly dominated by the total length of the duplex and neither by the B/apH ratio nor the position of the junctions (see Table 6). The effective energies, as well as their intramolecular and solvation free-energy components, for B and apH helices and for the different chimeras are very similar (the same result was obtained when GB/SA was used to compute the solvation term).


View this table:
[in this window]
[in a new window]
  		TABLE 6  Effective energy (<Eeffec>; bold), their solvation (<Gsol>; Roman), and intramolecular (<Eintra>; italics) contributions for the B, apH, and chimeric duplexes

 
The statistical analysis (see Methods and Eqs. 4 and 5) shows that the effective energy penalty due to the existence of two junctions in our B/apH helices (see Eq. 5) is only 3.4 ± 1.0 kcal/mol, that is, 1.7 kcal/mol per junction, which is <0.05% the total effective energy of the duplexes. Such a penalty is even reduced (0.7 kcal/mol x junction) when solvation free energies determined from GB/SA calculations are considered. In other words, the small destabilization (in terms of effective energy) required to insert an apH fragment into a B duplex appears to have little impact on the whole stability of a DNA duplex under physiological conditions.

The presence of syn-adenosines destabilizes the apH helix (~0.7 kcal/mol x step; see Cubero et al. (19Go) and Stolarski et al. (57Go), but this is mostly compensated by hydrogen bonding, which is expected to be more stable in the Hoogsteen pairing than in the Watson-Crick scheme (19Go,58Go,59Go). Thus, the analysis of all these trajectories confirms that on average the Hoogsteen H-bonds are around –1.2 ± 0.01 kcal/mol (see Table 7) more stable than the Watson-Crick ones, in good agreement with previous estimates (19Go). The Watson-Crick or Hoogsteen hydrogen bonds are equally stable in the junction or in the middle of the fragments (data not shown), and accordingly, the hydrogen-bond energy associated to a given junction is just the average between those of Watson-Crick and Hoogsteen pairs. The stacking energy of (A·T)-(T·A) dimers is the same irrespective of whether the hydrogen bond (the "·" here) pattern is Watson-Crick or Hoogsteen (Table 7). On the contrary, the stacking (T·A)-(A·T) is on average 3.8 kcal/mol more stable if the two dimers are hydrogen-bonded following the Hoogsteen scheme (Table 7). Despite the low twist values (Table 4), the stacking in the junction is quite favorable. Thus, for the (A·T)-(T·A) dimers there is only a loss of 0.5 kcal/mol with respect to B and apH stackings, whereas for the (T·A)-(A·T) dimers there is a loss of stability (1.3 kcal/mol) relative to apH stacking, but a gain (2.5 kcal/mol) with regard to B stacking. In summary, the d (AT)n DNA seems so flexible that the distortions needed to accommodate a junction between two different types of helices do not introduce any major destabilization in the structure. As noted in Supplementary Material, the conclusions on stability derived here from 2- to 5-ns trajectories remain valid when tested in 10-ns trajectories, or when more physiological temperatures (310 K) are considered.


View this table:
[in this window]
[in a new window]
  		TABLE 7  Average hydrogen-bond energy (kcal/mol) in Hoogsteen and Watson-Crick pairings in chimeras, B, and apH helices; stacking energy (kcal/mol) for the different steps and its intra- and intermolecular components

 
Previous calculations (19Go) suggested that whereas in terms of effective energy B and apH helices were equally stable, the former was favored by entropic reasons. As noted elsewhere (39Go,60Go,61Go), 2–5-ns simulations are in general too short to provide reliable measures of relative intramolecular entropy of different helices, which precludes the calculation of the total free energy using Eq. 1 and derived equations. However, for each oligonucleotide, comparison of the five estimates (B, apH, and the three chimeras) of the intramolecular entropy can be useful to identify general trends about the impact of chimeras on the entropy of the helix (19Go). It is also worth noting that entropy values obtained for selected cases from 10-ns trajectories are obviously larger, but in relative terms (apH, B, junctions) no significant change is found, confirming the validity of our conclusions (see Table 8). Chimeras have intramolecular entropies intermediate between those of B and apH helices (Table 8), and the increase in the content of apH helix in duplexes of constant length produces a linear (r2 ranging between 1.00 and 0.96) decrease in the entropy of the system, as shown in Fig. 5, thus reflecting the larger rigidity of apH helices (19Go). The fact that entropies of chimeras and pure helices (B and apH) can be fitted in the same regression line (see Fig. 5) demonstrates that the junctions does not introduce any major alteration in the global intramolecular entropy of the helix, and that no entropy penalty can be specifically assigned to the junctions.


View this table:
[in this window]
[in a new window]
  		TABLE 8  Entropy estimates (cal/Kmol) computed using Andricioaei & Karplus (in Roman) and Schlitter's (italics) methods for 2-ns simulations (values obtained for 10-ns trajectories are displayed for reference when available)

 

   CONCLUSIONS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
ACKNOWLEDGEMENTS
 REFERENCES
 
B/apH chimeras appear stable under simulation conditions. The global structural, dynamical, and molecular recognition properties of chimeras depend linearly on the percentage of apH and B helices. The structural transition from B to apH helices (and vice versa) is sharp (limited to one step) and happens without any dramatic alteration in the structure of the helical fragments. A certain local unwinding is the only remarkable conformational change detected at the junctions.

The fragments of B and apH helices in the chimera have a very strong memory of the structure, flexibility, and molecular interaction properties of the corresponding pure helices. Interestingly, such a memory does not depend on the length of the fragment, and even a very short apH fragment suffices to define local properties close to those of an apH helix.

B and apH helices possess similar effective energy. The B form is however entropically favored. The stability of the chimeras is very similar to that of the constituting helical fragments, without any dramatic lost of stability related to the presence of the junctions.

Overall, our simulations strongly suggest that, beside kinetic factors related to the anti-syn rotation of adenosine, the formation of short apH helices in the middle of long B-type duplexes might be possible if external factors (i.e., by interactions with proteins or drugs) reduce the flexibility of the duplex, thus compensating the entropic preference for the B form. This suggestion, in conjunction with the ability of small apH fragments to maintain a memory of their structural, dynamical, and reactive properties, support the possibility that apH fragments might play a greater than expected role in the control of gene expression.


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
REFERENCES
 
Supplementary material is available upon request to the corresponding authors.

We thank Prof. J. Subirana for many helpful discussions. We are indebted to X.-J. Lu and W. Olson for a copy of X3DNA, and help in the use of the code for a nonstandard structure. We also thank D. Bashford for a copy of his MEAD program.

We acknowledge the Centre de Supercomputació de Catalunya (CESCA) and the Spanish Ministry of Science and Technology (SAF2002-4282 and BIO2003-06848), and the Fundación BBVA for financial support.

Submitted on January 14, 2005; accepted for publication October 20, 2005.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Saenger, W. 1984. Principles of Nucleic Acid Structure. Springer-Verlag, New York, NY.

2. Bloomfield, V. A., D. N. Crothers, and I. Tinoco. 2000. Nucleic Acids: Structures, Properties and Functions. V. A. Bloomfield, D. N. Crothers, and I. Tinoco, editors. University Science Books, Sausalito, CA.

3. Blackburn, G. M., and M. J. Gait. 1990. Nucleic Acids in Chemistry and Biology. G. M. Blackburn and M. J. Gait, editors. IRL Press, Oxford, UK.

4. Robles, J., A. Grandas, E. Pedroso, F. J. Luque, R. Eritja, and M. Orozco. 2002. Nucleic acid triple helices: stability effects of nucleobase modifications. Curr. Org. Chem. 6:1333–1368.[CrossRef]

5. Gehring, K., J. L. Leroy, and M. Guéron. 1993. A tetrameric DNA structure with protonated cytosine-cytosine base pairs. Nature. 363:561–565.[CrossRef][Medline]

6. Gray, D. M., M. Vaughn, R. L. Ratliff, and F. N. Hayes. 1980. Circular dichroism spectra show that repeating dinucleotide DNAs may form helices in which every other base is looped out. Nucleic Acids Res. 8:3695–3707.[Abstract/Free Full Text]

7. Laughlan, G., A. I. Murchie, M. H. Norman, P. C. Moody, D. M. Lilley, and B. Luisi. 1994. The high-resolution crystal structure of a parallel-stranded guanine tetraplex. Science. 265:520–524.[Abstract/Free Full Text]

8. Wang, A. H., G. J. Quiqkey, F. J. Kolpak, J. L. Crawford, J. H. van Boom, G. van der Marel, and A. Rich. 1979. Molecular structure of a left-handed double helical DNA fragment at atomic resolution. Nature. 282:680–686.[CrossRef][Medline]

9. Cubero, E., F. J. Luque, and M. Orozco. 2001. Theoretical studies of d (A:T)-based parallel-stranded DNA duplexes. J. Am. Chem. Soc. 123:12018–12025.[CrossRef][Medline]

10. Cubero, E., A. Aviñó, B. G. de la Torre, M. Frieden, R. Eritja, F. J. Luque, C. González, and M. Orozco. 2002. Hoogsteen-based parallel-stranded duplexes of DNA. Effect of 8-amino-purine derivatives. J. Am. Chem. Soc. 124:3133–3142.[CrossRef][Medline]

11. Guzikevich-Guerstein, G., and Z. Shakked. 1996. A novel form of the DNA double helix imposed on the TATA-box by the TATA-binding protein. Nat. Struct. Biol. 3:32–37.[CrossRef][Medline]

12. Viswamitra, M. A., P. Kennard, P. G. Jones, G. M. Sheldrick, S. Salisbury, L. Falvello, and Z. Shakked. 1978. DNA double helical fragment at atomic resolution. Nature. 273:687–688.[CrossRef][Medline]

13. Viswamitra, M. A., Z. Shakked, P. G. Jones, G. M. Sheldrick, S. A. Salisbury, and O. Kennard. 1982. Structure of the deoxytetranucleotide d-pApTpApT and a sequence-dependent model for poly (dA-dT). Biopolymers. 21:513–533.[CrossRef][Medline]

14. Brahms, S., J. Brahms, and K. E. Van Holde. 1976. Nature of conformational changes in poly [d (A-T) · d (A-T)] in the premelting region. Proc. Natl. Acad. Sci. USA. 73:3453–3457.[Abstract/Free Full Text]

15. Davis, D. R., and R. L. Baldwin. 1963. X-ray studies on two synthetic DNA copolymers. J. Mol. Biol. 6:251–255.[Medline]

16. Radwan, M. M., and H. R. Wilson. 1982. Fibre and molecular structure of thymidylyl-3', 5'-deoxyadenosine. Int. J. Biol. Macromol. 4:145–149.[CrossRef]

17. Abrescia, N. G. A., A. Thompson, T. H. Dinh, and J. A. Subirana. 2002. Crystal structure of an antiparallel DNA fragment with Hoogsteen base pairing. Proc. Natl. Acad. Sci. USA. 99:2806–2811.[Abstract/Free Full Text]

18. Abrescia, N. G. A., C. González, C. Gouyette, and J. A. Subirana. 2004. X-ray and NMR studies of the DNA oligomer d (ATATAT): Hoogsteen base pairing in duplex DNA. Biochemistry. 43:4092–4100.[CrossRef][Medline]

19. Cubero, E., N. G. A. Abrescia, J. A. Subirana, F. J. Luque, and M. Orozco. 2003. Theoretical study of a new DNA structure: the antiparallel Hoogsteen duplex. J. Am. Chem. Soc. 125:14603–14612.[CrossRef][Medline]

20. Felsenfeld, G., D. R. Davis, and A. Rich. 1957. Formation of a three-stranded polynucleotide molecule. J. Am. Chem. Soc. 79:2023–2024.[CrossRef]

21. Pauling, L., and R. B. Corey. 1953. A proposed structure for the nucleic acids. Proc. Natl. Acad. Sci. USA. 39:84–97.[Free Full Text]

22. Aishima, J., R. K. Gitti, J. E. Noah, H. H. Gan, T. Schlick, and C. Wolberger. 2002. A Hoogsteen base pair embedded in undistorted B-DNA. Nucleic Acids Res. 30:5244–5252.[Abstract/Free Full Text]

23. Gilbert, D. E., G. A. van der Marel, J. H. van Boom, and H. Feigon. 1989. Unstable Hoogsteen base pairs adjacent to echinomycin binding sites within a DNA duplex. Proc. Natl. Acad. Sci. USA. 86:3006–3010.[Abstract/Free Full Text]

24. Patikoglou, G. A., J. L. Kim, L. Sun, S.-H. Yang, T. Kodadek, and S. K. Burley. 1999. TATA element recognition by the TATA box-binding protein has been conserved throughout evolution. Genes Dev. 13:3217–3230.[Abstract/Free Full Text]

25. Rice, P. A., S. Yang, K. Mizuuchi, and H. A. Nash. 1996. Crystal structure of an IHF-DNA complex: a protein-induced DNA U-turn. Cell. 87:1295–1306.[CrossRef][Medline]

26. Wang, A. H., G. Ughetto, G. J. Quigley, T. Hakoshima, G. A. van der Marel, J. H. van Boom, and A. Rich. 1984. The molecular structure of a DNA-triostin A complex. Science. 225:1115–1121.[Abstract/Free Full Text]

27. Moreau, J., L. Marcaud, F. Maschat, J. Kejzlarova-Lepesant, J. A. Lepesant, and K. Scherer. 1982. A + T-rich linkers define functional domains in eukaryotic DNA. Nature. 295:260–262.[CrossRef][Medline]

28. Oosumi, T., B. Garlick, and W. R. Belknap. 1995. Identification and characterization of putative transposable DNA elements in solanaceous plants and caenorhabditis elegans. Proc. Natl. Acad. Sci. USA. 92:8886–8890.[Abstract/Free Full Text]

29. Arnott, S., and D. W. L. Hukins. 1972. Optimized parameters for A-DNA and B-DNA. Biochem. Biophys. Res. Commun. 47:1504–1509.[CrossRef][Medline]

30. Shields, G. C., C. A. Laughton, and M. Orozco. 1997. Molecular dynamics simulations of the d (T·A·T) triple helix. J. Am. Chem. Soc. 119:7463–7469.[CrossRef]

31. Shields, G. C., C. A. Laughton, and M. Orozco. 1998. Molecular dynamics simulation of a PNA·DNA·PNA triple helix in aqueous solution. J. Am. Chem. Soc. 120:5895–5904.[CrossRef]

32. Soliva, R., C. A. Laughton, F. J. Luque, and M. Orozco. 1998. Molecular dynamics simulations in aqueous solution of triple helices containing d (G·C·C) trios. J. Am. Chem. Soc. 120:11226–11233.[CrossRef]

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

34. Ryckaert, J. P., G. Ciccotti, and H. J. C. Berendsen. 1977. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J. Comput. Phys. 23:327–341.[CrossRef]

35. Cheatham, T. E., P. Cieplak, and P. A. Kollman. 1999. A modified version of the Cornell et al. force field with improved sugar pucker phases and helical repeat. J. Biomol. Struct. Dyn. 16:845–862.[Medline]

36. Cornell, W. D., P. Cieplak, C. I. Bayly, I. R. Gould, K. M. Merz, D. M. Fergurson, D. C. Spellmeyer, T. Fox, J. W. Caldwell, and P. A. Kollman. 1995. A second generation force field for the simulation of proteins, nucleic acids, and organic molecules. J. Am. Chem. Soc. 117:5179–5197.[CrossRef]

37. Jorgensen, W. L., J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein. 1983. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79:926–935.[CrossRef]

38. Case, D. A., D. A. Pearlman, J. W. Caldwell, T. E. Cheatham, W. S. Ross, C. L. Simmerling, T. A. Darden, K. M. Merz, R. V. Stanton, A. L. Cheng, J. J. Vincent, M. Crowley, et al. 1999. AMBER 6.0. University of California, San Francisco, CA.

39. Orozco, M., A. Perez, A. Noy, and F. J. Luque. 2003. Theoretical methods for the simulation of nucleic acids. Chem. Soc. Rev. 32:350–364.[CrossRef][Medline]

40. Srinivasan, J., T. E. Cheatham, P. Cieplak, P. A. Kollman, and D. A. Case. 1998. Continuum solvent studies of the stability of DNA, RNA, and phosphoramidate-DNA helices. J. Am. Chem. Soc. 120:9401–9409.[CrossRef]

41. Bashford, D., and K. Gerwert. 1992. Electrostatic calculations of the pKa values of ionizable groups in bacteriorhodopsin. J. Mol. Biol. 224:473–486.[CrossRef][Medline]

42. Bashford, D. 1997. An object-oriented programming suite for electrostatic effects in biological molecules. In Scientific Computing in Object-Oriented Parallel Environments. Y. Ishikawa, R. Rodney, R. R. Oldehoeft, J. V. M. Reynders, and M. Tholburn, editors. Lecture Notes in Computer Science. Springer, Berlin, Germany. 1343:233–240.

43. Still, W. C., A. Tempcyck, R. C. Hawley, and T. Hendrickson. 1990. Semianalytical treatment of solvation for molecular mechanics and dynamics. J. Am. Chem. Soc. 112:6127–6129.[CrossRef]

44. Orozco, M., and F. J. Luque. 2000. Theoretical methods for the description of the solvent effect in biomolecular systems. Chem. Rev. 100:4187–4226.[CrossRef][Medline]

45. Andricioaei, I., and M. Karplus. 2001. On the calculation of entropy from covariance matrices of the atomic fluctuations. J. Chem. Phys. 115:6289–6292.[CrossRef]

46. Schlitter, J. 1993. Estimation of absolute and relative entropies of macromolecules using the covariance matrix. Chem. Phys. Lett. 215:617–621.[CrossRef]

47. Harris, S. A., E. Gavathiotis, M. S. Searle, M. Orozco, and C. A. Laughton. 2001. Cooperativity in drug-DNA recognition: a molecular dynamics study. J. Am. Chem. Soc. 123:12658–12663.[CrossRef][Medline]

48. Gelpí, J. L., S. G. Kalko, X. de la Cruz, X. Barril, J. Cirera, F. J. Luque, and M. Orozco. 2001. Classical molecular interaction potentials: improved setup procedure in molecular dynamics simulations of proteins. Proteins. 45:428–437.[CrossRef][Medline]

49. Sherer, E., S. A. Harris, R. Soliva, M. Orozco, and C. A. Laughton. 1999. Molecular dynamics studies of DNA A-tract structure and flexibility. J. Am. Chem. Soc. 121:5981–5991.[CrossRef]

50. Noy, A., A. Pérez, F. Lankas, F. J. Luque, and M. Orozco. 2004. Relative flexibility of DNA and RNA: a molecular dynamics study. J. Mol. Biol. 343:627–638.[CrossRef][Medline]

51. Pérez, A., A. Noy, F. Lankas, F. J. Luque, and M. Orozco. 2004. The relative flexibility of B-DNA and A-RNA duplexes: database analysis. Nucleic Acids Res. 32:6144–6151.[Abstract/Free Full Text]

52. Rueda, M., S. G. Kalko, F. J. Luque, and M. Orozco. 2003. The structure and dynamics of DNA in the gas phase. J. Am. Chem. Soc. 125:8007–8014.[CrossRef][Medline]

53. Pérez, A., J. R. Blas, M. Rueda, J. M. López-Bes, X. de la Cruz, and M. Orozco. 2005. Exploring the essential dynamics of B-DNA. J. Chem. Theory and Comput. 1:790–800.[CrossRef]

54. Lu, X.-J., and W. K. Olson. 2003. 3DNA: a software package for the analysis, rebuilding and visualization of three-dimensional nucleic acid structures. Nucleic Acids Res. 31:5108–5121.[Abstract/Free Full Text]

55. Lavery, R., and J. Sklenar. 1988. The definition of generalized helicoidal parameters and of axis curvature of irregular nucleic acids. J. Biomol. Struct. Dyn. 6:63–91.[Medline]

56. Drew, H. R., R. M. Wing, T. Tanako, C. Broka, S. Tanaka, K. Itakura, and R. Dickerson. 1981. Structure of a B-DNA dodecamer: conformation and dynamics. Proc. Natl. Acad. Sci. USA. 78:2179–2183.[Abstract/Free Full Text]

57. Stolarski, R., A. Pohorille, L. Dudycz, and D. Shugar. 1980. Comparison of theoretical and experimental approaches to determination of conformation of nucleosides about the glycosidic bond. Biochim. Biophys. Acta. 610:1–19.[Medline]

58. Hoogsteen, K. 1959. The structure of crystals containing a hydrogen-bonded complex of 1-methylthymine and 9-methyladenine. Acta Crystallogr. 12:822–823.[CrossRef]

59. Sponer, J., J. Leszczynski, and P. Hobza. 1996. Structures and energies of hydrogen-bonded DNA base pairs. A nonempirical study with inclusion of electron correlation. J. Phys. Chem. 100:1965–1974.[CrossRef]

60. Ponomarev, S. Y., K. M. Thayer, and D. L. Beveridge. 2004. Ion motions in molecular dynamics simulations on DNA. Proc. Natl. Acad. Sci. USA. 101:14771–14775.[Abstract/Free Full Text]

61. Rueda, M., E. Cubero, C. A. Laughton, and M. Orozco. 2004. Exploring the counterion atmosphere around DNA: what can be learned from molecular dynamics simulations? Biophys. J. 87:800–811.[Abstract/Free Full Text]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.105.059535v1
90/3/1000    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Cubero, E.
Right arrow Articles by Orozco, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Cubero, E.
Right arrow Articles by Orozco, M.

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2006 by the Biophysical Society.