1
|
Ortiz de Luzuriaga I, Sánchez-González Á, Synoradzki W, Lopez X, Gil A. Unravelling the binding affinity and selectivity of molybdenum(II) phenanthroline complexes with DNA G-quadruplexes by using linear-scaling DFT studies. The important role of ancillary ligands. Phys Chem Chem Phys 2022; 24:25918-25929. [PMID: 36260061 DOI: 10.1039/d2cp02241a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
We have used near linear-scaling density functional theory (LS-DFT) methods including dispersion, for the first time, to study the interaction of two isomers, equatorial (Eq) and axial (Ax), of the [Mo(η3-C3H5)Br(CO)2(phen)] metal complex with the DNA G-quadruplexes (GQ) to gain insight into its cytotoxicity. The LMKLL/DZDP level of calculation, which includes van der Waals contributions, with the SIESTA software was used to treat by means of first-principles computations the whole biological studied model system with ∼1000 atoms. Computed formation energies point to systems containing the Ax isomer as the most stable although the nearest system in energy containing the Eq isomer is only 7.5 kcal mol-1 above. On the other hand, the energy decomposition analysis (EDA) favours interaction energies for the systems containing the Eq isomer. However, when solvent effects are taken into account the systems containing the Ax isomer are again the most stable. This Ax isomer was found interacting by means of end-stacking with the GQ and surprisingly totally inside the non-canonical secondary structure, where all the ligands of the metal complex produce several weak interactions with the DNA structure. On the other hand, the Eq isomer prefers to interact from outside by means of intercalation in which the ancillary ligands also have some role in the interaction. Such features and comparison with the results regarding the interaction of the [Mo(η3-C3H5)Br(CO)2(phen)] metal complex with duplex DNA suggest that the [Mo(η3-C3H5)Br(CO)2(phen)] would have a higher affinity and eventual selectivity for non-canonical DNA GQ structures.
Collapse
Affiliation(s)
- Iker Ortiz de Luzuriaga
- CIC-nanoGUNE BRTA, Tolosa Hiribidea 76, E-20018, Donostia - San Sebastián, Euskadi, Spain.
- Polímero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, 20080, Donostia, Euskadi, Spain
| | - Ángel Sánchez-González
- BioISI - Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016, Lisboa, Portugal
| | - Wojciech Synoradzki
- BioISI - Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016, Lisboa, Portugal
| | - Xabier Lopez
- Polímero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, 20080, Donostia, Euskadi, Spain
- Donostia International Physics Center (DIPC), P. K. 1072, 20080, Donostia, Euskadi, Spain
| | - Adrià Gil
- CIC-nanoGUNE BRTA, Tolosa Hiribidea 76, E-20018, Donostia - San Sebastián, Euskadi, Spain.
- BioISI - Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, 1749-016, Lisboa, Portugal
- ARAID Foundation, Zaragoza, Spain
- Departamento de Química Inorgánica, Instituto de Síntesis Química y Catálisis Homogénea (ISQCH) - Consejo Superior de Investigaciones Científicas (CSIC). Universidad de Zaragoza, c/Pedro Cerbuna 12, 50009, Zaragoza, Spain
| |
Collapse
|
2
|
Chen J, Harper JB, Ho J. Improving the Accuracy of Quantum Mechanics/Molecular Mechanics (QM/MM) Models with Polarized Fragment Charges. J Chem Theory Comput 2022; 18:5607-5617. [PMID: 35952004 DOI: 10.1021/acs.jctc.2c00491] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This paper introduces an economical approach for improving the accuracy and convergence of quantum mechanics/molecular mechanics (QM/MM) models. The approach is tested on a series of neutral and charged amino acids embedded in a 160-water cluster, where their intramolecular proton transfer energies (neutral amino acid → zwitterionic amino acid) were previously obtained at the ωB97X-D/6-31G(d) level of theory. When the charges on the MM atoms were replaced with those obtained at the same QM level of theory used to treat the QM atoms, this significantly improved the accuracy and convergence of the QM/MM models. In particular, the QM/MM model converged to within 1.4 kcal mol-1 of directly calculated DFT energies for smaller (by as many as 20 waters) QM regions. The use of atomic charges obtained from the natural population analysis yielded the most significant improvement, while other charge schemes such as Mulliken, electrostatic potential, or CM5 led to poorer outcomes. It is further demonstrated that the QM atomic charges can be accurately estimated in a highly efficient manner using an iterative fragmentation approach based on the moving-domain QM/MM method. Similar observations were made when the approach was used to predict the barrier of an SN2 reaction. Thus, the use of QM-quality atomic charges on MM atoms represents a simple and easy-to-implement strategy for improving the accuracy of QM/MM models.
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Jason B Harper
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, New South Wales 2052, Australia
| |
Collapse
|
3
|
Abiri A, Lavigne M, Rezaei M, Nikzad S, Zare P, Mergny JL, Rahimi HR. Unlocking G-Quadruplexes as Antiviral Targets. Pharmacol Rev 2021; 73:897-923. [PMID: 34045305 DOI: 10.1124/pharmrev.120.000230] [Citation(s) in RCA: 53] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Guanine-rich DNA and RNA sequences can fold into noncanonical nucleic acid structures called G-quadruplexes (G4s). Since the discovery that these structures may act as scaffolds for the binding of specific ligands, G4s aroused the attention of a growing number of scientists. The versatile roles of G4 structures in viral replication, transcription, and translation suggest direct applications in therapy or diagnostics. G4-interacting molecules (proteins or small molecules) may also affect the balance between latent and lytic phases, and increasing evidence reveals that G4s are implicated in generally suppressing viral processes, such as replication, transcription, translation, or reverse transcription. In this review, we focus on the discovery of G4s in viruses and the role of G4 ligands in the antiviral drug discovery process. After assessing the role of viral G4s, we argue that host G4s participate in immune modulation, viral tumorigenesis, cellular pathways involved in virus maturation, and DNA integration of viral genomes, which can be potentially employed for antiviral therapeutics. Furthermore, we scrutinize the impediments and shortcomings in the process of studying G4 ligands and drug discovery. Finally, some unanswered questions regarding viral G4s are highlighted for prospective future projects. SIGNIFICANCE STATEMENT: G-quadruplexes (G4s) are noncanonical nucleic acid structures that have gained increasing recognition during the last few decades. First identified as relevant targets in oncology, their importance in virology is now increasingly clear. A number of G-quadruplex ligands are known: viral transcription and replication are the main targets of these ligands. Both viral and cellular G4s may be targeted; this review embraces the different aspects of G-quadruplexes in both host and viral contexts.
Collapse
Affiliation(s)
- Ardavan Abiri
- Department of Medicinal Chemistry, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (A.A., S.N.); Institut Pasteur, Department of Virology, UMR 3569 CNRS, Paris, France (M.L.); Faculty of Medicine, Kerman University of Medical Sciences, Kerman, Iran (M.R.); Dioscuri Center of Chromatin Biology and Epigenomics, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Warsaw, Poland (P.Z.); Faculty of Medicine, Cardinal Stefan Wyszyński University in Warsaw, Warsaw, Poland (P.Z.); Laboratoire d'Optique et Biosciences, Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, Palaiseau cedex, France (J.-L.M.); Pharmaceutics Research Center, Institute of Neuropharmacology, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.); and Department of Pharmacology and Toxicology, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.)
| | - Marc Lavigne
- Department of Medicinal Chemistry, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (A.A., S.N.); Institut Pasteur, Department of Virology, UMR 3569 CNRS, Paris, France (M.L.); Faculty of Medicine, Kerman University of Medical Sciences, Kerman, Iran (M.R.); Dioscuri Center of Chromatin Biology and Epigenomics, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Warsaw, Poland (P.Z.); Faculty of Medicine, Cardinal Stefan Wyszyński University in Warsaw, Warsaw, Poland (P.Z.); Laboratoire d'Optique et Biosciences, Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, Palaiseau cedex, France (J.-L.M.); Pharmaceutics Research Center, Institute of Neuropharmacology, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.); and Department of Pharmacology and Toxicology, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.)
| | - Masoud Rezaei
- Department of Medicinal Chemistry, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (A.A., S.N.); Institut Pasteur, Department of Virology, UMR 3569 CNRS, Paris, France (M.L.); Faculty of Medicine, Kerman University of Medical Sciences, Kerman, Iran (M.R.); Dioscuri Center of Chromatin Biology and Epigenomics, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Warsaw, Poland (P.Z.); Faculty of Medicine, Cardinal Stefan Wyszyński University in Warsaw, Warsaw, Poland (P.Z.); Laboratoire d'Optique et Biosciences, Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, Palaiseau cedex, France (J.-L.M.); Pharmaceutics Research Center, Institute of Neuropharmacology, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.); and Department of Pharmacology and Toxicology, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.)
| | - Sanaz Nikzad
- Department of Medicinal Chemistry, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (A.A., S.N.); Institut Pasteur, Department of Virology, UMR 3569 CNRS, Paris, France (M.L.); Faculty of Medicine, Kerman University of Medical Sciences, Kerman, Iran (M.R.); Dioscuri Center of Chromatin Biology and Epigenomics, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Warsaw, Poland (P.Z.); Faculty of Medicine, Cardinal Stefan Wyszyński University in Warsaw, Warsaw, Poland (P.Z.); Laboratoire d'Optique et Biosciences, Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, Palaiseau cedex, France (J.-L.M.); Pharmaceutics Research Center, Institute of Neuropharmacology, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.); and Department of Pharmacology and Toxicology, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.)
| | - Peyman Zare
- Department of Medicinal Chemistry, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (A.A., S.N.); Institut Pasteur, Department of Virology, UMR 3569 CNRS, Paris, France (M.L.); Faculty of Medicine, Kerman University of Medical Sciences, Kerman, Iran (M.R.); Dioscuri Center of Chromatin Biology and Epigenomics, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Warsaw, Poland (P.Z.); Faculty of Medicine, Cardinal Stefan Wyszyński University in Warsaw, Warsaw, Poland (P.Z.); Laboratoire d'Optique et Biosciences, Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, Palaiseau cedex, France (J.-L.M.); Pharmaceutics Research Center, Institute of Neuropharmacology, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.); and Department of Pharmacology and Toxicology, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.)
| | - Jean-Louis Mergny
- Department of Medicinal Chemistry, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (A.A., S.N.); Institut Pasteur, Department of Virology, UMR 3569 CNRS, Paris, France (M.L.); Faculty of Medicine, Kerman University of Medical Sciences, Kerman, Iran (M.R.); Dioscuri Center of Chromatin Biology and Epigenomics, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Warsaw, Poland (P.Z.); Faculty of Medicine, Cardinal Stefan Wyszyński University in Warsaw, Warsaw, Poland (P.Z.); Laboratoire d'Optique et Biosciences, Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, Palaiseau cedex, France (J.-L.M.); Pharmaceutics Research Center, Institute of Neuropharmacology, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.); and Department of Pharmacology and Toxicology, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.)
| | - Hamid-Reza Rahimi
- Department of Medicinal Chemistry, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (A.A., S.N.); Institut Pasteur, Department of Virology, UMR 3569 CNRS, Paris, France (M.L.); Faculty of Medicine, Kerman University of Medical Sciences, Kerman, Iran (M.R.); Dioscuri Center of Chromatin Biology and Epigenomics, Nencki Institute of Experimental Biology, Polish Academy of Sciences, Warsaw, Poland (P.Z.); Faculty of Medicine, Cardinal Stefan Wyszyński University in Warsaw, Warsaw, Poland (P.Z.); Laboratoire d'Optique et Biosciences, Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, Palaiseau cedex, France (J.-L.M.); Pharmaceutics Research Center, Institute of Neuropharmacology, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.); and Department of Pharmacology and Toxicology, Faculty of Pharmacy, Kerman University of Medical Sciences, Kerman, Iran (H.-R.R.)
| |
Collapse
|
4
|
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
|
5
|
Sutyak KB, Lee W, Zavalij PV, Gutierrez O, Davis JT. Templating and Catalyzing [2+2] Photocycloaddition in Solution Using a Dynamic G-Quadruplex. Angew Chem Int Ed Engl 2018; 57:17146-17150. [PMID: 30395701 DOI: 10.1002/anie.201811202] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2018] [Indexed: 11/08/2022]
Abstract
We describe a templating/covalent capture strategy that enables photochemical formation of 8 cyclobutanes in one noncovalent assembly. This process was characterized by experiment and quantum mechanical/molecular mechanics (ONIOM) calculations. Thus, KI and 16 units of 5'-cinnamate guanosine form a G-quadruplex where C=C π bonds in neighboring G4 -quartets are separated by 3.3 Å, enabling [2+2] photocycloaddition in solution. This reaction is high-yielding (>90 %), regio- and diastereoselective. Since all components are in dynamic equilibrium this photocycloaddition is catalytic in K+ .
Collapse
Affiliation(s)
- Keith B Sutyak
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| | - Wes Lee
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| | - Peter V Zavalij
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| | - Osvaldo Gutierrez
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| | - Jeffery T Davis
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| |
Collapse
|
6
|
Sutyak KB, Lee W, Zavalij PV, Gutierrez O, Davis JT. Templating and Catalyzing [2+2] Photocycloaddition in Solution Using a Dynamic G‐Quadruplex. Angew Chem Int Ed Engl 2018. [DOI: 10.1002/ange.201811202] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Keith B. Sutyak
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| | - Wes Lee
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| | - Peter V. Zavalij
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| | - Osvaldo Gutierrez
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| | - Jeffery T. Davis
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| |
Collapse
|
7
|
Bartas M, Brázda V, Karlický V, Červeň J, Pečinka P. Bioinformatics analyses and in vitro evidence for five and six stacked G-quadruplex forming sequences. Biochimie 2018; 150:70-75. [PMID: 29733879 DOI: 10.1016/j.biochi.2018.05.002] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2018] [Accepted: 05/02/2018] [Indexed: 02/08/2023]
Abstract
Quadruplexes are noncanonical DNA structures that arise in guanine rich loci and have important biological functions. Classically, quadruplexes contain four stacked intramolecular G-tetrads. Surprisingly, although some algorithms allow searching for longer than 4G tracts for quadruplex formation, these have not yet been systematically studied. Therefore, we analyzed the human genome for sequences that are predicted to adopt stacked intramolecular G-tetrads with greater than four stacks. The data provide evidence for numerous G-quadruplexes that contain five or six stacked intramolecular G-tetrads. These sequences are predominantly found in known gene regulatory regions. Electrophoretic mobility assays and circular dichroism spectroscopy indicate that these sequences form quadruplex structures in vitro under physiological conditions. The localization and in vitro stability of these G-quadruplexes indicate their potentially important roles in gene regulation and their potential for therapeutic applications.
Collapse
Affiliation(s)
- Martin Bartas
- Department of Biology and Ecology/Institute of Environmental Technologies, Faculty of Science, University of Ostrava, Chittussiho 10, Ostrava, CZ, 71000, Czech Republic
| | - Václav Brázda
- Department of Biophysical Chemistry and Molecular Oncology, Institute of Biophysics, The Czech Academy of Sciences, Brno, CZ, 61265, Czech Republic
| | - Václav Karlický
- Department of Physics, Faculty of Science, University of Ostrava, Chittussiho 10, Ostrava, CZ, 71000, Czech Republic; Global Change Research Institute CAS, V. V. I., Bělidla 986/4a, Brno, CZ, 60300, Czech Republic
| | - Jiří Červeň
- Department of Biology and Ecology/Institute of Environmental Technologies, Faculty of Science, University of Ostrava, Chittussiho 10, Ostrava, CZ, 71000, Czech Republic
| | - Petr Pečinka
- Department of Biology and Ecology/Institute of Environmental Technologies, Faculty of Science, University of Ostrava, Chittussiho 10, Ostrava, CZ, 71000, Czech Republic.
| |
Collapse
|
8
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
9
|
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
|
10
|
Wang J, Videla PE, Batista VS. Effects of aligned α-helix peptide dipoles on experimental electrostatic potentials. Protein Sci 2017; 26:1692-1697. [PMID: 28556371 DOI: 10.1002/pro.3204] [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: 03/24/2017] [Revised: 05/22/2017] [Accepted: 05/24/2017] [Indexed: 12/26/2022]
Abstract
Aligned protein α-helix dipoles have been implicated in protein function and structure. The recent breakthroughs in high-resolution electron microscopy (EM) of macromolecules makes it possible to explore fundamental aspects of structural biology at the detailed molecular level. The electrostatic potential (ESP) generated by aligned protein α-helix dipole should be observable in high-resolution EM maps despite the fact that the effect may be partially screened by induced electric fields. Here, we show that aligned backbone dipoles in protein α-helices account for long-range features in the protein ESP functions. Our results are consistent with experimental EM maps and density functional theory calculations, including direct Fourier summation for proper calculation of the ESP due to the nonlocal nature of the ESP function from aligned dipoles and other partial atomic charges.
Collapse
Affiliation(s)
- Jimin Wang
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, Connecticut, 06520-8114
| | - Pablo E Videla
- Department of Chemistry, Yale University, New Haven, Connecticut, 06520-8107
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, Connecticut, 06520-8107
| |
Collapse
|
11
|
Jin X, Zhang JZH, He X. Full QM Calculation of RNA Energy Using Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps Method. J Phys Chem A 2017; 121:2503-2514. [PMID: 28264557 DOI: 10.1021/acs.jpca.7b00859] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
In this study, the electrostatically embedded generalized molecular fractionation with conjugate caps (concaps) method (EE-GMFCC) was employed for efficient linear-scaling quantum mechanical (QM) calculation of total energies of RNAs. In the EE-GMFCC approach, the total energy of RNA is calculated by taking a proper combination of the QM energy of each nucleotide-centric fragment with large caps or small caps (termed EE-GMFCC-LC and EE-GMFCC-SC, respectively) deducted by the energies of concaps. The two-body QM interaction energy between non-neighboring ribonucleotides which are spatially in close contact are also taken into account for the energy calculation. Numerical studies were carried out to calculate the total energies of a number of RNAs using the EE-GMFCC-LC and EE-GMFCC-SC methods at levels of the Hartree-Fock (HF) method, density functional theory (DFT), and second-order many-body perturbation theory (MP2), respectively. The results show that the efficiency of the EE-GMFCC-SC method is about 3 times faster than the EE-GMFCC-LC method with minimal accuracy sacrifice. The EE-GMFCC-SC method is also applied for relative energy calculations of 20 different conformers of two RNA systems using HF and DFT, respectively. Both single-point and relative energy calculations demonstrate that the EE-GMFCC method has deviations from the full system results of only a few kcal/mol.
Collapse
Affiliation(s)
- Xinsheng Jin
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China
| | - John Z H Zhang
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062, China.,Department of Chemistry, New York University , New York, New York 10003, United States
| | - Xiao He
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062, China
| |
Collapse
|
12
|
Mirzaei M, Gülseren O, Hadipour N. DFT explorations of quadrupole coupling constants for planar 5-fluorouracil pairs. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.06.004] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
|
13
|
Askerka M, Ho J, Batista ER, Gascón JA, Batista VS. The MOD-QM/MM Method: Applications to Studies of Photosystem II and DNA G-Quadruplexes. Methods Enzymol 2016; 577:443-81. [PMID: 27498648 PMCID: PMC5304415 DOI: 10.1016/bs.mie.2016.05.021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/25/2024]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) hybrid methods are currently the most powerful computational tools for studies of structure/function relations and catalytic sites embedded in macrobiomolecules (eg, proteins and nucleic acids). QM/MM methodologies are highly efficient since they implement quantum chemistry methods for modeling only the portion of the system involving bond-breaking/forming processes (QM layer), as influenced by the surrounding molecular environment described in terms of molecular mechanics force fields (MM layer). Some of the limitations of QM/MM methods when polarization effects are not explicitly considered include the approximate treatment of electrostatic interactions between QM and MM layers. Here, we review recent advances in the development of computational protocols that allow for rigorous modeling of electrostatic interactions in biomacromolecules and structural refinement, beyond the common limitations of QM/MM hybrid methods. We focus on photosystem II (PSII) with emphasis on the description of the oxygen-evolving complex (OEC) and its high-resolution extended X-ray absorption fine structure spectra (EXAFS) in conjunction with Monte Carlo structural refinement. Furthermore, we review QM/MM structural refinement studies of DNA G4 quadruplexes with embedded monovalent cations and direct comparisons to NMR data.
Collapse
Affiliation(s)
- M Askerka
- Yale University, New Haven, CT, United States
| | - J Ho
- Yale University, New Haven, CT, United States
| | - E R Batista
- Los Alamos National Laboratory, Los Alamos, NM, United States
| | - J A Gascón
- University of Connecticut, Storrs, CT, United States
| | - V S Batista
- Yale University, New Haven, CT, United States.
| |
Collapse
|
14
|
Gkionis K, Kruse H, Šponer J. Derivation of Reliable Geometries in QM Calculations of DNA Structures: Explicit Solvent QM/MM and Restrained Implicit Solvent QM Optimizations of G-Quadruplexes. J Chem Theory Comput 2016; 12:2000-16. [DOI: 10.1021/acs.jctc.5b01025] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Konstantinos Gkionis
- Institute of Biophysics,
Academy of Sciences of the Czech Republic, Královopolská 135, 612
65 Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics,
Academy of Sciences of the Czech Republic, Královopolská 135, 612
65 Brno, 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
|
15
|
Jin X, Zhu T, Zhang JZH, He X. A systematic study on RNA NMR chemical shift calculation based on the automated fragmentation QM/MM approach. RSC Adv 2016. [DOI: 10.1039/c6ra22518g] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
1H, 13C and 15N NMR chemical shift calculations on RNAs were performed using the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach.
Collapse
Affiliation(s)
- Xinsheng Jin
- School of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
- China
| | - Tong Zhu
- School of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
- China
- NYU-ECNU Center for Computational Chemistry
| | - John Z. H. Zhang
- School of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
- China
- NYU-ECNU Center for Computational Chemistry
| | - Xiao He
- School of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
- China
- NYU-ECNU Center for Computational Chemistry
| |
Collapse
|
16
|
Sproviero EM, Newcomer MB, Gascón JA, Batista ER, Brudvig GW, Batista VS. The MoD-QM/MM methodology for structural refinement of photosystem II and other biological macromolecules. PHOTOSYNTHESIS RESEARCH 2009; 102:455-470. [PMID: 19633920 PMCID: PMC2954272 DOI: 10.1007/s11120-009-9467-6] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/18/2008] [Accepted: 06/25/2009] [Indexed: 05/28/2023]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) hybrid methods are currently the most powerful computational tools for studies of structure/function relations and structural refinement of macrobiomolecules (e.g., proteins and nucleic acids). These methods are highly efficient, since they implement quantum chemistry techniques for modeling only the small part of the system (QM layer) that undergoes chemical modifications, charge transfer, etc., under the influence of the surrounding environment. The rest of the system (MM layer) is described in terms of molecular mechanics force fields, assuming that its influence on the QM layer can be roughly decomposed in terms of electrostatic interactions and steric hindrance. Common limitations of QM/MM methods include inaccuracies in the MM force fields, when polarization effects are not explicitly considered, and the approximate treatment of electrostatic interactions at the boundaries between QM and MM layers. This article reviews recent advances in the development of computational protocols that allow for rigorous modeling of electrostatic interactions in extended systems beyond the common limitations of QM/MM hybrid methods. We focus on the moving-domain QM/MM (MoD-QM/MM) methodology that partitions the system into many molecular domains and obtains the electrostatic and structural properties of the whole system from an iterative self-consistent treatment of the constituent molecular fragments. We illustrate the MoD-QM/MM method as applied to the description of photosystem II as well as in conjunction with the application of spectroscopically constrained QM/MM optimization methods, based on high-resolution spectroscopic data (extended X-ray absorption fine structure spectra, and exchange coupling constants).
Collapse
Affiliation(s)
- Eduardo M. Sproviero
- Yale University, Department of Chemistry, P. O. Box 208107, New Haven Connecticut 06520-8107 USA
| | - Michael B. Newcomer
- Yale University, Department of Chemistry, P. O. Box 208107, New Haven Connecticut 06520-8107 USA
| | | | - Enrique R. Batista
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545
| | - Gary W. Brudvig
- Yale University, Department of Chemistry, P. O. Box 208107, New Haven Connecticut 06520-8107 USA
| | - Victor S. Batista
- Yale University, Department of Chemistry, P. O. Box 208107, New Haven Connecticut 06520-8107 USA
| |
Collapse
|