A rush to explore protein–ligand electrostatic interaction energy with Charger
- 22 September 2021
- journal article
- research article
- Published by International Union of Crystallography (IUCr) in Acta Crystallographica Section D Structural Biology
- Vol. 77 (10) , 1292-1304
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.Keywords
Funding Information
- MENSRT
This publication has 35 references indexed in Scilit:
- Charge Density and Electrostatic Interactions of Fidarestat, an Inhibitor of Human Aldose ReductaseJournal of the American Chemical Society, 2009
- Combining crystallographic information and an aspherical-atom data bank in the evaluation of the electrostatic interaction energy in an enzyme–substrate complex: influenza neuraminidase inhibitionActa Crystallographica Section D Structural Biology, 2009
- A Quantitative Model of Thermal Stabilization and Destabilization of Proteins by LigandsBiophysical Journal, 2008
- Kernel Energy Method: Application to DNABiochemistry, 2005
- A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculationsJournal of Computational Chemistry, 2003
- High-Density Miniaturized Thermal Shift Assays as a General Strategy for Drug DiscoverySLAS Discovery, 2001
- Roothaan-Hartree-Fock wave functions for atoms withZ≤54Physical Review A, 1993
- Benchmark values for two‐center Coulomb integrals over slater‐type orbitalsInternational Journal of Quantum Chemistry, 1993
- Testing aspherical atom refinements on small-molecule data setsActa Crystallographica Section A, 1978
- Roothaan-Hartree-Fock atomic wavefunctionsAtomic Data and Nuclear Data Tables, 1974