How Different are Electron-Rich and Electron-Deficient π Interactions?

Revealing the Role of N Heteroatoms in Noncovalent Aromatic Interactions by Ultrafast Intermolecular Coulombic Decay

Despite the widely recognized importance of noncovalent interactions involving aromatic rings in many fields, our understanding of the underlying forces and structural patterns, especially the impact of heteroaromaticity, is still incomplete. Here, we investigate the relaxation processes that follow inner-valence ionization in a range of molecular dimers involving various combinations of benzene, pyridine, and pyrimidine, which initiate an ultrafast intermolecular Coulombic decay process. Multiparticle coincidence momentum spectroscopy, combined with ab initio calculations, enables us to explore the principal orientations of these fundamental dimers and, thus, to elucidate the influence of N heteroatoms on the relative preference of the aromatic π-stacking, H-bonding, and CH-π interactions and their dependence on the number of nitrogen atoms in the rings. Our studies reveal a sensitive tool for the structural imaging of molecular complexes and provide a more complete understanding of the effects of N heteroatoms on the noncovalent aromatic interactions at the molecular level.

Hemiacetal-based dynamic systems: a new mechanistic insight

The formation of hemiacetals from pyrazine trifluoromethylketone as a model receptor and four simple alcohols was studied by using quantum chemical calculations and NMR spectroscopy. Free energy profiles for four types of mechanistic pathways were calculated and discussed with respect to kinetic and thermodynamic measurements. We show that hemiacetal formation is facilitated by an assisted proton transfer process via a pseudo eight-membered transition state which brings the theory and experiment into close agreement. Also, a newly proposed mechanistic pathway for hemiacetal formation via a five-membered transition state leading to zwitterionic intermediates is discussed. Direct proton transfer in a pseudo four-membered transition state can be ruled out due to the high energy of transition states with respect to other mechanistic pathways. We also show that in the case of hemiacetals, water and alcohol molecules cannot account sufficiently for the H-transfer process via six-membered transition states.

Parallelized Raman Difference Spectroscopy for the Investigation of Chemical Interactions

Raman spectroscopy provides an extremely high chemical selectivity. Raman difference spectroscopy is a technique to reveal even the smallest differences that occur due to weak interactions between substances and changes in the molecular structure. To enable parallelized and highly sensitive Raman difference spectroscopy in a microtiter-array, a diffractive optical element, a lens array, and a fiber bundle were integrated into a Raman spectroscopy setup in a unique fashion. The setup was evaluated with a microtiter-array containing pyridine-water complexes, and subwavenumber changes below the spectrometer's resolution could be resolved. The spectral changes were emphasized with two-dimensional correlation analysis. Density functional theory calculation and "atoms in molecule" analysis were performed to simulate the intermolecular long-range interactions between water and pyridine molecules and to get insight into the involved noncovalent interactions, respectively. It was found that by the addition of pyridine, the energy portion of hydrogen bonds to the total complexation energy between pyridine and water reduces. These results demonstrate the unique abilities of the new setup to investigate subtle changes due to biochemically important molecular interactions and opens new avenues to perform drug binding assays and to monitor highly parallelized chemical reactions.

Experimental and computational evidence for stabilising parallel, offset π[C(O)N(H)NC]⋯π(phenyl) interactions in acetohydrazide derivatives

Stabilising π[C(O)N(H)NC]⋯π(phenyl) interactions are described.

Cage-···Cage- Interaction: Boron Cluster-Based Noncovalent Bond and Its Applications in Solid-State Materials

Carboranes are boron-carbon clusters with important applications in the fields of materials, catalysis, pharmaceuticals, etc. However, the noncovalent interactions that could determine the solid-state structures and properties of such boron clusters have rarely been investigated. Herein, inspired by the coordinate bond in metallacarborane or ferrocene, the boron cluster-based noncovalent interaction (denoted as cage^-···cage^- interaction) between two nido-carborane clusters was successfully realized by using a pyridinium-based molecular barrier. The X-ray diffraction studies uncover that the cage^-···cage^- interaction has a contacting distance of 5.4-7.0 Å from centroid to centroid in the systems reported here. Theoretical calculations validate the formation of the noncovalent interaction and disclose its repulsive bonding nature that is overcome thanks to the positively charged pyridinium-based framework. Interestingly, such bulk crystalline materials containing the cage^-···cage^- interaction show relevant properties such as full-color absorption in the visible light range and important photothermal effect, which are absent for the control compound without carboranes. This study may offer fundamental insights into the boron cluster-based noncovalent interactions and open a new research avenue to rationally design boron cluster-based materials.

Anion-dependent structural variations and charge transport property analysis of 4′-(3-pyridyl)-4,2′:6′,4′′-terpyridinium salts

Anion-dependent structural variations and charge transport analysis of three 4′-(3-pyridyl)-4,2′:6′,4′′-terpyridinium salts are reported. They exhibit pronounced photosensing behavior when illuminated using visible light.

Recent Advances in the Fabrication and Application of Graphene Microfluidic Sensors

This review reports the progress of the recent development of graphene-based microfluidic sensors. The introduction of microfluidics technology provides an important possibility for the advance of graphene biosensor devices for a broad series of applications including clinical diagnosis, biological detection, health, and environment monitoring. Compared with traditional (optical, electrochemical, and biological) sensing systems, the combination of graphene and microfluidics produces many advantages, such as achieving miniaturization, decreasing the response time and consumption of chemicals, improving the reproducibility and sensitivity of devices. This article reviews the latest research progress of graphene microfluidic sensors in the fields of electrochemistry, optics, and biology. Here, the latest development trends of graphene-based microfluidic sensors as a new generation of detection tools in material preparation, device assembly, and chip materials are summarized. Special emphasis is placed on the working principles and applications of graphene-based microfluidic biosensors, especially in the detection of nucleic acid molecules, protein molecules, and bacterial cells. This article also discusses the challenges and prospects of graphene microfluidic biosensors.

Electrostatics does not dictate the slip-stacked arrangement of aromatic π-π interactions

Benzene dimer has long been an archetype for π-stacking. According to the Hunter–Sanders model, quadrupolar electrostatics favors an edge-to-face CH⋯π geometry but competes with London dispersion that favors cofacial π-stacking, with a compromise “slip-stacked” structure emerging as the minimum-energy geometry. This model is based on classical electrostatics, however, and neglects charge penetration. A fully quantum-mechanical analysis, presented here, demonstrates that electrostatics actually exerts very little influence on the conformational landscape of (C₆H₆)₂. Electrostatics also cannot explain the slip-stacked arrangement of C₆H₆⋯C₆F₆, where the sign of the quadrupolar interaction is reversed. Instead, the slip-stacked geometry emerges in both systems due to competition between dispersion and Pauli repulsion, with electrostatics as an ambivalent spectator. This revised interpretation helps to rationalize the persistence of offset π-stacking in larger polycyclic aromatic hydrocarbons and across the highly varied electrostatic environments that characterize π–π interactions in proteins.

According to the Hunter–Sanders model, geometries in π–π systems arise from competition between quadrupolar electrostatics (favoring an edge-to-face geometry) and London dispersion (favoring stacking), but this model misrepresents the molecular physics.

Interplay of π-stacking and inter-stacking interactions in two-component crystals of neutral closed-shell aromatic compounds: periodic DFT study

This paper bridges the gap between high-level ab initio computations of gas-phase models of 1 : 1 arene–arene complexes and calculations of the two-component (binary) organic crystals using atom–atom potentials. The studied crystals consist of electron-rich and electron-deficient compounds, which form infinite stacks (columns) of heterodimers. The sublimation enthalpy of crystals has been evaluated by DFT periodic calculations, while intermolecular interactions have been characterized by Bader analysis of the periodic electronic density. The consideration of aromatic compounds without a dipole moment makes it possible to reveal the contribution of quadrupole–quadrupole interactions to the π-stacking energy. These interactions are significant for heterodimers formed by arenes with more than 2 rings, with absolute values of the traceless quadrupole moment (Q_zz) larger than 10 D Å. The further aggregation of neighboring stacks is due to the C–H⋯F interactions in arene/perfluoroarene crystals. In crystals consisting of arene and an electron-deficient compound such as pyromellitic dianhydride, aggregation occurs due to the C–H⋯O interactions. The C–H⋯F and C–H⋯O inter-stacking interactions make the main contribution to the sublimation enthalpy, which exceeds 150 kJ mol⁻¹ for the two-component crystals formed by arenes with more than 2 rings.

The interplay of π-stacking and inter-stacking interactions in two-component organic crystals without conventional hydrogen bonds.

Influence of hydrogen bonds on edge-to-face interactions between pyridine molecules

Edge-to-face interactions between two pyridine molecules and the influence of simultaneous hydrogen bonding of one or both of the pyridines to water on those interactions were studied by analyzing data from ab initio calculations. The results show that the edge-to-face interactions of pyridine dimers that are hydrogen bonded to water are generally stronger than those of non-H-bonded pyridine dimers, especially when the donor pyridine forms a hydrogen bond. The binding energy of the most stable edge-to-face interacting H-bonded pyridine dimer is -5.05 kcal/mol, while that for the most stable edge-to-face interacting non-H-bonded pyridine dimer is -3.64 kcal/mol. The interaction energy data obtained in this study cannot be explained solely by the differences in electrostatic potential between pyridine and the pyridine-water dimer. However, the calculated cooperative effect can be predicted using electrostatic potential maps.

Chemical sensing with 2D materials

During the last decade, two-dimensional materials (2DMs) have attracted great attention due to their unique chemical and physical properties, which make them appealing platforms for diverse applications in sensing of gas, metal ions as well as relevant chemical entities.

Discovering the stacking landscape of a pyridine-pyridine system

Extremely extensive calculations of potential energy surfaces for the parallel-displaced configuration of pyridine dimer systems have been carried out using a dispersion-corrected density functional. Instead of focusing on stationary geometries these calculations provide much deeper insight into the "landscape" of the interaction energies of the particular systems-one can learn how the pyridine dimer stability changes along with various geometrical parameters. Other calculations such as natural bond orbital and energy decomposition have also been applied. The interplay of two significant factors, electrostatic forces and electron correlation effects, have been evaluated. The role of π···π interactions in the stacked pyridine systems has also been confirmed, and surprisingly, this happened to be true even for the geometries where the formation of C-H···π interactions might be proposed instead. The combination of many different methods has revealed the complexity of the stacking interactions. Apart from providing a "literal new look" into pyridine interaction patterns another picture has emerged. A stacking interaction in a pyridine dimer system is perceived as a combination of many different sources of the interaction energy, including orbital ones, and this is true for many different geometries.

Pharmacology and Toxicology of N-Benzylphenethylamine ("NBOMe") Hallucinogens

Serotonergic hallucinogens induce profound changes in perception and cognition. The characteristic effects of hallucinogens are mediated by 5-HT2A receptor activation. One class of hallucinogens are 2,5-dimethoxy-substituted phenethylamines, such as the so-called 2C-X compounds 2,5-dimethoxy-4-bromophenethylamine (2C-B) and 2,5-dimethoxy-4-iodophenethylamine (2C-I). Addition of an N-benzyl group to phenethylamine hallucinogens produces a marked increase in 5-HT2A-binding affinity and hallucinogenic potency. N-benzylphenethylamines ("NBOMes") such as N-(2-methoxybenzyl)-2,5-dimethoxy-4-iodophenethylamine (25I-NBOMe) show subnanomolar affinity for the 5-HT2A receptor and are reportedly highly potent in humans. Several NBOMEs have been available from online vendors since 2010, resulting in numerous cases of toxicity and multiple fatalities. This chapter reviews the structure-activity relationships, behavioral pharmacology, metabolism, and toxicity of members of the NBOMe hallucinogen class. Based on a review of 51 cases of NBOMe toxicity reported in the literature, it appears that rhabdomyolysis is a relatively common complication of severe NBOMe toxicity, an effect that may be linked to NBOMe-induced seizures, hyperthermia, and vasoconstriction.

The stacking interactions of bipyridine complexes: the influence of the metal ion type on the strength of interactions

The strength of the stacking interactions in the bipy complexes of nickel, palladium, and platinum, [M(CN)2 bipy]2 (M = Ni, Pd, Pt), was calculated using the ωB97xD/def2-TZVP method. The results show that for all considered geometries, interactions are the strongest for platinum, and weakest for nickel complexes, as a result of higher dispersion contributions of platinum over the palladium and nickel complexes. It was also shown that strength of interactions considerably rises with an increase of the stacking overlap area. As a consequence of the favorable electrostatic term, the strength of interactions also rises when metal atom and cyano ligands are involved in the overlap with bipy ligand. The strongest interaction was calculated in the platinum complex, for the geometry that has overlap of metal and cyano ligands with bipy ligand with an energy of -39.80 kcal mol(-1). The energies for similar geometries of palladium and nickel complexes are -34.60 and -32.45 kcal mol(-1). These energies, remarkably, exceed the strength of the stacking interactions between organic aromatic molecules. These results can be of importance in all systems with stacking interactions, from materials to biomolecules.

Four supramolecular isomers of dichloridobis(1,10-phenanthroline)cobalt(II): synthesis, structure characterization and isomerization

Crystal engineering can be described as the understanding of intermolecular interactions in the context of crystal packing and the utilization of such understanding to design new solids with desired physical and chemical properties. Free-energy differences between supramolecular isomers are generally small and minor changes in the crystallization conditions may result in the occurrence of new isomers. The study of supramolecular isomerism will help us to understand the mechanism of crystallization, a very central concept of crystal engineering. Two supramolecular isomers of dichloridobis(1,10-phenanthroline-κ(2)N,N')cobalt(II), [CoCl2(C12H8N2)2], i.e. (IA) (orthorhombic) and (IB) (monoclinic), and two supramolecular isomers of dichloridobis(1,10-phenanthroline-κ(2)N,N')cobalt(II) N,N-dimethylformamide monosolvate, [CoCl2(C12H8N2)2]·C3H7NO, i.e. (IIA) (orthorhombic) and (IIB) (monoclinic), were synthesized in dimethylformamide (DMF) and structurally characterized. Of these, (IA) and (IIA) have been prepared and structurally characterized previously [Li et al. (2007). Acta Cryst. E63, m1880-m1880; Cai et al. (2008). Acta Cryst. E64, m1328-m1329]. We found that the heating rate is a key factor for the crystallization of (IA) or (IB), while the temperature difference is responsible for the crystallization of (IIA) or (IIB). Based on the crystallization conditions, isomerization behaviour, the KPI (Kitajgorodskij packing index) values and the density data, (IB) and (IIA) are assigned as the thermodynamic and stable kinetic isomers, respectively, while (IA) and (IIB) are assigned as the metastable kinetic products. The 1,10-phenanthroline (phen) ligands interact with each other through offset face-to-face (OFF) π-π stacking in (IB) and (IIB), but by edge-to-face (EF) C-H...π interactions in (IA) and (IIA). Meanwhile, the DMF molecules in (IIB) connect to neighbouring [CoCl2(phen)2] units through two C-H...Cl hydrogen bonds, whereas there are no obvious interactions between DMF molecules and [CoCl2(phen)2] units in (IIA). Since OFF π-π stacking is generally stronger than EF C-H...π interactions for transition-metal complexes with nitrogen-containing aromatic ligands, (IIA) is among the uncommon examples that are stable and densely packed but that do not following Etter's intermolecular interaction hierarchy.

Fluorescence Quenching Nanoprobes Dedicated to In Vivo Photoacoustic Imaging and High-Efficient Tumor Therapy in Deep-Seated Tissue

Photoacoustic imaging (PAI) and photoacoustic (PA) therapy have promising applications for treating tumors. It is known that the utilization of high-absorption-coefficient probes can selectively enhance the PAI target contrast and PA tumor therapy efficiency in deep-seated tissue. Here, the design of a probe with the highest availability of optical-thermo conversion by using graphene oxide (GO) and dyes via π-π stacking interactions is reported. The GO serves as a base material for loading dyes and quenching dye fluorescence via fluorescence resonance energy transfer (FRET), with the one purpose of maximum of PA efficiency. Experiments verify that the designed fluorescence quenching nanoprobes can produce stronger PA signals than the sum of the separate signals generated in the dye and the GO. Potential applications of the fluorescence quenching nanoprobes are demonstrated, dedicating to enhance PA contrast of targets in deep-seated tissues and tumors in living mice. PA therapy efficiency both in vitro and in vivo by using the fluorescence quenching nanoprobes is found to be higher than with the commonly used PA therapy agents. Taken together, quenching dye fluorescence via FRET will provide a valid means for developing high-efficiency PA probes. Fluorescence quenching nanoprobes are likely to become a promising candidate for deep-seated tumor imaging and therapy.

Substituent Effects in π-Stacking of Histidine on Functionalized-SWNT and Graphene

Adsorptions of histidine on the functionalized (10,0) single-walled carbon nanotube (SWNT) and graphene were investigated using density function theory methods, M05-2x and DFT-D. The results show that the binding of the histidine ring to the functionalized SWNT is weaker than that to the pristine SWNT for both singlet and triplet complexes, regardless of the electron-donating (-OH, -NH₂) or electron-withdrawing (-COOH) character and their attached sites. The present decreased binding is opposite to the well-known enhanced binding in the substituted benzene dimers. Since the atoms of the histidine are distant from the substituent atoms by over 6Å, there would be no direct interaction between histidine and the substituent as in the case of the substituted benzene systems. The decreased binding can be mainly driven by the aromaticity of the functionalized SWNT. The nucleus-independent chemical shift (NICS) index analysis for the functionalized SWNTs in deed shows that local aromaticity of SWNT is decreased because of the electron redistribution induced by functional groups, and the π-π stacking between the histidine ring and functionalized-SWNT is therefore decreased as compared to the pristine SWNT. However, the above trend does not remain for the binding between the histidine and graphene. The binding of the histidine to the functionalized graphene with -OH and -NH₂ is just slightly weaker than that to the pristine graphene, while its binding to COOH-SWNT becomes a little bit stronger.

What are the preferred horizontal displacements of aromatic–aromatic interactions in proteins? Comparison with the calculated benzene–benzene potential energy surface

Stacking interactions of phenylalanine residues show preference for large offsets (3.5–5.0 Å), while the calculations show substantially strong interactions, of about −2.0 kcal mol⁻¹.

Influence of supramolecular structures in crystals on parallel stacking interactions between pyridine molecules

Parallel stacking interactions between pyridines in crystal structures and the influence of hydrogen bonding and supramolecular structures in crystals on the geometries of interactions were studied by analyzing data from the Cambridge Structural Database (CSD). In the CSD 66 contacts of pyridines have a parallel orientation of molecules and most of these pyridines simultaneously form hydrogen bonds (44 contacts). The geometries of stacked pyridines observed in crystal structures were compared with the geometries obtained by calculations and explained by supramolecular structures in crystals. The results show that the mean perpendicular distance (R) between pyridine rings with (3.48 Å) and without hydrogen bonds (3.62 Å) is larger than that calculated, because of the influence of supramolecular structures in crystals. The pyridines with hydrogen bonds show a pronounced preference for offsets of 1.25-1.75 Å, close to the position of the calculated minimum (1.80 Å). However, stacking interactions of pyridines without hydrogen bonds do not adopt values at or close to that of the calculated offset. This is because stacking interactions of pyridines without hydrogen bonds are less strong, and they are more susceptible to the influence of supramolecular structures in crystals. These results show that hydrogen bonding and supramolecular structures have an important influence on the geometries of stacked pyridines in crystals.

Energy component analysis of π interactions

Fundamental features of biomolecules, such as their structure, solvation, and crystal packing and even the docking of drugs, rely on noncovalent interactions. Theory can help elucidate the nature of these interactions, and energy component analysis reveals the contributions from the various intermolecular forces: electrostatics, London dispersion terms, induction (polarization), and short-range exchange-repulsion. Symmetry-adapted perturbation theory (SAPT) provides one method for this type of analysis. In this Account, we show several examples of how SAPT provides insight into the nature of noncovalent π-interactions. In cation-π interactions, the cation strongly polarizes electrons in π-orbitals, leading to substantially attractive induction terms. This polarization is so important that a cation and a benzene attract each other when placed in the same plane, even though a consideration of the electrostatic interactions alone would suggest otherwise. SAPT analysis can also support an understanding of substituent effects in π-π interactions. Trends in face-to-face sandwich benzene dimers cannot be understood solely in terms of electrostatic effects, especially for multiply substituted dimers, but SAPT analysis demonstrates the importance of London dispersion forces. Moreover, detailed SAPT studies also reveal the critical importance of charge penetration effects in π-stacking interactions. These effects arise in cases with substantial orbital overlap, such as in π-stacking in DNA or in crystal structures of π-conjugated materials. These charge penetration effects lead to attractive electrostatic terms where a simpler analysis based on atom-centered charges, electrostatic potential plots, or even distributed multipole analysis would incorrectly predict repulsive electrostatics. SAPT analysis of sandwich benzene, benzene-pyridine, and pyridine dimers indicates that dipole/induced-dipole terms present in benzene-pyridine but not in benzene dimer are relatively unimportant. In general, a nitrogen heteroatom contracts the electron density, reducing the magnitude of both the London dispersion and the exchange-repulsion terms, but with an overall net increase in attraction. Finally, using recent advances in SAPT algorithms, researchers can now perform SAPT computations on systems with 200 atoms or more. We discuss a recent study of the intercalation complex of proflavine with a trinucleotide duplex of DNA. Here, London dispersion forces are the strongest contributors to binding, as is typical for π-π interactions. However, the electrostatic terms are larger than usual on a fractional basis, which likely results from the positive charge on the intercalator and its location between two electron-rich base pairs. These cation-π interactions also increase the induction term beyond those of typical noncovalent π-interactions.