Abstract
An algorithm for finding an approximate global minimum of a funnel shaped function with many local minima is described. It is applied to compute the minimum energy docking position of a ligand with respect to a protein molecule. The method is based on the iterative use of a convex, general quadratic approximation that underestimates a set of local minima, where the error in the approximation is minimized in the L1 norm. The quadratic approximation is used to generate a reduced domain, which is assumed to contain the global minimum of the funnel shaped function. Additional local minima are computed in this reduced domain, and an improved approximation is computed. This process is iterated until a convergence tolerance is satisfied. The algorithm has been applied to find the global minimum of the energy function generated by the Docking Mesh Evaluator program. Results for three different protein docking examples are presented. Each of these energy functions has thousands of local minima. Convergence of the algorithm to an approximate global minimum is shown for all three examples.
Similar content being viewed by others
References
N. Baker, M. Holst, and F. Wang, “Adaptive multilevel finite element solution of the Poisson-Boltzmann equation II: Refinement at solvent accessible surfaces in biomolecular systems,” J. Comput. Chem., vol. 21, pp. 1343–1352, 2000.
K.A. Dill, A.T. Phillips, and J.B. Rosen, CGU: An algorithm for molecular structure prediction, in Large-scale optimization with applications, Part III (Minneapolis, MN, 1995), vol. 94 of IMA Vol. Math. Appl., Springer: New York, 1997, pp. 1–21.
P.E. Gill, W. Murray, M.A. Saunders, and M.H. Wright, “User's guide for NPSOL (version 4.0): A Fortran package for nonlinear programming,” Tech. Rep. SOL-86-2, Systems Optimization Laboratory, Stanford University, Stanford, CA, 1986.
D.S. Goodsell and A.J. Olson, “Automated docking of substrates to proteins by simulated annealing,” Proteins: Struct. Fun. Gen., vol. 8, pp. 195–202, 1990.
M. Holst, N. Baker, and F. Wang, “Adaptive multilevel finite element solution of the Poisson-Boltzmann equation I: Algorithms and examples,” J. Comput. Chem., vol. 21, pp. 1319–1342, 2000.
R. Hooke and T. A. Jeeves, “Direct search solution of numerical and statistical problems,” J. Assoc. Comput. Mach., vol. 8, pp. 212–229, 1961.
J. Janin, “Welcome to CAPRI: A critical assessment of predicted interactions,” Proteins: Struct. Func. Genet., vol. 47, no. 3, p. 257, 2002.
E.C. Meng, B.K. Shoichet, and I.D. Kuntz, “Automated docking with grid-based energy evaluation,” J. Comp. Chem., vol. 13, pp. 505–524, 1992.
J.C. Mitchell, J.B. Rosen, A.T. Phillips, and L.F. Ten Eyck, “Coupled optimization in protein docking,” in Proceedings of the Third Annual International Conference on Computational Molecular Biology, ACM Press, 1999, pp. 280–284.
J.B. Rosen, K.A. Dill, and A.T. Phillips, “Protein structure and energy landscape dependence on sequence using a continuous energy function,” J. Comp. Biol., vol. 4, pp. 227–239, 1997.
J.B. Rosen and R.F. Marcia, “Convex quadratic approximation,” Comput. Optim. Appl., vol. 28, pp. 173–184, 2004.
J.B. Rosen, H. Park, J. Glick, and L. Zhang, “Accurate solution to overdetermined linear equations with errors using L sb 1 norm minimization,” Comput. Optim. Appl., vol. 17, pp. 329–341, 2000.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Marcia, R.F., Mitchell, J.C. & Rosen, J.B. Iterative Convex Quadratic Approximation for Global Optimization in Protein Docking. Comput Optim Applic 32, 285–297 (2005). https://doi.org/10.1007/s10589-005-4799-4
Received:
Revised:
Accepted:
Issue Date:
DOI: https://doi.org/10.1007/s10589-005-4799-4