Locations of frequently opening regions on natural DNAs and their relation to functional loci

The shaping of a molecular linguist: How a career studying DNA energetics revealed the language of molecular communication

My personal and professional journeys have been far from predictable based on my early childhood. Owing to a range of serendipitous influences, I miraculously transitioned from a rebellious, apathetic teenage street urchin who did poorly in school to a highly motivated, disciplined, and ambitious academic honors student. I was the proverbial “late bloomer.” Ultimately, I earned my PhD in biophysical chemistry at Yale, followed by a postdoc fellowship at Berkeley. These two meccas of thermodynamics, coupled with my deep fascination with biology, instilled in me a passion to pursue an academic career focused on mapping the energy landscapes of biological systems. I viewed differential energetics as the language of molecular communication that would dictate and control biological structures, as well as modulate the modes of action associated with biological functions. I wanted to be a “molecular linguist.” For the next 50 years, my group and I used a combination of spectroscopic and calorimetric techniques to characterize the energy profiles of the polymorphic conformational space of DNA molecules, their differential ligand-binding properties, and the energy landscapes associated with mutagenic DNA damage recognition, repair, and replication. As elaborated below, the resultant energy databases have enabled the development of quantitative molecular biology through the rational design of primers, probes, and arrays for diagnostic, therapeutic, and molecular-profiling protocols, which collectively have contributed to a myriad of biomedical assays. Such profiling is further justified by yielding unique energy-based insights that complement and expand elegant, structure-based understandings of biological processes.

Relating DNA base-pairing in aqueous media to DNA polymerase fidelity

Controversy surrounds the perceived absence of a relationship between DNA polymerase fidelity (kinetic discrimination) and free energy changes determined from DNA melting studies (thermodynamic discrimination). Thermodynamic discrimination together with aqueous solvent effects can account for kinetic fidelities on the order of those observed experimentally.

Opening the Strands of Replication Origins-Still an Open Question

The local separation of duplex DNA strands (strand opening) is necessary for initiating basic transactions on DNA such as transcription, replication, and homologous recombination. Strand opening is commonly a stage at which these processes are regulated. Many different mechanisms are used to open the DNA duplex, the details of which are of great current interest. In this review, we focus on a few well-studied cases of DNA replication origin opening in bacteria. In particular, we discuss the opening of origins that support the theta (θ) mode of replication, which is used by all chromosomal origins and many extra-chromosomal elements such as plasmids and phages. Although the details of opening can vary among different origins, a common theme is binding of the initiator to multiple sites at the origin, causing stress that opens an adjacent and intrinsically unstable A+T rich region. The initiator stabilizes the opening by capturing one of the open strands. How the initiator binding energy is harnessed for strand opening remains to be understood.

Prediction of melting temperatures in fluorescence in situ hybridization (FISH) procedures using thermodynamic models

The thermodynamics and kinetics of DNA hybridization, i.e. the process of self-assembly of one, two or more complementary nucleic acid strands, has been studied for many years. The appearance of the nearest-neighbor model led to several theoretical and experimental papers on DNA thermodynamics that provide reasonably accurate thermodynamic information on nucleic acid duplexes and allow estimation of the melting temperature. Because there are no thermodynamic models specifically developed to predict the hybridization temperature of a probe used in a fluorescence in situ hybridization (FISH) procedure, the melting temperature is used as a reference, together with corrections for certain compounds that are used during FISH. However, the quantitative relation between melting and experimental FISH temperatures is poorly described. In this review, various models used to predict the melting temperature for rRNA targets, for DNA oligonucleotides and for nucleic acid mimics (chemically modified oligonucleotides), will be addressed in detail, together with a critical assessment of how this information should be used in FISH.

Analysis of mutational spectra by denaturing capillary electrophoresis

The point mutational spectrum over nearly any 75- to 250-bp DNA sequence isolated from cells, tissues or large populations may be discovered using denaturing capillary electrophoresis (DCE). A modification of the standard DCE method that uses cycling temperature (e.g., +/-5 degrees C), CyDCE, permits optimal resolution of mutant sequences using computer-defined target sequences without preliminary optimization experiments. The protocol consists of three steps: computer design of target sequence including polymerase chain reaction (PCR) primers, high-fidelity DNA amplification by PCR and mutant sequence separation by CyDCE and takes about 6 h. DCE and CyDCE have been used to define quantitative point mutational spectra relating to errors of DNA polymerases, human cells in development and carcinogenesis, common gene-disease associations and microbial populations. Detection limits are about 5 x 10(-3) (mutants copies/total copies) but can be as low as 10(-6) (mutants copies/total copies) when DCE is used in combination with fraction collection for mutant enrichment. No other technological approach for unknown mutant detection and enumeration offers the sensitivity, generality and efficiency of the approach described herein.

DnaA: controlling the initiation of bacterial DNA replication and more

Escherichia coli is a model system to study the mechanism of DNA replication and its regulation during the cell cycle. One regulatory pathway ensures that initiation of DNA replication from the chromosomal origin, oriC, is synchronous and occurs at the proper time in the bacterial cell cycle. A major player in this pathway is SeqA protein and involves its ability to bind preferentially to oriC when it is hemi-methylated. The second pathway modulates DnaA activity by stimulating the hydrolysis of ATP bound to DnaA protein. The regulatory inactivation of DnaA function involves an interaction with Hda protein and the beta dimer, which functions as a sliding clamp for the replicase, DNA polymerase III holoenzyme. The datA locus represents a third mechanism, which appears to influence the availability of DnaA protein in supporting rifampicin-resistant initiations.

Dissection of the bacteriophage Mu strong gyrase site (SGS): significance of the SGS right arm in Mu biology and DNA gyrase mechanism

The bacteriophage Mu strong gyrase site (SGS), required for efficient phage DNA replication, differs from other gyrase sites in the efficiency of gyrase binding coupled with a highly processive supercoiling activity. Genetic studies have implicated the right arm of the SGS as a key structural feature for promoting rapid Mu replication. Here, we show that deletion of the distal portion of the right arm abolishes efficient binding, cleavage, and supercoiling by DNA gyrase in vitro. DNase I footprinting analysis of the intact SGS revealed an adenylyl imidodiphosphate-dependent change in protection in the right arm, indicating that this arm likely forms the T segment that is passed through the cleaved G segment during the supercoiling reaction. Furthermore, in an SGS derivative with an altered right-arm sequence, the left arm showed these changes, suggesting that the selection of a T segment by gyrase is determined primarily by the sequences of the arms. Analysis of the sequences of the SGS and other gyrase sites suggests that the choice of T segment correlates with which arm possesses the more extensive set of phased anisotropic bending signals, with the Mu right arm possessing an unusually extended set of such signals. The implications of these observations for the structure of the gyrase-DNA complex and for the biological function of the Mu SGS are discussed.

Rapid diagnosis of clonal immunoglobulin heavy chain gene rearrangements in cutaneous B-cell lymphomas using the LightCycler-Polymerase Chain Reaction with DNA melting curve analysis

We have recently developed a novel Immunoglobulin heavy chain gene rearrangement (IgH-R) assay that combines polymerase chain reaction (PCR) amplification and analysis in the same closed capillary tube using the LightCycler System. IgH-R can be identified by DNA melting curve analysis within 40 minutes after DNA preparation and amplification. To test the clinical utility of this new IgH-R assay for rapidly diagnosing cutaneous B-cell lymphomas, we prospectively analyzed 44 formalin-fixed, paraffin-embedded tissues suspected of B-cell malignant lymphoma: skin (n = 31), lymph node (n = 7), stomach (n = 3), spleen (n = 1), colon (n = 1), and soft tissue (n = 1). We detected IgH-R in 12 DNA samples, including 8 skin biopsies, with the following diagnoses: B-cell chronic lymphocytic leukemia (n = 4), extranodal marginal zone B-cell lymphoma (n = 4), diffuse large B-cell lymphoma (n = 2), Burkitt lymphoma (n = 1), and precursor B-lymphoblastic lymphoma (n = 1). DNA melting curve analysis, compared with polyacrylamide gel electrophoresis, achieved a sensitivity equal to 92.3% and a specificity equal to 100%. There was a single false negative result because DNA melting curve analysis could not detect less than 10.0% clonal B-cells. We conclude that this new, rapid PCR assay for detecting IgH-R based on DNA melting curve analysis can be clinically useful for confirming the initial diagnosis of B-cell malignant lymphoma.

Re-replication from non-sequesterable origins generates three-nucleoid cells which divide asymmetrically

In rapidly growing Escherichia coli cells replication cycles overlap and initiation occurs at multiple replication origins (oriCs). All origins within a cell are initiated essentially in synchrony and only once per cell cycle. Immediate re-initiation of new origins is avoided by sequestration, a mechanism dependent on the SeqA protein and Dam methylation of GATC sites in oriC. Here, GATC sites in oriC were changed to GTTC. This reduced the sequestration to essentially the level found in SeqA-less cells. The mutant origins underwent re-initiation, showing that the GATC sites in oriC are required for sequestration. Each re-initiation eventually gave rise to a cell containing an extra nucleoid. The three-nucleoid cells displayed one asymmetrically placed FtsZ-ring and divided into a two-nucleoid cell and a one-nucleoid cell. The three nucleoid-cells thus divided into three daughters by two consecutive divisions. The results show that extra rounds of replication cause extra daughter cells to be formed prematurely. The fairly normal mutant growth rate and size distribution show, however, that premature rounds of replication, chromosome segregation, and cell division are flexibly accommodated by the existing cell cycle controls.

Cytosine methylation and CpG, TpG (CpA) and TpA frequencies

An analysis of dinucleotide frequencies was carried out on DNAs from insects and mammals, as well as on large DNA sequences from the genomes of Drosophila melanogaster, Anopheles gambiae, puffer fish (Takifugu rubripes), zebra fish (Danio rerio) and human. These organisms were chosen because Drosophila and Anopheles DNAs have an extremely low level of methylation, human DNA a high level and fish DNA a two-fold higher level compared to human. The results indicate that: (i) CpG deficiency and the corresponding TpG (CpA) excess show no correlation with the level of DNA methylation; indeed, genomes endowed with strikingly different levels of DNA methylation (such as those of Drosophila and human) exhibited similar TpG (CpA) levels; (ii) the correlation between GC levels of large (50 kb) DNA sequences and TpA or CpG shortage levels do not appear to be due to CpG methylation followed by deamination; (iii) CpG dinucleotides are more frequent in fishes than in mammals; interestingly, the monotreme Ornitorhinchus anatinus shows an intermediate CpG frequency. The implications of these results are discussed.

Non-random features of loop-size chromatin fragmentation

Upon isolation of DNA from normal eukaryotic cells by standard methods involving extensive proteolytic treatment, a rather homogeneous population of loop-size, double-stranded DNA fragments is regularly obtained. These DNA molecules can be efficiently end-labeled by the DNA polymerase I Klenow fragment, as well as by a 3'- to -5'-exonuclease-free Klenow enzyme, but not by terminal transferase (TdT) unless the ends have been filled up by Klenow, suggesting that dominantly 5' protruding termini are generated upon fragmentation. The filled-up termini were used for cloning the distal parts of the approximately 50 kb fragments. BLAST analysis of the sequence of several clones allowed us to determine the sequence of the non-cloned side of the breakpoints. Comparison of 25, 600 bp-long breakpoint sequences demonstrated prevalence of repetitive elements. Consensus motives characteristic of the breakpoint sequences have been identified. Several sequences exhibit peculiar computed conformational characteristics, with sharp transition or center of symmetry located exactly at the breakpoint. Our data collectively suggest that chromatin fragmentation involves nucleolytic cleavages at fragile/hypersensitive sites delimiting loop-size fragments in a non-random manner. Interestingly, the sequence characteristics of the breakpoints are reminiscent of certain breakpoint cluster regions frequently subject to gene rearrangements.

Rapid and accurate detection of monoclonal immunoglobulin heavy chain gene rearrangement by DNA melting curve analysis in the LightCycler System

The detection of immunoglobulin heavy chain gene rearrangement (IgH-R) is a standard tool for distinguishing polyclonal from monoclonal B-cell populations. Current DNA-based polymerase chain reactions (PCR) strategies can diagnose monoclonal IgH-R either by measuring the length of the amplicon or by detecting gel mobility variations owing to sequence-dependent conformational changes. However, amplification and analysis remain sequential operations usually requiring manual transfer. We have developed a novel PCR strategy for detecting monoclonal IgH-R that monitors fluorescence of the specific double-stranded DNA binding dye SYBR Green I during melting curve analysis using the LightCycler System. We compared polyacrylamide gel electrophoresis (PAGE) versus melting curve analysis in 130 clinical DNA samples from formalin-fixed, paraffin-embedded (FFPE) tissues (mostly skin biopsies) of 128 patients. The identical FR3 primers were used to amplify the IgH variable region for both analytic techniques. We detected IgH-R in 24 DNA samples from FFPE tissue of 22 patients. Melting curve analysis, compared to PAGE, revealed no false negative and no false positive results, yielding both sensitivity and specificity equal to 100%. We also compared Southern blot analysis versus melting curve analysis in 23 clinical DNA samples from fresh-frozen lymph nodes of 23 patients. We detected IgH-R by melting curve analysis in 7 DNA samples from fresh-frozen lymph nodes. Melting curve analysis, compared to Southern blot analysis, revealed sensitivity equal to 58.3% (7 of 12) and specificity equal to 100% (11 of 11). We conclude that continuous fluorescence monitoring of PCR products with DNA melting curve analysis can rapidly and reproducibly distinguish polyclonal from monoclonal B-cell populations.

The immediate early gene of canine herpesvirus is transcribed through early and late phases

The immediate early (IE) gene of canine herpesvirus (CHV), homologue of the infected cell protein 4 (ICP4) gene of herpes simplex virus 1, is transcribed as a 4.9kb mRNA during IE phase. The IE gene was further transcribed as a 4.8kb mRNA through early (E) and late (L) phases of productive infection. Transcription of the 4.8kb mRNA initiated from downstream of the TATA box in an intron which was spliced out during IE phase. The reverse transcription-polymerase chain reaction revealed that the IE promoter was turned off during L phase at a permissive temperature. We, thus, propose to redesignate the IE gene of CHV as CICP4 gene.

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.

The N-terminal methionine is a major determinant of the DNA binding specificity of MEF-2C

Members of the MEF-2 family of transcriptional regulators positively modulate the activity of basic helix-loop-helix proteins in both myogenic and neurogenic cell lineages. Previous work had shown that MEF-2C(2-117), a protein fragment comprising the dimerization and DNA-binding domains of MEF-2C but lacking the N-terminal methionine, bound to AT-rich DNA sequences with high affinity. MEF-2C(2-117) did not discriminate between different AT-rich sequences. We now report the in vitro DNA binding properties of a MEF-2C fragment containing the N-terminal methionine. Measurements of the apparent dissociation constants of the complexes of GG-MEF-2C(1-117) revealed that different AT-rich sequences are bound with different affinities; in particular MEF site containing DNA (CTATAAATAG) is bound preferentially to DNA containing a SRF site (CATAAATG). Strikingly, when the shorter AT run consisted of six alternating thymines and adenines, almost wild-type affinity was observed. Irrespective of the particular DNA sequence, all circular dichroism spectra of the DNA complexes of GG-MEF-2C(1-117) were superimposable and characterized by an identical maximal ellipticity at 269.5 nm, suggesting similar DNA conformations. Bending analysis by circular permutation assay revealed that on complex formation MEF-2C(2-117) induced cognate DNA to bend by 49 degrees, while heterologous DNA remained unbent. In the presence of the N-terminal methionine, however, all DNA sequences were bent by 70 degrees. The above results suggest an important function for the N-terminal methionine in properly orientating MEF-2C on the DNA.

Rapid, comprehensive analysis of human cytokine mRNA and its application to the study of acute renal allograft rejection

Cytokine mRNA analysis was performed on human renal allograft needle core biopsies by a PCR-based assay. The assay was specifically developed to be capable of simultaneous analysis of multiple interleukin transcripts (IL-1-IL-12), as well as those of other relevant cytokines, by one person in less than 1 day from cultured cells or directly from tissue samples. It was initially used on preparations containing known amounts of plasmid DNA encoding individual cytokine cDNA sequences, confirming that the sensitivity of this technique was both well defined and comparable for all target sequences tested. Analysis of human PBLs prior to stimulation, after polyclonal stimulation with PHA and after simultaneous treatment with PHA and MP or CyA, was also performed to show a proportional relationship between mRNA levels measured by PCR and protein release measured by ELISA (R2 = 0.86). This correlation was not adversely altered by pharmacologic immunosuppression by MP or CyA. Thus, this method of PCR primer design and usage was appropriate for the clinical study of cytokine mRNA levels during allograft rejection. Direct study of cytokine mRNA in allograft biopsy tissue showed that IL-2 was specifically and significantly (p = 0.006) elevated during ACR when compared to other causes of graft dysfunction. Transcripts from the IFN-gamma and IL-6 genes were also increased in ACR (p = 0.001 and 0.017, respectively), whereas increased IL-8 mRNA was correlated with irreversible loss of graft function (p = 0.02). TNF-alpha, IL-1 beta, and IL-10 gene transcripts were also detected during ACR, but were not quantitatively increased compared to other forms of graft injury (p > 0.2). We conclude that acute cellular rejection is associated with intragraft mRNA from the IL-2 gene. Other transcripts, including those from the IFN-gamma, IL-6, and IL-8 genes, are detected in increased amounts during this process. Messenger RNA from the TNF-alpha, IL-1 beta and IL-10 genes is also detected during ACR, but the presence of these transcripts is not exclusive to this process.

DNA spatial considerations in the arrangement of G/C and A/T blocks

Intensive computations have been carried out on eukaryotic and prokaryotic DNA sequences present in the GenBank database. These calculations were aimed at studying particular sequence patterns. Previously we have observed a trend in the relative positioning of DNA oligomers with respect to each other. Specifically, we have shown that there is a preference for A/T stretches to be inserted inside G/C blocks, partitioning the latter. This arrangement is preferred over having a longer G/C block with an A/T stretch next to it. That is (G/C)n(A/T)m(G/C)2 greater than (G/C)n + 2(A/T)m. Previously we have attributed this preferred pattern to nucleosome packaging of the DNA. Since the average length of the DNA involved in a single nucleosome formation is 200 bp, oligomers were scored in the analysis only if two identical ones occurred within that distance. Since, in addition, the total length of the oligomers studied is n + m + 2 less than or equal to 7, the counts of the longer oligomers were quite low. Here we have repeated the analysis on the newer, larger database, with that restriction removed. Our new analysis confirms and strengthens the older findings. Moreover, we have carried out DNA structural analysis of these trends, focusing on the twist and roll angular parameters. We have obtained these values from detailed structural computations on the Cray super computer. Some correlations between these values, DNA flexibility and the sequence trends are observed. We suggest that the positioning of the (A/T)m sequence stretch between the two (G/C) blocks affords greater flexibility in the spatial positioning of the G/C sequence elements.

The effect of temperature and oligonucleotide primer length on the specificity and efficiency of amplification by the polymerase chain reaction

The polymerase chain reaction (PCR) is most effectively performed using a thermostable DNA polymerase such as that isolated from Thermus aquaticus. Since temperature and oligonucleotide length are known to control the specificity of oligonucleotide hybridization, we have investigated the effect of oligonucleotide length, base composition, and the annealing temperature on the specificity and efficiency of amplification by the PCR. Generally, the specificity of PCR is controlled by the length of the oligonucleotide and/or the temperature of annealing of the primer to the template. An empirical relationship between oligonucleotide length and ability to support amplification was determined. This relationship allows for the design of specific oligonucleotide primers. A model is proposed which helps explain the observed dependence of PCR on annealing temperature and length of the primer.

Premelting thermal fluctuational base pair opening probability of poly(dA).poly(dT) as predicted by the modified self-consistent phonon theory

We employ a mean field, modified, self-consistent phonon theory to evaluate the single base-pair opening rate and the probability of a base pair in the amino proton exchangeable state for the homopolymer poly(dA).poly(dT) at temperatures below the helix-coil transition region. Our calculated premelting single base-pair opening probabilities are in general agreement with several available experimental estimates from imino proton exchange and formaldehyde-induced DNA melting measurements. These calculated opening probabilities, however, are in disagreement with the prediction of the helix-coil transition theory. Possible reasons for the differences are discussed, especially the possible different definition of a meaningful open state in the premelting region. The premelting open state of the modified self-consistent phonon approximation theory seems to be appropriate to describe a solvent-accessible open configuration that is sufficient to facilitate important chemical reactions such as imino proton exchange and formaldehyde reaction with the bases. This can be compared with the completely unstacked open state of the helix-coil transition theory originally defined in the helix-coil transition region. We propose that the amino proton exchangeable state is different from the open state associated with melting and only involves the breaking of the amino interbase H bond. The agreement between the calculated and experimentally estimated probability of a base pair in the amino proton exchangeable state seems to support this hypothesis.

Long range and symmetry considerations in the DNA

The common description of DNA is that of a polymer composed of overlapping dinucleotide base pairs. Conformational, thermodynamic and flexibility calculations of the DNA are frequently based on that description. As demonstrated by the good fit with experimental data, such as the free energy of opening the double stranded DNA and the correspondence between the DNA crystal and computed structures, this approach is to a large extent a faithful reproduction of the DNA. Here I show that longer range effects should be considered as well. Statistically strong longer-range sequence patterns are described. Such sequence preferences are clearly translated in the DNA into long range structural propagation. Long range effects have also been observed in some experimental studies, like gel mobility and patterns of DNA cleavage. Long range structural effects are the likely explanation of the effect of mutations at a distance from a protein binding site. They might also aid in understanding DNA looping.

Long-range structural effects in supercoiled DNA: statistical thermodynamics reveals a correlation between calculated cooperative melting and contextual influence on cruciform extrusion

We have previously described [K. M. Sullivan and D. M. J. Lilley (1986) Cell 47, 817-827] a set of sequences, called C-type inducing sequences, which cause cruciform extrusion by adjacent inverted repeats to occur by an abnormal kinetic pathway involving a large denatured region of DNA. In this paper we apply statistical thermodynamic DNA helix melting theory to these sequences. We find a marked correlation between the ability of sequences to confer C-type cruciform character experimentally and their calculated propensity to undergo cooperative melting, and no exceptions have been found. The correlations are both qualitative and quantitative. Thus the ColE1 flanking sequences behave as single melting units, while the DNA of the S-type plasmid pIRbke8 exhibits no propensity to melt in the region of the bke cruciform. The results of the calculations are also fully consistent with the following experimental observations: 1. the ability of the isolated colL and colR fragments of the ColE1 flanking sequences, as well as the short sequence col30, to confer C-type character; 2. C-type induction by an A + T rich Drosophila sequence; 3. low-temperature cruciform extrusion by an (AT)34 sequence; 4. the effect of changing sequences at a site 90 base pairs (bp) removed from the inverted repeat; 5. the effects of systematic deletion of the colL sequence; and 6. the effects of insertion of various sequences in between the colL sequence and the xke inverted repeat. These studies show that telestability effects on thermal denaturation as predicted from equilibrium helix melting theory of linear DNA molecules may explain all the features that are revealed by studying the extrusion of cruciforms in circular DNA molecules subjected to superhelical stress.

Integration of cellular-SV40 DNA and the thermal stability of the recombinant sites. Brief report

We calculated the predicted melting profiles of evolutionary variants of SV40 and compared the nucleotide sequences of large T antigen binding sites and recombinant sites for homology. Neither is critically important in selecting a recombinant site for SV40-SV40 and SV40-cellular DNA fragments.

An NMR study of polymorphous behaviour of the mismatched DNA octamer d(m5C-G-m5C-G-A-G-m5C-G) in solution. The B-duplex and hairpin forms

By means of one- and two-dimensional NMR spectroscopy the solution structures of the partly self-complementary octamer d(m5C-G-m5C-G-A-G-m5C-G) were investigated. It is shown that this DNA fragment, under conditions of high DNA concentration (8 mM DNA) and/or high ionic strength prefers to adopt a duplex structure. At low DNA concentration (0.4 mM DNA), the duplex exists in a 1:1 slow equilibrium with a monomeric hairpin form. Addition of salt destabilizes the hairpin structure in favour of the dimer. At high temperatures the hairpin form, as well as the dimer structure, exist in a fast equilibrium with the random-coil form. For the hairpin/random-coil equilibrium a Tm of 329 K and a delta H degree of -121 kJ.mol-1 were deduced. These thermodynamic parameters are independent of the DNA concentration, as is expected for a monomeric structure. For the dimer to coil transition a Tm of 359 K (1 M DNA) and a delta H degree of -285 kJ.mol duplex-1 were derived. The thermodynamic data of the hairpin-coil transition mutually agree with those recently reported for the hairpin to random coil equilibrium of the DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) [Orbons, L. P. M., van der Marel, G. A., van Boom, J. H. & Altona, C. (1987) J. Biomol. Struct. Dyns. 4, 939-963]. It is demonstrated that the dimer structure exhibits B-DNA characteristics, as is witnessed by the NOESY experiments and the analysis of the proton-proton coupling data. It is shown that the base-pair formation of the G x A mismatches is anti-anti. A comparison of 1H and 31P chemical-shift data of the title compound with those of a well-characterized B-DNA structure reveals large differences in the dm5C(3)-dG(4)-dA(5) part of the mismatched dimer structure. These differences apparently indicate some major local structural changes due to the incorporation of the G x A mismatches. Under the most extreme conditions used (i.e. up to 3 M NaCl or 75% CH3OH in the presence of 10 mM MgCl2) no Z-DNA structure was observed. It is shown that the structural features of the hairpin form of the title compound mimic those of the hairpin structure of d(m5C-G-m5C-G-T-G-m5C-G). An energy-minimized model of the hairpin form is given.

Influence of the base sequence on the conformational behaviour of DNA polynucleotides in solution

NMR studies were carried out on samples of the non-self-complementary tetramers d(C-A-C-A), d(T-G-T-G), d(G-A-G-A) and d(T-C-T-C) and of 1:1 mixtures of the complementary tetramers d(C-A-C-A) X d(T-G-T-G) and d(G-A-G-A) X d(T-C-T-C) at two DNA concentrations and of the self-complementary octamers d(C-A-C-A-T-G-T-G) and d(G-A-G-A-T-C-T-C). Assignments, based upon one-dimensional NOE and homonuclear-decoupling and two-dimensional correlated and NOE spectroscopies are given of the resonances of most of the base and sugar protons. Chemical shift vs temperature profiles, constructed for all samples, yielded insight into the temperature- and concentration-dependent conformational behaviour of the compounds and were used to obtain thermodynamic parameters pertaining to the stacked-single-strand----random-coil and duplex----random-coil equilibria. Vicinal proton-proton couplings were analyzed in terms of the conformation of the deoxyribose rings in the single-stranded tetramers and duplexed octamers. The NOE patterns, chemical shift profiles, imino-proton resonances and coupling data revealed that the compounds adopt B-DNA-like structures. The ratio duplexed/stacked-single-strand/random coil depends upon external conditions as well as upon base sequence. The thermodynamic data indicate that: in terms of single-helical stacking, the R-R steps (Tm 321-328 K) appear more stable than the Y-R or R-Y steps (Tm 308-316 K) and the Y-Y steps score least (Tm 290-300 K), and the duplexes consisting of alternating, d(Y-R)n, strands are more stable, in terms of delta H degrees, compared to the d(R-R)n X d(Y-Y)n duplexes. The analyses of the couplings demonstrated that the sugars of the single-stranded tetramers and duplexed octamers occur as a blend of N- and S-type conformers, with a preference for the S-type (C2'-endo) sugar conformation: upon duplex formation, no significant shift in the N-type/S-type ratio was observed. The fraction S-type sugar conformation of a given residue, %S, in the stacked-single strands was found to depend upon the nature of its own base and that of the adjacent residues: sugars in an R-R stretch display high values of %S (90-100), whereas those in Y-Y stretches show relatively low values (approximately equal to 65).(ABSTRACT TRUNCATED AT 400 WORDS)

Relationship between helix-coil transition and gene organization of ColE1 plasmid DNA. Differential scanning calorimetric and theoretical studies

Differential scanning calorimetry (DSC) was carried out to analyze the transition of helix to coil state of DNA, using ColE1 DNA molecules digested with EcoRI. The DSC curves showed multimodal transition, consisting of nine to 11 peaks over a temperature range, depending on the ionic strength of the DNA solution. These DSC curves were essentially in good agreement with the optical melting curves of ColE1 DNA. The theoretical melting profiles of ColE1 DNA were predicted from calculations based on the helix-coil transition theory and the nucleotide sequence of the DNA. These profiles resembled the DSC curves and made it possible to assign the peaks seen in the DSC curves to the helix-coil transition of particular regions of the nucleotide sequence of ColE1. The helix-coil transition of each of the small genes gave rise to a single peak in the DSC curve, while the helix-coil transition of large genes contributed to two or more peaks in the DSC curve. This multimodal transition within a single coding region might correspond to the melting of individual segments encoding the different domains of the proteins. The helix-coil transition at the specific sites including ori, the origin of replication of ColE1, was also found to occur in a particular temperature range. DSC, a simple method, is thus useful for analyzing the multimodal helix-coil transition of DNA, and for providing information on the genetic organization of DNA.

Influence of N6-methylation of residue A(5) on the conformational behaviour of d(C-C-G-A-A-T-T-C-G-G) in solution studied by 1H-NMR spectroscopy. 1. The duplex form

One- and two-dimensional NMR studies at 300 MHz and 500 MHz were carried out on the two oligonucleotides d(C-C-G-A-A-T-T-C-G-G) and d(C-C-G-A-m6A-T-T-C-G-G) in aqueous solution. NMR spectra were observed at 10 mM sample concentration over the temperature range 273-368 K. Assignments are given of the base, H1', H2', H2", H3' and of some H4' resonances, based upon a combination of two-dimensional correlation spectra (COSY) and two-dimensional nuclear Overhauser effect spectra (NOESY); imino-proton resonances were assigned with the aid of a two-dimensional NOE experiment. Chemical shift vs temperature profiles were constructed in order to gain insight into the influence of N6-methylation of residue A(5) on the temperature-dependent conformational behaviour of the decamer and to determine thermodynamic parameters for the duplex-to-coil transition. The NOESY spectra, the imino-proton spectra and the shift profiles of the two compounds, under conditions where each forms a B-DNA-type duplex, are very similar. This is taken to indicate that the influence of N6-methylation of residue A(5) on the local structure of the duplex must be small. However, the temperature dependence of the (non-)exchangeable proton resonances of the two compounds reveals that methylation slows down the duplex-single-strand exchange. Furthermore, a thermodynamic analysis of the two compounds indicates that N6-methylation slightly decreases the stability of the duplex relative to the monomeric forms (Tm is reduced from 332 K down to 325 K at 10 mM sample concentration). Proton-proton couplings were obtained by means of one-dimensional and two-dimensional NMR experiments and were used in a conformational analysis of the sugar ring of each residue of the two compounds in the duplex form. The analysis indicated that all sugar rings display conformational flexibility in the intact duplex: population S-type sugar conformation ranges from 70% to 100%. A more refined analysis of the sugar rings of the parent compound revealed a sequence-dependent variation of the sugar geometry. This variation does not follow well the trend predicted by the Calladine/Dickerson sigma 3-sum rule [Dickerson, R. E. (1983) J. Mol. Biol. 166, 419-441; Calladine, C. R. (1982) J. Mol. Biol. 161, 343-352]; moreover the actual variations appear to be smaller in solution than those expected on the basis of known X-ray structures.

Dissociation kinetics of 19 base paired oligonucleotide-DNA duplexes containing different single mismatched base pairs

The dissociation kinetics of 19 base paired oligonucleotide-DNA duplex containing a various single mismatched base pair are studied on dried agarose gels. The kinetics of the dissociation are first order under our experimental conditions. The incorporation of a single mismatched base pair destabilizes the DNA duplexes to some extent, the amount depending on the nature of the mismatched base pair. G-T and G-A mismatches slightly destabilize a duplex, while A-A, T-T, C-T and C-A mismatches significantly destabilize it. The activation energy for the overall dissociation processes for these oligonucleotide-DNA duplexes containing 19 base pairs is 52 +/- 2 Kcal mol-1 as determined from the slope of Arrhenius plot.

Conformational analysis of the octamer d(G-G-C-C-G-G-C-C) in aqueous solution. A one-dimensional and two-dimensional proton NMR study at 300 MHz and 500 MHz

Proton NMR studies in H2O and 2H2O were carried out on the self-complementary DNA octamer d(G-G-C-C-G-G-C-C) and compared with similar studies on the dimethylated analogue d(G-G-m5C-m5C-G-G-C-C) [Sanderson, M. R., Mellema, J.-R., van der Marel, G. A., Wille, G., van Boom, J. H. & Altona, C. (1983) Nucleic Acids Res. 11, 3333-3346]. Base, H1', H2' and H2" resonances were assigned by means of two-dimensional nuclear Overhauser spectroscopy (NOESY) and correlated spectroscopy (COSY) experiments. Chemical shift vs temperature profiles were used to analyze the temperature-dependent conformational behaviour and to extract thermodynamic parameters for the duplex-to-coil transition. Analysis of proton-proton couplings were used to discriminate between J1'2' and J1'2" and between the H2' and H2" resonances as well as to obtain conformational parameters of the sugar rings. From the NOESY and COSY experiments, the temperature profiles and the analysis of the coupling data it is concluded that: (a) the compound adopts a B-DNA-type helix in solution; (b) the sugar rings in the intact duplex display limited conformational freedom; (c) methylation of the cytosine bases at the C5 position has no measurable effect on the conformational behaviour of the octamer, nor on the conformation of the sugar rings; however, methylation does increase the stability of the duplex in aqueous solution under conditions of low salt concentration.

Correlation between the rate of productive transcription initiation and the strand-melting property of Escherichia coli promoters

"Opening potential"s of DNA segments of about ten base pairs in length were calculated for eleven promoters of Escherichia coli using the thermodynamic parameter values [Gotoh and Tagashira (1981) Biopolymers 20, 1043-1058] for the stabilities of ten kinds of nearest neighbor base pair doublets in a melting reaction. They were compared with the strength of each promoter determined experimentally with a "mixed transcription" system [Kajitani and Ishihama (1983) Nucleic Acids Res. 11, 3873-3888]. A positive correlation was found between the calculated opening potential in the -9 to +3 region (+1 denotes the nucleotide position at which the transcription starts) and the rate of open complex formation.

Stability distribution in the phage lambda-DNA double helix: a correlation between physical and genetic structure

Statistical analyses on the positional correlation of physical-stability and base-sequence distribution maps with genetic map are made for the whole DNA (48502 bases) of lambda-phage. The susceptibility to a double-helix unfolding perturbation and the fraction of the transient opening of a particular region of the double helix are adopted to define this physical stability. The principal features obtained are: A) The DNA double strand of protein coding regions is found to have homostabilizing propensity around a defined stability which is characteristic to each individual gene. B) The stability of the double helix in non-protein coding region fluctuates, on average over the whole region, more than that in protein coding region. C) Boundary regions of protein coding and non-protein coding regions are regions of high stability-fluctuation. Stability especially fluctuates at the protein-coding-region side of the boundary. Contrary to the quiet feature of the interior part of protein coding region rather noisy part exists at its edge. D) One frequently opening region coincides with the attaching site for the site specific recombination between phage and bacterial DNA. There are two possible ways to explain the noisy feature in the stability distribution in non-protein coding regions: 1) The region has been used as the locus of recombination as evolution took place. Thus DNAs which were homostabilized around a different value characteristic to each individual DNA, have been joined there many times, so that the noise has accumulated as a remnant of evolutional history; and/or 2) the base-composition homogenizing or double-helix homostabilizing mechanism does not work in unneeded region such as non-protein coding region or introns. Since corresponding characteristics have been found in our previous analyses on other viral and globin-gene DNAs, the rules mentioned above may be comprehensively extended to other DNAs.

Allowance for heterogeneous stacking in the DNA helix-coil transition theory

Melting profiles of eight DNA molecules with lengths ranging from 849 to 4362 bp have been measured in an SSC buffer where the melting is an equilibrium process up to complete strand separation. A theoretical analysis shows that the melting profiles depend on only eight invariants that are linear combinations of 10 original stacking parameters. As a result it is impossible to determine the 10 parameters from the melting profiles. The 8 variants have been determined by fitting theoretical profiles to experimental ones for two fragments. Then theoretical and experimental profiles are compared for those 6 fragments that were not used in the fitting procedure. This comparison demonstrates that allowance for heterogeneous stacking considerably improves the agreement between theory and experiment. The values of invariants have proved to be small. This confirms the previous conclusion that heterogeneous stacking interactions produce only small corrections to the major effect of the difference in the mean stabilities of AT and GC pairs. Some discrepancy between theory and experiment that remains after the allowance for heterogeneous stacking is probably due to even finer effects of long-range interactions.

Dual-mode computer processing for high resolution DNA thermal denaturation experiments

Two modes of data processing are appropriate in conducting high resolution thermal denaturation experiments (thermal increments of 0.05 degrees or closer). In the first mode, a general purpose microcomputer provides on-line services important to the control and monitoring of the initial experiment, including control of the spectrophotometer and heater, the recording of data, and the display of current hyperchromicities and approximate first derivatives. A subsequent microcomputer program then reads the recorded data files and carries out accurate calculations of derivative denaturation profiles and the estimate of the statistical error of the first derivative at each point. The data collection program handles three samples at a time and was designed to provide optimal results in thermal denaturation experiments with a single-beam spectrophotometer.

Viability of palindromic DNA is restored by deletions occurring at low but variable frequency in plasmids of Escherichia coli

Palindromic arrangements of 2 X 400 to 2 X 1517 bp of nucleotide sequences were created by in vitro manipulation of Escherichia coli plasmids. As a consequence of its method of formation, each palindrome possessed at its center the recognition site for a particular restriction endonuclease. Eight out of eight palindromes, having at their centers the recognition sites for BamHI, BglI, BglII, HindIII, PstI, SalI, XhoI, and Xma III, were shown to be inviable when transformed into E. coli. The smallest of these palindromes had a half-length of approx. 400 bp. The lethality of palindromic sequences for their carrier plasmids was circumvented at low frequencies by spontaneous in vivo deletion events which removed the centers of symmetry of the palindromes. The frequencies of such deletions were less than 1% in most cases, but varied significantly both with the palindromic sequence in question and with the surrounding nonpalindromic sequences of the carrier plasmid. We confirmed the viability of a plasmid with a 147-bp palindrome [Bergsma et al., Gene 20 (1982) 157-167] and found that this palindrome (derived from SV40) does not confer viability on the plasmids with long palindromes.

Assignment of non-exchangeable base proton and H1' resonances of a deoxyoctanucleoside heptaphosphate d(G-G-C*-C*-G-G-C-C) by using the nuclear Overhauser effect

The resonances of the non-exchangeable base protons and 1' protons of the octamer d(G-G-C*-C*-G-G-C-C), C* = m5dC, have been assigned by means of NOE difference NMR spectroscopy at 500 MHz. From the measured J1'2' and J1'2" it follows that the octamer at low temperature prefers to adopt a B-DNA double-helical conformation in solution, however, some residual conformational freedom is detected at the 3' terminus. From the chemical shift versus temperature profiles it is concluded that no major conformational change occurs below 60-65 degrees C where the duplex formation for residues (2) to (6) is essentially completed under the conditions used.

DNA fragments differing by single base-pair substitutions are separated in denaturing gradient gels: correspondence with melting theory

DNA fragments 536 base pairs long differing by single base-pair substitutions were clearly separated in denaturing gradient gel electrophoresis. Transversions as well as transitions were detected. The correspondence between the gradient gel measurements and the sequence-specific statistical mechanical theory of melting shows that mutations affecting final gradient penetration lie within the first cooperatively melting sequence. Fragments carrying substitutions in domains melting at a higher temperature reach final gel positions indistinguishable from wild type. The gradient data and the sites of substitution bracket the boundary between the first domain and its neighboring higher-melting domain within eight base pairs or fewer, in agreement with the calculated boundary. The correspondence between the gradient displacement of the mutants and the calculated change in helix stability permits substantial inference as to the type of substitution. Excision of the lowest melting domain allows recognition of mutants in the next ranking domain.

Correlation of thermodynamic and genetic properties in the Tn10 encoded TET gene control region

The thermal stability of the Tn10 encoded tetracycline resistance (TET) gene control region is investigated by melting studies using purified DNA restriction fragments containing various amounts of flanking sequences. In order to study the thermodynamic properties of this control region under conditions, where enough flanking DNA is present to mimic the situation in the chromosome, the five step melting process of a 1450-bp DNA fragment is analyzed. Because most of the sequence of this DNA is not known, the assignment of the melting transitions to segments of the DNA is done by an experimental method. This employs the preparation of subfragments from the 1450-bp DNA and comparison of their denaturation profiles with the one of the intact sequence. This approach results in the complete assignment of the five denaturation steps. Rather than from the ends, the unwinding starts from the TET gene control region in the middle of the 1450-bp sequence. A clear correlation between the thermodynamic and genetic properties of this DNA is observed. The regulatory sequence forms a small cooperative unit with the lowest stability in the entire fragment. The thermal denaturation of the TET repressor. TET operator complex reveals, that the TET repressor specifically recognizes the double stranded TET operator DNA and stabilizes this structure by 2.4 degrees C. This results is also discussed as an example of the possible action of denaturing or stabilizing proteins on this genetic control region.