1
|
Zhang Z, Mlýnský V, Krepl M, Šponer J, Stadlbauer P. Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations. J Chem Inf Model 2024; 64:3896-3911. [PMID: 38630447 PMCID: PMC11094737 DOI: 10.1021/acs.jcim.4c00227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2024] [Revised: 04/02/2024] [Accepted: 04/02/2024] [Indexed: 05/14/2024]
Abstract
Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.
Collapse
Affiliation(s)
- Zhengyue Zhang
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic
- CEITEC−Central
European Institute of Technology, Masaryk
University, Kamenice
5, Brno 625 00, Czech Republic
- National
Center for Biomolecular Research,
Faculty of Science, Masaryk University, Kamenice 5, Brno 625 00, Czech Republic
| | - Vojtěch Mlýnský
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic
| | - Miroslav Krepl
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic
| | - Petr Stadlbauer
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, Brno 61200, Czech Republic
| |
Collapse
|
2
|
Zhao S, Cieplak P, Duan Y, Luo R. Assessment of Amino Acid Electrostatic Parametrizations of the Polarizable Gaussian Multipole Model. J Chem Theory Comput 2024; 20:2098-2110. [PMID: 38394331 PMCID: PMC11060985 DOI: 10.1021/acs.jctc.3c01347] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2024]
Abstract
Accurate parametrization of amino acids is pivotal for the development of reliable force fields for molecular modeling of biomolecules such as proteins. This study aims to assess amino acid electrostatic parametrizations with the polarizable Gaussian Multipole (pGM) model by evaluating the performance of the pGM-perm (with atomic permanent dipoles) and pGM-ind (without atomic permanent dipoles) variants compared to the traditional RESP model. The 100-conf-combterm fitting strategy on tetrapeptides was adopted, in which (1) all peptide bond atoms (-CO-NH-) share identical set of parameters and (2) the total charges of the two terminal N-acetyl (ACE) and N-methylamide (NME) groups were set to neutral. The accuracy and transferability of electrostatic parameters across peptides with varying lengths and real-world examples were examined. The results demonstrate the enhanced performance of the pGM-perm model in accurately representing the electrostatic properties of amino acids. This insight underscores the potential of the pGM-perm model and the 100-conf-combterm strategy for the future development of the pGM force field.
Collapse
Affiliation(s)
- Shiji Zhao
- Nurix Therapeutics, Inc., 1700 Owens St. Suite 205, San Francisco, CA 94158, USA
| | - Piotr Cieplak
- SBP Medical Discovery Institute, 10901 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Yong Duan
- UC Davis Genome Center and Department of Biomedical Engineering, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Ray Luo
- Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine. Irvine, California 92697, United States
| |
Collapse
|
3
|
DeLuca M, Sensale S, Lin PA, Arya G. Prediction and Control in DNA Nanotechnology. ACS APPLIED BIO MATERIALS 2024; 7:626-645. [PMID: 36880799 DOI: 10.1021/acsabm.2c01045] [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] [Indexed: 03/08/2023]
Abstract
DNA nanotechnology is a rapidly developing field that uses DNA as a building material for nanoscale structures. Key to the field's development has been the ability to accurately describe the behavior of DNA nanostructures using simulations and other modeling techniques. In this Review, we present various aspects of prediction and control in DNA nanotechnology, including the various scales of molecular simulation, statistical mechanics, kinetic modeling, continuum mechanics, and other prediction methods. We also address the current uses of artificial intelligence and machine learning in DNA nanotechnology. We discuss how experiments and modeling are synergistically combined to provide control over device behavior, allowing scientists to design molecular structures and dynamic devices with confidence that they will function as intended. Finally, we identify processes and scenarios where DNA nanotechnology lacks sufficient prediction ability and suggest possible solutions to these weak areas.
Collapse
Affiliation(s)
- Marcello DeLuca
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Sebastian Sensale
- Department of Physics, Cleveland State University, Cleveland, Ohio 44115, United States
| | - Po-An Lin
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Gaurav Arya
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
4
|
Polêto MD, Lemkul JA. Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes. J Chem Inf Model 2023; 63:6851-6862. [PMID: 37847037 PMCID: PMC10841373 DOI: 10.1021/acs.jcim.3c01305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication and maintaining the genome, among others; thus, changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and telomeric repeat-containing RNA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on the GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in the electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle but important differences between DNA and RNA GQs with the same fold.
Collapse
Affiliation(s)
- Marcelo D Polêto
- Department of Biochemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
5
|
Gber TE, Louis H, Ngana OC, Amodu IO, Ekereke EE, Benjamin I, Adalikwu SA, Adeyinka A. Yttrium- and zirconium-decorated Mg 12O 12-X (X = Y, Zr) nanoclusters as sensors for diazomethane (CH 2N 2) gas. RSC Adv 2023; 13:25391-25407. [PMID: 37636506 PMCID: PMC10448449 DOI: 10.1039/d3ra02939e] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Accepted: 07/21/2023] [Indexed: 08/29/2023] Open
Abstract
Diazomethane (CH2N2) presents a notable hazard as a respiratory irritant, resulting in various adverse effects upon exposure. Consequently, there has been increasing concern in the field of environmental research to develop a sensor material that exhibits heightened sensitivity and conductivity for the detection and adsorption of this gas. Therefore, this study aims to provide a comprehensive analysis of the geometric structure of three systems: CH2N2@MgO (C1), CH2N2@YMgO (CY1), and CH2N2@ZrMgO (CZ1), in addition to pristine MgO nanocages. The investigation involves a theoretical analysis employing the DFT/ωB97XD method at the GenECP/6-311++G(d,p)/SDD level of theory. Notably, the examination of bond lengths within the MgO cage yielded specific values, including Mg15-O4 (1.896 Å), Mg19-O4 (1.952 Å), and Mg23-O4 (1.952 Å), thereby offering valuable insights into the structural properties and interactions with CH2N2 gas. Intriguingly, after the interaction, bond length variations were observed, with CH2N2@MgO exhibiting shorter bonds and CH2N2@YMgO showcasing longer bonds. Meanwhile, CH2N2@ZrMgO displayed shorter bonds, except for a longer bond in Mg19-O4, suggesting increased stability due to shorter bond distances. The study further investigated the electronic properties, revealing changes in the energy gap that influenced electrical conductivity and sensitivity. The energy gap increased for Zr@MgO, CH2N2@MgO, CH2N2@YMgO, and CH2N2@ZrMgO, indicating weak interactions on the MgO surface. Conversely, Y@MgO showed a decrease in energy, suggesting a strong interaction. The pure MgO surface exhibited the ability to donate and accept electrons, resulting in an energy gap of 4.799 eV. Surfaces decorated with yttrium and zirconium exhibited decreased energies of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO), as well as decreased energy gap, indicating increased conductivity and sensitivity. Notably, Zr@MgO had the highest energy gap before CH2N2 adsorption, but C1 exhibited a significantly higher energy gap after adsorption, implying increased conductivity and sensitivity. The study also examined the density of states, demonstrating significant variations in the electronic properties of MgO and its decorated surfaces due to CH2N2 adsorption. Moreover, various analysis techniques were employed, including natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM), and noncovalent interaction (NCI) analysis, which provided insights into bonding, charge density, and intermolecular interactions. The findings contribute to a deeper understanding of the sensing mechanisms of CH2N2 gas on nanocage surfaces, shedding light on adsorption energy, conductivity, and recovery time. These results hold significance for gas-sensing applications and provide a basis for further exploration and development in this field.
Collapse
Affiliation(s)
- Terkumbur E Gber
- Computational and Bio-Simulation Research Group, University of Calabar Calabar Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar Calabar Nigeria
| | - Hitler Louis
- Computational and Bio-Simulation Research Group, University of Calabar Calabar Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar Calabar Nigeria
- Centre for Herbal Pharmacology and Environmental Sustainability, Chettinad Hospital and Research Institute, Chettinad Academy of Research and Education Kelambakkam-603103 Tamil Nadu India
| | - Obinna C Ngana
- Department of Chemical Sciences, Federal University of Wukari Wukari Taraba State Nigeria
| | - Ismail O Amodu
- Computational and Bio-Simulation Research Group, University of Calabar Calabar Nigeria
- Department of Mathematics, Faculty of Physical Sciences, University of Calabar Calabar Nigeria
| | - Ernest E Ekereke
- Computational and Bio-Simulation Research Group, University of Calabar Calabar Nigeria
- Department of Mathematics, Faculty of Physical Sciences, University of Calabar Calabar Nigeria
| | - Innocent Benjamin
- Computational and Bio-Simulation Research Group, University of Calabar Calabar Nigeria
| | - Stephen A Adalikwu
- Computational and Bio-Simulation Research Group, University of Calabar Calabar Nigeria
| | - Adedapo Adeyinka
- Department of Chemical Sciences, University of Johannesburg South Africa
| |
Collapse
|
6
|
Polêto MD, Lemkul JA. Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.10.552810. [PMID: 37645825 PMCID: PMC10461924 DOI: 10.1101/2023.08.10.552810] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication; and maintaining the genome, among others, thus changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and TERRA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle, but important, differences between DNA and RNA GQs with the same fold.
Collapse
Affiliation(s)
- Marcelo D Polêto
- Department of Biochemistry, Virginia Tech, Blacksburg, VA 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, VA 24061, United States
- Center for Drug Discovery, Virginia Tech, Blacksburg, VA 24061, United States
| |
Collapse
|
7
|
Zhang Z, Šponer J, Bussi G, Mlýnský V, Šulc P, Simmons CR, Stephanopoulos N, Krepl M. Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations. J Chem Inf Model 2023; 63:2794-2809. [PMID: 37126365 PMCID: PMC10170514 DOI: 10.1021/acs.jcim.3c00358] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.
Collapse
Affiliation(s)
- Zhengyue Zhang
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- CEITEC─Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
- National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Petr Šulc
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Chad R Simmons
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Nicholas Stephanopoulos
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 241/27, 783 71 Olomouc, Czech Republic
| |
Collapse
|
8
|
Zhao S, Cieplak P, Duan Y, Luo R. Transferability of the Electrostatic Parameters of the Polarizable Gaussian Multipole Model. J Chem Theory Comput 2023; 19:924-941. [PMID: 36696564 PMCID: PMC10152989 DOI: 10.1021/acs.jctc.2c01048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Accuracy and transferability are the two highly desirable properties of molecular mechanical force fields. Compared with the extensively used point-charge additive force fields that apply fixed atom-centered point partial charges to model electrostatic interactions, polarizable force fields are thought to have the advantage of modeling the atomic polarization effects. Previous works have demonstrated the accuracy of the recently developed polarizable Gaussian multipole (pGM) models. In this work, we assessed the transferability of the electrostatic parameters of the pGM models with (pGM-perm) and without (pGM-ind) atomic permanent dipoles in terms of reproducing the electrostatic potentials surrounding molecules/oligomers absent from electrostatic parameterizations. Encouragingly, both the pGM-perm and pGM-ind models show significantly improved transferability than the additive model in the tests (1) from water monomer to water oligomer clusters; (2) across different conformations of amino acid dipeptides and tetrapeptides; (3) from amino acid tetrapeptides to longer polypeptides; and (4) from nucleobase monomers to Watson-Crick base pair dimers and tetramers. Furthermore, we demonstrated that the double-conformation fittings using amino acid tetrapeptides in the αR and β conformations can result in good transferability not only across different tetrapeptide conformations but also from tetrapeptides to polypeptides with lengths ranging from 1 to 20 repetitive residues for both the pGM-ind and pGM-perm models. In addition, the observation that the pGM-ind model has significantly better accuracy and transferability than the point-charge additive model, even though they have an identical number of parameters, strongly suggest the importance of intramolecular polarization effects. In summary, this and previous works together show that the pGM models possess both accuracy and transferability, which are expected to serve as foundations for the development of next-generation polarizable force fields for modeling various polarization-sensitive biological systems and processes.
Collapse
Affiliation(s)
- Shiji Zhao
- Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, Irvine, California 92697, United States
| | - Piotr Cieplak
- SBP Medical Discovery Institute, 10901 North Torrey Pines Road, La Jolla, California 92037, United States
| | - Yong Duan
- UC Davis Genome Center and Department of Biomedical Engineering, University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Ray Luo
- Departments of Molecular Biology and Biochemistry, Chemical and Biomolecular Engineering, Materials Science and Engineering, and Biomedical Engineering, University of California, Irvine, Irvine, California 92697, United States
| |
Collapse
|
9
|
Paloncýová M, Pykal M, Kührová P, Banáš P, Šponer J, Otyepka M. Computer Aided Development of Nucleic Acid Applications in Nanotechnologies. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2204408. [PMID: 36216589 DOI: 10.1002/smll.202204408] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2022] [Revised: 09/12/2022] [Indexed: 06/16/2023]
Abstract
Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.
Collapse
Affiliation(s)
- Markéta Paloncýová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Martin Pykal
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Petra Kührová
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Pavel Banáš
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
| | - Jiří Šponer
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, Brno, 612 65, Czech Republic
| | - Michal Otyepka
- Regional Center of Advanced Technologies and Materials, The Czech Advanced Technology and Research Institute (CATRIN), Palacký University Olomouc, Šlechtitelů 27, Olomouc, 779 00, Czech Republic
- IT4Innovations, VŠB - Technical University of Ostrava, 17. listopadu 2172/15, Ostrava-Poruba, 708 00, Czech Republic
| |
Collapse
|
10
|
Salsbury A, Michel HM, Lemkul JA. Ion-Dependent Conformational Plasticity of Telomeric G-Hairpins and G-Quadruplexes. ACS OMEGA 2022; 7:23368-23379. [PMID: 35847338 PMCID: PMC9280957 DOI: 10.1021/acsomega.2c01600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Telomeric DNA is guanine-rich and can adopt structures such as G-quadruplexes (GQs) and G-hairpins. Telomeric GQs influence genome stability and telomerase activity, making understanding of enzyme-GQ interactions and dynamics important for potential drug design. GQs have a characteristic tetrad core, which is connected by loop regions. Within this architecture are G-hairpins, fold-back motifs that are thought to represent the first intermediate in GQ folding. To better understand the relationship between G-hairpin motifs and GQs, we performed polarizable simulations of a two-tetrad telomeric GQ and an isolated SC11 telomeric G-hairpin. The telomeric GQ contains a G-triad, which functions as part of the tetrad core or linker regions, depending on local conformational change. This triad and another motif below the tetrad core frequently bound ions and may represent druggable sites. Further, we observed the unbiased formation of a G-triad and a G-tetrad in simulations of the SC11 G-hairpin and found that cations can be partially hydrated while facilitating the formation of these motifs. Finally, we demonstrated that K+ ions form specific interactions with guanine bases, while Na+ ions interact nonspecifically with bases in the structure. Together, these simulations provide new insights into the influence of ions on GQs, G-hairpins, and G-triad motifs.
Collapse
Affiliation(s)
- Alexa
M. Salsbury
- Department
of Biochemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Haley M. Michel
- Department
of Biochemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Justin A. Lemkul
- Department
of Biochemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Center
for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
11
|
Tucker MR, Piana S, Tan D, LeVine MV, Shaw DE. Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA-Protein Complexes. J Phys Chem B 2022; 126:4442-4457. [PMID: 35694853 PMCID: PMC9234960 DOI: 10.1021/acs.jpcb.1c10971] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
![]()
Although molecular
dynamics (MD) simulations have been used extensively
to study the structural dynamics of proteins, the role of MD simulation
in studies of nucleic acid based systems has been more limited. One
contributing factor to this disparity is the historically lower level
of accuracy of the physical models used in such simulations to describe
interactions involving nucleic acids. By modifying nonbonded and torsion
parameters of a force field from the Amber family of models, we recently
developed force field parameters for RNA that achieve a level of accuracy
comparable to that of state-of-the-art protein force fields. Here
we report force field parameters for DNA, which we developed by transferring
nonbonded parameters from our recently reported RNA force field and
making subsequent adjustments to torsion parameters. We have also
modified the backbone charges in both the RNA and DNA parameter sets
to make the treatment of electrostatics compatible with our recently
developed variant of the Amber protein and ion force field. We name
the force field resulting from the union of these three parameter
sets (the new DNA parameters, the revised RNA parameters, and the
existing protein and ion parameters) DES-Amber. Extensive
testing of DES-Amber indicates that it can describe the thermal stability
and conformational flexibility of single- and double-stranded DNA
systems with a level of accuracy comparable to or, especially for
disordered systems, exceeding that of state-of-the-art nucleic acid
force fields. Finally, we show that, in certain favorable cases, DES-Amber
can be used for long-timescale simulations of protein–nucleic
acid complexes.
Collapse
Affiliation(s)
| | - Stefano Piana
- D. E. Shaw Research, New York, New York 10036, United States
| | - Dazhi Tan
- D. E. Shaw Research, New York, New York 10036, United States
| | | | - David E Shaw
- D. E. Shaw Research, New York, New York 10036, United States.,Department of Biochemistry and Molecular Biophysics, Columbia University, New York, New York 10032, United States
| |
Collapse
|
12
|
Castelli M, Doria F, Freccero M, Colombo G, Moroni E. Studying the Dynamics of a Complex G-Quadruplex System: Insights into the Comparison of MD and NMR Data. J Chem Theory Comput 2022; 18:4515-4528. [PMID: 35666124 PMCID: PMC9281369 DOI: 10.1021/acs.jctc.2c00291] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Molecular dynamics
(MD) simulations are coming of age in the study
of nucleic acids, including specific tertiary structures such as G-quadruplexes.
While being precious for providing structural and dynamic information
inaccessible to experiments at the atomistic level of resolution,
MD simulations in this field may still be limited by several factors.
These include the force fields used, different models for ion parameters,
ionic strengths, and water models. We address various aspects of this
problem by analyzing and comparing microsecond-long atomistic simulations
of the G-quadruplex structure formed by the human immunodeficiency
virus long terminal repeat (HIV LTR)-III sequence for which nuclear
magnetic resonance (NMR) structures are available. The system is studied
in different conditions, systematically varying the ionic strengths,
ion numbers, and water models. We comparatively analyze the dynamic
behavior of the G-quadruplex motif in various conditions and assess
the ability of each simulation to satisfy the nuclear magnetic resonance
(NMR)-derived experimental constraints and structural parameters.
The conditions taking into account K+-ions to neutralize
the system charge, mimicking the intracellular ionic strength, and
using the four-atom water model are found to be the best in reproducing
the experimental NMR constraints and data. Our analysis also reveals
that in all of the simulated environments residues belonging to the
duplex moiety of HIV LTR-III exhibit the highest flexibility.
Collapse
Affiliation(s)
- Matteo Castelli
- Department of Chemistry, University of Pavia, V.le Taramelli 12, 27100 Pavia, Italy
| | - Filippo Doria
- Department of Chemistry, University of Pavia, V.le Taramelli 12, 27100 Pavia, Italy
| | - Mauro Freccero
- Department of Chemistry, University of Pavia, V.le Taramelli 12, 27100 Pavia, Italy
| | - Giorgio Colombo
- Department of Chemistry, University of Pavia, V.le Taramelli 12, 27100 Pavia, Italy.,Institute of Chemical Sciences and Technologies SCITEC-CNR, Via Mario Bianco, 9, 20131 Milano, Italy
| | - Elisabetta Moroni
- Institute of Chemical Sciences and Technologies SCITEC-CNR, Via Mario Bianco, 9, 20131 Milano, Italy
| |
Collapse
|
13
|
Sengul MY, MacKerell AD. Accurate Modeling of RNA Hairpins Through the Explicit Treatment of Electronic Polarizability with the Classical Drude Oscillator Force Field. JOURNAL OF COMPUTATIONAL BIOPHYSICS AND CHEMISTRY 2022; 21:461-471. [DOI: 10.1142/s2737416521420060] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Molecular dynamics (MD) simulations play a crucial role in modeling biomolecular systems in which the electrostatic interactions are critical in dictating the structural and dynamical properties. Thus, the treatment of the electrostatic interactions defined in the underlying force field (FF) strongly affects the simulation accuracy. Most FFs use fixed partial atomic charges to include electrostatic interactions, and therefore lack the electronic polarization response, representing an intrinsic limitation. To address this limitation, polarizable FFs have been developed that treat atomic polarizabilities explicitly. Here we present the application of the all-atom polarizable (Drude) and non-polarizable (CHARMM) nucleic acid FFs in RNA hairpin systems to investigate the impact of polarization on structural properties, dipole moment distributions, and cation interactions. Results show that the presence of polarizability in the FF significantly improves the stabilization of RNA hairpin structure. As expected, the distributions of dipole moments show more fluctuations when simulated using the polarizable FF, with the variation in dipoles contributing to the stabilization of the structures of the loop regions of the RNAs. Contact map analyses between the bases and cations show that the variation of the ion distribution around the entire hairpin is larger for the polarizable FF and the cations occupy the outer hydration shell to a greater extent. The presented results indicate the importance of the explicit treatment of electronic polarizability in molecular simulations of RNA, including in non-canonical regions.
Collapse
Affiliation(s)
- Mert Y. Sengul
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, USA
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, USA
| |
Collapse
|
14
|
Green AT, Pickard AJ, Li R, MacKerell AD, Bierbach U, Cho SS. Computational and Experimental Characterization of rDNA and rRNA G-Quadruplexes. J Phys Chem B 2022; 126:609-619. [PMID: 35026949 DOI: 10.1021/acs.jpcb.1c08340] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
DNA G-quadruplexes in human telomeres and gene promoters are being extensively studied for their role in controlling the growth of cancer cells. G-quadruplexes have been unambiguously shown to exist both in vitro and in vivo, including in the guanine (G)-rich DNA genes encoding pre-ribosomal RNA (pre-rRNA), which is transcribed in the cell's nucleolus. Recent studies strongly suggest that these DNA sequences ("rDNA"), and the transcribed rRNA, are a potential anticancer target through the inhibition of RNA polymerase I (Pol I) in ribosome biogenesis, but the structures of ribosomal G-quadruplexes at atomic resolution are unknown and very little biophysical characterization has been performed on them to date. In the present study, circular dichroism (CD) spectroscopy is used to show that two putative rDNA G-quadruplex sequences, NUC 19P and NUC 23P and their counterpart rRNAs, predominantly adopt parallel topologies, reminiscent of the analogous telomeric quadruplex structures. Based on this information, we modeled parallel topology atomistic structures of the putative ribosomal G-quadruplexes. We then validated and refined the modeled ribosomal G-quadruplex structures using all-atom molecular dynamics (MD) simulations with the CHARMM36 force field in the presence and absence of stabilizing K+. Motivated by preliminary MD simulations of the telomeric parallel G-quadruplex (TEL 24P) in which the K+ ion is expelled, we used updated CHARMM36 force field K+ parameters that were optimized, targeting the data from quantum mechanical calculations and the polarizable Drude model force field. In subsequent MD simulations with optimized CHARMM36 parameters, the K+ ions are predominantly in the G-quadruplex channel and the rDNA G-quadruplexes have more well-defined, predominantly parallel-topology structures as compared to rRNA. In addition, NUC 19P is more structured than NUC 23P, which contains extended loops. Results from this study set the structural foundation for understanding G-quadruplex functions and the design of novel chemotherapeutics against these nucleolar targets and can be readily extended to other DNA and RNA G-quadruplexes.
Collapse
Affiliation(s)
- Adam T Green
- Department of Physics, Wake Forest University, Winston-Salem, North Carolina 27109, United States
| | - Amanda J Pickard
- Department of Chemistry, Wake Forest University, Wake Downtown Campus, Winston-Salem, North Carolina 27101, United States
| | - Rongzhong Li
- Department of Physics, Wake Forest University, Winston-Salem, North Carolina 27109, United States.,Department of Computer Science, Wake Forest University, Winston-Salem, North Carolina 27109, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, United States
| | - Ulrich Bierbach
- Department of Chemistry, Wake Forest University, Wake Downtown Campus, Winston-Salem, North Carolina 27101, United States
| | - Samuel S Cho
- Department of Physics, Wake Forest University, Winston-Salem, North Carolina 27109, United States.,Department of Computer Science, Wake Forest University, Winston-Salem, North Carolina 27109, United States
| |
Collapse
|
15
|
Gber TE, Louis H, Owen AE, Etinwa BE, Benjamin I, Asogwa FC, Orosun MM, Eno EA. Heteroatoms (Si, B, N, and P) doped 2D monolayer MoS 2 for NH 3 gas detection. RSC Adv 2022; 12:25992-26010. [PMID: 36199611 PMCID: PMC9468912 DOI: 10.1039/d2ra04028j] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Accepted: 08/29/2022] [Indexed: 11/21/2022] Open
Abstract
2D transition metal dichalcogenide MoS2 monolayer quantum dots (MoS2-QD) and their doped boron (B@MoS2-QD), nitrogen (N@MoS2-QD), phosphorus (P@MoS2-QD), and silicon (Si@MoS2-QD) surfaces have been theoretically investigated using density functional theory (DFT) computation to understand their mechanistic sensing ability, such as conductivity, selectivity, and sensitivity toward NH3 gas. The results from electronic properties showed that P@MoS2-QD had the lowest energy gap, which indicated an increase in electrical conductivity and better adsorption behavior. By carrying out comparative adsorption studies using m062-X, ωB97XD, B3LYP, and PBE0 methods at the 6-311G++(d,p) level of theory, the most negative values were observed from ωB97XD for the P@MoS2-QD surface, signifying the preferred chemisorption surface for NH3 detection. The mechanistic studies provided in this study also indicate that the P@MoS2-QD dopant is a promising sensing material for monitoring ammonia gas in the real world. We hope this research work will provide informative knowledge for experimental researchers to realize the potential of MoS2 dopants, specifically the P@MoS2-QD surface, as a promising candidate for sensors to detect gas. 2D transition metal dichalcogenide MoS2 monolayer quantum dots (MoS2-QD) and their doped boron (B@MoS2-QD), nitrogen (N@MoS2-QD), phosphorus (P@MoS2-QD), and silicon (Si@MoS2-QD) counterparts are proposed as selective sensors for NH3 gas.![]()
Collapse
Affiliation(s)
- Terkumbur E. Gber
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar, Calabar, Nigeria
| | - Hitler Louis
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
- Department of Chemistry, Akwa-Ibom State University, Uyo, Nigeria
| | - Aniekan E. Owen
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
- Department of Chemistry, Akwa-Ibom State University, Uyo, Nigeria
| | - Benjamin E. Etinwa
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar, Calabar, Nigeria
| | - Innocent Benjamin
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar, Calabar, Nigeria
| | - Fredrick C. Asogwa
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar, Calabar, Nigeria
| | | | - Ededet A. Eno
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
- Department of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar, Calabar, Nigeria
| |
Collapse
|
16
|
Milovanović B, Petković M, Popov I, Etinski M. Water-Mediated Interactions Enhance Alkaline Earth Cation Chelation in Neighboring Cavities of a Cytosine Quartet in the DNA Quadruplex. J Phys Chem B 2021; 125:11996-12005. [PMID: 34694801 DOI: 10.1021/acs.jpcb.1c05598] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Larger Coulombic repulsion between divalent cations compared to the monovalent counterparts dictates the cation-cation distance in the central ion channel of quadruplexes. In this work, density functional theory and a continuum solvation model were employed to study bond energies of alkaline earth cations in adjacent cavities of the central ion channel. Four crystallized tetramolecular quadruplexes with various geometric constraints and structural motifs available in the Protein Data Bank were examined in order to understand how the cation binding affinities could be increased in aqueous solution. A cytosine quartet sandwiched between guanine quartets has a larger bond energy of the second alkaline earth cation in comparison with guanine and uracil quartets. Four highly conserved hydrogen-bonded water molecules in the center of the cytosine quartet are responsible for a higher electrostatic interaction with the cations in comparison with guanines' carbonyl groups. The reported findings are valuable for the design of synthetic quadruplexes templated with divalent cations for optoelectronic applications.
Collapse
Affiliation(s)
- Branislav Milovanović
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16 11000 Belgrade, Serbia
| | - Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16 11000 Belgrade, Serbia
| | - Igor Popov
- Institute for Multidisciplinary Research, University of Belgrade, Kneza Višeslava 1, 11030 Belgrade, Serbia.,Institut of Physics, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16 11000 Belgrade, Serbia
| |
Collapse
|
17
|
Li N, Gao Y, Qiu F, Zhu T. Benchmark Force Fields for the Molecular Dynamic Simulation of G-Quadruplexes. Molecules 2021; 26:5379. [PMID: 34500812 PMCID: PMC8434458 DOI: 10.3390/molecules26175379] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2021] [Revised: 08/30/2021] [Accepted: 08/31/2021] [Indexed: 11/29/2022] Open
Abstract
G-quadruplexes have drawn widespread attention for serving as a potential anti-cancer target and their application in material science. Molecular dynamics (MD) simulation is the key theoretical tool in the study of GQ's structure-function relationship. In this article, we systematically benchmarked the five force fields of parmbsc0, parmbsc1, OL15, AMOEBA, and Drude2017 on the MD simulation of G-quadruplex from four aspects: structural stability, central ion channel stability, description of Hoogsteen hydrogen bond network, and description of the main chain dihedral angle. The results show that the overall performance of the Drude force field is the best. Although there may be a certain over-polarization effect, it is still the best choice for the MD simulation of G-quadruplexes.
Collapse
Affiliation(s)
- Na Li
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200241, China;
| | - Ya Gao
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai 201620, China
| | - Feng Qiu
- Institute of Artificial Intelligence on Education, Shanghai Normal University, Shanghai 200234, China
| | - Tong Zhu
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200241, China;
- NYU-ECNU Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
- Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| |
Collapse
|
18
|
Sánchez-González Á, Bandeira NAG, Ortiz de Luzuriaga I, Martins FF, Elleuchi S, Jarraya K, Lanuza J, Lopez X, Calhorda MJ, Gil A. New Insights on the Interaction of Phenanthroline Based Ligands and Metal Complexes and Polyoxometalates with Duplex DNA and G-Quadruplexes. Molecules 2021; 26:4737. [PMID: 34443326 PMCID: PMC8397986 DOI: 10.3390/molecules26164737] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Revised: 07/28/2021] [Accepted: 07/29/2021] [Indexed: 12/24/2022] Open
Abstract
This work provides new insights from our team regarding advances in targeting canonical and non-canonical nucleic acid structures. This modality of medical treatment is used as a form of molecular medicine specifically against the growth of cancer cells. Nevertheless, because of increasing concerns about bacterial antibiotic resistance, this medical strategy is also being explored in this field. Up to three strategies for the use of DNA as target have been studied in our research lines during the last few years: (1) the intercalation of phenanthroline derivatives with duplex DNA; (2) the interaction of metal complexes containing phenanthroline with G-quadruplexes; and (3) the activity of Mo polyoxometalates and other Mo-oxo species as artificial phosphoesterases to catalyze the hydrolysis of phosphoester bonds in DNA. We demonstrate some promising computational results concerning the favorable interaction of these small molecules with DNA that could correspond to cytotoxic effects against tumoral cells and microorganisms. Therefore, our results open the door for the pharmaceutical and medical applications of the compounds we propose.
Collapse
Affiliation(s)
- Ángel Sánchez-González
- Faculdade de Ciências, BioISI-Biosystems and Integrative Sciences Institute, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal; (Á.S.-G.); (N.A.G.B.); (F.F.M.); (M.J.C.)
| | - Nuno A. G. Bandeira
- Faculdade de Ciências, BioISI-Biosystems and Integrative Sciences Institute, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal; (Á.S.-G.); (N.A.G.B.); (F.F.M.); (M.J.C.)
| | - Iker Ortiz de Luzuriaga
- CIC nanoGUNE BRTA, Tolosa Hiribidea 76, Euskadi, 20018 Donostia-San Sebastián, Spain;
- Polimero eta Material Aurreratuak, Fisika, Kimika eta Teknologia Saila, Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), Paseo Manuel de Lardizabal 3, 20018 Donostia-San Sebastián, Spain; (J.L.); (X.L.)
| | - Frederico F. Martins
- Faculdade de Ciências, BioISI-Biosystems and Integrative Sciences Institute, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal; (Á.S.-G.); (N.A.G.B.); (F.F.M.); (M.J.C.)
| | - Sawssen Elleuchi
- Laboratoire de Chimie Inorganique, LR17ES07, Faculté de Sciences de Sfax, Université de Sfax, Sfax 3000, Tunisia; (S.E.); (K.J.)
| | - Khaled Jarraya
- Laboratoire de Chimie Inorganique, LR17ES07, Faculté de Sciences de Sfax, Université de Sfax, Sfax 3000, Tunisia; (S.E.); (K.J.)
| | - Jose Lanuza
- Polimero eta Material Aurreratuak, Fisika, Kimika eta Teknologia Saila, Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), Paseo Manuel de Lardizabal 3, 20018 Donostia-San Sebastián, Spain; (J.L.); (X.L.)
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Xabier Lopez
- Polimero eta Material Aurreratuak, Fisika, Kimika eta Teknologia Saila, Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), Paseo Manuel de Lardizabal 3, 20018 Donostia-San Sebastián, Spain; (J.L.); (X.L.)
- Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 4, 20018 Donostia-San Sebastián, Spain
| | - Maria José Calhorda
- Faculdade de Ciências, BioISI-Biosystems and Integrative Sciences Institute, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal; (Á.S.-G.); (N.A.G.B.); (F.F.M.); (M.J.C.)
| | - Adrià Gil
- Faculdade de Ciências, BioISI-Biosystems and Integrative Sciences Institute, Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal; (Á.S.-G.); (N.A.G.B.); (F.F.M.); (M.J.C.)
- CIC nanoGUNE BRTA, Tolosa Hiribidea 76, Euskadi, 20018 Donostia-San Sebastián, Spain;
| |
Collapse
|
19
|
Salsbury AM, Lemkul JA. Cation competition and recruitment around the c-kit1 G-quadruplex using polarizable simulations. Biophys J 2021; 120:2249-2261. [PMID: 33794153 PMCID: PMC8390831 DOI: 10.1016/j.bpj.2021.03.022] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Revised: 02/22/2021] [Accepted: 03/25/2021] [Indexed: 11/24/2022] Open
Abstract
Nucleic acid-ion interactions are fundamentally important to the physical, energetic, and conformational properties of DNA and RNA. These interactions help fold and stabilize highly ordered secondary and tertiary structures, such as G-quadruplexes (GQs), which are functionally relevant in telomeres, replication initiation sites, and promoter sequences. The c-kit proto-oncogene encodes for a receptor tyrosine kinase and is linked to gastrointestinal stromal tumors, mast cell disease, and leukemia. This gene contains three unique GQ-forming sequences that have proposed antagonistic effects on gene expression. The dominant GQ, denoted c-kit1, has been shown to decrease expression of c-kit transcripts, making the c-kit1 GQ a promising drug target. Toward disease intervention, more information is needed regarding its conformational dynamics and ion binding properties. Therefore, we performed molecular dynamics simulations of the c-kit1 GQ with K+, Na+, Li+, and mixed salt solutions using the Drude-2017 polarizable force field. We evaluated GQ structure, ion sampling, core energetics, ion dehydration and binding, and ion competition and found that each analysis supported the known GQ-ion specificity trend (K+ > Na+ > Li+). We also found that K+ ions coordinate in the tetrad core antiprismatically, whereas Na+ and Li+ align coplanar to guanine tetrads, partially because of their attraction to surrounding water. Further, we showed that K+ occupancy is higher around the c-kit1 GQ and its nucleobases than Na+ and Li+, which tend to interact with backbone and sugar moieties. Finally, we showed that K+ binding to the c-kit1 GQ is faster and more frequent than Na+ and Li+. Such descriptions of GQ-ion dynamics suggest the rate of dehydration as the dominant factor for preference of K+ by DNA GQs and provide insight into noncanonical nucleic acids for which little experimental data exist.
Collapse
Affiliation(s)
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, Virginia; Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia.
| |
Collapse
|
20
|
Ortiz de Luzuriaga I, Lopez X, Gil A. Learning to Model G-Quadruplexes: Current Methods and Perspectives. Annu Rev Biophys 2021; 50:209-243. [PMID: 33561349 DOI: 10.1146/annurev-biophys-060320-091827] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
G-quadruplexes have raised considerable interest during the past years for the development of therapies against cancer. These noncanonical structures of DNA may be found in telomeres and/or oncogene promoters, and it has been observed that the stabilization of such G-quadruplexes may disturb tumor cell growth. Nevertheless, the mechanisms leading to folding and stabilization of these G-quadruplexes are still not well established, and they are the focus of much current work in this field. In seminal works, stabilization was observed to be produced by cations. However, subsequent studies showed that different kinds of small molecules, from planar and nonplanar organic molecules to square-planar and octahedral metal complexes, may also lead to the stabilization of G-quadruplexes. Thus, the comprehension and rationalization of the interaction of these small molecules with G-quadruplexes are also important topics of current interest in medical applications. To shed light on the questions arising from the literature on the formation of G-quadruplexes, their stabilization, and their interaction with small molecules, synergies between experimental studies and computational works are needed. In this review, we mainly focus on in silico approaches and provide a broad compilation of different leading studies carried out to date by different computational methods. We divide these methods into twomain categories: (a) classical methods, which allow for long-timescale molecular dynamics simulations and the corresponding analysis of dynamical information, and (b) quantum methods (semiempirical, quantum mechanics/molecular mechanics, and density functional theory methods), which allow for the explicit simulation of the electronic structure of the system but, in general, are not capable of being used in long-timescale molecular dynamics simulations and, therefore, give a more static picture of the relevant processes.
Collapse
Affiliation(s)
- Iker Ortiz de Luzuriaga
- CIC nanoGUNE BRTA, 20018 Donostia, Euskadi, Spain; .,Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Uniberstitatea, UPV/EHU, 20080 Donostia, Euskadi, Spain
| | - Xabier Lopez
- Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Uniberstitatea, UPV/EHU, 20080 Donostia, Euskadi, Spain.,Donostia International Physics Center, 20018 Donostia, Spain
| | - Adrià Gil
- CIC nanoGUNE BRTA, 20018 Donostia, Euskadi, Spain; .,BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal;
| |
Collapse
|
21
|
Stadlbauer P, Islam B, Otyepka M, Chen J, Monchaud D, Zhou J, Mergny JL, Šponer J. Insights into G-Quadruplex-Hemin Dynamics Using Atomistic Simulations: Implications for Reactivity and Folding. J Chem Theory Comput 2021; 17:1883-1899. [PMID: 33533244 DOI: 10.1021/acs.jctc.0c01176] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Guanine quadruplex nucleic acids (G4s) are involved in key biological processes such as replication or transcription. Beyond their biological relevance, G4s find applications as biotechnological tools since they readily bind hemin and enhance its peroxidase activity, creating a G4-DNAzyme. The biocatalytic properties of G4-DNAzymes have been thoroughly studied and used for biosensing purposes. Despite hundreds of applications and massive experimental efforts, the atomistic details of the reaction mechanism remain unclear. To help select between the different hypotheses currently under investigation, we use extended explicit-solvent molecular dynamics (MD) simulations to scrutinize the G4/hemin interaction. We find that besides the dominant conformation in which hemin is stacked atop the external G-quartets, hemin can also transiently bind to the loops and be brought to the external G-quartets through diverse delivery mechanisms. The simulations do not support the catalytic mechanism relying on a wobbling guanine. Similarly, the catalytic role of the iron-bound water molecule is not in line with our results; however, given the simulation limitations, this observation should be considered with some caution. The simulations rather suggest tentative mechanisms in which the external G-quartet itself could be responsible for the unique H2O2-promoted biocatalytic properties of the G4/hemin complexes. Once stacked atop a terminal G-quartet, hemin rotates about its vertical axis while readily sampling shifted geometries where the iron transiently contacts oxygen atoms of the adjacent G-quartet. This dynamics is not apparent from the ensemble-averaged structure. We also visualize transient interactions between the stacked hemin and the G4 loops. Finally, we investigated interactions between hemin and on-pathway folding intermediates of the parallel-stranded G4 fold. The simulations suggest that hemin drives the folding of parallel-stranded G4s from slip-stranded intermediates, acting as a G4 chaperone. Limitations of the MD technique are briefly discussed.
Collapse
Affiliation(s)
- Petr Stadlbauer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Barira Islam
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute, Palacký University, Křížkovského 8, 779 00 Olomouc, Czech Republic
| | - Jielin Chen
- State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry & Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - David Monchaud
- Institut de Chimie Moléculaire (ICMUB), CNRS UMR6302, UBFC, Dijon 21078, France
| | - Jun Zhou
- State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry & Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Jean-Louis Mergny
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,State Key Laboratory of Analytical Chemistry for Life Science, School of Chemistry & Chemical Engineering, Nanjing University, Nanjing 210023, China.,Laboratoire d'Optique et Biosciences, Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, 91128 Palaiseau cedex, France
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
22
|
Ratnasinghe BD, Salsbury AM, Lemkul JA. Ion Binding Properties and Dynamics of the bcl-2 G-Quadruplex Using a Polarizable Force Field. J Chem Inf Model 2020; 60:6476-6488. [PMID: 33264004 DOI: 10.1021/acs.jcim.0c01064] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
G-quadruplexes (GQs) are topologically diverse, highly thermostable noncanonical nucleic acid structures that form in guanine-rich sequences in DNA and RNA. GQs are implicated in transcriptional and translational regulation and genome maintenance, and deleterious alterations to their structures contribute to diseases such as cancer. The expression of the B-cell lymphoma 2 (Bcl-2) antiapoptotic protein, for example, is under transcriptional control of a GQ in the promoter of the bcl-2 gene. Modulation of the bcl-2 GQ by small molecules is of interest for chemotherapeutic development but doing so requires knowledge of the factors driving GQ folding and stabilization. To develop a greater understanding of the electrostatic properties of the bcl-2 promoter GQ, we performed molecular dynamics simulations using the Drude-2017 polarizable force field and compared relevant outcomes to the nonpolarizable CHARMM36 force field. Our simulation outcomes highlight the importance of dipole-dipole interactions in the bcl-2 GQ, particularly during the recruitment of a bulk K+ ion to the solvent-exposed face of the tetrad stem. We also predict and characterize an "electronegative pocket" at the tetrad-long loop junction that induces local backbone conformational change and may induce local conformational changes at cellular concentrations of K+. These outcomes suggest that moieties within the bcl-2 GQ can be targeted by small molecules to modulate bcl-2 GQ stability.
Collapse
Affiliation(s)
- Brian D Ratnasinghe
- Department of Biochemistry, Virginia Tech, 303 Engel Hall, 340 West Campus Dr., Blacksburg, Virginia 24061, United States
| | - Alexa M Salsbury
- Department of Biochemistry, Virginia Tech, 303 Engel Hall, 340 West Campus Dr., Blacksburg, Virginia 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, 303 Engel Hall, 340 West Campus Dr., Blacksburg, Virginia 24061, United States
| |
Collapse
|
23
|
Salsbury AM, Lemkul JA. Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics. Curr Opin Struct Biol 2020; 67:9-17. [PMID: 32950937 DOI: 10.1016/j.sbi.2020.08.003] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 08/04/2020] [Accepted: 08/13/2020] [Indexed: 01/24/2023]
Abstract
Nucleic acids play critical roles in carrying genetic information, participating in catalysis, and preserving chromosomal structure. Despite over a century of study, efforts to understand the dynamics and structure-function relationships of DNA and RNA at the atomic level are still ongoing. Molecular dynamics (MD) simulations augment experiments by providing atomistic resolution and quantitative relationships between structure and conformational energy. Steady advancements in computer hardware, software, and atomistic force fields (FFs) over 40 years have facilitated new discoveries. Here, we review nucleic acid FF development with emphasis on recent refinements that have improved descriptions of important nucleic acid properties. We then discuss several key examples of successes and challenges in modeling nucleic acid structure and dynamics using the latest FFs.
Collapse
Affiliation(s)
- Alexa M Salsbury
- Department of Biochemistry, Virginia Tech, Blacksburg, VA 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Tech, Blacksburg, VA 24061, United States; Center for Drug Discovery, Virginia Tech, Blacksburg, VA 24061, United States.
| |
Collapse
|
24
|
Islam B, Stadlbauer P, Vorlíčková M, Mergny JL, Otyepka M, Šponer J. Stability of Two-Quartet G-Quadruplexes and Their Dimers in Atomistic Simulations. J Chem Theory Comput 2020; 16:3447-3463. [PMID: 32163706 DOI: 10.1021/acs.jctc.9b01068] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
G-quadruplexes (GQs) are four-stranded noncanonical DNA and RNA architectures that can be formed by guanine-rich sequences. The stability of GQs increases with the number of G-quartets, and three G-quartets generally form stable GQs. However, the stability of two-quartet GQs is an open issue. To understand the intrinsic stability of two-quartet GQ stems, we have carried out a series of unbiased molecular dynamics (MD) simulations (505 μs in total) of two- and four-quartet DNA and RNA GQs, with attention paid mainly to parallel-stranded arrangements. We used AMBER DNA parmOL15 and RNA parmOL3 force fields and tested different ion and water models. Two-quartet parallel-stranded DNA GQs unfolded in all the simulations, while the equivalent RNA GQ was stable in most of the simulations. GQs composed of two stacked units of two-quartet GQs were stable for both DNA and RNA. The simulations suggest that a minimum of three quartets are needed to form an intrinsically stable all-anti parallel-stranded DNA GQ. Parallel two-quartet DNA GQ may exist if substantially stabilized by another molecule or structural element, including multimerization. On the other hand, we predict that isolated RNA two-quartet parallel GQs may form, albeit being weakly stable. We also show that ionic parameters and water models should be chosen with caution because some parameter combinations can cause spurious instability of GQ stems. Some in-so-far unnoticed limitations of force-field description of multiple ions inside the GQs are discussed, which compromise the capability of simulations to fully capture the effect of increase in the number of quartets on the GQ stability.
Collapse
Affiliation(s)
- Barira Islam
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michaela Vorlíčková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jean-Louis Mergny
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic.,Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
25
|
Salsbury AM, Dean TJ, Lemkul JA. Polarizable Molecular Dynamics Simulations of Two c-kit Oncogene Promoter G-Quadruplexes: Effect of Primary and Secondary Structure on Loop and Ion Sampling. J Chem Theory Comput 2020; 16:3430-3444. [PMID: 32307997 DOI: 10.1021/acs.jctc.0c00191] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
G-quadruplexes (GQs) are highly ordered nucleic acid structures that play fundamental roles in regulating gene expression and maintaining genomic stability. GQs are topologically diverse and enriched in promoter sequences of growth regulatory genes and proto-oncogenes, suggesting that they may serve as attractive targets for drug design at the level of transcription rather than inhibiting the activity of the protein products of these genes. The c-kit promoter contains three adjacent GQ-forming sequences that have proposed antagonistic effects on gene expression and thus are promising drug targets for diseases such as gastrointestinal stromal tumors, mast cell disease, and leukemia. Because GQ stability is influenced by primary structure, secondary structure, and ion interactions, a greater understanding of GQ structure, dynamics, and ion binding properties is needed to develop novel, GQ-targeting therapeutics. Here, we performed molecular dynamics simulations to systematically study the c-kit2 and c-kit* GQs, evaluating nonpolarizable and polarizable force fields (FFs) and examining the effects of base substitutions and cation type (K+, Na+, and Li+) on the dynamics of their isolated and linked structures. We found that the Drude polarizable FF outperformed the additive CHARMM36 FF in two- and three-tetrad GQs and solutions of KCl, NaCl, and LiCl. Drude simulations with different cations agreed with the known GQ stabilization preference (K+ > Na+ > Li+) and illustrated that tetrad core-ion coordination differs as a function of cation type. Finally, we showed that differences in primary and secondary structure influence loop sampling, ion binding, and core-ion energetics of GQs.
Collapse
Affiliation(s)
- Alexa M Salsbury
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Tanner J Dean
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Justin A Lemkul
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
26
|
Milovanović B, Stanojević A, Etinski M, Petković M. Intriguing Intermolecular Interplay in Guanine Quartet Complexes with Alkali and Alkaline Earth Cations. J Phys Chem B 2020; 124:3002-3014. [DOI: 10.1021/acs.jpcb.0c01165] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Branislav Milovanović
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Ana Stanojević
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| | - Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11158 Belgrade, Serbia
| |
Collapse
|
27
|
Lemkul JA. Same fold, different properties: polarizable molecular dynamics simulations of telomeric and TERRA G-quadruplexes. Nucleic Acids Res 2020; 48:561-575. [PMID: 31807754 PMCID: PMC6954416 DOI: 10.1093/nar/gkz1154] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 11/21/2019] [Accepted: 11/26/2019] [Indexed: 12/12/2022] Open
Abstract
DNA and RNA sequences rich in guanine can fold into noncanonical structures called G-quadruplexes (GQs), which exhibit a common stem structure of Hoogsteen hydrogen-bonded guanine tetrads and diverse loop structures. GQ sequence motifs are overrepresented in promoters, origins of replication, telomeres, and untranslated regions in mRNA, suggesting roles in modulating gene expression and preserving genomic integrity. Given these roles and unique aspects of different structures, GQs are attractive targets for drug design, but greater insight into GQ folding pathways and the interactions stabilizing them is required. Here, we performed molecular dynamics simulations to study two bimolecular GQs, a telomeric DNA GQ and the analogous telomeric repeat-containing RNA (TERRA) GQ. We applied the Drude polarizable force field, which we show outperforms the additive CHARMM36 force field in both ion retention and maintenance of the GQ folds. The polarizable simulations reveal that the GQs bind bulk K+ ions differently, and that the TERRA GQ accumulates more K+ ions, suggesting different ion interactions stabilize these structures. Nucleobase dipole moments vary as a function of position and also contribute to ion binding. Finally, we show that the TERRA GQ is more sensitive than the telomeric DNA GQ to water-mediated modulation of ion-induced dipole-dipole interactions.
Collapse
Affiliation(s)
- Justin A Lemkul
- Department of Biochemistry and Center for Drug Discovery, Virginia Tech, Blacksburg, VA 24061, USA
| |
Collapse
|
28
|
Molecular dynamics simulations of G-quadruplexes: The basic principles and their application to folding and ligand binding. ANNUAL REPORTS IN MEDICINAL CHEMISTRY 2020. [DOI: 10.1016/bs.armc.2020.04.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
29
|
Stadlbauer P, Kührová P, Vicherek L, Banáš P, Otyepka M, Trantírek L, Šponer J. Parallel G-triplexes and G-hairpins as potential transitory ensembles in the folding of parallel-stranded DNA G-Quadruplexes. Nucleic Acids Res 2019; 47:7276-7293. [PMID: 31318975 PMCID: PMC6698752 DOI: 10.1093/nar/gkz610] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2019] [Revised: 06/26/2019] [Accepted: 07/03/2019] [Indexed: 12/29/2022] Open
Abstract
Guanine quadruplexes (G4s) are non-canonical nucleic acids structures common in important genomic regions. Parallel-stranded G4 folds are the most abundant, but their folding mechanism is not fully understood. Recent research highlighted that G4 DNA molecules fold via kinetic partitioning mechanism dominated by competition amongst diverse long-living G4 folds. The role of other intermediate species such as parallel G-triplexes and G-hairpins in the folding process has been a matter of debate. Here, we use standard and enhanced-sampling molecular dynamics simulations (total length of ∼0.9 ms) to study these potential folding intermediates. We suggest that parallel G-triplex per se is rather an unstable species that is in local equilibrium with a broad ensemble of triplex-like structures. The equilibrium is shifted to well-structured G-triplex by stacked aromatic ligand and to a lesser extent by flanking duplexes or nucleotides. Next, we study propeller loop formation in GGGAGGGAGGG, GGGAGGG and GGGTTAGGG sequences. We identify multiple folding pathways from different unfolded and misfolded structures leading towards an ensemble of intermediates called cross-like structures (cross-hairpins), thus providing atomistic level of description of the single-molecule folding events. In summary, the parallel G-triplex is a possible, but not mandatory short-living (transitory) intermediate in the folding of parallel-stranded G4.
Collapse
Affiliation(s)
- Petr Stadlbauer
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacky University, Šlechtitelů 27, 771 46 Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacky University, Šlechtitelů 27, 771 46 Olomouc, Czech Republic.,Department of Physical Chemistry, Faculty of Science, Palacky University, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lukáš Vicherek
- Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacky University, Šlechtitelů 27, 771 46 Olomouc, Czech Republic.,Department of Physical Chemistry, Faculty of Science, Palacky University, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacky University, Šlechtitelů 27, 771 46 Olomouc, Czech Republic.,Department of Physical Chemistry, Faculty of Science, Palacky University, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lukáš Trantírek
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, 612 65 Brno, Czech Republic.,Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, v. v. i., Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacky University, Šlechtitelů 27, 771 46 Olomouc, Czech Republic.,Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
| |
Collapse
|
30
|
Jing Z, Liu C, Cheng SY, Qi R, Walker BD, Piquemal JP, Ren P. Polarizable Force Fields for Biomolecular Simulations: Recent Advances and Applications. Annu Rev Biophys 2019; 48:371-394. [PMID: 30916997 DOI: 10.1146/annurev-biophys-070317-033349] [Citation(s) in RCA: 216] [Impact Index Per Article: 43.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Realistic modeling of biomolecular systems requires an accurate treatment of electrostatics, including electronic polarization. Due to recent advances in physical models, simulation algorithms, and computing hardware, biomolecular simulations with advanced force fields at biologically relevant timescales are becoming increasingly promising. These advancements have not only led to new biophysical insights but also afforded opportunities to advance our understanding of fundamental intermolecular forces. This article describes the recent advances and applications, as well as future directions, of polarizable force fields in biomolecular simulations.
Collapse
Affiliation(s)
- Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Sara Y Cheng
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Rui Qi
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Brandon D Walker
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA; .,Sorbonne Université, CNRS, Laboratoire de Chimie Theórique, 75252 Paris CEDEX 05, France.,Institut Universitaire de France, 75005 Paris, France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA;
| |
Collapse
|
31
|
Andrałojć W, Małgowska M, Sarzyńska J, Pasternak K, Szpotkowski K, Kierzek R, Gdaniec Z. Unraveling the structural basis for the exceptional stability of RNA G-quadruplexes capped by a uridine tetrad at the 3' terminus. RNA (NEW YORK, N.Y.) 2019; 25:121-134. [PMID: 30341177 PMCID: PMC6298561 DOI: 10.1261/rna.068163.118] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2018] [Accepted: 10/16/2018] [Indexed: 05/24/2023]
Abstract
Uridine tetrads (U-tetrads) are a structural element encountered in RNA G-quadruplexes, for example, in the structures formed by the biologically relevant human telomeric repeat RNA. For these molecules, an unexpectedly strong stabilizing influence of a U-tetrad forming at the 3' terminus of a quadruplex was reported. Here we present the high-resolution solution NMR structure of the r(UGGUGGU)4 quadruplex which, in our opinion, provides an explanation for this stabilization. Our structure features a distinctive, abrupt chain reversal just prior to the 3' uridine tetrad. Similar "reversed U-tetrads" were already observed in the crystalline phase. However, our NMR structure coupled with extensive explicit solvent molecular dynamics (MD) simulations identifies some key features of this motif that up to now remained overlooked. These include the presence of an exceptionally stable 2'OH to phosphate hydrogen bond, as well as the formation of an additional K+ binding pocket in the quadruplex groove.
Collapse
Affiliation(s)
- Witold Andrałojć
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704 Poznan, Poland
| | - Magdalena Małgowska
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704 Poznan, Poland
| | - Joanna Sarzyńska
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704 Poznan, Poland
| | - Karol Pasternak
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704 Poznan, Poland
| | - Kamil Szpotkowski
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704 Poznan, Poland
| | - Ryszard Kierzek
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704 Poznan, Poland
| | - Zofia Gdaniec
- Institute of Bioorganic Chemistry, Polish Academy of Sciences, 61-704 Poznan, Poland
| |
Collapse
|
32
|
Salsbury AM, Lemkul JA. Molecular Dynamics Simulations of the c-kit1 Promoter G-Quadruplex: Importance of Electronic Polarization on Stability and Cooperative Ion Binding. J Phys Chem B 2018; 123:148-159. [DOI: 10.1021/acs.jpcb.8b11026] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
|
33
|
Islam B, Stadlbauer P, Krepl M, Havrila M, Haider S, Sponer J. Structural Dynamics of Lateral and Diagonal Loops of Human Telomeric G-Quadruplexes in Extended MD Simulations. J Chem Theory Comput 2018; 14:5011-5026. [DOI: 10.1021/acs.jctc.8b00543] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Barira Islam
- Institute of Biophysics
of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics
of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics
of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Marek Havrila
- Institute of Biophysics
of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Shozeb Haider
- UCL School of
Pharmacy, 29-39 Brunswick Square, London WC1N 1AX, U.K
| | - Jiri Sponer
- Institute of Biophysics
of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
34
|
Villani G. Quantum Mechanical Investigation of the G-Quadruplex Systems of Human Telomere. ACS OMEGA 2018; 3:9934-9944. [PMID: 31459122 PMCID: PMC6644616 DOI: 10.1021/acsomega.8b01678] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/17/2018] [Accepted: 08/09/2018] [Indexed: 05/17/2023]
Abstract
The three G-quadruplexes involved in the human telomere have been studied with an accurate quantum mechanical approach, and the possibility of reducing them to a simpler model has been tested. The similarities and the differences of these three systems are shown and discussed. Each system has been analyzed through different properties and compared to the others. In particular, we have considered: (1) the shape of the cavity and the atomic charges around it; (2) the electric field in and out of the cavity; (3) the stabilization energy due to the stacking of G-tetrads, to the H-bonds and to the ion interactions; and, finally, (4) to study the mechanism of the process of the ion inclusion in the cavity, the curves of potential energy due to the movement of the Na+ and K+ ions toward the cavity. The results suggest that a detailed study is essential in order to obtain the quantitative properties of these complex systems, but also that some qualitative behaviors can be schematized. Our study makes it clear that the entry of an ion in the cavity of these systems is a complex process, where it is possible to find stable structures with the ion out and in the cavity. Moreover, it is possible that more than one diabatic state is involved in this process.
Collapse
Affiliation(s)
- Giovanni Villani
- Istituto di Chimica dei Composti OrganoMetallici, ICCOM—CNR
(UOS Pisa), Area della Ricerca di Pisa, Via G. Moruzzi, 1, I-56124 Pisa, Italy
| |
Collapse
|
35
|
Gülbakan B, Barylyuk K, Schneider P, Pillong M, Schneider G, Zenobi R. Native Electrospray Ionization Mass Spectrometry Reveals Multiple Facets of Aptamer–Ligand Interactions: From Mechanism to Binding Constants. J Am Chem Soc 2018; 140:7486-7497. [DOI: 10.1021/jacs.7b13044] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Basri Gülbakan
- Department of Chemistry and Applied Bioscience, ETH Zürich, CH-8093 Zürich, Switzerland
- Hacettepe University Institute of Child Health, Ihsan Dogramaci Children’s Hospital, Sıhhiye Square, 06100 Ankara, Turkey
| | - Konstantin Barylyuk
- Department of Chemistry and Applied Bioscience, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Petra Schneider
- Department of Chemistry and Applied Bioscience, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Max Pillong
- Department of Chemistry and Applied Bioscience, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Gisbert Schneider
- Department of Chemistry and Applied Bioscience, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Renato Zenobi
- Department of Chemistry and Applied Bioscience, ETH Zürich, CH-8093 Zürich, Switzerland
| |
Collapse
|
36
|
Š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: 327] [Impact Index Per Article: 54.5] [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
|
37
|
Zhang C, Lu C, Jing Z, Wu C, Piquemal JP, Ponder JW, Ren P. AMOEBA Polarizable Atomic Multipole Force Field for Nucleic Acids. J Chem Theory Comput 2018; 14:2084-2108. [PMID: 29438622 PMCID: PMC5893433 DOI: 10.1021/acs.jctc.7b01169] [Citation(s) in RCA: 152] [Impact Index Per Article: 25.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, 3 J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.
Collapse
Affiliation(s)
- Changsheng Zhang
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Chao Lu
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Zhifeng Jing
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| | - Chuanjie Wu
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC, UMR7616 CNRS, Paris, France
| | - Jay W. Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, MO 63130, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
38
|
Yoo J, Aksimentiev A. New tricks for old dogs: improving the accuracy of biomolecular force fields by pair-specific corrections to non-bonded interactions. Phys Chem Chem Phys 2018; 20:8432-8449. [PMID: 29547221 DOI: 10.1039/c7cp08185e] [Citation(s) in RCA: 134] [Impact Index Per Article: 22.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
In contrast to ordinary polymers, the vast majority of biological macromolecules adopt highly ordered three-dimensional structures that define their functions. The key to folding of a biopolymer into a unique 3D structure or to an assembly of several biopolymers into a functional unit is a delicate balance between the attractive and repulsive forces that also makes such self-assembly reversible under physiological conditions. The all-atom molecular dynamics (MD) method has emerged as a powerful tool for studies of individual biomolecules and their functional assemblies, encompassing systems of ever increasing complexity. However, advances in parallel computing technology have outpaced the development of the underlying theoretical models-the molecular force fields, pushing the MD method into an untested territory. Recent tests of the MD method have found the most commonly used molecular force fields to be out of balance, overestimating attractive interactions between charged and hydrophobic groups, which can promote artificial aggregation in MD simulations of multi-component protein, nucleic acid, and lipid systems. One route towards improving the force fields is through the NBFIX corrections method, in which the intermolecular forces are calibrated against experimentally measured quantities such as osmotic pressure by making atom pair-specific adjustments to the non-bonded interactions. In this article, we review development of the NBFIX (Non-Bonded FIX) corrections to the AMBER and CHARMM force fields and discuss their implications for MD simulations of electrolyte solutions, dense DNA systems, Holliday junctions, protein folding, and lipid bilayer membranes.
Collapse
Affiliation(s)
- Jejoong Yoo
- Center for the Physics of Living Cells, Department of Physics and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA. and Center for Self-assembly and Complexity, Institute for Basic Science, Pohang, 37363, Republic of Korea
| | - Aleksei Aksimentiev
- Center for the Physics of Living Cells, Department of Physics and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA.
| |
Collapse
|
39
|
Siebenmorgen T, Zacharias M. Origin of Ion Specificity of Telomeric DNA G-Quadruplexes Investigated by Free-Energy Simulations. Biophys J 2017; 112:2280-2290. [PMID: 28591601 DOI: 10.1016/j.bpj.2017.04.036] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2017] [Revised: 04/12/2017] [Accepted: 04/24/2017] [Indexed: 11/29/2022] Open
Abstract
Telomeric DNA consists of tandem repeats of the sequence d(TTAGGG) that form G-quadruplex structures made of stacked guanines with monovalent cations bound at a central cavity. Although different ions can stabilize a G-quadruplex structure, the preferred bound ions are typically K+ or Na+. Several different strand-folding topologies have been reported for Q-quadruplexes formed from telomeric repeats depending on DNA length and ion solution condition. This suggests a possible dependence of the ion selectivity of the central pore on the folding topology of the quadruplex. Molecular dynamics free energy perturbation has been employed to systematically study the relative affinity of the central quadruplex pore for different cation types and the associated energetic and solvation contributions to ion selectivity. The calculations have been performed on two different common quadruplex folding topologies. For both topologies, the same ion selectivity was found with a preference for K+ followed by Rb+ and Na+, which agrees with the experimentally determined preference for most investigated quadruplexes. The selectivity is determined by a balance between attractive Coulomb interactions and loss of hydration but also modulated by van der Waals contributions. Specificity is mediated by the central guanines and no significant contribution of the nucleic acid backbone. The simulations indicate that different topologies might be stabilized by ions bound at the surface or alternative sites of the quadruplex because the ion specificity of the central pore does not depend on the strand folding topology.
Collapse
Affiliation(s)
- Till Siebenmorgen
- Physics Department, Technical University of Munich, Garching, Germany
| | - Martin Zacharias
- Physics Department, Technical University of Munich, Garching, Germany.
| |
Collapse
|
40
|
Havrila M, Stadlbauer P, Islam B, Otyepka M, Šponer J. Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes. J Chem Theory Comput 2017; 13:3911-3926. [PMID: 28657760 DOI: 10.1021/acs.jctc.7b00257] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]4 GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]4 simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na+ ions move inside the GQs in a concerted manner, while larger relocations of the K+ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K+ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.
Collapse
Affiliation(s)
- Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Barira Islam
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
41
|
Lemkul JA, MacKerell AD. Polarizable Force Field for DNA Based on the Classical Drude Oscillator: II. Microsecond Molecular Dynamics Simulations of Duplex DNA. J Chem Theory Comput 2017; 13:2072-2085. [PMID: 28398748 PMCID: PMC5485260 DOI: 10.1021/acs.jctc.7b00068] [Citation(s) in RCA: 61] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The structure and dynamics of DNA are governed by a sensitive balance between base stacking and pairing, hydration, and interactions with ions. Force-field models that include explicit representations of electronic polarization are capable of more accurately modeling the subtle details of these interactions versus commonly used additive force fields. In this work, we validate our recently refined polarizable force field for DNA based on the classical Drude oscillator model, in which electronic degrees of freedom are represented as negatively charged particles attached to their parent atoms via harmonic springs. The previous version of the force field, called Drude-2013, produced stable A- and B-DNA trajectories on the order of hundreds of nanoseconds, but deficiencies were identified that included weak base stacking ultimately leading to distortion of B-DNA duplexes and unstable Z-DNA. As a result of extensive refinement of base nonbonded terms and bonded parameters in the deoxyribofuranose sugar and phosphodiester backbone, we demonstrate that the new version of the Drude DNA force field is capable of simulating A- and B-forms of DNA on the microsecond time scale and the resulting conformational ensembles agree well with a broad set of experimental properties, including solution X-ray scattering profiles. In addition, simulations of Z-form duplex DNA in its crystal environment are stable on the order of 100 ns. The revised force field is to be called Drude-2017.
Collapse
Affiliation(s)
- Justin A. Lemkul
- 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
|
42
|
Lemkul JA, MacKerell AD. Polarizable Force Field for DNA Based on the Classical Drude Oscillator: I. Refinement Using Quantum Mechanical Base Stacking and Conformational Energetics. J Chem Theory Comput 2017; 13:2053-2071. [PMID: 28399366 DOI: 10.1021/acs.jctc.7b00067] [Citation(s) in RCA: 59] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Empirical force fields seek to relate the configuration of a set of atoms to its energy, thus yielding the forces governing its dynamics, using classical physics rather than more expensive quantum mechanical calculations that are computationally intractable for large systems. Most force fields used to simulate biomolecular systems use fixed atomic partial charges, neglecting the influence of electronic polarization, instead making use of a mean-field approximation that may not be transferable across environments. Recent hardware and software developments make polarizable simulations feasible, and to this end, polarizable force fields represent the next generation of molecular dynamics simulation technology. In this work, we describe the refinement of a polarizable force field for DNA based on the classical Drude oscillator model by targeting quantum mechanical interaction energies and conformational energy profiles of model compounds necessary to build a complete DNA force field. The parametrization strategy employed in the present work seeks to correct weak base stacking in A- and B-DNA and the unwinding of Z-DNA observed in the previous version of the force field, called Drude-2013. Refinement of base nonbonded terms and reparametrization of dihedral terms in the glycosidic linkage, deoxyribofuranose rings, and important backbone torsions resulted in improved agreement with quantum mechanical potential energy surfaces. Notably, we expand on previous efforts by explicitly including Z-DNA conformational energetics in the refinement.
Collapse
Affiliation(s)
- Justin A Lemkul
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| |
Collapse
|
43
|
Gresh N, Naseem-Khan S, Lagardère L, Piquemal JP, Sponer JE, Sponer J. Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li +, Na +, K +, and Rb + Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential. J Phys Chem B 2017; 121:3997-4014. [PMID: 28363025 DOI: 10.1021/acs.jpcb.7b01836] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K+ and Na+, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li+, Na+, K+, or Rb+) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K+ or two Rb+ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.
Collapse
Affiliation(s)
- Nohad Gresh
- Laboratoire de Chimie Théorique, Sorbonne Universités , UPMC, UMR7616 CNRS, 75006Paris, France
| | - Sehr Naseem-Khan
- Laboratoire de Chimie Théorique, Sorbonne Universités , UPMC, UMR7616 CNRS, 75006Paris, France
| | - Louis Lagardère
- Laboratoire de Chimie Théorique, Sorbonne Universités , UPMC, UMR7616 CNRS, 75006Paris, France
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Universités , UPMC, UMR7616 CNRS, 75006Paris, France.,Institut Universitaire de France, Paris Cedex 05, 75231, France.,Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas, 78712, United States
| | - Judit E Sponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovpolska 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiri Sponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovpolska 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
44
|
Yurenko YP, Novotný J, Marek R. Weak Supramolecular Interactions Governing Parallel and Antiparallel DNA Quadruplexes: Insights from Large-Scale Quantum Mechanics Analysis of Experimentally Derived Models. Chemistry 2017; 23:5573-5584. [PMID: 28225208 DOI: 10.1002/chem.201700236] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2017] [Indexed: 12/30/2022]
Abstract
The topology and energetics of guanine (G) quadruplexes is governed by supramolecular interactions within their strands. In this work, an extensive quantum mechanical (QM) study has been performed to analyze supramolecular interactions that shape the stems of (4+0) parallel (P) and (2+2) antiparallel (AP) quadruplex systems. The large-scale (≈400 atoms) models of P and AP were constructed from high-quality experimental structures. The results provide evidence that each of the P and AP structures is shaped by a distinct network of supramolecular interactions. Analysis of electron topological characteristics of hydrogen bonds in P and AP systems indicates that the P model benefits from stronger intratetrad hydrogen bonding. For intertetrad stacking interactions, both noncovalent interaction plot and energy decomposition analysis approaches suggest that the stem of the P quadruplex benefits more from stacking than that of the AP stem; the difference in energetic stabilization for the two topologies is about 10 %. Stronger hydrogen-bonding and stacking interactions in the stem of the P quadruplex, relative to those in the AP system, can be an important indicator to explain the experimental observations that guanine-rich oligonucleotides tend to form all-parallel stems with an all-anti orientation of nucleobases. However, in addition to intrinsic stabilization, partial desolvation effects, which affect the energetics and dynamics of the G-quadruplex folding process, call for further investigations.
Collapse
Affiliation(s)
- Yevgen P Yurenko
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| | - Jan Novotný
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| | - Radek Marek
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| |
Collapse
|
45
|
Zhang C, Bell D, Harger M, Ren P. Polarizable Multipole-Based Force Field for Aromatic Molecules and Nucleobases. J Chem Theory Comput 2017; 13:666-678. [PMID: 28030769 PMCID: PMC5312700 DOI: 10.1021/acs.jctc.6b00918] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Aromatic molecules with π electrons
are commonly involved
in chemical and biological recognitions. For example, nucleobases
play central roles in DNA/RNA structure and their interactions with
proteins. The delocalization of the π electrons is responsible
for the high polarizability of aromatic molecules. In this work, the
AMOEBA force field has been developed and applied to 5 regular nucleobases
and 12 aromatic molecules. The permanent electrostatic energy is expressed
as atomic multipole interactions between atom pairs, and many-body
polarization is accounted for by mutually induced atomic dipoles.
We have systematically investigated aromatic ring stacking and aromatic-water
interactions for nucleobases and aromatic molecules, as well as base–base
hydrogen-bonding pair interactions, all at various distances and orientations.
van der Waals parameters were determined by comparison to the quantum
mechanical interaction energy of these dimers and fine-tuned using
condensed phase simulation. By comparing to quantum mechanical calculations,
we show that the resulting classical potential is able to accurately
describe molecular polarizability, molecular vibrational frequency,
and dimer interaction energy of these aromatic systems. Condensed
phase properties, including hydration free energy, liquid density,
and heat of vaporization, are also in good overall agreement with
experimental values. The structures of benzene liquid phase and benzene-water
solution were also investigated by simulation and compared with experimental
and PDB structure derived statistical results.
Collapse
Affiliation(s)
- Changsheng Zhang
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| | - David Bell
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Matthew Harger
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| |
Collapse
|
46
|
Šponer J, Bussi G, Stadlbauer P, Kührová P, Banáš P, Islam B, Haider S, Neidle S, Otyepka M. Folding of guanine quadruplex molecules-funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies. Biochim Biophys Acta Gen Subj 2016; 1861:1246-1263. [PMID: 27979677 DOI: 10.1016/j.bbagen.2016.12.008] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2016] [Revised: 12/04/2016] [Accepted: 12/11/2016] [Indexed: 01/18/2023]
Abstract
BACKGROUND Guanine quadruplexes (GQs) play vital roles in many cellular processes and are of much interest as drug targets. In contrast to the availability of many structural studies, there is still limited knowledge on GQ folding. SCOPE OF REVIEW We review recent molecular dynamics (MD) simulation studies of the folding of GQs, with an emphasis paid to the human telomeric DNA GQ. We explain the basic principles and limitations of all types of MD methods used to study unfolding and folding in a way accessible to non-specialists. We discuss the potential role of G-hairpin, G-triplex and alternative GQ intermediates in the folding process. We argue that, in general, folding of GQs is fundamentally different from funneled folding of small fast-folding proteins, and can be best described by a kinetic partitioning (KP) mechanism. KP is a competition between at least two (but often many) well-separated and structurally different conformational ensembles. MAJOR CONCLUSIONS The KP mechanism is the only plausible way to explain experiments reporting long time-scales of GQ folding and the existence of long-lived sub-states. A significant part of the natural partitioning of the free energy landscape of GQs comes from the ability of the GQ-forming sequences to populate a large number of syn-anti patterns in their G-tracts. The extreme complexity of the KP of GQs typically prevents an appropriate description of the folding landscape using just a few order parameters or collective variables. GENERAL SIGNIFICANCE We reconcile available computational and experimental studies of GQ folding and formulate basic principles characterizing GQ folding landscapes. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, 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.
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, 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
| | - Petra Kührová
- 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
| | - Pavel Banáš
- 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
| | - Barira Islam
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Shozeb Haider
- UCL School of Pharmacy, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Stephen Neidle
- UCL School of Pharmacy, 29-39 Brunswick Square, London WC1N 1AX, UK
| | - Michal Otyepka
- 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
|
47
|
Stadlbauer P, Mazzanti L, Cragnolini T, Wales DJ, Derreumaux P, Pasquali S, Šponer J. Coarse-Grained Simulations Complemented by Atomistic Molecular Dynamics Provide New Insights into Folding and Unfolding of Human Telomeric G-Quadruplexes. J Chem Theory Comput 2016; 12:6077-6097. [DOI: 10.1021/acs.jctc.6b00667] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Petr Stadlbauer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Departments of Physical
Chemistry, Faculty of Science, Palacký University, 17. listopadu
1192/12, 771 46 Olomouc, Czech Republic
| | - Liuba Mazzanti
- Laboratoire
de Biochimie Théorique, IBPC, CNRS UPR9080, Université Sorbonne Paris Cite, Paris Diderot, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Tristan Cragnolini
- Department
of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, U.K
| | - David J. Wales
- Department
of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Philippe Derreumaux
- Laboratoire
de Biochimie Théorique, IBPC, CNRS UPR9080, Université Sorbonne Paris Cite, Paris Diderot, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Samuela Pasquali
- Laboratoire
de Biochimie Théorique, IBPC, CNRS UPR9080, Université Sorbonne Paris Cite, Paris Diderot, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Královopolská
135, 612 65 Brno, Czech Republic
| |
Collapse
|
48
|
Šponer J, Krepl M, Banáš P, Kührová P, Zgarbová M, Jurečka P, Havrila M, Otyepka M. How to understand atomistic molecular dynamics simulations of RNA and protein-RNA complexes? WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 8. [PMID: 27863061 DOI: 10.1002/wrna.1405] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2016] [Revised: 09/13/2016] [Accepted: 10/10/2016] [Indexed: 01/01/2023]
Abstract
We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| | - Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, Olomouc, Czech Republic
| |
Collapse
|
49
|
Ratcliff LE, Mohr S, Huhs G, Deutsch T, Masella M, Genovese L. Challenges in large scale quantum mechanical calculations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1290] [Citation(s) in RCA: 70] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Laura E. Ratcliff
- Argonne Leadership Computing Facility Argonne National Laboratory Lemon IL USA
| | - Stephan Mohr
- Department of Computer Applications in Science and Engineering Barcelona Supercomputing Center (BSC‐CNS) Barcelona Spain
| | - Georg Huhs
- Department of Computer Applications in Science and Engineering Barcelona Supercomputing Center (BSC‐CNS) Barcelona Spain
| | - Thierry Deutsch
- University Grenoble Alpes INAC‐MEM Grenoble France
- CEA, INAC‐MEM Grenoble France
| | - Michel Masella
- Laboratoire de Biologie Structurale et Radiologie, Service de Bioénergétique, Biologie Structurale et Mécanisme Institut de Biologie et de Technologie de Saclay, CEA Saclay Gif‐sur‐Yvette Cedex France
| | - Luigi Genovese
- University Grenoble Alpes INAC‐MEM Grenoble France
- CEA, INAC‐MEM Grenoble France
| |
Collapse
|
50
|
Vangaveti S, Ranganathan SV, Chen AA. Advances in RNA molecular dynamics: a simulator's guide to RNA force fields. WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 8. [DOI: 10.1002/wrna.1396] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/18/2016] [Revised: 08/17/2016] [Accepted: 08/25/2016] [Indexed: 01/04/2023]
Affiliation(s)
- Sweta Vangaveti
- The RNA Institute; University at Albany State University of New York; Albany NY USA
| | | | - Alan A. Chen
- The RNA Institute; University at Albany State University of New York; Albany NY USA
- Department of Chemistry; University at Albany State University of New York; Albany NY USA
| |
Collapse
|