1
|
Lu W, Onuchic JN, Di Pierro M. An associative memory Hamiltonian model for DNA and nucleosomes. PLoS Comput Biol 2023; 19:e1011013. [PMID: 36972316 PMCID: PMC10079229 DOI: 10.1371/journal.pcbi.1011013] [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: 10/15/2022] [Revised: 04/06/2023] [Accepted: 03/08/2023] [Indexed: 03/29/2023] Open
Abstract
A model for DNA and nucleosomes is introduced with the goal of studying chromosomes from a single base level all the way to higher-order chromatin structures. This model, dubbed the Widely Editable Chromatin Model (WEChroM), reproduces the complex mechanics of the double helix including its bending persistence length and twisting persistence length, and their respective temperature dependence. The WEChroM Hamiltonian is composed of chain connectivity, steric interactions, and associative memory terms representing all remaining interactions leading to the structure, dynamics, and mechanical characteristics of the B-DNA. Several applications of this model are discussed to demonstrate its applicability. WEChroM is used to investigate the behavior of circular DNA in the presence of positive and negative supercoiling. We show that it recapitulates the formation of plectonemes and of structural defects that relax mechanical stress. The model spontaneously manifests an asymmetric behavior with respect to positive or negative supercoiling, similar to what was previously observed in experiments. Additionally, we show that the associative memory Hamiltonian is also capable of reproducing the free energy of partial DNA unwrapping from nucleosomes. WEChroM is designed to emulate the continuously variable mechanical properties of the 10nm fiber and, by virtue of its simplicity, is ready to be scaled up to molecular systems large enough to investigate the structural ensembles of genes. WEChroM is implemented in the OpenMM simulation toolkits and is freely available for public use.
Collapse
Affiliation(s)
- Weiqi Lu
- Center for Theoretical Biological Physics, & Department of Physics and Astronomy, Rice University, Houston, Texas, United States of America
| | - José N. Onuchic
- Center for Theoretical Biological Physics, & Department of Physics and Astronomy, Rice University, Houston, Texas, United States of America
- Department of Chemistry, & Department of Biosciences, Rice University, Houston, Texas, United States of America
- * E-mail: (JNO); (MDP)
| | - Michele Di Pierro
- Department of Physics, Northeastern University, Boston, Massachusetts, United States of America
- Center for Theoretical Biological Physics, Northeastern University, Boston, Massachusetts, United States of America
- * E-mail: (JNO); (MDP)
| |
Collapse
|
2
|
Sun LZ, Qian JL, Cai P, Xu X. Mutual effects between single-stranded DNA conformation and Na +-Mg 2+ ion competition in mixed salt solutions. Phys Chem Chem Phys 2022; 24:20867-20881. [PMID: 36043348 DOI: 10.1039/d2cp02737b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The ion-dependence of single-stranded DNA (ssDNA) conformational changes has attracted growing attention because of its biological and technological importance. Although single-species ion effects have been extensively explored, it is challenging to study the ssDNA conformational properties under mixed monovalent/divalent ion conditions due to the complications of ssDNA flexibility and ion-ion competition. In this study, we apply Langevin dynamics simulations to investigate mixed Na+/Mg2+ ion-dependent ssDNA conformations. The ssDNA structure is described using a coarse-grained model, in which the phosphate, base, and sugar of each nucleotide are represented by three different beads. A novel improvement in our simulation model is that mixed-salt-related electrostatic interactions are computed via combining Manning counterion condensation (MCC) theory with the Monte Carlo tightly bound ion (MCTBI) model. Based on this MCC-MCTBI combination, we report new empirical functions to describe the ion-concentration-dependent and ssDNA conformation/structure-dependent electrostatic effects. The calculation results relating to the ion binding properties and the simulation results relating to the ssDNA conformational properties are validated against experimental results. In addition, our simulation results suggest a quantitative relationship between the ssDNA conformation and Na+-Mg2+ competition; this in turn reveals their mutual impact in the ion atmosphere.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.
| | - Jun-Lin Qian
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.
| | - Pinggen Cai
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.
| | - Xiaojun Xu
- Institute of Bioinformatics and Medical Engineering, Jiangsu University of Technology, Changzhou, 213001, China
| |
Collapse
|
3
|
Minhas V, Sun T, Mirzoev A, Korolev N, Lyubartsev AP, Nordenskiöld L. Modeling DNA Flexibility: Comparison of Force Fields from Atomistic to Multiscale Levels. J Phys Chem B 2019; 124:38-49. [DOI: 10.1021/acs.jpcb.9b09106] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Vishal Minhas
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Tiedong Sun
- School of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Alexander Mirzoev
- 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
|
4
|
Harrison RM, Romano F, Ouldridge TE, Louis AA, Doye JPK. Identifying Physical Causes of Apparent Enhanced Cyclization of Short DNA Molecules with a Coarse-Grained Model. J Chem Theory Comput 2019; 15:4660-4672. [PMID: 31282669 PMCID: PMC6694408 DOI: 10.1021/acs.jctc.9b00112] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
![]()
DNA
cyclization is a powerful technique to gain insight into the nature
of DNA bending. While the wormlike chain model provides a good description
of small to moderate bending fluctuations, it is expected to break
down for large bending. Recent cyclization experiments on strongly
bent shorter molecules indeed suggest enhanced flexibility over and
above that expected from the wormlike chain. Here, we use a coarse-grained
model of DNA to investigate the subtle thermodynamics of DNA cyclization
for molecules ranging from 30 to 210 base pairs. As the molecules
get shorter, we find increasing deviations between our computed equilibrium j-factor and the classic wormlike chain predictions of Shimada
and Yamakawa for a torsionally aligned looped molecule. These deviations
are due to sharp kinking, first at nicks, and only subsequently in
the body of the duplex. At the shortest lengths, substantial fraying
at the ends of duplex domains is the dominant method of relaxation.
We also estimate the dynamic j-factor measured in
recent FRET experiments. We find that the dynamic j-factor is systematically larger than its equilibrium counterpart—with
the deviation larger for shorter molecules—because not all
the stress present in the fully cyclized state is present in the transition
state. These observations are important for the interpretation of
recent cyclization experiments, suggesting that measured anomalously
high j-factors may not necessarily indicate non-WLC
behavior in the body of duplexes.
Collapse
Affiliation(s)
- Ryan M Harrison
- Physical & Theoretical Chemistry Laboratory, Department of Chemistry , University of Oxford , South Parks Road , Oxford OX1 3QZ , United Kingdom
| | - Flavio Romano
- Dipartimento di Scienze Molecolari e Nanosistemi , Universitá Ca' Foscari Venezia , I-30123 Venezia , Italy
| | - Thomas E Ouldridge
- Imperial College Centre for Synthetic Biology and Department of Bioengineering , Imperial College London , 180 Queen's Road , London SW7 2AZ , United Kingdom
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, Department of Physics , University of Oxford , 1 Keble Road , Oxford OX1 3NP , United Kingdom
| | - Jonathan P K Doye
- Physical & Theoretical Chemistry Laboratory, Department of Chemistry , University of Oxford , South Parks Road , Oxford OX1 3QZ , United Kingdom
| |
Collapse
|
5
|
Shimizu M, Takada S. Reconstruction of Atomistic Structures from Coarse-Grained Models for Protein-DNA Complexes. J Chem Theory Comput 2018; 14:1682-1694. [PMID: 29397721 DOI: 10.1021/acs.jctc.7b00954] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
While coarse-grained (CG) simulations have widely been used to accelerate structure sampling of large biomolecular complexes, they are unavoidably less accurate and thus the reconstruction of all-atom (AA) structures and the subsequent refinement is desirable. In this study we developed an efficient method to reconstruct AA structures from sampled CG protein-DNA complex models, which attempts to model the protein-DNA interface accurately. First we developed a method to reconstruct atomic details of DNA structures from a three-site per nucleotide CG model, which uses a DNA fragment library. Next, for the protein-DNA interface, we referred to the side chain orientations in the known structure of the target interface when available. The other parts are modeled by existing tools. We confirmed the accuracy of the protocol in various aspects including the structure deviation in the self-reproduction, the base pair reproducibility, atomic contacts at the protein-DNA interface, and feasibility of the posterior AA simulations.
Collapse
Affiliation(s)
- Masahiro Shimizu
- Department of Biophysics, Graduate School of Science , Kyoto University , Sakyo, Kyoto 606-8502 Japan
| | - Shoji Takada
- Department of Biophysics, Graduate School of Science , Kyoto University , Sakyo, Kyoto 606-8502 Japan
| |
Collapse
|
6
|
Yagyu H, Lee JY, Kim DN, Tabata O. Coarse-Grained Molecular Dynamics Model of Double-Stranded DNA for DNA Nanostructure Design. J Phys Chem B 2017; 121:5033-5039. [DOI: 10.1021/acs.jpcb.7b03931] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hiromasa Yagyu
- Department
of Mechanical Engineering, Kanto Gakuin University, 1-50-1 Mutsuura-higashi, Kanazawa-ku, Yokohama 236-8501, Japan
| | - Jae-Young Lee
- Department
of Mechanical and Aerospace Engineering, Seoul National University, 1 Gwanak-ro, Daehak-dong, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Do-Nyun Kim
- Department
of Mechanical and Aerospace Engineering, Seoul National University, 1 Gwanak-ro, Daehak-dong, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Osamu Tabata
- Department
of Micro Engineering, Kyoto University, Kyoto Daigaku-Katsura C3, Nishikyo-ku, Kyoto, 606-8501, Japan
| |
Collapse
|
7
|
Lyubartsev AP, Naômé A, Vercauteren DP, Laaksonen A. Systematic hierarchical coarse-graining with the inverse Monte Carlo method. J Chem Phys 2016; 143:243120. [PMID: 26723605 DOI: 10.1063/1.4934095] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
We outline our coarse-graining strategy for linking micro- and mesoscales of soft matter and biological systems. The method is based on effective pairwise interaction potentials obtained in detailed ab initio or classical atomistic Molecular Dynamics (MD) simulations, which can be used in simulations at less accurate level after scaling up the size. The effective potentials are obtained by applying the inverse Monte Carlo (IMC) method [A. P. Lyubartsev and A. Laaksonen, Phys. Rev. E 52(4), 3730-3737 (1995)] on a chosen subset of degrees of freedom described in terms of radial distribution functions. An in-house software package MagiC is developed to obtain the effective potentials for arbitrary molecular systems. In this work we compute effective potentials to model DNA-protein interactions (bacterial LiaR regulator bound to a 26 base pairs DNA fragment) at physiological salt concentration at a coarse-grained (CG) level. Normally the IMC CG pair-potentials are used directly as look-up tables but here we have fitted them to five Gaussians and a repulsive wall. Results show stable association between DNA and the model protein as well as similar position fluctuation profile.
Collapse
Affiliation(s)
- Alexander P Lyubartsev
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S 106 91 Stockholm, Sweden
| | - Aymeric Naômé
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S 106 91 Stockholm, Sweden
| | | | - Aatto Laaksonen
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, S 106 91 Stockholm, Sweden
| |
Collapse
|
8
|
Dai X, Yin Q, Wan G, Wang R, Shi X, Qiao Y. Effects of Concentrations on the Transdermal Permeation Enhancing Mechanisms of Borneol: A Coarse-Grained Molecular Dynamics Simulation on Mixed-Bilayer Membranes. Int J Mol Sci 2016; 17:E1349. [PMID: 27548141 PMCID: PMC5000745 DOI: 10.3390/ijms17081349] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2016] [Revised: 07/22/2016] [Accepted: 07/27/2016] [Indexed: 01/13/2023] Open
Abstract
Borneol is a natural permeation enhancer that is effective in drugs used in traditional clinical practices as well as in modern scientific research. However, its molecular mechanism is not fully understood. In this study, a mixed coarse-grained model of stratum corneum (SC) lipid bilayer comprised of Ceramide-N-sphingosine (CER NS) 24:0, cholesterol (CHOL) and free fatty acids (FFA) 24:0 (2:2:1) was used to examine the permeation enhancing mechanism of borneol on the model drug osthole. We found two different mechanisms that were dependent on concentrations levels of borneol. At low concentrations, the lipid system maintained a bilayer structure. The addition of borneol made the lipid bilayer loosen and improved drug permeation. The "pull" effect of borneol also improved drug permeation. However, for a strongly hydrophobic drug like osthole, the permeation enhancement of borneol was limited. When most borneol molecules permeated into bilayers and were located at the hydrophobic tail region, the spatial competition effect inhibited drug molecules from permeating deeper into the bilayer. At high concentrations, borneol led to the formation of water pores and long-lived reversed micelles. This improved the permeation of osthole and possibly other hydrophobic or hydrophilic drugs through the SC. Our simulation results were supported by Franz diffusion tests and transmission electron microscope (TEM) experiments.
Collapse
Affiliation(s)
- Xingxing Dai
- School of Traditional Chinese Medicine, Beijing University of Chinese Medicine, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
- Key Laboratory of TCM-information Engineer of State Administration of TCM, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
- Beijing Key Laboratory of Manufacturing Process Control and Quality Evaluation of Chinese Medicine, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
| | - Qianqian Yin
- School of Traditional Chinese Medicine, Beijing University of Chinese Medicine, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
- Key Laboratory of TCM-information Engineer of State Administration of TCM, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
- Beijing Key Laboratory of Manufacturing Process Control and Quality Evaluation of Chinese Medicine, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
| | - Guang Wan
- School of Traditional Chinese Medicine, Capital Medical University, No. 10 of Xitoutiao Outside Youanmen, Fengtai District, Beijing 100069, China.
| | - Ran Wang
- School of Traditional Chinese Medicine, Capital Medical University, No. 10 of Xitoutiao Outside Youanmen, Fengtai District, Beijing 100069, China.
| | - Xinyuan Shi
- School of Traditional Chinese Medicine, Beijing University of Chinese Medicine, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
- Key Laboratory of TCM-information Engineer of State Administration of TCM, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
- Beijing Key Laboratory of Manufacturing Process Control and Quality Evaluation of Chinese Medicine, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
| | - Yanjiang Qiao
- School of Traditional Chinese Medicine, Beijing University of Chinese Medicine, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
- Key Laboratory of TCM-information Engineer of State Administration of TCM, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
- Beijing Key Laboratory of Manufacturing Process Control and Quality Evaluation of Chinese Medicine, No. 6 of Zhonghuan South Road, Wangjing, Chaoyang District, Beijing 100102, China.
| |
Collapse
|
9
|
Dans PD, Walther J, Gómez H, Orozco M. Multiscale simulation of DNA. Curr Opin Struct Biol 2016; 37:29-45. [DOI: 10.1016/j.sbi.2015.11.011] [Citation(s) in RCA: 99] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2015] [Revised: 11/23/2015] [Accepted: 11/25/2015] [Indexed: 01/05/2023]
|