C-H...O hydrogen bonding in 4-phenyl-benzaldehyde: a comprehensive crystallographic, spectroscopic and computational study

The crystal structure of 4-phenyl-benzaldehyde reveals the presence of a dimer linked by the C=O and C9-H groups of adjacent molecules. In the liquid phase, the presence of C-H...O bonded forms is revealed by both vibrational and NMR spectroscopy. A DeltaH value of -8.2 +/- 0.5 kJ mol(-1) for the dimerisation equilibrium is established from the temperature-dependent intensities of the bands assigned to the carbonyl-stretching modes. The NMR data suggest the preferential engagement of the C(2,6)-H and C(10/12)/C(11)-H groups as hydrogen bond donors, instead of the C(9)-H group. While ab initio calculations for the isolated dimers are unable to corroborate these NMR results, the radial distribution functions obtained from molecular dynamics simulations show a preference for C(2,6)-H and C(10/12)/C(11)-H...O contacts relative to the C(9)-H...O ones.     

All-atom Molecular Dynamic Simulations and NMR Spectra Study on Intermolecular Interactionsof N,N-dimethylacetamide-Water System

N,N-dimethylacetamide (DMA) has been investigated extensively in studying models of peptide bonds. An all-atom MD simulation and the NMR spectra were performed to investigate the interactions in the DMA-water system. The radial distribution functions (RDFs) and the hydrogen-bonding network were used in MD simulations. There are strong hydrogen bonds and weak C-H￠ ￠ ￠O contacts in the mixtures, as shown by the analysis of the RDFs. The insight structures in the DMA-water mixtures can be classified into different regions by the analysis of the hydrogen-bonding network. Chemical shifts of the hydrogen atom of water molecule with concentration and temperatures are adopted to study the interactions in the mixtures. The results of NMR spectra show good agreement with the statistical results of hydrogen bonds in MD simulations.     

Hydrogen bond stabilization in 1,3-dimethylimidazolium methyl sulfate and 1-butyl-3-methylimidazolium hexafluorophosphate probed by high pressure: the role of charge-enhanced C-H...O interactions in the room-temperature ionic liquid

The hydrogen bonding structures of room-temperature ionic liquids 1,3-dimethylimidazolium methyl sulfate and 1-butyl-3-methylimidazolium hexafluorophosphate have been studied by infrared spectroscopy. High-pressure infrared spectral profiles and theoretical calculations allow us to make a vibrational assignment of these compounds. The imidazolium C-H bands of 1,3-dimethylimidazolium methyl sulfate display anomalous non-monotonic pressure-induced frequency shifts. This discontinuity in frequency shift is related to enhanced C-H...O hydrogen bonding. This behavior is in contrast with the trend of blue shifts in frequency for the methyl C-H stretching mode at ca. 2960 cm(-1). Our results indicated that the imidazolium C-H groups are more favorable sites for hydrogen bonding than the methyl C-H groups in the pure 1,3-dimethylimidazolium methyl sulfate. Nevertheless, both methyl C-H and imidazolium C-H groups are favorable sites for C-H...O hydrogen bonding in a dilute 1,3-dimethylimidazolium methyl sulfate/D(2)O mixture. Hydrogen bond-like C-H...F interactions were observed between PF(6)(-) and H atoms on the alkyl side chains and imidazolium ring for 1-butyl-3-methylimidazolium hexafluorophosphate.     

Improper hydrogen C-H...O bonds cause self-association of alpha, beta-enaminoketones containing fluorosubstituted alkyl groups.

Concentration- and temperature-dependent IR, NMR and dipole-moment studies on 4-N,N-dimethylamino-1,1,1-trifluoro-3-buten-2-one and two of its higher homologues showed that these compounds undergo reversible dimerization in nonpolar solvents. Antiparallel "closed" dimers are formed with a network of improper intermolecular C-H...O hydrogen bonds. Quantitative analysis of the 1H NMR data yielded delta H0 = -17.6 kJ mol-1 and delta S0 = -46.9 J deg-1 mol-1. The interactions observed are the strongest among those involving a C-H group reported so far. The complex described here is the first example of a cyclic complex stabilized by two improper C-H...O hydrogen bonds. The conclusions drawn from the solution and solid-state data were confirmed by ab initio calculations.     

All‐atom Molecular Dynamic Simulations Combined with the Chemical Shifts Study on the Weak Interactions of Ethanol‐water System

All‐atom molecular dynamics (MD) simulation combined with chemical shifts was performed to investigate the interactions over the entire concentration range of the ethanol (EtOH)‐water system. The results of the simulation were adopted to explain the NMR experiments by hydrogen bonding analysis. The strong hydrogen bonds and weak C–H···O contacts coexist in the mixtures through the analysis of the radial distribution functions. And the liquid structures in the whole concentration of EtOH‐water mixtures can be classified into three regions by the statistic analysis of the hydrogen‐bonding network in the MD simulations. Moreover, the chemical shifts of the hydrogen atom are in agreement with the statistical results of the average number hydrogen bonds in the MD simulations. Interestingly, the excess relative extent of η_rel^E calculated by the MD simulations and chemical shifts in the EtOH aqueous solutions shows the largest deviation at x_EtOH≈0.18. The excess properties present good agreement with the excess enthalpy in the concentration dependence.     

Neutron and X-ray diffraction and spectroscopic investigations of intramolecular [C-H...F-C] contacts in post-metallocene polyolefin catalysts: modeling weak attractive polymer-ligand interactions

A family of Group 4 post-metallocene catalysts, supported by fluorine-functionalized tridentate ligands with the fluorine substituent in the locality of the metal center, is described. For the first time, the contentious C-H...F-C interaction has been characterized by a neutron diffraction study, which has allowed the position of the hydrogen atoms to be accurately determined. The nature of the weak intramolecular C-H...F-C contacts in these complexes in solution and the solid state was probed by using multinuclear NMR spectroscopy in tandem with neutron and X-ray crystallography. Evidence is presented to demonstrate that the spectroscopic C-H...F-C coupling occurs "through-space" rather than "through-bond" or by MF coordination. The titanium catalysts exhibit excellent activities and high co-monomer incorporation in olefin polymerization. The observed intramolecular C-H...F-C interactions are important with regards to potential applications in polyolefin catalysis because they substantiate the proposed ortho-F...H(beta) ligand-(polymer chain) contacts derived from DFT calculations for the remarkable fluorinated phenoxyimine Group 4 catalysts. Compared with agostic and co-catalyst...metal contacts, weak attractive noncovalent interactions between a polymer chain and a judiciously designed "active" ligand is a new concept in polyolefin catalysis.     

Influence of N-H...O and C-H...O hydrogen bonds on the 17O NMR tensors in crystalline uracil: computational study

We report a computational study for the 17O NMR tensors (electric field gradient and chemical shielding tensors) in crystalline uracil. We found that N-H...O and C-H...O hydrogen bonds around the uracil molecule in the crystal lattice have quite different influences on the 17O NMR tensors for the two C=O groups. The computed 17O NMR tensors on O4, which is involved in two strong N-H...O hydrogen bonds, show remarkable sensitivity toward the choice of cluster model, whereas the 17O NMR tensors on O2, which is involved in two weak C-H...O hydrogen bonds, show much smaller improvement when the cluster model includes the C-H...O hydrogen bonds. Our results demonstrate that it is important to have accurate hydrogen atom positions in the molecular models used for 17O NMR tensor calculations. In the absence of low-temperature neutron diffraction data, an effective way to generate reliable hydrogen atom positions in the molecular cluster model is to employ partial geometry optimization for hydrogen atom positions using a cluster model that includes all neighboring hydrogen-bonded molecules. Using an optimized seven-molecule model (a total of 84 atoms), we were able to reproduce the experimental 17O NMR tensors to a reasonably good degree of accuracy. However, we also found that the accuracy for the calculated 17O NMR tensors at O2 is not as good as that found for the corresponding tensors at O4. In particular, at the B3LYP/6-311++G(d,p) level of theory, the individual 17O chemical shielding tensor components differ by less than 10 and 30 ppm from the experimental values for O4 and O2, respectively. For the 17O quadrupole coupling constant, the calculated values differ by 0.30 and 0.87 MHz from the experimental values for O4 and O2, respectively.     

All-atom Molecular Dynamics Simulationsand NMR Spectroscopy Study on Interactions and Structures in N-Glycylglycine Aqueous Solution

All-atom molecular dynamics (MD) simulation and the NMR spectra are used to investi-gate the interactions in N-glycylglycine aqueous solution. Different types of atoms exhibit different capability in forming hydrogen bonds by the radial distribution function analysis. Some typical dominant aggregates are found in different types of hydrogen bonds by the statistical hydrogen-bonding network. Moreover, temperature-dependent NMR are used to compare with the results of the MD simulations. The chemical shifts of the three hydrogen atoms all decrease with the temperature increasing which reveals that the hydrogen bonds are dominant in the glycylglycine aqueous solution. And the NMR results show agreement with the MD simulations. All-atom MD simulations and NMR spectra are successful in revealing the structures and interactions in the N-glycylglycine-water mixtures.     

Structure of liquid formic acid investigated by first principle and classical molecular dynamics simulations

The structure of liquid formic acid has been investigated by Car-Parrinello and classical molecular dynamics simulations, focusing on the characterization of the H-bond network and on the mutual arrangement of pairs of bonded molecules. In agreement with previous computational studies, two levels of H-bonded structures have been found. Small clusters, characterized by O-H...O bonds, are held together by weak C-H...O bonds to form large branched structures. From the ab initio simulation we infer the importance of cyclic H-bond dimer configurations, typical of the gas phase. Most of these dimer structures are however found to be embedded into H-bonded chains. When only O-H...O bonds are taken into account, linear H-bond chains are detected as basic structures of the liquid. More branched structures occur when C-H...O bonds are also considered. Regarding the arrangement of molecular pairs, we observed that O-H...O bonds favor the occurrence of configurations with parallel molecular planes, whereas no preferential orientation is observed for molecules forming C-H...O bonds.     

Methyl rotation barriers in proteins from 2H relaxation data. Implications for protein structure

Side-chain 2H and backbone 15N relaxation data have been collected at multiple temperatures in the samples of the SH3 domain from alpha-spectrin. Combined analyses of the data allowed for determination of the temperature-dependent correlation times tauf characterizing fast methyl motion. Molecular dynamics simulations confirmed that tauf are dominated by methyl rotation; the corresponding activation energies approximate methyl rotation barriers. For 33 methyl groups in the alpha-spectrin SH3 domain the average barrier height was thus determined to be 2.8 +/- 0.9 kcal/mol. This value is deemed representative of the "fluid" hydrophobic protein core where some barriers are increased and others are lowered because of the contacts with surrounding atoms, but there is no local order that could produce systematically higher (lower) barriers. For comparison, the MD simulation predicts the average barrier of 3.1 kcal/mol (calculated via the potential of mean force) or 3.4-3.5 kcal/mol (rigid barriers after appropriate averaging over multiple MD snapshots). The latter result prompted us to investigate rigid methyl rotation barriers in a series of NMR structures from the Protein Databank. In most cases the barriers proved to be higher than expected, 4-6 kcal/mol. To a certain degree, this is caused by tight packing of the side chains in the NMR structures and stems from the structure calculation procedure where the coordinates are first annealed toward the temperature of 0 K and then subjected to energy minimization. In several cases the barriers >10 kcal/mol are indicative of van der Waals violations. The notable exceptions are (i) the structures solved using the GROMOS force field where tight methyl packing is avoided (3.0-3.6 kcal/mol) and (ii) the structure solved by means of the dynamic ensemble refinement method (Lindorff-Larsen, K.; Best, R. B.; DePristo, M. A.; Dobson, C. M.; Vendruscolo, M. Nature 2005, 433, 128) (3.5 kcal/mol). These results demonstrate that methyl rotation barriers, derived from the experiments that are traditionally associated with studies of protein dynamics, can be also used in the context of structural work. This is particularly interesting in view of the recent efforts to incorporate dynamics data in the process of protein structure elucidation.     

Structures and Intermolecular Interactions in Dimethyl Sulfoxide-Water System Studied by All-atom Molecular Simulations

An all-atom dimethyl sulfoxide (DMSO) and water model have been used for molecular dy-namics simulation. The NMR and IR spectra are also performed to study the structures and interactions in the DMSO-water system. And there are traditional strong hydrogen bondsand weak C—H…O contacts existing in the mixtures according to the analysis of the radial distribution functions. The insight structures in the DMSO-water mixtures can be classified into different regions by the analysis of the hydrogen-bonding network. Interestingly, the molar fraction of DMSO 0.35 is found to be a special concentration by the network. It is the transitional region which is from the water rich region to the DMSO rich region. The sta-ble aggregates of (DMSO)_m·S=O…HW—OW·(H₂O)_n might play a key role in this region.Moreover, the simulation is compared with the chemical shifts in NMR and wavenumbers in IR with concentration dependence. And the statistical results of the average number hydrogen bonds in the MD simulations are in agreement with the experiment data in NMR and IR spectra.     

The unexpected and large enhancement of the dipole moment in the 3,4-bis(dimethylamino)-3-cyclobutene-1,2-dione (DMACB) molecule upon crystallization: a new role of the intermolecular CH...O interactions.

The molecular dipole moment of the 3,4-bis(dimethylamino)-3-cyclobutene-1,2-dione (DMACB) molecule and its enhancement in the crystal was evaluated by periodic RHF ab initio computations. A discrete boundary partitioning of the electronic density that allows an unambiguous partitioning of the molecular space in the condensed phase was adopted. The resulting molecular dipole in the crystal compares favorably with the experimental value obtained by a multipolar analysis of single-crystal X-ray diffraction data recorded at 20 K, using a fuzzy boundary partitioning of the derived pseudoatom densities. We show that a large and highly significant molecular dipole enhancement may occur upon crystallization, despite the lack of a strongly hydrogen bonded environment in the crystal. The 23 unique C-H...O interactions which are formed upon packing of the DMACB molecule induce an increase in the molecular dipole (over 75%) that is comparable to or greater than that found in systems which are characterized by the stronger O-H...O and N-H...O hydrogen bonds. The DMACB molecule constitutes an excellent system for the study of C-H...O interactions in the condensed phase, since no other kind of competing hydrogen bonds is present in its crystal. A simple and qualitative model for the matrix contribution to the DMACB molecular dipole enhancement in the crystal is proposed. The formation of several weak C-H...O bonds is found to yield a small (about 0.2 e) net flux of electronic charge flowing from the hydrogens of the methyl groups to the carbonyl oxygen atoms. Despite the limited increase of the intramolecular charge transfer upon crystallization, a large molecular dipole enhancement occurs because the centroids of the positive and negative induced charges are quite far apart. This work highlights a new and important role of the C-H...O bond, besides those already known in the literature.     

Direct observation of Calpha-Halpha...O=C hydrogen bonds in proteins by interresidue h3JCalphaC' scalar couplings

The role of C-H...O hydrogen bonds in the stabilization of biomolecules is increasingly being recognized from the evidence of close C-H...O contacts in crystal structures. However, relatively little is known about their strength. Here, we report the observation of NMR scalar couplings (h3JCalphaC') between the two carbons on each side of Calpha-Halpha...O=C H-bonds in proteins. These couplings give direct evidence of the correlation of the electronic wave functions in the donor and acceptor groups of Calpha-Halpha...O=C H-bonds. A long-range H(NCO)CA experiment or a selective long-range H(NCA)CO experiment was used for the detection of h3JCalphaC' correlations in the beta-sheet regions of the immunoglobulin binding domain of protein G. In total, six such correlations were detectable. These correspond to half of the Calpha-Halpha...O=C H-bonds of protein G with Halpha...O distances shorter than 2.5 A. The h3JCalphaC' couplings range from 0.2 to 0.3 Hz and are in good agreement with predicted average values based on DFT/FPT calculations. An anticorrelation is observed with the size of h3JNC' coupling constants across N-HN...O=C H-bonds, which share the same acceptor carbonyl oxygen.     

The molecular and crystal structure of tert-butyl Nalpha-tert-butoxycarbonyl-L-(S-trityl)cysteinate and the conformation-stabilizing function of weak intermolecular bonding

The title compound, C31H37NO4S [systematic name: (R)-tert-butyl-2-[(tert-butoxycarbonyl)amino]-3-(tritylsulfanyl)propanoate] is an L-cysteine derivative with three functions: NH2, COOH and SH, blocked by protecting groups tert-butoxycarbonyl, tert-butyl and trityl, respectively. The main chain of the molecule adopts the extended, nearly all-trans C5 conformation with the intramolecular N-H...O=C hydrogen bond. The urethane group is not involved in any intermolecular hydrogen bonding. Only weak intermolecular hydrogen bonds and hydrophobic contacts are observed in the crystal structure. These are C-H...O hydrogen bonds and CH/pi interactions with donor...acceptor distances, C...O ca. 3.5 A and C...C ca. 3.7 A, respectively. The first type of interaction links phenyl H-atoms and carbonyl groups. The second type of interaction is formed between a methyl group of the tert-butyl fragment and a trityl phenyl ring. The resulting molecular conformation in the crystal is very close to an ab initio minimum energy conformer of the isolated molecule. The extended C5 conformation of the main peptide chain is the same and there is slight discrepancy in the disposition of trityl phenyl rings. Their small dislocation creates the possibility of forming the entire network above of extensive, specific, weak intermolecular interactions; these constrain the molecule and permit it to retain the minimum energy C5 conformation of its main chain in the solid state. In contrast, in n-hexane solution, where such specific interactions cannot occur, only a small population of the molecules adopts the extended C5 conformation.     

The magnitude of [C-H...O] hydrogen bonding in molecular and supramolecular assemblies

Ab initio calculations at the MP2/6-311++G** level on model systems (N-methylpyridinium complexes of dimethyl ether and dimethyl phosphate anion) provide quantitative measures of the large stabilization energies that arise from [C-H...O] contacts in charged systems. These attractive interactions control (i) the self-assembly of bipyridinium-based catenanes and rotaxanes in solution, (ii) the self-organization of left-handed Z-DNA with alternating [dC-dG] sequences in the solid state, and (iii) the binding of pyridinium derivatives with single- and double-stranded DNA. Slightly attractive interactions occur between the donor ether and phosphate moieties and a neutral pyridine molecule in the gas phase. Electrostatic potential and solvation calculations demonstrate that [C-H...O] interactions which involve a cationic [C-H] donor are dominated by electrostatic terms.     

Crystal structure of the clathrate of 5,12-dihydro-5,12-o-benzenonaphthacene-1,4-dione with benzene (2:1).

In the crystal of the title compound, the quinone fragment of the host (1) participates in weak D-A interactions with a neighboring naphthalene moiety of 1 and with the guest molecule via fairly short intermolecular face-to-face contacts. An intermolecular electrostatic interaction between a quinone-oxygen atom and a neighboring quinone ring is suggested by short O...C contacts. There exist several C-H...O hydrogen bonds, one of which seems to contribute to the formation of a centrosymmetric dimer of 1.     

An investigation of weak CH...O hydrogen bonds in maltose anomers by a combination of calculation and experimental solid-state NMR spectroscopy

Two-dimensional (1)H-(13)C MAS-J-HMQC solid-state NMR spectra of the two anomeric forms of maltose at natural abundance are presented. The experimental (1)H chemical shifts of the CH and CH(2) protons are assigned using first-principles chemical shift calculations that employ a plane-wave pseudopotential approach. Further calculations show that the calculated change in the (1)H chemical shift when comparing the full crystal and an isolated molecule is a quantitative measure of intermolecular C-H...O weak hydrogen bonding. Notably, a clear correlation between a large chemical shift change (up to 2 ppm) and both a short H...O distance (<2.7 A) and a CHO bond angle greater than 130 degrees is observed, thus showing that directionality is important in C-H...O hydrogen bonding.     

Interactions of D2O with methane and fluoromethane surfaces

TOF-SIMS is used to investigate the interactions between D2O and hydrophobic molecules, such as CH4, CH3F, CH2F2, CHF3, and CF4, at cryogenic temperatures (15 K). By irradiation with a 1.5-keV He+ beam, the D(+)(D2O)n ions are ejected efficiently from the D2O nanoclusters physisorbed on the CF4 layer due to Coulomb explosion: the ion yields are by about two orders of magnitude higher than those from a thick D2O layer via the kinetic sputtering. The D(+)(D2O)n yields decrease on the CHnF(4-n) layer with increasing the number of the C-H group. This is because the Coulombic fission is quenched due to the delocalization of valence holes through the C-H...H-C and C-H...D2O contacts. A pure D2O film is hardly grown on the CH4 layer as a consequence of intermixing whereas the D2O molecules basically adsorb on the surfaces of fluoromethanes, suggesting the attractive (water-repellent) interactions in the C-H...D2O (C-F...D2O) contacts. The C-H...O bond behaves like a conventional O-H...O hydrogen bond as far as the collision-induced proton transfer reaction is concerned.     

Intrinsic conformational characteristics of alpha,alpha-diphenylglycine

Quantum mechanical calculations at the B3LYP/6-31+G(d,p) level have been used to investigate the intrinsic conformational preferences of alpha,alpha-diphenylglycine, a simple alpha,alpha-dialkylated amino acid bearing two phenyl substituents on the alpha-carbon, in both the gas phase and aqueous solution. Nine minimum energy conformations have been characterized for the N-acetyl-N'-methylamide derivative within a relative energy range of about 9 kcal/mol. The relative stability of these structures is largely influenced by specific backbone...side chain and side chain...side chain interactions that can be attractive (N-H...pi and C-H...pi) or repulsive (C=O...pi). On the other hand, comparison with the minimum energy conformations calculated for alpha-aminoisobutyric acid, in which the two phenyl substituents are replaced by methyl groups, revealed that the bulky aromatic rings of alpha,alpha-diphenylglycine induce strain in the internal geometry of the peptide. Finally, a set of force-field parameters for classical Molecular Mechanics calculations was developed for the investigated amino acid. Molecular Dynamics simulations in aqueous solutions have been carried out to validate the parameters obtained.     

特殊缔合体系TFE水溶液分子动力学模拟

三氟乙醇(TFE)水溶液是一类特殊的缔合体系. 采用分子动力学模拟方法结合核磁共振化学位移研究了TFE水溶液体系全浓度范围的氢键网络, 并对动力学模拟结果和核磁共振化学位移进行了比较. 从径向分布函数(RDF)发现, TFE水溶液中存在着强氢键, 而体系中的C—H…O弱相互作用较为明显, 也不能忽略. 氢键网络分析发现TFE 水溶液体系的氢键大致分为以下三个区域: 在水富集区域, 水分子倾向于自身缔合形成稳定的簇结构, 随着TFE 浓度的增加, 水的有序结构受到破坏, 水分子和TFE分子发生交叉缔合作用形成氢键; 在TFE富集区域, 水分子较少, TFE分子自身通过氢键形成多缔体结构. 此外, 分子动力学统计的平均氢键数的变化和文献报导的核磁共振化学位移变化趋势相同, 实验和理论的结果吻合较好.