Packing optimization for automated generation of complex system's initial configurations for molecular dynamics and docking

被引:624
作者
Martínez, JM
Martínez, L
机构
[1] Univ Estadual Campinas, UNICAMP, IMECC, Dept Appl Math, BR-13081970 Campinas, SP, Brazil
[2] Univ Estadual Campinas, Inst Chem, BR-13081970 Campinas, SP, Brazil
关键词
large-scale optimization; box constraints; molecular dynamics; docking;
D O I
10.1002/jcc.10216
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
Molecular Dynamics is a powerful methodology for the comprehension at molecular level of many chemical and biochemical systems. The theories and techniques developed for structural and thermodynamic analyses are well established, and many software packages are available. However, designing starting configurations for dynamics can be cumbersome. Easily generated regular lattices can be used when simple liquids or mixtures are studied. However, for complex mixtures, polymer solutions or solid adsorbed liquids (for example) this approach is inefficient, and it turns out to be very hard to obtain an adequate coordinate file. In this article, the problem of obtaining an adequate initial configuration is treated as a "packing" problem and solved by an optimization procedure. The initial configuration is chosen in such a way that the minimum distance between atoms of different molecules is greater than a fixed tolerance. The optimization uses a well-known algorithm for box-constrained minimization. Applications are given for biomolecule solvation, many-component mixtures, and interfaces. This approach can reduce the work of designing starting configurations from days or weeks to few minutes or hours, in an automated fashion. Packing optimization is also shown to be a powerful methodology for space search in docking of small ligands to proteins. This is demonstrated by docking of the thyroid hormone to its nuclear receptor. (C) 2003 Wiley Periodicals, Inc.
引用
收藏
页码:819 / 825
页数:7
相关论文
共 32 条
[1]   A critical comparison of equilibrium, non-equilibrium and boundary-driven molecular dynamics techniques for studying transport in microporous materials [J].
Arya, G ;
Chang, HC ;
Maginn, EJ .
JOURNAL OF CHEMICAL PHYSICS, 2001, 115 (17) :8112-8124
[2]  
Aste T., 2000, The pursuit of perfect packing
[3]  
Bielschowsky R. H., 1998, INVESTIGACION OPERAT, V7, P67
[4]   MOLECULAR MECHANICS SIMULATION OF PROTEIN LIGAND INTERACTIONS - BINDING OF THYROID-HORMONE ANALOGS TO PRE-ALBUMIN [J].
BLANEY, JM ;
WEINER, PK ;
DEARING, A ;
KOLLMAN, PA ;
JORGENSEN, EC ;
OATLEY, SJ ;
BURRIDGE, JM ;
BLAKE, CCF .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1982, 104 (23) :6424-6434
[5]   CHARMM - A PROGRAM FOR MACROMOLECULAR ENERGY, MINIMIZATION, AND DYNAMICS CALCULATIONS [J].
BROOKS, BR ;
BRUCCOLERI, RE ;
OLAFSON, BD ;
STATES, DJ ;
SWAMINATHAN, S ;
KARPLUS, M .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1983, 4 (02) :187-217
[6]  
Conn A. R., 2000, TRUST REGION METHODS, DOI [10.1137/1.9780898719857, DOI 10.1137/1.9780898719857]
[7]  
DENNIS JE, 1983, NUMERICAL METHODS UN
[8]   Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution [J].
Duan, Y ;
Kollman, PA .
SCIENCE, 1998, 282 (5389) :740-744
[9]   An active set Newton algorithm for large-scale nonlinear programs with box constraints [J].
Facchinei, F ;
Judice, J ;
Soares, J .
SIAM JOURNAL ON OPTIMIZATION, 1998, 8 (01) :158-186
[10]  
Frenkel D., 1996, UNDERSTANDING MOL SI