Abstract:A genetic algorithm (GA) combined with a tabu search (TA) has been applied as a minimization method to rake the appropriate associated sites for some biomolecular systems. In our docking procedure, surface complementarity and energetic complementarity of a ligand with its receptor have been considered separately in a two-stage docking method. The first stage was to find a set of potential associated sites mainly based on surface complementarity using a genetic algorithm combined with a tabu search. This step corresponds with the process of finding the potential binding sites where pharmacophores will bind. In the second stage, several hundreds of GA minimization steps were performed for each associated site derived from the first stage mainly based on the energetic complementarity. After calculations for both of the two stages, we can offer several solutions of associated sites for every complex. In this paper, seven biomolecular systems, including five bound complexes and two unbound complexes, were chosen from the Protein Data Bank (PDB) to test our method. The calculated results were very encouraging-the hybrid minimization algorithm successfully reaches the correct solutions near the best binded modes for these protein complexes. The docking results not only predict the bound complexes very well, but also get a relatively accurate complexed conformation for unbound systems. For the five bound complexes, the results show that surface complementarity is enough to find the precise binding modes, the top solution from the tabu list generally corresponds to the correct binding mode. For the two unbound complexes, due to the conformational changes upon binding, it seems more difficult to get their correct binding conformations. The predicted results show that the correct binding mode also corresponds to a relatively large surface complementarity score. In these two test cases, the correct solution can be found in the top several solutions from the tabu list. For unbound complexes, the interaction energy from energetic complementarity is very important, it can be used to filter these solutions from the surface complementarity. After the evaluation of the energetic complementarity, the conformations and orientations close to the crystallographically determined structures are resolved. In most cases, the smallest root mean square distance (r.m.s.d.) from the GA combined with TA solutions is in a relatively small region. Our program of automatic docking is really a universal one among the procedures used for the theoretical study of molecular recognition.

Conformational Analysis of Peptides Using Monte Carlo Simulations Combined with the Genetic Algorithm

Constraint Monte-Carlo Simulated Annealing: Application in Peptide Conformational Analysis in Solution

Automated Docking of Peptides and Proteins by Genetic Algorithm

First-principles molecular structure search with a genetic algorithm

Improving the efficiency of evolutionary de novo peptide design: strategies for probing configuration and parameter settings.

An Enhanced Comparative Molecular Field Analysis Method Using Genetic Algorithm

Protein Folding Simulations of the Hydrophobic–hydrophilic Model by Combining Tabu Search with Genetic Algorithms

Application of Simulated Annealing: Study on Relationship between Conformation of Linear Peptides and Their Cyclization Yields

GA-Novo: De Novo Peptide Sequencing via Tandem Mass Spectrometry using Genetic Algorithm

A Configurational Bias Monte Carlo Method for Linear and Cyclic Peptides

Computational peptide discovery with a genetic programming approach

A Simplified Systematic Search Method and Its Application in Peptide Conformational Analysis

Annealing Genetic Algorithm for Protein Folding Simulations in the 3D HP Model

Guided macro-mutation in a graded energy based genetic algorithm for protein structure prediction

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

Computational Methods for Studying Conformational Behaviors of Cyclic Peptides.

Gaussian process: an alternative approach for QSAM modeling of peptides

Conformational Analysis of Two Glycoproteins: A Monte Carlo Simulated Annealing Approach Using A Soft-Sphere Potential

ConfGen: A Conformational Search Method for Efficient Generation of Bioactive Conformers

Searching low-energy conformers of neutral and protonated di-, tri-, and tetra-glycine using first-principles accuracy assisted by the use of neural network potentials