1
|
Ligorio R, Grosskopf P, Dos Santos LHR, Krawczuk A. Unveiling GruPol: Predicting Electric and Electrostatic Properties of Macromolecules via the Building Block Approach. J Phys Chem B 2024; 128:7954-7965. [PMID: 38976348 PMCID: PMC11345817 DOI: 10.1021/acs.jpcb.4c03062] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2024] [Revised: 06/14/2024] [Accepted: 06/20/2024] [Indexed: 07/10/2024]
Abstract
Understanding electrostatics and electric properties of macromolecules is crucial in uncovering the intricacies of their behavior and functionality. The precise knowledge of these properties enhances our ability to manipulate and engineer macromolecules for diverse applications, spanning from drug design to materials science. Having that in mind, we present here the GruPol database approach to characterize and accurately predict dipole moments, static polarizabilities, and electrostatic potential of proteins and their subunits. The method involves partitioning of the electron density, calculated at the M06-HF/aug-cc-pVDZ level of theory, of small peptides into predefined building blocks that are averaged over the database. By manipulating and positioning these building blocks, GruPol enables the description of proteins assembled from over nearly 100 residual entries, allowing for efficient and precise computation of the above-mentioned properties across a broad range of proteins. The database enables the user to include solvent effects as well as define protonation states on the protein's backbone to account for pH variations. The precision of the proposed scheme is benchmarked against experimental data for myoglobin species.
Collapse
Affiliation(s)
- Raphael
F. Ligorio
- Institute
of Inorganic Chemistry, University of Goettingen, Tammannstrasse 4, D-37077 Göttingen, Germany
| | - Paul Grosskopf
- Institute
of Inorganic Chemistry, University of Goettingen, Tammannstrasse 4, D-37077 Göttingen, Germany
| | - Leonardo H. R. Dos Santos
- Departamento
de Química, Instituto de Ciências Exatas, Universidade Federal de Minas Gerais, Av. Antônio Carlos 6627, 31270-901 Belo Horizonte, MG, Brazil
| | - Anna Krawczuk
- Institute
of Inorganic Chemistry, University of Goettingen, Tammannstrasse 4, D-37077 Göttingen, Germany
| |
Collapse
|
2
|
Mocchetti E, Morette L, Mulliert G, Mathiot S, Guillot B, Dehez F, Chauvat F, Cassier-Chauvat C, Brochier-Armanet C, Didierjean C, Hecker A. Biochemical and Structural Characterization of Chi-Class Glutathione Transferases: A Snapshot on the Glutathione Transferase Encoded by sll0067 Gene in the Cyanobacterium Synechocystis sp. Strain PCC 6803. Biomolecules 2022; 12:biom12101466. [PMID: 36291676 PMCID: PMC9599700 DOI: 10.3390/biom12101466] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Revised: 09/30/2022] [Accepted: 10/05/2022] [Indexed: 11/24/2022] Open
Abstract
Glutathione transferases (GSTs) constitute a widespread superfamily of enzymes notably involved in detoxification processes and/or in specialized metabolism. In the cyanobacterium Synechocsytis sp. PCC 6803, SynGSTC1, a chi-class GST (GSTC), is thought to participate in the detoxification process of methylglyoxal, a toxic by-product of cellular metabolism. A comparative genomic analysis showed that GSTCs were present in all orders of cyanobacteria with the exception of the basal order Gloeobacterales. These enzymes were also detected in some marine and freshwater noncyanobacterial bacteria, probably as a result of horizontal gene transfer events. GSTCs were shorter of about 30 residues compared to most cytosolic GSTs and had a well-conserved SRAS motif in the active site (10SRAS13 in SynGSTC1). The crystal structure of SynGSTC1 in complex with glutathione adopted the canonical GST fold with a very open active site because the α4 and α5 helices were exceptionally short. A transferred multipolar electron-density analysis allowed a fine description of the solved structure. Unexpectedly, Ser10 did not have an electrostatic influence on glutathione as usually observed in serinyl-GSTs. The S10A variant was only slightly less efficient than the wild-type and molecular dynamics simulations suggested that S10 was a stabilizer of the protein backbone rather than an anchor site for glutathione.
Collapse
Affiliation(s)
- Eva Mocchetti
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
| | - Laura Morette
- Université de Lorraine, INRAE, IAM, F-54000 Nancy, France
| | | | | | - Benoît Guillot
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
| | - François Dehez
- Université de Lorraine, CNRS, LPCT, F-54000 Nancy, France
| | - Franck Chauvat
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), F-91190 Gif-sur-Yvette, France
| | - Corinne Cassier-Chauvat
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), F-91190 Gif-sur-Yvette, France
| | | | - Claude Didierjean
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
- Correspondence: (C.D.); (A.H.)
| | - Arnaud Hecker
- Université de Lorraine, INRAE, IAM, F-54000 Nancy, France
- Correspondence: (C.D.); (A.H.)
| |
Collapse
|
3
|
Jara-Cortés J, Matta CF, Hernández-Trujillo J. A fast approximate extension of the interacting quantum atoms energy decomposition to excited states. J Comput Chem 2022; 43:1068-1078. [PMID: 35470908 DOI: 10.1002/jcc.26863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Revised: 02/25/2022] [Accepted: 03/23/2022] [Indexed: 11/11/2022]
Abstract
An approach is developed for the fast calculation of the interacting quantum atoms energy decomposition (IQA) from the information contained in the first order reduced density matrix only. The proposed methodology utilizes an approximate exchange-correlation density from Density Matrix Functional Theory without the need to evaluate the correlation-exchange contribution directly. Instead, weight factors are estimated to decompose the exact Vxc into atomic and pairwise contributions. In this way, the sum of the IQA contributions recovers the energy obtained from the electronic structure calculation. This method can, hence, be applied to obtain atomic contributions in excited states on the same footing as in their ground states using any method that delivers the reduced first-order density matrix. In this way, one can locate chromophores from first principles quantum chemical calculations. Test calculations on the ground and excited states of a set of small molecules indicate that the scaled atomic contributions reproduce vertical electronic transition energies calculated exactly. This approach may be useful to extend the applicability of the IQA approach in the study of large photochemical systems especially when the calculations of the second order reduced density matrices is prohibitive or not possible.
Collapse
Affiliation(s)
- Jesús Jara-Cortés
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic, Mexico
| | - Chérif F Matta
- Department of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia, Canada
| | - Jesús Hernández-Trujillo
- Departamento de Física y Química Teórica, Facultad de Química, UNAM. Circuito Escolar, Ciudad Universitaria, Mexico City, Mexico
| |
Collapse
|
4
|
Lefrancois-Gagnon KM, Mawhinney RC. Toward universal substituent constants: Transferability of group descriptors from the quantum theory of atoms in molecules. J Comput Chem 2022; 43:265-278. [PMID: 34842294 DOI: 10.1002/jcc.26787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Revised: 09/29/2021] [Accepted: 11/06/2021] [Indexed: 11/06/2022]
Abstract
Traditionally, substituents are described not by their intrinsic properties, but by their effect elsewhere in a molecule. However, the quantum theory of atoms in molecules (QTAIM) provides a route to intrinsic substituent descriptors. Ideally, these descriptors would exhibit minimal change as the local environment changes, and hence be considered transferable. Whether this is true is an open question. Here, we evaluated the transferability of QTAIM functional group descriptors for 117 functional groups in a series of 17 substrates to determine whether descriptors obtained using hydrogen as substrate are transferable. The functional group volume has a strong, consistent, linear relationship throughout. All other hydrogen-based group descriptors exhibit a relatively strong linear relationship with those in carbon-based substrates and a reasonable linear relationship with those in non-carbon-based substrates. Outliers are readily interpretable in terms of substrate induced functional group geometry changes. As expected, directional descriptors lying along the substituent-substrate axis are the least conserved.
Collapse
|
5
|
Vuković V, Leduc T, Jelić-Matošević Z, Didierjean C, Favier F, Guillot B, Jelsch C. A rush to explore protein-ligand electrostatic interaction energy with Charger. Acta Crystallogr D Struct Biol 2021; 77:1292-1304. [PMID: 34605432 DOI: 10.1107/s2059798321008433] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 08/13/2021] [Indexed: 11/10/2022] Open
Abstract
The mutual penetration of electron densities between two interacting molecules complicates the computation of an accurate electrostatic interaction energy based on a pseudo-atom representation of electron densities. The numerical exact potential and multipole moment (nEP/MM) method is time-consuming since it performs a 3D integration to obtain the electrostatic energy at short interaction distances. Nguyen et al. [(2018), Acta Cryst. A74, 524-536] recently reported a fully analytical computation of the electrostatic interaction energy (aEP/MM). This method performs much faster than nEP/MM (up to two orders of magnitude) and remains highly accurate. A new program library, Charger, contains an implementation of the aEP/MM method. Charger has been incorporated into the MoProViewer software. Benchmark tests on a series of small molecules containing only C, H, N and O atoms show the efficiency of Charger in terms of execution time and accuracy. Charger is also powerful in a study of electrostatic symbiosis between a protein and a ligand. It determines reliable protein-ligand interaction energies even when both contain S atoms. It easily estimates the individual contribution of every residue to the total protein-ligand electrostatic binding energy. Glutathione transferase (GST) in complex with a benzophenone ligand was studied due to the availability of both structural and thermodynamic data. The resulting analysis highlights not only the residues that stabilize the ligand but also those that hinder ligand binding from an electrostatic point of view. This offers new perspectives in the search for mutations to improve the interaction between the two partners. A proposed mutation would improve ligand binding to GST by removing an electrostatic obstacle, rather than by the traditional increase in the number of favourable contacts.
Collapse
Affiliation(s)
- Vedran Vuković
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
| | - Theo Leduc
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
| | - Zoe Jelić-Matošević
- Department of Chemistry, Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia
| | | | | | - Benoît Guillot
- Université de Lorraine, CNRS, CRM2, F-54000 Nancy, France
| | | |
Collapse
|
6
|
Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease. Molecules 2021; 26:molecules26133872. [PMID: 34202892 PMCID: PMC8270314 DOI: 10.3390/molecules26133872] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Revised: 06/19/2021] [Accepted: 06/22/2021] [Indexed: 11/23/2022] Open
Abstract
Computational analysis of protein–ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions, such as in docking programs, all rely on equations that sum each type of protein–ligand interactions in order to predict the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions constitute one of the most important part of total interactions between macromolecules. Unlike dispersion forces, they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this study, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and darunavir) were analyzed by using charge densities from the transferable aspherical-atom University at Buffalo Databank (UBDB). Moreover, we analyzed the electrostatic interaction energy for an ensemble of structures, using molecular dynamic simulations to highlight the main features of electrostatic interactions important for binding affinity.
Collapse
|
7
|
Kleemiss F, Wieduwilt EK, Hupf E, Shi MW, Stewart SG, Jayatilaka D, Turner MJ, Sugimoto K, Nishibori E, Schirmeister T, Schmidt TC, Engels B, Grabowsky S. Similarities and Differences between Crystal and Enzyme Environmental Effects on the Electron Density of Drug Molecules. Chemistry 2021; 27:3407-3419. [PMID: 33090581 PMCID: PMC7898524 DOI: 10.1002/chem.202003978] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Indexed: 01/28/2023]
Abstract
The crystal interaction density is generally assumed to be a suitable measure of the polarization of a low-molecular weight ligand inside an enzyme, but this approximation has seldomly been tested and has never been quantified before. In this study, we compare the crystal interaction density and the interaction electrostatic potential for a model compound of loxistatin acid (E64c) with those inside cathepsin B, in solution, and in vacuum. We apply QM/MM calculations and experimental quantum crystallography to show that the crystal interaction density is indeed very similar to the enzyme interaction density. Less than 0.1 e are shifted between these two environments in total. However, this difference has non-negligible consequences for derived properties.
Collapse
Affiliation(s)
- Florian Kleemiss
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Department of Chemistry and BiochemistryUniversity of BernFreiestrasse 3, 3012 BernSwitzerland
| | - Erna K. Wieduwilt
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019Université de Lorraine & CNRSBoulevard Arago, 57078 MetzFrance
| | - Emanuel Hupf
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
| | - Ming W. Shi
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Scott G. Stewart
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Dylan Jayatilaka
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Michael J. Turner
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Kunihisa Sugimoto
- Japan Synchrotron Radiation Research InstituteSPring-81-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5198Japan
- Institute for Integrated Cell-Materials SciencesKyoto UniversityYoshida-Ushinomiya-cho, Sakyo-ku, Kyoto 606-8501Japan
| | - Eiji Nishibori
- Division of Physics, Faculty of Pure and Applied Sciences, Tsukuba Research Center for Energy Materials ScienceUniversity of TsukubaTsukubaJapan
| | - Tanja Schirmeister
- Institute of Pharmaceutical and Biomedical SciencesJohannes-Gutenberg University MainzStaudingerweg 5, 55128 MainzGermany
| | - Thomas C. Schmidt
- Institute for Physical and Theoretical ChemistryJulius-Maximilians-University WürzburgEmil-Fischer-Str. 42, 97074 WürzburgGermany
| | - Bernd Engels
- Institute for Physical and Theoretical ChemistryJulius-Maximilians-University WürzburgEmil-Fischer-Str. 42, 97074 WürzburgGermany
| | - Simon Grabowsky
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Department of Chemistry and BiochemistryUniversity of BernFreiestrasse 3, 3012 BernSwitzerland
| |
Collapse
|
8
|
Jabłuszewska A, Krawczuk A, Dos Santos LHR, Macchi P. Accurate Modelling of Group Electrostatic Potential and Distributed Polarizability in Dipeptides. Chemphyschem 2020; 21:2155-2165. [DOI: 10.1002/cphc.202000441] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2020] [Revised: 07/23/2020] [Indexed: 12/28/2022]
Affiliation(s)
- Angelika Jabłuszewska
- Faculty of Chemistry Jagiellonian University in Krakow Gronostajowa 2 30-387 Krakow Poland
| | - Anna Krawczuk
- Faculty of Chemistry Jagiellonian University in Krakow Gronostajowa 2 30-387 Krakow Poland
| | - Leonardo H. R. Dos Santos
- Departamento de Química Universidade Federal de Minas Gerais Av. Antônio Carlos 6627 31270-901 Belo Horizonte MG Brazil
| | - Piero Macchi
- Department of Chemistry, Materials and Chemical Engineering Polytechnics of Milan Via Mancinelli 7 20131 Milan Italy
| |
Collapse
|