Understanding the role of temperature change and the presence of NaCl salts on caffeine aggregation in aqueous solution: from structural and thermodynamics point of view.

To examine the molecular level understanding of temperature induced self-association of caffeine molecules in aqueous solution both in the presence and absence of salt NaCl, we have performed long MD simulations at a regime of temperatures ranging from 275 to 350 K with a temperature difference of 25 K. The calculations of different site-site radial distribution functions followed by coordination number analyses, the calculations of preferential interaction parameters, solvent accessible surface area, and cluster structure analyses show a depletion in the caffeine association propensity with increasing temperature. We have also observed the salting out effect of caffeine molecules in salt solution. The simultaneous presence of polar and nonpolar groups in a caffeine molecule leads to anisotropic hydration. Specifically, the hydration tendency of caffeine hydrophobic sites increases with increasing temperature, while hydrophilic sites tend to be less hydrated. This leads to a decrease in caffeine association. In accordance with some experimental studies on thermodynamics of caffeine association, we have also observed enthalpy driven association in pure water. But the presence of salt leads to entropy driven association specifically at higher temperature. This is due to the relatively stronger interactions of salt ions with caffeine at higher temperature.

[1]  Zukoski,et al.  Effect of Solvents on Interactions between Hydrophobic Self-Assembled Monolayers. , 1999, Journal of colloid and interface science.

[2]  Sandip Paul,et al.  Hydrophobic interactions in presence of osmolytes urea and trimethylamine-N-oxide. , 2011, The Journal of chemical physics.

[3]  Amalendu Chandra,et al.  Molecular dynamics simulations of aqueous NaCl and KCl solutions: Effects of ion concentration on the single-particle, pair, and collective dynamical properties of ions and water molecules , 2001 .

[4]  A. N. Veselkov,et al.  Hetero-association of caffeine and aromatic drugs and their competitive binding with a DNA oligomer , 2001, European Biophysics Journal.

[5]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

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

[7]  Paul E. Smith,et al.  Equilibrium dialysis data and the relationships between preferential interaction parameters for biological systems in terms of Kirkwood-Buff integrals. , 2006, The journal of physical chemistry. B.

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

[9]  S. Weerasinghe,et al.  Study of aggregate formation of caffeine in water by molecular dynamics simulation , 2011 .

[10]  S. Weerasinghe,et al.  Development of a molecular mechanics force field for caffeine to investigate the interactions of caffeine in different solvent media , 2010 .

[11]  S. L. Hansen,et al.  Thermodynamics of caffeine and the 1-1 caffeine-salicylate complex in water at 25.deg. , 1974 .

[12]  C. Vega,et al.  Radial distribution functions and densities for the SPC/E, TIP4P and TIP5P models for liquid water and ices Ih, Ic, II, III, IV, V, VI, VII, VIII, IX, XI and XII. , 2005, Physical chemistry chemical physics : PCCP.

[13]  N. Boden,et al.  Deuteron nuclear magnetic resonance spin echoes in solids , 1978 .

[14]  B. Berne,et al.  Temperature dependence of dimerization and dewetting of large-scale hydrophobes: a molecular dynamics study. , 2008, Journal of Physical Chemistry B.

[15]  D. Souza,et al.  Caffeine and an adenosine A2A receptor antagonist prevent memory impairment and synaptotoxicity in adult rats triggered by a convulsive episode in early life , 2010, Journal of neurochemistry.

[16]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[17]  T. Cheatham,et al.  Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations , 2008, The journal of physical chemistry. B.

[18]  José Mario Martínez,et al.  PACKMOL: A package for building initial configurations for molecular dynamics simulations , 2009, J. Comput. Chem..

[19]  M. Himmel,et al.  Caffeine and sugars interact in aqueous solutions: a simulation and NMR study. , 2012, The journal of physical chemistry. B.

[20]  A. Shestopalova Computer simulation of the association of caffeine and actinocin derivatives in aqueous solutions , 2006, Biofizika.

[21]  Sandip Paul,et al.  The effect of pressure on the hydration structure around hydrophobic solute: a molecular dynamics simulation study. , 2012, The Journal of chemical physics.

[22]  Alan K. Soper,et al.  The radial distribution functions of water and ice from 220 to 673 K and at pressures up to 400 MPa , 2000 .

[23]  Z. Tan,et al.  The standard molar enthalpy of formation, molar heat capacities, and thermal stability of anhydrous caffeine , 2007 .

[24]  B. Pettitt,et al.  Enthalpy-entropy contributions to the potential of mean force of nanoscopic hydrophobic solutes. , 2006, The journal of physical chemistry. B.

[25]  Paul E. Smith,et al.  A Kirkwood-Buff Derived Force Field for Aqueous Alkali Halides. , 2011, Journal of chemical theory and computation.

[26]  Weakly hydrated surfaces and the binding interactions of small biological solutes , 2012, European Biophysics Journal.

[27]  W. Pohle,et al.  On the interaction of caffeine with nucleic acids. IV. Studies of the caffeine-DNA interaction by infrared and ultraviolet linear dichroism, proton and deuteron nuclear magnetic resonance. , 1980, Biophysical chemistry.

[28]  M. Paulaitis,et al.  Effect of solute size and solute-water attractive interactions on hydration water structure around hydrophobic solutes. , 2001, Journal of the American Chemical Society.

[29]  J. Kirkwood,et al.  The Statistical Mechanical Theory of Solutions. I , 1951 .

[30]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

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

[32]  B. A. Patterson,et al.  Thermodynamics for self-association of caffeine in water: apparent molar volumes and apparent molar heat capacities of aqueous caffeine at temperatures from 278.15 to 393.15 K and at the pressure 0.35 MPa , 2002 .

[33]  P. Cremer,et al.  Interactions between macromolecules and ions: The Hofmeister series. , 2006, Current opinion in chemical biology.

[34]  Chandra,et al.  Effects of ion atmosphere on hydrogen-bond dynamics in aqueous electrolyte solutions , 2000, Physical review letters.

[35]  V. Crescenzi,et al.  Thermodynamics of caffeine aqueous solutions , 1976 .

[36]  A. Chandra Dynamical Behavior of Anion−Water and Water−Water Hydrogen Bonds in Aqueous Electrolyte Solutions: A Molecular Dynamics Study , 2003 .

[37]  Samantha Weerasinghe,et al.  Recent Applications of Kirkwood–Buff Theory to Biological Systems , 2007, Cell Biochemistry and Biophysics.

[38]  D. Chandler Interfaces and the driving force of hydrophobic assembly , 2005, Nature.

[39]  Chandler,et al.  Effect of environment on hydrogen bond dynamics in liquid water. , 1996, Physical review letters.

[40]  F. Stillinger Water Revisited , 1980, Science.

[41]  A. Cesàro,et al.  Thermodynamic properties of caffeine crystal forms , 1980 .

[42]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[43]  T. Head-Gordon,et al.  A New Solvent Model for Hydrophobic Association in Water I. Thermodynamics , 1995 .

[44]  Ricardo L. Mancera,et al.  Influence of Salt on Hydrophobic Effects: A Molecular Dynamics Study Using the Modified Hydration-Shell Hydrogen-Bond Model , 1999 .

[45]  P. Borer,et al.  1H‐ and 13C‐nmr studies on caffeine and its interaction with nucleic acids , 1980, Biopolymers.

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

[47]  David E. Smith,et al.  Free energy, entropy, and internal energy of hydrophobic interactions: Computer simulations , 1993 .

[48]  H. Chan,et al.  Configuration-dependent heat capacity of pairwise hydrophobic interactions. , 2001, Journal of the American Chemical Society.

[49]  Sandip Paul,et al.  Effects of dilute aqueous NaCl solution on caffeine aggregation. , 2013, The Journal of chemical physics.

[50]  A. Garcia,et al.  Enthalpy and entropy contributions to the pressure dependence of hydrophobic interactions , 2002 .