1
|
Klempahn S, Schiessel H, Blossey R. Chromatin remodelers: a concise introduction for biophysicists. Biophys Rev 2024; 16:357-363. [PMID: 39099840 PMCID: PMC11296983 DOI: 10.1007/s12551-024-01199-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2024] [Accepted: 05/16/2024] [Indexed: 08/06/2024] Open
Abstract
Chromatin remodelers are molecular motors that act on nucleosomes: they move them along DNA or (dis-)assemble them. Despite the fact that they perform essential regulatory functions in cells-their deregulation can contribute to the development of cancers and lead to cell death-chromatin remodelers have only received meager attention in the biophysics community so far. In this short text, we attempt to present the key features of this interesting class of enzymes obtained with different experimental and theoretical methods, thereby providing a concise introduction for biophysicists to further stimulate interest in their properties.
Collapse
Affiliation(s)
- Sophie Klempahn
- Cluster of Excellence Physics of Life, TUD Dresden University of Technology, 01307 Dresden, Germany
| | - Helmut Schiessel
- Cluster of Excellence Physics of Life, TUD Dresden University of Technology, 01307 Dresden, Germany
- Institut für Theoretische Physik, TUD Dresden University of Technology, 01069 Dresden, Germany
| | - Ralf Blossey
- University of Lille, Unité de Glycobiologie Structurale et Fonctionnelle (UGSF), CNRS UMR8576, 59000 Lille, France
| |
Collapse
|
2
|
Chakraborty D, Mondal B, Thirumalai D. Brewing COFFEE: A Sequence-Specific Coarse-Grained Energy Function for Simulations of DNA-Protein Complexes. J Chem Theory Comput 2024; 20:1398-1413. [PMID: 38241144 DOI: 10.1021/acs.jctc.3c00833] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2024]
Abstract
DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription and gene expression to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse-grained Force Field for Energy Estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the self-organized polymer model with side-chains for proteins and the three interaction site model for DNA in a modular fashion, without recalibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a data set of high-resolution crystal structures. The only parameter in COFFEE is the strength (λDNAPRO) of the DNA-protein contact potential. For an optimal choice of λDNAPRO, the crystallographic B-factors for DNA-protein complexes with varying sizes and topologies are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts scattering profiles that are in quantitative agreement with small-angle X-ray scattering experiments, as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which do not alter the balance of electrostatic interactions but affect chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.
Collapse
Affiliation(s)
- Debayan Chakraborty
- Department of Chemistry, The University of Texas at Austin, 105 E 24th Street, Stop A5300, Austin 78712, Texas, United States
| | - Balaka Mondal
- Department of Chemistry, The University of Texas at Austin, 105 E 24th Street, Stop A5300, Austin 78712, Texas, United States
| | - D Thirumalai
- Department of Chemistry, The University of Texas at Austin, 105 E 24th Street, Stop A5300, Austin 78712, Texas, United States
- Department of Physics, The University of Texas at Austin, 2515 Speedway, Austin 78712, Texas, United States
| |
Collapse
|
3
|
Li X, Chou T. Stochastic nucleosome disassembly mediated by remodelers and histone fragmentation. J Chem Phys 2023; 159:204107. [PMID: 38010331 PMCID: PMC10684310 DOI: 10.1063/5.0165136] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Accepted: 10/14/2023] [Indexed: 11/29/2023] Open
Abstract
We construct and analyze monomeric and multimeric models of the stochastic disassembly of a single nucleosome. Our monomeric model predicts the time needed for a number of histone-DNA contacts to spontaneously break, leading to dissociation of a non-fragmented histone from DNA. The dissociation process can be facilitated by DNA binding proteins or processing molecular motors that compete with histones for histone-DNA contact sites. Eigenvalue analysis of the corresponding master equation allows us to evaluate histone detachment times under both spontaneous detachment and protein-facilitated processes. We find that competitive DNA binding of remodeling proteins can significantly reduce the typical detachment time but only if these remodelers have DNA-binding affinities comparable to those of histone-DNA contact sites. In the presence of processive motors, the histone detachment rate is shown to be proportional to the product of the histone single-bond dissociation constant and the speed of motor protein procession. Our simple intact-histone model is then extended to allow for multimeric nucleosome kinetics that reveal additional pathways of disassembly. In addition to a dependence of complete disassembly times on subunit-DNA contact energies, we show how histone subunit concentrations in bulk solutions can mediate the disassembly process by rescuing partially disassembled nucleosomes. Moreover, our kinetic model predicts that remodeler binding can also bias certain pathways of nucleosome disassembly, with higher remodeler binding rates favoring intact-histone detachment.
Collapse
Affiliation(s)
- Xiangting Li
- Department of Computational Medicine, University of California, Los Angeles, California 90095-1766, USA
| | - Tom Chou
- Author to whom correspondence should be addressed:
| |
Collapse
|
4
|
Abstract
Nearly three-fourths of all eukaryotic DNA is occupied by nucleosomes, protein-DNA complexes comprising octameric histone core proteins and ∼150 base pairs of DNA. In addition to acting as a DNA compaction vehicle, the dynamics of nucleosomes regulate the DNA site accessibility for the nonhistone proteins, thereby controlling regulatory processes involved in determining the cell identity and cell fate. Here, we propose an analytical framework to analyze the role of nucleosome dynamics on the target search process of transcription factors through a simple discrete-state stochastic description of the search process. By considering the experimentally determined kinetic rates associated with protein and nucleosome dynamics as the only inputs, we estimate the target search time of a protein via first-passage probability calculations separately during nucleosome breathing and sliding dynamics. Although both the nucleosome dynamics permit transient access to the DNA sites that are otherwise occluded by the histone proteins, our result suggests substantial differences between the protein search mechanism on a nucleosome performing breathing and sliding dynamics. Furthermore, we identify the molecular factors that influence the search efficiency and demonstrate how these factors together portray a highly dynamic landscape of gene regulation. Our analytical results are validated using extensive Monte Carlo simulations.
Collapse
Affiliation(s)
- Sujeet Kumar Mishra
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Arnab Bhattacherjee
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| |
Collapse
|
5
|
Onufriev AV. Biologically relevant small variations of intra-cellular pH can have significant effect on stability of protein-DNA complexes, including the nucleosome. Front Mol Biosci 2023; 10:1067787. [PMID: 37143824 PMCID: PMC10151541 DOI: 10.3389/fmolb.2023.1067787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2022] [Accepted: 01/16/2023] [Indexed: 05/06/2023] Open
Abstract
Stability of a protein-ligand complex may be sensitive to pH of its environment. Here we explore, computationally, stability of a set of protein-nucleic acid complexes using fundamental thermodynamic linkage relationship. The nucleosome, as well as an essentially random selection of 20 protein complexes with DNA or RNA, are included in the analysis. An increase in intra-cellular/intra-nuclear pH destabilizes most complexes, including the nucleosome. We propose to quantify the effect by ΔΔG0.3-the change in the binding free energy due to pH increase of 0.3 units, corresponding to doubling of the H + activity; variations of pH of this amplitude can occur in living cells, including in the course of the cell cycle, and in cancer cells relative to normal ones. We suggest, based on relevant experimental findings, a threshold of biological significance of 1 2 k B T ( ∼ 0.3 k c a l / m o l ) for changes of stability of chromatin-related protein-DNA complexes: a change in the binding affinity above the threshold may have biological consequences. We find that for 70% of the examined complexes, Δ Δ G 0.3 > 1 2 k B T (for 10%, ΔΔG0.3 is between 3 and 4 k B T). Thus, small but relevant variations of intra-nuclear pH of 0.3 may have biological consequences for many protein-nucleic acid complexes. The binding affinity between the histone octamer and its DNA, which directly affects the DNA accessibility in the nucleosome, is predicted to be highly sensitive to intra-nuclear pH. A variation of 0.3 units results in ΔΔG0.3 ∼ 10k B T ( ∼ 6 k c a l / m o l ) ; for spontaneous unwrapping of 20 bp long entry/exit fragments of the nucleosomal DNA, ΔΔG0.3 = 2.2k B T; partial disassembly of the nucleosome into the tetrasome is characterized by ΔΔG0.3 = 5.2k B T. The predicted pH -induced modulations of the nucleosome stability are significant enough to suggest that they may have consequences relevant to the biological function of the nucleosome. Accessibility of the nucleosomal DNA is predicted to positively correlate with pH variations during the cell cycle; an increase in intra-cellular pH seen in cancer cells is predicted to lead to a more accessible nucleosomal DNA; a drop in pH associated with apoptosis is predicted to make nucleosomal DNA less accessible. We speculate that processes that depend on accessibility to the DNA in the nucleosomes, such as transcription or DNA replication, might become upregulated due to relatively small, but nevertheless realistic increases of intra-nuclear pH.
Collapse
Affiliation(s)
- Alexey V. Onufriev
- Department of Physics, Virginia Tech, Blacksburg, Blacksburg, VA, United States
- Department of Computer Science, Virginia Tech, Blacksburg, Blacksburg, VA, United States
- Center from Soft Matter and Biological Physics, Virginia Tech, Blacksburg, VA, United States
- *Correspondence: Alexey V. Onufriev,
| |
Collapse
|
6
|
Zuiddam M, Shakiba B, Schiessel H. Multiplexing mechanical and translational cues on genes. Biophys J 2022; 121:4311-4324. [PMID: 36230003 PMCID: PMC9703045 DOI: 10.1016/j.bpj.2022.10.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2022] [Revised: 07/06/2022] [Accepted: 10/07/2022] [Indexed: 12/14/2022] Open
Abstract
The genetic code gives precise instructions on how to translate codons into amino acids. Due to the degeneracy of the genetic code-18 out of 20 amino acids are encoded for by more than one codon-more information can be stored in a basepair sequence. Indeed, various types of additional information have been discussed in the literature, e.g., the positioning of nucleosomes along eukaryotic genomes and the modulation of the translating efficiency in ribosomes to influence cotranslational protein folding. The purpose of this study is to show that it is indeed possible to carry more than one additional layer of information on top of a gene. In particular, we show how much translation efficiency and nucleosome positioning can be adjusted simultaneously without changing the encoded protein. We achieve this by mapping genes on weighted graphs that contain all synonymous genes, and then finding shortest paths through these graphs. This enables us, for example, to readjust the disrupted translational efficiency profile after a gene has been introduced from one organism (e.g., human) into another (e.g., yeast) without greatly changing the nucleosome landscape intrinsically encoded by the DNA molecule.
Collapse
Affiliation(s)
- Martijn Zuiddam
- Institute Lorentz for Theoretical Physics, Leiden University, Leiden, the Netherlands
| | - Bahareh Shakiba
- Institute Lorentz for Theoretical Physics, Leiden University, Leiden, the Netherlands
| | - Helmut Schiessel
- Cluster of Excellence Physics of Life, TU Dresden, Dresden, Germany.
| |
Collapse
|
7
|
MacCarthy CM, Huertas J, Ortmeier C, Vom Bruch H, Tan DS, Reinke D, Sander A, Bergbrede T, Jauch R, Schöler HR, Cojocaru V. OCT4 interprets and enhances nucleosome flexibility. Nucleic Acids Res 2022; 50:10311-10327. [PMID: 36130732 PMCID: PMC9561370 DOI: 10.1093/nar/gkac755] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2022] [Revised: 07/28/2022] [Accepted: 08/24/2022] [Indexed: 01/06/2023] Open
Abstract
Pioneer transcription factors are proteins that induce cellular identity transitions by binding to inaccessible regions of DNA in nuclear chromatin. They contribute to chromatin opening and recruit other factors to regulatory DNA elements. The structural features and dynamics modulating their interaction with nucleosomes are still unresolved. From a combination of experiments and molecular simulations, we reveal here how the pioneer factor and master regulator of pluripotency, Oct4, interprets and enhances nucleosome structural flexibility. The magnitude of Oct4’s impact on nucleosome dynamics depends on the binding site position and the mobility of the unstructured tails of nucleosomal histone proteins. Oct4 uses both its DNA binding domains to propagate and stabilize open nucleosome conformations, one for specific sequence recognition and the other for nonspecific interactions with nearby regions of DNA. Our findings provide a structural basis for the versatility of transcription factors in engaging with nucleosomes and have implications for understanding how pioneer factors induce chromatin dynamics.
Collapse
Affiliation(s)
- Caitlin M MacCarthy
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany
| | - Jan Huertas
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany.,Yusuf Hamied Department of Chemistry, University of Cambridge, UK
| | - Claudia Ortmeier
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany
| | - Hermann Vom Bruch
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany
| | - Daisylyn Senna Tan
- School of Biomedical Sciences, Li Ka Shing Faculty of Medicine, University of Hong Kong, Hong Kong SAR, China
| | - Deike Reinke
- Max Planck Institute of Molecular Physiology, Dortmund, Germany
| | - Astrid Sander
- Max Planck Institute of Molecular Physiology, Dortmund, Germany
| | | | - Ralf Jauch
- School of Biomedical Sciences, Li Ka Shing Faculty of Medicine, University of Hong Kong, Hong Kong SAR, China
| | - Hans R Schöler
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany.,Medical Faculty, University of Münster, Germany
| | - Vlad Cojocaru
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany.,Computational Structural Biology Group, University of Utrecht, The Netherlands.,STAR-UBB Institute, Babeş-Bolyai University, Cluj-Napoca, Romania
| |
Collapse
|
8
|
Gurova K. Can aggressive cancers be identified by the "aggressiveness" of their chromatin? Bioessays 2022; 44:e2100212. [PMID: 35452144 DOI: 10.1002/bies.202100212] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2021] [Revised: 04/08/2022] [Accepted: 04/12/2022] [Indexed: 12/15/2022]
Abstract
Phenotypic plasticity is a crucial feature of aggressive cancer, providing the means for cancer progression. Stochastic changes in tumor cell transcriptional programs increase the chances of survival under any condition. I hypothesize that unstable chromatin permits stochastic transitions between transcriptional programs in aggressive cancers and supports non-genetic heterogeneity of tumor cells as a basis for their adaptability. I present a mechanistic model for unstable chromatin which includes destabilized nucleosomes, mobile chromatin fibers and random enhancer-promoter contacts, resulting in stochastic transcription. I suggest potential markers for "unsettled" chromatin in tumors associated with poor prognosis. Although many of the characteristics of unstable chromatin have been described, they were mostly used to explain changes in the transcription of individual genes. I discuss approaches to evaluate the role of unstable chromatin in non-genetic tumor cell heterogeneity and suggest using the degree of chromatin instability and transcriptional noise in tumor cells to predict cancer aggressiveness.
Collapse
Affiliation(s)
- Katerina Gurova
- Department of Cell Stress Biology, Roswell Park Comprehensive Cancer Center, Buffalo, New York, USA
| |
Collapse
|
9
|
Mines RC, Lipniacki T, Shen X. Slow nucleosome dynamics set the transcriptional speed limit and induce RNA polymerase II traffic jams and bursts. PLoS Comput Biol 2022; 18:e1009811. [PMID: 35143483 PMCID: PMC8865691 DOI: 10.1371/journal.pcbi.1009811] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Revised: 02/23/2022] [Accepted: 01/06/2022] [Indexed: 11/19/2022] Open
Abstract
Nucleosomes are recognized as key regulators of transcription. However, the relationship between slow nucleosome unwrapping dynamics and bulk transcriptional properties has not been thoroughly explored. Here, an agent-based model that we call the dynamic defect Totally Asymmetric Simple Exclusion Process (ddTASEP) was constructed to investigate the effects of nucleosome-induced pausing on transcriptional dynamics. Pausing due to slow nucleosome dynamics induced RNAPII convoy formation, which would cooperatively prevent nucleosome rebinding leading to bursts of transcription. The mean first passage time (MFPT) and the variance of first passage time (VFPT) were analytically expressed in terms of the nucleosome rate constants, allowing for the direct quantification of the effects of nucleosome-induced pausing on pioneering polymerase dynamics. The mean first passage elongation rate γ(hc, ho) is inversely proportional to the MFPT and can be considered to be a new axis of the ddTASEP phase diagram, orthogonal to the classical αβ-plane (where α and β are the initiation and termination rates). Subsequently, we showed that, for β = 1, there is a novel jamming transition in the αγ-plane that separates the ddTASEP dynamics into initiation-limited and nucleosome pausing-limited regions. We propose analytical estimates for the RNAPII density ρ, average elongation rate v, and transcription flux J and verified them numerically. We demonstrate that the intra-burst RNAPII waiting times tin follow the time-headway distribution of a max flux TASEP and that the average inter-burst interval tIBI¯ correlates with the index of dispersion De. In the limit γ→0, the average burst size reaches a maximum set by the closing rate hc. When α≪1, the burst sizes are geometrically distributed, allowing large bursts even while the average burst size NB¯ is small. Last, preliminary results on the relative effects of static and dynamic defects are presented to show that dynamic defects can induce equal or greater pausing than static bottle necks. To perform specific functions, cells must express specific genes by copying the information in DNA into RNA via transcription. Structural proteins called nucleosomes are spaced every 200 base pairs along the length of a strand of DNA and play a crucial function in the regulation of gene activity by tightly binding DNA strands and condensing them into heterochromatin, preventing transcription by RNA polymerase II (RNAPII). Even on active genes where nucleosomes are loosely attached to DNA strands, the wrapping and unwrapping of nucleosomes pause transcription as RNAPII passes by. Previous mathematical models of transcription have compared this biological process to traffic on a one lane highway without obstructions. In contrast, our proposed model simulates transcription like traffic in a grid system where nucleosomes can be thought of as pedestrians or other vehicles crossing the road at regularly spaced intersections. Just as side street traffic and pedestrian crossings can cause cars to form convoys and cause jams limiting the max speed in an area, nucleosomes can cause RNAPII to form convoys that lead to bursts of mRNA production and limit the average polymerase flux through the gene.
Collapse
Affiliation(s)
- Robert C. Mines
- Department of Biomedical Engineering, Duke University, Durham, North Carolina, United States of America
| | - Tomasz Lipniacki
- Department of Biosystems and Soft Matter, Institute of Fundamental Technological Research, Polish Academy of Sciences, Warsaw, Poland
- * E-mail: (TL); (XS)
| | - Xiling Shen
- Department of Biomedical Engineering, Duke University, Durham, North Carolina, United States of America
- Center for Genomic and Computational Biology, Duke University, Durham, North Carolina, United States of America
- Woo Center for Big Data and Precision Health, Duke University, Durham, North Carolina, United States of America
- * E-mail: (TL); (XS)
| |
Collapse
|
10
|
Dohnalová H, Lankaš F. Deciphering the mechanical properties of
B‐DNA
duplex. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1575] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Hana Dohnalová
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| | - Filip Lankaš
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| |
Collapse
|
11
|
Brandani GB, Tan C, Takada S. The kinetic landscape of nucleosome assembly: A coarse-grained molecular dynamics study. PLoS Comput Biol 2021; 17:e1009253. [PMID: 34314440 PMCID: PMC8345847 DOI: 10.1371/journal.pcbi.1009253] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 08/06/2021] [Accepted: 07/06/2021] [Indexed: 11/23/2022] Open
Abstract
The organization of nucleosomes along the Eukaryotic genome is maintained over time despite disruptive events such as replication. During this complex process, histones and DNA can form a variety of non-canonical nucleosome conformations, but their precise molecular details and roles during nucleosome assembly remain unclear. In this study, employing coarse-grained molecular dynamics simulations and Markov state modeling, we characterized the complete kinetics of nucleosome assembly. On the nucleosome-positioning 601 DNA sequence, we observe a rich transition network among various canonical and non-canonical tetrasome, hexasome, and nucleosome conformations. A low salt environment makes nucleosomes stable, but the kinetic landscape becomes more rugged, so that the system is more likely to be trapped in off-pathway partially assembled intermediates. Finally, we find that the co-operativity between DNA bending and histone association enables positioning sequence motifs to direct the assembly process, with potential implications for the dynamic organization of nucleosomes on real genomic sequences. Nucleosomes are biomolecular complexes formed by DNA wrapped around histone proteins. They represent the basic units of Eukaryotic chromosomes, compacting the genome so that it fits into the small nucleus, and regulating important biological processes such as gene expression. Nucleosomes are disassembled during disruptive events such as DNA replication, and re-assembled afterwards to preserve the correct organization of chromatin. However, the molecular details of nucleosome assembly are still not well understood. In particular, experiments found that histones and DNA may associate into a variety of non-canonical complexes, but their precise conformation and role during assembly remain unclear. In this study, we addressed these problems by performing extensive molecular dynamics simulations of nucleosomes undergoing assembly and disassembly. The simulations reveal many insights into the kinetics of assembly, the structure of non-canonical nucleosome intermediates, and the influence of salt concentration and DNA sequence on the assembly process.
Collapse
Affiliation(s)
- Giovanni B. Brandani
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
- * E-mail: (GBB); (ST)
| | - Cheng Tan
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
- * E-mail: (GBB); (ST)
| |
Collapse
|
12
|
Abstract
DNA dynamics can only be understood by taking into account its complex mechanical behavior at different length scales. At the micrometer level, the mechanical properties of single DNA molecules have been well-characterized by polymer models and are commonly quantified by a persistence length of 50 nm (~150 bp). However, at the base pair level (~3.4 Å), the dynamics of DNA involves complex molecular mechanisms that are still being deciphered. Here, we review recent single-molecule experiments and molecular dynamics simulations that are providing novel insights into DNA mechanics from such a molecular perspective. We first discuss recent findings on sequence-dependent DNA mechanical properties, including sequences that resist mechanical stress and sequences that can accommodate strong deformations. We then comment on the intricate effects of cytosine methylation and DNA mismatches on DNA mechanics. Finally, we review recently reported differences in the mechanical properties of DNA and double-stranded RNA, the other double-helical carrier of genetic information. A thorough examination of the recent single-molecule literature permits establishing a set of general 'rules' that reasonably explain the mechanics of nucleic acids at the base pair level. These simple rules offer an improved description of certain biological systems and might serve as valuable guidelines for future design of DNA and RNA nanostructures.
Collapse
|
13
|
Huertas J, Schöler HR, Cojocaru V. Histone tails cooperate to control the breathing of genomic nucleosomes. PLoS Comput Biol 2021; 17:e1009013. [PMID: 34081696 PMCID: PMC8174689 DOI: 10.1371/journal.pcbi.1009013] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2020] [Accepted: 04/28/2021] [Indexed: 11/19/2022] Open
Abstract
Genomic DNA is packaged in chromatin, a dynamic fiber variable in size and compaction. In chromatin, repeating nucleosome units wrap 145–147 DNA basepairs around histone proteins. Genetic and epigenetic regulation of genes relies on structural transitions in chromatin which are driven by intra- and inter-nucleosome dynamics and modulated by chemical modifications of the unstructured terminal tails of histones. Here we demonstrate how the interplay between histone H3 and H2A tails control ample nucleosome breathing motions. We monitored large openings of two genomic nucleosomes, and only moderate breathing of an engineered nucleosome in atomistic molecular simulations amounting to 24 μs. Transitions between open and closed nucleosome conformations were mediated by the displacement and changes in compaction of the two histone tails. These motions involved changes in the DNA interaction profiles of clusters of epigenetic regulatory aminoacids in the tails. Removing the histone tails resulted in a large increase of the amplitude of nucleosome breathing but did not change the sequence dependent pattern of the motions. Histone tail modulated nucleosome breathing is a key mechanism of chromatin dynamics with important implications for epigenetic regulation. In the cell, the DNA is packed in chromatin. Chromatin is a highly dynamic fiber structure made of arrays of nucleosomes with different degrees of compaction. Each nucleosome has 145–147 basepairs of DNA wrapped around a protein octamer made of four unique histone proteins. Each histone is present twice and has a structured part and one or two disordered terminal tails. The regulation of gene expression in the cell and during cellular transitions depends on dynamic changes in chromatin structure. Chromatin dynamics are modulated by intra and inter nucleosome motions and by posttranslational chemical modifications of the histone tails. Here we reveal how histone tails control the intra nucleosome dynamics at atomic resolution. From extensive sampling of nucleosome dynamics in atomistic molecular simulations, we show that genomic nucleosomes breath more extensively than engineered ones and we describe how two histone tails cooperate to control nucleosome breathing through interactions between clusters of positively charged residues and the DNA. Nucleosome conformations with different degrees of opening are associated with different conformations, positions, and DNA interaction patterns of the tails. With this mechanism, we contribute to the understanding of chromatin dynamics at atomic resolution.
Collapse
Affiliation(s)
- Jan Huertas
- In Silico Biomolecular Structure and Dynamics Group, Hubrecht Institute, Utrecht, The Netherlands
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany
- Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster, Germany
| | - Hans Robert Schöler
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany
- Medical Faculty, University of Münster, Münster, Germany
| | - Vlad Cojocaru
- In Silico Biomolecular Structure and Dynamics Group, Hubrecht Institute, Utrecht, The Netherlands
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany
- Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster, Germany
- * E-mail: ,
| |
Collapse
|
14
|
Huertas J, Cojocaru V. Breaths, Twists, and Turns of Atomistic Nucleosomes. J Mol Biol 2020; 433:166744. [PMID: 33309853 DOI: 10.1016/j.jmb.2020.166744] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Revised: 12/02/2020] [Accepted: 12/02/2020] [Indexed: 02/07/2023]
Abstract
Gene regulation programs establish cellular identity and rely on dynamic changes in the structural packaging of genomic DNA. The DNA is packaged in chromatin, which is formed from arrays of nucleosomes displaying different degree of compaction and different lengths of inter-nucleosomal linker DNA. The nucleosome represents the repetitive unit of chromatin and is formed by wrapping 145-147 basepairs of DNA around an octamer of histone proteins. Each of the four histones is present twice and has a structured core and intrinsically disordered terminal tails. Chromatin dynamics are triggered by inter- and intra-nucleosome motions that are controlled by the DNA sequence, the interactions between the histone core and the DNA, and the conformations, positions, and DNA interactions of the histone tails. Understanding chromatin dynamics requires studying all these features at the highest possible resolution. For this, molecular dynamics simulations can be used as a powerful complement or alternative to experimental approaches, from which it is often very challenging to characterize the structural features and atomic interactions controlling nucleosome motions. Molecular dynamics simulations can be performed at different resolutions, by coarse graining the molecular system with varying levels of details. Here we review the successes and the remaining challenges of the application of atomic resolution simulations to study the structure and dynamics of nucleosomes and their complexes with interacting partners.
Collapse
Affiliation(s)
- Jan Huertas
- In Silico Biomolecular Structure and Dynamics Group, Hubrecht Institute, Utrecht, the Netherlands; Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany; Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster, Germany
| | - Vlad Cojocaru
- In Silico Biomolecular Structure and Dynamics Group, Hubrecht Institute, Utrecht, the Netherlands; Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany; Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster, Germany.
| |
Collapse
|
15
|
Eck E, Liu J, Kazemzadeh-Atoufi M, Ghoreishi S, Blythe SA, Garcia HG. Quantitative dissection of transcription in development yields evidence for transcription-factor-driven chromatin accessibility. eLife 2020; 9:e56429. [PMID: 33074101 PMCID: PMC7738189 DOI: 10.7554/elife.56429] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2020] [Accepted: 10/16/2020] [Indexed: 12/28/2022] Open
Abstract
Thermodynamic models of gene regulation can predict transcriptional regulation in bacteria, but in eukaryotes, chromatin accessibility and energy expenditure may call for a different framework. Here, we systematically tested the predictive power of models of DNA accessibility based on the Monod-Wyman-Changeux (MWC) model of allostery, which posits that chromatin fluctuates between accessible and inaccessible states. We dissected the regulatory dynamics of hunchback by the activator Bicoid and the pioneer-like transcription factor Zelda in living Drosophila embryos and showed that no thermodynamic or non-equilibrium MWC model can recapitulate hunchback transcription. Therefore, we explored a model where DNA accessibility is not the result of thermal fluctuations but is catalyzed by Bicoid and Zelda, possibly through histone acetylation, and found that this model can predict hunchback dynamics. Thus, our theory-experiment dialogue uncovered potential molecular mechanisms of transcriptional regulatory dynamics, a key step toward reaching a predictive understanding of developmental decision-making.
Collapse
Affiliation(s)
- Elizabeth Eck
- Biophysics Graduate Group, University of California at BerkeleyBerkeleyUnited States
| | - Jonathan Liu
- Department of Physics, University of California at BerkeleyBerkeleyUnited States
| | | | - Sydney Ghoreishi
- Department of Molecular and Cell Biology, University of California at BerkeleyBerkeleyUnited States
| | - Shelby A Blythe
- Department of Molecular Biosciences, Northwestern UniversityEvanstonUnited States
| | - Hernan G Garcia
- Biophysics Graduate Group, University of California at BerkeleyBerkeleyUnited States
- Department of Physics, University of California at BerkeleyBerkeleyUnited States
- Department of Molecular and Cell Biology, University of California at BerkeleyBerkeleyUnited States
- Institute for Quantitative Biosciences-QB3, University of California at BerkeleyBerkeleyUnited States
| |
Collapse
|
16
|
Kono H, Ishida H. Nucleosome unwrapping and unstacking. Curr Opin Struct Biol 2020; 64:119-125. [DOI: 10.1016/j.sbi.2020.06.020] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2020] [Revised: 05/31/2020] [Accepted: 06/23/2020] [Indexed: 02/07/2023]
|
17
|
Neipel J, Brandani G, Schiessel H. Translational nucleosome positioning: A computational study. Phys Rev E 2020; 101:022405. [PMID: 32168683 DOI: 10.1103/physreve.101.022405] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2019] [Accepted: 11/25/2019] [Indexed: 01/26/2023]
Abstract
About three-quarters of eukaryotic DNA is wrapped into nucleosomes; DNA spools with a protein core. The affinity of a given DNA stretch to be incorporated into a nucleosome is known to depend on the base-pair sequence-dependent geometry and elasticity of the DNA double helix. This causes the rotational and translational positioning of nucleosomes. In this study we ask the question whether the latter can be predicted by a simple coarse-grained DNA model with sequence-dependent elasticity, the rigid base-pair model. Whereas this model is known to be rather robust in predicting rotational nucleosome positioning, we show that the translational positioning is a rather subtle effect that is dominated by the guanine-cytosine content dependence of entropy rather than energy. A correct qualitative prediction within the rigid base-pair framework can only be achieved by assuming that DNA elasticity effectively changes on complexation into the nucleosome complex. With that extra assumption we arrive at a model which gives an excellent quantitative agreement to experimental in vitro nucleosome maps, under the additional assumption that nucleosomes equilibrate their positions only locally.
Collapse
Affiliation(s)
- J Neipel
- Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany.,Faculty of Physics, Ludwig-Maximilians-Universität München, 80333 München, Germany.,Instituut-Lorentz, Universiteit Leiden, Postbus 9506, 2300 RA Leiden, The Netherlands
| | - G Brandani
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - H Schiessel
- Instituut-Lorentz, Universiteit Leiden, Postbus 9506, 2300 RA Leiden, The Netherlands
| |
Collapse
|
18
|
Huertas J, MacCarthy CM, Schöler HR, Cojocaru V. Nucleosomal DNA Dynamics Mediate Oct4 Pioneer Factor Binding. Biophys J 2020; 118:2280-2296. [PMID: 32027821 PMCID: PMC7202942 DOI: 10.1016/j.bpj.2019.12.038] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Revised: 12/04/2019] [Accepted: 12/23/2019] [Indexed: 01/16/2023] Open
Abstract
Transcription factor (TF) proteins bind to DNA to regulate gene expression. Normally, accessibility to DNA is required for their function. However, in the nucleus, the DNA is often inaccessible, wrapped around histone proteins in nucleosomes forming the chromatin. Pioneer TFs are thought to induce chromatin opening by recognizing their DNA binding sites on nucleosomes. For example, Oct4, a master regulator and inducer of stem cell pluripotency, binds to DNA in nucleosomes in a sequence-specific manner. Here, we reveal the structural dynamics of nucleosomes that mediate Oct4 binding from molecular dynamics simulations. Nucleosome flexibility and the amplitude of nucleosome motions such as breathing and twisting are enhanced in nucleosomes with multiple TF binding sites. Moreover, the regions around the binding sites display higher local structural flexibility. Probing different structures of Oct4-nucleosome complexes, we show that alternative configurations in which Oct4 recognizes partial binding sites display stable TF-DNA interactions similar to those observed in complexes with free DNA and compatible with the DNA curvature and DNA-histone interactions. Therefore, we propose a structural basis for nucleosome recognition by a pioneer TF that is essential for understanding how chromatin is unraveled during cell fate conversions.
Collapse
Affiliation(s)
- Jan Huertas
- In Silico Biomolecular Structure and Dynamics Group, Hubrecht Institute, Utrecht, the Netherlands; Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany; Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster, Germany; and
| | - Caitlin M MacCarthy
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany
| | - Hans R Schöler
- Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany; Medical Faculty, University of Münster, Münster, Germany
| | - Vlad Cojocaru
- In Silico Biomolecular Structure and Dynamics Group, Hubrecht Institute, Utrecht, the Netherlands; Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany; Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster, Germany; and.
| |
Collapse
|
19
|
Ensembles of Breathing Nucleosomes: A Computational Study. Biophys J 2019; 118:2297-2308. [PMID: 31882248 DOI: 10.1016/j.bpj.2019.11.3395] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2019] [Revised: 11/15/2019] [Accepted: 11/25/2019] [Indexed: 12/13/2022] Open
Abstract
About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.
Collapse
|
20
|
Onufriev AV, Schiessel H. The nucleosome: from structure to function through physics. Curr Opin Struct Biol 2019; 56:119-130. [DOI: 10.1016/j.sbi.2018.11.003] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2018] [Revised: 11/07/2018] [Accepted: 11/13/2018] [Indexed: 02/07/2023]
|
21
|
Abstract
The elasticity of the DNA double helix varies with the underlying base pair sequence. This allows one to put mechanical cues into sequences that in turn influence the packaging of DNA into nucleosomes, DNA-wrapped protein cylinders. Nucleosomes dictate a broad range of biological processes, ranging from gene regulation, recombination, and replication to chromosome condensation. Here we map base pair sequences onto graphs and use shortest paths algorithms to determine which DNA stretches are easiest or hardest to bend inside a nucleosome. We further demonstrate how genetic and mechanical information can be multiplexed by studying paths through graphs of synonymous codons. Using this method we find that nucleosomes can be placed by mechanical cues nearly everywhere on the genome of baker's yeast (Saccharomyces cerevisiae).
Collapse
Affiliation(s)
- Martijn Zuiddam
- Institute Lorentz for Theoretical Physics, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, the Netherlands
| | - Helmut Schiessel
- Institute Lorentz for Theoretical Physics, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, the Netherlands
| |
Collapse
|
22
|
Winogradoff D, Aksimentiev A. Molecular Mechanism of Spontaneous Nucleosome Unraveling. J Mol Biol 2019; 431:323-335. [PMID: 30468737 PMCID: PMC6331254 DOI: 10.1016/j.jmb.2018.11.013] [Citation(s) in RCA: 49] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2018] [Revised: 10/26/2018] [Accepted: 11/12/2018] [Indexed: 11/18/2022]
Abstract
Meters of DNA wrap around histone proteins to form nucleosomes and fit inside the micron-diameter nucleus. For the genetic information encoded in the DNA to become available for transcription, replication, and repair, the DNA-histone assembly must be disrupted. Experiment has indicated that the outer stretches of nucleosomal DNA "breathe" by spontaneously detaching from and reattaching to the histone core. Here, we report direct observation of spontaneous DNA breathing in atomistic molecular dynamics simulations, detailing a microscopic mechanism of the DNA breathing process. According to our simulations, the outer stretches of nucleosomal DNA detach in discrete steps involving 5 or 10 base pairs, with the detachment process being orchestrated by the motion of several conserved histone residues. The inner stretches of nucleosomal DNA are found to be more stably associated with the histone core by more abundant nonspecific DNA-protein contacts, providing a microscopic interpretation of nucleosome unraveling experiments. The CG content of nucleosomal DNA is found to anticorrelate with the extent of unwrapping, supporting the possibility that AT-rich segments may signal the start of transcription by forming less stable nucleosomes.
Collapse
Affiliation(s)
- David Winogradoff
- Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA; Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Aleksei Aksimentiev
- Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA; Department of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA; Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA.
| |
Collapse
|
23
|
Schiessel H. Telling Left from Right in Breathing Nucleosomes. Biophys J 2018; 115:749-750. [PMID: 30119832 DOI: 10.1016/j.bpj.2018.07.026] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2018] [Accepted: 07/26/2018] [Indexed: 11/16/2022] Open
Affiliation(s)
- Helmut Schiessel
- Institute Lorentz for Theoretical Physics, Leiden University, Leiden, the Netherlands.
| |
Collapse
|
24
|
Local DNA Sequence Controls Asymmetry of DNA Unwrapping from Nucleosome Core Particles. Biophys J 2018; 115:773-781. [PMID: 30072033 DOI: 10.1016/j.bpj.2018.07.009] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2018] [Revised: 06/18/2018] [Accepted: 07/05/2018] [Indexed: 01/06/2023] Open
Abstract
DNA is tightly wrapped around histone proteins in nucleosome core particles (NCPs) yet must become accessible for processing in the cell. This accessibility, a key component of transcription regulation, is influenced by the properties of both the histone proteins and the DNA itself. Small angle x-ray scattering with contrast variation is used to examine how sequence variations affect DNA unwrapping from NCPs at different salt concentrations. Salt destabilizes NCPs, populating multiple unwrapped states as many possible unwrapping pathways are explored by the complexes. We apply coarse-grained Monte Carlo methods to generate realistic sequence-dependent unwrapped structures for the nucleosomal DNA with thermal variations. An ensemble optimization method is employed to determine the composition of the overall ensemble as electrostatic interactions are weakened. Interesting DNA-sequence-dependent differences are revealed in the unwrapping paths and equilibrium constants. These differences are correlated with specific features within the nucleic acid sequences.
Collapse
|
25
|
Blossey R, Schiessel H. The Latest Twists in Chromatin Remodeling. Biophys J 2018; 114:2255-2261. [PMID: 29310890 DOI: 10.1016/j.bpj.2017.12.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2017] [Revised: 12/09/2017] [Accepted: 12/12/2017] [Indexed: 02/03/2023] Open
Abstract
In its most restrictive interpretation, the notion of chromatin remodeling refers to the action of chromatin-remodeling enzymes on nucleosomes with the aim of displacing and removing them from the chromatin fiber (the effective polymer formed by a DNA molecule and proteins). This local modification of the fiber structure can have consequences for the initiation and repression of the transcription process, and when the remodeling process spreads along the fiber, it also results in long-range effects essential for fiber condensation. There are three regulatory levels of relevance that can be distinguished for this process: the intrinsic sequence preference of the histone octamer, which rules the positioning of the nucleosome along the DNA, notably in relation to the genetic information coded in DNA; the recognition or selection of nucleosomal substrates by remodeling complexes; and, finally, the motor action on the nucleosome exerted by the chromatin remodeler. Recent work has been able to provide crucial insights at each of these three levels that add new twists to this exciting and unfinished story, which we highlight in this perspective.
Collapse
Affiliation(s)
- Ralf Blossey
- University of Lille 1, Unité de Glycobiologie Structurale et Fonctionnelle, CNRS UMR8576, Lille, France.
| | - Helmut Schiessel
- Institute Lorentz for Theoretical Physics, Leiden University, Leiden, the Netherlands
| |
Collapse
|