Inclusion of loss of translational and rotational freedom in theoretical estimates of free energies of binding. Application to a complex of benzene and mutant T4 lysozyme

被引:299
作者
Hermans, J
Wang, L
机构
[1] Dept. of Biochemistry and Biophysics, School of Medicine, University of North Carolina, Chapel Hill
[2] Dept. of Pharmaceutical Chemistry, University of California, San Francisco
关键词
D O I
10.1021/ja963568+
中图分类号
O6 [化学];
学科分类号
0703 ;
摘要
We present the first complete treatment for calculating theoretical estimates of free energies of formation of macromolecule-ligand complexes with molecular dynamics simulations, as the free energy for transforming the ligand into a non-interacting state by gradually diminishing the forces between macromolecule (plus solvent) and ligand. The calculations become possible due to the introduction of a specially designed potential (''molecular tweezers'') which restrains the spatial position and orientation of the ligand molecule and is gradually applied as the transformation proceeds from complexed to non-interacting components. The binding of benzene to a mutant T4 lysozyme (Morton et al. Biochemistry 1995, 34, 8564-8575) has been used as a test case. The simulations reproduce the value of the free energy of binding (-5.19 kcal/mol if the standard state of benzene is a 1 M aqueous solution) within the sum of experimental and statistical error. Another series of such simulations with rigid protein models provides an estimate of the dependence of the free energy of binding on the protein conformation. The free energy of binding is found to decrease in the series: energy-minimized ligand-free protein (-3.5 kcal/mol), energy-minimized ligand-containing protein (-6.3 kcal/mol), and crystal structure (-8.5 kcal/mol). The free energy of binding to a series of snapshots from a protein-ligand dynamics trajectory varies between -7 and -9 kcal/mol. The ''cratic'' free energy contribution, which corresponds to the loss of translational and rotational freedom of the ligand molecule, was estimated at 7 kcal/mol. It has proved possible to decompose this into translational and rotational components and, from these free energies,estimate the remaining freedom of the benzene in the binding pocket, at 0.6 Angstrom for positional range and 10-15 degrees for angular range, in excellent agreement with the motion observed in a dynamics trajectory.
引用
收藏
页码:2707 / 2714
页数:8
相关论文
共 40 条
  • [1] FREE-ENERGY CALCULATIONS BY COMPUTER-SIMULATION
    BASH, PA
    SINGH, UC
    LANGRIDGE, R
    KOLLMAN, PA
    [J]. SCIENCE, 1987, 236 (4801) : 564 - 568
  • [2] MOLECULAR-DYNAMICS WITH COUPLING TO AN EXTERNAL BATH
    BERENDSEN, HJC
    POSTMA, JPM
    VANGUNSTEREN, WF
    DINOLA, A
    HAAK, JR
    [J]. JOURNAL OF CHEMICAL PHYSICS, 1984, 81 (08) : 3684 - 3690
  • [3] FREE-ENERGY SIMULATIONS - THE MEANING OF THE INDIVIDUAL CONTRIBUTIONS FROM A COMPONENT ANALYSIS
    BORESCH, S
    ARCHONTIS, G
    KARPLUS, M
    [J]. PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS, 1994, 20 (01) : 25 - 33
  • [4] Denti TZM, 1996, J AM CHEM SOC, V118, P6044
  • [5] A CAVITY-CONTAINING MUTANT OF T4 LYSOZYME IS STABILIZED BY BURIED BENZENE
    ERIKSSON, AE
    BAASE, WA
    WOZNIAK, JA
    MATTHEWS, BW
    [J]. NATURE, 1992, 355 (6358) : 371 - 373
  • [6] BROWNIAN DYNAMICS WITH HYDRODYNAMIC INTERACTIONS
    ERMAK, DL
    MCCAMMON, JA
    [J]. JOURNAL OF CHEMICAL PHYSICS, 1978, 69 (04) : 1352 - 1360
  • [7] DETERMINATION OF THE RELATIVE BINDING FREE-ENERGIES OF PEPTIDE INHIBITORS TO THE HIV-1 PROTEASE
    FERGUSON, DM
    RADMER, RJ
    KOLLMAN, PA
    [J]. JOURNAL OF MEDICINAL CHEMISTRY, 1991, 34 (08) : 2654 - 2659
  • [8] FERRO DR, 1980, J MOL BIOL, V136, P1
  • [9] THE PRICE OF LOST FREEDOM - ENTROPY OF BIMOLECULAR COMPLEX-FORMATION
    FINKELSTEIN, AV
    JANIN, J
    [J]. PROTEIN ENGINEERING, 1989, 3 (01): : 1 - 3
  • [10] GILSON MK, 1997, IN PRESS BIOPHYS J