A Gibbs free energy correlation for automated docking of carbohydrates

Thermodynamic information can be inferred from static atomic configurations. To model the thermodynamics of carbohydrate binding to proteins accurately, a large binding data set has been assembled from the literature. The data set contains information from 262 unique protein-carbohydrate crystal structures for which experimental binding information is known. Hydrogen atoms were added to the structures and training conformations were generated with the automated docking program AutoDock 3.06, resulting in a training set of 225,920 all-atom conformations. In all, 288 formulations of the AutoDock 3.0 free energy model were trained against the data set, testing each of four alternate methods of computing the van der Waals, solvation, and hydrogen-bonding energetic components. The van der Waals parameters from AutoDock I produced the lowest errors, and an entropic model derived from statistical mechanics produced the only models with five physically and statistically significant coefficients. Eight models predict the Gibbs free energy of binding with an error of less than 40% of the error of any similar models previously published.