Insights into affinity and specificity in the complexes of α-lytic protease and its inhibitor proteins: binding free energy from molecular dynamics simulation

被引:17
作者
Deng, Nan-Jie [1 ]
Cieplak, Piotr [2 ]
机构
[1] Univ Sci & Technol China, Sch Life Sci, Hefei 230027, Anhui, Peoples R China
[2] Burnham Inst Med Res, La Jolla, CA 92037 USA
关键词
ELECTROSTATIC INTERACTIONS; CONTINUUM SOLVENT; EGLIN-C; DECOMPOSITION; ASSOCIATION; MECHANICS; SOLVATION; STABILITY; RNA;
D O I
10.1039/b820961h
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
We report the binding free energy calculation and its decomposition for the complexes of alpha-lytic protease and its protein inhibitors using molecular dynamics simulation. Standard mechanism serine protease inhibitors eglin C and OMTKY3 are known to have strong binding affinity for many serine proteases. Their binding loops have significant similarities, including a common P1 Leu as the main anchor in the binding interface. However, recent experiments demonstrate that the two inhibitors have vastly different affinity towards alpha-lytic protease (ALP), a bacterial serine protease. OMTKY3 inhibits the enzyme much more weakly (by similar to 10(6) times) than eglin C. Moreover, a variant of OMTKY3 with five mutations, OMTKY3M, has been shown to inhibit 10(4) times more strongly than the wild-type inhibitor. The underlying mechanisms for the unusually large difference in binding affinities and the effect of mutation are not well understood. Here we use molecular dynamics simulation with molecular mechanics-Poisson Boltzmann/surface area method (MM-PB/SA) to investigate quantitatively the binding specificity. The calculated absolute binding free energies correctly differentiate the thermodynamic stabilities of these protein complexes, but the magnitudes of the binding affinities are systematically overestimated. Analysis of the binding free energy components provides insights into the molecular mechanism of binding specificity. The large Delta Delta G(bind) between eglin C and wild type OMTKY3 towards ALP is mainly attributable to the stronger nonpolar interactions in the ALP-eglin C complex, arising from a higher degree of structural complementarity. Here the electrostatic interaction contributes to a lesser extent. The enhanced inhibition in the penta-mutant OMTKY3M over its wild type is entirely due to an overall improvement in the solvent-mediated electrostatic interactions in the ALP-OMTKY3M complex. The results suggest that for these protein-complexes and similar enzyme-inhibitor systems (1) the binding is driven by nonpolar interactions, opposed by overall electrostatic and solute entropy contributions; (2) binding specificity can be tuned by improving the complementarity in electrostatics between two associating proteins. Binding free energy decomposition into contributions from individual protein residues provides additional detailed information on the structural determinants and subtle conformational changes responsible for the binding specificity.
引用
收藏
页码:4968 / 4981
页数:14
相关论文
共 37 条
[1]   Binding free energies and free energy components from molecular dynamics and Poisson-Boltzmann calculations. Application to amino acid recognition by aspartyl-tRNA synthetase [J].
Archontis, G ;
Simonson, T ;
Karplus, M .
JOURNAL OF MOLECULAR BIOLOGY, 2001, 306 (02) :307-327
[2]   Small-molecule inhibitors of protein-protein interactions: Progressing towards the dream [J].
Arkin, MR ;
Wells, JA .
NATURE REVIEWS DRUG DISCOVERY, 2004, 3 (04) :301-317
[3]  
Bamberg P., 1991, COURSE MATH STUDENTS
[4]   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
[5]   ZDOCK: An initial-stage protein-docking algorithm [J].
Chen, R ;
Li, L ;
Weng, ZP .
PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS, 2003, 52 (01) :80-87
[6]   Computer simulation of protein-protein interactions [J].
Elcock, AH ;
Sept, D ;
McCammon, JA .
JOURNAL OF PHYSICAL CHEMISTRY B, 2001, 105 (08) :1504-1518
[7]   A SMOOTH PARTICLE MESH EWALD METHOD [J].
ESSMANN, U ;
PERERA, L ;
BERKOWITZ, ML ;
DARDEN, T ;
LEE, H ;
PEDERSEN, LG .
JOURNAL OF CHEMICAL PHYSICS, 1995, 103 (19) :8577-8593
[8]   REFINED STRUCTURE OF ALPHA-LYTIC PROTEASE AT 1.7-A RESOLUTION - ANALYSIS OF HYDROGEN-BONDING AND SOLVENT STRUCTURE [J].
FUJINAGA, M ;
DELBAERE, LTJ ;
BRAYER, GD ;
JAMES, MNG .
JOURNAL OF MOLECULAR BIOLOGY, 1985, 184 (03) :479-502
[9]   Entropy-enthalpy compensation in solvation and ligand binding revisited [J].
Gallicchio, E ;
Kubo, MM ;
Levy, RM .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1998, 120 (18) :4526-4527
[10]   The statistical-thermodynamic basis for computation of binding affinities: A critical review [J].
Gilson, MK ;
Given, JA ;
Bush, BL ;
McCammon, JA .
BIOPHYSICAL JOURNAL, 1997, 72 (03) :1047-1069