Comparison of analyses of DNA curvature.

Popular programs for characterizing DNA structure include Curves 5.1 (Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 63-91, 1988; Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 655-67, 1989) and Freehelix98 (Dickerson, R. E., Nucleic Acids Res. 26, 1906-1926, 1998), along with the more recent 3DNA (X. J. Lu, Z. Shakked and W. K. Olson., J. Mol. Biol. 300, 819-840 (2000). Given input of structural coordinates, all of these programs return values of the local helical parameters, such as roll, tilt, twist, etc. The first two programs also provide characterization of global curvature. Madbend (Strahs, D. and Schlick, T., J. Mol. Biol. 301, 643-663, 2000), a program that computes global curvature from local roll, tilt, and twist parameters, can be applied to the output of all three structural programs. We have compared the curvature predicted by the three programs with and without the use of Madbend. Global bend magnitudes and directions as well as values of helical kinks were calculated for four high-resolution DNA structures and four model DNA helices. Global curvature determined by Curves 5.1 without Madbend was found to differ from values obtained using Freehelix98 with or without Madbend or 3DNA and Curves 5.1 with Madbend. Using model helices, this difference was attributed the fact that Curves 5.1 is the only program sensitive to changes in axial displacement, such as shift and slide. Madbend produced robust values of bend magnitude and direction, and displayed little sensitivity to axis displacement or the source of local helical parameters. Madbend also appears to be the method of choice for bending comparisons of high-resolution structures with results from cyclization kinetics, a method that measures DNA curvature as a vectorial sum of local roll and tilt angles.

Kinetic consequences of covalent linkage of DNA binding polyamides.

Polyamides composed of N-methylpyrrole (Py) and N-methylimidazole (Im) subunits can bind in the minor groove of DNA at predetermined sequences with subnanomolar affinity and high specificity. Covalent linkage of polymer subunits using a gamma-aminobutyric acid linker has been shown to increase both the affinity and specificity of polyamides. Using a fluorescence detected stopped-flow assay, we have studied the differences in association and dissociation kinetics of a series of polyamides representing unlinked, hairpin and cyclic analogues of the four ring polyamide ImPyPyPy-beta-Dp. Whereas the large differences seen in the equilibrium association constants between the unlinked and covalently linked polyamides are primarily due to higher association rate constants, discrimination between matched and mismatched sites by each polyamide can be ascribed in large part to differences in their dissociation rate constants. The consequences of this kinetic behavior for future design are discussed.

RNA folding pathways.

A general overview of the questions and problems in RNA folding is presented. Topics include the differences in the folding problems/questions that apply to RNA versus proteins, methods for determination of final structures, folding versus unfolding, resolution of space and time, and conformational switching.

Global structure and mechanical properties of a 10-bp nucleosome positioning motif.

The method of DNA cyclization kinetics reveals special properties of the TATAAACGCC sequence motif found in DNA sequences that have high affinity for core histones. Replacement of 30 bp of generic DNA by three 10-bp repeats of the motif in small cyclization constructs increases cyclization rates by two orders of magnitude. We document a 13 degrees bend in the motif and characterize the direction of curvature. The bending force constant is smaller by nearly 2-fold and there is a 35% decrease in the twist modulus, relative to generic DNA. These features are the likely source of the high affinity for bending around core histones to form nucleosomes. Our results establish a protocol for determination of the ensemble-averaged global solution structure and mechanical properties of any approximately 10-bp DNA sequence element of interest, providing information complementary to that from NMR and crystallographic structural studies.

On the kinetics of distamycin binding to its target sites on duplex DNA.

Distamycin A is a well known polyamide antibiotic that can bind in the minor groove of duplex DNA primarily at AT-rich sequences both as a monomer or as a side-by-side antiparallel dimer. The association phase of the distamycin binding reaction has not been studied in either of its binding modes, because of the lack of an adequate UV or CD signal at the low concentrations needed to monitor the fast bimolecular reaction. We report a significant increase in fluorescence amplitude, accompanied by a small red shift, on binding distamycin to its specific target sites. This signal can be used to monitor drug binding in steady-state and time-resolved processes. Distamycin shows extremely fast association with the 1:1 binding site, with a bimolecular rate of 7 x 10(7) M(-1) small middle dots(-1) and also fairly rapid dissociation ( approximately 3 s(-1)). When DNA is in excess, there is a slow component in the association reaction whose rate decreases strongly with increasing DNA concentration. Binding of the drug to the 2:1 site occurs in two distinct steps: fast, sequential binding of each drug molecule to the DNA with a bimolecular rate comparable to that at the 1:1 site, followed by a slow ( approximately 4 s(-1)) equilibration to the final population. Dissociation from the 2:1 site is approximately 40-fold slower than from the 1:1 site. This study provides the groundwork for analysis of the binding kinetics of longer polyamides and covalently linked polyamides that have recently been shown to inhibit transcription in vivo.

Interaction of calicheamicin gamma1(I) and its related carbohydrates with DNA-protein complexes.

We report studies of the contribution of DNA structure, holding the sequence constant, to the affinity of calicheamicin gamma(1)(I) and its aryltetrasaccharide moiety for DNA. We used polynucleotide chains as models of known protein-binding sequences [the catabolite activator protein (CAP) consensus sequence, AP-1 and cAMP response element (CRE) sites] in their free and protein-bound forms. The proteins were selected to provide examples in which the minor-groove binding site for the carbohydrate is (CAP) or is not (GCN4) covered by the protein. Additionally, peptides related to the GCN4 and CREB families, which have different bending effects on their DNA-binding sites, were used. We observe that proteins of the CREB class, which induce a tendency to bend toward the minor groove at the center of the site, inhibit drug-cleavage sites located at the center of the free AP-1 or CRE DNA sites. In the case of GCN4, which does not induce DNA bending, there is no effect on calicheamicin cleavage of the CRE site, but we observe a GCN4-induced rearrangement of the cutting pattern in the AP-1 site. This effect may arise from either a subtle local conformational rearrangement not accompanied by bending or a localized reduction in DNA flexibility. Whereas GCN4 binding is not inhibited by the calicheamicin aryltetrasaccharide, binding of CAP to its DNA target is significantly inhibited, and calicheamicin cutting of DNA at the center of the CAP-DNA complex site is strongly reduced by protein binding. This result probably reflects steric inhibition of drug binding by the protein.

Characterization of the solution conformations of unbound and Tat peptide-bound forms of HIV-1 TAR RNA.

Basic peptides from the carboxy terminus of the HIV-1 Tat protein bind to the apical stem-loop region of TAR RNA with high affinity and moderate specificity. The conformations of the unbound and 24 residue Tat peptide (Tfr24)-bound forms of TAR RNA have been characterized by NMR spectroscopy. The unbound form of TAR exists in major and minor forms having different trinucleotide bulge conformations. A specific TAR RNA conformational change is observed upon complex formation with Tfr24, consisting of coaxial stacking of helical stems and base triple formation. A U23-A27-U38 base triple is proposed based on exchangeable proton NMR data, where U23 forms a base pair with A27 in the major groove. No evidence for base triple formation was found for Tat peptides in which lysine residues are extensively substituted for arginine.

Decreased imino proton exchange and base-pair opening in the IHF-DNA complex measured by NMR.

Integration Host Factor, IHF, is an E. coli DNA binding protein that imposes a substantial bend on DNA. Previous footprinting studies and bending assays have characterized several recognition sequences in the bacterial and lambda phage genome as unique in the way they are bound by IHF. We have chosen one of the lambda phage sites, H1, for study because it presents a small yet sequence-specific substrate for NMR analysis of the complex. A 19 base-pair duplex, H19, corresponding to the recognition sequence at the H1 site was constructed by isotopically labeling one of the strands with 15N. (1H, 15N) heteronuclear NMR experiments aided in assigning the imino proton resonances of the DNA alone and in complex with IHF. The NMR results are consistent with a mode of binding observed in the recent crystal structure of IHF bound to another of its sites from the lambda phage genome. Additionally, the dramatic change that IHF imposes on the imino proton chemical shifts is indicative of a severe deviation from canonical B-DNA structure. In order to understand the dynamic properties of the DNA in the complex with IHF, the exchange rates of the imino protons with the solvent have been measured for H19 with and without IHF bound. A drastic reduction in exchange is observed for the imino protons in the IHF bound DNA. In the DNA-protein complex, groups of adjacent base-pair exchange at the same rate, and appear to close more slowly than the rate of imino proton exchange with bulk water, since their exchange rate is independent of catalyst concentration. We infer that segments of the double helix as large as 6 bp open in a cooperative process, and remain open much longer than is typical for opening fluctuations in naked duplex DNA. We discuss these results in terms of the specific protein-DNA contacts observed in the crystal structure.

An in vitro transcription assay for probing drug-DNA interactions at individual drug sites.

An in vitro transcription assay of drug-DNA interactions has been described and is based largely on the stable lac UV5-initiated transcription complex. This system utilizes a synchronized population of radiolabeled nascent RNA 10 nucleotides long. Reaction of this initiated transcription complex with drug and subsequent elongation of the nascent RNA by Escherichia coli RNA polymerase, reveals blockages at drug binding sites. From these blockages it is possible to obtain four features of the drug-DNA interaction: the sequence of preferred drug binding sites, the relative drug occupancy at each binding site, the drug dissociation rate at each site, and the probability of drug-induced termination of transcription at each site. The unidirectional transcription assay has been extended to a two-promoter, counter-directed system, which yields a bidirectional transcription footprint of drug sites.

Characterization of covalent adriamycin-DNA adducts.

Adriamycin is a popular antineoplastic agent whose ability to form covalent adducts with DNA has been correlated to cellular apoptosis (programmed cell death) in tumor models. We have isolated and purified this adduct formed under oxido-reductive (Fenton) conditions in Tris buffer. We show by homo- and heteronuclear NMR spectroscopy that the covalent Adriamycin-DNA adduct is structurally equivalent to that resulting from direct reaction with formaldehyde. Covalent linkage of the drug to one of the DNA strands confers remarkable stability to the duplex, indicated by a 162-fold reduction in the rate of strand displacement compared with the complex with noncovalently bound drug. Glyceraldehyde also engenders covalent Adriamycin-DNA complexes, providing a possible relevant biological context for in vivo adduct formation.

The solution structure of an RNA loop-loop complex: the ColE1 inverted loop sequence.

BACKGROUND: Replication of the ColE1 plasmid of Escherichia coli is regulated by the interaction of sense and antisense plasmid-encoded transcripts. The antisense RNA I negatively regulates the replication of the plasmid by duplex formation with complementary RNA II. The interaction is initiated by the formation of a double helix between seven-nucleotide loops from each RNA and is stabilized by binding of the RNA one modulator (ROM) protein. The ROM protein is thought to recognize a specific RNA structure, regardless of sequence. RESULTS: The solution structure of a loop-loop complex between model RNA hairpins that resemble RNA I and RNA II has been determined by nuclear magnetic resonance spectroscopy. The model hairpins have loop sequences inverted 5' to 3' relative to the wild-type sequence and were chosen because of their complex's slow dissociation in comparison to the wild type. The complex has continuous stacking from the 3'-side of one stem helix through the loop-loop helix to the other stem helix. One residue from each hairpin has a unique phosphodiester bond which bridges and narrows the major groove. These bridging phosphates are in close proximity to the phosphate groups of the adjacent bases, forming unique structural motifs called phosphate clusters. The purine residue at the 3'-end of the loop-loop helix of one RNA stacks on a purine residue on the 5'-side of the other RNA stem, and there are strong cross-strand stacking interactions between guanine bases in the stem helices adjacent to the loops. CONCLUSIONS: Unique distortions, such as the strong bend and the phosphate clusters flanking the major groove of the loop-loop helix, provide an attractive nonsequence-specific structural feature for recognition by the ROM protein. The structure provides a basis for rationalizing the sequence dependence of the stability of loop-loop interaction.

Simultaneous determination of helical unwinding angles and intrinsic association constants in ligand-DNA complexes: the interaction between DNA and calichearubicin B.

We present a helical unwinding assay for reversibly binding DNA ligands that uses closed circular DNA, topoisomerase I (Topo I), and two-dimensional agarose gel electrophoresis. Serially diluted Topo I relaxation reactions at constant DNA/ligand ratio are performed, and the resulting apparent unwinding of the closed circular DNA is used to calculate both ligand unwinding angle (phi) and intrinsic association constant (Ka). Mathematical treatment of apparent unwinding is formally analogous to that of apparent extinction coefficient data for optical binding titrations. Extrapolation to infinite DNA concentration yields the true unwinding angle of a given ligand and its association constant under Topo I relaxation conditions. Thus this assay delivers simultaneous structural and thermodynamic information describing the ligand-DNA complex. The utility of this assay has been demonstrated by using calichearubicin B (CRB), a synthetic hybrid molecule containing the anthraquinone chromophore of (DA) and the carbohydrate domain of calicheamicin gamma1I. The unwinding angle for CRB calculated by this method is -5. 3 +/- 0.5 degrees. Its Ka value is 0.20 x 10(6) M-1. For comparison, the unwinding angles of ethidium bromide and DA have been independently calculated, and the results are in agreement with canonical values for these compounds. Although a stronger binder to selected sites, CRB is a less potent unwinder than its parent compound DA. The assay requires only small amounts of ligand and offers an attractive option for analysis of DNA binding by synthetic and natural compounds.

DNA-binding domains of Fos and Jun do not induce DNA curvature: an investigation with solution and gel methods.

We demonstrate the use of a DNA minicircle competition binding assay, together with DNA cyclization kinetics and gel-phasing methods, to show that the DNA-binding domains (dbd) of the heterodimeric leucine zipper protein Fos-Jun do not bend the AP-1 target site. Our DNA constructs contain an AP-1 site phased by 1-4 helical turns against an A-tract-directed bend. Competition binding experiments reveal that (dbd)Fos-Jun has a slight preference for binding to linear over circular AP-1 DNAs, independent of whether the site faces in or out on the circle. This result suggests that (dbd)Fos-Jun slightly stiffens rather than bends its DNA target site. A single A-tract bend replacing the AP-1 site is readily detected by its effect on cyclization kinetics, in contrast to the observations for Fos-Jun bound at the AP-1 locus. In contrast, comparative electrophoresis reveals that Fos-Jun-DNA complexes, in which the A-tract bend is positioned close (1-2 helical turns) to the AP-1 site, show phase-dependent variations in gel mobilities that are comparable with those observed when a single A-tract bend replaces the AP-1 site. Whereas gel mobility variations of Fos-Jun-DNA complexes decrease linearly with increasing Mg2+ contained in the gel, the solution binding preference of (dbd)Fos-Jun for linear over circular DNAs is independent of Mg2+ concentration. Hence, gel mobility variations of Fos-Jun-DNA complexes are not indicative of (dbd)Fos-Jun-induced DNA bending (upper limit 5 degrees) in the low salt conditions of gel electrophoresis. Instead, we propose that the gel anomalies depend on the steric relationship of the leucine zipper region with respect to a DNA bend.

Characterization of the ATF/CREB site and its complex with GCN4.

We have studied DNA minicircles containing the ATF/CREB binding site for GCN4 by using a combination of cyclization kinetics experiments and Monte Carlo simulations. Cyclization rates were determined with and without GCN4 for DNA constructs containing the ATF/CREB site separated from a phased A-tract multimer bend by a variable length phasing adaptor. The cyclization results show that GCN4 binding does not significantly change the conformation of the ATF/CREB site, which is intrinsically slightly bent toward the major groove. Monte Carlo simulations quantitate the ATF/CREB site structure as an 8 degrees bend toward the major groove in a coordinate frame near the center of the site. The ATF/CREB site is underwound by 53 degrees relative to the related AP-1 site DNA. The effect of GCN4 binding can be modeled either as a decrease in the local flexibility, corresponding to an estimated 60% increase in the persistence length for the 10-bp binding site, or possibly as a small decrease (1 degrees) in intrinsic bend angle. Our results agree with recent electrophoretic and crystallographic studies and demonstrate that cyclization and simulation can characterize subtle changes in DNA structure and flexibility.

Measurement of the DNA bend angle induced by the catabolite activator protein using Monte Carlo simulation of cyclization kinetics.

A Monte Carlo simulation method for studying DNA cyclization (or ring-closure) has been extended to the case of protein-induced bending, and its application to experimental data has been demonstrated. Estimates for the geometric parameters describing the DNA bend induced by the catabolite activator protein (CAP or CRP) were obtained which correctly predict experimental DNA cyclization probabilities (J factors), determined for a set of 11 150 to 166 bp DNA restriction fragments bearing A tracts phased against CAP binding sites. We find that simulation of out-of-phase molecules is difficult and time consuming, requiring the geometric parameters to be optimized individually rather than globally. A wedge angle model for DNA bending was found to make reasonable predictions for the free DNA. The bend angle in the CAP-DNA complex is estimated to be 85 to 90 degrees, in agreement with estimates from gel electrophoresis and X-ray co-crystal structures. Since the DNA is found to have a pre-existing bend of 15 degrees, the change in bend angle induced by CAP is 70 to 75 degrees, in a agreement with an estimate from topological measurements. We find evidence for slight (approximately 10 degrees) unwinding by CAP. The persistence length and helical repeat of the unbound portion of the DNA are in accord with literature-cited values, but the best-fit DNA torsional modulus C is found to be 1.7 (+/- 0.2) x 10(-19) erg. cm, versus literature estimates and best-fit values for the free DNA of 2.0 x 10(-19) to 3.4 x 10(-19) erg.com. Simulations using this low value of C predict that cyclization of molecules with out-of-phase bends proceeds via undertwisting or overtwisting of the DNA between the bends, so as to align the bends, rather than through conformations with substantial writhe. We present experiments on the topoisomers formed by cyclization with CAP which support this conclusion, and thereby rationalize the surprising result that cyclization can actually be enhanced by out-of-phase bends if the twist required to align the bends improves the torsional alignment of the ends. The relationship between the present work and previous studies on DNA bending by CAP is discussed, and recommendations are given for the efficient application of the cyclization/simulation approach to DNA bending.

Measurement of diffusion constants for nucleic acids by NMR.

Pulsed field-gradient NMR experiments can be used to measure the diffusion constants of nucleic acids. The diffusion constants measured in this way for double-helical DNAs of defined length agree well both with theory and with measurements done using other techniques. When applied to RNAs, this experiment easily distinguishes duplex RNAs from RNA hairpins and thus it can solve one of the perennial problems faced by RNA spectroscopists, i.e. assessing whether their samples are monomeric or not.