Letizia Tavagnacco

Molecular Dynamics Simulation Studies of Caffeine Aggregation in Aqueous Solution

Molecular dynamics simulations were carried out on a system of eight independent caffeine molecules in a periodic box of water at 300 K, representing a solution near the solubility limit for caffeine at room temperature, using a newly developed CHARMM-type force field for caffeine in water. Simulations were also conducted for single caffeine molecules in water using two different water models (TIP3P and TIP4P). Water was found to structure in a complex fashion around the planar caffeine molecules, which was not sensitive to the water model used. As expected, extensive aggregation of the caffeine molecules was observed, with the molecules stacking their flat faces against one another like coins, with their methylene groups staggered to avoid steric clashes. A dynamic equilibrum was observed between large n-mers, including stacks with all eight solute molecules, and smaller clusters, with the calculated osmotic coefficient being in acceptable agreement with the experimental value. The insensitivity of the results to water model and the congruence with experimental thermodynamic data suggest that the observed stacking interactions are a realistic representation of the actual association mechanism in aqueous caffeine solutions.

Sugar-binding sites on the surface of the carbohydrate-binding module of CBH I from Trichoderma reesei.

Molecular dynamics simulations were carried out for a system consisting of the carbohydrate-binding module (CBM) of the cellulase CBH I from Trichoderma reesei (Hypocrea jecorina) in a concentrated solution of β-D-glucopyranose, to determine whether there is any tendency for the sugar molecules to bind to the CBM. In spite of the general tendency of glucose to behave as an osmolyte, a marked tendency for the sugar molecules to bind to the protein was observed. However, the glucose molecules tended to bind only to specific sites on the protein. As expected, the hydrophobic face of the sugar molecules, comprising the axial H1, H3, and H5 aliphatic protons, tended to adhere to the flat faces of the three tyrosine side chains on the planar binding surface of the CBM. However, a significant tendency to bind to a groove-like feature on the upper surface of the CBM was also observed. These results would not be inconsistent with a model of the mechanism for this globular domain in which the cellodextrin chain being removed from the surface of crystalline cellulose passes over the upper surface of the CBM, presumably then available for hydrolysis in the active site tunnel of this processive cellulase.

Weakly hydrated surfaces and the binding interactions of small biological solutes

Extended planar hydrophobic surfaces, such as are found in the side chains of the amino acids histidine, phenylalanine, tyrosine, and tryptophan, exhibit an affinity for the weakly hydrated faces of glucopyranose. In addition, molecular species such as these, including indole, caffeine, and imidazole, exhibit a weak tendency to pair together by hydrophobic stacking in aqueous solution. These interactions can be partially understood in terms of recent models for the hydration of extended hydrophobic faces and should provide insight into the architecture of sugar-binding sites in proteins.

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

The solvation of caffeine in water is examined with neutron diffraction experiments at 353 K. The experimental data, obtained by taking advantage of isotopic H/D substitution in water, were analyzed by empirical potential structure refinement (EPSR) in order to extract partial structure factors and site-site radial distribution functions. In parallel, molecular dynamics (MD) simulations were carried out to interpret the data and gain insight into the intermolecular interactions in the solutions and the solvation process. The results obtained with the two approaches evidence differences in the individual radial distribution functions, although both confirm the presence of caffeine stacks at this temperature. The two approaches point to different accessibility of water to the caffeine sites due to different stacking configurations.

Myelography Iodinated Contrast Media. I. Unraveling the Atropisomerism Properties in Solution

The present work reports a thorough conformational analysis of iodinated contrast media: iomeprol, iopamidol (the worlds most utilized contrast agent), and iopromide. Its main aim is the understanding of the complex structural features of these atropisomeric molecules, characterized by the presence of many conformers with hindered rotations, and of the role of atropisomerism in the physicochemical properties of their aqueous solutions. The problem was tackled by using an extensive analysis of (13)C NMR data on the solutions of whole molecules and of simple precursors in addition to FT-IR investigation and molecular simulations. This analysis demonstrated that out of the many possible atropisomers, only a few are significantly populated, and their relative population is provided. The conformational analysis also indicated that the presence of a sterically hindered amidic bond, allowing a significant population of cis forms (E in iopromide and exo in iomeprol), may be the basis for an increased thermodynamic solubility of concentrated solutions of iomeprol.

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.

Myelography iodinated contrast media. 2. Conformational versatility of iopamidol in the solid state

The phenomenon of polymorphism is of great relevance in pharmaceutics, since different polymorphs have different physicochemical properties, e.g., solubility, hence, bioavailability. Coupling diffractometric and spectroscopic experiments with thermodynamic analysis and computational work opens to a methodological approach which provides information on both structure and dynamics in the solid as well as in solution. The present work reports on the conformational changes in crystalline iopamidol, which is characterized by atropisomerism, a phenomenon that influences both the solution properties and the distinct crystal phases. The conformation of iopamidol is discussed for three different crystal phases. In the anhydrous and monohydrate crystal forms, iopamidol molecules display a syn conformation of the long branches stemming out from the triiodobenzene ring, while in the pentahydrate phase the anti conformation is found. IR and Raman spectroscopic studies carried out on the three crystal forms, jointly with quantum chemical computations, revealed that the markedly different spectral features can be specifically attributed to the different molecular conformations. Our results on the conformational versatility of iopamidol in different crystalline phases, linking structural and spectroscopic evidence for the solution state and the solid forms, provide a definite protocol for grasping the signals that can be taken as conformational markers. This is the first step for understanding the crystallization mechanism occurring in supersaturated solution of iopamidol molecules.

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.