Dataset Information

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA.

ABSTRACT: We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.

SUBMITTER: Chakraborty D

PROVIDER: S-EPMC6423546 | biostudies-literature | 2018 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA.

Chakraborty Debayan D Hori Naoto N Thirumalai D D

Journal of chemical theory and computation 20180626 7

We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inv ...[more]

PMID: 29870236

Similar Datasets

Project description:BackgroundTelomeres are protective cap structures at the ends of the linear eukaryotic chromosomes, which provide stability to the genome by shielding from degradation and chromosome fusions. The cap consists of telomere-specific proteins binding to the respective single- and double-stranded parts of the telomeric sequence. In addition to the nucleation of the chromatin structure the telomere-binding proteins are involved in the regulation of the telomere length. However, the telomeric sequences are highly diverged among yeast species. During the evolution this high rate of divergency presents a challenge for the sequence recognition of the telomere-binding proteins.ResultsWe found that the Saccharomyces castellii protein Rap1, a negative regulator of telomere length, binds a 12-mer minimal binding site (MBS) within the double-stranded telomeric DNA. The sequence specificity is dependent on the interaction with two 5 nucleotide motifs, having a 6 nucleotide centre-to-centre spacing. The isolated DNA-binding domain binds the same MBS and retains the same motif binding characteristics as the full-length Rap1 protein. However, it shows some deviations in the degree of sequence-specific dependence in some nucleotide positions. Intriguingly, the positions of most importance for the sequence-specific binding of the full-length Rap1 protein coincide with 3 of the 4 nucleotides utilized by the 3' overhang binding protein Cdc13. These nucleotides are very well conserved within the otherwise highly divergent telomeric sequences of yeasts.ConclusionsRap1 and Cdc13 are two very distinct types of DNA-binding proteins with highly separate functions. They interact with the double-stranded vs. the single-stranded telomeric DNA via significantly different types of DNA-binding domain structures. However, we show that they are dependent on coinciding nucleotide positions for their sequence-specific binding to telomeric sequences. Thus, we conclude that during the molecular evolution they act together to preserve a core sequence of the telomeric DNA.

Dataset Information

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA.

Publications

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets