Equilibrium distributions of topological states in circular DNA: Interplay of supercoiling and knotting

Two variables define the topological state of closed double-stranded DNA: the knot type, K, and ΔLk, the linking number difference from relaxed DNA. The equilibrium distribution of probabilities of these states, P(ΔLk, K), is related to two conditional distributions: P(ΔLk|K), the distribution of ΔLk for a particular K, and P(K|ΔLk) and also to two simple distributions: P(ΔLk), the distribution of ΔLk irrespective of K, and P(K). We explored the relationships between these distributions. P(ΔLk, K), P(ΔLk), and P(K|ΔLk) were calculated from the simulated distributions of P(ΔLk|K) and of P(K). The calculated distributions agreed with previous experimental and theoretical results and greatly advanced on them. Our major focus was on P(K|ΔLk), the distribution of knot types for a particular value of ΔLk, which had not been evaluated previously. We found that unknotted circular DNA is not the most probable state beyond small values of ΔLk. Highly chiral knotted DNA has a lower free energy because it has less torsional deformation. Surprisingly, even at |ΔLk| > 12, only one or two knot types dominate the P(K|ΔLk) distribution despite the huge number of knots of comparable complexity. A large fraction of the knots found belong to the small family of torus knots. The relationship between supercoiling and knotting in vivo is discussed.

RNA polymerase sigma factor determines start-site selection but is not required for upstream promoter element activation on heteroduplex (bubble) templates

Sequence-selective transcription by bacterial RNA polymerase (RNAP) requires σ factor that participates in both promoter recognition and DNA melting. RNAP lacking σ (core enzyme) will initiate RNA synthesis from duplex ends, nicks, gaps, and single-stranded regions. We have used DNA templates containing short regions of heteroduplex (bubbles) to compare initiation in the presence and absence of various σ factors. Using bubble templates containing the σD-dependent flagellin promoter, with or without its associated upstream promoter (UP) element, we demonstrate that UP element stimulation occurs efficiently even in the absence of σ. This supports a model in which the UP element acts primarily through the α subunit of core enzyme to increase the initial association of RNAP with the promoter. Core and holoenzyme do differ substantially in the template positions chosen for initiation: σD restricts initiation to sites 8–9 nucleotides downstream of the conserved −10 element. Remarkably, σA also has a dramatic effect on start-site selection even though the σA holoenzyme is inactive on the corresponding homoduplexes. The start sites chosen by the σA holoenzyme are located 8 nucleotides downstream of sequences on the nontemplate strand that resemble the conserved −10 hexamer recognized by σA. Thus, σA appears to recognize the −10 region even in a single-stranded state. We propose that in addition to its described roles in promoter recognition and start-site melting, σ also localizes the transcription start site.

Torsional directed walks, entropic elasticity, and DNA twist stiffness

DNA and other biopolymers differ from classical polymers because of their torsional stiffness. This property changes the statistical character of their conformations under tension from a classical random walk to a problem we call the “torsional directed walk.” Motivated by a recent experiment on single lambda-DNA molecules [Strick, T. R., Allemand, J.-F., Bensimon, D., Bensimon, A. & Croquette, V. (1996) Science 271, 1835–1837], we formulate the torsional directed walk problem and solve it analytically in the appropriate force regime. Our technique affords a direct physical determination of the microscopic twist stiffness C and twist-stretch coupling D relevant for DNA functionality. The theory quantitatively fits existing experimental data for relative extension as a function of overtwist over a wide range of applied force; fitting to the experimental data yields the numerical values C = 120 nm and D = 50 nm. Future experiments will refine these values. We also predict that the phenomenon of reduction of effective twist stiffness by bend fluctuations should be testable in future single-molecule experiments, and we give its analytic form.

Uracil-DNA glycosylase–DNA substrate and product structures: Conformational strain promotes catalytic efficiency by coupled stereoelectronic effects

Enzymatic transformations of macromolecular substrates such as DNA repair enzyme/DNA transformations are commonly interpreted primarily by active-site functional-group chemistry that ignores their extensive interfaces. Yet human uracil–DNA glycosylase (UDG), an archetypical enzyme that initiates DNA base-excision repair, efficiently excises the damaged base uracil resulting from cytosine deamination even when active-site functional groups are deleted by mutagenesis. The 1.8-Å resolution substrate analogue and 2.0-Å resolution cleaved product cocrystal structures of UDG bound to double-stranded DNA suggest enzyme–DNA substrate-binding energy from the macromolecular interface is funneled into catalytic power at the active site. The architecturally stabilized closing of UDG enforces distortions of the uracil and deoxyribose in the flipped-out nucleotide substrate that are relieved by glycosylic bond cleavage in the product complex. This experimentally defined substrate stereochemistry implies the enzyme alters the orientation of three orthogonal electron orbitals to favor electron transpositions for glycosylic bond cleavage. By revealing the coupling of this anomeric effect to a delocalization of the glycosylic bond electrons into the uracil aromatic system, this structurally implicated mechanism resolves apparent paradoxes concerning the transpositions of electrons among orthogonal orbitals and the retention of catalytic efficiency despite mutational removal of active-site functional groups. These UDG/DNA structures and their implied dissociative excision chemistry suggest biology favors a chemistry for base-excision repair initiation that optimizes pathway coordination by product binding to avoid the release of cytotoxic and mutagenic intermediates. Similar excision chemistry may apply to other biological reaction pathways requiring the coordination of complex multistep chemical transformations.

A coupled complex of T4 DNA replication helicase (gp41) and polymerase (gp43) can perform rapid and processive DNA strand-displacement synthesis

We have developed a coupled helicase–polymerase DNA unwinding assay and have used it to monitor the rate of double-stranded DNA unwinding catalyzed by the phage T4 DNA replication helicase (gp41). This procedure can be used to follow helicase activity in subpopulations in systems in which the unwinding-synthesis reaction is not synchronized on all the substrate-template molecules. We show that T4 replication helicase (gp41) and polymerase (gp43) can be assembled onto a loading site located near the end of a long double-stranded DNA template in the presence of a macromolecular crowding agent, and that this coupled “two-protein” system can carry out ATP-dependent strand displacement DNA synthesis at physiological rates (400 to 500 bp per sec) and with high processivity in the absence of other T4 DNA replication proteins. These results suggest that a direct helicase–polymerase interaction may be central to fast and processive double-stranded DNA replication, and lead us to reconsider the roles of the other replication proteins in processivity control.

A synthetic cryIC gene, encoding a Bacillus thuringiensis δ-endotoxin, confers Spodoptera resistance in alfalfa and tobacco

Spodoptera species, representing widespread polyphagous insect pests, are resistant to Bacillus thuringiensis δ-endotoxins used thus far as insecticides in transgenic plants. Here we describe the chemical synthesis of a cryIC gene by a novel template directed ligation–PCR method. This simple and economical method to construct large synthetic genes can be used when routine resynthesis of genes is required. Chemically phosphorylated adjacent oligonucleotides of the gene to be synthesized are assembled and ligated on a single-stranded, partially homologous template derived from a wild-type gene (cryIC in our case) by a thermostable Pfu DNA ligase using repeated cycles of melting, annealing, and ligation. The resulting synthetic DNA strands are selectively amplified by PCR with short specific flanking primers that are complementary only to the new synthetic DNA. Optimized expression of the synthetic cryIC gene in alfalfa and tobacco results in the production of 0.01–0.2% of total soluble proteins as CryIC toxin and provides protection against the Egyptian cotton leafworm (Spodoptera littoralis) and the beet armyworm (Spodoptera exigua). To facilitate selection and breeding of Spodoptera-resistant plants, the cryIC gene was linked to a pat gene, conferring resistance to the herbicide BASTA.

t(11;22)(q23;q11.2) in acute myeloid leukemia of infant twins fuses MLL with hCDCrel, a cell division cycle gene in the genomic region of deletion in DiGeorge and velocardiofacial syndromes

We examined the MLL genomic translocation breakpoint in acute myeloid leukemia of infant twins. Southern blot analysis in both cases showed two identical MLL gene rearrangements indicating chromosomal translocation. The rearrangements were detectable in the second twin before signs of clinical disease and the intensity relative to the normal fragment indicated that the translocation was not constitutional. Fluorescence in situ hybridization with an MLL-specific probe and karyotype analyses suggested t(11;22)(q23;q11.2) disrupting MLL. Known 5′ sequence from MLL but unknown 3′ sequence from chromosome band 22q11.2 formed the breakpoint junction on the der(11) chromosome. We used panhandle variant PCR to clone the translocation breakpoint. By ligating a single-stranded oligonucleotide that was homologous to known 5′ MLL genomic sequence to the 5′ ends of BamHI-digested DNA through a bridging oligonucleotide, we formed the stem–loop template for panhandle variant PCR which yielded products of 3.9 kb. The MLL genomic breakpoint was in intron 7. The sequence of the partner DNA from band 22q11.2 was identical to the hCDCrel (human cell division cycle related) gene that maps to the region commonly deleted in DiGeorge and velocardiofacial syndromes. Both MLL and hCDCrel contained homologous CT, TTTGTG, and GAA sequences within a few base pairs of their respective breakpoints, which may have been important in uniting these two genes by translocation. Reverse transcriptase-PCR amplified an in-frame fusion of MLL exon 7 to hCDCrel exon 3, indicating that an MLL-hCDCrel chimeric mRNA had been transcribed. Panhandle variant PCR is a powerful strategy for cloning translocation breakpoints where the partner gene is undetermined. This application of the method identified a region of chromosome band 22q11.2 involved in both leukemia and a constitutional disorder.

Efficient translesion replication in the absence of Escherichia coli Umu proteins and 3′–5′ exonuclease proofreading function

Translesion replication (TR) past a cyclobutane pyrimidine dimer in Escherichia coli normally requires the UmuD′2C complex, RecA protein, and DNA polymerase III holoenzyme (pol III). However, we find that efficient TR can occur in the absence of the Umu proteins if the 3′–5′ exonuclease proofreading activity of the pol III ɛ-subunit also is disabled. TR was measured in isogenic uvrA6 ΔumuDC strains carrying the dominant negative dnaQ allele, mutD5, or ΔdnaQ spq-2 mutations by transfecting them with single-stranded M13-based vectors containing a specifically located cis-syn T–T dimer. As expected, little TR was observed in the ΔumuDC dnaQ+ strain. Surprisingly, 26% TR occurred in UV-irradiated ΔumuDC mutD5 cells, one-half the frequency found in a uvrA6 umuDC+mutD5 strain. lexA3 (Ind−) derivatives of the strains showed that this TR was contingent on two inducible functions, one LexA-dependent, responsible for ≈70% of the TR, and another LexA-independent, responsible for the remaining ≈30%. Curiously, the ΔumuDC ΔdnaQ spq-2 strain exhibited only the LexA-independent level of TR. The cause of this result appears to be the spq-2 allele, a dnaE mutation required for viability in ΔdnaQ strains, since introduction of spq-2 into the ΔumuDC mutD5 strain also reduces the frequency of TR to the LexA-independent level. The molecular mechanism responsible for the LexA-independent TR is unknown but may be related to the UVM phenomenon [Palejwala, V. A., Wang, G. E., Murphy, H. S. & Humayun, M. Z. (1995) J. Bacteriol. 177, 6041–6048]. LexA-dependent TR does not result from the induction of pol II, since TR in the ΔumuDC mutD5 strain is unchanged by introduction of a ΔpolB mutation.

Immobilization of oligodeoxyribonucleotides with multiple anchors to microchips

Facile modification of oligodeoxyribonucleotides is required for efficient immobilization to a pre-activated glass surface. This report presents an oligodeoxyribonucleotide which contains a hairpin stem–loop structure with multiple phosphorothioate moieties in the loop. These moieties are used to anchor the oligo to glass slides that are pre-activated with bromoacetamidopropylsilane. The efficiency of the attachment reaction was improved by increasing the number of phosphorothioates in the loop, as shown in the remarkable enhancement of template hybridization and single base extension through catalysis by DNA polymerase. The loop and stem presumably serve as lateral spacers between neighboring oligodeoxyribonucleotides and as a linker arm between the glass surface and the single-stranded sequence of interest. The oligodeoxyribonucleotides of this hairpin stem–loop architecture with multiple phosphorothioate moieties have broad application in DNA chip-based gene analysis.

Thermostable and site-specific DNA binding of the gene product ORF56 from the Sulfolobus islandicus plasmid pRN1, a putative archael plasmid copy control protein

There is still a lack of information on the specific characteristics of DNA-binding proteins from hyperthermophiles. Here we report on the product of the gene orf56 from plasmid pRN1 of the acidophilic and thermophilic archaeon Sulfolobus islandicus. orf56 has not been characterised yet but low sequence similarily to several eubacterial plasmid-encoded genes suggests that this 6.5 kDa protein is a sequence-specific DNA-binding protein. The DNA-binding properties of ORF56, expressed in Escherichia coli, have been investigated by EMSA experiments and by fluorescence anisotropy measurements. Recombinant ORF56 binds to double-stranded DNA, specifically to an inverted repeat located within the promoter of orf56. Binding to this site could down-regulate transcription of the orf56 gene and also of the overlapping orf904 gene, encoding the putative initiator protein of plasmid replication. By gel filtration and chemical crosslinking we have shown that ORF56 is a dimeric protein. Stoichiometric fluorescence anisotropy titrations further indicate that ORF56 binds as a tetramer to the inverted repeat of its target binding site. CD spectroscopy points to a significant increase in ordered secondary structure of ORF56 upon binding DNA. ORF56 binds without apparent cooperativity to its target DNA with a dissociation constant in the nanomolar range. Quantitative analysis of binding isotherms performed at various salt concentrations and at different temperatures indicates that approximately seven ions are released upon complex formation and that complex formation is accompanied by a change in heat capacity of –6.2 kJ/mol. Furthermore, recombinant ORF56 proved to be highly thermostable and is able to bind DNA up to 85°C.

Protein–RNA interactions: a structural analysis

A detailed computational analysis of 32 protein–RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein–double-stranded DNA and protein–single-stranded DNA complexes. The interface properties of the protein–RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein–RNA and protein–DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein–RNA complexes, backbone contacts were more dominant in the protein–DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.

Efficient construction of long DNA duplexes with internal non-Watson–Crick motifs and modifications

We have developed a semi-synthetic approach for preparing long stretches of DNA (>100 bp) containing internal chemical modifications and/or non-Watson–Crick structural motifs which relies on splint-free, cell-free DNA ligations and recycling of side-products by non-PCR thermal cycling. A double-stranded DNA PCR fragment containing a polylinker in its middle is digested with two restriction enzymes and a small insert (∼20 bp) containing the modification or non-Watson–Crick motif of interest is introduced into the middle. Incorrect products are recycled to starting materials by digestion with appropriate restriction enzymes, while the correct product is resistant to digestion since it does not contain these restriction sites. This semi-synthetic approach offers several advantages over DNA splint-mediated ligations, including fewer steps, substantially higher yields (∼60% overall yield) and ease of use. This method has numerous potential applications, including the introduction of modifications such as fluorophores and cross-linking agents into DNA, controlling the shape of DNA on a large scale and the study of non-sequence-specific nucleic acid–protein interactions.

A nicked duplex decamer DNA with a PEG6 tether

A dumbbell double-stranded DNA decamer tethered with a hexaethylene glycol linker moiety (DDSDPEG), with a nick in the centre of one strand, has been synthesised. The standard NMR methods, E.COSY, TOCSY, NOESY and HMQC, were used to measure 1H, 31P and T1 spectral parameters. Molecular modelling using rMD-simulated annealing was used to compute the structure. Scalar couplings and dipolar contacts show that the molecule adopts a right-handed B-DNA helix in 38 mM phosphate buffer at pH 7. Its high melting temperature confirms the good base stacking and stability of the duplex. This is partly attributed to the presence of the PEG6 linker at both ends of the duplex that restricts the dynamics of the stem pentamers and thus stabilises the oligonucleotide. The inspection of the global parameters shows that the linker does not distort the B-DNA geometry. The computed structure suggests that the presence of the nick is not disturbing the overall tertiary structure, base pair geometry or duplex base pairing to a substantial extent. The nick has, however, a noticeable impact on the local geometry at the nick site, indicated clearly by NMR analysis and reflected in the conformational parameters of the computed structure. The 1H spectra also show much sharper resonances in the presence of K+ indicating that conformational heterogeneity of DDSDPEG is reduced in the presence of potassium as compared to sodium or caesium ions. At the same time the 1H resonances have longer T1 times. This parameter is suggested as a sensitive gauge of stabilisation.

Herpes simplex virus DNA packaging sequences adopt novel structures that are specifically recognized by a component of the cleavage and packaging machinery

The product of the herpes simplex virus type 1 UL28 gene is essential for cleavage of concatemeric viral DNA into genome-length units and packaging of this DNA into viral procapsids. To address the role of UL28 in this process, purified UL28 protein was assayed for the ability to recognize conserved herpesvirus DNA packaging sequences. We report that DNA fragments containing the pac1 DNA packaging motif can be induced by heat treatment to adopt novel DNA conformations that migrate faster than the corresponding duplex in nondenaturing gels. Surprisingly, these novel DNA structures are high-affinity substrates for UL28 protein binding, whereas double-stranded DNA of identical sequence composition is not recognized by UL28 protein. We demonstrate that only one strand of the pac1 motif is responsible for the formation of novel DNA structures that are bound tightly and specifically by UL28 protein. To determine the relevance of the observed UL28 protein–pac1 interaction to the cleavage and packaging process, we have analyzed the binding affinity of UL28 protein for pac1 mutants previously shown to be deficient in cleavage and packaging in vivo. Each of the pac1 mutants exhibited a decrease in DNA binding by UL28 protein that correlated directly with the reported reduction in cleavage and packaging efficiency, thereby supporting a role for the UL28 protein–pac1 interaction in vivo. These data therefore suggest that the formation of novel DNA structures by the pac1 motif confers added specificity on recognition of DNA packaging sequences by the UL28-encoded component of the herpesvirus cleavage and packaging machinery.

The length of telomeric G-rich strand 3′-overhang measured by oligonucleotide ligation assay

A typical G-rich telomeric DNA strand, which runs 5′→3′ toward the chromosome ends, protrudes by several nucleotides in lower eukaryotes. In human chromosomes long G-rich 3′-overhangs have been found. Apart from the standard G-rich tail, several non-canonical terminal structures have been proposed. However, the mechanism of long-tail formation, the presence and the role of these structures in telomere maintenance or shortening are not completely understood. In a search for a simple method to accurately measure the 3′-overhang we have established a protocol based on the ligation of telomeric oligonucleotide hybridized to non-denatured DNA under stringent conditions (oligonucleotide ligation assay with telomeric repeat oligonucleotide). This method enabled us to detect a large proportion of G-rich single-stranded telomeric DNA that was as short as 24 nt. Nevertheless, we showed G-tails longer than 400 nt. In all tested cells the lengths ranging from 108 to 270 nt represented only 37% of the whole molecule population, while 56–62% were <90 nt. Our protocol provides a simple and sensitive method for measuring the length of naturally occurring unpaired repeated DNA.