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.

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 (DNANCP) 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 DNANCP 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 DNANCP 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 DNANCP 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 DNANCP 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.

Histone deposition promotes recombination-dependent replication at arrested forks.

Replication stress poses a serious threat to genome stability. Recombination-Dependent-Replication (RDR) promotes DNA synthesis resumption from arrested forks. Despite the identification of chromatin restoration pathways after DNA repair, crosstalk coupling RDR and chromatin assembly is largely unexplored. The fission yeast Chromatin Assembly Factor-1, CAF-1, is known to promote RDR. Here, we addressed the contribution of histone deposition to RDR. We expressed a mutated histone, H3-H113D, to genetically alter replication-dependent chromatin assembly by destabilizing (H3-H4)2 tetramer. We established that DNA synthesis-dependent histone deposition, by CAF-1 and Asf1, promotes RDR by preventing Rqh1-mediated disassembly of joint-molecules. The recombination factor Rad52 promotes CAF-1 binding to sites of recombination-dependent DNA synthesis, indicating that histone deposition occurs downstream Rad52. Histone deposition and Rqh1 activity act synergistically to promote cell resistance to camptothecin, a topoisomerase I inhibitor that induces replication stress. Moreover, histone deposition favors non conservative recombination events occurring spontaneously in the absence of Rqh1, indicating that the stabilization of joint-molecules by histone deposition also occurs independently of Rqh1 activity. These results indicate that histone deposition plays an active role in promoting RDR, a benefit counterbalanced by stabilizing at-risk joint-molecules for genome stability.

Structural Explorations of NCp7-Nucleic Acid Complexes Give Keys to Decipher the Binding Process.

A comprehensive view of all the structural aspects related to NCp7 is essential to understand how this protein, crucial in many steps of the HIV-1 cycle, binds and anneals nucleic acids (NAs), mainly thanks to two zinc fingers, ZF1 and ZF2. Here, we inspected the structural properties of the available experimental models of NCp7 bound to either DNA or RNA molecules, or free of ligand. Our analyses included the characterization of the relative positioning of ZF1 and ZF2, accessibility measurements and the exhaustive, quantitative mapping of the contacts between amino acids and nucleotides by a recent tessellation method, VLDM. This approach unveiled the intimate connection between NA binding process and the conformations explored by the free protein. It also provided new insights into the functional specializations of ZF1 and ZF2. The larger accessibility of ZF2 in free NCp7 and the consistency of the ZF2/NA interface in different models and conditions give ZF2 the lead of the binding process. ZF1 contributes to stabilize the complexes through various organizations of the ZF1/NA interface. This work outcome is a global binding scheme of NCp7 to DNA and RNA, and an example of how protein-NA complexes are stabilized.

Holding the Nucleosome Together: A Quantitative Description of the DNA-Histone Interface in Solution.

The nucleosome is the fundamental unit of eukaryotic genome packaging in the chromatin. In this complex, the DNA wraps around eight histone proteins to form a superhelical double helix. The resulting bending, stronger than anything observed in free DNA, raises the question of how such a distortion is stabilized by the proteic and solvent environments. In this work, the DNA-histone interface in solution was exhaustively analyzed from nucleosome structures generated by molecular dynamics. An original Voronoi tessellation technique, measuring the topology of interacting elements without any empirical or subjective adjustment, was used to characterize the interface in terms of contact area and occurrence. Our results revealed an interface more robust than previously known, combining extensive, long-lived nonelectrostatic and electrostatic interactions between DNA and both structured and unstructured histone regions. Cation accumulation makes the proximity of juxtaposed DNA gyres in the superhelix possible by shielding the strong electrostatic repulsion of the charged phosphate groups. Overall, this study provides new insights on the nucleosome cohesion, explaining how DNA distortions can be maintained in a nucleoprotein complex.

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.

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(31)P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (Dinter) 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 Dinterare 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.

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

We investigated how the intrinsic sequence-dependent properties probed via the phosphate linkages (BI ↔ BII equilibrium) influence the preferred shape of free DNA, and how this affects the nucleosome formation. First, this exploits NMR solution studies of four B-DNA dodecamers that together cover 39 base pairs of the 5' half of the sequence 601, of special interest for nucleosome formation. The results validate our previous prediction of a systematic, general sequence effect on the intrinsic backbone BII propensities. NMR provides new evidence that the backbone behavior is intimately coupled to the minor groove width. Second, application of the backbone behavior predictions to the full sequence 601 and other relevant sequences demonstrates that alternation of intrinsic low and high BII propensities, coupled to intrinsic narrow and wide minor grooves, largely coincides with the sinusoidal variations of the DNA minor groove width observed in crystallographic structures of the nucleosome. This correspondence is much poorer with low affinity sequences. Overall, the results indicate that nucleosome formation involves an indirect readout process implicating pre-existing DNA minor groove conformations. It also illustrates how the prediction of the intrinsic structural DNA behavior offers a powerful framework to gain explanatory insight on how proteins read DNA.

Simulating DNA by molecular dynamics: aims, methods, and validation.

The structure and dynamics of the B-DNA double helix involves subtle sequence-dependent effects which are decisive for its function, but difficult to characterize. These structural and dynamic effects can be addressed by simulations of DNA sequences in explicit solvent. Here, we present and discuss the state-of-art of B-DNA molecular dynamics simulations with the major force fields in use today. We explain why a critical analysis of the MD trajectories is required to assess their reliability, and estimate the value and limitations of these models. Overall, simulations of DNA bear great promise towards deciphering the structural and physical subtleties of this biopolymer, where much remains to be understood.

Mg2+ in the major groove modulates B-DNA structure and dynamics.

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.

How cations can assist DNase I in DNA binding and hydrolysis.

DNase I requires Ca²+ and Mg²+ for hydrolyzing double-stranded DNA. However, the number and the location of DNase I ion-binding sites remain unclear, as well as the role of these counter-ions. Using molecular dynamics simulations, we show that bovine pancreatic (bp) DNase I contains four ion-binding pockets. Two of them strongly bind Ca²+ while the other two sites coordinate Mg²+. These theoretical results are strongly supported by revisiting crystallographic structures that contain bpDNase I. One Ca²+ stabilizes the functional DNase I structure. The presence of Mg²+ in close vicinity to the catalytic pocket of bpDNase I reinforces the idea of a cation-assisted hydrolytic mechanism. Importantly, Poisson-Boltzmann-type electrostatic potential calculations demonstrate that the divalent cations collectively control the electrostatic fit between bpDNase I and DNA. These results improve our understanding of the essential role of cations in the biological function of bpDNase I. The high degree of conservation of the amino acids involved in the identified cation-binding sites across DNase I and DNase I-like proteins from various species suggests that our findings generally apply to all DNase I-DNA interactions.

Understanding the sequence-dependence of DNA groove dimensions: implications for DNA interactions.

BACKGROUND: The B-DNA major and minor groove dimensions are crucial for DNA-protein interactions. It has long been thought that the groove dimensions depend on the DNA sequence, however this relationship has remained elusive. Here, our aim is to elucidate how the DNA sequence intrinsically shapes the grooves. METHODOLOGY/PRINCIPAL FINDINGS: The present study is based on the analysis of datasets of free and protein-bound DNA crystal structures, and from a compilation of NMR (31)P chemical shifts measured on free DNA in solution on a broad range of representative sequences. The (31)P chemical shifts can be interpreted in terms of the BI↔BII backbone conformations and dynamics. The grooves width and depth of free and protein-bound DNA are found to be clearly related to the BI/BII backbone conformational states. The DNA propensity to undergo BI↔BII backbone transitions is highly sequence-dependent and can be quantified at the dinucleotide level. This dual relationship, between DNA sequence and backbone behavior on one hand, and backbone behavior and groove dimensions on the other hand, allows to decipher the link between DNA sequence and groove dimensions. It also firmly establishes that proteins take advantage of the intrinsic DNA groove properties. CONCLUSIONS/SIGNIFICANCE: The study provides a general framework explaining how the DNA sequence shapes the groove dimensions in free and protein-bound DNA, with far-reaching implications for DNA-protein indirect readout in both specific and non specific interactions.

Multi-resolution approach for interactively locating functionally linked ion binding sites by steering small molecules into electrostatic potential maps using a haptic device.

Metal ions drive important parts of biology, yet it remains experimentally challenging to locate their binding sites. Here we present an innovative computational approach. We use interactive steering of charged ions or small molecules in an electrostatic potential map in order to identify potential binding sites. The user interacts with a haptic device and experiences tactile feedback related to the strength of binding at a given site. The potential field is the first level of resolution used in this model. Any type of potential field can be used, implicitly taking into account various conditions such as ionic strength, dielectric constants or the presence of a membrane. At a second level, we represent the accessibility of all binding sites by modelling the shape of the target macromolecule via non-bonded van der Waals interactions between its static atomic or coarse-grained structure and the probe molecule(s). The third independent level concerns the representation of the molecular probe itself. Ion selectivity can be assessed by using multiple interacting ions as probes. This method was successfully applied to the DNase I enzyme, where we recently identified two new cation binding sites by computationally expensive extended molecular dynamics simulations.

Intrinsic flexibility of B-DNA: the experimental TRX scale.

B-DNA flexibility, crucial for DNA-protein recognition, is sequence dependent. Free DNA in solution would in principle be the best reference state to uncover the relation between base sequences and their intrinsic flexibility; however, this has long been hampered by a lack of suitable experimental data. We investigated this relationship by compiling and analyzing a large dataset of NMR (31)P chemical shifts in solution. These measurements reflect the BI <--> BII equilibrium in DNA, intimately correlated to helicoidal descriptors of the curvature, winding and groove dimensions. Comparing the ten complementary DNA dinucleotide steps indicates that some steps are much more flexible than others. This malleability is primarily controlled at the dinucleotide level, modulated by the tetranucleotide environment. Our analyses provide an experimental scale called TRX that quantifies the intrinsic flexibility of the ten dinucleotide steps in terms of Twist, Roll, and X-disp (base pair displacement). Applying the TRX scale to DNA sequences optimized for nucleosome formation reveals a 10 base-pair periodic alternation of stiff and flexible regions. Thus, DNA flexibility captured by the TRX scale is relevant to nucleosome formation, suggesting that this scale may be of general interest to better understand protein-DNA recognition.

DNA structures from phosphate chemical shifts.

For B-DNA, the strong linear correlation observed by nuclear magnetic resonance (NMR) between the (31)P chemical shifts (deltaP) and three recurrent internucleotide distances demonstrates the tight coupling between phosphate motions and helicoidal parameters. It allows to translate deltaP into distance restraints directly exploitable in structural refinement. It even provides a new method for refining DNA oligomers with restraints exclusively inferred from deltaP. Combined with molecular dynamics in explicit solvent, these restraints lead to a structural and dynamical view of the DNA as detailed as that obtained with conventional and more extensive restraints. Tests with the Jun-Fos oligomer show that this deltaP-based strategy can provide a simple and straightforward method to capture DNA properties in solution, from routine NMR experiments on unlabeled samples.

Sequence-dependent DNA flexibility mediates DNase I cleavage.

Understanding the preference of nonspecific proteins for certain DNA structural features requires an accurate description of the properties of free DNA, especially regarding their possible predisposition to adopt a conformation that favors the formation of a complex. Exploiting previous exhaustive NMR studies performed on free DNA oligomers, we investigated the molecular basis of DNase I sensitivity under conditions where DNase I binding limits the probability of cleavage. We showed that cleavage intensity was correlated with adjacent 3' phosphate linkage flexibility, monitored by (31)P chemical shifts. Examining NMR-refined DNA structures highlighted that sequence-dependent flexible phosphates were associated with large minor groove variations that may promote the affinity of DNase I, according to relevant DNA-protein complexes. In sum, this work demonstrates that specificity in DNA-DNase I interaction is mediated by DNA flexibility, which influences the induced-fit transitions required to form productive complexes.

Ets-1 p51 and p42 isoforms differentially modulate Stromelysin-1 promoter according to induced DNA bend orientation.

The Stromelysin-1 gene promoter contains a palindrome of two Ets-binding sites (EBS) that bind the p51 and p42 isoforms of the human Ets-1-transcription factor. A previous study established that full gene transactivation is associated with a ternary complex consisting of two p51 bound to the two EBS on the promoter. p42, only able to bind one of the two EBS, induces only very weak activity. Here, we investigate the mechanism by which the Stromelysin-1 promoter discriminates between p51 and p42. The differential stoichiometry of the two Ets-1 isoforms arises from the Stromelysin-1 EBS palindrome. The ternary complex requires the presence of two inhibitory domains flanking the DNA-binding domain and the ability to form an intramolecular autoinhibition module. Most importantly, the p51-ternary and the p42-binary complexes induce DNA curvatures with opposite orientations. These results establish that differential DNA bending, via p51 and p42 differential binding, is correlated with the Stromelysin-1 promoter activation process.

An optimized extended DNA kappa B site that enhances plasmid DNA nuclear import and gene expression.

BACKGROUND: The nuclear factor kappa B (NF kappaB) transcription factor, which shuttles between the cytoplasm and the nucleus under specific conditions, is a suitable intracellular target to increase the nuclear import of plasmid DNA. We report the design of an optimized and extended NF kappaB DNA binding sequence that promotes an efficient plasmid nuclear import. METHODS: On the basis of structural studies, the 5'-CTGGGGACTTTCCAGCTGGGGACTTTCCAGCTGGGGACTTTCCAGG-3' segment (termed 3NF) comprising three 10-bp kappaB sites (GGGACTTTCC) separated by a 5-bp optimized spacer (AGCTG) was selected for its capacity to ensure the best structural fit with NF kappaB and to fix simultaneously three proteins. Plasmids encoding luciferase and bearing this sequence (3NF-plasmids) were constructed and their nuclear import and gene expression efficiencies compared with that of plasmids containing classical kappaB motifs. RESULTS: A high luciferase expression was associated with plasmids containing one (p3NF-luc) or two (p3NF-luc-3NF) 3NF sequences. In situ hybridization experiments and quantitative measurement of the number of plasmid copies demonstrated that the nuclear delivery of 3NF-plasmids was more efficient than that of 3NF-free plasmids. Cross-linked immunoprecipitation showed that 3NF-plasmids were recognized by NF kappaB inside cells upon transfection. The nuclear delivery was inhibited with BAY 11-7085, an inhibitor of NF kappaB activation. Finally, p3NF-luc-3NF, the most efficient construct for in vitro transfection, had a long-lived luciferase expression in vivo. CONCLUSIONS: The results obtained in the present study demonstrate the NF kappaB-mediated nuclear delivery of 3NF-plasmids. Due to its high affinity for fixing several NF kappaB, the 3NF sequence is a very promising helper for a nonviral gene delivery system.

Importance of accurate DNA structures in solution: the Jun-Fos model.

Understanding the recognition of DNA sequences by proteins requires an accurate description of the structural dynamics of free DNA, especially regarding indirect readout. This involves subtle sequence-dependent effects that are difficult to characterize in solution. To progress in this area, we applied NMR and extensive simulations to a DNA sequence relevant to the Jun-Fos system. The backbone and base behaviors demonstrate that unrestrained simulations with major force fields (Parm98, Parmbsc0, and CHARMM27) are not reliable enough for in silico predictions of detailed DNA structures. More realistic structures required molecular dynamics simulations supplemented by NMR restraints. A new methodological element involved restraints inferred from the phosphate chemical shifts and from the phosphate dynamics. This provided a detailed and dynamic view of the intrinsic properties of the free DNA sequence that can be related to its recognition, by comparison with a relevant DNA-protein complex. We show how to exploit the relationship between phosphate motions and helicoidal descriptors for structure determination toward an accurate description of DNA structures and dynamics in solution.