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.

Forty-one conformations of alpha- and beta-d-galactopyranose were geometry optimized using the B3LYP density functional and 6-311++G** basis set. Full geometry optimization was performed on different ring geometries and different hydroxymethyl rotamers (gg/gt/tg). Analytically derived Hessians were used to calculate zero point energy, enthalpy, and entropy. The lowest energy and free-energy conformation found is the alpha-gg-(4)C(1)-c chair conformation, which is of lower electronic and free energy than the lowest energy alpha-d-glucopyranose conformer because of favorable hydrogen-bonding interactions. The in vacuo calculations showed considerable ( approximately 2.2kcal/mol) energetic preference for the alpha over the beta anomer for galactopyranose in both the (4)C(1) and (1)C(4) chair conformations. Results are compared to glucopyranose and mannopyranose calculations in vacuo. Boat and skew-boat forms were found that remained stable upon gradient optimization, although many starting conformations moved to other boat forms upon optimization. As with glucopyranose and mannopyranose, the orientation and interaction of the hydroxyl groups make the most significant contributions to the conformation-energy relationship in vacuo.

No abstract

The result of DFT geometry optimization of 68 unique a-maltose conformers at the B3LYP/6-311++G** level of theory is described. Particular attention is paid to the hydroxyl group rotational positions and their influence on the glycosidic bond dihedral angles. The orientation of lone pair electrons across the bridging hydrogen bonds are implicated in directing the glycosidic dihedral angles with for example, conformers gg-gg and gt-gt, having different minimum energy conformations for the clockwise (c) and the reverse clockwise (r) forms. Conformers tg-gg, gg-tg, tg-tg, gt-gg, and gg-gt were studied, to understand the intermediate glycosidic bond conformations. The conformation, tg-gg-c, was found to be the lowest energy structure. When the hydroxyl groups on each glucose residue were made to point in opposite directions, i.e., c/r and r/c, the optimized structures were found to have high relative energies. Several optimized 'kink' structures were found around (/ H , w H )~(-40°, -40°), the lowest relative energy conformation being~3 kcal/mol. ''Kink'' conformations are observed in crystalline CA-10 and CA-14mers. Bandflip conformations, also observed in X-ray structures of CA-26 fragments, were studied with the lowest energy amaltose conformations~4.0 kcal/mol above the global energy minimum. Several trends in geometry resulting from hydroxyl rotamer directions are described.Keywords Conformation Á DFT Á B3LYP/6-311++G** Á a-Maltose Á Band flip Á Kinks Á Exocyclic hydroxyls Electronic supplementary material The online version of this article (

Molecular dynamics simulations were carried out on several systems of caffeine interacting with simple sugars. These included a single caffeine molecule in a 3 molal solution of α-D-glucopyranose, at a caffeine concentration of 0.083 molal; a single caffeine in a 3 molal solution of β-D-glucopyranose, and a single caffeine molecule in a 1.08 molal solution of sucrose (table sugar). Parallel Nuclear Magnetic Resonance titration experiments were carried out on the same solutions under similar conditions. Consistent with previous thermodynamic experiments, the sugars were found to have an affinity for the caffeine molecules in both the simulations and experiments, and that the binding in these complexes occurs by face-to-face stacking of the hydrophobic triad of protons of the pyranose rings against the caffeine face, rather than by hydrogen bonding. For the disaccharide, the binding occurs via stacking of the glucose ring against the caffeine, with a lesser affinity for the fructose observed. These findings are consistent with the association being driven by hydrophobic hydration, and are similar to the previously observed binding of glucose rings to various other planar molecules, including indole, serotonin, and phenol.

No abstract

Flea beetles in several genera are known to possess male-specific sesquiterpenes, at least some of which serve as aggregation pheromones that attract both sexes. In continuing research on the chemical ecology of Phyllotreta flea beetles, six new male-specific sesquiterpenes were identified, one from P. striolata (hydroxyketone 9) and five from P. pusilla (aldehydes 10-12 and 14 and alcohol 13); both species are crop pests. The minute amounts from beetles provided mass spectra and chromatographic data but were insufficient for complete structure determination. However, it was discovered that the new compounds could all be produced by applying organic reactions to previously identified flea beetle sesquiterpenes, and the resulting, larger amounts of material permitted definitive structure analysis by NMR. Molecular modeling was used in conjunction with NMR to define relative configurations of several newly created stereogenic centers. The absolute configurations of natural 9-14 were established by chiral gas chromatography/mass spectrometry. In electrophysiological tests (GC-EAD) conducted with P. striolata, compound 9 was detected with high sensitivity by the beetle antennae, which is consistent with a pheromonal function. The research opens new possibilities for using behavioral chemicals to monitor or manage these pest species.

No abstract