Enhanced docking with the mining minima optimizer: acceleration and side-chain flexibility

The ligand-protein docking algorithm based on the Mining Minima method has been substantially enhanced. First, the basic algorithm is accelerated by: (1) adaptively determining the extent of each energy well to help avoid previously discovered energy minima; (2) biasing the search away from ligand positions at the surface of the receptor to prevent the ligand from staying at the surface when large sampling regions are used; (3) quickly testing multiple different ligand positions and orientations for each ligand conformation; and (4) tuning the source code to increase computational efficiency. These changes markedly shorten the time needed to discover an accurate result, especially when large sampling regions are used. The algorithm now also allows user-selected receptor sidechains to be treated as mobile during the docking procedure. The energies associated with the mobile side chains are computed as if they belonged to the ligand, except that atoms at the boundary between side chains and the rigid backbone are treated specially. This new capability is tested for several well-known ligand/protein systems, and preliminary application to an enzyme whose substrate is unknown--the recently solved hypothetical protein YecO (HI0319) from Haemophilus influenzae--indicates that side-chains relaxations allow candidate substrates of various sizes to be accommodated.