Design of immobile nucleic acid junctions

Nucleic acids that interact to generate structures in which three or more double helices emanate from a single point are said to form a junction. Such structures arise naturally as intermediates in DNA replication and recombination. It has been proposed that stable junctions can be created by synthesizing sets of oligonucleotides of defined sequence that can associate by maximizing Watson-Crick complementarity (Seeman N. C., 1981, Biomolecular Stereodynamics. Adenine Press, New York. 1: 269-278; Seeman, N. C., 1982, J. Theor. Biol. 99:237-247.) To make it possible to design molecules that will form junctions of specific architecture, we present here an efficient algorithm for generating nucleic acid sequences that optimize two fundamental properties: fidelity and stability. Fidelity refers to the relative probability of forming the junction complex relative to all alternative paired structures. Calculations are described that permit approximate prediction of the melting curves for junction complexes.     

Macromolecular design, nucleic acid junctions, and crystal formation

The simplest form of macromolecular design involves the ligation of nucleic acids. Recent results on the concatenation of nucleic acid junctions show that these molecules can act as fairly rigid macromolecular valence clusters on the nanometer scale. These clusters can be joined to form closed stick figures in which each edge is double helical DNA or RNA and each vertex is a nucleic acid junction. The geometrical criteria for forming discrete-closed and periodic structures from these components are established. The helicity of each edge limits the possible structures that can be formed. The formation of a periodic array from nucleic acid junction building blocks is compared with the crystallization of molecular systems. This comparison leads to a new interpretation of the nature of order in the solid state for molecular crystals. The suggestion is made that the structure of a solid molecular system described by the fewest unique orthogonal (Fourier) components is the one which will be entropically favored, since it contains the least information. This is the crystalline state, with a small number of molecules per asymmetric unit. The free energy from the proposed entropic driving force responsible for this behavior is available, in principle, to correct small deviations from ideality in forming covalent crystals from nucleic acid junction components, as well as in non-bonded molecular systems. Nucleic acid junction periodic arrays provide an appropriate vehicle with which to test this interpretation.     

Macromolecular Design,Nucleic Acid Junctions,and Crystal Formation

Abstract

The simplest form of macromolecular design involves the ligation of nucleic acids. Recent results on the concatenation of nucleic acid junctions show that these molecules can act as fairly rigid macromolecular valence clusters on the nanometer scale. These clusters can be joined to form closed stick figures in which each edge is double helical DNA or RNA and each vertex is a nucleic acid junction. The geometrical criteria for forming discrete-closed and periodic structures from these components are established. The helicity of each edge limits the possible structures that can be formed.

The formation of a periodic array from nucleic acid junction building blocks is compared with the crystallization of molecular systems. This comparison leads to a new interpretation of the nature of order in the solid state for molecular crystals. The suggestion is made that the structure of a solid molecular system described by the fewest unique orthogonal (Fourier) components is the one which will be entropically favored, since it contains the least information. This is the crystalline state, with a small number of molecules per asymmetric unit. The free energy from the proposed entropie driving force responsible for this behavior is available, in principle, to correct small deviations from ideality in forming covalent crystals from nucleic acid junction components, as well as in non-bonded molecular systems. Nucleic acid junction periodic arrays provide an appropriate vehicle with which to test this interpretation.     

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.     

Nucleic acid structure and recognition

We review the global structures adopted by branched nucleic acids, including three- and four-way helical junctions in DNA and RNA. We find that some general folding principles emerge. First, all the structures exhibit a tendency to undergo pairwise coaxial helical stacking when permitted by the local stereochemistry of strand exchange. Second, metal ions generally play an important role in facilitating folding of branched nucleic acids. These principles can be applied to functionally important branched nucleic acids, such as the Holliday DNA junction of genetic recombination, and the hammerhead ribozyme in RNA.     

Computational approaches toward the design of pools for the in vitro selection of complex aptamers

It is well known that using random RNA/DNA sequences for SELEX experiments will generally yield low-complexity structures. Early experimental results suggest that having a structurally diverse library, which, for instance, includes high-order junctions, may prove useful in finding new functional motifs. Here, we develop two computational methods to generate sequences that exhibit higher structural complexity and can be used to increase the overall structural diversity of initial pools for in vitro selection experiments. Random Filtering selectively increases the number of five-way junctions in RNA/DNA pools, and Genetic Filtering designs RNA/DNA pools to a specified structure distribution, whether uniform or otherwise. We show that using our computationally designed DNA pool greatly improves access to highly complex sequence structures for SELEX experiments (without losing our ability to select for common one-way and two-way junction sequences).     

The design of single-stranded nucleic acid knots

A general strategy is described for the synthesis of single-stranded nucleic acid knots. Control of nucleic acid sequence is used to direct the formation of secondary structures that produce the target topology. The key feature of the strategy is the equation of a half-turn of double helical DNA or RNA with a node in a knot. By forming nodes from complementary DNA sequences, it appears possible to direct the assembly of any simple knot. Stabilization of individual nodes may be achieved by constructing them from long regions containing both B-DNA and Z-DNA. Control over the braiding of DNA that acts as a link between node-forming domains can be realized by condensing the nodes into well-defined DNA structures, such as extended domains of linear duplex, branched junctions, antijunctions or mesojunctions. Further topological control may be derived from the pairing of linker regions to complementary single-stranded molecules, thereby preventing them from braiding in an undesirable fashion.     

Analysis of branched nucleic acid structure using comparative gel electrophoresis

Electrophoresis in polyacrylamide gels provides a simple yet powerful means of analyzing the relative disposition of helical arms in branched nucleic acids. The electrophoretic mobility of DNA or RNA with a central discontinuity is determined by the angle subtended between the arms radiating from the branchpoint. In a multi-helical branchpoint, comparative gel electrophoresis can provide a relative measure of all the inter-helical angles and thus the shape and symmetry of the molecule. Using the long-short arm approach, the electrophoretic mobility of all the species with two helical arms that are longer than all others is compared. This can be done as a function of conditions, allowing the analysis of ion-dependent folding of branched DNA and RNA species. Notable successes for the technique include the four-way (Holliday) junction in DNA and helical junctions in functionally significant RNA species such as ribozymes. Many of these structures have subsequently been proved correct by crystallography or other methods, up to 10 years later in the case of the Holliday junction. Just as important, the technique has not failed to date. Comparative gel electrophoresis can provide a window on both fast and slow conformational equilibria such as conformer exchange in four-way DNA junctions. But perhaps the biggest test of the approach has been to deduce the structures of complexes of four-way DNA junctions with proteins. Two recent crystallographic structures show that the global structures were correctly deduced by electrophoresis, proving the worth of the method even in these rather complex systems. Comparative gel electrophoresis is a robust method for the analysis of branched nucleic acids and their complexes.     

DNA junctions, antijunctions, and mesojunctions.

Antijunctions and mesojunctions are new classes of multistranded DNA complexes. They represent a generalization of DNA branched junction complexes, such as the Holliday recombination intermediate. Each strand of a conventional branched junction participates in two different double helices, and this is also true for mesojunctions and antijunctions. The helix axes of conventional branched junction complexes may be drawn to converge at a point, but this convergence occurs for lines drawn perpendicular to the helix axes of antijunctions. Mesojunctions are complexes that mix these features of junctions and antijunctions. Antijunction complexes require an even number of strands. We have synthesized the mesojunction containing three strands, the two mesojunctions containing four strands, and the antijunction containing four strands; we compare them with branched junctions containing three and four strands, derived by permutations of the same sequences. Each double helix is designed to contain 1.5 turns of DNA. A tendency to oligomerize makes it difficult to capture antijunctions and mesojunctions in stable discrete complexes, in contrast to conventional branched junctions. For both three-strand and four-strand complexes, Tm is highest for conventional branched junctions. Ferguson analysis indicates similarities in the occluded surface area of junctions, antijunctions, and one four-strand mesojunction, but the other four-strand mesojunction has a much lower apparent surface area. Hydroxyl radical cleavage patterns suggest that the four-strand antijunction and the low-surface-area four-strand mesojunction form stacking domains, analogous to the behavior of conventional branched junctions. These new structures are related to replicational and recombinational intermediates and to single-stranded nucleic acid knots.     

Holliday junction binding and processing by the RuvA protein of Mycoplasma pneumoniae.

The RuvA, RuvB and RuvC proteins of Escherichia coli act together to process Holliday junctions formed during recombination and DNA repair. RuvA has a well-defined DNA binding surface that is sculptured specifically to accommodate a Holliday junction and allow subsequent loading of RuvB and RuvC. A negatively charged pin projecting from the centre limits binding of linear duplex DNA. The amino-acid sequences forming the pin are highly conserved. However, in certain Mycoplasma and Ureaplasma species the structure is extended by four amino acids and two acidic residues forming a crucial charge barrier are missing. We investigated the significance of these differences by analysing RuvA from Mycoplasma pneumoniae. Gel retardation and surface plasmon resonance assays revealed that this protein binds Holliday junctions and other branched DNA structures in a manner similar to E. coli RuvA. Significantly, it binds duplex DNA more readily. However it does not support branch migration mediated by E. coli RuvB and when bound to junction DNA is unable to provide a platform for stable binding of E. coli RuvC. It also fails to restore radiation resistance to an E. coli ruvA mutant. The data presented suggest that the modified pin region retains the ability to promote junction-specific DNA binding, but acts as a physical obstacle to linear duplex DNA rather than as a charge barrier. They also indicate that such an obstacle may interfere with the binding of a resolvase. Mycoplasma species may therefore process Holliday junctions via uncoupled branch migration and resolution reactions.     

Specific detection of DNA and RNA targets using a novel isothermal nucleic acid amplification assay based on the formation of a three-way junction structure

The formation of DNA three-way junction (3WJ) structures has been utilised to develop a novel isothermal nucleic acid amplification assay (SMART) for the detection of specific DNA or RNA targets. The assay consists of two oligonucleotide probes that hybridise to a specific target sequence and, only then, to each other forming a 3WJ structure. One probe (template for the RNA signal) contains a non-functional single-stranded T7 RNA polymerase promoter sequence. This promoter sequence is made double-stranded (hence functional) by DNA polymerase, allowing T7 RNA polymerase to generate a target-dependent RNA signal which is measured by an enzyme-linked oligosorbent assay (ELOSA). The sequence of the RNA signal is always the same, regardless of the original target sequence. The SMART assay was successfully tested in model systems with several single-stranded synthetic targets, both DNA and RNA. The assay could also detect specific target sequences in both genomic DNA and total RNA from Escherichia coli. It was also possible to generate signal from E.coli samples without prior extraction of nucleic acid, showing that for some targets, sample purification may not be required. The assay is simple to perform and easily adaptable to different targets.     

The Archaeal Topoisomerase Reverse Gyrase Is a Helix-destabilizing Protein That Unwinds Four-way DNA Junctions

Four-way junctions are non-B DNA structures that originate as intermediates of recombination and repair (Holliday junctions) or from the intrastrand annealing of palindromic sequences (cruciforms). These structures have important functional roles but may also severely interfere with DNA replication and other genetic processes; therefore, they are targeted by regulatory and architectural proteins, and dedicated pathways exist for their removal. Although it is well known that resolution of Holliday junctions occurs either by recombinases or by specialized helicases, less is known on the mechanisms dealing with secondary structures in nucleic acids. Reverse gyrase is a DNA topoisomerase, specific to microorganisms living at high temperatures, which comprises a type IA topoisomerase fused to an SF2 helicase-like module and catalyzes ATP hydrolysis-dependent DNA positive supercoiling. Reverse gyrase is likely involved in regulation of DNA structure and stability and might also participate in the cell response to DNA damage. By applying FRET technology to multiplex fluorophore gel imaging, we show here that reverse gyrase induces unwinding of synthetic four-way junctions as well as forked DNA substrates, following a mechanism independent of both the ATPase and the strand-cutting activity of the enzyme. The reaction requires high temperature and saturating protein concentrations. Our results suggest that reverse gyrase works like an ATP-independent helix-destabilizing protein specific for branched DNA structures. The results are discussed in light of reverse gyrase function and their general relevance for protein-mediated unwinding of complex DNA structures.     

Assembly and characterization of five-arm and six-arm DNA branched junctions.

DNA branched junctions have been constructed that contain either five arms or six arms surrounding a branch point. These junctions are not as stable as junctions containing three or four arms; unlike the smaller junctions, they cannot be shown to migrate as a single band on native gels when each of their arms contains eight nucleotide pairs. However, they can be stabilized if their arms contain 16 nucleotide pairs. Ferguson analysis of these junctions in combination with three-arm and four-arm junctions indicates a linear increase in friction constant as the number of arms increases, with the four-arm junction migrating anomalously. The five-arm junction does not appear to have any unusual stacking structure, and all strands show similar responses to hydroxyl radical autofootprinting analysis. By contrast, one strand of the six-arm junction shows virtually no protection from hydroxyl radicals, suggesting that it is the helical strand of a preferred stacking domain. Both junctions are susceptible to digestion by T4 endonuclease VII, which resolves Holliday junctions. However, the putative helical strand of the six-arm junction shows markedly reduced cleavage, supporting the notion that its structure is largely found in a helical conformation. Branched DNA molecules can be assembled into structures whose helix axes form multiply connected objects and networks. The ability to construct five-arm and six-arm junctions vastly increases the number of structures and networks that can be built from branched DNA components. Icosahedral deltahedra and 11 networks with 432 symmetry, constructed from Platonic and Archimedean solids, are among the structures whose construction is feasible, now that these junctions can be made.     

Cruciform-resolvase interactions in supercoiled DNA

T4 endonuclease VII, which cleaves Holliday-like junctions in DNA, specifically cleaves short inverted repeats in supercoiled plasmids. These sequences are subject to site-specific cleavage by single-strand-specific nucleases, and cruciform formation has been suggested as an explanation for this observation. This proposal is greatly strengthened by the present data, since a formal analogy between cruciform structures and Holliday junctions exists. Resolution of a variety of unrelated cruciform sequences demonstrates that the cleavage process results in a linear molecule with hairpin ends and single ligatable nicks at positions corresponding to the stem-base of the cruciform. In two examples mapped in detail, the cleavages are exclusively introduced at two or three nucleotides from the end of the symmetric sequence at the 5' side on each strand. These studies demonstrate the potential of endonuclease VII as a probe of cruciform structure and the utility of short cruciform structures as Holliday junction models.     

Modeling protein–nucleic acid complexes with extremely large conformational changes using Flex-LZerD

Proteins and nucleic acids are key components in many processes in living cells, and interactions between proteins and nucleic acids are often crucial pathway components. In many cases, large flexibility of proteins as they interact with nucleic acids is key to their function. To understand the mechanisms of these processes, it is necessary to consider the 3D atomic structures of such protein–nucleic acid complexes. When such structures are not yet experimentally determined, protein docking can be used to computationally generate useful structure models. However, such docking has long had the limitation that the consideration of flexibility is usually limited to small movements or to small structures. We previously developed a method of flexible protein docking which could model ordered proteins which undergo large-scale conformational changes, which we also showed was compatible with nucleic acids. Here, we elaborate on the ability of that pipeline, Flex-LZerD, to model specifically interactions between proteins and nucleic acids, and demonstrate that Flex-LZerD can model more interactions and types of conformational change than previously shown.     

Evaluation of polymer sequence fragment data using graph theory

Much of recent work to determine primary structures of nucleic acids and proteins employs the “fragmentation” or “overlap” stratagem. Typically, a preparation of a given polymer with unknown sequence is purified and then subjected to an enzyme known to cut the polymer at certain specific sites. The quantities and sequences of the resulting fragments are determined. For RNA primary sequences, pancreatic ribonuclease and T1 ribonuclease are ordinarily used as fragmenting enzymes. A technique is described for evaluating such fragment data. It has the following properties: It is easily determined whether or not the fragment data is inconsistent. It is always possible to determine the first and last nucleotides of the unknown sequence from the data of two limit digests. Consistent data from two limit digests can always be fitted into a convenient conceptual framework developed within the theory of graphs. In most cases, partial digest information can be used to modify the framework constructed from two limit digests, as such information is obtained. An efficient analysis of all fragment data in this conceptual framework can always be made. One can detect inconsistencies and can generate the entire list of polymer sequences consistent with the fragment data.     

The ligation and flexibility of four-arm DNA junctions

Four-arm DNA branched junctions are stable analogs of Holliday recombination intermediates, constructed from oligonucleotides. The conformational flexibility of junctions can be estimated by ligating them together and determining the set of closed macrocyclic products that are obtained among the linked units. We have performed a series of these experiments, using pairs of sticky ends that flank each of the six angles of a four-arm junction. In every case, the ligated junctions are separated by 20 nucleotide pairs, about two turns of DNA. All expected short linear products, starting with dimers, are observed for all ligations. All ligations result in a macrocyclic series that begins with trimers. Thus, over the time scale of these reactions, the arms of this junction can form angles as low as 60°. The response of this junction to torsional stress has been tested in a companion experiment. A smaller version of this same four-arm junction has been oligomerized so that successive junctions are separated by 16 nucleotide pairs, approximately 1.5 turns of DNA. If junctions were as rigid as linear duplex DNA, this system would not be expected to form macrocycles until the continuous chain approaches the Shore–Baldwin limit, ca. 160 base pairs. However, macrocyclic closure is observed in a regular ligation ladder, starting from tetramers. Model building suggests that the most likely explanation for the observed closure is that the junction adopts two different conformations, which bend the continuous strand toward opposite grooves. The junction structures formed by these ligations represent fluctuations from equilibrium structures.     

Model for the interaction of DNA junctions and resolving enzymes.

Four-way DNA junctions are thought to be important intermediates in a number of recombination processes. Resolution of these junctions occurs by cleavage of two strands of DNA to generate two duplex molecules. The interaction between DNA junctions and resolving enzymes appears to be largely structure-specific, reflecting a molecular recognition on a significant scale. We propose a working model for this interaction that takes account of the present state of knowledge of the structure of the DNA junction, and the substrate requirements of the enzymes. We note that three different enzymes introduce cleavages at phosphodiester bonds that are presented on one side of the molecule, suggesting that the enzymes selectively interact with this face of the junction. By forcing a junction of constant sequence to adopt one or other of the two possible antiparallel isomers, we show that the junction is cleaved in such a way as to suggest a constant mode of interaction with the protein that is dependent on structure rather than sequence. We propose that the feature that is recognized is a mutual inclination of two DNA helices at approximately 120 degrees. We show that a number of DNA substrates that contain similar inclined helices, such as a three-way junction, bulged duplexes and a duplex that is curved because of repeated runs of oligoadenine sequences, are each cleaved by phage T4 endonuclease VII. This mode of DNA-protein interaction could be significant in either recombination or DNA repair processes.     

Diversity of Tn4001 Transposition Products: the Flanking IS256 Elements Can Form Tandem Dimers and IS Circles

We show that both flanking IS256 elements carried by transposon Tn4001 are capable of generating head-to-tail tandem copies and free circular forms, implying that both are active. Our results suggest that the tandem structures arise from dimeric copies of the donor or vector plasmid present in the population by a mechanism in which an IS256 belonging to one Tn4001 copy attacks an IS256 end carried by the second Tn4001 copy. The resulting structures carry abutted left (inverted left repeat [IRL]) and right (inverted right repeat [IRR]) IS256 ends. Examination of the junction sequence suggested that it may form a relatively good promoter capable of driving transposase synthesis in Escherichia coli. This behavior resembles that of an increasing number of bacterial insertion sequences which generate integrative junctions as part of the transposition cycle. Sequence analysis of the IRL-IRR junctions demonstrated that attack of one end by the other is largely oriented (IRL attacks IRR). Our experiments also defined the functional tips of IS256 as the tips predicted from sequence alignments, confirming that the terminal 4 bp at each end are indeed different. The appearance of these multiple plasmid and transposon forms indicates that care should be exercised when Tn4001 is used in transposition mutagenesis. This is especially true when it is used with naturally transformable hosts, such as Streptococcus pneumoniae, in which reconstitution of the donor plasmid may select for higher-order multimers.     

Role of the septate junction in the regulation of paracellular transepithelial flow

A comparison of the distribution of septate junctions in invertebrate epithelia and tight junctions in vertebrate systems suggests that these structures may be functionally analogous. This proposition is supported by the internal design of each junction which constitutes a serial arrangement of structures crossing the intercellular space between cells to effectively provide resistance to the paracellular flow of water and small molecules. We have tested the validity of such an analogy by examining whether the osmotic sensitivity of the septate junctions of planarian epidermis follow the rather striking pattern observed for the junctions of very tight vertebrate epithelia (e.g. toad urinary bladder). It has been found that the septate junctions in this system respond in similar fashion to their vertebrate counterparts, blistering with accumulated fluid when the medium outside the epidermis is made hypertonic with small, water-soluble molecules. We conclude that the two types of junction probably are functionally analogous and that, in each case, this rectified structural response to transepithelial osmotic gradients may be indicative of the role of such structures in the transport function of epithelia.