1
|
Filliaux S, Sun Z, Lyubchenko YL. Nanoscale Structure, Interactions, and Dynamics of Centromere Nucleosomes. Biomacromolecules 2024. [PMID: 38959412 DOI: 10.1021/acs.biomac.3c01440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2024]
Abstract
Centromeres are specific segments of chromosomes comprising two types of nucleosomes: canonical nucleosomes containing an octamer of H2A, H2B, H3, and H4 histones and CENP-A nucleosomes in which H3 is replaced with its analogue CENP-A. This modification leads to a difference in DNA wrapping (∼121 bp), considerably less than 147 bp in canonical nucleosomes. We used atomic force microscopy (AFM) and high-speed AFM (HS-AFM) to characterize nanoscale features and dynamics for both types of nucleosomes. For both nucleosomes, spontaneous asymmetric unwrapping of DNA was observed, and this process occurs via a transient state with ∼100 bp DNA wrapped around the core, followed by a rapid dissociation of DNA. Additionally, HS-AFM revealed higher stability of CENP-A nucleosomes compared with H3 nucleosomes in which dissociation of the histone core occurs prior to the nucleosome dissociation. These results help elucidate the differences between these nucleosomes and the potential biological necessity for CENP-A nucleosomes.
Collapse
Affiliation(s)
- Shaun Filliaux
- Department of Pharmaceutical Sciences, University of Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| | - Zhiqiang Sun
- Department of Pharmaceutical Sciences, University of Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| | - Yuri L Lyubchenko
- Department of Pharmaceutical Sciences, University of Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| |
Collapse
|
2
|
Khatua P, Tang PK, Moulick AG, Patel R, Manandhar A, Loverde SM. Sequence Dependence in Nucleosome Dynamics. J Phys Chem B 2024; 128:3090-3101. [PMID: 38530903 PMCID: PMC11181342 DOI: 10.1021/acs.jpcb.3c07363] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/28/2024]
Abstract
The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 μs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.
Collapse
Affiliation(s)
- Prabir Khatua
- Department of Chemistry, College of Staten Island, The City University of New York, 2800 Victory Boulevard, Staten Island, New York, 10314, United States
- Present Address: Department of Chemistry, GITAM School of Science, GITAM (Deemed to be University), Bengaluru, 562163, INDIA
| | - Phu K. Tang
- Department of Chemistry, College of Staten Island, The City University of New York, 2800 Victory Boulevard, Staten Island, New York, 10314, United States
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY, 10016, United States
- Present Address: Flatiron Institute, New York, NY, 10010, United States
- These authors contributed equally: Phu K. Tang, Abhik Ghosh Moulick, and Rutika Patel
| | - Abhik Ghosh Moulick
- Department of Chemistry, College of Staten Island, The City University of New York, 2800 Victory Boulevard, Staten Island, New York, 10314, United States
- These authors contributed equally: Phu K. Tang, Abhik Ghosh Moulick, and Rutika Patel
| | - Rutika Patel
- Department of Chemistry, College of Staten Island, The City University of New York, 2800 Victory Boulevard, Staten Island, New York, 10314, United States
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY, 10016, United States
- These authors contributed equally: Phu K. Tang, Abhik Ghosh Moulick, and Rutika Patel
| | - Anjela Manandhar
- Department of Chemistry, College of Staten Island, The City University of New York, 2800 Victory Boulevard, Staten Island, New York, 10314, United States
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY, 10016, United States
- Present Address: Takeda Pharmaceuticals, Boston, MA, 02139, United States
| | - Sharon M. Loverde
- Department of Chemistry, College of Staten Island, The City University of New York, 2800 Victory Boulevard, Staten Island, New York, 10314, United States
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY, 10016, United States
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, United States
- Ph.D. Program in Physics, The Graduate Center of the City University of New York, New York, NY, 10016, United States
| |
Collapse
|
3
|
Kapoor U, Kim YC, Mittal J. Coarse-Grained Models to Study Protein-DNA Interactions and Liquid-Liquid Phase Separation. J Chem Theory Comput 2024; 20:1717-1731. [PMID: 37988476 PMCID: PMC10911113 DOI: 10.1021/acs.jctc.3c00525] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2023] [Revised: 10/20/2023] [Accepted: 10/27/2023] [Indexed: 11/23/2023]
Abstract
Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define nonbonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), which is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond time scale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micrometer-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.
Collapse
Affiliation(s)
- Utkarsh Kapoor
- Artie
McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas 78743, United States
| | - Young C. Kim
- Center
for Materials Physics and Technology, Naval
Research Laboratory, Washington, District of Columbia 20375, United States
| | - Jeetain Mittal
- Artie
McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas 78743, United States
- Department
of Chemistry, Texas A&M University, College Station, Texas 78743, United States
- Interdisciplinary
Graduate Program in Genetics in Genomics, Texas A&M University, College
Station, Texas 78743, United States
| |
Collapse
|
4
|
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
|
5
|
Lin X, Zhang B. Explicit ion modeling predicts physicochemical interactions for chromatin organization. eLife 2024; 12:RP90073. [PMID: 38289342 PMCID: PMC10945522 DOI: 10.7554/elife.90073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/01/2023] Open
Abstract
Molecular mechanisms that dictate chromatin organization in vivo are under active investigation, and the extent to which intrinsic interactions contribute to this process remains debatable. A central quantity for evaluating their contribution is the strength of nucleosome-nucleosome binding, which previous experiments have estimated to range from 2 to 14 kBT. We introduce an explicit ion model to dramatically enhance the accuracy of residue-level coarse-grained modeling approaches across a wide range of ionic concentrations. This model allows for de novo predictions of chromatin organization and remains computationally efficient, enabling large-scale conformational sampling for free energy calculations. It reproduces the energetics of protein-DNA binding and unwinding of single nucleosomal DNA, and resolves the differential impact of mono- and divalent ions on chromatin conformations. Moreover, we showed that the model can reconcile various experiments on quantifying nucleosomal interactions, providing an explanation for the large discrepancy between existing estimations. We predict the interaction strength at physiological conditions to be 9 kBT, a value that is nonetheless sensitive to DNA linker length and the presence of linker histones. Our study strongly supports the contribution of physicochemical interactions to the phase behavior of chromatin aggregates and chromatin organization inside the nucleus.
Collapse
Affiliation(s)
- Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of TechnologyCambridgeUnited States
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of TechnologyCambridgeUnited States
| |
Collapse
|
6
|
Mondal A, Kolomeisky AB. Why Are Nucleosome Breathing Dynamics Asymmetric? J Phys Chem Lett 2024; 15:422-431. [PMID: 38180351 DOI: 10.1021/acs.jpclett.3c03339] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2024]
Abstract
In eukaryotic cells, DNA is bound to nucleosomes, but DNA segments occasionally unbind in the process known as nucleosome breathing. Although DNA can unwrap simultaneously from both ends of the nucleosome (symmetric breathing), experiments indicate that DNA prefers to dissociate from only one end (asymmetric breathing). However, the molecular origin of the asymmetry is not understood. We developed a new theoretical approach that gives microscopic explanations of asymmetric breathing. It is based on a stochastic description that leads to a comprehensive evaluation of dynamics by using effective free-energy landscapes. It is shown that asymmetric breathing follows the kinetically preferred pathways. In addition, it is also found that asymmetric breathing leads to a faster target search by transcription factors. Theoretical predictions, supported by computer simulations, agree with experiments. It is proposed that nature utilizes the symmetry of nucleosome breathing to achieve a better dynamic accessibility of chromatin for more efficient genetic regulation.
Collapse
Affiliation(s)
- Anupam Mondal
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | - Anatoly B Kolomeisky
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
7
|
Zhao H, Wu H, Guseman A, Abeykoon D, Camara CM, Dalal Y, Fushman D, Papoian GA. The role of cryptic ancestral symmetry in histone folding mechanisms across Eukarya and Archaea. PLoS Comput Biol 2024; 20:e1011721. [PMID: 38181064 PMCID: PMC10796010 DOI: 10.1371/journal.pcbi.1011721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Revised: 01/18/2024] [Accepted: 11/28/2023] [Indexed: 01/07/2024] Open
Abstract
Histones compact and store DNA in both Eukarya and Archaea, forming heterodimers in Eukarya and homodimers in Archaea. Despite this, the folding mechanism of histones across species remains unclear. Our study addresses this gap by investigating 11 types of histone and histone-like proteins across humans, Drosophila, and Archaea through multiscale molecular dynamics (MD) simulations, complemented by NMR and circular dichroism experiments. We confirm and elaborate on the widely applied "folding upon binding" mechanism of histone dimeric proteins and report a new alternative conformation, namely, the inverted non-native dimer, which may be a thermodynamically metastable configuration. Protein sequence analysis indicated that the inverted conformation arises from the hidden ancestral head-tail sequence symmetry underlying all histone proteins, which is congruent with the previously proposed histone evolution hypotheses. Finally, to explore the potential formations of homodimers in Eukarya, we utilized MD-based AWSEM and AI-based AlphaFold-Multimer models to predict their structures and conducted extensive all-atom MD simulations to examine their respective structural stabilities. Our results suggest that eukaryotic histones may also form stable homodimers, whereas their disordered tails bring significant structural asymmetry and tip the balance towards the formation of commonly observed heterotypic dimers.
Collapse
Affiliation(s)
- Haiqing Zhao
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, Maryland, United States of America
- Laboratory of Receptor Biology and Gene Expression, National Cancer Institute, National Institutes of Health, Bethesda, Maryland, United States of America
| | - Hao Wu
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, Maryland, United States of America
| | - Alex Guseman
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland, United States of America
| | - Dulith Abeykoon
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland, United States of America
| | - Christina M. Camara
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland, United States of America
| | - Yamini Dalal
- Laboratory of Receptor Biology and Gene Expression, National Cancer Institute, National Institutes of Health, Bethesda, Maryland, United States of America
| | - David Fushman
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, Maryland, United States of America
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland, United States of America
| | - Garegin A. Papoian
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, Maryland, United States of America
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland, United States of America
| |
Collapse
|
8
|
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
|
9
|
Yang H, Shi X. The Free Energy of Nucleosomal DNA Based on the Landau Model and Topology. Biomolecules 2023; 13:1686. [PMID: 38136559 PMCID: PMC10741420 DOI: 10.3390/biom13121686] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2023] [Revised: 10/30/2023] [Accepted: 11/09/2023] [Indexed: 12/24/2023] Open
Abstract
The free energy of nucleosomal DNA plays a key role in the formation of nucleosomes in eukaryotes. Some work on the free energy of nucleosomal DNA have been carried out in experiments. However, the relationships between the free energy of nucleosomal DNA and its conformation, especially its topology, remain unclear in theory. By combining the Landau theory, the Hopfion model and experimental data, we find that the free energy of nucleosomal DNA is at the lower level. With the help of the energy minimum principle, we conclude that nucleosomal DNA stays in a stable state. Moreover, we discover that small perturbations on nucleosomal DNA have little effect on its free energy. This implies that nucleosomal DNA has a certain redundancy in order to stay stable. This explains why nucleosomal DNA will not change significantly due to small perturbations.
Collapse
Affiliation(s)
| | - Xuguang Shi
- College of Science, Beijing Forestry University, Beijing 100083, China;
| |
Collapse
|
10
|
Mondal A, Kolomeisky AB. Role of Nucleosome Sliding in the Protein Target Search for Covered DNA Sites. J Phys Chem Lett 2023; 14:7073-7082. [PMID: 37527481 DOI: 10.1021/acs.jpclett.3c01704] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/03/2023]
Abstract
Associations of transcription factors (TFs) with specific sites on DNA initiate major cellular processes. But DNA in eukaryotic cells is covered by nucleosomes which prevent TFs from binding. However, nucleosome structures on DNA are not static and exhibit breathing and sliding. We develop a theoretical framework to investigate the effect of nucleosome sliding on a protein target search. By analysis of a discrete-state stochastic model of nucleosome sliding, search dynamics are explicitly evaluated. It is found that for long sliding lengths the target search dynamics are faster for normal TFs that cannot enter the nucleosomal DNA. But for more realistic short sliding lengths, the so-called pioneer TFs, which can invade nucleosomal DNA, locate specific sites faster. It is also suggested that nucleosome breathing, which is a faster process, has a stronger effect on protein search dynamics than that of nucleosome sliding. Theoretical arguments to explain these observations are presented.
Collapse
Affiliation(s)
- Anupam Mondal
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | - Anatoly B Kolomeisky
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
11
|
Chakraborty D, Mondal B, Thirumalai D. Brewing COFFEE: A sequence-specific coarse-grained energy function for simulations of DNA-protein complexes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.07.544064. [PMID: 37333386 PMCID: PMC10274755 DOI: 10.1101/2023.06.07.544064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/20/2023]
Abstract
DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription, 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 re-calibrating 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 dataset of high-resolution crystal structures. The only parameter in COFFEE is the strength (λ D N A P R O ) of the DNA-protein contact potential. For an optimal choice of λ D N A P R O , 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 the scattering profiles that are in quantitative agreement with SAXS 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 does not alter the balance of electrostatic interactions, but affects 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 St, Stop A5300, Austin TX 78712, USA
| | - Balaka Mondal
- Department of Chemistry, The University of Texas at Austin, 105 E 24th St, Stop A5300, Austin TX 78712, USA
| | - D Thirumalai
- Department of Chemistry, The University of Texas at Austin, 105 E 24th St, Stop A5300, Austin TX 78712, USA
- Department of Physics, The University of Texas at Austin, 2515 Speedway,Austin TX 78712, USA
| |
Collapse
|
12
|
Kapoor U, Kim YC, Mittal J. A coarse-grained DNA model to study protein-DNA interactions and liquid-liquid phase separation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.05.19.541513. [PMID: 37292850 PMCID: PMC10245785 DOI: 10.1101/2023.05.19.541513] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define non-bonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), that is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond timescale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micron-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.
Collapse
Affiliation(s)
- Utkarsh Kapoor
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas 78743, United States
| | - Young C. Kim
- Center for Materials Physics and Technology, Naval Research Laboratory, Washington, District of Columbia
| | - Jeetain Mittal
- Artie McFerrin Department of Chemical Engineering, Texas A&M University, College Station, Texas 78743, United States
- Department of Chemistry, Texas A&M University, College Station, Texas 78743, United States
- Interdisciplinary Graduate Program in Genetics in Genomics, Texas A&M University, College Station, Texas 78743, United States
| |
Collapse
|
13
|
Mondal A, Felipe C, Kolomeisky AB. Nucleosome Breathing Facilitates the Search for Hidden DNA Sites by Pioneer Transcription Factors. J Phys Chem Lett 2023; 14:4096-4103. [PMID: 37125729 DOI: 10.1021/acs.jpclett.3c00529] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Transfer of genetic information starts with transcription factors (TFs) binding to specific sites on DNA. But in living cells, DNA is mostly covered by nucleosomes. There are proteins, known as pioneer TFs, that can efficiently reach the DNA sites hidden by nucleosomes, although the underlying mechanisms are not understood. Using the recently proposed idea of interaction-compensation mechanism, we develop a stochastic model for the target search on DNA with nucleosome breathing. It is found that nucleosome breathing can significantly accelerate the search by pioneer TFs in comparison to situations without breathing. We argue that this is the result of the interaction-compensation mechanism that allows proteins to enter the inner nucleosome region through the outer DNA segment. It is suggested that nature optimized pioneer TFs to take advantage of nucleosome breathing. The presented theoretical picture provides a possible microscopic explanation for the successful invasion of nucleosome-buried genes.
Collapse
Affiliation(s)
- Anupam Mondal
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | - Cayke Felipe
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, United States
| | - Anatoly B Kolomeisky
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, United States
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
14
|
Li S, Wei T, Panchenko AR. Histone variant H2A.Z modulates nucleosome dynamics to promote DNA accessibility. Nat Commun 2023; 14:769. [PMID: 36765119 PMCID: PMC9918499 DOI: 10.1038/s41467-023-36465-5] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2022] [Accepted: 02/02/2023] [Indexed: 02/12/2023] Open
Abstract
Nucleosomes, containing histone variants H2A.Z, are important for gene transcription initiation and termination, chromosome segregation and DNA double-strand break repair, among other functions. However, the underlying mechanisms of how H2A.Z influences nucleosome stability, dynamics and DNA accessibility are not well understood, as experimental and computational evidence remains inconclusive. Our modeling efforts of human nucleosome stability and dynamics, along with comparisons with experimental data show that the incorporation of H2A.Z results in a substantial decrease of the energy barrier for DNA unwrapping. This leads to the spontaneous DNA unwrapping of about forty base pairs from both ends, nucleosome gapping and increased histone plasticity, which otherwise is not observed for canonical nucleosomes. We demonstrate that both N- and C-terminal tails of H2A.Z play major roles in these events, whereas the H3.3 variant exerts a negligible impact in modulating the DNA end unwrapping. In summary, our results indicate that H2A.Z deposition makes nucleosomes more mobile and DNA more accessible to transcriptional machinery and other chromatin components.
Collapse
Affiliation(s)
- Shuxiang Li
- Department of Pathology and Molecular Medicine, Queen's University, Kingston, ON, Canada
| | - Tiejun Wei
- Department of Pathology and Molecular Medicine, Queen's University, Kingston, ON, Canada
| | - Anna R Panchenko
- Department of Pathology and Molecular Medicine, Queen's University, Kingston, ON, Canada. .,Department of Biology and Molecular Sciences, Queen's University, Kingston, ON, Canada. .,School of Computing, Queen's University, Kingston, ON, Canada. .,Ontario Institute of Cancer Research, Toronto, Canada.
| |
Collapse
|
15
|
Meng Y, Liu R, Wang L, Li F, Tian Y, Lu H. Binding affinity and conformational change predictions for a series of inhibitors with RuBisCO in a carbon dioxide gas and water environment by multiple computational methods. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121478] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/19/2023]
|
16
|
Zhang H, Huo QY, Gao YQ. DNA Sequence-Dependent Binding of Linker Histone gH1 Regulates Nucleosome Conformations. J Phys Chem B 2022; 126:6771-6779. [PMID: 36062461 DOI: 10.1021/acs.jpcb.2c03785] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Sequence-dependent binding between DNA and proteins in chromatin is an essential part of gene expression. Linker histone H1 is an important protein in the regulation of chromatin compartmentalization and compaction, and its binding with the nucleosome is sensitive to the DNA sequence. Although the interactions of H1 and DNA have been widely investigated, the mechanism of nucleosome conformation changes induced by the DNA-sequence-dependent binding with gH1 (globular H1.0) remains largely unclear at the atomic level. In the present molecular dynamics simulations, both linker and dyad DNAs were mutated to investigate the conformational changes of the nucleosome induced by the sequence-dependent binding of gH1 based on the on-dyad binding mode. Our results indicate that gH1 is insensitive to the DNA sequence of the dyad DNA but presents an apparent preference to linker DNA with an AT-rich sequence. Moreover, this specific binding induces the entry/exit region of a nucleosome to a tight conformation and regulates the accessibility of core histones. Considering that the entry/exit region of the nucleosome is a crucial binding site for many functional proteins related to gene expression, the conformational change at this region could represent an important gene regulation signal.
Collapse
Affiliation(s)
- Hong Zhang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Qin Yuan Huo
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China
| |
Collapse
|
17
|
Liu S, Lin X, Zhang B. Chromatin fiber breaks into clutches under tension and crowding. Nucleic Acids Res 2022; 50:9738-9747. [PMID: 36029149 PMCID: PMC9508854 DOI: 10.1093/nar/gkac725] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2022] [Revised: 07/08/2022] [Accepted: 08/12/2022] [Indexed: 11/18/2022] Open
Abstract
The arrangement of nucleosomes inside chromatin is of extensive interest. While in vitro experiments have revealed the formation of 30 nm fibers, most in vivo studies have failed to confirm their presence in cell nuclei. To reconcile the diverging experimental findings, we characterized chromatin organization using a residue-level coarse-grained model. The computed force–extension curve matches well with measurements from single-molecule experiments. Notably, we found that a dodeca-nucleosome in the two-helix zigzag conformation breaks into structures with nucleosome clutches and a mix of trimers and tetramers under tension. Such unfolded configurations can also be stabilized through trans interactions with other chromatin chains. Our study suggests that unfolding from chromatin fibers could contribute to the irregularity of in vivo chromatin configurations. We further revealed that chromatin segments with fibril or clutch structures engaged in distinct binding modes and discussed the implications of these inter-chain interactions for a potential sol–gel phase transition.
Collapse
Affiliation(s)
- Shuming Liu
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| |
Collapse
|
18
|
Ishida H, Kono H. Free Energy Landscape of H2A-H2B Displacement From Nucleosome. J Mol Biol 2022; 434:167707. [PMID: 35777463 DOI: 10.1016/j.jmb.2022.167707] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 06/11/2022] [Accepted: 06/23/2022] [Indexed: 12/14/2022]
Abstract
Nucleosome reconstitution plays an important role in many cellular functions. As an initial step, H2A-H2B dimer displacement, which is accompanied by disruption of many of the interactions within the nucleosome, should occur. To understand how H2A-H2B dimer displacement occurs, an adaptively biased molecular dynamics (ABMD) simulation was carried out to generate a variety of displacements of the H2A-H2B dimer from the fully wrapped to partially unwrapped nucleosome structures. With regards to these structures, the free energy landscape of the dimer displacement was investigated using umbrella sampling simulations. We found that the main contributors to the free energy were the docking domain of H2A and the C-terminal of H4. There were various paths for the dimer displacement which were dependent on the extent of nucleosomal DNA wrapping, suggesting that modulation of the intra-nucleosomal interaction by external factors such as histone chaperons could control the path for the H2A-H2B dimer displacement. Key residues which contributed to the free energy have also been reported to be involved in the mutations and posttranslational modifications (PTMs) which are important for assembling and/or reassembling the nucleosome at the molecular level and are found in cancer cells at the phenotypic level. Our results give insight into how the H2A-H2B dimer displacement proceeds along various paths according to different interactions within the nucleosome.
Collapse
Affiliation(s)
- Hisashi Ishida
- Institute for Quantum Life Science, National Institutes for Quantum Science and Technology, 619-0215 Kizugawa, Kyoto, Japan.
| | - Hidetoshi Kono
- Institute for Quantum Life Science, National Institutes for Quantum Science and Technology, 619-0215 Kizugawa, Kyoto, Japan
| |
Collapse
|
19
|
Huertas J, Woods EJ, Collepardo-Guevara R. Multiscale modelling of chromatin organisation: Resolving nucleosomes at near-atomistic resolution inside genes. Curr Opin Cell Biol 2022; 75:102067. [DOI: 10.1016/j.ceb.2022.02.001] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Revised: 01/24/2022] [Accepted: 02/04/2022] [Indexed: 12/15/2022]
|
20
|
Release of linker histone from the nucleosome driven by polyelectrolyte competition with a disordered protein. Nat Chem 2022; 14:224-231. [PMID: 34992286 DOI: 10.1038/s41557-021-00839-3] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2020] [Accepted: 10/19/2021] [Indexed: 12/13/2022]
Abstract
Highly charged intrinsically disordered proteins are essential regulators of chromatin structure and transcriptional activity. Here we identify a surprising mechanism of molecular competition that relies on the pronounced dynamical disorder present in these polyelectrolytes and their complexes. The highly positively charged human linker histone H1.0 (H1) binds to nucleosomes with ultrahigh affinity, implying residence times incompatible with efficient biological regulation. However, we show that the disordered regions of H1 retain their large-amplitude dynamics when bound to the nucleosome, which enables the highly negatively charged and disordered histone chaperone prothymosin α to efficiently invade the H1-nucleosome complex and displace H1 via a competitive substitution mechanism, vastly accelerating H1 dissociation. By integrating experiments and simulations, we establish a molecular model that rationalizes the remarkable kinetics of this process structurally and dynamically. Given the abundance of polyelectrolyte sequences in the nuclear proteome, this mechanism is likely to be widespread in cellular regulation.
Collapse
|
21
|
Shino G, Takada S. Modeling DNA Opening in the Eukaryotic Transcription Initiation Complexes via Coarse-Grained Models. Front Mol Biosci 2021; 8:772486. [PMID: 34869598 PMCID: PMC8636136 DOI: 10.3389/fmolb.2021.772486] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 10/07/2021] [Indexed: 01/24/2023] Open
Abstract
Recently, the molecular mechanisms of transcription initiation have been intensively studied. Especially, the cryo-electron microscopy revealed atomic structure details in key states in the eukaryotic transcription initiation. Yet, the dynamic processes of the promoter DNA opening in the pre-initiation complex remain obscured. In this study, based on the three cryo-electron microscopic yeast structures for the closed, open, and initially transcribing complexes, we performed multiscale molecular dynamics (MD) simulations to model structures and dynamic processes of DNA opening. Combining coarse-grained and all-atom MD simulations, we first obtained the atomic model for the DNA bubble in the open complexes. Then, in the MD simulation from the open to the initially transcribing complexes, we found a previously unidentified intermediate state which is formed by the bottleneck in the fork loop 1 of Pol II: The loop opening triggered the escape from the intermediate, serving as a gatekeeper of the promoter DNA opening. In the initially transcribing complex, the non-template DNA strand passes a groove made of the protrusion, the lobe, and the fork of Rpb2 subunit of Pol II, in which several positively charged and highly conserved residues exhibit key interactions to the non-template DNA strand. The back-mapped all-atom models provided further insights on atomistic interactions such as hydrogen bonding and can be used for future simulations.
Collapse
Affiliation(s)
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| |
Collapse
|
22
|
Galaz-Davison P, Román EA, Ramírez-Sarmiento CA. The N-terminal domain of RfaH plays an active role in protein fold-switching. PLoS Comput Biol 2021; 17:e1008882. [PMID: 34478435 PMCID: PMC8454952 DOI: 10.1371/journal.pcbi.1008882] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Revised: 09/21/2021] [Accepted: 08/07/2021] [Indexed: 11/19/2022] Open
Abstract
The bacterial elongation factor RfaH promotes the expression of virulence factors by specifically binding to RNA polymerases (RNAP) paused at a DNA signal. This behavior is unlike that of its paralog NusG, the major representative of the protein family to which RfaH belongs. Both proteins have an N-terminal domain (NTD) bearing an RNAP binding site, yet NusG C-terminal domain (CTD) is folded as a β-barrel while RfaH CTD is forming an α-hairpin blocking such site. Upon recognition of the specific DNA exposed by RNAP, RfaH is activated via interdomain dissociation and complete CTD structural rearrangement into a β-barrel structurally identical to NusG CTD. Although RfaH transformation has been extensively characterized computationally, little attention has been given to the role of the NTD in the fold-switching process, as its structure remains unchanged. Here, we used Associative Water-mediated Structure and Energy Model (AWSEM) molecular dynamics to characterize the transformation of RfaH, spotlighting the sequence-dependent effects of NTD on CTD fold stabilization. Umbrella sampling simulations guided by native contacts recapitulate the thermodynamic equilibrium experimentally observed for RfaH and its isolated CTD. Temperature refolding simulations of full-length RfaH show a high success towards α-folded CTD, whereas the NTD interferes with βCTD folding, becoming trapped in a β-barrel intermediate. Meanwhile, NusG CTD refolding is unaffected by the presence of RfaH NTD, showing that these NTD-CTD interactions are encoded in RfaH sequence. Altogether, these results suggest that the NTD of RfaH favors the α-folded RfaH by specifically orienting the αCTD upon interdomain binding and by favoring β-barrel rupture into an intermediate from which fold-switching proceeds. Proteins commonly adopt a single three-dimensional structure that is required for biological function. Nevertheless, proteins are not isolated in the cell, and the presence of binding partners can give rise to alternate structural configurations. Metamorphic proteins represent an extreme case of the latter, by folding into at least two well-defined configurations that are both structurally and functionally different. For RfaH, a virulence factor in enterobacteria, two distinct folds are found: an autoinhibited state in which its two protein domains strongly interact, and an active state in which these domains dissociate due to a specific DNA signal on RNA polymerases. This activation is accompanied by the refolding of the C-terminal domain (CTD) from an α-helical structure to a β-barrel. Our work employs computational simulations to explore the role of the N-terminal domain (NTD) in regulating the metamorphic behavior of RfaH, determining that this domain has a major part in orienting and binding to the CTD in its α-helical fold, and in stabilizing an intermediate state instead of the fully folded β-barrel. These results suggest that the NTD not only participates in stabilizing the autoinhibited state, but also aids in fold-switching back to it after active RfaH is released from RNA polymerase.
Collapse
Affiliation(s)
- Pablo Galaz-Davison
- Institute for Biological and Medical Engineering, Schools of Engineering, Medicine and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile
- ANID–Millennium Science Initiative Program–Millennium Institute for Integrative Biology (iBio), Santiago, Chile
| | - Ernesto A. Román
- Instituto de Química y Fisicoquímica Biológicas (UBA-CONICET), Ciudad Autónoma de Buenos Aires, Argentina
- Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Autónoma de Buenos Aires, Argentina
| | - César A. Ramírez-Sarmiento
- Institute for Biological and Medical Engineering, Schools of Engineering, Medicine and Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile
- ANID–Millennium Science Initiative Program–Millennium Institute for Integrative Biology (iBio), Santiago, Chile
- * E-mail:
| |
Collapse
|
23
|
Ma YW, Lin TY, Tsai MY. Fibril Surface-Dependent Amyloid Precursors Revealed by Coarse-Grained Molecular Dynamics Simulation. Front Mol Biosci 2021; 8:719320. [PMID: 34422910 PMCID: PMC8378332 DOI: 10.3389/fmolb.2021.719320] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Accepted: 07/26/2021] [Indexed: 01/05/2023] Open
Abstract
Amyloid peptides are known to self-assemble into larger aggregates that are linked to the pathogenesis of many neurodegenerative disorders. In contrast to primary nucleation, recent experimental and theoretical studies have shown that many toxic oligomeric species are generated through secondary processes on a pre-existing fibrillar surface. Nucleation, for example, can also occur along the surface of a pre-existing fibril—secondary nucleation—as opposed to the primary one. However, explicit pathways are still not clear. In this study, we use molecular dynamics simulation to explore the free energy landscape of a free Abeta monomer binding to an existing fibrillar surface. We specifically look into several potential Abeta structural precursors that might precede some secondary events, including elongation and secondary nucleation. We find that the overall process of surface-dependent events can be described at least by the following three stages: 1. Free diffusion 2. Downhill guiding 3. Dock and lock. And we show that the outcome of adding a new monomer onto a pre-existing fibril is pathway-dependent, which leads to different secondary processes. To understand structural details, we have identified several monomeric amyloid precursors over the fibrillar surfaces and characterize their heterogeneity using a probability contact map analysis. Using the frustration analysis (a bioinformatics tool), we show that surface heterogeneity correlates with the energy frustration of specific local residues that form binding sites on the fibrillar structure. We further investigate the helical twisting of protofilaments of different sizes and observe a length dependence on the filament twisting. This work presents a comprehensive survey over the properties of fibril growth using a combination of several openMM-based platforms, including the GPU-enabled openAWSEM package for coarse-grained modeling, MDTraj for trajectory analysis, and pyEMMA for free energy calculation. This combined approach makes long-timescale simulation for aggregation systems as well as all-in-one analysis feasible. We show that this protocol allows us to explore fibril stability, surface binding affinity/heterogeneity, as well as fibrillar twisting. All these properties are important for understanding the molecular mechanism of surface-catalyzed secondary processes of fibril growth.
Collapse
Affiliation(s)
- Yuan-Wei Ma
- Department of Chemistry, Tamkang University, New Taipei City, Taiwan
| | - Tong-You Lin
- Department of Chemistry, Tamkang University, New Taipei City, Taiwan
| | - Min-Yeh Tsai
- Department of Chemistry, Tamkang University, New Taipei City, Taiwan
| |
Collapse
|
24
|
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
|
25
|
Torsional stress can regulate the unwrapping of two outer half superhelical turns of nucleosomal DNA. Proc Natl Acad Sci U S A 2021; 118:2020452118. [PMID: 33558240 DOI: 10.1073/pnas.2020452118] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Torsional stress has a significant impact on the structure and stability of the nucleosome. RNA polymerase imposes torsional stress on the DNA in chromatin and unwraps the DNA from the nucleosome to access the genetic information encoded in the DNA. To understand how the torsional stress affects the stability of the nucleosome, we examined the unwrapping of two half superhelical turns of nucleosomal DNA from either end of the DNA under torsional stress with all-atom molecular dynamics simulations. The free energies for unwrapping the DNA indicate that positive stress that overtwists DNA facilitates a large-scale asymmetric unwrapping of the DNA without a large extension of the DNA. During the unwrapping, one end of the DNA was dissociated from H3 and H2A-H2B, while the other end of the DNA stably remained wrapped. The detailed analysis indicates that this asymmetric dissociation is facilitated by the geometry and bendability of the DNA under positive stress. The geometry stabilized the interaction between the major groove of the twisted DNA and the H3 αN-helix, and the straightened DNA destabilized the interaction with H2A-H2B. Under negative stress, the DNA became more bendable and flexible, which facilitated the binding of the unwrapped DNA to the octamer in a stable state. Consequently, we conclude that the torsional stress has a significant impact on the affinity of the DNA and the octamer through the inherent nature of the DNA and can change the accessibility of regulatory proteins.
Collapse
|
26
|
Lin X, Leicher R, Liu S, Zhang B. Cooperative DNA looping by PRC2 complexes. Nucleic Acids Res 2021; 49:6238-6248. [PMID: 34057467 PMCID: PMC8216278 DOI: 10.1093/nar/gkab441] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2020] [Revised: 04/25/2021] [Accepted: 05/07/2021] [Indexed: 01/14/2023] Open
Abstract
Polycomb repressive complex 2 (PRC2) is an essential protein complex that silences gene expression via post-translational modifications of chromatin. This paper combined homology modeling, atomistic and coarse-grained molecular dynamics simulations, and single-molecule force spectroscopy experiments to characterize both its full-length structure and PRC2-DNA interactions. Using free energy calculations with a newly parameterized protein-DNA force field, we studied a total of three potential PRC2 conformations and their impact on DNA binding and bending. Consistent with cryo-EM studies, we found that EZH2, a core subunit of PRC2, provides the primary interface for DNA binding, and its curved surface can induce DNA bending. Our simulations also predicted the C2 domain of the SUZ12 subunit to contact DNA. Multiple PRC2 complexes bind with DNA cooperatively via allosteric communication through the DNA, leading to a hairpin-like looped configuration. Single-molecule experiments support PRC2-mediated DNA looping and the role of AEBP2 in regulating such loop formation. The impact of AEBP2 can be partly understood from its association with the C2 domain, blocking C2 from DNA binding. Our study suggests that accessory proteins may regulate the genomic location of PRC2 by interfering with its DNA interactions.
Collapse
Affiliation(s)
- Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Rachel Leicher
- Laboratory of Nanoscale Biophysics and Biochemistry, The Rockefeller University, New York, NY 10065, USA.,Tri-Institutional PhD Program in Chemical Biology, New York, NY 10065, USA
| | - Shixin Liu
- Laboratory of Nanoscale Biophysics and Biochemistry, The Rockefeller University, New York, NY 10065, USA
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| |
Collapse
|
27
|
Lin X, Qi Y, Latham AP, Zhang B. Multiscale modeling of genome organization with maximum entropy optimization. J Chem Phys 2021; 155:010901. [PMID: 34241389 PMCID: PMC8253599 DOI: 10.1063/5.0044150] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Accepted: 04/28/2021] [Indexed: 12/15/2022] Open
Abstract
Three-dimensional (3D) organization of the human genome plays an essential role in all DNA-templated processes, including gene transcription, gene regulation, and DNA replication. Computational modeling can be an effective way of building high-resolution genome structures and improving our understanding of these molecular processes. However, it faces significant challenges as the human genome consists of over 6 × 109 base pairs, a system size that exceeds the capacity of traditional modeling approaches. In this perspective, we review the progress that has been made in modeling the human genome. Coarse-grained models parameterized to reproduce experimental data via the maximum entropy optimization algorithm serve as effective means to study genome organization at various length scales. They have provided insight into the principles of whole-genome organization and enabled de novo predictions of chromosome structures from epigenetic modifications. Applications of these models at a near-atomistic resolution further revealed physicochemical interactions that drive the phase separation of disordered proteins and dictate chromatin stability in situ. We conclude with an outlook on the opportunities and challenges in studying chromosome dynamics.
Collapse
Affiliation(s)
- Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Yifeng Qi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Andrew P Latham
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
28
|
Reddy G, Thirumalai D. Asymmetry in histone rotation in forced unwrapping and force quench rewrapping in a nucleosome. Nucleic Acids Res 2021; 49:4907-4918. [PMID: 33877361 PMCID: PMC8136794 DOI: 10.1093/nar/gkab263] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Revised: 02/25/2021] [Accepted: 03/30/2021] [Indexed: 01/07/2023] Open
Abstract
Single molecule pulling experiments have shown that DNA in the nucleosomes unwraps in two stages from the histone protein core (HPC). The first stage, attributed to the rupture of the outer DNA turn, occurs between 3 and 5 pNs, and is reversible. The inner DNA turn ruptures irreversibly at forces between 9 and 15 pNs (or higher) in the second stage. Molecular simulations using the Self-Organized Polymer model capture the experimental findings. The unwrapping of the outer DNA turn is independent of the pulling direction. The rupture of the DNA inner turn depends on the pulling direction and involves overcoming substantial energetic (most likely electrostatic in origin) and kinetic barriers. They arise because the mechanical force has to generate sufficient torque to rotate the HPC by 180°. On the other hand, during the rewrapping process, HPC rotation is stochastic, with force playing no role. The assembly of the outer DNA wrap upon force quench nearly coincides with the unwrapping process, confirming the reversibility of the outer turn rupture. The asymmetry in HPC rotation during unwrapping and rewrapping explains the observed hysteresis in the stretch-release cycles in experiments. We propose experiments to test the prediction that HPC rotation produces kinetic barriers in the unwrapping process.
Collapse
Affiliation(s)
- Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012, Karnataka, India
| | - D Thirumalai
- Department of Chemistry, The University of Texas, Austin, TX 78712, USA
| |
Collapse
|
29
|
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
|
30
|
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
|
31
|
Farr SE, Woods EJ, Joseph JA, Garaizar A, Collepardo-Guevara R. Nucleosome plasticity is a critical element of chromatin liquid-liquid phase separation and multivalent nucleosome interactions. Nat Commun 2021; 12:2883. [PMID: 34001913 PMCID: PMC8129070 DOI: 10.1038/s41467-021-23090-3] [Citation(s) in RCA: 71] [Impact Index Per Article: 23.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 03/23/2021] [Indexed: 12/19/2022] Open
Abstract
Liquid-liquid phase separation (LLPS) is an important mechanism that helps explain the membraneless compartmentalization of the nucleus. Because chromatin compaction and LLPS are collective phenomena, linking their modulation to the physicochemical features of nucleosomes is challenging. Here, we develop an advanced multiscale chromatin model-integrating atomistic representations, a chemically-specific coarse-grained model, and a minimal model-to resolve individual nucleosomes within sub-Mb chromatin domains and phase-separated systems. To overcome the difficulty of sampling chromatin at high resolution, we devise a transferable enhanced-sampling Debye-length replica-exchange molecular dynamics approach. We find that nucleosome thermal fluctuations become significant at physiological salt concentrations and destabilize the 30-nm fiber. Our simulations show that nucleosome breathing favors stochastic folding of chromatin and promotes LLPS by simultaneously boosting the transient nature and heterogeneity of nucleosome-nucleosome contacts, and the effective nucleosome valency. Our work puts forward the intrinsic plasticity of nucleosomes as a key element in the liquid-like behavior of nucleosomes within chromatin, and the regulation of chromatin LLPS.
Collapse
Affiliation(s)
- Stephen E Farr
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge, UK
| | - Esmae J Woods
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge, UK
| | - Jerelle A Joseph
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge, UK
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK
- Department of Genetics, University of Cambridge, Cambridge, UK
| | - Adiran Garaizar
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge, UK
| | - Rosana Collepardo-Guevara
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, Cambridge, UK.
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK.
- Department of Genetics, University of Cambridge, Cambridge, UK.
| |
Collapse
|
32
|
Histone dynamics mediate DNA unwrapping and sliding in nucleosomes. Nat Commun 2021; 12:2387. [PMID: 33888707 PMCID: PMC8062685 DOI: 10.1038/s41467-021-22636-9] [Citation(s) in RCA: 52] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Accepted: 03/23/2021] [Indexed: 02/06/2023] Open
Abstract
Nucleosomes are elementary building blocks of chromatin in eukaryotes. They tightly wrap ∼147 DNA base pairs around an octamer of histone proteins. How nucleosome structural dynamics affect genome functioning is not completely clear. Here we report all-atom molecular dynamics simulations of nucleosome core particles at a timescale of 15 microseconds. At this timescale, functional modes of nucleosome dynamics such as spontaneous nucleosomal DNA breathing, unwrapping, twisting, and sliding were observed. We identified atomistic mechanisms of these processes by analyzing the accompanying structural rearrangements of the histone octamer and histone-DNA contacts. Octamer dynamics and plasticity were found to enable DNA unwrapping and sliding. Through multi-scale modeling, we showed that nucleosomal DNA dynamics contribute to significant conformational variability of the chromatin fiber at the supranucleosomal level. Our study further supports mechanistic coupling between fine details of histone dynamics and chromatin functioning, provides a framework for understanding the effects of various chromatin modifications.
Collapse
|
33
|
Procyk J, Poppleton E, Šulc P. Coarse-grained nucleic acid-protein model for hybrid nanotechnology. SOFT MATTER 2021; 17:3586-3593. [PMID: 33398312 DOI: 10.1039/d0sm01639j] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The emerging field of hybrid DNA-protein nanotechnology brings with it the potential for many novel materials which combine the addressability of DNA nanotechnology with the versatility of protein interactions. However, the design and computational study of these hybrid structures is difficult due to the system sizes involved. To aid in the design and in silico analysis process, we introduce here a coarse-grained DNA/RNA-protein model that extends the oxDNA/oxRNA models of DNA/RNA with a coarse-grained model of proteins based on an anisotropic network model representation. Fully equipped with analysis scripts and visualization, our model aims to facilitate hybrid nanomaterial design towards eventual experimental realization, as well as enabling study of biological complexes. We further demonstrate its usage by simulating DNA-protein nanocage, DNA wrapped around histones, and a nascent RNA in polymerase.
Collapse
Affiliation(s)
- Jonah Procyk
- School of Molecular Sciences and Center for Molecular Design and Biomimetics, The Biodesign Institute, Arizona State University, 1001 South McAllister Avenue, Tempe, Arizona 85281, USA.
| | | | | |
Collapse
|
34
|
Mukherjee A, Saurabh S, Olive E, Jang YH, Lansac Y. Protamine Binding Site on DNA: Molecular Dynamics Simulations and Free Energy Calculations with Full Atomistic Details. J Phys Chem B 2021; 125:3032-3044. [DOI: 10.1021/acs.jpcb.0c09166] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Arnab Mukherjee
- GREMAN, CNRS UMR 7347, Université de Tours, 37200 Tours, France
| | - Suman Saurabh
- GREMAN, CNRS UMR 7347, Université de Tours, 37200 Tours, France
| | - Enrick Olive
- GREMAN, CNRS UMR 7347, Université de Tours, 37200 Tours, France
| | - Yun Hee Jang
- Department of Energy Science and Engineering, DGIST, Daegu 42988, Korea
| | - Yves Lansac
- GREMAN, CNRS UMR 7347, Université de Tours, 37200 Tours, France
- Department of Energy Science and Engineering, DGIST, Daegu 42988, Korea
- Laboratoire de Physique des Solides, CNRS UMR 8502, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
35
|
Winogradoff D, Li P, Joshi H, Quednau L, Maffeo C, Aksimentiev A. Chiral Systems Made from DNA. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2021; 8:2003113. [PMID: 33717850 PMCID: PMC7927625 DOI: 10.1002/advs.202003113] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2020] [Revised: 10/13/2020] [Indexed: 05/05/2023]
Abstract
The very chemical structure of DNA that enables biological heredity and evolution has non-trivial implications for the self-organization of DNA molecules into larger assemblies and provides limitless opportunities for building functional nanostructures. This progress report discusses the natural organization of DNA into chiral structures and recent advances in creating synthetic chiral systems using DNA as a building material. How nucleic acid chirality naturally comes into play in a diverse array of situations is considered first, at length scales ranging from an individual nucleotide to entire chromosomes. Thereafter, chiral liquid crystal phases formed by dense DNA mixtures are discussed, including the ongoing efforts to understand their origins. The report then summarizes recent efforts directed toward building chiral structures, and other structures of complex topology, using the principle of DNA self-assembly. Discussed last are existing and proposed functional man-made nanostructures designed to either probe or harness DNA's chirality, from plasmonics and spintronics to biosensing.
Collapse
Affiliation(s)
- David Winogradoff
- Center for the Physics of Living CellsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
- Department of PhysicsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
| | - Pin‐Yi Li
- Department of PhysicsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
| | - Himanshu Joshi
- Department of PhysicsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
| | - Lauren Quednau
- Center for the Physics of Living CellsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
| | - Christopher Maffeo
- Center for the Physics of Living CellsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
- Department of PhysicsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
- Beckman Institute for Advanced Science and TechnologyUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
| | - Aleksei Aksimentiev
- Center for the Physics of Living CellsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
- Department of PhysicsUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
- Beckman Institute for Advanced Science and TechnologyUniversity of Illinois at Urbana–ChampaignUrbanaILUSA
| |
Collapse
|
36
|
Wu H, Dalal Y, Papoian GA. Binding Dynamics of Disordered Linker Histone H1 with a Nucleosomal Particle. J Mol Biol 2021; 433:166881. [PMID: 33617899 DOI: 10.1016/j.jmb.2021.166881] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 02/11/2021] [Accepted: 02/12/2021] [Indexed: 01/30/2023]
Abstract
Linker histone H1 is an essential regulatory protein for many critical biological processes, such as eukaryotic chromatin packaging and gene expression. Mis-regulation of H1s is commonly observed in tumor cells, where the balance between different H1 subtypes has been shown to alter the cancer phenotype. Consisting of a rigid globular domain and two highly charged terminal domains, H1 can bind to multiple sites on a nucleosomal particle to alter chromatin hierarchical condensation levels. In particular, the disordered H1 amino- and carboxyl-terminal domains (NTD/CTD) are believed to enhance this binding affinity, but their detailed dynamics and functions remain unclear. In this work, we used a coarse-grained computational model, AWSEM-DNA, to simulate the H1.0b-nucleosome complex, namely chromatosome. Our results demonstrate that H1 disordered domains restrict the dynamics and conformation of both globular H1 and linker DNA arms, resulting in a more compact and rigid chromatosome particle. Furthermore, we identified regions of H1 disordered domains that are tightly tethered to DNA near the entry-exit site. Overall, our study elucidates at near-atomic resolution the way the disordered linker histone H1 modulates nucleosome's structural preferences and conformational dynamics.
Collapse
Affiliation(s)
- Hao Wu
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, United States
| | - Yamini Dalal
- Laboratory of Receptor Biology and Gene Expression, National Cancer Institute, National Institutes of Health, Bethesda, MD 20892, United States.
| | - Garegin A Papoian
- Biophysics Program, Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742, United States; Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, United States.
| |
Collapse
|
37
|
Ding X, Lin X, Zhang B. Stability and folding pathways of tetra-nucleosome from six-dimensional free energy surface. Nat Commun 2021; 12:1091. [PMID: 33597548 PMCID: PMC7889939 DOI: 10.1038/s41467-021-21377-z] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2020] [Accepted: 01/22/2021] [Indexed: 01/01/2023] Open
Abstract
The three-dimensional organization of chromatin is expected to play critical roles in regulating genome functions. High-resolution characterization of its structure and dynamics could improve our understanding of gene regulation mechanisms but has remained challenging. Using a near-atomistic model that preserves the chemical specificity of protein-DNA interactions at residue and base-pair resolution, we studied the stability and folding pathways of a tetra-nucleosome. Dynamical simulations performed with an advanced sampling technique uncovered multiple pathways that connect open chromatin configurations with the zigzag crystal structure. Intermediate states along the simulated folding pathways resemble chromatin configurations reported from in situ experiments. We further determined a six-dimensional free energy surface as a function of the inter-nucleosome distances via a deep learning approach. The zigzag structure can indeed be seen as the global minimum of the surface. However, it is not favored by a significant amount relative to the partially unfolded, in situ configurations. Chemical perturbations such as histone H4 tail acetylation and thermal fluctuations can further tilt the energetic balance to stabilize intermediate states. Our study provides insight into the connection between various reported chromatin configurations and has implications on the in situ relevance of the 30 nm fiber. The three-dimensional organization of chromatin plays critical roles in regulating genome function. Here the authors apply a near atomistic model to study the structure and dynamics of the chromatin folding unit - the tetra-nucleosome - to provide insight into how chromatin folds.
Collapse
Affiliation(s)
- Xinqiang Ding
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA.
| |
Collapse
|
38
|
Lu W, Bueno C, Schafer NP, Moller J, Jin S, Chen X, Chen M, Gu X, Davtyan A, de Pablo JJ, Wolynes PG. OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations. PLoS Comput Biol 2021; 17:e1008308. [PMID: 33577557 PMCID: PMC7906472 DOI: 10.1371/journal.pcbi.1008308] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Revised: 02/25/2021] [Accepted: 01/09/2021] [Indexed: 01/28/2023] Open
Abstract
We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology. The cell’s most important pieces of machinery are large complexes of proteins often along with nucleic acids. From the ribosome, to CRISPR-Cas9, to transcription factors and DNA-wrangling proteins like the SMC-Kleisins, these complexes allow organisms to replicate and enable cells to respond to environmental cues. Computer simulation is a key technology that can be used to connect physical theories with biological reality. Unfortunately, the time and length scales accessible to molecular simulation have not kept pace with our ambition to study the cell’s molecular factories. Many simulation codes also unfortunately remain effectively locked away from the user community who need to modify them as more of the underlying physics is learned. In this paper, we present OpenAWSEM and Open3SPN2, two new easy-to-use and easy to modify implementations of efficient and accurate coarse-grained protein and DNA simulation forcefields that can now be run hundreds of times faster than before, thereby making studies of large biomolecular machines more facile.
Collapse
Affiliation(s)
- Wei Lu
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Department of Physics, Rice University, Houston, Texas, United States of America
| | - Carlos Bueno
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Department of Chemistry, Rice University, Houston, Texas, United States of America
| | - Nicholas P. Schafer
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Department of Chemistry, Rice University, Houston, Texas, United States of America
- Schafer Science, LLC, Houston, Texas United States of America
| | - Joshua Moller
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois, United States of America
- Argonne National Laboratory, Lemont, Illinois, United States of America
| | - Shikai Jin
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Department of Biosciences, Rice University, Houston, Texas, United States of America
| | - Xun Chen
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Department of Chemistry, Rice University, Houston, Texas, United States of America
| | - Mingchen Chen
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
| | - Xinyu Gu
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Department of Chemistry, Rice University, Houston, Texas, United States of America
| | - Aram Davtyan
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
| | - Juan J. de Pablo
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois, United States of America
- Argonne National Laboratory, Lemont, Illinois, United States of America
| | - Peter G. Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, Texas, United States of America
- Department of Chemistry, Rice University, Houston, Texas, United States of America
- Department of Physics, Rice University, Houston, Texas, United States of America
- Department of Biosciences, Rice University, Houston, Texas, United States of America
- * E-mail:
| |
Collapse
|
39
|
Kenzaki H, Takada S. Linker DNA Length is a Key to Tri-nucleosome Folding. J Mol Biol 2020; 433:166792. [PMID: 33383034 DOI: 10.1016/j.jmb.2020.166792] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Revised: 12/07/2020] [Accepted: 12/22/2020] [Indexed: 01/18/2023]
Abstract
The folding of a nucleosome array has long been one of the fundamental and unsolved problems in chromatin biology. In this study, we address how nucleosome array folding depends on the length of linker DNA. We performed molecular dynamics simulations of a tri-nucleosome, a minimal model of chromatin folding, with various linker lengths (LLs) ranging from 20 to 40 base pairs (bps). We found that the tri-nucleosome folding strongly depends on LLs, and classified the structure ensemble into five classes, named from trinuc-1 to trinuc-5. As a function of LL, the different classes appear, on average, every 2 bps with a period of 10 bps, and are characterized by distinct inter-nucleosome interactions. The trinuc-1 conformation corresponds to LL ~ 10n, where n is an integer, and is stabilized by the tight packing between the first and the third nucleosomes, consistent with a zigzag fiber form. Structures of the other four classes are more diverse and distributed continuously in the space of possible configurations. Histone-DNA electrostatic interactions in the tri-nucleosome are further analyzed.
Collapse
Affiliation(s)
- Hiroo Kenzaki
- Information Systems Division, Head Office for Information Systems and Cybersecurity, RIKEN, Hirosawa, Wako, Saitama 351-0198, Japan.
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kitashirakawa, Sakyo, Kyoto 606-8502, Japan.
| |
Collapse
|
40
|
Leicher R, Ge EJ, Lin X, Reynolds MJ, Xie W, Walz T, Zhang B, Muir TW, Liu S. Single-molecule and in silico dissection of the interaction between Polycomb repressive complex 2 and chromatin. Proc Natl Acad Sci U S A 2020; 117:30465-30475. [PMID: 33208532 PMCID: PMC7720148 DOI: 10.1073/pnas.2003395117] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Polycomb repressive complex 2 (PRC2) installs and spreads repressive histone methylation marks on eukaryotic chromosomes. Because of the key roles that PRC2 plays in development and disease, how this epigenetic machinery interacts with DNA and nucleosomes is of major interest. Nonetheless, the mechanism by which PRC2 engages with native-like chromatin remains incompletely understood. In this work, we employ single-molecule force spectroscopy and molecular dynamics simulations to dissect the behavior of PRC2 on polynucleosome arrays. Our results reveal an unexpectedly diverse repertoire of PRC2 binding configurations on chromatin. Besides reproducing known binding modes in which PRC2 interacts with bare DNA, mononucleosomes, and adjacent nucleosome pairs, our data also provide direct evidence that PRC2 can bridge pairs of distal nucleosomes. In particular, the "1-3" bridging mode, in which PRC2 engages two nucleosomes separated by one spacer nucleosome, is a preferred low-energy configuration. Moreover, we show that the distribution and stability of different PRC2-chromatin interaction modes are modulated by accessory subunits, oncogenic histone mutations, and the methylation state of chromatin. Overall, these findings have implications for the mechanism by which PRC2 spreads histone modifications and compacts chromatin. The experimental and computational platforms developed here provide a framework for understanding the molecular basis of epigenetic maintenance mediated by Polycomb-group proteins.
Collapse
Affiliation(s)
- Rachel Leicher
- Laboratory of Nanoscale Biophysics and Biochemistry, The Rockefeller University, New York, NY 10065
- Tri-Institutional PhD Program in Chemical Biology, New York, NY 10065
| | - Eva J Ge
- Department of Chemistry, Princeton University, Princeton, NJ 08544
| | - Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Matthew J Reynolds
- Laboratory of Structural Biophysics and Mechanobiology, The Rockefeller University, New York, NY 10065
| | - Wenjun Xie
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Thomas Walz
- Laboratory of Molecular Electron Microscopy, The Rockefeller University, New York, NY 10065
| | - Bin Zhang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139;
| | - Tom W Muir
- Department of Chemistry, Princeton University, Princeton, NJ 08544
| | - Shixin Liu
- Laboratory of Nanoscale Biophysics and Biochemistry, The Rockefeller University, New York, NY 10065;
| |
Collapse
|
41
|
Mondal A, Bhattacherjee A. Mechanism of Dynamic Binding of Replication Protein A to ssDNA. J Chem Inf Model 2020; 60:5057-5069. [PMID: 32990435 DOI: 10.1021/acs.jcim.0c00564] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Replication protein A (RPA) serves as a hub protein inside eukaryotic cells, where it coordinates crucial DNA metabolic processes and activates the DNA-damage response system. A characteristic feature of its action is to associate with single-stranded DNA (ssDNA) intermediates before handing them over to downstream proteins. The length of ssDNA intermediates differs for different pathways. This means that RPA must have mechanisms for selective processing of ssDNA intermediates based on their length, the knowledge of which is fundamental to elucidate when and how DNA repair and replication processes are symphonized. By employing extensive molecular dynamics simulations, we investigated the mechanism of binding of RPA to ssDNA of different lengths. We show that the binding involves dynamic equilibrium with a stable intermediate, the population of which increases with the length of ssDNA. The vital underlying factors are decoded through collective variable principal component analysis. It suggests a differently orchestrated set of interactions that define the action of RPA based on the length of ssDNA intermediates. We further estimated the association kinetics that matches excellently well with previous experimental studies and probed the diffusion mechanism of RPA to ssDNA. RPA diffuses on short ssDNA through progressive "bulge" formation. With long ssDNA, we observed a conformational change in ssDNA coupled with its binding to RPA in a cooperative fashion. This unanticipated binding mechanism successfully explains how the "short-lived", long ssDNA intermediates are processed quickly in vivo. This study thus reveals the molecular basis of several recent experimental observations related to RPA binding to ssDNA and provides novel insights into the RPA functioning in DNA repair and replication.
Collapse
Affiliation(s)
- Anupam Mondal
- 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
|
42
|
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]
|
43
|
Abstract
While recent experiments revealed that some pioneer transcription factors (TFs) can bind to their target DNA sequences inside a nucleosome, the binding dynamics of their target recognitions are poorly understood. Here we used the latest coarse-grained models and molecular dynamics simulations to study the nucleosome-binding procedure of the two pioneer TFs, Sox2 and Oct4. In the simulations for a strongly positioning nucleosome, Sox2 selected its target DNA sequence only when the target was exposed. Otherwise, Sox2 entropically bound to the dyad region nonspecifically. In contrast, Oct4 plastically bound on the nucleosome mainly in two ways. First, the two POU domains of Oct4 separately bound to the two parallel gyres of the nucleosomal DNA, supporting the previous experimental results of the partial motif recognition. Second, the POUS domain of Oct4 favored binding on the acidic patch of histones. Then, simulating the TFs binding to a genomic nucleosome, the LIN28B nucleosome, we found that the recognition of a pseudo motif by Sox2 induced the local DNA bending and shifted the population of the rotational position of the nucleosomal DNA. The redistributed DNA phase, in turn, changed the accessibility of a distant TF binding site, which consequently affected the binding probability of a second Sox2 or Oct4. These results revealed a nucleosomal DNA-mediated allosteric mechanism, through which one TF binding event can change the global conformation, and effectively regulate the binding of another TF at distant sites. Our simulations provide insights into the binding mechanism of single and multiple TFs on the nucleosome.
Collapse
|
44
|
Peng J, Yuan C, Hua X, Zhang Z. Molecular mechanism of histone variant H2A.B on stability and assembly of nucleosome and chromatin structures. Epigenetics Chromatin 2020; 13:28. [PMID: 32664941 PMCID: PMC7362417 DOI: 10.1186/s13072-020-00351-x] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2019] [Accepted: 07/09/2020] [Indexed: 01/23/2023] Open
Abstract
BACKGROUND H2A.B, the most divergent histone variant of H2A, can significantly modulate nucleosome and chromatin structures. However, the related structural details and the underlying mechanism remain elusive to date. In this work, we built atomic models of the H2A.B-containing nucleosome core particle (NCP), chromatosome, and chromatin fiber. Multiscale modeling including all-atom molecular dynamics and coarse-grained simulations were then carried out for these systems. RESULTS It is found that sequence differences at the C-terminal tail, the docking domain, and the L2 loop, between H2A.B and H2A are directly responsible for the DNA unwrapping in the H2A.B NCP, whereas the N-terminus of H2A.B may somewhat compensate for the aforementioned unwrapping effect. The assembly of the H2A.B NCP is more difficult than that of the H2A NCP. H2A.B may also modulate the interactions of H1 with both the NCP and the linker DNA and could further affect the higher-order structure of the chromatin fiber. CONCLUSIONS The results agree with the experimental results and may shed new light on the biological function of H2A.B. Multiscale modeling may be a valuable tool for investigating structure and dynamics of the nucleosome and the chromatin induced by various histone variants.
Collapse
Affiliation(s)
- Junhui Peng
- MOE Key Laboratory for Membraneless Organelles & Cellular Dynamics, National Science Center for Physical Sciences at Microscale, School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, 230026, People's Republic of China.,Laboratory of Evolutionary Genetics and Genomics, The Rockefeller University, New York, NY, 10065, USA
| | - Chuang Yuan
- MOE Key Laboratory for Membraneless Organelles & Cellular Dynamics, National Science Center for Physical Sciences at Microscale, School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, 230026, People's Republic of China
| | - Xinfan Hua
- MOE Key Laboratory for Membraneless Organelles & Cellular Dynamics, National Science Center for Physical Sciences at Microscale, School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, 230026, People's Republic of China
| | - Zhiyong Zhang
- MOE Key Laboratory for Membraneless Organelles & Cellular Dynamics, National Science Center for Physical Sciences at Microscale, School of Life Sciences, University of Science and Technology of China, Hefei, Anhui, 230026, People's Republic of China.
| |
Collapse
|
45
|
Jin S, Chen M, Chen X, Bueno C, Lu W, Schafer NP, Lin X, Onuchic JN, Wolynes PG. Protein Structure Prediction in CASP13 Using AWSEM-Suite. J Chem Theory Comput 2020; 16:3977-3988. [PMID: 32396727 DOI: 10.1021/acs.jctc.0c00188] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Recently several techniques have emerged that significantly enhance the quality of predictions of protein tertiary structures. In this study, we describe the performance of AWSEM-Suite, an algorithm that incorporates template-based modeling and coevolutionary restraints with a realistic coarse-grained force field, AWSEM. With its roots in neural networks, AWSEM contains both physical and bioinformatical energies that have been optimized using energy landscape theory. AWSEM-Suite participated in CASP13 as a server predictor and generated reliable predictions for most targets. AWSEM-Suite ranked eighth in both the free-modeling category and the hard-to-model category and in one case provided the best submitted prediction. Here we critically discuss the prediction performance of AWSEM-Suite using several examples from different categories in CASP13. Structure prediction tests on these selected targets, two of them being hard-to-model targets, show that AWSEM-Suite can achieve high-resolution structure prediction after incorporating both template guidances and coevolutionary restraints even when homology is weak. For targets with reliable templates (template-easy category), introducing coevolutionary restraints sometimes damages the overall quality of the predictions. Free energy profile analyses demonstrate, however, that the incorporations of both of these evolutionarily informed terms effectively increase the funneling of the landscape toward native-like structures while still allowing sufficient flexibility to correct for discrepancies between the correct target structure and the provided guidance. In contrast to other predictors that are exclusively oriented toward structure prediction, the connection of AWSEM-Suite to a statistical mechanical basis and affiliated molecular dynamics and importance sampling simulations makes it suitable for functional explorations.
Collapse
Affiliation(s)
| | | | - Xun Chen
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | | | - Wei Lu
- Department of Physics, Rice University, Houston, Texas 77005, United States
| | | | - Xingcheng Lin
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - José N Onuchic
- Department of Chemistry, Rice University, Houston, Texas 77005, United States.,Department of Physics, Rice University, Houston, Texas 77005, United States
| | - Peter G Wolynes
- Department of Chemistry, Rice University, Houston, Texas 77005, United States.,Department of Physics, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
46
|
Moller J, de Pablo JJ. Bottom-Up Meets Top-Down: The Crossroads of Multiscale Chromatin Modeling. Biophys J 2020; 118:2057-2065. [PMID: 32320675 DOI: 10.1016/j.bpj.2020.03.014] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2019] [Revised: 03/04/2020] [Accepted: 03/20/2020] [Indexed: 01/19/2023] Open
Abstract
Chromatin can be viewed as a hierarchically structured fiber that regulates gene expression. It consists of a complex network of DNA and proteins whose characteristic dynamical modes facilitate compaction and rearrangement in the cell nucleus. These modes stem from chromatin's fundamental unit, the nucleosome, and their effects are propagated across length scales. Understanding the effects of nucleosome dynamics on the chromatin fiber, primarily through post-translational modifications that occur on the histones, is of central importance to epigenetics. Within the last decade, imaging and chromosome conformation capture techniques have revealed a number of structural and statistical features of the packaged chromatin fiber at a hitherto unavailable level of resolution. Such experiments have led to increased efforts to develop polymer models that aim to reproduce, explain, and predict the contact probability scaling and density heterogeneity. At nanometer scales, available models have focused on the role of the nucleosome and epigenetic marks on local chromatin structure. At micrometer scales, existing models have sought to explain scaling laws and density heterogeneity. Less work, however, has been done to reconcile these two approaches: bottom-up and top-down models of chromatin. In this perspective, we highlight the multiscale simulation models that are driving toward an understanding of chromatin structure and function, from the nanometer to the micron scale, and we highlight areas of opportunity and some of the prospects for new frameworks that bridge these two scales. Taken together, experimental and modeling advances over the last few years have established a robust platform for the study of chromatin fiber structure and dynamics, which will be of considerable use to the chromatin community in developing an understanding of the interplay between epigenomic regulation and molecular structure.
Collapse
Affiliation(s)
- Joshua Moller
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois
| | - Juan J de Pablo
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois; Material Science Division, Argonne National Laboratory, Lemont, Illinois.
| |
Collapse
|
47
|
Chromatin Compaction Leads to a Preference for Peripheral Heterochromatin. Biophys J 2020; 118:1479-1488. [PMID: 32097622 DOI: 10.1016/j.bpj.2020.01.034] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2019] [Revised: 01/23/2020] [Accepted: 01/24/2020] [Indexed: 12/14/2022] Open
Abstract
A layer of dense heterochromatin is found at the periphery of the nucleus. Because this peripheral heterochromatin functions as a repressive phase, mechanisms that relocate genes to the periphery play an important role in regulating transcription. Using Monte Carlo simulations, we show that an interaction that attracts euchromatin and heterochromatin equally to the nuclear envelope will still preferentially locate heterochromatin to the nuclear periphery. This observation considerably broadens the class of possible interactions that result in peripheral positioning to include boundary interactions that either weakly attract all chromatin or strongly bind to a randomly chosen 0.05% of nucleosomes. The key distinguishing feature of heterochromatin is its high chromatin density with respect to euchromatin. In our model, this densification is caused by heterochromatin protein 1's preferential binding to histone H3 tails with a methylated lysine at the ninth residue, a hallmark of heterochromatin. We find that a global rearrangement of chromatin to place heterochromatin at the nuclear periphery can be accomplished by attaching a small subset of loci, even if these loci are uncorrelated with heterochromatin. Hence, factors that densify chromatin determine which genomic regions condense to form peripheral heterochromatin.
Collapse
|
48
|
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
|
49
|
Takada S, Brandani GB, Tan C. Nucleosomes as allosteric scaffolds for genetic regulation. Curr Opin Struct Biol 2020; 62:93-101. [PMID: 31901887 DOI: 10.1016/j.sbi.2019.11.013] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Revised: 11/13/2019] [Accepted: 11/14/2019] [Indexed: 12/11/2022]
Abstract
Nucleosomes are stable yet highly dynamic complexes exhibiting diverse types of motions, such as sliding, DNA unwrapping, and disassembly, encoding a landscape with a large number of metastable states. In this review, describing recent studies on these nucleosome structure changes, we propose that the nucleosome can be viewed as an ideal allosteric scaffold: regulated by effector molecules such as transcription factors and chromatin remodelers, the nucleosome controls the downstream gene activity. Binding of transcription factors to the nucleosome can enhance DNA unwrapping or slide the DNA, altering either the binding or the unbinding of other transcription factors to nearby sites. ATP-dependent chromatin remodelers induce a series of DNA deformations, which allosterically propagate throughout the nucleosome to induce DNA sliding or histone exchange.
Collapse
Affiliation(s)
- Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwake, Sakyo Kyoto, 606-8502, Japan.
| | - Giovanni B Brandani
- Department of Biophysics, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwake, Sakyo Kyoto, 606-8502, Japan
| | - Cheng Tan
- RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo, Kobe, 650-0047 Japan
| |
Collapse
|
50
|
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
|