Free-energy analysis of the molecular binding into lipid membrane with the method of energy representation

Journal of Chemical Physics

Volumn 128, Issue 19, 2008, Pages

  Matubayasi, Nobuyuki ^a,b   Shinoda, Wataru ^c   Nakahara, Masaru ^a  

^a KYOTO UNIVERSITY   (Japan)

^b JAPAN SCIENCE AND TECHNOLOGY AGENCY   (Japan)

^c NATIONAL INSTITUTE OF ADVANCED INDUSTRIAL SCIENCE AND TECHNOLOGY AIST   (Japan)

Author keywords

[No Author keywords available]

Indexed keywords

CHEMICAL BONDS; COMPUTER SIMULATION; FREE ENERGY; HYDROPHOBICITY; MOLECULAR DYNAMICS; SOLVATION;

ETHYLBENZENE SOLUTES; KIRKWOOD-BUFF THEORY; MOLECULAR SIMULATION; STATISTICAL-MECHANICAL TREATMENT;

LIPID BILAYERS;

CARBON MONOXIDE; DIMYRISTOYLPHOSPHATIDYLCHOLINE; LIPID;

ARTICLE; BIOPHYSICS; CHEMICAL MODEL; CHEMICAL STRUCTURE; CHEMISTRY; COMPUTER SIMULATION; ENTROPY; MEMBRANE; METHODOLOGY; MICELLE; PHYSICAL CHEMISTRY; SOLUBILITY; SOLUTION AND SOLUBILITY; STATISTICAL MODEL; THERMODYNAMICS;

BIOPHYSICS; CARBON MONOXIDE; CHEMISTRY, PHYSICAL; COMPUTER SIMULATION; DIMYRISTOYLPHOSPHATIDYLCHOLINE; ENTROPY; LIPIDS; MEMBRANES; MICELLES; MODELS, CHEMICAL; MODELS, MOLECULAR; MODELS, STATISTICAL; SOLUBILITY; SOLUTIONS; THERMODYNAMICS;

EID: 44349122912     PISSN: 00219606     EISSN: None     Source Type: Journal    
DOI: 10.1063/1.2919117     Document Type: Article

References (89)

1. The Structure of Biological Membranes, 2nd ed., edited by, P. L. Yeagle, (CRC, Boca Raton, FL, 2005).
2. T. Hanai and D. A. Haydon, J. Theor. Biol. 11, 370 (1966). 0022-5193
3. Y. Katz and J. M. Diamond, J. Membr. Biol. 0022-2631 10.1007/BF01870175 17, 101 (1974).
4. A. Finkelstein, J. Gen. Physiol. 0022-1295 10.1085/jgp.68.2.127 68, 127 (1976).
5. H. Kamaya, S. Kaneshina, and I. Ueda, Biochim. Biophys. Acta 646, 135 (1981). 0006-3002
6. S. A. Simon, R. V. McDaniel, and T. J. McIntosh, J. Phys. Chem. 0022-3654 10.1021/j100397a048 86, 1449 (1982).
7. L. R. De Young and K. A. Dill, Biochemistry 0006-2960 10.1021/bi00414a050 27, 5281 (1988).
8. Y. Xu and P. Tang, Biochim. Biophys. Acta 1323, 154 (1997). 0006-3002
9. P. Tang, B. Yan, and Y. Xu, Biophys. J. 72, 1676 (1997). 0006-3495
10. W. H. J. Vaes, E. U. Ramos, H. J. M. Verhaar, C. J. Cramer, and J. L. M. Hermens, Chem. Res. Toxicol. 11, 847 (1998). 0893-228X
11. P. Westh and C. Trandum, J. Phys. Chem. B 1089-5647 10.1021/jp001540l 104, 11334 (2000).
12. H. Patel, T. W. Schultz, and M. T. D. Cronin, J. Mol. Struct.: THEOCHEM 0166-1280 593, 9 (2002).
13. T. -X. Xiang and B. D. Anderson, Adv. Drug Deliv. Rev. 58, 1357 (2006). 0169-409X
14. S. -J. Marrink and H. J. C. Berendsen, J. Phys. Chem. 0022-3654 10.1021/j100066a040 98, 4155 (1994).
15. D. Bassolino-Klimas, H. E. Alper, and T. R. Stouch, J. Am. Chem. Soc. 0002-7863 10.1021/ja00119a028 117, 4118 (1995).
16. S. J. Marrink and H. J. C. Berendsen, J. Phys. Chem. 0022-3654 10.1021/jp952956f 100, 16729 (1996).
17. A. Pohorille and M. A. Wilson, J. Chem. Phys. 0021-9606 10.1063/1.471030 104, 3760 (1996).
18. M. L. Huertas, V. Cruz, J. J. L. Cascales, A. U. Acua, and J. G. de la Torre, Biophys. J. 71, 1428 (1996). 0006-3495
19. B. Roux, Biophys. J. 77, 139 (1999). 0006-3495
20. P. Jedlovszky and M. Mezei, J. Am. Chem. Soc. 0002-7863 10.1021/ja000156z 122, 5125 (2000).
21. L. Koubi, M. Tarek, M. L. Klein, and D. Scharf, Biophys. J. 78, 800 (2000). 0006-3495
22. T. -X. Xiang and B. D. Anderson, Biophys. J. 82, 2052 (2002). 0006-3495
23. D. Bemporad, J. W. Essex, and C. Luttmann, J. Phys. Chem. B 1089-5647 10.1021/jp035260s 108, 4875 (2004).
24. W. Shinoda, M. Mikami, T. Baba, and M. Hato, J. Phys. Chem. B 1089-5647 10.1021/jp035998+ 108, 9346 (2004).
25. J. L. MacCallum and D. P. Tieleman, J. Am. Chem. Soc. 128, 125 (2006). 0002-7863
26. K. E. Norman and H. Nymeyer, Biophys. J. 91, 2046 (2006). 0006-3495
27. S. Vemparala, L. Saiz, R. G. Eckenhoff, and M. L. Klein, Biophys. J. 0006-3495 10.1529/biophysj.106.085324 91, 2815 (2006).
28. L. Li, H. Davande, D. Bedrov, and G. D. Smith, J. Phys. Chem. B 111, 4067 (2007). 1520-6106
29. N. Matubayasi and M. Nakahara, J. Chem. Phys. 0021-9606 10.1063/1.1309013 113, 6070 (2000).
30. N. Matubayasi and M. Nakahara, J. Chem. Phys. 0021-9606 10.1063/1.1495850 117, 3605 (2002); N. Matubayasi and M. Nakahara, J. Chem. Phys. 0021-9606 10.1063/1.1533752 118, 2446 (2003).
31. N. Matubayasi and M. Nakahara, J. Chem. Phys. 0021-9606 10.1063/1.1613938 119, 9686 (2003).
32. H. Takahashi, N. Matubayasi, M. Nakahara, and T. Nitta, J. Chem. Phys. 0021-9606 10.1063/1.1774981 121, 3989 (2004).
33. N. Matubayasi, K. K. Liang, and M. Nakahara, J. Chem. Phys. 0021-9606 10.1063/1.2186324 124, 154908 (2006).
34. S. E. Feller, Curr. Opin. Colloid Interface Sci. 1359-0294 10.1016/S1359-0294(00)00058-3 5, 217 (2000).
35. H. L. Scott, Curr. Opin. Struct. Biol. 0959-440X 10.1016/S0959-440X(02) 00353-6 12, 495 (2002).
36. R. W. Pastor, R. M. Venable, and S. E. Feller, Acc. Chem. Res. 0001-4842 10.1021/ar0100529 35, 438 (2002).
37. L. Saiz and M. L. Klein, Acc. Chem. Res. 0001-4842 10.1021/ar010167c 35, 482 (2002).
38. P. van der Ploeg and H. J. C. Berendsen, J. Chem. Phys. 0021-9606 10.1063/1.443321 76, 3271 (1982).
39. H. Heller, M. Schaefer, and K. Schulten, J. Phys. Chem. 0022-3654 10.1021/j100133a034 97, 8343 (1993).
40. K. Tu, D. J. Tobias, and M. L. Klein, Biophys. J. 69, 2558 (1995). 0006-3495
41. S. E. Feller, R. M. Venable, and R. W. Pastor, Langmuir 0743-7463 10.1021/la970746j 13, 6555 (1997).
42. W. Shinoda, N. Namiki, and S. Okazaki, J. Chem. Phys. 0021-9606 10.1063/1.473592 106, 5731 (1997).
43. K. Tu, M. Tarek, M. L. Klein, and D. Scharf, Biophys. J. 75, 2123 (1998). 0006-3495
44. S. E. Feller and R. W. Pastor, J. Chem. Phys. 0021-9606 10.1063/1.479313 111, 1281 (1999).
45. S. A. Pandit, D. Bostick, and M. L. Berkowitz, Biophys. J. 84, 3743 (2003). 0006-3495
46. J. N. Sachs, H. Nanda, H. I. Petrache, and T. B. Woolf, Biophys. J. 0006-3495 10.1529/biophysj.103.035816 86, 3772 (2004).
47. M. Patra, E. Salonen, E. Terama, I. Vattulainen, R. Faller, B. W. Lee, J. Holopainen, and M. Karttunen, Biophys. J. 90, 1121 (2006). 0006-3495
48. U. R. Pedersen, G. H. Peters, and P. Westh, Biophys. Chem. 125, 104 (2007). 0301-4622
49. M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford University Press, Oxford, 1987).
50. The drawback of the approach is that the solvation free energy is evaluated from an approximate functional.
51. A. Ben-Naim, Solvation Thermodynamics (Plenum, New York, 1987).
52. The one-to-one association is implied in Eq.. This is because both the membrane and solute are assumed to be dilute in the present formulation. When more than one solute molecule is involved in the membrane solution, it suffices to extend the development straightforwardly by viewing one of the solute molecules as the "solute" species and the others as the "solvent" species constituting the "mixed solvent" system.
53. The extension of the formulation is straightforward when the system contains cosolvents.
54. The condition of a fixed (total) chemical potential of the solute is realized, for example, when the counterpart of the equilibrium of the solute insertion is common to the membrane solution and neat solvent.
55. When S (0) (Km + Vm) is much larger than nA, Kx is simply close to Cw /S (0) and is determined essentially by the solubility in the neat solvent. Indeed, when the affinity of the solute to the membrane is strong enough, the mole fraction in the membrane is at the order of 1 and the order of Kx is given by the solubility in neat water for a sparingly soluble solute. In this case, Kx (or Δ Gx) is not very informative for the solute-membrane correlation.
56. J. G. Kirkwood and F. P. Buff, J. Chem. Phys. 0021-9606 10.1063/1.1748352 19, 774 (1951).
57. I. L. Shulgin and E. Ruckenstein, J. Chem. Phys. 0021-9606 10.1063/1.2011388 123, 054909 (2005).
58. In Eq., one membrane aggregate is treated as a single particle. The configuration fluctuates within the aggregate, and the coordinates for the intra-aggregate motion are incorporated in the definition of x.
59. Actually, an arbitrary position in the solute molecule can be taken as the molecular center. The Gsm /A value calculated from Eq. or is independent of the choice of the molecular center of the solute.
60. N. Matubayasi and R. M. Levy, J. Phys. Chem. 0022-3654 10.1021/jp951618b 100, 2681 (1996).
61. A. D. MacKerell, Jr., D. Bashford, M. Bellott, R. L. Dunbrack, Jr., J. D. Evanseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, D. Joseph-McCarthy, L. Kuchnir, K. Kuczera, F. T. K. Lau, C. Mattos, S. Michnick, T. Ngo, D. T. Nguyen, B. Prodhom, W. E. Reiher III, B. Roux, M. Schlenkrich, J. C. Smith, R. Stote, J. Straub, M. Watanabe, J. Wiórkiewicz-Kuczera, D. Yin, and M. Karplus, J. Phys. Chem. B 1089-5647 10.1021/jp973084f 102, 3586 (1998).
62. J. B. Klauda, B. R. Brooks, A. D. MacKerell, Jr., R. M. Venable, and R. W. Pastor, J. Phys. Chem. B 1089-5647 10.1021/jp0468096 109, 5300 (2005).
63. W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 0021-9606 10.1063/1.445869 79, 926 (1983).
64. N. Kučerka, Y. Liu, N. Chu, H. I. Petrache, S. Tristram-Nagle, and J. F. Nagle, Biophys. J. 0006-3495 10.1529/biophysj.104.056606 88, 2626 (2005).
65. W. Shinoda and M. L. Klein, manuscript in preparation. The equilibration scheme adopted is summarized as follows. First, 128 DMPC molecules were placed in lattice form with the conformation at the crystal phase. Then 3405 water molecules were added, and the system was equilibrated through an 80 ns molecular dynamics simulation in the N Pn AT ensemble with the membrane area corresponding to the experimental value (Ref.). Finally, the number of water molecules was increased to 3886 and the equilibration was further performed for 6 ns.
66. G. J. Martyna, M. E. Tuckerman, D. J. Tobias, and M. L. Klein, Mol. Phys. 0026-8976 10.1080/00268979650027054 87, 1117 (1996).
67. B. R. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus, J. Comput. Chem. 0192-8651 10.1002/jcc.540040211 4, 187 (1983).
68. W. Rawicz, K. C. Olbrich, T. McIntosh, D. Needham, and E. Evans, Biophys. J. 79, 328 (2000). 0006-3495
69. The area compressibility modulus is given by kB T 〈A〉 / 〈 δ A2 〉, where 〈A〉 is the average of the total area of the membrane and 〈 δ A2 〉 is its mean square deviation (Ref.). When the experimental value of 234 dyne/cm is used for the modulus (Ref.) and the membrane area is considered to be ∼60 Å2 per lipid, 〈 δ A2 〉 is ∼80 Å2 with 128 DMPC molecules and is larger than the area expected to be occupied by the solute treated in the present work. According to the membrane simulations involving alcoholic or anesthetic solute (Refs.), the membrane area changes by a few percent when the mole ratio of the solute to the phospholipid is ∼0.1. In our simulation, the mole ratio is 1/128 (a single solute molecule in 128 DMPC molecules). The area expansion or contraction is thus negligible upon introduction of CO, CO2, benzene, or ethylbenzene as the solute; the area can be safely taken to be identical between the solution and reference solvent systems.
70. M. Parrinello and A. Rahman, Phys. Rev. Lett. 0031-9007 10.1103/PhysRevLett.45.1196 45, 1196 (1980).
71. G. J. Martyna, D. J. Tobias, and M. L. Klein, J. Chem. Phys. 0021-9606 10.1063/1.467468 101, 4177 (1994).
72. G. J. Martyna, M. L. Klein, and M. Tuckerman, J. Chem. Phys. 0021-9606 10.1063/1.463940 97, 2635 (1992).
73. W. Shinoda and M. Mikami, J. Comput. Chem. 0192-8651 10.1002/jcc.10249 24, 920 (2003).
74. For all the solutes and the regions of solute position, ρ k,0 e and kl,0 e can be constructed from a single trajectory of the membrane system without the solute since the solute insertion as a test particle does not affect the trajectory. This is useful to reduce the computational demand when a multiple set of solvations is to be treated against a common (mixed) solvent system. A similar remark applies to the calculation in neat water.
75. H. C. Andersen, J. Comput. Phys. 0021-9991 10.1016/0021-9991(83)90014-1 52, 24 (1983).
76. The free energies presented in Ref. are different from those in this work. According to the notation in Eq. (2) of Ref., Jedlovszky and Mezei showed A′ in their paper while we show μ′.
77. The cavity-insertion Widom method is highly efficient when it works. When the solute size is large, however, an appropriate cavity is not found in the molecular simulation system and the method does not work. For the DMPC-water system, benzene and ethylbenzene are too large and are beyond the applicability of the cavity-insertion Widom method. In the approximate procedure described in Secs., the overlapping configuration produced upon the solute insertion provides nonzero ρ k,0 e and kl,0 e at large energy coordinates and contributes to the solvation free energy as the excluded-volume effect.
78. When the solvation free energy is given by Δμ (z) as a continuous function of z and the region of interest is identified as A
79. The contributions to Kx from regions V-VI are 0.09%, 0.7%, 0.04%, and 0.001% in weight for CO, CO2, benzene, and ethylbenzene, respectively.
80. The discussion concerning Δμ is based on the approximation expressed as Eqs.. The arguments concerning 〈u〉 are exact within the statistical error, on the other hand, for the potential functions employed.
81. N. Matubayasi and M. Nakahara, J. Chem. Phys. 0021-9606 10.1063/1.1849165 122, 074509 (2005).
82. As shown in Refs., the density dependence is steep for the contribution from the repulsive component of the solute-solvent interaction to the solvation free energy, while it is milder for the contribution from the attractive component. Correspondingly, when the water density decreases in the interfacial region of the membrane, the reduction of the repulsive effect is stronger and the solvation free energy can be negative.
83. The "excluded-volume" region can be identified over the energy coordinate as > c, where c is rather arbitrary and is typically 10-20 kcal/mol. When the integral over in Eqs. is restricted to the domain corresponding to the excluded-volume region, the integral value is larger in the order of bulk water, region IV, and region II.
84. N. Matubayasi, L. H. Reed, and R. M. Levy, J. Phys. Chem. 0022-3654 10.1021/j100092a040 98, 10640 (1994).
85. The change in the solvent-solvent interaction energy contributes to the thermodynamic enthalpy and energy. The excess partial molar enthalpy is the sum of the change in the intramolecular energy, 〈u〉, the solvent reorganization term, and the pressure-volume term (Ref.).
86. Handbook of Surface and Colloid Chemistry, 2nd ed., edited by, K. S. Birdi, (CRC, Boca Raton, FL, 2003).
87. F. Figueirido, G. S. Del Buono, and R. M. Levy, J. Chem. Phys. 0021-9606 10.1063/1.470721 103, 6133 (1995).
88. G. Hummer, L. R. Pratt, and A. E. García, J. Phys. Chem. 0022-3654 10.1021/jp951011v 100, 1206 (1996).
89. P. H. Hünenberger and J. A. McCammon, J. Chem. Phys. 0021-9606 10.1063/1.477873 110, 1856 (1999).