Structural insight into the sequence dependence of nucleosome positioning

DrugReAlign: a multisource prompt framework for drug repurposing based on large language models

Drug repurposing is a promising approach in the field of drug discovery owing to its efficiency and cost-effectiveness. Most current drug repurposing models rely on specific datasets for training, which limits their predictive accuracy and scope. The number of both market-approved and experimental drugs is vast, forming an extensive molecular space. Due to limitations in parameter size and data volume, traditional drug-target interaction (DTI) prediction models struggle to generalize well within such a broad space. In contrast, large language models (LLMs), with their vast parameter sizes and extensive training data, demonstrate certain advantages in drug repurposing tasks. In our research, we introduce a novel drug repurposing framework, DrugReAlign, based on LLMs and multi-source prompt techniques, designed to fully exploit the potential of existing drugs efficiently. Leveraging LLMs, the DrugReAlign framework acquires general knowledge about targets and drugs from extensive human knowledge bases, overcoming the data availability limitations of traditional approaches. Furthermore, we collected target summaries and target-drug space interaction data from databases as multi-source prompts, substantially improving LLM performance in drug repurposing. We validated the efficiency and reliability of the proposed framework through molecular docking and DTI datasets. Significantly, our findings suggest a direct correlation between the accuracy of LLMs' target analysis and the quality of prediction outcomes. These findings signify that the proposed framework holds the promise of inaugurating a new paradigm in drug repurposing.

Collaboration through chromatin: motors of transcription and chromatin structure

Packaging of the eukaryotic genome into chromatin places fundamental physical constraints on transcription. Clarifying how transcription operates within these constraints is essential to understand how eukaryotic gene expression programs are established and maintained. Here we review what is known about the mechanisms of transcription on chromatin templates. Current models indicate that transcription through chromatin is accomplished by the combination of an inherent nucleosome disrupting activity of RNA polymerase and the action of ATP-dependent chromatin remodeling motors. Collaboration between these two types of molecular motors is proposed to occur at all stages of transcription through diverse mechanisms. Further investigation of how these two motors combine their basic activities is essential to clarify the interdependent relationship between genome structure and transcription.

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.

Hydroxyl-radical footprinting combined with molecular modeling identifies unique features of DNA conformation and nucleosome positioning

Nucleosomes are the most abundant protein–DNA complexes in eukaryotes that provide compaction of genomic DNA and are implicated in regulation of transcription, DNA replication and repair. The details of DNA positioning on the nucleosome and the DNA conformation can provide key regulatory signals. Hydroxyl-radical footprinting (HRF) of protein–DNA complexes is a chemical technique that probes nucleosome organization in solution with a high precision unattainable by other methods. In this work we propose an integrative modeling method for constructing high-resolution atomistic models of nucleosomes based on HRF experiments. Our method precisely identifies DNA positioning on nucleosome by combining HRF data for both DNA strands with the pseudo-symmetry constraints. We performed high-resolution HRF for Saccharomyces cerevisiae centromeric nucleosome of unknown structure and characterized it using our integrative modeling approach. Our model provides the basis for further understanding the cooperative engagement and interplay between Cse4p protein and the A-tracts important for centromere function.

Contributions of Sequence to the Higher-Order Structures of DNA

One of the critical unanswered questions in genome biophysics is how the primary sequence of DNA bases influences the global properties of very-long-chain molecules. The local sequence-dependent features of DNA found in high-resolution structures introduce irregularities in the disposition of adjacent residues that facilitate the specific binding of proteins and modulate the global folding and interactions of double helices with hundreds of basepairs. These features also determine the positions of nucleosomes on DNA and the lengths of the interspersed DNA linkers. Like the patterns of basepair association within DNA, the arrangements of nucleosomes in chromatin modulate the properties of longer polymers. The intrachromosomal loops detected in genomic studies contain hundreds of nucleosomes, and given that the simulated configurations of chromatin depend on the lengths of linker DNA, the formation of these loops may reflect sequence-dependent information encoded within the positioning of the nucleosomes. With knowledge of the positions of nucleosomes on a given genome, methods are now at hand to estimate the looping propensities of chromatin in terms of the spacing of nucleosomes and to make a direct connection between the DNA base sequence and larger-scale chromatin folding.

Yeast HMO1: Linker Histone Reinvented

Eukaryotic genomes are packaged in chromatin. The higher-order organization of nucleosome core particles is controlled by the association of the intervening linker DNA with either the linker histone H1 or high mobility group box (HMGB) proteins. While H1 is thought to stabilize the nucleosome by preventing DNA unwrapping, the DNA bending imposed by HMGB may propagate to the nucleosome to destabilize chromatin. For metazoan H1, chromatin compaction requires its lysine-rich C-terminal domain, a domain that is buried between globular domains in the previously characterized yeast Saccharomyces cerevisiae linker histone Hho1p. Here, we discuss the functions of S. cerevisiae HMO1, an HMGB family protein unique in containing a terminal lysine-rich domain and in stabilizing genomic DNA. On ribosomal DNA (rDNA) and genes encoding ribosomal proteins, HMO1 appears to exert its role primarily by stabilizing nucleosome-free regions or "fragile" nucleosomes. During replication, HMO1 likewise appears to ensure low nucleosome density at DNA junctions associated with the DNA damage response or the need for topoisomerases to resolve catenanes. Notably, HMO1 shares with the mammalian linker histone H1 the ability to stabilize chromatin, as evidenced by the absence of HMO1 creating a more dynamic chromatin environment that is more sensitive to nuclease digestion and in which chromatin-remodeling events associated with DNA double-strand break repair occur faster; such chromatin stabilization requires the lysine-rich extension of HMO1. Thus, HMO1 appears to have evolved a unique linker histone-like function involving the ability to stabilize both conventional nucleosome arrays as well as DNA regions characterized by low nucleosome density or the presence of noncanonical nucleosomes.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

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.

Nucleosome dynamics: Sequence matters

About three quarter of all eukaryotic DNA is wrapped around protein cylinders, forming nucleosomes. Even though the histone proteins that make up the core of nucleosomes are highly conserved in evolution, nucleosomes can be very different from each other due to posttranslational modifications of the histones. Another crucial factor in making nucleosomes unique has so far been underappreciated: the sequence of their DNA. This review provides an overview of the experimental and theoretical progress that increasingly points to the importance of the nucleosomal base pair sequence. Specifically, we discuss the role of the underlying base pair sequence in nucleosome positioning, sliding, breathing, force-induced unwrapping, dissociation and partial assembly and also how the sequence can influence higher-order structures. A new view emerges: the physical properties of nucleosomes, especially their dynamical properties, are determined to a large extent by the mechanical properties of their DNA, which in turn depends on DNA sequence.

Role of transcription factor-mediated nucleosome disassembly in PHO5 gene expression

Studying nucleosome dynamics in promoter regions is crucial for understanding gene regulation. Nucleosomes regulate gene expression by sterically occluding transcription factors (TFs) and other non–histone proteins accessing genomic DNA. How the binding competition between nucleosomes and TFs leads to transcriptionally compatible promoter states is an open question. Here, we present a computational study of the nucleosome dynamics and organization in the promoter region of PHO5 gene in Saccharomyces cerevisiae. Introducing a model for nucleosome kinetics that takes into account ATP-dependent remodeling activity, DNA sequence effects, and kinetics of TFs (Pho4p), we compute the probability of obtaining different “promoter states” having different nucleosome configurations. Comparing our results with experimental data, we argue that the presence of local remodeling activity (LRA) as opposed to basal remodeling activity (BRA) is crucial in determining transcriptionally active promoter states. By modulating the LRA and Pho4p binding rate, we obtain different mRNA distributions—Poisson, bimodal, and long-tail. Through this work we explain many features of the PHO5 promoter such as sequence-dependent TF accessibility and the role of correlated dynamics between nucleosomes and TFs in opening/coverage of the TATA box. We also obtain possible ranges for TF binding rates and the magnitude of LRA.

X-ray structure of the MMTV-A nucleosome core

The conformation of DNA bound in nucleosomes depends on the DNA sequence. Questions such as how nucleosomes are positioned and how they potentially bind sequence-dependent nuclear factors require near-atomic resolution structures of the nucleosome core containing different DNA sequences; despite this, only the DNA for two similar α-satellite sequences and a sequence (601) selected in vitro have been visualized bound in the nucleosome core. Here we report the 2.6-Å resolution X-ray structure of a nucleosome core particle containing the DNA sequence of nucleosome A of the 3'-LTR of the mouse mammary tumor virus (147 bp MMTV-A). To our knowledge, this is the first nucleosome core particle structure containing a promoter sequence and crystallized from Mg(2+) ions. It reveals sequence-dependent DNA conformations not seen previously, including kinking into the DNA major groove.

Finding optimal interaction interface alignments between biological complexes

Motivation: Biological molecules perform their functions through interactions with other molecules. Structure alignment of interaction interfaces between biological complexes is an indispensable step in detecting their structural similarities, which are keys to understanding their evolutionary histories and functions. Although various structure alignment methods have been developed to successfully access the similarities of protein structures or certain types of interaction interfaces, existing alignment tools cannot directly align arbitrary types of interfaces formed by protein, DNA or RNA molecules. Specifically, they require a ‘blackbox preprocessing’ to standardize interface types and chain identifiers. Yet their performance is limited and sometimes unsatisfactory.

Results: Here we introduce a novel method, PROSTA-inter, that automatically determines and aligns interaction interfaces between two arbitrary types of complex structures. Our method uses sequentially remote fragments to search for the optimal superimposition. The optimal residue matching problem is then formulated as a maximum weighted bipartite matching problem to detect the optimal sequence order-independent alignment. Benchmark evaluation on all non-redundant protein–DNA complexes in PDB shows significant performance improvement of our method over TM-align and iAlign (with the ‘blackbox preprocessing’). Two case studies where our method discovers, for the first time, structural similarities between two pairs of functionally related protein–DNA complexes are presented. We further demonstrate the power of our method on detecting structural similarities between a protein–protein complex and a protein–RNA complex, which is biologically known as a protein–RNA mimicry case.

Availability and implementation: The PROSTA-inter web-server is publicly available at http://www.cbrc.kaust.edu.sa/prosta/.

Contact:xin.gao@kaust.edu.sa

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.

Stereochemical control of nucleosome targeting by platinum-intercalator antitumor agents

Platinum-based anticancer drugs act therapeutically by forming DNA adducts, but suffer from severe toxicity and resistance problems, which have not been overcome in spite of decades of research. And yet defined chromatin targets have generally not been considered in the drug development process. Here we designed novel platinum-intercalator species to target a highly deformed DNA site near the nucleosome center. Between two seemingly similar structural isomers, we find a striking difference in DNA site selectivity in vitro, which comes about from stereochemical constraints that limit the reactivity of the trans isomer to special DNA sequence elements while still allowing the cis isomer to efficiently form adducts at internal sites in the nucleosome core. This gives the potential for controlling nucleosome site targeting in vivo, which would engender sensitivity to epigenetic distinctions and in particular cell type/status-dependent differences in nucleosome positioning. Moreover, while both compounds yield very similar DNA-adduct structures and display antitumor cell activity rivalling that of cisplatin, the cis isomer, relative to the trans, has a much more rapid cytotoxic effect and distinct impact on cell function. The novel stereochemical principles for controlling DNA site selectivity we discovered could aid in the design of improved site discriminating agents.

Repair efficiency of clustered abasic sites by APE1 in nucleosome core particles is sequence and position dependent

DNA sequence context and nucleosomal positioning guide the repair efficiency of clustered abasic sites by APE1 enzyme.

Reading sequence-directed computational nucleosome maps

Recently developed latest version of the sequence-directed single-base resolution nucleosome mapping reveals existence of strong nucleosomes and chromatin columnar structures (columns). Broad application of this simple technique for further studies of chromatin and chromosome structure requires some basic understanding as to how it works and what information it affords. The paper provides such an introduction to the method. The oscillating maps of singular nucleosomes, of short and long oligonucleosome columns, are explained, as well as maps of chromatin on satellite DNA and occurrences of counter-phase (antiparallel) nucleosome neighbors.

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.

DNA damage may drive nucleosomal reorganization to facilitate damage detection

One issue in genome maintenance is how DNA repair proteins find lesions at rates that seem to exceed diffusion-limited search rates. We propose a phenomenon where DNA damage induces nucleosomal rearrangements which move lesions to potential rendezvous points in the chromatin structure. These rendezvous points are the dyad and the linker DNA between histones, positions in the chromatin which are more likely to be accessible by repair proteins engaged in a random search. The feasibility of this mechanism is tested by considering the statistical mechanics of DNA containing a single lesion wrapped onto the nucleosome. We consider lesions which make the DNA either more flexible or more rigid by modeling the lesion as either a decrease or an increase in the bending energy. We include this energy in a partition function model of nucleosome breathing. Our results indicate that the steady state for a breathing nucleosome will most likely position the lesion at the dyad or in the linker, depending on the energy of the lesion. A role for DNA binding proteins and chromatin remodelers is suggested based on their ability to alter the mechanical properties of the DNA and DNA-histone binding, respectively. We speculate that these positions around the nucleosome potentially serve as rendezvous points where DNA lesions may be encountered by repair proteins which may be sterically hindered from searching the rest of the nucleosomal DNA. The strength of the repositioning is strongly dependent on the structural details of the DNA lesion and the wrapping and breathing of the nucleosome. A more sophisticated evaluation of this proposed mechanism will require detailed information about breathing dynamics, the structure of partially wrapped nucleosomes, and the structural properties of damaged DNA.

Biodamage via shock waves initiated by irradiation with ions

Radiation damage following the ionising radiation of tissue has different scenarios and mechanisms depending on the projectiles or radiation modality. We investigate the radiation damage effects due to shock waves produced by ions. We analyse the strength of the shock wave capable of directly producing DNA strand breaks and, depending on the ion's linear energy transfer, estimate the radius from the ion's path, within which DNA damage by the shock wave mechanism is dominant. At much smaller values of linear energy transfer, the shock waves turn out to be instrumental in propagating reactive species formed close to the ion's path to large distances, successfully competing with diffusion.

Combined micrococcal nuclease and exonuclease III digestion reveals precise positions of the nucleosome core/linker junctions: implications for high-resolution nucleosome mapping

Micrococcal nuclease (MNase) is extensively used in genome-wide mapping of nucleosomes but its preference for AT-rich DNA leads to errors in establishing precise positions of nucleosomes. Here, we show that the MNase digestion of nucleosomes assembled on a strong nucleosome positioning sequence, Widom's clone 601, releases nucleosome cores whose sizes are strongly affected by the linker DNA sequence. Our experiments produced nucleosomal DNA sizes varying between 147 and 155 bp, with positions of the MNase cuts reflecting positions of the A⋅T pairs rather than the nucleosome core/linker junctions determined by X-ray crystallography. Extent of chromatosomal DNA protection by linker histone H1 also depends on the linker DNA sequence. Remarkably, we found that a combined treatment with MNase and exonuclease III (exoIII) overcomes MNase sequence preference producing nucleosomal DNA trimmed symmetrically and precisely at the core/linker junctions regardless of the underlying DNA sequence. We propose that combined MNase/exoIII digestion can be applied to in situ chromatin for unbiased genome-wide mapping of nucleosome positions that is not influenced by DNA sequences at the core/linker junctions. The same approach can be also used for the precise mapping of the extent of linker DNA protection by H1 and other protein factors associated with nucleosome linkers.

ATP-Dependent Chromatin Remodeling

In the eukaryotic nucleus, processes of DNA metabolism such as transcription, DNA replication, and repair occur in the context of DNA packaged into nucleosomes and higher order chromatin structures. In order to overcome the barrier presented by chromatin structures to the protein machinery carrying out these processes, the cell relies on a class of enzymes called chromatin remodeling complexes which catalyze ATP-dependent restructuring and repositioning of nucleosomes. Chromatin remodelers are large multi-subunit complexes which all share a common SF2 helicase ATPase domain in their catalytic subunit, and are classified into four different families-SWI/SNF, ISWI, CHD, INO80-based on the arrangement of other domains in their catalytic subunit as well as their non-catalytic subunit composition. A large body of structural, biochemical, and biophysical evidence suggests chromatin remodelers operate as histone octamer-anchored directional DNA translocases in order to disrupt DNA-histone interactions and catalyze nucleosome sliding. Remodeling mechanisms are family-specific and depend on factors such as how the enzyme engages with nucleosomal and linker DNA, features of DNA loop intermediates, specificity for mono- or oligonucleosomal substrates, and ability to remove histones and exchange histone variants. Ultimately, the biological function of chromatin remodelers and their genomic targeting in vivo is regulated by each complex's subunit composition, association with chromatin modifiers and histone chaperones, and affinity for chromatin signals such as histone posttranslational modifications.

Dps biomineralizing proteins: multifunctional architects of nature

Dps proteins are the structural relatives of bacterioferritins and ferritins ubiquitously present in the bacterial and archaeal kingdoms. The ball-shaped enzymes play important roles in the detoxification of ROS (reactive oxygen species), in iron scavenging to prevent Fenton reactions and in the mechanical protection of DNA. Detoxification of ROS and iron chaperoning represent the most archetypical functions of dodecameric Dps enzymes. Recent crystallographic studies of these dodecameric complexes have unravelled species-dependent mechanisms of iron uptake into the hollow spheres. Subsequent functions in iron oxidation at ferroxidase centres are highly conserved among bacteria. Final nucleation of iron as iron oxide nanoparticles has been demonstrated to originate at acidic residues located on the inner surface. Some Dps enzymes are also implicated in newly observed catalytic functions related to the formation of molecules playing roles in bacterium–host cell communication. Most recently, Dps complexes are attracting attention in semiconductor science as biomimetic tools for the technical production of the smallest metal-based quantum nanodots used in nanotechnological approaches, such as memory storage or solar cell development.

The mechanics behind DNA sequence-dependent properties of the nucleosome

Chromatin organization and composition impart sophisticated regulatory features critical to eukaryotic genomic function. Although DNA sequence-dependent histone octamer binding is important for nucleosome activity, many aspects of this phenomenon have remained elusive. We studied nucleosome structure and stability with diverse DNA sequences, including Widom 601 derivatives with the highest known octamer affinities, to establish a simple model behind the mechanics of sequence dependency. This uncovers the unique but unexpected role of TA dinucleotides and a propensity for G|C-rich sequence elements to conform energetically favourably at most locations around the histone octamer, which rationalizes G|C% as the most predictive factor for nucleosome occupancy in vivo. In addition, our findings reveal dominant constraints on double helix conformation by H3–H4 relative to H2A–H2B binding and DNA sequence context-dependency underlying nucleosome structure, positioning and stability. This provides a basis for improved prediction of nucleosomal properties and the design of tailored DNA constructs for chromatin investigations.

On the sequence-directed nature of human gene mutation: the role of genomic architecture and the local DNA sequence environment in mediating gene mutations underlying human inherited disease

Different types of human gene mutation may vary in size, from structural variants (SVs) to single base-pair substitutions, but what they all have in common is that their nature, size and location are often determined either by specific characteristics of the local DNA sequence environment or by higher order features of the genomic architecture. The human genome is now recognized to contain "pervasive architectural flaws" in that certain DNA sequences are inherently mutation prone by virtue of their base composition, sequence repetitivity and/or epigenetic modification. Here, we explore how the nature, location and frequency of different types of mutation causing inherited disease are shaped in large part, and often in remarkably predictable ways, by the local DNA sequence environment. The mutability of a given gene or genomic region may also be influenced indirectly by a variety of noncanonical (non-B) secondary structures whose formation is facilitated by the underlying DNA sequence. Since these non-B DNA structures can interfere with subsequent DNA replication and repair and may serve to increase mutation frequencies in generalized fashion (i.e., both in the context of subtle mutations and SVs), they have the potential to serve as a unifying concept in studies of mutational mechanisms underlying human inherited disease.

Specific DNA structural attributes modulate platinum anticancer drug site selection and cross-link generation

Heavy metal compounds have toxic and medicinal potential through capacity to form strong specific bonds with macromolecules, and the interaction of platinum drugs at the major groove nitrogen atom of guanine bases primarily underlies their therapeutic activity. By crystallographic analysis of transition metal-and in particular platinum compound-DNA site selectivity in the nucleosome core, we establish that steric accessibility, which is controlled by specific structural parameters of the double helix, modulates initial guanine-metal bond formation. Moreover, DNA conformational features can be linked to both similarities and distinctions in platinum drug adduct formation between the naked and nucleosomal DNA states. Notably, structures that facilitate initial platinum-guanine bond formation can oppose cross-link generation, rationalizing the occurrence of long-lived therapeutically ineffective monofunctional adducts. These findings illuminate DNA structure-dependent reactivity and provide a novel framework for understanding metal-double helix interactions, which should facilitate the development of improved chromatin-targeting medicinal agents.

Modeling DNA-bending in the nucleosome: role of AA periodicity

This paper uses atomistic molecular mechanics within the framework of the JUMNA model to study the bending properties of DNA segments, with emphasis on understanding the role of the 10 bp periodicity associated with AA repeats that has been found to dominate in nucleosomal DNA. The calculations impose a bending potential on 18 bp segments that is consistent with nucleosome structures (i.e., radius of curvature of 4.1 nm), and then determine the energies of the minimum energy structures for different values of the rotational register (a measure of the direction of bending of the DNA) subject to forces derived from the Amber force field (parm99bsc0). The results show that sequences that contain the 10 bp repeats but are otherwise random have a narrow distribution of rotational register values that minimize the energy such that it is possible to combine several minimized structures to give the 147 bp nearly planar loop structure of the nucleosome. The rotational register values that lead to minimum bending energy with 10 bp AA repeats have a narrower minor groove, which points toward the histone interior at the positions of the AA repeats, which is a result that matches the experiments. The calculations also show that these sequences have a relatively flat potential energy landscape for bending to a 4.1 nm radius of curvature. Random sequences that do not have the 10 bp AA repeats have less stable bent structures, and a flat rotational register distribution, such that low energy nearly planar loops are less likely.

Working the kinks out of nucleosomal DNA

Condensation of DNA in the nucleosome takes advantage of its double-helical architecture. The DNA deforms at sites where the base pairs face the histone octamer. The largest so-called kink-and-slide deformations occur in the vicinity of arginines that penetrate the minor groove. Nucleosome structures formed from the 601 positioning sequence differ subtly from those incorporating an AT-rich human α-satellite DNA. Restraints imposed by the histone arginines on the displacement of base pairs can modulate the sequence-dependent deformability of DNA and potentially contribute to the unique features of the different nucleosomes. Steric barriers mimicking constraints found in the nucleosome induce the simulated large-scale rearrangement of canonical B DNA to kink-and-slide states. The pathway to these states shows nonharmonic behavior consistent with bending profiles inferred from AFM measurements.

Nucleosome structural studies

Chromatin plays a fundamental role in eukaryotic genomic regulation, and the increasing awareness of the importance of epigenetic processes in human health and disease emphasizes the need for understanding the structure and function of the nucleosome. Recent advances in chromatin structural studies, including the first structures of nucleosomes containing the Widom 601 sequence and the structure of a chromatin protein-nucleosome assembly, have provided new insight into stretching of nucleosomal DNA, nucleosome positioning, binding of metal ions, drugs and therapeutic candidates to nucleosomes, and nucleosome recognition by nuclear proteins. These discoveries ensure promising future prospects for unravelling structural attributes of chromatin.

The effects of histone H4 tail acetylations on cation-induced chromatin folding and self-association

Understanding the molecular mechanisms behind regulation of chromatin folding through covalent modifications of the histone N-terminal tails is hampered by a lack of accessible chromatin containing precisely modified histones. We study the internal folding and intermolecular self-association of a chromatin system consisting of saturated 12-mer nucleosome arrays containing various combinations of completely acetylated lysines at positions 5, 8, 12 and 16 of histone H4, induced by the cations Na⁺, K⁺, Mg²⁺, Ca²⁺, cobalt-hexammine³⁺, spermidine³⁺ and spermine⁴⁺. Histones were prepared using a novel semi-synthetic approach with native chemical ligation. Acetylation of H4-K16, but not its glutamine mutation, drastically reduces cation-induced folding of the array. Neither acetylations nor mutations of all the sites K5, K8 and K12 can induce a similar degree of array unfolding. The ubiquitous K⁺, (as well as Rb⁺ and Cs⁺) showed an unfolding effect on unmodified arrays almost similar to that of H4-K16 acetylation. We propose that K⁺ (and Rb⁺/Cs⁺) binding to a site on the H2B histone (R96-L99) disrupts H4K16 ε-amino group binding to this specific site, thereby deranging H4 tail-mediated nucleosome–nucleosome stacking and that a similar mechanism operates in the case of H4-K16 acetylation. Inter-array self-association follows electrostatic behavior and is largely insensitive to the position or nature of the H4 tail charge modification.

Crystal structures of nucleosome core particles containing the '601' strong positioning sequence

Nucleosome positioning plays a key role in genomic regulation by defining histone-DNA context and by modulating access to specific sites. Moreover, the histone-DNA register influences the double-helix structure, which in turn can affect the association of small molecules and protein factors. Analysis of genomic and synthetic DNA has revealed sequence motifs that direct nucleosome positioning in vitro; thus, establishing the basis for the DNA sequence dependence of positioning would shed light on the mechanics of the double helix and its contribution to chromatin structure in vivo. However, acquisition of well-diffracting nucleosome core particle (NCP) crystals is extremely dependent on the DNA fragment used for assembly, and all previous NCP crystal structures have been based on human α-satellite sequences. Here, we describe the crystal structures of Xenopus NCPs containing one of the strongest known histone octamer binding and positioning sequences, the so-called '601' DNA. Two distinct 145-bp 601 crystal forms display the same histone-DNA register, which coincides with the occurrence of DNA stretching-overtwisting in both halves of the particle around five double-helical turns from the nucleosome center, giving the DNA an 'effective length' of 147 bp. As we have found previously with stretching around two turns from the nucleosome center for a centromere-based sequence, the terminal stretching observed in the 601 constructs is associated with extreme kinking into the minor groove at purine-purine (pyrimidine-pyrimidine) dinucleotide steps. In other contexts, these step types display an overall nonflexible behavior, which raises the possibility that DNA stretching in the nucleosome or extreme distortions in general have unique sequence dependency characteristics. Our findings indicate that DNA stretching is an intrinsically predisposed site-specific property of the nucleosome and suggest how NCP crystal structures with diverse DNA sequences can be obtained.

Perturbations in nucleosome structure from heavy metal association

Heavy metals have the potential to engage in strong bonding interactions and can thus function in essential as well as toxic or therapeutic capacities. We conducted crystallographic analyses of heavy cation binding to the nucleosome core particle and found that Co²⁺ and Ni²⁺ preferentially associate with the DNA major groove, in a sequence- and conformation-dependent manner. Conversely, Rb⁺ and Cs⁺ are found to bind only opportunistically to minor groove elements of the DNA, in particular at narrow AT dinucleotide sites. Furthermore, relative to Mn²⁺ the aggressive coordination of Co²⁺ and Ni²⁺ to guanine bases is observed to induce a shift in histone–DNA register around the nucleosome center by stabilizing DNA stretching over one region accompanied by expulsion of two bases at an opposing location. These ‘softer’ transition metals also associate with multiple histone protein sites, including inter-nucleosomal cross-linking, and display a proclivity for coordination to histidine. Sustained binding and the ability to induce structural perturbations at specific locations in the nucleosome may contribute to genetic and epigenetic mechanisms of carcinogenesis mediated by Co²⁺ and Ni²⁺.