A genetic algorithm for flexible molecular overlay and pharmacophore elucidation

Summary A genetic algorithm (GA) has been developed for the superimposition of sets of flexible molecules. Molecules are represented by a chromosome that encodes angles of rotation about flexible bonds and mappings between hydrogen-bond donor proton, acceptor lone pair and ring centre features in pairs of molecules. The molecule with the smallest number of features in the data set is used as a template, onto which the remaining molecules are fitted with the objective of maximising structural equivalences. The fitness function of the GA is a weighted combination of: (i) the number and the similarity of the features that have been overlaid in this way; (ii) the volume integral of the overlay; and (iii) the van der Waals energy of the molecular conformations defined by the torsion angles encoded in the chromosomes. The algorithm has been applied to a number of pharmacophore elucidation problems, i.e., angiotensin II receptor antagonists, Leu-enkephalin and a hybrid morphine molecule, 5-HT_1D agonists, benzodiazepine receptor ligands, 5-HT₃ antagonists, dopamine D₂ antagonists, dopamine reuptake blockers and FKBP12 ligands. The resulting pharmacophores are generated rapidly and are in good agreement with those derived from alternative means.     

简化的系统搜索法及其在构象分析中的应用

多肽的构象研究除了可以用X-射线晶体学及二维核磁共振等实验方法外,理论计算的方法有系统搜索法、蒙特卡洛方法、距离几何方法、分子动力学方法及能量极小化等.系统搜索法具有构象空间搜索彻底的特点,相对来说找到系统整体极小值的可能性较大,但由于其计算量较大,对于多肽及蛋白质很难实现.我们对系统搜索法进行了简化,不是同时旋转所有的二面角,而是成对地进行,以迭代的方法达到收敛,最后得到可能的构象.     

Do active site conformations of small ligands correspond to low free-energy solution structures?

We compare the low free energy structures of ten small, polar ligands in solution to their conformations in their respective receptor active sites. The solution conformations are generated by a systematic search and the free energies of representative structures are computed with a continuum solvation model. Based on the values of torsion angles, we find little similarity between low energy solution structures of small ligands and their active site conformations. However, in nine out of ten cases, the positions of 'anchor points' (key atoms responsible for tight binding) in the lowest energy solution structures are very similar to the positions of these atoms in the active site conformations. A metric that more closely captures the essentials of binding supports the basic premise underlying pharmacophore mapping, namely that active site conformations of small flexible ligands correspond to their low energy structures in solution. This work supports the efforts of building pharmacophore models based on the information present in solution structures of small isolated ligands.     

HUNTER: A conformational search program for acyclic to polycyclic molecules with special emphasis on stereochemistry

A new conformational search program, HUNTER, connected with the force fields MMP2 and MM3(92) is presented. The program accepts all types of molecules with most different substructures, considers stereochemical facts, and covers conformational space efficiently and completely. The most important facilities are an automated analysis of the stereochemistry including topographical facts, a separate perturbation of the acyclic and cyclic parts of the molecule using modified corner flapping, and an incremental rotation around single bonds with fixed flap and rotation angles, respectively; an exclusion of high energy structures by simulated annealing; the choice of the conformer lowest in energy, which is new as an initial structure for the next sampling run; and the use of a reduced set of dihedral angles to define a conformation. A specifically devised graphic interface, SERVANT, is used to feed in and control all informations necessary for a program run and to visualize the results. Most of the parameters are user-defined and thereby allow a flexible search, including a search for the most stable diastereomer. The efficiency of the different parameter sets was tested in calculation with cycloundecane ( 12 ), (Z)-oct-3-ene ( 13 ), and sipholenol-A monoacetate ( 14 ). The best performance regarding the number of different low-energy conformers was achieved with 60° ( 14 ) and 90° flaps ( 12 ), respectively, including substituent correction for the cyclic parts, and with 105° ( 14 ) and 120° rotations ( 13 ), respectively, for the acyclic parts. In comparison to the stochastic search routine implemented in MM3(92), HUNTER performed two ( 12 ) to six ( 14 ) times better. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1264–1281, 1997     

Effect of electronic interactions between double bonds on the conformational flexibility of 1,4-cyclohexadiene

The electronic structure of 1,4-cyclohexadiene for various angles between the double bond planes has been calculated by the AM1 method. The effects of through-bond and through-space interactions, which result in flattening and unflattening, respectively, are oppositely directed. These effects are specified for various bending angles. In addition to steric factors, electronic factors affect the conformational flexibility of 1,4-cyclohexadiene.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1676–1677, September, 1994.     

An analytic ring closure condition for geometrical algorithm to search the conformational space

The recently reported geometrical algorithm to search the conformational space (GASCOS) scans conformational space exhaustively using an internal coordinate tree search. Using only geometrical operations and a set of criteria for eliminating chemically unreasonable atomic arrangements, the algorithm generates starting geometries for optimizations by molecular mechanics or by molecular orbital procedures. Up until now GASCOS has been used for linear structures, but an extension to cyclic structures is reported here.     

Raman Optical Activity of N-Acetyl-L-Cysteine in Water and in Methanol: The “Clusters-in-a-Liquid” Model and ab Initio Molecular Dynamics Simulations

Raman and Raman Optical Activity (ROA) spectra of N-acetyl-L-cysteine (NALC), a flexible chiral molecule, were measured in water and in methanol to evaluate the solvent effects. Two different solvation approaches, that is, the DFT based “clusters-in-a-liquid” solvent model and the ab initio molecular dynamics (AIMD) simulations, were applied to simulate the Raman and ROA spectra. Systematic conformational searches were carried out using a recently developed conformational searching tool, CREST, with the inclusion of polarizable continuum model of water and of methanol. The CREST candidates of NALC and the NALC-solvent complexes were re-optimized and their Raman and ROA simulations were done at the B3LYP−D3BJ/def2-TZVP and the B3LYP-aug-cc-pVDZ//cc-pVTZ levels. Also, AIMD simulations, which includes some anharmonic effects and all intermolecular interactions in solution, were performed. By empirically weighting the computed Raman and ROA spectra of each conformer, good agreements with the experimental data were achieved with both approaches, while AIMD offered some improvements in the carbonyl and in the low wavenumber regions over the static DFT approach. The pros and cons of these two different approaches for accounting the solvent effects on Raman and ROA of this flexible chiral system will also be discussed.     

Flexible ligand docking using a genetic algorithm

Summary Two computational techniques have been developed to explore the orientational and conformational space of a flexible ligand within an enzyme. Both methods use the Genetic Algorithm (GA) to generate conformationally flexible ligands in conjunction with algorithms from the DOCK suite of programs to characterize the receptor site. The methods are applied to three enzyme-ligand complexes: dihydrofolate reductase-methotrexate, thymidylate synthase-phenolpthalein and HIV protease-thioketal haloperidol. Conformations and orientations close to the crystallographically determined structures are obtained, as well as alternative structures with low energy. The potential for the GA method to screen a database of compounds is also examined. A collection of ligands is evaluated simultaneously, rather than docking the ligands individually into the enzyme.Abbreviations GA genetic algorithm; dhfr, dihydrofolate reductase - mtx methotrexate - ts thymidylate synthase - fen phenolphalein - HIV human immune deficiency virus - hivp HIV protease - thk thioketal haloperidol