Analysis and Application of Potential Energy Smoothing and Search Methods for Global Optimization

Abstract

Global energy optimization of a molecular system is difficult due to the well-known "multiple minimum" problem. The rugged potential energy surface (PES) characteristic of multidimensional systems can be transformed reversibly using potential smoothing to generate a new surface that is easier to search for favorable configurations. The diffusion equation method (DEM) is one example of a potential smoothing algorithm. Potential smoothing as implemented in DEM is intuitively appealing and has certain appropriate statistical mechanical properties, but often fails to identify the global minimum even for relatively small problems. In the present paper, extensions to DEM capable of correcting its empirical behavior are systematically investigated. Two types of local search (LS) procedures are applied during the reversing schedule from the smooth deformed PES to the undeformed surface. Changes needed to generate smoothable versions of standard molecular mechanics force fields such as AMBER/OPLS and MM2 are also described. The resulting methods are applied in an attempt to determine the global energy minimum for a variety of systems in different coordinate representations. The problems studied include argon clusters, cycloheptadecane, capped polyalanine, and the docking of R-helices. Depending on the specific problem, potential smoothing and search (PSS) is performed in Cartesian, torsional, or rigid body space. For example, PSS finds a very low energy structure for cycloheptadecane with much greater efficiency than a search restricted to the undeformed potential surface. It is shown that potential smoothing is characterized by three salient features. As the level of smoothing is increased, unique minima merge into a common basin, crossings can occur in the relative energies of a pair of minima, and the spatial locations of minima are shifted due to the averaging effects of smoothing. Local search procedures improve the ability of smoothing methods to locate global minima because they facilitate the post facto correction of errors due to energy crossings that may have occurred at higher levels of smoothing. PSS methods should serve as useful tools for global energy optimization on a variety of difficult problems of practical interest.

Structural modeling of mutant α-glucosidases resulting in a processing/transport defect in Pompe disease

“…The root-mean-square gradient value was set at 0.05 kcal/mol Å . The structural models of mutant species incorporating the amino acid substitutions were constructed using molecular modeling software, TINKER, developed by Ponder et al [13][14][15][16][17] Calculation of the numbers of atoms influenced by amino acid substitutions…”

Evolution of the potential energy surface with size for Lennard-Jones clusters

“…14 For the above reasons LJ 38 and LJ 75 provide relatively difficult test cases for any unbiased global optimization algorithm. It is only relatively recently that algorithms have begun to find the truncated octahedron, [15][16][17]38,[45][46][47][48][49][50] and only one unbiased method 17 and one method that involves seeding 38 have reported finding the Marks decahedron.…”

Structural study on mutant α-l-iduronidases: insight into mucopolysaccharidosis type I

“…Furthermore, we examined the distributions and degrees of threedimensional structural changes caused by representative mutations by means of color imaging. (Kundrot et al 1991;Dudek and Ponder 1995;Kong and Ponder 1997;Pappu et al 1998;Ren and Ponder 2003). As a template, the homology model of human IDUA (PDB: 1Y24) (Rempel et al 2005) was used, and energy minimization was performed.…”