1
|
Ferreira FHDC, Farrell NP, Costa LAS. Spermine and spermidine SI-PPCs: Molecular dynamics reveals enhanced biomolecular interactions. Int J Biol Macromol 2024; 278:134654. [PMID: 39128748 DOI: 10.1016/j.ijbiomac.2024.134654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2024] [Revised: 08/07/2024] [Accepted: 08/08/2024] [Indexed: 08/13/2024]
Abstract
In this paper the effects on the interaction of highly positively charged substitution-inert platinum polynuclear complexes (SI-PPCs) with negatively charged DNA and heparin are examined and compared by theoretical chemistry methods. Electrostatic and hydrogen bonding interactions contribute to the overall effects on the biomolecule. Root Mean Square (RMS) deviation, Solvent Accessible Surface, RMS fluctuation, and interaction analysis all confirm similar effects on both biomolecules, dictated predominantly by the total positive charge and total number of hydrogen bonds formed. Especially, changes in structural parameters suggesting condensation and reduction of available surface area will reduce or prevent normal protein recognition and may thus potentially inhibit biological mechanisms related to apoptosis (DNA) or reduced vascularization viability (HEP). Thermodynamic analyses supported these findings with favourable interaction energies. The comparison of DNA and heparin confirms the general intersectionality between the two biomolecules and confirms the intrinsic dual-nature function of this chemotype. The distinction between the two-limiting mode of actions (HS or DNA-centred) could reflect an intriguing balance between extracellular (GAG) and intracellular (DNA) binding and affinities. The results underline the need to fully understand GAG-small molecule interactions and their contribution to drug pharmacology and related therapeutic modalities. This report contributes to that understanding.
Collapse
Affiliation(s)
- Frederico Henrique do C Ferreira
- NEQC - Núcleo de Estudos em Química Computacional, Departamento de Química, ICE, Universidade Federal de Juiz de Fora, Juiz de Fora, Minas Gerais 36036-900, Brazil
| | - Nicholas P Farrell
- Department of Chemistry, Virginia Commonwealth University, Richmond, VA 23284-2006, USA
| | - Luiz Antônio S Costa
- NEQC - Núcleo de Estudos em Química Computacional, Departamento de Química, ICE, Universidade Federal de Juiz de Fora, Juiz de Fora, Minas Gerais 36036-900, Brazil.
| |
Collapse
|
2
|
Yuan S, Han X, Zhang J, Xie Z, Fan C, Xiao Y, Gao YQ, Yang YI. Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms Using Molecular Configuration Transformer. J Phys Chem A 2024; 128:4378-4390. [PMID: 38759697 DOI: 10.1021/acs.jpca.4c01267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2024]
Abstract
Theoretical studies on chemical reaction mechanisms have been crucial in organic chemistry. Traditionally, calculating the manually constructed molecular conformations of transition states for chemical reactions using quantum chemical calculations is the most commonly used method. However, this way is heavily dependent on individual experience and chemical intuition. In our previous study, we proposed a research paradigm that used enhanced sampling in molecular dynamics simulations to study chemical reactions. This approach can directly simulate the entire process of a chemical reaction. However, the computational speed limited the use of high-precision potential energy functions for simulations. To address this issue, we presented a scheme for training high-precision force fields for molecular modeling using a previously developed graph-neural-network-based molecular model, molecular configuration transformer. This potential energy function allowed for highly accurate simulations at a low computational cost, leading to more precise calculations of the mechanism of chemical reactions. We applied this approach to study a Claisen rearrangement reaction and a carbonyl insertion reaction catalyzed by manganese.
Collapse
Affiliation(s)
- Sihao Yuan
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Xu Han
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jun Zhang
- Changping Laboratory, Beijing 102200, China
| | - Zhaoxin Xie
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Cheng Fan
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Yunlong Xiao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
- Changping Laboratory, Beijing 102200, China
| | - Yi Isaac Yang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| |
Collapse
|
3
|
Patra P, Gao YQ. Structural and dynamical aspect of DNA motif sequence specific binding of AP-1 transcription factor. J Chem Phys 2024; 160:115103. [PMID: 38506297 DOI: 10.1063/5.0196508] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2024] [Accepted: 02/26/2024] [Indexed: 03/21/2024] Open
Abstract
Activator protein-1 (AP-1) comprises one of the largest and most evolutionary conserved families of ubiquitous eukaryotic transcription factors that act as a pioneer factor. Diversity in DNA binding interaction of AP-1 through a conserved basic-zipper (bZIP) domain directs in-depth understanding of how AP-1 achieves its DNA binding selectivity and consequently gene regulation specificity. Here, we address the structural and dynamical aspects of the DNA target recognition process of AP-1 using microsecond-long atomistic simulations based on the structure of the human AP-1 FosB/JunD bZIP-DNA complex. Our results show the unique role of DNA shape features in selective base specific interactions, characteristic ion population, and solvation properties of DNA grooves to form the motif sequence specific AP-1-DNA complex. The TpG step at the two terminals of the AP-1 site plays an important role in the structural adjustment of DNA by modifying the helical twist in the AP-1 bound state. We addressed the role of intrinsic motion of the bZIP domain in terms of opening and closing gripper motions of DNA binding helices, in target site recognition and binding of AP-1 factors. Our observations suggest that binding to the cognate motif in DNA is mainly accompanied with the precise adjustment of closing gripper motion of DNA binding helices of the bZIP domain.
Collapse
Affiliation(s)
- Piya Patra
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518107 Shenzhen, China
| | - Yi Qin Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518107 Shenzhen, China
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China
- Biomedical Pioneering Innovation Center, Peking University, 100871 Beijing, China
- Changping Laboratory, Beijing 102200, China
| |
Collapse
|
4
|
Xue Y, Wang J, He Y, Patra P, Gao YQ. Multi-scale gene regulation mechanism: Spatiotemporal transmission of genetic information. Curr Opin Struct Biol 2022; 77:102487. [PMID: 36274420 DOI: 10.1016/j.sbi.2022.102487] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2022] [Revised: 09/09/2022] [Accepted: 09/18/2022] [Indexed: 12/14/2022]
Abstract
Gene expression is regulated by many factors, including transcription factors, chromatin three-dimensional topology, modifications of DNA and histone proteins, and non-coding RNAs. The execution of these complex mechanisms requires an effectively coordinated regulation system. In this review, we emphasize that the multi-scale heterogeneous DNA sequence plays a fundamental and important role for gene expression activity and usage of different means of epigenetic regulation. We illustrate here that the chromatin structure organization provides a stage for spatiotemporal regulation between different genes or gene modules and to realize their downstream functional cooperation. Such a perspective expands our understanding of the central dogma: In addition to one-dimensional sequence information, inter-gene interactions can also be transferred from DNA and RNA to protein levels.
Collapse
Affiliation(s)
- Yue Xue
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jingyao Wang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yueying He
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Piya Patra
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Yi Qin Gao
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China; Biomedical Pioneering Innovation Center (BIOPIC), Peking University, Beijing 100871, China; Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China.
| |
Collapse
|
5
|
Wang SD, Zhang RB, Eriksson LA. Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base. Nucleic Acids Res 2022; 50:9072-9082. [PMID: 35979954 PMCID: PMC9458442 DOI: 10.1093/nar/gkac691] [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: 05/23/2022] [Revised: 07/15/2022] [Accepted: 07/30/2022] [Indexed: 12/24/2022] Open
Abstract
The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol-1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol-1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson-Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.
Collapse
Affiliation(s)
- Shu-dong Wang
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, South Street No. 5, Zhongguancun, Haidan District, 100081 Beijing, China
| | - Ru-bo Zhang
- Correspondence may also be addressed to Ru-bo Zhang.
| | - Leif A Eriksson
- To whom correspondence should be addressed. Tel: +46 31 786 9117;
| |
Collapse
|
6
|
Satange R, Kao SH, Chien CM, Chou SH, Lin CC, Neidle S, Hou MH. Staggered intercalation of DNA duplexes with base-pair modulation by two distinct drug molecules induces asymmetric backbone twisting and structure polymorphism. Nucleic Acids Res 2022; 50:8867-8881. [PMID: 35871296 PMCID: PMC9410880 DOI: 10.1093/nar/gkac629] [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/26/2022] [Revised: 06/24/2022] [Accepted: 07/22/2022] [Indexed: 12/12/2022] Open
Abstract
The use of multiple drugs simultaneously targeting DNA is a promising strategy in cancer therapy for potentially overcoming single drug resistance. In support of this concept, we report that a combination of actinomycin D (ActD) and echinomycin (Echi), can interact in novel ways with native and mismatched DNA sequences, distinct from the structural effects produced by either drug alone. Changes in the former with GpC and CpG steps separated by a A:G or G:A mismatch or in a native DNA with canonical G:C and C:G base pairs, result in significant asymmetric backbone twists through staggered intercalation and base pair modulations. A wobble or Watson–Crick base pair at the two drug-binding interfaces can result in a single-stranded ‘chair-shaped’ DNA duplex with a straight helical axis. However, a novel sugar-edged hydrogen bonding geometry in the G:A mismatch leads to a ‘curved-shaped’ duplex. Two non-canonical G:C Hoogsteen base pairings produce a sharply kinked duplex in different forms and a four-way junction-like superstructure, respectively. Therefore, single base pair modulations on the two drug-binding interfaces could significantly affect global DNA structure. These structures thus provide a rationale for atypical DNA recognition via multiple DNA intercalators and a structural basis for the drugs’ potential synergetic use.
Collapse
Affiliation(s)
- Roshan Satange
- Institute of Genomics and Bioinformatics, National Chung Hsing University , Taichung 402, Taiwan
- Ph.D. Program in Medical Biotechnology, National Chung Hsing University , Taichung 402, Taiwan
| | - Shih-Hao Kao
- Institute of Biotechnology, National Chung Hsing University , Taichung 402, Taiwan
| | - Ching-Ming Chien
- Institute of Genomics and Bioinformatics, National Chung Hsing University , Taichung 402, Taiwan
| | - Shan-Ho Chou
- Institute of Biochemistry, National Chung Hsing University , Taichung 402, Taiwan
| | - Chi-Chien Lin
- Institute of Biomedical Science, National Chung Hsing University , Taichung 402, Taiwan
| | - Stephen Neidle
- The School of Pharmacy, University College London , London WC1N 1AX, United Kingdom
| | - Ming-Hon Hou
- Institute of Genomics and Bioinformatics, National Chung Hsing University , Taichung 402, Taiwan
- Ph.D. Program in Medical Biotechnology, National Chung Hsing University , Taichung 402, Taiwan
- Institute of Biotechnology, National Chung Hsing University , Taichung 402, Taiwan
| |
Collapse
|
7
|
Mondal M, Gao YQ. Microscopic Insight into pH-Dependent Conformational Dynamics and Noncanonical Base Pairing in Telomeric i-Motif DNA. J Phys Chem Lett 2022; 13:5109-5115. [PMID: 35657602 DOI: 10.1021/acs.jpclett.2c00640] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Gene regulatory functions of noncanonical i-motif DNA are associated with dynamic i-motif formation in the cellular environment and pH variation. With atomistic simulations, we show the dramatic influence of solvent pH on the conformational dynamics of biologically relevant telomeric i-motif DNA coupled with protonation of cytosine bases in different conformations. We rationalized the pH-dependent dynamics and conformational variability of the i-motif in terms of base pairing and specific loop motions. The human telomeric i-motif is found to acquire various metastable folded conformations at pH values near the pKa of cytosine with the formation of a noncanonical C:C W:W trans base pair along with the hemiprotonated C:C+ pairs in the i-motif core. pH-dependent dynamics and the local solvent structure of i-motif DNA imply that the presence of a cosolvent or molecular crowding can promote i-motif formation in vivo by changing the conformational fluctuations and hydration state of the structure.
Collapse
Affiliation(s)
- Manas Mondal
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518107 Shenzhen, China
| | - Yi Qin Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518107 Shenzhen, China
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China
- Biomedical Pioneering Innovation Center, Peking University, 100871 Beijing, China
- Beijing Advanced Innovation Center for Genomics, Peking University, 100871 Beijing, China
| |
Collapse
|
8
|
Wang SD, Eriksson LA, Zhang RB. Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond. J Chem Inf Model 2022; 62:386-398. [PMID: 34994562 PMCID: PMC8790752 DOI: 10.1021/acs.jcim.1c01169] [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] [Indexed: 11/28/2022]
Abstract
![]()
Damaged or mismatched
DNA bases are normally thought to be able
to flip out of the helical stack, providing enzymes with access to
the faulty genetic information otherwise hidden inside the helix.
Thymine glycol (Tg) is one of the most common products of nucleic
acid damage. However, the static and dynamic structures of DNA duplexes
affected by 5R-Tg epimers are still not clearly understood, including
the ability of these to undergo spontaneous base flipping. Structural
effects of the 5R-Tg epimers on the duplex DNA are herein studied
using molecular dynamics together with reliable DFT based calculations.
In comparison with the corresponding intact DNA, the cis-5R,6S-Tg epimer base causes little perturbation to the duplex DNA,
and a barrier of 4.9 kcal mol–1 is obtained by meta-eABF
for cis-5R,6S-Tg base flipping out of the duplex
DNA, comparable to the 5.4 kcal mol–1 obtained for
the corresponding thymine flipping in intact DNA. For the trans-5R,6R-Tg epimer, three stable local structures were
identified, of which the most stable disrupts the Watson–Crick
hydrogen-bonded G5/C20 base pair, leading to conformational distortion
of the duplex. Interestingly, the relative barrier height of the 5R-Tg
flipping is only 1.0 kcal mol–1 for one of these trans-5R,6R-Tg epimers. Water bridge interactions were identified
to be essential for 5R-Tg flipping. The study clearly demonstrates
the occurrence of partial trans-5R,6R-Tg epimer flipping
in solution.
Collapse
Affiliation(s)
- Shu Dong Wang
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, South Street no 5, Zhongguancun, Haidian District, 100081 Beijing, China
| | - Leif A Eriksson
- Department of Chemistry and Molecular Biology, University of Gothenburg, Medicinaregatan 9c, 405 30 Göteborg, Sweden
| | - Ru Bo Zhang
- School of Chemistry and Chemical Engineering, Beijing Institute of Technology, South Street no 5, Zhongguancun, Haidian District, 100081 Beijing, China
| |
Collapse
|