A Competition between Relative Stability and Binding Energy in Caffeine Phenyl-Glucose Aggregates: Implications in Biological Mechanisms.

Hydrogen bonds and stacking interactions are pivotal in biological mechanisms, although their proper characterisation within a molecular complex remains a difficult task. We used quantum mechanical calculations to characterise the complex between caffeine and phenyl-ß-D-glucopyranoside, in which several functional groups of the sugar derivative compete with each other to attract caffeine. Calculations at different levels of theory (M06-2X/6-311++G(d,p) and B3LYP-ED=GD3BJ/def2TZVP) agree to predict several structures similar in stability (relative energy) but with different affinity (binding energy). These computational results were experimentally verified by laser infrared spectroscopy, through which the caffeine·phenyl-ß-D-glucopyranoside complex was identified in an isolated environment, produced under supersonic expansion conditions. The experimental observations correlate with the computational results. Caffeine shows intermolecular interaction preferences that combine both hydrogen bonding and stacking interactions. This dual behaviour had already been observed with phenol, and now with phenyl-ß-D-glucopyranoside, it is confirmed and maximised. In fact, the size of the complex's counterparts affects the maximisation of the intermolecular bond strength because of the conformational adaptability given by the stacking interaction. Comparison with the binding of caffeine within the orthosteric site of the A2A adenosine receptor shows that the more strongly bound caffeine·phenyl-ß-D-glucopyranoside conformer mimics the interactions occurring within the receptor.

Exploring the Influence of Intermolecular Interactions in Prebiotic Chemistry Using Laser Spectroscopy and Calculations.

One of the most fascinating questions in chemistry is why nature chose CGAT as the alphabet of life. Very likely, such selection was the result of multiple factors and a long period of refinement. Here, we explore how the intermolecular interactions influenced such process, by characterizing the formation of dimers between adenine, theobromine and 4-aminopyrimidine. Using a combination of mass-resolved excitation spectroscopy and DFT calculations, we determined the structure of adenine-theobromine and 4-aminopyrimidine-theobromine dimers. The binding energy of these dimers is very close to the canonical adenine-thymine nucleobases. Likewise, the dimers are able to adopt Watson-Crick conformations. These findings seem to indicate that there were many options available to build the first versions of the informational polymers, which also had to compete with other molecules, such as 4-aminopyrimidine, which does not have a valid attaching point for a saccharide. For some reason, nature did not select the most strongly-bonded partners or if it did, such proto-bases were later replaced by the nowadays canonical CGAT.

Bond Length Alternation and Internal Dynamics in Model Aromatic Substituents of Lignin.

Broadband microwave spectra were recorded over the 2-18 GHz frequency range for a series of four model aromatic components of lignin; namely, guaiacol (ortho-methoxy phenol, G), syringol (2,6-dimethoxy phenol, S), 4-methyl guaiacol (MG), and 4-vinyl guaiacol (VG), under jet-cooled conditions in the gas phase. Using a combination of 13 C isotopic data and electronic structure calculations, distortions of the phenyl ring by the substituents on the ring are identified. In all four molecules, the rC(1)-C(6) bond between the two substituted C-atoms lengthens, leading to clear bond alternation that reflects an increase in the phenyl ring resonance structure with double bonds at rC(1)-C(2) , rC(3)-C(4) and rC(5)-C(6) . Syringol, with its symmetric methoxy substituents, possesses a microwave spectrum with tunneling doublets in the a-type transitions associated with H-atom tunneling. These splittings were fit to determine a barrier to hindered rotation of the OH group of 1975 cm-1 , a value nearly 50 % greater than that in phenol, due to the presence of the intramolecular OH⋅⋅⋅OCH3 H-bonds at the two equivalent planar geometries. In 4-methyl guaiacol, methyl rotor splittings are observed and used to confirm and refine an earlier measurement of the three-fold barrier V3 =67 cm-1 . Finally, 4-vinyl guaiacol shows transitions due to two conformers differing in the relative orientations of the vinyl and OH groups.

The Shapes of Sulfonamides: A Rotational Spectroscopy Study.

Benzenesulfonamides are a class of molecules of extreme interest in the biochemical field because many of them are active against a variety of diseases. In this work, the pharmacophoric group benzensulfonamide, its derivatives para-toluensulfonamide and ortho-toluensulfonamide, and the bioactive molecule sulfanilamide, were investigated using rotational spectroscopy to determine their conformations and the influence of different substituents on their structures. For all species, the hyperfine structure due to the 14N atom was analyzed, and this provided crucial information for the unambiguous identification of the observed conformation of all molecules. In addition, for ortho-toluensulfonamide, the vibration-rotation hyperfine structure related to the methyl torsion was analyzed, and the methyl group rotation barrier was determined. For benzensulfonamide, partial rS and r0 structures were established from the experimental rotational constants of the parent and two deuterated isotopic species. In all compounds except ortho-toluensulfonamide, the amino group of the sulfonamide group lies perpendicular to the benzene plane with the aminic hydrogens eclipsing the oxygen atoms. In ortho-toluensulfonamide, where weak attractive interactions occur between the nitrogen lone pair and the methyl hydrogen atoms, the amino group lies in a gauche orientation, retaining the eclipsed configuration with respect to the SO2 frame. A comparison of the geometrical arrangements found in the PDB database allowed us to understand that the bioactive conformations are different from those found in isolated conditions. The conformations within the receptor are reached with an energy cost, which is balanced by the interactions established in the receptor.

Characterizing hydrogen and tetrel bonds in clusters of CO2 with carboxylic acids.

The interaction between carbon dioxide and planar carboxylic acids has been investigated through the analysis of the microwave spectrum of the acrylic acid·CO2 complex and quantum chemical modeling of the R-COOH·(CO2)1,16 clusters, where R = H, CH2CH. As regards the 1 : 1 compounds, two species, involving the s-cis and s-trans conformers of acrylic acid were observed. For both of them, a similar bidentate interaction arises between the carbonyl group of CO2 and the carboxylic group of the organic acid, leading to the formation of a planar six-membered ring. The binding energy is estimated to be De ≃ 21 kJ mol-1, 1/3 being the energy contributions of the tetrel to hydrogen bonds, respectively. In the 1 : 16 clusters, the ring arrangement is broken, allowing for the interaction of the acid with several CO2 molecules. The CO2 molecules completely surround formic acid, whereas, in the case of acrylic acid, they tend to avoid the allyl chain.

Characterizing the lone pairπ-hole interaction in complexes of ammonia with perfluorinated arenes.

When hydrogen is completely replaced by fluorine, arenes become prone to forming a lone pairπ-hole non-covalent bond with ligands presenting electron rich regions. Such a species is ammonia, which confirms this behavior engaging its lone pair as the electron donor counterpart in the 1 : 1 adducts with hexafluorobenzene and pentafluoropyridine. In this work, the geometrical parameters of the interaction have been unambiguously identified through the detection, by means of Fourier transform microwave spectroscopy, of the rotational spectra of both normal species and their 15NH3 isotopologues. An accurate analysis of the experimental data, including internal dynamics effects, endorsed by quantum chemical calculations, both with topological analysis and energy decomposition method, extended to the hydrogenated arenes and their water complexes, proved the ability of ammonia to create a stronger and more flexible lone pairπ-hole interaction than water. Interestingly, the higher binding energies of the ammonia lone pairπ-hole interactions correspond to larger intermolecular distances.

The Role of Non-Covalent Interactions on Cluster Formation: Pentamer, Hexamers and Heptamer of Difluoromethane.

The role of non-covalent interactions (NCIs) has broadened with the inclusion of new types of interactions and a plethora of weak donor/acceptor partners. This work illustrates the potential of chirped-pulse Fourier transform microwave technique, which has revolutionized the field of rotational spectroscopy. In particular, it has been exploited to reveal the role of NCIs' in the molecular self-aggregation of difluoromethane where a pentamer, two hexamers and a heptamer were detected. The development of a new automated assignment program and a sophisticated computational screening protocol was essential for identifying the homoclusters in conditions of spectral congestion. The major role of dispersion forces leads to less directional interactions and more distorted structures than those found in polar clusters, although a detailed analysis demonstrates that the dominant interaction energy is the pairwise interaction. The tetramer cluster is identified as a structural unit in larger clusters, representing the maximum expression of bond between dimers.

Exploration of the theobromine-water dimer: comparison with DNA microhydration.

Understanding the molecular basis of the appearance of life on Earth is an exciting research field. Many factors may have influenced the election of the molecules used by living beings and evolution may have modified those original compounds. In an attempt to understand the role played by intermolecular interactions in the election of CGAT as the alphabet of life, we present here a thorough experimental and computational study on the interaction of theobromine with water. Theobromine is a xanthine derivative, structurally related to the nucleobases, and also present in many living beings. The experimental results demonstrate that the most stable isomer of theobromine-water was formed and detected in supersonic expansions. This isomer very well resembles the structure of the dimers between nucleobases and water, offering similar values of binding energy. A comparison between the results obtained for theobromine-water with those reported in the literature for monohydrates of nucleobases is also offered.

Chlorination and tautomerism: a computational and UPS/XPS study of 2-hydroxypyridine ⇌ 2-pyridone equilibrium.

The prototropic tautomeric equilibrium in 2-hydroxypyridine serves as a prototype model for the study of nucleobases' behaviour. The position of such an equilibrium in parent and chlorine monosubstituted 2-hydroxypyridine compounds in the gas phase was determined using synchrotron based techniques. The lactim tautomer is dominant for the 5- and 6-substituted compounds, whereas the parent, 3- and 4-substituted isomers have comparable populations for both tautomers. Information was obtained by measuring valence band and core level photoemission spectra at the chlorine L-edge and carbon, nitrogen, and oxygen K-edges. The effect of chlorine on the core ionization potentials of the atoms in the heterocycle was evaluated and reasonable agreement with a simple model was obtained. Basic considerations of resonance structures correctly predicts the tautomeric equilibrium for the 5- and 6-substituted compounds. The vibrationally resolved structure of the low energy portion of the valence band photoionization spectra is assigned based on quantum-chemical calculations of the neutral and charged species followed by simulation of the vibronic structure. It is shown that the first ionization occurs from a π orbital of similar shape for both tautomers. In addition, the highly distinctive vibronic structure observed just above the first ionization of the lactim, for three of the five species investigated, is assigned to the second ionization of the lactam.

Structure, Dynamics, and Accurate Laboratory Rotational Frequencies of the Acrylonitrile-Methanol Complex.

The hydrogen-bonded complex between acrylonitrile (CH2═CHCN) and methanol has been characterized spectroscopically in the millimeter wave range (59.6-74.4 GHz) using a free jet absorption millimeter wave spectrometer. Precise values of the rotational and centrifugal distortion constants were obtained from the measured frequencies of the complex of acrylonitrile with CH3OH and CD3OD. The analysis of the splittings of the rotational lines due to the hindered internal rotation of the methanol methyl group led to the determination of a V3 value of 221.9(7) and 218(5) cm-1 for the complexes of CH3OH and CD3OD, respectively, and these values are about 40% lower than that of free methanol. The structure of the observed conformation is in agreement with the global minimum determined at the MP2/aug-cc-pVTZ level of calculation, and the counterpoise corrected intermolecular binding energy, obtained at the same theoretical level, is De = 26.3 kJ mol-1.

Exploring Caffeine-Phenol Interactions by the Inseparable Duet of Experimental and Theoretical Data.

Intermolecular interactions are difficult to model, especially in systems formed by multiple interactions. Such is the case of caffeine-phenol. Structural data has been extracted by using mass-resolved excitation spectroscopy and double resonance techniques. Then the predictions of seven different computational methods have been explored to discover structural and energetic discrepancies between them that may even result in different assignments of the system. The results presented herein highlight the difficulty of constructing functionals to model systems with several competing interactions, and raise awareness of problems with assignments of complex systems with limited experimental information that rely exclusively on energetic data.

Bond Length Alternation Observed Experimentally: The Case of 1H-Indazole.

Bond length alternation is a chemical phenomenon in benzene rings fused to other rings, which has been mainly predicted theoretically. Its physical origin is still not clear and has generated discussion. Here, by using a strategy that combines microwave spectroscopy, custom-made synthesis and high-level ab initio calculations, we demonstrate that this phenomenon is clearly observed in the prototype indazole molecule isolated in the gas phase. The 1H-indazole conformer was detected by rotational spectroscopy, and its 17 isotopologues resulting from single and double heavy atom substitution (13 C and 15 N) were also unambiguously observed. Several experimental structures were determined and, in particular, the most useful semi-experimental equilibrium structure (re SE ), allowed determination of the heavy atom bond lengths to milli-Ångstrom precision. The experimentally determined bond length alternation is estimated to correspond to 60:40 contributions from the two resonant forms of 1H-indazole.

Atmospherically relevant acrolein-water complexes: spectroscopic evidence of aldehyde hydration and oxygen atom exchange.

Direct spectroscopic evidence of a reaction occurring between acrolein and water and involving the exchange of an oxygen atom has been obtained by characterizing the non-covalently bound water complexes and their isotopic forms, via rotational spectroscopy. The experimental geometries of the binary and ternary water complexes have been determined, and other stationary points on the reaction path have been characterized using ab initio quantum chemical methods at the MP2/6-311++G(d,p) level. These results can enhance the understanding of the water-mediated atmospherically important reactions involving acrolein.

A General Treatment to Study Molecular Complexes Stabilized by Hydrogen-, Halogen-, and Carbon-Bond Networks: Experiment and Theory of (CH2 F2 )n ⋠⋠⋠(H2 O)m.

Rotational spectra of several difluoromethane-water adducts have been observed using two broadband chirped-pulse Fourier-transform microwave (CP-FTMW) spectrometers. The experimental structures of (CH2 F2 )⋅⋅⋅(H2 O)2 , (CH2 F2 )2 ⋅⋅⋅(H2 O), (CH2 F2 )⋅⋅⋅(H2 O)3 , and (CH2 F2 )2 ⋅⋅⋅(H2 O)2 were unambiguously identified with the aid of 18 isotopic substituted species. A subtle competition between hydrogen, halogen, and carbon bonds is observed and a detailed analysis was performed on the complex network of non-covalent interactions which stabilize each cluster. The study shows that the combination of stabilizing contact networks is able to reinforce the interaction strength through a cooperative effect, which can lead to large stable oligomers.

Effects of Chlorination on the Tautomeric Equilibrium of 2-Hydroxypyridine: Experiment and Theory.

The effects of halogenation on the tautomeric and conformational equilibria of the model system 2-hydroxypyridine/2-pyridone have been investigated through chlorine substitution at positions 3, 4, 5, and 6. In the gas phase, the lactim syn-periplanar tautomer (OHs ) was the predominant species for all compounds over the lactam form (C=O) and the less abundant anti-periplanar lactim (OHa ). However, the population of the three species was shown to be dependent on the position of the chlorine substitution. Chlorination in position 5 or 6 strongly stabilizes the OHs tautomer, whereas the C=O form has a significant population when the ring is chlorinated in positions 3 or 4. Overall, the OHa form is the least favourable form, although the 3-substitution favours the population of this tautomer. In addition, the C=O tautomer is strongly stabilized in the solvent, which makes it the dominant form in some substituted species. This study has been performed by means of rotational spectroscopy in the gas phase and/or theoretical calculations in the isolated phase and in solution. Both the OHs and C=O forms of 5-chloro-2-hydroxypyridine and the OHs form of 6-chloro-2-hydroxypyridine were experimentally observed. All transitions displayed a complex nuclear hyperfine structure owing to the presence of the chlorine and nitrogen nuclei. For all species, a full quadrupolar hyperfine analysis has been performed. This has provided crucial information for the unambiguous identification of tautomers.

Regarding the torsional flexibility of the dihydrolipoic acid's pharmacophore: 1,3-propanedithiol.

The conformational space of antioxidant dihydrolipoic acid has been explored through the investigation of its pharmacophore, 1,3-propanedithiol. Five of the possible 25 non-equivalent isomers (namely: gGGg', gGGg, g'AGg, gAGg and g'AGg') were observed in the 59.6-74.4 GHz frequency region using free-jet absorption rotational spectroscopy. Furthermore, for three of them, the 34S mono-substituted isotopologues were also detected in natural abundance. Theoretical simulations show that the balance of steric and electronic intramolecular interactions arises on a shallow conformational potential energy surface and suggest that in polar solvents the flexibility of the dithiol chain is greater than that in the isolated phase.

Effects of Fluorine Substitution on the Microsolvation of Aromatic Azines: The Microwave Spectrum of 3-Fluoropyridine-Water.

The rotational spectra of the parent species and of four water isotopologues (-DOH, -HOD, -DOD, and -H(18)OH) of the adduct 3-fluoropyridine-water have been investigated using pulsed supersonic-jet Fourier-transform microwave spectroscopy. From the rotational constants, the structure of the adduct was deduced where the water is linked to the aromatic ring through an intermolecular O-H···N hydrogen bond with a bond distance of 1.9961(5) Å and an O-H···N angle of 156.8(1)°. The shape of the complex is such that the water oxygen is in the plane of the aromatic ring, on the opposite side of fluorine, and forming a C-H···N weak hydrogen bond with the adjacent aromatic hydrogen.

How CO2 Interacts with Carboxylic Acids: A Rotational Study of Formic Acid-CO2.

The rotational spectra of the 1:1 formic acid-carbon dioxide molecular complex and of its monodeuterated isotopologues are analysed in the 6.5-18.5 and 59.6-74.4 GHz frequency ranges using a pulsed jet Fourier transform microwave spectrometer and a free-jet absorption millimetre wave spectrometer, respectively. Precise values of the rotational and quartic centrifugal distortion constants are obtained from the measured frequencies, and quadrupole coupling constants are determined from the deuterium hyperfine splittings. Structural parameters are estimated from the moments of inertia and their differences among isotopologues: the complex has a planar structure with the two subunits held together by a HC(O)OH⋅⋅⋅O=C=O (2.075 Å) and a HC(OH)O⋅⋅⋅CO2 (2.877 Å) interactions. The ab initio intermolecular binding energy, obtained at the counterpoise corrected MP2/aug-cc-pVTZ level of calculation, is De =17 kJ mol(-1).

Millimeter Wave Spectrum of the Weakly Bound Complex CH2═CHCN·H2O: Structure, Dynamics, and Implications for Astronomical Search.

The weakly bound 1:1 complex between acrylonitrile (CH2═CHCN) and water has been characterized spectroscopically in the millimeter wave range (59.6-74.4 GHz) using a Free Jet Absorption Millimeter Wave spectrometer. Precise values of the rotational and quartic centrifugal distortion constants have been obtained from the measured frequencies of the normal and isotopically substituted water moiety (DOH, DOD, H(18)OH). Structural parameters have been estimated from the rotational constants and their differences among isotopologues: the complex has a planar structure with the two subunits held together by a O-H···N (2.331(3) Å) and a C-H···O (2.508(4) Å) interaction. The ab initio intermolecular binding energy, obtained at the counterpoise corrected MP2/aug-cc-pVTZ level of calculation, is De = 24.4 kJ mol(-1).

Accurate spectroscopy of polycyclic aromatic compounds: from the rotational spectrum of fluoren-9-one in the millimeter wave region to its infrared spectrum.

The rotational spectrum of fluoren-9-one, a small oxygenated polycyclic aromatic hydrocarbon, has been recorded and assigned in the 52-74.4 GHz region. The determined small negative value of the inertia defect (-0.3 u Å(2)) has been explained in terms of vibrational-rotational coupling constants calculated at the B3LYP/cc-pVTZ level of theory. Vibrational anharmonic analysis together with second-order vibrational perturbation theory approximation was applied both to fluorenone and its reduced form, fluorene, to predict the mid- and near-infrared spectra. The data presented here give precise indication on the fluorenone ground state structure, allow for an accurate spectral characterization in the millimeter wave and infrared regions, and hopefully will facilitate extensive radio astronomical searches with large radio telescopes.