DNA stress and strain, in silico, in vitro and in vivo

R-Loop Physiology and Pathology: A Brief Review

Physiological and pathological roles for R-loop structures continue to be discovered, and studies suggest that R-loops could contribute to human disease. R-loops are nucleic acid structures characterized by a DNA:RNA hybrid and displaced single-stranded DNA that occur in connection with transcription. R-loops form naturally and have been shown to be important for a number of physiological processes such as mitochondrial replication initiation, class switch recombination, DNA repair, modulating DNA topology, and regulation of gene expression. However, subsets of R-loops or persistent R-loops lead to DNA breaks, chromosome rearrangement, and genome instability. In addition, R-loops have been linked to human diseases, specifically neurological disorders and cancer. Of the large amount of research produced recently on R-loops, this review covers evidence for R-loop involvement in normal cellular physiology and pathophysiology, as well as describing factors that contribute to R-loop regulation.

Emerging roles for R-loop structures in the management of topological stress

R-loop structures are a prevalent class of alternative non-B DNA structures that form during transcription upon invasion of the DNA template by the nascent RNA. R-loops form universally in the genomes of organisms ranging from bacteriophages, bacteria, and yeasts to plants and animals, including mammals. A growing body of work has linked these structures to both physiological and pathological processes, in particular to genome instability. The rising interest in R-loops is placing new emphasis on understanding the fundamental physicochemical forces driving their formation and stability. Pioneering work in Escherichia coli revealed that DNA topology, in particular negative DNA superhelicity, plays a key role in driving R-loops. A clear role for DNA sequence was later uncovered. Here, we review and synthesize available evidence on the roles of DNA sequence and DNA topology in controlling R-loop formation and stability. Factoring in recent developments in R-loop modeling and single-molecule profiling, we propose a coherent model accounting for the interplay between DNA sequence and DNA topology in driving R-loop structure formation. This model reveals R-loops in a new light as powerful and reversible topological stress relievers, an insight that significantly expands the repertoire of R-loops' potential biological roles under both normal and aberrant conditions.

Nucleosome stability measured in situ by automated quantitative imaging

Current approaches have limitations in providing insight into the functional properties of particular nucleosomes in their native molecular environment. Here we describe a simple and powerful method involving elution of histones using intercalators or salt, to assess stability features dependent on DNA superhelicity and relying mainly on electrostatic interactions, respectively, and measurement of the fraction of histones remaining chromatin-bound in the individual nuclei using histone type- or posttranslational modification- (PTM-) specific antibodies and automated, quantitative imaging. The method has been validated in H3K4me3 ChIP-seq experiments, by the quantitative assessment of chromatin loop relaxation required for nucleosomal destabilization, and by comparative analyses of the intercalator and salt induced release from the nucleosomes of different histones. The accuracy of the assay allowed us to observe examples of strict association between nucleosome stability and PTMs across cell types, differentiation state and throughout the cell-cycle in close to native chromatin context, and resolve ambiguities regarding the destabilizing effect of H2A.X phosphorylation. The advantages of the in situ measuring scenario are demonstrated via the marked effect of DNA nicking on histone eviction that underscores the powerful potential of topological relaxation in the epigenetic regulation of DNA accessibility.

The competing mini-dumbbell mechanism: new insights into CCTG repeat expansion

CCTG repeat expansions in intron 1 of the cellular nucleic acid-binding protein gene are associated with myotonic dystrophy type 2. Recently, we have reported a novel mini-dumbbell (MDB) structure formed by two CCTG or TTTA repeats, which potentially has a critical role in repeat expansions. Here we present a mechanism, called the competing MDB mechanism, to explain how the formation of MDB can lead to efficient mismatch repair (MMR) escape and thus CCTG repeat expansions during DNA replication. In a long tract of CCTG repeats, two competing MDBs can be formed in any segment of three repeats. Fast exchange between these MDBs will make the commonly occupied repeat behave like a mini-loop. Further participations of the 5'- or 3'-flanking repeat in forming competing MDBs will make the mini-loop shift in the 5'- or 3'-direction, thereby providing a pathway for the mini-loop to escape from MMR. To avoid the complications due to the formation of hairpin conformers in longer CCTG repeats, we made use of TTTA repeats as model sequences to demonstrate the formation of competing MDBs and shifting of mini-loop in a long tract of repeating sequence.

Unusual structures of CCTG repeats and their participation in repeat expansion

CCTG repeat expansion in intron 1 of the cellular nucleic acid-binding protein (CNBP) gene has been identified to be the genetic cause of myotonic dystrophy type 2 (DM2). Yet the underlying reasons for the genetic instability in CCTG repeats remain elusive. In recent years, CCTG repeats have been found to form various types of unusual secondary structures including mini-dumbbell (MDB), hairpin and dumbbell, revealing that there is a high structural diversity in CCTG repeats intrinsically. Upon strand slippage, the formation of unusual structures in the nascent strand during DNA replication has been proposed to be the culprit of CCTG repeat expansions. On the one hand, the thermodynamic stability, size, and conformational dynamics of these unusual structures affect the propensity of strand slippage. On the other hand, these structural properties determine whether the unusual structure can successfully escape from DNA repair. In this short overview, we first summarize the recent advances in elucidating the solution structures of CCTG repeats. We then discuss the potential pathways by which these unusual structures bring about variable sizes of repeat expansion, high strand slippage propensity and efficient repair escape.

Potential non-B DNA regions in the human genome are associated with higher rates of nucleotide mutation and expression variation

While individual non-B DNA structures have been shown to impact gene expression, their broad regulatory role remains elusive. We utilized genomic variants and expression quantitative trait loci (eQTL) data to analyze genome-wide variation propensities of potential non-B DNA regions and their relation to gene expression. Independent of genomic location, these regions were enriched in nucleotide variants. Our results are consistent with previously observed mutagenic properties of these regions and counter a previous study concluding that G-quadruplex regions have a reduced frequency of variants. While such mutagenicity might undermine functionality of these elements, we identified in potential non-B DNA regions a signature of negative selection. Yet, we found a depletion of eQTL-associated variants in potential non-B DNA regions, opposite to what might be expected from their proposed regulatory role. However, we also observed that genes downstream of potential non-B DNA regions showed higher expression variation between individuals. This coupling between mutagenicity and tolerance for expression variability of downstream genes may be a result of evolutionary adaptation, which allows reconciling mutagenicity of non-B DNA structures with their location in functionally important regions and their potential regulatory role.

Integration of syntactic and semantic properties of the DNA code reveals chromosomes as thermodynamic machines converting energy into information

Understanding genetic regulation is a problem of fundamental importance. Recent studies have made it increasingly evident that, whereas the cellular genetic regulation system embodies multiple disparate elements engaged in numerous interactions, the central issue is the genuine function of the DNA molecule as information carrier. Compelling evidence suggests that the DNA, in addition to the digital information of the linear genetic code (the semantics), encodes equally important continuous, or analog, information that specifies the structural dynamics and configuration (the syntax) of the polymer. These two DNA information types are intrinsically coupled in the primary sequence organisation, and this coupling is directly relevant to regulation of the genetic function. In this review, we emphasise the critical need of holistic integration of the DNA information as a prerequisite for understanding the organisational complexity of the genetic regulation system.

Cellular MYCro economics: Balancing MYC function with MYC expression

The expression levels of the MYC oncoprotein have long been recognized to be associated with the outputs of major cellular processes including proliferation, cell growth, apoptosis, differentiation, and metabolism. Therefore, to understand how MYC operates, it is important to define quantitatively the relationship between MYC input and expression output for its targets as well as the higher-order relationships between the expression levels of subnetwork components and the flow of information and materials through those networks. Two different views of MYC are considered, first as a molecular microeconomic manager orchestrating specific positive and negative responses at individual promoters in collaboration with other transcription and chromatin components, and second, as a macroeconomic czar imposing an overarching rule onto all active genes. In either case, c-myc promoter output requires multiple inputs and exploits diverse mechanisms to tune expression to the appropriate levels relative to the thresholds of expression that separate health and disease.

DNA thermodynamics shape chromosome organization and topology

How much information is encoded in the DNA sequence of an organism? We argue that the informational, mechanical and topological properties of DNA are interdependent and act together to specify the primary characteristics of genetic organization and chromatin structures. Superhelicity generated in vivo, in part by the action of DNA translocases, can be transmitted to topologically sensitive regions encoded by less stable DNA sequences.

The genome-wide distribution of non-B DNA motifs is shaped by operon structure and suggests the transcriptional importance of non-B DNA structures in Escherichia coli

Although the right-handed double helical B-form DNA is most common under physiological conditions, DNA is dynamic and can adopt a number of alternative structures, such as the four-stranded G-quadruplex, left-handed Z-DNA, cruciform and others. Active transcription necessitates strand separation and can induce such non-canonical forms at susceptible genomic sequences. Therefore, it has been speculated that these non-B DNA motifs can play regulatory roles in gene transcription. Such conjecture has been supported in higher eukaryotes by direct studies of several individual genes, as well as a number of large-scale analyses. However, the role of non-B DNA structures in many lower organisms, in particular proteobacteria, remains poorly understood and incompletely documented. In this study, we performed the first comprehensive study of the occurrence of B DNA-non-B DNA transition-susceptible sites (non-B DNA motifs) within the context of the operon structure of the Escherichia coli genome. We compared the distributions of non-B DNA motifs in the regulatory regions of operons with those from internal regions. We found an enrichment of some non-B DNA motifs in regulatory regions, and we show that this enrichment cannot be simply explained by base composition bias in these regions. We also showed that the distribution of several non-B DNA motifs within intergenic regions separating divergently oriented operons differs from the distribution found between convergent ones. In particular, we found a strong enrichment of cruciforms in the termination region of operons; this enrichment was observed for operons with Rho-dependent, as well as Rho-independent terminators. Finally, a preference for some non-B DNA motifs was observed near transcription factor-binding sites. Overall, the conspicuous enrichment of transition-susceptible sites in these specific regulatory regions suggests that non-B DNA structures may have roles in the transcriptional regulation of specific operons within the E. coli genome.

Thermodynamics of writhe in DNA minicircles from molecular dynamics simulations

DNA supercoiling plays a role in genetic control by imposing torsional stress. This can induce writhe, which changes the global shape of the DNA. We have used atomistic molecular dynamics simulations to partition the free energy changes driving the writhing and unwrithing transitions in supercoiled minicircle DNA. The calculations show that while writhing is energetically driven, the unwrithing transition occurs because the circular state has a higher configurational entropy than the plectoneme. Writhing improves the van der Waals interactions between stacked bases, but can be suppressed by electrostatic repulsion within the negatively charged backbone strands in low salt conditions where electrostatic screening is poor. The free energy difference between circular and plectonemic DNA is determined by such a delicate balance of opposing thermodynamic terms that any perturbation in the environment, such as a change in salt concentration, can be sufficient to convert between these two states. This switchable behavior provides a mechanism for supercoiled DNA to store and communicate biological information physically as well as chemically.

DNA thermodynamic stability and supercoil dynamics determine the gene expression program during the bacterial growth cycle

The chromosomal DNA polymer constituting the cellular genetic material is primarily a device for coding information. Whilst the gene sequences comprise the digital (discontinuous) linear code, physiological alterations of the DNA superhelical density generate in addition analog (continuous) three-dimensional information essential for regulation of both chromosome compaction and gene expression. Insight into the relationship between the DNA analog information and the digital linear code is of fundamental importance for understanding genetic regulation. Our previous study in the model organism Escherichia coli suggested that the chromosomal gene order and a spatiotemporal gradient of DNA superhelicity associated with DNA replication determine the growth phase-dependent gene transcription. In this study we reveal a general gradient of DNA thermodynamic stability correlated with the polarity of chromosomal replication and manifest in the spatiotemporal pattern of gene transcription during the bacterial growth cycle. Furthermore, by integrating the physical and dynamic features of the transcribed sequences with their functional content we identify spatiotemporal domains of gene expression encompassing different functions. We thus provide both an insight into the organisational principle of the bacterial growth program and a novel holistic methodology for exploring chromosomal dynamics.

The condensed chromatin fiber: an allosteric chemo-mechanical machine for signal transduction and genome processing

Allostery is a key concept of molecular biology which refers to the control of an enzyme activity by an effector molecule binding the enzyme at another site rather than the active site (allos = other in Greek). We revisit here allostery in the context of chromatin and argue that allosteric principles underlie and explain the functional architecture required for spacetime coordination of gene expression at all scales from DNA to the whole chromosome. We further suggest that this functional architecture is provided by the chromatin fiber itself. The structural, mechanical and topological features of the chromatin fiber endow chromosomes with a tunable signal transduction from specific (or nonspecific) effectors to specific (or nonspecific) active sites. Mechanical constraints can travel along the fiber all the better since the fiber is more compact and regular, which speaks in favor of the actual existence of the (so-called 30 nm) chromatin fiber. Chromatin fiber allostery reconciles both the physical and biochemical approaches of chromatin. We illustrate this view with two supporting specific examples. Moreover, from a methodological point of view, we suggest that the notion of chromatin fiber allostery is particularly relevant for systemic approaches. Finally we discuss the evolutionary power of allostery in the context of chromatin and its relation to modularity.