1
|
Zhu JY, Liu Q, Jiang XN, Zheng XH, Wang L, Hao Q, Wang CS. From bonds to interactions: comprehensive molecular characterization via polarizable bond-dipole approach. Phys Chem Chem Phys 2023; 25:29867-29880. [PMID: 37888898 DOI: 10.1039/d3cp04060g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2023]
Abstract
Accurately characterizing molecular interactions stands as a pivotal requirement for ensuring the reliability of molecular dynamics simulations. In line with our bond-dipole-based interaction model proposed by Gao et al. [X.-C. Gao, Q. Hao and C.-S. Wang, J. Chem. Theory Comput., 2017, 13, 2730-2741.], we have implemented an efficient and concise approach to compute electrostatic potential. This methodology capitalizes on the polarizable nature of chemical bond dipoles, resulting in a model of remarkable simplicity. In this study, we have revised the polarizable bond-dipole-based force field (PBFF) through the meticulous curation of quantum chemical data sets. These data sets encompass a comprehensive collection of 40 000 conformations, including those of water, methylamine, methanol, and N-methylacetamide. Additionally, we incorporate 520 hydrogen-bonded dimers into our data sets. In pursuit of enhanced accuracy in molecular dynamics simulations and a more faithful representation of potential energy landscapes, we undertook the re-optimization of the nonbonded parameters within the PBFF framework. Concurrently, we intricately fine-tuned the bonded parameters. The results of our comprehensive evaluation denote that this newly optimized force field method adeptly and efficiently computes structural characteristics, harmonic frequencies, and interaction energies. Overall, this study provides further validation for the applicability of PBFF in molecular dynamics simulations.
Collapse
Affiliation(s)
- Jia-Yi Zhu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, China.
| | - Qi Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, China.
| | - Xiao-Nan Jiang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, China.
| | - Xiao-Han Zheng
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, China.
| | - Lei Wang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, China.
| | - Qiang Hao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, China.
| | - Chang-Sheng Wang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, China.
| |
Collapse
|
2
|
Weng J, Yang S, Shen J, Liu H, Xu Y, Hao D, Wang S. Molecular dynamics simulation reveals DNA-specific recognition mechanism via c-Myb in pseudo-palindromic consensus of mim-1 promoter. J Zhejiang Univ Sci B 2023; 24:883-895. [PMID: 37752090 PMCID: PMC10522569 DOI: 10.1631/jzus.b2200634] [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: 12/07/2022] [Accepted: 02/07/2023] [Indexed: 09/28/2023]
Abstract
This study aims to gain insight into the DNA-specific recognition mechanism of c-Myb transcription factor during the regulation of cell early differentiation and proliferation. Therefore, we chose the chicken myeloid gene, mitochondrial import protein 1 (mim-1), as a target to study the binding specificity between potential dual-Myb-binding sites. The c-Myb-binding site in mim-1 is a pseudo-palindromic sequence AACGGTT, which contains two AACNG consensuses. Simulation studies in different biological scenarios revealed that c-Myb binding with mim-1 in the forward strand (complex F) ismore stable than that inthereverse strand (complex R). The principal component analysis (PCA) dynamics trajectory analyses suggested an opening motion of the recognition helices of R2 and R3 (R2R3), resulting in the dissociation of DNA from c-Myb in complex R at 330 K, triggered by the reduced electrostatic potential on the surface of R2R3. Furthermore, the DNA confirmation and hydrogen-bond interaction analyses indicated that the major groove width of DNA increased in complex R, which affected on the hydrogen-bond formation ability between R2R3 and DNA, and directly resulted in the dissociation of DNA from R2R3. The steered molecular dynamics (SMD) simulation studies also suggested that the electrostatic potential, major groove width, and hydrogen bonds made major contribution to the DNA-specific recognition. In vitro trials confirmed the simulation results that c-Myb specifically bound to mim-1 in the forward strand. This study indicates that the three-dimensional (3D) structure features play an important role in the DNA-specific recognition mechanism by c-Myb besides the AACNG consensuses, which is beneficial to understanding the cell early differentiation and proliferation regulated by c-Myb, as well as the prediction of novel c-Myb-binding motifs in tumorigenesis.
Collapse
Affiliation(s)
- Jinru Weng
- Stomatology Hospital, School of Stomatology, Zhejiang University School of Medicine, Zhejiang Provincial Clinical Research Center for Oral Diseases, Key Laboratory of Oral Biomedical Research of Zhejiang Province, Cancer Center of Zhejiang University, Engineering Research Center of Oral Biomaterials and Devices of Zhejiang Province, Hangzhou 310000, China
| | - Shuo Yang
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, Jilin University, Changchun 130021, China
| | - Jinkang Shen
- Stomatology Hospital, School of Stomatology, Zhejiang University School of Medicine, Zhejiang Provincial Clinical Research Center for Oral Diseases, Key Laboratory of Oral Biomedical Research of Zhejiang Province, Cancer Center of Zhejiang University, Engineering Research Center of Oral Biomaterials and Devices of Zhejiang Province, Hangzhou 310000, China
| | - Hongsen Liu
- Stomatology Hospital, School of Stomatology, Zhejiang University School of Medicine, Zhejiang Provincial Clinical Research Center for Oral Diseases, Key Laboratory of Oral Biomedical Research of Zhejiang Province, Cancer Center of Zhejiang University, Engineering Research Center of Oral Biomaterials and Devices of Zhejiang Province, Hangzhou 310000, China
| | - Yuzi Xu
- Stomatology Hospital, School of Stomatology, Zhejiang University School of Medicine, Zhejiang Provincial Clinical Research Center for Oral Diseases, Key Laboratory of Oral Biomedical Research of Zhejiang Province, Cancer Center of Zhejiang University, Engineering Research Center of Oral Biomaterials and Devices of Zhejiang Province, Hangzhou 310000, China
| | - Dongyun Hao
- Institute of Agricultural Biotechnology, Jilin Academy of Agricultural Sciences (JAAS), Changchun 130033, China.
| | - Shan Wang
- Stomatology Hospital, School of Stomatology, Zhejiang University School of Medicine, Zhejiang Provincial Clinical Research Center for Oral Diseases, Key Laboratory of Oral Biomedical Research of Zhejiang Province, Cancer Center of Zhejiang University, Engineering Research Center of Oral Biomaterials and Devices of Zhejiang Province, Hangzhou 310000, China.
| |
Collapse
|
3
|
Hanke M, Hansen N, Tomm E, Grundmeier G, Keller A. Time-Dependent DNA Origami Denaturation by Guanidinium Chloride, Guanidinium Sulfate, and Guanidinium Thiocyanate. Int J Mol Sci 2022; 23:ijms23158547. [PMID: 35955680 PMCID: PMC9368935 DOI: 10.3390/ijms23158547] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 07/29/2022] [Accepted: 07/29/2022] [Indexed: 11/16/2022] Open
Abstract
Guanidinium (Gdm) undergoes interactions with both hydrophilic and hydrophobic groups and, thus, is a highly potent denaturant of biomolecular structure. However, our molecular understanding of the interaction of Gdm with proteins and DNA is still rather limited. Here, we investigated the denaturation of DNA origami nanostructures by three Gdm salts, i.e., guanidinium chloride (GdmCl), guanidinium sulfate (Gdm2SO4), and guanidinium thiocyanate (GdmSCN), at different temperatures and in dependence of incubation time. Using DNA origami nanostructures as sensors that translate small molecular transitions into nanostructural changes, the denaturing effects of the Gdm salts were directly visualized by atomic force microscopy. GdmSCN was the most potent DNA denaturant, which caused complete DNA origami denaturation at 50 °C already at a concentration of 2 M. Under such harsh conditions, denaturation occurred within the first 15 min of Gdm exposure, whereas much slower kinetics were observed for the more weakly denaturing salt Gdm2SO4 at 25 °C. Lastly, we observed a novel non-monotonous temperature dependence of DNA origami denaturation in Gdm2SO4 with the fraction of intact nanostructures having an intermediate minimum at about 40 °C. Our results, thus, provide further insights into the highly complex Gdm–DNA interaction and underscore the importance of the counteranion species.
Collapse
|
4
|
Hanke M, Dornbusch D, Hadlich C, Rossberg A, Hansen N, Grundmeier G, Tsushima S, Keller A, Fahmy K. Anion-specific structure and stability of guanidinium-bound DNA origami. Comput Struct Biotechnol J 2022; 20:2611-2623. [PMID: 35685373 PMCID: PMC9163702 DOI: 10.1016/j.csbj.2022.05.037] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Revised: 05/19/2022] [Accepted: 05/19/2022] [Indexed: 12/13/2022] Open
Abstract
While the folding of DNA into rationally designed DNA origami nanostructures has been studied extensively with the aim of increasing structural diversity and introducing functionality, the fundamental physical and chemical properties of these nanostructures remain largely elusive. Here, we investigate the correlation between atomistic, molecular, nanoscopic, and thermodynamic properties of DNA origami triangles. Using guanidinium (Gdm) as a DNA-stabilizing but potentially also denaturing cation, we explore the dependence of DNA origami stability on the identity of the accompanying anions. The statistical analyses of atomic force microscopy (AFM) images and circular dichroism (CD) spectra reveals that sulfate and chloride exert stabilizing and destabilizing effects, respectively, already below the global melting temperature of the DNA origami triangles. We identify structural transitions during thermal denaturation and show that heat capacity changes ΔCp determine the temperature sensitivity of structural damage. The different hydration shells of the anions and their potential to form Gdm+ ion pairs in concentrated salt solutions modulate ΔCp by altered wetting properties of hydrophobic DNA surface regions as shown by molecular dynamics simulations. The underlying structural changes on the molecular scale become amplified by the large number of structurally coupled DNA segments and thereby find nanoscopic correlations in AFM images.
Collapse
|
5
|
Pisetsky DS, Shaffer R, Armstrong DD, Spencer DM. The interaction of anti-DNA antibodies with DNA antigen: Evidence for hysteresis for high avidity binding. Clin Immunol 2021; 231:108848. [PMID: 34492381 DOI: 10.1016/j.clim.2021.108848] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Revised: 07/11/2021] [Accepted: 09/02/2021] [Indexed: 10/20/2022]
Abstract
Antibodies to DNA (anti-DNA) are the serological hallmark of systemic lupus erythematosus. Previous studies have indicated that the phosphodiester backbone is the main antigenic target, with electrostatic interactions important for high avidity. To define further these interactions, the effects of ionic strength on anti-DNA binding of SLE plasmas were assessed in association and dissociation assays by ELISA. As these studies demonstrated, increasing ionic strength to a concentration of 1000 mM NaCl reduced antibody binding although the extent of the reduction varied among samples. In dissociation assays, differences among plasmas were also observed. For one of the plasmas, binding to DNA displayed resistance to dissociation by increasing ionic strength even though these concentrations limited binding in association assays. Time course studies showed a gradual change in binding interactions. These studies indicate that anti-DNA binding can involve both electrostatic and non-electrostatic interactions, with binding in some plasmas showing evidence of hysteresis.
Collapse
Affiliation(s)
- David S Pisetsky
- Medical Research Service, Durham Veterans Administration Medical Center, Durham, NC 27705, USA; Division of Rheumatology and Immunology, Duke University Medical Center, Durham, NC 27710, USA.
| | | | | | - Diane M Spencer
- Division of Rheumatology and Immunology, Duke University Medical Center, Durham, NC 27710, USA
| |
Collapse
|
6
|
Sarkar S, Singh PC. Sequence specific hydrogen bond of DNA with denaturants affects its stability: Spectroscopic and simulation studies. Biochim Biophys Acta Gen Subj 2020; 1865:129735. [PMID: 32946929 DOI: 10.1016/j.bbagen.2020.129735] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2020] [Revised: 08/21/2020] [Accepted: 09/13/2020] [Indexed: 12/31/2022]
Abstract
BACKGROUND Several different small molecules have been used to target the DNA helix in order to treat the diseases caused by its mutation. Guanidinium(Gdm+) and urea based drugs have been used for the diseases related to central nervous system, also as the anti-inflammatory and chemotherapeutic agent. However, the role of Gdm+ and urea in the stabilization/destabilization of DNA is not well understood. METHODS Spectroscopic techniques along with molecular dynamics (MD) simulation have been performed on different sequences of DNA in the presence of guanidinium chloride (GdmCl) and urea to decode the binding of denaturants with DNA and the role of hydrogen bond with the different regions of DNA in its stability/destability. RESULTS AND CONCLUSION Our study reveals that, Gdm+ of GdmCl and urea both intrudes into the groove region of DNA along with the interaction with its phosphate backbone. However, interaction of Gdm+ and urea with the nucleobases in the groove region is different. Gdm+ forms the intra-strand hydrogen bond with the central region of the both sequences of DNA whereas inter-strand hydrogen bond along with water assisted hydrogen bond takes place in the case of urea. The intra-strand hydrogen bond formation capability of Gdm+ with the nucleobases in the minor groove of DNA decreases its groove width which probably causes the stabilization of B-DNA in GdmCl. In contrast, the propensity of the formation of inter-strand hydrogen bond of urea with the nucleobases in the groove region of DNA without affecting the groove width destabilizes B-DNA as compared to GdmCl. This study depicts that the opposite effect of GdmCl and urea on the stability is a general property of B-DNA. However, the extent of stabilization/destabilization of DNA in Gdm+ and urea depend on its sequence probably due to the difference in the intra/inter-strand hydrogen bonding with different bases present in both the sequences of DNA. GENERAL SIGNIFICANCE The information obtained from this study will be useful for the designing of Gdm+ based drug molecule which can target the DNA more specifically and selectively.
Collapse
Affiliation(s)
- Sunipa Sarkar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 70032, India
| | - Prashant Chandra Singh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 70032, India.
| |
Collapse
|
7
|
Sarkar S, Rajdev P, Singh PC. Hydrogen bonding of ionic liquids in the groove region of DNA controls the extent of its stabilization: synthesis, spectroscopic and simulation studies. Phys Chem Chem Phys 2020; 22:15582-15591. [PMID: 32613973 DOI: 10.1039/d0cp01548b] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Ionic liquids (ILs) have been extensively used for stabilization and long-term DNA storage. However, molecular level understanding of the role of the hydrogen bond of DNA with ILs in its stabilization is still inadequate. Two ILs, namely, 1,1,3,3-tetramethylguanidinium acetate (TMG) and 2,2-diethyl-1,1,3,3-tetramethylguanidinium acetate (DETMG), have been synthesized, of which TMG has a hydrogen bonding N-H group whereas DETMG does not contain any hydrogen bonding site. It has been found that both TMG and DETMG cations interact in the groove region of DNA; however, their mode of interaction is distinctly different, which causes the stabilization of DNA in the presence of TMG, whereas the effect is opposite in the case of DETMG. It is apparent from the data that only the accommodation of ILs in the groove region is not enough for the stabilization of DNA. MD simulation and spectroscopic studies combinedly indicate that the hydrogen bonding capability of the TMG cation enhances the hydrogen bonding between the Watson-Crick base pairs of DNA, resulting in its stabilization. In contrast, the bigger size as well as the absence of the hydrogen bonding site of the DETMG cation perturbs the minor groove width and base pair step parameters of DNA during its intrusion into the minor groove, which decreases the hydrogen bond between the Watson-Crick base pairs of DNA, leading to destabilization.
Collapse
Affiliation(s)
- Sunipa Sarkar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Kolkata 700032, India.
| | | | | |
Collapse
|