Periodic ab initio approach for the cooperative effect of CH/pi interaction in crystals: relative energy of CH/pi and hydrogen-bonding interactions

The bonding property of the CH/pi interaction in organic crystals has been investigated by the means of a periodic ab initio method. The energy of the CH(sp(2))/pi interaction in crystals, estimated with periodic RHF/6-21G*, showed a reasonable attractive CH(sp(2))/pi interaction owing to a cooperative effect, whereas the results calculated with RHF/cc-pVDZ indicate a negligibly small or repulsive interaction. The relative contribution of the CH(sp(2))/pi interaction to the column packing energy was found to be roughly half of the energy of a conventional hydrogen bond. The calculation of the charge distributions on the aromatic rings participating in the CH(sp(2))/pi interaction in crystals revealed that the atoms were more ionic than those in the gas phase. These theoretical calculations suggest a hydrogen-bonding characteristic for the CH(sp(2))/pi interaction in crystals, which does not occur in solution nor gas phase. We present computational evidence of the existence of the cooperative effect of CH(sp(2))/pi interaction in crystals.     

CH/pi interactions in methane clusters with polycyclic aromatic hydrocarbons

Geometries and interaction energies for methane clusters with naphthalene and pyrene were studied. Estimated CCSD(T) interaction energies for the clusters at the basis set limit were -1.92 and -2.50 kcal mol(-1), respectively. Dispersion is mainly responsible for the attraction. Electrostatic interaction is very small. Although the benzene-methane cluster prefers a monodentate structure, in which a C-H bond of the methane points toward the benzene, the methane clusters with the polycyclic aromatic hydrocarbons do not prefer monodentate structures. In the benzene-methane cluster, the weak electrostatic interaction stabilizes the monodentate structure. On the other hand the dispersion interaction controls the orientation of methane in the naphthalene and pyrene clusters. The dispersion interactions in these clusters are significantly larger than those in the benzene-methane cluster. The methane prefers the orientation which is suitable for stabilization by dispersion. Hydrogen atoms of the methane locate above the centers of hexagonal rings of the polycyclic aromatic hydrocarbons in the stable structures. The structures have a small steric repulsion and this positions them only a short distance from the aromatic plane. The large dispersion contribution to the attraction shows that interactions between carbon atoms are mainly responsible for the attraction, and that hydrogen atoms are not important for the attraction. This shows that the interactions between the methane and polycyclic aromatic hydrocarbons are not pi-hydrogen bonds.     

Saturated hydrocarbon-benzene complexes: theoretical study of cooperative CH/pi interactions

High-level ab initio calculations at the CCSD(T)/aug-cc-pVTZ//MP2/aug(d,p)-6-311G(d,p) level were employed to investigate the cooperative CH/pi effects between the pi face of benzene and several modeled saturated hydrocarbons, propane, isobutane, cyclopropane, cyclobutane, cyclopentane, cyclohexane, cyclopentane, cyclooctane, and bicyclo[2.2.2]octane. In all cases, multiple C-H groups (2-4) are found to interact with the pi face of benzene, with one C-H group pointing close to the center of the benzene ring. The geometries of these complexes are governed predominantly by electrostatic interaction between the interacting systems. The calculated interaction energies (10-14 kJ mol(-1)) are 2-3 times larger than that of the prototypical methane-benzene complex. The trends of geometries, interaction energies, binding properties, as well as electron-density topological properties were analyzed. The calculated interaction energies correlate well with the polarizabilities of the hydrocarbons. AIM analysis confirms the hydrogen-bonded nature of the CH/pi interactions. Significant changes in proton chemical shift and stretching frequency (blue shift) are predicted for the ring C-H bond in these complexes.     

CH/pi interactions in DNA and proteins. A theoretical study

A systematic study of the CH/pi interactions of methane with the purine and pyrimidine bases of nucleic acids and with the lateral chains of the four natural aromatic amino acids has been carried out for the first time. The MPWB1K/6-31+G(d,p) method has shown to be adequate for the study of these weak interactions in which dispersion forces play a main role. It has been shown that two different kinds of clusters exist, depending on whether one or two CH bonds point to the aromatic system. The latter one, which we have called bifurcated, is usually more stable. With regard to aromatic amino acids, our calculations agree with experimental data in the fact that tryptophan leads to the strongest interaction, while hystidine leads to the weakest one. In the case of nucleic acid bases, the differences in binding energies are not large. This is specially true for thymine and uracil, showing that these two bases have a similar acceptor character in CH/pi interactions.     

Noncovalent interactions between cytosine and SWCNT: curvature dependence of complexes via pi...pi stacking and cooperative CH...pi/NH...pi

Fragments of C24H12, adapted from a variety of armchair [(n,n), (n = 5, 7, and 8)] and zigzag [(m,0) (m = 8, 10, and 12)] single-walled carbon nanotube (SWCNT), are used to model corresponding SWCNTs with different diameters and electronic structures. The parallel binding mainly through pi...pi stacking interaction, as well as the perpendicular binding via cooperative NH...pi and CH...pi between cytosine and the fragments of SWCNT have been extensively investigated with a GGA type of DFT, PW91LYP/6-311++G(d,p). The eclipsed tangential (ET) conformation with respect to the six-membered ring of cytosine and the central ring of SWCNT fragments is less stable than the slipped tangential (ST) conformation for the given fragment; perpendicular conformations with NH2 and CH ends have higher negative binding energy than those with NH and CH ends. At PW91LYP/6-311++G(d,p) level, two tangential complexes are less bound than perpendicular complexes. However, as electron correlation is treated with MP2/6-311G(d,p) for PW91LYP/6-311++G(d,p) optimized complexes, it turns out there is an opposite trend that two tangential complexes become more stable than three perpendicular complexes. This result implies that electron correlation, a primary source to dispersion energy, has more significant contributions to the pi...pi stacking complexes than to the complexes via cooperative NH...pi and CH...pi interactions. In addition, it was found for the first time that binding energies for two tangential complexes become more negative with increasing nanotube diameter, while those for three perpendicular complexes have a weaker dependence on the curvature; i.e., binding energies are slightly less and less negative. The performance of a novel hybrid DFT, MPWB1K, was also discussed.     

Crystal engineering for topochemical polymerization of muconic esters using halogen-halogen and CH/pi interactions as weak intermolecular interactions

We now report the molecular and crystal structure design of muconic ester derivatives on the basis of crystal engineering using halogen-halogen contacts and CH/pi interactions. The solid-state photoreaction pathway of the dibenzyl (Z,Z)-muconates as the 1,3-diene dicarboxylic acid monomers depends on the structure of the ester groups. The substitution of a halogen atom for the aromatic hydrogen of a benzyl group induces topochemical polymerization to produce stereoregular polymers in a crystalline form, whereas the unsubstituted benzyl derivative isomerizes to yield the corresponding E,E isomer under similar conditions. The topochemical polymerization process is directly confirmed by the fact that the single-crystal structures before and after the polymerization are very similar to each other. From the crystal structure analysis for a series of substituted benzyl (Z,Z)- and (E,E)-muconates, it has been revealed that the planar diene moieties are closely packed to form a columnar structure in the crystals. The stacking of the polymerizable monomers is characterized by a stacking distance of 4.9-5.2 A along the columns. This structure is supported by a halogen-halogen interaction between the chlorine or bromine atoms introduced at the p position of the benzyl groups in addition to an aromatic stacking due to the CH/pi interaction between the benzylic methylene hydrogens and aromatic rings. The design of a monomer packing corresponds to the type and position of the introduced halogen atom and also the polymorphs. To make a stacking distance of 5 A using both halogen-halogen and CH/pi interactions as supramolecular synthons is important for the molecular design of muconic ester derivatives appropriate for topochemical polymerization.     

CH/pi interactions involving aromatic amino acids: refinement of the CHARMM tryptophan force field

High-level ab initio calculations have been carried out to study weak CH/pi interactions and as a check of the CHARMM force field for aromatic amino acids. Comparisons with published data indicate that the MP2/cc-pVTZ level of theory is suitable for calculations of CH/pi interaction, including the T-shape benzene dimer. This level of theory was, therefore, applied to investigate CH/pi interactions between ethene or cis-2-butene and benzene in a variety of orientations. In addition, complexes between ethene and a series of model compounds (toluene, methylindole and p-cresol) representing the aromatic amino acids were studied motivated by the presence of CH/pi interactions in biological systems. Ab initio binding energies were compared to the binding energies obtained with the CHARMM22 force field. In the majority of orientations, CHARMM22 reproduces the preferred binding modes, with excellent agreement for the benzene dimer. Small discrepancies found in the calculations involving methylindole along with a survey of published thermodynamic data for the aromatic amino acids prompted additional optimization of the tryptophan force field. Partial atomic charges, Lennard-Jones parameters, and force constants were improved to obtain better intra- and intermolecular properties, with significant improvements obtained in the reproduction of experimental heats of sublimation for indole and free energies of aqueous solvation for methylindole.     

Rational design of CH/pi interaction sites in a basic resolving agent

A novel synthetic basic resolving agent, cis-1-aminobenz[f]indan-2-ol (ABI), was rationally designed by introducing effective CH/pi interaction sites to cis-1-aminoindan-2-ol (AI), whose chiral recognition ability has been reported from our laboratory. ABI was applicable to a wide variety of racemic arylalkanoic acids and showed moderate to excellent chiral recognition ability, which was obviously higher than that of AI. The fundamental and important role of CH/pi interactions, such as tunable CH(sp(2))/pi and CH(sp(3))/pi interactions, in the chiral recognition by ABI was revealed by X-ray crystallographic study.     

The role of CH/pi interaction in the stabilization of less-soluble diastereomeric salt crystals

Enantiopure 2-naphthylglycolic acid (NGA) and cis-1-aminobenz[f]indan-2-ol (ABI) were rationally designed as new resolving agents on the model of mandelic acid (MA) and cis-1-aminoindan-2-ol (AI), respectively. As expected, NGA and ABI showed superior chiral recognition ability to racemates, compared with MA and AI. In order to clarify any factors governing the chiral recognition abilities of NGA and ABI, the crystal structures of their less- and more-soluble diastereomeric salts were determined by X-ray crystallographic analyses and revealed that CH/pi interactions play an intrinsic role in chiral recognitions. A theoretical investigation was also performed with the periodic ab initio method by using the X-ray crystal structures of the less-soluble salt crystals with AI and ABI to find the unique properties of CH/pi interaction in the crystalline state, which largely contributed to the stabilization of the crystals.     

Multicoefficient extrapolated density functional theory studies of pi...pi interactions: the benzene dimer

We report tests of new (2005) and established (1999-2003) multilevel methods against essentially converged benchmark results for nonbonded interactions in benzene dimers. We found that the newly developed multicoefficient extrapolated density functional theory (DFT) methods (which combine DFT with correlated wave function methods) give better performance than multilevel methods such as G3SX, G3SX(MP3), and CBS-QB3 that are based purely on wave function theory (WFT); furthermore, they have a lower computational cost. We conclude that our empirical approach for combining WFT methods with DFT methods is a very efficient and effective way for describing not only covalent interactions (as shown previously) but also nonbonded interactions.     

Calorimetric measurement of the CH/pi interaction involved in the molecular recognition of saccharides by aromatic compounds

Can a benzene molecule differentiate between two isomeric carbohydrates? It is generally accepted that two factors govern molecular recognition: complementarity and preorganization. Preorganization requires the presence of cavities for positioning the host's groups of complementary nature to those of the guest. This study shows that, in fact, groups should be complementary to recognize each other (for the case presented here, it is controlled by the CH/pi interaction) but preorganization is not essential. Since weak interactions have their origin in dispersion forces, they also have impact on the enthalpic term of the free energy, so it was considered that their participation can be demonstrated by measuring the energy involved. For recognition to happen, two conditions must be satisfied: specificity and associated stabilizing energy. In this study we evaluated the heat of dissolution of different carbohydrates such as methyl 2,3,4,6-tetra-O-methyl-alpha-d-mannopyranoside and methyl 2,3,4,6-tetra-O-methyl-beta-d-galactopyranoside using different aromatic solvents. The solvation enthalpies in benzene were -78.8 +/- 3.9 and -88.7 +/- 5.5 kJ mol(-1) for each carbohydrate, respectively; and these values yielded a CH/pi energy of interaction of 9.9 kJ mol(-1). In addition, NMR studies of the effect of the addition of benzene to chloroform solutions of the two carbohydrates showed that benzene specifically interacts with the hydrogen atoms of the pyranose ring at positions 3, 4, and 5 located on the alpha face of the methyl-beta-galactoside, so it is, in fact, able to recognize it. Thus, the interactions between carbohydrates and the aromatic residues of proteins occur in the absence of the confinement generated by the protein structure. By experimentally measuring the energy associated with this interaction and comparing it to theoretical calculations, it was also possible to unequivocally determine the existence of CH/pi interactions between carbohydrates and proteins.     

Electron delocalization in the metallabenzenes: a computational analysis of ring currents

Many metallabenzene complexes appear to exhibit an enhanced thermodynamic stability which has been attributed to the concept of aromaticity. Analysis of the ring currents induced by a magnetic field, either by direct visualization or by considering nuclear or nucleus-independent chemical shielding values (NMR or NICS), have become useful theoretical tools to characterize the aromaticity of many molecules involving the main group elements. We have analyzed 21 metallabenzenes using variations of these techniques, which take account of the large core and metal orbital contributions which often lead to transition-metal-containing systems exhibiting anomalous shielding values. Analysis of individual orbital contributions to both the ring currents and chemical shielding values based upon the ipsocentric and CSGT (continuous set of gauge transformations) methods has shown that complexes such as the 18 electron Ir or Rh(C 5H 5)(PH 3) 2Cl 2 molecules should be classed as aromatic, whereas the 16 electron complexes such as Os or Ru(C 5H 5)(PH 3) 2Cl 2 should not, despite having the same occupancy of pi-MOs. The differences can be directly attributed to the HOMO/LUMO b 2 in-plane (d xy ) molecular orbital, which, when unoccupied, is available to disrupt the delocalized currents typical of aromatic systems. A range of Pd and Pt metallabenzenes with cyclopentadienyl and phosphine ligands is also discussed as having aromatic and nonaromatic character, respectively.     

A single CH/pi weak hydrogen bond governs stability and the photocycle of the photoactive yellow protein

Importance of the CH/pi interaction on the structure and function of the photoactive yellow protein (PYP) was substantiated. Focusing on the phenyl ring of Phe6 adjacent to the alkyl chain of Lys123, the mutants for these amino acid residues were characterized. The results demonstrated that the mutants lacking the pi-electron at position 6 or the alkyl chain at position 123 show substantial malfunction. This is a clear example that single CH/pi weak interaction plays a crucial role in the normal action of the protein.     

Magnitude of CH/O interactions between carbohydrate and water

The interaction energy potentials for six orientations of the fucose–water complex were calculated for evaluating the magnitude of the CH/O interactions in the complex. The calculations show that the C–H bonds of the nonpolar surface of fucose prefer to have contact with the oxygen atom of the water. The substantial attraction exists between the C–H bonds of fucose and water. The interaction energy calculated for the fucose–water complex at the MP2/aug-cc-pVTZ level is −2.55 kcal/mol. The CH/O interactions in the fucose–water complex are significantly larger than those in the cyclohexane–water complex (−1.13 kcal/mol), which shows that the oxygen atoms of fucose enhance the CH/O interactions. The electrostatic and dispersion interactions are responsible for the attraction in the CH/O interactions in the fucose–water complex, while the electrostatic contributions to the attraction in the CH/O interactions in the cyclohexane–water complex is small. The DFT-SAPT calculations also show that the electrostatic interactions are responsible for the larger attraction in the fucose–water complex. These results suggest that the nature of the CH/O interactions between carbohydrate and water is significantly different from that of the CH/O interactions between saturated hydrocarbon and water.     

Nature and physical origin of CH/pi interaction: significant difference from conventional hydrogen bonds

Recently reported high-level ab initio calculations and gas phase spectroscopic measurements show that the nature of CH/pi interactions is considerably different from conventional hydrogen bonds, although the CH/pi interactions were often regarded as the weakest class of hydrogen bonds. The major source of attraction in the CH/pi interaction is the dispersion interaction and the electrostatic contribution is small, while the electrostatic interaction is mainly responsible for the attraction in the conventional hydrogen bonds. The nature of the "typical" CH/pi interactions is similar to that of van der Waals interactions, if some exceptional "activated" CH/pi interactions of highly acidic C-H bonds are excluded. Shifts of C-H vibrational frequencies and electronic spectra also support the similarity. The hydrogen bond is important in controlling structures of molecular assemblies, since the hydrogen bond is sufficiently strong and directional due to the large electrostatic contribution. On the other hand, the directionality of the "typical" CH/pi interaction is very weak. Although the "typical" CH/pi interaction is often regarded as an important interaction in controlling the structures of molecular assemblies as in the cases of conventional hydrogen bonds, the importance of the "typical" CH/pi interactions is questionable.     

Cooperation of multiple CH...pi interactions to stabilize polymers in aromatic nanochannels as indicated by 2D solid state NMR

Advanced 2D solid state NMR techniques reveal weak intermolecular interactions that cooperatively sustain nanostructures of high molecular mass aliphatic polymers entrapped as guests in channels formed by an aromatic host.     

Magnitude and directionality of the interaction energy of the aliphatic CH/pi interaction: significant difference from hydrogen bond

The CCSD(T) level interaction energies of CH/pi complexes at the basis set limit were estimated. The estimated interaction energies of the benzene complexes with CH(4), CH(3)CH(3), CH(2)CH(2), CHCH, CH(3)NH(2), CH(3)OH, CH(3)OCH(3), CH(3)F, CH(3)Cl, CH(3)ClNH(2), CH(3)ClOH, CH(2)Cl(2), CH(2)FCl, CH(2)F(2), CHCl(3), and CH(3)F(3) are -1.45, -1.82, -2.06, -2.83, -1.94, -1.98, -2.06, -2.31, -2.99, -3.57, -3.71, -4.54, -3.88, -3.22, -5.64, and -4.18 kcal/mol, respectively. Dispersion is the major source of attraction, even if substituents are attached to the carbon atom of the C-H bond. The dispersion interaction between benzene and chlorine atoms, which is not the CH/pi interaction, is the cause of the very large interaction energy of the CHCl(3) complex. Activated CH/pi interaction (acetylene and substituted methanes with two or three electron-withdrawing groups) is not very weak. The nature of the activated CH/pi interaction may be similar to the hydrogen bond. On the other hand, the nature of other typical (nonactivated) CH/pi interactions is completely different from that of the hydrogen bond. The typical CH/pi interaction is significantly weaker than the hydrogen bond. Dispersion interaction is mainly responsible for the attraction in the CH/pi interaction, whereas electrostatic interaction is the major source of attraction in the hydrogen bond. The orientation dependence of the interaction energy of the typical CH/pi interaction energy is very small, whereas the hydrogen bond has strong directionality. The weak directionality suggests that the hydrogen atom of the interacting C-H bond is not essential for the attraction and that the typical CH/pi interaction does not play critical roles in determining the molecular orientation in molecular assemblies.     

Intramolecular rearrangement for regioselective complexation by intramolecular CH/pi interaction in a hydrophobic cavity of a ruthenium coordination sphere

A Ru(II) complex with a hydrophobic cavity formed from two 1-naphthoylamide groups was prepared. Its reactions with beta-diketones gave beta-diketonato complexes in which hydrophobic pi-pi or CH/pi interactions were confirmed by NMR spectroscopy and X-ray crystallography. In the case of the asymmetric beta-diketone benzoylacetone, an isomer with a CH/pi interaction was afforded as the sole product owing to thermodynamic control. The reaction was found to involve a novel intramolecular rearrangement from the phenyl-included isomer to the methyl-included one without rupture of Ru-beta-diketonato coordination bonds (activation energy 52 kJ mol(-1)). This indicates that CH/pi interactions can be more favored thermodynamically than pi-pi interactions in a suitable hydrophobic environment.