1
|
Reilley DJ, Arraf Z, Alexandrova AN. Contrasting Effects of Inhibitors Li + and Be 2+ on Catalytic Cycle of Glycogen Synthase Kinase-3β. J Phys Chem B 2021; 125:9480-9489. [PMID: 34404214 DOI: 10.1021/acs.jpcb.1c05099] [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
Ionic lithium shows rare effectiveness for treating bipolar disorder and is a potential drug for neurodegenerative diseases. Unfortunately, lithium suffers from significant drawbacks, mainly a narrow therapeutic window. Among the targets of lithium, glycogen synthase kinase 3β (GSK-3β) may be responsible for its therapeutic effects. The development of alternative, selective inhibitors of this kinase could prevent lithium side effects, but such efforts have met little success so far. An atomistic understanding of Li+ inhibition and the GSK-3β phosphorylation reaction would therefore facilitate the development of new drugs. In this study, we use extensive sampling of catalytic states with our mixed quantum-classical dynamics method QM/DMD and binding affinities from a competitive metal affinity (CMA) approach to expand the atomistic picture of Li+ GSK-3β inhibition. We compare Li+ action with Be2+ and find our results in agreement with in vitro kinetics studies. Ultimately, our simulations show that Li+ inhibition is driven by decreasing the phosphorylation reaction rate, rather than reducing catalytic turnover through tight binding to different GSK-3β states like Be2+ inhibition. The effect of these metals derive from electrostatic differences and especially their smaller atomic radii compared to the native Mg2+ and thus provide insight for the development of GSK-3β inhibitors based on other paradigms.
Collapse
Affiliation(s)
- David J Reilley
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Zaher Arraf
- Department of Education in Technology and Science, Technion-Israel Institute of Technology, 32000, Haifa, Israel
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States.,California NanoSystems Institute, Los Angeles California 90095-1569, United States
| |
Collapse
|
2
|
Abstract
QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.
Collapse
Affiliation(s)
- Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Tanmoy Pal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Luke Xie
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
3
|
Reilley DJ, Fuller JT, Nechay MR, Victor M, Li W, Ruberry JD, Mujika JI, Lopez X, Alexandrova AN. Toxic and Physiological Metal Uptake and Release by Human Serum Transferrin. Biophys J 2020; 118:2979-2988. [PMID: 32497515 PMCID: PMC7300305 DOI: 10.1016/j.bpj.2020.05.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2020] [Revised: 04/03/2020] [Accepted: 04/29/2020] [Indexed: 10/24/2022] Open
Abstract
An atomistic understanding of metal transport in the human body is critical to anticipate the side effects of metal-based therapeutics and holds promise for new drugs and drug delivery designs. Human serum transferrin (hTF) is a central part of the transport processes because of its ubiquitous ferrying of physiological Fe(III) and other transition metals to tightly controlled parts of the body. There is an atomistic mechanism for the uptake process with Fe(III), but not for the release process, or for other metals. This study provides initial insight into these processes for a range of transition metals-Ti(IV), Co(III), Fe(III), Ga(III), Cr(III), Fe(II), Zn(II)-through fully atomistic, extensive quantum mechanical/discrete molecular dynamics sampling and provides, to our knowledge, a new technique we developed to calculate relative binding affinities between metal cations and the protein. It identifies protonation of Tyr188 as a trigger for metal release rather than protonation of Lys206 or Lys296. The study identifies the difficulty of metal release from hTF as potentially related to cytotoxicity. Simulations identify a few critical interactions that stabilize the metal binding site in a flexible, nuanced manner.
Collapse
Affiliation(s)
- David J Reilley
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California
| | - Jack T Fuller
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California
| | - Michael R Nechay
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California
| | - Marie Victor
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California; Institut Lumire Matire, Villeurbanne, France
| | - Wei Li
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California; Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts
| | - Josiah D Ruberry
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California
| | - Jon I Mujika
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU) and Donostia, International Physics Center, Donostia, Euskadi, Spain
| | - Xabier Lopez
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU) and Donostia, International Physics Center, Donostia, Euskadi, Spain
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California; California NanoSystems Institute, University of California, Los Angeles, Los Angeles, California.
| |
Collapse
|
4
|
Reilley DJ, Hennefarth MR, Alexandrova AN. The Case for Enzymatic Competitive Metal Affinity Methods. ACS Catal 2020; 10:2298-2307. [PMID: 34012720 PMCID: PMC8130888 DOI: 10.1021/acscatal.9b04831] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- David J Reilley
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, CA 90095-1569, USA
| | - Matthew R Hennefarth
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, CA 90095-1569, USA
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, CA 90095-1569, USA
- California NanoSystems Institute, University of California, Los Angeles, 570 Westwood Plaza, Los Angeles, California 90095-1569, USA
| |
Collapse
|
5
|
Quantitative prediction of electronic absorption spectra of copper(II)-bioligand systems: Validation and applications. J Inorg Biochem 2019; 204:110953. [PMID: 31816442 DOI: 10.1016/j.jinorgbio.2019.110953] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2019] [Revised: 11/21/2019] [Accepted: 11/25/2019] [Indexed: 02/06/2023]
Abstract
The visible region of the electronic absorption spectra of Cu(II) complexes was studied by time-dependent density functional theory (TD-DFT). The performance of twelve functionals in the prediction of absorption maxima (λmax) was tested on eleven compounds with different geometry, donors and charge. The ranking of the functionals for λmax was determined in terms of mean absolute percent deviation (MAPD) and standard deviation (SD) and it is as follows: BHandHLYP > M06 ≫ CAM-B3LYP ≫ MPW1PW91 ~ B1LYP ~ BLYP > HSE06 ~ B3LYP > B3P86 ~ ω-B97x-D ≫ TPSSh ≫ M06-2X (MAPD) and BHandHLYP > M06 ~ HSE06 > ω-B97x-D ~ CAM-B3LYP ~ MPW1PW91 > B1LYP ~ B3LYP > B3P86 > BLYP ≫ TPSSh ≫ M06-2X (SD). With BHandHLYP functional the MAPD is 3.1% and SD is 2.3%, while with M06 the MAPD is 3.7% and SD is 3.7%. The protocol validated in the first step of the study was applied to: i) calculate the number of transitions in the spectra and relate them to the geometry of Cu(II) species; ii) determine the coordination of axial water(s); iii) predict the electronic spectra of the systems where Cu(II) is bound to human serum albumin (HSA) and to the regions 94-97 and 108-112 of prion protein (PrP). The results indicate that the proposed computational protocol allows a successful prediction of the electronic spectra of Cu(II) species and to relate an experimental spectrum to a specific structure.
Collapse
|
6
|
Reilley DJ, Popov KI, Dokholyan NV, Alexandrova AN. Uncovered Dynamic Coupling Resolves the Ambiguous Mechanism of Phenylalanine Hydroxylase Oxygen Binding. J Phys Chem B 2019; 123:4534-4539. [PMID: 31038957 DOI: 10.1021/acs.jpcb.9b02893] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Phenylalanine hydroxylase (PAH) is an iron enzyme catalyzing the oxidation of l-Phe to l-Tyr during phenylalanine catabolism. Dysfunction of PAH leads to the debilitating condition phenylketonuria (PKU), which prompted research into the structure and function of PAH over the last 50 years. Despite intensive study, there is no consensus on the atomistic details of the mechanism of O2 binding and splitting by wild-type (WT) PAH and how it varies with PKU-inducing mutations, Arg158Gln and Glu280Lys. We studied structures involved in a proposed mechanism for the WT and mutants using extensive mixed quantum-classical molecular dynamics simulations. Simulations reveal a previously unobserved dynamic coupling between the active site and the mutation sites, suggesting how they can affect the catalytic performance of PAH. Furthermore, the effect of the coupling on the PAH structure agrees with and expands our understanding of the experimentally observed differences in activity between the WT and mutants.
Collapse
Affiliation(s)
- David J Reilley
- Department of Chemistry and Biochemistry , University of California, Los Angeles , Los Angeles , California 90095-1569 , United States
| | - Konstantin I Popov
- Department of Biochemistry and Biophysics , University of North Carolina at Chapel Hill , Chapel Hill , North Carolina 27599 , United States
| | - Nikolay V Dokholyan
- Department of Biochemistry and Biophysics , University of North Carolina at Chapel Hill , Chapel Hill , North Carolina 27599 , United States.,Department of Pharmacology, Department of Biochemistry & Molecular Biology , Penn State University College of Medicine , Hershey , Pennsylvania 17033 , United States
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry , University of California, Los Angeles , Los Angeles , California 90095-1569 , United States.,California NanoSystems Institute, Los Angeles , California 90095-1569 , United States
| |
Collapse
|
7
|
Bignon E, Rizza S, Filomeni G, Papaleo E. Use of Computational Biochemistry for Elucidating Molecular Mechanisms of Nitric Oxide Synthase. Comput Struct Biotechnol J 2019; 17:415-429. [PMID: 30996821 PMCID: PMC6451115 DOI: 10.1016/j.csbj.2019.03.011] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 03/17/2019] [Accepted: 03/21/2019] [Indexed: 12/25/2022] Open
Abstract
Nitric oxide (NO) is an essential signaling molecule in the regulation of multiple cellular processes. It is endogenously synthesized by NO synthase (NOS) as the product of L-arginine oxidation to L-citrulline, requiring NADPH, molecular oxygen, and a pterin cofactor. Two NOS isoforms are constitutively present in cells, nNOS and eNOS, and a third is inducible (iNOS). Despite their biological relevance, the details of their complex structural features and reactivity mechanisms are still unclear. In this review, we summarized the contribution of computational biochemistry to research on NOS molecular mechanisms. We described in detail its use in studying aspects of structure, dynamics and reactivity. We also focus on the numerous outstanding questions in the field that could benefit from more extensive computational investigations.
Collapse
Affiliation(s)
- Emmanuelle Bignon
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Salvatore Rizza
- Redox Signaling and Oxidative Stress Group, Cell Stress and Survival Unit, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Giuseppe Filomeni
- Redox Signaling and Oxidative Stress Group, Cell Stress and Survival Unit, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark.,Department of Biology, University of Rome Tor Vergata, Rome, Italy
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark.,Translational Disease Systems Biology, Faculty of Health and Medical Sciences, Novo Nordisk Foundation Center for Protein Research University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
8
|
Fuller J, Wilson TR, Eberhart ME, Alexandrova AN. Charge Density in Enzyme Active Site as a Descriptor of Electrostatic Preorganization. J Chem Inf Model 2019; 59:2367-2373. [DOI: 10.1021/acs.jcim.8b00958] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jack Fuller
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Tim R. Wilson
- Molecular Theory Group, Colorado School of Mines, Golden, Colorado 80401, United States
| | - Mark E. Eberhart
- Molecular Theory Group, Colorado School of Mines, Golden, Colorado 80401, United States
| | - Anastassia N. Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
- California
NanoSystems
Institute, Los Angeles, California 90095, United States
| |
Collapse
|
9
|
A High Coordination of Cross-Links Is Beneficial for the Strength of Cross-Linked Fibers. Biomimetics (Basel) 2019; 4:biomimetics4010012. [PMID: 31105198 PMCID: PMC6477605 DOI: 10.3390/biomimetics4010012] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2018] [Revised: 01/18/2019] [Accepted: 01/22/2019] [Indexed: 01/23/2023] Open
Abstract
The influence of the coordination of (reversible) cross-links on the mechanical properties of aligned fiber bundles is investigated. Two polymeric systems containing cross-links of different coordination (two- and three-fold coordination) but having the same binding energy are investigated. In particular, the response to loading of these systems is compared. Mechanical parameters (strength, stiffness and work-to-fracture) are obtained by computational loading tests. The influence of coordination is studied for simple test systems with pre-defined topologies that maximize strength as well as for more realistic fiber bundles containing nine chains. The results show that a higher coordination of cross-links has a beneficial effect on the strength and the stiffness of the systems, while the work-to-fracture was found larger for the system having a smaller coordination of cross-links. It can be concluded that controlling the coordination of cross-links is a versatile tool to specifically tailor the mechanical properties of polymeric structures.
Collapse
|
10
|
Simm GN, Vaucher AC, Reiher M. Exploration of Reaction Pathways and Chemical Transformation Networks. J Phys Chem A 2018; 123:385-399. [DOI: 10.1021/acs.jpca.8b10007] [Citation(s) in RCA: 103] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Gregor N. Simm
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Alain C. Vaucher
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
11
|
Morgenstern A, Jaszai M, Eberhart ME, Alexandrova AN. Quantified electrostatic preorganization in enzymes using the geometry of the electron charge density. Chem Sci 2017; 8:5010-5018. [PMID: 28970888 PMCID: PMC5612031 DOI: 10.1039/c7sc01301a] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2017] [Accepted: 04/20/2017] [Indexed: 11/21/2022] Open
Abstract
Electrostatic preorganization is thought to be a principle factor responsible for the impressive catalytic capabilities of enzymes. The full protein structure is believed to facilitate catalysis by exerting a highly specific electrostatic field on the active site. Computationally determining the extent of electrostatic preorganization is a challenging process. We propose using the topology and geometry of the electron charge density in the enzyme's active site to asses the effects of electrostatic preorganization. In support of this approach we study the convergence of features of the charge density as the size of the active site model increases in Histone Deacetylase 8. The magnitude of charge density at critical points and most Bader atomic charges are found to converge quickly as more of the protein is included in the simulation. The exact position of critical points however, is found to converge more slowly and be strongly influenced by the protein residues that are further away from the active site. We conjecture that the positions of critical points are affected through perturbations to the wavefunctions in the active site caused by dipole moments from amino acid residues throughout the protein. We further hypothesize that electrostatic preorganization, from the point of view of charge density, can not be easily understood through the charges on atoms or the nature of the bonding interactions, but through the relative positions of critical points that are known to correlate with reactivity and reaction barriers.
Collapse
Affiliation(s)
| | - Matthew Jaszai
- Molecular Theory Group , Colorado School of Mines , USA .
| | | | | |
Collapse
|