Definitions and analysis of DNA Holliday junction geometry

Structural Motifs at the Telomeres and Their Role in Regulatory Pathways

Telomeres are specialized structures, found at the ends of linear chromosomes in eukaryotic cells, that play a crucial role in maintaining the stability and integrity of genomes. They are composed of repetitive DNA sequences, ssDNA overhangs, and several associated proteins. The length of telomeres is linked to cellular aging in humans, and deficiencies in their maintenance are associated with various diseases. Key structural motifs at the telomeres serve to protect vulnerable chromosomal ends. Telomeric DNA also has the ability to form diverse complex DNA higher-order structures, including T-loops, D-loops, R-loops, G-loops, G-quadruplexes, and i-motifs, in the complementary C-rich strand. While many essential proteins at telomeres have been identified, the intricacies of their interactions and structural details are still not fully understood. This Perspective highlights recent advancements in comprehending the structures associated with human telomeres. It emphasizes the significance of telomeres, explores various telomeric structural motifs, and delves into the structural biology surrounding telomeres and telomerase. Furthermore, telomeric loops, their topologies, and the associated proteins that contribute to the safeguarding of telomeres are discussed.

Molecular Dynamic Studies of Dye-Dye and Dye-DNA Interactions Governing Excitonic Coupling in Squaraine Aggregates Templated by DNA Holliday Junctions

Dye molecules, arranged in an aggregate, can display excitonic delocalization. The use of DNA scaffolding to control aggregate configurations and delocalization is of research interest. Here, we applied Molecular Dynamics (MD) to gain an insight on how dye-DNA interactions affect excitonic coupling between two squaraine (SQ) dyes covalently attached to a DNA Holliday junction (HJ). We studied two types of dimer configurations, i.e., adjacent and transverse, which differed in points of dye covalent attachments to DNA. Three structurally different SQ dyes with similar hydrophobicity were chosen to investigate the sensitivity of excitonic coupling to dye placement. Each dimer configuration was initialized in parallel and antiparallel arrangements in the DNA HJ. The MD results, validated by experimental measurements, suggested that the adjacent dimer promotes stronger excitonic coupling and less dye-DNA interaction than the transverse dimer. Additionally, we found that SQ dyes with specific functional groups (i.e., substituents) facilitate a closer degree of aggregate packing via hydrophobic effects, leading to a stronger excitonic coupling. This work advances a fundamental understanding of the impacts of dye-DNA interactions on aggregate orientation and excitonic coupling.

Molecular dynamics simulations of cyanine dimers attached to DNA Holliday junctions

Dye aggregates and their excitonic properties are of interest for their applications to organic photovoltaics, non-linear optics, and quantum information systems. DNA scaffolding has been shown to be effective at promoting the aggregation of dyes in a controllable manner. Specifically, isolated DNA Holliday junctions have been used to achieve strongly coupled cyanine dye dimers. However, the structural properties of the dimers and the DNA, as well as the role of Holliday junction isomerization are not fully understood. To study the dynamics of cyanine dimers in DNA, molecular dynamics simulations were carried out for adjacent and transverse dimers attached to Holliday junctions in two different isomers. It was found that dyes attached to adjacent strands in the junction exhibit stronger dye-DNA interactions and larger inter-dye separations compared to transversely attached dimers, as well as end-to-end arrangements. Transverse dimers exhibit lower inter-dye separations and more stacked configurations. Furthermore, differences in Holliday junction isomer are analyzed and compared to dye orientations. For transverse dyes exhibiting the smaller inter-dye separations, excitonic couplings were calculated and shown to be in agreement with experiment. Our results suggested that dye attachment locations on DNA Holliday junctions affect dye-DNA interactions, dye dynamics, and resultant dye orientations which can guide the design of DNA-templated cyanine dimers with desired properties.

Molecular dynamics simulations reveal dye attachment and DNA Holliday junction isomer effects on dye dimer orientations and excitonic couplings. These simulations can guide synthesis and experiments of dye-DNA structures for excitonic applications.

High-throughput techniques enable advances in the roles of DNA and RNA secondary structures in transcriptional and post-transcriptional gene regulation

The most stable structure of DNA is the canonical right-handed double helix termed B DNA. However, certain environments and sequence motifs favor alternative conformations, termed non-canonical secondary structures. The roles of DNA and RNA secondary structures in transcriptional regulation remain incompletely understood. However, advances in high-throughput assays have enabled genome wide characterization of some secondary structures. Here, we describe their regulatory functions in promoters and 3’UTRs, providing insights into key mechanisms through which they regulate gene expression. We discuss their implication in human disease, and how advances in molecular technologies and emerging high-throughput experimental methods could provide additional insights.

The influence of Holliday junction sequence and dynamics on DNA crystal self-assembly

The programmable synthesis of rationally engineered crystal architectures for the precise arrangement of molecular species is a foundational goal in nanotechnology, and DNA has become one of the most prominent molecules for the construction of these materials. In particular, branched DNA junctions have been used as the central building block for the assembly of 3D lattices. Here, crystallography is used to probe the effect of all 36 immobile Holliday junction sequences on self-assembling DNA crystals. Contrary to the established paradigm in the field, most junctions yield crystals, with some enhancing the resolution or resulting in unique crystal symmetries. Unexpectedly, even the sequence adjacent to the junction has a significant effect on the crystal assemblies. Six of the immobile junction sequences are completely resistant to crystallization and thus deemed “fatal,” and molecular dynamics simulations reveal that these junctions invariably lack two discrete ion binding sites that are pivotal for crystal formation. The structures and dynamics detailed here could be used to inform future designs of both crystals and DNA nanostructures more broadly, and have potential implications for the molecular engineering of applied nanoelectronics, nanophotonics, and catalysis within the crystalline context.

Engineered crystal architectures from DNA have become a foundational goal for nanotechnological precise arrangement. Here, the authors systematically investigate the structures of 36 immobile Holliday junction sequences and identify the features allowing the crystallisation of most of them, while 6 are considered fatal.

Computational investigation of the impact of core sequence on immobile DNA four-way junction structure and dynamics

Immobile four-way junctions (4WJs) are core structural motifs employed in the design of programmed DNA assemblies. Understanding the impact of sequence on their equilibrium structure and flexibility is important to informing the design of complex DNA architectures. While core junction sequence is known to impact the preferences for the two possible isomeric states that junctions reside in, previous investigations have not quantified these preferences based on molecular-level interactions. Here, we use all-atom molecular dynamics simulations to investigate base-pair level structure and dynamics of four-way junctions, using the canonical Seeman J1 junction as a reference. Comparison of J1 with equivalent single-crossover topologies and isolated nicked duplexes reveal conformational impact of the double-crossover motif. We additionally contrast J1 with a second junction core sequence termed J24, with equal thermodynamic preference for each isomeric configuration. Analyses of the base-pair degrees of freedom for each system, free energy calculations, and reduced-coordinate sampling of the 4WJ isomers reveal the significant impact base sequence has on local structure, isomer bias, and global junction dynamics.

Holliday Junctions Formed from Human Telomeric DNA

Cells have evolved inherent mechanisms, like homologous recombination (HR), to repair damaged DNA. However, repairs at telomeres can lead to genomic instability, often associated with cancer. While most rapidly dividing cells employ telomerase, the others maintain telomere length through HR-dependent alternative lengthening of telomeres (ALT) pathways. Here we describe the crystal structures of Holliday junction intermediates of the HR-dependent ALT mechanism. Using an extended human telomeric repeat, we also report the crystal structure of two Holliday junctions in close proximity, which associate together through strand exchange to form a hemicatenated double Holliday junction. Our combined structural results demonstrate that ACC nucleotides in the C-rich lagging strand (5'-CTAACCCTAA-3') at the telomere repeat sequence constitute a conserved structural feature that constrains crossover geometry and is a preferred site for Holliday junction formation in telomeres.

Structure of the Holliday junction: applications beyond recombination

The Holliday junction (HJ) is an essential element in recombination and related mechanisms. The structure of this four-stranded DNA assembly, which is now well-defined alone and in complex with proteins, has led to its applications in areas well outside of molecular recombination, including nanotechnology and biophysics. This minireview explores some interesting recent research on the HJ, as it has been adapted to design regular two- or three-dimensional lattices for crystal engineering, and more complex systems through DNA origami. In addition, the sequence dependence of the structure is discussed in terms how it can be applied to characterize the geometries and energies of various noncovalent interactions, including halogen bonds in oxidatively damaged (halogenated) bases and hydrogen bonds associated with the epigenetic 5-hydroxylmethylcytosine base.

Tree-hierarchy of DNA and distribution of Holliday junctions

We define a DNA as a sequence of [Formula: see text]'s and embed it on a path of Cayley tree. Using group representation of the Cayley tree, we give a hierarchy of a countable set of DNAs each of which 'lives' on the same Cayley tree. This hierarchy has property that each vertex of the Cayley tree belongs only to one of DNA. Then we give a model (energy, Hamiltonian) of this set of DNAs by an analogue of Ising model with three spin values (considered as DNA base pairs) on a set of admissible configurations. To study thermodynamic properties of the model of DNAs we describe corresponding translation invariant Gibbs measures (TIGM) of the model on the Cayley tree of order two. We show that there is a critical temperature [Formula: see text] such that (i) if temperature [Formula: see text] then there exists unique TIGM; (ii) if [Formula: see text] then there are two TIGMs; (iii) if [Formula: see text] then there are three TIGMs. Each such measure describes a phase of the set of DNAs. We use these results to study distributions of Holliday junctions and branches of DNAs. In case of very high and very low temperatures we give stationary distributions and typical configurations of the Holliday junctions.

Effect of Hydroxymethylcytosine on the Structure and Stability of Holliday Junctions

5-Hydroxymethylcytosine (5hmC) is an epigenetic marker that has recently been shown to promote homologous recombination (HR). In this study, we determine the effects of 5hmC on the structure, thermodynamics, and conformational dynamics of the Holliday junction (the four-stranded DNA intermediate associated with HR) in its native stacked-X form. The hydroxymethyl and the control methyl substituents are placed in the context of an amphimorphic GxCC trinucleotide core sequence (where xC is C, 5hmC, or the methylated 5mC), which is part of a sequence also recognized by endonuclease G to promote HR. The hydroxymethyl group of the 5hmC junction adopts two distinct rotational conformations, with an in-base-plane form being dominant over the competing out-of-plane rotamer that has typically been seen in duplex structures. The in-plane rotamer is seen to be stabilized by a more stable intramolecular hydrogen bond to the junction backbone. Stabilizing hydrogen bonds (H-bonds) formed by the hydroxyl substituent in 5hmC or from a bridging water in the 5mC structure provide approximately 1.5-2 kcal/mol per interaction of stability to the junction, which is mostly offset by entropy compensation, thereby leaving the overall stability of the G5hmCC and G5mCC constructs similar to that of the GCC core. Thus, both methyl and hydroxymethyl modifications are accommodated without disrupting the structure or stability of the Holliday junction. Both 5hmC and 5mC are shown to open the structure to make the junction core more accessible. The overall consequences of incorporating 5hmC into a DNA junction are thus discussed in the context of the specificity in protein recognition of the hydroxymethyl substituent through direct and indirect readout mechanisms.

Holliday Junction Thermodynamics and Structure: Coarse-Grained Simulations and Experiments

Holliday junctions play a central role in genetic recombination, DNA repair and other cellular processes. We combine simulations and experiments to evaluate the ability of the 3SPN.2 model, a coarse-grained representation designed to mimic B-DNA, to predict the properties of DNA Holliday junctions. The model reproduces many experimentally determined aspects of junction structure and stability, including the temperature dependence of melting on salt concentration, the bias between open and stacked conformations, the relative populations of conformers at high salt concentration, and the inter-duplex angle (IDA) between arms. We also obtain a close correspondence between the junction structure evaluated by all-atom and coarse-grained simulations. We predict that, for salt concentrations at physiological and higher levels, the populations of the stacked conformers are independent of salt concentration, and directly observe proposed tetrahedral intermediate sub-states implicated in conformational transitions. Our findings demonstrate that the 3SPN.2 model captures junction properties that are inaccessible to all-atom studies, opening the possibility to simulate complex aspects of junction behavior.

Probing the Ion Binding Site in a DNA Holliday Junction Using Förster Resonance Energy Transfer (FRET)

Holliday Junctions are critical DNA intermediates central to double strand break repair and homologous recombination. The junctions can adopt two general forms: open and stacked-X, which are induced by protein or ion binding. In this work, fluorescence spectroscopy, metal ion luminescence and thermodynamic measurements are used to elucidate the ion binding site and the mechanism of junction conformational change. Förster resonance energy transfer measurements of end-labeled junctions monitored junction conformation and ion binding affinity, and reported higher affinities for multi-valent ions. Thermodynamic measurements provided evidence for two classes of binding sites. The higher affinity ion-binding interaction is an enthalpy driven process with an apparent stoichiometry of 2.1 ± 0.2. As revealed by Eu(3+) luminescence, this binding class is homogeneous, and results in slight dehydration of the ion with one direct coordination site to the junction. Luminescence resonance energy transfer experiments confirmed the presence of two ions and indicated they are 6-7 Å apart. These findings are in good agreement with previous molecular dynamics simulations, which identified two symmetrical regions of high ion density in the center of stacked junctions. These results support a model in which site-specific binding of two ions in close proximity is required for folding of DNA Holliday junctions into the stacked-X conformation.

Holliday Junctions Are Associated with Transposable Element Sequences in the Human Genome

Holliday junctions (HJs) constitute important intermediate structures for many cell functions such as DNA recombination and DNA repair. They derive from a 10-nt degenerate sequence, with a 3-nt core motif. In this study, we explored the human genome whether the HJ degenerate sequence associates with transposable elements (TEs) and mainly with those of the active and inactive ALU, LINE, SVA and HERV families. We identified six different forms of the HJ sequence motif, and we located the genomic coordinates of sequences containing both HJs and TEs. From 2982 total HJs, a significant number of 1319 TE-associated HJs were found, with a median distribution of 1 per 2.4 Mb. The HJs with higher GC content were observed more frequently at the genome. A high percentage of HJs were associated with all main TE families, with specificity for particular active or inactive elements: DNA elements and the retroelements ALUs, LINEs and HERVs up to 41.94%, 72.72%, 42.94% and 84.5%, respectively. Phylogenetic analysis revealed that HJs occur in both active and inactive TEs. Furthermore, the TE-associated HJs were almost exclusively found within a distance less than 1 Mb from human genes, while only 23 were not associated with any genes. This is the first report associating human HJs, with mobile elements. Our data pinpoint that particular HJ forms show preference for specific active retrotransposon families of ALUs and LINEs, suggesting that retrotransposon-incorporated HJs may relocate or replicate in the genome through retrotransposition, contributing to recombination, genome plasticity and DNA repair.

Rational design of DNA-actuated enzyme nanoreactors guided by single molecule analysis

The control of enzymatic reactions using nanoscale DNA devices offers a powerful application of DNA nanotechnology uniquely derived from actuation. However, previous characterization of enzymatic reaction rates using bulk biochemical assays reported suboptimal function of DNA devices such as tweezers. To gain mechanistic insight into this deficiency and to identify design rules to improve their function, here we exploit the synergy of single molecule imaging and computational modeling to characterize the three-dimensional structures and catalytic functions of DNA tweezer-actuated nanoreactors. Our analysis revealed two important deficiencies--incomplete closure upon actuation and conformational heterogeneity. Upon rational redesign of the Holliday junctions located at their hinge and arms, we found that the DNA tweezers could be more completely and uniformly closed. A novel single molecule enzyme assay was developed to demonstrate that our design improvements yield significant, independent enhancements in the fraction of active enzyme nanoreactors and their individual substrate turnover frequencies. The sequence-level design strategies explored here may aid more broadly in improving the performance of DNA-based nanodevices including biological and chemical sensors.

Computing Nonequilibrium Conformational Dynamics of Structured Nucleic Acid Assemblies

Synthetic nucleic acids can be programmed to form precise three-dimensional structures on the nanometer-scale. These thermodynamically stable complexes can serve as structural scaffolds to spatially organize functional molecules including multiple enzymes, chromophores, and force-sensing elements with internal dynamics that include substrate reaction-diffusion, excitonic energy transfer, and force-displacement response that often depend critically on both the local and global conformational dynamics of the nucleic acid assembly. However, high molecular weight assemblies exhibit long time-scale and large length-scale motions that cannot easily be sampled using all-atom computational procedures such as molecular dynamics. As an alternative, here we present a computational framework to compute the overdamped conformational dynamics of structured nucleic acid assemblies and apply it to a DNA-based tweezer, a nine-layer DNA origami ring, and a pointer-shaped DNA origami object, which consist of 204, 3,600, and over 7,000 basepairs, respectively. The framework employs a mechanical finite element model for the DNA nanostructure combined with an implicit solvent model to either simulate the Brownian dynamics of the assembly or alternatively compute its Brownian modes. Computational results are compared with an all-atom molecular dynamics simulation of the DNA-based tweezer. Several hundred microseconds of Brownian dynamics are simulated for the nine-layer ring origami object to reveal its long time-scale conformational dynamics, and the first ten Brownian modes of the pointer-shaped structure are predicted.

Applications of halogen bonding in solution

Halogen bonding is the noncovalent interaction where the halogen atom acts as an electrophile towards Lewis bases. Known for more than 200 years, only recently it has attracted interest in the context of solution-phase applications, especially during the last decade which was marked by the introduction of multitopic systems. In addition, the small yet rich collection of halogen-bond donor moieties that appeared in this period is shown to be versatile enough as to be applied in virtually any solvent system. This review covers the applications of halogen bonding in solution during the past ten years in a semi-comprehensive way. Emphasis is made on molecular recognition, catalytic applications and anion binding and transport. Medicinal applications are addressed as well with key examples. Focussing on the major differences observed for halogen bonding, as compared to the ubiquitous hydrogen bonding, it aims to contribute to the design of future solution-phase applications.

Lattice-free prediction of three-dimensional structure of programmed DNA assemblies

DNA can be programmed to self-assemble into high molecular weight 3D assemblies with precise nanometer-scale structural features. Although numerous sequence design strategies exist to realize these assemblies in solution, there is currently no computational framework to predict their 3D structures on the basis of programmed underlying multi-way junction topologies constrained by DNA duplexes. Here, we introduce such an approach and apply it to assemblies designed using the canonical immobile four-way junction. The procedure is used to predict the 3D structure of high molecular weight planar and spherical ring-like origami objects, a tile-based sheet-like ribbon, and a 3D crystalline tensegrity motif, in quantitative agreement with experiments. Our framework provides a new approach to predict programmed nucleic acid 3D structure on the basis of prescribed secondary structure motifs, with possible application to the design of such assemblies for use in biomolecular and materials science.

DNA may be used to fabricate functional nanostructures with various possible geometries, but first being able to predict these structures is a challenging task. Here, the authors use coarse-grained modelling to predict the shape of artificial DNA nanostructures in solution.

Making the bend: DNA tertiary structure and protein-DNA interactions

DNA structure functions as an overlapping code to the DNA sequence. Rapid progress in understanding the role of DNA structure in gene regulation, DNA damage recognition and genome stability has been made. The three dimensional structure of both proteins and DNA plays a crucial role for their specific interaction, and proteins can recognise the chemical signature of DNA sequence ("base readout") as well as the intrinsic DNA structure ("shape recognition"). These recognition mechanisms do not exist in isolation but, depending on the individual interaction partners, are combined to various extents. Driving force for the interaction between protein and DNA remain the unique thermodynamics of each individual DNA-protein pair. In this review we focus on the structures and conformations adopted by DNA, both influenced by and influencing the specific interaction with the corresponding protein binding partner, as well as their underlying thermodynamics.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

Crystallographic and spectroscopic studies of d(CCGGTACCGG)

The decanucleotide sequence d(CCGGTACCGG) crystallizes as a four-way junction at low cobalt ion concentrations (i.e., 1 mM). When the cobalt concentration in the crystallization solution is increased to 5 mM, the sequence crystallizes as resolved B-DNA duplexes. Gel retardation studies of the decamer show both a faint slow-moving band and a much thicker fast-moving band at low cobalt ion concentrations, and only the intense fast-moving band at higher ion concentration. Circular dichroism (CD) spectroscopy of the decamer indicates a structural transition as the cobalt ion concentration in the solution is increased, probably from B-type to A-type DNA. These studies revealed that the oligomer sequence has several conformations and structures accessible to it, in a manner dependent on sequence, ion concentration, and DNA concentration. [Supplementary materials are available for this article. Go to the publisher's online edition of Nucleosides, Nucleotides & Nucleic Acids for the following free supplemental resources(s): Supplementary Figures 1, 2, and 3.].

Measurements of DNA-loop formation via Cre-mediated recombination

The Cre-recombination system has become an important tool for genetic manipulation of higher organisms and a model for site-specific DNA-recombination mechanisms employed by the λ-Int superfamily of recombinases. We report a novel quantitative approach for characterizing the probability of DNA-loop formation in solution using time-dependent ensemble Förster resonance energy transfer measurements of intra- and inter-molecular Cre-recombination kinetics. Our method uses an innovative technique for incorporating multiple covalent modifications at specific sites in covalently closed DNA. Because the mechanism of Cre recombinase does not conform to a simple kinetic scheme, we employ numerical methods to extract rate constants for fundamental steps that pertain to Cre-mediated loop closure. Cre recombination does not require accessory proteins, DNA supercoiling or particular metal-ion cofactors and is thus a highly flexible system for quantitatively analyzing DNA-loop formation in vitro and in vivo.

Molecular dynamics of a DNA Holliday junction: the inverted repeat sequence d(CCGGTACCGG)₄

All-atom molecular dynamics (MD) computer simulations have been applied successfully to duplex DNA structures in solution for some years and found to give close accord with observed results. However, the MD force fields have generally not been parameterized against unusual DNA structures, and their use to obtain dynamical models for this class of systems needs to be independently validated. The four-way junction (4WJ), or Holliday junction, is a dynamic DNA structure involved in central cellular processes of homologous replication and double strand break repair. Two conformations are observed in solution: a planar open-X form (OPN) with a mobile center and four duplex arms, and an immobile stacked-X (STX) form with two continuous strands and two crossover strands, stabilized by high salt conditions. To characterize the accuracy of MD modeling on 4WJ, we report a set of explicit solvent MD simulations of ∼100 ns on the repeat sequence d(CCGGTACCGG)(4) starting from the STX structure (PDB code 1dcw), and an OPN structure built for the same sequence. All 4WJ MD simulations converged to a stable STX structure in close accord with the crystal structure. Our MD beginning in the OPN form converts to the STX form spontaneously at both high and low salt conditions, providing a model for the conformational transition. Thus, these simulations provide a successful account of the dynamical structure of the STX form of d(CCGGTACCGG)(4) in solution, and provide new, to our knowledge, information on the conformational stability of the junction and distribution of counterions in the junction interior.

Structure of d(CGGGTACCCG)4 as a four-way Holliday junction

The crystal structure of the decamer sequence d(CGGGTACCCG)(4) as a four-way Holliday junction has been determined at 2.35 Å resolution. The sequence was designed in order to understand the principles that govern the relationship between sequence and branching structure. It crystallized as a four-way junction structure with an overall geometry similar to those of previously determined Holliday junction structures.

The sequence d(CGGCGGCCGC) self-assembles into a two dimensional rhombic DNA lattice

We report here the crystal structure of the partially self-complementary decameric sequence d(CGGCGGCCGC), which self assembles to form a four-way junction with sticky ends. Each junction binds to four others through Watson-Crick base pairing at the sticky ends to form a rhombic structure. The rhombuses bind to each other and form two dimensional tiles. The tiles stack to form the crystal. The crystal diffracted in the space group P1 to a resolution of 2.5Å. The junction has the anti-parallel stacked-X conformation like other junction structures, though the formation of the rhombic net noticeably alters the details of the junction geometry.

HU binding to a DNA four-way junction probed by Förster resonance energy transfer

The Escherichia coli protein HU is a non-sequence-specific DNA-binding protein that interacts with DNA primarily through electrostatic interactions. In addition to nonspecific binding to linear DNA, HU has been shown to bind with nanomolar affinity to discontinuous DNA substrates, such as repair and recombination intermediates. This work specifically examines the HU-four-way junction (4WJ) interaction using fluorescence spectroscopic methods. The conformation of the junction in the presence of different counterions was investigated by Förster resonance energy transfer (FRET) measurements, which revealed an ion-type conformational dependence, where Na(+) yields the most stacked conformation followed by K(+) and Mg(2+). HU binding induces a greater degree of stacking in the Na(+)-stabilized and Mg(2+)-stabilized junctions but not the K(+)-stabilized junction, which is attributed to differences in the size of the ionic radii and potential differences in ion binding sites. Interestingly, junction conformation modulates binding affinity, where HU exhibits the lowest affinity for the Mg(2+)-stabilized form (24 μM(-1)), which is the least stacked conformation. Protein binding to a mixed population of open and stacked forms of the junction leads to nearly complete formation of a protein-stabilized stacked-X junction. These results strongly support a model in which HU binds to and stabilizes the stacked-X conformation.

Non-B DNA structure-induced genetic instability and evolution

Repetitive DNA motifs are abundant in the genomes of various species and have the capacity to adopt non-canonical (i.e., non-B) DNA structures. Several non-B DNA structures, including cruciforms, slipped structures, triplexes, G-quadruplexes, and Z-DNA, have been shown to cause mutations, such as deletions, expansions, and translocations in both prokaryotes and eukaryotes. Their distributions in genomes are not random and often co-localize with sites of chromosomal breakage associated with genetic diseases. Current genome-wide sequence analyses suggest that the genomic instabilities induced by non-B DNA structure-forming sequences not only result in predisposition to disease, but also contribute to rapid evolutionary changes, particularly in genes associated with development and regulatory functions. In this review, we describe the occurrence of non-B DNA-forming sequences in various species, the classes of genes enriched in non-B DNA-forming sequences, and recent mechanistic studies on DNA structure-induced genomic instability to highlight their importance in genomes.

A rare nucleotide base tautomer in the structure of an asymmetric DNA junction

The single-crystal structure of a DNA Holliday junction assembled from four unique sequences shows a structure that conforms to the general features of models derived from similar constructs in solution. The structure is a compact stacked-X form junction with two sets of stacked B-DNA-type arms that coaxially stack to form semicontinuous duplexes interrupted only by the crossing of the junction. These semicontinuous helices are related by a right-handed rotation angle of 56.5 degrees, which is nearly identical to the 60 degree angle in the solution model but differs from the more shallow value of approximately 40 degrees for previous crystal structures of symmetric junctions that self-assemble from single identical inverted-repeat sequences. This supports the model in which the unique set of intramolecular interactions at the trinucleotide core of the crossing strands, which are not present in the current asymmetric junction, affects both the stability and geometry of the symmetric junctions. An unexpected result, however, is that a highly wobbled A.T base pair, which is ascribed here to a rare enol tautomer form of the thymine, was observed at the end of a CCCC/GGGG sequence within the stacked B-DNA arms of this 1.9 A resolution structure. We suggest that the junction itself is not responsible for this unusual conformation but served as a vehicle for the study of this CG-rich sequence as a B-DNA duplex, mimicking the form that would be present in a replication complex. The existence of this unusual base lends credence to and defines a sequence context for the "rare tautomer hypothesis" as a mechanism for inducing transition mutations during DNA replication.

Assembly pathway analysis of DNA nanostructures and the construction of parallel motifs

We present a system for analyzing the assembly pathway of DNA nanostructures. This enables the identification, explanation, and avoidance of obstacles to proper structure formation. Potential problems include strand end-pinning and misfolding caused by the structural bias of nominally flexible junctions. We have used this system to guide the construction of parallel motifs that had previously, for unknown reasons, resisted assembly.

Functional architecture of T7 RNA polymerase transcription complexes

Bacteriophage T7 RNA polymerase is the best-characterized member of a widespread family of single-subunit RNA polymerases. Crystal structures of T7 RNA polymerase initiation and elongation complexes have provided a wealth of detailed information on RNA polymerase interactions with the promoter and transcription bubble, but the absence of DNA downstream of the melted region of the template in the initiation complex structure, and the absence of DNA upstream of the transcription bubble in the elongation complex structure means that our picture of the functional architecture of T7 RNA polymerase transcription complexes remains incomplete. Here, we use the site-specifically tethered chemical nucleases and functional characterization of directed T7 RNAP mutants to both reveal the architecture of the duplex DNA that flanks the transcription bubble in the T7 RNAP initiation and elongation complexes, and to define the function of the interactions made by these duplex elements. We find that downstream duplex interactions made with a cluster of lysine residues (K711/K713/K714) are present during both elongation and initiation, where they contribute to stabilizing a bend in the downstream DNA that is important for promoter opening. The upstream DNA in the elongation complex is also found to be sharply bent at the upstream edge of the transcription bubble, thereby allowing formation of upstream duplex:polymerase interactions that contribute to elongation complex stability.

Directing macromolecular conformation through halogen bonds

The halogen bond, a noncovalent interaction involving polarizable chlorine, bromine, or iodine molecular substituents, is now being exploited to control the assembly of small molecules in the design of supramolecular complexes and new materials. We demonstrate that a halogen bond formed between a brominated uracil and phosphate oxygen can be engineered to direct the conformation of a biological molecule, in this case to define the conformational isomer of a four-stranded DNA junction when placed in direct competition against a classic hydrogen bond. As a result, this bromine interaction is estimated to be approximately 2-5 kcal/mol stronger than the analogous hydrogen bond in this environment, depending on the geometry of the halogen bond. This study helps to establish halogen bonding as a potential tool for the rational design and construction of molecular materials with DNA and other biological macromolecules.

The stacked-X DNA Holliday junction and protein recognition

The crystal structure of the four-stranded DNA Holliday junction has now been determined in the presence and absence of junction binding proteins, with the extended open-X form of the junction seen in all protein complexes, but the more compact stacked-X structure observed in free DNA. The structures of the stacked-X junction were crystallized because of an unexpected sequence dependence on the stability of this structure. Inverted repeat sequences that contain the general motif NCC or ANC favor formation of stacked-X junctions, with the junction cross-over occurring between the first two positions of the trinucleotides. This review focuses on the sequence dependent structure of the stacked-X junction and how it may play a role in structural recognition by a class of dimeric junction resolving enzymes that themselves show no direct sequence recognition.

Solution formation of Holliday junctions in inverted-repeat DNA sequences

The structure of Holliday junctions has now been well characterized at the atomic level through single-crystal X-ray diffraction in symmetric (inverted-repeat) DNA sequences. At issue, however, is whether the formation of these four-stranded complexes in solution is truly sequence dependent in the manner proposed or is an artifact of the crystallization process and, therefore, has no relevance to the behavior of this central intermediate in homologous recombination and recombination-dependent cellular processes. Here, we apply analytical ultracentrifugation to demonstrate that the sequence d(CCGGTACCGG), which crystallizes in the stacked-X form of the junction, assembles into four-stranded junctions in solution in a manner that is dependent on the DNA and cation concentrations, with an equilibrium established between the junction and duplex forms at 100-200 microM DNA duplex. In contrast, the sequence d(CCGCTAGCGG), which has been crystallized as B-DNA, is seen to adopt only the double-helical form at all DNA and salt concentrations that were tested. Thus, the ACC trinucleotide core is now shown to be important for the formation of Holliday junctions in both crystals and in solution and can be estimated to contribute approximately -4 kcal/mol to stabilizing this recombination intermediate in inverted-repeat sequences.

Conformational model of the Holliday junction transition deduced from molecular dynamics simulations

Homologous recombination plays a key role in the restart of stalled replication forks and in the generation of genetic diversity. During this process, two homologous DNA molecules undergo strand exchange to form a four-way DNA (Holliday) junction. In the presence of metal ions, the Holliday junction folds into the stacked-X structure that has two alternative conformers. Experiments have revealed the spontaneous transitions between these conformers, but their detailed pathways are not known. Here, we report a series of molecular dynamics simulations of the Holliday junction at physiological and elevated (400 K) temperatures. The simulations reveal new tetrahedral intermediates and suggest a schematic framework for conformer transitions. The tetrahedral intermediates bear resemblance to the junction conformation in complex with a junction-resolving enzyme, T7 endonuclease I, and indeed, one intermediate forms a stable complex with the enzyme as demonstrated in one simulation. We also describe free energy minima for various states of the Holliday junction system, which arise during conformer transitions. The results show that magnesium ions stabilize the stacked-X form and destabilize the open and tetrahedral intermediates. Overall, our study provides a detailed dynamic model of the Holliday junction undergoing a conformer transition.

Influence of minor groove substituents on the structure of DNA Holliday junctions

The inosine-containing sequence d(CCIGTACm(5)CGG) is shown to crystallize as a four-stranded DNA junction. This structure is nearly identical to the antiparallel junction formed by the parent d(CCGGTACm(5)()CGG) sequence [Vargason, J. M., and Ho, P. S. (2002) J. Biol. Chem. 277, 21041-21049] in terms of its conformational geometry, and inter- and intramolecular interactions within the DNA and between the DNA and solvent, even though the 2-amino group in the minor groove of the important G(3).m(5)C(8) base pair of the junction core trinucleotide (italicized) has been removed. In contrast, the analogous 2,6-diaminopurine sequence d(CCDGTACTGG) crystallizes as resolved duplex DNAs, just like its parent sequence d(CCAGTACTGG) [Hays, F. A., Vargason, J. M., and Ho, P. S. (2003) Biochemistry 42, 9586-9597]. These results demonstrate that it is not the presence or absence of the 2-amino group in the minor groove of the R(3).Y(8) base pair that specifies whether a sequence forms a junction, but the positions of the extracyclic amino and keto groups in the major groove. Finally, the study shows that the arms of the junction can accommodate perturbations to the B-DNA conformation of the stacked duplex arms associated with the loss of the 2-amino substituent, and that two hydrogen bonding interactions from the C(7) and Y(8) pyrimidine nucleotides to phosphate oxygens of the junction crossover specify the geometry of the Holliday junction.