1
|
Delgado JM, Nagy PR, Varma S. Polarizable AMOEBA Model for Simulating Mg 2+·Protein·Nucleotide Complexes. J Chem Inf Model 2024; 64:378-392. [PMID: 38051630 PMCID: PMC11345861 DOI: 10.1021/acs.jcim.3c01513] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
Molecular mechanics (MM) simulations have the potential to provide detailed insights into the mechanisms of enzymes that utilize nucleotides as cofactors. In most cases, the activities of these enzymes also require the binding of divalent cations to catalytic sites. However, modeling divalent cations in MM simulations has been challenging. The inclusion of explicit polarization was considered promising, but despite improvements over nonpolarizable force fields and despite the inclusion of "Nonbonded-fix (NB-fix)" corrections, errors in interaction energies of divalent cations with proteins remain large. Importantly, the application of these models fails to reproduce the experimental structural data on Mg2+·Protein·ATP complexes. Focusing on these complexes, here we provide a systematic assessment of the polarizable AMOEBA model and recommend critical changes that substantially improve its predictive performance. Our key results are as follows. We first show that our recent revision of the AMOEBA protein model (AMOEBABIO18-HFC), which contains high field corrections (HFCs) to induced dipoles, dramatically improves Mg2+-protein interaction energies, reducing the mean absolute error (MAE) from 17 to 10 kcal/mol. This further supports the general applicability of AMOEBABIO18-HFC. The inclusion of many-body NB-fix corrections further reduces MAE to 6 kcal/mol, which amounts to less than 2% error. The errors are estimated with respect to vdW-inclusive density functional theory that we benchmark against CCSD(T) calculations and experiments. We also present a new model of ATP with revised polarization parameters to better capture its high field response, as well as new vdW and dihedral parameters. The ATP model accurately predicts experimental Mg2+-ATP binding free energy in the aqueous phase and provides new insights into how Mg2+ associates with ATP. Finally, we show that molecular dynamics (MD) simulations of Mg2+·Kinase·ATP complexes carried out with these improvements lead to a better agreement in global and local catalytic site structures between MD and X-ray crystallography.
Collapse
Affiliation(s)
- Julian M Delgado
- Department of Molecular Biosciences, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., Budapest H-1111, Hungary
- HUN-REN-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., Budapest H-1111, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., Budapest H-1111, Hungary
| | - Sameer Varma
- Department of Molecular Biosciences, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
- Department of Physics, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| |
Collapse
|
2
|
Szabó PB, Csóka J, Kállay M, Nagy PR. Linear-Scaling Local Natural Orbital CCSD(T) Approach for Open-Shell Systems: Algorithms, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2023; 19:8166-8188. [PMID: 37921429 PMCID: PMC10687875 DOI: 10.1021/acs.jctc.3c00881] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Revised: 10/05/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023]
Abstract
The extension of the highly optimized local natural orbital (LNO) coupled cluster (CC) with single-, double-, and perturbative triple excitations [LNO-CCSD(T)] method is presented for high-spin open-shell molecules based on restricted open-shell references. The techniques enabling the outstanding efficiency of the closed-shell LNO-CCSD(T) variant are adopted, including the iteration- and redundancy-free second-order Møller-Plesset and (T) formulations as well as the integral-direct, memory- and disk use-economic, and OpenMP-parallel algorithms. For large molecules, the efficiency of our open-shell LNO-CCSD(T) method approaches that of its closed-shell parent method due to the application of restricted orbital sets for demanding integral transformations and a novel approximation for higher-order long-range spin-polarization effects. The accuracy of open-shell LNO-CCSD(T) is extensively tested for radicals and reactions thereof, ionization processes, as well as spin-state splittings, and transition-metal compounds. At the size range where the canonical CCSD(T) reference is accessible (up to 20-30 atoms), the average open-shell LNO-CCSD(T) correlation energies are found to be 99.9 to 99.95% accurate, which translates into average absolute deviations of a few tenths of kcal/mol in the investigated energy differences already with the default settings. For more extensive molecules, the local errors may grow, but they can be estimated and decreased via affordable systematic convergence studies. This enables the accurate modeling of large systems with complex electronic structures, as illustrated on open-shell organic radicals and transition-metal complexes of up to 179 atoms as well as on challenging biochemical systems, including up to 601 atoms and 11,000 basis functions. While the protein models involve difficulties for local approximations, such as the spin states of a bounded iron ion or an extremely delocalized singly occupied orbital, the corresponding single-node LNO-CCSD(T) computations were feasible in a matter of days with 10s to 100 GB of memory use. Therefore, the new LNO-CCSD(T) implementation enables highly accurate computations for open-shell systems of unprecedented size and complexity with widely accessible hardware.
Collapse
Affiliation(s)
- P. Bernát Szabó
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - József Csóka
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| |
Collapse
|
3
|
Belle R, Kamps JJAG, Poater J, Kumar K, Pieters BJGE, Salah E, Claridge TDW, Paton RS, Bickelhaupt FM, Kawamura A, Schofield CJ, Mecinović J. Reading and erasing of the phosphonium analogue of trimethyllysine by epigenetic proteins. Commun Chem 2022; 5:10.1038/s42004-022-00640-4. [PMID: 36071790 PMCID: PMC7613515 DOI: 10.1038/s42004-022-00640-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Accepted: 02/03/2022] [Indexed: 01/27/2023] Open
Abstract
N ε-Methylation of lysine residues in histones plays an essential role in the regulation of eukaryotic transcription. The 'highest' methylation mark, N ε-trimethyllysine, is specifically recognised by N ε-trimethyllysine binding 'reader' domains, and undergoes demethylation, as catalysed by 2-oxoglutarate dependent JmjC oxygenases. We report studies on the recognition of the closest positively charged N ε-trimethyllysine analogue, i.e. its trimethylphosphonium derivative (KPme3), by N ε-trimethyllysine histone binding proteins and Nε-trimethyllysine demethylases. Calorimetric and computational studies with histone binding proteins reveal that H3KP4me3 binds more tightly than the natural H3K4me3 substrate, though the relative differences in binding affinity vary. Studies with JmjC demethylases show that some, but not all, of them can accept the phosphonium analogue of their natural substrates and that the methylation state selectivity can be changed by substitution of nitrogen for phosphorus. The combined results reveal that very subtle changes, e.g. substitution of nitrogen for phosphorus, can substantially affect interactions between ligand and reader domains / demethylases, knowledge that we hope will inspire the development of highly selective small molecules modulating their activity.
Collapse
Affiliation(s)
- Roman Belle
- Department of Chemistry and the Ineos Oxford Institute for Antimicrobial Research, Chemistry Research Laboratory, University of Oxford, 12 Mansfield Road, Oxford, OX1 3TA UK
- Chemistry—School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU UK
| | - Jos J. A. G. Kamps
- Department of Chemistry and the Ineos Oxford Institute for Antimicrobial Research, Chemistry Research Laboratory, University of Oxford, 12 Mansfield Road, Oxford, OX1 3TA UK
- Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Jordi Poater
- ICREA, Pg. Lluís Companys 23, 08010 Barcelona, Spain
- Departament de Química Inorgànica i Orgànica & IQTCUB, Universitat de Barcelona, Martí i Franquès 1-11, 08028 Barcelona, Spain
| | - Kiran Kumar
- Department of Chemistry and the Ineos Oxford Institute for Antimicrobial Research, Chemistry Research Laboratory, University of Oxford, 12 Mansfield Road, Oxford, OX1 3TA UK
| | - Bas J. G. E. Pieters
- Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Eidarus Salah
- Department of Chemistry and the Ineos Oxford Institute for Antimicrobial Research, Chemistry Research Laboratory, University of Oxford, 12 Mansfield Road, Oxford, OX1 3TA UK
| | - Timothy D. W. Claridge
- Department of Chemistry and the Ineos Oxford Institute for Antimicrobial Research, Chemistry Research Laboratory, University of Oxford, 12 Mansfield Road, Oxford, OX1 3TA UK
| | - Robert S. Paton
- Department of Chemistry and the Ineos Oxford Institute for Antimicrobial Research, Chemistry Research Laboratory, University of Oxford, 12 Mansfield Road, Oxford, OX1 3TA UK
| | - F. Matthias Bickelhaupt
- Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
- Department of Theoretical Chemistry, Amsterdam Center for Multiscale Modeling, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Akane Kawamura
- Department of Chemistry and the Ineos Oxford Institute for Antimicrobial Research, Chemistry Research Laboratory, University of Oxford, 12 Mansfield Road, Oxford, OX1 3TA UK
- Chemistry—School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU UK
| | - Christopher J. Schofield
- Department of Chemistry and the Ineos Oxford Institute for Antimicrobial Research, Chemistry Research Laboratory, University of Oxford, 12 Mansfield Road, Oxford, OX1 3TA UK
| | - Jasmin Mecinović
- Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense, Denmark
| |
Collapse
|
4
|
Delgado JA, Wineman-Fisher V, Pandit S, Varma S. Inclusion of High-Field Target Data in AMOEBA's Calibration Improves Predictions of Protein-Ion Interactions. J Chem Inf Model 2022; 62:4713-4726. [PMID: 36173398 DOI: 10.1021/acs.jcim.2c00758] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The reliability of molecular mechanics simulations to predict effects of ion binding to proteins depends on their ability to simultaneously describe ion-protein, ion-water, and protein-water interactions. Force fields (FFs) to describe protein-water and ion-water interactions have been constructed carefully and have also been refined routinely to improve accuracy. Descriptions for ion-protein interactions have also been refined, although in an a posteriori manner through the use of "nonbonded-fix (NB-fix)" approaches in which parameters from default Lennard-Jones mixing rules are replaced with those optimized against some reference data. However, even after NB-fix corrections, there remains a significant need for improvement. This is also true for polarizable FFs that include self-consistent inducible moments. Our recent studies on the polarizable AMOEBA FF suggested that the problem associated with modeling ion-protein interactions could be alleviated by recalibrating polarization models of cation-coordinating functional groups so that they respond better to the high electric fields present near ions. Here, we present such a recalibration of carbonyls, carboxylates, and hydroxyls in the AMOEBA protein FF and report that it does improve predictions substantially─mean absolute errors in Na+-protein and K+-protein interaction energies decrease from 8.7 to 5.3 and 9.6 to 6.3 kcal/mol, respectively. Errors are computed with respect to estimates from van der Waals-inclusive density functional theory benchmarked against high-level quantum mechanical calculations and experiments. While recalibration does improve ion-protein interaction energies, they still remain underestimated, suggesting that further improvements can be made in a systematic manner through modifications in classical formalism. Nevertheless, we show that by applying our many-body NB-fix correction to Lennard-Jones components, these errors are further reduced to 2.7 and 2.6 kcal/mol, respectively, for Na+ and K+ ions. Finally, we show that the recalibrated AMOEBA protein FF retains its intrinsic reliability in predicting protein structure and dynamics in the condensed phase.
Collapse
Affiliation(s)
- Julián A Delgado
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| | - Vered Wineman-Fisher
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| | - Sagar Pandit
- Department of Physics, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| | - Sameer Varma
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States.,Department of Physics, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, United States
| |
Collapse
|
5
|
Theoretical Study of Stability and Electronic Characteristics in Various Complexes of Psoralen as an Anticancer Drug in Gas Phase, Water and CCl4 Solutions. Chem Res Chin Univ 2022. [DOI: 10.1007/s40242-022-1475-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|