A new method is presented for docking molecular fragments to a rigid protein with evaluation of the binding energy. Polar fragments are docked with at least one hydrogen bond with the protein while apolar fragments are positioned in the hydrophobic pockets. The electrostatic contribution to the binding energy, which consists of screened intermolecular energy and protein and fragment desolvation terms, is evaluated efficiently by a numerical approach based on the continuum dielectric approximation. The latter is also used to predetermine the hydrophobic pockets of the protein by rolling a low dielectric sphere over the protein surface and calculating the electrostatic desolvation of the protein and van der Waals interaction energy. The method was implemented in the program SEED (solvation energy for exhaustive docking). The SEED continuum electrostatic approach has been successfully validated by a comparison with finite difference solutions of the Poisson equation for more than 2,500 complexes of small molecules with thrombin and the monomer of HIV-1 aspartic proteinase. The fragments docked by SEED in the active site of thrombin reproduce the structural features of the interaction patterns between known inhibitors and thrombin. Moreover, the combinatorial connection of these fragments yields a number of compounds that are very similar to potent inhibitors of thrombin. Proteins 1999;37:88-105.

No abstract

A new computational approach for the efficient docking of flexible ligands in a rigid protein is presented. It exploits the binding modes of functional groups determined by an exhaustive search with solvation. The search in ligand conformational space is performed by a genetic algorithm whose scoring function approximates steric effects and intermolecular hydrogen bonds. Ligand conformations generated by the genetic algorithm are docked in the protein binding site by optimizing the fit of their fragments to optimal positions of chemically related functional groups. We show that the use of optimal binding modes of molecular fragments allows to dock known inhibitors with about ten rotatable bonds in the active site of the uncomplexed and complexed conformations of thrombin and HIV-1 protease.

An improved version of the fragment-based flexible ligand docking approach SEED-FFLD is tested on inhibitors of human immunodeficiency virus type 1 protease, human alpha-thrombin and the estrogen receptor beta. The docking results indicate that it is possible to correctly reproduce the binding mode of inhibitors with more than ten rotatable bonds if the strain in their covalent geometry upon binding is not large. A high degree of convergence towards a unique binding mode in multiple runs of the genetic algorithm is proposed as a necessary condition for successful docking.

A new method is presented to quantitatively estimate and graphically display the propensity of nonpolar groups to bind at the surface of proteins. It is based on the calculation of the binding energy, i.e., van der Waals interaction plus protein electrostatic desolvation, of a nonpolar probe sphere rolled over the protein surface, and on the color coding of this quantity on a smooth molecular surface (hydrophobicity map). The method is validated on ten protein-ligand complexes and is shown to distinguish precisely where polar and nonpolar groups preferentially bind. Comparisons with existing approaches, like the display of the electrostatic potential or the curvature, illustrate the advantages and the better predictive power of the present method. Hydrophobicity maps will play an important role in the characterization of binding sites for the large number of proteins emerging from the genome projects and structure modeling approaches. Proteins 1999;37:565-575. 1999 Wiley-Liss, Inc.