Cresset has combined the analysis of ligand electrostatics with that of protein electrostatics to generate a visual and numerical assessment of the Electrostatic Complementarity™ (EC) of ligand protein complexes. The Cresset EC method is computationally inexpensive and can be applied to large data sets. We have analysed both the visual and numerical components by applying the method to a range of literature sets, showing that it correlates with and can predict reported bioactivity differences – see our paper in J. Med Chem. EC analysis has also been applied to kinase selectivity prediction, streptavidin mutant analysis and selection of ligand protonation states in protein context. The results of this work have shown that assessing EC can be a powerful tool for analysis and optimization of electrostatic protein-ligand interactions, making it possible to quantify electrostatically driven SAR and predict electrostatic target selectivity. The Cresset protein-ligand EC is calculated from comparison of protein and ligand electrostatic potential (ESP) values at all vertices of a generated ligand or protein solvent accessible surface (SAS) [1,2]. Perfect electrostatic complementarity would mean that at each vertex the ligand ESP value would be paired with a protein ESP value of the same magnitude with reverse sign. Below is a visualization of both ESP and EC of the biotin-streptavidin complex on both ligand and protein SAS, showing how a close matching of positive and negative electrostatic potential areas of protein and ligand leads to a good EC.

Below, we calculate an EC score that provides an approximate correction for some desolvation effects and allows local visualization of EC on a protein or ligand solvent-accessible surface where the integral is over the ligand SAS, ESPL and ESPP, the ligand and protein ESP values, and MAX(ESPL,ESPP,k) the protein or ligand ESP value with the larger deviation from zero, or a constant k if that is larger. Ligand and protein ESP values are capped to the maximum ESP values observed for water molecules to approximately correct for desolvation. EC scores range from 1 to -1, corresponding to a perfect EC or complete electrostatic clash, respectively. As solvent-exposed portions of the ligand contribute less information about the electrostatic complementarity of protein-ligand complexes, regions of the ligand SAS that are more than 3 Å away from any protein atom are scaled down by a distance-dependent factor.

The Cresset EC method has been applied to a reported XIAP-BIR3 data set, demonstrating that this method can detect and quantify electrostatic differences in XIAP ligands that cause changes in bioactivity.
References
- Xu, D.; Li, H.; Zhang, Y. Protein Depth Calculation and the Use for Improving Accuracy of Protein Fold Recognition. J Comput Biol 2013, 20 (10), 805–816. https://doi.org/10.1089/cmb.2013.0071.
- Xu, D.; Zhang, Y. Generating Triangulated Macromolecular Surfaces by Euclidean Distance Transform. PLoS ONE 2009, 4 (12), e8140. https://doi.org/10.1371/journal.pone.0008140.
