1
|
Wang R, Lin Y, Sun Y, Zhao B, Chen L. Insight into the molecular recognition of human and polar bear pregnane X receptor by three organic pollutants using molecular docking and molecular dynamics simulations. ENVIRONMENT INTERNATIONAL 2024; 190:108926. [PMID: 39098090 DOI: 10.1016/j.envint.2024.108926] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2024] [Revised: 07/30/2024] [Accepted: 07/30/2024] [Indexed: 08/06/2024]
Abstract
Pregnane X receptor (PXR) is a heterologous biosensor that is involved in the metabolic pathway of environmental pollutants, regulating the transcription of genes involved in biotransformation. There are significant differences in the selectivity and specificity of organic pollutants (OPs) toward polar bear PXR (pbPXR) and human PXR (hPXR), but the detailed dynamical characteristics of their interactions are unclear. Homology Modeling, molecular docking, molecular dynamics simulation, and free energy calculation were used to analyze the recognition of pbPXR and hPXR by three OPs: BPA, chlordane and toxaphene. Comparing interaction patterns along with binding free energy of pbPXR and hPXR with these three OPs revealed that although pbPXR and hPXR interact similar with these three OPs, these OPs have different effects on the internal dynamics of pbPXR and hPXR. This results in significant alterations in the interaction of key residues near Leu209, Met243, Phe288, Met323, and His407 with OPs, thereby influencing their binding energy. Non-polar interactions, especially van der Waals interactions, were found to be the dominating factors in interacting of these OPs with PXRs. The region surrounding these key residues facilitates hydrophobic contacts with PXR, which are crucial for the selective activation of PXRs in different species by these three OPs. These findings are of significant guidance in understanding the impacts of environmental endocrine disruptors on different organisms.
Collapse
Affiliation(s)
- Ruige Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, China
| | - Yaqi Lin
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, China
| | - Ying Sun
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, China
| | - Bing Zhao
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, China; Heilongjiang Provincial Key Laboratory of Surface Active Agent and Auxiliary, Qiqihar University, Qiqihar, 161006, China
| | - Lin Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, China.
| |
Collapse
|
2
|
Bu F, Chen L, Sun Y, Zhao B, Wang R. Insight into the Binding Interaction between PEDCs and hERRγ Utilizing Molecular Docking and Molecular Dynamics Simulations. Molecules 2024; 29:3256. [PMID: 39064835 PMCID: PMC11278984 DOI: 10.3390/molecules29143256] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2024] [Revised: 07/03/2024] [Accepted: 07/03/2024] [Indexed: 07/28/2024] Open
Abstract
Phenolic environmental endocrine-disrupting chemicals (PEDCs) are persistent EDCs that are widely found in food packaging materials and environmental media and seriously threaten human health and ecological security. Human estrogen-related receptor γ (hERRγ) has been proposed as a mediator for the low-dose effects of many environmental PEDCs; however, the atomic-level descriptions of dynamical structural features and interactions of hERRγ and PEDCs are still unclarified. Herein, how three PEDCs, 4-(1-methylpropyl)phenol (4-sec-butylphenol), 5,6,7,8-tetrahydro-2-naphthol (tetrahydro-2-napthol), and 2,2-bis(4-hydroxy-3,5-dimethoxyphenyl)propane (BP(2,2)(Me)), interact with hERRγ to produce its estrogenic disruption effects was studied. Molecular docking and multiple molecular dynamics (MD) simulations were first conducted to distinguish the detailed interaction pattern of hERRγ with PEDCs. These binding structures revealed that residues around Leu271, Leu309, Leu345, and Phe435 are important when binding with PEDCs. Furthermore, the binding energies of PEDCs with hERRγ were also characterized using the molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methods, and the results showed that the interactions of CH-π, π-π, and hydrogen bonds are the major contributors for hERRγ binding to these three PEDCs. What is striking is that the methoxide groups of BP(2,2)(Me), as hydrophobic groups, can help to reduce the binding energy of PEDCs binding with hERRγ. These results provide important guidance for further understanding the influence of PEDCs on human health problems.
Collapse
Affiliation(s)
- Fanqiang Bu
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, China; (F.B.); (Y.S.); (B.Z.)
- State Key Laboratory of Organic-Inorganic Composites, Beijing Laboratory of Biomedical Materials, Beijing University of Chemical Technology, Beijing 100029, China
| | - Lin Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, China; (F.B.); (Y.S.); (B.Z.)
| | - Ying Sun
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, China; (F.B.); (Y.S.); (B.Z.)
| | - Bing Zhao
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, China; (F.B.); (Y.S.); (B.Z.)
- Heilongjiang Provincial Key Laboratory of Surface Active Agent and Auxiliary, Qiqihar University, Qiqihar 161006, China
| | - Ruige Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, China; (F.B.); (Y.S.); (B.Z.)
| |
Collapse
|
3
|
Frezza E, Laage D, Duboué-Dijon E. Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences. J Phys Chem Lett 2024; 15:4351-4358. [PMID: 38619551 DOI: 10.1021/acs.jpclett.4c00629] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/16/2024]
Abstract
Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.
Collapse
Affiliation(s)
- Elisa Frezza
- Université Paris Cité, CNRS, CiTCoM, Paris 75006, France
| | - Damien Laage
- PASTEUR, Department of Chemistry, École Normale Supérieure-PSL, Sorbonne Université, CNRS, Paris 75005, France
| | - Elise Duboué-Dijon
- Université Paris Cité, CNRS, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, Paris 75005, France
| |
Collapse
|
4
|
Rinaldi S, Moroni E, Rozza R, Magistrato A. Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation. J Chem Theory Comput 2024; 20:993-1018. [PMID: 38287883 DOI: 10.1021/acs.jctc.3c01239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2024]
Abstract
Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.
Collapse
Affiliation(s)
- Silvia Rinaldi
- National Research Council of Italy (CNR) - Institute of Chemistry of OrganoMetallic Compounds (ICCOM), c/o Area di Ricerca CNR di Firenze Via Madonna del Piano 10, 50019 Sesto Fiorentino, Florence, Italy
| | - Elisabetta Moroni
- National Research Council of Italy (CNR) - Institute of Chemical Sciences and Technologies (SCITEC), via Mario Bianco 9, 20131 Milano, Italy
| | - Riccardo Rozza
- National Research Council of Italy (CNR) - Institute of Material Foundry (IOM) c/o International School for Advanced Studies (SISSA), Via Bonomea, 265, 34136 Trieste, Italy
| | - Alessandra Magistrato
- National Research Council of Italy (CNR) - Institute of Material Foundry (IOM) c/o International School for Advanced Studies (SISSA), Via Bonomea, 265, 34136 Trieste, Italy
| |
Collapse
|
5
|
Chen J, Wang W, Sun H, He W. Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters. Mini Rev Med Chem 2024; 24:1323-1333. [PMID: 38265367 DOI: 10.2174/0113895575252165231122095555] [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: 03/06/2023] [Revised: 09/05/2023] [Accepted: 10/16/2023] [Indexed: 01/25/2024]
Abstract
Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan-250357, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan-250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan-250357, China
| | - Weikai He
- School of Science, Shandong Jiaotong University, Jinan-250357, China
| |
Collapse
|
6
|
Sun J, Liang S, Liu X, Zhang S, Li M, Zhang Q, Chen J. Insights into the selectivity of a brain-penetrant CDK4/6 vs CDK1/2 inhibitor for glioblastoma used in multiple replica molecular dynamics simulations. J Biomol Struct Dyn 2023:1-20. [PMID: 38112295 DOI: 10.1080/07391102.2023.2294175] [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: 07/27/2023] [Accepted: 11/23/2023] [Indexed: 12/21/2023]
Abstract
Cyclin dependent kinases (CDKs) play an important role in cell cycle regulation and their dysfunction is associated with many cancers. That is why CDKs have been attractive targets for the treatment of cancer. Glioblastoma is a cancer caused by the aberrant expression of CDK4/6, so exploring the mechanism of the selection of CDK4/6 toward inhibitors relative to the other family members CDK1/2 is essential. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), molecular mechanics Poisson-Boltzmann/Generalized Born surface area (MM-PB/GBSA) and other methods were integrated to decipher the selectively binding mechanism of the inhibitor N1J to CDK4/6 and CDK1/2. Molecular electrostatic potential (MESP) analysis provides an explanation for the N1J selectivity. Residue-based free energy decomposition reveals that most of the hot residues are located at the same location of CDKs proteins, but the different types of residues in different proteins cause changes in binding energy, which is considered as a potential developmental direction to improve the selectivity of inhibitors to CDK4/6. These results provide insights into the source of inhibitor and CDK4/6 selectivity for the future development of more selective inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Jiahao Sun
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shanshan Liang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
7
|
Drici N. The influence of the hydrogen-bond network on the structure and dynamics of the RAPRKKG heptapeptide and its mutants. J Mol Graph Model 2023; 125:108598. [PMID: 37586130 DOI: 10.1016/j.jmgm.2023.108598] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2023] [Revised: 08/02/2023] [Accepted: 08/08/2023] [Indexed: 08/18/2023]
Abstract
The structural behaviour of the RAPRKKG heptapeptide after individual or multiple mutations was inspected through molecular dynamics simulation. The nature of the mutations provided information on the flexibility of the heptapeptide and on how water molecules establish hydrogen bonds with it. The structural behaviour of the wild-type and the mutated structures were measured through the analysis of protein‒protein and protein‒solvent hydrogen bonds. The conformational behaviours of the different structures were analysed through free energy landscape analysis. The flexibility characteristics of the mutants seem to depend on the reorganization of water molecules and their static or dynamic behaviour around amino acid side chains.
Collapse
Affiliation(s)
- Nedjoua Drici
- University of Mostaganem, Abdelhamid Ibn Badis, Faculty of Exact Sciences and Informatics, Chemin des cretes ex INES, Mostaganem, 27000, Algeria; Laboratoire de Chimie Physique Macromoleculaire LCPM, University of Oran1 Ahmed benbella, Oran, 31000, Algeria.
| |
Collapse
|
8
|
Sun J, Liu X, Zhang S, Li M, Zhang Q, Chen J. Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study. Phys Chem Chem Phys 2023; 25:28479-28496. [PMID: 37846774 DOI: 10.1039/d3cp03392a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continues to spread globally, and rapid viral evolution and the emergence of new variants pose challenges to pandemic control. During infection, the spike protein of SARS-CoV-2 interacts with the human ACE2 protein via its receptor binding domain (RBD), and it is known that engineered forms of ACE2 can compete with wild-type (WT) ACE2 for binding to inhibit infection. Here, we conducted multiple replica molecular dynamics (MRMD) simulations to study the mechanisms of the engineered ACE2 variants 3N39 and 3N94 and provide directions for optimization. Our findings reveal that engineered ACE2 is notably more efficacious in systems that show weaker binding to WT ACE2 (i.e., WT and BA.1 RBD), but also faces immune escape as the virus evolves. Moreover, by modifying residue types near the binding interface, engineered ACE2 alters the electrostatic potential distribution and reconfigures the hydrogen bonding network, which results in modified binding to the RBD. However, this structural rearrangement does not occur in all RBD variants. In addition, we identified potentially engineerable beneficial residues and potentially engineerable detrimental residues in both ACE2 and RBD. Functional conservation can thus enable the optimization of these residues and improve the binding competitiveness of engineered ACE2, which therefore provides additional immune escape prevention. Finally, we conclude that these findings have implications for understanding the mechanisms responsible for engineered ACE2 and can help us to develop engineered ACE2 proteins that show superior performance.
Collapse
Affiliation(s)
- Jiahao Sun
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| |
Collapse
|
9
|
Yang F, Wang Y, Yan D, Liu Z, Wei B, Chen J, He W. Binding Mechanism of Inhibitors to Heat Shock Protein 90 Investigated by Multiple Independent Molecular Dynamics Simulations and Prediction of Binding Free Energy. Molecules 2023; 28:4792. [PMID: 37375347 DOI: 10.3390/molecules28124792] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Revised: 06/09/2023] [Accepted: 06/12/2023] [Indexed: 06/29/2023] Open
Abstract
The heat shock protein (HSP90) has been an import target of drug design in the treatment of human disease. An exploration of the conformational changes in HSP90 can provide useful information for the development of efficient inhibitors targeting HSP90. In this work, multiple independent all-atom molecular dynamics (AAMD) simulations followed by calculations of the molecular mechanics generalized Born surface area (MM-GBSA) were performed to explore the binding mechanism of three inhibitors (W8Y, W8V, and W8S) to HSP90. The dynamics analyses verified that the presence of inhibitors impacts the structural flexibility, correlated movements, and dynamics behavior of HSP90. The results of the MM-GBSA calculations suggest that the selection of GB models and empirical parameters has important influences on the predicted results and verify that van der Waals interactions are the main forces that determine inhibitor-HSP90 binding. The contributions of separate residues to the inhibitor-HSP90 binding process indicate that hydrogen-bonding interactions (HBIs) and hydrophobic interactions play important roles in HSP90-inhibitor identifications. Moreover, residues L34, N37, D40, A41, D79, I82, G83, M84, F124, and T171 are recognized as hot spots of inhibitor-HSP90 binding and provide significant target sites of for the design of drugs related to HSP90. This study aims to contribute to the development of efficient inhibitors that target HSP90 by providing an energy-based and theoretical foundation.
Collapse
Affiliation(s)
- Fen Yang
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
| | - Yiwen Wang
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Aeronautics, Shandong Jiaotong University, Jinan 250357, China
| | - Dongliang Yan
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Zhongtao Liu
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
| | - Benzheng Wei
- Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao 266112, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Weikai He
- School of Aeronautics, Shandong Jiaotong University, Jinan 250357, China
| |
Collapse
|
10
|
Wang Y, Yang F, Yan D, Zeng Y, Wei B, Chen J, He W. Identification Mechanism of BACE1 on Inhibitors Probed by Using Multiple Separate Molecular Dynamics Simulations and Comparative Calculations of Binding Free Energies. Molecules 2023; 28:4773. [PMID: 37375328 DOI: 10.3390/molecules28124773] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Revised: 06/12/2023] [Accepted: 06/12/2023] [Indexed: 06/29/2023] Open
Abstract
β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward the treatment of AD.
Collapse
Affiliation(s)
- Yiwen Wang
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Aeronautics, Shandong Jiaotong University, Jinan 250357, China
| | - Fen Yang
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
| | - Dongliang Yan
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Yalin Zeng
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
| | - Benzheng Wei
- Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao 266112, China
| | - Jianzhong Chen
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Weikai He
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Aeronautics, Shandong Jiaotong University, Jinan 250357, China
| |
Collapse
|
11
|
Mlotkowski AJ, Schlegel HB, Chow CS. Calculated p Ka Values for a Series of Aza- and Deaza-Modified Nucleobases. J Phys Chem A 2023; 127:3526-3534. [PMID: 37037184 PMCID: PMC10123669 DOI: 10.1021/acs.jpca.3c01358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
Abstract
A variety of synthetic modified nucleobases have been used to investigate the structure and function of RNA and DNA or act as enzyme inhibitors. A set of these modifications involves the addition or removal of a nitrogen atom in the ring. These aza and deaza modifications have garnered interest as useful biochemical tools, but information on some of their physical characteristics is lacking. In this study, the B3LYP density functional with the 6-31+G(d,p) basis set and an implicit-explicit solvent model was used to perform ab initio quantum mechanical studies to estimate pKa values of aza- and deaza-modified nucleobases. A comparison between theoretical and known experimental pKa values was carried out, and adjustment factors were applied to 57 pKa values in the purine and pyrimidine data sets.
Collapse
Affiliation(s)
- Alan J Mlotkowski
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - H Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Christine S Chow
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| |
Collapse
|
12
|
Chen L, Huang X, Li Y, Zhao B, Liang M, Wang R. Structural and energetic basis of interaction between human estrogen-related receptor γ and environmental endocrine disruptors from multiple molecular dynamics simulations and free energy predictions. JOURNAL OF HAZARDOUS MATERIALS 2023; 443:130174. [PMID: 36265380 DOI: 10.1016/j.jhazmat.2022.130174] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Revised: 10/04/2022] [Accepted: 10/09/2022] [Indexed: 06/16/2023]
Abstract
Environmental endocrine disruptors (EEDs), a class of molecules that are widespread in our environment, may adversely affect the endocrine system. Exploring the interactions between these compounds and their potential targets is important for assessing their role in the organism. Focused on the human estrogen-related receptor γ (hERRγ) with BPA, BPB, HPTE, BPE, BP(2,2)(Et), and BP(2,2)(MeO) complexes, respectively, we groped for the mechanisms of conformational changes and interactions of hERRγ when binding to these six EEDs by combining multiple molecular dynamics (MD) simulations with energy prediction (MM-PBSA and solvated interaction energy (SIE)). Dynamics analysis results revealed these six EEDs have different effects on the internal dynamics of hERRγ, resulting in significant changes in the interaction of key residues around Leu268, Val313, Leu345, and Phe435 with EEDs, and thus affected its binding energy with these EEDs. The energy calculations further demonstrated that van der Waals interactions are critical for these EEDs binding to hERRγ. These results present detailed molecular insight into the interaction features between EEDs and hERRγ and help guide the search for safer alternatives to BPA.
Collapse
Affiliation(s)
- Lin Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China.
| | - Xu Huang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China
| | - Yufei Li
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China
| | - Bing Zhao
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China; Heilongjiang Provincial Key Laboratory of Surface Active Agent and Auxiliary, Qiqihar 161006, PR China
| | - Min Liang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China
| | - Ruige Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar 161006, PR China
| |
Collapse
|
13
|
Wang L, Lu D, Wang Y, Xu X, Zhong P, Yang Z. Binding selectivity-dependent molecular mechanism of inhibitors towards CDK2 and CDK6 investigated by multiple short molecular dynamics and free energy landscapes. J Enzyme Inhib Med Chem 2023; 38:84-99. [DOI: 10.1080/14756366.2022.2135511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Affiliation(s)
- Lifei Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Dan Lu
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| | - Yan Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Xiaoyan Xu
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Peihua Zhong
- College of Computer Information and Engineering, Jiangxi Agriculture University, Nanchang, PR China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| |
Collapse
|
14
|
Wang R, Sun H, Chen W, Zhao B, Chen L. Molecular basis of ssDNA recognition by RBM45 protein of neurodegenerative disease from multiple molecular dynamics simulations and energy predictions. J Mol Graph Model 2023; 118:108377. [PMID: 36435031 DOI: 10.1016/j.jmgm.2022.108377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Revised: 11/09/2022] [Accepted: 11/14/2022] [Indexed: 11/21/2022]
Abstract
Neurodegenerative diseases (NDD) are a group of cognitive and behavioral disorders characterized by progressive loss of neuronal structure and function. As the population ages, the incidence is getting higher and higher, but there is currently no effective treatment. The details of RNA/DNA recognition by the RNA-binding protein RBM45 closely related to neurodegenerative diseases through its two tandem RNA-recognition domains at its N-terminus have important implications for structure-based drug discovery against degenerative diseases. To explore the key characteristics of interaction between ssDNA and RBM45, we performed multiple molecular dynamics (MD) simulations along with MM-PBSA energy prediction on the complexes of wild type (WT) and three mutant RBM45s (K100A, F124A/Y165A, and F29A/F70A/F124A/Y165A) with ssDNA, respectively. The findings suggest that these mutated residues of RBM45 modify the interaction of their surrounding residues with ssDNA, thereby affecting RBM45 protein binding to ssDNA. In contrast with WT RBM45 protein, variations in van der Waals and electrostatic interactions with ssDNA caused by these three RBM45 mutants are critical to affect binding between them. In addition, energy analysis showed that RBM45 is a specific ssDNA-binding protein. The results of our work provide valuable theoretical guidance for the design effective drugs of NDD.
Collapse
Affiliation(s)
- Ruige Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China
| | - Han Sun
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China
| | - Wei Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China.
| | - Bing Zhao
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China; Heilongjiang Provincial Key Laboratory of Surface Active Agent and Auxiliary, Qiqihar, 161006, PR China
| | - Lin Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, 161006, PR China.
| |
Collapse
|
15
|
Wang R, Zheng Q. Multiple Molecular Dynamics Simulations and Energy Analysis Unravel the Dynamic Properties and Binding Mechanism of Mutants HIV-1 Protease with DRV and CA-p2. Microbiol Spectr 2022; 10:e0074821. [PMID: 35319278 PMCID: PMC9045218 DOI: 10.1128/spectrum.00748-21] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Accepted: 02/10/2022] [Indexed: 11/20/2022] Open
Abstract
PRS17, a variant of human immunodeficiency virus type I protease (HIV-1 PR), has 17 mutated residues showing high levels of multidrug resistance. To describe the effects of these mutated residues on the dynamic properties and the binding mechanism of PR with substrate and inhibitor, focused on six systems (two complexes of WT PR and PRS17 with inhibitor Darunavir (DRV), two complexes of WT PR and PRS17 with substrate analogue CA-p2, two unligand WT PR and PRS17), we performed multiple molecular dynamics (MD) simulations combined with MM-PBSA and solvated interaction energy (SIE) methods. For both the unligand PRs and ligand-PR complexes, the results from simulations revealed 17 mutated residues alter the flap-flap distance, the distance from flap regions to catalytic sites, and the curling degree of the flap tips. These mutated residues changed the flexibility of the flap region in PR, and thus affected its binding energy with DRV and CA-p2, resulting in differences in sensitivity. Hydrophobic cavity makes an important contribution to the binding of PR and ligands. And most noticeable of all, the binding of the guanidine group in CA-p2 and Arg8' of PRS17 is useful for increasing their binding ability. These results have important guidance for the further design of drugs against multidrug resistant PR. IMPORTANCE Developing effective anti-HIV inhibitors is the current requirement to cope with the emergence of the resistance of mutants. Compared with the experiments, MD simulations along with energy calculations help reduce the time and cost of designing new inhibitors. Based on our simulation results, we propose two factors that may help design effective inhibitors against HIV-1 PR: (i) importance of hydrophobic cavity, and (ii) introduction of polar groups similar to the guanidine group.
Collapse
Affiliation(s)
- Ruige Wang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun, People's Republic of China
| | - Qingchuan Zheng
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun, People's Republic of China
| |
Collapse
|
16
|
Dutta N, Deb I, Sarzynska J, Lahiri A. Data-informed reparameterization of modified RNA and the effect of explicit water models: application to pseudouridine and derivatives. J Comput Aided Mol Des 2022; 36:205-224. [PMID: 35338419 PMCID: PMC8956458 DOI: 10.1007/s10822-022-00447-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Accepted: 03/04/2022] [Indexed: 11/29/2022]
Abstract
Pseudouridine is one of the most abundant post-transcriptional modifications in RNA. We have previously shown that the FF99-derived parameters for pseudouridine and some of its naturally occurring derivatives in the AMBER distribution either alone or in combination with the revised γ torsion parameters (parmbsc0) failed to reproduce their conformational characteristics observed experimentally (Deb et al. in J Chem Inf Model 54:1129–1142, 2014; Deb et al. in J Comput Chem 37:1576–1588, 2016; Dutta et al. in J Chem Inf Model 60:4995–5002, 2020). However, the application of the recommended bsc0 correction did lead to an improvement in the description not only of the distribution in the γ torsional space but also of the sugar pucker distributions. In an earlier study, we examined the transferability of the revised glycosidic torsion parameters (χIDRP) for Ψ to its derivatives. We noticed that although these parameters in combination with the AMBER FF99-derived parameters and the revised γ torsional parameters resulted in conformational properties of these residues that were in better agreement with experimental observations, the sugar pucker distributions were still not reproduced accurately. Here we report a new set of partial atomic charges for pseudouridine, 1-methylpseudouridine, 3-methylpseudouridine and 2′-O-methylpseudouridine and a new set of glycosidic torsional parameters (χND) based on chosen glycosidic torsional profiles that most closely corresponded to the NMR data for conformational propensities and studied their effect on the conformational distributions using REMD simulations at the individual nucleoside level. We have also studied the effect of the choice of water model on the conformational characteristics of these modified nucleosides. Our observations suggest that the current revised set of parameters and partial atomic charges describe the sugar pucker distributions for these residues more accurately and that the choice of a suitable water model is important for the accurate description of their conformational properties. We have further validated the revised sets of parameters by studying the effect of substitution of uridine with pseudouridine within single stranded RNA oligonucleotides on their conformational and hydration characteristics.
Collapse
Affiliation(s)
- Nivedita Dutta
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92, Acharya Prafulla Chandra Road, Kolkata, West Bengal, 700009, India
| | - Indrajit Deb
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92, Acharya Prafulla Chandra Road, Kolkata, West Bengal, 700009, India
| | - Joanna Sarzynska
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, Noskowskiego 12/14, 61-704, Poznan, Poland
| | - Ansuman Lahiri
- Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, 92, Acharya Prafulla Chandra Road, Kolkata, West Bengal, 700009, India.
| |
Collapse
|
17
|
Jones EL, Mlotkowski AJ, Hebert SP, Schlegel HB, Chow CS. Calculations of p Ka Values for a Series of Naturally Occurring Modified Nucleobases. J Phys Chem A 2022; 126:1518-1529. [PMID: 35201779 DOI: 10.1021/acs.jpca.1c10905] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Modified nucleobases are found in functionally important regions of RNA and are often responsible for essential structural roles. Many of these nucleobase modifications are dynamically regulated in nature, with each modification having a different biological role in RNA. Despite the high abundance of modifications, many of their characteristics are still poorly understood. One important property of a nucleobase is its pKa value, which has been widely studied for unmodified nucleobases, but not for the modified versions. In this study, the pKa values of modified nucleobases were determined by performing ab initio quantum mechanical calculations using a B3LYP density functional with the 6-31+G(d,p) basis set and a combination of implicit-explicit solvation systems. This method, which was previously employed to determine the pKa values of unmodified nucleobases, is applicable to a variety of modified nucleobases. Comparisons of the pKa values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.
Collapse
Affiliation(s)
- Evan L Jones
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Alan J Mlotkowski
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Sebastien P Hebert
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - H Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Christine S Chow
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| |
Collapse
|
18
|
Liang S, Liu X, Zhang S, Li M, Zhang Q, Chen J. Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations. Phys Chem Chem Phys 2022; 24:1743-1759. [PMID: 34985081 DOI: 10.1039/d1cp04361g] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (Mpro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, Mpro has been an attractive target of efficient drug design for anti-SARS-CoV-2 treatment. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), and the molecular mechanics-generalized Born surface area (MM-GBSA) method were integrated together to decipher the binding mechanism of four inhibitors masitinib, O6K, FJC and GQU to Mpro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of Mpro along with dihedral angle changes of key residues. The analysis of FELs unveils that the stability in the relative orientation and geometric position of inhibitors to Mpro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-Mpro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against Mpro. The hot spot residues including H41, M49, F140, N142, G143, C145, H163, H164, M165, E166 and Q189 identified by computational alanine scanning are considered as reliable targets of clinically available inhibitors inhibiting the activities of Mpro.
Collapse
Affiliation(s)
- Shanshan Liang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| |
Collapse
|
19
|
Chen J, Li N, Wang X, Chen J, Zhang JZH, Zhu T. Molecular mechanism related to the binding of fluorophores to Mango-II revealed by multiple-replica molecular dynamics simulations. Phys Chem Chem Phys 2021; 23:10636-10649. [PMID: 33904542 DOI: 10.1039/d0cp06438f] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Recently, RNA aptamers activating small-molecule fluorophores have been successfully applied to tag and track RNAs in vivo. It is of significance to investigate the molecular mechanism of the fluorophore-RNA aptamer bindings at the atomic level to seek a possible pathway to enhance the fluorescence efficiency of fluorophores. In this work, multiple replica molecular dynamics (MRMD) simulations, essential dynamics (ED) analysis, and hierarchical clustering analysis were coupled to probe the effect of A22U mutation on the binding of two fluorophores, TO1-Biotin (TO1) and TO3-Biotin (TO3), to the Mango-II RNA aptamer (Mango-II). ED analysis reveals that A22U induces alterations in the binding pocket and sites of TO1 and TO3 to the Mango-II, which in turn tunes the fluorophore-RNA interface and changes the interactions of TO1 and TO3 with separate nucleotides of Mango-II. Dynamics analyses also uncover that A22U exerts the opposite impact on the molecular surface areas of the Mango-II and sugar puckers of nucleotides 22 and 23 in Mango-II complexed with TO1 and TO3. Moreover, the calculations of binding free energies suggest that A22U strengthens the binding ability of TO1 to the mutated Mango-II but weakens TO3 to the mutated Mango-II when compared with WT. These findings imply that point mutation in nucleotides possibly tune the fluorescence of fluorophores binding to RNA aptamers, providing a possible scheme to enhance the fluorescence of fluorophores.
Collapse
Affiliation(s)
- Junxiao Chen
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and School of Chemistry and Pharmaceutical Engineering, Qilu University of Technology, Jinan, 250353, People's Republic of China
| | - Na Li
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China.
| | - Xingyu Wang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, People's Republic of China.
| | - John Z H Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China
| | - Tong Zhu
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China and Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| |
Collapse
|
20
|
Wang L, Wang Y, Yang Z, Xu S, Li H. Binding Selectivity of Inhibitors toward Bromodomains BAZ2A and BAZ2B Uncovered by Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations. ACS OMEGA 2021; 6:12036-12049. [PMID: 34056358 PMCID: PMC8154142 DOI: 10.1021/acsomega.1c00687] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Accepted: 04/19/2021] [Indexed: 06/12/2023]
Abstract
Two Bromodomain-Containing proteins BAZ2A and BAZ2B are responsible for remodeling chromatin and regulating noncoding RNAs. As for our current studies, integration of multiple short molecular dynamics simulations (MSMDSs) with molecular mechanics generalized Born surface area (MM-GBSA) method is adopted for insights into binding selectivity of three small molecules D8Q, D9T and UO1 to BAZ2A against BAZ2B. The calculations of MM-GBSA unveil that selectivity of inhibitors toward BAZ2A and BAZ2B highly depends on the enthalpy changes and the details uncover that D8Q has better selectivity toward BAZ2A than BAZ2B, D9T more favorably bind to BAZ2B than BAZ2A, and UO1 does not show obvious selectivity toward these two proteins. The analysis of interaction network between residues and inhibitors indicates that seven residues are mainly responsible for the selectivity of D8Q, six residues for D9T and four residues provide significant contributions to associations of UO1 with two proteins. Moreover the analysis of interaction network not only reveals warm spots of inhibitor bindings to BAZ2A and BAZ2B but also unveils that common residue pairs, including (W1816, W1887), (P1817, P1888), (F1818, F1889), (V1822, V1893), (N1823, N1894),(L1826, L1897), (V1827, V1898), (F1872, F1943), (N1873, N1944) and (V1879, I1950) belonging to (BAZ2A, BAZ2B), induce mainly binding differences of inhibitors to BAZ2A and BAZ2B. Hence, insights from our current studies offer useful dynamics information relating with conformational alterations and structure-affinity relationship at atomistic levels for novel therapeutic strategies toward BAZ2A and BAZ2B.
Collapse
Affiliation(s)
- Lifei Wang
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Yan Wang
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Zhiyong Yang
- Department
of Physics, Jiangxi Agricultural University, 1101 Zhimin Road, Economic and Technological
Development Zone, Nanchang, Jiangxi Province 330045, China
| | - Shuobo Xu
- School
of Information Science and Electrical Engineering, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Hongyun Li
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| |
Collapse
|
21
|
Yan FF, Gao F. Comparison of the binding characteristics of SARS-CoV and SARS-CoV-2 RBDs to ACE2 at different temperatures by MD simulations. Brief Bioinform 2021; 22:1122-1136. [PMID: 33611368 PMCID: PMC7929385 DOI: 10.1093/bib/bbab044] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2020] [Revised: 01/14/2021] [Accepted: 01/28/2021] [Indexed: 12/12/2022] Open
Abstract
Temperature plays a significant role in the survival and transmission of SARS-CoV (severe acute respiratory syndrome coronavirus) and SARS-CoV-2. To reveal the binding differences of SARS-CoV and SARS-CoV-2 receptor-binding domains (RBDs) to angiotensin-converting enzyme 2 (ACE2) at different temperatures at atomic level, 20 molecular dynamics simulations were carried out for SARS-CoV and SARS-CoV-2 RBD-ACE2 complexes at five selected temperatures, i.e. 200, 250, 273, 300 and 350 K. The analyses on structural flexibility and conformational distribution indicated that the structure of the SARS-CoV-2 RBD was more stable than that of the SARS-CoV RBD at all investigated temperatures. Then, molecular mechanics Poisson-Boltzmann surface area and solvated interaction energy approaches were combined to estimate the differences in binding affinity of SARS-CoV and SARS-CoV-2 RBDs to ACE2; it is found that the binding ability of ACE2 to the SARS-CoV-2 RBD was stronger than that to the SARS-CoV RBD at five temperatures, and the main reason for promoting such binding differences is electrostatic and polar interactions between RBDs and ACE2. Finally, the hotspot residues facilitating the binding of SARS-CoV and SARS-CoV-2 RBDs to ACE2, the key differential residues contributing to the difference in binding and the interaction mechanism of differential residues that exist at all investigated temperatures were analyzed and compared in depth. The current work would provide a molecular basis for better understanding of the high infectiousness of SARS-CoV-2 and offer better theoretical guidance for the design of inhibitors targeting infectious diseases caused by SARS-CoV-2.
Collapse
Affiliation(s)
- Fang-Fang Yan
- Department of Physics, School of Science, Tianjin University, Tianjin 300072, China
| | - Feng Gao
- Department of Physics, School of Science, Tianjin University, Tianjin 300072, China
- Frontiers Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin 300072, China
- SynBio Research Platform, Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin 300072, China
| |
Collapse
|
22
|
Wang Y, Wu S, Wang L, Yang Z, Zhao J, Zhang L. Binding selectivity of inhibitors toward the first over the second bromodomain of BRD4: theoretical insights from free energy calculations and multiple short molecular dynamics simulations. RSC Adv 2020; 11:745-759. [PMID: 35423696 PMCID: PMC8693360 DOI: 10.1039/d0ra09469b] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2020] [Accepted: 12/17/2020] [Indexed: 12/12/2022] Open
Abstract
Bromodomain-containing protein 4 (BRD4) plays an important role in mediating gene transcription involved in cancers and non-cancer diseases such as acute heart failure and inflammatory diseases. In this work, multiple short molecular dynamics (MSMD) simulations are integrated with a molecular mechanics generalized Born surface area (MM-GBSA) approach to decipher binding selectivity of three inhibitors 8NS, 82Y, and 837 toward two domains BD1 and BD2 of BRD4. The results demonstrate that the enthalpy effects play critical roles in selectivity identification of inhibitors toward BD1 and BD2, determining that 8NS has better selectivity toward BD2 than BD1, while 82Y and 837 more favorably bind to BD1 than BD2. A residue-based free-energy decomposition method was used to calculate an inhibitor-residue interaction spectrum and unveil contributions of separate residues to binding selectivity. The results identify six common residues, containing (P82, P375), (V87, V380), (L92, L385), (L94, L387), (N140, N433), and (I146, V439) individually belonging to (BD1, BD2) of BRD4, and yield a considerable binding difference of inhibitors to BD1 and BD2, suggesting that these residues play key roles in binding selectivity of inhibitors toward BD1 and BD2 of BRD4. Therefore, these results provide useful dynamics information and a structure affinity relationship for the development of highly selective inhibitors targeting BD1 and BD2 of BRD4.
Collapse
Affiliation(s)
- Yan Wang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Shiliang Wu
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Lifei Wang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University Nanchang 330045 China
| | - Juan Zhao
- School of Science, Shandong Jiaotong University Jinan 250357 China
| | - Lulu Zhang
- School of Science, Shandong Jiaotong University Jinan 250357 China
| |
Collapse
|
23
|
Wang R, Zheng Q. Multiple Molecular Dynamics Simulations of the Inhibitor GRL-02031 Complex with Wild Type and Mutant HIV-1 Protease Reveal the Binding and Drug-Resistance Mechanism. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:13817-13832. [PMID: 33175558 DOI: 10.1021/acs.langmuir.0c02151] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Human immunodeficiency virus type 1 (HIV-1) protease is regarded as a fascinating target for drug development against HIV infection. However, mutations causing drug resistance severely limit the efficiency of the recently marketed drugs in the treatment of HIV replication. To elucidate the binding mechanism of HIV-1 protease with promising inhibitor GRL-02031 and further to probe the resistance mechanism associated with mutations (I47V, L76V, V82A, and N88D) to the inhibitor, we applied multiple molecular dynamics (MMD) simulations along with energy analysis by the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methodology on specific HIV-1 protease with GRL-0231 complexes. On the basis of detail analysis of the simulations, we revealed key characteristics that constitute the drug resistance of four mutation HIV-1 proteases toward GRL-02031: substitution of the side chain in these four mutation residues leads to a change in the distances between the flaps and catalytic sites, thereby reducing the affinity for GRL-02031 with these four mutation proteases, even though the L76V and N88D residues cannot directly contact GRL-02031. The results of energy analysis according to the MM-PBSA and SIE methods further indicated that hydrophobic interaction was considered to be the prime driving force for inhibitor GRL-02031 binding to protease and the decrease in van der Waals interactions between inhibitor GRL-02031 and mutant proteases as the primary cause of the drug resistance. Analyses of the hydrogen bonds and atomic interactions further provided detailed explanations for the resistance of these four mutation proteases toward inhibitor GRL-02031. The present study provides potential guidance on the structure-based inhibitors' design targeting HIV-1 protease.
Collapse
Affiliation(s)
- Ruige Wang
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Qingchuan Zheng
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, P. R. China
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, College of Life Science, Jilin University, Changchun 130023, People's Republic of China
| |
Collapse
|
24
|
Xu J, Tojo S, Fujitsuka M, Kawai K. Dynamics of Single‐Stranded RNA Looping Probed and Photoregulated by Sulfonated Pyrene. ChemistrySelect 2020. [DOI: 10.1002/slct.202002231] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Jie Xu
- The Institute of Scientific and Industrial Research (SANKEN)Osaka University Mihogaoka 8–1 Ibaraki Osaka 567-0047 Japan
| | - Sachiko Tojo
- The Institute of Scientific and Industrial Research (SANKEN)Osaka University Mihogaoka 8–1 Ibaraki Osaka 567-0047 Japan
| | - Mamoru Fujitsuka
- The Institute of Scientific and Industrial Research (SANKEN)Osaka University Mihogaoka 8–1 Ibaraki Osaka 567-0047 Japan
| | - Kiyohiko Kawai
- The Institute of Scientific and Industrial Research (SANKEN)Osaka University Mihogaoka 8–1 Ibaraki Osaka 567-0047 Japan
| |
Collapse
|
25
|
Wang LF, Wang Y, Yang ZY, Zhao J, Sun HB, Wu SL. Revealing binding selectivity of inhibitors toward bromodomain-containing proteins 2 and 4 using multiple short molecular dynamics simulations and free energy analyses. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2020; 31:373-398. [PMID: 32496901 DOI: 10.1080/1062936x.2020.1748107] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Accepted: 03/24/2020] [Indexed: 06/11/2023]
Abstract
Emerging evidences indicate bromodomain-containing proteins 2 and 4 (BRD2 and BRD4) play critical roles in cancers, inflammations, cardiovascular diseases and other pathologies. Multiple short molecular dynamics (MSMD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were applied to investigate the binding selectivity of three inhibitors 87D, 88M and 89G towards BRD2 over BRD4. The root-mean-square fluctuation (RMSF) analysis indicates that the structural flexibility of BRD4 is stronger than that of BRD2. Moreover the calculated distances between the Cα atoms in the centres of the ZA_loop and BC_loop of BRD4 are also bigger than that of BRD2. The rank of binding free energies calculated using MM-GBSA method agrees well with that determined by experimental data. The results show that 87D can bind more favourably to BRD2 than BRD4, while 88M has better selectivity on BRD4 over BRD2. Residue-based free-energy decomposition method was utilized to estimate the inhibitor-residue interaction spectrum and the results not only identify the hot interaction spots of inhibitors with BRD2 and BRD4, but also demonstrate that several common residues, including (W370, W374), (P371, P375), (V376, V380) and (L381, L385) belonging to (BRD2, BRD4), generate significant binding difference of inhibitors to BRD2 and BRD4.
Collapse
Affiliation(s)
- L F Wang
- School of Science, Shandong Jiaotong University , Jinan, China
| | - Y Wang
- School of Science, Shandong Jiaotong University , Jinan, China
| | - Z Y Yang
- Department of Physics, Jiangxi Agricultural University , Nanchang, China
| | - J Zhao
- School of Science, Shandong Jiaotong University , Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University , Jinan, China
| | - S L Wu
- School of Science, Shandong Jiaotong University , Jinan, China
| |
Collapse
|
26
|
Chen J, Liu X, Zhang S, Chen J, Sun H, Zhang L, Zhang Q. Molecular mechanism with regard to the binding selectivity of inhibitors toward FABP5 and FABP7 explored by multiple short molecular dynamics simulations and free energy analyses. Phys Chem Chem Phys 2020; 22:2262-2275. [DOI: 10.1039/c9cp05704h] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Recently, fatty acid binding proteins 5 and 7 (FABP5 and FABP7) have been regarded as the prospective targets for clinically treating multiple diseases related to FABPs.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Xinguo Liu
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| | - Shaolong Zhang
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| | - Junxiao Chen
- School of Chemistry and Pharmaceutical Engineering
- Qilu University of Technology
- Jinan
- People's Republic of China
| | - Haibo Sun
- School of Science
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Lin Zhang
- School of Construction Machinery
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Qinggang Zhang
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| |
Collapse
|
27
|
Computational and NMR studies of RNA duplexes with an internal pseudouridine-adenosine base pair. Sci Rep 2019; 9:16278. [PMID: 31700156 PMCID: PMC6838189 DOI: 10.1038/s41598-019-52637-0] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2019] [Accepted: 10/03/2019] [Indexed: 01/28/2023] Open
Abstract
Pseudouridine (Ψ) is the most common chemical modification present in RNA. In general, Ψ increases the thermodynamic stability of RNA. However, the degree of stabilization depends on the sequence and structural context. To explain experimentally observed sequence dependence of the effect of Ψ on the thermodynamic stability of RNA duplexes, we investigated the structure, dynamics and hydration of RNA duplexes with an internal Ψ-A base pair in different nearest-neighbor sequence contexts. The structures of two RNA duplexes containing 5′-GΨC/3′-CAG and 5′-CΨG/3′-GAC motifs were determined using NMR spectroscopy. To gain insight into the effect of Ψ on duplex dynamics and hydration, we performed molecular dynamics (MD) simulations of RNA duplexes with 5′-GΨC/3′-CAG, 5′-CΨG/3′-GAC, 5′-AΨU/3′-UAA and 5′-UΨA/3′-AAU motifs and their unmodified counterparts. Our results showed a subtle impact from Ψ modification on the structure and dynamics of the RNA duplexes studied. The MD simulations confirmed the change in hydration pattern when U is replaced with Ψ. Quantum chemical calculations showed that the replacement of U with Ψ affected the intrinsic stacking energies at the base pair steps depending on the sequence context. The calculated intrinsic stacking energies help to explain the experimentally observed sequence dependent changes in the duplex stability from Ψ modification.
Collapse
|
28
|
Chen J, Wang J, Yin B, Pang L, Wang W, Zhu W. Molecular Mechanism of Binding Selectivity of Inhibitors toward BACE1 and BACE2 Revealed by Multiple Short Molecular Dynamics Simulations and Free-Energy Predictions. ACS Chem Neurosci 2019; 10:4303-4318. [PMID: 31545898 DOI: 10.1021/acschemneuro.9b00348] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The β-amyloid cleaving enzymes 1 and 2 (BACE1 and BACE2) have been regarded as the prospective targets for clinically treating Alzheimer's disease (AD) in the last two decades. Thus, insight into the binding differences of inhibitors to BACE1 and BACE2 is of significance for designing highly selective inhibitors toward the two proteins. In this work, multiple short molecular dynamics (MSMD) simulations are coupled with the molecular mechanics generalized Born surface area (MM-GBSA) method to probe the binding selectivity of three inhibitors DBO, CS9, and SC7 on BACE1 over BACE2. The results show that the entropy effect plays a key role in selectivity identification of inhibitors toward BACE1 and BACE2, which determines that DBO has better selectivity toward BACE2 over BACE1, while CS9 and CS7 can more favorably bind to BACE1 than BACE2. The hierarchical clustering analysis based on energetic contributions of residues suggests that BACE1 and BACE2 share the common hot interaction spots. The residue-based free-energy decomposition method was applied to compute the inhibitor-residue interaction spectrum, and the results recognize four common binding subpockets corresponding to the different groups of inhibitors, which can be used as efficient targets for designing highly selective inhibitors toward BACE1 and BACE2. Therefore, these results provide a useful molecular basis and dynamics information for development of highly selective inhibitors targeting BACE1 and BACE2.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Jinan Wang
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Baohua Yin
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan 250355, China
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Weiliang Zhu
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| |
Collapse
|
29
|
Krepl M, Vögele J, Kruse H, Duchardt-Ferner E, Wöhnert J, Sponer J. An intricate balance of hydrogen bonding, ion atmosphere and dynamics facilitates a seamless uracil to cytosine substitution in the U-turn of the neomycin-sensing riboswitch. Nucleic Acids Res 2018; 46:6528-6543. [PMID: 29893898 PMCID: PMC6061696 DOI: 10.1093/nar/gky490] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2018] [Accepted: 05/23/2018] [Indexed: 01/04/2023] Open
Abstract
The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.
Collapse
Affiliation(s)
- Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jennifer Vögele
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Elke Duchardt-Ferner
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Jens Wöhnert
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Jiri Sponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
30
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 336] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|
31
|
Bjerregaard N, Andreasen PA, Dupont DM. Expected and unexpected features of protein-binding RNA aptamers. WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 7:744-757. [PMID: 27173731 DOI: 10.1002/wrna.1360] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Received: 11/05/2015] [Revised: 04/03/2016] [Accepted: 04/04/2016] [Indexed: 12/28/2022]
Abstract
RNA molecules with high affinity to specific proteins can be isolated from libraries of up to 1016 different RNA sequences by systematic evolution of ligands by exponential enrichment (SELEX). These so-called protein-binding RNA aptamers are often interesting, e.g., as modulators of protein function for therapeutic use, for probing the conformations of proteins, for studies of basic aspects of nucleic acid-protein interactions, etc. Studies on the interactions between RNA aptamers and proteins display a number of expected and unexpected features, including the chemical nature of the interacting RNA-protein surfaces, the conformation of protein-bound aptamer versus free aptamer, the conformation of aptamer-bound protein versus free protein, and the effects of aptamers on protein function. Here, we review current insights into the details of RNA aptamer-protein interactions. WIREs RNA 2016, 7:744-757. doi: 10.1002/wrna.1360 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Nils Bjerregaard
- Department of Molecular Biology and Genetics, Aarhus University, Aarhus, Denmark
| | - Peter A Andreasen
- Department of Molecular Biology and Genetics, Aarhus University, Aarhus, Denmark
| | - Daniel M Dupont
- Department of Molecular Biology and Genetics, Aarhus University, Aarhus, Denmark.
| |
Collapse
|
32
|
Voltz K, Léonard J, Touceda PT, Conyard J, Chaker Z, Dejaegere A, Godet J, Mély Y, Haacke S, Stote RH. Quantitative sampling of conformational heterogeneity of a DNA hairpin using molecular dynamics simulations and ultrafast fluorescence spectroscopy. Nucleic Acids Res 2016; 44:3408-19. [PMID: 26896800 PMCID: PMC4838372 DOI: 10.1093/nar/gkw077] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2015] [Accepted: 01/29/2016] [Indexed: 01/22/2023] Open
Abstract
Molecular dynamics (MD) simulations and time resolved fluorescence (TRF) spectroscopy were combined to quantitatively describe the conformational landscape of the DNA primary binding sequence (PBS) of the HIV-1 genome, a short hairpin targeted by retroviral nucleocapsid proteins implicated in the viral reverse transcription. Three 2-aminopurine (2AP) labeled PBS constructs were studied. For each variant, the complete distribution of fluorescence lifetimes covering 5 orders of magnitude in timescale was measured and the populations of conformers experimentally observed to undergo static quenching were quantified. A binary quantification permitted the comparison of populations from experimental lifetime amplitudes to populations of aromatically stacked 2AP conformers obtained from simulation. Both populations agreed well, supporting the general assumption that quenching of 2AP fluorescence results from pi-stacking interactions with neighboring nucleobases and demonstrating the success of the proposed methodology for the combined analysis of TRF and MD data. Cluster analysis of the latter further identified predominant conformations that were consistent with the fluorescence decay times and amplitudes, providing a structure-based rationalization for the wide range of fluorescence lifetimes. Finally, the simulations provided evidence of local structural perturbations induced by 2AP. The approach presented is a general tool to investigate fine structural heterogeneity in nucleic acid and nucleoprotein assemblies.
Collapse
Affiliation(s)
- Karine Voltz
- Department of Integrative Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC), Institut National de Santé et de Recherche Médicale (INSERM) U964, Centre National de Recherche Scientifique (CNRS) UMR 7104, Université de Strasbourg, 67404 Illkirch, France
| | - Jérémie Léonard
- Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, CNRS UMR 7504, Université de Strasbourg, Strasbourg, France
| | - Patricia Tourón Touceda
- Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, CNRS UMR 7504, Université de Strasbourg, Strasbourg, France
| | - Jamie Conyard
- Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, CNRS UMR 7504, Université de Strasbourg, Strasbourg, France
| | - Ziyad Chaker
- Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, CNRS UMR 7504, Université de Strasbourg, Strasbourg, France
| | - Annick Dejaegere
- Department of Integrative Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC), Institut National de Santé et de Recherche Médicale (INSERM) U964, Centre National de Recherche Scientifique (CNRS) UMR 7104, Université de Strasbourg, 67404 Illkirch, France
| | - Julien Godet
- Laboratoire de Biophotonique et Pharmacologie, UMR 7213 CNRS, Faculté de Pharmacie, Université de Strasbourg, Illkirch Cedex, France
| | - Yves Mély
- Laboratoire de Biophotonique et Pharmacologie, UMR 7213 CNRS, Faculté de Pharmacie, Université de Strasbourg, Illkirch Cedex, France
| | - Stefan Haacke
- Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, CNRS UMR 7504, Université de Strasbourg, Strasbourg, France
| | - Roland H Stote
- Department of Integrative Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC), Institut National de Santé et de Recherche Médicale (INSERM) U964, Centre National de Recherche Scientifique (CNRS) UMR 7104, Université de Strasbourg, 67404 Illkirch, France
| |
Collapse
|
33
|
Auffinger P, Cheatham TE, Vaiana AC. Spontaneous Formation of KCl Aggregates in Biomolecular Simulations: A Force Field Issue? J Chem Theory Comput 2015; 3:1851-9. [PMID: 26627627 DOI: 10.1021/ct700143s] [Citation(s) in RCA: 152] [Impact Index Per Article: 16.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Realistic all-atom simulation of biological systems requires accurate modeling of both the biomolecules and their ionic environment. Recently, ion nucleation phenomena leading to the rapid growth of KCl or NaCl clusters in the vicinity of biomolecular systems have been reported. To better understand this phenomenon, molecular dynamics simulations of KCl aqueous solutions at three (1.0, 0.25, and 0.10 M) concentrations were performed. Two popular water models (TIP3P and SPC/E) and two Lennard-Jones parameter sets (AMBER and Dang) were combined to produce a total of 80 ns of molecular dynamics trajectories. Results suggest that the use of the Dang cation Lennard-Jones parameters instead of those adopted by the AMBER force-field produces a more accurate description of the ionic solution. In the later case, formation of salt aggregates is probably indicative of an artifact resulting from misbalanced force-field parameters. Because similar results were obtained with two different water parameter sets, the simulations exclude a water model dependency in the formation of anomalous ionic clusters. Overall, the results strongly suggest that for accurate modeling of ions in biomolecular systems, great care should be taken in choosing balanced ionic parameters even when using the most popular force-fields. These results invite a reexamination of older data obtained using available force-fields and a thorough check of the quality of current parameters sets by performing simulations at finite (>0.25 M) instead of minimal salt conditions.
Collapse
Affiliation(s)
- Pascal Auffinger
- Architecture et Réactivité de l'ARN, Université Louis Pasteur de Strasbourg, CNRS, IBMC, 15 rue René Descartes, 67084 Strasbourg, France, and Department of Medical Chemistry, Pharmaceutical Chemistry and Pharmaceutics and Bioengineering, University of Utah, Salt Lake City, Utah 84112
| | - Thomas E Cheatham
- Architecture et Réactivité de l'ARN, Université Louis Pasteur de Strasbourg, CNRS, IBMC, 15 rue René Descartes, 67084 Strasbourg, France, and Department of Medical Chemistry, Pharmaceutical Chemistry and Pharmaceutics and Bioengineering, University of Utah, Salt Lake City, Utah 84112
| | - Andrea C Vaiana
- Architecture et Réactivité de l'ARN, Université Louis Pasteur de Strasbourg, CNRS, IBMC, 15 rue René Descartes, 67084 Strasbourg, France, and Department of Medical Chemistry, Pharmaceutical Chemistry and Pharmaceutics and Bioengineering, University of Utah, Salt Lake City, Utah 84112
| |
Collapse
|
34
|
Yu H, Mu Y, Nordenskiöld L, Stock G. Influence of Nitroxide Spin Labels on RNA Structure: A Molecular Dynamics Simulation Study. J Chem Theory Comput 2015; 4:1781-7. [PMID: 26620180 DOI: 10.1021/ct800266e] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Pulsed electron double resonance (PELDOR) experiments on oligonucleotides provide a distance ruler that allows the measurement of nanometer distances accurately. The technique requires attachment of nitroxide spin labels to the nucleotides, which may possibly perturb its conformation. To study to what extent nitroxide spin labels may affect RNA structure, all-atom molecular dynamics simulations in explicit solvent are performed for six double-labeled RNA duplexes. A new parametrization of the force field for the nitroxide spin label is developed, which leads to intramolecular distances that are in good agreement with experimental results. Comparison of the results for spin-labeled and unlabeled RNA reveals that the conformational effect of the spin label depends significantly on whether the spin label is attached to the major or the minor groove of RNA. While major-groove spin labeling may to some extent affect the conformation of nearby base pairs, minor-groove spin labeling has the advantage of mostly preserving the RNA conformation.
Collapse
Affiliation(s)
- Hang Yu
- School of Biological Sciences, Nanyang Technological University, Singapore 637551, and Institute of Physical and Theoretical Chemistry, J. W. Goethe University, D-60438 Frankfurt, Germany
| | - Yuguang Mu
- School of Biological Sciences, Nanyang Technological University, Singapore 637551, and Institute of Physical and Theoretical Chemistry, J. W. Goethe University, D-60438 Frankfurt, Germany
| | - Lars Nordenskiöld
- School of Biological Sciences, Nanyang Technological University, Singapore 637551, and Institute of Physical and Theoretical Chemistry, J. W. Goethe University, D-60438 Frankfurt, Germany
| | - Gerhard Stock
- School of Biological Sciences, Nanyang Technological University, Singapore 637551, and Institute of Physical and Theoretical Chemistry, J. W. Goethe University, D-60438 Frankfurt, Germany
| |
Collapse
|
35
|
MD Simulations of tRNA and Aminoacyl-tRNA Synthetases: Dynamics, Folding, Binding, and Allostery. Int J Mol Sci 2015; 16:15872-902. [PMID: 26184179 PMCID: PMC4519929 DOI: 10.3390/ijms160715872] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2015] [Revised: 07/07/2015] [Accepted: 07/08/2015] [Indexed: 12/21/2022] Open
Abstract
While tRNA and aminoacyl-tRNA synthetases are classes of biomolecules that have been extensively studied for decades, the finer details of how they carry out their fundamental biological functions in protein synthesis remain a challenge. Recent molecular dynamics (MD) simulations are verifying experimental observations and providing new insight that cannot be addressed from experiments alone. Throughout the review, we briefly discuss important historical events to provide a context for how far the field has progressed over the past few decades. We then review the background of tRNA molecules, aminoacyl-tRNA synthetases, and current state of the art MD simulation techniques for those who may be unfamiliar with any of those fields. Recent MD simulations of tRNA dynamics and folding and of aminoacyl-tRNA synthetase dynamics and mechanistic characterizations are discussed. We highlight the recent successes and discuss how important questions can be addressed using current MD simulations techniques. We also outline several natural next steps for computational studies of AARS:tRNA complexes.
Collapse
|
36
|
Galindo-Murillo R, Roe DR, Cheatham TE. Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC). Biochim Biophys Acta Gen Subj 2014; 1850:1041-1058. [PMID: 25219455 DOI: 10.1016/j.bbagen.2014.09.007] [Citation(s) in RCA: 119] [Impact Index Per Article: 11.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2014] [Revised: 09/01/2014] [Accepted: 09/03/2014] [Indexed: 12/12/2022]
Abstract
BACKGROUND The structure and dynamics of DNA are critically related to its function. Molecular dynamics simulations augment experiment by providing detailed information about the atomic motions. However, to date the simulations have not been long enough for convergence of the dynamics and structural properties of DNA. METHODS Molecular dynamics simulations performed with AMBER using the ff99SB force field with the parmbsc0 modifications, including ensembles of independent simulations, were compared to long timescale molecular dynamics performed with the specialized Anton MD engine on the B-DNA structure d(GCACGAACGAACGAACGC). To assess convergence, the decay of the average RMSD values over longer and longer time intervals was evaluated in addition to assessing convergence of the dynamics via the Kullback-Leibler divergence of principal component projection histograms. RESULTS These molecular dynamics simulations-including one of the longest simulations of DNA published to date at ~44μs-surprisingly suggest that the structure and dynamics of the DNA helix, neglecting the terminal base pairs, are essentially fully converged on the ~1-5μs timescale. CONCLUSIONS We can now reproducibly converge the structure and dynamics of B-DNA helices, omitting the terminal base pairs, on the μs time scale with both the AMBER and CHARMM C36 nucleic acid force fields. Results from independent ensembles of simulations starting from different initial conditions, when aggregated, match the results from long timescale simulations on the specialized Anton MD engine. GENERAL SIGNIFICANCE With access to large-scale GPU resources or the specialized MD engine "Anton" it is possible for a variety of molecular systems to reproducibly and reliably converge the conformational ensemble of sampled structures. This article is part of a Special Issue entitled: Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, L.S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112, USA
| | - Daniel R Roe
- Department of Medicinal Chemistry, L.S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112, USA
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, L.S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112, USA.
| |
Collapse
|
37
|
Abstract
We investigate the role of water molecules in 89 protein–RNA complexes taken from the Protein Data Bank. Those with tRNA and single-stranded RNA are less hydrated than with duplex or ribosomal proteins. Protein–RNA interfaces are hydrated less than protein–DNA interfaces, but more than protein–protein interfaces. Majority of the waters at protein–RNA interfaces makes multiple H-bonds; however, a fraction do not make any. Those making H-bonds have preferences for the polar groups of RNA than its partner protein. The spatial distribution of waters makes interfaces with ribosomal proteins and single-stranded RNA relatively ‘dry’ than interfaces with tRNA and duplex RNA. In contrast to protein–DNA interfaces, mainly due to the presence of the 2′OH, the ribose in protein–RNA interfaces is hydrated more than the phosphate or the bases. The minor groove in protein–RNA interfaces is hydrated more than the major groove, while in protein–DNA interfaces it is reverse. The strands make the highest number of water-mediated H-bonds per unit interface area followed by the helices and the non-regular structures. The preserved waters at protein–RNA interfaces make higher number of H-bonds than the other waters. Preserved waters contribute toward the affinity in protein–RNA recognition and should be carefully treated while engineering protein–RNA interfaces.
Collapse
Affiliation(s)
- Amita Barik
- Department of Biotechnology, Indian Institute of Technology Kharagpur, Kharagpur-721302, India
| | - Ranjit Prasad Bahadur
- Department of Biotechnology, Indian Institute of Technology Kharagpur, Kharagpur-721302, India
| |
Collapse
|
38
|
Zhang X, Walker RC, Phizicky EM, Mathews DH. Influence of Sequence and Covalent Modifications on Yeast tRNA Dynamics. J Chem Theory Comput 2014; 10:3473-3483. [PMID: 25136272 PMCID: PMC4132867 DOI: 10.1021/ct500107y] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2014] [Indexed: 12/25/2022]
Abstract
![]()
Modified nucleotides are prevalent
in tRNA. Experimental studies
reveal that these covalent modifications play an important role in
tuning tRNA function. In this study, molecular dynamics (MD) simulations
were used to investigate how modifications alter tRNA dynamics. The
X-ray crystal structures of tRNA(Asp), tRNA(Phe), and tRNA(iMet),
both with and without modifications, were used as initial structures
for 333 ns explicit solvent MD simulations with AMBER. For each tRNA
molecule, three independent trajectory calculations were performed,
giving an aggregate of 6 μs of total MD across six molecules.
The global root-mean-square deviations (RMSD) of atomic positions
show that modifications only introduce significant rigidity to the
global structure of tRNA(Phe). Interestingly, RMSDs of the anticodon
stem-loop (ASL) suggest that modified tRNA has a more rigid structure
compared to the unmodified tRNA in this domain. The anticodon RMSDs
of the modified tRNAs, however, are higher than those of corresponding
unmodified tRNAs. These findings suggest that the rigidity of the
anticodon stem-loop is finely tuned by modifications, where rigidity
in the anticodon arm is essential for tRNA translocation in the ribosome,
and flexibility of the anticodon is important for codon recognition.
Sugar pucker and water residence time of pseudouridines in modified
tRNAs and corresponding uridines in unmodified tRNAs were assessed,
and the results reinforce that pseudouridine favors the 3′-endo
conformation and has a higher tendency to interact with water. Principal
component analysis (PCA) was used to examine correlated motions in
tRNA. Additionally, covariance overlaps of PCAs were compared for
trajectories of the same molecule and between trajectories of modified
and unmodified tRNAs. The comparison suggests that modifications alter
the correlated motions. For the anticodon bases, the extent of stacking
was compared between modified and unmodified molecules, and only unmodified
tRNA(Asp) has significantly higher percentage of stacking time. Overall,
the simulations reveal that the effect of covalent modification on
tRNA dynamics is not simple, with modifications increasing flexibility
in some regions of the structure and increasing rigidity in other
regions.
Collapse
Affiliation(s)
- Xiaoju Zhang
- Department of Biochemistry and Biophysics and Center for RNA Biology, University of Rochester Medical Center , Rochester, New York 14642, United States
| | - Ross C Walker
- San Diego Supercomputer Center, University of California San Diego , La Jolla, California 92093, United States ; Department of Chemistry and Biochemistry, University of California San Diego , La Jolla, California 92093, United States
| | - Eric M Phizicky
- Department of Biochemistry and Biophysics and Center for RNA Biology, University of Rochester Medical Center , Rochester, New York 14642, United States
| | - David H Mathews
- Department of Biochemistry and Biophysics and Center for RNA Biology, University of Rochester Medical Center , Rochester, New York 14642, United States
| |
Collapse
|
39
|
Yoon J, Lin JC, Hyeon C, Thirumalai D. Dynamical Transition and Heterogeneous Hydration Dynamics in RNA. J Phys Chem B 2014; 118:7910-9. [DOI: 10.1021/jp500643u] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Jeseong Yoon
- Korea Institute for Advanced Study, 130-722 Seoul, Korea
| | - Jong-Chin Lin
- Department
of Chemistry and Biochemistry, and Biophysics
Program, Institute for Physical Sciences and Technology, University of Maryland, College
Park, Maryland 20742, United States
| | | | - D. Thirumalai
- Department
of Chemistry and Biochemistry, and Biophysics
Program, Institute for Physical Sciences and Technology, University of Maryland, College
Park, Maryland 20742, United States
| |
Collapse
|
40
|
Negureanu L, Salsbury FR. Destabilization of the MutSα's protein-protein interface due to binding to the DNA adduct induced by anticancer agent carboplatin via molecular dynamics simulations. J Mol Model 2013; 19:4969-89. [PMID: 24061854 DOI: 10.1007/s00894-013-1998-2] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2013] [Accepted: 09/05/2013] [Indexed: 12/22/2022]
Abstract
DNA mismatch repair (MMR) proteins maintain genetic integrity in all organisms by recognizing and repairing DNA errors. Such alteration of hereditary information can lead to various diseases, including cancer. Besides their role in DNA repair, MMR proteins detect and initiate cellular responses to certain type of DNA damage. Its response to the damaged DNA has made the human MMR pathway a useful target for anticancer agents such as carboplatin. This study indicates that strong, specific interactions at the interface of MutSα in response to the mismatched DNA recognition are replaced by weak, non-specific interactions in response to the damaged DNA recognition. Data suggest a severe impairment of the dimerization of MutSα in response to the damaged DNA recognition. While the core of MutSα is preserved in response to the damaged DNA recognition, the loss of contact surface and the rearrangement of contacts at the protein interface suggest a different packing in response to the damaged DNA recognition. Coupled in response to the mismatched DNA recognition, interaction energies, hydrogen bonds, salt bridges, and solvent accessible surface areas at the interface of MutSα and within the subunits are uncoupled or asynchronously coupled in response to the damaged DNA recognition. These pieces of evidence suggest that the loss of a synchronous mode of response in the MutSα's surveillance for DNA errors would possibly be one of the mechanism(s) of signaling the MMR-dependent programed cell death much wanted in anticancer therapies. The analysis was drawn from dynamics simulations.
Collapse
|
41
|
Negureanu L, Salsbury FR. Non-specificity and synergy at the binding site of the carboplatin-induced DNA adduct via molecular dynamics simulations of the MutSα-DNA recognition complex. J Biomol Struct Dyn 2013; 32:969-92. [PMID: 23799640 DOI: 10.1080/07391102.2013.799437] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
MutSα is the most abundant mismatch-binding factor of human DNA mismatch repair (MMR) proteins. MMR maintains genetic stability by recognizing and repairing DNA defects. Failure to accomplish their function may lead to cancer. In addition, MutSα recognizes at least some types of DNA damage making it a target for anticancer agents. Here, complementing scarce experimental data, we report unique hydrogen-bonding motifs associated with the recognition of the carboplatin induced DNA damage by MutSα. These data predict that carboplatin and cisplatin induced damaging DNA adducts are recognized by MutSα in a similar manner. Our simulations also indicate that loss of base pairing at the damage site results in (1) non-specific binding and (2) changes in the atomic flexibility at the lesion site and beyond. To further quantify alterations at MutSα-DNA interface in response to damage recognition, non-bonding interactions and salt bridges were investigated. These data indicate (1) possible different packing and (2) disruption of the salt bridges at the MutSα-DNA interface in the damaged complex. These findings (1) underscore the general observation of disruptions at the MutSα-DNA interface and (2) highlight the nature of the anticancer effect of the carboplatin agent. The analysis was carried out from atomistic simulations.
Collapse
|
42
|
Jiang J, Sakakibara Y, Chow CS. Helix 69: A Multitasking RNA Motif as a Novel Drug Target. Isr J Chem 2013. [DOI: 10.1002/ijch.201300012] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
43
|
Suga K, Tanabe T, Umakoshi H. Heterogeneous cationic liposomes modified with 3β-{N-[(N',N'-dimethylamino)ethyl]carbamoyl}cholesterol can induce partial conformational changes in messenger RNA and regulate translation in an Escherichia coli cell-free translation system. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2013; 29:1899-1907. [PMID: 23323854 DOI: 10.1021/la3050576] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
The effect of cationic liposomes (CLs) on messenger RNA(mRNA) conformation and translation was studied, focusing on membrane heterogeneity. CLs, composed of 1,2-dioleoyl-sn-glycerol-3-phosphocholine/1,2-dioleoyl-3-timethylammonium propane (DOPC/DOTAP) and DOPC/3β-{N-[(N',N'-dimethylamino)ethyl]carbamoyl}cholesterol (DOPC/DC-Ch), inhibited mRNA translation in an Escherichia coli cell-free translation system. Analysis of the membrane fluidity and polarity indicated a heterogeneous DOPC/DC-Ch (70/30) membrane, while other CLs exhibited homogeneous disordered membranes. mRNA adsorbed onto DOPC/DC-Ch liposomes showed translational activity, while DOPC/DOTAP liposomes inhibited mRNA translation in proportion to its adsorption onto membranes. Dehydration of DOPC/DOTAP (70/30) and DOPC/DC-Ch (70/30) was observed in the presence of mRNA but not in the case of zwitterionic DOPC liposomes, indicating that mRNA binds in regions between the phosphate [-PO(2)(-)-] and carbonyl [-C=O-] moieties of lipids. UV resonance Raman spectroscopy suggests that adenine, cytosine, and guanine interact with DOPC/DOTAP (70/30) and DOPC/DC-Ch (70/30) but not with DOPC. Circular dichroism indicates that DOPC/DOTAP (70/30) extensively denatured the mRNA. In contrast, heterogeneous DOPC/DC-Ch (70/30) induced partial conformational changes but maintained the translational activity of mRNA.
Collapse
Affiliation(s)
- Keishi Suga
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama-cho, Toyonaka, Osaka 560-8531, Japan
| | | | | |
Collapse
|
44
|
Golebiowski J, Topin J, Charlier L, Briand L. Interaction between odorants and proteins involved in the perception of smell: the case of odorant-binding proteins probed by molecular modelling and biophysical data. FLAVOUR FRAG J 2012. [DOI: 10.1002/ffj.3121] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Jérôme Golebiowski
- Institut de Chimie de Nice, UMR7272 CNRS; Université de Nice Sophia Antipolis; 06108; Nice; France
| | - Jérémie Topin
- Institut de Chimie de Nice, UMR7272 CNRS; Université de Nice Sophia Antipolis; 06108; Nice; France
| | - Landry Charlier
- Institut de Chimie de Nice, UMR7272 CNRS; Université de Nice Sophia Antipolis; 06108; Nice; France
| | - Loïc Briand
- Centre des Sciences du Goût et de l'Alimentation, INRA UMR1324, CNRS UMR6265; Université de Bourgogne; 21000; Dijon; France
| |
Collapse
|
45
|
Giegé R, Jühling F, Pütz J, Stadler P, Sauter C, Florentz C. Structure of transfer RNAs: similarity and variability. WILEY INTERDISCIPLINARY REVIEWS-RNA 2011; 3:37-61. [DOI: 10.1002/wrna.103] [Citation(s) in RCA: 112] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
46
|
Sumita M, Jiang J, SantaLucia J, Chow CS. Comparison of solution conformations and stabilities of modified helix 69 rRNA analogs from bacteria and human. Biopolymers 2011; 97:94-106. [PMID: 21858779 DOI: 10.1002/bip.21706] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2011] [Accepted: 07/31/2011] [Indexed: 12/20/2022]
Abstract
The helix 69 (H69) region of the large subunit (28S) ribosomal RNA (rRNA) of Homo sapiens contains five pseudouridine (Ψ) residues out of 19 total nucleotides, three of which are highly conserved. In this study, the effects of this abundant modified nucleotide on the structure and stability of H69 were compared with those of uridine in double-stranded (stem) regions. These results were compared with previous hairpin (stem plus single-stranded loop) studies to understand the contributions of the loop sequences to H69 structure and stability. The role of a loop nucleotide substitution from an A in bacteria (position 1918 in Escherichia coli 23S rRNA) to a G in eukaryotes (position 3734 in H. sapiens 28S rRNA) was examined. Thermodynamic parameters for the duplex RNAs were obtained through UV melting studies, and differences in the modified and unmodified RNA structures were examined by circular dichroism spectroscopy. The overall folded structure of human H69 appears to be similar to the bacterial RNA, consistent with the idea that ribosome structure and function are highly conserved; however, our results reveal subtle differences in structure and stability between the bacterial and human H69 RNAs in both the stem and loop regions. These findings may be significant with respect to H69 as a potential drug target site.
Collapse
Affiliation(s)
- Minako Sumita
- Department of Chemistry, Wayne State University, Detroit, MI 48202, USA
| | | | | | | |
Collapse
|
47
|
Denning EJ, Priyakumar UD, Nilsson L, MacKerell AD. Impact of 2'-hydroxyl sampling on the conformational properties of RNA: update of the CHARMM all-atom additive force field for RNA. J Comput Chem 2011; 32:1929-43. [PMID: 21469161 PMCID: PMC3082605 DOI: 10.1002/jcc.21777] [Citation(s) in RCA: 288] [Impact Index Per Article: 22.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2010] [Revised: 01/24/2011] [Accepted: 01/30/2011] [Indexed: 01/02/2023]
Abstract
Here, we present an update of the CHARMM27 all-atom additive force field for nucleic acids that improves the treatment of RNA molecules. The original CHARMM27 force field parameters exhibit enhanced Watson-Crick base pair opening which is not consistent with experiment, whereas analysis of molecular dynamics (MD) simulations show the 2'-hydroxyl moiety to almost exclusively sample the O3' orientation. Quantum mechanical (QM) studies of RNA related model compounds indicate the energy minimum associated with the O3' orientation to be too favorable, consistent with the MD results. Optimization of the dihedral parameters dictating the energy of the 2'-hydroxyl proton targeting the QM data yielded several parameter sets, which sample both the base and O3' orientations of the 2'-hydroxyl to varying degrees. Selection of the final dihedral parameters was based on reproduction of hydration behavior as related to a survey of crystallographic data and better agreement with experimental NMR J-coupling values. Application of the model, designated CHARMM36, to a collection of canonical and noncanonical RNA molecules reveals overall improved agreement with a range of experimental observables as compared to CHARMM27. The results also indicate the sensitivity of the conformational heterogeneity of RNA to the orientation of the 2'-hydroxyl moiety and support a model whereby the 2'-hydroxyl can enhance the probability of conformational transitions in RNA.
Collapse
Affiliation(s)
- Elizabeth J. Denning
- Department of Pharmaceutical Sciences, School of Pharmacy, University
of Maryland, Baltimore, MD 21201
| | - U. Deva Priyakumar
- Department of Pharmaceutical Sciences, School of Pharmacy, University
of Maryland, Baltimore, MD 21201
| | - Lennart Nilsson
- Department of Pharmaceutical Sciences, School of Pharmacy, University
of Maryland, Baltimore, MD 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University
of Maryland, Baltimore, MD 21201
| |
Collapse
|
48
|
Wu G, Yu AT, Kantartzis A, Yu YT. Functions and mechanisms of spliceosomal small nuclear RNA pseudouridylation. WILEY INTERDISCIPLINARY REVIEWS. RNA 2011; 2:571-81. [PMID: 21957045 PMCID: PMC4161978 DOI: 10.1002/wrna.77] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Pseudouridines are the most abundant and highly conserved modified nucleotides identified in spliceosomal small nuclear RNAs (snRNAs). Most pseudouridines are also clustered in functionally important regions of spliceosomal snRNAs. Experiments carried out in several independent experimental systems show that the pseudouridines in spliceosomal snRNAs are functionally important for pre-messenger RNA (mRNA) splicing. Experimental data also indicate that spliceosomal snRNA pseudouridylation can be catalyzed by both RNA-dependent (box H/ACA Ribonucleoproteins) and RNA-independent (protein-only enzymes) mechanisms.
Collapse
Affiliation(s)
- Guowei Wu
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, NY, USA
| | - Andrew T. Yu
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, NY, USA
| | - Athena Kantartzis
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, NY, USA
| | - Yi-Tao Yu
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, Rochester, NY, USA
| |
Collapse
|
49
|
Blaby IK, Majumder M, Chatterjee K, Jana S, Grosjean H, de Crécy-Lagard V, Gupta R. Pseudouridine formation in archaeal RNAs: The case of Haloferax volcanii. RNA (NEW YORK, N.Y.) 2011; 17:1367-80. [PMID: 21628430 PMCID: PMC3138572 DOI: 10.1261/rna.2712811] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Pseudouridine (Ψ), the isomer of uridine, is commonly found at various positions of noncoding RNAs of all organisms. Ψ residues are formed by a number of single- or multisite specific Ψ synthases, which generally act as stand-alone proteins. In addition, in Eukarya and Archaea, specific ribonucleoprotein complexes, each containing a distinct box H/ACA guide RNA and four core proteins, can produce Ψ at many sites of different cellular RNAs. Cbf5 is the core Ψ synthase in these complexes. Using Haloferax volcanii as an archaeal model organism, we show that, contrary to eukaryotes, the Cbf5 homolog (HVO_2493) is not essential in this archaeon. The Cbf5-deleted strain of H. volcanii completely lacks Ψ at positions 1940, 1942, 2605, and 2591 (Escherichia coli positions 1915, 1917, 2572, and 2586) of its 23S rRNA, and contains reduced steady-state levels of some box H/ACA RNAs. Archaeal Cbf5 is known to have tRNA Ψ55 synthase activity in vitro but we could not confirm this activity in vivo in H. volcanii. Conversely, the Pus10 (previously PsuX) homolog (HVO_1979), which can produce tRNA Ψ55, as well as Ψ54 in vitro, is shown here to be essential in H. volcanii, whereas the corresponding tRNA Ψ55 synthases, Pus4 and TruB, are not essential in yeast and E. coli, respectively. Finally, we demonstrate that HVO_1852, the TruA/Pus3 homolog, is responsible for the pseudouridylation of position 39 in H. volcanii tRNAs and that the corresponding gene is not essential.
Collapse
Affiliation(s)
- Ian K. Blaby
- Department of Microbiology & Cell Science, University of Florida, Gainesville, Florida 32611-0700, USA
| | - Mrinmoyee Majumder
- Department of Biochemistry and Molecular Biology, Southern Illinois University, Carbondale, Illinois 62901-4413, USA
| | - Kunal Chatterjee
- Department of Biochemistry and Molecular Biology, Southern Illinois University, Carbondale, Illinois 62901-4413, USA
| | - Sujata Jana
- Department of Biochemistry and Molecular Biology, Southern Illinois University, Carbondale, Illinois 62901-4413, USA
| | - Henri Grosjean
- Université Paris 11, IGM, CNRS, UMR 8621, Orsay, F 91405, France
| | - Valérie de Crécy-Lagard
- Department of Microbiology & Cell Science, University of Florida, Gainesville, Florida 32611-0700, USA
- Corresponding authors.E-mail .E-mail .
| | - Ramesh Gupta
- Department of Biochemistry and Molecular Biology, Southern Illinois University, Carbondale, Illinois 62901-4413, USA
- Corresponding authors.E-mail .E-mail .
| |
Collapse
|
50
|
Karijolich J, Yu YT. Converting nonsense codons into sense codons by targeted pseudouridylation. Nature 2011; 474:395-8. [PMID: 21677757 PMCID: PMC3381908 DOI: 10.1038/nature10165] [Citation(s) in RCA: 270] [Impact Index Per Article: 20.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2011] [Accepted: 04/28/2011] [Indexed: 11/09/2022]
Abstract
All three translation termination codons, or nonsense codons, contain a uridine residue at the first position of the codon. Here, we demonstrate that pseudouridylation (conversion of uridine into pseudouridine (Ψ), ref. 4) of nonsense codons suppresses translation termination both in vitro and in vivo. In vivo targeting of nonsense codons is accomplished by the expression of an H/ACA RNA capable of directing the isomerization of uridine to Ψ within the nonsense codon. Thus, targeted pseudouridylation represents a novel approach for promoting nonsense suppression in vivo. Remarkably, we also show that pseudouridylated nonsense codons code for amino acids with similar properties. Specifically, ΨAA and ΨAG code for serine and threonine, whereas ΨGA codes for tyrosine and phenylalanine, thus suggesting a new mode of decoding. Our results also suggest that RNA modification, as a naturally occurring mechanism, may offer a new way to expand the genetic code.
Collapse
Affiliation(s)
- John Karijolich
- Department of Biochemistry and Biophysics, University of Rochester Medical Center, 601 Elmwood Avenue, Rochester, New York 14642, USA
| | | |
Collapse
|