A molecular model for RecA-promoted strand exchange via parallel triple-stranded helices

Modeling the Homologous Recombination Process: Methods, Successes and Challenges

Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.

Weaving DNA strands: structural insight on ATP hydrolysis in RecA-induced homologous recombination

Homologous recombination is a fundamental process in all living organisms that allows the faithful repair of DNA double strand breaks, through the exchange of DNA strands between homologous regions of the genome. Results of three decades of investigation and recent fruitful observations have unveiled key elements of the reaction mechanism, which proceeds along nucleofilaments of recombinase proteins of the RecA family. Yet, one essential aspect of homologous recombination has largely been overlooked when deciphering the mechanism: while ATP is hydrolyzed in large quantity during the process, how exactly hydrolysis influences the DNA strand exchange reaction at the structural level remains to be elucidated. In this study, we build on a previous geometrical approach that studied the RecA filament variability without bound DNA to examine the putative implication of ATP hydrolysis on the structure, position, and interactions of up to three DNA strands within the RecA nucleofilament. Simulation results on modeled intermediates in the ATP cycle bring important clues about how local distortions in the DNA strand geometries resulting from ATP hydrolysis can aid sequence recognition by promoting local melting of already formed DNA heteroduplex and transient reverse strand exchange in a weaving type of mechanism.

Mechanism of the formation of the RecA-ssDNA nucleoprotein filament structure: a coarse-grained approach

In prokaryotes, the RecA protein catalyzes the repair and strand exchange of double-stranded DNA. RecA binds to single-stranded DNA (ssDNA) and forms a presynaptic complex in which the protein polymerizes around the ssDNA to form a right-handed helical nucleoprotein filament structure. In the present work, the mechanism for the formation of the RecA-ssDNA filament structure is modeled using coarse-grained molecular dynamics simulations. Information from the X-ray structure was used to model the protein itself but not its interactions; the interactions between the protein and the ssDNA were modeled solely by electrostatic, aromatic, and repulsive energies. For the present study, the monomeric, dimeric, and trimeric units of RecA and 4, 8, and 11 NT-long ssDNA, respectively, were studied. Our results indicate that monomeric RecA is not sufficient for nucleoprotein filament formation; rather, dimeric RecA is the elementary binding unit, with higher multimeric units of RecA facilitating filament formation. Our results reveal that loop region flexibility at the primary binding site of RecA is essential for it to bind the incoming ssDNA, that the aromatic residues present in the loop region play an important role in ssDNA binding, and that ATP may play a role in guiding the ssDNA by changing the electrostatic potential of the RecA protein.

Integrating multi-scale data on homologous recombination into a new recognition mechanism based on simulations of the RecA-ssDNA/dsDNA structure

RecA protein is the prototypical recombinase. Members of the recombinase family can accurately repair double strand breaks in DNA. They also provide crucial links between pairs of sister chromatids in eukaryotic meiosis. A very broad outline of how these proteins align homologous sequences and promote DNA strand exchange has long been known, as are the crystal structures of the RecA-DNA pre- and postsynaptic complexes; however, little is known about the homology searching conformations and the details of how DNA in bacterial genomes is rapidly searched until homologous alignment is achieved. By integrating a physical model of recognition to new modeling work based on docking exploration and molecular dynamics simulation, we present a detailed structure/function model of homology recognition that reconciles extremely quick searching with the efficient and stringent formation of stable strand exchange products and which is consistent with a vast body of previously unexplained experimental results.

Targeting duplex DNA with chimeric α,β-triplex-forming oligonucleotides

Triplex-directed DNA recognition is strictly limited by polypurine sequences. In an attempt to address this problem with synthetic biology tools, we designed a panel of short chimeric α,β-triplex-forming oligonucleotides (TFOs) and studied their interaction with fluorescently labelled duplex hairpins using various techniques. The hybridization of hairpin with an array of chimeric probes suggests that recognition of double-stranded DNA follows complicated rules combining reversed Hoogsteen and non-canonical homologous hydrogen bonding. In the presence of magnesium ions, chimeric TFOs are able to form highly stable α,β-triplexes, as indicated by native gel-electrophoresis, on-array thermal denaturation and fluorescence-quenching experiments. CD spectra of chimeric triplexes exhibited features typically observed for anti-parallel purine triplexes with a GA or GT third strand. The high potential of chimeric α,β-TFOs in targeting double-stranded DNA was demonstrated in the EcoRI endonuclease protection assay. In this paper, we report, for the first time, the recognition of base pair inversions in a duplex by chimeric TFOs containing α-thymidine and α-deoxyguanosine.

RecA homology search is promoted by mechanical stress along the scanned duplex DNA

A RecA–single-stranded DNA (RecA–ssDNA) filament searches a genome for sequence homology by rapidly binding and unbinding double-stranded DNA (dsDNA) until homology is found. We demonstrate that pulling on the opposite termini (3′ and 5′) of one of the two DNA strands in a dsDNA molecule stabilizes the normally unstable binding of that dsDNA to non-homologous RecA–ssDNA filaments, whereas pulling on the two 3′, the two 5′, or all four termini does not. We propose that the ‘outgoing’ strand in the dsDNA is extended by strong DNA–protein contacts, whereas the ‘complementary’ strand is extended by the tension on the base pairs that connect the ‘complementary’ strand to the ‘outgoing’ strand. The stress resulting from different levels of tension on its constitutive strands causes rapid dsDNA unbinding unless sufficient homology is present.

Optimizing the design of oligonucleotides for homology directed gene targeting

Background

Gene targeting depends on the ability of cells to use homologous recombination to integrate exogenous DNA into their own genome. A robust mechanistic model of homologous recombination is necessary to fully exploit gene targeting for therapeutic benefit.

Methodology/Principal Findings

In this work, our recently developed numerical simulation model for homology search is employed to develop rules for the design of oligonucleotides used in gene targeting. A Metropolis Monte-Carlo algorithm is used to predict the pairing dynamics of an oligonucleotide with the target double-stranded DNA. The model calculates the base-alignment between a long, target double-stranded DNA and a probe nucleoprotein filament comprised of homologous recombination proteins (Rad51 or RecA) polymerized on a single strand DNA. In this study, we considered different sizes of oligonucleotides containing 1 or 3 base heterologies with the target; different positions on the probe were tested to investigate the effect of the mismatch position on the pairing dynamics and stability. We show that the optimal design is a compromise between the mean time to reach a perfect alignment between the two molecules and the stability of the complex.

Conclusion and Significance

A single heterology can be placed anywhere without significantly affecting the stability of the triplex. In the case of three consecutive heterologies, our modeling recommends using long oligonucleotides (at least 35 bases) in which the heterologous sequences are positioned at an intermediate position. Oligonucleotides should not contain more than 10% consecutive heterologies to guarantee a stable pairing with the target dsDNA. Theoretical modeling cannot replace experiments, but we believe that our model can considerably accelerate optimization of oligonucleotides for gene therapy by predicting their pairing dynamics with the target dsDNA.

Study of force induced melting of dsDNA as a function of length and conformation

We measure the constant force required to melt double-stranded (ds) DNA as a function of length for lengths from 12 to 100,000 base pairs, where the force is applied to the 3'3' or 5'5' ends of the dsDNA. Molecules with 32 base pairs or fewer melt before overstretching. For these short molecules, the melting force is independent of the ends to which the force is applied and the shear force as a function of length is well described by de Gennes theory with a de Gennes length of less than 10 bp. Molecules with lengths of 500 base pairs or more overstretch before melting. For these long molecules, the melting force depends on the ends to which the force is applied. The melting force as a function of length increases even when the length exceeds 1000 bp, where the length dependence is inconsistent with de Gennes theory. Finally, we expand de Gennes melting theory to 3'5' pulling and compare the predictions with experimental results.

Modeling the early stage of DNA sequence recognition within RecA nucleoprotein filaments

Homologous recombination is a fundamental process enabling the repair of double-strand breaks with a high degree of fidelity. In prokaryotes, it is carried out by RecA nucleofilaments formed on single-stranded DNA (ssDNA). These filaments incorporate genomic sequences that are homologous to the ssDNA and exchange the homologous strands. Due to the highly dynamic character of this process and its rapid propagation along the filament, the sequence recognition and strand exchange mechanism remains unknown at the structural level. The recently published structure of the RecA/DNA filament active for recombination (Chen et al., Mechanism of homologous recombination from the RecA-ssDNA/dsDNA structure, Nature 2008, 453, 489) provides a starting point for new exploration of the system. Here, we investigate the possible geometries of association of the early encounter complex between RecA/ssDNA filament and double-stranded DNA (dsDNA). Due to the huge size of the system and its dense packing, we use a reduced representation for protein and DNA together with state-of-the-art molecular modeling methods, including systematic docking and virtual reality simulations. The results indicate that it is possible for the double-stranded DNA to access the RecA-bound ssDNA while initially retaining its Watson–Crick pairing. They emphasize the importance of RecA L2 loop mobility for both recognition and strand exchange.

Homologous recombination in real time: DNA strand exchange by RecA

Homologous recombination, the exchange of strands between different DNA molecules, is essential for proper maintenance and accurate duplication of the genome. Using magnetic tweezers, we monitor RecA-driven homologous recombination of individual DNA molecules in real time. We resolve several key aspects of DNA structure during and after strand exchange. Changes in DNA length and twist yield helical parameters for the protein-bound three-stranded structure in conditions in which ATP was not hydrolyzed. When strand exchange was completed under ATP hydrolysis conditions that allow protein dissociation, a "D wrap" structure formed. During homologous recombination, strand invasion at one end and RecA dissociation at the other end occurred at the same rate, and our single-molecule analysis indicated that a region of only about 80 bp is actively involved in the synapsis at any time during the entire reaction involving a long ( approximately 1 kb) region of homology.

RecA-DNA filament topology: the overlooked alternative of an unconventional syn-syn duplex intermediate

The helical filaments of RecA protein mediate strand exchange for homologous recombination, but the paths of the interacting DNAs have yet to be determined. Although this interaction is commonly limited to three strands, it is reasoned here that the intrinsic symmetry relationships of quadruplex topology are superior in explaining a range of observations. In particular, this topology suggests the potential of post-exchange base pairing in the unorthodox configuration of syn-syn glycosidic bonds between the nucleotide bases and the pentose rings in the sugar-phosphate backbone, which would transiently be stabilized by the external scaffolding of the RecA protein filament.

Direct evaluation of a kinetic model for RecA-mediated DNA-strand exchange: the importance of nucleic acid dynamics and entropy during homologous genetic recombination

The Escherichia coli RecA protein is the prototype of a class of proteins that play central roles in genomic repair and recombination in all organisms. The unresolved mechanistic strategy by which RecA aligns a single strand of DNA with a duplex DNA and mediates a DNA strand switch is central to understanding homologous recombination. We explored the mechanism of RecA-mediated DNA-strand exchange using oligonucleotide substrates with the intrinsic fluorophore 6-methylisoxanthopterin. Pre-steady-state spectrofluorometric analysis elucidated the earliest transient intermediates formed during recombination and delineated the mechanistic strategy by which RecA facilitates this process. The structural features of the first detectable intermediate and the energetic characteristics of its formation were consistent with interactions between a few bases of the single-stranded DNA and the minor groove of a locally melted or stretched duplex DNA. Further analysis revealed RecA to be an unusual enzyme in that entropic rather than enthalpic contributions dominate its catalytic function, and no unambiguously active role for the protein was detected in the earliest molecular events of recombination. The data best support the conclusion that the mechanistic strategy of RecA likely relies on intrinsic DNA dynamics.

Origins of sequence selectivity in homologous genetic recombination: insights from rapid kinetic probing of RecA-mediated DNA strand exchange

Despite intense effort over the past 30 years, the molecular determinants of sequence selectivity in RecA-mediated homologous recombination have remained elusive. Here, we describe when and how sequence homology is recognized between DNA strands during recombination in the context of a kinetic model for RecA-mediated DNA strand exchange. We characterized the transient intermediates of the reaction using pre-steady-state kinetic analysis of strand exchange using oligonucleotide substrates containing a single fluorescent G analog. We observed that the reaction system was sensitive to heterology between the DNA substrates; however, such a "heterology effect" was not manifest when functional groups were added to or removed from the edges of the base-pairs facing the minor groove of the substrate duplex. Hence, RecA-mediated recombination must occur without the involvement of a triple helix, even as a transient intermediate in the process. The fastest detectable reaction phase was accelerated when the structure or stability of the substrate duplex was perturbed by internal mismatches or the replacement of G.C by I.C base-pairs. These findings indicate that the sequence specificity in recombination is achieved by Watson-Crick pairing in the context of base-pair dynamics inherent to the extended DNA structure bound by RecA during strand exchange.

Construction and evaluation of a kinetic scheme for RecA-mediated DNA strand exchange

The Escherichia coli RecA protein is the prototype of a class of proteins playing a central role in genomic repair and recombination in all organisms. The unresolved mechanistic strategy by which RecA aligns a single strand of DNA with a duplex DNA and mediates a DNA strand switch is central to understanding its recombinational activities. Toward a molecular-level understanding of RecA-mediated DNA strand exchange, we explored its mechanism using oligonucleotide substrates and the intrinsic fluorescence of 6-methylisoxanthopterin (6MI). Steady- and presteady-state spectrofluorometric data demonstrate that the reaction proceeds via a sequential four-step mechanism comprising a rapid, bimolecular association step followed by three slower unimolecular steps. Previous authors have proposed multistep mechanisms involving two or three steps. Careful analysis of the differences among the experimental systems revealed a previously undiscovered intermediate (N1) whose formation may be crucial in the kinetic discrimination of homologous and heterologous sequences. This observation has important implications for probing the fastest events in DNA strand exchange using 6MI to further elucidate the molecular mechanisms of recombination and recombinational repair.

Electromagnetic field effect on separation of nucleotide sequences and unwinding of a double helix during DNA replication

Our previous pre-clinic experimental results have showed that the bacterium infection can be suppressed and the epithelialization can be enhanced by the externally applied rectangular pulsed electrical current stimulation (RPECS). The results are clinically significant for patients, especially for those difficult patients whose skin wounds need long periods to heal. However, the results also raise questions: How does the RPECS accelerate the epithelium cell proliferation? What is the relationship among the bacterium infection, the epithelialization and the RPECS? To answer these questions, we have previously modeled mitosis and cytokinesis mechanisms for animal cells and amitosis for bacteria at a cellular level and in a view of physics. In this paper, we model the separation of nucleotide sequences and the unwinding of a double helix during DNA replication at a molecular level and also in wild types of cells. Firstly, we define a new concept of nucleotide (NT) electromagnetic field (EMF) box (sequence) which carries genetic information: The continuous NT EMF boxes compose a nucleotide strand. Then, we hypothesize the symmetry, repulsion and attraction of NT EMF boxes: If a pair of NT EMF boxes are (quasi) mirror or complementary symmetric about a plane (curve) or point, they repulse or attract from each other because there is a repulsive or attractive EMF force between them. Our models suggest, the repulsive EMF force from children DNA strands simultaneously separates the children DNA strands, splits the hydrogen bonds of parental base pairs, and unwinds the parental double helix while DNA polymerases are synchronously synthesizing the new children DNA strands. To understand the mechanism of epithelialization enhanced with the externally applied RPECS at a molecular level, we hypothesize that the normal separation of nucleotide sequences and unwinding of a double helix during DNA replication could be suppressed in the bacteria but not in the epithelium cells because: (a) the spontaneous EMF in the epithelium could be 1000 times stronger than that in bacteria; (b) the epithelium cells have one more non-conducting envelope (nuclear membrane) to protect the normal separation and unwinding; (c) based on our previous experimental data, the RPECS amount received by the bacteria are three times as much as the amount the epithelium cells receive. Therefore, the epithelium cellular proliferation may be directly, as well as indirectly (e.g., somatic reflex) accelerated by the RPECS.

Numerical investigation of sequence dependence in homologous recognition: evidence for homology traps

During the initial phase of RecA-mediated recombination, known as the search for homology, a single-stranded DNA coated by RecA protein and a homologous double-stranded DNA have to perfectly align and pair. We designed a model for the homology search between short molecules, and performed Monte Carlo Metropolis computer simulations of the process. The central features of our model are 1), the assumption that duplex DNA longitudinal thermal fluctuations are instrumental in the binding; and 2), the explicit consideration of the nucleotide sequence. According to our results, recognition undergoes a first slow nucleation step over a few basepairs, followed by a quick extension of the pairing to adjacent bases. The formation of the three-stranded complex tends to be curbed by heterologies but also by another possible obstacle: the presence of partially homologous stretches, such as mono- or polynucleotide repeats. Actually, repeated sequences are observed to trap the molecules in unproductive configurations. We investigate the dependence of the phenomenon on various energy parameters. This mechanism of homology trapping could have a strong biological relevance in the light of the genomic instability experimentally known to be triggered by repeated sequences.

Probing protein mechanics: residue-level properties and their use in defining domains

It is becoming clear that, in addition to structural properties, the mechanical properties of proteins can play an important role in their biological activity. It nevertheless remains difficult to probe these properties experimentally. Whereas single-molecule experiments give access to overall mechanical behavior, notably the impact of end-to-end stretching, it is currently impossible to directly obtain data on more local properties. We propose a theoretical method for probing the mechanical properties of protein structures at the single-amino acid level. This approach can be applied to both all-atom and simplified protein representations. The probing leads to force constants for local deformations and to deformation vectors indicating the paths of least mechanical resistance. It also reveals the mechanical coupling that exists between residues. Results obtained for a variety of proteins show that the calculated force constants vary over a wide range. An analysis of the induced deformations provides information that is distinct from that obtained with measures of atomic fluctuations and is more easily linked to residue-level properties than normal mode analyses or dynamic trajectories. It is also shown that the mechanical information obtained by residue-level probing opens a new route for defining so-called dynamical domains within protein structures.

Direct evaluation of a mechanism for activation of the RecA nucleoprotein filament

The RecA protein of Escherichia coli controls the SOS response for DNA damage tolerance and plays a crucial role in recombinational DNA repair. The formation of a RecA.ATP.ssDNA complex initiates all RecA activities, and yet this process is not understood at the molecular level. An analysis of RecA.DNA interactions was performed using both a mutant RecA protein containing a tryptophan (Trp) reporter and oligodeoxyribonucleotides (ODNs) containing a fluorescent guanine analogue, 6-methylisoxanthopterin (6MI). Experiments using fluorescent ODNs allowed structurally distinct nucleoprotein filaments, formed in the absence and presence of ATPgammaS (a slowly hydrolyzed analogue of ATP), to be differentiated directly. Stopped-flow spectrofluorometry, combined with presteady-state kinetic analyses, revealed unexpected differences in the rates of RecA.ODN and RecA.ATPgammaS.ODN complex assembly. This is the first demonstration that such intrinsically fluorescent synthetic DNAs can be used to characterize definitively the real-time assembly and activation of RecA.ssDNA complexes. Surprisingly, the ssDNA binding event is almost 50-fold slower in the presence of the activating ATPgammaS cofactor. Furthermore, a combination of time-dependent emission changes from 6MI and Trp allowed the first direct chemical test of whether an inactive filament can isomerize to the active state. The results revealed that, unlike the hexameric motor proteins, the inactive RecA filament cannot directly convert to the active state upon ATPgammaS binding. These results have implications for understanding how a coincidence of functions--an ATP-communicated signal-like activity and an ATP-driven motorlike activity--are resolved within a single protein molecule.

The synaptic complex of RecA protein participates in hybridization and inverse strand exchange reactions

RecA protein catalyzes strand exchange between homologous single-stranded and double-stranded DNAs. In the presence of ATPgammaS, the post-strand exchange synaptic complex is a stable end product that can be studied. Here we ask whether such complexes can hybridize to or exchange with DNA, 2'-OMe RNA, PNA, or LNA oligonucleotides. Using a gel mobility shift assay, we show that the displaced strand of a 45 bp synaptic complex can hybridize to complementary oligonucleotides with different backbones to form a four-stranded (double D-loop) joint that survives removal of the RecA protein. This hybridization reaction, which confirms the single-stranded character of the displaced strand in a synaptic complex, might initiate recombination-dependent DNA replication if it occurs in vivo. We also show that either strand of the heteroduplex in a 30 bp synaptic complex can be replaced with a homologous DNA oligonucleotide in a strand exchange reaction that is mediated by the RecA filament. Consistent with the important role that deoxyribose plays in strand exchange, oligonucleotides with non-DNA backbones did not participate in this reaction. The hybridization and strand exchange reactions reported here demonstrate that short synaptic complexes are dynamic structures even in the presence of ATPgammaS.

Elucidating a key intermediate in homologous DNA strand exchange: structural characterization of the RecA-triple-stranded DNA complex using fluorescence resonance energy transfer

The RecA protein of Escherichia coli plays essential roles in homologous recombination and restarting stalled DNA replication forks. In vitro, the protein mediates DNA strand exchange between single-stranded (ssDNA) and homologous double-stranded DNA (dsDNA) molecules that serves as a model system for the in vivo processes. To date, no high-resolution structure of the key intermediate, comprised of three DNA strands simultaneously bound to a RecA filament (RecA-tsDNA complex), has been reported. We present a systematic characterization of the helical geometries of the three DNA strands of the RecA-tsDNA complex using fluorescence resonance energy transfer (FRET) under physiologically relevant solution conditions. FRET donor and acceptor dyes were used to label different DNA strands, and the interfluorophore distances were inferred from energy transfer efficiencies measured as a function of the base-pair separation between the two dyes. The energy transfer efficiencies were first measured on a control RecA-dsDNA complex, and the calculated helical parameters (h approximately 5 A, Omega(h) approximately 20 degrees ) were consistent with structural conclusions derived from electron microscopy (EM) and other classic biochemical methods. Measurements of the helical parameters for the RecA-tsDNA complex revealed that all three DNA strands adopt extended and unwound conformations similar to those of RecA-bound dsDNA. The structural data are consistent with the hypothesis that this complex is a late, post-strand-exchange intermediate with the outgoing strand shifted by about three base-pairs with respect to its registry with the incoming and complementary strands. Furthermore, the bases of the incoming and complementary strands are displaced away from the helix axis toward the minor groove of the heteroduplex, and the bases of the outgoing strand lie in the major groove of the heteroduplex. We present a model for the strand exchange intermediate in which homologous contacts preceding strand exchange arise in the minor groove of the substrate dsDNA.

DNA exhibits multi-stranded binding recognition on glass microarrays

In the course of exploring the hybridization properties of glass DNA microarrays, multi-stranded DNA structures were observed that could not be accounted for by classical Watson-Crick base pairing. Non-denatured double-stranded DNA array elements were shown to hybridize to single-stranded (ss)DNA probes. Similarly, ssDNA array elements were shown to bind duplex DNA probes. This led to a series of experiments demonstrating the formation of multi-stranded DNA structures on the surface of microarrays. These structures were observed with a number of heterogeneous sequences, including both purine and pyrimidine bases, with shared sequence identity between the ssDNA and one of the duplex strands. Furthermore, we observed a strong binding preference near the ends of duplexes containing a 3'-homologous strand. We suggest that such binding interactions on cationic solid surfaces could serve as a model for a number of biological processes mediated through multi-stranded DNA.

Evidence for a DNA triplex in a recombination-like motif: I. Recognition of Watson-Crick base pairs by natural bases in a high-stability triplex

Data are presented on a triplex type with two parallel homologous strands for which triplex formation is almost as strong as duplex formation at least for some sequences and even at pH 7 and 0.2 M NaCl. The evidence mainly rests upon comparing thermodynamic properties of similar systems. A paperclip oligonucleotide d(A12C4T12C4A12) with two linkers C4 obviously can form a triplex with parallel back-folded adenine strand regions, because the single melting transition of this complex splits in two transitions by introducing mismatches only in the third strand region. Respectively, a hairpin duplex d(A12C4T12) and a single strand d(A12) form a triplex as a 1:1 complex in which the second adenine strand is parallel oriented to the homologous one in the Watson-Crick paired duplex. In this system the melting temperature T(m) of the triplex is practically the same as that of the duplex d(A12)-d(T12), at least within a complex concentration range of 0.2-4.0 microM. The melting behaviour of complexes between triplex stabilizing ligand BePI and the system hairpin duplex plus single strand supports the triplex model. Non-denaturing gel electrophoresis suggests the existence of a triplex for a system in which five of the twelve A-T*A base triads are substituted by C-G*C base triads. The recognition between any substituted Watson-Crick base pair (X-Y) in the hairpin duplex d(A4XA7C4T7YT4) and the correspondingly replaced base (Z) in the third strand d(A4ZA7) is mutually selective. All triplexes with matching base substitutions (Z = X) have nearly the same stability (T(m) values from 29 to 33.5 degrees C), whereas triplexes with non-matching substitutions (Z not equal X) show a clearly reduced stability (T(m) values from 15 to 22 degrees C) at 2microM equimolar oligonucleotide concentration. Most nucleic acid triple helices hitherto known are limited to homopurine-homopyrimidine sequences in the target duplex. A stable triplex formation is demonstrated for inhomogeneous sequences tolerating at least 50% pyrimidine content in the homologous strands. On the basis of the surprisingly similar thermodynamic parameters for duplex and triplex, and of the fact that this triplex type seems to be more stable than many other natural DNA triplexes known, and on the basis of semiempirical and molecule mechanical calculations, we postulate bridging interactions of the third strand with the two other strands in the triplex according to the recombination motif. This triplex, denoted by us 'recombination-like form', tolerates heterogeneous base sequences.

Protein-free parallel triple-stranded DNA complex formation

A 14 nt DNA sequence 5'-AGAATGTGGCAAAG-3' from the zinc finger repeat of the human KRAB zinc finger protein gene ZNF91 bearing the intercalator 2-methoxy,6-chloro,9-amino acridine (Acr) attached to the sugar-phosphate backbone in various positions has been shown to form a specific triple helix (triplex) with a 16 bp hairpin (intramolecular) or a two-stranded (intermolecular) duplex having the identical sequence in the same (parallel) orientation. Intramolecular targets with the identical sequence in the antiparallel orientation and a non-specific target sequence were tested as controls. Apparent binding constants for formation of the triplex were determined by quantitating electrophoretic band shifts. Binding of the single-stranded oligonucleotide probe sequence to the target led to an increase in the fluorescence anisotropy of acridine. The parallel orientation of the two identical sequence segments was confirmed by measurement of fluorescence resonance energy transfer between the acridine on the 5'-end of the probe strand as donor and BODIPY-Texas Red on the 3'-amino group of either strand of the target duplex as acceptor. There was full protection from OsO(4)-bipyridine modification of thymines in the probe strand of the triplex, in accordance with the presumed triplex formation, which excluded displacement of the homologous duplex strand by the probe-intercalator conjugate. The implications of these results for the existence of protein-independent parallel triplexes are discussed.

Force-induced melting of the DNA double helix 1. Thermodynamic analysis

The highly cooperative elongation of a single B-DNA molecule to almost twice its contour length upon application of a stretching force is interpreted as force-induced DNA melting. This interpretation is based on the similarity between experimental and calculated stretching profiles, when the force-dependent free energy of melting is obtained directly from the experimental force versus extension curves of double- and single-stranded DNA. The high cooperativity of the overstretching transition is consistent with a melting interpretation. The ability of nicked DNA to withstand forces greater than that at the transition midpoint is explained as a result of the one-dimensional nature of the melting transition, which leads to alternating zones of melted and unmelted DNA even substantially above the melting midpoint. We discuss the relationship between force-induced melting and the B-to-S transition suggested by other authors. The recently measured effect on T7 DNA polymerase activity of the force applied to a ssDNA template is interpreted in terms of preferential stabilization of dsDNA by weak forces approximately equal to 7 pN.

Comb-type cationic copolymer expedites DNA strand exchange while stabilizing DNA duplex

The accelerating effect of cationic substances on the DNA strand exchange reaction between a 20 bp DNA duplex and its complementary single strand was studied. A polycationic comb-type copolymer, that consists of a poly(L-lysine) backbone and a dextran graft chain (PLL-g-Dex) and known to stabilize triplex DNA, expedites the strand exchange reaction under physiological relevant conditions. Electrostatically a small excess of the copolymer let to a 300-1500-fold increase in the DNA strand exchange while large excess of spermine or cetyltrimethylammonium bromide, a cationic detergent known to promote markedly hybridization of complementary DNA strands, shows only a slight effect. The efficacy of the copolymer was not affected by a 10 mM Mg2+ concentration. Notably the copolymer promotes the strand exchange reaction while it stabilizes double-stranded DNA. The stabilization of strand exchange intermediates consisting of the parent duplex and the single strand by the copolymer is believed to be responsible for the observed acceleration behavior.

Evidence for a four-strand exchange catalyzed by the RecA protein

Strand exchange between two duplexes is usually initiated as a three-strand event that requires the presence of a single-stranded overhang or gap in one of the two molecules. Here we show that the RecA protein can catalyze a four-strand exchange. Specifically, it can recombine short hairpin substrates with homologous stems provided that one of the hairpins possesses a chimeric DNA/RNA backbone. This four-strand exchange reaction goes to completion in the presence of ATPgammaS and releases a stable heteroduplex upon removal of the RecA protein. Under identical conditions, strand exchange between two DNA hairpins is incomplete and generates a nascent heteroduplex that rapidly dissociates when the RecA protein is denatured. Since presynaptic filament formation does not appear to melt either type of hairpin, we propose that exchange occurs between homologously aligned duplexes that are extended and unwound within a RecA filament. The first reaction provides a mechanism for gene targeting by chimeric double-hairpin oligonucleotides while the second reaction explains the ability of the RecA protein to transiently align double-stranded DNA molecules.

Binding of SSB and RecA protein to DNA-containing stem loop structures: SSB ensures the polarity of RecA polymerization on single-stranded DNA

Single-stranded DNA-binding proteins play an important role in homologous pairing and strand exchange promoted by the class of RecA-like proteins. It is presumed that SSB facilitates binding of RecA on to ssDNA by melting secondary structure, but direct physical evidence for the disruption of secondary structure by either SSB or RecA is still lacking. Using a series of oligonucleotides with increasing amounts of secondary structure, we show that stem loop structures impede contiguous binding of RecA and affect the rate of ATP hydrolysis. The electrophoretic mobility shift of a ternary complex of SSB-DNA-RecA and a binary complex of SSB-DNA are similar; however, the mechanism remains obscure. Binding of RecA on to stem loop is rapid in the presence of SSB or ATPgammaS and renders the complex resistant to cleavage by HaeIII, to higher amounts of competitor DNA or low temperature. The elongation of RecA filament in a 5' to 3' direction is halted at the proximal end of the stem. Consequently, RecA nucleates at the loop and cooperative binding propagates the RecA filament down the stem causing its disruption. The pattern of modification of thymine residues in the loop region indicates that RecA monomer is the minimum binding unit. Together, these results suggest that SSB plays a novel role in ensuring the directionality of RecA polymerization across stem loop in ssDNA. These observations have fundamental implications on the role of SSB in multiple aspects of cellular DNA metabolism.

Modeling DNA deformations

Recent developments have been made in modeling double-helical DNA at four levels of three-dimensional structure: the all-atom level, whereby an oligonucleotide duplex is surrounded by a shroud of solvent molecules; the base-pair level, with explicit backbone atoms; the mesoscopic level, that is, a few hundred base pairs, with the local duplex conformation described by knowledge-based harmonic energy functions; and the scale of several thousand nucleotides, with the duplex described as an ideal elastic rod. Predictions of the sequence-dependent bending and twisting of the double helix, as well as solvent- and force-induced B-->A and over-stretching conformational transitions, are compared with experimental data. These subtle conformational changes are critical to the functioning of the double helix, including its packaging in the close confines of the cell, the mutual fit of DNA and protein in nucleoprotein complexes, and the effective recognition of base pairs in recombination and transcription.