Hydration of Caffeine at High Temperature by Neutron Scattering and Simulation Studies

Salt Effects on Caffeine across Concentration Regimes

Salts affect the solvation thermodynamics of molecules of all sizes; the Hofmeister series is a prime example in which different ions lead to salting-in or salting-out of aqueous proteins. Early work of Tanford led to the discovery that the solvation of molecular surface motifs is proportional to the solvent accessible surface area (SASA), and later studies have shown that the proportionality constant varies with the salt concentration and type. Using multiscale computer simulations combined with vapor-pressure osmometry on caffeine-salt solutions, we reveal that this SASA description captures a rich set of molecular driving forces in tertiary solutions at changing solute and osmolyte concentrations. Central to the theoretical work is a new potential energy function that depends on the instantaneous surface area, salt type, and concentration. Used in, e.g., Monte Carlo simulations, this allows for a highly efficient exploration of many-body interactions and the resulting thermodynamics at elevated solute and salt concentrations.

Elucidating the hydrotropism behaviour of aqueous caffeine and sodium benzoate solution through NMR and neutron total scattering analysis

Hydrotropism is a convenient way to increase the solubility of drugs by up to several orders of magnitude, and even though it has been researched for decades with both experimental and simulation methods, its mechanism is still unknown. Here, we use caffeine/sodium benzoate (CAF-SB) as model system to explore the behaviour of caffeine solubility enhancement in water through NMR spectroscopy and neutron total scattering. 1H NMR shows strong interaction between caffeine and sodium benzoate in water. Neutron total scattering combined with empirical potential structure refinement, a systematic method to study the solution structure, reveals π-stacking between caffeine and the benzoate anion as well as Coulombic interactions with the sodium cation. The strongest hydrogen bond interaction in the system is between benzoate and water, which help dissolve CAF-SB complex and increase the solubility of CAF in water. Besides, the stronger interaction between CAF and water and the distortion of water structure are further mechanisms of the CAF solubility enhancement. It is likely that the variety of mechanisms for hydrotropism shown in this system can be found for other hydrotropes, and NMR spectroscopy and neutron total scattering can be used as complementary techniques to generate a holistic picture of hydrotropic solutions.

Using Prenucleation Aggregation of Caffeine-Benzoic Acid as a Rapid Indication of Co-crystallization from Solutions

Co-crystal design is a convenient way to remedy the poor biopharmaceutical properties of drugs. Most studies focus on experimental co-crystal screening or computational prediction, but hardly any work has been done toward fast, efficient, and reliable prediction of solution crystallization for co-crystal formation. Here, we study the caffeine-benzoic acid co-crystal system, due to its reported difficulty to crystallize from the solution phase. With this work, we investigate whether there is a link between prenucleation aggregation in solution and co-crystal formation and how to harness this for crystallization prediction. ¹H and ¹³C NMR spectroscopy is used to study the prenucleation interaction between caffeine and benzoic acid in methanol, acetone, and acetonitrile as examples of common solvents. In this system, crystallization from methanol leads to no co-crystallization, from acetone to concomitant crystallization of co-crystal and caffeine, and from acetonitrile to pure co-crystal formation from solution. Strong heteromeric dimers were found to exist in all three solvents. Ternary phase diagrams were defined and a solution-accessible co-crystal region was found for all solvents. For this system, the prenucleation clusters found in solution could be linked to the crystallization of the co-crystal. Crystallization from DMSO did not yield the co-crystal and there were no detectable prenucleation aggregates. NMR spectroscopy to probe dimers in solution can thus be used as a fast, reliable, and promising tool to predict co-crystallization from specific solvents and to screen for suitable solvents for manufacturing and scale-up.

Anion-cation contrast of small molecule solvation in salt solutions

The contributions from anions and cations from salt are inseparable in their perturbation of molecular systems by experimental and computational methods, rendering it difficult to dissect the effects exerted by the anions and cations individually. Here we investigate the solvation of a small molecule, caffeine, and its perturbation by monovalent salts from various parts of the Hofmeister series. Using molecular dynamics and the energy-representation theory of solvation, we estimate the solvation free energy of caffeine and decompose it into the contributions from anions, cations, and water. We also decompose the contributions arising from the solute-solvent and solute-ions interactions and that from excluded volume, enabling us to pin-point the mechanism of salt. Anions and cations revealed high contrast in their perturbation of caffeine solvation, with the cations salting-in caffeine via binding to the polar ketone groups, while the anions were found to be salting-out via perturbations of water. In agreement with previous findings, the perturbation by salt is mostly anion dependent, with the magnitude of the excluded-volume effect found to be the governing mechanism. The free-energy decomposition as conducted in the present work can be useful to understand ion-specific effects and the associated Hofmeister series.

Interaction of Caffeine with Model Lipid Membranes

Caffeine is not only a widely consumed active stimulant, but it is also a model molecule commonly used in pharmaceutical sciences. In this work, by performing quartz-crystal microbalance and neutron reflectometry experiments we investigate the interaction of caffeine molecules with a model lipid membrane. We determined that caffeine molecules are not able to spontaneously partition from an aqueous environment, enriched in caffeine, into a bilayer. Caffeine could be however included in solid-supported lipid bilayers if present with lipids during self-assembly. In this case, thanks to surface-sensitive techniques, we determined that caffeine molecules are preferentially located in the hydrophobic region of the membrane. These results are highly relevant for the development of new drug delivery vectors, as well as for a deeper understanding of the membrane permeation role of purine molecules.

Hydration of Carboxyl Groups: A Route toward Molecular Recognition?

On Earth, water plays an active role in cellular life, over several scales of distance and time. At a nanoscale, water drives macromolecular conformation through hydrophobic forces and at short times acts as a proton donor/acceptor providing charge carriers for signal transmission. At longer times and larger distances, water controls osmosis, transport, and protein mobility. Neutron diffraction experiments augmented by computer simulation, show that the three-dimensional shape of the hydration shell of carboxyl and carboxylate groups belonging to different molecules is characteristic of each molecule. Different hydration shells identify and distinguish specific sites with the same chemical structure. This experimental evidence suggests an active role of water also in controlling, modulating, and mediating chemical reactions involving carboxyl and carboxylate groups.

Ligand binding to G-quadruplex DNA: new insights from ultraviolet resonance Raman spectroscopy

G-Quadruplexes (G4s) are noncanonical nucleic acid structures involved in the regulation of several biological processes of many organisms. The rational design of G4-targeting molecules developed as potential anticancer and antiviral therapeutics is a complex problem intrinsically due to the structural polymorphism of these peculiar DNA structures. The aim of the present work is to show how Ultraviolet Resonance Raman (UVRR) spectroscopy can complement other techniques in providing valuable information about ligand/G4 interactions in solution. Here, the binding of BRACO-19 and Pyridostatin - two of the most potent ligands - to selected biologically relevant G4s was investigated by polarized UVRR scattering at 266 nm. The results give new insights into the binding mode of these ligands to G4s having different sequences and topologies by performing an accurate analysis of peaks assigned to specific groups and their changes upon binding. Indeed, the UVRR data not only show that BRACO-19 and Pyridostatin interact with different G4 sites, but also shed light on the ligand and G4 chemical groups really involved in the interaction. In addition, UVRR results complemented by circular dichroism data clearly indicate that the binding mode of a ligand can also depend on the conformation(s) of the target G4. Overall, these findings demonstrate the utility of using UVRR spectroscopy in the investigation of G4s and G4-ligand interactions in solution.

Role of Water in Sucrose, Lactose, and Sucralose Taste: The Sweeter, The Wetter?

Natural sugars combine energy supply and, except a few cases, a pleasant taste. On the other hand, exaggerated consumption may impact population health. This has busted the research for the synthesis of increasingly cheaper artificial sweeteners, with low energy content and intense taste. Here, we suggest that studies of the hydration properties of three disaccharides, namely, the natural sucrose and lactose and the artificial sucralose, may explain the difference by orders of magnitude among their sweetness. This is done by analyzing via Monte Carlo simulations the neutron diffraction differential cross sections of aqueous solutions of the three sugars and their isotopes. Our results show that the strength of the sugar-water hydrogen bond interaction is one of the factors influencing sweetness, another being the number of water molecules within the first neighboring shell of the sugar whether bonded or not.

Solubility of Caffeine in Supercritical CO2: A Molecular Dynamics Simulation Study

The extraction of caffeine from green tea leaves and cocoa beans is a common industrial process for the production of decaffeinated beverages and pharmaceuticals. The choice of the solvent critically determines the yield of this extraction process. Being an environmentally benign and recyclable solvent, supercritical carbon dioxide (scCO₂) has emerged as the most desirable green solvent for caffeine extraction. The present study investigates the solvation properties of caffeine in scCO₂ at two different temperatures (318 and 350 K) using molecular dynamics simulations. Unlike in water, the caffeine molecules in scCO₂ do not aggregate to form clusters due to relatively stronger caffeine-CO₂ interactions. A well-structured scCO₂ solvent shell envelops each caffeine molecule as a result of strong electron-donor-acceptor (EDA) and hydrogen-bonding interactions between these two species. Upon heating, although marginal site-specific changes in the distribution of nearest CO₂ around caffeine are observed, the overall distribution is retained. At a higher temperature, the caffeine-CO₂ hydrogen-bonding interactions are weakened, while their EDA interactions become relatively stronger. The results underscore the importance of the interplay of these interactions in determining stable solvent structures and solubility of caffeine in scCO₂.

Pre-nucleation aggregation based on solvent microheterogeneity

The microheterogeneous region of aqueous acetonitrile leads to preferred localisation and aggregation of caffeine and theophylline on the interface.

Hydrogen Bond Length as a Key To Understanding Sweetness

Neutron diffraction experiments have been performed to investigate and compare the structure of the hydration shell of three monosaccharides, namely, fructose, glucose, and mannose. It is found that despite their differences with respect to many thermodynamical quantities, bioprotective properties against environmental stresses, and taste, the influence of these monosaccharides on the bulk water solvent structure is virtually identical. Conversely, these sugars interact with the neighboring water molecules by forming H bonds of different length and strength. Interestingly, the sweetness of these monosaccharides, along with that of the disaccharide trehalose, is correlated with the length of these H bonds. This suggests that the small differences in stereochemistry between the different sugars determine a relevant change in polarity, which has a fundamental impact on the behavior of these molecules in vivo.

Molecular Dynamics and Neutron Scattering Studies of Mixed Solutions of Caffeine and Pyridine in Water

Insight into the molecular interactions of homotactic and heterotactic association of caffeine and pyridine in aqueous solution is given on the basis of both experimental and simulation studies. Caffeine is about 5 times more soluble in a 3 m aqueous pyridine solution than it is in pure water (an increase from ∼0.1 m to 0.5 m). At this elevated concentration the system becomes suitable for neutron scattering study. Caffeine-pyridine interactions were studied by neutron scattering and molecular dynamics simulations, allowing a detailed characterization of the spatial and orientational structure of the solution. It was found that while pyridine-caffeine interactions are not as strong as caffeine-caffeine interactions, the pyridine-caffeine interactions still significantly disrupted caffeine-caffeine stacking. The alteration of the caffeine-caffeine stacking, occasioned by the presence of pyridine molecules in solution and the consequent formation of heterotactic interactions, leads to the experimentally detected increase in caffeine solubility.

Stacking and Branching in Self-Aggregation of Caffeine in Aqueous Solution: From the Supramolecular to Atomic Scale Clustering

The dynamical and structural properties of caffeine solutions at the solubility limit have been investigated as a function of temperature by means of MD simulations, static and dynamic light scattering, and small angle neutron scattering experiments. A clear picture unambiguously supported by both experiment and simulation emerges: caffeine self-aggregation promotes the formation of two distinct types of clusters: linear aggregates of stacked molecules, formed by 2-14 caffeine molecules depending on the thermodynamic conditions and disordered branched aggregates with a size in the range 1000-3000 Å. While the first type of association is well-known to occur under room temperature conditions for both caffeine and other purine systems, such as nucleotides, the presence of the supramolecular aggregates has not been reported previously. MD simulations indicate that branched structures are formed by caffeine molecules in a T-shaped arrangement. An increase of the solubility limit (higher temperature but also higher concentration) broadens the distribution of cluster sizes, promoting the formation of stacked aggregates composed by a larger number of caffeine molecules. Surprisingly, the effect on the branched aggregates is rather limited. Their internal structure and size do not change considerably in the range of solubility limits investigated.

Atomic scale insights into urea–peptide interactions in solution

Investigations on the β-turn forming peptide, GPG, suggest that urea denatures proteins by replacing water molecules and subsequently weakening the peptide bonds as a possible mechanism of protein denaturation by urea.

The solvation structure of alprazolam

Alprazolam is a benzodiazepine that is commonly prescribed for the treatment of anxiety and other related disorders.

Stacking of purines in water: the role of dipolar interactions in caffeine

Concentration dependence of the NCE and the dephasing time show that caffeine molecules aggregate at 80 °C by planar stacking with a relevant contribution of dipole interactions.