1
|
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
|
2
|
Borges-Araújo L, Patmanidis I, Singh AP, Santos LHS, Sieradzan AK, Vanni S, Czaplewski C, Pantano S, Shinoda W, Monticelli L, Liwo A, Marrink SJ, Souza PCT. Pragmatic Coarse-Graining of Proteins: Models and Applications. J Chem Theory Comput 2023; 19:7112-7135. [PMID: 37788237 DOI: 10.1021/acs.jctc.3c00733] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Ilias Patmanidis
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Akhil P Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
| | - Lucianna H S Santos
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
- Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d'Azur, Inserm, CNRS, 06560 Valbonne, France
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Sergio Pantano
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita, Okayama 700-8530, Japan
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| |
Collapse
|
3
|
Klein F, Soñora M, Helene Santos L, Nazareno Frigini E, Ballesteros-Casallas A, Rodrigo Machado M, Pantano S. The SIRAH force field: A suite for simulations of complex biological systems at the coarse-grained and multiscale levels. J Struct Biol 2023; 215:107985. [PMID: 37331570 DOI: 10.1016/j.jsb.2023.107985] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 05/18/2023] [Accepted: 06/13/2023] [Indexed: 06/20/2023]
Abstract
The different combinations of molecular dynamics simulations with coarse-grained representations have acquired considerable popularity among the scientific community. Especially in biocomputing, the significant speedup granted by simplified molecular models opened the possibility of increasing the diversity and complexity of macromolecular systems, providing realistic insights on large assemblies for more extended time windows. However, a holistic view of biological ensembles' structural and dynamic features requires a self-consistent force field, namely, a set of equations and parameters that describe the intra and intermolecular interactions among moieties of diverse chemical nature (i.e., nucleic and amino acids, lipids, solvent, ions, etc.). Nevertheless, examples of such force fields are scarce in the literature at the fully atomistic and coarse-grained levels. Moreover, the number of force fields capable of handling simultaneously different scales is restricted to a handful. Among those, the SIRAH force field, developed in our group, furnishes a set of topologies and tools that facilitate the setting up and running of molecular dynamics simulations at the coarse-grained and multiscale levels. SIRAH uses the same classical pairwise Hamiltonian function implemented in the most popular molecular dynamics software. In particular, it runs natively in AMBER and Gromacs engines, and porting it to other simulation packages is straightforward. This review describes the underlying philosophy behind the development of SIRAH over the years and across families of biological molecules, discussing current limitations and future implementations.
Collapse
Affiliation(s)
- Florencia Klein
- Laboratoire de Biochimie Théorique, UPR9080, CNRS, Paris, France
| | - Martín Soñora
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay
| | | | - Ezequiel Nazareno Frigini
- Instituto Multidisciplinario de Investigaciones Biológicas de San Luis (IMIBIO-SL), Universidad Nacional de San Luis - CONICET, San Luis, Argentina
| | - Andrés Ballesteros-Casallas
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay; Area Bioinformática, DETEMA, Facultad de Química, Universidad de la República, General Flores 2124, Montevideo, 11600, Uruguay
| | | | - Sergio Pantano
- Institut Pasteur de Montevideo, Mataojo 2020, 11400, Montevideo, Uruguay; Area Bioinformática, DETEMA, Facultad de Química, Universidad de la República, General Flores 2124, Montevideo, 11600, Uruguay.
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
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]
|
7
|
Bae S, Kim JS. Potential of Mean Force for DNA Wrapping Around a Cationic Nanoparticle. J Chem Theory Comput 2021; 17:7952-7961. [PMID: 34792353 DOI: 10.1021/acs.jctc.1c00797] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Sharp bending and wrapping of DNA around proteins and nanoparticles (NPs) has been of extensive research interest. Here, we present the potential of mean force (PMF) for wrapping a DNA double helix around a cationic NP using coarse-grained models of a double-stranded DNA and a cationic NP. Starting from a NP wrapped around by DNA, the PMF was calculated along the distance between the center of the NP and one end of the DNA molecule. A relationship between the distance and the extent of DNA wrapping is used to calculate the PMF as a function of DNA wrapping around a NP. In particular, the PMF was compared for two DNA sequences of (AT)25/(AT)25 and (AC)25/(GT)25, for which the persistence lengths are different by ∼10 nm. The simulation results provide solid evidence of the thermodynamic preference for complex formation of a cationic NP with more flexible DNA over the less flexible DNA. Furthermore, we estimated the elastic energy of DNA bending, which was in good order-of-magnitude agreement with the theoretical prediction of elastic rods. This work suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnologies, in which the position and dynamics of NPs are regulated on large-scale DNA structures, or the structural transformation of DNA is triggered by the sequence-dependent binding of NPs.
Collapse
Affiliation(s)
- Sehui Bae
- Department of Chemistry and Nanoscience, Ewha Womans University, Seoul 03760, Republic of Korea
| | - Jun Soo Kim
- Department of Chemistry and Nanoscience, Ewha Womans University, Seoul 03760, Republic of Korea
| |
Collapse
|
8
|
Caceres-Delpiano J, Wang LP, Essex JW. The automated optimisation of a coarse-grained force field using free energy data. Phys Chem Chem Phys 2021; 23:24842-24851. [PMID: 34723311 PMCID: PMC8579472 DOI: 10.1039/d0cp05041e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2020] [Accepted: 10/18/2021] [Indexed: 11/21/2022]
Abstract
Atomistic models provide a detailed representation of molecular systems, but are sometimes inadequate for simulations of large systems over long timescales. Coarse-grained models enable accelerated simulations by reducing the number of degrees of freedom, at the cost of reduced accuracy. New optimisation processes to parameterise these models could improve their quality and range of applicability. We present an automated approach for the optimisation of coarse-grained force fields, by reproducing free energy data derived from atomistic molecular simulations. To illustrate the approach, we implemented hydration free energy gradients as a new target for force field optimisation in ForceBalance and applied it successfully to optimise the un-charged side-chains and the protein backbone in the SIRAH protein coarse-grain force field. The optimised parameters closely reproduced hydration free energies of atomistic models and gave improved agreement with experiment.
Collapse
Affiliation(s)
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, California 95616, USA.
| | - Jonathan W Essex
- School of Chemistry, University of Southampton, Southapton, S017 1BJ, UK.
| |
Collapse
|
9
|
Liwo A, Czaplewski C, Sieradzan AK, Lipska AG, Samsonov SA, Murarka RK. Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems. Biomolecules 2021; 11:1347. [PMID: 34572559 PMCID: PMC8465211 DOI: 10.3390/biom11091347] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Revised: 09/03/2021] [Accepted: 09/09/2021] [Indexed: 12/16/2022] Open
Abstract
Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.
Collapse
Affiliation(s)
- Adam Liwo
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Adam K. Sieradzan
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Agnieszka G. Lipska
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Sergey A. Samsonov
- Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland; (C.C.); (A.K.S.); (A.G.L.); (S.A.S.)
| | - Rajesh K. Murarka
- Department of Chemistry, Indian Institute of Science Education and Research Bhopal, Bhopal Bypass Road, Bhopal 462066, MP, India;
| |
Collapse
|
10
|
Multiscale simulations of large complexes in conjunction with cryo-EM analysis. Curr Opin Struct Biol 2021; 72:27-32. [PMID: 34399155 DOI: 10.1016/j.sbi.2021.07.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2021] [Revised: 07/15/2021] [Accepted: 07/16/2021] [Indexed: 11/21/2022]
Abstract
The cellular environment is highly crowded with most proteins and RNA/DNA forming homomeric and heteromeric complexes. Essential questions regarding how these complexes switch between functional, rest, and abnormal states with regulators or modifications remain challenging and complicated. Here, we review the recent progress integrating cryoelectron microscopy and multiscale molecular modeling to understand the dynamics and function-related mechanism in protein-RNA/DNA complexes, protein-protein complexes/assemblies, and membrane protein complexes. One future direction of multiscale simulations will be to interpret the large complex multibody regulation in assembly-induced function enhancement in conjunction with advanced atomic resolution structural-biology techniques and specialized computing architectures.
Collapse
|
11
|
Bae S, Oh I, Yoo J, Kim JS. Effect of DNA Flexibility on Complex Formation of a Cationic Nanoparticle with Double-Stranded DNA. ACS OMEGA 2021; 6:18728-18736. [PMID: 34337212 PMCID: PMC8319935 DOI: 10.1021/acsomega.1c01709] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 07/05/2021] [Indexed: 06/13/2023]
Abstract
We present extensive molecular dynamics simulations of a cationic nanoparticle and a double-stranded DNA molecule to discuss the effect of DNA flexibility on the complex formation of a cationic nanoparticle with double-stranded DNA. Martini coarse-grained models were employed to describe double-stranded DNA molecules with two different flexibilities and cationic nanoparticles with three different electric charges. As the electric charge of a cationic nanoparticle increases, the degree of DNA bending increases, eventually leading to the wrapping of DNA around the nanoparticle at high electric charges. However, a small increase in the persistence length of DNA by 10 nm requires a cationic nanoparticle with a markedly increased electric charge to bend and wrap DNA around. Thus, a more flexible DNA molecule bends and wraps around a cationic nanoparticle with an intermediate electric charge, whereas a less flexible DNA molecule binds to a nanoparticle with the same electric charge without notable bending. This work provides solid evidence that a small difference in DNA flexibility (as small as 10 nm in persistence length) has a substantial influence on the complex formation of DNA with proteins from a biological perspective and suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnology as a new tool to manipulate the structure of DNA molecules mediated by nanoparticle binding.
Collapse
Affiliation(s)
- Sehui Bae
- Department
of Chemistry and Nanoscience, Ewha Womans
University, Seoul 03760, Republic of Korea
| | - Inrok Oh
- LG
Chem Ltd., LG Science Park, Seoul 07796, Republic of Korea
| | - Jejoong Yoo
- Department
of Physics, Sungkyunkwan University, Suwon 16419, Republic of Korea
| | - Jun Soo Kim
- Department
of Chemistry and Nanoscience, Ewha Womans
University, Seoul 03760, Republic of Korea
| |
Collapse
|
12
|
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
|
13
|
Zhang X, Sundram S, Oppelstrup T, Kokkila-Schumacher SIL, Carpenter TS, Ingólfsson HI, Streitz FH, Lightstone FC, Glosli JN. ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field. J Chem Phys 2020; 153:045103. [PMID: 32752727 DOI: 10.1063/5.0014500] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
We have implemented the Martini force field within Lawrence Livermore National Laboratory's molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein-lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns. The ddcMD code has been made open source and is available on GitHub at https://github.com/LLNL/ddcMD.
Collapse
Affiliation(s)
- Xiaohua Zhang
- Physical and Life Sciences (PLS) Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Shiv Sundram
- Computing Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Tomas Oppelstrup
- Physical and Life Sciences (PLS) Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | | | - Timothy S Carpenter
- Physical and Life Sciences (PLS) Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Helgi I Ingólfsson
- Physical and Life Sciences (PLS) Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Frederick H Streitz
- Physical and Life Sciences (PLS) Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Felice C Lightstone
- Physical and Life Sciences (PLS) Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - James N Glosli
- Physical and Life Sciences (PLS) Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|
14
|
Wilson KA, Wang L, MacDermott-Opeskin H, O'Mara ML. The Fats of Life: Using Computational Chemistry to Characterise the Eukaryotic Cell Membrane. Aust J Chem 2020. [DOI: 10.1071/ch19353] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Our current knowledge of the structural dynamics and complexity of lipid bilayers is still developing. Computational techniques, especially molecular dynamics simulations, have increased our understanding significantly as they allow us to model functions that cannot currently be experimentally resolved. Here we review available computational tools and techniques, the role of the major lipid species, insights gained into lipid bilayer structure and function from molecular dynamics simulations, and recent progress towards the computational modelling of the physiological complexity of eukaryotic lipid bilayers.
Collapse
|
15
|
Ricci M, Roscioni OM, Querciagrossa L, Zannoni C. MOLC. A reversible coarse grained approach using anisotropic beads for the modelling of organic functional materials. Phys Chem Chem Phys 2019; 21:26195-26211. [PMID: 31755499 DOI: 10.1039/c9cp04120f] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
We describe the development and implementation of a coarse grained (CG) modelling approach where complex organic molecules, and particularly the π-conjugated ones often employed in organic electronics, are modelled in terms of connected sets of attractive-repulsive biaxial Gay-Berne ellipsoidal beads. The CG model is aimed at reproducing realistically large scale morphologies (e.g. up to 100 nm thick films) for the materials involved, while being able to generate, with a back-mapping procedure, atomistic coordinates suitable, with limited effort, to be applied for charge transport calculations. Detailed methodology and an application to the common hole transporter material α-NPD are provided.
Collapse
Affiliation(s)
- Matteo Ricci
- Dipartimento di Chimica Industriale "Toso Montanari" and INSTM, Università di Bologna, Viale Risorgimento 4, IT-40136 Bologna, Italy.
| | | | | | | |
Collapse
|
16
|
Sun T, Mirzoev A, Minhas V, Korolev N, Lyubartsev AP, Nordenskiöld L. A multiscale analysis of DNA phase separation: from atomistic to mesoscale level. Nucleic Acids Res 2019; 47:5550-5562. [PMID: 31106383 PMCID: PMC6582353 DOI: 10.1093/nar/gkz377] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2019] [Revised: 04/26/2019] [Accepted: 05/09/2019] [Indexed: 12/17/2022] Open
Abstract
DNA condensation and phase separation is of utmost importance for DNA packing in vivo with important applications in medicine, biotechnology and polymer physics. The presence of hexagonally ordered DNA is observed in virus capsids, sperm heads and in dinoflagellates. Rigorous modelling of this process in all-atom MD simulations is presently difficult to achieve due to size and time scale limitations. We used a hierarchical approach for systematic multiscale coarse-grained (CG) simulations of DNA phase separation induced by the three-valent cobalt(III)-hexammine (CoHex3+). Solvent-mediated effective potentials for a CG model of DNA were extracted from all-atom MD simulations. Simulations of several hundred 100-bp-long CG DNA oligonucleotides in the presence of explicit CoHex3+ ions demonstrated aggregation to a liquid crystalline hexagonally ordered phase. Following further coarse-graining and extraction of effective potentials, we conducted modelling at mesoscale level. In agreement with electron microscopy observations, simulations of an 10.2-kb-long DNA molecule showed phase separation to either a toroid or a fibre with distinct hexagonal DNA packing. The mechanism of toroid formation is analysed in detail. The approach used here is based only on the underlying all-atom force field and uses no adjustable parameters and may be generalised to modelling chromatin up to chromosome size.
Collapse
Affiliation(s)
- Tiedong Sun
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Alexander Mirzoev
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Vishal Minhas
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Nikolay Korolev
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Alexander P Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm University, 10691 Stockholm, Sweden
| | - Lars Nordenskiöld
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| |
Collapse
|
17
|
Sanbonmatsu KY. Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR. Curr Opin Struct Biol 2019; 55:104-113. [PMID: 31125796 DOI: 10.1016/j.sbi.2019.03.004] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Accepted: 03/01/2019] [Indexed: 12/11/2022]
Abstract
Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.
Collapse
|
18
|
Pant P, Jayaram B. C5' omitted DNA enhances bendability and protein binding. Biochem Biophys Res Commun 2019; 514:979-984. [PMID: 31092333 DOI: 10.1016/j.bbrc.2019.05.051] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2019] [Accepted: 05/06/2019] [Indexed: 12/24/2022]
Abstract
Protein-DNA interactions are of great biological importance. The specificity and strength of these intimate contacts are crucial in the proper functioning of a cell, wherein the role of DNA dynamic bendability has been a matter of discussion. We relate DNA bendability to protein binding by introducing some simple modifications in the DNA structure. We removed C5' carbon in first modified structure and the second has an additional carbon between C3' and 3'-OH, hereby pronounced as C(-) and C(+) nucleic acids respectively. We observed that C(+) nucleic acid retains B-DNA duplex as seen by means of 500 ns long molecular dynamics (MD) simulations, structural and energetic calculations, while C(-) nucleic acid attains a highly bend structure. We transferred these observations to a protein-DNA system in order to monitor as to what extent the bendability enhances the protein binding. The energetics of binding is explored by performing 100 ns long MD simulations on control and modified DNA-protein complexes followed by running MM-PBSA/GBSA calculations on the resultant structures. It is observed that C(+) nucleic acid has protein binding in close correspondence to the control system (∼-14 kcal/mol) due to their relatable structure, while the C(-) nucleic acid displayed high binding to the protein (∼-18 kcal/mol). DelPhi based calculations reveal that the high binding could be the result of enhanced electrostatic interactions caused by exposed bases in the bend structure for protein recognition. Such modified oligonucleotides, due to their improved binding to protein and resistance to nuclease degradation, have a great therapeutic value.
Collapse
Affiliation(s)
- Pradeep Pant
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, 110016, India; Supercomputing Facility for Bioinformatics & Computational Biology, Hauz Khas, New Delhi, 110016, India
| | - B Jayaram
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, 110016, India; Supercomputing Facility for Bioinformatics & Computational Biology, Hauz Khas, New Delhi, 110016, India; Kusuma School of Biological Sciences, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, 110016, India.
| |
Collapse
|
19
|
Marrink SJ, Corradi V, Souza PC, Ingólfsson HI, Tieleman DP, Sansom MS. Computational Modeling of Realistic Cell Membranes. Chem Rev 2019; 119:6184-6226. [PMID: 30623647 PMCID: PMC6509646 DOI: 10.1021/acs.chemrev.8b00460] [Citation(s) in RCA: 410] [Impact Index Per Article: 82.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2018] [Indexed: 12/15/2022]
Abstract
Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.
Collapse
Affiliation(s)
- Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Paulo C.T. Souza
- Groningen
Biomolecular Sciences and Biotechnology Institute & Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Helgi I. Ingólfsson
- Biosciences
and Biotechnology Division, Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94550, United States
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Mark S.P. Sansom
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, U.K.
| |
Collapse
|
20
|
Machado MR, Barrera EE, Klein F, Sóñora M, Silva S, Pantano S. The SIRAH 2.0 Force Field: Altius, Fortius, Citius. J Chem Theory Comput 2019; 15:2719-2733. [PMID: 30810317 DOI: 10.1021/acs.jctc.9b00006] [Citation(s) in RCA: 86] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
A new version of the coarse-grained (CG) SIRAH force field for proteins has been developed. Modifications to bonded and non-bonded interactions on the existing molecular topologies significantly ameliorate the structural description and flexibility of a non-redundant set of proteins. The SIRAH 2.0 force field has also been ported to the popular simulation package AMBER, which along with the former implementation in GROMACS expands significantly the potential range of users and performance of this CG force field on CPU/GPU codes. As a non-trivial example of its application, we undertook the structural and dynamical analysis of the most abundant and conserved calcium-binding protein, calmodulin (CaM). CaM is composed of two calcium-binding motifs called EF-hands, which in the presence of calcium specifically recognize a cognate peptide by embracing it. CG simulations of CaM bound to four calcium ions in the presence or absence of a binding peptide (holo and apo forms, respectively) resulted in good and stable ion coordination. The simulation of the holo form starting from an experimental structure sampled near-native conformations, retrieving quasi-atomistic precision. Removing the binding peptide enabled the EF-hands to perform large reciprocal movements, comparable to those observed in NMR structures. On the other hand, the isolated peptide starting from the helical conformation experienced spontaneous unfolding, in agreement with previous experimental data. However, repositioning the peptide in the neighborhood of one EF-hand not only prevented the peptide from unfolding but also drove CaM to a fully bound conformation, with both EF-hands embracing the cognate peptide, resembling the experimental holo structure. Therefore, SIRAH 2.0 shows the capacity to handle a number of structurally and dynamically challenging situations, including metal ion coordination, unbiased conformational sampling, and specific protein-peptide recognition.
Collapse
Affiliation(s)
- Matías R Machado
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Exequiel E Barrera
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Florencia Klein
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Martín Sóñora
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Steffano Silva
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| | - Sergio Pantano
- Biomolecular Simulations Group , Institut Pasteur de Montevideo , Mataojo 2020 , CP 11400 Montevideo , Uruguay
| |
Collapse
|
21
|
Sulpizi M, Faller R, Pantano S. Multiscale modeling on biological systems. Biochem Biophys Res Commun 2018; 498:263. [PMID: 29496451 DOI: 10.1016/j.bbrc.2018.02.179] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Marialore Sulpizi
- Department of Physics, Johannes Gutenberg University, Mainz, Germany
| | - Roland Faller
- Department of Chemical Engineering, University of California-Davis, Davis, CA, USA
| | | |
Collapse
|