1
|
Pirojsirikul T, Lee VS, Nimmanpipug P. Unraveling Bacterial Single-Stranded Sequence Specificities: Insights from Molecular Dynamics and MMPBSA Analysis of Oligonucleotide Probes. Mol Biotechnol 2024; 66:582-591. [PMID: 38374320 DOI: 10.1007/s12033-024-01082-0] [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: 10/03/2022] [Accepted: 01/10/2024] [Indexed: 02/21/2024]
Abstract
We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.
Collapse
Affiliation(s)
- Teerapong Pirojsirikul
- Division of Physical Science, Faculty of Science, Prince of Songkla University, Songkhla, 90110, Thailand.
| | - Vannajan Sanghiran Lee
- Department of Chemistry, Center of Theoretical and Computational Physics, Faculty of Science, University of Malaya, 50603, Kuala Lumpur, Malaysia
| | - Piyarat Nimmanpipug
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, 50200, Thailand
| |
Collapse
|
2
|
Ishida H, Kono H. Free Energy Landscape of H2A-H2B Displacement From Nucleosome. J Mol Biol 2022; 434:167707. [PMID: 35777463 DOI: 10.1016/j.jmb.2022.167707] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 06/11/2022] [Accepted: 06/23/2022] [Indexed: 12/14/2022]
Abstract
Nucleosome reconstitution plays an important role in many cellular functions. As an initial step, H2A-H2B dimer displacement, which is accompanied by disruption of many of the interactions within the nucleosome, should occur. To understand how H2A-H2B dimer displacement occurs, an adaptively biased molecular dynamics (ABMD) simulation was carried out to generate a variety of displacements of the H2A-H2B dimer from the fully wrapped to partially unwrapped nucleosome structures. With regards to these structures, the free energy landscape of the dimer displacement was investigated using umbrella sampling simulations. We found that the main contributors to the free energy were the docking domain of H2A and the C-terminal of H4. There were various paths for the dimer displacement which were dependent on the extent of nucleosomal DNA wrapping, suggesting that modulation of the intra-nucleosomal interaction by external factors such as histone chaperons could control the path for the H2A-H2B dimer displacement. Key residues which contributed to the free energy have also been reported to be involved in the mutations and posttranslational modifications (PTMs) which are important for assembling and/or reassembling the nucleosome at the molecular level and are found in cancer cells at the phenotypic level. Our results give insight into how the H2A-H2B dimer displacement proceeds along various paths according to different interactions within the nucleosome.
Collapse
Affiliation(s)
- Hisashi Ishida
- Institute for Quantum Life Science, National Institutes for Quantum Science and Technology, 619-0215 Kizugawa, Kyoto, Japan.
| | - Hidetoshi Kono
- Institute for Quantum Life Science, National Institutes for Quantum Science and Technology, 619-0215 Kizugawa, Kyoto, Japan
| |
Collapse
|
3
|
Abstract
Mismatched base pairs alter the flexibility and intrinsic curvature of DNA. The role of such DNA features is not fully understood in the mismatch repair pathway. MutS/DNA complexes exhibit DNA bending, PHE intercalation, and changes of base-pair parameters near the mismatch. Recently, we have shown that base-pair opening in the absence of MutS can discriminate mismatches from canonical base pairs better than DNA bending. However, DNA bending in the absence of MutS was found to be rather challenging to describe correctly. Here, we present a computational study on the DNA bending of canonical and G/T mismatched DNAs. Five types of geometric parameters covering template-based bending toward the experimental DNA structure, global, and local geometry parameters were employed in biased molecular dynamics in the absence of MutS. None of these parameters showed higher discrimination than the base-pair opening. Only roll could induce a sharply localized bending of DNA as observed in the experimental MutS/DNA structure. Further, we demonstrated that the intercalation of benzene mimicking PHE decreases the energetic cost of DNA bending without any effect on mismatch discrimination.
Collapse
Affiliation(s)
- Tomáš Bouchal
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,CEITEC─Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Ivo Durník
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,CEITEC─Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Petr Kulhánek
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,CEITEC─Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
4
|
Torsional stress can regulate the unwrapping of two outer half superhelical turns of nucleosomal DNA. Proc Natl Acad Sci U S A 2021; 118:2020452118. [PMID: 33558240 DOI: 10.1073/pnas.2020452118] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Torsional stress has a significant impact on the structure and stability of the nucleosome. RNA polymerase imposes torsional stress on the DNA in chromatin and unwraps the DNA from the nucleosome to access the genetic information encoded in the DNA. To understand how the torsional stress affects the stability of the nucleosome, we examined the unwrapping of two half superhelical turns of nucleosomal DNA from either end of the DNA under torsional stress with all-atom molecular dynamics simulations. The free energies for unwrapping the DNA indicate that positive stress that overtwists DNA facilitates a large-scale asymmetric unwrapping of the DNA without a large extension of the DNA. During the unwrapping, one end of the DNA was dissociated from H3 and H2A-H2B, while the other end of the DNA stably remained wrapped. The detailed analysis indicates that this asymmetric dissociation is facilitated by the geometry and bendability of the DNA under positive stress. The geometry stabilized the interaction between the major groove of the twisted DNA and the H3 αN-helix, and the straightened DNA destabilized the interaction with H2A-H2B. Under negative stress, the DNA became more bendable and flexible, which facilitated the binding of the unwrapped DNA to the octamer in a stable state. Consequently, we conclude that the torsional stress has a significant impact on the affinity of the DNA and the octamer through the inherent nature of the DNA and can change the accessibility of regulatory proteins.
Collapse
|
5
|
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
|
6
|
Wang E, Sun H, Wang J, Wang Z, Liu H, Zhang JZH, Hou T. End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design. Chem Rev 2019; 119:9478-9508. [DOI: 10.1021/acs.chemrev.9b00055] [Citation(s) in RCA: 578] [Impact Index Per Article: 115.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Ercheng Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Huiyong Sun
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Junmei Wang
- Department of Pharmaceutical Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Zhe Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Hui Liu
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - John Z. H. Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU−ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai 200122, China
- Department of Chemistry, New York University, New York, New York 10003, United States
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Tingjun Hou
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| |
Collapse
|
7
|
Maffeo C, Chou HY, Aksimentiev A. Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations. ADVANCED THEORY AND SIMULATIONS 2019; 2:1800191. [PMID: 31728433 PMCID: PMC6855400 DOI: 10.1002/adts.201800191] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2018] [Indexed: 12/15/2022]
Abstract
Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.
Collapse
Affiliation(s)
- Christopher Maffeo
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| | - Han-Yi Chou
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| | - Aleksei Aksimentiev
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| |
Collapse
|
8
|
Kono H, Sakuraba S, Ishida H. Free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA. PLoS Comput Biol 2018; 14:e1006024. [PMID: 29505570 PMCID: PMC5854429 DOI: 10.1371/journal.pcbi.1006024] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2017] [Revised: 03/15/2018] [Accepted: 02/01/2018] [Indexed: 11/19/2022] Open
Abstract
The eukaryotic genome is packaged into a nucleus in the form of chromatin. The fundamental structural unit of chromatin is a protein-DNA complex, the nucleosome, where 146 or 147 base pairs of DNA wrap 1.75 times around a histone core. To function in cellular processes, however, nucleosomal DNA must be unwrapped. Although this unwrapping has been experimentally investigated, details of the process at an atomic level are not yet well understood. Here, we used molecular dynamics simulation with an enhanced sampling method to calculate the free energy profiles for unwrapping the outer superhelical turn of nucleosomal DNA. A free energy change of about 11.5 kcal/mol for the unwrapping agrees well with values obtained in single molecule experiments. This simulation revealed a variety of conformational states, indicating there are many potential paths to outer superhelicdal turn unwrapping, but the dominant path is likely asymmetric. At one end of the DNA, the first five bps unwrap, after which a second five bps unwrap at the same end with no increase in free energy. The unwrapping then starts at the other end of the DNA, where 10 bps are unwrapped. During further unwrapping of 15 bps, the unwrapping advances at one of the ends, after which the other end of the DNA unwraps to complete the unwrapping of the outer superhelical turn. These results provide insight into the construction, disruption, and repositioning of nucleosomes, which are continuously ongoing during cellular processes.
Collapse
Affiliation(s)
- Hidetoshi Kono
- Molecular Modeling and Simulation Group, Department of Quantum Beam Life Science, National Institutes for Quantum and Radiological Science and Technology, Umemidai, Kizugawa, Kyoto, Japan
- * E-mail:
| | - Shun Sakuraba
- Graduate School of Frontier Sciences, The University of Tokyo, Kashiwanoha, Kashiwa, Chiba, Japan
| | - Hisashi Ishida
- Molecular Modeling and Simulation Group, Department of Quantum Beam Life Science, National Institutes for Quantum and Radiological Science and Technology, Umemidai, Kizugawa, Kyoto, Japan
| |
Collapse
|
9
|
Wang C, Greene D, Xiao L, Qi R, Luo R. Recent Developments and Applications of the MMPBSA Method. Front Mol Biosci 2018; 4:87. [PMID: 29367919 PMCID: PMC5768160 DOI: 10.3389/fmolb.2017.00087] [Citation(s) in RCA: 332] [Impact Index Per Article: 55.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2017] [Accepted: 11/30/2017] [Indexed: 12/23/2022] Open
Abstract
The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.
Collapse
Affiliation(s)
- Changhao Wang
- Chemical and Materials Physics Graduate Program, University of California, Irvine, Irvine, CA, United States
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Physics and Astronomy, University of California, Irvine, Irvine, CA, United States
| | - D'Artagnan Greene
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
| | - Li Xiao
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, United States
| | - Ruxi Qi
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
| | - Ray Luo
- Chemical and Materials Physics Graduate Program, University of California, Irvine, Irvine, CA, United States
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, United States
- Department of Chemical Engineering and Materials Science, University of California, Irvine, Irvine, CA, United States
| |
Collapse
|
10
|
Ishida H, Kono H. H4 Tails Potentially Produce the Diversity in the Orientation of Two Nucleosomes. Biophys J 2017; 113:978-990. [PMID: 28877499 DOI: 10.1016/j.bpj.2017.07.015] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2017] [Revised: 06/29/2017] [Accepted: 07/20/2017] [Indexed: 11/17/2022] Open
Abstract
Histone tails play an important role in internucleosomal interaction and chromatin compaction. To understand how the H4 tails are involved in the internucleosomal interaction, an adaptively biased molecular dynamics simulation of 63 models of two stacked nucleosomes, each with the H4 tails in different locations, was carried out. This simulation generated a variety of orientations of the separated nucleosomes depending on the formation of the H4 tail bridge between the H4 tails and the DNA of the neighboring nucleosomes. For the models that showed distinctive orientations of the two nucleosomes, the free energies of the separation of the nucleosomes were further investigated using umbrella sampling simulations. The attractive force between the nucleosomes was estimated from the free energies; the force when two H4 tail bridges formed varied from 36 to 63 pN, depending on the formation of the H4 tail-bridge and the interfacial interaction, whereas the force reduced to 15-18 pN after either one of the H4 tail bridges had broken, regardless of the conformation of the H4 tail. Additional simulations of the nucleosomes show that when the H4 tail was truncated, the force between the nucleosomes became repulsive (from-3 to -7 pN). We concluded that the H4 tails potentially produce the diversity in the orientation of the two nucleosomes, which would contribute to the polymorphism of the chromatin structure.
Collapse
Affiliation(s)
- Hisashi Ishida
- Molecular Modeling and Simulation Group, Department of Quantum Beam Life Science, National Institutes for Quantum and Radiological Science and Technology, Kizugawa, Kyoto, Japan.
| | - Hidetoshi Kono
- Molecular Modeling and Simulation Group, Department of Quantum Beam Life Science, National Institutes for Quantum and Radiological Science and Technology, Kizugawa, Kyoto, Japan
| |
Collapse
|
11
|
Lakhani B, Thayer KM, Hingorani MM, Beveridge DL. Evolutionary Covariance Combined with Molecular Dynamics Predicts a Framework for Allostery in the MutS DNA Mismatch Repair Protein. J Phys Chem B 2017; 121:2049-2061. [PMID: 28135092 PMCID: PMC5346969 DOI: 10.1021/acs.jpcb.6b11976] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
![]()
Mismatch
repair (MMR) is an essential, evolutionarily conserved
pathway that maintains genome stability by correcting base-pairing
errors in DNA. Here we examine the sequence and structure of MutS
MMR protein to decipher the amino acid framework underlying its two
key activities—recognizing mismatches in DNA and using ATP
to initiate repair. Statistical coupling analysis (SCA) identified
a network (sector) of coevolved amino acids in the MutS protein family.
The potential functional significance of this SCA sector was assessed
by performing molecular dynamics (MD) simulations for alanine mutants
of the top 5% of 160 residues in the distribution, and control nonsector
residues. The effects on three independent metrics were monitored:
(i) MutS domain conformational dynamics, (ii) hydrogen bonding between
MutS and DNA/ATP, and (iii) relative ATP binding free energy. Each
measure revealed that sector residues contribute more substantively
to MutS structure–function than nonsector residues. Notably,
sector mutations disrupted MutS contacts with DNA and/or ATP from
a distance via contiguous pathways and correlated motions, supporting
the idea that SCA can identify amino acid networks underlying allosteric
communication. The combined SCA/MD approach yielded novel, experimentally
testable hypotheses for unknown roles of many residues distributed
across MutS, including some implicated in Lynch cancer syndrome.
Collapse
Affiliation(s)
- Bharat Lakhani
- Molecular Biology and Biochemistry Department, ‡Molecular Biophysics Program, §Chemistry Department, and ∥Computer Science Department, Wesleyan University , Middletown, Connecticut 06459, United States
| | - Kelly M Thayer
- Molecular Biology and Biochemistry Department, ‡Molecular Biophysics Program, §Chemistry Department, and ∥Computer Science Department, Wesleyan University , Middletown, Connecticut 06459, United States
| | - Manju M Hingorani
- Molecular Biology and Biochemistry Department, ‡Molecular Biophysics Program, §Chemistry Department, and ∥Computer Science Department, Wesleyan University , Middletown, Connecticut 06459, United States
| | - David L Beveridge
- Molecular Biology and Biochemistry Department, ‡Molecular Biophysics Program, §Chemistry Department, and ∥Computer Science Department, Wesleyan University , Middletown, Connecticut 06459, United States
| |
Collapse
|