Intrinsic DNA Curvature Research Articles

Bending of Canonical and G/T Mismatched DNAs.

Mismatched base pairs alter the flexibility and intrinsic curvature of DNA. The role of such DNA features is not fully understood in the mismatch repair pathway. MutS/DNA complexes exhibit DNA bending, PHE intercalation, and changes of base-pair parameters near the mismatch. Recently, we have shown that base-pair opening in the absence of MutS can discriminate mismatches from canonical base pairs better than DNA bending. However, DNA bending in the absence of MutS was found to be rather challenging to describe correctly. Here, we present a computational study on the DNA bending of canonical and G/T mismatched DNAs. Five types of geometric parameters covering template-based bending toward the experimental DNA structure, global, and local geometry parameters were employed in biased molecular dynamics in the absence of MutS. None of these parameters showed higher discrimination than the base-pair opening. Only roll could induce a sharply localized bending of DNA as observed in the experimental MutS/DNA structure. Further, we demonstrated that the intercalation of benzene mimicking PHE decreases the energetic cost of DNA bending without any effect on mismatch discrimination.

Reduced intrinsic DNA curvature leads to increased mutation rate

BackgroundMutation rates vary across the genome. Many trans factors that influence mutation rates have been identified, as have specific sequence motifs at the 1–7-bp scale, but cis elements remain poorly characterized. The lack of understanding regarding why different sequences have different mutation rates hampers our ability to identify positive selection in evolution and to identify driver mutations in tumorigenesis.ResultsHere, we use a combination of synthetic genes and sequences of thousands of isolated yeast colonies to show that intrinsic DNA curvature is a major cis determinant of mutation rate. Mutation rate negatively correlates with DNA curvature within genes, and a 10% decrease in curvature results in a 70% increase in mutation rate. Consistently, both yeast and humans accumulate mutations in regions with small curvature. We further show that this effect is due to differences in the intrinsic mutation rate, likely due to differences in mutagen sensitivity and not due to differences in the local activity of DNA repair.ConclusionsOur study establishes a framework for understanding the cis properties of DNA sequence in modulating the local mutation rate and identifies a novel causal source of non-uniform mutation rates across the genome.

Intrinsic DNA curvature in trypanosomes

BackgroundTrypanosoma cruzi and Trypanosoma brucei are protozoan parasites causing Chagas disease and African sleeping sickness, displaying unique features of cellular and molecular biology. Remarkably, no canonical signals for RNA polymerase II promoters, which drive protein coding genes transcription, have been identified so far. The secondary structure of DNA has long been recognized as a signal in biological processes and more recently, its involvement in transcription initiation in Leishmania was proposed. In order to study whether this feature is conserved in trypanosomatids, we undertook a genome wide search for intrinsic DNA curvature in T. cruzi and T. brucei.ResultsUsing a region integrated intrinsic curvature (RIIC) scoring that we previously developed, a non-random distribution of sequence-dependent curvature was observed. High RIIC scores were found to be significantly correlated with transcription start sites in T. cruzi, which have been mapped in divergent switch regions, whereas in T. brucei, the high RIIC scores correlated with sites that have been involved not only in RNA polymerase II initiation but also in termination. In addition, we observed regions with high RIIC score presenting in-phase tracts of Adenines, in the subtelomeric regions of the T. brucei chromosomes that harbor the variable surface glycoproteins genes.ConclusionsIn both T. cruzi and T. brucei genomes, a link between DNA conformational signals and gene expression was found. High sequence dependent curvature is associated with transcriptional regulation regions. High intrinsic curvature also occurs at the T. brucei chromosome subtelomeric regions where the recombination processes involved in the evasion of the immune host system take place. These findings underscore the relevance of indirect DNA readout in these ancient eukaryotes.

Modulation of extracytoplasmic function (ECF) sigma factor promoter selectivity by spacer region sequence.

The ability of bacteria to adapt to stress depends on the conditional expression of specific sets of genes. Bacillus subtilis encodes seven extracytoplasmic function (ECF) sigma (σ) factors that regulate functions important for survival under conditions eliciting cell envelope stress. Of these, four have been studied in detail: σM, σW, σX and σV. These four σ factors recognize overlapping sets of promoters, although the sequences that determine this overlapping recognition are incompletely understood. A major role in promoter selectivity has been ascribed to the core −10 and −35 promoter elements. Here, we demonstrate that a homopolymeric T-tract motif, proximal to the −35 element, functions in combination with the core promoter sequences to determine selectivity for ECF sigma factors. This motif is most critical for promoter activation by σV, and contributes variably to activation by σM, σX and σW. We propose that this motif, which is a feature of the deduced promoter consensus for a subset of ECF σ factors from many species, imparts intrinsic DNA curvature to influence promoter activity. The differential effect of this region among ECF σ factors thereby provides a mechanism to modulate the nature and extent of regulon overlap.

Quaternion modeling of the helical path for analysis of the shape of the DNA molecule

The three­dimensional shape of a DNA molecule is a key property influencing its functional specificity and the nature of its molecular interactions. The characteristic shape into which a DNA molecule folds under certain conditions is a manifestation of its micromechanical and structural features, which are sequence­dependent. DNA shape­related properties can there fore be determined in a predictable manner. A number of models have been designed to describe intrinsic DNA curvature, incorporating a set of helical parameters which can be applied to operative three­dimensional reconstruction of the DNA structures. Alternatively, desired base pair parameters can be computed based on publicly available information about atomic DNA structures. Further, taking the base pairs as rigid bodies, their relative location in space can be estimated based on these parameters. Matrices are a common method to implement any rigid body transformations and are widely used in the modeling of DNA structures. Quaternions are the more straightforward and robust alternative for matrices. Unit quaternions can represent only a rotation, whereas dual quaternions combine rotation and translation into a single state. In the present guide, the algebra of unit and dual quaternions is applied for the first time to modeling of the DNA helical path, based on conformational parameters of the base pair steps. Although dual quaternions are preferable for modeling of DNA structure in detail, the use of unit quaternions is sufficient to predict the DNA trajectory and all calculations of DNA shape features. In order to analyze DNA shape and chain sta ­ tistics, and predict the micromechanical properties of DNA molecules based on coordinates of the helical path, the widely used as well as original algorithms for computing DNA curvature, radius of gyration, persistence length and phasing of DNA bends are described. Taken together, these algorithms will be useful both in the in silico analysis of relatively short DNA fragments as well as in topological mapping of whole genomes.

Conserved Curvature of RNA Polymerase I Core Promoter Beyond rRNA Genes: The Case of the Tritryps

In trypanosomatids, the RNA polymerase I (RNAPI)-dependent promoters controlling the ribosomal RNA (rRNA) genes have been well identified. Although the RNAPI transcription machinery recognizes the DNA conformation instead of the DNA sequence of promoters, no conformational study has been reported for these promoters. Here we present the in silico analysis of the intrinsic DNA curvature of the rRNA gene core promoters in Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major. We found that, in spite of the absence of sequence conservation, these promoters hold conformational properties similar to other eukaryotic rRNA promoters. Our results also indicated that the intrinsic DNA curvature pattern is conserved within the Leishmania genus and also among strains of T. cruzi and T. brucei. Furthermore, we analyzed the impact of point mutations on the intrinsic curvature and their impact on the promoter activity. Furthermore, we found that the core promoters of protein-coding genes transcribed by RNAPI in T. brucei show the same conserved conformational characteristics. Overall, our results indicate that DNA intrinsic curvature of the rRNA gene core promoters is conserved in these ancient eukaryotes and such conserved curvature might be a requirement of RNAPI machinery for transcription of not only rRNA genes but also protein-coding genes.

The Role of DNA Shape in Nucleosome Formation and Positioning

Nucleosomes provide the basic unit of compaction in eukaryotic genomes and the mechanisms that dictate their position along a DNA sequence are of central importance to epigenetics. In this work, we use molecular models of DNA and proteins to explore various aspects of nucleosome positioning. First, we show how the affinity of DNA for histones is encoded in its sequence-dependent shape, including subtle deviations from ideal straight B-DNA and local variations in minor groove width. Using high-precision simulations of the free energy of nucleosome complexes, we also demonstrate that, depending on DNA's intrinsic curvature, histone binding can be dominated by bending interactions or electrostatic interactions. Secondly, we examine the sequence-dependent mechanism by which DNA slides to its optimal position around the histone proteins. More generally, the results presented here explain how sequence, manifested as the shape of the DNA molecule, dominates molecular recognition in the problem of nucleosome formation.

DNA shape dominates sequence affinity in nucleosome formation.

Nucleosomes provide the basic unit of compaction in eukaryotic genomes, and the mechanisms that dictate their position at specific locations along a DNA sequence are of central importance to genetics. In this Letter, we employ molecular models of DNA and proteins to elucidate various aspects of nucleosome positioning. In particular, we show how DNA's histone affinity is encoded in its sequence-dependent shape, including subtle deviations from the ideal straight B-DNA form and local variations of minor groove width. By relying on high-precision simulations of the free energy of nucleosome complexes, we also demonstrate that, depending on DNA's intrinsic curvature, histone binding can be dominated by bending interactions or electrostatic interactions. More generally, the results presented here explain how sequence, manifested as the shape of the DNA molecule, dominates molecular recognition in the problem of nucleosome positioning.

A functional variant in the 5'-flanking region of the chicken serotonin transporter gene is associated with increased body weight and locomotor activity.

In this study, we identified a polymorphism in the 5'-flanking region of the chicken serotonin transporter (5-HTT) gene. Sequencing analysis revealed that in comparison with the wild-type variant (W), a deleted variant (D) is generated by deletion of four nucleotides (5'-AATT-3') and a single nucleotide change (A→T). Using a polyacrylamide gel electrophoresis system, we found that the 360-bp DNA fragment containing the W variant with the wild-type sequence 5'-AATTAATT-3' shows intrinsic DNA curvature while the 356-bp fragment containing the D variant lacking the four base pairs AATT is not curved. Quantitative real-time RT-PCR and ELISA demonstrated that the expression of 5-HTT in D/D chickens was higher than that in W/W and W/D chickens. In addition, transient transfection experiments with chloramphenicol acetyltransferase reporter gene constructs revealed increased 5-HTT promoter activity mediated by the D variant and a silencer activity of the W variant. Interestingly, females and males with D/D genotype showed significant greater increase in body weight from 6 weeks and 16 weeks of age, respectively, and higher body mass index. Moreover, we found that D/D chickens of both genders were physically more active than W/W and W/D chickens.

Investigating the Interplay between Nucleoid-Associated Proteins, DNA Curvature, and CRISPR Elements Using Comparative Genomics

Many prokaryotic and eukaryotic genomes feature a characteristic periodic signal in distribution of short runs of A or T (A-tracts) phased with the DNA helical period of ∼10–11 bp. Such periodic spacing of A-tracts has been associated with intrinsic DNA curvature. In eukaryotes, this periodicity is a major component of the nucleosome positioning signal but its physiological role in prokaryotes is not clear. One hypothesis centers on possible role of intrinsic DNA bends in nucleoid compaction. We use comparative genomics to investigate possible relationship between the A-tract periodicity and nucleoid-associated proteins in prokaryotes. We found that genomes with DNA-bridging proteins tend to exhibit stronger A-tract periodicity, presumably indicative of more prevalent intrinsic DNA curvature. A weaker relationship was detected for nucleoid-associated proteins that do not form DNA bridges. We consider these results an indication that intrinsic DNA curvature acts collaboratively with DNA-bridging proteins in maintaining the compact structure of the nucleoid, and that previously observed differences among prokaryotic genomes in terms DNA curvature-related sequence periodicity may reflect differences in nucleoid organization. We subsequently investigated the relationship between A-tract periodicity and presence of CRISPR elements and we found that genomes with CRISPR tend to have stronger A-tract periodicity. This result is consistent with our earlier hypothesis that extensive A-tract periodicity could help protect the chromosome against integration of prophages, possibly due to its role in compaction of the nucleoid.

Assessing Diversity of DNA Structure-Related Sequence Features in Prokaryotic Genomes

Prokaryotic genomes are diverse in terms of their nucleotide and oligonucleotide composition as well as presence of various sequence features that can affect physical properties of the DNA molecule. We present a survey of local sequence patterns which have a potential to promote non-canonical DNA conformations (i.e. different from standard B-DNA double helix) and interpret the results in terms of relationships with organisms' habitats, phylogenetic classifications, and other characteristics. Our present work differs from earlier similar surveys not only by investigating a wider range of sequence patterns in a large number of genomes but also by using a more realistic null model to assess significant deviations. Our results show that simple sequence repeats and Z-DNA-promoting patterns are generally suppressed in prokaryotic genomes, whereas palindromes and inverted repeats are over-represented. Representation of patterns that promote Z-DNA and intrinsic DNA curvature increases with increasing optimal growth temperature (OGT), and decreases with increasing oxygen requirement. Additionally, representations of close direct repeats, palindromes and inverted repeats exhibit clear negative trends with increasing OGT. The observed relationships with environmental characteristics, particularly OGT, suggest possible evolutionary scenarios of structural adaptation of DNA to particular environmental niches.

Evaluation of DNA bending models in their capacity to predict electrophoretic migration anomalies of satellite DNA sequences

DNA containing a sequence that generates a local curvature exhibits a pronounced retardation in electrophoretic mobility. Various theoretical models have been proposed to explain relationship between DNA structural features and migration anomaly. Here, we studied the capacity of 15 static wedge-bending models to predict electrophoretic behavior of 69 satellite monomers derived from four divergent families. All monomers exhibited retarded mobility in PAGE corresponding to retardation factors ranging 1.02-1.54. The curvature varied both within and across the groups and correlated with the number, position, and lengths of A-tracts. Two dinucleotide models provided strong correlation between gel mobility and curvature prediction; two trinucleotide models were satisfactory while remaining dinucleotide models provided intermediate results with reliable prediction for subsets of sequences only. In some cases, similarly shaped molecules exhibited relatively large differences in mobility and vice versa. Generally less accurate predictions were obtained in groups containing less homogeneous sequences possessing distinct structural features. In conclusion, relatively universal theoretical models were identified suitable for the analysis of natural sequences known to harbor relatively moderate curvature. These models could be potentially applied to genome wide studies. However, in silico predictions should be viewed in context of experimental measurement of intrinsic DNA curvature.

Genomic analysis of sequence-dependent DNA curvature in Leishmania.

Leishmania major is a flagellated protozoan parasite of medical importance. Like other members of the Trypanosomatidae family, it possesses unique mechanisms of gene expression such as constitutive polycistronic transcription of directional gene clusters, gene amplification, mRNA trans-splicing, and extensive editing of mitochondrial transcripts. The molecular signals underlying most of these processes remain under investigation. In order to investigate the role of DNA secondary structure signals in gene expression, we carried out a genome-wide in silico analysis of the intrinsic DNA curvature. The L. major genome revealed a lower frequency of high intrinsic curvature regions as well as inter- and intra- chromosomal distribution heterogeneity, when compared to prokaryotic and eukaryotic organisms. Using a novel method aimed at detecting region-integrated intrinsic curvature (RIIC), high DNA curvature was found to be associated with regions implicated in transcription initiation. Those include divergent strand-switch regions between directional gene clusters and regions linked to markers of active transcription initiation such as acetylated H3 histone, TRF4 and SNAP50. These findings suggest a role for DNA curvature in transcription initiation in Leishmania supporting the relevance of DNA secondary structures signals.

Modeling of DNA curvature: Comment on “Sequence-dependent collective properties of DNAs and their role in biological systems” by Pasquale De Santis and Anita Scipioni

Since the pioneering works of Trifonov [1,2] the idea of DNA sequence dependent intrinsic properties, which are involved or even govern many biological processes, has become widely accepted. The review of De Santis and Scipioni [3] is devoted to such superstructural DNA features, which they call “sequence-dependent collective properties of DNA”, mainly having in mind DNA curvature and bendability. The authors of this broad and detailed review do not share Trifonov’s opinion [4,5] that these are rather different properties that require separate models and study approaches. Moreover, De Santis and Scipioni even do not mention that there are researchers clearly distinguishing between intrinsic DNA curvature that is primarily the property of naked DNA and DNA bendability particularly evident in the case of sequence dependent DNA–histone octamer association. The great advantage of the review is in the consistent attitude, which I deeply appreciate and share the assurance of the authors that in this case it is the only correct scientific approach. It is in description of the physical model that explores different aspects of DNA collective properties instead of application by some means irrelevant models related to signal processing, machine learning and data mining. (In the text of the review these approaches are called “the bioinformatic approaches”.) Because of the space limitations I would like to comment on the DNA curvature model, exclusively. De Santis and Scipioni use the term curvature from the classical mathematical point of view. Though, we want to measure DNA curvature associated with a certain base pair. I do not see any sense in usage of the regular mathematical definition of the trajectory line curvature, even De Santis and Scipioni appeal to the Landau and Lifshitz [6] authorities. The DNA ‘axis’, i.e. the broken line passing through the base pair centers, or DNA strand, i.e. the broken line connecting the phosphates, cannot strictly be described in terms of curvature, unless represented as smooth curves. DNA molecule may be presented as a line in space only in the frame of the certain model, which should be clearly stated. In our model [7,8] we proposed to describe the overall DNA curvature of a relatively small DNA fragment by a curvature of an arc approximating to the path of the axis of the given DNA fragment. Overall DNA shape is a meaningful notion only if the length of the DNA molecule is substantially larger than its diameter. This is, only fragments longer than 20 bp, with elongated shape, could have that shape curved. To predict the mentioned above DNA axis from the DNA sequence we used the nearest-neighbor approximation of the wedge model for the DNA

Label-free, atomic force microscopy-based mapping of DNA intrinsic curvature for the nanoscale comparative analysis of bent duplexes

We propose a method for the characterization of the local intrinsic curvature of adsorbed DNA molecules. It relies on a novel statistical chain descriptor, namely the ensemble averaged product of curvatures for two nanosized segments, symmetrically placed on the contour of atomic force microscopy imaged chains. We demonstrate by theoretical arguments and experimental investigation of representative samples that the fine mapping of the average product along the molecular backbone generates a characteristic pattern of variation that effectively highlights all pairs of DNA tracts with large intrinsic curvature. The centrosymmetric character of the chain descriptor enables targetting strands with unknown orientation. This overcomes a remarkable limitation of the current experimental strategies that estimate curvature maps solely from the trajectories of end-labeled molecules or palindromes. As a consequence our approach paves the way for a reliable, unbiased, label-free comparative analysis of bent duplexes, aimed to detect local conformational changes of physical or biological relevance in large sample numbers. Notably, such an assay is virtually inaccessible to the automated intrinsic curvature computation algorithms proposed so far. We foresee several challenging applications, including the validation of DNA adsorption and bending models by experiments and the discrimination of specimens for genetic screening purposes.

Interference of intrinsic curvature of DNA by DNA-intercalating agents

It has been demonstrated in our studies that the intrinsic curvature of DNA can be easily interrupted by low concentrations of chloroquine and ethidium bromide. In addition, the changes of DNA curvature caused by varying the concentration of these two DNA intercalators can be readily verified through using an atomic force microscope.

Intrinsic curvature in duplex DNA inhibits Human Topoisomerase I

Human Topoisomerase I (hTopo I) have been known as a potential target for cancer therapy. A series of duplex DNA with different intrinsic curvatures have been designed as inhibitors to hTopo I. The activities of hTopo I on relaxing supercoiled plasmid pUC 19 are apparently diminished in the presence of the curved DNA. More potent inhibitions and smaller IC(50) are achieved by duplex DNA with higher curvatures. EMSA indicates that hTopo I can recognize the curved DNA through binding interactions. Our studies demonstrate that the activity of hTopo I can be modulated by the intrinsic curvature of linear DNA and provide a new avenue to design curved DNA as hTopo I inhibitors with high therapeutic efficiency and low toxicity.

PerPlot & PerScan: tools for analysis of DNA curvature-related periodicity in genomic nucleotide sequences

BackgroundPeriodic spacing of short adenine or thymine runs phased with DNA helical period of ~10.5 bp is associated with intrinsic DNA curvature and deformability, which play important roles in DNA-protein interactions and in the organization of chromosomes in both eukaryotes and prokaryotes. Local differences in DNA sequence periodicity have been linked to differences in gene expression in some organisms. Despite the significance of these periodic patterns, there are virtually no publicly accessible tools for their analysis.ResultsWe present novel tools suitable for assessments of DNA curvature-related sequence periodicity in nucleotide sequences at the genome scale. Utility of the present software is demonstrated on a comparison of sequence periodicities in the genomes of Haemophilus influenzae, Methanocaldococcus jannaschii, Saccharomyces cerevisiae, and Arabidopsis thaliana. The software can be accessed through a web interface and the programs are also available for download.ConclusionsThe present software is suitable for comparing DNA curvature-related sequence periodicity among different genomes as well as for analysis of intrachromosomal heterogeneity of the sequence periodicity. It provides a quick and convenient way to detect anomalous regions of chromosomes that could have unusual structural and functional properties and/or distinct evolutionary history.

Contribution of nonlocal interactions to DNA elasticity

A nonlocal harmonic elastic rod model is proposed to describe the elastic behavior of short DNA molecules. We show that the nonlocal interactions contribute to effective bending energy of the molecule and affect its apparent persistence length. It is also shown that the anomalous behavior which has been observed in all-atom molecular dynamic simulations [A. K. Mazur, Biophys. J. 134, 4507 (2006)] can be a consequence of both nonlocal interactions between DNA base pairs and the intrinsic curvature of DNA.

Computational Analysis of Looping of a Large Family of Highly Bent DNA by LacI

Sequence-dependent intrinsic curvature of DNA influences looping by regulatory proteins such as LacI and NtrC. Curvature can enhance stability and control shape, as observed in LacI loops formed with three designed sequences with operators bracketing an A-tract bend. We explore geometric, topological, and energetic effects of curvature with an analysis of a family of highly bent sequences, using the elastic rod model from previous work. A unifying straight-helical-straight representation uses two phasing parameters to describe sequences composed of two straight segments that flank a common helically supercoiled segment. We exercise the rod model over this two-dimensional space of phasing parameters to evaluate looping behaviors. This design space is found to comprise two subspaces that prefer parallel versus anti-parallel binding topologies. The energetic cost of looping varies from 4 to 12 kT. Molecules can be designed to yield distinct binding topologies as well as hyperstable or hypostable loops and potentially loops that can switch conformations. Loop switching could be a mechanism for control of gene expression. Model predictions for linking numbers and sizes of LacI-DNA loops can be tested using multiple experimental approaches, which coupled with theory could address whether proteins or DNA provide the observed flexibility of protein-DNA loops.