1
|
Portillo-Ledesma S, Li Z, Schlick T. Genome modeling: From chromatin fibers to genes. Curr Opin Struct Biol 2023; 78:102506. [PMID: 36577295 PMCID: PMC9908845 DOI: 10.1016/j.sbi.2022.102506] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2022] [Revised: 11/01/2022] [Accepted: 11/06/2022] [Indexed: 12/27/2022]
Abstract
The intricacies of the 3D hierarchical organization of the genome have been approached by many creative modeling studies. The specific model/simulation technique combination defines and restricts the system and phenomena that can be investigated. We present the latest modeling developments and studies of the genome, involving models ranging from nucleosome systems and small polynucleosome arrays to chromatin fibers in the kb-range, chromosomes, and whole genomes, while emphasizing gene folding from first principles. Clever combinations allow the exploration of many interesting phenomena involved in gene regulation, such as nucleosome structure and dynamics, nucleosome-nucleosome stacking, polynucleosome array folding, protein regulation of chromatin architecture, mechanisms of gene folding, loop formation, compartmentalization, and structural transitions at the chromosome and genome levels. Gene-level modeling with full details on nucleosome positions, epigenetic factors, and protein binding, in particular, can in principle be scaled up to model chromosomes and cells to study fundamental biological regulation.
Collapse
Affiliation(s)
- Stephanie Portillo-Ledesma
- Department of Chemistry, New York University, 100 Washington Square East, Silver Building, New York, 10003, NY, USA
| | - Zilong Li
- Department of Chemistry, New York University, 100 Washington Square East, Silver Building, New York, 10003, NY, USA
| | - Tamar Schlick
- Department of Chemistry, New York University, 100 Washington Square East, Silver Building, New York, 10003, NY, USA; Courant Institute of Mathematical Sciences, New York University, 251 Mercer St., New York, 10012, NY, USA; New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Room 340, Geography Building, 3663 North Zhongshan Road, Shanghai, 200122, China; Simons Center for Computational Physical Chemistry, 24 Waverly Place, Silver Building, New York University, New York, 10003, NY, USA.
| |
Collapse
|
2
|
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
|
3
|
Extended ensemble simulations of a SARS-CoV-2 nsp1-5'-UTR complex. PLoS Comput Biol 2022; 18:e1009804. [PMID: 35045069 PMCID: PMC8803185 DOI: 10.1371/journal.pcbi.1009804] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2021] [Revised: 01/31/2022] [Accepted: 01/04/2022] [Indexed: 11/19/2022] Open
Abstract
Nonstructural protein 1 (nsp1) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a 180-residue protein that blocks translation of host mRNAs in SARS-CoV-2-infected cells. Although it is known that SARS-CoV-2’s own RNA evades nsp1’s host translation shutoff, the molecular mechanism underlying the evasion was poorly understood. We performed an extended ensemble molecular dynamics simulation to investigate the mechanism of the viral RNA evasion. Simulation results suggested that the stem loop structure of the SARS-CoV-2 RNA 5’-untranslated region (SL1) binds to both nsp1’s N-terminal globular region and intrinsically disordered region. The consistency of the results was assessed by modeling nsp1-40S ribosome structure based on reported nsp1 experiments, including the X-ray crystallographic structure analysis, the cryo-EM electron density map, and cross-linking experiments. The SL1 binding region predicted from the simulation was open to the solvent, yet the ribosome could interact with SL1. Cluster analysis of the binding mode and detailed analysis of the binding poses suggest residues Arg124, Lys47, Arg43, and Asn126 may be involved in the SL1 recognition mechanism, consistent with the existing mutational analysis. The pandemic of COVID-19 is still rampant all over the world as of 2021 June. SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), the causative pathogen of COVID-19, encodes a protein called nsp1 (nonstructural protein 1), which modulates and hijacks the ribosome of the infected host cells. With nsp1, infected human cells selectively translate SARS-CoV-2’s RNA, which increases the virus reproduction efficiency while evading the host immunity. Though it has been known that nsp1 recognizes characteristic stem-loop structure at 5’-end of SARS-CoV-2’s RNA (called SL1), the molecular mechanism underlying the recognition has been poorly understood. We investigated the mechanism of selective translation using the all-atom molecular dynamics simulation of nsp1-SL1 complex. Our simulation results suggest that the binding between nsp1 and SL1 is multi-modal. The results also imply that both the N-terminal globular part and the C-terminal flexible tail of nsp1 are involved in the binding. The residues involved in nsp1-SL1 binding coincides with the known mutant analyses of SARS-CoV-1 and SARS-CoV-2, as well as experimental evidence about nsp1-ribosome interactions.
Collapse
|
4
|
Ohtomo H, Kurita JI, Sakuraba S, Li Z, Arimura Y, Wakamori M, Tsunaka Y, Umehara T, Kurumizaka H, Kono H, Nishimura Y. The N-terminal Tails of Histones H2A and H2B Adopt Two Distinct Conformations in the Nucleosome with Contact and Reduced Contact to DNA. J Mol Biol 2021; 433:167110. [PMID: 34153285 DOI: 10.1016/j.jmb.2021.167110] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2021] [Revised: 06/11/2021] [Accepted: 06/11/2021] [Indexed: 11/28/2022]
Abstract
The nucleosome comprises two histone dimers of H2A-H2B and one histone tetramer of (H3-H4)2, wrapped around by ~145 bp of DNA. Detailed core structures of nucleosomes have been established by X-ray and cryo-EM, however, histone tails have not been visualized. Here, we have examined the dynamic structures of the H2A and H2B tails in 145-bp and 193-bp nucleosomes using NMR, and have compared them with those of the H2A and H2B tail peptides unbound and bound to DNA. Whereas the H2A C-tail adopts a single but different conformation in both nucleosomes, the N-tails of H2A and H2B adopt two distinct conformations in each nucleosome. To clarify these conformations, we conducted molecular dynamics (MD) simulations, which suggest that the H2A N-tail can locate stably in either the major or minor grooves of nucleosomal DNA. While the H2B N-tail, which sticks out between two DNA gyres in the nucleosome, was considered to adopt two different orientations, one toward the entry/exit side and one on the opposite side. Then, the H2A N-tail minor groove conformation was obtained in the H2B opposite side and the H2B N-tail interacts with DNA similarly in both sides, though more varied conformations are obtained in the entry/exit side. Collectively, the NMR findings and MD simulations suggest that the minor groove conformer of the H2A N-tail is likely to contact DNA more strongly than the major groove conformer, and the H2A N-tail reduces contact with DNA in the major groove when the H2B N-tail is located in the entry/exit side.
Collapse
Affiliation(s)
- Hideaki Ohtomo
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Jun-Ichi Kurita
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Shun Sakuraba
- Institute for Quantum Life Science (iQLS), National Institutes for Quantum and Radiological Science and Technology (QST), 8-1-7 Umemidai, Kizugawa, Kyoto 619-0215, Japan
| | - Zhenhai Li
- Institute for Quantum Life Science (iQLS), National Institutes for Quantum and Radiological Science and Technology (QST), 8-1-7 Umemidai, Kizugawa, Kyoto 619-0215, Japan
| | - Yasuhiro Arimura
- Laboratory of Chromatin Structure and Function, Institute for Quantitative Biosciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan
| | - Masatoshi Wakamori
- Laboratory for Epigenetics Drug Discovery, RIKEN Center for Biosystems Dynamics Research (BDR), 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Yasuo Tsunaka
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Takashi Umehara
- Laboratory for Epigenetics Drug Discovery, RIKEN Center for Biosystems Dynamics Research (BDR), 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, 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
| | - Hidetoshi Kono
- Institute for Quantum Life Science (iQLS), National Institutes for Quantum and Radiological Science and Technology (QST), 8-1-7 Umemidai, Kizugawa, Kyoto 619-0215, Japan
| | - Yoshifumi Nishimura
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan; Graduate School of Integrated Sciences for Life, Hiroshima University, 1-4-4 Kagamiyama, Higashi-Hiroshima 739-8258, Japan.
| |
Collapse
|
5
|
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: 61] [Impact Index Per Article: 20.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
|
6
|
Huertas J, Cojocaru V. Breaths, Twists, and Turns of Atomistic Nucleosomes. J Mol Biol 2020; 433:166744. [PMID: 33309853 DOI: 10.1016/j.jmb.2020.166744] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Revised: 12/02/2020] [Accepted: 12/02/2020] [Indexed: 02/07/2023]
Abstract
Gene regulation programs establish cellular identity and rely on dynamic changes in the structural packaging of genomic DNA. The DNA is packaged in chromatin, which is formed from arrays of nucleosomes displaying different degree of compaction and different lengths of inter-nucleosomal linker DNA. The nucleosome represents the repetitive unit of chromatin and is formed by wrapping 145-147 basepairs of DNA around an octamer of histone proteins. Each of the four histones is present twice and has a structured core and intrinsically disordered terminal tails. Chromatin dynamics are triggered by inter- and intra-nucleosome motions that are controlled by the DNA sequence, the interactions between the histone core and the DNA, and the conformations, positions, and DNA interactions of the histone tails. Understanding chromatin dynamics requires studying all these features at the highest possible resolution. For this, molecular dynamics simulations can be used as a powerful complement or alternative to experimental approaches, from which it is often very challenging to characterize the structural features and atomic interactions controlling nucleosome motions. Molecular dynamics simulations can be performed at different resolutions, by coarse graining the molecular system with varying levels of details. Here we review the successes and the remaining challenges of the application of atomic resolution simulations to study the structure and dynamics of nucleosomes and their complexes with interacting partners.
Collapse
Affiliation(s)
- Jan Huertas
- In Silico Biomolecular Structure and Dynamics Group, Hubrecht Institute, Utrecht, the Netherlands; Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany; Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster, Germany
| | - Vlad Cojocaru
- In Silico Biomolecular Structure and Dynamics Group, Hubrecht Institute, Utrecht, the Netherlands; Department of Cellular and Developmental Biology, Max Planck Institute for Molecular Biomedicine, Münster, Germany; Center for Multiscale Theory and Computation, Westfälische Wilhelms University, Münster, Germany.
| |
Collapse
|
7
|
Kono H, Sakuraba S, Ishida H. Free energy profile for unwrapping outer superhelical turn of CENP-A nucleosome. Biophys Physicobiol 2019; 16:337-343. [PMID: 31984189 PMCID: PMC6975924 DOI: 10.2142/biophysico.16.0_337] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Accepted: 08/05/2019] [Indexed: 01/19/2023] Open
Abstract
Eukaryotic genome is packaged in a nucleus in the form of chromatin. The fundamental structural unit of the chromatin is the protein-DNA complex, nucleosome, where DNA of about 150 bp is wrapped around a histone core almost twice. In cellular processes such as gene expression, DNA repair and duplication, the nucleosomal DNA has to be unwrapped. Histone proteins have their variants, indicating there are a variety of constitutions of nucleosomes. These different constitutions are observed in different cellular processes. To investigate differences among nucleosomes, we calculated free energy profiles for unwrapping the outer superhelical turn of CENP-A nucleosome and compared them with those of the canonical nucleosome. The free energy profiles for CENP-A nucleosome suggest that CENP-A nucleosome is the most stable when 16 to 22 bps are unwrapped in total whereas the canonical nucleosome is the most stable when it is fully wrapped. This indicates that the flexible conformation of CENP-A nucleosome is ready to provide binding sites for the structural integrity of the centromere.
Collapse
Affiliation(s)
- Hidetoshi Kono
- Molecular Modelling and Simulation Group, Institute for Quantum Life Science (iQLS), National Institutes for Quantum and Radiological Science and Technology, Kizugawa, Kyoto 619-0215 Japan
| | - Shun Sakuraba
- Molecular Modelling and Simulation Group, Institute for Quantum Life Science (iQLS), National Institutes for Quantum and Radiological Science and Technology, Kizugawa, Kyoto 619-0215 Japan
| | - Hisashi Ishida
- Molecular Modelling and Simulation Group, Institute for Quantum Life Science (iQLS), National Institutes for Quantum and Radiological Science and Technology, Kizugawa, Kyoto 619-0215 Japan
| |
Collapse
|
8
|
Li D, Ji B. Protein conformational transitions coupling with ligand interactions: Simulations from molecules to medicine. MEDICINE IN NOVEL TECHNOLOGY AND DEVICES 2019. [DOI: 10.1016/j.medntd.2019.100026] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
|