Pasquale De Santis

A Theoretical Model for the Prediction of Sequence-Dependent Nucleosome Thermodynamic Stability

A theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence is advanced. The model is based on a statistical mechanical approach, which allows the calculation of the canonical ensemble free energy involved in the competitive nucleosome reconstitution. It is based on the hypothesis that nucleosome stability mainly depends on the bending and twisting elastic energy to transform the DNA intrinsic superstructure into the nucleosomal structure. The ensemble average free energy is calculated starting from the intrinsic curvature, obtained by integrating the dinucleotide step deviations from the canonical B-DNA and expressed in terms of a Fourier series, in the framework of first-order elasticity. The sequence-dependent DNA flexibility is evaluated from the differential double helix thermodynamic stability. A large number of free-energy experimental data, obtained in different laboratories by competitive nucleosome reconstitution assays, are successfully compared to the theoretical results. They support the hypothesis that the stacking energies are the major factor in DNA rigidity and could be a measure of DNA stiffness. A dual role of DNA intrinsic curvature and flexibility emerges in the determination of nucleosome stability. The difference between the experimental and theoretical (elastic) nucleosome-reconstitution free energy for the whole pool of investigated DNAs suggests a significant role for the curvature-dependent DNA hydration and counterion interactions, which appear to destabilize nucleosomes in highly curved DNAs. This model represents an attempt to clarify the main features of the nucleosome thermodynamic stability in terms of physical-chemical parameters and suggests that in molecular systems with a large degree of complexity, the average molecular properties dominate over the local features, as in a statistical ensemble.

Mapping the intrinsic curvature and flexibility along the DNA chain

The energy of DNA deformation plays a crucial and active role in its packaging and its function in the cell. Considerable effort has gone into developing methodologies capable of evaluating the local sequence-directed curvature and flexibility of a DNA chain. These studies thus far have focused on DNA constructs expressly tailored either with anomalous flexibility or curvature tracts. Here we demonstrate that these two structural properties can be mapped also along the chain of a “natural” DNA with any sequence on the basis of its scanning force microscope (SFM) images. To know the orientation of the sequence of the investigated DNA molecules in their SFM images, we prepared a palindromic dimer of the long DNA molecule under study. The palindromic symmetry also acted as an internal gauge of the statistical significance of the analysis carried out on the SFM images of the dimer molecules. It was found that although the curvature modulus is not efficient in separating static and dynamic contributions to the curvature of the population of molecules, the curvature taken with its direction (its sign in two dimensions) permits the direct separation of the intrinsic curvature from the flexibility contributions. The sequence-dependent flexibility seems to vary monotonically with the chains intrinsic curvature; the chain rigidity was found to modulate as its local thermodynamic stability and does not correlate with the dinucleotide chain rigidities evaluation made from x-ray data by other authors.

Sequence-Dependent DNA Curvature and Flexibility from Scanning Force Microscopy Images

This paper reports a study of the sequence-dependent DNA curvature and flexibility based on scanning force microscopy (SFM) images. We used a palindromic dimer of a 1878-bp pBR322 fragment and collected a large pool of SFM images. The curvature of each imaged chain was measured in modulus and direction. It was found that the ensemble curvature modulus does not allow the separation of static and dynamic contributions to the curvature, whereas the curvature, when its direction in the two dimensions is taken into account, permits the direct separation of the intrinsic curvature contributions static and dynamic contributions. The palindromic symmetry also acted as an internal gauge of the validity of the SFM images statistical analysis. DNA static curvature resulted in good agreement with the predicted sequence-dependent intrinsic curvature. Furthermore, DNA sequence-dependent flexibility was found to correlate with the occurrence of A.T-rich dinucleotide steps along the chain and, in general, with the normalized basepair stacking energy distribution.

The main role of the sequence-dependent DNA elasticity in determining the free energy of nucleosome formation on telomeric DNAs.

Using a competitive reconstitution assay, we measured the free energy spent in nucleosome formation of eight telomeric DNAs, differing in sequence and/or in length. The obtained values are in satisfactorily good agreement with those derived from a theoretical model that allows the calculation of the free energy of nucleosome formation on the basis of sequence-dependent DNA elasticity, using a statistical thermodynamic approach. Both theoretical and experimental evaluations show that telomeres are characterized by the highest free energies of nucleosome formation among all the DNA sequences so far studied. The free energy of nucleosome formation varies according to the different telomeric sequences and the length of the fragments. Theoretical analysis and experimental mapping by lambda exonuclease show that telomeric nucleosomes occupy multiple positions spaced every telomeric repeat. Sequence-dependent DNA elasticity appears as the main determinant of the stability of telomeric nucleosomes and their multiple translational positioning.

Prediction of nucleosome positioning in genomes: limits and perspectives of physical and bioinformatic approaches.

Abstract Nucleosomes, the fundamental repeating subunits of all eukaryotic chromatin, are responsible for packaging DNA into chromosomes inside the cell nucleus and controlling gene expression. While it has been well established that nucleosomes exhibit higher affinity for select DNA sequences, until recently it was unclear whether such preferences exerted a significant, genome-wide effect on nucleosome positioning in vivo. For this reason, an increasing interest is arising on a wide-ranging series of experimental and computational analyses capable of predicting the nucleosome positioning along genomes. Toward this goal, we propose a theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence. Based on a statistical mechanical approach, the model allows the calculation of the sequence-dependent canonical ensemble free energy involved in nucleosome formation. The theoretical free energies were evaluated for 90 single nucleosome DNA tracts and successfully compared with those obtained with nucleosome competitive reconstitution. These results, obtained for single nucleosomes, could in principle allow the calculation of the intrinsic affinity of nucleosomes along DNA sequences virtually opening the possibility of predicting the nucleosome positioning along genomes on physical basis. The theoretical nucleosome distribution was compared and validated with that of yeast and human genome experimentally determined. The results interpret on a physical basis the experimental nucleosome positioning and are comparable with those obtained adopting models based on the identification of some recurrent sequence features obtained from the statistical analysis of a very large pool of nucleosomal DNA sequences provided by the positioning maps of genomes.

Recognition of the DNA sequence by an inorganic crystal surface

The sequence-dependent curvature is generally recognized as an important and biologically relevant property of DNA because it is involved in the formation and stability of association complexes with proteins. When a DNA tract, intrinsically curved for the periodical recurrence on the same strand of A-tracts phased with the B-DNA periodicity, is deposited on a flat surface, it exposes to that surface either a T- or an A-rich face. The surface of a freshly cleaved mica crystal recognizes those two faces and preferentially interacts with the former one. Statistical analysis of scanning force microscopy (SFM) images provides evidence of this recognition between an inorganic crystal surface and nanoscale structures of double-stranded DNA. This finding could open the way toward the use of the sequence-dependent adhesion to specific crystal faces for nanotechnological purposes.

Dual role of sequence-dependent DNA curvature in nucleosome stability: the critical test of highly bent Crithidia fasciculata DNA tract

In spite of the knowledge of the nucleosome molecular structure, the role of DNA intrinsic curvature in determining nucleosome stabilization is still an open question. In this paper, we describe a general model that allows the prediction of the nucleosome stability, tested on 83 different DNA sequences, in surprising good agreement with the experimental data, carried out in ours as well as in many other laboratories. The model is based on the dual role of DNA curvature in nucleosome thermodynamic stabilization. A critical test is the evaluation of the nucleosome free energy relative to a Crithidia fasciculata kinetoplast DNA fragment, which represents the most curved DNA found so far in biological systems and, therefore, is generally believed to form a highly stable nucleosome.

Geometrical, conformational and topological restraints in regular nucleosome compaction in chromatin.

The folding of the nucleosome array into a chromatin fiber modulates DNA accessibility and is therefore an important factor for the control of gene expression. The statistical analysis of the nucleosome repeat length in chromatin fibers reveals the presence of a ten-fold periodicity suggesting the existence of orientational constraints of the nucleosome units that provide the geometrical conditions of helical conformations. Recently, the elucidation of the x-ray crystal structure of a nucleosome tetramer array and the interpretation of electron microscopy images of reconstituted nucleosome arrays suggested two different architectures of the chromatin fiber. We approached the problem by integrating the experimental findings with geometrical, conformational and topological restraints, under the hypothesis of the minimum distortion of the nucleosome and linker DNA structures. We show that the excluded volume at linker crossing and the torsional energy limit the possible close packing of the nucleosomes in the chromatin fiber. In particular, the torsional energy of the chromatin fiber appears crucial in determining the kind of nucleosome packing for short nucleosome repeat lengths as in telomeres and yeast chromatin.

Sequence-Dependent DNA Dynamics by Scanning Force Microscopy Time-Resolved Imaging

Scanning force microscopy was used to study in fluid the conformational fluctuations of two double-stranded DNA molecules resulting from differently cut pBR322 circular DNAs. A new approach was conceived to monitor the thermodynamic equilibrium of the chain dynamics on different scale lengths. This method made it possible to demonstrate that both the observed DNA molecules were allowed to equilibrate only on their local small-scale dynamics during the time of the experiment. This capability of monitoring the length scale and the time scale of the equilibration processes in the dynamics of a DNA chain is relevant to give an insight in the thermodynamics of the DNA binding with proteins and synthetic ligands. It was also shown that the small-scale equilibration of the DNA chain during surface-restricted dynamics is enough to allow a valid measurement of the local sequence-dependent curvature.

A statistical thermodynamic approach for predicting the sequence-dependent nucleosome positioning along genomes

Nucleosomes are the fundamental repeating unit of chromatin and constitute the structural building blocks of the eukaryotic genome. The distribution of nucleosomes along the genome is a significant aspect of chromatin structure and influences gene regulation through modulation of DNA accessibility. For this reason, an increasing interest is arising in models capable of predicting the nucleosome positioning along genomes. Toward this goal, we propose a theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence. The model, based on a statistical mechanical approach allows the calculation of the canonical ensemble free energy involved in nucleosome formation. The theoretical free energies were evaluated for about one hundred nucleosome DNA tracts and successfully compared with those obtained in different laboratories with nucleosome competitive reconstitution (correlation coefficient equal to 0.92). We extended these results to the nucleosome positioning along genomes. To test our model, the theoretical nucleosome distribution was compared with that of yeast genome experimentally determined. The results are comparable with those obtained by different authors adopting models based on identifying some recurrent sequence features obtained from the statistical analysis of a very large pool of nucleosomal DNA sequences provided by the positioning maps of genomes.