BACKGROUND
experimental high-throughput screening processes that drive structure-guided drug design efforts are effective tools for the identification of candidate molecular ligands that may tightly bind a target protein; however, such an approach often proves to be a costly endeavor, in terms of both time and financial expense, one that can potentially be alleviated with reliable in silico protein-ligand binding affinity models to assist in winnowing the search space  <cit> .
a diverse array of computational approaches to model binding affinity have been described in the literature, each of which focuses on a unique combination of physicochemical properties and interactions: x-score  <cit> , lig-score  <cit> , drugscore  <cit> , sfcscore  <cit> , autodock <dig>  <cit> , itscore  <cit> , and phoenix  <cit>  are just a few examples of such predictive tools.
here we describe our development of a model for predicting protein-ligand binding energy that relies on delaunay tessellation, a computational geometry technique  <cit> , for the purpose of objectively capturing nearest neighbor atomic four-body interactions in the structures of macromolecular complexes .
first, we compute the propensities for occurrence of all atomic quadruplet interactions by applying the tessellation procedure to atomic coordinates for a diverse cross-section of over  <dig> high-resolution macromolecular crystal structures, and the data collected are used in generating an atomic four-body potential.
tasked with distinguishing native structures as having global energy minima relative to decoys, our knowledge-based potential performs well compared to several related atomic energy functions  <cit> ; however this work constitutes substantial research outside the immediate focus of this study, and accordingly it will be reported elsewhere.
next, we apply our atomic potential to a separate dataset of three hundred diverse protein-ligand complexes, each selected for having both a solved high-resolution crystal structure and a known dissociation constant , the latter quantity being useful for determining the gibbs free energy of binding .
two thirds of the complexes are randomly selected to train our predictive model of binding affinity: in each case, the entire complex is tessellated and then scored using the four-body potential, as is the structure of the isolated protein without its bound ligand, and we derive an empirical linear function of the difference between these scores to predict Δg values.
the remaining one hundred complexes are then used to validate the capability of the trained linear model to predict binding energies for new protein-ligand complexes.
the steps taken to develop our model formed the basis of a recently published companion study  <cit> , and here we begin by carefully outlining those details below, since they lay the foundation for the next stage of the work to be presented.
in particular, the model is subsequently applied to the prediction of binding affinities for an independent, diverse test set of three hundred enzyme-inhibitor complexes for which high-resolution crystal structures, as well as experimentally determined inhibition constants , are available.
also, model performance is comprehensively benchmarked against a number of related methods from the literature.
methods
datasets
high-resolution  crystallographic structures for  <dig> macromolecular complexes , culled using the pisces server  <cit>  and having protein chains that share low  sequence identity, were selected to develop the four-body statistical potential.
dataset diversity is also reflected in the fact that the complexes consist of both single chain and multimeric proteins, many of which have bound ligands in the form of either small molecules or peptides.
each complex has a coordinate file deposited in the protein data bank   <cit> , and following the removal of all hydrogen atoms and water molecules, delaunay tessellation is applied to each structure file by using all the remaining atomic coordinates.
in order to train and validate our model for predicting binding affinity, we selected another diverse set of three hundred protein-ligand complexes  from the binding moad  <cit>  database.
the binding moad is a repository for all protein-ligand complexes that have high-resolution  structures deposited in the pdb, and if available, published experimental binding energy data.
focusing specifically on a non-redundant subset of the binding moad, both to ensure diversity of complexes as well as to minimize bias due to over-represented structures, we identified three hundred complexes having both pdb coordinate files as well as experimental dissociation constants .
the pdb accession codes and kd values for the protein-ligand complexes are tabulated in additional file  <dig>  as is the identity of the subset  into which each is randomly placed.
software and performance measurements
we use the qhull software package  <cit>  to carry out the atomic delaunay tessellations, matlab  service pack 1) to produce graphical depictions of the tessellations, and the ucsf chimera software package  <cit>  to generate all other molecular visualizations in this study.
codes to perform all data formatting and analyses tasks are written in the perl programming language.
given the dissociation constant  for a protein-ligand complex, the standard gibbs free energy of binding  can be determined using
Δg=rt ln= <dig> ×ln,
where r =  <dig>  ×10- <dig> kcal k- <dig> mol- <dig> is the gas constant and t = 298° k is the absolute temperature.
we evaluate the agreement between known  and predicted  binding energies by reporting the pearson's correlation coefficient
r=∑∑2 ∑ <dig>
the standard error of the predictions
se= ∑2n-2=1nn-2 n ∑yi2-∑yi2-n ∑xiyi-∑xi∑yi2n ∑xi2-∑xi <dig>
and the equation of the fitted regression line.
RESULTS
four-body statistical potential
to generate our knowledge-based potential, a six-letter alphabet  is used for labelling all atoms .
the qhull software uses the 3-dimensional  coordinates of atoms in a pdb file to generate a delaunay tessellation of the structure, a space-filling convex hull formed by hundreds of solid, non-overlapping, irregular tetrahedra whose vertices are the 3d atomic points.
each atom serves as a vertex, with most being shared by numerous adjacent tetrahedra, and every tetrahedral simplex objectively identifies a quadruplet of nearest neighbor atoms at its four vertices.
to ensure this is indeed the case, we eliminate all edges longer than 8Å immediately upon tessellation, which is in agreement with related research in this arena at the atomic  <cit>  and residue  <cit>  levels of coarse-graining.
the combined total number of tetrahedra remaining for analysis after tessellating the  <dig> pdb coordinate files is provided in table  <dig>  as are the total number of atoms of each type as well as their relative frequencies.
without regards to the ordering of a quadruplet of atoms , and allowing for the repeated occurrence of atom types in any given quadruplet , there are  <dig> possible types of atomic quadruplets that can be enumerated based on the use of a 6-letter atomic alphabet .
for each quadruplet , we define fijkl as the observed proportion of all tetrahedral simplices obtained by tessellating all  <dig> structures to have those four types of atoms at the vertices; similarly, we let pijkl represent the rate expected by chance, which is based on relative frequencies of the six atom types in the structures  and calculated using a multinomial background distribution given by
fijkl
pijkl
sijkl
fijkl
pijkl
sijkl
pijkl=4!∏n=16tn!
∏n=16antn,where  ∑n=16an=1 and  ∑n=16tn= <dig>
in eq.
, an is the relative frequency of atom type n, while tn counts how many times atom type n appears in the quadruplet .
as a consequence of the inverted boltzmann principle  <cit> , the score sijkl = log  is proportional to the energy of quadruplet atomic interaction, and the set of  <dig> scored atomic quadruplets defines our four-body statistical potential .
topological scores
in order to develop our predictive model, the four-body potential is applied to the dataset of three hundred protein-ligand complexes compiled from the binding moad in the following manner.
for each complex, the atomic coordinates  in the pdb file are tessellated , each tetrahedron in the tessellation is scored using table  <dig> according to the four atoms at its vertices, and a normalized topological score  is calculated to be the sum of all the tetrahedral scores divided by the number of tetrahedra in the tessellation, a quantity that can be summarized compactly by the equation
q=1n ∑sijkl
next, atomic coordinates of the ligand are removed from the pdb file of the complex, and the procedure is repeated to compute q for the isolated protein .
lastly, we define the topological score difference
Δq=qcomplex-qprotein
for the complex.
in the next section, we compare computed Δq quantities with known Δg values for these complexes in order to develop a model for predicting binding energy.
an important underlying assumption in this formulation is that ligand size is small enough so that tetrahedra formed at the interface with the protein dominate purely internal atomic interactions within the ligand.
the calculated q values, for structures of the three hundred protein-ligand complexes, as well as the isolated proteins without their bound ligands, are tabulated in additional file  <dig>
predictive model of binding energy
a comparison of the calculated Δq quantities for our training set of two hundred randomly selected complexes with their experimental Δg values  reveals a correlation coefficient of r =  <dig> .
however, the Δq values are not uniform in sign, and they are on a significantly smaller scale relative to the standard gibbs free energy of binding  data; hence, they cannot be used directly as a representation of predicted Δg values .
both issues related to Δq values for the training data are addressed with an empirically derived linear function
Δgcalc=1/ <dig> ×Δq- <dig> ,
resulting in negative Δgcalc values in each case that also scale similarly to Δgexp.
owing to Δgcalc arising from a simple linear transformation of the Δq values, Δgcalc and Δgexp also display a correlation of r =  <dig>   with a fitted regression line of y =  <dig> x .
turning next to the validation set of one hundred complexes, we obtain Δgcalc values from their computed Δq quantities by utilizing the linear model given in eq.
that we empirically derived from the training data.
the predicted Δgcalc and known Δgexp values for these complexes again display a correlation of r =  <dig>   with a fitted regression line of y =  <dig> x -  <dig> , and a scatter plot of the validation data is superimposed over that of the training data in figure  <dig>  tabulated in additional file  <dig> are Δgexp and Δgcalc values for all three hundred protein-ligand complexes.
discussion
enzyme-inhibitor binding affinity prediction
in order to test the utility of our model through a practical application, we predict binding affinities for a diverse dataset of three hundred enzyme-inhibitor complexes , independent of those protein-ligand complexes used for training and validation, which are annotated with their respective experimental inhibition constants  in the non-redundant binding moad.
analogous to eq.
, we obtain the standard gibbs free energy of binding for each complex with the equation
Δgexp=rt ln= <dig> ×ln,
where r =  <dig>  ×10- <dig> kcal k- <dig> mol- <dig> is the gas constant and t = 298° k is the absolute temperature.
tabulated in additional file  <dig> are the pdb accession codes of the high-resolution  crystallographic structures, as well as the ki and Δgexp values, corresponding to these enzyme-inhibitor complexes.
next, we use the atomic coordinates  provided by the pdb structure file for each complex to generate a delaunay tessellation , from which we obtain a normalized topological score  by employing eq.
in conjunction with our atomic four-body statistical potential .
in a similar fashion, we generate a normalized topological score for the isolated protein without the bound inhibitor , by tessellating a modified version of the pdb file that excludes the atomic coordinates for the inhibitor.
lastly, we calculate the difference  between these normalized topological scores according to eq.
, which is subsequently used by our model in eq.
to yield a prediction for the enzyme-inhibitor binding affinity .
all normalized topological score and calculated binding affinity data are also tabulated in additional file  <dig>
for this dataset of three hundred enzyme-inhibitor complexes, the calculated Δq values and the experimental binding affinity Δgexp data display a correlation of r =  <dig> ; likewise, as discussed previously, the correlation between Δgcalc and Δgexp is similarly given by r =  <dig> , in this case with a calculated standard error for the predictions of se =  <dig>  kcal/mol .
comparisons to related methods
in the same way that our predictive model of protein-ligand binding affinity is evaluated on a test set of three hundred enzyme-inhibitor complexes as described in the previous section, other related methods similarly use test sets of complexes to validate their models.
hence, to directly compare our performance to that of other methods, binding affinity predictions are generated using our approach for complexes that form their test sets.
starting with x-score, wang et al.
<cit>  report predictions with their model on a test set of ten complexes that reflect a correlation of r =  <dig>  between experimental and predicted binding affinity , with a fitted regression line of y =  <dig> x +  <dig> .
on the identical dataset, predictions obtained with our model yield a correlation of r =  <dig>  and fitted regression line of y =  <dig> x -  <dig> , results that signify a clear improvement over those of x-score .
pkd
experimental Δg data for the complexes are derived from the experimental pkd values.
turning next to itscore, we discover that huang et al.
<cit>  utilize a benchmarking test set consisting of one hundred protein-ligand complexes, originally constructed by wang et al.
<cit> , to compare their scoring function and  <dig> other methods by ranking the respective pearson's correlation coefficients  between experimental and predicted binding affinities.
the test set is diverse, consisting of  <dig> different proteins as well as binding affinities that span nearly nine orders of magnitude.
by generating binding affinity predictions for these one hundred complexes with our model and calculating their correlation with the experimental data, we can subsequently determine our ranking among these  <dig> related approaches: itscore  <cit> , x-score  <cit> , dfire  <cit> , drugscorecsd  <cit> , drugscorepdb  <cit> , cerius2/plp  <cit> , sybyl/g-score  <cit> , sybyl/d-score  <cit> , sybyl/chemscore  <cit> , cerius2/pmf  <cit> , dock/ff  <cit> , cerius2/ludi  <cit> , cerius2/lig-score  <cit> , sybyl/f-score  <cit> , and autodock  <cit> .
the results of our predictions are summarized in figure  <dig>  which provides a scatter plot of calculated versus experimental binding energies for this dataset of one hundred complexes.
the plot reflects a correlation of r =  <dig>  , surpassing all of the other methods  and validating the reliability of our approach.
CONCLUSIONS
delaunay tessellation of atomic coordinates in a diverse dataset of macromolecular structures objectively identifies four-body atomic interactions, providing the raw data for developing a knowledge-based atomic four-body statistical contact potential.
this potential is used to score a separate diverse set of three hundred protein-ligand complexes with known binding affinities, as well as to score the isolated proteins without their bound ligands, based on their respective structure tessellations.
initially, the difference  between scores calculated for an entire complex and for its isolated protein is considered as a predictor of binding affinity; however, since these Δq do not scale as binding free energy values, two hundred randomly selected protein-ligand complexes from this set are used to empirically derive a linear function of Δq as a model for calculating the binding energy.
for this training set, we observe a correlation of r =  <dig>  between calculated and experimental binding energies, with a standard error of se =  <dig>  kcal/mol and a regression line of y =  <dig> x. validation of this model with the remaining one hundred complexes that were held out yields performance measures of r =  <dig>  and se =  <dig>  kcal/mol.
in an application of the method, our model is then used to predict binding energies for an independent and diverse test set of three hundred enzyme-inhibitor complexes, producing results that are consistent with those based on the training and validation data.
finally, we utilize a diverse test set of one hundred protein-ligand complexes to benchmark the binding energy predictions made with our model, and our correlation between calculated and experimental binding energies for this dataset surpasses those of all  <dig> related methods to which it is compared.
a key advantage with our approach is the ability to generate rapid predictions, typically under one second per complex.
competing interests
the author declares that they have no competing interests.
authors' contributions
mm conceived of the study, implemented the methods, analyzed the data, and wrote the manuscript.
supplementary material
additional file 1
pdb accession codes for the  <dig> macromolecular structures used to derive the atomic four-body statistical potential.
http://proteins.gmu.edu/automute/additional file  <dig> txt
click here for file
additional file 2
three hundred protein-ligand complexes used to train and validate the model.
http://proteins.gmu.edu/automute/additional file  <dig> txt
click here for file
additional file 3
independent set of three hundred enzyme-inhibitor complexes used to test the model.
http://proteins.gmu.edu/automute/additional file  <dig> txt
click here for file
acknowledgements
thanks to researchers affiliated with the binding moad repository for creating and making publicly available their centralized database of macromolecular complexes.
declarations
publication of this article was funded in part by the george mason university libraries open access publishing fund.
this article has been published as part of bmc structural biology volume  <dig> supplement  <dig>  2013: selected articles from the computational structural bioinformatics workshop  <dig>  the full contents of the supplement are available online at http://www.biomedcentral.com/bmcstructbiol/supplements/13/s <dig>

 experimental high-throughput screening processes that drive structure-guided drug design efforts are effective tools for the identification of candidate molecular ligands that may tightly bind a target protein; however, such an approach often proves to be a costly endeavor, in terms of both time and financial expense, one that can potentially be alleviated with reliable in silico protein-ligand binding affinity models to assist in winnowing the search space. a diverse array of computational approaches to model binding affinity have been described in the literature, each of which focuses on a unique combination of physicochemical properties and interactions: x-score, lig-score, drugscore, sfcscore, autodock, itscore, and phoenix are just a few examples of such predictive tools. here we describe our development of a model for predicting protein-ligand binding energy that relies on delaunay tessellation, a computational geometry technique, for the purpose of objectively capturing nearest neighbor atomic four-body interactions in the structures of macromolecular complexes. first, we compute the propensities for occurrence of all atomic quadruplet interactions by applying the tessellation procedure to atomic coordinates for a diverse cross-section of over high-resolution macromolecular crystal structures, and the data collected are used in generating an atomic four-body potential. tasked with distinguishing native structures as having global energy minima relative to decoys, our knowledge-based potential performs well compared to several related atomic energy functions; however this work constitutes substantial research outside the immediate focus of this study, and accordingly it will be reported elsewhere. next, we apply our atomic potential to a separate dataset of three hundred diverse protein-ligand complexes, each selected for having both a solved high-resolution crystal structure and a known dissociation constant, the latter quantity being useful for determining the gibbs free energy of binding. two thirds of the complexes are randomly selected to train our predictive model of binding affinity: in each case, the entire complex is tessellated and then scored using the four-body potential, as is the structure of the isolated protein without its bound ligand, and we derive an empirical linear function of the difference between these scores to predict Δg values. the remaining one hundred complexes are then used to validate the capability of the trained linear model to predict binding energies for new protein-ligand complexes. the steps taken to develop our model formed the basis of a recently published companion study, and here we begin by carefully outlining those details below, since they lay the foundation for the next stage of the work to be presented. in particular, the model is subsequently applied to the prediction of binding affinities for an independent, diverse test set of three hundred enzyme-inhibitor complexes for which high-resolution crystal structures, as well as experimentally determined inhibition constants, are available. also, model performance is comprehensively benchmarked against a number of related methods from the literature. methods datasets high-resolution crystallographic structures for macromolecular complexes, culled using the pisces server and having protein chains that share low sequence identity, were selected to develop the four-body statistical potential. dataset diversity is also reflected in the fact that the complexes consist of both single chain and multimeric proteins, many of which have bound ligands in the form of either small molecules or peptides. each complex has a coordinate file deposited in the protein data bank, and following the removal of all hydrogen atoms and water molecules, delaunay tessellation is applied to each structure file by using all the remaining atomic coordinates. in order to train and validate our model for predicting binding affinity, we selected another diverse set of three hundred protein-ligand complexes from the binding moad database. the binding moad is a repository for all protein-ligand complexes that have high-resolution structures deposited in the pdb, and if available, published experimental binding energy data. focusing specifically on a non-redundant subset of the binding moad, both to ensure diversity of complexes as well as to minimize bias due to over-represented structures, we identified three hundred complexes having both pdb coordinate files as well as experimental dissociation constants. the pdb accession codes and kd values for the protein-ligand complexes are tabulated in additional file as is the identity of the subset into which each is randomly placed. software and performance measurements we use the qhull software package to carry out the atomic delaunay tessellations, matlab service pack 1) to produce graphical depictions of the tessellations, and the ucsf chimera software package to generate all other molecular visualizations in this study. codes to perform all data formatting and analyses tasks are written in the perl programming language. given the dissociation constant for a protein-ligand complex, the standard gibbs free energy of binding can be determined using Δg=rt ln= ×ln, where r = ×10- kcal k- mol- is the gas constant and t = 298° k is the absolute temperature. we evaluate the agreement between known and predicted binding energies by reporting the pearson's correlation coefficient r=∑∑2 ∑ the standard error of the predictions se= ∑2n-2=1nn-2 n ∑yi2-∑yi2-n ∑xiyi-∑xi∑yi2n ∑xi2-∑xi and the equation of the fitted regression line.  four-body statistical potential to generate our knowledge-based potential, a six-letter alphabet is used for labelling all atoms. the qhull software uses the 3-dimensional coordinates of atoms in a pdb file to generate a delaunay tessellation of the structure, a space-filling convex hull formed by hundreds of solid, non-overlapping, irregular tetrahedra whose vertices are the 3d atomic points. each atom serves as a vertex, with most being shared by numerous adjacent tetrahedra, and every tetrahedral simplex objectively identifies a quadruplet of nearest neighbor atoms at its four vertices. to ensure this is indeed the case, we eliminate all edges longer than 8Å immediately upon tessellation, which is in agreement with related research in this arena at the atomic and residue levels of coarse-graining. the combined total number of tetrahedra remaining for analysis after tessellating the pdb coordinate files is provided in table as are the total number of atoms of each type as well as their relative frequencies. without regards to the ordering of a quadruplet of atoms, and allowing for the repeated occurrence of atom types in any given quadruplet, there are possible types of atomic quadruplets that can be enumerated based on the use of a 6-letter atomic alphabet. for each quadruplet, we define fijkl as the observed proportion of all tetrahedral simplices obtained by tessellating all structures to have those four types of atoms at the vertices; similarly, we let pijkl represent the rate expected by chance, which is based on relative frequencies of the six atom types in the structures and calculated using a multinomial background distribution given by fijkl pijkl sijkl fijkl pijkl sijkl pijkl=4!∏n=16tn! ∏n=16antn,where  ∑n=16an=1 and  ∑n=16tn= in eq., an is the relative frequency of atom type n, while tn counts how many times atom type n appears in the quadruplet. as a consequence of the inverted boltzmann principle, the score sijkl = log is proportional to the energy of quadruplet atomic interaction, and the set of scored atomic quadruplets defines our four-body statistical potential. topological scores in order to develop our predictive model, the four-body potential is applied to the dataset of three hundred protein-ligand complexes compiled from the binding moad in the following manner. for each complex, the atomic coordinates in the pdb file are tessellated, each tetrahedron in the tessellation is scored using table according to the four atoms at its vertices, and a normalized topological score is calculated to be the sum of all the tetrahedral scores divided by the number of tetrahedra in the tessellation, a quantity that can be summarized compactly by the equation q=1n ∑sijkl next, atomic coordinates of the ligand are removed from the pdb file of the complex, and the procedure is repeated to compute q for the isolated protein. lastly, we define the topological score difference Δq=qcomplex-qprotein for the complex. in the next section, we compare computed Δq quantities with known Δg values for these complexes in order to develop a model for predicting binding energy. an important underlying assumption in this formulation is that ligand size is small enough so that tetrahedra formed at the interface with the protein dominate purely internal atomic interactions within the ligand. the calculated q values, for structures of the three hundred protein-ligand complexes, as well as the isolated proteins without their bound ligands, are tabulated in additional file predictive model of binding energy a comparison of the calculated Δq quantities for our training set of two hundred randomly selected complexes with their experimental Δg values reveals a correlation coefficient of r =. however, the Δq values are not uniform in sign, and they are on a significantly smaller scale relative to the standard gibbs free energy of binding data; hence, they cannot be used directly as a representation of predicted Δg values. both issues related to Δq values for the training data are addressed with an empirically derived linear function Δgcalc=1/ ×Δq-, resulting in negative Δgcalc values in each case that also scale similarly to Δgexp. owing to Δgcalc arising from a simple linear transformation of the Δq values, Δgcalc and Δgexp also display a correlation of r = with a fitted regression line of y = x. turning next to the validation set of one hundred complexes, we obtain Δgcalc values from their computed Δq quantities by utilizing the linear model given in eq. that we empirically derived from the training data. the predicted Δgcalc and known Δgexp values for these complexes again display a correlation of r = with a fitted regression line of y = x -, and a scatter plot of the validation data is superimposed over that of the training data in figure tabulated in additional file are Δgexp and Δgcalc values for all three hundred protein-ligand complexes. discussion enzyme-inhibitor binding affinity prediction in order to test the utility of our model through a practical application, we predict binding affinities for a diverse dataset of three hundred enzyme-inhibitor complexes, independent of those protein-ligand complexes used for training and validation, which are annotated with their respective experimental inhibition constants in the non-redundant binding moad. analogous to eq., we obtain the standard gibbs free energy of binding for each complex with the equation Δgexp=rt ln= ×ln, where r = ×10- kcal k- mol- is the gas constant and t = 298° k is the absolute temperature. tabulated in additional file are the pdb accession codes of the high-resolution crystallographic structures, as well as the ki and Δgexp values, corresponding to these enzyme-inhibitor complexes. next, we use the atomic coordinates provided by the pdb structure file for each complex to generate a delaunay tessellation, from which we obtain a normalized topological score by employing eq. in conjunction with our atomic four-body statistical potential. in a similar fashion, we generate a normalized topological score for the isolated protein without the bound inhibitor, by tessellating a modified version of the pdb file that excludes the atomic coordinates for the inhibitor. lastly, we calculate the difference between these normalized topological scores according to eq., which is subsequently used by our model in eq. to yield a prediction for the enzyme-inhibitor binding affinity. all normalized topological score and calculated binding affinity data are also tabulated in additional file for this dataset of three hundred enzyme-inhibitor complexes, the calculated Δq values and the experimental binding affinity Δgexp data display a correlation of r =; likewise, as discussed previously, the correlation between Δgcalc and Δgexp is similarly given by r =, in this case with a calculated standard error for the predictions of se = kcal/mol. comparisons to related methods in the same way that our predictive model of protein-ligand binding affinity is evaluated on a test set of three hundred enzyme-inhibitor complexes as described in the previous section, other related methods similarly use test sets of complexes to validate their models. hence, to directly compare our performance to that of other methods, binding affinity predictions are generated using our approach for complexes that form their test sets. starting with x-score, wang et al. report predictions with their model on a test set of ten complexes that reflect a correlation of r = between experimental and predicted binding affinity, with a fitted regression line of y = x +. on the identical dataset, predictions obtained with our model yield a correlation of r = and fitted regression line of y = x -, results that signify a clear improvement over those of x-score. pkd experimental Δg data for the complexes are derived from the experimental pkd values. turning next to itscore, we discover that huang et al. utilize a benchmarking test set consisting of one hundred protein-ligand complexes, originally constructed by wang et al., to compare their scoring function and other methods by ranking the respective pearson's correlation coefficients between experimental and predicted binding affinities. the test set is diverse, consisting of different proteins as well as binding affinities that span nearly nine orders of magnitude. by generating binding affinity predictions for these one hundred complexes with our model and calculating their correlation with the experimental data, we can subsequently determine our ranking among these related approaches: itscore, x-score, dfire, drugscorecsd, drugscorepdb, cerius2/plp, sybyl/g-score, sybyl/d-score, sybyl/chemscore, cerius2/pmf, dock/ff, cerius2/ludi, cerius2/lig-score, sybyl/f-score, and autodock. the results of our predictions are summarized in figure which provides a scatter plot of calculated versus experimental binding energies for this dataset of one hundred complexes. the plot reflects a correlation of r =, surpassing all of the other methods and validating the reliability of our approach.  delaunay tessellation of atomic coordinates in a diverse dataset of macromolecular structures objectively identifies four-body atomic interactions, providing the raw data for developing a knowledge-based atomic four-body statistical contact potential. this potential is used to score a separate diverse set of three hundred protein-ligand complexes with known binding affinities, as well as to score the isolated proteins without their bound ligands, based on their respective structure tessellations. initially, the difference between scores calculated for an entire complex and for its isolated protein is considered as a predictor of binding affinity; however, since these Δq do not scale as binding free energy values, two hundred randomly selected protein-ligand complexes from this set are used to empirically derive a linear function of Δq as a model for calculating the binding energy. for this training set, we observe a correlation of r = between calculated and experimental binding energies, with a standard error of se = kcal/mol and a regression line of y = x. validation of this model with the remaining one hundred complexes that were held out yields performance measures of r = and se = kcal/mol. in an application of the method, our model is then used to predict binding energies for an independent and diverse test set of three hundred enzyme-inhibitor complexes, producing results that are consistent with those based on the training and validation data. finally, we utilize a diverse test set of one hundred protein-ligand complexes to benchmark the binding energy predictions made with our model, and our correlation between calculated and experimental binding energies for this dataset surpasses those of all related methods to which it is compared. a key advantage with our approach is the ability to generate rapid predictions, typically under one second per complex. competing interests the author declares that they have no competing interests. authors' contributions mm conceived of the study, implemented the methods, analyzed the data, and wrote the manuscript. supplementary material additional file 1 pdb accession codes for the macromolecular structures used to derive the atomic four-body statistical potential. http://proteins.gmu.edu/automute/additional file txt click here for file additional file 2 three hundred protein-ligand complexes used to train and validate the model. http://proteins.gmu.edu/automute/additional file txt click here for file additional file 3 independent set of three hundred enzyme-inhibitor complexes used to test the model. http://proteins.gmu.edu/automute/additional file txt click here for file acknowledgements thanks to researchers affiliated with the binding moad repository for creating and making publicly available their centralized database of macromolecular complexes. declarations publication of this article was funded in part by the george mason university libraries open access publishing fund. this article has been published as part of bmc structural biology volume supplement 2013: selected articles from the computational structural bioinformatics workshop the full contents of the supplement are available online at http://www.biomedcentral.com/bmcstructbiol/supplements/13/s 
two-thirds of these complexes are randomly selected to train a predictive model of binding affinity as follows: two tessellations are generated in each case, one for the entire complex and another strictly for the isolated protein without its bound ligand, and a topological score is computed for each tessellation with the four-body potential. next, we identify a second, independent set of three hundred protein-ligand complexes, having both high-resolution structures and known dissociation constants. here, we begin by applying the computational geometry technique of delaunay tessellation to each set of atomic coordinates for over diverse macromolecular structures, for the purpose of deriving a four-body statistical potential that serves as a topological scoring function. applying our model to an independent test set of high-resolution structures for three hundred diverse enzyme-inhibitor complexes, each with an experimentally known inhibition constant, also yields a correlation of r = between experimental and calculated binding energies. lastly, we generate predictions with our model on a diverse test set of one hundred protein-ligand complexes previously used to benchmark related methods, and our correlation of r = between the calculated and experimental binding energies for this dataset exceeds those of the other approaches. models that are capable of reliably predicting binding affinities for protein-ligand complexes play an important role the field of structure-guided drug design. to validate the method, we similarly generated two tessellations for each of the remaining protein-ligand complexes, computed their topological scores and the difference between the two scores for each complex, and applied the previously derived linear transformation of this topological score difference to predict binding affinities. predicted protein-ligand binding affinity is then based on an empirically derived linear function of the difference between both topological scores, one that appropriately scales the value of this difference.
 models that are capable of reliably predicting binding affinities for protein-ligand complexes play an important role the field of structure-guided drug design. methods here, we begin by applying the computational geometry technique of delaunay tessellation to each set of atomic coordinates for over diverse macromolecular structures, for the purpose of deriving a four-body statistical potential that serves as a topological scoring function. next, we identify a second, independent set of three hundred protein-ligand complexes, having both high-resolution structures and known dissociation constants. two-thirds of these complexes are randomly selected to train a predictive model of binding affinity as follows: two tessellations are generated in each case, one for the entire complex and another strictly for the isolated protein without its bound ligand, and a topological score is computed for each tessellation with the four-body potential. predicted protein-ligand binding affinity is then based on an empirically derived linear function of the difference between both topological scores, one that appropriately scales the value of this difference.  a comparison between experimental and calculated binding affinity values over the two hundred complexes reveals a pearson's correlation coefficient of r = with a standard error of se = kcal/mol. to validate the method, we similarly generated two tessellations for each of the remaining protein-ligand complexes, computed their topological scores and the difference between the two scores for each complex, and applied the previously derived linear transformation of this topological score difference to predict binding affinities. for these one hundred complexes, we again observe a correlation of r = between known and calculated binding affinities. applying our model to an independent test set of high-resolution structures for three hundred diverse enzyme-inhibitor complexes, each with an experimentally known inhibition constant, also yields a correlation of r = between experimental and calculated binding energies.  lastly, we generate predictions with our model on a diverse test set of one hundred protein-ligand complexes previously used to benchmark related methods, and our correlation of r = between the calculated and experimental binding energies for this dataset exceeds those of the other approaches. compared with these related prediction methods, our approach stands out based on salient features that include the reliability of our model, combined with the rapidity of the generated predictions, which are less than one second for an average sized complex. computational structural bioinformatics workshop philadelphia, pa, usa 