{"id":704,"date":"2023-03-15T13:43:30","date_gmt":"2023-03-15T06:43:30","guid":{"rendered":"https:\/\/conf.icgbio.ru\/bgrs98\/?page_id=704"},"modified":"2023-04-11T14:33:45","modified_gmt":"2023-04-11T07:33:45","slug":"054_fine-structure-of-the-profile-of-bending-stiffness-energy-of-nucleosomal-dna","status":"publish","type":"page","link":"https:\/\/conf.icgbio.ru\/bgrs98\/abstracts\/abstract-list\/054_fine-structure-of-the-profile-of-bending-stiffness-energy-of-nucleosomal-dna\/","title":{"rendered":"FINE STRUCTURE OF THE PROFILE OF BENDING STIFFNESS ENERGY OF NUCLEOSOMAL DNA"},"content":{"rendered":"

LEVITSKY V.G.<\/a><\/p>\n

Laboratory of Theoretical Genetics, Institute of Cytology and Genetics, Siberian Branch of the Russian Academy of Sciences, 10 Lavrentiev Ave., Novosibirsk, 630090, Russia; e-mail:\u00a0levitsky@bionet.nsc.ru<\/a><\/p>\n

Keywords<\/a>: bending stiffness energy, nucleosomal DNA<\/p>\n

 <\/p>\n

Nucleosome is an elementary unit of chromatin packaging and are able to play an active role in regulation of gene expression [1]. A minimal nucleosome (core) comprises eight histones: two dimers H2A-H2B and tetramer (H3-H4)2<\/sub>. Nucleosomal DNA wrapped around the core has a shape of solenoid containing 1.8 coils that require 145 bp. [2]. A histone H1 attached to the core-DNA complex completes the second coil on its surface. Histone H1 binds 20 additional base pairs and contributes to further chromatin compactization. Thus, a nucleosome binds about 165 bp in a more compact state of chromatin packaging. Analysis of nucleosomal DNA has so far demonstrated the existence of a weak context signal providing for precise nucleosome positioning [3-5]. According to recent X-ray structure studies, the DNA has detectable bends on the surface of nucleosome [2]. Hence, the ability of DNA to adopt certain conformation can enhance considerably its binding to nucleosome. That is why it is especially interesting to find the DNA conformational properties that are significant for nucleosome formation.<\/p>\n

The fine structure of the profiles of bending stiffness energy of nucleosomal DNA was studied in this work. The approach based on employment of DNA persistent lengths for dinucleotide pairs [6] was used. The persistent length of DNA characterizes the length within which the DNA axis remains detectably constant. In case of dinucleotides, the parameters of persistent length [6] reflect the local contribution of a dinucleotide in the persistent length of more extended DNA region. We calculated the energy of DNA bending stiffness\u00a0E\u00a0<\/b>as described in [6]:<\/p>\n\n\n\n
\"\"<\/a><\/td>\n(1)<\/b><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

where\u00a0R\u00a0<\/b>= 4.3 nm is the radius of curvature of nucleosomal DNA;\u00a0k,<\/b>\u00a0Boltzmann constant;\u00a0T,<\/b>\u00a0absolute temperature;\u00a0h<\/b>\u00a0= 0.34 nm, the distance between two adjacent base pairs;\u00a0W,<\/b>\u00a0the length (bp) of a region to calculate the energy;\u00a0pi,<\/b>\u00a0the dinucleotide parameters of persistent length (Table 1) that reflect the linear dependence between the DNA persistent length and its melting temperature [6]; and\u00a0i,<\/b>\u00a0position of a dinucleotide in the DNA region considered. Equation\u00a0(1)<\/b>\u00a0estimates the work necessary to wrap the chain containing\u00a0W-1<\/b>\u00a0dinucleotide links around the nucleosome. Thus, the smaller is the value\u00a0E,<\/b>\u00a0the more flexible is the DNA region of length\u00a0W.<\/b><\/p>\n

Equation\u00a0(1)<\/b>\u00a0was used to analyze 143 nucleotide sequences of nucleosome binding sites, described in the database [7]. These sites are located in various regions of the eukaryotic genome: 3′, 5′, and inner regions of the protein-encoding genes; rRNA and tRNA genes; satellite and centromeric regions of chromosomes. The site nucleotide sequences 185 bp in length (with left and right flanks of 92 bp relative to the center of the site) were extracted form the EMBL database using the codes described in [7]. Selection of such sequence length, exceeding the size of nucleosomal DNA (145 bp), is substantiated by the necessity to take into consideration possible experimental errors in determining the position of the site’s center [9]. The profile of DNA bending stiffness energy\u00a0E(X)<\/b>\u00a0was calculated by Equation\u00a0(1)<\/b>\u00a0for each sequence analyzed at all the possible\u00a0L-W+1<\/b>\u00a0positions of the scanning window, where\u00a0L<\/b>\u00a0is the length of the sequence;\u00a0W,<\/b>\u00a0size of the window; and\u00a0X,<\/b>\u00a0position of the center.<\/p>\n

The profiles of DNA bending stiffness of all nucleosome sites [7] obtained for both orientations of nucleotide sequences were summed pointwise from 2 to 150 at window lengths being multiple of two. Thus, the symmetrized profiles\u00a0E(X),<\/b>\u00a0mean for the sample, were determined at each length of the window (see examples in Fig. 1). Note that the profiles contain maximums and minimums, where the value\u00a0E(X)<\/b>\u00a0differs considerably from the means averaged over the profile. Note also the dependence of the position and markedness of the extremums on the window length. The following analysis was carried out to reveal most significant minimums and maximums of the profiles\u00a0E(X).<\/b>\u00a0The mean value\u00a0Emean<\/sub><\/b> and dispersion \"\"\u00a0were determined for the profile-forming parameters\u00a0{\u00c5(\u00d5)}<\/b>\u00a0at each window size. These parameters were used to construct the approximation of\u00a0{\u00c5(\u00d5)}<\/b>\u00a0distribution to the normal distribution\u00a0N(Emean<\/sub>,\"\"<\/b>\u00a0).<\/b>\u00a0The test on normality of distributions was positive in 72 of 75 window sizes considered; thus, the energy distribution along the profiles were considered normal. Minimum\u00a0Emin<\/b>\u00a0and maximum\u00a0Emax<\/b>\u00a0within the profile\u00a0E(X)<\/b>\u00a0were determined at each window size as well as the corresponding positions of the window on the DNA. Then the distribution\u00a0N(Emean<\/sub>, \"\"<\/b>\u00a0)<\/b>\u00a0was used to estimate the probabilities\u00a0\u00d01(\u00c5<Emin)<\/b>\u00a0and\u00a0\u00d02(\u00c5>Emax)<\/b>\u00a0to obtain accidentally the energy\u00a0E<\/b>\u00a0that is lower or higher than the minimal value\u00a0Emin<\/b>\u00a0or maximal value\u00a0Emax\u00a0<\/b>observed for this profile, respectively, were estimated. Depencences of the probabilities\u00a0P1<\/b>\u00a0and\u00a0P2<\/b>\u00a0on the window length are shown in Fig. 2. Eight local minimums for each of the probabilities\u00a0P1<\/b>\u00a0and\u00a0P2<\/b>, corresponding to different lengths of the window\u00a0W<\/b>, were revealed.<\/p>\n

It is evident from Fig 2. that the probability\u00a0P1<\/b>\u00a0is minimal for the windows 8, 30, 34, 40, 50, 80, 116, and 122 bp long, the lowest minimums being for the windows of 8 and 166 bp long. The position of the center of nucleosome site was taken as the origin. The centers of the regions with an increased flexibility at a 8-bp window are located at positions -21 and +21 (Fig. 1). Calculations demonstrate that these regions [-25; -17] and [+17; +25] are most flexible in the nucleosomal DNA (Fig. 1). The regions with increased flexibility revealed at the window sizes of 30, 34, 40, and 50 bp are located within one region [-26; +26]. Note that the regions revealed at a 8-bp window correspond to the ends of the region [-26; +26]. In turn, the latter region is a constituent of the regions with increased flexibility that are revealed at a longer window. Thus, the hierarchy of the regions of different lengths with increased DNA flexibility is observed within a nucleosome site.<\/p>\n

As for the probability\u00a0P2,<\/b>\u00a0it is minimal at window lengths of 6, 66, 70, 84, 88, 96, 110, and 148 bp. (Fig. 2). The least probabilities\u00a0P2<\/b>\u00a0were found for the windows of 6 and 110 bp. Regions of an increased stiffness appear at a window of 6 bp; their centers are located at positions -43 and +43 (Fig. 1). According to the calculations, it is these regions [-46; +40] and [+40; +46] that are the most stiff within the nucleosome site (Fig. 1). At a window of 60 bp, two regions with an increased stiffness were found: they are located at the right and left ends of the site at positions [-92; -26] and [+26; +92]. If the window length is increased over 66 bp, two regions of the increased DNA stiffness, also located at the ends of the site, are coming closer and overlap when the window length exceeds 96 bp. Thus, the hierarchy of the regions of different lengths and locations with increased DNA stiffness is also observed in this case.<\/p>\n

Returning to the profile constructed for the 8-bp window (Fig. 1), note the regions of increased flexibility [-56; +64] and [+56; +64] in addition to the most flexible regions [-17; -25] and [+17; +25]. Similarly, the regions of high stiffness [-80; -72] and [+72; +80] are also evident in addition to the regions of the highest stiffness [-48; -40] and [+40; +48]. Thus, we are observing that the regions with nonaccidentally high and low flexibility are located in a contrasting vicinity, which is likely to provide the necessary conditions for DNA bending and precise positioning of the nucleosome. Analysis of the profile at a window of 8 bp in length also demonstrates that the central region about 20 bp long fail to manifest extremal bending stiffness. At longer windows, note the region of DNA flexibility [-25; +25] (50-bp window) and stiffness [-92; -26] and [+26; +92] (66 bp window).<\/p>\n

Spatial location on the surface of nucleosome of the above-described most stiff and most flexible regions revealed at 8-bp window are shown in Fig. 3. Broken line indicate the region of DNA flexibility [-25; +25] revealed at 50-bp window; solid line, the regions of DNA stiffness [-92; -26] and [+26; +92] revealed at 66-bp window.<\/p>\n

The conclusions of our study are consistent with the data on analysis of nucleosomes by other methods. As for the central region of the site about 20 bp in length, recent studies of nucleosomal DNA by X-ray analysis [2] and multiple alignment [4] confirm that this region is not bent. Significant periodicities of the DNA nucleosomal context were revealed in [4] in other regions of the site. Discussing the shape of bending stiffness profile of nucleosomal DNA (Fig. 1; 8-bp window), note also that it fits the shapes of the profiles of nucleosomal DNA obtained using DNA geometric parameters [5]. In that work, the authors suggest that the flexibility pattern of nucleosomal DNA consists of two bent elements 50-60 bp each that are separated by a region of smaller curvature 20-30 bp long.<\/p>\n

The analysis performed allowed us to obtain the additional information on structural peculiarities of nucleosomal DNA. The lengths and locations of the regions within the nucleosome site that correspond to the most significant minimums and maximums of DNA bending stiffness were determined, and their hierarchy was demonstrated.<\/p>\n

 <\/p>\n

Table 1. Dinucleotide parameters of DNA persistent length [6]<\/p>\n

\n

 <\/p>\n\n\n\n\n\n\n\n\n\n\n\n\n\n
\n

Dinucleotide pair<\/p>\n<\/td>\n

\n

\u00f0 (nm)<\/p>\n<\/td>\n<\/tr>\n

\n

TA<\/p>\n<\/td>\n

\n

11<\/p>\n<\/td>\n<\/tr>\n

\n

AT<\/p>\n<\/td>\n

\n

41<\/p>\n<\/td>\n<\/tr>\n

\n

AA\/TT<\/p>\n<\/td>\n

\n

37<\/p>\n<\/td>\n<\/tr>\n

\n

CA\/TG<\/p>\n<\/td>\n

\n

38<\/p>\n<\/td>\n<\/tr>\n

\n

AG\/CT<\/p>\n<\/td>\n

\n

43<\/p>\n<\/td>\n<\/tr>\n

\n

GA\/TC<\/p>\n<\/td>\n

\n

84<\/p>\n<\/td>\n<\/tr>\n

\n

AC\/GT<\/p>\n<\/td>\n

\n

100<\/p>\n<\/td>\n<\/tr>\n

\n

CG<\/p>\n<\/td>\n

\n

63<\/p>\n<\/td>\n<\/tr>\n

\n

GC<\/p>\n<\/td>\n

\n

156<\/p>\n<\/td>\n<\/tr>\n

\n

GG\/CC<\/p>\n<\/td>\n

\n

83<\/p>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

 <\/p>\n<\/div>\n

 <\/p>\n

\"\"<\/a><\/p>\n

Fig. 1. The profiles\u00a0E(X)\u00a0<\/b>of DNA bending stiffness at window lengths of 8 and 50 bp:\u00a0X,\u00a0<\/b>position of the window relative to the center of the site;\u00a0E,\u00a0<\/b>in the units of persistent length.<\/p>\n

\"\"<\/a><\/p>\n

Fig. 2. The probabilities\u00a0\u00d01(\u00c5<Emin)\u00a0<\/b>and\u00a0\u00d02(E>Emax)<\/b>\u00a0representing the significance of extremums of profiles\u00a0E(X)<\/b>\u00a0depending on the window length\u00a0W.<\/b><\/p>\n

 <\/p>\n

\"\"<\/a><\/p>\n

Fig. 3. Spatial location of the revealed regions contrasting in DNA flexibility on the histone octamer: (1 and 2) the most stiff regions and (3) the regions with increased stiffness; (4 and 5) the most flexible regions and (6) the regions with increased flexibility.<\/p>\n

 <\/p>\n

This work was supported by grants from the Russian Foundation for Basic Research (No.97-04-49740, 97-07-90309, 96-04-50006, 98-04-49479, 98-07-90126); Russian Ministry of Science and Technologies; Russian Human Genome Project; Russian Ministry of High Education; Siberian Department of RAS (Programms for support of reseach of young scientists and Programm of Integration projects); National Institutes of Health, U.S.A. (No.5-R01-RR-04026-08)<\/p>\n

References<\/b><\/p>\n

    \n
  1. F. Thoma,. \u0166Nucleosome positioning (Review)\u0167 Biochim. Biophys. Acta\u00a01130<\/b>, 1-19 (1992).<\/li>\n
  2. K. Luger, A.W. Mader, R.K. Richmond, D.F. Sargent, and T.J. Richmond, \u0166Crystal structure of the nucleosome core particle at 2.8 A resolution\u0167 Nature\u00a0389<\/b>, 251-260 (1997).<\/li>\n
  3. S. Satchwell, H.R. Drew, and A.A. Travers, \u0166Sequence periodicities in chicken nucleosome core DNA.\u0167 J. Mol. Biol.\u00a0191<\/b>, 659-675 (1986).<\/li>\n
  4. I. Ioshikhes, E.N. Trifonov, et al., \u0166Nucleosomal DNA sequence database\u0167 J. Mol. Biol.\u00a0262<\/b>, 129-139. (1996).<\/li>\n
  5. D.J. Fitzgerald, et al., \u0166Conserved patterns of bending in satellite and nucleosome DNA\u0167 J. Biol. Chem.\u00a0269<\/b>, 21303-21314 (1994).<\/li>\n
  6. A.V. Sivolob and S.N. Kharpunov, \u0166Translational positioning of nucleosomes on DNA: the role of sequence-dependent isotropic DNA bending stiffness\u0167 J. Mol. Biol.\u00a0247<\/b>, 918-931 (1995).<\/li>\n
  7. I. Ioshikhes and E.N. Trifonov, \u0166Nucleosome DNA sequence pattern revealed by multiple alignment of experimentally mapped sequences\u0167 Nucl. Acids. Res.\u00a021<\/b>, 4857-4859 (1993)<\/li>\n
  8. The profiles of DNA bending stiffness used in the work were built using the program Craft available at\u00a0http:\/\/wwwmgs.bionet.nsc.ru\/Programs\/CRAFT\/<\/a><\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"

    LEVITSKY V.G. Laboratory of Theoretical Genetics, Institute of Cytology and Genetics, Siberian Branch of the Russian Academy of Sciences, 10 Lavrentiev Ave., Novosibirsk, 630090, Russia; e-mail:\u00a0levitsky@bionet.nsc.ru Keywords: bending stiffness energy, nucleosomal DNA   Nucleosome is an elementary unit of chromatin … Continue reading →<\/span><\/a><\/p>\n","protected":false},"author":13,"featured_media":0,"parent":97,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":[],"_links":{"self":[{"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages\/704"}],"collection":[{"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/users\/13"}],"replies":[{"embeddable":true,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/comments?post=704"}],"version-history":[{"count":5,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages\/704\/revisions"}],"predecessor-version":[{"id":1391,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages\/704\/revisions\/1391"}],"up":[{"embeddable":true,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages\/97"}],"wp:attachment":[{"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/media?parent=704"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}