DNA structure and gene regulation

Effect of loop sequence on unzipping of short DNA hairpins

The dependence of stability on the sequence of a DNA hairpin has been investigated through atomistic simulations. For this, a sequence of 16 bases of a hairpin, which consists of a loop of four bases and a stem of six base pairs, has been considered. We have taken eight different sequences, where the first five base pairs were kept fixed in all sequences, whereas the loop sequence and the identity of the duplex base pair closing the loop have been varied. For these hairpin structures, force-induced melting (unzipping) studies were carried out to investigate the effect of the variables on the stability of hairpin. The temperature at which half of the base pairs are open is termed the melting temperature. We defined the unzipping force F_{h} (half of the base pairs are open) and showed that it may not provide the effect of closing the base pair or loop sequence on the stability of the DNA hairpin. In order to have a better understanding of the stability of a DNA hairpin, the closing base pair or hairpin loop must be open. This requires complete opening of the stem. We defined a force F_{c} at which all base pairs of the stem are open, and we showed that the F_{c} gives better understanding of DNA hairpin stability.

Uncertainty in protein-ligand binding constants: asymmetric confidence intervals versus standard errors

Equilibrium binding constants (K_b) between chemical compounds and target proteins or between interacting proteins provide a quantitative understanding of biological interaction mechanisms. Reported uncertainties of measured experimental parameters are critical for decision-making in many scientific areas, e.g., in lead compound discovery processes and in comparing computational predictions with experimental results. Uncertainties in measured K_b values are commonly represented by a symmetric normal distribution, often quoted in terms of the experimental value plus-minus the standard deviation. However, in general, the distributions of measured K_b (and equivalent K_d) values and the corresponding free energy change ΔG_b are all asymmetric to varying degree. Here, using a simulation approach, we illustrate the effect of asymmetric K_b distributions within the realm of isothermal titration calorimetry (ITC) experiments. Further we illustrate the known, but perhaps not widely appreciated, fact that when distributions of any of K_b, K_d and ΔG_b are transformed into each other, their degree of asymmetry is changed. Consequently, we recommend that a more accurate way of expressing the uncertainties of K_b, K_d, and ΔG_b values is to consistently report 95% confidence intervals, in line with other authors' suggestions. The ways to obtain such error ranges are discussed in detail and exemplified for a binding reaction obtained by ITC.

Anomalous Separation of Small Y-Chromosomal DNA Fragments on Microchip Electrophoresis

We investigated an anomalous DNA separation where two DNA fragments from the human Y-chromosome sY638 (64 bp) and sY592 (65 bp), with only one base pair difference, were separated. This result is abnormal since in a previous study, we found that 5 bp was the minimum difference between two DNA fragments that the microchip electrophoresis system can separate. The formation of a mini-loop in the structure of the DNA fragment of sY638 (64 bp) was strongly expected to be the reason. To investigate this, we synthesized three modified DNA fragments for sY638 (64 bp), and the modifications were in two expected locations for possible mini-loop formation. Later, the separation between sY592 (65 bp) and the three modified fragments of sY638 (64 bp) was not possible. Thus, we conclude that the formation of a mini-loop in the structure of the DNA is the reason behind this anomalous separation.

Structure-specific nucleic acid recognition by L-motifs and their diverse roles in expression and regulation of the genome

The high-mobility group (HMG) domain containing proteins regulate transcription, DNA replication and recombination. They adopt L-shaped folds and are structure-specific DNA binding motifs. Here, I define the L-motif super-family that consists of DNA-binding HMG-box proteins and the L-motif of the histone mRNA binding domain of stem-loop binding protein (SLBP). The SLBP L-motif and HMG-box domains adopt similar L-shaped folds with three α-helices and two or three small hydrophobic cores that stabilize the overall fold, but have very different and distinct modes of nucleic acid recognition. A comparison of the structure, dynamics, protein-protein and nucleic acid interactions, and regulation by PTMs of the SLBP and the HMG-box L-motifs reveals the versatile and diverse modes by which L-motifs utilize their surfaces for structure-specific recognition of nucleic acids to regulate gene expression.

Microarray-based STR genotyping using RecA-mediated ligation

We describe a novel assay capable of accurately determining the length of short tandem repeat (STR) alleles. STR genotyping is achieved utilizing RecA-mediated ligation (RML), which combines the high fidelity of RecA-mediated homology searching with allele-specific ligation. RecA catalyzes the pairing of synthetic oligonucleotides with one strand of a double-stranded DNA target, in this case a PCR amplicon. Ligation occurs only when two adjacent oligonucleotides are base paired to the STR region without any overlap or gap. RecA activity is required to overcome the inherent difficulty of annealing repeated sequences in register. This assay is capable of determining STR genotypes of human samples, is easily adapted to high throughput or automated systems and can have widespread utility in diagnostic and forensic applications.

Bent helical structure in kinetoplast DNA

We have investigated the unusual physical properties of a restriction fragment of Leishmania tarentolae kinetoplast DNA. A gel-purified fragment comprising slightly more than half of a minicircle was determined by Maxam-Gilbert sequence determination to be 490 base pairs (bp) in length. This fragment has dramatically anomalous electrophoretic behavior; it has an apparent size of 450 bp on a 1% agarose gel but migrates as 1,380 bp on a 12% polyacrylamide gel. However, in gel filtration on Sephacryl S-500, the fragment elutes with an apparent size of 375 bp. Finally, it behaves anomalously in electric dichroism experiments. Field-free rotational relaxation times from transient electric dichroism studies are highly sensitive to effective molecular dimensions. The rotational relaxation time of the kinetoplast fragment is smaller than that of a 309-bp control fragment from pBR322. Because rigorous control experiments rule out the possibility that this fragment is modified, these anomalous properties must be dictated by the sequence itself. Fragment behavior indicates that it has an unusually compact configuration; we propose that this molecule contains a region of systematically bent B-DNA. This model accounts for the fragment's difficulty in snaking through the pores of a polyacrylamide gel, its ease in diffusing into Sephacryl beads, and its smaller rotational relaxation time. Bending of this molecule may be caused by periodicities in the DNA sequence.

Predicting mutation frequencies in stem-loop structures of derepressed genes: implications for evolution

This work provides evidence that, during transcription, the mutability (propensity to mutate) of a base in a DNA secondary structure depends both on the stability of the structure and on the extent to which the base is unpaired. Zuker's DNA folding computer program reveals the most probable stem-loop structures (SLSs) and negative energies of folding (-DeltaG) for any given nucleotide sequence. We developed an interfacing program that calculates (i) the percentage of folds in which each base is unpaired during transcription; and (ii) the mutability index (MI) for each base, expressed as an absolute value and defined as -follows: MI = (% total folds in which the base is unpaired) x (highest -DeltaG of all folds in which it is unpaired). Thus, MIs predict the relative mutation or reversion frequencies of unpaired bases in SLSs. MIs for 16 mutable bases in auxotrophs, selected during starvation in derepressed genes, are compared with 70 background mutations in lacI and ebgR that were not derepressed during mutant selection. All the results are consistent with the location of known mutable bases in SLSs. Specific conclusions are: (i) Of 16 mutable bases in transcribing genes, 87% have higher MIs than the average base of the sequence analysed, compared with 50% for the 70 background mutations. (ii) In 15 of the mutable bases of transcribing genes, the correlation between MIs and relative mutation frequencies determined experimentally is good. There is no correlation for 35 mutable bases in the lacI gene. (iii) In derepressed auxotrophs, 100% of the codons containing the mutable bases are within one codon's length of a stem, compared with 53% for the background mutable bases in lacI. (iv) The data suggest that environmental stressors may cause as well as select mutations in derepressed genes. The implications of these results for evolution are discussed.

Single-strand-specific nucleases

Single-strand-specific nucleases are multifunctional enzymes and widespread in distribution. Their ability to act selectively on single-stranded nucleic acids and single-stranded regions in double-stranded nucleic acids has led to their extensive application as probes for the structural determination of nucleic acids. Intracellularly, they have been implicated in recombination, repair and replication, whereas extracellular enzymes have a role in nutrition. Although more than 30 single-strand-specific nucleases from various sources have been isolated till now, only a few enzymes (S1 nuclease from Aspergillus oryzae, P1 nuclease from Penicillium citrinum and nucleases from Alteromonas espejiana, Neurospora crassa, Ustilago maydis and mung bean) have been characterized to a significant extent. Recently, some of these enzymes have been cloned, their crystal structures solved and their interactions with different substrates have been established. The detection, purification, characteristics, structure-function correlations, biological role and applications of single-strand-specific nucleases are reviewed.

(CGA)(4): parallel, anti-parallel, right-handed and left-handed homoduplexes of a trinucleotide repeat DNA

Conformational properties of microsatellite DNA regions are the probable reason of their expansions in genomes which lead to serious genetic diseases in some cases. Using CD spectroscopy, UV absorption spectroscopy and polyacrylamide gel electrophoresis, we study in this paper conformational properties of (CGA)(4) and compare them with those of (CAG)(4) - a related repeat, connected with Huntington's disease. We show that (CGA)(4) can adopt several distinct conformations in solution. Around neutral pH it forms a parallel-stranded homoduplex containing C(+).C, G.G, and A.A base pairs. Under the same conditions (CAG)(4) forms a hairpin. At slightly alkaline pH values and low ionic strength, (CGA)(4) also folded into a hairpin which transformed into a bimolecular anti-parallel homoduplex at increasing salt concentrations. The duplex easily isomerized into left-handed Z-DNA, implying that the mismatched adenines between G.C pairs facilitate rather than hinder the B-Z transition. No similar changes took place with (CAG)(4). Thus, the conformational repertoire of (CGA)(4) includes parallel, anti-parallel, right-handed, and left-handed homoduplexes. In contrast, (CAG)(4) invariably adopts only a single conformation, namely the very stable hairpin.

Genes and the physics of the DNA double-helix

The processing of the genetic information stored in the double-helical DNA implies the separation of the two strands, the physics of which is described by the helix-coil transition model. Is there a relationship between genetic maps and DNA physical stability maps that plot the sequence-specific propensity for the thermal disruption of the double-helix? Here, with appropriate methodological formulations, such maps are derived for a large set of sequences, including complete genomes. The superposition of the two maps leads to a contrasted picture with correlations ranging between two extremes: from almost perfect (with the genes precisely delineated as stable regions) to more or less complete unrelatedness. The simplest explanation for the results is that the observed striking correlations correspond to the relics of a primeval organisation of the genetic message, with the physics of DNA playing a role in the delimitation of coding regions. In order to trace the evolutionary fate of this signal further, a detailed study of the yeast complete genome is performed. In this study, the superposition of the genetic and physical stability maps is examined in the light of information concerning gene duplication. On the basis of this analysis it is concluded that the 'signature' associated with the supposed archaic signal is in the process of being erased, most probably because the underlying feature is no longer under selective pressure. There are many evolutionary implications for the results presented and for their proposed interpretations, notably concerning models of mutational dynamics in relation to erasure processes.

Melting studies of short DNA hairpins: influence of loop sequence and adjoining base pair identity on hairpin thermodynamic stability

Spectroscopic and calorimetric melting studies of 28 DNA hairpins were performed. These hairpins form by intramolecular folding of 16 base self-complementary DNA oligomer sequences. Sequence design dictated that the hairpin structures have a six base pair duplex linked by a four base loop and that the first five base pairs in the stem are the same in every molecule. Only loop sequence and identity of the duplex base pair closing the loop vary for the set of hairpins. For these DNA samples, melting studies were carried out to investigate effects of the variables on hairpin stability. Stability of the 28 oligomers was ascertained from their temperature-induced melting transitions in buffered 115 mM Na(+) solvent, monitored by ultraviolet absorbance and differential scanning calorimetry (DSC). Experiments revealed the melting temperatures of these molecules range from 32.4 to 60.5 degrees C and are concentration independent over strand concentrations of 0.5 to 260 microM; thus, as expected for hairpins, the melting transitions are apparently unimolecular. Model independent thermodynamic transition parameters, DeltaH(cal), DeltaS(cal), and DeltaG(cal), were determined from DSC measurements. Model dependent transition parameters, DeltaH(vH), DeltaS(vH), and DeltaG(vH) were estimated from a van't Hoff (two-state) analysis of optical melting transitions. Results of these studies reveal a significant sequence dependence to DNA hairpin stability. Thermodynamic parameters evaluated by either procedure reveal the transition enthalpy, DeltaH(cal) (DeltaH(vH)) can differ by as much as 20 kcal/mol depending on sequence. Similarly, values of the transition entropy DeltaS(cal) (DeltaS(vH)) can differ by as much as 60 cal/Kmol (eu) for different molecules. Differences in free energies DeltaG(cal) (DeltaG(vH)) are as large as 4 kcal/mol for hairpins with different sequences. Comparisons between the model independent calorimetric values and the thermodynamic parameters evaluated assuming a two-state model reveal that 10 of the 28 hairpins display non-two-state melting behavior. The database of sequence-dependent melting free energies obtained for the hairpins was employed to extract a set of n-n (nearest-neighbor) sequence dependent loop parameters that were able to reproduce the input data within error (with only two exceptions). Surprisingly, this suggests that the thermodynamic stability of the DNA hairpins can in large part be reasonably represented in terms of sums of appropriate nearest-neighbor loop sequence parameters.

Circular dichroism spectroscopy analysis of conformational transitions of a 54 base pair DNA duplex composed of alternating CGCGCG and TATATA blocks

CD spectroscopy was used to analyze conformational properties of a self-complementary 54-mer DNA composed of alternating (CG)3 and (TA)3 hexamers. NaCl induced the Z-form in poly(GC), but the 54-mer remained the B-form under the same conditions. The B-Z transition was induced only after the addition of NiCl2. However, the Z-form was adopted by the whole molecule, not by the (CG)3 blocks alone. Two orders of magnitude higher concentrations of NiCl2 were required to induce the Z-form in poly(AT). The Z-form was also induced in poly(GC) by CsF that switched poly(AT) into the X-form, which seems to be a solution counterpart of D-DNA. Under these conditions the CD spectrum of the 54-mer was consistent with the (TA)3 blocks being in the X-form and the (CG)3 blocks in the B-form. At high concentrations of ethanol or trifluoroethanol, poly(AT) was an A-form, while poly(GC) adopted either Z-form, A-form or Z'-form. At the high trifluoroethanol concentrations the 54-mer cooperatively switched into a conformation whose CD spectrum was most consistent with the A-form in the (TA)3 blocks and the Z'-form in the (CG)3 blocks. This suggests that the base pairs are tilted in the Z'-form as in the A-form. The present article illustrates that CD spectroscopy can provide interesting pieces of information about conformational isomerizations and coexistence of multiple conformations in DNA molecules containing blocks of different simple sequence repeats.

Prokaryotic promoters in biotechnology

The basic properties of prokaryotic promoters and the promotor region are described with special emphasis on promoters that are found in Escherichia coli and Bacillus subtilis. Promoters recognized by major and minor forms of RNA polymerase holoenzymes are compared for their specificities and differences. Both natural and hybrid promoters that have been constructed for purposes of efficient and regulated transcription are discussed in terms of their utility. Since promoter regions contain sequences that are recognized not only by RNA polymerase but by positive and negative regulatory factors that regulate expression from promoters, the functions and properties of these promoter regions are also described. The current utility and the future prospects of the prokaryotic promoters in expressing heterologous genes for biotechnology purposes are discussed.

NMR determination of the conformational and drug binding properties of the DNA heptamer d(GpCpGpApApGpC) in aqueous solution

1D and 2D NMR spectroscopy (500/600 MHz) has been used to investigate the equilibrium conformational states of the deoxyheptanucleotide 5'-d(GpCpGpApApGpC), as well as its complexation with the phenanthridinium drug ethidium bromide (EB). Quantitative determination (reaction constants and thermodynamic parameters) of the conformational equilibrium of the heptamer in solution and its complexation with EB was based on analysis of the dependence of proton chemical shifts on concentration (at two temperatures, 298 and 308 K) and on temperature (in the range 278-353 K). The experimental results were analysed in terms of a model of the dynamic equilibrium between single-stranded, hairpin and bulged dimer forms of the deoxyheptanucleotide and its complexes with EB. Calculation of the relative amounts of the different complexes reveals important features of the dynamic equilibrium as a function of both temperature and the ratio of the drug and heptamer concentrations. The quantitative analysis also provides the limiting proton chemical shifts of EB in each complex which have been used to determine the most favourable structures of the intercalated complexes of EB with the (GC) sites of both the hairpin and dimer forms of the heptanucleotide.

Structure and dynamics of the DNA hairpins formed by tandemly repeated CTG triplets associated with myotonic dystrophy

Anomalous expansion of the DNA triplet (CTG)n causes myotonic dystrophy. Structural studies have been carried out on (CTG)n repeats in an attempt to better understand the molecular mechanism of repeat expansion. NMR and gel electrophoretic studies demonstrate the presence of hairpin structures for (CTG)5 and (CTG)6 in solution. The monomeric hairpin structure remains invariant over a wide range of salt concentrations (10-200 mM NaCl), DNA concentrations (micromolar to millimolar in DNA strand) and pH (6.0-7.5). The (CTG)n hairpin contains three bases in the loop when n is odd and four bases when n is even. For both odd and even n the stacking and pairing in the stem remain the same, i.e, two hydrogen bond T.T pairs stack with the neighboring G.C pairs. All the nucleotides in (CTG)5 and (CTG)6 adopt C2'-endo, anti conformations. Full-relaxation matrix analysis has been performed to derive the NOE distance constraints from NOESY experiments at seven different mixing times (25, 50, 75, 100, 125, 200 and 500 ms). NOESY-derived distance constraints were subsequently used in restrained molecular dynamics simulations to obtain a family of structures consistent with the NMR data. The theoretical order parameters are computed for H5-H6(cytosines) and H2'-H2" dipolar correlations for both (CTG)5 and (CTG)6 by employing the Lipari-Szabo formalism. Experimental data show that the cytosine in the loop of the (CTG)5 hairpin is slightly more flexible than those in the stem. The cytosine in the loop of the (CTG)6 hairpin is extremely flexible, implying that the dynamics of the four base loop is intrinsically different from that of the three base loop.

Thermodynamic investigation of the association of ethidium, propidium and bis-ethidium to DNA hairpins

We have used a combination of calorimetric and spectroscopic techniques to investigate the association of the bis-intercalator ethidium homodimer (bis-ethidium) to short DNA hairpins with sequences: d(GCGCT5GCGC) and d(CGCGT5CGCG). The helix-coil transition of each hairpin, investigated by UV and calorimetric melting protocol, takes place in monomolecular two-state transitions with characteristic enthalpies of approximately 37 kcal mol-1 for disrupting the four dG-dC base pairs of the hairpin stems. Deconvolution of the bis-ethidium-hairpin calorimetric titration curves indicate that each hairpin contains two distinct binding sites for the ligand: a high affinity site in the stem (Kb approximately 10(7)) that accommodates one bis-ethidium molecule and a lower affinity site (Kb approximately 10(6)) located probably at the loop that accommodates two bis-ethidium molecules. The overall stoichiometries of three ligands per hairpin are in agreement with those obtained in continuous variation experiments using visible spectroscopy. The interaction of bis-ethidium for each type of sites results in enthalpy driven reactions, with average binding enthalpies, delta Hb, of -13.1 and -12.1 kcal mol-1 for the stem and loop sites, respectively. Comparison to the thermodynamic profiles of ethidium and propidium binding reveals that the bis-ethidium binding to the stem site of each hairpin has a more favorable free energy term of -1.4 kcal mol-1 and more favorable enthalpy of -4.2 kcal mol-1. These suggest that only one phenanthridine ring of bis-ethidium intercalates in the stem, while the second planar ring is exposed to solvent or weakly associated to the surface of DNA.

Thermodynamics of DNA hairpins: contribution of loop size to hairpin stability and ethidium binding

A combination of calorimetric and spectroscopic techniques was used to evaluate the thermodynamic behavior of a set of DNA hairpins with the sequence d(GCGCTnGCGC), where n = 3, 5 and 7, and the interaction of each hairpin with ethidium. All three hairpins melt in two-state monomolecular transitions, with tm's ranging from 79.1 degrees C (T3) to 57.5 degrees C (T7), and transition enthalpies of approximately 38.5 kcal mol-1. Standard thermodynamic profiles at 20 degrees C reveal that the lower stability of the T5 and T7 hairpins corresponds to a delta G degree term of +0.5 kcal mol-1 per thymine residue, due to the entropic ordering of the thymine loops and uptake of counterions. Deconvolution of the ethidium-hairpin calorimetric titration curves indicate two sets of binding sites that correspond to one ligand in the stem with binding affinity, Kb, of approximately 1.8 x 10(6) M-1, and two ligands in the loops with Kb of approximately 4.3 x 10(4) M-1. However, the binding enthalpy, delta Hb, ranges from -8.6 (T3) to -11.6 kcal mol-1 (T7) for the stem site, and -6.6 (T3) to -12.7 kcal mol-1 (T7) for the loop site. Relative to the T3 hairpin, we obtained an overall thermodynamic contribution (per dT residue) of delta delta Hb = delta(T delta Sb) = -0.7(5) kcal mol-1 for the stem sites and delta delta Hb = delta(T delta Sb) = -1.5 kcal mol-1 for the loop sites. Therefore, the induced structural perturbations of ethidium binding results in a differential compensation of favorable stacking interactions with the unfavorable ordering of the ligands.

Studies of DNA dumbbells. II. Construction and characterization of DNA dumbbells with a 16 base-pair duplex stem and Tn end loops (n = 2, 3, 4, 6, 8, 10, 14)

The preparation and characterization of DNA dumbbells that contain the 16 base-pair duplex sequences 5'G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C3' (set 1) and 5'G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C3' (set 2) are reported. The dumbbells of set 1 have the duplex stem nucleated on both ends by Tn (n = 2, 3, 4, 6, 8, 10, and 14) loops. The dumbbells of set 2 have Tn (n = 2, 4, 8, 10) end loops. For the molecules of set 1, effects of end loop size on the electrophoretic mobility, CD and UV absorbance spectra, and cleavage by restriction enzymes, were investigated. Effects of loop size on the CD spectra and restriction enzyme cleavage of the molecules of set 2 were also examined. Optical melting curves of the molecules of set 1 were collected as a function of sodium ion concentration from 30 to 120 mM. These investigations revealed that as loop size decreases, the electrophoretic mobilities, rates of enzyme cleavage, and optical melting temperatures increase. For end loops with at least three T's the observed increases are inversely proportional to loop size. The behavior of the dumbbell with T2 end loops departs from this linear dependence and is anomalous in every experimental context. For molecules with end loops comprised of at least four T's CD spectra were virtually indistinguishable. However, these spectra differed considerably from the CD spectrum of the T2-looped molecule. The CD spectrum of the dumbbell with T3 end loops displayed features common to the dumbbells with larger loops and T2 end loops. Thermodynamic evidence that the terminal G.C base pairs (bps) nucleating the T2 end loops were intact was obtained from a comparison of the melting temperature of this molecule with that of a DNA dumbbell containing the 14 central bps of the set 1 duplex sequence linked instead by end loops comprised of the four base sequence, C-T-T-C. The tm of this latter molecule was determined to be 9 degrees C less than that of the former dumbbell assumed to contain a 16-bp stem and T2 end loops.(ABSTRACT TRUNCATED AT 400 WORDS)

The solution structure of a DNA hairpin containing a loop of three thymidines determined by nuclear magnetic resonance and molecular mechanics

We have determined by two-dimensional nuclear magnetic resonance studies and molecular mechanics calculations the three-dimensional solution structure of a 21 residue oligonucleotide capable of forming a hairpin structure with a loop of three thymidine residues. This structure is in equilibrium with a duplex form. At 33 degrees C, low ionic strength and in the presence of MgCl2 the hairpin form dominates in solution. Six Watson-Crick base pairs are formed topped by the loop structure. The residues 1-3 and 18-21 are not complementary and form dangling ends. Distance constraints have been derived from nuclear Overhauser enhancement measurements. These, together with molecular mechanics calculations, have been used to determine the structure. We do not observe stacking of thymidine residues either over the 3' or the 5' end of the stem.

rbcR [correction of rcbR], a gene coding for a member of the LysR family of transcriptional regulators, is located upstream of the expressed set of ribulose 1,5-bisphosphate carboxylase/oxygenase genes in the photosynthetic bacterium Chromatium vinosum

An open reading frame, rbcR, was identified 226 bp upstream of rbcAB, i.e., the ribulose 1,5-bisphosphate carboxylase genes expressed in the phototrophic purple bacterium Chromatium vinosum. Several features reveal that rbcR encodes a member of the LysR family of transcriptional regulators, in which an anomalous content of lysine and arginine residues (Lys/Arg anomaly) was found. The expression of rbcR in Escherichia coli as a protein fused to the N-terminal region of beta-galactosidase led to reduced expression of rbcAB. Thus, rbcR is likely to encode a trans-acting transcriptional regulator of rbcAB expression in C. vinosum.

Coupling of sequential transitions in a DNA double hairpin: energetics, ion binding, and hydration

In an effort to evaluate the relative contributions of sequence, ion binding, and hydration to the thermodynamic stability of nucleic acids, we have investigated the melting behavior of a double hairpin and that of its component single hairpins. Temperature-dependent UV absorption and differential scanning calorimetry techniques have been used to characterize the helix-coil transitions of three deoxyoligonucleotides: d(GTACT5GTAC), d(GCGCT5GCGC), and d(GCGCT5GCGCGTACT5GTAC). The first two oligomers melt with transition temperatures equal to 28 and 69 degrees C, respectively, in 10 mM dibasic sodium phosphate at pH 7.0. The Tm's are independent of strand concentration, strongly indicating the presence of single-stranded hairpin structures at low temperatures. The third oligomer, with a sequence corresponding to the joined sequences of the first two oligomers, melts with two apparently independent monomolecular transitions with Tm's of 41 and 69 degrees C. These transitions correspond to the melting of a double hairpin. In the salt range of 10-100 mM in NaCl, we obtain average enthalpies of 24 and 38 kcal/mol for the transitions in the single-hairpin molecules. Each transition in the double hairpin has an enthalpy of 32 kcal/mol. In addition, dtm/d log [Na+] for the transitions are 4.1 and 4.7 degrees C for the single hairpins and 12.6 and 11.2 degrees C for each transition in the double hairpin. The differential ion binding parameter between the double hairpin and that of the sum of single hairpins is roughly equal to 1.1 mol of Na+ ions/mol of double hairpin and is consistent with an increase in the electrostatic behavior of the stem phosphates of this molecule.

Analysis of melting transitions of the DNA hairpins formed from the oligomer sequences d[GGATAC(X)4GTATCC] (X = A, T, G, C)

Optical melting transitions of the short DNA hairpins formed from the self-complementary DNA oligomers d[GGATACX4GTATCC] where X = A, T, G, or C measured in 100 mM NaCl are presented. A significant dependence of the melting transitions on loop sequence is observed and transition temperatures, tm, of the hairpins vary from 58.3 degrees C for the T4 loop hairpin to 55.3 degrees C for the A4 loop. A nearest-neighbor sequence-dependent theoretical algorithm for calculating melting curves of DNA hairpins is presented and employed to analyze the experimental melting transitions. Experimental melting curves were fit by adjustment of a single theoretical parameter, Fend(n), the weighting function for a hairpin loop comprised of n single-strand bases. Empirically determined values of Fend(n) provide an evaluation of the free-energy of hairpin loop formation and stability. Effects of heterogeneous nearest-neighbor sequence interactions in the duplex stem on hairpin loop formation were investigated by evaluating Fend(n) in individual fitting procedures using two of the published sets of nearest-neighbor stacking interactions in DNA evaluated in 100 mM NaCl and given by Wartell and Benight, 1985. In all cases, evaluated values of Fend(n) were obtained that provided exact theoretical predictions of the experimental transitions. Results of the evaluations indicate: (1) Evaluated free-energies of hairpin loop formation are only slightly dependent on loop sequences examined. At the transition temperature, Tm, the free-energy of forming a loop of four bases is approximately equal for T4, G4, or C4 loops and varies from 3.9 to 4.8 kcal/mole depending on the set of nearest-neighbor interactions employed in the evaluations. This result suggests, in light of the observed differences in stability between the T4, G4, and C4 loop hairpins, that sequence-dependent interactions between base residues of the loop are most likely not the source of the enhanced stability of a T4 loop.(ABSTRACT TRUNCATED AT 400 WORDS)

B to Z transitions of the short DNA hairpins formed from the oligomer sequences: d[(CG)3X4(CG)3] (X = A, T, G, C)

Circular Dichroism (CD) spectra were collected as a function of sodium perchlorate concentration [NaClO4] for the set of DNA hairpins formed from the oligomer sequences d[(CG)3X4(CG)3] where X = A, T, G or C. Over the range in salt concentration from 0 to 4.0 M NaClO4, the CD spectra invert in a manner characteristic of the B to Z transition. A factor analysis routine is described and employed to determine the least number of basis spectra required to fit the measured spectra of each hairpin over the entire salt range examined. In every case, linear combinations of only two sub-spectra fit the experimental spectra of the hairpins with greater than 98% accuracy, indicating the spectrally monitored structural transitions are two-state. From the relative weights of the individual sub-spectra, B-Z transition curves are constructed. The transitions are analyzed in terms of a simple two-state equilibrium model which yields an evaluation of the transition free-energy, delta GB-Z, as a function of NaClO4 concentration. At 1.0 M NaClO4 and 21 degrees C, delta GB-Z = 5.4, 4.9, 3.6 and 2.3 kcal/mole for the G4, T4, A4 and C4 loop hairpins, respectively.

Thermodynamic stability of the 5' dangling-ended DNA hairpins formed from sequences 5'-(XY)2GGATAC(T)4GTATCC-3', where X, Y = A, T, G, C

Expressions for the partition function Q (T) of DNA hairpins are presented. Calculations of Q (T), in conjunction with our previously reported numerically exact algorithm [T. M. Paner, M. Amaratunga, M. J. Doktycz, and A. S. Benight (1990) Biopolymers, 29, 1715-1734], yield a numerical method to evaluate the temperature dependence of the transition enthalpy, entropy, and free energy of a DNA hairpin directly from its optical melting curve. No prior assumptions that the short hairpins melt in a two-state manner are required. This method is then applied in a systematic manner to investigate the stability of the six basepair duplex stem 5'-GGATAC-3' having four-base dangling single-strand ends with the sequences (XY)2, where X, Y = A, T, G, C, on the 5' end and a T4 loop on the 3' end. Results show that all dangling ends of the sample set stabilize the hairpin against melting. Increases in transition temperatures as great as 4.0 degrees C above the blunt-ended control hairpin were observed. The hierarchy of the hairpin transition temperatures is dictated by the identity of the first base of the dangling end adjoining the duplex in the order: purine greater than T greater than C. Calculated melting curves of every hairpin were fit to experimental curves by adjustment of a single parameter in the numerically exact theoretical algorithm. Exact fits were obtained in all cases. Experimental melting curves were also calculated assuming a two-state melting process. Equally accurate fits of all dangling-ended hairpin melting curves were obtained with the two-state model calculation. This was not the case for the melting curve of the blunt-ended hairpin, indicating the presence of a four-base dangling-end drives hairpin melting to a two-state process. Q (T) was calculated as a function of temperature for each hairpin using the theoretical parameters that provided calculated curves in exact agreement with the experimentally obtained optical melting curves. From Q (T), the temperature dependence of the transition enthalpy delta H, entropy delta S, and free energy delta G were calculated for every hairpin providing a quantitative assessment of the effects of dangling ends on hairpin thermodynamics. Comparisons of our results are made with those of the Breslauer group [M. Senior, R. A. Jones, and K. J. Breslauer (1988) Biochemistry 27, 3879-3885] on the T2 5' dangling-ended d(GC)3 duplexes.(ABSTRACT TRUNCATED AT 400 WORDS)

Inverted repeat sequences can influence the melting transitions of linear DNAs

The influence of inverted repeat sequences on the melting transitions of linear DNAs has been examined. Derivative melting curves (DMC) of a 514 base pair (bp) DNA, seven subfragments of this DNA, and four other DNAs have been compared to predictions of DNA melting theory. The 514-bp DNA contains three inverted repeat sequences that can form cruciform structures in supercoiled DNA. We refer to these sequences as c-inverted repeats. Previous work showed that the DMC of this DNA, unlike a number of other DNAs, is not accurately predicted by DNA melting theory. Since the theoretical model does not include hairpin-like structures, it was suggested that hairpin or cruciform formation in these inverted repeats may be responsible for this discrepancy. Our results support this hypothesis. Predicted DMCs are in good agreement with DNAs with no inverted repeats, or inverted repeats not evident in supercoiled DNA. Differences between the theoretical and experimental Tm's are less than or equal to 0.3 degrees C. DNA molecules that contain one or more of the three c-inverted repeats are not as accurately predicted. Experimental Tm values are lower than predicted values by 0.7-3.8 degrees C. It is concluded that some inverted repeat sequences can form hairpin-like structures during the melting of linear DNAs. These structures appear to lower overall DNA stability.

Hairpin formation of d(CGCG-TA-CGCG), d(CGCG-TG-CGCG) and their cytosine methylated analogs

Hairpin formations of decamers d(CGCG-TA-CGCG), d(CGCG-TG-CGCG), and their m5dC analogs are evidenced by the existence of biphasic absorbance melting profiles in which the lower transition temperature increases with increasing oligomer concentration, whereas the higher melting temperature is concentration independent. The corresponding temperature dependent CD intensity at 285 nm exhibits a maximum around 55 degrees C. These observations are consistent with the interpretation that the lower temperature transition corresponds to the duplex to hairpin transformation while the melting of hairpins into single strands constitutes the higher temperature transition. The CD spectrum of the hairpin conformation appears to be characterized by a couplet with nearly equal positive and negative intensities at 285 and 255 nm, respectively, while a significantly smaller intensity at 285 nm is apparent for the duplex form. The hairpin conformation is suspected to contain a two-nucleotide loop. Titrations with NaCl further suggest that, in contrast to the TA sequence, the TG sequence with wobble base pairing favors Z formation under high salt conditions.

Conformation and dynamics of a left-handed Z-DNA hairpin: studies of d(CGCGCGTTTTCGCGCG) in solution

The physical properties of the DNA oligomer d(CGCGCGTTTTCGCGCG) in solvents containing 4 M NaClO4 and 0.1 M NaCl were investigated by proton NMR, optical melting, and circular dichroism spectroscopy. Results of these investigations are as follows: (i) The DNA hexadecamer exists as a unimolecular hairpin in either high or low salt. (ii) In high salt the stem region of the hairpin is in the left-handed Z conformation. (iii) In either high or low salt, the duplex stem of the hairpin is stabilized against melting by approximately 40 degrees C compared to the linear core duplex. The added stability of the hairpin is entropic in origin. (iv) In high salt, as the temperature is elevated, the equilibrium structure of the duplex stem of the hairpin shifts from the Z to the B conformation before melting. (v) In low salt, when the DNA duplex exists in the B conformation, attachment of a T4 single-strand loop to one end only slightly decreases (by 14%) the correlation time of the CH5-CH6 interproton vector. In high salt, when the DNA duplex exists in the Z conformation, the correlation time of the CH5-CH6 interproton vector decreases by 51%. Since these viscosity-corrected correlation times are taken to be indicators of duplex motions on the nanosecond time scale, this result directly suggests a larger amplitude of these motions is present in the duplex stem of the hairpin when it exists in the Z conformation.

Immunochemistry of DNA

Since the first reports of anti-DNA antibodies in sera of patients with systemic lupus erythematosus (SLE) in 1957, studies of nucleic acid immunochemistry have grown in two directions. One has been the analysis of the specificity, the nature and the origins of these autoantibodies. The second has been exploration of anti-nucleic acid antibodies that can be induced experimentally, their specificities, and their application as biochemical reagents. Although the properties of autoantibodies and experimentally induced antibodies differ in certain respects, these two lines of research are complementary and provide important information for each other. For example, the production of autoantibodies by adjuvant-stimulated B cells yields a background that has to be considered in evaluating the specificity of weak responses to experimental nucleic acid immunogens: in turn, the possibilities and limitations of experimental immunization should be considered in evaluating possible stimuli for autoantibody production. Several aspects of nucleic acid immunochemistry have been described and evaluated in previous reviews. Following some general statements of historical perspective, this review will emphasize questions addressed and findings of about the last five years.

DNA conformational variations in the in vitro torsionally strained Ig kappa light chain gene localize on consensus sequences

We have analyzed the localization and the dependence upon superhelical density of the DNA sites which modify their conformation under torsional strain in a mouse Ig L kappa gene. The conformational variations occur on DNA sites which have been defined as protein interaction sites and consensus sequence motifs: the 5'-upstream regulatory decanucleotides, the TATA sequence, the consensus heptanucleotides of the J recombinational sequences.

The B- to Z-DNA equilibrium in vivo is perturbed by biological processes

Right-handed B and left-handed Z conformations coexist in equilibrium in portions of plasmids in Escherichia coli. The equilibria are influenced by the length of the sequences that undergo the structural transitions and are perturbed by biological processes. The composite results of three types of determinations indicate a supercoil density of -0.025 in vivo. The coexistence of alternative DNA conformations in living cells implies the potential of these structures or their transitions for important functions in genetic regulatory processes.

Influence of loop residues on the relative stabilities of DNA hairpin structures

We have determined the relative stabilities and melting behaviors of DNA hairpin structures as a function of the nonbonded residues in the loop. The specific family of hairpin structures we investigated in this work is formed by the 16-mer sequence d[CGAACG(X)4CGTTCG], where X is deoxyadenosine, deoxycytidine, deoxyguanosine, or deoxythymidine. As shown below, this 16-mer can fold back on itself to form a family of DNA hairpin structures that possess a common hexameric stem duplex and a nonbonded loop of 4 nucleotides. For the hairpin structures investigated in this work, we varied the loop composition from all purine residues to all pyrimidine residues. (Formula: see text). We thermodynamically characterized the relative stabilities and melting profiles of these hairpin structures by a combination of spectroscopic and calorimetric techniques. To establish a thermodynamic "baseline," we also conducted parallel studies on the isolated hexameric duplex d[CGAACG).(CG-TTCG)], which corresponds to the common stem duplex present in each hairpin structure. Our spectroscopic and calorimetric data reveal the following: (i) The hairpin structure with four dT residues in the loop exhibits the highest melting temperature, while the corresponding hairpin structure with four dA residues in the loop exhibits the lowest melting temperature. (ii) The free energy data at 25 degrees C reveal the following order of DNA hairpin stability for the four structures studied here: T loop greater than C loop greater than G loop greater than A loop. In other words, the pyrimidine-looped hairpins of four residues are more stable than the purine-looped hairpins. (iii) The loop-dependent order of hairpin stability is paralleled by a similar trend in the calorimetrically determined transition enthalpies for hairpin disruption. Thus, the enhanced stability of the pyrimidine-looped hairpin structures relative to purine-looped hairpin structures is enthalpic in origin. To develop insight into the molecular basis for the thermodynamic differences, proton NMR spectroscopy was used to probe for structural disparities between the most stable hairpin structure (T loop) and the least stable hairpin structure (A loop). Two-dimensional nuclear Overhauser enhancement spectroscopy revealed connectivities between the residues in the stem duplexes of both hairpin structures that are consistent with B-form DNA. In addition, the nonbonded residues in both the T and A loops exhibited the same connectivity patterns. However, on the 5' side of the stem-loop junction, the T-loop residue exhibited a connectivity with the adjacent base pair of the stem duplex that is not observed for the corresponding A-loop residue. This difference in connectivities at the stem-loop junction may provide a structural basis for our observation that the T-looped hairpin structure is more stable than the corresponding A-looped hairpin structure.

Oligodeoxynucleotide folding in solution: loop size and stability of B-hairpins

The secondary structures of the synthetic DNA fragments d(CGCGCGTTTTTCGCGCG) (T5), d(CGCGCGAAAAACGCGCG) (A5), d(CGCGCGTACGCGCG) (TA), and d(CGCGCGATCGCGCG) (AT) were investigated in a combined electrophoretic and spectroscopic study. All the oligomers exist, at low temperature and over a wide range of ionic strength (0.5-100 mM salt) and of nucleotide concentration [0.1-2.0 mM (phosphate)], as a mixture of two slowly interconverting species, identified as the dimeric duplex and the monomeric hairpin structure. The thermodynamic parameters for hairpin denaturation of T5, A5, TA, and AT and for duplex denaturation of d(CGCGCG) show that (a) the hairpins are more stable than the reference hexamer duplex at all accessible nucleotide concentrations; (b) the loop contributes favorably to the enthalpy change of hairpin denaturation in the four DNA fragments; (c) the base composition of the loop (A vs T) and the size of the loop (A5/T5 vs TA/AT) do not appreciably influence the enthalpic contents of the hairpins; (d) hairpins TA and AT, with two AT bases intervening in the CG self-complementary part of the molecule, exhibit a markedly higher thermal stability than hairpins T5 and A5, which is entropic in origin. These findings are consistent with the presence of two-residue loops in the tetradecamers TA and AT.

The B-Z conformational transition in folded oligodeoxynucleotides: loop size and stability of Z-hairpins

The capacity to assume a left-handed conformation and the thermodynamics of loop formation in concentrated aqueous NaClO4 have been investigated for the following palindromic sequences: d-(CGCGCGAAAAACGCGCG) (A5), d(CGCGCGTTTTTCGCGCG) (T5), d(CGCGCGTACGCGCG) (TA), and d(CGCGCGATCGCGCG) (AT). The results show that (a) each oligomer assumes a Z conformation upon exposure to increasing NaClO4 concentrations; the salt concentration at the transition midpoint is 1.8 M for both A5 and T5 and 3 and 3.5 M for TA and AT, respectively; (b) in high salt the four oligomers exist, over a wide range of nucleotide concentrations (up to 10(-3) M) and of temperature (greater than 0 degrees C), as unimolecular hairpin structures; (c) hairpins TA and AT exhibit, in buffer A, a lower thermal stability with respect to A5 and T5 (delta T about 16 degrees C), contrary to what is observed at low ionic strength; (d) on hairpin formation, the enthalpic term is about -52 kcal/mol for the two 17-mers and -38 kcal/mol for the two 14-mers, while the change in entropy is found to be around -150 eu for A5 and T5 and -115 eu for TA and AT. This thermodynamic picture suggests that a two-residue loop for TA and AT, found at low ionic strength [see preceding paper (Xodo, L.E., Manzini, G., Quadrifoglio, F., van der Marel, G.A., & van Boom, J.H. (1988) Biochemistry (preceding paper in this issue)], is substituted by a longer one including two additional residues from a missing dC.dG base pairing at the top of the stem.

Stability of RNA stem-loop structure and distribution of non-random structure in the human immunodeficiency virus (HIV-I)

The stability of potential RNA stem-loop structures in human immunodeficiency virus isolates, HTLV-III and ARV, has been calculated, and the relevance to the local significant secondary structures in the sequence has been tested statistically using a Monte Carlo simulation method. Potentially significant structures exist in the 5'non-coding region, the boundary regions between the protein coding frames, and the 3' non-coding region. The locally optimal secondary structure occurring in the 5' terminal region has been assessed using different overlapping segment sizes and the Monte Carlo method. The results show that the most favorable structure for the 5' mRNA leader sequence of HIV has two stem-loops folded at nucleotides 5-104 in the R region (stem-loop I, 5-54 and stem-loop II, 58-104). A large fluctuation of segment score of the local optimal secondary structure also occurs in the boundary between the exterior glycosylated protein or outer membrane protein and transmembrane protein coding region. This finding is surprising since no RNA signals or RNA processing are expected to occur at this site. In addition, regions of the genome predicted to have significantly more open structure at the RNA level correlate closely with hypervariable sites found in these viral genomes. The possible importance of local secondary structure to the biological function of the human immunodeficiency virus genome is discussed.

Nonlinear waves in DNA and regulation of transcription

The transmission of conformational changes along the DNA double helix is interpreted as a propagation of nonlinear solitary waves. Experimental and theoretical data in favour of this approach are presented. The possible role of nonlinear waves in regulation of transcription is discussed.

Left-handed DNA in vivo

Left-handed DNA is shown to exist and elicit a biological response in Escherichia coli. A plasmid encoding the gene for a temperature-sensitive Eco RI methylase (MEco RI) was cotransformed with different plasmids containing inserts that had varying capacities to form left-handed helices or cruciforms with a target Eco RI site in the center or at the ends of the inserts. Inhibition of methylation in vivo was found for the stable inserts with the longest left-handed (presumably Z) helices. In vitro methylation with the purified MEco RI agreed with the results in vivo. Supercoil-induced changes in the structure of the primary helix in vitro provided confirmation that left-handed helices were responsible for this behavior. The presence in vivo of left-handed inserts elicits specific deletions and plasmid incompatibilities in certain instances.

Unusual protonated structure in the homopurine.homopyrimidine tract of supercoiled and linearized plasmids recognized by chemical probes

Plasmid pEJ4, which is a derivative of pUC19 containing an insert with 60-bp-long homopurine.homopyrimidine tract from sea urchin P. miliaris histone gene spacer, was studied by chemical probes of the DNA structure osmium tetroxide and glyoxal. The former probe reacts with pyrimidine bases, while the latter forms a stable product only with guanine residues. These probes can thus be applied as specific probes for the homopyrimidine and homopurine strands, respectively. At pH 6.0 the site-specific modification of the homopurine.homopyrimidine tract by both probes was observed at native superhelical density of the plasmid. In the linear plasmid under the same conditions this modification was absent; it appeared, however, at more acid pH values. In supercoiled DNA the hypersensitivity of the homopurine.homopyrimidine tract to osmium tetroxide did not substantially change when pH was decreased from 6.0 to 4.0. Changes in NaCl concentration at pH 4.5 did not influence the hypersensitivity to osmium tetroxide; at pH 6.0 this hypersensitivity decreased with increasing NaCl concentration. These results thus show that the chemical probes recognize an unusual protonated structure containing unpaired bases or non-Watson-Crick base pairs. At pH 5.6 the site-specific modification occurred at or near to the middle of the homopurine.homopyrimidine tract, suggesting that a hairpin may be involved in the unusual structure under the given conditions. From the models suggested so far for the unusual structure of homopurine.homopyrimidine tracts our results fit best the protonated triplex H form suggest by V.I. Lyamichev, S.M. Mirkin and M.D. Frank-Kamenetskii, J. Biomol. Struct. Dyn. 3,667 (1986).

Conformational and model-building studies of the hairpin form of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G)

The hairpin form of the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G) was studied by means of NMR spectroscopy. In a companion study it is shown that the hairpin form of this DNA fragment consists of a structure with a stem of three Watson-Crick-type base pairs and a loop consisting of only two nucleotides. The non-exchangeable proton resonances were assigned by means of two-dimensional correlation spectroscopy and two-dimensional nuclear Overhauser effect spectroscopy. Proton-proton coupling constants were used for the conformational analysis of the deoxyribose ring and for some of the backbone torsion angles. From the two-dimensional NMR spectra and the coupling-constant analysis it is concluded that: (i) the stem of the hairpin exhibits B-DNA characteristics; (ii) the sugar rings are not conformationally pure, but display a certain amount of conformational flexibility; (iii) the stacking interaction in the stem of the hairpin is elongated from the 3'-side in a more or less regular fashion with the two loop nucleotides; (iv) at the 5'-side of the stem a stacking discontinuity occurs between the stem and the loop; (v) at the 5'-side of the stem the loop is closed by means of a sharp backbone turn which involves unusual gamma' and beta+ torsion angles in residue dG(6). The NMR results led to the construction of a hairpin-loop model which was energy-minimized by means of a molecular-mechanics program. The results clearly show that a DNA hairpin-loop structure in which the loop consists of only two nucleotides bridging the minor groove in a straightforward fashion, (i) causes no undue steric strain, and (ii) involves well-known conformational principles throughout the course of the backbone. The hairpin form of the title compound is compared with the hairpin form of d(A-T-C-C-T-A-T4-T-A-G-G-A-T), in which the central -T4- part forms a loop of four nucleotides. Both models display similarities as far as stacking interactions are concerned.