Are buckyballs hydrophobic?

Buckyballs exhibit two seemingly opposing characters. On one hand, they are known to be insoluble in water with all potential chemical properties of being hydrophobic. On the other hand, their pairwise effective interaction in water includes a repulsive solvent-induced contribution. We perform molecular dynamics simulations of the association process of two C60 fullerenes in water at different temperatures in order to reconcile these contradicting observations. For comparison, the simulations were also performed in a nonpolar solvent, and the results were further contrasted with those obtained previously for the association of graphene sheets. Considering the association in water, we find small magnitudes for the enthalpy and entropy changes with small positive slopes as a function of temperature, implying an almost negligible change in the heat capacity at constant pressure. These findings sharply contradict the behavior of typical hydrophobic interactions. The reason for these abnormalities, as well as for the repulsive nature of the solvent-induced interactions, is the shape of the contact state that supports the existence of a distinct type of interfacial waters located between two convex surfaces of two buckyballs. These interfacial waters are characterized by smaller entropy and lower density and form a smaller number of hydrogen bonds with surrounding waters compared with those of the interfacial waters around the dissociated solutes. Thus, upon bringing the C60 fullerenes into contact, the changes associated with the liberation of the latter to bulk waters are opposed by the concomitant conversion of the latter also to the distinct waters between the two spherical solutes. We argue that although the effective pair interaction is not hydrophobic, the solvation properties are hydrophobic. In the hydration free energy of a single solute, there is no contact state. Furthermore, in the solubility at the macroscopic scale, the relative number of these distinct waters around a large aggregate is smaller, and hence, they are not predicted to influence much the solvation properties. Therefore, buckyballs can serve as an example in which hydrophobic interaction cannot be deducted from hydrophobic solvation.

[1]  R. Zangi,et al.  Ribose 2'-Hydroxyl Groups Stabilize RNA Hairpin Structures Containing GCUAA Pentaloop. , 2013, Journal of chemical theory and computation.

[2]  Brad A. Bauer,et al.  Free energetics of carbon nanotube association in pure and aqueous ionic solutions. , 2012, The journal of physical chemistry. B.

[3]  S. Garde,et al.  Hydrophobicity of proteins and interfaces: insights from density fluctuations. , 2011, Annual review of chemical and biomolecular engineering.

[4]  R. Zangi Driving force for hydrophobic interaction at different length scales. , 2011, The journal of physical chemistry. B.

[5]  B. Farouk,et al.  Molecular dynamics simulations of the interactions and dispersion of carbon nanotubes in polyethylene oxide/water systems , 2011 .

[6]  G. Graziano On the pairwise hydrophobic interaction of fullerene , 2010 .

[7]  A. Liwo,et al.  Potential of mean force of association of large hydrophobic particles: toward the nanoscale limit. , 2010, The journal of physical chemistry. B.

[8]  B. Berne,et al.  Temperature dependence of dimerization and dewetting of large-scale hydrophobes: a molecular dynamics study. , 2008, The journal of physical chemistry. B.

[9]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[10]  B. Berne,et al.  Effect of ions on the hydrophobic interaction between two plates. , 2007, Journal of the American Chemical Society.

[11]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[12]  Jason R. Brown,et al.  Single-molecule force spectroscopy measurements of "hydrophobic bond" between tethered hexadecane molecules. , 2006, The journal of physical chemistry. B.

[13]  Dmitry Bedrov,et al.  Water-induced interactions between carbon nanoparticles. , 2006, The journal of physical chemistry. B.

[14]  D. Bedrov,et al.  A molecular-dynamics simulation study of solvent-induced repulsion between C60 fullerenes in water. , 2005, The Journal of chemical physics.

[15]  M. Sasai,et al.  Fluctuating hydration structure around nanometer-size hydrophobic solutes. I. Caging and drying around C60 and C60H60 spheres. , 2005, The journal of physical chemistry. B.

[16]  J. Engberts,et al.  Physisorption of hydroxide ions from aqueous solution to a hydrophobic surface. , 2005, Journal of the American Chemical Society.

[17]  B Montgomery Pettitt,et al.  On the mechanism of hydrophobic association of nanoscopic solutes. , 2005, Journal of the American Chemical Society.

[18]  D. Bedrov,et al.  Repulsive solvent-induced interaction between C60 fullerenes in water. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[19]  Greg L. Hura,et al.  Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. , 2004, The Journal of chemical physics.

[20]  D. Paschek Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: a detailed structural and energetical analysis based on molecular dynamics simulations. , 2004, The Journal of chemical physics.

[21]  D. Paschek Temperature dependence of the hydrophobic hydration and interaction of simple solutes: an examination of five popular water models. , 2003, The Journal of chemical physics.

[22]  B. Berne,et al.  Dewetting-induced collapse of hydrophobic particles , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[23]  Petros Koumoutsakos,et al.  On the Water−Carbon Interaction for Use in Molecular Dynamics Simulations of Graphite and Carbon Nanotubes , 2003 .

[24]  Lukas D. Schuler,et al.  On the Choice of Dihedral Angle Potential Energy Functions for n-Alkanes , 2000 .

[25]  Schuster,et al.  Synthesis and characterization of water-soluble amino fullerene derivatives , 2000, Organic letters.

[26]  D. Chandler,et al.  Hydrophobicity at Small and Large Length Scales , 1999 .

[27]  E. E. Braudo,et al.  Heat capacities of aqueous solutions of amino acid and dipeptide derivatives of fullerene , 1998 .

[28]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[29]  R. Abseher,et al.  The Temperature-Dependence of Hydrophobic Association in Water. Pair versus Bulk Hydrophobic Interactions , 1997 .

[30]  García,et al.  Origin of Entropy Convergence in Hydrophobic Hydration and Protein Folding. , 1996, Physical review letters.

[31]  Anders Wallqvist,et al.  Computer Simulation of Hydrophobic Hydration Forces on Stacked Plates at Short Range , 1995 .

[32]  B. Berne,et al.  Molecular Dynamics Study of the Dependence of Water Solvation Free Energy on Solute Curvature and Surface Area , 1995 .

[33]  M. Costas,et al.  Thermodynamics of the hydrophobic effect in surfactant solutions: Micellization and adsorption , 1995 .

[34]  Y. Guissani,et al.  A computer simulation study of the temperature dependence of the hydrophobic hydration , 1993 .

[35]  J. Engberts,et al.  Hydrophobic Effects. Opinions and Facts , 1993 .

[36]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[37]  Rodney S. Ruoff,et al.  Solubility of fullerene (C60) in a variety of solvents , 1993 .

[38]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[39]  K A Dill,et al.  The meaning of hydrophobicity. , 1990, Science.

[40]  K. Dill Dominant forces in protein folding. , 1990, Biochemistry.

[41]  R. L. Baldwin,et al.  Temperature dependence of the hydrophobic interaction in protein folding. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[42]  S. Gill,et al.  Anomalous heat capacity of hydrophobic solvation , 1985 .

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

[44]  S. Gill,et al.  An equation of state describing hydrophobic interactions. , 1976, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Frank H. Stillinger,et al.  Structure in aqueous solutions of nonpolar solutes from the standpoint of scaled-particle theory , 1973 .

[46]  Harold A. Scheraga,et al.  Structure of Water and Hydrophobic Bonding in Proteins. II. Model for the Thermodynamic Properties of Aqueous Solutions of Hydrocarbons , 1962 .

[47]  Henry S. Frank,et al.  Free Volume and Entropy in Condensed Systems III. Entropy in Binary Liquid Mixtures; Partial Molal Entropy in Dilute Solutions; Structure and Thermodynamics in Aqueous Electrolytes , 1945 .

[48]  J. Edsall Apparent Molal Heat Capacities of Amino Acids and Other Organic Compounds , 1935 .

[49]  D. Heymann Solubility of C60 in alcohols and alkanes , 1996 .

[50]  Norbert Muller,et al.  Search for a realistic view of hydrophobic effects , 1990 .

[51]  P. Privalov,et al.  Stability of protein structure and hydrophobic interaction. , 1988, Advances in protein chemistry.

[52]  Alfons Geiger,et al.  Molecular dynamics study of the hydration of Lennard‐Jones solutes , 1979 .

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

[54]  D. D. Eley On the solubility of gases. Part II.—A comparison of organic solvents with water , 1939 .

[55]  D. D. Eley On the solubility of gases. Part I.—The inert gases in water , 1939 .