Next Article in Journal
Dose Titration of Solid Dosage Forms via FDM 3D-Printed Mini-Tablets
Previous Article in Journal
Compression Density as an Alternative to Identify an Optimal Moisture Content for High Shear Wet Granulation as an Initial Step for Spheronisation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Sodium Salicylate on Self-Aggregation and Caffeine Solubility in Water—A New Hypothesis from Experimental and Computational Data

1
Faculty of Sciences, University of Novi Sad, Trg Dositeja Obradovića 3, 21000 Novi Sad, Serbia
2
Faculty of Sport and Physical Education, University of Novi Sad, Lovćenska 16, 21000 Novi Sad, Serbia
3
Faculty of Chemistry and Technology, University of Split, Ruđera Boškovića 35, HR-21000 Split, Croatia
4
Department of Chemistry, Faculty of Science, University of Kragujevac, Radoja Domanovića 12, 34000 Kragujevac, Serbia
*
Author to whom correspondence should be addressed.
Pharmaceutics 2022, 14(11), 2304; https://doi.org/10.3390/pharmaceutics14112304
Submission received: 9 September 2022 / Revised: 14 October 2022 / Accepted: 21 October 2022 / Published: 26 October 2022
(This article belongs to the Section Pharmacokinetics and Pharmacodynamics)

Abstract

:
The present study analyzed experimental data from volumetric and viscosimetric measurements and computational simulations to understand caffeine hydration and aggregation properties in 0.1 mol∙kg−1 of sodium salicylate aqueous solution. Sodium salicylate reduces the bitter taste and increases the solubility of caffeine in water, which is the main reason for their combination in food products. The results noted in volumetric and viscosimetric measurements indicate that sodium salicylate promotes the self-aggregation of caffeine in water. After self-aggregation, the hydration number of caffeine significantly increases. Molecular simulations have allowed us to hypothesize how salicylate increases caffeine solubility. At the molecular level, relocating salicylate moiety from the parallel stacking (π–π) aromatic complex with caffeine and its hydration could be the main reason for increasing the solubility of caffeine in water. The presented study provides clear guidelines on the choice of additives to increase caffeine’s solubility in aqueous media. The choice of salicylate as an additive to increase the solubility of caffeine is very important because caffeine and salicylate are found in combination in a large number of formulations.

1. Introduction

Caffeine (1,3,7-trimethylxanthine, see Figure 1) belongs to the group of xanthine alkaloids together with theophylline and theobromine [1]. As an ingredient in coffee, black and green tea, and soft and energy drinks, caffeine is the most widely used psychoactive substance [2]. Pure caffeine is a solid substance, odorless and slightly bitter-tasting. The important physiological effects of caffeine are stimulation of the central nervous system, gastric acid secretion, and blood pressure increase in the short term of action [3]. The solubility of caffeine in water is relatively low (approximately 16 mg·mL−1 at room temperature), which is one of the crucial problems, especially in preparations consumed or stored at low temperatures [4]. The hydrophobicity of pharmacologically active drugs and active ingredients of supplements is a major problem in drug design, food science, delivery, drug effect, and bioavailability [5]. The possibility of transport and achieving optimal concentration in biological fluids is limited due to their low solubility in water.
It is considered that caffeine has limited solubility in water due to the self-association and aggregation of caffeine molecules by hydrophobic interactions [6,7,8,9]. A common solution for increasing solubility and limiting aggregation is to add some biocompatible molecules, excipients, or hydrotropes. Ahmad et al. investigated different salts’ (Na2SO4, NaCl, NaClO4, and NaSCN) effects on caffeine solubility and self-association, of which NaClO4 and NaSCN increase the solubility of caffeine in water [7]. The same effects were investigated by Shumilin et al. using sugars [5]. Mejri et al. used sucrose and β-cyclodextrin to determine hydration and self-association of caffeine molecules in an aqueous solution [10]. Cui explores interactions of a system consisting of caffeine, riboflavin, and water [11]. Vraneš et al. analyzed caffeine hydration and aggregation properties in aqueous solutions and compared these results to caffeine properties in the presence of adenosine triphosphate (ATP) molecules. They concluded that ATP negatively affects caffeine due to self-aggregation, leading to reduced bioavailability in the biological systems [12].
In earlier works, interactions between caffeine and benzoic acid derivatives (sodium benzoate, sodium salicylate, salicylic acid, acetylsalicylic acid, etc.) [13,14,15] were investigated, but how these benzoic acid derivatives affect the better solubility of caffeine has not been investigated so far. The mechanism by which these compounds increase the solubility of caffeine in water has remained unknown. The influence of sodium salicylate on the solubility of caffeine has not yet been experimentally determined. Sodium salicylate (Figure 2) (SS) is a non-steroidal anti-inflammatory drug, and it belongs to the derivates of salicylic acid and is used as an analgesic, antipyretic, antioxidant, and antimicrobial drug [16]. In spite of that, salicylate is most frequently used in the food industry as a dietary supplement and preservative. It can also be found in food as a naturally occurring compound [17,18,19].
Ankita et al. analyzed physicochemical properties of sodium salicylate in aqueous solution and aqueous glucose/sucrose solutions at temperatures from 293.15 to 318.15 K. They concluded that sodium salicylate is a typical structure-maker (positive values of B-coefficient which decrease with increasing temperature) [20]. Drugs containing a combination of caffeine and acetylsalicylic acid have been available for many years and are widely used clinically [21]. It has been proven that the addition of caffeine to drugs containing acetylsalicylic acid leads to increased activity of the drug; for example, the onset of analgesia is accelerated, and pain relief is more effective compared to drugs that do not contain caffeine in their composition [22]. Caffeine in combination with SS can be found in various drugs or preparations, such as Algopirin (used for acute low back pain), Acetaminophen (used for temporary relief of pain), Excedrin (used for Migraine headache), Ephedrine-Caffeine-Acetylsalicylate stack (drug combination used in weight loss), etc. Today, there is an increasing number of scientific papers that research the influence of sodium salicylate molecules on caffeine, both in drug formulations and aqueous solutions. For example, one of the recent works by Kim et al. talks about the influence of caffeine on the absorption rate of aspirin [23]. Therefore, our goal in this paper is to discuss how adding a typical structure-maker in 0.1 mol·kg−1 of SS affects caffeine hydration and aggregation properties through hydration parameters obtained from volumetric and viscosimetric measurements, solubility determination, and computational simulations.

2. Materials and Methods

2.1. Materials

Caffeine (anhydrous, CAS number: 58-08-2, mass fraction purity ≥ 99%), chloroform (CAS number: 67-66-3, mass fraction purity ≥ 99%), and sodium salicylate (CAS number: 54-21-7, mass fraction purity ≥ 99,5%), presented in Table S1, were procured from Sigma Aldrich (St. Louis, MI, USA), and were used as-received without further purification.

2.2. Apparatus and Procedure

The 0.1 mol∙kg−1 SS aqueous solution was used as a solvent to prepare a set of caffeine solutions for density and viscosity measurements up to caffeine molality of m = 0.10584 mol·kg−1. Experimental measurements of solution densities and viscosities were performed at five temperatures: T = 293.15, 298.15, 303.15, 308.15, and 313.15 K. A detailed description of the apparatuses used is provided in the Supplementary Materials.
The determination of caffeine solubility in pure water was performed at five temperatures: T = 293.15, 298.15, 303.15, 308.15, and 313.15 K. Solubility of caffeine in pure water was determined using the gravimetric method already described in the work of Romero and Oviedo [24]. A double-layer glass flask with a volume of 60 cm3 was used for solubility determination. It was held in a constant temperature bath for the necessary time to reach equilibrium at each temperature. The temperature was controlled with a Lauda E-100 circulator. The flask was placed on an automatic mixer with a magnetic core. Caffeine saturated solutions were prepared by adding an appropriate amount of caffeine to 40 cm3 of ultrapure water and stirring the flasks containing the samples with a temperature control of ±0.01 K. Before transferring the samples, the syringe with a needle was immersed in the flask, and the stirrer was turned off for 2 h to allow the undissolved solute to settle.
The samples of 5 cm3 were withdrawn using a syringe. The needle was quickly removed, and a filter (33 mm diameter sterile Millex-HA syringe filter with a 0.45 µm pore size mixed cellulose esters membrane) was placed on the syringe to prevent undissolved caffeine from causing inaccurate results. The samples were transferred to 5 cm3 pre-weighed glass flasks and their mass was measured with an accuracy of ±1 × 10−5 g in the lower range. The samples were placed in a vacuum oven at 383.15 K and evaporated to dryness to recover the caffeine. The mass of the caffeine was determined gravimetrically. Solubility measurements were performed at temperatures between 293.15 and 313.15 K. Each obtained value represents the average of at least five independent measurements. The uncertainty in the mass fraction solubility is ±0.001. Reproducibility was found to be better than ±0.002.
The determination of caffeine solubility in 0.1 mol∙kg−1 SS aqueous solutions was performed at the same temperature range, glass flask, and conditions as in the previously described method. However, in the presence of sodium salicylate, it is necessary to separate the caffeine. The separation of caffeine and sodium salicylate molecules was performed by liquid–liquid extraction with chloroform. Caffeine dissolves well in chloroform, while the solubility of sodium salicylate is negligible. Due to that, the procedure has been slightly modified. After taking a 5 cm3 sample from the double-layer glass flask, the sample was transferred to a glass separating funnel, and 10 cm3 of chloroform was added. Extraction with chloroform was repeated three times for each sample. After each extraction, the bottom chloroform layer was dropped into a previously measured balloon. Chloroform was removed using a rotary evaporator. The sample was further dried for 2 h at 383.15 K until a constant mass was reached. The extracted mass of caffeine corresponds to the mass of caffeine that can be dissolved in a sampled volume of 0.1 mol·kg−1 of SS aqueous solution at a given temperature.
The interactions of caffeine in pure water and in SS aqueous solution were theoretically investigated with molecular dynamics (MD) simulations using the Yasara Structure (version 16.3.8.) [25] packages by opting for the use of AMBER14 as a force field. For AMBER14, the TIP3P water model was used. These MD simulations were accomplished in NPT canonical ensemble at a constant pressure of 1 × 105 Pa and a temperature of 298.15 K using periodic boundary conditions. Pressure and temperature control was performed using a Nose–Hoover barostat, while the “cut-off” radius used for all Lenard–Jones interactions was 10 Å. To imitate physiological conditions, counter-ions were added to neutralize the system. Na or Cl was added in replacement of water to give a total NaCl concentration of 0.9%. pH was maintained at 7.4. The particle-mash Ewald (PME) summation was applied to correct the effect of the cut-off radius with respect to Coulomb interactions. The initial phase of the MD simulation involved establishing equilibrium, which lasted an average of about 10 ns, and these data were not analyzed later. The system after equilibration was analyzed by examining the trajectory from which the radial distribution functions (RDF) were examined. All simulations lasted a minimum of 60 ns.

3. Results and Discussion

3.1. Solubility Results

Experimental data for the solubility of caffeine in 0.1 mol∙kg1 SS aqueous solutions and the solubility of caffeine in pure water were obtained experimentally and presented in Table 1. As can be seen, the solubility of caffeine increased with temperature in both cases. In 0.1 mol∙kg1 SS aqueous solutions, caffeine had a 2–2.7 times higher solubility than caffeine in pure water, depending on the temperature.
Based on the modified Van ’t Hoff equation, the apparent standard dissolution enthalpy (ΔsolH°), entropy (ΔsolS°), and Gibbs energy (ΔsolG°) change of solutes in water and 0.1 mol·kg1 SS aqueous solutions were calculated from the results of experimental solubility data [26,27]:
ln x 1 1 / T 1 / T hm p = Δ sol H o R .
where x1 is the mole fraction of caffeine in solutions, T is the absolute temperature (K), and R is the universal gas constant with a value of 8.314 J·K1·mol1. Thm represents the harmonic temperature, defined as:
T hm = n / j 1 n 1 T j ,
where n equals the number of temperature points, Tj is the experimental temperature, and the obtained value is 302.99 K. ΔsolH° is defined as the apparent standard mole dissolution enthalpy change of solute dissolved in water and 0.1 mol∙kg1 SS aqueous solution, obtained from the slope of the fitted line curves, lnx1 versus (1/T − 1/Thm) (Figure 3).
From the intercept obtained from the modified Van ’t Hoff plot, the ΔsolG° and ΔsolS° can be calculated by the following equations:
Δ sol G o = R T hm I n t e r c e p t ,
Δ sol S o = Δ sol H o Δ sol G o T hm .
Determination of the relative contribution of enthalpy (ζH) and the relative contribution of entropy (ζTS) were calculated using the following expressions [26,27]:
ς H = Δ sol H o Δ sol H o + T hm Δ sol S o ,
ς T S = T hm Δ sol S o Δ sol H o + T hm Δ sol S o .
The obtained apparent thermodynamic properties’ values of ΔsolG°, ΔsolH°, and Δsol, as well as of ζH and ζTS for caffeine in water and caffeine in SS aqueous solution, are presented in Table 2.
Table 2 shows that ΔsolG° and ΔsolH° values were positive for both examined solutions, caffeine in aqueous solution and caffeine in SS aqueous solution. The dissolution processes are endothermic, which indicates that dissolving caffeine in water and aqueous sodium salicylate solution proceeds with absorption of heat. A positive value of the enthalpy of dissolution indicates that the heat energy released due to the hydration of compounds is less than the energy required to break down the interactions between caffeine molecules. Based on these results, we can assume that in the presence of sodium salicylate, caffeine self-aggregates are formed in an aqueous solution, similar to that in pure water. Lower enthalpy values in sodium salicylate solution may be due to the hydration of sodium and salicylate ions. The negative ΔsolS° value for both examined solutions denotes that dissolution of caffeine in water proceeds with an increase in the order of the system. The formation of caffeine self-aggregates is undoubtedly one of the explanations for the negative values of dissolution entropy.
On the other hand, the influence of caffeine molecules on the structural organization of water molecules, i.e., its structure-making/breaking properties, will be examined based on volumetric and viscometric measurements. The calculated values of the relative contribution of enthalpy and entropy imply that the main contributor to ΔsolG°, in caffeine dissolution in aqueous solutions and SS aqueous solutions, is enthalpy. Moreover, the endothermic enthalpy of the dissolution process is in accordance with the obtained solubility increase with rising temperature.

3.2. Volumetric Properties

For the determination of caffeine volumetric properties in 0.1 mol·kg−1 SS aqueous solution, data obtained from experimental measurements of solution densities were collected and provided in Table S2. Densities of solutions increased with a molality of caffeine at all temperatures (Figure S1).
The obtained solution densities were used to calculate the apparent molar volume ( V ϕ ) of caffeine, from the equation presented in the Supplementary Materials. The calculated values are shown in Table S2. As shown in Figure 4, the caffeine values in 0.1 mol∙kg−1 SS aqueous solutions increased with caffeine molality at all temperatures.
The V ϕ values were fitted with the Masson’s equation modified for non-electrolytes [28]:
V ϕ = V ϕ o + S v m ,
and the plots of V ϕ versus m presented in Figure 3. The V ϕ o value represents the apparent molar volume of solute at infinite dilution, and Sv is the experimental slope, which provides information about solute–solute interactions. The obtained values of V ϕ o and Sv are shown in Table 3 with the standard deviation (σ) and regression coefficient (R2).
The positive Sv values indicate the existence of strong caffeine–caffeine interactions. These values are significantly higher than pure water and 0.1 mol∙kg−1 ATP aqueous solution (see Table S3). As caffeine self-aggregation occurs in pure water and 0.1 mol∙kg−1 ATP water solution [12], higher values of the Sv coefficient indicate that the process of self-aggregation in the presence of SS is even more pronounced. At first glance, this is unusual. The self-aggregation of molecules usually reduces their solubility in water, while our results showed an increase in solubility of caffeine in the presence of sodium salicylate.
Furthermore, the obtained V ϕ o values of caffeine in SS aqueous solutions were positive and increased with rising temperature (Figure S2), but they were significantly lower than in pure water or in the presence of ATP molecules (see Table S3). This means that adding caffeine to a SS aqueous solution will lead to a significantly smaller increase in volume than in pure water or in the presence of ATP. Lower values of apparent molar volumes may also be a consequence of enhanced self-aggregation of caffeine molecules in the presence of sodium salicylate and/or strong interactions between caffeine molecules and salicylate anion. It is not possible to draw a conclusion based on volumetric measurements alone. That is why we need viscosity measurements and computer simulations.
The values of V ϕ o were fitted as a function of the temperature using the equation of the second order:
V ϕ o = a 0 + a 1 T + a 2 T 2 .
Calculated coefficients a0, a1, and a2 of caffeine in 0.1 mol∙kg−1 SS aqueous solutions together with regression coefficients (R2) are presented in Table S4. Based on the calculated coefficients, the limiting apparent molar expansibility, E ϕ o , of caffeine in SS aqueous solutions was calculated and presented in Table 4, together with values in pure water:
E ϕ o = V ϕ o T p = a 1 + 2 a 2 T .
The values of E ϕ o were positive in the whole temperature range and significantly higher at lower temperatures than in pure water (Table 4 and Figure S3). Higher E ϕ o values indicate a faster release of water molecules from the hydration sphere. However, E ϕ o values decreased much faster with increasing temperature in the presence of SS compared to pure water. The Hepler’s coefficient defines the rate of this change and can be calculated with the following equation [29]:
E ϕ o T p = 2 V ϕ o T 2 p = 2 a 2 .
The obtained negative values for Hepler’s coefficient, presented in Table 4, suggest caffeine structure-breaking properties in sodium salicylate aqueous solutions. However, the efficiency of Hepler’s coefficient to estimate the structure-making/breaking properties in water systems where significant self-aggregation of molecules occurs is quite questionable.
From the available data for V ϕ o of caffeine in water and SS aqueous solutions, the limiting apparent transfer molar volumes ( Δ t r V ϕ o ) can be calculated with the following equation [12]:
Δ tr V ϕ o = V ϕ o ( caffeine + 0.1 mol · kg 1 aqueous   SS   solution ) V ϕ o ( caffeine + water ) .
Values of V ϕ o for the aqueous solution of caffeine are presented in Table S3. The obtained values are presented in Table 5. The negative values of Δ t r V ϕ o indicate a more pronounced self-aggregation of caffeine in the presence of SS or a strong interaction between caffeine and SS, which was discussed earlier in this paper.
According to Birch et al. [30], an organization of water molecules around the solute regulates transport in the taste epithelium, binding to the receptor and finally inducing a taste response. In this study, to define the taste quality of studied solutes, the apparent specific volume at infinite dilution, V ϕ o (results in Table 5), was calculated by dividing caffeine V ϕ o values with caffeine molar mass (M = 194.19 g·mol−1). It was found that the apparent specific volume at infinite dilution is a parameter of taste quality in the following order: salty (0.1–0.3 cm3·g−1) < sour (0.3–0.5 cm3·g−1) < sweet (0.5–0.7 cm3·g−1) < bitter (0.7–0.9 cm3·g−1). From Table 5, it can be seen that the calculated values of V ϕ o 0.69 cm3·g−1 are on the border between the sweet and bitter taste of caffeine in SS aqueous solutions at all examined temperatures. On the other hand, Vraneš et al. calculated values for caffeine in pure water, 0.73–0.75 cm3·g−1, pointing out the bitter taste of caffeine at all examined temperatures (283.15–313.15 K) [12]. Thus, the presence of SS in water reduces the bitter taste of caffeine. These data are essential due to the joint presence of caffeine and benzoic acid derivatives in food products.

3.3. Viscosimetric Properties

The viscosity measurements of caffeine in 0.1 mol∙kg−1 SS aqueous solutions were measured in the temperature range from T = 293.15 to 313.15 K up to molality of caffeine m = 0.1058 mol·kg−1 at the atmospheric pressure (p = 1 × 105 Pa). Data obtained from experimental caffeine measurements in 0.1 mol∙kg−1 SS aqueous solution viscosities are presented in Table S5, together with the values for aqueous solutions of caffeine [12]. Table S5 shows viscosity values for 0.1 mol∙kg−1 SS in water. Viscosity values of 0.1 mol∙kg−1 SS in water were higher than viscosities of pure water.
The viscosity values for caffeine in 0.1 mol∙kg−1 SS aqueous solution were higher than the values in pure water. In both systems, viscosity decreased with the temperature increase. One of the most common criteria for the determination of water organization around solute is the viscosity B-coefficient obtained from the Jones–Dole’s equation [31]:
η η o 1 = B c ,
where ηo is the viscosity of the pure solvent, η is the viscosity of the caffeine in SS aqueous solutions, and c is concentration. Following the Jones–Dole’s equation, the reduced viscosity (η/ηo) dependence on the concentration was linear (Figure 5). Calculated B-coefficients are presented in Table 6.
The relation between the B-coefficient from the Jones–Dole’s equation and the change in the structural order of water molecules after adding a solute is well-known. If values of the B-coefficient are positive, that indicates the structure-making properties of the solute, pointing out an increase in the local order of the water molecules. Negative B-coefficient values indicate structure-breaking properties, weakening hydrogen bonds between water molecules near the solute [32]. From the calculated results given in Table 6, B-coefficients for caffeine in SS aqueous solutions had positive values and decreased with the rising temperature. Thus, both conditions (B > 0 and dB/dT < 0) were satisfied, and the caffeine molecule in 0.1 mol∙kg−1 sodium salicylate water solution can be declared a structure-maker. Table 6 also includes the B-coefficients of caffeine in aqueous solutions tested by Vraneš et al. [12]. From Table 6, it can be seen that B-coefficients in pure water were positive and increased with the rising temperature, suggesting that caffeine is an atypical structure-maker. The values of the B-coefficient in the presence of sodium salicylate were higher compared to pure water, especially at lower temperatures. At higher temperatures (313.15 K), these values were similar in both systems. The obtained trends indicate that the influence of caffeine on the local arrangement of water molecules at higher temperatures is similar and that the presence of SS makes a significant difference only at lower temperatures. After analyzing the results of volumetric and viscometric measurements, we can conclude that the presence of SS promotes the self-aggregation of caffeine in water. The formation of self-aggregates contributes to the structural arrangement of water. However, it remains unclear what interactions occur between caffeine molecules and sodium salicylate and whether salicylate anion is incorporated into aggregates.

3.4. Thermodynamic Activation Parameters for Viscous Flow

Feakins et al. [33] suggested calculation of the transition-state treatment of relative viscosity using viscosity obtained data:
Δ μ 1 o = R T ln η o V ¯ 1 o h N A
where Δ μ 1 o   is free energy of activation of viscous flow per mole of the pure water, N A is Avogadro’s number, V ¯ 1 o is partial molar volumes of water at infinite dilution, h is the Planck constant, ηo is the viscosity of the solvent, R is the universal gas constant, and T is the absolute temperature. The following equation calculates the free energy of activation of viscous flow per mole of solute:
Δ μ 2 o = Δ μ 1 o + R T V ¯ 1 o B + V ¯ 2 o V ¯ 1 o ,
where B is the viscosity coefficient, and V ¯ 2 o is the partial molar volume at infinite dilution of solute. The calculated values of V ¯ 1 o , V ¯ 2 o , Δ μ 1 o , and Δ μ 2 o are presented in Table 7.
The values of Δ μ 2 o were positive and larger than Δ μ 1 o at all investigated temperatures, indicating that the interactions between caffeine and solvent molecules (water and SS) are more pronounced in the ground state than in the transition state. The transition state formation occurs with the breaking and deformation of the intermolecular bonds between solvent molecules [34,35]. The positive values of the term ( Δ μ 2 o Δ μ 1 o   ) were described by Glasstone et al. as “loss of the structure-breaking contribution special to liquid water” [36]. In other words, the positive values of a mentioned term are typical for structure-making compounds. These data agree well with the analyzed results of the Jones-Dole’s equation’s B-coefficient.
The equation for determination of the entropy of activation for the viscous flow ( Δ S 2 o ) of the caffeine in 0.1 mol∙kg−1 SS aqueous solutions is:
Δ S 2 o = Δ μ 2 o T ,
while the activation enthalpy ( Δ H 2 o ) can be calculated with the equation:
Δ H 2 o = Δ μ 2 o + T Δ S 2 o .
The values of Δ S 2 o and Δ H 2 o were 239.66 J·mol−1·K−1 and 162.33 kJ·mol−1 for caffeine in SS aqueous solutions. These results support that formation of the transition state of caffeine in SS aqueous solutions is associated with bond-breaking and a decrease in order [34,35].

3.5. Hydration Number

The results of volumetric and viscosimetric measurements and higher solubility of caffeine indicate more significant caffeine interactions with water molecules in the presence of SS. Therefore, it is necessary to calculate the hydration number of caffeine in 0.1 mol∙kg−1 sodium salicylate water solution. The hydration number of caffeine is calculated from the volumetric properties using the following equation [37]:
h n = V ϕ o e V e o V b o ,
where V ϕ o e is the electrostriction partial molar volume, V e o is the molar volume of electrostricted water, and V b o is the molar volume of bulk water. The literature values of ( V e o V b o ) are presented in Table S6. The mathematical approach for calculating the hydration number from the volumetric and viscosimetric parameters is described in the Supplementary Materials. The hydration numbers calculated in both described ways for caffeine in 0.1 mol∙kg−1 SS aqueous solution are presented in Table 8, along with the data for pure water [12].
From Table 8, it can be seen that the hydration numbers of caffeine in 0.1 mol∙kg−1 SS aqueous solutions are in good correlation. The values of the hydration number decreased with temperature, which is in accordance with the results obtained from expansibility. In our previous study, we calculated the hydration numbers of caffeine in pure water (Table 8), and the values obtained were significantly lower than those obtained in the presence of SS. However, these values were obtained at significantly lower concentrations of caffeine (about 0.06 mol∙kg−1), at which self-aggregation does not occur. To better understand the hydrating properties of caffeine at higher concentrations and the role of SS in increasing caffeine solubility, we will look for the answer in computer simulations.

3.6. Computational Study

Computational studies were used to confirm the results obtained from volumetric, viscometric, and solubility measurements. MD simulations were performed at caffeine molality m = 0.25 mol∙kg−1 for both tested cases at a temperature 298.15 K (250 molecules of caffeine, 100 molecules of SS (0.1 mol∙kg−1), and 55,555 molecules of water, while in the case of caffeine in pure water, 250 molecules of caffeine and 55,555 molecules of water were used).
It is a well-known scientific fact that caffeine molecules show a tendency towards self-aggregation in water solutions at higher concentrations [5,38]. As shown in Figure 6a, caffeine molecules tend to interact with each other, forming self-aggregates at molality m = 0.25 mol∙kg−1. Experimental results indicate that the presence of sodium salicylates increases the self-aggregation of caffeine molecules in water. Therefore, the radial distributions of caffeine molecules around the caffeine in pure water and the presence of SS were calculated. Figure 7 shows the obtained results from RDFs of the caffeine–caffeine center of mass. The first and second peaks of this distribution function appeared at 3.65 and 7.1 Å, indicating the existence of caffeine’s first and second coordination shells around caffeine.
The intensity of both peaks of radial function increased after the addition of SS to water, indicating increasing probabilities of self-aggregation between caffeine molecules. The large number of caffeine molecules in the coordination shells of caffeine confirmed that the presence of SS increased the self-aggregation of caffeine in water. Caffeine molecules interact through π–π interactions, as shown in Figure 6b. In Figure 6b, it can also be seen that salicylate molecules are not involved in caffeine self-aggregates. To confirm the absence of SS molecules from caffeine aggregates, RDF caffeine-salicylate was calculated, considering the center of masses of caffeine and salicylate. The calculated values were g(r) = 0 in the whole distance range from 0 to 10 Å, which confirms the absence of any interactions between these particles in equilibrium and at selected concentrations (Figure 8).
To further analyze the hydration of caffeine molecules in the presence and absence of SS, calculations of the radial distribution function (RDFs) were applied. Figure 9 shows the RDFs of water oxygen atoms (Ow) around selected atomic sites of caffeine: O1, O2, N4, and H1. Hydrophilic centers in the caffeine molecule have been selected because they can form the strongest interactions with water molecules. As hydrogen bond acceptors, oxygen atoms from carbonyl groups, O1 and O2, as well as N4 atom from the imidazole ring were selected. The most acidic hydrogen atom (H1) was also considered the potential proton-donor for forming the H-bond with water molecules. Figure 9a–d show the values of g(r) in pure water and in the presence of sodium salicylate for all selected caffeine atoms. The obtained results indicate that sodium salicylate reduces the number of water molecules in the hydration spheres of O1 and H1 atoms. In the hydration sphere of the imidazole nitrogen atom (N4), the number of water molecules increased in the presence of SS, while in the case of carbonyl oxygen (O2), it did not change.
According to the procedure published elsewhere [12], the hydration number of caffeine in the presence and absence of SS was calculated to quantify the obtained results. The calculated value of the hydration number of caffeine in pure water at 298.15 K and caffeine molality m = 0.25 mol∙kg−1 was hn = 5.25. In 0.1 mol∙kg−1 sodium salicylate aqueous solutions, at the same molality, the hydration number of caffeine had a lower value, hn = 4.24. This value of hydration number in the presence of SS is in good agreement with the experimentally obtained results (Table 8). In our previous work [12], the hydration number of caffeine at a molality of 0.06 mol∙kg−1 (unfortunately, these data are missing in that paper) was calculated using an identical computer simulations procedure, and a value of hn = 0.82 at T = 298.15 K was obtained. At a molality of 0.06 mol∙kg−1, self-aggregation of caffeine in water was not observed. Thus, computer simulations show that the caffeine hydration number increased significantly during self-aggregation (from 0.82 to 5.25). The main reason for the increase in hydration number is that caffeine molecules form aggregates primarily through π–π interactions (see Figure 6b). As a result, caffeine molecules’ hydrophobic surfaces in self-aggregate are reduced [38], while hydrophilic centers have more probability for interaction with water molecules.
The presence of SS reduces the hydration number of caffeine (from 5.25 to 4.24), which is in concordance with the results of other studies. Sodium salicylate and many other polar molecules, such as sodium chloride [38], ATP [12], and sucrose [5], have a dehydrating effect on caffeine molecules. The mentioned molecules have more pronounced hydration abilities than caffeine, and by binding water molecules to themselves, they reduce the number of water molecules in the hydration shells of caffeine.
However, the main question remains open: how does sodium salicylate increase the solubility of caffeine in water? If we compare the results of the molecular dynamics simulations for the system NaCl + caffeine with our results, we can conclude that they are very similar [38]. In the presence of NaCl, caffeine self-aggregation in water increases, and RDFs peaks occur at the exact distances from the caffeine–caffeine center of mass. The addition of NaCl reduces the hydration number of water molecules around the hydrophilic centers of caffeine. These results were also obtained in this paper. Furthermore, both NaCl and SS have structure-making properties to water molecules.
The work of Ahmad et al. [7] showed that the presence of a cosmotrope (structure-makers, such as NaCl and Na2SO4) increases self-aggregation and decreases caffeine solubility, which is a consequence of the salting-out effect of these salts. The salting-out effect on caffeine molecules is also shown by some neutral molecules, such as sucrose [5]. Sucrose has structure-making properties, increasing the self-aggregation of caffeine and reducing its solubility.
However, SS has structure-maker properties but increases the solubility of caffeine.
The cosmotropic behavior of some additives does not appear to be a “driving force” to promote caffeine self-aggregation in water. Combining the Kirkwood–Buff solution theory with the isodesmic model of caffeine association, Shimizu [39] concluded that changes in the water structure upon the addition of cosmotropic or chaotropic additives have a negligible effect on the self-aggregation of caffeine. In the same paper, Shimizu concludes that the real driving force for caffeine self-aggregation is the interaction between additives and caffeine. Self-aggregation increases if additives are removed from caffeine molecules and decreases if additives stick around caffeine. In the mentioned manuscript, the author emphasizes that a new theory has important consequences for predicting the influence of additives on the solubility of macromolecules, but he did not elaborate on that influence.
Therefore, we need to analyze what interactions SS can form with the caffeine molecule, increasing its solubility, unlike sodium chloride and sucrose. The logical answer is π–π interactions because both salicylate and caffeine are aromatic compounds. However, we have already discussed (Figure 6b) that salicylate anions are not included in caffeine self-aggregates. Therefore, we set up a new MD simulation with a significantly lower concentration of caffeine (0.06 mol∙kg−1). At this concentration, self-association in pure water does not occur. The concentration of SS was 0.1 mol∙kg−1 and the number of water molecules was 5555. Figure 10 shows a snapshot of the distribution of caffeine and SS molecules after equilibrium (30 ns). It can be seen from the figure that the presence of SS initiated the self-aggregation of caffeine, but that not all molecules are included in the aggregates. However, Figure 10 also shows π–π interactions between caffeine and salicylate molecules. We calculated the radial distribution function concerning the center of masses to better understand the spatial relationship of salicylate and caffeine molecules (see Figure 8). The peak maximum appeared at a distance of 3.86 Å, which is very close to the value of the first peak of the radial distribution of the caffeine–caffeine center of mass (3.65 Å). This means that a molecule of caffeine displaces salicylate during aggregation and takes its place next to another molecule of caffeine. Finally, we can assume the mechanism of salicylate increasing caffeine’s solubility.
After adding caffeine to the SS solution (Figure 11a), a parallel stacking (π–π) aromatic complex between aromatic rings was formed (Figure 11b). Caffeine-salicylate complexes provide better solubility of caffeine monomers and allow local caffeine supersaturation formation (Figure 10c). After contact with two or more complexes, the salicylate anion was released, forming a π–π interaction between caffeine molecules (Figure 11d). The driving force for this process is the good solvation of the salicylate anion by water molecules and the tendency of caffeine molecules to aggregate with each other. We calculated the hydration number for salicylate anions not bonded with caffeine (see Figure 6b), and the average value was 5.47. Hydration of the salicylate ion is an exothermic process. Additional heat release due to hydration of the salicylate anion is a reason why the enthalpy of dissolution is lower in the presence of SS than in pure water. Releasing salicylate molecules from the complex with caffeine and its hydration after that is probably the main reason for increasing the solubility of caffeine in water.
Therefore, for an additive to increase caffeine’s solubility in water through forming caffeine self-aggregates, it needs to form π–π interactions with caffeine and have pronounced hydrating properties. This hypothesis needs to be tested in the presence of several additives. It is necessary to choose additives that can and cannot form π–π interactions with caffeine and additives with different hydrating properties, i.e., structure-makers and structure-breakers.

4. Conclusions

The objective of this study was to investigate the solubility of caffeine in the presence of SS both experimentally and computationally. The solubility of caffeine was 2–2.7 times higher in 0.1 mol∙kg−1 SS aqueous solutions than in pure water, and values increased with the rising temperature. Obtained results from volumetric and viscosimetric measurements indicated that SS increased caffeine self-aggregation in water and significantly reduced the bitterness of caffeine molecule in aqueous solution. The hydration number of caffeine increased after self-aggregation. The increase in caffeine solubility in the presence of salicylates is not a consequence of the structure-making properties of the salicylate anion but its interactions with caffeine molecules. The mechanism by which salicylate affects the solubility of caffeine in an aqueous solution has been hypothesized. The better solubility of caffeine monomers and the ability to achieve local supersaturation of solution are probably consequences of forming complexes with the salicylate anion through π–π interactions. Release and hydration of the salicylate anion after self-aggregation of caffeine is the main reason for the lower enthalpy dissolution of caffeine in the presence of SS. It remains an open question whether self-aggregation of caffeine always leads to a decrease in solubility. The authors know that it is not easy to answer these questions, so the plan is to start a series of experiments, for example, self-aggregation of caffeine in the presence of the ionic liquid based on benzoic acid derivates. One of them is the well-known ionic liquid 1-butyl-3-methylimidazolium salicylate, where we can assume that caffeine molecules can interact through pi–pi interactions with a cation and anion, or just one of them. The research conducted in this work is of great importance for future formulations in pharmacy, which will contain a combination of caffeine and sodium salicylate or only one of them.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pharmaceutics14112304/s1, Figure S1: Variation of density (d) for caffeine in SS aqueous solutions as a function of molality, at different temperatures, (T) = (■) 293.15, () 298.15, () 303.15, () 308.15, and () 313.15 K; Figure S2: Variation of the caffeine apparent molar volume at infinite dilution ( where V ϕ o ) with temperature in (■) SS aqueous solutions; Figure S3: Variation of the limiting apparent molar expansibility ( where E ϕ o ) of caffeine () in SS aqueous solutions with temperature; Table S1: Provenance, purity, and structure of the used chemicals; Table S2: Molality (m), density (d), and apparent molar volume (Vφ) of caffeine in 0.1 mol∙kg−1 SS aqueous solutions in the temperature range from T = 283.15 to 313.15 K at p = 1·105 Pa; Table S3: Fitting parameter obtained using the Masson′s equation modified for non-electrolytes of caffeine in aqueous solutions and in 0.1 ATP aqueous solutions in the temperature range from T = 293.15 to 313.15 K; Table S4: Fitting coefficients a0, a1, and a2 of caffeine in 0.1 mol∙kg−1 sodium salicylate aqueous solutions, with the regression coefficient (R2); Table S5: Viscosity of caffeine in aqueous solution [1] and in 0.1 mol∙kg−1 SS aqueous solutions in the temperature range from T = 283.15 to 313.15 K at p = 1·105 Pa; Table S6: The literature values of ( V e o   V b o   ) in the temperature range from T = 293.15 to 313.15 K. References [12,40,41,42] are cited in the Supplementary Materials.

Author Contributions

M.V., conceptualization, supervision, writing—original draft; T.T.B., investigation, writing—original draft; P.D., funding acquisition, writing– review and editing; T.T., funding acquisition, writing—review and editing; N.J., supervision, writing—review and editing; R.T., supervision, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Provincial Secretariat for Higher Education and Scientific Research (142-451-2594).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge financial support of the Ministry of Education, Science and Technological Development of the Republic of Serbia (Grant No. 451-03-68/2022-14/200125), and the Secretariat for Higher Education and Scientific Research of the Autonomous Province of Vojvodina, Serbia (Grant No. 142-451-2545/2021-01/2). The authors are grateful to Aleksandar Tot (KTH Royal Institute of Technology: Stockholm) for help during the interpretation of the results from the computational study (molecular dynamic simulations).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Evattand, D.P.; Griffiths, R.R. Caffeine. In Encyclopedia of Human Nutrition; Elsevier Inc.: Amsterdam, The Netherlands, 2013; pp. 221–227. [Google Scholar]
  2. Ioannidis, K.; Chamberlain, S.R.; Müller, U. Ostracising caffeine from the pharmacological arsenal for attention-deficit hyperactivity disorder—Was this a correct decision? A literature review. J. Psychopharmacol. 2014, 28, 830–836. [Google Scholar] [CrossRef] [PubMed]
  3. Ashihara, H.; Crozier, A. Caffeine: A well known but little mentioned compound in plant science. Trends Plant Sci. 2001, 6, 407–413. [Google Scholar] [CrossRef]
  4. Oestreich, S.J. Caffeine: Characterization and properties. In Encyclopedia of Food and Health; Elsevier Inc.: Amsterdam, The Netherlands, 2016; pp. 556–572. [Google Scholar]
  5. Shumilin, I.; Allolio, C.; Harries, D. How Sugars Modify Caffeine Self-Association and Solubility: Resolving a Mechanism of Selective Hydrotropy. J. Am. Chem. Soc. 2019, 141, 18056–18063. [Google Scholar] [CrossRef] [PubMed]
  6. Cesaro, A.; Russo, E.; Crescenzi, V. Thermodynamics of caffeine aqueous solutions. J. Phys. Chem. 1976, 80, 335–339. [Google Scholar] [CrossRef]
  7. Al-Maaieh, A.; Flanagan, D.R. Salt effects on caffeine solubility, distribution, and self-association. J. Pharm. Sci. 2002, 91, 1000–1008. [Google Scholar] [CrossRef]
  8. Lilley, T.H.; Linsdell, H.; Maestre, A. Association of caffeine in water and in aqueous solutions of sucrose. J. Chem. Soc. Faraday Trans. 1992, 88, 2865–2870. [Google Scholar] [CrossRef]
  9. Guttman, D.; Higuchi, T. Reversible association of caffeine and of some caffeine homologs in aqueous solution. J. Am. Pharm. Assoc. 1957, 46, 4–10. [Google Scholar] [CrossRef]
  10. Mejri, M.; BenSouissi, A.; Aroulmoji, V.; Rogé, B. Hydration and self-association of caffeine molecules in aqueous solution: Comparative effects of sucrose and β-cyclodextrin. Spectrochim. Acta A Mol. Biomol. Spectrosc. 2009, 73, 6–10. [Google Scholar] [CrossRef]
  11. Cui, Y. Parallel stacking of caffeine with riboflavin in aqueous solutions: The potential mechanism for hydrotropic solubilization of riboflavin. Int. J. Pharm. 2010, 397, 36–43. [Google Scholar] [CrossRef]
  12. Vraneš, M.; Panić, J.; Tot, A.; Gadžurić, S.; Podlipnik, Č.; Bešter-Rogač, M. How the presence of ATP affect caffeine hydration and self-aggregation? J. Mol. Liq. 2020, 318, 113885. [Google Scholar] [CrossRef]
  13. Higuchi, T.; Zuck, D.A. Investigation of Some Complexes Formed in Solution by Caffeine. J. Am. Pharm. Assoc. 1953, 42, 138–145. [Google Scholar] [CrossRef]
  14. Baranovskii, S.F.; Bolotin, P.A. Association of riboflavin, caffeine, and sodium salicylate in aqueous solution. J. Appl. Spectrosc. 2007, 74, 211–218. [Google Scholar] [CrossRef] [Green Version]
  15. Zeinalipour-Yazdi, C.D. A DFT study of the interaction of aspirin, paracetamol and caffeine with one water molecule. J. Mol. Model. 2022, 28, 285. [Google Scholar] [CrossRef]
  16. Muthuselvi, C.; Abirami, A.; Poornalatha, P.; Seethalakshmi, A. Growth and Characterization of Sodium Salicylate Single Crystal. Pharm. Chem. J. 2018, 5, 46–55. [Google Scholar]
  17. Malakar, S.; Gibson, P.R.; Barrett, J.S.; Muir, J.G. Naturally occurring dietary salicylates: A closer look at common Australian foods. J. Food Compos. Anal. 2017, 57, 31–39. [Google Scholar] [CrossRef]
  18. Kęszycka, P.K.; Szkop, M.; Gajewska, D. Overall Content of Salicylic Acid and Salicylates in Food Available on the European Market. J. Agric. Food Chem. 2017, 65, 11085–11091. [Google Scholar] [CrossRef]
  19. Scotter, M.J.; Roberts, D.P.T.; Wilson, L.A.; Howart, F.A.C.; Davis, J.; Mansell, N. Free salicylic acid and acetyl salicylic acid content of foods using gas chromatography-mass spectrometry. Food Chem. 2007, 105, 273–279. [Google Scholar] [CrossRef]
  20. Ankita; Nain, A.K. Solute-solute and solute-solvent interactions of drug sodium salicylate in aqueous-glucose/sucrose solutions at temperatures from 293.15 to 318.15 K: A physicochemical study. J. Mol. Liq. 2020, 298, 112006. [Google Scholar] [CrossRef]
  21. Weiser, T.W. Treatments, Mechanism and Adverse Reactions of Anesthetics and Analgesic; Academic Press: Cambridge, MA, USA, 2001. [Google Scholar]
  22. Weiser, T.; Weigmann, H. Effect of Caffeine on the Bioavailability and Pharmacokinetics of an Acetylsalicylic Acid-Paracetamol Combination: Results of a Phase I Study. Adv. Ther. 2019, 36, 597–607. [Google Scholar] [CrossRef] [Green Version]
  23. Kim, J.K.; Choi, M.S.; Yoo, H.H.; Kim, D.H. The Intake of Coffee Increases the Absorption of Aspirin in Mice by Modifying Gut Microbiome. Pharmaceutics 2022, 14, 746. [Google Scholar] [CrossRef]
  24. Romero, C.M.; Oviedo, C.A. Effect of temperature on the solubility of α-amino acids and α,ω-amino acids in water. J. Solut. Chem. 2013, 42, 1355–1362. [Google Scholar] [CrossRef]
  25. Krieger, E.; Koraimann, G.; Vriend, G. Increasing the precision of comparative models with YASARA NOVA—A self-parameterizing force field. Proteins 2002, 47, 393–402. [Google Scholar] [CrossRef] [PubMed]
  26. Sha, J.; Ma, T.; Zhao, R.; Zhang, P.; Sun, R.; Jiang, G.; Wan, Y.; He, H.; Yao, X.; Li, Y.; et al. The dissolution behaviour and apparent thermodynamic analysis of doxifluridine in twelve pure solvents at various temperatures. J. Chem. Thermodyn. 2020, 144, 106073. [Google Scholar] [CrossRef]
  27. Shakeel, F.; Imran, M.; Haq, N.; Alshehri, S.; Anwer, M.K. Synthesis, characterization and solubility determination of 6-Phenyl-pyridazin-3 (2 H)-one in different pharmaceutical solvents. Molecules 2019, 24, 3404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Masson, D.O. XXVIII. Solute molecular volumes in relation to solvation and ionization. Lond. Edinb. Dubl. Phil. Mag. 1929, 8, 218–235. [Google Scholar] [CrossRef]
  29. Hepler, L.G. Thermal expansion and structure in water and aqueous solutions. Can. J. Chem. 1969, 47, 4613–4617. [Google Scholar] [CrossRef]
  30. Birch, G.G.; Parke, S.; Siertsema, R.; Westwell, J.M. Importance of molar volumes and related parameters in sweet taste chemoreception. Pure App. Chem. 1997, 69, 685–692. [Google Scholar] [CrossRef] [Green Version]
  31. Jones, G.; Dole, M. The viscosity of aqueous solutions of strong electrolytes with special reference to barium chloride. J. Am. Chem. Soc. 1929, 51, 2950–2964. [Google Scholar] [CrossRef]
  32. Marcus, Y. Effect of ions on the structure of water: Structure making and breaking. Chem. Rev. 2009, 109, 1346–1370. [Google Scholar] [CrossRef]
  33. Feakins, D.; Freemantle, D.J.; Lawrence, K.G. Transition state of the relative viscosity of electrolytic solutions. J. Chem. Soc. Faraday Trans. 1974, 70, 795–806. [Google Scholar] [CrossRef]
  34. Shekaari, H.; Zafarani-Moattar, M.T.; Mirheydari, S.N. Thermodynamic study of aspirin in the presence of ionic liquid, 1-hexyl-3-methylimidazolium bromide in acetonitrile at T = (288.15 to 318.15) K. J. Mol. Liq. 2015, 209, 138–148. [Google Scholar] [CrossRef]
  35. Shekaari, H.; Zafarani-Moattar, M.T.; Mirheydari, S.N. Thermodynamic properties of 1-butyl-3-methylimidazolium salicylate as an active pharmaceutical ingredient ionic liquid (API-IL) in aqueous solutions of glycine and L-alanine at T = (288.15–318.15) K. Thermochim. Acta 2016, 637, 51–68. [Google Scholar] [CrossRef]
  36. Glasstone, S.; Laidler, K.J.; Eyring, H. The Theory of Rate Processes: The Kinetics of Chemical Reactions, Viscosity, Diffusion and Electrochemical Phenomena; McGraw-Hill: New York, NY, USA, 1941; p. 611. [Google Scholar]
  37. Millero, F.J.; Surdo, A.L.; Shin, C. The Apparent Molal Volumes and Adiabatic Compressibilities of Aqueous Amino Acids at 25 °C. J. Phys. Chem. 1978, 82, 784–792. [Google Scholar] [CrossRef]
  38. Sharma, B.; Paul, S. Effects of dilute aqueous NaCl solution on caffeine aggregation. J. Chem. Phys. 2013, 139, 194504. [Google Scholar] [CrossRef]
  39. Shimizu, S. Caffeine dimerization: Effects of sugar, salts, and water structure. Food Funct. 2015, 6, 3228–3235. [Google Scholar] [CrossRef]
  40. Páez, M.; Figueredo, S.; Pérez, D.; Vergara, M.; Lans, E. Volumetric, viscometric and molecular simulation studies of glycine in aqueous sodium sulphate solutions at different temperatures. J. Mol. Liq. 2018, 266, 718–726. [Google Scholar] [CrossRef]
  41. Begard, J.; Jadot, R. Solvation in aqueous N,N-dimethylformamide. J. Chem. Thermodyn. 1975, 7, 1185–1188. [Google Scholar] [CrossRef]
  42. Thomas, D.G. Transport characteristics of suspension: VIII. A Note on the viscosity of Newtonian suspension of uniform spherical particles. J. Colloid Sci. 1965, 20, 267–271. [Google Scholar] [CrossRef]
Figure 1. Structure of caffeine molecule and atom-numbering scheme.
Figure 1. Structure of caffeine molecule and atom-numbering scheme.
Pharmaceutics 14 02304 g001
Figure 2. Structure of sodium salicylate molecule and atom-numbering scheme.
Figure 2. Structure of sodium salicylate molecule and atom-numbering scheme.
Pharmaceutics 14 02304 g002
Figure 3. The plot of the natural logarithm of the mole fraction of (■) caffeine in a saturated 0.1 mol∙kg−1 SS aqueous solution and () caffeine in a saturated aqueous solution, lnx1 versus 1/T − 1/Thm.
Figure 3. The plot of the natural logarithm of the mole fraction of (■) caffeine in a saturated 0.1 mol∙kg−1 SS aqueous solution and () caffeine in a saturated aqueous solution, lnx1 versus 1/T − 1/Thm.
Pharmaceutics 14 02304 g003
Figure 4. The dependence of the apparent molar volume ( V ϕ ) on caffeine molality (m) in caffeine in 0.1 mol∙kg−1 SS aqueous solutions at different temperatures (T) = (■) 293.15, () 298.15, () 303.15, () 308.15, and () 313.15 K.
Figure 4. The dependence of the apparent molar volume ( V ϕ ) on caffeine molality (m) in caffeine in 0.1 mol∙kg−1 SS aqueous solutions at different temperatures (T) = (■) 293.15, () 298.15, () 303.15, () 308.15, and () 313.15 K.
Pharmaceutics 14 02304 g004
Figure 5. Plot of reduced viscosity (η/ηo − 1) versus concentration (c) of the caffeine in SS aqueous solutions, at different temperatures, T = (■) 293.15, () 298.15, () 303.15, () 308.15, and () 313.15 K.
Figure 5. Plot of reduced viscosity (η/ηo − 1) versus concentration (c) of the caffeine in SS aqueous solutions, at different temperatures, T = (■) 293.15, () 298.15, () 303.15, () 308.15, and () 313.15 K.
Pharmaceutics 14 02304 g005
Figure 6. Visual representation obtained from the Yasara Structure program of caffeine organization in aqueous solution (a) and 0.1 mol∙kg−1 SS aqueous solution (b), and blue dashed lines represent π–π interactions between caffeine molecules.
Figure 6. Visual representation obtained from the Yasara Structure program of caffeine organization in aqueous solution (a) and 0.1 mol∙kg−1 SS aqueous solution (b), and blue dashed lines represent π–π interactions between caffeine molecules.
Pharmaceutics 14 02304 g006
Figure 7. RDFs for the caffeine–caffeine center of mass (black line—aqueous solution, red line—0.1 mol∙kg−1 SS aqueous solution). The lines represent the radial distribution functions obtained from MD simulations’ trajectories.
Figure 7. RDFs for the caffeine–caffeine center of mass (black line—aqueous solution, red line—0.1 mol∙kg−1 SS aqueous solution). The lines represent the radial distribution functions obtained from MD simulations’ trajectories.
Pharmaceutics 14 02304 g007
Figure 8. RDFs for caffeine-sodium salicylate center of mass. The lines represent the radial distribution functions obtained from MD simulations’ trajectories. The black line represents RDF for caffeine-sodium salicylate center of mass at 0.06 mol∙kg−1 and the red line represents RDF for caffeine-sodium salicylate center of mass at 0.25 mol∙kg−1.
Figure 8. RDFs for caffeine-sodium salicylate center of mass. The lines represent the radial distribution functions obtained from MD simulations’ trajectories. The black line represents RDF for caffeine-sodium salicylate center of mass at 0.06 mol∙kg−1 and the red line represents RDF for caffeine-sodium salicylate center of mass at 0.25 mol∙kg−1.
Pharmaceutics 14 02304 g008
Figure 9. Radial distribution functions of water oxygen atom around selected atoms of caffeine in aqueous solution (black lines) and the presence of SS (red lines): (a) O1-Ow, (b) O2-Ow, (c) N4-Ow, and (d) H1-Ow.
Figure 9. Radial distribution functions of water oxygen atom around selected atoms of caffeine in aqueous solution (black lines) and the presence of SS (red lines): (a) O1-Ow, (b) O2-Ow, (c) N4-Ow, and (d) H1-Ow.
Pharmaceutics 14 02304 g009aPharmaceutics 14 02304 g009b
Figure 10. Visual representation obtained from the Yasara Structure program of caffeine organization in 0.1 mol∙kg−1 SS aqueous solution, and blue dashed lines represent π–π interactions. The molal concentration of caffeine is 0.06 mol∙kg−1.
Figure 10. Visual representation obtained from the Yasara Structure program of caffeine organization in 0.1 mol∙kg−1 SS aqueous solution, and blue dashed lines represent π–π interactions. The molal concentration of caffeine is 0.06 mol∙kg−1.
Pharmaceutics 14 02304 g010
Figure 11. Visual representation of the potential role of salicylates in increasing the solubility of caffeine: (a) adding caffeine in SS solution, (b) forming π–π interactions between caffeine and salicylate anions, (c) caffeine and salicylate complexes that provide better solubility of caffeine monomers, and (d) forming caffeine–caffeine complexes through π–π interactions with salicylate anion release.
Figure 11. Visual representation of the potential role of salicylates in increasing the solubility of caffeine: (a) adding caffeine in SS solution, (b) forming π–π interactions between caffeine and salicylate anions, (c) caffeine and salicylate complexes that provide better solubility of caffeine monomers, and (d) forming caffeine–caffeine complexes through π–π interactions with salicylate anion release.
Pharmaceutics 14 02304 g011
Table 1. Determined solubility of caffeine in water and caffeine in 0.1 mol∙kg1 SS aqueous solutions in the temperature range from T = 291.15 to 313.15 K.
Table 1. Determined solubility of caffeine in water and caffeine in 0.1 mol∙kg1 SS aqueous solutions in the temperature range from T = 291.15 to 313.15 K.
T/(K)Caffeine Solubility/(g (1000 gwater)−1)
Caffeine in WaterCaffeine in SS Aqueous Solution
293.1516.3343.83
298.1520.7149.26
303.1525.2955.73
308.1534.8368.62
313.1543.9089.22
Table 2. The obtained apparent thermodynamic values of Δsol, Δsol, and Δsol, as well as ζH and ζTS for caffeine in water and SS aqueous solutions.
Table 2. The obtained apparent thermodynamic values of Δsol, Δsol, and Δsol, as well as ζH and ζTS for caffeine in water and SS aqueous solutions.
Caffeine in Aqueous SolutionCaffeine in 0.1 mol∙kg−1 SS Aqueous Solution
ΔsolG°/(kJ·mol−1)13.3013.28
ΔsolH°/(kJ·mol−1)0.410.13
ΔsolS°/(J·K−1·mol−1)−0.0425−0.0434
ζH0.0310.009
ζTS0.960.99
Table 3. Fitting parameters obtained using the Masson′s equation modified for non-electrolytes of caffeine in SS aqueous solutions in the temperature range from T = 293.15 to 313.15 K, with the standard deviations (σ) and regression coefficients (R2).
Table 3. Fitting parameters obtained using the Masson′s equation modified for non-electrolytes of caffeine in SS aqueous solutions in the temperature range from T = 293.15 to 313.15 K, with the standard deviations (σ) and regression coefficients (R2).
T/(K) V ϕ o ( cm 3 · mol 1 ) Sv/(cm3·kg·mol−2)σ (cm3∙mol−1)R2
293.15131.4675.460.5030.9497
298.15132.4773.660.2340.9882
303.15133.6272.560.3040.9796
308.15134.1574.280.3750.9706
313.15134.3382.780.2370.9903
Table 4. Values of the limiting apparent molar expansibilities, E ϕ o , and Hepler’s coefficients of caffeine in SS aqueous solutions in the temperature range from T = 283.15 to 313.15 K.
Table 4. Values of the limiting apparent molar expansibilities, E ϕ o , and Hepler’s coefficients of caffeine in SS aqueous solutions in the temperature range from T = 283.15 to 313.15 K.
T/(K)293.15298.15303.15308.15313.15 E ϕ o T p /(cm3·mol−1·K−2)
E ϕ o /(cm3·mol−1·K−1)
m(SS)/(mol·kg−1)
0.00000.16900.16090.15280.14470.1366−0.0016 [12]
0.10000.27910.21370.14840.08310.0178−0.0130
Table 5. Limiting apparent transfer molar volumes, Δ tr v ϕ o , and the apparent specific volume at infinite dilution, v ϕ o , of caffeine in 0.1 mol∙kg−1 SS aqueous solutions in the temperature range from T = 293.15 to 313.15 K.
Table 5. Limiting apparent transfer molar volumes, Δ tr v ϕ o , and the apparent specific volume at infinite dilution, v ϕ o , of caffeine in 0.1 mol∙kg−1 SS aqueous solutions in the temperature range from T = 293.15 to 313.15 K.
T/(K) Δ tr v ϕ o /(cm3·mol−1) v ϕ o /(cm3·g−1)
293.15−11.470.68
298.15−11.060.68
303.15−10.830.69
308.15−11.050.69
313.15−11.580.69
Table 6. Values of viscosity B-coefficients obtained from the Jones-Dole’s equation for caffeine in aqueous solution [12] and in 0.1 mol∙kg−1 SS aqueous solutions in the temperature range from T = 293.15 to 313.15 K.
Table 6. Values of viscosity B-coefficients obtained from the Jones-Dole’s equation for caffeine in aqueous solution [12] and in 0.1 mol∙kg−1 SS aqueous solutions in the temperature range from T = 293.15 to 313.15 K.
T/(K)293.15298.15303.15308.15313.15
m(SS)/(mol·kg−1)B/(dm3·mol−1)
0.10000.5010.4900.4610.4490.436
0.0000 [11]0.4140.4210.4290.4320.435
Table 7. Thermodynamic parameters V ¯ 1 o , V ¯ 2 o , Δ μ 1 o , and Δ μ 2 o of viscous flow values of caffeine in 0.1 mol∙kg−1 SS aqueous solution in the temperature range from T = 293.15 to 313.15 K.
Table 7. Thermodynamic parameters V ¯ 1 o , V ¯ 2 o , Δ μ 1 o , and Δ μ 2 o of viscous flow values of caffeine in 0.1 mol∙kg−1 SS aqueous solution in the temperature range from T = 293.15 to 313.15 K.
ParametersT/(K)
293.15298.15303.15308.15313.15
V ¯ 2 o 131.46132.47133.62134.15134.33
V ¯ 1 o 18.1818.2118.2318.2618.30
Δ μ 1 o 9.429.289.159.038.92
Δ μ 2 o 91.7791.6288.7888.3187.43
Table 8. Values of caffeine hydration numbers gathered from volumetry and viscosimetric measurements in aqueous solution and 0.1 mol∙kg−1 SS aqueous solutions in the temperature range from T = 283.15 to 313.15 K.
Table 8. Values of caffeine hydration numbers gathered from volumetry and viscosimetric measurements in aqueous solution and 0.1 mol∙kg−1 SS aqueous solutions in the temperature range from T = 283.15 to 313.15 K.
T/(K)hn (Volumetry)hn (Viscosimetry)hn (for Caffeine in Water, Volumetry [11])
293.154.234.170.94
298.153.673.280.68
303.153.042.580.38
308.152.612.170.16
313.152.231.39−0.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vraneš, M.; Borović, T.T.; Drid, P.; Trivić, T.; Tomaš, R.; Janković, N. Influence of Sodium Salicylate on Self-Aggregation and Caffeine Solubility in Water—A New Hypothesis from Experimental and Computational Data. Pharmaceutics 2022, 14, 2304. https://doi.org/10.3390/pharmaceutics14112304

AMA Style

Vraneš M, Borović TT, Drid P, Trivić T, Tomaš R, Janković N. Influence of Sodium Salicylate on Self-Aggregation and Caffeine Solubility in Water—A New Hypothesis from Experimental and Computational Data. Pharmaceutics. 2022; 14(11):2304. https://doi.org/10.3390/pharmaceutics14112304

Chicago/Turabian Style

Vraneš, Milan, Teona Teodora Borović, Patrik Drid, Tatjana Trivić, Renato Tomaš, and Nenad Janković. 2022. "Influence of Sodium Salicylate on Self-Aggregation and Caffeine Solubility in Water—A New Hypothesis from Experimental and Computational Data" Pharmaceutics 14, no. 11: 2304. https://doi.org/10.3390/pharmaceutics14112304

APA Style

Vraneš, M., Borović, T. T., Drid, P., Trivić, T., Tomaš, R., & Janković, N. (2022). Influence of Sodium Salicylate on Self-Aggregation and Caffeine Solubility in Water—A New Hypothesis from Experimental and Computational Data. Pharmaceutics, 14(11), 2304. https://doi.org/10.3390/pharmaceutics14112304

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop