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
|
Luo S, Xiong D, Zhao X, Duan L. An Attempt of Seeking Favorable Binding Free Energy Prediction Schemes Considering the Entropic Effect on Fis-DNA Binding. J Phys Chem B 2023; 127:1312-1324. [PMID: 36735878 DOI: 10.1021/acs.jpcb.2c07811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Protein-DNA binding mechanisms in a complex manner are essential for understanding many biological processes. Over the past decades, numerous experiments and calculations have analyzed the specificity of protein-DNA binding. However, the accuracy of binding free energy prediction for multi-base DNA systems still needs to be improved. Fis is a DNA-binding protein that regulates various transcription and recombination reactions. In the present work, we tested several methods of predict binding free energy based on this system to find a favorable prediction scheme and explore the binding mechanism of Fis protein and DNA. Two solvent models (explicit and implicit solvent models) were chosen for the dynamics process, and the predicted binding free energy was more accurate under the explicit solvent model. When different Poisson-Boltzmann/Generalized Born (PB/GB) models were tested for DNA force fields (BSC1 and OL15), it was found that the binding free energy predicted by the selected OL15 force field performed better and the correlation between predicted and experimental values was improved with the increasing interior dielectric constant (Dk). Finally, using Dk = 8, the GBOBC1 model combined with interaction entropy (IE), which was calculated for entropic contribution (GBOBC1_IE_8), was screened out for the binding free energy prediction and analysis of the Fis-DNA system, and the validity of the method was further verified by testing the Cren7-DNA system. By performing conformational analysis of the minor groove, it was found that mutation of the DNA central sequence A/T to C/G and deletion of the guanine 2-amino group would change the minor groove width and thus affect the formation of the major groove, altering the interaction and atomic contact between the protein and the major groove, thus changing the binding affinity of Fis and DNA. Hopefully, the series of tests in this work can shed some light on the related studies of protein and DNA systems.
Collapse
Affiliation(s)
- Song Luo
- School of Physics and Electronics, Shandong Normal University, Jinan, Shandong250014, China
| | - Danyang Xiong
- School of Physics and Electronics, Shandong Normal University, Jinan, Shandong250014, China
| | - Xiaoyu Zhao
- School of Physics and Electronics, Shandong Normal University, Jinan, Shandong250014, China
| | - Lili Duan
- School of Physics and Electronics, Shandong Normal University, Jinan, Shandong250014, China
| |
Collapse
|
3
|
Ameen F, Siddiqui S, Jahan I, Nayeem SM, Rehman SU, Tabish M. Studying the interaction of scopolamine with calf-thymus DNA: An in-vitro and in-silico approach and genotoxicity. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2022; 265:120391. [PMID: 34571375 DOI: 10.1016/j.saa.2021.120391] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Revised: 09/03/2021] [Accepted: 09/07/2021] [Indexed: 06/13/2023]
Abstract
Scopolamine is used to treat various CNS disorder like urinary incontinence, motion sickness, spasmic movements. Despite its pharmaceutical properties, its interaction with DNA is not yet reported. In this article, the interaction between scopolamine and ct-DNA is reported using a combination of biophysical techniques. UV-visible and steady-state fluorescence spectroscopy were used to study interaction and complex formation. Competitive displacement assays and potassium iodide quenching confirmed the mode of binding between scopolamine and DNA. Structural changes induced in the ct-DNA in the presence of scopolamine were evaluated by CD spectroscopy. The plasmid nicking and NBT assay confirmed the genotoxic effect of scopolamine. In-silico study by molecular docking and molecular dynamics simulation revealed the mode of interaction, major stabilizing forces as well as the nucleotide sequences to which the scopolamine binds.
Collapse
Affiliation(s)
- Faisal Ameen
- Department of Biochemistry, Faculty of Life Sciences, A.M. University, Aligarh, U.P. 202002, India
| | - Sharmin Siddiqui
- Department of Biochemistry, Faculty of Life Sciences, A.M. University, Aligarh, U.P. 202002, India
| | - Ishrat Jahan
- Department of Chemistry, Faculty of Science, A.M. University, Aligarh, U.P. 202002, India
| | - Shahid M Nayeem
- Department of Chemistry, Faculty of Science, A.M. University, Aligarh, U.P. 202002, India
| | - Sayeed Ur Rehman
- Department of Biochemistry, School of Chemical and Life Sciences, Jamia Hamdard, New Delhi 110062, India
| | - Mohammad Tabish
- Department of Biochemistry, Faculty of Life Sciences, A.M. University, Aligarh, U.P. 202002, India.
| |
Collapse
|
4
|
Maddah M, Shahabi M, Peyvandi K. How Does DcAFP, a Plant Antifreeze Protein, Control Ice Inhibition through the Kelvin Effect? Ind Eng Chem Res 2021. [DOI: 10.1021/acs.iecr.1c02559] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mina Maddah
- Department of Chemistry, K.N. Toosi University of Technology, 1969764499 Tehran, Iran
- Super Computing Institute, University of Tehran, 1417935840 Tehran, Iran
| | - Maryam Shahabi
- Faculty of Chemical, Petroleum and Gas Engineering, Semnan University, 3513119111 Semnan, Iran
| | - Kiana Peyvandi
- Faculty of Chemical, Petroleum and Gas Engineering, Semnan University, 3513119111 Semnan, Iran
| |
Collapse
|
5
|
Ghadari R, Mohsenzadeh E. Effect of COF Presence on DNA Molecular Interactions: A QM/MM and MD Simulations Study. ChemistrySelect 2021. [DOI: 10.1002/slct.202102157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Rahim Ghadari
- Computational Chemistry Laboratory Department of Organic and Biochemistry Faculty of Chemistry University of Tabriz Tabriz Iran, P.O. 5166616471
| | - Enayat Mohsenzadeh
- Computational Chemistry Laboratory Department of Organic and Biochemistry Faculty of Chemistry University of Tabriz Tabriz Iran, P.O. 5166616471
| |
Collapse
|
6
|
Maddah M, Karami L. An atomistic investigation on the interaction of distamycin A and its derivative with the telomeric G-Quadruplex as anticancer agents by molecular dynamics simulation. Arch Biochem Biophys 2021; 701:108797. [PMID: 33607110 DOI: 10.1016/j.abb.2021.108797] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Revised: 02/05/2021] [Accepted: 02/06/2021] [Indexed: 10/22/2022]
Abstract
Human telomerase that activates within cancer cells has a telomeric sequence at the 3' end. Each factor that stabilizes the G-quadruplex in guanine-rich telomeric sequences can inhibit the regular telomerase activity. Therefore, the telomeric G-quadruplex is known as a promising target in cancer treatment. In this work, we studied the binding of positively charged distamycin A and its uncharged derivative to the G-quadruplex in a solution environment by Molecular Dynamics (MD) simulation. The binding mechanism and subtle conformational changes were investigated as a result of the ligand attachment. Moreover, binding free energy and clustering analysis describe the stability and flexibility of G-quadruplexes upon ligand binding. Structural analyses displayed that the favorable binding of both ligands imposes significant stability and rigidity in G-quadruplex conformation compared to free G-quadruplex, especially charged distamycin. Hydration pattern and ion distribution were different for free G-quadruplex and both of the ligand complexes. Energy decomposition reveals the electrostatic effect on the stability of G-quadruplex. The radial distribution function displayed the solvent shell and ion moving away from the groove. The hydrogen bond played an essential role in the binding of both ligands, especially for the charged derivative. van der Waals interaction is the only factor that is more important in binding uncharged distamycin into G-quadruplex than the charged one. The calculated ΔGbind showed the stability of both ligands within grooves and good agreement with the experimental binding free energy data. Finally, the results suggest that ligand modification improves the binding mode toward stabilizing G-quadruplexes.
Collapse
Affiliation(s)
- Mina Maddah
- Depatment of Chemistry, K. N. Toosi University of Technology, Tehran, Iran; Super Computing Institute, University of Tehran, Tehran, Iran.
| | - Leila Karami
- Department of Cell and Molecular Biology, Faculty of Biological Sciences, Kharazmi University, Tehran, Iran.
| |
Collapse
|
7
|
Liu F, Zhang Y, Wang H, Li L, Zhao W, Shen JW, Liang L. Study on the adsorption orientation of DNA on two-dimensional MoS2 surface via molecular dynamics simulation: A vertical orientation phenomenon. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2019.110546] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
8
|
Jawad B, Poudel L, Podgornik R, Steinmetz NF, Ching WY. Molecular mechanism and binding free energy of doxorubicin intercalation in DNA. Phys Chem Chem Phys 2019; 21:3877-3893. [PMID: 30702122 DOI: 10.1039/c8cp06776g] [Citation(s) in RCA: 54] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of ΔGele, ΔGvDW, ΔGpolar, ΔGnon-polar, TΔStotal, ΔGdeform, ΔGcon, and ΔGion. To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (ΔGdeform), the energy cost from translational and rotational entropic contributions (TΔStran+rot), the total electrostatic interactions (ΔGpolar-PB/GB + ΔGele) of incorporation, the intramolecular electrostatic interactions (ΔGele) and electrostatic polar solvation interactions (ΔGpolar-PB/GB) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (ΔGvDW), the non-polar solvation interaction (ΔGnon-polar), the vibrational entropic contribution (TΔSvib), and the standard concentration dependent free energies of DOX (ΔGcon) and the ionic solution (ΔGion). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX-DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. We have calculated the total binding free energy (BFE) (ΔGt-sim) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX-DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX-DNA1 and DOX-DNA2 in the 1 : 1 complex is -9.1 and -5.1 kcal mol-1 respectively. The same quantities from the modified MM-PBSA method are -12.74 and -8.35 kcal mol-1 respectively. The value of the total BFE ΔGt-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of -7.7 ± 0.3 kcal mol-1.
Collapse
Affiliation(s)
- Bahaa Jawad
- Department of Physics and Astronomy, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA.
| | | | | | | | | |
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
|
Chen Y, Wang X, Zhang K, Bao Y. Intercalation of cytosine into Eu3+-doped hydrocalumite and their fluorescent responses. JOURNAL OF THE IRANIAN CHEMICAL SOCIETY 2017. [DOI: 10.1007/s13738-017-1176-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|