基于DNA变形能的核小体定位预测方法研究进展

真核细胞中,作为染色质基本结构单元的核小体参与调控基因的转录、DNA复制、重组以及RNA剪接等诸多生物学过程。阐明核小体定位机制并准确预测核小体在染色体上的位置对解读染色质结构与功能有重要生物学意义。在过去30多年时间里,研究人员发展了多种预测核小体位置的方法。最理想的方法应考虑DNA序列、组蛋白修饰和染色质重塑等影响核小体定位的诸多因素,然而现实中,捕捉主要因素的模型也往往具有很高的鲁棒性和实用价值。DNA序列偏好性是在全基因组尺度上影响核小体定位的最重要因素之一,因此基于DNA序列的核小体定位预测方法也最常见。这种方法可大致分为两类,即基于DNA序列信息的生物信息学模型和基于DNA变形能的生物物理学模型。本文重点介绍生物物理学模型近些年取得的主要进展。     

The Molecular Evolution of Nucleosome Positioning Through Sequence-Dependent Deformation of the DNA Polymer

Abstract

The computational prediction of nucleosome positioning from DNA sequence now allows for in silico investigation of the molecular evolution of biophysical properties of the DNA molecule responsible for primary chromatin organization in the genome. To discern what signal components driving nucleosome positioning in the yeast genome are potentially targeted by natural selection, we compare the performance of various models predictive of nucleosome positioning within the context of a simple statistical test, the repositioned mutation test. We demonstrate that while nucleosome occupancy is driven largely by translational exclusion in response to AT content, there is also a strong signature of evolutionary conservation of regular patterns within nucleosomal DNA sequence related to the structural organization of the nucleosome core (e.g., 10-bp dinucleotide periodicity). We also use computer simulations to investigate hypothetical coding and regulatory constraints on the ability of sequence properties affecting nucleosome formation to adaptively evolve. Our results demonstrate that natural selection may act independently on different DNA sequence properties responsible for local chromatin organization. Furthermore, at least with respect to the deformation energy of the DNA molecule in the nucleosome, the presence of the genetic code has greatly restricted the ability of sequences to evolve the dynamic nucleosome organization typically observed in promoter regions.     

Structure-based Analysis of DNA Sequence Patterns Guiding Nucleosome Positioning in vitro

Abstract

Recent studies of genome-wide nucleosomal organization suggest that the DNA sequence is one of the major determinants of nucleosome positioning. Although the search for underlying patterns encoded in nucleosomal DNA has been going on for about 30 years, our knowledge of these patterns still remains limited. Based on our evaluations of DNA deformation energy, we developed new scoring functions to predict nucleosome positioning. There are three principal differences between our approach and earlier studies: (i) we assume that the length of nucleosomal DNA varies from 146 to 147 bp; (ii) we consider the anisotropic flexibility of pyrimidine-purine (YR) dimeric steps in the context of their neighbors (e.g., YYRR versus RYRY); (iii) we postulate that alternating AT-rich and GC-rich motifs reflect sequence-dependent interactions between histone arginines and DNA in the minor groove. Using these functions, we analyzed 20 nucleosome positions mapped in vitro at single nucleotide resolution (including clones 601, 603, 605, the pGUB plasmid, chicken β-globin and three 5S rDNA genes). We predicted 15 of the 20 positions with 1-bp precision, and two positions with 2-bp precision. The predicted position of the ‘601’ nucleosome (i.e., the optimum of the computed score) deviates from the experimentally determined unique position by no more than 1 bp—an accuracy exceeding that of earlier predictions.

Our analysis reveals a clear heterogeneity of the nucleosomal sequences which can be divided into two groups based on the positioning ‘rules’ they follow. The sequences of one group are enriched by highly deformable YR/YYRR motifs at the minor-groove bending sites SHL ±3.5 and ±5.5, which is similar to the α-satellite sequence used in most crystallized nucleosomes. Apparently, the positioning of these nucleosomes is determined by the interactions between histones H2A/H2B and the terminal parts of nucleosomal DNA. In the other group (that includes the ‘601’ clone) the same YR/YYRR motifs occur predominantly at the sites SHL ±1.5. The interaction between the H3/H4 tetramer and the central part of the nucleosomal DNA is likely to be responsible for the positioning of nucleosomes of this group, and the DNA trajectory in these nucleosomes may differ in detail from the published structures.

Thus, from the stereochemical perspective, the in vitro nucleosomes studied here follow either an X-ray-like pattern (with strong deformations in the terminal parts of nucleosomal DNA), or an alternative pattern (with the deformations occurring predominantly in the central part of the nucleosomal DNA). The results presented here may be useful for genome-wide classification of nucleosomes, linking together structural and thermodynamic characteristics of nucleosomes with the underlying DNA sequence patterns guiding their positions.     

体外装配核小体过程组蛋白浓度依赖的动力学模型

为探索组蛋白浓度对核小体体外装配的影响,本研究表达纯化了4种组蛋白,通过控制实验反应体系中组蛋白的浓度,利用盐透析法在体外装配了核小体,检测分析了组蛋白浓度与核小体组装效率的关系。以此实验数据为基础,提出了核小体组装过程组蛋白浓度依赖性的动力学模型。实验结果发现,反应体系中组蛋白浓度与核小体生成量呈典型的线性关系。依据动力学理论模型,进行线性回归拟合,回归系数达到0.963;经计算601 DNA序列组装核小体的反应速率常数k为1.49×10^-5mL·h·μg^-1。CS1序列验证动力学模型的线性回归相关系数为0.989,反应速率常数为1.52×10^-5mL·h·μg^-1。该实验方法及动力学模型中反应速率常数k可用于评价相同长度的DNA序列组装核小体的能力、组蛋白与其突变体以及组蛋白变体之间形成核小体结构能力的差异。该动力学模型的建立为理解核小体装配、核小体定位、染色质结构等相关问题提供了理论指导。     

Nucleosome Positioning,Nucleosome Spacing and the Nucleosome Code

Abstract

Nucleosome positioning has been the subject of intense study for many years. The properties of micrococcal nuclease, the enzyme central to these studies, are discussed. The various methods used to determine nucleosome positions in vitro and in vivo are reviewed critically. These include the traditional low resolution method of indirect end-labelling, high resolution methods such as primer extension, monomer extension and nucleosome sequencing, and the high throughput methods for genome-wide analysis (microarray hybridisation and parallel sequencing). It is established that low resolution mapping yields an averaged chromatin structure, whereas high resolution mapping reveals the weighted superposition of all the chromatin states in a cell population. Mapping studies suggest that yeast DNA contains information specifying the positions of nucleosomes and that this code is made use of by the cell. It is proposed that the positioning code facilitates nucleosome spacing by encoding information for multiple alternative overlapping nucleosomal arrays. Such a code might facilitate the shunting of nucleosomes from one array to another by ATP-dependent chromatin remodelling machines.     

Nucleosome DNA Bendability Matrix (C. elegans)

Abstract

An original signal extraction procedure is applied to database of 146 base nucleosome core DNA sequences from C. elegans (S. M. Johnson et al. Genome Research 16, 1505–1516, 2006). The positional preferences of various dinucleotides within the 10.4 base nucleosome DNA repeat are calculated, resulting in derivation of the nucleosome DNA bendability matrix of 16x10 elements. A simplified one-line presentation of the matrix (“consensus” repeat) is…A(TTTCCGGAAA)T…. All 6 chromosomes of C. elegans conform to the bendability pattern. The strongest affinity to their respective positions is displayed by dinucleotides AT and CG, separated within the repeat by 5 bases. The derived pattern makes a basis for sequence-directed mapping of nucleosome positions in the genome of C. elegans. As the first complete matrix of bendability available the pattern may serve for iterative calculations of the species-specific matrices of bendability applicable to other genomic sequences.     

Sequence-directed Mapping of Nucleosome Positions

Abstract

A nucleosome DNA sequence probe is designed that combines recently derived RR/YY counter-phase and AA/TT in-phase periodical patterns. A simple nucleosome mapping procedure is introduced for prediction of the nucleosome positions in the sequence of interest, to serve as a guide for experimental studies of the chromatin structure.     

Three Sequence Rules for Chromatin

Abstract

Extensive DNA sequence analysis of three eukaryotes, S. cerevisiae, C. elegans, and D. melanogaster, reveals two different AA/TT periodical patterns associated with the nucleosome positioning. The first pattern is the counter-phase oscillation of AA and TT dinucleotides, which has been frequently considered as the nucleosome DNA pattern. This represents the sequence rule I for chromatin structure. The second pattern is the in-phase oscillation of the AA and TT dinucleotides with the same nucleosome DNA period, 10.4 bases. This pattern apparently corresponds to curved DNA, that also participates in the nucleosome formation, and represents the sequence rule II for chromatin. The positional correlations of AA and TT dinucleotides also indicate that the nucleosomes are separated by specific linker sizes (preferably 8, 18,…bases), dictated by the steric exclusion rules. Thus, the sequence positions of the neighboring nucleosomes are correlated, and this represents the sequence rule III.     

基于核小体位置预测的酵母进化印迹研究

在利用核小体定位实验数据训练支持向量机(SVM)对任意酵母DNA序列的核小体形成能力进行预测的过程中,发现染色质结构对基因组DNA分子进化过程有着显著影响。我们观察到核小体DNA比连接DNA的平均预测准确率低15%,这种普遍存在的局部预测准确率差异性代表了酵母核小体定位的进化印迹(Evolutionary footprint),它揭示了核小体组织在基因组的整个进化过程中所具有的保守性。     

5-Bromouracil disrupts nucleosome positioning by inducing A-form-like DNA conformation in yeast cells

5-Bromodeoxyuridine (BrdU) modulates expression of particular genes associated with cellular differentiation and senescence. Our previous studies have suggested an involvement of chromatin structure in this phenomenon. Here, we examined the effect of 5-bromouracil on nucleosome positioning in vivo using TALS plasmid in yeast cells. This plasmid can stably and precisely be assembled nucleosomes aided by the α2 repressor complex bound to its α2 operator. Insertion of AT-rich sequences into a site near the operator destabilized nucleosome positioning dependent on their length and sequences. Addition of BrdU almost completely disrupted nucleosome positioning through specific AT-tracts. The effective AT-rich sequences migrated faster on polyacrylamide gel electrophoresis, and their mobility was further accelerated by substitution of thymine with 5-bromouracil. Since this property is indicative of a rigid conformation of DNA, our results suggest that 5-bromouracil disrupts nucleosome positioning by inducing A-form-like DNA.     

Histone hyperacetylation facilitates chromatin remodelling in a Drosophila embryo cell-free system

We found that Drosophila embryo extract contains a protein activity (or activities) that can destabilize nucleosomes, resulting in increased sensitivity to DNase I, release of nucleosomal supercoiling, high levels of conformational flexibility of DNA and more diffuse micrococcal nuclease digestion patterns. Incorporation of histone H1 did not significantly affect this nucleosome remodelling. Remodelling occurs more efficiently in hyperacetylated chromatin. It was shown previously that hyperacetylated chromatin, reconstituted in a Drosophila embryo cell-free system, exhibits increased DNase I sensitivity and a high degree of conformational flexibility of DNA. The present data suggest that the more diffuse structure of acetylated chromatin is a result of chromatin remodelling by protein activities in the Drosophila embryo extract. Received: 4 November 1998 / Accepted: 10 May 1999     

Prediction of nucleosome occupancy in Saccharomyces cerevisiae using position-correlation scoring function

Knowledge of the detailed organization of nucleosomes across genomes and the mechanisms of nucleosome positioning is critical for the understanding of gene regulation and expression. In the present work, the bias of 4-mer frequency in nucleosome and linker sequences of the S. cerevisiae genome was analyzed statistically. A novel position-correlation scoring function algorithm based on the bias of 4-mer frequency in linker sequences was presented to distinguish nucleosome vs linker sequences. Five-fold cross-validation demonstrated that the algorithm achieved a good performance with mean area under the receiver operator characteristics curve of 0.981. Next, the algorithm was used to predict nucleosome occupancy throughout the S. cerevisiae genome and relatively high correlation coefficients with experiment maps of nucleosome positioning were obtained. Besides, the distinct nucleosome depleted regions in the vicinity of regulatory sites were confirmed. The results suggest that intrinsic DNA sequence preferences in linker regions have a significant impact on the nucleosome occupancy.     

Brownian dynamics simulation of nucleosome formation and disruption under stretching

Using a Brownian dynamics simulation, we numerically studied the interaction of DNA with histone and proposed an octamer-rotation model to describe the process of nucleosome formation. Nucleosome disruption under stretching was also simulated. The theoretical curves of extension versus time as well as of force versus extension are consistent with previous experimental results.