Parallel and antiparallel Holliday junctions differ in structure and stability

Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations

Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.

Effect of DNA supercoiling on the geometry of holliday junctions

Unusual DNA conformations including cruciforms play an important role in gene regulation and various DNA transactions. Cruciforms are also the models for Holliday junctions, the transient DNA conformations critically involved in DNA homologous and site-specific recombination, repair, and replication. Although the conformations of immobile Holliday junctions in linear DNA molecules have been analyzed with the use of various techniques, the role of DNA supercoiling has not been studied systematically. We utilized atomic force microscopy (AFM) to visualize cruciform geometry in plasmid DNA with different superhelical densities at various ionic conditions. Both folded and unfolded conformations of the cruciform were identified, and the data showed that DNA supercoiling shifts the equilibrium between folded and unfolded conformations of the cruciform toward the folded one. In topoisomers with low superhelical density, the population of the folded conformation is 50-80%, depending upon the ionic strength of the buffer and a type of cation added, whereas in the sample with high superhelical density, this population is as high as 98-100%. The time-lapse studies in aqueous solutions allowed us to observe the conformational transition of the cruciform directly. The time-dependent dynamics of the cruciform correlates with the structural changes revealed by the ensemble-averaged analysis of dry samples. Altogether, the data obtained show directly that DNA supercoiling is the major factor determining the Holliday junction conformation.

The regulatory region of the human desmocollin 3 promoter forms a DNA four-way junction

Adhesion between desmosomal junctions is mediated by structural proteins of the cadherin family, viz. three desmocollins (DSC) and three desmogleins (DSG). Promoter and primer extension analysis of human DSC3 showed a TATA-less sequence initiating transcription via a cluster of sites upstream of the coding region. Deletion analysis of 1 kb of the promoter showed that expression is regulated between --303 and --203 bp upstream of the start-site of translation. Tertiary structure analysis of this cis-active region (cis 1) revealed a potential DNA 4-way junction which is notably G/C-rich in sequence. PAGE analysis of this region identified four differently migrating forms of the DNA. Structure-specific cleavage of the DNA with bacteriophage T7 endonuclease I showed the slowest migrating form to be either an extended/cruciform or stacked-X 4-way junction. DNA-binding, gel retardation assays of the cis 1 region showed distinct DNA-protein complexes and by competition experiments and using purified junction DNA we show that one of these complexes bound with both sequence and structure specificity to the 4-way junction DNA.

Thermodynamic properties of an intramolecular DNA four-way junction

We have investigated the thermodynamic properties of two homologous DNA four-way junctions, J4 and J4M, based on 46-mer linear DNA molecules. J4 and J4M have the same base sequence with the only difference that the latter contains an uncharged methylene-acetal linkage, -O3'-CH2-O5', instead of the phosphodiester linkage, -O3'-PO2-O5'-, between the residues T18 and C19. The comparison of the thermal unfolding of the J4 junction and J4M junction serves to investigate the effect of the uncharged methylene-acetal linkage on the stability of the junction. Our analysis is based on CD, UV absorbance spectroscopy, DSC, and chemical footprinting. The aim is to characterize in detail the structure and stability of the junctions. As demonstrated before by NMR, in the presence of 5 mM MgCl2 +/- 50 mM NaCl, both J4 and J4M form a complete four-way junction. This is now evidenced by protection from OsO4 cleavage (chemical footprinting). We can assume that full base pairing occurs throughout the arms even at the center of the junction. CD spectra suggest that the helices within the junctions adopt the regular B-DNA conformation. Almost identical melting temperatures and unfolding enthalpies are obtained for J4 and J4M both by UV and DSC. Furthermore, the Van't Hoff enthalpy (DeltaHVH) derived from UV melting equals the calorimetric enthalpy (DeltaHcal), which means that the melting process of the structures proceeds in a two-state manner. All results taken together support the conclusion that there are no major conformational and energetic differences between J4 and J4M. The inclusion of the uncharged methylene-acetal group into the junction has no effect on its stability.

DNA-sequence asymmetry directs the alignment of recombination sites in the FLP synaptic complex

The FLP recombinase promotes site-specific recombination in the 2 micrometer circle of Saccharomyces cerevisiae. FLP recognizes a 48 bp target site (FLP recombination target, or FRT) consisting of three 13 bp protein binding sites, or symmetry elements, flanking an 8 bp spacer region. Efficient recombination also occurs with DNA substrates that have minimal FRT sites, consisting only of the spacer and two surrounding 13 bp symmetry elements arranged in inverse orientation; thus, the wild-type spacer sequence is the main asymmetric feature of the minimal recombination site. FLP carries out recombination with many minimal target sites bearing symmetric or asymmetric mutant spacer sequences; however, the overall directionality of recombination defined in terms of inversion or excision of a DNA domain is determined by spacer-sequence asymmetry. In order to evaluate the potential influence of spacer-sequence asymmetry on structures formed during early steps in recombination, we used electron microscopy to investigate the structure of the FLP synaptic complex, which is the intermediate protein-DNA complex involved in site pairing and strand exchange. Using linear substrate DNAs that have minimal FRTs with wild-type spacer sequences, we find that 85 to 90% of the FLP synaptic complexes examined contain the two FRTs aligned in parallel. This strong preference for parallel site alignment stands in contrast with prevailing models for lambda integrase-class recombination systems, which postulate antiparallel site alignment, and results from biophysical studies on synthetic, immobile four-way DNA junctions. Our results show that the strong preference for parallel alignment can be attributed to conformational preferences of Holliday junctions present in the synaptosome.

Sequence dependence of branch migratory minima

The Holliday junction is a central intermediate in the process of genetic recombination. The position of its branch-point can relocate through an isomerization known as branch migration. This migration occurs because the branch-point is flanked by homologous symmetry. All attempts at modeling the kinetics of branch migration have relied on the assumption that branch migration minima are sequence-independent. We have tested that assumption here, using a competition assay based on symmetric immobile branched junctions; these are junctions that cannot undergo branch migration, despite the fact that they are flanked by homology. The assay used is predicated on the non-association of strands displaced in the assay; we have tested this assumption, and have performed our experiments under conditions where we know that it is true. We have measured the free energy of relocating a branched junction from a fixed non-homologous sequence to all possible dimeric symmetric sequences. We find that the assumption of sequence-independence is often valid, but that it is not universally true. We find that the flanking sequences can have a marked effect on the free energy measured, both for extensions of symmetry and for reversals of flanking nucleotides. We have varied the temperature in our experiments, and have derived both enthalpies and entropies for the different sequences. The entropies are largely unfavorable, whereas the enthalpies are largely favorable; regardless of the signs of these quantities, we see that this is another system where enthalpy-entropy compensation is operative.

Structure of the Holliday junction intermediate in Cre-loxP site-specific recombination

We have determined the X-ray crystal structures of two DNA Holliday junctions (HJs) bound by Cre recombinase. The HJ is a four-way branched structure that occurs as an intermediate in genetic recombination pathways, including site-specific recombination by the lambda-integrase family. Cre recombinase is an integrase family member that recombines 34 bp loxP sites in the absence of accessory proteins or auxiliary DNA sequences. The 2.7 A structure of Cre recombinase bound to an immobile HJ and the 2.5 A structure of Cre recombinase bound to a symmetric, nicked HJ reveal a nearly planar, twofold-symmetric DNA intermediate that shares features with both the stacked-X and the square conformations of the HJ that exist in the unbound state. The structures support a protein-mediated crossover isomerization of the junction that acts as the switch responsible for activation and deactivation of recombinase active sites. In this model, a subtle isomerization of the Cre recombinase-HJ quaternary structure dictates which strands are cleaved during resolution of the junction via a mechanism that involves neither branch migration nor helical restacking.

Structure and dynamics of supercoil-stabilized DNA cruciforms

Understanding DNA function requires knowledge of the structure of local, sequence-dependent conformations that can be dramatically different from the B-form helix. One alternative DNA conformation is the cruciform, which has been shown to have a critical role in the initiation of DNA replication and the regulation of transcription in certain systems. In addition, cruciforms provide a model system for structural studies of Holliday junctions, intermediates in homologous DNA recombination. Cruciforms are not thermodynamically stable in linear DNA due to branch point migration, which makes their study using many biophysical techniques problematic. Atomic Force Microscopy (AFM) was applied to visualize cruciforms in negatively supercoiled plasmid DNA. Cruciforms are seen as clear-cut extrusions on the DNA filament with the lengths of the arms consistent with the size of the hairpins expected from a 106 bp inverted repeat. The cruciform exists in two different conformations, an extended one with the angle of ca. 180 degrees between the hairpin arms and a compact, X-type conformation, with acute angles between the hairpin arms and the main DNA strands. The ratio of molecules with the different conformations of cruciforms depends on ionic conditions. In the presence of high salt or Mg cations, a compact, X-type conformation is highly preferable. Remarkably, the X-conformation was highly mobile allowing the cruciform arms to adopt a parallel orientation. The structure observed is consistent with a model of the Holliday junction with a parallel orientation of the exchanging strands.

Effect of DNA topology on Holliday junction resolution by Escherichia coli RuvC and bacteriophage T7 endonuclease I

Holliday junctions are key intermediates in homologous genetic recombination. Their resolution requires specialised nucleases that nick pairs of strands at the junction point, leading to the separation of mature recombinants. Resolution occurs in either of two orientations, according to which DNA strands are cut. We show that DNA topology can determine the efficiency and outcome of a recombination reaction. Using two Holliday junction resolvases, Escherichia coli RuvC protein and T7 endonuclease I, we observed that supercoiled figure-8 DNA molecules containing Holliday junctions were resolved with a specific orientation bias, and that this bias was reversed by the presence of a topological tether (catenation). In contrast, when all topological constraints were removed by restriction digestion, the recombination intermediates were resolved equally in the two orientations. These results show that topological constraints affecting Holliday junction structure influence the orientation of resolution by cellular resolvases.

Action of site-specific recombinases XerC and XerD on tethered Holliday junctions

In Xer site-specific recombination, two related recombinases, XerC and XerD, mediate the formation of recombinant products using Holliday junction-containing DNA molecules as reaction intermediates. Each recombinase catalyses the exchange of one pair of specific strands. By using synthetic Holliday junction-containing recombination substrates in which two of the four arms are tethered in an antiparallel configuration by a nine thymine oligonucleotide, we show that XerD catalyses efficient strand exchange only when its substrate strands are 'crossed'. XerC also catalyses very efficient strand exchange when its substrate strands are 'crossed', though it also appears to be able to mediate strand exchange when its substrate strands are 'continuous'. By using chemical probes of Holliday junction structure in the presence and absence of bound recombinases, we show that recombinase binding induces unstacking of the bases in the centre of the recombination site, indicating that the junction branch point is positioned there and is distorted as a consequence of recombinase binding.

Three-way and four-way junctions in DNA: a conformational viewpoint

DNA junctions are potential intermediates in various important genetic processes, including mutagenesis and recombination. The quantity of research carried out in this area is rapidly increasing. Examples of three-way and four-way junctions are now relatively well characterized and a few common properties have been recognized, of which the most important is the tendency of junctions to fold into one or more coaxially stacked helical conformations or cross-over structures.

Relative stabilities of DNA three-way, four-way and five-way junctions (multi-helix junction loops): unpaired nucleotides can be stabilizing or destabilizing

Competition binding and UV melting studies of a DNA model system consisting of three, four or five mutually complementary oligonucleotides demonstrate that unpaired bases at the branch point stabilize three- and five-way junction loops but destabilize four-way junctions. The inclusion of unpaired nucleotides permits the assembly of five-way DNA junction complexes (5WJ) having as few as seven basepairs per arm from five mutually complementary oligonucleotides. Previous work showed that 5WJ, having eight basepairs per arm but lacking unpaired bases, could not be assembled [Wang, Y.L., Mueller, J.E., Kemper, B. and Seeman, N.C. (1991) Biochemistry, 30, 5667-5674]. Competition binding experiments demonstrate that four-way junctions (4WJ) are more stable than three-way junctions (3WJ), when no unpaired bases are included at the branch point, but less stable when unpaired bases are present at the junction. 5WJ complexes are in all cases less stable than 4WJ or 3WJ complexes. UV melting curves confirm the relative stabilities of these junctions. These results provide qualitative guidelines for improving the way in which multi-helix junction loops are handled in secondary structure prediction programs, especially for single-stranded nucleic acids having primary sequences that can form alternative structures comprising different types of junctions.

A model for repair of radiation-induced DNA double-strand breaks in the extreme radiophile Deinococcus radiodurans

The bacterium Deinococcus (formerly Micrococcus) radiodurans and other members of the eubacterial family Deinococaceae are extremely resistant to ionizing radiation and many other agents that damage DNA. Stationary phase D. radiodurans exposed to 1.0-1.5 Mrad gamma-irradiation sustains > 120 DNA double-strand breaks (dsbs) per chromosome; these dsbs are mended over a period of hours with 100% survival and virtually no mutagenesis. This contrasts with nearly all other organisms in which just a few ionizing radiation induced-dsbs per chromosome are lethal. In this article we present an hypothesis that resistance of D. radiodurans to ionizing radiation and its ability to mend radiation-induced dsbs are due to a special form of redundancy wherein chromosomes exist in pairs, linked to each other by thousands of four-stranded (Holliday) junctions. Thus, a dsb is not a lethal event because the identical undamaged duplex is nearby, providing an accurate repair template. As addressed in this article, much of what is known about D. radiodurans suggests that it is particularly suited for this proposed novel form of DNA repair.

Symmetric immobile DNA branched junctions

Branch migration is an isomerization of Holliday recombination intermediates that arises from their homologous (2-fold) sequence symmetry. This isomerization relocates the branch point in an apparently random fashion and thereby complicates the study of the physical and structural properties of these structures. For the past decade, these properties have been studied in low-symmetry immobile junctions, whose sequence asymmetry eliminates branch migration. The asymmetric findings of many of these studies suggest the need for a system combining both immobility and symmetry. Double-crossover DNA molecules have been used to create molecules with both these properties. Immobility is achieved by flanking one crossover with a symmetric junction and the other crossover with an asymmetric junction. Close torsional coupling between the two junctions renders the symmetric junction immobile. These molecules will enable the characterization of thermodynamic, structural, dynamic, liganding, and substrate properties of symmetric branched DNA molecules in a sequence-specific fashion.

The structure of branched DNA species

Branched DNA molecules provide a challenging set of structural problems. Operationally we define branched DNA species as molecules in which double helical segments are interrupted by abrupt discontinuities, and we draw together a number of different kinds of structure in the class, including helical junctions of different orders, and base bulges (Fig. 1).

DNA double-crossover molecules

DNA molecules containing two crossover sites between helical domains have been suggested as intermediates in recombination processes involving double-strand breaks. We have modeled these double-crossover structures in an oligonucleotide system. Whereas the relative orientations of the helical domains must be specified in designing these molecules, there are two broad classes of the molecules, the parallel, DP, and antiparallel, DA, molecules. The distance between crossover points must be specified as multiples of half-turns, in order to avoid torsional stress in this system; hence, there are two further subdivisions, those double-crossover molecules separated by odd, O, and even, E, numbers of half-turns. In addition, the parallel molecules with odd numbers of half-turns between crossovers must be divided into those with an excess major or wide-groove separation, W, or those with an excess minor- or narrow-groove separation, N. We have constructed models of all five of these classes, DAE, DAO, DPE, DPOW, and DPON. DPE molecules containing 1 and 2 helical turns between crossovers have been constructed; the DAE molecule contains 1 turn between crossovers, and the DAO, DPOW, and DPON molecules contain 1.5 helical turns between crossovers. None of the parallel molecules is well-behaved; the molecules either dissociate or form multimers when visualized on native polyacrylamide gels. In contrast, antiparallel molecules form single bands when assayed in this fashion. Hydroxyl radical autofootprinting analysis of these molecules reveals protection at expected sites of crossover and of occlusion, suggesting that all the complexes contain linear helix axes that are roughly coplanar between crossovers. However, the DPOW molecule and the DPE molecule with 2 turns between crossovers show decreased protection in the portion between crossovers, suggesting that their helices may bow in response to charge repulsion. We conclude that the helices between parallel double crossovers must be shielded from each other or distorted from linearity if they are to participate in recombination. We have analyzed the possibilities of branch migration and crossover isomerization in double-crossover molecules. Parallel molecules need no sequence symmetry beyond homology to branch migrate, but the sequence symmetry requirements for antiparallel molecules restrict migration to directly repetitive segments that iterate the sequence between crossovers. Crossover isomerization appears to be a very complex process in parallel double-crossover molecules, suggesting that it may be catalyzed by topoisomerases if it occurs within the cell.

Thermodynamics of DNA branching

Branched DNA molecules arise transiently as intermediates in genetic recombination or on extrusion of cruciforms from covalent circular DNA duplexes that contain palindromic sequences. The free energy of these structures relative to normal DNA duplexes is of interest both physically and biologically. Oligonucleotide complexes that can form stable branched structures, DNA junctions, have made it possible to model normally unstable branched states of DNA such as Holliday recombinational intermediates. We present here an evaluation of the free energy of creating four-arm branch points in duplex DNA, using a system of two complementary junctions and four DNA duplexes formed from different combinations of the same set of eight 16-mer strands. The thermodynamics of formation of each branched structure from the matching pair of intact duplexes have been estimated in two experiments. In the first, labeled strands are allowed to partition between duplexes and junctions in a competition assay on polyacrylamide gels. In the second, the heats of forming branched or linear molecules from the component strands have been determined by titration microcalorimetry at several temperatures. Taken together these measurements allow us to determine the standard thermodynamic parameters for the process of creating a branch in an otherwise normal DNA duplex. The free energy for reacting two 16-mer duplexes to yield a four-arm junction in which the branch site is incapable of migrating is + 1.1 (+/- 0.4) kcal mol-1 (at 18 degrees C, 10 mM-Mg2+). Analysis of the distribution of duplex and tetramer products by electrophoresis confirms that the free energy difference between the four duplexes and two junctions is small at this temperature. The associated enthalpy change at 18 degrees C is +27.1 (+/- 1.3) kcal mol-1, while the entropy is +89 (+/- 30) cal K-1 mol-1. The free energy for branching is temperature dependent, with a large unfavorable enthalpy change compensated by a favorable entropy term. Since forming one four-stranded complex from two duplexes should be an entropically unfavorable process, branch formation is likely to be accompanied by significant changes in hydration and ion binding. A significant apparent delta Cp is also observed for the formation of one mole of junction, +0.97 (+/-0.05) kcal deg-1 mol-1.