Local vs. global search strategies in evolutionary GRID-based conformational sampling u00026; docking

摘要

Conformational sampling, the computational prediction of the experimental geometries of small proteins (folding) or of protein-ligand complexes (docking), is often cited as one of the most challenging multimodal optimization problems. Due to the extreme ruggedness of the energy landscape as a function of geometry, sampling heuristics must rely on an appropriate trade-off between global and local searching efforts. A previously reported ldquoplanetary strategyrdquo, a generalization of the classical island model used to deploy a hybrid genetic algorithm on computer grids, has shown a good ability to quickly discover low-energy geometries of small proteins and sugars, and sometimes even pinpoint their native structures-although not reproducibly. The procedure focused on broad exploration and used a tabu strategy to avoid revisiting the neighborhood of known solutions, at the risk of ldquoburyingrdquo important minima in overhastily set tabu areas. The strategy reported here, termed ldquodivide-and-conquer planetary modelrdquo couples this global search procedure to a local search tool. Grid nodes are now shared between global and local exploration tasks. The phase space is cut into ldquocellsrdquo corresponding to a specified sampling width for each of the N degrees of freedom. Global search locates cells containing low-energy geometries. Local searches pinpoint even deeper minima within a cell. Sampling width controls the important trade-off between the number of cells and the local search effort needed to reproducibly sample each cell. The probability to submit a cell to local search depends on the energy of the most stable geometry found within. Local searches are allotted limited resources and are not expected to converge. However, as long as they manage to discover some deeper local minima, the explored cell remains eligible for further local search, now relying on the improved energy level to enhance chances to be picked again. This competition prevents the syst-nem to waste too much effort in fruitless local searches. Eventually, after a limited number of local searches, a cell will be ldquoclosedrdquo and used - first as ldquoseedrdquo, later as tabu zone-to bias future global searches. Technical details and some folding and docking results will be discussed.