1
|
Kunze T, Dreßler C, Lauer C, Paul W, Sebastiani D. Reverse Mapping of Coarse Grained Polyglutamine Conformations from PRIME20 Sampling. Chemphyschem 2024; 25:e202300521. [PMID: 38314956 DOI: 10.1002/cphc.202300521] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Revised: 02/01/2024] [Accepted: 02/02/2024] [Indexed: 02/07/2024]
Abstract
An inverse coarse-graining protocol is presented for generating and validating atomistic structures of large (bio-) molecules from conformations obtained via a coarse-grained sampling method. Specifically, the protocol is implemented and tested based on the (coarse-grained) PRIME20 protein model (P20/SAMC), and the resulting all-atom conformations are simulated using conventional biomolecular force fields. The phase space sampling at the coarse-grained level is performed with a stochastical approximation Monte Carlo approach. The method is applied to a series of polypeptides, specifically dimers of polyglutamine with varying chain length in aqueous solution. The majority (>70 %) of the conformations obtained from the coarse-grained peptide model can successfully be mapped back to atomistic structures that remain conformationally stable during 10 ns of molecular dynamics simulations. This work can be seen as the first step towards the overarching goal of improving our understanding of protein aggregation phenomena through simulation methods.
Collapse
Affiliation(s)
- Thomas Kunze
- Faculty of Natural Sciences II, Martin-Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120, Halle, Germany
| | - Christian Dreßler
- Institut für Physik, Ilmenau University of Technology, Weimarer Straße 32, 98693, Ilmenau, Germany
| | - Christian Lauer
- Faculty of Natural Sciences II, Martin-Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120, Halle, Germany
| | - Wolfgang Paul
- Faculty of Natural Sciences II, Martin-Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120, Halle, Germany
| | - Daniel Sebastiani
- Faculty of Natural Sciences II, Martin-Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120, Halle, Germany
| |
Collapse
|
2
|
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
|
3
|
Nishimura M, Fujii T, Tanaka H, Maehara K, Morishima K, Shimizu M, Kobayashi Y, Nozawa K, Takizawa Y, Sugiyama M, Ohkawa Y, Kurumizaka H. Genome-wide mapping and cryo-EM structural analyses of the overlapping tri-nucleosome composed of hexasome-hexasome-octasome moieties. Commun Biol 2024; 7:61. [PMID: 38191828 PMCID: PMC10774305 DOI: 10.1038/s42003-023-05694-1] [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] [Received: 03/26/2023] [Accepted: 12/11/2023] [Indexed: 01/10/2024] Open
Abstract
The nucleosome is a fundamental unit of chromatin in which about 150 base pairs of DNA are wrapped around a histone octamer. The overlapping di-nucleosome has been proposed as a product of chromatin remodeling around the transcription start site, and previously found as a chromatin unit, in which about 250 base pairs of DNA continuously bind to the histone core composed of a hexamer and an octamer. In the present study, our genome-wide analysis of human cells suggests another higher nucleosome stacking structure, the overlapping tri-nucleosome, which wraps about 300-350 base-pairs of DNA in the region downstream of certain transcription start sites of actively transcribed genes. We determine the cryo-electron microscopy (cryo-EM) structure of the overlapping tri-nucleosome, in which three subnucleosome moieties, hexasome, hexasome, and octasome, are associated by short connecting DNA segments. Small angle X-ray scattering and coarse-grained molecular dynamics simulation analyses reveal that the cryo-EM structure of the overlapping tri-nucleosome may reflect its structure in solution. Our findings suggest that nucleosome stacking structures composed of hexasome and octasome moieties may be formed by nucleosome remodeling factors around transcription start sites for gene regulation.
Collapse
Affiliation(s)
- Masahiro Nishimura
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo, 113-0032, Japan
- Epigenetics and Stem Cell Biology Laboratory, National Institute of Environmental Health Sciences, 111 TW, Alexander Drive, Research Triangle Park, NC, 27707, USA
| | - Takeru Fujii
- Division of Transcriptomics, Medical Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi, Higashi, Fukuoka, 812-0054, Japan
| | - Hiroki Tanaka
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo, 113-0032, Japan
- Department of Structural Virology, National Center for Global Health and Medicine, 1-21-1 Toyama, Shinjuku-ku, Tokyo, 162-8655, Japan
| | - Kazumitsu Maehara
- Division of Transcriptomics, Medical Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi, Higashi, Fukuoka, 812-0054, Japan
| | - Ken Morishima
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, Kumatori, Sennan-gun, Osaka, 590-0494, Japan
| | - Masahiro Shimizu
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, Kumatori, Sennan-gun, Osaka, 590-0494, Japan
| | - Yuki Kobayashi
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo, 113-0032, Japan
| | - Kayo Nozawa
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo, 113-0032, Japan
- School of Life Science and Technology, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa, 226-8501, Japan
| | - Yoshimasa Takizawa
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo, 113-0032, Japan
| | - Masaaki Sugiyama
- Institute for Integrated Radiation and Nuclear Science, Kyoto University, Kumatori, Sennan-gun, Osaka, 590-0494, Japan
| | - Yasuyuki Ohkawa
- Division of Transcriptomics, Medical Institute of Bioregulation, Kyushu University, 3-1-1 Maidashi, Higashi, Fukuoka, 812-0054, Japan.
| | - Hitoshi Kurumizaka
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo, 113-0032, Japan.
| |
Collapse
|
4
|
Kim S. Backmapping with Mapping and Isomeric Information. J Phys Chem B 2023. [PMID: 38049145 DOI: 10.1021/acs.jpcb.3c05593] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/06/2023]
Abstract
I present a powerful and flexible backmapping tool named Multiscale Simulation Tool (mstool) that converts a coarse-grained (CG) system into all-atom (AA) resolution and only requires AA to CG mapping and isomeric information (cis/trans/dihedral/chiral). The backmapping procedure includes two simple steps: (a) AA atoms are randomly placed near the corresponding CG beads according to the provided mapping scheme. (b) Energy minimization is performed with two modifications in the AA force field (FF). First, nonbonded interactions are replaced with cosine functions to ensure the numerical stability. Second, additional torsions are imposed to maintain the molecules' isomeric properties. To test the simplicity and robustness of the tool, I backmapped multiple membrane and protein CG structures into AA resolution, including a four-bead CG lipid model (resolution increased by a factor of 34) without using intermediate resolution. The tool is freely available at github.com/ksy141/mstool.
Collapse
Affiliation(s)
- Siyoung Kim
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637 United States
| |
Collapse
|
5
|
Laurent H, Hughes MDG, Walko M, Brockwell DJ, Mahmoudi N, Youngs TGA, Headen TF, Dougan L. Visualization of Self-Assembly and Hydration of a β-Hairpin through Integrated Small and Wide-Angle Neutron Scattering. Biomacromolecules 2023; 24:4869-4879. [PMID: 37874935 PMCID: PMC10646990 DOI: 10.1021/acs.biomac.3c00583] [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: 06/14/2023] [Revised: 10/03/2023] [Indexed: 10/26/2023]
Abstract
Fundamental understanding of the structure and assembly of nanoscale building blocks is crucial for the development of novel biomaterials with defined architectures and function. However, accessing self-consistent structural information across multiple length scales is challenging. This limits opportunities to exploit atomic scale interactions to achieve emergent macroscale properties. In this work we present an integrative small- and wide-angle neutron scattering approach coupled with computational modeling to reveal the multiscale structure of hierarchically self-assembled β hairpins in aqueous solution across 4 orders of magnitude in length scale from 0.1 Å to 300 nm. Our results demonstrate the power of this self-consistent cross-length scale approach and allows us to model both the large-scale self-assembly and small-scale hairpin hydration of the model β hairpin CLN025. Using this combination of techniques, we map the hydrophobic/hydrophilic character of this model self-assembled biomolecular surface with atomic resolution. These results have important implications for the multiscale investigation of aqueous peptides and proteins, for the prediction of ligand binding and molecular associations for drug design, and for understanding the self-assembly of peptides and proteins for functional biomaterials.
Collapse
Affiliation(s)
- Harrison Laurent
- School
of Physics and Astronomy, University of
Leeds, Leeds, United Kingdom, LS2
9JT
| | - Matt D. G. Hughes
- School
of Physics and Astronomy, University of
Leeds, Leeds, United Kingdom, LS2
9JT
- Astbury
Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom LS2
9JT
| | - Martin Walko
- School
of Chemistry, University of Leeds, Leeds, United
Kingdom, LS2 9JT
| | - David J. Brockwell
- Astbury
Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom LS2
9JT
| | - Najet Mahmoudi
- ISIS
Neutron and Muon Source, Rutherford Appleton
Laboratory, Harwell Oxford, Didcot, United Kingdom, OX11 0QX
| | - Tristan G. A. Youngs
- ISIS
Neutron and Muon Source, Rutherford Appleton
Laboratory, Harwell Oxford, Didcot, United Kingdom, OX11 0QX
| | - Thomas F. Headen
- ISIS
Neutron and Muon Source, Rutherford Appleton
Laboratory, Harwell Oxford, Didcot, United Kingdom, OX11 0QX
| | - Lorna Dougan
- School
of Physics and Astronomy, University of
Leeds, Leeds, United Kingdom, LS2
9JT
- Astbury
Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom LS2
9JT
| |
Collapse
|
6
|
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
|
7
|
NARall: a novel tool for reconstruction of the all-atom structure of nucleic acids from heavily coarse-grained model. CHEMICAL PAPERS 2022. [DOI: 10.1007/s11696-022-02634-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
AbstractNucleic acids are one of the most important cellular components. These molecules have been studied both experimentally and theoretically. As all-atom simulations are still limited to short time scales, coarse-grain modeling allows to study of those molecules on a longer time scale. Nucleic-Acid united RESidue (NARES-2P) is a low-resolution coarse-grained model with two centers of interaction per repeating unit. It has been successfully applied to study DNA self-assembly and telomeric properties. This force field enables the study of nucleic acids Behavior on a long time scale but lacks atomistic details. In this article, we present new software to reconstruct atomistic details from the NARES-2P model. It has been applied to RNA pseudoknot, nucleic acid four-way junction, G-quadruplex and DNA duplex converted to NARES-2P model and back. Moreover, it was applied to DNA structure folded and self-assembled with NARES-2P.
Collapse
|
8
|
Brandani GB, Gopi S, Yamauchi M, Takada S. Molecular dynamics simulations for the study of chromatin biology. Curr Opin Struct Biol 2022; 77:102485. [PMID: 36274422 DOI: 10.1016/j.sbi.2022.102485] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2022] [Revised: 09/05/2022] [Accepted: 09/18/2022] [Indexed: 12/14/2022]
Abstract
The organization of Eukaryotic DNA into chromatin has profound implications for the processing of genetic information. In the past years, molecular dynamics (MD) simulations proved to be a powerful tool to investigate the mechanistic basis of chromatin biology. We review recent all-atom and coarse-grained MD studies revealing how the structure and dynamics of chromatin underlie its biological functions. We describe the latest method developments; the structural fluctuations of nucleosomes and the various factors affecting them; the organization of chromatin fibers, with particular emphasis on its liquid-like character; the interactions and dynamics of transcription factors on chromatin; and how chromatin organization is modulated by molecular motors acting on DNA.
Collapse
Affiliation(s)
- Giovanni B Brandani
- Department of Biophysics, Graduate School of Science, Kyoto University, Japan.
| | - Soundhararajan Gopi
- Department of Biophysics, Graduate School of Science, Kyoto University, Japan
| | - Masataka Yamauchi
- Department of Biophysics, Graduate School of Science, Kyoto University, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Japan
| |
Collapse
|
9
|
Tan C, Jung J, Kobayashi C, Torre DUL, Takada S, Sugita Y. Implementation of residue-level coarse-grained models in GENESIS for large-scale molecular dynamics simulations. PLoS Comput Biol 2022; 18:e1009578. [PMID: 35381009 PMCID: PMC9012402 DOI: 10.1371/journal.pcbi.1009578] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 04/15/2022] [Accepted: 03/26/2022] [Indexed: 12/28/2022] Open
Abstract
Residue-level coarse-grained (CG) models have become one of the most popular tools in biomolecular simulations in the trade-off between modeling accuracy and computational efficiency. To investigate large-scale biological phenomena in molecular dynamics (MD) simulations with CG models, unified treatments of proteins and nucleic acids, as well as efficient parallel computations, are indispensable. In the GENESIS MD software, we implement several residue-level CG models, covering structure-based and context-based potentials for both well-folded biomolecules and intrinsically disordered regions. An amino acid residue in protein is represented as a single CG particle centered at the Cα atom position, while a nucleotide in RNA or DNA is modeled with three beads. Then, a single CG particle represents around ten heavy atoms in both proteins and nucleic acids. The input data in CG MD simulations are treated as GROMACS-style input files generated from a newly developed toolbox, GENESIS-CG-tool. To optimize the performance in CG MD simulations, we utilize multiple neighbor lists, each of which is attached to a different nonbonded interaction potential in the cell-linked list method. We found that random number generations for Gaussian distributions in the Langevin thermostat are one of the bottlenecks in CG MD simulations. Therefore, we parallelize the computations with message-passing-interface (MPI) to improve the performance on PC clusters or supercomputers. We simulate Herpes simplex virus (HSV) type 2 B-capsid and chromatin models containing more than 1,000 nucleosomes in GENESIS as examples of large-scale biomolecular simulations with residue-level CG models. This framework extends accessible spatial and temporal scales by multi-scale simulations to study biologically relevant phenomena, such as genome-scale chromatin folding or phase-separated membrane-less condensations.
Collapse
Affiliation(s)
- Cheng Tan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | - Jaewoon Jung
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
| | - Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | - Diego Ugarte La Torre
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Yuji Sugita
- Computational Biophysics Research Team, RIKEN Center for Computational Science, Kobe, Hyogo, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo, Japan
- * E-mail:
| |
Collapse
|
10
|
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
|
11
|
Brandani GB, Tan C, Takada S. The kinetic landscape of nucleosome assembly: A coarse-grained molecular dynamics study. PLoS Comput Biol 2021; 17:e1009253. [PMID: 34314440 PMCID: PMC8345847 DOI: 10.1371/journal.pcbi.1009253] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 08/06/2021] [Accepted: 07/06/2021] [Indexed: 11/23/2022] Open
Abstract
The organization of nucleosomes along the Eukaryotic genome is maintained over time despite disruptive events such as replication. During this complex process, histones and DNA can form a variety of non-canonical nucleosome conformations, but their precise molecular details and roles during nucleosome assembly remain unclear. In this study, employing coarse-grained molecular dynamics simulations and Markov state modeling, we characterized the complete kinetics of nucleosome assembly. On the nucleosome-positioning 601 DNA sequence, we observe a rich transition network among various canonical and non-canonical tetrasome, hexasome, and nucleosome conformations. A low salt environment makes nucleosomes stable, but the kinetic landscape becomes more rugged, so that the system is more likely to be trapped in off-pathway partially assembled intermediates. Finally, we find that the co-operativity between DNA bending and histone association enables positioning sequence motifs to direct the assembly process, with potential implications for the dynamic organization of nucleosomes on real genomic sequences. Nucleosomes are biomolecular complexes formed by DNA wrapped around histone proteins. They represent the basic units of Eukaryotic chromosomes, compacting the genome so that it fits into the small nucleus, and regulating important biological processes such as gene expression. Nucleosomes are disassembled during disruptive events such as DNA replication, and re-assembled afterwards to preserve the correct organization of chromatin. However, the molecular details of nucleosome assembly are still not well understood. In particular, experiments found that histones and DNA may associate into a variety of non-canonical complexes, but their precise conformation and role during assembly remain unclear. In this study, we addressed these problems by performing extensive molecular dynamics simulations of nucleosomes undergoing assembly and disassembly. The simulations reveal many insights into the kinetics of assembly, the structure of non-canonical nucleosome intermediates, and the influence of salt concentration and DNA sequence on the assembly process.
Collapse
Affiliation(s)
- Giovanni B. Brandani
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
- * E-mail: (GBB); (ST)
| | - Cheng Tan
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
- * E-mail: (GBB); (ST)
| |
Collapse
|
12
|
Koide H, Kodera N, Bisht S, Takada S, Terakawa T. Modeling of DNA binding to the condensin hinge domain using molecular dynamics simulations guided by atomic force microscopy. PLoS Comput Biol 2021; 17:e1009265. [PMID: 34329301 PMCID: PMC8357123 DOI: 10.1371/journal.pcbi.1009265] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Revised: 08/11/2021] [Accepted: 07/10/2021] [Indexed: 11/19/2022] Open
Abstract
The condensin protein complex compacts chromatin during mitosis using its DNA-loop extrusion activity. Previous studies proposed scrunching and loop-capture models as molecular mechanisms for the loop extrusion process, both of which assume the binding of double-strand (ds) DNA to the hinge domain formed at the interface of the condensin subunits Smc2 and Smc4. However, how the hinge domain contacts dsDNA has remained unknown. Here, we conducted atomic force microscopy imaging of the budding yeast condensin holo-complex and used this data as basis for coarse-grained molecular dynamics simulations to model the hinge structure in a transient open conformation. We then simulated the dsDNA binding to open and closed hinge conformations, predicting that dsDNA binds to the outside surface when closed and to the outside and inside surfaces when open. Our simulations also suggested that the hinge can close around dsDNA bound to the inside surface. Based on these simulation results, we speculate that the conformational change of the hinge domain might be essential for the dsDNA binding regulation and play roles in condensin-mediated DNA-loop extrusion.
Collapse
Affiliation(s)
- Hiroki Koide
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Noriyuki Kodera
- Nano Life Science Institute (WPI-NanoLSI), Kanazawa University, Kanazawa, Japan
| | - Shveta Bisht
- Cell Biology and Biophysics Unit, Structural and Computational Unit, European Molecular Biology Laboratory (EMBL), Heidelberg, Germany
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
| | - Tsuyoshi Terakawa
- Department of Biophysics, Graduate School of Science, Kyoto University, Kyoto, Japan
- PREST, Japan Science and Technology Agency (JST), Kawaguchi, Japan
| |
Collapse
|
13
|
Thaler S, Praprotnik M, Zavadlav J. Back-mapping augmented adaptive resolution simulation. J Chem Phys 2020; 153:164118. [PMID: 33138420 DOI: 10.1063/5.0025728] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Concurrent multiscale techniques such as Adaptive Resolution Scheme (AdResS) can offer ample computational advantages over conventional atomistic (AT) molecular dynamics simulations. However, they typically rely on aphysical hybrid regions to maintain numerical stability when high-resolution degrees of freedom (DOFs) are randomly re-inserted at the resolution interface. We propose an Energy Minimized AT (DOF) Insertion (EMATI) method that uses an informed rather than random AT DOF insertion to tackle the root cause of the issue, i.e., overlapping AT potentials. EMATI enables us to directly couple AT and coarse-grained resolutions without any modifications of the interaction potentials. We exemplify AdResS-EMATI in a system of liquid butane and show that it yields improved structural and thermodynamic properties at the interface compared to competing AdResS approaches. Furthermore, our approach extends the applicability of the AdResS without a hybrid region to systems for which force capping is inadequate.
Collapse
Affiliation(s)
- S Thaler
- Professorship of Multiscale Modeling of Fluid Materials, Department of Mechanical Engineering, Technical University of Munich, Munich, Germany
| | - M Praprotnik
- Laboratory for Molecular Modeling, National Institute of Chemistry, SI-1001 Ljubljana, Slovenia
| | - J Zavadlav
- Professorship of Multiscale Modeling of Fluid Materials, Department of Mechanical Engineering, Technical University of Munich, Munich, Germany
| |
Collapse
|
14
|
Tokuhisa A, Kanada R, Chiba S, Terayama K, Isaka Y, Ma B, Kamiya N, Okuno Y. Coarse-Grained Diffraction Template Matching Model to Retrieve Multiconformational Models for Biomolecule Structures from Noisy Diffraction Patterns. J Chem Inf Model 2020; 60:2803-2818. [PMID: 32469517 DOI: 10.1021/acs.jcim.0c00131] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Biomolecular imaging using X-ray free-electron lasers (XFELs) has been successfully applied to serial femtosecond crystallography. However, the application of single-particle analysis for structure determination using XFELs with 100 nm or smaller biomolecules has two practical problems: the incomplete diffraction data sets for reconstructing 3D assembled structures and the heterogeneous conformational states of samples. A new diffraction template matching method is thus presented here to retrieve a plausible 3D structural model based on single noisy target diffraction patterns, assuming candidate structures. Two concepts are introduced here: prompt candidate diffraction, generated by enhanced sampled coarse-grain (CG) candidate structures, and efficient molecular orientation searching for matching based on Bayesian optimization. A CG model-based diffraction-matching protocol is proposed that achieves a 100-fold speed increase compared to exhaustive diffraction matching using an all-atom model. The conditions that enable multiconformational analysis were also investigated by simulated diffraction data for various conformational states of chromatin and ribosomes. The proposed method can enable multiconformational analysis, with a structural resolution of at least 20 Å for 270-800 Å flexible biomolecules, in experimental single-particle structure analyses that employ XFELs.
Collapse
Affiliation(s)
- Atsushi Tokuhisa
- RIKEN Cluster for Science and Technology Hub, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN Center for Computational Science, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN Medical Sciences Innovation Hub Program, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Ryo Kanada
- RIKEN Cluster for Science and Technology Hub, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Shuntaro Chiba
- RIKEN Medical Sciences Innovation Hub Program, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Kei Terayama
- RIKEN Medical Sciences Innovation Hub Program, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan.,RIKEN Center for Advanced Intelligence Project, 1-4-1 Nihombashi, Chuo-ku, Tokyo 103-0027, Japan.,Graduate School of Medicine, Kyoto University, Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| | - Yuta Isaka
- RIKEN Cluster for Science and Technology Hub, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Center for Cluster Development and Coordination, Foundation for Biomedical Research and Innovation at Kobe, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Biao Ma
- RIKEN Cluster for Science and Technology Hub, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Center for Cluster Development and Coordination, Foundation for Biomedical Research and Innovation at Kobe, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Narutoshi Kamiya
- Graduate School of Simulation Studies, University of Hyogo, 7-1-28, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yasushi Okuno
- RIKEN Cluster for Science and Technology Hub, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN Medical Sciences Innovation Hub Program, 1-7-22, Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan.,Graduate School of Medicine, Kyoto University, Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan.,Center for Cluster Development and Coordination, Foundation for Biomedical Research and Innovation at Kobe, 6-3-5, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
15
|
Roel-Touris J, Bonvin AM. Coarse-grained (hybrid) integrative modeling of biomolecular interactions. Comput Struct Biotechnol J 2020; 18:1182-1190. [PMID: 32514329 PMCID: PMC7264466 DOI: 10.1016/j.csbj.2020.05.002] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2020] [Revised: 04/23/2020] [Accepted: 05/06/2020] [Indexed: 12/23/2022] Open
Abstract
The computational modeling field has vastly evolved over the past decades. The early developments of simplified protein systems represented a stepping stone towards establishing more efficient approaches to sample intricated conformational landscapes. Downscaling the level of resolution of biomolecules to coarser representations allows for studying protein structure, dynamics and interactions that are not accessible by classical atomistic approaches. The combination of different resolutions, namely hybrid modeling, has also been proved as an alternative when mixed levels of details are required. In this review, we provide an overview of coarse-grained/hybrid models focusing on their applicability in the modeling of biomolecular interactions. We give a detailed list of ready-to-use modeling software for studying biomolecular interactions allowing various levels of coarse-graining and provide examples of complexes determined by integrative coarse-grained/hybrid approaches in combination with experimental information.
Collapse
|
16
|
An Y, Deshmukh SA. Machine learning approach for accurate backmapping of coarse-grained models to all-atom models. Chem Commun (Camb) 2020; 56:9312-9315. [PMID: 32667366 DOI: 10.1039/d0cc02651d] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Four different machine learning (ML) regression models: artificial neural network, k-nearest neighbors, Gaussian process regression and random forest were built to backmap coarse-grained models to all-atom models.
Collapse
Affiliation(s)
- Yaxin An
- Department of Chemical Engineering
- Virginia Tech
- Blacksburg
- USA
| | | |
Collapse
|
17
|
Badaczewska-Dawid AE, Kolinski A, Kmiecik S. Computational reconstruction of atomistic protein structures from coarse-grained models. Comput Struct Biotechnol J 2019; 18:162-176. [PMID: 31969975 PMCID: PMC6961067 DOI: 10.1016/j.csbj.2019.12.007] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Accepted: 12/10/2019] [Indexed: 01/02/2023] Open
Abstract
Three-dimensional protein structures, whether determined experimentally or theoretically, are often too low resolution. In this mini-review, we outline the computational methods for protein structure reconstruction from incomplete coarse-grained to all atomistic models. Typical reconstruction schemes can be divided into four major steps. Usually, the first step is reconstruction of the protein backbone chain starting from the C-alpha trace. This is followed by side-chains rebuilding based on protein backbone geometry. Subsequently, hydrogen atoms can be reconstructed. Finally, the resulting all-atom models may require structure optimization. Many methods are available to perform each of these tasks. We discuss the available tools and their potential applications in integrative modeling pipelines that can transfer coarse-grained information from computational predictions, or experiment, to all atomistic structures.
Collapse
Affiliation(s)
| | | | - Sebastian Kmiecik
- Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
18
|
Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications. Int J Mol Sci 2019; 20:ijms20153774. [PMID: 31375023 PMCID: PMC6696403 DOI: 10.3390/ijms20153774] [Citation(s) in RCA: 60] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2019] [Revised: 07/28/2019] [Accepted: 07/30/2019] [Indexed: 12/23/2022] Open
Abstract
Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.
Collapse
|
19
|
Computational approaches to macromolecular interactions in the cell. Curr Opin Struct Biol 2019; 55:59-65. [PMID: 30999240 DOI: 10.1016/j.sbi.2019.03.012] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2018] [Accepted: 03/08/2019] [Indexed: 12/15/2022]
Abstract
Structural modeling of a cell is an evolving strategic direction in computational structural biology. It takes advantage of new powerful modeling techniques, deeper understanding of fundamental principles of molecular structure and assembly, and rapid growth of the amount of structural data generated by experimental techniques. Key modeling approaches to principal types of macromolecular assemblies in a cell already exist. The main challenge, along with the further development of these modeling approaches, is putting them together in a consistent, unified whole cell model. This opinion piece addresses the fundamental aspects of modeling macromolecular assemblies in a cell, and the state-of-the-art in modeling of the principal types of such assemblies.
Collapse
|
20
|
Peng J, Yuan C, Ma R, Zhang Z. Backmapping from Multiresolution Coarse-Grained Models to Atomic Structures of Large Biomolecules by Restrained Molecular Dynamics Simulations Using Bayesian Inference. J Chem Theory Comput 2019; 15:3344-3353. [PMID: 30908042 DOI: 10.1021/acs.jctc.9b00062] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Coarse-grained (CG) simulations have allowed access to larger length scales and longer time scales in the study of the dynamic processes of large biomolecules than all-atom (AA) molecular dynamics (MD) simulations. Backmapping from CG models to AA structures has long been studied because it enables us to gain detailed structure insights from CG simulations. Many methods first construct an AA structure from the CG model by fragments, random placement, or geometrical rules and subsequently optimize the solution via energy minimization, simulated annealing or position-restrained simulations. However, such methods may only work well on residue-level CG models and cannot consider the deviations of CG models. In this work, we describe, to the best of our knowledge, a new backmapping method based on Bayesian inference and restrained MD simulations. Restraints with log harmonic energy terms are defined according to the target CG model using the Bayesian inference in which the CG deviations can be estimated. From an initial AA structure obtained from either high-resolution experiments or homology modeling, a MD simulation with the aforementioned restraints is performed to obtain a final AA structure that is a backmapping of the target CG model. The method was validated using multiresolution CG models of the soluble extracellular region of the human epidermal growth factor receptor and was further applied to construct AA structures from CG simulations of the nucleosome core particle. The results demonstrate that our method can generate accurate AA structures of different types of biomolecules from multiple CG models with either residue-level resolution or much lower resolution than one-site-per-residue.
Collapse
Affiliation(s)
- Junhui Peng
- Hefei National Laboratory for Physical Science at Microscale and School of Life Sciences , University of Science and Technology of China , Hefei , Anhui 230026 , People's Republic of China
| | - Chuang Yuan
- Hefei National Laboratory for Physical Science at Microscale and School of Life Sciences , University of Science and Technology of China , Hefei , Anhui 230026 , People's Republic of China
| | - Rongsheng Ma
- Hefei National Laboratory for Physical Science at Microscale and School of Life Sciences , University of Science and Technology of China , Hefei , Anhui 230026 , People's Republic of China
| | - Zhiyong Zhang
- Hefei National Laboratory for Physical Science at Microscale and School of Life Sciences , University of Science and Technology of China , Hefei , Anhui 230026 , People's Republic of China
| |
Collapse
|
21
|
Advances in coarse-grained modeling of macromolecular complexes. Curr Opin Struct Biol 2018; 52:119-126. [PMID: 30508766 DOI: 10.1016/j.sbi.2018.11.005] [Citation(s) in RCA: 83] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2018] [Revised: 11/05/2018] [Accepted: 11/17/2018] [Indexed: 01/12/2023]
Abstract
Recent progress in coarse-grained (CG) molecular modeling and simulation has facilitated an influx of computational studies on biological macromolecules and their complexes. Given the large separation of length-scales and time-scales that dictate macromolecular biophysics, CG modeling and simulation are well-suited to bridge the microscopic and mesoscopic or macroscopic details observed from all-atom molecular simulations and experiments, respectively. In this review, we first summarize recent innovations in the development of CG models, which broadly include structure-based, knowledge-based, and dynamics-based approaches. We then discuss recent applications of different classes of CG models to explore various macromolecular complexes. Finally, we conclude with an outlook for the future in this ever-growing field of biomolecular modeling.
Collapse
|