1
|
Kladova OA, Tyugashev TE, Mikushina ES, Soloviev NO, Kuznetsov NA, Novopashina DS, Kuznetsova AA. Human Polβ Natural Polymorphic Variants G118V and R149I Affects Substate Binding and Catalysis. Int J Mol Sci 2023; 24:ijms24065892. [PMID: 36982964 PMCID: PMC10051265 DOI: 10.3390/ijms24065892] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2023] [Revised: 03/17/2023] [Accepted: 03/17/2023] [Indexed: 03/30/2023] Open
Abstract
DNA polymerase β (Polβ) expression is essential for the cell's response to DNA damage that occurs during natural cellular processes. Polβ is considered the main reparative DNA polymerase, whose role is to fill the DNA gaps arising in the base excision repair pathway. Mutations in Polβ can lead to cancer, neurodegenerative diseases, or premature aging. Many single-nucleotide polymorphisms have been identified in the POLB gene, but the consequences of these polymorphisms are not always clear. It is known that some polymorphic variants in the Polβ sequence reduce the efficiency of DNA repair, thereby raising the frequency of mutations in the genome. In the current work, we studied two polymorphic variants (G118V and R149I separately) of human Polβ that affect its DNA-binding region. It was found that each amino acid substitution alters Polβ's affinity for gapped DNA. Each polymorphic variant also weakens its binding affinity for dATP. The G118V variant was found to greatly affect Polβ's ability to fill gapped DNA and slowed the catalytic rate as compared to the wild-type enzyme. Thus, these polymorphic variants seem to decrease the ability of Polβ to maintain base excision repair efficiency.
Collapse
Affiliation(s)
- Olga A Kladova
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
| | - Timofey E Tyugashev
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
| | - Elena S Mikushina
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
| | - Nikita O Soloviev
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
| | - Nikita A Kuznetsov
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
- Department of Natural Sciences, Novosibirsk State University, 630090 Novosibirsk, Russia
| | - Daria S Novopashina
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
| | - Aleksandra A Kuznetsova
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of Russian Academy of Sciences, 630090 Novosibirsk, Russia
| |
Collapse
|
2
|
Schackert F, Biedermann J, Abdolvand S, Minniberger S, Song C, Plested AJR, Carloni P, Sun H. Mechanism of Calcium Permeation in a Glutamate Receptor Ion Channel. J Chem Inf Model 2023; 63:1293-1300. [PMID: 36758214 PMCID: PMC9976283 DOI: 10.1021/acs.jcim.2c01494] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) are neurotransmitter-activated cation channels ubiquitously expressed in vertebrate brains. The regulation of calcium flux through the channel pore by RNA-editing is linked to synaptic plasticity while excessive calcium influx poses a risk for neurodegeneration. Unfortunately, the molecular mechanisms underlying this key process are mostly unknown. Here, we investigated calcium conduction in calcium-permeable AMPAR using Molecular Dynamics (MD) simulations with recently introduced multisite force-field parameters for Ca2+. Our calculations are consistent with experiment and explain the distinct calcium permeability in different RNA-edited forms of GluA2. For one of the identified metal binding sites, multiscale Quantum Mechanics/Molecular Mechanics (QM/MM) simulations further validated the results from MD and revealed small but reproducible charge transfer between the metal ion and its first solvation shell. In addition, the ion occupancy derived from MD simulations independently reproduced the Ca2+ binding profile in an X-ray structure of an NaK channel mimicking the AMPAR selectivity filter. This integrated study comprising X-ray crystallography, multisite MD, and multiscale QM/MM simulations provides unprecedented insights into Ca2+ permeation mechanisms in AMPARs, and paves the way for studying other biological processes in which Ca2+ plays a pivotal role.
Collapse
Affiliation(s)
- Florian
Karl Schackert
- Computational
Biomedicine (IAS-5/INM-9), Forschungszentrum
Jülich GmbH, 52428 Jülich, Germany,Department
of Physics, RWTH Aachen University, 52062 Aachen, Germany
| | - Johann Biedermann
- Institute
of Biology, Cellular Biophysics, Humboldt
Universität zu Berlin, 10115 Berlin, Germany,Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany
| | - Saeid Abdolvand
- Institute
of Biology, Cellular Biophysics, Humboldt
Universität zu Berlin, 10115 Berlin, Germany,Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany
| | - Sonja Minniberger
- Institute
of Biology, Cellular Biophysics, Humboldt
Universität zu Berlin, 10115 Berlin, Germany,Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany
| | - Chen Song
- Center
for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China,Peking-Tsinghua
Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China
| | - Andrew J. R. Plested
- Institute
of Biology, Cellular Biophysics, Humboldt
Universität zu Berlin, 10115 Berlin, Germany,Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany
| | - Paolo Carloni
- Computational
Biomedicine (IAS-5/INM-9), Forschungszentrum
Jülich GmbH, 52428 Jülich, Germany,Department
of Physics, RWTH Aachen University, 52062 Aachen, Germany,
| | - Han Sun
- Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany,Institute
of Chemistry, TU Berlin, Straße des 17 Juni 135, 10623 Berlin, Germany,
| |
Collapse
|
3
|
Rahnamoun A, O'Hearn KA, Kaymak MC, Li Z, Merz KM, Aktulga HM. A Polarizable Cationic Dummy Metal Ion Model. J Phys Chem Lett 2022; 13:5334-5340. [PMID: 35675715 DOI: 10.1021/acs.jpclett.2c01279] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
A novel locally polarizable multisite model based on the original cation dummy atom (CDA) model is described for molecular dynamics simulations of ions in condensed phases. Polarization effects are introduced by the electronegativity equalization model (EEM) method where charges on the metal ion and its dummy atoms can fluctuate to respond to the environment. This model includes explicit polarization and ion-induced interactions and can be coupled with nonpolarizable or polarizable water models, making it more transferable to simpler force fields. This approach allows us to enhance the original fixed charge CDA model where the charge distribution cannot adapt to the local solvent structure. To illustrate the new CDApol model, we examined properties of the Zn2+, Al3+, and Zr4+ ions in aqueous solution. The polarizable model and Lennard-Jones parameters were refined for octahedrally coordinated Zn2+, Al3+, and Zr4+ CDAs to reproduce thermodynamic and geometrical properties. Using this locally polarizable model, we were able to obtain the experimental hydration free energy, ion-oxygen distance, and coordination number coupled with the standard 12-6 Lennard-Jones model. This model can be used in myriad additional applications where local polarization and charge transfer is important.
Collapse
Affiliation(s)
- Ali Rahnamoun
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 South Shaw Lane, East Lansing, Michigan 48824-1322, United States
- Department of Computer Science and Engineering, Michigan State University, 428 South Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Kurt A O'Hearn
- Department of Computer Science and Engineering, Michigan State University, 428 South Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Mehmet Cagri Kaymak
- Department of Computer Science and Engineering, Michigan State University, 428 South Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Zhen Li
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 South Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Kenneth M Merz
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 South Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Hasan Metin Aktulga
- Department of Computer Science and Engineering, Michigan State University, 428 South Shaw Lane, East Lansing, Michigan 48824-1322, United States
| |
Collapse
|
4
|
Piskorz TK, Martí-Centelles V, Young TA, Lusby PJ, Duarte F. Computational Modeling of Supramolecular Metallo-organic Cages-Challenges and Opportunities. ACS Catal 2022; 12:5806-5826. [PMID: 35633896 PMCID: PMC9127791 DOI: 10.1021/acscatal.2c00837] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Revised: 04/09/2022] [Indexed: 01/18/2023]
Abstract
![]()
Self-assembled
metallo-organic
cages have emerged as promising
biomimetic platforms that can encapsulate whole substrates akin to
an enzyme active site. Extensive experimental work has enabled access
to a variety of structures, with a few notable examples showing catalytic
behavior. However, computational investigations of metallo-organic
cages are scarce, not least due to the challenges associated with
their modeling and the lack of accurate and efficient protocols to
evaluate these systems. In this review, we discuss key molecular principles
governing the design of functional metallo-organic cages, from the
assembly of building blocks through binding and catalysis. For each
of these processes, computational protocols will be reviewed, considering
their inherent strengths and weaknesses. We will demonstrate that
while each approach may have its own specific pitfalls, they can be
a powerful tool for rationalizing experimental observables and to
guide synthetic efforts. To illustrate this point, we present several
examples where modeling has helped to elucidate fundamental principles
behind molecular recognition and reactivity. We highlight the importance
of combining computational and experimental efforts to speed up supramolecular
catalyst design while reducing time and resources.
Collapse
Affiliation(s)
- Tomasz K. Piskorz
- Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford OX1 3TA, United Kingdom
| | - Vicente Martí-Centelles
- Instituto Interuniversitario de Investigación de Reconocimiento Molecular y Desarrollo Tecnológico (IDM), Universitat Politècnica de València, Universitat de València, Valencia 46022, Spain
| | - Tom A. Young
- Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford OX1 3TA, United Kingdom
| | - Paul J. Lusby
- EaStCHEM School of Chemistry, University of Edinburgh, Joseph Black Building, David Brewster Road, Edinburgh, Scotland EH9 3FJ, United Kingdom
| | - Fernanda Duarte
- Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford OX1 3TA, United Kingdom
| |
Collapse
|
5
|
Modeling glycosaminoglycan–protein complexes. Curr Opin Struct Biol 2022; 73:102332. [DOI: 10.1016/j.sbi.2022.102332] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Revised: 12/20/2021] [Accepted: 01/06/2022] [Indexed: 12/23/2022]
|
6
|
Tanimoto S, Itoh SG, Okumura H. "Bucket brigade" using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2. Biophys J 2021; 120:3615-3627. [PMID: 34339634 PMCID: PMC8324383 DOI: 10.1016/j.bpj.2021.07.026] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Revised: 06/07/2021] [Accepted: 07/27/2021] [Indexed: 01/18/2023] Open
Abstract
The RNA-dependent RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a promising drug target for coronavirus disease 2019 (COVID-19) because it plays the most important role in the replication of the RNA genome. Nucleotide analogs such as remdesivir and favipiravir are thought to interfere with the RNA replication by RdRp. More specifically, they are expected to compete with nucleoside triphosphates, such as ATP. However, the process in which these drug molecules and nucleoside triphosphates are taken up by RdRp remains unknown. In this study, we performed all-atom molecular dynamics simulations to clarify the recognition mechanism of RdRp for these drug molecules and ATP that were at a distance. The ligand recognition ability of RdRp decreased in the order of remdesivir, favipiravir, and ATP. We also identified six recognition paths. Three of them were commonly found in all ligands, and the remaining three paths were ligand-dependent ones. In the common two paths, it was observed that the multiple lysine residues of RdRp carried the ligands to the binding site like a "bucket brigade." In the remaining common path, the ligands directly reached the binding site. Our findings contribute to the understanding of the efficient ligand recognition by RdRp at the atomic level.
Collapse
Affiliation(s)
- Shoichi Tanimoto
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi, Japan
| | - Satoru G Itoh
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi, Japan; Exploratory Research Center on Life and Living Systems, National Institutes of Natural Sciences, Okazaki, Aichi, Japan; Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Aichi, Japan
| | - Hisashi Okumura
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki, Aichi, Japan; Exploratory Research Center on Life and Living Systems, National Institutes of Natural Sciences, Okazaki, Aichi, Japan; Department of Structural Molecular Science, SOKENDAI (The Graduate University for Advanced Studies), Okazaki, Aichi, Japan.
| |
Collapse
|
7
|
Peng J, Zhang Y, Jiang Y, Zhang H. Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References. J Chem Inf Model 2021; 61:2981-2997. [PMID: 34080414 DOI: 10.1021/acs.jcim.1c00281] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg2+ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl2 solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg2+ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl2 solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.
Collapse
Affiliation(s)
- Jiarong Peng
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
8
|
Spinello A, Borišek J, Pavlin M, Janoš P, Magistrato A. Computing Metal-Binding Proteins for Therapeutic Benefit. ChemMedChem 2021; 16:2034-2049. [PMID: 33740297 DOI: 10.1002/cmdc.202100109] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Indexed: 01/18/2023]
Abstract
Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.
Collapse
Affiliation(s)
- Angelo Spinello
- National Research Council of Italy (CNR)-, Institute of Materials (IOM) c/o International School for Advanced Studies (SISSA), via Bonomea 265, 34136, Trieste, Italy
| | - Jure Borišek
- National Institute of Chemistry Institution Hajdrihova ulica 19, 1000, Ljubljana, Slovenia
| | - Matic Pavlin
- Laboratory of Microsensor Structures and Electronics Faculty of Electrical Engineering, University of Ljubljana Tržaška cesta 25, 1000, Ljubljana, Slovenia
| | - Pavel Janoš
- National Research Council of Italy (CNR)-, Institute of Materials (IOM) c/o International School for Advanced Studies (SISSA), via Bonomea 265, 34136, Trieste, Italy
| | - Alessandra Magistrato
- National Research Council of Italy (CNR)-, Institute of Materials (IOM) c/o International School for Advanced Studies (SISSA), via Bonomea 265, 34136, Trieste, Italy
| |
Collapse
|
9
|
Zhang P, Han J, Cieplak P, Cheung MS. Determining the atomic charge of calcium ion requires the information of its coordination geometry in an EF-hand motif. J Chem Phys 2021; 154:124104. [PMID: 33810667 DOI: 10.1063/5.0037517] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
It is challenging to parameterize the force field for calcium ions (Ca2+) in calcium-binding proteins because of their unique coordination chemistry that involves the surrounding atoms required for stability. In this work, we observed a wide variation in Ca2+ binding loop conformations of the Ca2+-binding protein calmodulin, which adopts the most populated ternary structures determined from the molecular dynamics simulations, followed by ab initio quantum mechanical (QM) calculations on all 12 amino acids in the loop that coordinate Ca2+ in aqueous solution. Ca2+ charges were derived by fitting to the electrostatic potential in the context of a classical or polarizable force field (PFF). We discovered that the atomic radius of Ca2+ in conventional force fields is too large for the QM calculation to capture the variation in the coordination geometry of Ca2+ in its ionic form, leading to unphysical charges. Specifically, we found that the fitted atomic charges of Ca2+ in the context of PFF depend on the coordinating geometry of electronegative atoms from the amino acids in the loop. Although nearby water molecules do not influence the atomic charge of Ca2+, they are crucial for compensating for the coordination of Ca2+ due to the conformational flexibility in the EF-hand loop. Our method advances the development of force fields for metal ions and protein binding sites in dynamic environments.
Collapse
Affiliation(s)
- Pengzhi Zhang
- Department of Physics, University of Houston, Houston, Texas 77204, USA
| | - Jaebeom Han
- Department of Physics, University of Houston, Houston, Texas 77204, USA
| | - Piotr Cieplak
- Sanford Burnham Prebys Medical Discovery Institute, La Jolla, California 92037, USA
| | - Margaret S Cheung
- Department of Physics, University of Houston, Houston, Texas 77204, USA
| |
Collapse
|
10
|
Evaluating the Performance of a Non-Bonded Cu 2+ Model Including Jahn-Teller Effect into the Binding of Tyrosinase Inhibitors. Int J Mol Sci 2020; 21:ijms21134783. [PMID: 32640730 PMCID: PMC7369908 DOI: 10.3390/ijms21134783] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Revised: 06/08/2020] [Accepted: 06/09/2020] [Indexed: 02/08/2023] Open
Abstract
Tyrosinase (TYR) is a metalloenzyme classified as a type-3 copper protein, which is involved in the synthesis of melanin through a catalytic process beginning with the conversion of the amino acid l-Tyrosine (l-Tyr) to l-3,4-dihydroxyphenylalanine (l-DOPA). It plays an important role in the mechanism of melanogenesis in various organisms including mammals, plants, and fungi. Herein, we used a combination of computational molecular modeling techniques including molecular dynamic (MD) simulations and the linear interaction energy (LIE) model to evaluate the binding free energy of a set of analogs of kojic acid (KA) in complex with TYR. For the MD simulations, we used a dummy model including the description of the Jahn–Teller effect for Cu2+ ions in the active site of this enzyme. Our results show that the LIE model predicts the TYR binding affinities of the inhibitor in close agreement to experimental results. Overall, we demonstrate that the classical model provides a suitable description of the main interactions between analogs of KA and Cu2+ ions in the active site of TYR.
Collapse
|
11
|
Cannon JF. Novel phosphorylation-dependent regulation in an unstructured protein. Proteins 2019; 88:366-384. [PMID: 31512287 DOI: 10.1002/prot.25812] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2019] [Revised: 07/15/2019] [Accepted: 09/04/2019] [Indexed: 12/15/2022]
Abstract
This work explores how phosphorylation of an unstructured protein region in inhibitor-2 (I2) regulates protein phosphatase-1 (PP1) enzyme activity using molecular dynamics (MD). Free I2 is largely unstructured; however, when bound to PP1, three segments adopt a stable structure. In particular, an I2 helix (i-helix) blocks the PP1 active site and inhibits phosphatase activity. I2 phosphorylation in the PP1-I2 complex activates phosphatase activity without I2 dissociation. The I2 Thr74 regulatory phosphorylation site is in an unstructured domain in PP1-I2. PP1-I2 MD demonstrated that I2 phosphorylation promotes early steps of PP1-I2 activation in explicit solvent models. Moreover, phosphorylation-dependent activation occurred in PP1-I2 complexes derived from I2 orthologs with diverse sequences from human, yeast, worm, and protozoa. This system allowed exploration of features of the 73-residue unstructured human I2 domain critical for phosphorylation-dependent activation. These studies revealed that components of I2 unstructured domain are strategically positioned for phosphorylation responsiveness including a transient α-helix. There was no evidence that electrostatic interactions of I2 phosphothreonine74 influenced PP1-I2 activation. Instead, phosphorylation altered the conformation of residues around Thr74. Phosphorylation uncurled the distance between I2 residues Glu71 to Tyr76 to promote PP1-I2 activation, whereas reduced distances reduced activation. This I2 residue Glu71 to Tyr76 distance distribution, independently from Thr74 phosphorylation, controls I2 i-helix displacement from the PP1 active site leading to PP1-I2 activation.
Collapse
Affiliation(s)
- John F Cannon
- Department of Molecular Microbiology and Immunology, University of Missouri, Columbia, Missouri
| |
Collapse
|
12
|
Maffeo C, Chou HY, Aksimentiev A. Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations. ADVANCED THEORY AND SIMULATIONS 2019; 2:1800191. [PMID: 31728433 PMCID: PMC6855400 DOI: 10.1002/adts.201800191] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2018] [Indexed: 12/15/2022]
Abstract
Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.
Collapse
Affiliation(s)
- Christopher Maffeo
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| | - Han-Yi Chou
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| | - Aleksei Aksimentiev
- Department of Physics, Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign,1110 W Green St, Urbana, IL 61801, USA
| |
Collapse
|
13
|
Zuo Z, Liu J. Assessing the Performance of the Nonbonded Mg 2+ Models in a Two-Metal-Dependent Ribonuclease. J Chem Inf Model 2018; 59:399-408. [PMID: 30521334 DOI: 10.1021/acs.jcim.8b00627] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Magnesium ions (Mg2+), abundant in living cells, are essential for biomolecular structure, dynamics, and function. The biological importance of Mg2+ has motivated continuous development and improvement of various Mg2+ models for molecular dynamics (MD) simulations during the last decades. There are four types of nonbonded Mg2+ models: the point charge models based on a 12-6 or 12-6-4 type Lennard-Jones (LJ) potential, and the multisite models based on a 12-6 or 12-6-4 LJ potential. Here, we systematically assessed the performance of these four types of nonbonded Mg2+ models (21 models in total) in terms of maintaining a challenging intermediate state configuration captured in the structure of a prototypical two-metal-ion RNase H complex with an RNA/DNA hybrid. Our data demonstrate that the 12-6-4 multisite models, which account for charge-induced dipole interactions, perform the best in reproducing all the unique coordination modes in this intermediate state and maintaining the correct carboxylate denticity. Our benchmark work provides a useful guideline for MD simulations and structural refinement of Mg2+-containing biomolecular systems.
Collapse
Affiliation(s)
- Zhicheng Zuo
- Department of Pharmaceutical Sciences , University of North Texas System College of Pharmacy, University of North Texas Health Science Center , Fort Worth , Texas 76107 , United States
| | - Jin Liu
- Department of Pharmaceutical Sciences , University of North Texas System College of Pharmacy, University of North Texas Health Science Center , Fort Worth , Texas 76107 , United States
| |
Collapse
|
14
|
Rinaldi S, Assimon VA, Young ZT, Morra G, Shao H, Taylor IR, Gestwicki JE, Colombo G. A Local Allosteric Network in Heat Shock Protein 70 (Hsp70) Links Inhibitor Binding to Enzyme Activity and Distal Protein-Protein Interactions. ACS Chem Biol 2018; 13:3142-3152. [PMID: 30372610 DOI: 10.1021/acschembio.8b00712] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Allosteric inhibitors can be more difficult to optimize without an understanding of how their binding influences the conformational motions of the target. Here, we used an integrated computational and experimental approach to probe the molecular mechanism of an allosteric inhibitor of heat shock protein 70 (Hsp70). The anticancer compound, MKT-077, is known to bind a conserved site in members of the Hsp70 family, which favors the ADP-bound state and interferes with a protein-protein interaction (PPI) at long range. However, the binding site does not overlap with either the nucleotide-binding cleft or the PPI contact surface, so its mechanism is unclear. To this end, we modeled Hsp70's internal dynamics and studied how MKT-077 alters local sampling of its allosteric states. The results pointed to a set of concerted motions between five loops in Hsp70's nucleotide-binding domain (NBD), surrounding the MKT-077 binding site. To test this prediction, we mutated key residues and monitored chaperone activities in vitro. Together, the results indicate that MKT-077 interacts with loop222 to favor a pseudo-ADP bound conformer of Hsp70's NBD, even when ATP is present. We used this knowledge to synthesize an analog of MKT-077 that would better prevent motions of loop222 and confirmed that it had improved antiproliferative activity in breast cancer cells. These results provide an example of how to unlock and leverage the complex mechanisms of allosteric inhibitors.
Collapse
Affiliation(s)
- Silvia Rinaldi
- Istituto di Chimica del Riconoscimento Molecolare, CNR Via Mario Bianco, 9 20131 Milano, Italy
| | - Victoria A. Assimon
- Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, California 94158, United States
| | - Zapporah T. Young
- Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, California 94158, United States
| | - Giulia Morra
- Istituto di Chimica del Riconoscimento Molecolare, CNR Via Mario Bianco, 9 20131 Milano, Italy
| | - Hao Shao
- Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, California 94158, United States
| | - Isabelle R. Taylor
- Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, California 94158, United States
| | - Jason E. Gestwicki
- Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, California 94158, United States
| | - Giorgio Colombo
- Istituto di Chimica del Riconoscimento Molecolare, CNR Via Mario Bianco, 9 20131 Milano, Italy
- Department of Chemistry, University of Pavia, V.le Taramelli, 12 27100, Pavia, Italy
| |
Collapse
|
15
|
Sala D, Musiani F, Rosato A. Application of Molecular Dynamics to the Investigation of Metalloproteins Involved in Metal Homeostasis. Eur J Inorg Chem 2018. [DOI: 10.1002/ejic.201800602] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Davide Sala
- Magnetic Resonance Center (CERM); University of Florence; Via Luigi Sacconi 6 50019 Sesto Fiorentino Italy
| | - Francesco Musiani
- Laboratory of Bioinorganic Chemistry; Department of Pharmacy and Biotechnology; University of Bologna; Viale Giuseppe Fanin 40, I 40127 Bologna Italy
| | - Antonio Rosato
- Magnetic Resonance Center (CERM); University of Florence; Via Luigi Sacconi 6 50019 Sesto Fiorentino Italy
- Consorzio Interuniversitario di Risonanze Magnetiche di Metallo Proteine; Via Luigi Sacconi 6 50019 Sesto Fiorentino Italy
- Department of Chemistry; University of Florence; Via della Lastruccia 3 50019 Sesto Fiorentino Italy
| |
Collapse
|
16
|
Esposito C, Wiedmer L, Caflisch A. In Silico Identification of JMJD3 Demethylase Inhibitors. J Chem Inf Model 2018; 58:2151-2163. [DOI: 10.1021/acs.jcim.8b00539] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- C. Esposito
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - L. Wiedmer
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - A. Caflisch
- Department of Biochemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| |
Collapse
|
17
|
Š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: 326] [Impact Index Per Article: 54.3] [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
|
18
|
Riniker S. Fixed-Charge Atomistic Force Fields for Molecular Dynamics Simulations in the Condensed Phase: An Overview. J Chem Inf Model 2018; 58:565-578. [DOI: 10.1021/acs.jcim.8b00042] [Citation(s) in RCA: 106] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
19
|
Liao Q, Pabis A, Strodel B, Kamerlin SCL. Extending the Nonbonded Cationic Dummy Model to Account for Ion-Induced Dipole Interactions. J Phys Chem Lett 2017; 8:5408-5414. [PMID: 29022713 PMCID: PMC5672556 DOI: 10.1021/acs.jpclett.7b02358] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2017] [Accepted: 10/12/2017] [Indexed: 05/28/2023]
Abstract
Modeling metalloproteins often requires classical molecular dynamics (MD) simulations in order to capture their relevant motions, which in turn necessitates reliable descriptions of the metal centers involved. One of the most successful approaches to date is provided by the "cationic dummy model", where the positive charge of the metal ion is transferred toward dummy particles that are bonded to the central metal ion in a predefined coordination geometry. While this approach allows for ligand exchange, and captures the correct electrostatics as demonstrated for different divalent metal ions, current dummy models neglect ion-induced dipole interactions. In the present work, we resolve this weakness by taking advantage of the recently introduced 12-6-4 type Lennard-Jones potential to include ion-induced dipole interactions. We revise our previous dummy model for Mg2+ and demonstrate that the resulting model can simultaneously reproduce the experimental solvation free energy and metal-ligand distances without the need for artificial restraints or bonds. As ion-induced dipole interactions become particularly important for highly charged metal ions, we develop dummy models for the biologically relevant ions Al3+, Fe3+, and Cr3+. Finally, the effectiveness of our new models is demonstrated in MD simulations of several diverse (and highly challenging to simulate) metalloproteins.
Collapse
Affiliation(s)
- Qinghua Liao
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC
Box 596, Uppsala 75124, Sweden
| | - Anna Pabis
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC
Box 596, Uppsala 75124, Sweden
| | - Birgit Strodel
- Institute
of Complex Systems: Structural Biochemistry (ICS-6), Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
- Institute
of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, 40204 Düsseldorf, Germany
| | - Shina Caroline Lynn Kamerlin
- Science
for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC
Box 596, Uppsala 75124, Sweden
| |
Collapse
|
20
|
Biswal D, Kusalik PG. Molecular simulations of self-assembly processes in metal-organic frameworks: Model dependence. J Chem Phys 2017; 147:044702. [PMID: 28764378 DOI: 10.1063/1.4994700] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Molecular simulation is a powerful tool for investigating microscopic behavior in various chemical systems, where the use of suitable models is critical to successfully reproduce the structural and dynamic properties of the real systems of interest. In this context, molecular dynamics simulation studies of self-assembly processes in metal-organic frameworks (MOFs), a well-known class of porous materials with interesting chemical and physical properties, are relatively challenging, where a reasonably accurate representation of metal-ligand interactions is anticipated to play an important role. In the current study, we both investigate the performance of some existing models and introduce and test new models to help explore the self-assembly in an archetypal Zn-carboxylate MOF system. To this end, the behavior of six different Zn-ion models, three solvent models, and two ligand models was examined and validated against key experimental structural parameters. To explore longer time scale ordering events during MOF self-assembly via explicit solvent simulations, it is necessary to identify a suitable combination of simplified model components representing metal ions, organic ligands, and solvent molecules. It was observed that an extended cationic dummy atom (ECDA) Zn-ion model combined with an all-atom carboxylate ligand model and a simple dipolar solvent model can reproduce characteristic experimental structures for the archetypal MOF system. The successful use of these models in extensive sets of molecular simulations, which provide key insights into the self-assembly mechanism of this archetypal MOF system occurring during the early stages of this process, has been very recently reported.
Collapse
Affiliation(s)
- Debasmita Biswal
- Department of Chemistry, University of Calgary, 2500 University Dr. NW, Calgary, Alberta T2N 1N4, Canada
| | - Peter G Kusalik
- Department of Chemistry, University of Calgary, 2500 University Dr. NW, Calgary, Alberta T2N 1N4, Canada
| |
Collapse
|
21
|
Masetti M, Musiani F, Bernetti M, Falchi F, Cavalli A, Ciurli S, Recanatini M. Development of a multisite model for Ni(II) ion in solution from thermodynamic and kinetic data. J Comput Chem 2017; 38:1834-1843. [PMID: 28558120 DOI: 10.1002/jcc.24827] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2017] [Revised: 04/12/2017] [Accepted: 04/14/2017] [Indexed: 12/14/2022]
Abstract
Force-field parameters are developed for a multisite model of Ni(II) ions to be used in molecular dynamics simulations combined to enhanced sampling methods. The performances of two charge-partitioning schemes are validated by taking into account structural, thermodynamic, and kinetic observables. One of the two models, featuring partial charges on the dummy atoms only, matches both Ni(II) free energy of solvation and water exchange rates. Such model is particularly suited to study complexation events at a fully dynamic description. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Matteo Masetti
- Laboratory of Computational Medicinal Chemistry, Department of Pharmacy and Biotechnology, Alma Mater Studiorum - Università di Bologna, Via Belmeloro 6, Bologna, I-40126, Italy
| | - Francesco Musiani
- Laboratory of Bioinorganic Chemistry, Department of Pharmacy and Biotechnology, Alma Mater Studiorum - Università di Bologna, Viale G. Fanin 40, Bologna, I-40127, Italy
| | - Mattia Bernetti
- Laboratory of Computational Medicinal Chemistry, Department of Pharmacy and Biotechnology, Alma Mater Studiorum - Università di Bologna, Via Belmeloro 6, Bologna, I-40126, Italy
- Compunet, Istituto Italiano di Tecnologia, Via Morego 30, Genova, I-16163, Italy
| | - Federico Falchi
- Compunet, Istituto Italiano di Tecnologia, Via Morego 30, Genova, I-16163, Italy
| | - Andrea Cavalli
- Laboratory of Computational Medicinal Chemistry, Department of Pharmacy and Biotechnology, Alma Mater Studiorum - Università di Bologna, Via Belmeloro 6, Bologna, I-40126, Italy
- Compunet, Istituto Italiano di Tecnologia, Via Morego 30, Genova, I-16163, Italy
| | - Stefano Ciurli
- Laboratory of Bioinorganic Chemistry, Department of Pharmacy and Biotechnology, Alma Mater Studiorum - Università di Bologna, Viale G. Fanin 40, Bologna, I-40127, Italy
| | - Maurizio Recanatini
- Laboratory of Computational Medicinal Chemistry, Department of Pharmacy and Biotechnology, Alma Mater Studiorum - Università di Bologna, Via Belmeloro 6, Bologna, I-40126, Italy
| |
Collapse
|
22
|
Abstract
Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.
Collapse
Affiliation(s)
| | - Kenneth M. Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, and Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
23
|
Chen J. Functional roles of magnesium binding to extracellular signal-regulated kinase 2 explored by molecular dynamics simulations and principal component analysis. J Biomol Struct Dyn 2017; 36:351-361. [PMID: 28030988 DOI: 10.1080/07391102.2016.1277783] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Molecular dynamics (MD) simulations coupled with principal component (PC) analysis were carried out to study functional roles of Mg2+ binding to extracellular signal-regulated kinase 2 (ERK2). The results suggest that Mg2+ binding heavily decreases eigenvalue of the first principal component and totally inhibits motion strength of ERK2, which favors stabilization of ERK2 structure. Binding free energy predictions indicate that Mg2+ binding produces an important effect on binding ability of adenosine triphosphate (ATP) to ERK2 and strengthens the ATP binding. The calculations of residue-based free energy decomposition show that lack of Mg2+ weakens interactions between the hydrophobic rings of ATP and five residues I29, V37, A50, L105, and L154. Hydrogen bond analyses also prove that Mg2+ binding increases occupancies of hydrogen bonds formed between ATP and residues K52, Q103, D104, and M106. We expect that this study can provide a significant theoretical hint for designs of anticancer drugs targeting ERK2.
Collapse
Affiliation(s)
- Jianzhong Chen
- a School of Science, Shandong Jiaotong University , Jinan 250357 , China
| |
Collapse
|
24
|
Casalino L, Palermo G, Abdurakhmonova N, Rothlisberger U, Magistrato A. Development of Site-Specific Mg(2+)-RNA Force Field Parameters: A Dream or Reality? Guidelines from Combined Molecular Dynamics and Quantum Mechanics Simulations. J Chem Theory Comput 2016; 13:340-352. [PMID: 28001405 DOI: 10.1021/acs.jctc.6b00905] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The vital contribution of Mg2+ ions to RNA biology is challenging to dissect at the experimental level. This calls for the integrative support of atomistic simulations, which at the classical level are plagued by limited accuracy. Indeed, force fields intrinsically neglect nontrivial electronic effects that Mg2+ exerts on its surrounding ligands in varying RNA coordination environments. Here, we present a combined computational study based on classical molecular dynamics (MD) and Density Functional Theory (DFT) calculations, aimed at characterizing (i) the performance of five Mg2+ force field (FF) models in RNA systems and (ii) how charge transfer and polarization affect the binding of Mg2+ ions in different coordination motifs. As a result, a total of ∼2.5 μs MD simulations (100/200 ns for each run) for two prototypical Mg2+-dependent ribozymes showed remarkable differences in terms of populations of inner-sphere coordination site types. Most importantly, complementary DFT calculations unveiled that differences in charge transfer and polarization among recurrent Mg2+-RNA coordination motifs are surprisingly small. In particular, the charge of the Mg2+ ions substantially remains constant through different coordination sites, suggesting that the common philosophy of developing site-specific Mg2+ ion parameters is not in line with the physical origin of the Mg2+-RNA MD simulations inaccuracies. Overall, this study constitutes a guideline for an adept use of current Mg2+ models and provides novel insights for the rational development of next-generation Mg2+ FFs to be employed for atomistic simulations of RNA.
Collapse
Affiliation(s)
- Lorenzo Casalino
- International School for Advanced Studies (SISSA) , Trieste, Italy
| | - Giulia Palermo
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne , CH-1015 Lausanne, Switzerland
| | - Nodira Abdurakhmonova
- International School for Advanced Studies (SISSA) , Trieste, Italy.,Università degli Studi di Trieste , Trieste, Italy
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne , CH-1015 Lausanne, Switzerland
| | - Alessandra Magistrato
- CNR-IOM-Democritos National Simulation Center c/o SISSA , via Bonomea 265, Trieste, Italy
| |
Collapse
|
25
|
Chen J. Clarifying binding difference of ATP and ADP to extracellular signal-regulated kinase 2 by using molecular dynamics simulations. Chem Biol Drug Des 2016; 89:548-558. [PMID: 27696674 DOI: 10.1111/cbdd.12877] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2016] [Revised: 09/03/2016] [Accepted: 09/26/2016] [Indexed: 12/21/2022]
Abstract
Extracellular signal-regulated kinase 2 is a promising target for designs and development of anticancer drugs. Molecular dynamics simulations and molecular mechanics Poisson-Boltzmann method were applied to study binding difference of ADP and ATP to extracellular signal-regulated kinase 2. The results prove that the binding ability of ATP to extracellular signal-regulated kinase 2 is stronger than that of ADP. Principal component analysis performed by using molecular dynamics trajectories suggests that binding of ADP and ATP to extracellular signal-regulated kinase 2 change motion directions of two helices α1 and α2. Residue-based free energy decomposition method was adopted to calculate contributions of separate residues to associations of ADP and ATP with extracellular signal-regulated kinase 2. The results show that ADP and ATP produce strong CH-π interactions with five residues Ile29, Val37, Ala50, Leu105, and Leu154. In addition, five hydrogen bonding interactions of ADP and ATP with residues Lys52, Gln103, Asp104, and Met106 also stabilize bindings of ADP and ATP to extracellular signal-regulated kinase 2. Overall, the CH-π interactions of ATP with five residues Ile29, Val37, Ala50, Leu105, and Leu154 are stronger than ADP. This study is expected to contribute a significant theoretical hint for designs of anticancer drugs targeting extracellular signal-regulated kinase 2.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
26
|
Casalino L, Magistrato A. Structural, dynamical and catalytic interplay between Mg2+ ions and RNA. Vices and virtues of atomistic simulations. Inorganica Chim Acta 2016. [DOI: 10.1016/j.ica.2016.02.011] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
27
|
Yoon H, Warshel A. The control of the discrimination between dNTP and rNTP in DNA and RNA polymerase. Proteins 2016; 84:1616-1624. [PMID: 27480935 DOI: 10.1002/prot.25104] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2016] [Revised: 07/06/2016] [Accepted: 07/11/2016] [Indexed: 12/11/2022]
Abstract
Understanding the origin of discrimination between rNTP and dNTP by DNA/RNA polymerases is important both for gaining fundamental knowledge on the corresponding systems and for advancing the design of specific drugs. This work explores the nature of this discrimination by systematic calculations of the transition state (TS) binding energy in RB69 DNA polymerase (gp43) and T7 RNA polymerase. The calculations reproduce the observed trend, in particular when they included the water contribution obtained by the water flooding approach. Our detailed study confirms the idea that the discrimination is due to the steric interaction between the 2'OH and Tyr416 in DNA polymerase, while the electrostatic interaction is the source of the discrimination in RNA polymerase. Proteins 2016; 84:1616-1624. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Hanwool Yoon
- Department of Chemistry, University of Southern California, Los Angeles, California, 90089-1062
| | - Arieh Warshel
- Department of Chemistry, University of Southern California, Los Angeles, California, 90089-1062.
| |
Collapse
|
28
|
Bergonzo C, Hall KB, Cheatham TE. Divalent Ion Dependent Conformational Changes in an RNA Stem-Loop Observed by Molecular Dynamics. J Chem Theory Comput 2016; 12:3382-9. [PMID: 27294370 PMCID: PMC4944181 DOI: 10.1021/acs.jctc.6b00173] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
We compare the performance of five
magnesium (Mg2+)
ion models in simulations of an RNA stem loop which has an experimentally
determined divalent ion dependent conformational shift. We show that
despite their differences in parametrization and resulting van der
Waals terms, including differences in the functional form of the nonbonded
potential, when the RNA adopts its folded conformation, all models
behave similarly across ten independent microsecond length simulations
with each ion model. However, when the entire structure ensemble is
accounted for, chelation of Mg2+ to RNA is seen in three
of the five models, most egregiously and likely artifactual in simulations
using a 12-6-4 model for the Lennard-Jones potential. Despite the
simple nature of the fixed point-charge and van der Waals sphere models
employed, and with the exception of the likely oversampled directed
chelation of the 12-6-4 potential models, RNA–Mg2+ interactions via first shell water molecules are surprisingly well
described by modern parameters, allowing us to observe the spontaneous
conformational shift from Mg2+ free RNA to Mg2+ associated RNA structure in unrestrained molecular dynamics simulations.
Collapse
Affiliation(s)
- Christina Bergonzo
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah , Salt Lake City, Utah 84112, United States
| | - Kathleen B Hall
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine , St. Louis, Missouri 63110, United States
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, College of Pharmacy, University of Utah , Salt Lake City, Utah 84112, United States
| |
Collapse
|
29
|
Jiang Y, Zhang H, Tan T. Rational Design of Methodology-Independent Metal Parameters Using a Nonbonded Dummy Model. J Chem Theory Comput 2016; 12:3250-60. [DOI: 10.1021/acs.jctc.6b00223] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yang Jiang
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| | - Haiyang Zhang
- Department
of Biological Science and Engineering, School of Chemistry and Biological
Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Tianwei Tan
- Beijing
Key Lab of Bioprocess, College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
30
|
Shaneh A, Purisima EO, Salavati R, Sulea T. Structural Studies of Trypanosoma brucei RNA Editing Ligases and Their Binding Partner Proteins. Biochemistry 2016; 55:2319-31. [PMID: 27031688 DOI: 10.1021/acs.biochem.5b01257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To study the mechanism of ligating nicked RNA strands, we conducted molecular dynamics simulations of Trypanosoma brucei RNA editing ligases L1 and L2 complexed with double-stranded RNA (dsRNA) fragments. In each resulting model, a Mg(2+) ion coordinates the 5'-PO4 of the nicked nucleotide and the 3'-OH of the terminal nucleotide for a nucleophilic reaction consistent with the postulated step 3 chemistry of the ligation mechanism. Moreover, coordination of the 3'-OH to the Mg(2+) ion may lower its pKa, thereby rendering it a more effective nucleophile as an oxyanion. Thus, Mg(2+) may play a twofold role: bringing the reactants into the proximity of each other and activating the nucleophile. We also conducted solvated interaction energy calculations to explore whether ligation specificities can be correlated to ligase-dsRNA binding affinity changes. The calculated dsRNA binding affinities are stronger for both L1 and L2 when the terminal nucleotide is changed from cytosine to guanine, in line with their experimentally measured ligation specificities. Because the ligation mechanism is also influenced by interactions of the ligase with partner proteins from the editosome subcomplex, we also modeled the structure of the RNA-bound L2 in complex with the oligonucleotide binding (OB) domain of largest editosome interacting protein A1. The resulting L2-dsRNA-A1 model, which is consistent with mutagenesis and binding data recorded to date, provides the first atomic-level glimpse of plausible interactions around the RNA ligation site in the presence of an OB domain presented in-trans to a nucleic acid ligase.
Collapse
Affiliation(s)
- Alireza Shaneh
- Institute of Parasitology, McGill University , 21111 Lakeshore Road, Sainte-Anne-de-Bellevue, H9X 3V9 Quebec, Canada.,McGill Centre for Bioinformatics, McGill University , Bellini Building, 3649 Promenade Sir William Osler, Montreal, H3G 0B1 Quebec, Canada
| | - Enrico O Purisima
- Human Health Therapeutics, National Research Council Canada , 6100 Royalmount Avenue, Montreal, H4P 2R2 Quebec, Canada.,Department of Biochemistry, McGill University , McIntyre Medical Building, 3655 Promenade Sir William Osler, Montreal, H3G 1Y6 Quebec, Canada.,McGill Centre for Bioinformatics, McGill University , Bellini Building, 3649 Promenade Sir William Osler, Montreal, H3G 0B1 Quebec, Canada
| | - Reza Salavati
- Institute of Parasitology, McGill University , 21111 Lakeshore Road, Sainte-Anne-de-Bellevue, H9X 3V9 Quebec, Canada.,Department of Biochemistry, McGill University , McIntyre Medical Building, 3655 Promenade Sir William Osler, Montreal, H3G 1Y6 Quebec, Canada.,McGill Centre for Bioinformatics, McGill University , Bellini Building, 3649 Promenade Sir William Osler, Montreal, H3G 0B1 Quebec, Canada
| | - Traian Sulea
- Institute of Parasitology, McGill University , 21111 Lakeshore Road, Sainte-Anne-de-Bellevue, H9X 3V9 Quebec, Canada.,Human Health Therapeutics, National Research Council Canada , 6100 Royalmount Avenue, Montreal, H4P 2R2 Quebec, Canada
| |
Collapse
|
31
|
Komuro Y, Re S, Kobayashi C, Muneyuki E, Sugita Y. CHARMM Force-Fields with Modified Polyphosphate Parameters Allow Stable Simulation of the ATP-Bound Structure of Ca(2+)-ATPase. J Chem Theory Comput 2015; 10:4133-42. [PMID: 26588553 DOI: 10.1021/ct5004143] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Adenosine triphosphate (ATP) is an indispensable energy source in cells. In a wide variety of biological phenomena like glycolysis, muscle contraction/relaxation, and active ion transport, chemical energy released from ATP hydrolysis is converted to mechanical forces to bring about large-scale conformational changes in proteins. Investigation of structure-function relationships in these proteins by molecular dynamics (MD) simulations requires modeling of ATP in solution and ATP bound to proteins with accurate force-field parameters. In this study, we derived new force-field parameters for the triphosphate moiety of ATP based on the high-precision quantum calculations of methyl triphosphate. We tested our new parameters on membrane-embedded sarcoplasmic reticulum Ca(2+)-ATPase and four soluble proteins. The ATP-bound structure of Ca(2+)-ATPase remains stable during MD simulations, contrary to the outcome in shorter simulations using original parameters. Similar results were obtained with the four ATP-bound soluble proteins. The new force-field parameters were also tested by investigating the range of conformations sampled during replica-exchange MD simulations of ATP in explicit water. Modified parameters allowed a much wider range of conformational sampling compared with the bias toward extended forms with original parameters. A diverse range of structures agrees with the broad distribution of ATP conformations in proteins deposited in the Protein Data Bank. These simulations suggest that the modified parameters will be useful in studies of ATP in solution and of the many ATP-utilizing proteins.
Collapse
Affiliation(s)
- Yasuaki Komuro
- Graduate School of Science and Engineering, Chuo University , 1-13-27, Kasuga, Bunkyo-ku, Tokyo 112-8551, Japan.,RIKEN Theoretical Molecular Science Laboratory , 2-1, Hirosawa, Wako-shi, Saitama 351-0198, Japan.,RIKEN Advanced Institute for Computational Science, International Medical Device Alliance (IMDA) 6F , 1-6-5 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Suyong Re
- RIKEN Theoretical Molecular Science Laboratory , 2-1, Hirosawa, Wako-shi, Saitama 351-0198, Japan
| | - Chigusa Kobayashi
- RIKEN Advanced Institute for Computational Science, International Medical Device Alliance (IMDA) 6F , 1-6-5 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Eiro Muneyuki
- Graduate School of Science and Engineering, Chuo University , 1-13-27, Kasuga, Bunkyo-ku, Tokyo 112-8551, Japan
| | - Yuji Sugita
- RIKEN Theoretical Molecular Science Laboratory , 2-1, Hirosawa, Wako-shi, Saitama 351-0198, Japan.,RIKEN Advanced Institute for Computational Science, International Medical Device Alliance (IMDA) 6F , 1-6-5 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN Quantitative Biology Center, International Medical Device Alliance (IMDA) 6F , 1-6-5 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN iTHES , 2-1, Hirosawa, Wako-shi, Saitama 351-0198, Japan
| |
Collapse
|
32
|
Jiang Y, Zhang H, Feng W, Tan T. Refined Dummy Atom Model of Mg2+ by Simple Parameter Screening Strategy with Revised Experimental Solvation Free Energy. J Chem Inf Model 2015; 55:2575-86. [DOI: 10.1021/acs.jcim.5b00286] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Yang Jiang
- Beijing
Key Lab of Bioprocess,
College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| | - Haiyang Zhang
- Beijing
Key Lab of Bioprocess,
College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| | - Wei Feng
- Beijing
Key Lab of Bioprocess,
College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| | - Tianwei Tan
- Beijing
Key Lab of Bioprocess,
College of Life Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
33
|
Liao Q, Kamerlin SCL, Strodel B. Development and Application of a Nonbonded Cu(2+) Model That Includes the Jahn-Teller Effect. J Phys Chem Lett 2015; 6:2657-62. [PMID: 26167255 PMCID: PMC4493862 DOI: 10.1021/acs.jpclett.5b01122] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2015] [Accepted: 06/18/2015] [Indexed: 05/26/2023]
Abstract
Metal ions are both ubiquitous to and crucial in biology. In classical simulations, they are typically described as simple van der Waals spheres, making it difficult to provide reliable force field descriptions for them. An alternative is given by nonbonded dummy models, in which the central metal atom is surrounded by dummy particles that each carry a partial charge. While such dummy models already exist for other metal ions, none is available yet for Cu(2+) because of the challenge to reproduce the Jahn-Teller distortion. This challenge is addressed in the current study, where, for the first time, a dummy model including a Jahn-Teller effect is developed for Cu(2+). We successfully validate its usefulness by studying metal binding in two biological systems: the amyloid-β peptide and the mixed-metal enzyme superoxide dismutase. We believe that our parameters will be of significant value for the computational study of Cu(2+)-dependent biological systems using classical models.
Collapse
Affiliation(s)
- Qinghua Liao
- Institute of Complex
Systems: Structural Biochemistry, Forschungszentrum Jülich, 52425 Jülich, Germany
| | | | - Birgit Strodel
- Institute of Complex
Systems: Structural Biochemistry, Forschungszentrum Jülich, 52425 Jülich, Germany
- Institute
of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Universitätsstrasse 1, 40225 Düsseldorf, Germany
| |
Collapse
|
34
|
Sgrignani J, Magistrato A. QM/MM MD Simulations on the Enzymatic Pathway of the Human Flap Endonuclease (hFEN1) Elucidating Common Cleavage Pathways to RNase H Enzymes. ACS Catal 2015. [DOI: 10.1021/acscatal.5b00178] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Jacopo Sgrignani
- Institute of Research in Biomedicine (IRB), Via Vincenzo Vela, 6500 Bellinzona, Switzerland
| | - Alessandra Magistrato
- CNR-IOM-Democritos
National Simulation Center c/o International School for Advanced Studies (SISSA/ISAS), Via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
35
|
Panteva MT, GiambaȈsu GM, York DM. Comparison of structural, thermodynamic, kinetic and mass transport properties of Mg(2+) ion models commonly used in biomolecular simulations. J Comput Chem 2015; 36:970-82. [PMID: 25736394 PMCID: PMC4409555 DOI: 10.1002/jcc.23881] [Citation(s) in RCA: 66] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2014] [Revised: 01/17/2015] [Accepted: 02/08/2015] [Indexed: 01/09/2023]
Abstract
The prevalence of Mg(2+) ions in biology and their essential role in nucleic acid structure and function has motivated the development of various Mg(2+) ion models for use in molecular simulations. Currently, the most widely used models in biomolecular simulations represent a nonbonded metal ion as an ion-centered point charge surrounded by a nonelectrostatic pairwise potential that takes into account dispersion interactions and exchange effects that give rise to the ion's excluded volume. One strategy toward developing improved models for biomolecular simulations is to first identify a Mg(2+) model that is consistent with the simulation force fields that closely reproduces a range of properties in aqueous solution, and then, in a second step, balance the ion-water and ion-solute interactions by tuning parameters in a pairwise fashion where necessary. The present work addresses the first step in which we compare 17 different nonbonded single-site Mg(2+) ion models with respect to their ability to simultaneously reproduce structural, thermodynamic, kinetic and mass transport properties in aqueous solution. None of the models based on a 12-6 nonelectrostatic nonbonded potential was able to reproduce the experimental radial distribution function, solvation free energy, exchange barrier and diffusion constant. The models based on a 12-6-4 potential offered improvement, and one model in particular, in conjunction with the SPC/E water model, performed exceptionally well for all properties. The results reported here establish useful benchmark calculations for Mg(2+) ion models that provide insight into the origin of the behavior in aqueous solution, and may aid in the development of next-generation models that target specific binding sites in biomolecules.
Collapse
Affiliation(s)
- Maria T. Panteva
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ 08854-8076, USA
| | - George M. GiambaȈsu
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ 08854-8076, USA
| | - Darrin M. York
- Center for Integrative Proteomics Research, BioMaPS Institute and Department of Chemistry & Chemical Biology, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ 08854-8076, USA
| |
Collapse
|
36
|
Bora RP, Mills MJL, Frushicheva MP, Warshel A. On the Challenge of Exploring the Evolutionary Trajectory from Phosphotriesterase to Arylesterase Using Computer Simulations. J Phys Chem B 2015; 119:3434-45. [DOI: 10.1021/jp5124025] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Ram Prasad Bora
- Department
of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, United States
| | - Matthew J. L. Mills
- Deconstruction Division, Joint BioEnergy Institute, Emeryville, California, United States & Biomass Science and Conversion Technology Department, Sandia National Laboratories, Livermore, California 94550, United States
| | - Maria P. Frushicheva
- Department
of Chemical Engineering, Massachusetts Institute of Technology, 25 Ames
Street, Cambridge, Massachusetts 02142, United States
| | - Arieh Warshel
- Department
of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089, United States
| |
Collapse
|
37
|
Zhang R, Bhattacharjee A, Field MJ, Salahub DR. Multiple proton relay routes in the reaction mechanism of RNAP II: Assessing the effect of structural model. Proteins 2014; 83:268-81. [DOI: 10.1002/prot.24732] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2014] [Revised: 11/05/2014] [Accepted: 11/06/2014] [Indexed: 12/22/2022]
Affiliation(s)
- Rui Zhang
- Department of Chemistry; Centre for Molecular Simulation, Institute for Quantum Science and Technology, University of Calgary; Calgary Canada
| | - Anirban Bhattacharjee
- Department of Chemistry; Centre for Molecular Simulation, Institute for Quantum Science and Technology, University of Calgary; Calgary Canada
| | - Martin J. Field
- DYNAMOP; Institut de Biologie Structurale, Jean-Pierre Ebel; Grenoble France
| | - Dennis R. Salahub
- Department of Chemistry; Centre for Molecular Simulation, Institute for Quantum Science and Technology, University of Calgary; Calgary Canada
| |
Collapse
|
38
|
Miller BR, Parish CA, Wu EY. Molecular dynamics study of the opening mechanism for DNA polymerase I. PLoS Comput Biol 2014; 10:e1003961. [PMID: 25474643 PMCID: PMC4256020 DOI: 10.1371/journal.pcbi.1003961] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2014] [Accepted: 10/02/2014] [Indexed: 11/18/2022] Open
Abstract
During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme:DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.
Collapse
Affiliation(s)
- Bill R. Miller
- Department of Biology, University of Richmond, Richmond, Virginia, United States of America
- Department of Chemistry, University of Richmond, Richmond, Virginia, United States of America
| | - Carol A. Parish
- Department of Chemistry, University of Richmond, Richmond, Virginia, United States of America
| | - Eugene Y. Wu
- Department of Biology, University of Richmond, Richmond, Virginia, United States of America
| |
Collapse
|
39
|
Vidossich P, Magistrato A. QM/MM molecular dynamics studies of metal binding proteins. Biomolecules 2014; 4:616-45. [PMID: 25006697 PMCID: PMC4192665 DOI: 10.3390/biom4030616] [Citation(s) in RCA: 64] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2014] [Revised: 06/05/2014] [Accepted: 06/06/2014] [Indexed: 11/16/2022] Open
Abstract
Mixed quantum-classical (quantum mechanical/molecular mechanical (QM/MM)) simulations have strongly contributed to providing insights into the understanding of several structural and mechanistic aspects of biological molecules. They played a particularly important role in metal binding proteins, where the electronic effects of transition metals have to be explicitly taken into account for the correct representation of the underlying biochemical process. In this review, after a brief description of the basic concepts of the QM/MM method, we provide an overview of its capabilities using selected examples taken from our work. Specifically, we will focus on heme peroxidases, metallo-β-lactamases, α-synuclein and ligase ribozymes to show how this approach is capable of describing the catalytic and/or structural role played by transition (Fe, Zn or Cu) and main group (Mg) metals. Applications will reveal how metal ions influence the formation and reduction of high redox intermediates in catalytic cycles and enhance drug metabolism, amyloidogenic aggregate formation and nucleic acid synthesis. In turn, it will become manifest that the protein frame directs and modulates the properties and reactivity of the metal ions.
Collapse
Affiliation(s)
- Pietro Vidossich
- Department of Chemistry, Autonomous University of Barcelona, 08193 Cerdanyola del Vallés, Spain.
| | - Alessandra Magistrato
- CNR-IOM-Democritos National Simulation Center c/o, International School for Advanced Studies (SISSA/ISAS), via Bonomea 265, 34165 Trieste, Italy.
| |
Collapse
|
40
|
Fujiwara SI, Sawada K, Amisaki T. Molecular dynamics study on conformational differences between dGMP and 8-oxo-dGMP: Effects of metal ions. J Mol Graph Model 2014; 51:158-67. [PMID: 24929814 DOI: 10.1016/j.jmgm.2014.05.007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2013] [Revised: 05/09/2014] [Accepted: 05/22/2014] [Indexed: 01/05/2023]
Abstract
The modified nucleotide base 7,8-dihydro-8-oxo-guanine (8-oxo-G) is one of the major sources of spontaneous mutagenesis. Nucleotide-sanitizing enzymes, such as the MutT homolog-1 (MTH1) and nudix-type motif 5 (NUDT5), selectively remove 8-oxo-G from the cellular pool of nucleotides. Previous studies showed that, although the syn conformation generally predominates in purine nucleotides with a bulky substituent at the 8-position, 8-oxo-dGMP binds to both MTH1 and NUDT5 in the anti conformation. This study was initiated to investigate the possibility that 8-oxo-dGMP itself may adopt the anti conformation. Molecular dynamics simulations of mononucleotides (dGMP, 8-oxo-dGMP) in aqueous solution were performed. 8-oxo-dGMP adopted the anti conformation as well as the syn conformation, and the proportion of adopting the anti conformation increased in the presence of metal ions. When 8-oxo-dGMP was in the anti conformation, a metal ion was located between the oxygen atom of phosphate and the oxygen atom at the 8-position of 8-oxo-G. The types of stable anti conformations of 8-oxo-dGMP differed, depending on the ionic radii and charges of coexisting ions. These data suggested a role for metal ions, other than as cofactors for the hydrolysis of the di- and tri-phosphate forms of mononucleotides; that the metal ions help retain the anti conformation of the N-glycosidic torsion angle of 8-oxo-dGMP to promote the binding between the 8-oxo-G deoxynucleotide and the nucleotide-sanitizing enzymes.
Collapse
Affiliation(s)
- Shin-Ichi Fujiwara
- Department of Biological Regulation, Faculty of Medicine, Tottori University, 86 Nishi-cho, Yonago 683-8503, Japan.
| | - Kenichiro Sawada
- Department of Biological Regulation, Faculty of Medicine, Tottori University, 86 Nishi-cho, Yonago 683-8503, Japan
| | - Takashi Amisaki
- Department of Biological Regulation, Faculty of Medicine, Tottori University, 86 Nishi-cho, Yonago 683-8503, Japan
| |
Collapse
|
41
|
Duarte F, Bauer P, Barrozo A, Amrein BA, Purg M, Aqvist J, Kamerlin SCL. Force field independent metal parameters using a nonbonded dummy model. J Phys Chem B 2014; 118:4351-62. [PMID: 24670003 PMCID: PMC4180081 DOI: 10.1021/jp501737x] [Citation(s) in RCA: 116] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
![]()
The cationic dummy atom approach
provides a powerful nonbonded
description for a range of alkaline-earth and transition-metal centers,
capturing both structural and electrostatic effects. In this work
we refine existing literature parameters for octahedrally coordinated
Mn2+, Zn2+, Mg2+, and Ca2+, as well as providing new parameters for Ni2+, Co2+, and Fe2+. In all the cases, we are able to reproduce
both M2+–O distances and experimental solvation
free energies, which has not been achieved to date for transition
metals using any other model. The parameters have also been tested
using two different water models and show consistent performance.
Therefore, our parameters are easily transferable to any force field
that describes nonbonded interactions using Coulomb and Lennard-Jones
potentials. Finally, we demonstrate the stability of our parameters
in both the human and Escherichia coli variants of
the enzyme glyoxalase I as showcase systems, as both enzymes are active
with a range of transition metals. The parameters presented in this
work provide a valuable resource for the molecular simulation community,
as they extend the range of metal ions that can be studied using classical
approaches, while also providing a starting point for subsequent parametrization
of new metal centers.
Collapse
Affiliation(s)
- Fernanda Duarte
- Department of Cell and Molecular Biology, Uppsala University , BMC Box 596, S-751 24 Uppsala, Sweden
| | | | | | | | | | | | | |
Collapse
|
42
|
Duarte F, Amrein BA, Kamerlin SCL. Modeling catalytic promiscuity in the alkaline phosphatase superfamily. Phys Chem Chem Phys 2013; 15:11160-77. [PMID: 23728154 PMCID: PMC3693508 DOI: 10.1039/c3cp51179k] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2013] [Accepted: 05/02/2013] [Indexed: 12/19/2022]
Abstract
In recent years, it has become increasingly clear that promiscuity plays a key role in the evolution of new enzyme function. This finding has helped to elucidate fundamental aspects of molecular evolution. While there has been extensive experimental work on enzyme promiscuity, computational modeling of the chemical details of such promiscuity has traditionally fallen behind the advances in experimental studies, not least due to the nearly prohibitive computational cost involved in examining multiple substrates with multiple potential mechanisms and binding modes in atomic detail with a reasonable degree of accuracy. However, recent advances in both computational methodologies and power have allowed us to reach a stage in the field where we can start to overcome this problem, and molecular simulations can now provide accurate and efficient descriptions of complex biological systems with substantially less computational cost. This has led to significant advances in our understanding of enzyme function and evolution in a broader sense. Here, we will discuss currently available computational approaches that can allow us to probe the underlying molecular basis for enzyme specificity and selectivity, discussing the inherent strengths and weaknesses of each approach. As a case study, we will discuss recent computational work on different members of the alkaline phosphatase superfamily (AP) using a range of different approaches, showing the complementary insights they have provided. We have selected this particular superfamily, as it poses a number of significant challenges for theory, ranging from the complexity of the actual reaction mechanisms involved to the reliable modeling of the catalytic metal centers, as well as the very large system sizes. We will demonstrate that, through current advances in methodologies, computational tools can provide significant insight into the molecular basis for catalytic promiscuity, and, therefore, in turn, the mechanisms of protein functional evolution.
Collapse
Affiliation(s)
- Fernanda Duarte
- Uppsala University, Science for Life Laboratory (SciLifeLab), Cell and Molecular Biology, Uppsala, Sweden. ; ;
| | - Beat Anton Amrein
- Uppsala University, Science for Life Laboratory (SciLifeLab), Cell and Molecular Biology, Uppsala, Sweden. ; ;
| | | |
Collapse
|
43
|
Saxena A, Sept D. Multisite Ion Models That Improve Coordination and Free Energy Calculations in Molecular Dynamics Simulations. J Chem Theory Comput 2013; 9:3538-42. [DOI: 10.1021/ct400177g] [Citation(s) in RCA: 71] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Akansha Saxena
- Department
of Biomedical Engineering and Center for Computational Medicine and
Bioinformatics, University of Michigan,
Ann Arbor, Michigan 48109, United States
- Department of Physics, Applied Physics and Astronomy, Rensselaer Polytechnic Institute, Troy, New York 12180,
United States
| | - David Sept
- Department
of Biomedical Engineering and Center for Computational Medicine and
Bioinformatics, University of Michigan,
Ann Arbor, Michigan 48109, United States
| |
Collapse
|
44
|
Maláč K, Barvík I. Complex between Human RNase HI and the phosphonate-DNA/RNA duplex: Molecular dynamics study. J Mol Graph Model 2013; 44:81-90. [DOI: 10.1016/j.jmgm.2013.05.004] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2012] [Revised: 04/30/2013] [Accepted: 05/05/2013] [Indexed: 10/26/2022]
|
45
|
Maláč K, Barvík I. Substrate recognition by norovirus polymerase: microsecond molecular dynamics study. J Comput Aided Mol Des 2013; 27:373-88. [PMID: 23619980 DOI: 10.1007/s10822-013-9652-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2012] [Accepted: 04/18/2013] [Indexed: 01/02/2023]
Abstract
Molecular dynamics simulations of complexes between Norwalk virus RNA dependent RNA polymerase and its natural CTP and 2dCTP (both containing the O5'-C5'-C4'-O4' sequence of atoms bridging the triphosphate and sugar moiety) or modified coCTP (C5'-O5'-C4'-O4'), cocCTP (C5'-O5'-C4'-C4'') substrates were produced by means of CUDA programmable graphical processing units and the ACEMD software package. It enabled us to gain microsecond MD trajectories clearly showing that similar nucleoside triphosphates can bind surprisingly differently into the active site of the Norwalk virus RNA dependent RNA polymerase. It corresponds to their different modes of action (CTP-substrate, 2dCTP-poor substrate, coCTP-chain terminator, cocCTP-inhibitor). Moreover, extremely rare events-as repetitive pervasion of Arg182 into a potentially reaction promoting arrangement-were captured.
Collapse
Affiliation(s)
- Kamil Maláč
- Faculty of Mathematics and Physics, Institute of Physics, Charles University, Ke Karlovu 5, Prague 2, 121 16, Czech Republic
| | | |
Collapse
|
46
|
Musiani F, Ippoliti E, Micheletti C, Carloni P, Ciurli S. Conformational fluctuations of UreG, an intrinsically disordered enzyme. Biochemistry 2013; 52:2949-54. [PMID: 23560717 DOI: 10.1021/bi4001744] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
UreG proteins are small GTP binding (G) proteins that catalyze the hydrolysis of GTP necessary for the maturation of urease, a virulence factor in bacterial pathogenesis. UreG proteins are the first documented cases of intrinsically disordered enzymes. The comprehension of the dynamics of folding-unfolding events occurring in this protein could shed light on the enzymatic mechanism of UreG. Here, we used the recently developed replica exchange with solute tempering (REST2) computational methodology to explore the conformational space of UreG from Helicobacter pylori (HpUreG) and to identify its structural fluctuations. The same simulation and analysis protocol has been applied to HypB from Methanocaldococcus jannaschii (MjHypB), which is closely related to UreG in both sequence and function, even though it is not intrinsically disordered. A comparison of the two systems reveals that both HpUreG and MjHypB feature a substantial rigidity of the protein regions involved in catalysis, justifying its residual catalytic activity. On the other hand, HpUreG tends to unfold more than MjHypB in portions involved in protein-protein interactions with metallochaperones necessary for the formation of multiprotein complexes known to be involved in urease activation.
Collapse
Affiliation(s)
- Francesco Musiani
- Laboratory of Bioinorganic Chemistry, Department of Pharmacy and Biotechnology, University of Bologna, Bologna, Italy
| | | | | | | | | |
Collapse
|
47
|
Álvarez R, Puebla P, Díaz JF, Bento AC, García-Navas R, de la Iglesia-Vicente J, Mollinedo F, Andreu JM, Medarde M, Peláez R. Endowing Indole-Based Tubulin Inhibitors with an Anchor for Derivatization: Highly Potent 3-Substituted Indolephenstatins and Indoleisocombretastatins. J Med Chem 2013; 56:2813-27. [DOI: 10.1021/jm3015603] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Raquel Álvarez
- Laboratorio de Química
Orgánica y Farmacéutica, CIETUS and IBSAL, Facultad
de Farmacia, Universidad de Salamanca,
Campus Miguel de Unamuno, E-37007 Salamanca, Spain
| | - Pilar Puebla
- Laboratorio de Química
Orgánica y Farmacéutica, CIETUS and IBSAL, Facultad
de Farmacia, Universidad de Salamanca,
Campus Miguel de Unamuno, E-37007 Salamanca, Spain
| | | | - Ana C. Bento
- Instituto de Biología Molecular
y Celular del Cáncer, Centro de Investigación del Cáncer, CSIC-Universidad de Salamanca, Campus Miguel de Unamuno,
E-37007 Salamanca, Spain
| | - Rósula García-Navas
- Instituto de Biología Molecular
y Celular del Cáncer, Centro de Investigación del Cáncer, CSIC-Universidad de Salamanca, Campus Miguel de Unamuno,
E-37007 Salamanca, Spain
| | - Janis de la Iglesia-Vicente
- Instituto de Biología Molecular
y Celular del Cáncer, Centro de Investigación del Cáncer, CSIC-Universidad de Salamanca, Campus Miguel de Unamuno,
E-37007 Salamanca, Spain
| | - Faustino Mollinedo
- Instituto de Biología Molecular
y Celular del Cáncer, Centro de Investigación del Cáncer, CSIC-Universidad de Salamanca, Campus Miguel de Unamuno,
E-37007 Salamanca, Spain
| | | | - Manuel Medarde
- Laboratorio de Química
Orgánica y Farmacéutica, CIETUS and IBSAL, Facultad
de Farmacia, Universidad de Salamanca,
Campus Miguel de Unamuno, E-37007 Salamanca, Spain
| | - Rafael Peláez
- Laboratorio de Química
Orgánica y Farmacéutica, CIETUS and IBSAL, Facultad
de Farmacia, Universidad de Salamanca,
Campus Miguel de Unamuno, E-37007 Salamanca, Spain
| |
Collapse
|
48
|
Lu SY, Huang ZM, Huang WK, Liu XY, Chen YY, Shi T, Zhang J. How calcium inhibits the magnesium-dependent kinase gsk3β: a molecular simulation study. Proteins 2013. [PMID: 23184735 DOI: 10.1002/prot.24221] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Glycogen synthase kinase 3β (GSK3β) is a ubiquitous serine/threonine kinase that plays a pivotal role in many biological processes. GSK3β catalyzes the transfer of γ-phosphate of ATP to the unique substrate Ser/Thr residues with the assistance of two natural activating cofactors Mg(2+). Interestingly, the biological observation reveals that a non-native Ca(2+) ion can inhibit the GSK3β catalytic activity. Here, the inhibitory mechanism of GSK3β by the displacement of native Mg(2+) at site 1 by Ca(2+) was investigated by means of 80 ns comparative molecular dynamics (MD) simulations of the GSK3β···Mg(2+)-2/ATP/Mg(2+) -1 and GSK3β···Mg(2+)-2/ATP/Ca(2+)-1 systems. MD simulation results revealed that using the AMBER point charge model force field for Mg(2+) was more appropriate in the reproduction of the active site architectural characteristics of GSK3β than using the magnesium-cationic dummy atom model force field. Compared with the native Mg(2+) bound system, the misalignment of the critical triphosphate moiety of ATP, the erroneous coordination environments around the Mg(2+) ion at site 2, and the rupture of the key hydrogen bond between the invariant Lys85 and the ATP O(β2) atom in the Ca(2+) substituted system were observed in the MD simulation due to the Ca(2+) ion in active site in order to achieve its preferred sevenfold coordination geometry, which adequately abolish the enzymatic activity. The obtained results are valuable in understanding the possible mechanism by why Ca(2+) inhibits the GSK3β activity and also provide insights into the mechanism of Ca(2+) inhibition in other structurally related protein kinases.
Collapse
Affiliation(s)
- Shao-Yong Lu
- Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai JiaoTong University, School of Medicine, Shanghai 200025, China
| | | | | | | | | | | | | |
Collapse
|
49
|
Abstract
Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.
Collapse
|
50
|
How phosphorylation activates the protein phosphatase-1 • inhibitor-2 complex. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2013; 1834:71-86. [DOI: 10.1016/j.bbapap.2012.09.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2012] [Revised: 08/28/2012] [Accepted: 09/05/2012] [Indexed: 11/16/2022]
|