A molecular dynamics study of thermodynamic and structural aspects of the hydration of cavities in proteins

The structure and activity of a protein molecule are strongly influenced by the extent of hydration of its cavities. This is, in turn, related to the free energy change on transfer of a water molecule from bulk solvent into a cavity. Such free energy changes have been calculated for two cavities in a sulfate‐binding protein. One of these cavities contains a crystallo graphically observed water molecule while the other does not. Thermodynamic integration and perturbation methods were used to calculate free energies of hydration for each of the cavities from molecular dynamics simulations of two separate events: the removal of a water molecule from pure water, and the introduction of a water molecule into each protein cavity. From the simulations for the pure water system, the excess chemical potential of water was computed to be −6.4 ± 0.4 kcal/mol, in accord with experiment and with other recent theoretical calculations. For the protein cavity containing an experimentally observed water molecule, the free energy change on hydrating it with one water molecule was calculated as −10.0 ± 1.3 kcal/mol, indicating the high probability that this cavity is occupied by a water molecule. By contrast, for the cavity in which no water molecules were experimentally observed, the free energy change on hydrating it with one water molecule was calculated as 0.2 ± 1.5 kcal/mol, indicating its low occupancy by water. The agreement of these results with experiment suggests that thermodynamic simulation methods may become useful for the prediction and analysis of internal hydration in proteins.

[1]  A. Joachimiak,et al.  Crystal structure of trp represser/operator complex at atomic resolution , 1988, Nature.

[2]  J. Mccammon,et al.  Dynamics of Proteins and Nucleic Acids , 2018 .

[3]  M Ptashne,et al.  Recognition of a DNA operator by the repressor of phage 434: a view at high resolution , 1988, Science.

[4]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[5]  Rebecca C. Wade.  Solvation of the active site of cytochrome P450-cam , 1990, J. Comput. Aided Mol. Des..

[6]  M. Mezei.  Excess free energy of different water models computed by Monte Carlo methods , 1982 .

[7]  B Pullman,et al.  New Paths in the Molecular Orbital Approach to Solvation of Biological Molecules , 1974, Quarterly Reviews of Biophysics.

[8]  J. Goodfellow.  Cooperative effects in water-biomolecule crystal systems. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[9]  P. A. Bash,et al.  Free energy calculations by computer simulation. , 1987, Science.

[10]  Arieh Ben-Naim,et al.  Solvation thermodynamics of nonionic solutes , 1984 .

[11]  P M Cullis,et al.  Affinities of amino acid side chains for solvent water. , 1981, Biochemistry.

[12]  B C Finzel,et al.  Three-dimensional structure of an antibody-antigen complex. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Terry P. Lybrand,et al.  Hydration of chloride and bromide anions: determination of relative free energy by computer simulation , 1985 .

[14]  L. Sieker,et al.  Water structure in a protein crystal: rubredoxin at 1.2 A resolution. , 1978, Journal of molecular biology.

[15]  Charles H. Bennett,et al.  Mass tensor molecular dynamics , 1975 .

[16]  R. Sharon,et al.  Accurate simulation of protein dynamics in solution. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[17]  T. Poulos,et al.  High-resolution crystal structure of cytochrome P450cam. , 1987, Journal of molecular biology.

[18]  M Karplus,et al.  Solvent effects on protein motion and protein effects on solvent motion. Dynamics of the active site region of lysozyme. , 1989, Journal of molecular biology.

[19]  M. Teeter,et al.  Water structure of a hydrophobic protein at atomic resolution: Pentagon rings of water molecules in crystals of crambin. , 1981, Proceedings of the National Academy of Sciences of the United States of America.

[20]  J. Hermans,et al.  The Free Energy of Xenon Binding to Myoglobin from Molecular Dynamics Simulation , 1986 .

[21]  A. Kossiakoff,et al.  The Neutron Structure of the Hydrophobic Plant Protein Crambin , 1982 .

[22]  Charles Tanford,et al.  The hydrophobic effect , 1980 .

[23]  W. Richards.  Ab Initio Methods and the Study of Molecular Hydration , 1979 .

[24]  Computer simulation of aqueous biomolecular systems. , 1985, Journal of biomolecular structure & dynamics.

[25]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[26]  J. Goodfellow,et al.  The Structure and Dynamics of Water in Globular Proteins , 1982 .

[27]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[28]  W. Kauzmann.  Some factors in the interpretation of protein denaturation. , 1959, Advances in protein chemistry.

[29]  W. Saenger.  Structure and dynamics of water surrounding biomolecules. , 1987, Annual review of biophysics and biophysical chemistry.

[30]  R. Wade,et al.  New hydrogen-bond potentials for use in determining energetically favorable binding sites on molecules of known structure. , 1989, Journal of medicinal chemistry.

[31]  R W Hockney,et al.  Computer Simulation Using Particles , 1966 .

[32]  George D. Rose,et al.  Prediction of chain turns in globular proteins on a hydrophobic basis , 1978, Nature.

[33]  William L. Jorgensen,et al.  Free energy of TIP4P water and the free energies of hydration of CH4 and Cl- from statistical perturbation theory , 1989 .

[34]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[35]  Hydration of cavities in proteins : a molecular dynamics approach , 1990 .

[36]  T. Straatsma,et al.  Free energy of hydrophobic hydration: A molecular dynamics study of noble gases in water , 1986 .

[37]  Felix Franks,et al.  Water:A Comprehensive Treatise , 1972 .

[38]  A. D. McLachlan,et al.  Solvation energy in protein folding and binding , 1986, Nature.

[39]  G. Otting,et al.  Studies of protein hydration in aqueous solution by direct NMR observation of individual protein-bound water molecules , 1989 .

[40]  F A Quiocho,et al.  The 2 A resolution structure of the sulfate-binding protein involved in active transport in Salmonella typhimurium. , 1988, Journal of molecular biology.

[41]  T. Straatsma,et al.  Free energy of ionic hydration: Analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations , 1988 .

[42]  Charles L. Brooks,et al.  Thermodynamics of aqueous solvation: Solution properties of alcohols and alkanes , 1987 .

[43]  P. Artymiuk,et al.  X-ray studies of water in crystals of lysozyme. , 1983, Journal of molecular biology.

[44]  H. Mckenzie,et al.  Water and proteins. II. The location and dynamics of water in protein systems and its relation to their stability and properties. , 1983, Advances in biophysics.

[45]  Ab initio calculation of the free energy of liquid water , 1978 .

[46]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[47]  Marçal De Oliveira Neto,et al.  Rapid Location of the Preferred Interaction Sites between Small Polar Molecules and Macromolecules. I. Binding of Water to the Component Units of Nucleic Acids , 1986 .

[48]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[49]  E. Baker,et al.  Hydrogen bonding in globular proteins. , 1984, Progress in biophysics and molecular biology.

[50]  P. Goodford.  A computational procedure for determining energetically favorable binding sites on biologically important macromolecules. , 1985, Journal of medicinal chemistry.

[51]  P. Nagy,et al.  Orientation and structure-building role of the water molecules bound at the contact surface of the dihydrofolate reductase-methotrexate complex , 1988, J. Comput. Aided Mol. Des..

[52]  J. Bolin,et al.  Crystal structures of Escherichia coli and Lactobacillus casei dihydrofolate reductase refined at 1.7 A resolution. I. General features and binding of methotrexate. , 1982, The Journal of biological chemistry.

[53]  B Honig,et al.  Internal cavities and buried waters in globular proteins. , 1986, Biochemistry.

[54]  W. L. Jorgensen,et al.  Monte Carlo simulation of differences in free energies of hydration , 1985 .

[55]  H. Berendsen,et al.  A consistent empirical potential for water–protein interactions , 1984 .

[56]  B. Widom,et al.  Some Topics in the Theory of Fluids , 1963 .

[57]  K. Nakanishi,et al.  Monte Carlo studies on the hydrophobic hydration in dilute aqueous solutions of nonpolar molecules , 1979 .

[58]  J. Hermans,et al.  Excess free energy of liquids from molecular dynamics simulations. Application to water models. , 1988, Journal of the American Chemical Society.

[59]  J. Andrew McCammon,et al.  Mass and step length optimization for the calculation of equilibrium properties by molecular dynamics simulation , 1990 .