Victor B. Zhurkin

Modeling DNA deformations.

Recent developments have been made in modeling double-helical DNA at four levels of three-dimensional structure: the all-atom level, whereby an oligonucleotide duplex is surrounded by a shroud of solvent molecules; the base-pair level, with explicit backbone atoms; the mesoscopic level, that is, a few hundred base pairs, with the local duplex conformation described by knowledge-based harmonic energy functions; and the scale of several thousand nucleotides, with the duplex described as an ideal elastic rod. Predictions of the sequence-dependent bending and twisting of the double helix, as well as solvent- and force-induced B-->A and over-stretching conformational transitions, are compared with experimental data. These subtle conformational changes are critical to the functioning of the double helix, including its packaging in the close confines of the cell, the mutual fit of DNA and protein in nucleoprotein complexes, and the effective recognition of base pairs in recombination and transcription.

Overall structure and sugar dynamics of a DNA dodecamer from homo- and heteronuclear dipolar couplings and 31P chemical shift anisotropy.

The solution structure of d(CGCGAATTCGCG)2 has been determined on the basis of an exceptionally large set of residual dipolar couplings. In addition to the heteronuclear 13C-1H and 15N-1H and qualitative homonuclear 1H-1H dipolar couplings, previously measured in bicelle medium, more than 300 quantitative 1H-1H and 22 31P-1H dipolar restraints were obtained in liquid crystalline Pf1 medium, and 22 31P chemical shift anisotropy restraints. High quality DNA structures can be obtained solely on the basis of these new restraints, and these structures are in close agreement with those calculated previously on the basis of 13C-1H and 15N-1H dipolar couplings. In the newly calculated structures, 31P-1H dipolar and 3JsubH3′Psub couplings and 31P CSA data restrain the phosphodiester backbone torsion angles. The final structure represents a quite regular B-form helix with a modest bending of ∼10°, which is essentially independent of whether or not electrostatic terms are used in the calculation. Combined, the number of homo- and heteronuclear dipolar couplings significantly exceeds the number of degrees of freedom in the system. Results indicate that the dipolar coupling data cannot be fit by a single structure, but are compatible with the presence of rapid equilibria between C2′-endo and C3′-endo deoxyribose puckers (sugar switching). The C2′-H2′/H2′′ dipolar couplings in B-form DNA are particularly sensitive to sugar pucker and yield the largest discrepancies when fit to a single structure. To resolve these discrepancies, we suggest a simplified dipolar coupling analysis that yields N/S equilibria for the ribose sugar puckers, which are in good agreement with previous analyses of NMR JHH couplings, with a population of the minor C3′-endo form higher for pyrimidines than for purines.

Gal repressosome contains an antiparallel DNA loop.

Gal repressosome assembly and repression of the gal operon in Escherichia coli occurs when two dimeric GalR proteins and the histone-like HU protein bind to cognate sites causing DNA looping. Structure-based genetic analysis defined the GalR surfaces interacting to form a stacked, V-shaped, tetrameric structure. Stereochemical models of the four possible DNA loops compatible with the GalR tetramer configuration were constructed using the sequence-dependent structural parameters of the interoperator DNA and conformation changes caused by GalR and asymmetric HU binding. Evaluation of their DNA elastic energies gave unambiguous preference to a loop structure in which the two gal operators adopt an antiparallel orientation causing undertwisting of DNA.

A-tract clusters may facilitate DNA packaging in bacterial nucleoid

Molecular mechanisms of bacterial chromosome packaging are still unclear, as bacteria lack nucleosomes or other apparent basic elements of DNA compaction. Among the factors facilitating DNA condensation may be a propensity of the DNA molecule for folding due to its intrinsic curvature. As suggested previously, the sequence correlations in genome reflect such a propensity [Trifonov and Sussman (1980) Proc. Natl Acad. Sci. USA, 77, 3816–3820]. To further elaborate this concept, we analyzed positioning of A-tracts (the sequence motifs introducing the most pronounced DNA curvature) in the Escherichia coli genome. First, we observed that the A-tracts are over-represented and distributed ‘quasi-regularly’ throughout the genome, including both the coding and intergenic sequences. Second, there is a 10–12 bp periodicity in the A-tract positioning indicating that the A-tracts are phased with respect to the DNA helical repeat. Third, the phased A-tracts are organized in ∼100 bp long clusters. The latter feature was revealed with the help of a novel approach based on the Fourier series expansion of the A-tract distance autocorrelation function. Since the A-tracts introduce local bends of the DNA duplex and these bends accumulate when properly phased, the observed clusters would facilitate DNA looping. Also, such clusters may serve as binding sites for the nucleoid-associated proteins that have affinities for curved DNA (such as HU, H-NS, Hfq and CbpA). Therefore, we suggest that the ∼100 bp long clusters of the phased A-tracts constitute the ‘structural code’ for DNA compaction by providing the long-range intrinsic curvature and increasing stability of the DNA complexes with architectural proteins.

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.

Distinctive sequence patterns in metazoan and yeast nucleosomes: implications for linker histone binding to AT-rich and methylated DNA.

Linker histones (LHs) bind to the DNA entry/exit points of nucleosomes and demonstrate preference for AT-rich DNA, although the recognized sequence patterns remain unknown. These patterns are expected to be more pronounced in metazoan nucleosomes with abundant LHs, compared to yeast nucleosomes with few LHs. To test this hypothesis, we compared the nucleosome core particle (NCP) sequences from chicken, Drosophila and yeast, extending them by the flanking sequences extracted from the genomes. We found that the known ∼10-bp periodic oscillation of AT-rich elements goes beyond the ends of yeast nucleosomes, but is distorted in metazoan sequences where the ‘out-of-phase’ AT-peaks appear at the NCP ends. The observed difference is likely to be associated with sequence-specific LH binding. We therefore propose a new structural model for LH binding to metazoan nucleosomes, postulating that the highly conserved nonpolar ‘wing’ region of the LH globular domain (tetrapeptide GVGA) recognizes AT-rich fragments through hydrophobic interactions with the thymine methyl groups. These interactions lead to DNA bending at the NCP ends and formation of a ‘stem-like’ structure. The same mechanism accounts for the high affinity of LH to methylated DNA—a feature critical for stabilization of the higher-order structure of chromatin and for repression of transcription.

Mitogen-activated protein kinase phosphatase 2: a novel transcription target of p53 in apoptosis.

The p53 tumor suppressor plays critical roles in diverse cellular responses such as cell cycle arrest, senescence, and apoptosis through transcriptional control of its target genes. Identification and characterization of new p53 target genes will advance our understanding of how p53 exerts its multiple regulatory functions. In this article, we show that mitogen-activated protein kinase phosphatase 2 (MKP2) is a novel transcription target of p53 in mediating apoptosis. Moreover, we identify a 10-bp perfect palindrome motif (CTGGCGCCAG) in the MKP2 promoter as a new binding site for p53 to activate the MKP2 gene. This GC-rich palindrome is completely different from the consensus p53 binding sequence. Induction of MKP2 is highly responsive to oxidative stress in a p53-dependent manner. Interestingly, the p53-dependent induction of MKP2 is prominent only in the cellular response to stimuli leading to apoptosis but not to cell cycle arrest. In response to oxidative stress, MKP2 is not only required for p53-mediated apoptosis, but ectopic MKP2 expression can also enhance apoptotic responses even independent of p53. These data suggest that p53 regulates distinct genes via different binding mechanisms and that MKP2 is an essential target of p53 in signaling apoptosis.

Sequence-dependent Kink-and-Slide Deformations of Nucleosomal DNA Facilitated by Histone Arginines Bound in the Minor Groove

Abstract In addition to bending and twisting deformabilities, the lateral displacements of the DNA axis (Kink-and-Slide) play an important role in DNA wrapping around the histone core (M. Y. Tolstorukov, A. V. Colasanti, D. M. McCandlish, W. K. Olson, V. B. Zhurkin, J. Mol. Biol. 371, 725-738 (2007)). Here, we show that these Kink-and-Slide deformations are likely to be stabilized by the arginine residues of histones interacting with the minor groove of DNA. The arginines are positioned asymmetrically in the minor groove, being closer to one strand. The asymmetric arginine-DNA interactions facilitate lateral displacement of base pairs across the DNA grooves, thus leading to a stepwise accumulation of the superhelical pitch of nucleosomal DNA. To understand the sequence dependence of such Kink-and-Slide deformations, we performed all-atom calculations of DNA hexamers with the YR and RY steps in the center. We found that when the unrestrained DNA deformations are allowed, the YR steps tend to bend into the major groove, and RY steps bend into the minor groove. However, when the nucleosomal Kink-and-Slide deformation is considered, the YR steps prove to be more favorable for bending into the minor groove. Overall, the Kink-and-Slide deformation energy of DNA increases in the order TA < CA < CG < GC < AC < AT. We propose a simple stereochemical model accounting for this sequence dependence. Our results agree with experimental data indicating that the TA step most frequently occurs in the minor-groove kink positions in the most stable nucleosomes. Our computations demonstrate that the Kink-and-Slide distortion is accompanied by the BI to BII transition. This fact, together with irregularities in the two-dimensional (Roll, Slide) energy contour maps, suggest that the Kink-and-Slide deformations represent a nonharmonic behavior of the duplex. This explains the difference between the two estimates of the DNA deformation energy in nucleosome - the earlier one made using knowledge-based elastic energy functions, and the current one based on all-atom calculations. Our findings are useful for refining the score functions for the prediction of nucleosome positioning. In addition, the reverse bending behavior of the YR and RY steps revealed under the Kink-and-Slide constraint is important for understanding the molecular mechanisms of binding transcription factors (such as p53) to DNA exposed on the surface of nucleosome.

nuScore: a web-interface for nucleosome positioning predictions.

UNLABELLED Sequence-directed mapping of nucleosome positions is of major biological interest. Here, we present a web-interface for estimation of the affinity of the histone core to DNA and prediction of nucleosome arrangement on a given sequence. Our approach is based on assessment of the energy cost of imposing the deformations required to wrap DNA around the histone surface. The interface allows the user to specify a number of options such as selecting from several structural templates for threading calculations and adding random sequences to the analysis. AVAILABILITY The nuScore interface is freely available for use at http://compbio.med.harvard.edu/nuScore. CONTACT peter_park@harvard.edu; tolstorukov@gmail.com SUPPLEMENTARY INFORMATION The site contains user manual, description of the methodology and examples.

Working the kinks out of nucleosomal DNA.

Condensation of DNA in the nucleosome takes advantage of its double-helical architecture. The DNA deforms at sites where the base pairs face the histone octamer. The largest so-called kink-and-slide deformations occur in the vicinity of arginines that penetrate the minor groove. Nucleosome structures formed from the 601 positioning sequence differ subtly from those incorporating an AT-rich human α-satellite DNA. Restraints imposed by the histone arginines on the displacement of base pairs can modulate the sequence-dependent deformability of DNA and potentially contribute to the unique features of the different nucleosomes. Steric barriers mimicking constraints found in the nucleosome induce the simulated large-scale rearrangement of canonical B DNA to kink-and-slide states. The pathway to these states shows nonharmonic behavior consistent with bending profiles inferred from AFM measurements.