PMC1474062.xml

<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Archiving and Interchange DTD v2.3 20070202//EN" "archivearticle.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" article-type="research-article"><?properties open_access?><front><journal-meta><journal-id journal-id-type="nlm-ta">Nucleic Acids Res</journal-id><journal-id journal-id-type="publisher-id">Nucleic Acids Research</journal-id><journal-title>Nucleic Acids Research</journal-title><issn pub-type="ppub">0305-1048</issn><issn pub-type="epub">1362-4962</issn><publisher><publisher-name>Oxford University Press</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">16738142</article-id><article-id pub-id-type="pmc">PMC1474062</article-id><article-id pub-id-type="doi">10.1093/nar/gkl397</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Structural analysis of hyperperiodic DNA from <italic>Caenorhabditis elegans</italic></article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Moreno-Herrero</surname><given-names>Fernando</given-names></name></contrib><contrib contrib-type="author"><name><surname>Seidel</surname><given-names>Ralf</given-names></name></contrib><contrib contrib-type="author"><name><surname>Johnson</surname><given-names>Steven M.</given-names></name><xref ref-type="aff" rid="au1">1</xref><xref ref-type="aff" rid="au2">2</xref></contrib><contrib contrib-type="author"><name><surname>Fire</surname><given-names>Andrew</given-names></name><xref ref-type="aff" rid="au1">1</xref><xref ref-type="aff" rid="au2">2</xref></contrib><contrib contrib-type="author"><name><surname>Dekker</surname><given-names>Nynke H.</given-names></name><xref ref-type="corresp" rid="cor1">*</xref></contrib><aff><institution>Kavli Institute of Nanoscience, Faculty of Applied Sciences, Delft University of Technology</institution><addr-line>Lorentzweg 1, 2628 CJ Delft, The Netherlands</addr-line></aff><aff id="au1"><sup>1</sup><institution>Department of Pathology, Stanford University School of Medicine</institution><addr-line>300 Pasteur Drive, Room L235, Stanford, CA 94305-5324, USA</addr-line></aff><aff id="au2"><sup>2</sup><institution>Department of Genetics, Stanford University School of Medicine</institution><addr-line>300 Pasteur Drive, Room L235, Stanford, CA 94305-5324, USA</addr-line></aff></contrib-group><author-notes><corresp id="cor1"><sup>*</sup>To whom correspondence should be addressed. Tel: +31 0 15 278 3219; Fax: +31 0 15 278 1202; Email: <email>Nynke.Dekker@mb.tn.tudelft.nl</email></corresp><fn><p>Present address: Ralf Seidel, BioTechnological Center, University of Technology Dresden, Tatzberg 47-51, D-01307 Dresden, Germany</p></fn></author-notes><!--For NAR: both ppub and collection dates generated for PMC processing 1/27/05 beck--><pub-date pub-type="collection"><year>2006</year></pub-date><pub-date pub-type="ppub"><year>2006</year></pub-date><pub-date pub-type="epub"><day>31</day><month>5</month><year>2006</year></pub-date><volume>34</volume><issue>10</issue><fpage>3057</fpage><lpage>3066</lpage><history><date date-type="received"><day>11</day><month>4</month><year>2006</year></date><date date-type="rev-recd"><day>11</day><month>5</month><year>2006</year></date><date date-type="accepted"><day>11</day><month>5</month><year>2006</year></date></history><copyright-statement>&#x000a9; 2006 The Author(s)</copyright-statement><copyright-year>2006</copyright-year><license license-type="openaccess"><p>This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by-nc/2.0/uk/"/>) which permits unrestricted non-commerical use, distribution, and reproduction in any medium, provided the original work is properly cited.</p></license><abstract><p>Several bioinformatics studies have identified an unexpected but remarkably prevalent &#x0223c;10 bp periodicity of AA/TT dinucleotides (hyperperiodicity) in certain regions of the <italic>Caenorhabditis elegans</italic> genome. Although the relevant <italic>C.elegans</italic> DNA segments share certain sequence characteristics with bent DNAs from other sources (e.g. trypanosome mitochondria), the nematode sequences exhibit a much more extensive and defined hyperperiodicity. Given the presence of hyperperiodic structures in a number of critical <italic>C.elegans</italic> genes, the physical characteristics of hyperperiodic DNA are of considerable interest. In this work, we demonstrate that several hyperperiodic DNA segments from <italic>C.elegans</italic> exhibit structural anomalies using high-resolution atomic force microscopy (AFM) and gel electrophoresis. Our quantitative analysis of AFM images reveals that hyperperiodic DNA adopts a significantly smaller mean square end-to-end distance, hence a more compact coil structure, compared with non-periodic DNA of similar length. While molecules remain capable of adopting both bent and straight (rod-like) configurations, indicating that their flexibility is still retained, examination of the local curvatures along the DNA contour length reveals that the decreased mean square end-to-end distance can be attributed to the presence of long-scale intrinsic bending in hyperperiodic DNA. Such bending is not detected in non-periodic DNA. Similar studies of shorter, nucleosome-length DNAs that survived micrococcal nuclease digestion show that sequence hyperperiodicity in short segments can likewise induce strong intrinsic bending. It appears, therefore, that regions of the <italic>C.elegans</italic> genome display a significant correlation between DNA sequence and unusual mechanical properties.</p></abstract></article-meta></front><body><sec><title>INTRODUCTION</title><p>DNA structures are frequently modeled as flexible, isotropic rods both in solution and in their ability to contribute to native biological structures such as chromatin (<xref ref-type="bibr" rid="b1">1</xref>). Despite the utility of such models, it has been clear for some time that the sequence of a DNA molecule can contribute very specifically to its 3D structure, both <italic>in vitro</italic> (<xref ref-type="bibr" rid="b2">2</xref>,<xref ref-type="bibr" rid="b3">3</xref>) and <italic>in vivo</italic> (<xref ref-type="bibr" rid="b4">4</xref>). Particularly striking in their non-linear structure, certain sequences can introduce a stereo-specific bend in DNA. The most prominent sequence contributing to these structural anomalies are strings of three, four or more consecutive A or T residues (A-tracts) (<xref ref-type="bibr" rid="b5">5</xref>,<xref ref-type="bibr" rid="b6">6</xref>), although other nucleotide segments can also induce structural anomalies (<xref ref-type="bibr" rid="b7">7</xref>). A-tract bending has been investigated using X-ray crystallography, gel electrophoresis, electric birefringence, circular dichroism and NMR [for reviews see (<xref ref-type="bibr" rid="b8">8</xref>,<xref ref-type="bibr" rid="b9">9</xref>)]. DNA curvature of non-genomic DNA has also been investigated with electron microscopy (<xref ref-type="bibr" rid="b10">10</xref>,<xref ref-type="bibr" rid="b11">11</xref>) and atomic force microscopy (AFM) (<xref ref-type="bibr" rid="b12">12</xref>&#x02013;<xref ref-type="bibr" rid="b14">14</xref>).</p><p>Many of the unusual DNA sequences that have been the subject of biophysical and modeling studies have involved relatively short (10&#x02013;80 bases) bent regions either examined in isolation (<xref ref-type="bibr" rid="b15">15</xref>) or examined in the context of a longer fragment (<xref ref-type="bibr" rid="b13">13</xref>). This type of DNA structure is common in many genome sequences and in numerous cases has been associated with biological function, either in facilitating a specific protein interaction at the site of the bend as in the case of nucleosomes (<xref ref-type="bibr" rid="b16">16</xref>) or in facilitating a long-range interaction of two sites in the genome that otherwise might not be capable of interacting (<xref ref-type="bibr" rid="b17">17</xref>).</p><p>Recently, it has become clear that certain eukaryotic genomes include sequences with a more extended anomalous character in which groups of A<italic><sub>n</sub></italic>/T<italic><sub>n</sub></italic> residues are present in a periodic and highly non-random manner over a range of several hundred nucleotides (<xref ref-type="bibr" rid="b18">18</xref>,<xref ref-type="bibr" rid="b19">19</xref>). In the <italic>Caenorhabditis elegans</italic> genome, a remarkable periodicity of &#x0223c;10 bp AA/TT dinucleotides has been described elsewhere (<xref ref-type="bibr" rid="b20">20</xref>&#x02013;<xref ref-type="bibr" rid="b24">24</xref>). One proposed role for such sequences has been to limit the capability of nucleosomes to slide into the kind of ordered arrays needed for heterochromatin formation (<xref ref-type="bibr" rid="b19">19</xref>). Although such long DNAs can be readily modeled using the extensive data available for shorter bent segments, it is by no means clear that such models continue to provide an accurate approximation as the length of the Periodic A<italic><sub>n</sub></italic>/T<italic><sub>n</sub></italic> Cluster (PATC) increases. At extremes of possibility, it could be imagined that such DNAs would form complex coils and knots in solution or, alternatively, might behave as relatively straight segments perhaps with unusual microscopic flexibility properties.</p><p>In order to understand the genesis and role for extended PATCs in eukaryotic genomes, we first need to define their physical properties. Genomes for <italic>C.elegans</italic> and other nematode species have been shown to exhibit an unusual high frequency of extended PATCs (<xref ref-type="bibr" rid="b20">20</xref>,<xref ref-type="bibr" rid="b23">23</xref>,<xref ref-type="bibr" rid="b24">24</xref>), particularly in regions of the genome that are transcribed in germline tissues. The <italic>C.elegans</italic> genome has the additional advantage that it is among the most extensively characterized both in structure and function, since the sequence is known and individual genetic studies, genome expression and phenotype analysis have been done.</p><p>In this work, we have used gel electrophoresis and AFM to carry out an extensive analysis of several extended PATCs from the genome of the nematode <italic>C.elegans</italic>. We find that the hyperperiodic DNA shows a significant and readily discernible degree of bending compared with non-periodic DNAs, adopting correlated angles over long distances, and that increased flexibility alone cannot explain the mechanical properties of this DNA.</p></sec><sec sec-type="materials|methods"><title>MATERIALS AND METHODS</title><sec><title>AFM constructs</title><p>For AFM imaging, we fabricated a pseudo-palindrome construct of 2885 bp using a method similar to that described in Ref. (<xref ref-type="bibr" rid="b14">14</xref>). This construct contained two segments of 923 bp of hyperperiodic &#x003a9;4 fragments from <italic>C.elegans</italic> in an inverse orientation each located at the ends of the construct (for description of the fragments of <italic>C.elegans</italic> used, see Results). We refer to this construct as the periodic palindrome construct (PPC). Our AFM analysis focused on the first 953 bp at both ends, which included the hyperperiodic sequences. We analyzed a fragment of DNA slightly longer than the pure hyperperiodic fragment because the very ends were not considered in the analysis. The extra non-periodic 979 bp, located in the middle part of the construct, was chosen to be non-palindromic, although it was also part of the <italic>C.elegans</italic> genome. In doing so, we avoided the synthesis of a real palindromic construct, which could generate two independent double-stranded DNA molecules of length equal to half of the original one. The scheme of the PPC, including the sequence of one of the hyperperiodic DNA fragments (capital letters), is shown in <xref ref-type="fig" rid="fig1">Figure 1a</xref>. PPC fabrication and long-term stability was verified using gel electrophoresis.</p><p>We also analyzed two short DNAs, which contained fragments &#x003a9;6 and &#x003a9;7 (for description of these fragments see Results). Fragments &#x003a9;7 and &#x003a9;6 were amplified by cloning into the pCR4Blunt-TOPO vector (Promega), yielding two short DNAs (SQ100 and SQ322) that were analyzed with gel electrophoresis and AFM. SQ100 (185 bp) contained the non-periodic fragment &#x003a9;7 (<xref ref-type="fig" rid="fig1">Figure 1b</xref>) and SQ322 (182 bp) contained the hyperperiodic fragment &#x003a9;6 (<xref ref-type="fig" rid="fig1">Figure 1c</xref>). Capital letters in <xref ref-type="fig" rid="fig1">Figure 1b and c</xref> denote the &#x003a9;7 and &#x003a9;6 sequence components. Sequences are specified excluding the single-stranded overhangs.</p><p>All sequences in <xref ref-type="fig" rid="fig1">Figure 1</xref> are displayed in a 10 bp arrangement to facilitate visualization of same-residue coincidence every helical turn. To assist this, two or more consecutive Ts are colored in red and two or more consecutive As in green. The sequences in <xref ref-type="fig" rid="fig1">Figure 1a and c</xref>, corresponding to the hyperperiodic fragment in the PPC and the SQ322 construct fragment, respectively, demonstrate a propensity for an AA/TT to be present every 10 bp. In contrast, the degree of periodicity in SQ100 (<xref ref-type="fig" rid="fig1">Figure 1b</xref>) is negligible.</p></sec><sec><title>Gel electrophoresis</title><p>Both 1D and 2D electrophoresis techniques were used to identify anomalous migration behavior of fragments of DNA from <italic>C.elegans</italic>. 2D gels were cast in 4% NuSieve agarose gels using 1&#x000d7; TAE Buffer in the absence of ethidium bromide. For DNA samples run with an internal standard, it was important to avoid overloading of DNA and consequent smearing of the first dimension during electrophoresis; this was achieved by constructing single-tooth combs that were relatively wide (2.5&#x02013;5 mm) rather than the narrower combs used in other 2D gel systems. The first &#x02018;vertical&#x02019; direction in all 2D gels was run at 53 V for &#x0223c;3 h in a cold room at &#x0223c;4&#x000b0;C. At this voltage, heating by the electrophoretic current was observed to be minimal. Each gel was then soaked for &#x0223c;2 h in 1&#x000d7; TAE Buffer with ethidium bromide (0.3 &#x000b5;g/ml) at room temperature. Next, the second dimension (horizontal direction in the 2D gels) was run with the gel box at room temperature at 80 V for 2 h. Temperatures are somewhat raised by the higher voltage of the second dimension electrophoresis, and were estimated to be &#x0223c;37&#x000b0;C.</p><p>Similar analysis using 1D electrophoresis was carried out by running two separate gels on each sample. One gel was cast with 0.3 &#x000b5;g/ml ethidium bromide present and run at 60 V for 6 h at room temperature, while the other was cast with no intercalator present and run at 30 V for 15 h in a cold room at &#x0223c;4&#x000b0;C and then stained with ethidium bromide (2.0 &#x000b5;g/ml) in 1&#x000d7; TAE for 1 h. Although the two gel systems are comparable in resolution, the paired 1D gels offer the ability to analyze greater numbers of samples, while the 2D gels allow more sensitive detection of modest electrophoretic anomalies.</p></sec><sec><title>Sample preparation and AFM imaging</title><p>DNA samples for AFM imaging were prepared by depositing 5 ng of DNA, diluted in 6 &#x000b5;l of 10 mM Tris&#x02013;HCl (pH 8.0) and then supplemented with 15 mM MgCl<sub>2</sub>, onto freshly cleaved mica. After 60 s, the mica was washed with MilliQ-filtered water and blown dry in a gentle stream of nitrogen gas. Samples were imaged at 1.95 nm per pixel resolution in air at room temperature and humidity using tapping mode. We used an AFM from Nanotec (Nanotec Electr&#x000f3;nica, Madrid, Spain) and also a Nanoscope IIIa (Digital Instruments, Santa Barbara, CA). We utilized tapping mode SuperSharpSilicon tips, type SSS-NCH-8 (Nanosensors, Neuch&#x000e2;tel, Switzerland). Image processing, including plane subtraction and flattening, was done using WSxM freeware (<ext-link ext-link-type="uri" xlink:href="www.nanotec.es"/>). Traces were obtained using the tracing routine described in (P.A. Wiggins, T. van der Heijden, F. Moreno-Herrero, A. Spakowitz, R. Phillips, J. Widom, C. Dekker and P.C. Nelson, manuscript submitted).</p></sec><sec><title>Calculation of periodicity value and predictions of curvature</title><p>Calculations of predicted radius of curvature for DNA segments and of overall periodicity used several available algorithms (<xref ref-type="bibr" rid="b5">5</xref>,<xref ref-type="bibr" rid="b7">7</xref>,<xref ref-type="bibr" rid="b19">19</xref>) (see also Supplementary Data). Source code for local analysis and display of these parameters is available on request.</p></sec></sec><sec><title>RESULTS</title><sec><title>Periodic DNA segments from <italic>C.elegans</italic></title><p>We chose a number of hyperperiodic DNA fragments from <italic>C.elegans</italic> for analysis based on the relevance of the proteins that they encode and on the periodicity value provided by the PATC algorithm described in Ref. (<xref ref-type="bibr" rid="b19">19</xref>). This algorithm is by no means unique and we note that many different algorithms identify a strongly periodic character to <italic>C.elegans</italic> DNA (data not shown) and could have been used to select specific fragments. For reference, the PATC algorithm returns a numerical score with greater numbers indicating more highly periodic character. A score of &#x0003e;95 is reached at &#x0003c;10<sup>&#x02212;6</sup> frequency in a random sequence with base composition equivalent to <italic>C.elegans</italic>. There is no segment in the entire unique genomes of <italic>Drosophila melanogaster</italic> (euchromatin) or <italic>Schizosaccharomyces pombe</italic> (two genomes with comparable AT content to <italic>C.elegans</italic>) with a periodicity score &#x0003e;200. Of the natural bent DNAs that have been analyzed previously on a biophysical level, only the mitochondrial DNA of <italic>Crithidia</italic> shows a substantial signal in the PATC assay (a score of 235 maintained over a distance of 175 bp). Supplementary Figure S2 shows the periodicity values as a function of sequence for the DNAs used in this work.</p><p>Fragment &#x003a9;4 is a 932 bp segment comprising the fourth intron of <italic>C.elegans</italic> gene F54C4.1 on chromosome III. The F54C4.1 gene is essential, encodes mitochondrial ribosomal protein L40, and is likely expressed in both somatic and germline tissue of the animal. The intron 4 fragment &#x003a9;4 was chosen from among the many hyperperiodic segments in the <italic>C.elegans</italic> genome based on its localization in a well-defined essential gene. Informatic analysis of fragment &#x003a9;4 using the PATC algorithm (<xref ref-type="bibr" rid="b19">19</xref>) or standard bending algorithms (<xref ref-type="bibr" rid="b5">5</xref>,<xref ref-type="bibr" rid="b7">7</xref>) indicated a degree of periodicity that is characteristic of a subset of germline-expressed <italic>C.elegans</italic> genes but which has not been seen outside of nematode DNA. Using the PATC algorithm, an average score of 433 for the &#x003a9;4 fragment was observed. This magnitude of score has never been observed in random DNA or in the genomes of non-nematode species. Of sequences in the &#x003a9;4 fragment 89.7% have PATC scores above 95 (the cutoff for 10<sup>&#x02212;6</sup> occurrence in random DNA).</p><p>Several other fragments have also been analyzed in this work. Fragment &#x003a9;6 is a strongly periodic 147 bp segment comprising the third intron of <italic>C.elegans</italic> gene ZK792.7 on chromosome IV. The gene is predicted to encode a protein with metallophosphoesterase activity. The periodicity value of &#x003a9;6 is 122 with 76.9% of the sequences with PATC scores above 95. Fragment &#x003a9;6 was originally isolated among a set of random clones of Microccocal nuclease treated <italic>C.elegans</italic> lysates, consistent with the possibility that this fragment represents a single nucleosome core (<xref ref-type="bibr" rid="b25">25</xref>). We stress, however, that no further <italic>in vivo</italic> analysis has been carried out on this segment of the chromosome.</p><p>As controls, we have examined the behavior of several DNAs which had no unusual periodic character. These include lambda DNA, pUC plasmid backbones and &#x003a9;7, a distinct <italic>C.elegans</italic> genome fragment. Fragment &#x003a9;7 was 150 bp and was also derived from Micrococcal nuclease digestion of <italic>C.elegans</italic> lysates, but did not show distinct periodicity. &#x003a9;7 is an exonic fragment from gene Y38F1A.3. <xref ref-type="table" rid="tbl1">Table 1</xref> summarizes the DNAs employed in this study including the origin of the DNA, the periodicity value and the derived AFM construct.</p><p>From the periodicities above one would expect unusually bent DNA structures for the fragments &#x003a9;4 and &#x003a9;6, but not for &#x003a9;7. We tested this prediction by gel electrophoresis, which exactly confirmed the expectations (<xref ref-type="fig" rid="fig2">Figure 2</xref>). Each of the hyperperiodic fragments (&#x003a9;4 and &#x003a9;6) was retarded compared to size markers, with retardation relieved at higher temperatures and with the addition of an intercalating agent (ethidium bromide) during electrophoresis. In contrast, the non-periodic fragment &#x003a9;7 did not show any retardation. Interestingly, our observations were made using agarose gel electrophoresis, which has not been seen to distinguish DNA structures in previous studies (<xref ref-type="bibr" rid="b3">3</xref>). Similar migration anomalies were observed in other two different hyperperiodic fragments from <italic>C.elegans</italic> (data not shown).</p></sec><sec><title>Single-molecule AFM imaging</title><p>After the electrophoresis experiments indicated the presence of unusual DNA structures, we analyzed our samples with quantitative AFM analysis. The different DNA samples were imaged at high resolution and traced using a custom-written tracing routine. Typically, measured DNA profiles displayed a width of 4&#x02013;7 nm at half height. We imaged the PPC (2885 bp), the control DNA (2743 bp) and the short <italic>C.elegans</italic> fragments (182 and 185 bp). <xref ref-type="fig" rid="fig3">Figure 3a</xref> shows an image of the PPC, where two DNA molecules are visible. The average measured contour length <italic>L</italic><sub>0</sub> was 1032 &#x000b1; 27 nm, yielding a rise per bp of 0.36 &#x000b1; 0.01 nm. To illustrate the analysis, we plotted in <xref ref-type="fig" rid="fig3">Figure 3b</xref> an example of a single trace. The tracing routine accurately generates XY pairs separated by &#x02113; = 2.5 nm that follow the contour of the molecule. We found this length to be the lowest meaningful length that faithfully followed the contour given the pixel size of the images, 1.95 nm, and the high resolution obtained by using ultra-sharp tips. A magnification of the trace illustrates the three quantities deduced that are important for the analysis that follows (<xref ref-type="fig" rid="fig3">Figure 3b</xref>, inset). These are the contour length spacing <italic>L</italic> (in the example, <italic>L</italic> equals 3 &#x000d7; 2.5 = 7.5 nm), and the end-to-end distance <italic>R</italic>. Clearly, a given contour length spacing <italic>L</italic> can be associated with several values of <italic>R</italic> (three of them are plotted in the inset). The inset also shows the definition of curvature, <italic>C</italic>. The curvature of a 10 nm long segment is defined as the inverse of the radius <italic>r</italic> of the circle that encompasses three points separated by 5 nm (showing in red in the inset in <xref ref-type="fig" rid="fig3">Figure 3b</xref>).</p><p>High-resolution AFM images of the short DNA samples, SQ100 and SQ322, are shown in <xref ref-type="fig" rid="fig3">Figure 3c and d</xref>. <xref ref-type="fig" rid="fig3">Figure 3c</xref> corresponds to the non-periodic sequence and <xref ref-type="fig" rid="fig3">Figure 3d</xref> to the hyperperiodic one. By eye, one can already notice that <xref ref-type="fig" rid="fig3">Figure 3d</xref> displays more bent molecules. This was confirmed in the subsequent numerical and analytical analysis (see below). Insets represent histograms of contour lengths for each sample. They show that our tracing routine precisely deduces the contour length with a standard deviation of 3 nm. For samples SQ100 and SQ322, 692 and 653 molecules were analyzed, respectively.</p></sec><sec><title>Single molecule AFM-based quantitative analysis of hyperperiodic <italic>C.elegans</italic> DNA</title><p>A quantitative analysis of the traces of hyperperiodic DNA showed systematic deviations from the worm-like-chain (WLC) model consistent with extra bending induced by sequence periodicity. As described by Rivetti <italic>et al</italic>. (<xref ref-type="bibr" rid="b26">26</xref>) DNA equilibrates on a mica surface when adsorption is promoted by divalent ions such as Mg<sup>2+</sup>. The behavior of DNA molecules under these conditions can be predicted by the WLC model (<xref ref-type="bibr" rid="b26">26</xref>). The WLC describes a mean trajectory of the polymer on the mica with the following equation:
<disp-formula id="e1"><label>1</label><mml:math id="M1"><mml:mrow><mml:mrow><mml:mo>&#x03008;</mml:mo><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mo>&#x03009;</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>4</mml:mn><mml:mi>P</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mi>L</mml:mi><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:mi>P</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>&#x02212;</mml:mo><mml:mrow><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mrow><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mn>2</mml:mn><mml:mi>P</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:mrow></mml:mrow></mml:msup><mml:mo>&#x02212;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
where <italic>P</italic> is the persistence length of the polymer, <italic>L</italic> is the contour length spacing and <inline-formula><mml:math id="M2"><mml:mrow><mml:mrow><mml:mo>&#x03008;</mml:mo><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mo>&#x03009;</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> is the mean-squared separation between points along the chains located at distances <italic>s</italic> and <italic>s+L</italic> from the end of the molecule. The average is computed over <italic>s</italic> and over all observed contours. The averaging is done up to the upper limit in the contour length of the molecule <italic>L</italic><sub>0</sub> (for the PPC, <italic>L</italic><sub>0</sub> equals 325 nm) minus the contour length spacing <italic>L</italic> analyzed, thus 0 &#x0003c; <italic>s</italic> &#x0003c; <italic>L</italic><sub>0</sub> &#x02212; <italic>L</italic>.</p><p>Using this analysis, the control DNA and the non-periodic short fragment SQ100 displayed exactly the behavior predicted by the WLC model; the hyperperiodic DNA, however, did not. <xref ref-type="fig" rid="fig4">Figure 4a</xref> shows <inline-formula><mml:math id="M3"><mml:mrow><mml:mrow><mml:mo>&#x03008;</mml:mo><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mo>&#x03009;</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> versus <italic>L</italic> for the Control and PPC samples. We calculated <inline-formula><mml:math id="M4"><mml:mrow><mml:mrow><mml:mo>&#x03008;</mml:mo><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mo>&#x03009;</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> up to a contour length spacing of <italic>L</italic> = 200 nm. This is because the length of the segment studied in the PPC was 325 nm, necessarily restricting one to values of <italic>L</italic> &#x0003c; <italic>L</italic><sub>0</sub> to obtain a statistically reliable calculation for all <italic>L</italic>. For consistency, we also used an upper limit for <italic>L</italic> of 200 nm for the control DNA. A fit using <xref ref-type="disp-formula" rid="e1">Equation 1</xref> was done for the control DNA (<xref ref-type="fig" rid="fig4">Figure 4a</xref>, black line). The fit yielded a value for the persistence length of 51 nm, in agreement with published results (<xref ref-type="bibr" rid="b26">26</xref>,<xref ref-type="bibr" rid="b27">27</xref>). These data were consistent with the hypothesis that our sample preparation yields contours that reflect equilibrium 2D chain conformations. We repeated the same type of analysis on the PPC, tracing 325 nm (953 bp) from both ends. The use of a palindrome construct is very convenient because it avoids ambiguity in the sequence analyzed: the traced fragments always included the hyperperiodic Fragment &#x003a9;4 from <italic>C.elegans</italic> independently of which end was selected. In this way, 190 fragments coming from three sets of experiments done on different days, mica surfaces and preparations were traced and analyzed. Remarkably, the same analysis done on the PPC showed significantly smaller <inline-formula><mml:math id="M5"><mml:mrow><mml:mrow><mml:mo>&#x03008;</mml:mo><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mo>&#x03009;</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> than predicted by the WLC model (<xref ref-type="fig" rid="fig4">Figure 4a</xref>), which is indicative of extra bending. To guide the eye we included plots of <xref ref-type="disp-formula" rid="e1">Equation 1</xref> using <italic>P</italic> = 20, 30, 40, 50 and 60 nm (solid lines in <xref ref-type="fig" rid="fig4">Figure 4a</xref>). The PPC data curves mostly overlaid the curve of <italic>P</italic> = 30 nm. It should be noted, though, that the WLC model does not necessarily hold for intrinsically curved DNA.</p><p>We conducted a similar analysis using samples SQ100 and SQ322. In this case due to the shorter length of the molecules we plotted <inline-formula><mml:math id="M6"><mml:mrow><mml:mrow><mml:mo>&#x03008;</mml:mo><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mo>&#x03009;</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> up to a contour length spacing of <italic>L</italic> = 47.5 nm (<xref ref-type="fig" rid="fig4">Figure 4b</xref>). This required tracing many more molecules to have comparable statistics with those for the PPC sample. We traced 692 and 653 molecules for SQ100 and SQ322, respectively. The analysis showed that SQ322 contours had smaller mean squared end-to-end distances than those of non-periodic SQ100 (<xref ref-type="fig" rid="fig4">Figure 4b</xref>). To stress the significance of the result, we included in <xref ref-type="fig" rid="fig4">Figure 4b</xref> plots of <xref ref-type="disp-formula" rid="e1">Equation 1</xref> with <italic>P</italic> = 20, 30, 40, 50 and 60 nm to guide the eye. Whereas SQ100 data lay as expected between the curves of <italic>P</italic> = 50 and 60 nm, SQ322 data overlaid the curve of <italic>P</italic> = 25 nm.</p><p>Summarizing the results from these different DNA samples (<xref ref-type="fig" rid="fig4">Figure 4a and b</xref>), we find consistently that hyperperiodic DNA results in a shorter <inline-formula><mml:math id="M7"><mml:mrow><mml:mrow><mml:mo>&#x03008;</mml:mo><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mo>&#x03009;</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> than predicted by the WLC at the two length scales studied: hundreds of bp and tens of bp. This can be most easily explained by increased bending of the DNA. Furthermore, the plot for <inline-formula><mml:math id="M8"><mml:mrow><mml:mrow><mml:mo>&#x03008;</mml:mo><mml:mrow><mml:msubsup><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow><mml:mo>&#x03009;</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> deviates from the WLC at lower persistence length. This suggests that the increase in bending is not due to an increase in random bending/flexibility, which would exhibit WLC behavior, but rather due to non-random bending.</p><p>The presence of increased bending in hyperperiodic <italic>C.elegans</italic> DNA compared to non-periodic DNA was further supported by examining distributions of absolute curvatures &#x02223;<italic>C</italic>&#x02223;. We calculated curvatures from traces at a fixed spacing <italic>L</italic> (see inset in <xref ref-type="fig" rid="fig3">Figure 3b</xref> as an example, for <italic>L</italic> = 5 nm). To do this, we fitted a circle to three points separated by <italic>L</italic> along the contour of the DNA. The curvature (<italic>C</italic>) at this spacing is defined as the inverse of the radius of that circle and its value was assigned to the middle point of the three contour points. Hence to calculate <italic>C<sub>i</sub></italic> at a given position <italic>i</italic>, for a given <italic>L</italic>, we used DNA stretches of <italic>2</italic> &#x000d7; <italic>L</italic> nm. Then we moved to the next point <italic>i</italic> + &#x02113; along the contour, which is &#x02113; = 2.5 nm (7 bp) apart from the previous, and we obtained a new value of <italic>C<sub>i+</sub></italic><sub>&#x02113;</sub>. We arbitrarily assigned positive curvature if by progressing along the DNA the contour turned to the right, and we assigned negative curvature if it turned to the left. As expected, we found Gaussian distributions of <italic>C</italic> centered at zero for all type of molecules studied, proving no preferential curvature orientations due to adsorption (data not shown).</p><p>Absolute curvatures &#x02223;<italic>C</italic>&#x02223; distributions were calculated at different values of <italic>L</italic>. <xref ref-type="fig" rid="fig4">Figure 4c</xref> shows absolute curvature distributions for the PPC and the control non-periodic DNA for <italic>L</italic> = 5 nm. We analyzed 100 control DNA molecules, which yielded &#x0223c;12&#x02009;000 curvature values, and 190 PPC molecules, which yielded &#x0223c;24&#x02009;000 curvature values. <xref ref-type="fig" rid="fig4">Figure 4d</xref> shows similar absolute curvature distributions for the short DNA fragments studied, SQ100 (non-periodic) and SQ322 (highly-periodic). We analyzed 613 SQ100 and 567 SQ322 molecules, each yielding &#x0223c;12&#x02009;000 curvature values. Distributions were normalized to have an area of one, and error bars denote normalized square root of number of points in the bin. All absolute curvature distributions were fit with a Gaussian function centered at curvature zero. From these fits we extracted the standard deviation &#x003c3;<italic><sub>C</sub></italic>. The values &#x003c3;<italic><sub>C</sub></italic> of the Control, PPC, SQ100 and SQ322 samples from curvature histograms calculated using a spacing of <italic>L</italic> = 5 and 25 nm are quoted in <xref ref-type="table" rid="tbl2">Table 2</xref> (the latter only for long DNAs). Systematically, we found that DNA with hyperperiodic sequences had larger &#x003c3;<italic><sub>C</sub></italic> than non-periodic DNA. This was found to be independent of the spacing <italic>L</italic> used to calculate curvatures up to 50 nm. The inset in <xref ref-type="fig" rid="fig4">Figure 4c</xref> proves this by showing &#x02223;<italic>C</italic>&#x02223; distributions for PPC and control samples calculated at a larger spacing of <italic>L</italic> = 25 nm. Spacings larger than 50 nm, the persistence length of the DNA, did not generate statistically significant curvature distributions.</p></sec><sec><title>Intrinsic curvature in hyperperiodic <italic>C.elegans</italic> DNA</title><p>In principle, a broad &#x02223;<italic>C</italic>&#x02223; distribution could result from either increased flexibility or intrinsic curvature (Supplementary Data). In order to determine which of the two phenomena dominates in our hyperperiodic DNA, we performed two calculations: (i) we calculated the distance over which the handedness of the curvature persists, and (ii) we calculated an average curvature per single imaged AFM molecule and considered its absolute value &#x02223;&#x0003c;<italic>C</italic>&#x0003e;<sub>molecule</sub>&#x02223;.</p><p>One way to measure the persistence of handedness is to plot a correlation function <italic>CC<sub>L,D</sub></italic> of the form:
<disp-formula id="e2"><label>2</label><mml:math id="M9"><mml:mrow><mml:mi>C</mml:mi><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>,</mml:mo><mml:mi>D</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msubsup><mml:mo>&#x02211;</mml:mo><mml:mi>s</mml:mi><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>&#x02212;</mml:mo><mml:mi>D</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>D</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow><mml:mrow><mml:msubsup><mml:mo>&#x02211;</mml:mo><mml:mi>s</mml:mi><mml:mrow><mml:msub><mml:mi>L</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>&#x02212;</mml:mo><mml:mi>D</mml:mi></mml:mrow></mml:msubsup><mml:mrow><mml:mrow><mml:mo>|</mml:mo><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>|</mml:mo></mml:mrow><mml:mrow><mml:mo>|</mml:mo><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mrow><mml:mi>L</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>+</mml:mo><mml:mi>D</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>|</mml:mo></mml:mrow></mml:mrow></mml:mrow></mml:mfrac><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
where <italic>C<sub>L,s</sub></italic> is the curvature calculated at spacing <italic>L</italic> at position <italic>s</italic> and <italic>L</italic><sub>0</sub> is the total contour length of the molecule. <italic>C<sub>L,s</sub></italic> can therefore be positive or negative as defined above. The function <italic>CC<sub>L,D</sub></italic> involves addition over all the molecules considered in this study, 100 for the control DNA and 190 for the PPC. To examine the persistence of handedness, renormalization was included to compare the two datasets. For a given spacing <italic>L</italic>, <italic>CC<sub>L,D</sub></italic> is a function that depends on the correlation spacing <italic>D</italic>. <xref ref-type="fig" rid="fig5">Figure 5a</xref> shows a plot of <italic>CC<sub>L,D</sub></italic> for <italic>L</italic> = 5, 15 and 25 nm for the PPC (dashed lines) and the control DNA (solid lines). Systematically, a significantly larger correlation distance was found for the hyperperiodic DNA than for the non-periodic DNA. As expected, at large correlation spacing <italic>D</italic>, the correlation function approached zero for non-periodic DNA, showing that the &#x02018;memory&#x02019; of handedness is lost very quickly. Interestingly, for hyperperiodic DNA the correlation function approaches a non-zero value, for large <italic>D</italic>. Qualitatively similar correlation functions were predicted from theoretical calculations for chains with long-range correlations (<xref ref-type="bibr" rid="b28">28</xref>). Therefore, these data are strongly indicative of the presence of an intrinsic bend rather than increased flexibility.</p><p>A complementary approach to the correlation distance consisted of calculating an average curvature per molecule and taking its absolute value &#x02223;&#x0003c;<italic>C</italic>&#x0003e;<sub>molecule</sub>&#x02223;. The distributions of absolute curvature per molecule are shown in <xref ref-type="fig" rid="fig5">Figure 5b and c</xref>. The statistics in these experiments is poorer (resulting in larger error bars) than in the &#x02223;<italic>C</italic>&#x02223; distributions (shown in <xref ref-type="fig" rid="fig4">Figure 4c and d</xref>), as per definition the total number of values for the histogram equals the number of molecules imaged. <xref ref-type="fig" rid="fig5">Figure 5b</xref> shows &#x02223;&#x0003c;<italic>C</italic>&#x0003e;<sub>molecule</sub>&#x02223; distributions for the PPC and the non-periodic control DNA, and <xref ref-type="fig" rid="fig5">Figure 5c</xref> shows &#x02223;&#x0003c;<italic>C</italic>&#x0003e;<sub>molecule</sub>&#x02223; for the SQ100 and SQ322 fragments. Distributions were calculated using a spacing of <italic>L</italic> = 5 nm to calculate the curvature values, except the inset in <xref ref-type="fig" rid="fig5">Figure 5b</xref>, which used <italic>L</italic> = 25 nm. The distributions were fit with Gaussian functions. As expected, the distributions for both non-periodic DNA molecules had a maximum at curvature zero. However, the distributions for the two hyperperiodic DNA samples displayed an offset, &#x02223;&#x0003c;<italic>C</italic>&#x0003e;<sub>0,molecule</sub> <italic><sub>L=5</sub></italic> <sub>nm</sub>&#x02223; of &#x0223c;0.01 nm<sup>&#x02212;1</sup>. This implies that a given molecule has the tendency to bend either more towards the right than towards the left or vice versa. This observation can only be explained by the presence of an intrinsic curvature, as the distributions for the &#x02018;random&#x02019; control DNAs had a maximum at around zero. We note that the effect is more pronounced for the long DNA samples, as here more curvature values contribute to the mean curvature, which can directly be seen from the width of the distributions &#x003c3;<sub>&#x0003c;</sub><italic><sub>C</sub></italic><sub>&#x0003e;molecule</sub> (<xref ref-type="table" rid="tbl2">Table 2</xref>). Similar conclusions can be drawn from the analysis done at <italic>L</italic> = 25 nm (inset in <xref ref-type="fig" rid="fig5">Figure 5b</xref>), showing that the detection of an intrinsic bend is independent of the spacing <italic>L</italic> used.</p></sec></sec><sec><title>DISCUSSION</title><p>We report anomalous gel migration of several fragments of the nematode <italic>C.elegans</italic> genome with unusual periodicity of AA/TT dinucleotides. AFM allowed detailed characterization of the bending of these DNA structures showing that PATCs generate intrinsic bending. Despite the presence of an intrinsic bend, the AFM images also showed that hyperperiodic DNA can still adopt a wide range of bending configurations (<xref ref-type="fig" rid="fig3">Figure 3d</xref>), suggesting that flexibility is retained.</p><p>In order to be able to extract quantitative information related to bending of DNA, we first verified that our DNA contours reflect 2D equilibrium states and that the average behavior can be accurately described by the WLC model. For this, we used as a control the plasmid pGem-3Z, which had no unusual periodic character. We observed that the average trajectory of adsorbed control DNA molecules fulfilled several criteria: (i) the distribution of curvatures was Gaussian (<xref ref-type="fig" rid="fig4">Figure 4c</xref>, black squares), and (ii) the mean square end-to-end distance followed <xref ref-type="disp-formula" rid="e1">Equation 1</xref> with the accepted value for the persistence length (<xref ref-type="fig" rid="fig4">Figure 4a</xref>, black squares). Thus, we conclude that the WLC describes accurately the average trajectory of our molecules and that these trajectories reflect equilibrium states of the molecules on the mica surface. The fact that the control DNA molecules are equilibrated on the surface is crucial for allowing us to extend the same analysis to DNAs with hyperperiodic sequences. Plots of mean squared end-to-end distances of hyperperiodic DNA samples, PPC (<xref ref-type="fig" rid="fig4">Figure 4a</xref>, colored symbols) and SQ322 (<xref ref-type="fig" rid="fig4">Figure 4b</xref>, red circles), showed that AA/TT periodicity induces extra bending in the DNA structure. This was demonstrated by significant deviations of the hyperperiodic DNA from the predictions of the WLC model using a persistence length of 50 nm.</p><p>Subsequent analysis to determine the source of this deviation took into account distributions of absolute values of curvatures &#x02223;<italic>C</italic>&#x02223;. As a check for the validity of this analysis, we found that curvature distributions had a maximum at zero, showing equal probability for the observation of positive and negative curvatures. The value &#x003c3;<sub>C</sub> (<xref ref-type="table" rid="tbl2">Table 2</xref>) obtained from Gaussian fits to the &#x02223;<italic>C</italic>&#x02223; distributions indicated how often large curvatures were present in the trace or the degree of bendability. <xref ref-type="fig" rid="fig4">Figure 4c</xref> and d showed that hyperperiodic DNA (PPC and SQ322) accommodates high curvatures more readily than non-periodic DNA (Control and SQ100). However, it is not possible from the value of &#x003c3;<sub>C</sub> alone to attribute the increased presence of high curvatures to increased flexibility or to the presence of intrinsic bends: indeed either intrinsic curvature or unusual flexibility will widen the <italic>C</italic> distribution (Supplementary Data). Nevertheless, it is accurate to say that if the DNA does not have an intrinsic bend, then the value of &#x003c3;<sub>C</sub> directly reflects the degree of flexibility.</p><p>One characteristic feature of the <italic>C.elegans</italic> hyperperiodic sequences analyzed in this work is the distance over which substantial periodicity is maintained. In the AFM images, this corresponds to what appears to be a considerable distance over which the handedness of the curvature persists (<xref ref-type="fig" rid="fig3">Figure 3a</xref>). We quantified the distance <italic>D</italic> over which handedness was correlated using <xref ref-type="disp-formula" rid="e2">Equation 2</xref>. These correlation plots were computed for the PPC and the Control DNA (<xref ref-type="fig" rid="fig5">Figure 5a</xref>). It was observed that at a given distance <italic>D</italic> the probability of finding a curvature with the same sign was significantly higher in the hyperperiodic DNA than in non-periodic DNA. The curvature correlation increased with larger values of <italic>L</italic>. This indicates that AA/TT periodicity induces a long-scale effect in the structure of DNA and provides a strong argument for bending as opposed to flexibility, since flexibility would definitely not increase the correlation distance. In principle, the process of adsorption itself could cause switching of handedness as the DNA might flip over at certain points, leading to an underestimation of the correlation distance <italic>D</italic>, compared to the 3D structure in solution.</p><p>A complementary approach to analyze the presence or absence of long-scale correlated intrinsic bends consists of calculating the mean curvatures per molecule (<xref ref-type="fig" rid="fig5">Figure 5b and c</xref>). This analysis showed a maximum at zero for the non-periodic DNA samples (control and SQ100), in agreement with purely random bending. A clear peak, however, was found for both of the hyperperiodic DNA samples (PPC and SQ322), suggesting the presence of a large-scale intrinsic bend with an average curvature of &#x0223c;0.01 nm<sup>&#x02212;1</sup> (at <italic>L</italic> = 5 nm). It should be noted that we underestimate the bending of the DNA, because it is not clear how the bent regions absorb to the mica, resulting in a mixture of regions with positive and negative curvature. The absence of intrinsic bending for the non-periodic SQ100 DNA suggests that its widened curvature distribution (<xref ref-type="fig" rid="fig4">Figure 4d</xref>) could be due to an increase in flexibility compared to the non-periodic Control (<xref ref-type="fig" rid="fig4">Figure 4c</xref>). For hyperperiodic DNA, however, it was not possible to address flexibility, as flexibility is convoluted by intrinsic bending in the curvature distributions. However, the variety of shapes displayed in the AFM images (<xref ref-type="fig" rid="fig3">Figure 3d</xref>) indicated that hyperperiodic DNA seems to be flexible enough to adopt straight and strongly bent conformations.</p><p>DNA sequences from numerous organisms have been shown to contain regions with periodic character corresponding to the helical repeat of 10&#x02013;10.5 bp. Despite a number of cases in which structural roles have been assigned to relatively short periodic regions (<xref ref-type="bibr" rid="b17">17</xref>), the functions of extended periodic regions in DNA are unknown. Of the numerous &#x02018;model&#x02019; genomes that have been subject to intensive analysis, <italic>C.elegans</italic> has (at least by some measures) the strongest periodicity signal (<xref ref-type="bibr" rid="b19">19</xref>,<xref ref-type="bibr" rid="b24">24</xref>). Given that 1&#x02013;2% of the base pairs in the <italic>C.elegans</italic> genome appear constrained by their contribution to periodic character (<xref ref-type="bibr" rid="b19">19</xref>), it seems likely that these sequences have substantial roles in the evolution and/or function of the genome. One role that has been considered by virtually every investigator who has studied these sequences has been to constrain or otherwise influence nucleosome positioning. Despite these proposals, it remains a possibility that other extended surfaces in the nucleus are implicated in the observed long-range sequence relationships.</p><p>In summary, we analyzed molecules with a greater degree and extent of periodicity than any previously analyzed DNA. We used gel electrophoresis and AFM to investigate the effect of long PATCs in the structure of the genome of the nematode <italic>C.elegans</italic>. Gel electrophoresis showed that hyperperiodic DNA migrates more slowly than non-periodic DNA. We quantified trajectories of equilibrated molecules on a mica surface obtained with high-resolution AFM, and proved that the PATCs exhibit large-scale intrinsic bending. This explained the anomalous gel migration. Finally, we studied two nucleosome-length DNA fragments that survived vigorous Micrococcal nuclease digestion <italic>in situ</italic>, and likewise found the presence of an intrinsic bend for the fragment containing hyperperiodic sequences. The fragment without intrinsic bending exhibited increased flexibility. This is consistent with the idea that bendability of the DNA might be a relevant factor in promoting binding to nucleosomes.</p></sec><sec><title>SUPPLEMENTARY DATA</title><p>Supplementary Data are available at NAR Online.</p></sec>
<sec sec-type="supplementary-material">
<title>Supplementary Material</title>
<supplementary-material id="PMC_1" content-type="local-data">
<caption>
<title>[Supplementary Material]</title>
</caption>
<media mimetype="text" mime-subtype="html" xlink:href="nar_34_10_3057__index.html"/>
<media xlink:role="associated-file" mimetype="application" mime-subtype="pdf" xlink:href="nar_34_10_3057__1.pdf"/>
</supplementary-material>
</sec>
</body><back><ack><p>We are grateful to Jamie Fleenor and Chaya Krishna for expert technical assistance, and to Javier Lopez-Molina and Michael McCaffery for early advice on protocols. We thank Fabian Czerwinski for the preliminary AFM experiments. F.M.-H. is supported by a postdoctoral contract from the Dutch research program NanoNed. S.M.J. is supported by a postdoctoral fellowship PF-05-121-01-DDC from the American Cancer Society. The laboratory of N.H.D. acknowledges financial support from FOM (Dutch Foundation for Fundamental Research on Matter) and NWO (Netherlands Organization for Scientific Research). The laboratory of A.F. acknowledges financial support from NIH (GM37706). Funding to pay the Open Access publication charges for this article was provided by FOM, NWO and NIH.</p><p><italic>Conflict of interest statement</italic>. None declared.</p></ack><ref-list><title>REFERENCES</title><ref id="b1"><label>1</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Luger</surname><given-names>K.</given-names></name><name><surname>Mader</surname><given-names>A.W.</given-names></name><name><surname>Richmond</surname><given-names>R.K.</given-names></name><name><surname>Sargent</surname><given-names>D.F.</given-names></name><name><surname>Richmond</surname><given-names>T.J.</given-names></name></person-group><article-title>Crystal structure of the nucleosome core particle at 2.8 A resolution</article-title><source>Nature</source><year>1997</year><volume>389</volume><fpage>251</fpage><lpage>260</lpage><pub-id pub-id-type="pmid">9305837</pub-id></citation></ref><ref id="b2"><label>2</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Trifonov</surname><given-names>E.N.</given-names></name><name><surname>Sussman</surname><given-names>J.L.</given-names></name></person-group><article-title>The pitch of chromatin DNA is reflected in its nucleotide sequence</article-title><source>Proc. Natl Acad. Sci. USA</source><year>1980</year><volume>77</volume><fpage>3816</fpage><lpage>3820</lpage><pub-id pub-id-type="pmid">6933438</pub-id></citation></ref><ref id="b3"><label>3</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marini</surname><given-names>J.C.</given-names></name><name><surname>Levene</surname><given-names>S.D.</given-names></name><name><surname>Crothers</surname><given-names>D.M.</given-names></name><name><surname>Englund</surname><given-names>P.T.</given-names></name></person-group><article-title>Bent helical structure in kinetoplast DNA</article-title><source>Proc. Natl Acad. Sci. USA</source><year>1982</year><volume>79</volume><fpage>7664</fpage><lpage>7668</lpage><pub-id pub-id-type="pmid">16593261</pub-id></citation></ref><ref id="b4"><label>4</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shen</surname><given-names>C.H.</given-names></name><name><surname>Clark</surname><given-names>D.J.</given-names></name></person-group><article-title>DNA sequence plays a major role in determining nucleosome positions in yeast CUP1 chromatin</article-title><source>J. Biol. Chem.</source><year>2001</year><volume>276</volume><fpage>35209</fpage><lpage>35216</lpage><pub-id pub-id-type="pmid">11461917</pub-id></citation></ref><ref id="b5"><label>5</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Koo</surname><given-names>H.S.</given-names></name><name><surname>Wu</surname><given-names>H.M.</given-names></name><name><surname>Crothers</surname><given-names>D.M.</given-names></name></person-group><article-title>DNA bending at adenine. thymine tracts</article-title><source>Nature</source><year>1986</year><volume>320</volume><fpage>501</fpage><lpage>506</lpage><pub-id pub-id-type="pmid">3960133</pub-id></citation></ref><ref id="b6"><label>6</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ulanovsky</surname><given-names>L.E.</given-names></name><name><surname>Trifonov</surname><given-names>E.N.</given-names></name></person-group><article-title>Estimation of wedge components in curved DNA</article-title><source>Nature</source><year>1987</year><volume>326</volume><fpage>720</fpage><lpage>722</lpage><pub-id pub-id-type="pmid">3561514</pub-id></citation></ref><ref id="b7"><label>7</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bolshoy</surname><given-names>A.</given-names></name><name><surname>McNamara</surname><given-names>P.</given-names></name><name><surname>Harrington</surname><given-names>R.E.</given-names></name><name><surname>Trifonov</surname><given-names>E.N.</given-names></name></person-group><article-title>Curved DNA without A-A: experimental estimation of all 16 DNA wedge angles</article-title><source>Proc. Natl Acad. Sci. USA</source><year>1991</year><volume>88</volume><fpage>2312</fpage><lpage>2316</lpage><pub-id pub-id-type="pmid">2006170</pub-id></citation></ref><ref id="b8"><label>8</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hagerman</surname><given-names>P.J.</given-names></name></person-group><article-title>Straightening out the bends in curved DNA</article-title><source>Biochim. Biophys. Acta</source><year>1992</year><volume>1131</volume><fpage>125</fpage><lpage>132</lpage><pub-id pub-id-type="pmid">1610891</pub-id></citation></ref><ref id="b9"><label>9</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dickerson</surname><given-names>R.E.</given-names></name><name><surname>Goodsell</surname><given-names>D.</given-names></name><name><surname>Kopka</surname><given-names>M.L.</given-names></name></person-group><article-title>MPD and DNA bending in crystals and in solution</article-title><source>J. Mol. Biol.</source><year>1996</year><volume>256</volume><fpage>108</fpage><lpage>125</lpage><pub-id pub-id-type="pmid">8609604</pub-id></citation></ref><ref id="b10"><label>10</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Griffith</surname><given-names>J.</given-names></name><name><surname>Bleyman</surname><given-names>M.</given-names></name><name><surname>Rauch</surname><given-names>C.A.</given-names></name><name><surname>Kitchin</surname><given-names>P.A.</given-names></name><name><surname>Englund</surname><given-names>P.T.</given-names></name></person-group><article-title>Visualization of the bent helix in kinetoplast DNA by electron microscopy</article-title><source>Cell</source><year>1986</year><volume>46</volume><fpage>717</fpage><lpage>724</lpage><pub-id pub-id-type="pmid">3742596</pub-id></citation></ref><ref id="b11"><label>11</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Muzard</surname><given-names>G.</given-names></name><name><surname>Theveny</surname><given-names>B.</given-names></name><name><surname>Revet</surname><given-names>B.</given-names></name></person-group><article-title>Electron microscopy mapping of pBR322 DNA curvature. Comparison with theoretical models</article-title><source>EMBO J.</source><year>1990</year><volume>9</volume><fpage>1289</fpage><lpage>1298</lpage><pub-id pub-id-type="pmid">2323339</pub-id></citation></ref><ref id="b12"><label>12</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hansma</surname><given-names>H.G.</given-names></name><name><surname>Browne</surname><given-names>K.A.</given-names></name><name><surname>Bezanilla</surname><given-names>M.</given-names></name><name><surname>Bruice</surname><given-names>T.C.</given-names></name></person-group><article-title>Bending and straightening of DNA induced by the same ligand: characterization with the atomic force microscope</article-title><source>Biochemistry</source><year>1994</year><volume>33</volume><fpage>8436</fpage><lpage>8441</lpage><pub-id pub-id-type="pmid">8031778</pub-id></citation></ref><ref id="b13"><label>13</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rivetti</surname><given-names>C.</given-names></name><name><surname>Walker</surname><given-names>C.</given-names></name><name><surname>Bustamante</surname><given-names>C.</given-names></name></person-group><article-title>Polymer chain statistics and conformational analysis of DNA molecules with bends or sections of different flexibility</article-title><source>J. Mol. Biol.</source><year>1998</year><volume>280</volume><fpage>41</fpage><lpage>59</lpage><pub-id pub-id-type="pmid">9653030</pub-id></citation></ref><ref id="b14"><label>14</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zuccheri</surname><given-names>G.</given-names></name><name><surname>Scipioni</surname><given-names>A.</given-names></name><name><surname>Cavaliere</surname><given-names>V.</given-names></name><name><surname>Gargiulo</surname><given-names>G.</given-names></name><name><surname>De Santis</surname><given-names>P.</given-names></name><name><surname>Samori</surname><given-names>B.</given-names></name></person-group><article-title>Mapping the intrinsic curvature and flexibility along the DNA chain</article-title><source>Proc. Natl Acad. Sci. USA</source><year>2001</year><volume>98</volume><fpage>3074</fpage><lpage>3079</lpage><pub-id pub-id-type="pmid">11248034</pub-id></citation></ref><ref id="b15"><label>15</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nelson</surname><given-names>H.C.</given-names></name><name><surname>Finch</surname><given-names>J.T.</given-names></name><name><surname>Luisi</surname><given-names>B.F.</given-names></name><name><surname>Klug</surname><given-names>A.</given-names></name></person-group><article-title>The structure of an oligo(dA).oligo(dT) tract and its biological implications</article-title><source>Nature</source><year>1987</year><volume>330</volume><fpage>221</fpage><lpage>226</lpage><pub-id pub-id-type="pmid">3670410</pub-id></citation></ref><ref id="b16"><label>16</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lowary</surname><given-names>P.T.</given-names></name><name><surname>Widom</surname><given-names>J.</given-names></name></person-group><article-title>New DNA sequence rules for high affinity binding to histone octamer and sequence-directed nucleosome positioning</article-title><source>J. Mol. Biol.</source><year>1998</year><volume>276</volume><fpage>19</fpage><lpage>42</lpage><pub-id pub-id-type="pmid">9514715</pub-id></citation></ref><ref id="b17"><label>17</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schleif</surname><given-names>R.</given-names></name></person-group><article-title>DNA looping</article-title><source>Annu. Rev. Biochem.</source><year>1992</year><volume>61</volume><fpage>199</fpage><lpage>223</lpage><pub-id pub-id-type="pmid">1497310</pub-id></citation></ref><ref id="b18"><label>18</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>J.P.</given-names></name><name><surname>Widom</surname><given-names>J.</given-names></name></person-group><article-title>Improved alignment of nucleosome DNA sequences using a mixture model</article-title><source>Nucleic Acids Res.</source><year>2005</year><volume>33</volume><fpage>6743</fpage><lpage>6755</lpage><pub-id pub-id-type="pmid">16339114</pub-id></citation></ref><ref id="b19"><label>19</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fire</surname><given-names>A.</given-names></name><name><surname>Alcazar</surname><given-names>R.</given-names></name><name><surname>Tan</surname><given-names>F.</given-names></name></person-group><article-title>Unusual DNA structures associated with germline genetic activity in <italic>Caenorhabditis elegans</italic></article-title><source>Genetics</source><year>2006</year><comment>[Epub 2006 Apr 28]</comment></citation></ref><ref id="b20"><label>20</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>VanWye</surname><given-names>J.D.</given-names></name><name><surname>Bronson</surname><given-names>E.C.</given-names></name><name><surname>Anderson</surname><given-names>J.N.</given-names></name></person-group><article-title>Species-specific patterns of DNA bending and sequence</article-title><source>Nucleic Acids Res.</source><year>1991</year><volume>19</volume><fpage>5253</fpage><lpage>5261</lpage><pub-id pub-id-type="pmid">1923808</pub-id></citation></ref><ref id="b21"><label>21</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kumar</surname><given-names>L.</given-names></name><name><surname>Futschik</surname><given-names>M.</given-names></name><name><surname>Herzel</surname><given-names>H.</given-names></name></person-group><article-title>DNA motifs and sequence periodicities</article-title><source>In Silico Biol.</source><year>2006</year><volume>6</volume><fpage>0008</fpage><pub-id pub-id-type="pmid">16789915</pub-id></citation></ref><ref id="b22"><label>22</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cohanim</surname><given-names>A.B.</given-names></name><name><surname>Kashi</surname><given-names>Y.</given-names></name><name><surname>Trifonov</surname><given-names>E.N.</given-names></name></person-group><article-title>Three sequence rules for chromatin</article-title><source>J. Biomol. Struct. Dyn.</source><year>2006</year><volume>23</volume><fpage>559</fpage><lpage>566</lpage><pub-id pub-id-type="pmid">16494506</pub-id></citation></ref><ref id="b23"><label>23</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fukushima</surname><given-names>A.</given-names></name><name><surname>Ikemura</surname><given-names>T.</given-names></name><name><surname>Kinouchi</surname><given-names>M.</given-names></name><name><surname>Oshima</surname><given-names>T.</given-names></name><name><surname>Kudo</surname><given-names>Y.</given-names></name><name><surname>Mori</surname><given-names>H.</given-names></name><name><surname>Kanaya</surname><given-names>S.</given-names></name></person-group><article-title>Periodicity in prokaryotic and eukaryotic genomes identified by power spectrum analysis</article-title><source>Gene</source><year>2002</year><volume>300</volume><fpage>203</fpage><lpage>211</lpage><pub-id pub-id-type="pmid">12468102</pub-id></citation></ref><ref id="b24"><label>24</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Widom</surname><given-names>J.</given-names></name></person-group><article-title>Short-range order in two eukaryotic genomes: relation to chromosome structure</article-title><source>J. Mol. Biol.</source><year>1996</year><volume>259</volume><fpage>579</fpage><lpage>588</lpage><pub-id pub-id-type="pmid">8683566</pub-id></citation></ref><ref id="b25"><label>25</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Satchwell</surname><given-names>S.C.</given-names></name><name><surname>Drew</surname><given-names>H.R.</given-names></name><name><surname>Travers</surname><given-names>A.A.</given-names></name></person-group><article-title>Sequence periodicities in chicken nucleosome core DNA</article-title><source>J. Mol. Biol.</source><year>1986</year><volume>191</volume><fpage>659</fpage><lpage>675</lpage><pub-id pub-id-type="pmid">3806678</pub-id></citation></ref><ref id="b26"><label>26</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rivetti</surname><given-names>C.</given-names></name><name><surname>Guthold</surname><given-names>M.</given-names></name><name><surname>Bustamante</surname><given-names>C.</given-names></name></person-group><article-title>Scanning force microscopy of DNA deposited onto mica: equilibration versus kinetic trapping studied by statistical polymer chain analysis</article-title><source>J. Mol. Biol.</source><year>1996</year><volume>264</volume><fpage>919</fpage><lpage>932</lpage><pub-id pub-id-type="pmid">9000621</pub-id></citation></ref><ref id="b27"><label>27</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bustamante</surname><given-names>C.</given-names></name><name><surname>Marko</surname><given-names>J.F.</given-names></name><name><surname>Siggia</surname><given-names>E.D.</given-names></name><name><surname>Smith</surname><given-names>S.</given-names></name></person-group><article-title>Entropic elasticity of lambda-phage DNA</article-title><source>Science</source><year>1994</year><volume>265</volume><fpage>1599</fpage><lpage>1600</lpage><pub-id pub-id-type="pmid">8079175</pub-id></citation></ref><ref id="b28"><label>28</label><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vaillant</surname><given-names>C.</given-names></name><name><surname>Audit</surname><given-names>B.</given-names></name><name><surname>Arneodo</surname><given-names>A.</given-names></name></person-group><article-title>Thermodynamics of DNA loops with long-range correlated structural disorder</article-title><source>Phys Rev. Lett.</source><year>2005</year><volume>95</volume><fpage>068101</fpage><pub-id pub-id-type="pmid">16090995</pub-id></citation></ref></ref-list><sec sec-type="display-objects"><title>Figures and Tables</title><fig id="fig1" position="float"><label>Figure 1</label><caption><p>Sequences of DNA used in the AFM experiments. (<bold>a</bold>) Schematic representation of the PPC (2885) and sequence of the region analyzed (953 bp) that contains the hyperperiodic fragment &#x003a9;4 from <italic>C.elegans</italic> (923 bp). In the PPC the hyperperiodic region is located at both ends of the molecule, hence tracing of the molecule is independent of the particular end chosen as a starting point. Capital case letters denote the hyperperiodic fragment &#x003a9;4 from <italic>C.elegans</italic>. (<bold>b</bold>) Sequence of SQ100 fragment (185 bp, non-periodic) that contains the fragment &#x003a9;7 and (<bold>c</bold>) sequence of SQ322 (182 bp, hyperperiodic) that contains the fragment &#x003a9;6 from <italic>C.elegans</italic>. Two or more consecutive As or Ts were colored in green and red, respectively. The sequences were plotted in columns of 10 bases to highlight &#x0223c;10 bp periodicities. Capital letters in (a) denote the hyperperiodic fragment &#x003a9;4 from <italic>C.elegans</italic> and in (b and c) denote the &#x003a9;7 and &#x003a9;6 fragments, respectively.</p></caption><graphic xlink:href="gkl397f1"/></fig><fig id="fig2" position="float"><label>Figure 2</label><caption><p>Migration of hyperperiodic and control DNA in agarose gels. (<bold>a</bold>) Two-dimensional 4% NuSieve agarose gel electrophoresis of lambda DNA fragments (cut with NciI and PvuII) and (<bold>b</bold>) lambda fragments plus &#x003a9;4 fragments cut with the same restriction enzymes. Gels were run first at 4&#x000b0;C in the absence of intercalator (Ethidium Bromide, vertical) and then at 37&#x000b0;C in the presence of intercalator in the perpendicular direction (horizontal). The hyperperiodic &#x003a9;4 fragments (indicated by arrows) migrated more slowly than one would expect on the basis of its length. Similar anomalies were observed when using 2% agarose gels (data not shown). Gel electrophoresis on fragments SQ100 (185 bp, non-periodic, lane 3) and SQ322 (182 bp, hyperperiodic, lane 4) run at 37&#x000b0;C with EtBr (<bold>c</bold>) and at 4&#x000b0;C without EtBr (<bold>d</bold>). Lanes 1, 2, 5 and 6 correspond to 100 and 25 bp DNA ladders. At low temperature, the hyperperiodic fragment migrates more slowly than expected.</p></caption><graphic xlink:href="gkl397f2"/></fig><fig id="fig3" position="float"><label>Figure 3</label><caption><p>High-resolution AFM images of DNA samples. (<bold>a</bold>) Two representative DNA molecules of the PPC adsorbed on a mica surface. (<bold>b</bold>) Illustration of the data analysis and tracing of the DNA molecules. A segmented line made of points separated by &#x02113; = 2.5 nm follows the contour of the molecule. The important parameters we employ are the contour length spacing <italic>L</italic>, the end-to-end distance <italic>R</italic> and the curvature <italic>C</italic> defined as the inverse of the curvature radius <italic>r</italic> that embraces three contour points, <italic>C =</italic> 1<italic>/r</italic>. In the inset in (<bold>b</bold>), these parameters are shown together with an example of three values of <italic>R</italic> for the same contour length spacing <italic>L</italic>, in this case <italic>L</italic> = 3 &#x000d7; &#x02113; = 7.5 nm. We also show an example for the calculation of <italic>r</italic> for three points separated by <italic>L</italic> = 2 &#x000d7; &#x02113; = 5 nm (red arc). (<bold>c</bold>) AFM image of sample SQ100 and histogram of contour lengths <italic>L</italic><sub>0</sub> [inset in (<bold>c</bold>)]. (<bold>d</bold>) AFM image of sample SQ322 and histogram of contour lengths <italic>L</italic><sub>0</sub> [inset in (d)].</p></caption><graphic xlink:href="gkl397f3"/></fig><fig id="fig4" position="float"><label>Figure 4</label><caption><p>Analysis of unusual character in hyperperiodic DNA from <italic>C.elegans</italic> using AFM. (<bold>a</bold> and <bold>b</bold>) Plots of &#x0003c;<italic>R</italic><sup>2</sup>&#x0003e; versus contour length spacing <italic>L</italic> of the ends of the PPC, control DNA, SQ100 and SQ322 samples. Different sets of data in (<bold>a</bold>) correspond to different experiments, mica surfaces and sample preparations. In total, 190 PPC fragments were studied. The control DNA follows the prediction given by the worm-like-chain model with a persistence length of 51 nm. The hyperperiodic DNA regions from the <italic>C.elegans</italic> genome clearly show a very different behavior, characterized by lower &#x0003c;<italic>R<sup>2</sup></italic>&#x0003e;. To guide the eye, predictions from WLC with different values for the persistence length are included (solid grey). The five curves included are for <italic>P</italic> = 20, 30, 40, 50 and 60 nm. The slope of the curves increases as <italic>P</italic> increases. Bending of the DNA, is evidenced by lower &#x0003c;<italic>R<sup>2</sup></italic>&#x0003e; than predicted by WLC. (<bold>c</bold>) Normalized histograms of absolute curvature |<italic>C</italic>| for the PPC and the control DNA sample for <italic>L</italic> = 5 nm and for <italic>L</italic> = 25 nm [inset in (c)]. (<bold>d</bold>) Normalized absolute curvature |<italic>C</italic>| histograms for SQ100 and SQ322 for <italic>L</italic> = 5 nm. Distributions were fit to a Gaussian function centered at zero. The hyperperiodic DNA consistently showed larger standard deviations, independently of the spacing <italic>L</italic> [inset in (c)]. Bins and associated errors (square root of number of points in the bin) were normalized by the area of the histogram. We note that occasionally one cannot distinguish the error bars because they are smaller than the points.</p></caption><graphic xlink:href="gkl397f4"/></fig><fig id="fig5" position="float"><label>Figure 5</label><caption><p>Measurement of the intrinsic curvature in hyperperiodic <italic>C.elegans</italic> DNA. (<bold>a</bold>) Curvature correlation function (<xref ref-type="disp-formula" rid="e2">Equation 2</xref> in the main text) for PPC and Control DNA calculated for curvatures measured at contour length separations of <italic>L</italic> = 5, 15 and 25 nm. A long-range persistence of handedness in hyperperiodic DNA is evident from the graph, indicating the presence of intrinsic bending. (<bold>b</bold>) Normalized histograms of absolute mean curvatures per molecule |&#x0003c;<italic>C</italic>&#x0003e;<sub>molecule</sub>| for the PPC and the control DNA sample for <italic>L</italic> = 5 nm and for <italic>L</italic> = 25 nm [inset in (b)]. (<bold>c</bold>) |&#x0003c;<italic>C</italic>&#x0003e;<sub>molecule</sub>| histograms for SQ100 and SQ322 for <italic>L</italic> = 5 nm. Hyperperiodic DNA showed a preferred mean curvature per molecule of &#x0223c;0.01 nm<sup>&#x02212;1</sup> at <italic>L</italic> = 5 nm, consistent with the presence of intrinsic bending. Bins and associated errors (square root of number of points in the bin) were normalized by the area of the histogram.</p></caption><graphic xlink:href="gkl397f5"/></fig><table-wrap id="tbl1" position="float"><label>Table 1</label><caption><p>Relevant DNAs used in this study</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left" colspan="1" rowspan="1">DNA (size)</th><th align="left" colspan="1" rowspan="1">Gene in <italic>C.elegans</italic></th><th align="left" colspan="1" rowspan="1">Periodicity value Ave, (Frac&#x0003e;95 cutoff)</th><th align="left" colspan="1" rowspan="1">AFM sample (size)</th></tr></thead><tbody><tr><td align="left" colspan="1" rowspan="1">Fragment &#x003a9;4 (932 bp)</td><td align="left" colspan="1" rowspan="1">F54C4.1</td><td align="left" colspan="1" rowspan="1">433, (89.7%)</td><td align="left" colspan="1" rowspan="1">PPC (2885 bp)</td></tr><tr><td align="left" colspan="1" rowspan="1">Fragment &#x003a9;6 (147 bp)</td><td align="left" colspan="1" rowspan="1">ZK792.7</td><td align="left" colspan="1" rowspan="1">122, (76.9%)</td><td align="left" colspan="1" rowspan="1">SQ322 (182 bp)</td></tr><tr><td align="left" colspan="1" rowspan="1">Fragment &#x003a9;7 (150 bp)</td><td align="left" colspan="1" rowspan="1">Y38F1A.3</td><td align="left" colspan="1" rowspan="1">0.93, (0%)</td><td align="left" colspan="1" rowspan="1">SQ100 (185 bp)</td></tr><tr><td align="left" colspan="1" rowspan="1">Control Pgem-3Z (2743 bp)</td><td align="left" colspan="1" rowspan="1">&#x02014;</td><td align="left" colspan="1" rowspan="1">3.21, (0%)</td><td align="left" colspan="1" rowspan="1">Control (2743 bp)</td></tr></tbody></table><table-wrap-foot><fn><p>Periodicity values were calculated using the algorithm described in Ref. (<xref ref-type="bibr" rid="b19">19</xref>).</p></fn></table-wrap-foot></table-wrap><table-wrap id="tbl2" position="float"><label>Table 2</label><caption><p>Standard deviation &#x003c3;<italic><sub>C</sub></italic> of curvature distributions and standard deviation &#x003c3;<sub>&#x0003c;<italic>C</italic>&#x0003e;molecule</sub> and intrinsic curvature &#x0003c;<italic>C</italic>&#x0003e;<sub>0,molecule</sub> of mean curvature distributions</p></caption><table frame="hsides" rules="groups"><thead><tr><th colspan="1" rowspan="1"/><th align="left" colspan="1" rowspan="1">Control</th><th align="left" colspan="1" rowspan="1">PPC</th><th align="left" colspan="1" rowspan="1">SQ100</th><th align="left" colspan="1" rowspan="1">SQ322</th></tr></thead><tbody><tr><td align="left" colspan="1" rowspan="1">&#x003c3;<sub><italic>C,L</italic>=5 nm</sub>(nm<sup>&#x02212;1</sup>)</td><td align="left" colspan="1" rowspan="1">0.0547 &#x000b1; 0.0007</td><td align="left" colspan="1" rowspan="1">0.0618 &#x000b1; 0.0005</td><td align="left" colspan="1" rowspan="1">0.0667 &#x000b1; 0.0004</td><td align="left" colspan="1" rowspan="1">0.0737 &#x000b1; 0.0004</td></tr><tr><td align="left" colspan="1" rowspan="1">&#x003c3;<sub><italic>C, L</italic>=25 nm</sub> (nm<sup>&#x02212;1</sup>)</td><td align="left" colspan="1" rowspan="1">0.0241 &#x000b1; 0.0003</td><td align="left" colspan="1" rowspan="1">0.0298 &#x000b1; 0.0003</td><td align="left" colspan="1" rowspan="1">&#x02014;</td><td align="left" colspan="1" rowspan="1">&#x02014;</td></tr><tr><td align="left" colspan="1" rowspan="1">&#x003c3;<sub>&#x0003c;C&#x0003e;molecule, L=5 nm</sub> (nm<sup>&#x02212;1</sup>)</td><td align="left" colspan="1" rowspan="1">0.0077 &#x000b1; 0.0003</td><td align="left" colspan="1" rowspan="1">0.0072 &#x000b1; 0.0008</td><td align="left" colspan="1" rowspan="1">0.0163 &#x000b1; 0.0005</td><td align="left" colspan="1" rowspan="1">0.024 &#x000b1; 0.002</td></tr><tr><td align="left" colspan="1" rowspan="1">&#x0003c;C&#x0003e;<sub>0,molecule, L=5 nm</sub> (nm<sup>&#x02212;1</sup>)</td><td align="left" colspan="1" rowspan="1">0</td><td align="left" colspan="1" rowspan="1">&#x000b1;0.010 &#x000b1; 0.001</td><td align="left" colspan="1" rowspan="1">0</td><td align="left" colspan="1" rowspan="1">&#x000b1;0.010 &#x000b1; 0.002</td></tr><tr><td align="left" colspan="1" rowspan="1">&#x003c3;<sub>&#x0003c;C&#x0003e;molecule, L=25 nm</sub> (nm<sup>&#x02212;1</sup>)</td><td align="left" colspan="1" rowspan="1">0.007 &#x000b1; 0.001</td><td align="left" colspan="1" rowspan="1">0.008 &#x000b1; 0.001</td><td align="left" colspan="1" rowspan="1">&#x02014;</td><td align="left" colspan="1" rowspan="1">&#x02014;</td></tr><tr><td align="left" colspan="1" rowspan="1">&#x0003c;C&#x0003e;<sub>0,molecule, L=25 nm</sub> (nm<sup>&#x02212;1</sup>)</td><td align="left" colspan="1" rowspan="1">0</td><td align="left" colspan="1" rowspan="1">&#x000b1;0.013 &#x000b1; 0.001</td><td align="left" colspan="1" rowspan="1">&#x02014;</td><td align="left" colspan="1" rowspan="1">&#x02014;</td></tr></tbody></table><table-wrap-foot><fn><p>Values correspond to Gaussian fits to distributions shown in <xref ref-type="fig" rid="fig4">Figures 4c</xref>, <xref ref-type="fig" rid="fig4">4d</xref>, <xref ref-type="fig" rid="fig5">5b</xref> and <xref ref-type="fig" rid="fig5">5c</xref>. The segment length was 5 and 25 nm, where indicated. Distributions for the short DNAs were not calculated at <italic>L</italic> = 25 nm.</p></fn></table-wrap-foot></table-wrap></sec></back></article>