Are there reliable DFT approaches for 13C NMR chemical shift predictions of fullerene C60 derivatives?

The relationships between experimental and theoretical ¹³C NMR chemical shifts of a pristine fullerene C₆₀, monoadducts from [2 + n] cycloaddition (n = 1–3), and one [2 + 1] bis‐adduct are systematically analyzed for the first time by using diverse quantum‐chemical levels of theory. These levels involved B3LYP, B3PW91, B97‐2, mPW1PW91, PBE1PBE, and X3LYP hybrid functionals combined with 3‐21G, 6‐31G, 6‐31G(d), 6‐31G(d,p), 6‐31G(d,2p), LanL2DZ, and SDDAll basis sets. X3LYP/6‐31G approach is determined to have the lowest deviations from the ¹³C NMR experimental data compared to the other methods for all the fullerene compounds (mean absolute error value is 0.856 ppm and root mean squared error value is 1.197 ppm). The highest deviations are characteristic for α (sp² C2/C5/C8/C10) and β (sp² C6/C7/C11/C12) carbon atoms relative to a functionalization site and for those (sp³ C1/C9) directly attached with a side fragment in the [2 + n] monoadducts (n = 1–3). A probable reason of such deviation is that the approaches do not take into account a contribution of paramagnetic ring currents to ¹³C NMR chemical shifts. The results will be useful in design of novel fullerene derivatives and in performing unambiguous ¹³C NMR chemical shift assignments with modern quantum chemistry calculations.     

Synthesis and study on 2,4,6,8-tetraaryl- 3,7-diazabicyclo[3.3.1]nonanes and their derivatives

1-Carbethoxy-2,4,6,8-tetraaryl-3,7-diazabicyclo[3.3.1] nonan-9-ones (1, 2) were synthesized and their 1H and 13C NMR data are reported. Chemical shifts and spectral assignments for 2,4,6,8-tetrakis(4-chlorophenyl)-3,7-diazabicyclo[3.3.1]nonan-9-one (3), 2,4,6,8-tetraphenyl-3-thia-7-azabicyclo[3.3.1]nonan-9-one (4) and 2,4,6,8-tetraaryl-3,7-diazabicyclo[3.3.1]nonanes (5-7) are also included.     

Prediction of the ~(13)C NMR Chemical Shifts of 9,10-Dihydrophenanthrene Analogues by the GIAO Method

After the geometry optimizations at the B3LYP/6-31+G(d,p) level, the NMR calculations of a series of 9,10-dihydrophenanthrene analogues have been carried out by GIAO method at the HF/6-31+G(d) level. The calculated 13C NMR chemical shifts are in agreement with the observed values. By a series of linear correlation equations (δpred = a + bδcal.c) of the 13C chemical shifts, accurate prediction of 13C chemical shifts was achieved for the new 9,10- dihydrophenanthrene compound, for which the predicted 13C NMR chemical shifts are in quite good agreement with the experimental values. The linear correlation between δpred and δexptl is excellent, and the square of correlation coefficient, r2, is up to 0.9973. The maximum absolute difference between δpred and δexptl, Δδ, is 4.5 ppm, and the rms error between δpred and δexptl is 2.55 ppm. In the meantime, according to the theoretical predicted result, we could confirm that the new 9,10-dihydrophenanthrene analogue is erianthridin (2,7-dihydroxy-3,4-dimethoxy-9,10-dihydro-phenanthrene).     

The impact of the pi-electron conjugation on (15)N, (13)C and (1)H NMR chemical shifts in push-pull benzothiazolium salts. Experimental and theoretical study

The (15)N as well as (13)C and (1)H chemical shifts of eight push-pull benzothiazolium iodides with various pi-conjugated chains between dimethylamino group and benzothiazolium moiety have been determined by NMR spectroscopy at the natural-abundance level of all nuclei in DMSO-d(6) solution. In general, the quaternary benzothiazolium nitrogen is more shielded [delta((15)N-3) vary between - 241.3 and - 201.9 ppm] with respect to parent 3-methylbenzothiazolium iodide [delta((15)N-3) = - 183.8 ppm], depending on the length and constitution of the pi-conjugated bridge. A larger variation in (15)N chemical shifts is observed on dimethylamino nitrogen, which covers the range of - 323.3 to - 257.2 ppm. The effect of pi-conjugation degree has a less pronounced influence on (13)C and (1)H chemical shifts. Experimental data are interpreted by means of density functional theory (DFT) calculations. Reasonable agreement between theoretical and experimental (15)N NMR chemical shifts was found, particularly when performing calculations with hybrid exchange-correlation functionals. A better accord with experiment is achieved by utilizing a polarizable continuum model (PCM) along with an explicit treatment of hydrogen-bonding between the solute and the water present in dimethylsulfoxide (DMSO). Finally, (13)C and (1)H NMR spectra were computed and analysed in order to compare them with available experimental data.     

Natural abundance nitrogen-15 nuclear magnetic resonance spectral studies on selected donors

The natural abundance 15N-NMR chemical shifts of selected aliphatic amines, 2-substituted pyridine type compounds, bialicyclic tertiary amines have been measured as a function of the nature of the solvent. In the case of cyclic aliphatic amines, like piperidine, morpholine, piperazine, thiomorpholine, the nitrogen is more shielded in concentrated solution compared to that in dilute solution whereas in the hydrogen bonding and protonating solvents there is a prominent deshielding. 2-Substituted pyridines studied can be further divided into four sub groups. The site of hydrogen bonding and protonation in 2-amino, 2-hydroxy and 2-mercapto pyridines have been conclusively proved from the 15N-NMR chemical shifts and the well-known tautomeric forms of the above compounds. Similarly in the case of 2-(2-thienyl)pyridine and 2-(3-thienyl)pyridine, the site of donation has been proved as the nitrogen of the pyridine ring in both the compounds. In a similar manner, the site of hydrogen bonding and protonation in two individual compounds 2-anilinopyridine and 2-(2-pyridyl)benzimidazole have also been established. Among the bialicyclic amines, 1,2-diazabicyclo[2.2.2]octane (DABCO) behaved differently from the other two compounds. In both 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU) and 1,5-diazabicyclo[4.3.0]non-5-ene (DBN), it was possible to show that N1-nitrogen in both the compounds is the site of donation. The effect of the second donor site on the 15N-NMR chemical shift, the site of donation in the selected compounds and some typical compounds reported in literature have been presented and discussed.     

Comparison of various density functional methods for distinguishing stereoisomers based on computed (1)H or (13)C NMR chemical shifts using diastereomeric penam beta-lactams as a test set

Full (1)H and (13)C NMR chemical shift assignments were made for two sets of penam beta-lactams: namely, the diastereomeric (2S, 5S, 6S)-, (2S, 5R, 6R)-, (2S, 5S, 6R)-, and (2S, 5R, 6S)-methyl 6-(1,3-dioxoisoindolin-2-yl)-3,3-dimethyl-7-oxo-4-thia-1-aza-bicyclo[3.2.0]heptane-2-carboxylates (1-4) and (2S, 5R, 6R)-, (2S, 5S, 6R)-, and (2S, 5R, 6S)-6-(1,3-dioxoisoindolin-2-yl)-3,3-dimethyl-7-oxo-4-thia-1-aza-bicyclo[3.2.0]heptane-2-carboxylic acids (6-8). Each penam was then modeled as a family of conformers obtained from Monte Carlo searches using the AMBER* force field followed by IEFPCM/B3LYP/6-31G(d) geometry optimization of each conformer using chloroform solvation. (1)H and (13)C chemical shifts for each conformer were computed at the WP04, WC04, B3LYP, and PBE1 density functional levels as Boltzmann averages of IEFPCM/B3LYP/6-311 + G(2d,p) energies over each family. Comparisons between experimental and theoretical chemical shift data were made using the total absolute error (|Deltadelta| (T)) criterion. For the (1)H shift data, all methods were sufficiently accurate to identify the proper stereoisomers. Computed (13)C shifts were not always successful in identifying the correct stereoisomer, regardless of which DFT method was used. The relative ability of each theoretical approach to discriminate among stereoisomers on the basis of proton shifts was also evaluated.     

Oxoperoxo vanadium(V) complexes of L-lactic acid: density functional theory study of structure and NMR chemical shifts

Various combinations of density functionals and pseudopotentials with associated valence basis-sets are compared for reproducing the known solid-state structure of [V 2O 2(OO) 2 l-lact 2] (2-) cis . Gas-phase optimizations at the B3LYP/SBKJC level have been found to provide a structure that is close to that seen in the solid state by X-ray diffraction. Although this may result in part from error compensation, this optimized structure allowed satisfactory reproduction of solution multinuclear NMR chemical shifts of the complex in all-electron DFT-IGLO calculations (UDFT-IGLO-PW91 level), suggesting that it is probably close to that found in solution. This combination of approaches has subsequently been used to optimize the structures of the vanadium oxoperoxo complexes [V 2O 3(OO) l-lact 2] (2-) cis , [V 2O 3(OO) l-lact 2] (2-) trans , and [VO(OO)( l-lact)(H 2O)] (-) cis . The (1)H, (13)C, (51)V, and (17)O NMR chemical shifts for these complexes have been calculated and compared with the experimental solution chemical shifts. Excellent agreement is seen with the (13)C chemical shifts, while somewhat inferior agreement is found for (1)H shifts. The (51)V and (17)O chemical shifts of the dioxo vanadium centers are well reproduced, with differences between theoretical and experimental shifts ranging from 22.9 to 35.6 ppm and from 25.1 to 43.7 ppm, respectively. Inferior agreement is found for oxoperoxo vanadium centers, with differences varying from 137.3 to 175.0 ppm for (51)V shifts and from 148.7 to 167.0 ppm for (17)O(oxo) shifts. The larger errors are likely to be due to overestimated peroxo O-O distances. The chosen methodology is able to predict and analyze a number of interesting structural features for vanadium(V) oxoperoxocomplexes of alpha-hydroxycarboxylic acids.     

Quantum Mechanical Calculation of 13C NMR Chemical Shifts in a Series of Isomeric Fucobiosides with the Account for Conformational Equilibrium

Chemical shift calculations for isomeric fucobiosides were performed by means of quantum mechanics, taking into account the conformational equilibrium. It was found during our previous investigations that generally two conformers are present for the studied compounds. As was expected, the major differences of the chemical shift values between these conformers were found for the atoms in the area near the glycosidic linkage. Taking into account the relative weight of each conformer, we succeeded in obtaining a good coincidence of calculated and experimental chemical shifts for the compounds for which experimental shifts were known. The deviation was not more than 2 ppm.     

1H, 13C and 15N NMR studies on adducts formation of rhodium(II) tetraacylates with some azoles in CDCl3 solution

Adduct formations of rhodium(II) tetraacetate and tetratrifluoroacetate with some 1H-imidazoles, oxazoles, thiazoles, 1H-pyrazoles and isoxazole have been investigated by the use of 1H, 13C, 15N NMR and electronic absorption spectroscopy (VIS) in the visible range. Azoles tend to form axial adducts containing rhodium(II) tetraacylates bonded via nitrogen atom. Bulky substituents close to the nitrogen atom prevent the Rh--N bond formation, and in several cases switch over the binding site to the oxygen or sulphur atoms. The (15)N adduct formation shift Deltadelta(15N) (Deltadelta = delta(adduct) - delta(ligand)) varied from ca - 40 to - 70 ppm for the nitrogen atom involved in complexation, and of a few parts per million only, from ca - 6 to 3 ppm, for the non-bonded nitrogen atom within the same molecule. The Deltadelta(1H) values do not exceed one ppm; Deltadelta(13C) ranges from - 1 to 6 ppm. Various complexation modes have been proved by electronic absorption spectroscopy in the visible region (VIS). For comparison purposes, some adducts of pyridine, thiophene and furan derivatives have been measured as well. The experimental findings were compared with calculated chemical shifts, obtained by means of DFT B3LYP method, using 6-311 + G(2d,p), 6-31(d)/LanL2DZ and 6-311G(d,p) basis set.     

Computational NMR of natural products: On the way to super large molecules exemplified with alasmontamine A

¹H and ¹³C nuclear magnetic resonance (NMR) chemical shifts of a tetrakis monoterpene indole alkaloid alasmontamine A with a molecular formula of C₈₄H₉₁N₈O₁₂ have been calculated at the PBE0/pcSseg-2//pcseg-2 level of theory on M06-2X/aug-cc-pVDZ geometry. In the course of the preliminary conformational search, six true minimum energy conformers were identified that can contribute to the actual conformation of this huge alkaloid. Calculated chemical shifts generally demonstrated a good agreement with available experimental data characterized with a corrected mean absolute error of 0.10 ppm for the range of about 7 ppm for protons and 1.1 ppm for the range of about 160 ppm for carbons.     

Thermodynamics and Spectroscopy of Halogen- and Hydrogen-Bonded Complexes of Haloforms with Aromatic and Aliphatic Amines

Similarities and differences of halogen and hydrogen bonding were explored via UV–Vis and ¹H NMR measurements, X-ray crystallography and computational analysis of the associations of CHX₃ (X=I, Br, Cl) with aromatic (tetramethyl-p-phenylenediamine) and aliphatic (4-diazabicyclo[2,2,2]octane) amines. When the polarization of haloforms was taken into account, the strengths of these complexes followed the same correlation with the electrostatic potentials on the surfaces of the interacting atoms. However, their spectral properties were quite distinct. While the halogen-bonded complexes showed new intense absorption bands in the UV–Vis spectra, the absorptions of their hydrogen-bonded analogues were close to the superposition of the absorption of reactants. Additionally, halogen bonding led to a shift in the NMR signal of haloform protons to lower ppm values compared with the individual haloforms, whereas hydrogen bonding of CHX₃ with aliphatic amines resulted in a shift in the opposite direction. The effects of hydrogen bonding with aromatic amines on the NMR spectra of haloforms were ambivalent. Titration of all CHX₃ with these nucleophiles produced consistent shifts in their protons’ signals to lower ppm values, whereas calculations of these pairs produced multiple hydrogen-bonded minima with similar structures and energies, but opposite directions of the NMR signals’ shifts. Experimental and computational data were used for the evaluation of formation constants of some halogen- and hydrogen-bonded complexes between haloforms and amines co-existing in solutions.     

Elucidating the Origin of Conformational Energy Differences in Substituted 1,3-Dioxanes: A Combined Theoretical and Experimental Study

13C NMR spectroscopy, ab initio quantum mechanics, and molecular mechanics have been used to investigate the trans-4-(trifluoromethyl)-2,2,6-trimethyl-1,3-dioxane chair/twist-boat equilibrium. The molecular mechanics calculations were based upon the MM3 and AMBER force fields. A 6-31G basis set was used for the ab initio calculations, and MP2 correlation corrections were applied. Both the ab initio and AMBER molecular mechanics calculations are consistent with the (13)C NMR chemical shift differences for the trans-4-(trifluoromethyl)-2,2,6-trimethyl-1,3-dioxane conformers. The predicted chair to twist-boat equilibrium suggested by the MM3 calculations is not consistent with the experimental data. These results support the suggestion by Howard et al. (Howard, A. E.; Cieplak, P.; Kollman, P. A. J. Comput.Chem. 1995, 16, 243-261) on the critical role of electrostatic interactions in determining the chair/twist-boat equilibrium.     

1,5] Sigmatropic hydrogen shifts in cyclic 1,3-dienes

Density functional calculations have been carried out for [1,5] hydrogen shifts in 1,3-cycloalkadienes (cyclohexadiene, cycloheptadiene, and cyclooctadiene). The complexity of the potential surfaces of these reactions was found to increase with ring size. For 1,3-cyclohexadiene a single transition structure for the [1,5] hydrogen shift was located, which connects the two enantiomeric conformers. For 1,3-cycloheptadiene two enantiomeric transition structures for the [1,5] hydrogen shift were located, which interconnect three conformers of the diene, a pair of enantiomeric conformers and a third achiral conformer. Finally for 1,3-cyclooctadiene two diastereomeric transition structures were found in addition to six conformers (three pairs of enantiomeric conformers) of the diene. Calculated activation energies for the [1,5] hydrogen shifts were found to be in qualitative agreement with experiment. Variation in these energies are attributed to strain energies present in either the diene or the transition structure.     

The prediction of (1)H chemical shifts in amines: a semiempirical and ab initio investigation

Twenty one conformationally fixed amines and their N,N-dimethyl derivatives were obtained commercially or synthesized. These included cis and trans 4-t-butyl cyclohexylamine, 2-exo and 2-endo norbornylamine, 2-adamantylamine, 4-phenylpiperidine, 1-napthylamine and tetrahydro-1-napthylamine. The (1)H NMR spectra of these amines were measured in CDCl(3) solution, assigned and the (1)H chemical shifts given. This data was used to investigate the effect of the amino group on the (1)H chemical shifts in these molecules. These effects were analyzed using the CHARGE model. This calculates the electric field and steric effects of the amino group for protons more than three bonds removed, together with functions for the calculation of two-bond and three-bond effects. The rotational isomerism about the C--N bond of the amino group was investigated by ab initio calculations of the potential energy surface (PES) about this bond at the HF/3-21G level. The resulting conformers were then minimized at the B3LYP/6-311 + + G (d,p) level. These geometries were then used to calculate the (1)H chemical shifts in the above compounds by CHARGE and the ab initio gauge-invariant atomic orbital (GIAO) method at the B3LYP/6-311 + + G(d,p) level and the shifts were compared with those observed. The compounds investigated gave 170 (1)H chemical shifts ranging from 0.60 to 8.2 ppm. The rms errors (obs.-calc.) were ca 0.1 ppm (CHARGE) and ca 0.2 ppm (GIAO). Large deviations of ca 1.0 ppm were observed for the NH protons in the GIAO calculations. The complex spectra of alkyl and aryl amines can thus be successfully predicted by both ab initio and semiempirical methods except for the NH protons, for which the ab initio calculations are not sufficiently accurate.     

Spherical shape complementarity as an overriding motif in the molecular recognition of noncharged organic guests by p-sulfonatocalix[4]arene: complexation of bicyclic azoalkanes

[reaction: see text] The complexation of p-sulfonatocalix[4]arene (CX4) with 2,3-diazabicyclo[2.2.1]hept-2-ene (1), 2,3-diazabicyclo[2.2.2]oct-2-ene (2), 2,3-diazabicyclo[3.2.2]non-2-ene (3), 1-methyl-4-isopropyl-2,3-diazabicyclo[2.2.2]oct-2-ene (4), and 1-phenyl-2,3-diazabicyclo[2.2.2]oct-2-ene (5) was studied in D2O at pD 7.4 by 1H NMR spectroscopy. The formation of deep inclusion complexes was indicated by large upfield 1H NMR shifts of the guest protons (up to 2.6 ppm), which were also used to assign, in combination with 2D ROESY spectra, a preferential inclusion of the isopropyl group of 4 and a dominant inclusion of the azo bicyclic residue for 5. The bicyclic azoalkanes 1-3 showed exceptionally high binding constants on the order of 1000 M(-1), 1-2 orders of magnitude larger than for previously investigated noncharged organic guest molecules. The strong binding was attributed to the spherical shape complementarity between the guest and the conical cavity offered by CX4. Interestingly, although the derivatives 4 and 5 are more hydrophobic, they showed a 2-3 times weaker binding, which was again attributed to the deviation from spherical shape in these bridgehead-substituted derivatives. The preferential inclusion of the hydrophilic but spherical bicyclic residue of 5 rather than the hydrophobic aromatic phenyl group provides a unique observation in aqueous host-guest chemistry and corroborates the pronounced spherical shape affinity of CX4.     

13-methyl-2,6-dithia[7]metacyclophane: a useful molecule to connect VT NMR results and structure with calculations

Ground state energies (DFT) and 1H and 13C NMR chemical shifts are calculated for the conformers of 13-methyl-2,6-dithia[7]metacyclophane (1), and the results are compared with X-ray structural data and variable-temperature NMR data, including the determination of the activation barrier. Calculations predict the correct low energy conformer with good agreement with chemical shifts, bond distances, and angles. VT NMR data for the 10-tert-butyl-substituted derivative 2 indicate that it undergoes the same conformational equilibria as 1. This paper should enhance the confidence that organic chemists have in calculations to satisfactorily predict conformer energies.     

Cation-pi interactions of a thiocarbonyl group and a carbonyl group with a pyridinium nucleus

Attractive interactions between a thiocarbonyl group and a pyridinium nucleus, and between a carbonyl group and a pyridinium nucleus have been proven by (1)H and (13)C NMR studies, UV-vis spectral analyses, and X-ray crystallographic analyses of nicotinic amides 1 and 3, and pyridinium salts 2 and 4. Comparison of the Deltadelta values, which are the differences in the chemical shifts with reference compounds 5 or 6, showed that the absolute Deltadelta values of 2 and 4 are much larger than those of 1 and 3. In the UV-vis spectra, the n-->pi absorption of the C=S group of 2a exhibited a significant blue shift in CHCl(3). X-ray crystallographic analysis of 1-4 clearly showed that the C=S group of 2a and the C=O group of 4 are very close to the pyridinium moiety compared to the case of 1 and 3. In addition, the X-ray crystal packing structure of 2a showed the C=S group is sandwiched between two pyridinium rings. These experimental results strongly suggested the existence of attractive (C=S)...Py(+) and (C=O)...Py(+) interactions in solution and in crystal. The optimized geometries of 1 and 2 calculated at the HF/6-311G level are in good agreement with their X-ray geometries. MP2/6-311G calculations for the model systems of pyridinium salts 2 and 4 predicted that the electrostatic and induction energies are the major source of the attractive interactions. Since the larger contribution of electrostatic and induction interactions are characteristic features of cation-pi interactions, the (C=S)...Py(+) and (C=O)...Py(+) interactions would be classified as a cation-pi interaction.     

Determination of absolute configuration using density functional theory calculations of optical rotation and electronic circular dichroism: chiral alkenes

In principle, the absolute configuration (AC) of a chiral molecule can be deduced from its optical rotation (OR) and/or its electronic circular dichroism (ECD). In practice, this requires reliable methodologies for predicting OR and ECD. The recent application of ab initio time-dependent density functional theory (TDDFT) to the calculation of transparent spectral region OR and ECD has greatly enhanced the reliability with which these phenomena can be predicted. TDDFT calculations of OR and ECD are being increasingly utilized in determining ACs. Nevertheless, such calculations are not perfect, and as a result, ACs determined are not 100% reliable. In this paper, we examine the reliability of the TDDFT methods in the case of chiral alkenes. Sodium d line specific rotations, [alpha]D, are predicted for 26 conformationally rigid alkenes of known AC, ranging in size from 5 to 20 C atoms, and with [alpha]D values in the range of 0-500. The mean absolute deviation of predicted [alpha]D values from experimental values is 28.7. With one exception, beta-pinene, the signs of [alpha]D are correctly predicted. Errors in calculated [alpha]D values are approximately random. Our results define a "zone of indeterminacy" within which calculated [alpha]D values cannot be used to determine ACs with >95% confidence. TDDFT ECD spectra are predicted for eight of the alkenes and compared to experimental spectra. Agreement ranges from modestly good to poor, leading to the conclusion that TDDFT calculations of ECD spectra are not yet of sufficient accuracy to routinely provide highly reliable ACs. TDDFT OR calculations for two conformationally flexible alkenes, 3-tert-butylcyclohexene and trans-4-carene, are also reported. For the former, predicted rotations are incorrect in sign over the range 589-365 nm. It is possible that the AC of this molecule has been incorrectly assigned.     

Quantifying weak hydrogen bonding in uracil and 4-cyano-4'-ethynylbiphenyl: a combined computational and experimental investigation of NMR chemical shifts in the solid state

Weak hydrogen bonding in uracil and 4-cyano-4'-ethynylbiphenyl, for which single-crystal diffraction structures reveal close CH...O=C and C[triple bond]CH...N[triple bond]C distances, is investigated in a study that combines the experimental determination of 1H, 13C, and 15N chemical shifts by magic-angle spinning (MAS) solid-state NMR with first-principles calculations using plane-wave basis sets. An optimized synthetic route, including the isolation and characterization of intermediates, to 4-cyano-4'-ethynylbiphenyl at natural abundance and with 13C[triple bond]13CH and 15N[triple bond]C labeling is described. The difference in chemical shifts calculated, on the one hand, for the full crystal structure and, on the other hand, for an isolated molecule depends on both intermolecular hydrogen bonding interactions and aromatic ring current effects. In this study, the two effects are separated computationally by, first, determining the difference in chemical shift between that calculated for a plane (uracil) or an isolated chain (4-cyano-4'-ethynylbiphenyl) and that calculated for an isolated molecule and by, second, calculating intraplane or intrachain nucleus-independent chemical shifts that quantify the ring current effects caused by neighboring molecules. For uracil, isolated molecule to plane changes in the 1H chemical shift of 2.0 and 2.2 ppm are determined for the CH protons involved in CH...O weak hydrogen bonding; this compares to changes of 5.1 and 5.4 ppm for the NH protons involved in conventional NH...O hydrogen bonding. A comparison of CH bond lengths for geometrically relaxed uracil molecules in the crystal structure and for geometrically relaxed isolated molecules reveals differences of no more than 0.002 A, which corresponds to changes in the calculated 1H chemical shifts of at most 0.1 ppm. For the C[triple bond]CH...N[triple bond]C weak hydrogen bonds in 4-cyano-4'-ethynylbiphenyl, the calculated molecule to chain changes are of similar magnitude but opposite sign for the donor 13C and acceptor 15N nuclei. In uracil and 4-cyano-4'-ethynylbiphenyl, the CH hydrogen-bonding donors are sp2 and sp hybridized, respectively; a comparison of the calculated changes in 1H chemical shift with those for the sp3 hybridized CH donors in maltose (Yates et al. J. Am. Chem. Soc. 2005, 127, 10216) reveals no marked dependence on hybridization for weak hydrogen-bonding strength.     

Conformational analysis of jet-cooled tryptophan analogs and Histamine using the MM3(94) force field: Comparison with experiment

MM3(94) has been used to predict the conformers of nitrogen-containing aromatic heterocycles with polar aliphatic sidechains. Computations were done for cases in which experimental gas-phase rotational constants have been determined and include histamine and analogs of tryptophan. The agreement with experiment for the tryptophan analogs is better than earlier MM2(87) computations but still not complete. A fairly good match can be made to experimental rotational constants of four histamine conformers, but other conformers are also predicted that may not be important experimentally. A comparison can be made with ab initio calculations undertaken for histamine. Similar structures were generally predicted, but there were significant discrepancies with MM3 in relative conformer energies. © 1995 by John Wiley & Sons, Inc.