A MULTIPLE-START MONTE-CARLO DOCKING METHOD

被引:148
作者
HART, TN
READ, RJ
机构
[1] Department of Medical Microbiology and Infectious Diseases, University of Alberta, Edmonton, Alberta
来源
PROTEINS-STRUCTURE FUNCTION AND GENETICS | 1992年 / 13卷 / 03期
关键词
MOLECULAR DOCKING; MONTE-CARLO; SIMULATED ANNEALING; RATIONAL DRUG-DESIGN; DIHYDROFOLATE REDUCTASE; PROTEINASE INHIBITORS;
D O I
10.1002/prot.340130304
中图分类号
Q5 [生物化学]; Q7 [分子生物学];
学科分类号
071010 ; 081704 ;
摘要
We present a method to search for possible binding modes of molecular fragments at a specific site of a potential drug target of known structure. Our method is based on a Monte Carlo (MC) algorithm applied to the translational and rotational degrees of freedom of the probe fragment. Starting from a randomly generated initial configuration, favorable binding modes are generated using a two-step process. An MC run is first performed in which the energy in the Metropolis algorithm is substituted by a score function that measures the average distance of the probe to the target surface. This has the effect of making buried probes move toward the target surface and also allows enhanced sampling of deep pockets. In a second MC run, a pairwise atom potential function is used, and the temperature parameter is slowly lowered during the run (Simulated Annealing). We repeat this procedure starting from a large number of different randomly generated initial configurations in order to find all energetically favorable docking modes in a specified region around the target. We test this method using two inhibitor-receptor systems: Streptomyces griseus proteinase B in complex with the third domain of the ovomucoid inhibitor from turkey, and dihydrofolate reductase from E. coli in complex with methotrexate. The method could consistently reproduce the complex found in the crystal structure searching from random initial positions in cubes ranging from 25 angstrom to 50 angstrom about the binding site. In the case of SGPB, we were also successful in docking to the native structure. In addition, we were successful in docking small probes in a search that included the entire protein surface.
引用
收藏
页码:206 / 222
页数:17
相关论文
共 29 条
[11]   USING SHAPE COMPLEMENTARITY AS AN INITIAL SCREEN IN DESIGNING LIGANDS FOR A RECEPTOR-BINDING SITE OF KNOWN 3-DIMENSIONAL STRUCTURE [J].
DESJARLAIS, RL ;
SHERIDAN, RP ;
SEIBEL, GL ;
DIXON, JS ;
KUNTZ, ID ;
VENKATARAGHAVAN, R .
JOURNAL OF MEDICINAL CHEMISTRY, 1988, 31 (04) :722-729
[12]   DRUG DESIGN BY THE METHOD OF RECEPTOR FIT [J].
GOODFORD, PJ .
JOURNAL OF MEDICINAL CHEMISTRY, 1984, 27 (05) :557-564
[14]   AUTOMATED DOCKING OF SUBSTRATES TO PROTEINS BY SIMULATED ANNEALING [J].
GOODSELL, DS ;
OLSON, AJ .
PROTEINS-STRUCTURE FUNCTION AND GENETICS, 1990, 8 (03) :195-202
[15]   CONSISTENT FORCE-FIELD STUDIES OF INTER-MOLECULAR FORCES IN HYDROGEN-BONDED CRYSTALS .2. BENCHMARK FOR THE OBJECTIVE COMPARISON OF ALTERNATIVE FORCE-FIELDS [J].
HAGLER, AT ;
LIFSON, S ;
DAUBER, P .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1979, 101 (18) :5122-5130
[16]   ENERGY FUNCTIONS FOR PEPTIDES AND PROTEINS .1. DERIVATION OF A CONSISTENT FORCE-FIELD INCLUDING HYDROGEN-BOND FROM AMIDE CRYSTALS [J].
HAGLER, AT ;
HULER, E ;
LIFSON, S .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1974, 96 (17) :5319-5327
[17]   CONSISTENT FORCE-FIELD STUDIES OF INTER-MOLECULAR FORCES IN HYDROGEN-BONDED CRYSTALS .3. C=O...H-O HYDROGEN-BOND AND THE ANALYSIS OF THE ENERGETICS AND PACKING OF CARBOXYLIC-ACIDS [J].
HAGLER, AT ;
DAUBER, P ;
LIFSON, S .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1979, 101 (18) :5131-5141
[18]   PROTEIN CRYSTALLOGRAPHY AND COMPUTER-GRAPHICS TOWARD RATIONAL DRUG DESIGN [J].
HOL, WGJ .
ANGEWANDTE CHEMIE-INTERNATIONAL EDITION IN ENGLISH, 1986, 25 (09) :767-778
[19]   OPTIMIZATION BY SIMULATED ANNEALING [J].
KIRKPATRICK, S ;
GELATT, CD ;
VECCHI, MP .
SCIENCE, 1983, 220 (4598) :671-680
[20]   A GEOMETRIC APPROACH TO MACROMOLECULE-LIGAND INTERACTIONS [J].
KUNTZ, ID ;
BLANEY, JM ;
OATLEY, SJ ;
LANGRIDGE, R ;
FERRIN, TE .
JOURNAL OF MOLECULAR BIOLOGY, 1982, 161 (02) :269-288