(1)H NMR spectra of alcohols in hydrogen bonding solvents: DFT/GIAO calculations of chemical shifts

An initial investigation of accuracy required for the identification of small molecules in complex samples using quantum chemical calculated NMR chemical shifts

The majority of primary and secondary metabolites in nature have yet to be identified, representing a major challenge for metabolomics studies that currently require reference libraries from analyses of authentic compounds. Using currently available analytical methods, complete chemical characterization of metabolomes is infeasible for both technical and economic reasons. For example, unambiguous identification of metabolites is limited by the availability of authentic chemical standards, which, for the majority of molecules, do not exist. Computationally predicted or calculated data are a viable solution to expand the currently limited metabolite reference libraries, if such methods are shown to be sufficiently accurate. For example, determining nuclear magnetic resonance (NMR) spectroscopy spectra in silico has shown promise in the identification and delineation of metabolite structures. Many researchers have been taking advantage of density functional theory (DFT), a computationally inexpensive yet reputable method for the prediction of carbon and proton NMR spectra of metabolites. However, such methods are expected to have some error in predicted 13C and 1H NMR spectra with respect to experimentally measured values. This leads us to the question-what accuracy is required in predicted 13C and 1H NMR chemical shifts for confident metabolite identification? Using the set of 11,716 small molecules found in the Human Metabolome Database (HMDB), we simulated both experimental and theoretical NMR chemical shift databases. We investigated the level of accuracy required for identification of metabolites in simulated pure and impure samples by matching predicted chemical shifts to experimental data. We found 90% or more of molecules in simulated pure samples can be successfully identified when errors of 1H and 13C chemical shifts in water are below 0.6 and 7.1 ppm, respectively, and below 0.5 and 4.6 ppm in chloroform solvation, respectively. In simulated complex mixtures, as the complexity of the mixture increased, greater accuracy of the calculated chemical shifts was required, as expected. However, if the number of molecules in the mixture is known, e.g., when NMR is combined with MS and sample complexity is low, the likelihood of confident molecular identification increased by 90%.

Fluorophore spectroscopy in aqueous glycerol solution: the interactions of glycerol with the fluorophore

A common perception exists that glycerol provides an inert-like environment modifying viscosity and index of refraction by its various concentrations in aqueous solution. Said perception is herein challenged by investigating the effects of the glycerol environment on the spectroscopic properties of fluorescein, as a representative fluorophore, using steady-state and time-resolved techniques and computational chemistry. Results strongly suggest that the fluorescence quantum yield, measured fluorescence lifetime (FLT), natural lifetime and calculated fluorescence lifetime are all highly sensitive to the presence of glycerol. Glycerol was found to impact both the ground and first excited states of fluorescein, quenching and modifying both absorption and emission spectra, affecting the fundamental electrical dipoles of the ground and first excited singlet states, and lowering FLT and quantum yield. Furthermore, the Stern-Volmer, Lippert-Mataga, Perrin and Strickler-Berg relations indicate that glycerol acts upon fluorescein in aqueous solution as a quencher and alters the fluorescein geometry. Predictions made by computational chemistry impressively correspond to experimental results, both indicating changes in the properties of fluorescein at around 35% v/v aqueous glycerol, a clear indication that glycerol is not an innocent medium. This study proposes the Strickler-Berg relation as a means of detecting non-negligible effects of a hosting medium on its host fluorophore. These new insights on the molecular structures, the interactions between glycerol and its host fluorophore, and the effects of one on the other may be essential for understanding fundamental phenomena in chemistry and related fields.

3,4-Ethylenedioxythiophene (EDOT) End-Group Functionalized Poly-ε-caprolactone (PCL): Self-Assembly in Organic Solvents and Its Coincidentally Observed Peculiar Behavior in Thin Film and Protonated Media

End-group functionalization of homopolymers is a valuable way to produce high-fidelity nanostructured and functional soft materials when the structures obtained have the capacity for self-assembly (SA) encoded in their structural details. Herein, an end-functionalized PCL with a π-conjugated EDOT moiety, (EDOT-PCL), designed exclusively from hydrophobic domains, as a functional "hydrophobic amphiphile", was synthesized in the bulk ROP of ε-caprolactone. The experimental results obtained by spectroscopic methods, including NMR, UV-vis, and fluorescence, using DLS and by AFM, confirm that in solvents with extremely different polarities (chloroform and acetonitrile), EDOT-PCL presents an interaction- and structure-based bias, which is strong and selective enough to exert control over supramolecular packing, both in dispersions and in the film state. This leads to the diversity of SA structures, including spheroidal, straight, and helical rods, as well as orthorhombic single crystals, with solvent-dependent shapes and sizes, confirming that EDOT-PCL behaves as a "block-molecule". According to the results from AFM imaging, an unexpected transformation of micelle-type nanostructures into single 2D lamellar crystals, through breakout crystallization, took place by simple acetonitrile evaporation during the formation of the film on the mica support at room temperature. Moreover, EDOT-PCL's propensity for spontaneous oxidant-free oligomerization in acidic media was proposed as a presumptive answer for the unexpected appearance of blue color during its dissolution in CDCl3 at a high concentration. FT-IR, UV-vis, and fluorescence techniques were used to support this claim. Besides being intriguing and unforeseen, the experimental findings concerning EDOT-PCL have raised new and interesting questions that deserve to be addressed in future research.

Nonrotational Mechanism of Polarization in Alcohols

Chemical polarity governs various mechanical, chemical, and thermodynamic properties of dielectrics. Polar liquids have been amply studied, yet the basic mechanisms underpinning their dielectric properties remain not fully understood, as standard models following Debye's phenomenological approach do not account for quantum effects and cannot aptly reproduce the full dc-up-to-THz spectral range. Here, using the illustrative case of monohydric alcohols, we show that deep tunneling and the consequent intermolecular separation of excess protons and "proton-holes" in the polar liquids govern their static and dynamic dielectric properties on the same footing. We performed systematic ultrabroadband (0-10 THz) spectroscopy experiments with monohydric alcohols of different (0.4-1.6 nm) molecular lengths and show that the finite lifetime of molecular species and the proton-hole correlation length are the two principle parameters responsible for the dielectric response of all the studied alcohols across the entire frequency range. Our results demonstrate that a quantum nonrotational intermolecular mechanism drives the polarization in alcohols while the rotational mechanism of molecular polarization plays a secondary role, manifesting itself in the sub-terahertz region only.

Clusters of Hydroxyl-Functionalized Cations Stabilized by Cooperative Hydrogen Bonds: The Role of Polarizability and Alkyl Chain Length

We explore quantum chemical calculations for studying clusters of hydroxyl-functionalized cations kinetically stabilized by hydrogen bonding despite strongly repulsive electrostatic forces. In a comprehensive study, we calculate clusters of ammonium, piperidinium, pyrrolidinium, imidazolium, pyridinium, and imidazolium cations, which are prominent constituents of ionic liquids. All cations are decorated with hydroxy-alkyl chains allowing H-bond formation between ions of like charge. The cluster topologies comprise linear and cyclic clusters up to the size of hexamers. The ring structures exhibit cooperative hydrogen bonds opposing the repulsive Coulomb forces and leading to kinetic stability of the clusters. We discuss the importance of hydrogen bonding and dispersion forces for the stability of the differently sized clusters. We find the largest clusters when hydrogen bonding is maximized in cyclic topologies and dispersion interaction is properly taken into account. The kinetic stability of the clusters with short-chained cations is studied for the different types of cations ranging from hard to polarizable or exhibiting additional functional groups such as the acidic C(2)-H position in the imidazolium-based cation. Increasing the alkyl chain length, the cation effect diminishes and the kinetic stability is exclusively governed by the alkyl chain tether increasing the distance between the positively charged rings of the cations. With adding the counterion tetrafluoroborate (BF4-) to the cationic clusters, the binding energies immediately switch from strongly positive to strongly negative. In the neutral clusters, the OH functional groups of the cations can interact either with other cations or with the anions. The hexamer cluster with the cyclic H-bond motive and "released" anions is almost as stable as the hexamer built by H-bonded ion pairs exclusively, which is in accord with recent IR spectra of similar ionic liquids detecting both types of hydrogen bonding. For the cationic and neutral clusters, we discuss geometric and spectroscopic properties as sensitive probes of opposite- and like-charge interaction. Finally, we show that NMR proton chemical shifts and deuteron quadrupole coupling constants can be related to each other, allowing to predict properties which are not easily accessible by experiment.

Monodisperse oligo(δ-valerolactones) and oligo(ε-caprolactones) with docosyl (C22) end-groups

Two different families of monodisperse oligoesters with α-hydroxyl-ω-docosyl (C₂₂) terminal groups [oligo(δ-valerolactone) and oligo(ϵ-caprolactone)] were isolated by flash column chromatography (FCC).

Ascorbigen A-NMR identification

The connectivities of all atoms in ascorbigen A, an important metabolite, were determined unambiguously for the first time. The connectivity between carbon atoms was established by 2D INADEQUATE, and one-bond ¹³ C-¹³ C coupling constants were determined for all pairs of directly connected carbon atoms except for two strongly coupled carbon pairs. The ¹³ C-¹³ C coupling in one of the pairs was proved by a modification of standard INADEQUATE; however, the signals from the other pair were too weak to be observed. The connectivity within the two strongly coupled C-C pairs was confirmed by a combination of COSY and gHSQC; the latter experiment also identified all C-H bonds. The proton nuclear magnetic resonance (¹ H NMR) spectra in dry dimethyl sulfoxide allowed identification and assignment of the signals due to NH and OH protons. The derived structure, 3-((1H-indol-3-yl)methyl)-3,3a,6-trihydroxytetrahydrofuro[3,2-b]furan-2(5H)-one, agrees with the structure suggested for ascorbigen A in 1966. The density functional theory (DFT) calculations showed that among 16 possible stereoisomers, only two complied with the almost zero value of the measured ³ J(H6-H6a). Of the two stereoisomers, 3S,3aS,6S,6aR and 3R,3aR,6R,6aS, the latter was excluded on synthetic grounds. The nuclear Overhauser effect measurements unveiled close proximity between H2' proton of the indole and the H6a proton of the tetrahydrofuro[3,2-b]furan part. Detailed structural interpretation of the measured NMR parameters by means of DFT NMR was hampered by rotational flexibility of the indole and tetrahydrofuro[3,2-b]furan parts and inadequacy of Polarizable Continuum Model (PCM) solvent model.

Cost-effective density functional theory (DFT) calculations of equilibrium isotopic fractionation in large organic molecules

The application of stable isotopes to address a wide range of biochemical, microbiological and environmental problems is hindered by the experimental difficulty and the computational cost of determining equilibrium isotopic fractionations (EIF) of large organic molecules. Here, we evaluate the factors that impact the accuracy of computed EIFs and develop a framework for cost-effective and accurate computation of EIFs by density functional theory (DFT). We generated two benchmark databases of experimentally determined EIFs, one for H isotopes and another for the isotopes of the heavy atoms C, N and O. The accuracy of several DFT exchange-correlation functionals in calculating EIFs was then evaluated by comparing the computational results to these experimental datasets. We find that with the def2-TZVP basis set, O3LYP had the lowest mean absolute deviation (21‰ and 3.9‰ for the isotopic fractionation of H and the heavier atoms, respectively), but the GGA/meta-GGA functionals τHCTH_D3BJ, τHCTH and HCTH have similar performances (22‰ and 4.1‰, respectively, for τHCTH_D3BJ). Leveraging the good performance of computationally efficient functionals, we provide a robust, practical, experimentally validated framework for using DFT to accurately predict EIFs of large organic molecules, including uncertainty estimates.

Solvent-Dependent Structures of Natural Products Based on the Combined Use of DFT Calculations and 1H-NMR Chemical Shifts

Detailed solvent and temperature effects on the experimental 1H-NMR chemical shifts of the natural products chrysophanol (1), emodin (2), and physcion (3) are reported for the investigation of hydrogen bonding, solvation and conformation effects in solution. Very small chemical shift of │Δδ│ < 0.3 ppm and temperature coefficients │Δδ/ΔΤ│ ≤ 2.1 ppb/K were observed in DMSO-d6, acetone-d6 and CDCl3 for the C(1)-OH and C(8)-OH groups which demonstrate that they are involved in a strong intramolecular hydrogen bond. On the contrary, large chemical shift differences of 5.23 ppm at 298 K and Δδ/ΔΤ values in the range of -5.3 to -19.1 ppb/K between DMSO-d6 and CDCl3 were observed for the C(3)-OH group which demonstrate that the solvation state of the hydroxyl proton is a key factor in determining the value of the chemical shift. DFT calculated 1H-NMR chemical shifts, using various functionals and basis sets, the conductor-like polarizable continuum model, and discrete solute-solvent hydrogen bond interactions, were found to be in very good agreement with the experimental 1H-NMR chemical shifts even with computationally less demanding level of theory. The 1H-NMR chemical shifts of the OH groups which participate in intramolecular hydrogen bond are dependent on the conformational state of substituents and, thus, can be used as molecular sensors in conformational analysis. When the X-ray structures of chrysophanol (1), emodin (2), and physcion (3) were used as input geometries, the DFT-calculated 1H-NMR chemical shifts were shown to strongly deviate from the experimental chemical shifts and no functional dependence could be obtained. Comparison of the most important intramolecular data of the DFT calculated and the X-ray structures demonstrate significant differences for distances involving hydrogen atoms, most notably the intramolecular hydrogen bond O-H and C-H bond lengths which deviate by 0.152 tο 0.132 Å and 0.133 to 0.100 Å, respectively, in the two structural methods. Further differences were observed in the conformation of -OH, -CH3, and -OCH3 substituents.

An automated framework for NMR chemical shift calculations of small organic molecules

When using nuclear magnetic resonance (NMR) to assist in chemical identification in complex samples, researchers commonly rely on databases for chemical shift spectra. However, authentic standards are typically depended upon to build libraries experimentally. Considering complex biological samples, such as blood and soil, the entirety of NMR spectra required for all possible compounds would be infeasible to ascertain due to limitations of available standards and experimental processing time. As an alternative, we introduce the in silico Chemical Library Engine (ISiCLE) NMR chemical shift module to accurately and automatically calculate NMR chemical shifts of small organic molecules through use of quantum chemical calculations. ISiCLE performs density functional theory (DFT)-based calculations for predicting chemical properties—specifically NMR chemical shifts in this manuscript—via the open source, high-performance computational chemistry software, NWChem. ISiCLE calculates the NMR chemical shifts of sets of molecules using any available combination of DFT method, solvent, and NMR-active nuclei, using both user-selected reference compounds and/or linear regression methods. Calculated NMR chemical shifts are provided to the user for each molecule, along with comparisons with respect to a number of metrics commonly used in the literature. Here, we demonstrate ISiCLE using a set of 312 molecules, ranging in size up to 90 carbon atoms. For each, calculation of NMR chemical shifts have been performed with 8 different levels of DFT theory, and with solvation effects using the implicit solvent Conductor-like Screening Model. The DFT method dependence of the calculated chemical shifts have been systematically investigated through benchmarking and subsequently compared to experimental data available in the literature. Furthermore, ISiCLE has been applied to a set of 80 methylcyclohexane conformers, combined via Boltzmann weighting and compared to experimental values. We demonstrate that our protocol shows promise in the automation of chemical shift calculations and, ultimately, the expansion of chemical shift libraries.

Evaluation of the Factors Impacting the Accuracy of 13C NMR Chemical Shift Predictions using Density Functional Theory-The Advantage of Long-Range Corrected Functionals

The various factors influencing the accuracy of ¹³C NMR calculations using density functional theory (DFT), including the basis set, exchange-correlation (XC) functional, and isotropic shielding calculation method, are evaluated. A wide selection of XC functionals (over 70) were considered, and it was found that long-range corrected functionals offer a significant improvement over the other classes of functionals. Based on a thorough study, it is recommended that for calculating NMR chemical shifts (δ) one should use the CSGT method, the COSMO solvation model, and the LC-TPSSTPSS exchange-correlation functional in conjunction with the cc-pVTZ basis set. A selection of problems in natural product identification are considered in light of the newly recommended level of theory.

Association of symmetrical alkane diols with pyridine: DFT/GIAO calculation of 1 H NMR chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of the free diol and of its 1 : 1 and 1 : 2 hydrogen-bonded complexes with pyridine have been computed for five symmetrical alkane diols on the basis of density functional theory, by applying the gauge-including atomic orbital method to geometry-optimized conformers. For certain conformers, intramolecular OH···OH interactions, evidenced by high NMR OH proton shifts, are further enhanced on going from the free diol to the corresponding 1 : 1 diol/pyridine complex. This is confirmed by atoms-in-molecules and non-covalent interaction plots. The computed OH and CH proton shifts for the diol and the two complexes correlate well with values obtained by analysing data from the NMR titration of the diols in benzene against pyridine. Shift values for the diols in neat pyridine are calculated by weighting the shifts of the various protons in the three forms (free diol, 1 : 1 and 1 : 2 diol/pyridine complexes) according to the experimentally determined association constants. The results are in good agreement with those observed, and after empirical scaling, the root mean square difference is 0.18 ppm. Copyright © 2016 John Wiley & Sons, Ltd.