Conformation simulated annealing

Other methods which are applied to conformational analysis and to generating multiple conformations and which can be regarded as random or stochastic techniques, since they explore the conformational space in a non-deterministic fashion, arc genetic algorithms (GA) [137, 1381 simulation methods, such as molecular dynamics (MD) and Monte Carlo (MC) simulations 1139], as well as simulated annealing [140], All of those approaches and their application to generate ensembles of conformations arc discussed in Chapter II, Section 7.2 in the Handbook. 

Quenched dynamics is a combination of high temperature molecular dynamics and energy minimization. This process determines the energy distribution of con formational families produced during molecular dynamics trajectories. To provide a better estimate of conformations, you should combine quenched dynamics with simulated annealing. 

Je next introduce the basic algorithms and then describe some of the mmy variants upon lem. We then discuss two methods called evolutionary algorithms and simulated anneal-ig, which are generic methods for locating the globally optimal solution. Finally, we discuss jme of the ways in which one might cinalyse the data from a conformational malysis in rder to identify a representative set of conformations. 

A particularly important application of molecular dynamics, often in conjunction with the simulated annealing method, is in the refinement of X-ray and NMR data to determine the three-dimensional structures of large biological molecules such as proteins. The aim of such refinement is to determine the conformation (or conformations) that best explain the experimental data. A modified form of molecular dynamics called restrained moleculai dynarrdcs is usually used in which additional terms, called penalty functions, are added tc the potential energy function. These extra terms have the effect of penalising conformations... 

There are similar algorithms, also called simulated annealing, that are Monte Carlo algorithms in which the choice conformations obey a Gaussian distribution centered on the lowest-energy value found thus far. The standard deviation of this distribution decreases over the course of the simulation. 

Model optimization is a further refinement of the secondary and tertiary structure. At a minimum, a molecular mechanics energy minimization is done. Often, molecular dynamics or simulated annealing are used. These are frequently chosen to search the region of conformational space relatively close to the starting structure. For marginal cases, this step is very important and larger simulations should be run. 

Reduce stress on molecules caused by a simulation at elevated temperatures. The cooling process, called simulated annealing, takes new, high energy conformational states toward stable conform ations. 

Quenched dynamics can trap structures in local minima. To prevent this problem, you can cool the system slowly to room temperature or some appropriate lower temperature. Then run room temperature molecular dynamics simulations to search for conformations that have lower energies, closer to the starting structure. Cooling a structure slowly is called simulated annealing. 

To overcome the limitations of the database search methods, conformational search methods were developed [95,96,109]. There are many such methods, exploiting different protein representations, objective function tenns, and optimization or enumeration algorithms. The search algorithms include the minimum perturbation method [97], molecular dynamics simulations [92,110,111], genetic algorithms [112], Monte Carlo and simulated annealing [113,114], multiple copy simultaneous search [115-117], self-consistent field optimization [118], and an enumeration based on the graph theory [119]. 

Restrained MD (rMD) is followed by the use of MD in explicit solvent, i.e. the conformation as determined above is taken into a box containing many solvent molecules around the molecule. Subsequently, simulated annealing (SA) and energy minimizahon steps are performed to draw the molecule into the global energy minimum. An MD run (the so-called trajectory) over at least 150ps to Ins is followed and a mean structure is calculated from such a trajectory. The con-formahon must be stable under this condihon even when the experimental constraints are removed. 

Refinement of Conformations by Computational Methods 243 9.4.2.6 Simulated Annealing (SA)... 

Roitberg, A. Elber, R., Modeling side chains in peptides and proteins application of the locally enhanced sampling technique and the simulated annealing methods to find minimum energy conformations, J. Chem. Phys. 1991, 95, 9277-9287... 

---