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.     

Plasticity and Epigenetic Inheritance of Centromere-specific Histone H3 (CENP-A)-containing Nucleosome Positioning in the Fission Yeast

Nucleosomes containing the specific histone H3 variant CENP-A mark the centromere locus on each chromatin and initiate kinetochore assembly. For the common type of regional centromeres, little is known in molecular detail of centromeric chromatin organization, its propagation through cell division, and how distinct organization patterns may facilitate kinetochore assembly. Here, we show that in the fission yeast S. pombe, a relatively small number of CENP-A/Cnp1 nucleosomes are found within the centromeric core and that their positioning relative to underlying DNA varies among genetically homogenous cells. Consistent with the flexible positioning of Cnp1 nucleosomes, a large portion of the endogenous centromere is dispensable for its essential activity in mediating chromosome segregation. We present biochemical evidence that Cnp1 occupancy directly correlates with silencing of the underlying reporter genes. Furthermore, using a newly developed pedigree analysis assay, we demonstrated the epigenetic inheritance of Cnp1 positioning and quantified the rate of occasional repositioning of Cnp1 nucleosomes throughout cell generations. Together, our results reveal the plasticity and the epigenetically inheritable nature of centromeric chromatin organization.     

Beyond the Nucleosome: Nucleosome-Protein Interactions and Higher Order Chromatin Structure

The regulation of chromatin biology ultimately depends on the manipulation of its smallest subunit, the nucleosome. The proteins that bind and operate on the nucleosome do so, while their substrate is part of a polymer embedded in the dense nuclear environment. Their molecular interactions must in some way be tuned to deal with this complexity. Due to the rapid increase in the number of high-resolution structures of nucleosome-protein complexes and the increasing understanding of the cellular chromatin structure, it is starting to become clearer how chromatin factors operate in this complex environment. In this review, we analyze the current literature on the interplay between nucleosome-protein interactions and higher-order chromatin structure. We examine in what way nucleosomes-protein interactions can affect and can be affected by chromatin organization at the oligonucleosomal level. In addition, we review the characteristics of nucleosome-protein interactions that can cause phase separation of chromatin. Throughout, we hope to illustrate the exciting challenges in characterizing nucleosome-protein interactions beyond the nucleosome.     

蓝氏贾第虫核小体的初步研究

最原始的真核生物蓝氏贾第虫虽已有五种组蛋白,但以微球菌核酸酶水解却得不到规则的DNA片段。鉴于贾第虫的特殊进化地位,探讨其是否具有核小体结构,对于研究核小体的起源和进化具有重要的意义。采用改进的染色质铺展技术制备核小体并进行透射电镜观察。结果表明蓝氏贾第虫已有了直径约l0nm的核小体结构。作者认为核小体的形成可追溯到真核生物形成的初期甚至更早,核小体的完善则在真核生物形成之后。     

基于位点保守性参量和位置权重矩阵预测酵母转录因子结合位点

目的:改进转录因子结合位点的理论预测方法。方法:构建转录因子结合位点位置权重矩阵,以转录因子结合位点每一位置的碱基保守性指数Mi为参量,利用位置权重打分函数算法(PWMSA)对酵母五种转录因子结合位点进行预测。结果:利用self-consistency和cross-validation两种方法对此算法进行检验,均获得了较高的预测成功率,结果表明5种转录因子结合位点的预测成功率均超过80%。结论:与已有的三种预测转录因子结合位点的软件进行比较,PWMSA算法明显优于其他三种算法,核苷酸水平上的关联系数和结合位点水平上的关联系数分别提高了0.25和0.22。     

Histone H3 Lysine 14 (H3K14) Acetylation Facilitates DNA Repair in a Positioned Nucleosome by Stabilizing the Binding of the Chromatin Remodeler RSC (Remodels Structure of Chromatin)

Histone H3 acetylation is induced by UV damage in yeast and may play an important role in regulating the repair of UV photolesions in nucleosome-loaded genomic loci. However, it remains elusive how H3 acetylation facilitates repair. We generated a strongly positioned nucleosome containing homogeneously acetylated H3 at Lys-14 (H3K14ac) and investigated possible mechanisms by which H3K14 acetylation modulates repair. We show that H3K14ac does not alter nucleosome unfolding dynamics or enhance the repair of UV-induced cyclobutane pyrimidine dimers by UV photolyase. Importantly, however, nucleosomes with H3K14ac have a higher affinity for purified chromatin remodeling complex RSC (Remodels the Structure of Chromatin) and show greater cyclobutane pyrimidine dimer repair compared with unacetylated nucleosomes. Our study indicates that, by anchoring RSC, H3K14 acetylation plays an important role in the unfolding of strongly positioned nucleosomes during repair of UV damage.     

Substituent Effect on the Absolute Stereochemistry of the Asymmetric Reduction of Fluorine-Containing β-Diketones by Bakers' Yeast

Fluorine-containing β-diketones (la-d) were reduced by free bakers' yeast (FBY) and immobilized baker yeast (IMBY) in water, to give optically-active fluorinated β-hydroxyketones (2a-d). It was found that the reaction is highly regioselective, and that the stereochemistry of the reduction is controlled by the R substituent.     

Role of Direct Interactions between the Histone H4 Tail and the H2A Core in Long Range Nucleosome Contacts

In eukaryotic nuclei the majority of genomic DNA is believed to exist in higher order chromatin structures. Nonetheless, the nature of direct, long range nucleosome interactions that contribute to these structures is poorly understood. To determine whether these interactions are directly mediated by contacts between the histone H4 amino-terminal tail and the acidic patch of the H2A/H2B interface, as previously demonstrated for short range nucleosomal interactions, we have characterized the extent and effect of disulfide cross-linking between residues in histones contained in different strands of nucleosomal arrays. We show that in 208-12 5 S rDNA and 601-177-12 nucleosomal array systems, direct interactions between histones H4-V21C and H2A-E64C can be captured. This interaction depends on the extent of initial cross-strand association but does not require these specific residues, because interactions with residues flanking H4-V21C can also be captured. Additionally, we find that trapping H2A-H4 intra-array interactions antagonizes the ability of these arrays to undergo intermolecular self-association.     

Structural Mechanics of DNA Wrapping in the Nucleosome

Experimental X-ray crystal structures and a database of calculated structural parameters of DNA octamers were used in combination to analyse the mechanics of DNA bending in the nucleosome core complex. The 1kx5 X-ray crystal structure of the nucleosome core complex was used to determine the relationship between local structure at the base-step level and the global superhelical conformation observed for nucleosome-bound DNA. The superhelix is characterised by a large curvature (597°) in one plane and very little curvature (10°) in the orthogonal plane. Analysis of the curvature at the level of 10-step segments shows that there is a uniform curvature of 30° per helical turn throughout most of the structure but that there are two sharper kinks of 50° at ± 2 helical turns from the central dyad base pair. The curvature is due almost entirely to the base-step parameter roll. There are large periodic variations in roll, which are in phase with the helical twist and account for 500° of the total curvature. Although variations in the other base-step parameters perturb the local path of the DNA, they make minimal contributions to the total curvature. This implies that DNA bending in the nucleosome is achieved using the roll-slide-twist degree of freedom previously identified as the major degree of freedom in naked DNA oligomers. The energetics of bending into a nucleosome-bound conformation were therefore analysed using a database of structural parameters that we have previously developed for naked DNA oligomers. The minimum energy roll, the roll flexibility force constant and the maximum and minimum accessible roll values were obtained for each base step in the relevant octanucleotide context to account for the effects of conformational coupling that vary with sequence context. The distribution of base-step roll values and corresponding strain energy required to bend DNA into the nucleosome-bound conformation defined by the 1kx5 structure were obtained by applying a constant bending moment. When a single bending moment was applied to the entire sequence, the local details of the calculated structure did not match the experiment. However, when local 10-step bending moments were applied separately, the calculated structure showed excellent agreement with experiment. This implies that the protein applies variable bending forces along the DNA to maintain the superhelical path required for nucleosome wrapping. In particular, the 50° kinks are constraints imposed by the protein rather than a feature of the 1kx5 DNA sequence. The kinks coincide with a relatively flexible region of the sequence, and this is probably a prerequisite for high-affinity nucleosome binding, but the bending strain energy is significantly higher at these points than for the rest of the sequence. In the most rigid regions of the sequence, a higher strain energy is also required to achieve the standard 30° curvature per helical turn. We conclude that matching of the DNA sequence to the local roll periodicity required to achieve bending, together with the increased flexibility required at the kinks, determines the sequence selectivity of DNA wrapping in the nucleosome.     

Absolute stereochemistry of 2-hydroxyglutaric acid present in methanopterin

The absolute stereochemistry of the 2-hydroxyglutaric acid present in methanopterin and its α-glutamyl derivative, sarcinapterin, has been determined to be (S). Thus, the absolute stereochemistry of the 2-hydroxyglutaric acid present in the side chain of these modified folates is the same as the L -glutamic acid present in the side chain of folates. Chirality 9:678–680, 1997. © 1997 Wiley-Liss, Inc.     

植物抗逆蛋白（LEA3）22-氨基酸耐盐结构域在酵母细胞中的鉴定

大豆PM2蛋白属LEA3 （late embryogenesis abundant 3）蛋白。本文构建了编码PM2全长及含22-氨基酸结构域多肽（PM2、PM2A、PM2B和 PM2C）的酵母表达载体。转化酵母得到四种重组菌。SDS-PAGE电泳和ESI-MS/MS或MALDI-TOF/TOF质谱结果表明,重组菌可表达目标多肽。测定对照菌及重组菌在无胁迫、高盐（1.5 mol·L^-1 NaCl）和高渗透（2 mol·L^-1山梨糖）下的生长曲线。结果表明,在高盐胁迫下,四种重组菌胁迫后的恢复明显好于对照菌。多肽对高盐耐受能力的大小为:PM2C>PM2B≈PM2A≈PM2。证明大豆PM2蛋白的表达可提高酵母耐盐性,且22-氨基酸基序为PM2蛋白的耐盐结构域。结合前文在大肠杆菌中的结果,为“LEA蛋白可以类似机制参与原核和真核生物耐盐保护作用”的假说提供实验支持。然而,在高山梨糖胁迫下,对照菌和酵母重组菌的生长情况无明显差异。     

MeCP2 preferentially binds to methylated linker DNA in the absence of the terminal tail of histone H3 and independently of histone acetylation

Methyl CpG binding protein 2 (MeCP2) is a basic protein that contains a DNA methyl binding domain. The mechanism by which the highly positive charge of MeCP2 and its ability to bind methylated DNA contribute to the specificity of its binding to chromatin has long remained elusive. In this paper, we show that MeCP2 binds to nucleosomes in a very similar way to linker histones both in vitro and in vivo. However, its binding specificity strongly depends on DNA methylation. We also observed that as with linker histones, this binding is independent of the core histone H3 N-terminal tail and is not affected by histone acetylation.     

Predicting the Binding Patterns of Hub Proteins: A Study Using Yeast Protein Interaction Networks

Background

Protein-protein interactions are critical to elucidating the role played by individual proteins in important biological pathways. Of particular interest are hub proteins that can interact with large numbers of partners and often play essential roles in cellular control. Depending on the number of binding sites, protein hubs can be classified at a structural level as singlish-interface hubs (SIH) with one or two binding sites, or multiple-interface hubs (MIH) with three or more binding sites. In terms of kinetics, hub proteins can be classified as date hubs (i.e., interact with different partners at different times or locations) or party hubs (i.e., simultaneously interact with multiple partners).

Methodology

Our approach works in 3 phases: Phase I classifies if a protein is likely to bind with another protein. Phase II determines if a protein-binding (PB) protein is a hub. Phase III classifies PB proteins as singlish-interface versus multiple-interface hubs and date versus party hubs. At each stage, we use sequence-based predictors trained using several standard machine learning techniques.

Conclusions

Our method is able to predict whether a protein is a protein-binding protein with an accuracy of 94% and a correlation coefficient of 0.87; identify hubs from non-hubs with 100% accuracy for 30% of the data; distinguish date hubs/party hubs with 69% accuracy and area under ROC curve of 0.68; and SIH/MIH with 89% accuracy and area under ROC curve of 0.84. Because our method is based on sequence information alone, it can be used even in settings where reliable protein-protein interaction data or structures of protein-protein complexes are unavailable to obtain useful insights into the functional and evolutionary characteristics of proteins and their interactions.

Availability

We provide a web server for our three-phase approach: http://hybsvm.gdcb.iastate.edu.