Validation of Poisson-Boltzmann Electrostatic Potential Fields in 3D QSAR: A CoMFA Study on Multiple Datasets

A CoMFA 3D QSAR (Quantitative Structure-Activity Relationships) study was performed on five diverse data sets that have been previously published or are otherwise generally available to assess the value of the ZAP electrostatic potential field (calculated from a solution of the Poisson-Boltzmann Equation) as compared to the Coulombic field as implemented in standard CoMFA. The ZAP field showed an average improvement of 0.114 in q 2 (leave-one-out cross-validation) as compared to the Coulombic field, and the steric/ZAP combination field showed an improvement of 0.128 in q 2 as compared to the Both (steric and electrostatic) field combination of CoMFA. However, the enhancement is not monolithic - some of the data sets in this study showed a much larger preference for the ZAP field than others. The ZAP field appears to be more physically realistic than the Coulombic field and should be of significant value in QSAR studies of some data sets.

[1]  Tudor I. Oprea,et al.  Ligand-based identification of environmental estrogens. , 1996, Chemical research in toxicology.

[2]  R. Cramer,et al.  Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. , 1988, Journal of the American Chemical Society.

[3]  A. Doweyko,et al.  The hypothetical active site lattice. An approach to modelling active sites from data on inhibitor molecules. , 1988, Journal of medicinal chemistry.

[4]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

[5]  Chris L. Waller,et al.  Polarizability Fields for Use in Three-Dimensional Quantitative Structure-Activity Relationship (3D-QSAR) , 2001, J. Chem. Inf. Comput. Sci..

[6]  Glen E Kellogg,et al.  3D QSAR in modern drug design. , 2003, EXS.

[7]  Tudor I. Oprea,et al.  Three-dimensional quantitative structure-activity relationship of human immunodeficiency virus (I) protease inhibitors. 2. Predictive power using limited exploration of alternate binding modes. , 1994, Journal of medicinal chemistry.

[8]  Simon K. Kearsley,et al.  An alternative method for the alignment of molecular structures: Maximizing electrostatic and steric overlap , 1990 .

[9]  Use of electron densities in comparative molecular field analysis (CoMFA): O?H bond dissociation energies in phenols , 1999 .

[10]  B. Honig,et al.  Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis , 1988, Proteins.

[11]  G. Cruciani,et al.  Generating Optimal Linear PLS Estimations (GOLPE): An Advanced Chemometric Tool for Handling 3D‐QSAR Problems , 1993 .

[12]  A. Newman,et al.  Three-dimensional quantitative structure-activity relationships of mazindol analogues at the dopamine transporter. , 2002, Journal of Medicinal Chemistry.

[13]  J. Andrew Grant,et al.  A smooth permittivity function for Poisson–Boltzmann solvation methods , 2001, J. Comput. Chem..

[14]  Tudor I. Oprea,et al.  Three-dimensional QSAR of human immunodeficiency virus (I) protease inhibitors. 1. A CoMFA study employing experimentally-determined alignment rules. , 1993, Journal of medicinal chemistry.

[15]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[16]  G. Klebe,et al.  Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity. , 1994, Journal of medicinal chemistry.

[17]  Lemont B. Kier,et al.  E-state fields: Applications to 3D QSAR , 1996, J. Comput. Aided Mol. Des..

[18]  Glen Eugene Kellogg,et al.  HINT: A new method of empirical hydrophobic field calculation for CoMFA , 1991, J. Comput. Aided Mol. Des..

[19]  Alexander Tropsha,et al.  The effect of region size on CoMFA analyses , 1999 .