NMR studies of DNA support the role of pre-existing minor groove variations in nucleosome indirect readout

Nucleosome assembly and disassembly pathways in vitro

Structural fluctuations of nucleosomes modulate the access to internal DNA in eukaryotic cells; clearly characterisation of this fundamental process is crucial to understanding gene regulation. Here we apply PhAST (Photochemical Analysis of Structural Transitions) to monitor at a base pair level, structural alterations induced all along the DNA upon histone binding or release. By offering the first reliable, detailed comparison of nucleosome assembly and disassembly in vitro, we reveal similarities and differences between the two processes. We identify multiple, sequential intermediate states characterised by specific PhAST signals whose localisation and amplitude reflect asymmetries of DNA/histone interactions with respect to the nucleosome pseudo dyad. These asymmetries involve not only the DNA extremities but also regions close to the pseudo dyad. Localisations of asymmetries develop in a consistent manner during both assembly and disassembly processes; they primarily reflect the DNA sequence effect on the efficiency of DNA-histone binding. More unexpectedly, the amplitude component of PhAST signals not only evolves as a function of intermediate states but does so differently between assembly and disassembly pathways. Our observation of differences between assembly and disassembly opens up new avenues to define the role of the DNA sequence in processes underlying the regulation of gene expression. Overall, we provide new insights into how the intrinsic properties of DNA are integrated into a holistic mechanism that controls chromatin structure.

Kinetics and thermodynamics of BI-BII interconversion altered by T:G mismatches in DNA

T:G mismatches in DNA result in humans primarily from deamination of methylated CpG sites. They are repaired by redundant systems, such as thymine DNA glycosylase (TDG) and methyl-binding domain enzyme (MBD4), and maintenance of these sites has been implicated in epigenetic processes. The process by which these enzymes identify a canonical DNA base in the incorrect basepairing context remains a mystery. However, the conserved contacts of the repair enzymes with the DNA backbone suggests a role for protein-phosphate interaction in the recognition and repair processes. We have used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion, and for this work have focused on alterations to the activation barriers to interconversion and the effect of a mismatch compared with canonical DNA. We have found that alterations to the ΔG of interconversion for T:G basepairs are remarkably similar to U:G basepairs in the form of stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA, suggesting a universality of this result for TDG substrates. Likewise, we see perturbations to the free energy (∼1 kcal/mol) and enthalpy (2-5 kcal/mol) of activation for the BI-BII interconversion localized to the phosphates flanking the mismatch. Overall our results strongly suggest that the perturbed backbone energetics in T:G basepairs play a significant role in the recognition process of DNA repair enzymes.

Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N6-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.

Slow motions in A·T rich DNA sequence

In free B-DNA, slow (microsecond-to-millisecond) motions that involve equilibrium between Watson-Crick (WC) and Hoogsteen (HG) base-pairing expand the DNA dynamic repertoire that could mediate DNA-protein assemblies. R_1ρ relaxation dispersion NMR methods are powerful tools to capture such slow conformational exchanges in solution using ¹³C/¹⁵ N labelled DNA. Here, these approaches were applied to a dodecamer containing a TTAAA element that was assumed to facilitate nucleosome formation. NMR data and inferred exchange parameters assign HG base pairs as the minor, transient conformers specifically observed in three successive A·T base pairs forming the TAA·TTA segment. The abundance of these HG A·T base pairs can be up to 1.2% which is high compared to what has previously been observed. Data analyses support a scenario in which the three adenines undergo non-simultaneous motions despite their spatial proximity, thus optimising the probability of having one HG base pair in the TAA·TTA segment. Finally, revisiting previous NMR data on H2 resonance linewidths on the basis of our results promotes the idea of there being a special propensity of A·T base pairs in TAA·TTA tracts to adopt HG pairing. In summary, this study provides an example of a DNA functional element submitted to slow conformational exchange. More generally, it strengthens the importance of the role of the DNA sequence in modulating its dynamics, over a nano- to milli-second time scale.

A dynamic view of DNA structure within the nucleosome: Biological implications

Most of eukaryotic cellular DNA is packed in nucleosome core particles (NCPs), in which the DNA (DNA_NCP) is wrapped around histones. The influence of this organization on the intrinsic local dynamics of DNA is largely unknown, in particular because capturing such information from experiments remains notoriously challenging. Given the importance of dynamical properties in DNA functions, we addressed this issue using CHARMM36 MD simulations of a nucleosome containing the NCP positioning 601 sequence and four related free dodecamers. Comparison between DNA_NCP and free DNA reveals a limited impact of the dense DNA-histone interface on correlated motions of dinucleotide constituents and on fluctuations of inter base pair parameters. A characteristic feature intimately associated with the DNA_NCP super-helical path is a set of structural periodicities that includes a marked alternation of regions enriched in backbone BI and BII conformers. This observation led to uncover a convincing correspondence between the sequence effect on BI/BII propensities in both DNA_NCP and free DNA, strengthening the idea that the histone preference for particular DNA sequences relies on those intrinsic structural properties. These results offer for the first time a detailed view of the DNA dynamical behavior within NCP. They show in particular that the DNA_NCP dynamics is substantial enough to preserve the ability to structurally adjust to external proteins, for instance remodelers. Also, fresh structural arguments highlight the relevance of relationships between DNA sequence and structural properties for NCP formation. Overall, our work offers a more rational framework to approach the functional, biological roles of NCP.

Unspinning chromatin: Revealing the dynamic nucleosome landscape by NMR

NMR is an essential technique for obtaining information at atomic resolution on the structure, motions and interactions of biomolecules. Here, we review the contribution of NMR to our understanding of the fundamental unit of chromatin: the nucleosome. Nucleosomes compact the genome by wrapping the DNA around a protein core, the histone octamer, thereby protecting genomic integrity. Crucially, the imposed barrier also allows strict regulation of gene expression, DNA replication and DNA repair processes through an intricate system of histone and DNA modifications and a wide range of interactions between nucleosomes and chromatin factors. In this review, we describe how NMR has contributed to deciphering the molecular basis of nucleosome function. Starting from pioneering studies in the 1960s using natural abundance NMR studies, we focus on the progress in sample preparation and NMR methodology that has allowed high-resolution studies on the nucleosome and its subunits. We summarize the results and approaches of state-of-the-art NMR studies on nucleosomal DNA, histone complexes, nucleosomes and nucleosomal arrays. These studies highlight the particular strength of NMR in studying nucleosome dynamics and nucleosome-protein interactions. Finally, we look ahead to exciting new possibilities that will be afforded by on-going developments in solution and solid-state NMR. By increasing both the depth and breadth of nucleosome NMR studies, it will be possible to offer a unique perspective on the dynamic landscape of nucleosomes and its interacting proteins.

Bacterial pathogen gene regulation: a DNA-structure-centred view of a protein-dominated domain

The mechanisms used by bacterial pathogens to regulate the expression of their genes, especially their virulence genes, have been the subject of intense investigation for several decades. Whole genome sequencing projects, together with more targeted studies, have identified hundreds of DNA-binding proteins that contribute to the patterns of gene expression observed during infection as well as providing important insights into the nature of the gene products whose expression is being controlled by these proteins. Themes that have emerged include the importance of horizontal gene transfer to the evolution of pathogens, the need to impose regulatory discipline upon these imported genes and the important roles played by factors normally associated with the organization of genome architecture as regulatory principles in the control of virulence gene expression. Among these architectural elements is the structure of DNA itself, its variable nature at a topological rather than just at a base-sequence level and its ability to play an active (as well as a passive) part in the gene regulation process.

DNA supercoiling is a fundamental regulatory principle in the control of bacterial gene expression

Although it has become routine to consider DNA in terms of its role as a carrier of genetic information, it is also an important contributor to the control of gene expression. This regulatory principle arises from its structural properties. DNA is maintained in an underwound state in most bacterial cells and this has important implications both for DNA storage in the nucleoid and for the expression of genetic information. Underwinding of the DNA through reduction in its linking number potentially imparts energy to the duplex that is available to drive DNA transactions, such as transcription, replication and recombination. The topological state of DNA also influences its affinity for some DNA binding proteins, especially in DNA sequences that have a high A + T base content. The underwinding of DNA by the ATP-dependent topoisomerase DNA gyrase creates a continuum between metabolic flux, DNA topology and gene expression that underpins the global response of the genome to changes in the intracellular and external environments. These connections describe a fundamental and generalised mechanism affecting global gene expression that underlies the specific control of transcription operating through conventional transcription factors. This mechanism also provides a basal level of control for genes acquired by horizontal DNA transfer, assisting microbial evolution, including the evolution of pathogenic bacteria.

DNA supercoiling is a fundamental regulatory principle in the control of bacterial gene expression

Although it has become routine to consider DNA in terms of its role as a carrier of genetic information, it is also an important contributor to the control of gene expression. This regulatory principle arises from its structural properties. DNA is maintained in an underwound state in most bacterial cells and this has important implications both for DNA storage in the nucleoid and for the expression of genetic information. Underwinding of the DNA through reduction in its linking number potentially imparts energy to the duplex that is available to drive DNA transactions, such as transcription, replication and recombination. The topological state of DNA also influences its affinity for some DNA binding proteins, especially in DNA sequences that have a high A + T base content. The underwinding of DNA by the ATP-dependent topoisomerase DNA gyrase creates a continuum between metabolic flux, DNA topology and gene expression that underpins the global response of the genome to changes in the intracellular and external environments. These connections describe a fundamental and generalised mechanism affecting global gene expression that underlies the specific control of transcription operating through conventional transcription factors. This mechanism also provides a basal level of control for genes acquired by horizontal DNA transfer, assisting microbial evolution, including the evolution of pathogenic bacteria.

High-resolution biophysical analysis of the dynamics of nucleosome formation

We describe a biophysical approach that enables changes in the structure of DNA to be followed during nucleosome formation in in vitro reconstitution with either the canonical "Widom" sequence or a judiciously mutated sequence. The rapid non-perturbing photochemical analysis presented here provides 'snapshots' of the DNA configuration at any given moment in time during nucleosome formation under a very broad range of reaction conditions. Changes in DNA photochemical reactivity upon protein binding are interpreted as being mainly induced by alterations in individual base pair roll angles. The results strengthen the importance of the role of an initial (H3/H4)2 histone tetramer-DNA interaction and highlight the modulation of this early event by the DNA sequence. (H3/H4)2 binding precedes and dictates subsequent H2A/H2B-DNA interactions, which are less affected by the DNA sequence, leading to the final octameric nucleosome. Overall, our results provide a novel, exciting way to investigate those biophysical properties of DNA that constitute a crucial component in nucleosome formation and stabilization.

Flanking A·T basepairs destabilize the B(∗) conformation of DNA A-tracts

Capillary electrophoresis has been used to characterize the interaction of monovalent cations with 26-basepair DNA oligomers containing A-tracts embedded in flanking sequences with different basepair compositions. A 26-basepair random-sequence oligomer was used as the reference; lithium and tetrabutylammonium (TBA(+)) ions were used as the probe ions. The free solution mobilities of the A-tract and random-sequence oligomers were identical in solutions containing <∼ 100 mM cation. At higher cation concentrations, the A-tract oligomers migrated faster than the reference oligomer in TBA(+) and slower than the reference in Li(+). Hence, cations of different sizes can interact very differently with DNA A-tracts. The increased mobilities observed in TBA(+) suggest that the large hydrophobic TBA(+) ions are preferentially excluded from the vicinity of the A-tract minor groove, increasing the effective net charge of the A-tract oligomers and increasing the mobility. By contrast, Li(+) ions decrease the mobility of A-tract oligomers because of the preferential localization of Li(+) ions in the narrow A-tract minor groove. Embedding the A-tracts in AT-rich flanking sequences markedly alters preferential interactions of monovalent cations with the B(∗) conformation. Hence, A-tracts embedded in genomic DNA may or may not interact preferentially with monovalent cations, depending on the relative number of A · T basepairs in the flanking sequences.

The intrinsic mechanics of B-DNA in solution characterized by NMR

Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the ³¹P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (D_inter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while D_inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics

The accurate prediction of the structure and dynamics of DNA remains a major challenge in computational biology due to the dearth of precise experimental information on DNA free in solution and limitations in the DNA force-fields underpinning the simulations. A new generation of force-fields has been developed to better represent the sequence-dependent B-DNA intrinsic mechanics, in particular with respect to the BI ↔ BII backbone equilibrium, which is essential to understand the B-DNA properties. Here, the performance of MD simulations with the newly updated force-fields Parmbsc0εζOLI and CHARMM36 was tested against a large ensemble of recent NMR data collected on four DNA dodecamers involved in nucleosome positioning. We find impressive progress towards a coherent, realistic representation of B-DNA in solution, despite residual shortcomings. This improved representation allows new and deeper interpretation of the experimental observables, including regarding the behavior of facing phosphate groups in complementary dinucleotides, and their modulation by the sequence. It also provides the opportunity to extensively revisit and refine the coupling between backbone states and inter base pair parameters, which emerges as a common theme across all the complementary dinucleotides. In sum, the global agreement between simulations and experiment reveals new aspects of intrinsic DNA mechanics, a key component of DNA-protein recognition.