Previous Article |
Table of Contents
| Next Article
BIOPHYSICS
Cooperative water filling of a nonpolar protein cavity observed by high-pressure crystallography and simulation
, 

, 

*Department of Physics, Cornell University, Ithaca, NY 14853;
Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Building 5, Bethesda, MD 20892-0520; and
Institute of Molecular Biology, Howard Hughes Medical Institute and Department of Physics, University of Oregon, Eugene, OR 97403
Contributed by Brian W. Matthews, September 20, 2005
| Abstract |
|---|
|
|
|---|
Formation of a water-expelling nonpolar core is the paradigm of protein folding and stability. Although experiment largely confirms this picture, water buried in "hydrophobic" cavities is required for the function of some proteins. Hydration of the protein core has also been suggested as the mechanism of pressure-induced unfolding. We therefore are led to ask whether even the most nonpolar protein core is truly hydrophobic (i.e., water-repelling). To answer this question we probed the hydration of an
160-Å3, highly hydrophobic cavity created by mutation in T4 lysozyme by using high-pressure crystallography and molecular dynamics simulation. We show that application of modest pressure causes approximately four water molecules to enter the cavity while the protein itself remains essentially unchanged. The highly cooperative filling is primarily due to a small change in bulk water activity, which implies that changing solvent conditions or, equivalently, cavity polarity can dramatically affect interior hydration of proteins and thereby influence both protein activity and folding.
hydrophobic effect | T4 lysozyme
Additional progress will require direct measurement of the free energy needed to insert water molecules into a nonpolar protein cavity. Such an experiment requires sufficient structural resolution to unambiguously locate water molecules, the ability to resolve small changes in structure and occupancy, and the means to bias the occupancy while minimally perturbing the protein. Protein crystallography can detect changes in cavity hydration directly, but the temperature or chemical perturbations useful in solution thermodynamic measurements frequently damage crystals. We have found that high pressures of up to several hundred megapascals generally do not damage crystals (15, 16). High-pressure crystallography makes it possible to shift the equilibrium to interesting, unexplored protein states.
The
160-Å3 cavity containing mutant L99A of T4 lysozyme studied here was originally produced to probe the stabilizing interactions between buried nonpolar residues (2, 17). The Leu-99
Ala mutation creates a large cavity that destabilizes the folded protein because of lost interactions between side chains and the reduced free-energy cost of exposing alanine (instead of leucine) to water (17). L99A T4 lysozyme has also been the subject of ligand- and noble-gasbinding studies (1820) and is believed to be entirely empty under ambient conditions. Under pressure, one may expect a large cavity, similar to that produced by the L99A mutation, to collapse. Indeed, this was the original impetus to study this mutant under pressure. Instead, we find this cavity to be remarkably rigid and to fill with water at
150 MPa. It provides an excellent system in which to examine the interactions of water with surrounding protein.
| Methods |
|---|
|
|
|---|
2.1 Å by beryllium powder diffraction. Table 1 lists refinement statistics and unit-cell parameters at each pressure.
|
calc for each pair of data sets. Fo,hp and Fo,0.1 MPa are the experimentally measured structure factors at high pressure and 0.1 MPa.
calc are phases calculated from the ambient pressure structure. Observed structure factors are scaled to calculated structure factors to put them on an absolute scale. Because the unit cell decreases in size by as much as 1 Å as pressure increases, the highest-resolution data yield only noise in the maps. We limited the resolution used in the maps to 4 Å for this reason. B factors of atoms surrounding the cavity are among the lowest in the model, excluding the possibility of partial cavity-wall collapse. Maps were masked by using cavity descriptions that were generated with the program VOIDOO (22) using a probe radius of 1.2 Å and standard optimized potentials for liquid simulations (OPLS) atomic radii. The masked maps were integrated by using the program MAPMAN (23) and finally averaged over all data sets.
Molecular Dynamics Simulations. Molecular dynamics (MD) simulations of T4 lysozyme mutant L99A were performed with the AMBER 6.0 code (University of California, San Francisco) and parm94 force field (24). The initial structure (PDB ID code 1L90
[PDB]
) was solvated by 4,964 three-point transferable intermolecular potential (TIP3P) water molecules (25) and 18 sodium and 27 chloride ions in an
52 x 54 x 63-Å3 rectangular box. After energy minimization, equilibrium MD simulations were performed at constant temperature (300 K) and pressure (0.1 and 200 MPa) with periodic boundary conditions, particle-mesh Ewald summation, and a 1-fs time step. After 54 ps of MD simulations at 0.1 MPa with an empty cavity, water molecules were transferred from the surrounding bulk into the cavity to create states with occupancy n = 1, 2, 3, 4, and 5. Those configurations were equilibrated for 50 to 300 ps, followed by production runs of 1- to 1.7-ns duration. Structures were saved for analysis every 0.5 ps. For of each of the 18 runs (three pressures, n = 0, 1,..., 5), the MD structures averaged over the final 0.25 ns deviated by <0.7 Å from the x-ray structure (rms deviation of all
carbons) and <0.55 Å for the N- and C-terminal domains, respectively, alone (residues 1060 and 80160).
Thermodynamics of Water Filling. The thermodynamics of cavity filling is determined effectively by a "grand-canonical partition function" for water molecules inside the cavity (17). The occupancy probabilities P(N) of having exactly N water molecules in the cavity are related through
![]() | [1] |
P(0) is determined by normalization,
. V is a spherical volume with a radius of6Åcovering the cavity whose center is the instantaneous positional average of selected C
atoms of residues lining the cavity.
and µ' are the number density and excess chemical potential of bulk water at ambient conditions, respectively. The canonical average
...
N is over the Boltzmann factor of the potential-energy change
U = UN+1 UN of randomly inserting a water molecule into V, already occupied by N water molecules.
U includes interactions with protein and solvent. We evaluate the canonical average by using Bennett's method of overlapping histograms for the water-insertion and -removal energies, pins(u) and prem(u), respectively (17). The ratio pins(u)/prem(u) is equal to exp(
U) times the average in Eq. 1. For the chemical potential of the reference bulk fluid at ambient conditions, we use the previously calculated µ'=25.3 kJ/mol and
= 33.33 nm3 (10). At higher pressures, we calculated the equation of state
(p, T = 298 K) of three-point transferable intermolecular potential (TIP3P) water from constant-pressure simulations. Integration then gives the change of the total bulk chemical potential,
. At p1 = 100 and 200 MPa, we obtain
µ = 1.82 and 3.64 kJ/mol, respectively, where exp(
µ) = [
2exp(
µ'2)]/[
1exp(
µ'1)]. The average occupancy is
and can be calculated directly at the simulation pressures. At intermediate pressures we use perturbation theory at
200 MPa.
| Results |
|---|
|
|
|---|
/
N-terminal domain (residues 1058) relative to the C-terminal domain (residues 80162) by
0.25 Å and deformation of one end of helix C (at approximately residues 7680) by
0.25 Å toward the main cavity. Despite these changes and other smaller deformations (
0.15 Å) in the C-terminal domain, the cavity volume decreases by <3% up to 200 MPa. This essentially pressure-independent protein structure is an ideal system in which we can study interior hydration as a function of water activity, modulated here by pressure. As the pressure increases to 200 MPa, the principal change in the system is the increase in chemical potential of water,
µwat = 1.4 kBT (1 kBT
2.5 kJ/mol at room temperature).
To locate possible water in the cavity, we constructed the obser ved electron-density difference maps (Fo,hp Fo,0.1 MPa)
c,0.1 MPa (see Methods) at 100, 150, and 200 MPa. Maps averaged in real space over data at constant pressure are shown in Fig. 1. The maps were integrated over the cavity volume to obtain the average occupancies (Fig. 2). The principal feature in the maps is an increased electron density in the cavity, particularly at 150 and 200 MPa. At 100 MPa only a small peak is observed, corresponding to an average cavity occupancy of
0.5 water molecules (Fig. 2). This peak rapidly swells into a broad distribution at 150 MPa, increasing further to an average of approximately two water molecules in the cavity at 200 MPa. The average occupancy at 200 MPa together with the electron-density difference maps indicate a cooperative transition into a variety of configurations of at least two, and possibly as many as four, water molecules in the cavity.
To explore the energetics and dynamics of interior water, we performed MD simulations of solvated T4 lysozyme at pressures of 0.1, 100, and 200 MPa, with cavity occupancies between zero and five water molecules. By constructing what is effectively a grand-canonical partition function of cavity water, we obtain the pressure-dependent equilibrium occupancy distributions. As shown in Figs. 2 and 3, the calculations indicate a sharp transition from a predominantly empty cavity at pressures below
100 MPa to a cavity cooperatively filled by approximately four water molecules at pressures above
200 MPa. The calculated electron densities (Fig. 1) and average occupancies (Fig. 2) are in good agreement with the experiments after shifting the chemical potential of the bulk water phase (or, equivalently, the average interaction energy of water with the protein cavity) by
0.4 kBT (1 kJ/mol).
|
During two simulations a water molecule escaped the cavity, one after
0.5 ns in the occupancy N = 1 simulation and one in the N = 5 simulation after 0.9 ps, both at atmospheric pressure (105 Pa). In each case, water escaped through a transient opening between the side chains of Phe-114, Ser-116, Asn-132, and Leu-133. The escaping water molecule transiently occupied the crystallographic water site Wat-196. Water escape is coupled to protein dynamics, in particular, motions of Phe-114. Indeed, the Phe-114 side-chain crystallographic B factors, a measure of atomic fluctuations, are almost twice as large as those of the four other phenylalanines in the L99A crystal structure. The water-escape route from our simulations corroborates earlier NMR data that suggested a pathway for binding of nonpolar molecules near the F and G helices (19).
|
µwat
1.4 kBT at 293 K between 0.1 and 200 MPa) is sufficient to induce
50% filling of the cavity, on average. Such biphasic, environment-dependent behavior finds support in the coexistence of liquid and vapor water under ambient conditions (27). Consequently, filling under near ambient conditions could also be induced by lowering the average interaction energy of water with the surrounding protein through changes in the electrostatic environment of the cavity. The free energy of dissolving water into oil is remarkably high and largely an entropic penalty because of the formation of a water-sized cavity in oil. What contributes to the much more favorable free energy of transferring a water molecule into the protein? Our MD simulations suggest that a significant part of the free energy comes from hydrogen bonds between multiple water molecules in the cavity, explaining the cooperativity of the transition. This cooperativity is reflected in the average potential energy of water in the cavity, which decreases from 23 kJ/mol (N = 1) to 60 kJ/mol (N = 4) per molecule, or approximately 10 kJ/mol per waterwater hydrogen bond. Van der Waals interactions with the surrounding medium account for an additional approximately 10 kJ/mol per water molecule. Other electrostatic interactions with the protein should also contribute, as will entropy from the static-free volume of the preexisting cavity. The subtle thermodynamics of filling show that a simple description of the interior as hydrophobic fails to capture all of its relevant features.
|
Pressure-induced hydration may help explain why many proteins undergo large changes in activity under pressure (32). Interactions between water and protein hydrophobic surfaces is particularly important for molecular interfaces in ligand or substrate binding and complex formation (33) and in multidomain protein folding (30). Changes in the interior hydration also affect enzymatic function. Weakly screened electrostatic interactions in the protein interior provide a strong coupling between charge sites and internal water. Consequently, changes in protein conformation and charge states during enzyme turnover directly influence the free energy of water in the low-polarizability protein interior, with implications for protein function. Cytochrome P450 (11) and bacteriorhodopsin (12) are examples in which changes in redox states and amino acid protonation determine the presence or absence of functional water.
| Acknowledgements |
|---|
| Footnotes |
|---|
Conflict of interest statement: No conflicts declared.
Abbreviation: MD, molecular dynamics.
Data deposition: The atomic coordinates and structure factors have been deposited in the Protein Data Bank, www.pdb.org (PDB ID codes 2B6T [PDB] , 2B6W, 2B6X, 2B6Y, 2B6Z, 2B70, 2B72, 2B73, 2B74, and 2B75).
To whom correspondence may be addressed. E-mail: smg26{at}cornell.edu (regarding high-pressure crystallography), gerhard.hummer{at}nih.gov (regarding simulation), or brian{at}uoxray.uoregon.edu.
© 2005 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
R. Roth, D. Gillespie, W. Nonner, and R. E. Eisenberg Bubbles, Gating, and Anesthetics in Ion Channels Biophys. J., June 1, 2008; 94(11): 4282 - 4298. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Qvist, M. Davidovic, D. Hamelberg, and B. Halle From the Cover: A dry ligand-binding cavity in a solvated protein PNAS, April 29, 2008; 105(17): 6296 - 6301. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. Schlessman, C. Abe, A. Gittis, D. A. Karp, M. A. Dolan, and B. Garcia-Moreno E. Crystallographic Study of Hydration of an Internal Cavity in Engineered Proteins with Buried Polar or Ionizable Groups Biophys. J., April 15, 2008; 94(8): 3208 - 3216. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Giovambattista, C. F. Lopez, P. J. Rossky, and P. G. Debenedetti Hydrophobicity of protein surfaces: Separating geometry from chemistry PNAS, February 19, 2008; 105(7): 2274 - 2279. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ebbinghaus, S. J. Kim, M. Heyden, X. Yu, U. Heugen, M. Gruebele, D. M. Leitner, and M. Havenith An extended dynamical hydration shell around proteins PNAS, December 26, 2007; 104(52): 20749 - 20752. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Imai, S. Ohyama, A. Kovalenko, and F. Hirata Theoretical study of the partial molar volume change associated with the pressure-induced structural transition of ubiquitin Protein Sci., September 1, 2007; 16(9): 1927 - 1933. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Young, R. Abel, B. Kim, B. J. Berne, and R. A. Friesner Motifs for molecular recognition exploiting hydrophobic enclosure in protein-ligand binding PNAS, January 16, 2007; 104(3): 808 - 813. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. L. Quillin, P. T. Wingfield, and B. W. Matthews Determination of solvent content in cavities in IL-1beta using experimentally phased electron density PNAS, December 26, 2006; 103(52): 19749 - 19753. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. J. Roche, F. Guo, and J. M. Friedman Molecular Level Probing of Preferential Hydration and Its Modulation by Osmolytes through the Use of Pyranine Complexed to Hemoglobin J. Biol. Chem., December 15, 2006; 281(50): 38757 - 38768. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||