Abstract
We performed molecular dynamics simulations of the collapse of a two-domain protein, the BphC enzyme, into a globular structure to examine how water molecules mediate hydrophobic collapse of proteins. In the interdomain region, liquid water persists with a density 10to 15% lower than in the bulk, even at small domain separations. Water depletion and hydrophobic collapse occur on a nanosecond time scale, which is two orders of magnitude slower than that found in the collapse of idealized paraffin-like plates. When the electrostatic protein-water forces are turned off, a dewetting transition occurs in the interdomain region and the collapse speeds up by more than an order of magnitude. When attractive van der Waals forces are turned off as well, the dewetting in the interdomain region is more profound, and the collapse is even faster.
References
33
Referenced
454
-
K. Lum, D. Chandler, J. D. Weeks, J. Phys. Chem. B103, 4570 (1999).
(
10.1021/jp984327m
) / J. Phys. Chem. B (1999) -
D. M. Huang, D. Chandler, J. Phys. Chem. B106, 2047 (2002).
(
10.1021/jp013289v
) / J. Phys. Chem. B (2002) -
X. Huang, C. J. Margulis, B. J. Berne, Proc. Natl. Acad. Sci. U.S.A.100, 11953 (2003).
(
10.1073/pnas.1934837100
) / Proc. Natl. Acad. Sci. U.S.A. (2003) -
F. H. Stillinger, J. Solution Chem.2, 141 (1973).
(
10.1007/BF00651970
) / J. Solution Chem. (1973) -
D. Chandler, Nature417, 491 (2002).
(
10.1038/417491a
) / Nature (2002) -
P. Ball, Nature423, 25 (2003).
(
10.1038/423025a
) / Nature (2003) -
C. Pangali, M. Rao, B. J. Berne, J. Chem. Phys.71, 2982 (1979).
(
10.1063/1.438702
) / J. Chem. Phys. (1979) 10.1038/35102535
-
R. Steitzet al., Langmuir19, 2409 (2003).
(
10.1021/la026731p
) / Langmuir (2003) -
T. R. Jensenet al., Phys. Rev. Lett.90(8), 86101 (2003).
(
10.1103/PhysRevLett.90.086101
) / Phys. Rev. Lett. (2003) -
A. Wallqvist, B. J. Berne, J. Phys. Chem.99, 2893 (1995).
(
10.1021/j100009a053
) / J. Phys. Chem. (1995) -
A. Wallqvist, B. J. Berne, J. Phys. Chem.99, 2885 (1995).
(
10.1021/j100009a052
) / J. Phys. Chem. (1995) -
K. Leung, A. Luzar, J. Chem. Phys.113, 5845 (2000).
(
10.1063/1.1290479
) / J. Chem. Phys. (2000) -
K. Leung, A. Luzar, D. Bratko, Phys. Rev. Lett.90, 65502 (2003).
(
10.1103/PhysRevLett.90.065502
) / Phys. Rev. Lett. (2003) -
M. E. Paulaitis, L. R. Pratt, Adv. Protein Chem.62, 283 (2002).
(
10.1016/S0065-3233(02)62011-X
) / Adv. Protein Chem. (2002) -
L. R. Pratt, Annu. Rev. Phys. Chem.53, 409 (2002).
(
10.1146/annurev.physchem.53.090401.093500
) / Annu. Rev. Phys. Chem. (2002) -
G. Hummer, S. Garde, A. E. Garcia, L. R. Pratt, Chem. Phys.258, 349 (2000).
(
10.1016/S0301-0104(00)00115-4
) / Chem. Phys. (2000) -
K. Lum, A. Luzar, Phys. Rev. E.56, R6283 (1997).
(
10.1103/PhysRevE.56.R6283
) / Phys. Rev. E. (1997) -
D. M. Huang, D. Chandler, Proc. Natl. Acad. Sci. U.S.A.97, 8324 (2000).
(
10.1073/pnas.120176397
) / Proc. Natl. Acad. Sci. U.S.A. (2000) -
P. R. ten Wolde, D. Chandler, Proc. Natl. Acad. Sci. U.S.A.99, 6539 (2002).
(
10.1073/pnas.052153299
) / Proc. Natl. Acad. Sci. U.S.A. (2002) -
K. Lum, D. Chandler, Int. J. Thermophys.19, 845 (1998).
(
10.1023/A:1022643006849
) / Int. J. Thermophys. (1998) -
R. D. Mountain, D. Thirumalai, J. Am. Chem. Soc.125, 1950 (2003).
(
10.1021/ja020496f
) / J. Am. Chem. Soc. (2003) -
R. Zhou, B. D. Silverman, A. Royyuru, P. Athma, Proteins52, 561 (2003).
(
10.1002/prot.10419
) / Proteins (2003) - The starting structure is taken from the crystal structure deposited in PDB (entry 1dhy.pdb). The interdomain distance of the crystal structure is increased by D along the direction of two domain centers of geometry to create “gaps” between the two domains (to make room for water). Various distances ranging from D = 2.5 Å to 20.0 Å are studied. The resulting protein configurations are then solvated in a water box with water layers at least 8 Å from the protein surfaces. Figure 1 shows one such solvated configuration. Eight Na + counterions are added to make the system electrically neutral. The solvated protein systems have up to 42 000 atoms (the actual size varies for different interdomain distances). The GROMACS simulation package is used here for this large system because of its fast speed ( 25 ). Each NPT MD simulation [( 26 ) 1 atm and 300 K] is run in parallel with 8 to 16 processors on IBM SP2-Power3-375 MHz clusters. The OPLSAA force field is used for the protein ( 27 ) and a simple point charge (SPC) water model is used for the explicit solvent. For the long-range electrostatic interactions the particle-mesh Ewald method is used. For the van der Waals interactions a typical 10 Å cutoff is used. A time step of 2.0 fs is used with bond lengths constrained. A standard equilibration procedure which includes a conjugate gradient minimization and a 100-ps MD simulation with position restraints is followed to equilibrate each solvated system. The final configurations from equilibration are then used for data collection. We have modified GROMACS to allow us to selectively turn off certain terms in the force field in particular the protein-water electrostatic potential the attractive part of the protein-water van der Waals potential the domain-domain electrostatic potential and the attractive part of the domain-domain van der Waals potential. This approach allows us to make a sensitivity analysis of hydrophobic collapse and to compare the protein collapse with what is already known about the hydrophobic collapse of the plates.
-
E. Lindahl, B. Hess, D. van der Spoel, J. Mol. Model.7, 306 (2001).
(
10.1007/s008940100045
) / J. Mol. Model. (2001) - One caveat is that NPT MD introduces fictitious dynamics from both the temperature and pressure control. New fictitious degrees of freedom are introduced with corresponding kinetic energies fictitious masses and fictitious frequencies. These influence the time scales and can affect the rate processes observed. Nevertheless the thermodynamic equilibrium water densities inside the interplate or interdomain region and the water depletion or dewetting phenomena should stay the same. Our observations about kinetics are borne out in careful studies.
-
W. L. Jorgensen, D. Maxwell, J. Tirado-Rives, J. Am. Chem. Soc.118, 11225 (1996).
(
10.1021/ja9621760
) / J. Am. Chem. Soc. (1996) - The original total interdomain distance is about 18.8 Å so after the enlargement the total interdomain distance is about 22.8 Å. For simplicity we will refer to it as interdomain gap distance D = 4 Å in the following.
- It is tricky to define the water density inside the domain gap region particularly for cases with very small D values because the fluctuation in the number of water molecules can be fairly large compared with the total number of water molecules inside the region. In general it is very difficult to directly measure the volume inside the domain gap region because of the irregularity of the protein domain structures; thus we used a comparative way to measure the water density.
-
Y. Cheng, P. J. Rossky, Nature392, 696 (1998).
(
10.1038/33653
) / Nature (1998) -
H. S. Ashbaugh, M. E. Paulaitis, J. Am. Chem. Soc.123, 10721 (2001).
(
10.1021/ja016324k
) / J. Am. Chem. Soc. (2001) - Single-letter abbreviations for the amino acid residues are as follows: A Ala; C Cys; D Asp; E Glu; F Phe; G Gly; H His; I Ile; K Lys; L Leu; M Met; N Asn; P Pro; Q Gln; R Arg; S Ser; T Thr; V Val; W Trp; and Y Tyr.
- This work was partially supported by a grant to B.J.B. from NIH (GM4330) and a Shared University Research (SUR) grant from the IBM Corporation to Columbia University.
Dates
Type | When |
---|---|
Created | 20 years, 11 months ago (Sept. 9, 2004, 5:17 p.m.) |
Deposited | 1 year, 7 months ago (Jan. 9, 2024, 9:05 p.m.) |
Indexed | 1 month ago (July 25, 2025, 5:55 a.m.) |
Issued | 20 years, 11 months ago (Sept. 10, 2004) |
Published | 20 years, 11 months ago (Sept. 10, 2004) |
Published Print | 20 years, 11 months ago (Sept. 10, 2004) |
@article{Zhou_2004, title={Hydrophobic Collapse in Multidomain Protein Folding}, volume={305}, ISSN={1095-9203}, url={http://dx.doi.org/10.1126/science.1101176}, DOI={10.1126/science.1101176}, number={5690}, journal={Science}, publisher={American Association for the Advancement of Science (AAAS)}, author={Zhou, Ruhong and Huang, Xuhui and Margulis, Claudio J. and Berne, Bruce J.}, year={2004}, month=sep, pages={1605–1609} }