Hydrophobic collapse in multidomain protein folding

Science

Volumn 305, Issue 5690, 2004, Pages 1605-1609

  Zhou, Ruhong ^a,b   Huang, Xuhui ^b   Margulis, Claudio J ^b   Berne, Bruce J ^a,b  

^a IBM T J WATSON RESEARCH CENTER   (United States)

^b Columbia University ^*  (United States)

Author keywords

[No Author keywords available]

Indexed keywords

COMPUTER SIMULATION; DEWATERING; ELECTROSTATICS; HYDROPHOBICITY; MOLECULAR DYNAMICS; PARAFFINS; VAN DER WAALS FORCES;

MULTIDOMAIN PROTEIN FOLDING; PARAFFIN-LIKE PLATES;

PROTEINS;

GLOBULAR PROTEIN; PARAFFIN; POLYCHLORINATED BIPHENYL; PROTEIN BPHC; UNCLASSIFIED DRUG; WATER; 2,3 DIHYDROXYBIPHENYL OXYGENASE; 2,3-DIHYDROXYBIPHENYL OXYGENASE; DIOXYGENASE; OXYGENASE;

BIOLOGY;

ARTICLE; HYDROPHOBICITY; MOLECULAR DYNAMICS; PRIORITY JOURNAL; PROTEIN AGGREGATION; PROTEIN DOMAIN; PROTEIN FOLDING; SOLVATION; CHEMICAL STRUCTURE; CHEMISTRY; COMPUTER SIMULATION; ELECTRICITY; KINETICS; PROTEIN CONFORMATION; PROTEIN TERTIARY STRUCTURE; SURFACE PROPERTY;

COMPUTER SIMULATION; DIOXYGENASES; ELECTROSTATICS; HYDROPHOBICITY; KINETICS; MODELS, MOLECULAR; OXYGENASES; PROTEIN CONFORMATION; PROTEIN FOLDING; PROTEIN STRUCTURE, TERTIARY; SURFACE PROPERTIES; WATER;

EID: 4644236471     PISSN: 00368075     EISSN: None     Source Type: Journal    
DOI: 10.1126/science.1101176     Document Type: Article

References (33)

1. K. Lum, D. Chandler, J. D. Weeks, J. Phys. Chem. B 103, 4570 (1999).
2. D. M. Huang, D. Chandler J. Phys. Chem. B 106, 2047 (2002).
3. X. Huang, C. J. Margulis, B. J. Berne, Proc. Natl. Acad. Sci. U.S.A. 100, 11953 (2003).
4. F. H. Stillinger, J. Solution Chem. 2, 141 (1973).
5. D. Chandler, Nature 417, 491 (2002).
6. P. Ball, Nature 423, 25 (2003).
7. C. Pangali, M. Rao, B. J. Berne, J. Chem. Phys. 71, 2982 (1979).
8. G. Hummer, J. C. Rasaiah, J. P. Noworyta, Nature 414, 188 (2001).
9. R. Steitz et al., Langmuir 19, 2409 (2003).
10. T. R. Jensen et al., Phys. Rev. Lett. 90(8), 86101 (2003).
11. A. Wallqvist, B. J. Berne, J. Phys. Chem. 99, 2893 (1995).
12. A. Wallqvist B. J. Berne, J. Phys. Chem. 99, 2885 (1995).
13. K. Leung, A. Luzar, J. Chem. Phys. 113, 5845 (2000).
14. K. Leung, A. Luzar, D. Bratko, Phys. Rev. Lett. 90, 65502 (2003).
15. M. E. Paulaitis, L. R. Pratt, Adv. Protein Chem. 62, 283 (2002).
16. L. R. Pratt, Annu. Rev. Phys. Chem. 53, 409 (2002).
17. G. Hummer, S. Garde, A. E. Garcia, L. R. Pratt, Chem. Phys. 258, 349 (2000).
18. K. Lum, A. Luzar, Phys. Rev. E. 56, R6283 (1997).
19. D. M. Huang, D. Chandler, Proc. Natl. Acad. Sci. U.S.A. 97, 8324 (2000).
20. P. R. ten Wolde, D. Chandler, Proc. Natl. Acad. Sci. U.S.A. 99, 6539 (2002).
21. K. Lum, D. Chandler, Int. J. Thermophys. 19, 845 (1998).
22. R. D. Mountain, D. Thirumalai, J. Am. Chem. Soc. 125, 1950 (2003).
23. R. Zhou, B. D. Silverman, A. Royyuru, P. Athma, Proteins 52, 561 (2003).
24. + 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 5P2-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.
25. E. Lindahl, B. Hess, D. van der Spoel, J. Mol. Model. 7, 306 (2001).
26. 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.
27. W. L. Jorgensen, D. Maxwell, J. Tirado-Rives, J. Am. Chem. Soc. 118, 11225 (1996).
28. The original total interdomain distance is about 1B.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.
29. 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.
30. Y. Cheng, P. J. Rossky, Nature 392, 696 (1998).
31. H. S. Ashbaugh, M. E. Paulaitis, J. Am. Chem. Soc. 123, 10721 (2001).
32. 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.
33. 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.