664 research outputs found
Protein Structure Prediction Using Basin-Hopping
Associative memory Hamiltonian structure prediction potentials are not overly
rugged, thereby suggesting their landscapes are like those of actual proteins.
In the present contribution we show how basin-hopping global optimization can
identify low-lying minima for the corresponding mildly frustrated energy
landscapes. For small systems the basin-hopping algorithm succeeds in locating
both lower minima and conformations closer to the experimental structure than
does molecular dynamics with simulated annealing. For large systems the
efficiency of basin-hopping decreases for our initial implementation, where the
steps consist of random perturbations to the Cartesian coordinates. We
implemented umbrella sampling using basin-hopping to further confirm when the
global minima are reached. We have also improved the energy surface by
employing bioinformatic techniques for reducing the roughness or variance of
the energy surface. Finally, the basin-hopping calculations have guided
improvements in the excluded volume of the Hamiltonian, producing better
structures. These results suggest a novel and transferable optimization scheme
for future energy function development
Thermodynamics and the Global Optimization of Lennard-Jones clusters
Theoretical design of global optimization algorithms can profitably utilize
recent statistical mechanical treatments of potential energy surfaces (PES's).
Here we analyze the basin-hopping algorithm to explain its success in locating
the global minima of Lennard-Jones (LJ) clusters, even those such as \LJ{38}
for which the PES has a multiple-funnel topography, where trapping in local
minima with different morphologies is expected. We find that a key factor in
overcoming trapping is the transformation applied to the PES which broadens the
thermodynamic transitions. The global minimum then has a significant
probability of occupation at temperatures where the free energy barriers
between funnels are surmountable.Comment: 13 pages, 13 figures, revte
The Energy Landscape, Folding Pathways and the Kinetics of a Knotted Protein
The folding pathway and rate coefficients of the folding of a knotted protein
are calculated for a potential energy function with minimal energetic
frustration. A kinetic transition network is constructed using the discrete
path sampling approach, and the resulting potential energy surface is
visualized by constructing disconnectivity graphs. Owing to topological
constraints, the low-lying portion of the landscape consists of three distinct
regions, corresponding to the native knotted state and to configurations where
either the N- or C-terminus is not yet folded into the knot. The fastest
folding pathways from denatured states exhibit early formation of the
N-terminus portion of the knot and a rate-determining step where the C-terminus
is incorporated. The low-lying minima with the N-terminus knotted and the
C-terminus free therefore constitute an off-pathway intermediate for this
model. The insertion of both the N- and C-termini into the knot occur late in
the folding process, creating large energy barriers that are the rate limiting
steps in the folding process. When compared to other protein folding proteins
of a similar length, this system folds over six orders of magnitude more
slowly.Comment: 19 page
The Activation-Relaxation Technique : ART nouveau and kinetic ART
The evolution of many systems is dominated by rare activated events that occur on timescale ranging from nanoseconds to the hour or more. For such systems, simulations must leave aside the full thermal description to focus specifically on mechanisms that generate a configurational change. We present here the activation relaxation technique (ART), an open-ended saddle point search algorithm, and a series of recent improvements to ART nouveau and kinetic ART, an ART-based on-the-fly off-lattice self-learning kinetic Monte Carlo method
Unbiased Global Optimization of Lennard-Jones Clusters for N <= 201 by Conformational Space Annealing Method
We apply the conformational space annealing (CSA) method to the Lennard-Jones
clusters and find all known lowest energy configurations up to 201 atoms,
without using extra information of the problem such as the structures of the
known global energy minima. In addition, the robustness of the algorithm with
respect to the randomness of initial conditions of the problem is demonstrated
by ten successful independent runs up to 183 atoms. Our results indicate that
the CSA method is a general and yet efficient global optimization algorithm
applicable to many systems.Comment: revtex, 4 pages, 2 figures. Physical Review Letters, in pres
Recommended from our members
Sphere Encapsulated Monte Carlo: Obtaining Minimum Energy Configurations of Large Aromatic Systems.
We introduce a simple global optimization approach that is able to find minimum energy configurations of clusters containing aromatic molecules. The translational and rotational perturbations required in Monte Carlo-based methods often lead to unrealistic configurations within which two or more molecular rings intersect, causing many of the computational steps to be rejected and the optimization process to be inefficient. Here we develop a modification of the basin-hopping global optimization procedure tailored to tackle problems with intersecting molecular rings. Termed the Sphere Encapsulated Monte Carlo (SEMC) method, this method introduces sphere-based rearrangement and minimization steps at each iteration, and its performance is shown through the exploration of potential energy landscapes of polycyclic aromatic hydrocarbon (PAH) clusters, systems of interest in combustion and astrophysics research. The SEMC method provides clusters that are accurate to 5% mean difference of the minimum energy at a 10-fold speed up compared to previous work using advanced molecular dynamics simulations. Importantly, the SEMC method captures key structural characteristics and molecular size partitioning trends as measured by the molecular radial distances and coordination numbers. The advantages of the SEMC method are further highlighted in its application to previously unstudied heterogeneous PAH clusters
A global optimization approach for searching low energy conformations of proteins
De novo protein structure prediction and understanding the protein folding mechanism is an outstanding challenge of Biological Physics. Relying on the thermodynamic hypothesis of protein folding it is expected that the native state of a protein can be found out if the global minimum of the free energy surface is found. To understand the energy landscape or the free energy surface is challenging. The structure and dynamics of proteins are the manifestations of the underlying potential energy surface. Here the potential energy function stands on a framework of all-atom representation and uses purely physics-based interactions. For the solvated proteins the effective free energy is defined as an implicit solvation model which includes the solvation free energy, along with a standard all-atom biomolecular forcefield. A major challenge is to search for the global minimum on this effective free energy surface. In this work the Minima Hopping Algorithm (MHOP) to find global minima on potential energy surfaces has been used for protein structure prediction or in general finding the lowest energy conformations of proteins. Here proteins have been studied both in vacuo and in the aqueous medium. For short peptides starting from a completely extended conformation we could find conformational minima which are very close to the experimentally observed structures
- …