1
|
Strelnikov IA, Kovaleva NA, Klinov AP, Zubova EA. C-B-A Test of DNA Force Fields. ACS OMEGA 2023; 8:10253-10265. [PMID: 36969447 PMCID: PMC10034787 DOI: 10.1021/acsomega.2c07781] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 02/24/2023] [Indexed: 06/18/2023]
Abstract
The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.
Collapse
|
2
|
Nasou AG, Pantatosaki E, Papadopoulos GK. A Simulation Study of the Effect of Naturally Occurring Point Mutations on the SRY-DNA Complex. J Phys Chem B 2022; 126:8921-8930. [PMID: 36315187 DOI: 10.1021/acs.jpcb.2c04852] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Molecular dynamics (MD) simulations were conducted in order to investigate the effect of the naturally occurring point mutations of the transcription factor (TF) sex-determining region Y (SRY) on the structure and dynamics of the SRY-DNA complex. The normal SRY, along with the two mutants I13T and G40R, comprising point mutations on the SRY chain, which have been clinically identified in patients with sex developmental disorders, were modeled as DNA complexes. Our modeling work aims at elucidating atomic-level structural determinants of the aberrant SRY-DNA complexation by means of μs-long MD. The results suggest that the observed disorders brought about by the G40R-DNA and I13T-DNA may arise predominantly from the destabilization of the complex being in accord with in vitro assays found elsewhere and from modifications of the DNA bending as revealed in this study. Comparative potential of mean force computations, over a sequence of short separation distances for the three complexes, verified a higher stability of the normal SRY-DNA. Examining the way the SRY mutations modulate the SRY-DNA complex dynamics at the microscopic level is important also toward elucidating molecular determinants of function for proteins capable of binding to DNA.
Collapse
Affiliation(s)
- Angeliki-Georgia Nasou
- School of Chemical Engineering, National Technical University of Athens, 157 80 Athens, Greece
| | - Evangelia Pantatosaki
- School of Chemical Engineering, National Technical University of Athens, 157 80 Athens, Greece
| | - George K Papadopoulos
- School of Chemical Engineering, National Technical University of Athens, 157 80 Athens, Greece
| |
Collapse
|
3
|
Bosire R, Fadel L, Mocsár G, Nánási P, Sen P, Sharma AK, Naseem MU, Kovács A, Kugel J, Kroemer G, Vámosi G, Szabó G. Doxorubicin impacts chromatin binding of HMGB1, Histone H1 and retinoic acid receptor. Sci Rep 2022; 12:8087. [PMID: 35577872 PMCID: PMC9110345 DOI: 10.1038/s41598-022-11994-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2021] [Accepted: 05/03/2022] [Indexed: 11/10/2022] Open
Abstract
Doxorubicin (Dox), a widely used anticancer DNA-binding drug, affects chromatin in multiple ways, and these effects contribute to both its efficacy and its dose-limiting side effects, especially cardiotoxicity. Here, we studied the effects of Dox on the chromatin binding of the architectural proteins high mobility group B1 (HMGB1) and the linker histone H1, and the transcription factor retinoic acid receptor (RARα) by fluorescence recovery after photobleaching (FRAP) and fluorescence correlation spectroscopy (FCS) in live cells. At lower doses, Dox increased the binding of HMGB1 to DNA while decreasing the binding of the linker histone H1. At higher doses that correspond to the peak plasma concentrations achieved during chemotherapy, Dox reduced the binding of HMGB1 as well. This biphasic effect is interpreted in terms of a hierarchy of competition between the ligands involved and Dox-induced local conformational changes of nucleosome-free DNA. Combined, FRAP and FCS mobility data suggest that Dox decreases the overall binding of RARα to DNA, an effect that was only partially overcome by agonist binding. The intertwined interactions described are likely to contribute to both the effects and side effects of Dox.
Collapse
Affiliation(s)
- Rosevalentine Bosire
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.,Doctoral School of Molecular Cell and Immune Biology, University of Debrecen, Debrecen, Hungary
| | - Lina Fadel
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.,Doctoral School of Molecular Medicine, University of Debrecen, Debrecen, Hungary
| | - Gábor Mocsár
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Péter Nánási
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Pialy Sen
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.,Doctoral School of Molecular Medicine, University of Debrecen, Debrecen, Hungary
| | - Anshu Kumar Sharma
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Muhammad Umair Naseem
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.,Doctoral School of Molecular Medicine, University of Debrecen, Debrecen, Hungary
| | - Attila Kovács
- Department of Radiation Therapy, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Jennifer Kugel
- Department of Biochemistry, University of Colorado, Boulder, USA
| | - Guido Kroemer
- Centre de Recherche Des Cordeliers, Equipe Labellisée Par La Ligue Contre Le Cancer, Université de Paris, Sorbonne Université, Inserm U1138, Institut Universitaire de France, Paris, France.,Metabolomics and Cell Biology Platforms, Institut Gustave Roussy, Villejuif, France.,Pôle de Biologie, Hôpital Européen Georges Pompidou, AP-HP, Paris, France
| | - György Vámosi
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.
| | - Gábor Szabó
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.
| |
Collapse
|
4
|
Singh RK, Mukherjee A. Molecular Mechanism of Dual Intercalation in Sac7d–DNA Complexation. J Phys Chem B 2022; 126:1682-1690. [DOI: 10.1021/acs.jpcb.1c09355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Reman Kumar Singh
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
- Department of Chemistry, Indian Institute of Technology, Bombay 400076, India
| | - Arnab Mukherjee
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| |
Collapse
|
5
|
Chowdhury R, Sai Sreyas Adury V, Vijay A, Singh RK, Mukherjee A. Atomistic De-novo Inhibitor Generation-Guided Drug Repurposing for SARS-CoV-2 Spike Protein with Free-Energy Validation by Well-Tempered Metadynamics. Chem Asian J 2021; 16:1634-1642. [PMID: 33949124 PMCID: PMC8207131 DOI: 10.1002/asia.202100268] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2021] [Revised: 04/20/2021] [Indexed: 12/13/2022]
Abstract
Computational drug design is increasingly becoming important with new and unforeseen diseases like COVID-19. In this study, we present a new computational de novo drug design and repurposing method and applied it to find plausible drug candidates for the receptor binding domain (RBD) of SARS-CoV-2 (COVID-19). Our study comprises three steps: atom-by-atom generation of new molecules around a receptor, structural similarity mapping to existing approved and investigational drugs, and validation of their binding strengths to the viral spike proteins based on rigorous all-atom, explicit-water well-tempered metadynamics free energy calculations. By choosing the receptor binding domain of the viral spike protein, we showed that some of our new molecules and some of the repurposable drugs have stronger binding to RBD than hACE2. To validate our approach, we also calculated the free energy of hACE2 and RBD, and found it to be in an excellent agreement with experiments. These pool of drugs will allow strategic repurposing against COVID-19 for a particular prevailing conditions.
Collapse
Affiliation(s)
- Rituparno Chowdhury
- Department of ChemistryIndian Institute of Science Education and Research PuneMaharashtraIndia
- Department of Chemical SciencesIndian Institute of Science Education and Research KolkataWest BengalIndia
| | | | - Amal Vijay
- Department of ChemistryIndian Institute of Science Education and Research PuneMaharashtraIndia
| | - Reman K. Singh
- Department of ChemistryIndian Institute of Science Education and Research PuneMaharashtraIndia
| | - Arnab Mukherjee
- Department of ChemistryIndian Institute of Science Education and Research PuneMaharashtraIndia
| |
Collapse
|
6
|
Singh RK, Mukherjee A. Molecular Mechanism of the Intercalation of the SOX-4 Protein into DNA Inducing Bends and Kinks. J Phys Chem B 2021; 125:3752-3762. [PMID: 33848164 DOI: 10.1021/acs.jpcb.0c11496] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
DNA-protein interactions regulate several biophysical functions, yet the mechanism of only a few is investigated in molecular detail. An important example is the intercalation of transcription factor proteins into DNA that produce bent and kinked DNA. Here, we have studied the molecular mechanism of the intercalation of a transcription factor SOX4 into DNA with a goal to understand the sequence of molecular events that precede the bending and kinking of the DNA. Our long well-tempered metadynamics and molecular dynamics (MD) simulations show that the protein primarily binds to the backbone of DNA and rotates around it to form an intercalative native state. We show that although there are multiple pathways for intercalation, the deintercalation pathway matches with the most probable intercalation pathway. In both cases, bending and kinking happen simultaneously, driven by the onset of the intercalation of the amino acid.
Collapse
Affiliation(s)
- Reman Kumar Singh
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| | - Arnab Mukherjee
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| |
Collapse
|
7
|
Narwade N, Patel S, Alam A, Chattopadhyay S, Mittal S, Kulkarni A. Mapping of scaffold/matrix attachment regions in human genome: a data mining exercise. Nucleic Acids Res 2019; 47:7247-7261. [PMID: 31265077 PMCID: PMC6698742 DOI: 10.1093/nar/gkz562] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Revised: 06/08/2019] [Accepted: 06/27/2019] [Indexed: 11/14/2022] Open
Abstract
Scaffold/matrix attachment regions (S/MARs) are DNA elements that serve to compartmentalize the chromatin into structural and functional domains. These elements are involved in control of gene expression which governs the phenotype and also plays role in disease biology. Therefore, genome-wide understanding of these elements holds great therapeutic promise. Several attempts have been made toward identification of S/MARs in genomes of various organisms including human. However, a comprehensive genome-wide map of human S/MARs is yet not available. Toward this objective, ChIP-Seq data of 14 S/MAR binding proteins were analyzed and the binding site coordinates of these proteins were used to prepare a non-redundant S/MAR dataset of human genome. Along with co-ordinate (location) details of S/MARs, the dataset also revealed details of S/MAR features, namely, length, inter-SMAR length (the chromatin loop size), nucleotide repeats, motif abundance, chromosomal distribution and genomic context. S/MARs identified in present study and their subsequent analysis also suggests that these elements act as hotspots for integration of retroviruses. Therefore, these data will help toward better understanding of genome functioning and designing effective anti-viral therapeutics. In order to facilitate user friendly browsing and retrieval of the data obtained in present study, a web interface, MARome (http://bioinfo.net.in/MARome), has been developed.
Collapse
Affiliation(s)
- Nitin Narwade
- Bioinformatics Centre, Savitribai Phule Pune University, Pune - 411 007, Maharashtra, India
| | - Sonal Patel
- Chromatin and Disease Biology Lab, National Centre for Cell Science, Pune - 411 007, Maharashtra, India
| | - Aftab Alam
- Chromatin and Disease Biology Lab, National Centre for Cell Science, Pune - 411 007, Maharashtra, India
| | - Samit Chattopadhyay
- Chromatin and Disease Biology Lab, National Centre for Cell Science, Pune - 411 007, Maharashtra, India.,Indian Institute of Chemical Biology, 4, Raja S.C. Mullick Road, Jadavpur, Kolkata - 700 032, West Bengal, India
| | - Smriti Mittal
- Department of Biotechnology, Savitribai Phule Pune University, Pune - 411 007, Maharashtra, India
| | - Abhijeet Kulkarni
- Bioinformatics Centre, Savitribai Phule Pune University, Pune - 411 007, Maharashtra, India
| |
Collapse
|
8
|
Intercalation of small molecules into DNA in chromatin is primarily controlled by superhelical constraint. PLoS One 2019; 14:e0224936. [PMID: 31747397 PMCID: PMC6867626 DOI: 10.1371/journal.pone.0224936] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Accepted: 09/30/2019] [Indexed: 12/22/2022] Open
Abstract
The restricted access of regulatory factors to their binding sites on DNA wrapped around the nucleosomes is generally interpreted in terms of molecular shielding exerted by nucleosomal structure and internucleosomal interactions. Binding of proteins to DNA often includes intercalation of hydrophobic amino acids into the DNA. To assess the role of constrained superhelicity in limiting these interactions, we studied the binding of small molecule intercalators to chromatin in close to native conditions by laser scanning cytometry. We demonstrate that the nucleosome-constrained superhelical configuration of DNA is the main barrier to intercalation. As a result, intercalating compounds are virtually excluded from the nucleosome-occupied regions of the chromatin. Binding of intercalators to extranucleosomal regions is limited to a smaller degree, in line with the existence of net supercoiling in the regions comprising linker and nucleosome free DNA. Its relaxation by inducing as few as a single nick per ~50 kb increases intercalation in the entire chromatin loop, demonstrating the possibility for long-distance effects of regulatory potential.
Collapse
|