Research progress in cation-π interactions

Cation-π interaction is a potent intermolecular interaction between a cation and an aromatic system,which has been viewed as a new kind of binding force,as being compared with the classical interactions(e.g. hydrogen bonding,electrostatic and hydrophobic interactions). Cation-π interactions have been observed in a wide range of biological contexts. In this paper,we present an overview of the typical cation-π interactions in biological systems,the experimental and theoretical investigations on cation-π interactions,as well as the research results on cation-π interactions in our group.     

The open-close mechanism of M2 channel protein in influenza A virus:A computational study on the hydrogen bonds and cation-π interactions among His37 and Trp41

The M2 protein from influenza A virus is a tetrameric ion channel. It was reported that the permeation of the ion channel is correlated with the hydrogen bond network among His37 residues and the cation-π interactions between His37 and Trp41. In the present study,the hydrogen bonding network of 4-methyl-imidazoles was built to mimic the hydrogen bonds between His37 residues,and the cation-π interactions between 4-methyl-imidazolium and indole systems were selected to represent the interac-tions between His37 and Trp41. Then,quantum chemistry calculations at the MP2/6-311G level were carried out to explore the properties of the hydrogen bonds and the cation-π interactions. The calcula-tion results indicate that the binding strength of the N-H···N hydrogen bond between imidazole rings is up to -6.22 kcal·mol-1,and the binding strength of the strongest cation-π interaction is up to -18.8 kcal·mol-1(T-shaped interaction) or -12.3 kcal·mol-1(parallel stacking interaction). Thus,the calcu-lated binding energies indicate that it is possible to control the permeation of the M2 ion channel through the hydrogen bond network and the cation-π interactions by altering the pH values.     

Investigation of the Substituent Effects on π-Type Pnicogen Bond Interaction

Intermolecular interactions between PH_2Cl and Ar–R(R = H,OH,NH_2,CH_3,Br,Cl,F,CN,NO_2) were calculated by using MP2/aug-cc-p VDZ quantum chemical method.It has been shown from our calculations that the aromatic rings with electron-withdrawing groups represent much weaker binding affinities than those with electron-donating groups.The charge-transfer interaction between PH_2Cl and Ar–R plays an important role in the formation of pnicogen bond complexes,as revealed by NBO analysis.Nevertheless,AIM analysis shows that the nature of the interactions between PH_2Cl and Ar–R is electrostatic,and the interaction energies of the complexes are correlated positively with the electron densities in the bond critical points(BCPs).RDG/ELF graphical analyses were performed to visualize the positions and strengths of the pnicogen bonding,as well as the spatial change of the electron localization upon the formation of complexes.The π-type halogen bond was also calculated,and it has been revealed that the π-type pnicogen bond systems are more stable than the halogen bond ones.     

Theoretical Study of Pnicogen Bonding Interactions between ECl_3(P,As) and Different Electron Donors

The pnicogen bond interaction between different electron donors(anion, π-electron, heteroatom) and ECl_3(E = As, P) was calculated by the method of MP_2/aug-cc-p VTZ. It has been indicated that the pnicogen bonds of complex formed by the anion and ECl_3 are more stable than that by the neutral electron donor, in which the pnicogen bonds of complex formed by NH_3 and ECl_3 are the most stable, and that by H_2S and ECl_3 is the least stable. The nature of pnicogen bond interaction is the closed shell interaction by AIM analysis, and BCP electron density is positively correlated to the complex interaction energy. RDG and DDF graphical analyses are performed to visualize the nature of pnicogen bond interaction from different donors, the position and strength of the pnicogen bond interaction, as well as the rearrangement of electron density after the formation of pnicogen bond system.     

Molecular and quantum mechanical studies on the monomer recognition of a highly-regular β-helical antifreeze protein

The possible interaction models for an antifreeze protein from Tenebrio molitar (TmAFP) have been systematically studied using the methods of molecular mechanics, molecular dynamics and quantum chemistry. It is hoped that these approaches would provide insights into the nature of interaction between protein monomers through sampling a number of interaction possibilities and evaluating their interaction energies between two monomers in the course of recognition. The results derived from the molecular mechanics indicate that monomer's β-sheets would be involved in interaction area and the side chains on two β-faces can match each other at the two-dimensional level. The results from molecular mechanics and ONIOM methods show that the strongest interaction energy could be gained through the formation of H-bonds when the two β-sheets are involved in the interaction model. Furthermore, the calculation of DFT and analysis of van der Waals bond charge density confirm further that recognition between the two     

Does the molecular structure of CaH2 affect the dihydrogen bonding in CaH2 … HY (Y = CH3, C2H3, C2H, CN, and NC) complexes? A quantum chemistry study using MP2 and B3LYP methods

Second-order M ller-Plesset(MP2) and density functional theory(DFT) calculations have been carried out in order to investigate the structures and properties of dihydrogen-bonded CaH 2 HY(Y = CH 3,C 2 H 3,C 2 H,CN,and NC) complexes.Our calculations revealed two possible structures for CaH 2 in CaH 2 HY complexes:linear(I) and bent(II).The bond lengths,interaction energies,and strengths for H H interactions obtained by both MP2 and B3LYP methods are quite close to each other.It was found that the interaction energy decreases with increasing electron density at the Ca-H bond critical point.Atom-in-molecule(AIM) results show that for all of Ca-H H-Y interactions considered here,the Laplacian of the electron density at the H H bond critical point is positive,indicating the electrostatic nature of these Ca-H H-Y dihydrogen bonded systems.     

Interaction specific binding hotspots in Endonuclease colicin-immunity protein complex from MD simulations

The binding of Endonuclease colicin 9 (E9) by Immunity protein 9 (Im9) was found to involve some hotspots from helix III of Im9 on protein-protein interface that contribute the dominant binding energy to the complex.In the current work,MD simulations of the WT and three hotspot mutants (D51A,Y54A and Y55A of Im9) of the E9-Im9 complexes were carried out to investigate specific interaction mechanisms of these three hotspot residues.The changes of binding energy between the WT and mutants of the complex were computed by the MM/PBSA method using a polarized force field and were in excellent agreement with experiment values,verifying that these three residues were indeed hotspots of the binding complex.Energy decomposition analysis revealed that binding by D51 to E9 was dominated by electrostatic interaction due to the presence of the carboxyl group of Asp51 which hydrogen bonds to K89.For binding by hotspots Y54 and Y55,van der Waals interaction from the aromatic side chain of tyrosine provided the dominant interaction.For comparison,calculation by using the standard (nonpolarizable) AMBER99SB force field produced binding energy changes from these mutations in opposite direction to the experimental observation.Dynamic hydrogen bond analysis showed that conformations sampled from MD simulation in the standard AMBER force field were distorted from the native state and they disrupted the inter-protein hydrogen bond network of the protein-protein complex.The current work further demonstrated that electrostatic polarization plays a critical role in modulating protein-protein binding.     

Halogen Bonding： An AIM Analysis of the Weak Interactions

A series of complexes formed between halogen-containing molecules and ammonia have been investigated by means of the atoms in molecules （AIM） approach to gain a deeper insight into halogen bonding. The existence of the halogen bond critical points （XBCP） and the values of the electron density （Pb） and Laplacian of electron density （V2pb） at the XBCP reveal the closed-shell interactions in these complexes. Integrated atomic properties such as charge, energy, polarization moment, volume of the halogen bond donor atoms, and the corresponding changes （△） upon complexation have been calculated. The present calculations have demonstrated that the halogen bond represents different AIM properties as compared to the well-documented hydrogen bond. Both the electron density and the Laplacian of electron density at the XBCP have been shown to correlate well with the interaction energy, which indicates that the topological parameters at the XBCP can be treated as a good measure of the halogen bond strength In addition, an excellent linear relationship between the interatomic distance d（X…N） and the logarithm of Pb has been established.     

Cation-π Interaction between the Aromatic Organic Counterion and DTAB Micelle in Mixed Solvents

The cation-π interaction between the aromatic organic counterion potassium hydrogen phthalate (KHP) and DTAB micelle in aqueous mixture of EG was investigated, using the techniques of conductivity measurements, UV absorption spectrum and NMR spectrum. The conductivity and UV spectrum studies were with respect to the effect of KHP on DTAB and that of DTAB micelle on KHP, respectively. According to the chemical shift changes of the aromatic ring and the surfactant methylene protons, it can be assumed that KHP penetrated into DTAB micelle with its carboxylic group protruding out of the micellar surface. And the strength of the interaction became weaker with the content of EG in the mixed solvent increasing.     

EFFECT OF pH ON THE ADSORPTION OF p-AMINOBENZOIC ACID ON POLYSTYRENE-BASED ADSORBENTS

In the present study the adsorptive properties of p-aminobenzoic acid with hypercrosslinked and multi-functional polymeric adsorbents at different solution pHs were systematically investigated in accordance with the particular physicochemical characteristics of the aromatic amphoteric compound involving both Lewis acid and Lewis base functional groups. It was found that the equilibrium adsorption data of the three polymeric adsorbents fitted well in the Langmuir and Freundlich isotherm equations. Studies at various pH levels indicate that the capacity of the adsorbents for adsorption of the ionic forms of adsorbate is less than that for the corresponding neutral species. At pH 3.78, the adsorption capacities of the three adsorbents are the highest. Whereas the adsorption property of multi-functional polymeric adsorbent NJ-99 is the largest, which may be attributed to the strong hydrogen-bonding interaction between the amino groups on the resin and the carboxyl group of p-aminobenzoic acid. The trend of the adsorption capacities of the three adsorbents towards p-aminobenzoic acid with the solution pH is in accord with the dissociation curve of the neutral molecular p-aminobenzoic acid. The adsorption forces include π-π interaction, hydrogen-bonding interaction and electrostatic attraction or repulsion when there exist the molecular and ionic adsorbing species at different pHs in aqueous solution.     

Tuning anion‐π interaction via halogen substituent effects in cyanuric acids and its derivatives

Several anion‐π complexes of isocyanuric acid, thioderivatives and their halogen substituted derivatives with chloride anion have been studied. The geometric and energetic features, charges transfer from chloride anion to the aromatic rings and “atoms‐in‐molecules” analysis are performed and discussed for these complexes. The results show that the strength of the anion‐π interaction between cyanuric derivatives and chloride anion can be tuned by halogen‐substituting. The localized molecular orbital energy decomposition analyses shows that, in the total interaction, exchange and electrostatic energies are the dominant stabilizing forces, and the polarization energies also make a favorable contribution. Finally, solvent effect significantly weakens the anion‐π interaction between the isocyanuric derivatives and chloride anion, especially in polar solvents. © 2015 Wiley Periodicals, Inc.     

Polyether Natural Product Inspired Cascade Cyclizations: Autocatalysis on π‐Acidic Aromatic Surfaces

Anion‐π catalysis functions by stabilizing anionic transition states on aromatic π surfaces, thus providing a new approach to molecular transformation. The delocalized nature of anion–π interactions suggests that they serve best in stabilizing long‐distance charge displacements. Aiming therefore for an anionic cascade reaction that is as charismatic as the steroid cyclization is for conventional cation‐π biocatalysis, reported here is the anion‐π‐catalyzed epoxide‐opening ether cyclizations of oligomers. Only on π‐acidic aromatic surfaces having a positive quadrupole moment, such as hexafluorobenzene to naphthalenediimides, do these polyether cascade cyclizations proceed with exceptionally high autocatalysis (rate enhancements k_auto/k_cat >10⁴ m ^?1). This distinctive characteristic adds complexity to reaction mechanisms (Goldilocks‐type substrate concentration dependence, entropy‐centered substrate destabilization) and opens intriguing perspectives for future developments.     

Enhanced Aerogen–π Interaction by a Cation–π Force

The interaction between a noble gas atom and an aromatic π‐electron system, which mainly originates from the London dispersion force, is very weak and has not attracted enough attention yet. Herein, we reported a type of notably enhanced aerogen–π interaction between cation–π systems and noble gas atoms. The binding strength of a divalent cation–π system with a xenon atom is comparable to a moderate hydrogen bond (up to ca. 7 kcal mol^?1), whereas krypton and argon atoms produce slightly weaker interactions. Energy‐decomposition analysis reveals that the induction interaction is responsible for the stabilization of divalent cation–π?Xe species besides the dispersion interaction. Our results might be helpful to increase the understanding of some unsolved mysteries of aerogens.     

The Prominent Enhancing Effect of the Cation–π Interaction on the Halogen–Hydride Halogen Bond in M1⋅⋅⋅C6H5X⋅⋅⋅HM2

We designed M¹???C₆H₅X???HM² (M¹=Li⁺, Na⁺; X=Cl, Br; M²=Li, Na, BeH, MgH) complexes to enhance halogen–hydride halogen bonding with a cation–π interaction. The interaction strength has been estimated mainly in terms of the binding distance and the interaction energy. The results show that halogen–hydride halogen bonding is strengthened greatly by a cation–π interaction. The interaction energy in the triads is two to six times as much as that in the dyads. The largest interaction energy is ?8.31 kcal mol^?1 for the halogen bond in the Li⁺???C₆H₅Br???HNa complex. The nature of the cation, the halogen donor, and the metal hydride influence the nature of the halogen bond. The enhancement effect of Li⁺ on the halogen bond is larger than that of Na⁺. The halogen bond in the Cl donor has a greater enhancement than that in the Br one. The metal hydride imposes its effect in the order HBeH<HMgH<HNa<HLi for the Cl complex and HBeH<HMgH<HLi<HNa for the Br complex. The large cooperative energy indicates that there is a strong interplay between the halogen–hydride halogen bonding and the cation–π interaction. Natural bond orbital and energy decomposition analyses indicate that the electrostatic interaction plays a dominate role in enhancing halogen bonding by a cation–π interaction.     

Understanding the Fundamental Role of π/π, σ/σ, and σ/π Dispersion Interactions in Shaping Carbon‐Based Materials

Noncovalent interactions involving aromatic rings, such as π‐stacking and CH/π interactions, are central to many areas of modern chemistry. However, recent studies proved that aromaticity is not required for stacking interactions, since similar interaction energies were computed for several aromatic and aliphatic dimers. Herein, the nature and origin of π/π, σ/σ, and σ/π dispersion interactions has been investigated by using dispersion‐corrected density functional theory, energy decomposition analysis, and the recently developed noncovalent interaction (NCI) method. Our analysis shows that π/π and σ/σ stacking interactions are equally important for the benzene and cyclohexane dimers, explaining why both compounds have similar boiling points. Also, similar dispersion forces are found in the benzene???methane and cyclohexane???methane complexes. However, for systems larger than naphthalene, there are enhanced stacking interactions in the aromatic dimers adopting a parallel‐displaced configuration compared to the analogous saturated systems. Although dispersion plays a decisive role in stabilizing all the complexes, the origin of the π/π, σ/σ, and σ/π interactions is different. The NCI method reveals that the dispersion interactions between the hydrogen atoms are responsible for the surprisingly strong aliphatic interactions. Moreover, whereas σ/σ and σ/π interactions are local, the π/π stacking are inherently delocalized, which give rise to a non‐additive effect. These new types of dispersion interactions between saturated groups can be exploited in the rational design of novel carbon materials.     

Cooperative and Diminutive Effects of Pnicogen Bonds and Cation–π Interactions

The interplay between pnicogen bonds and cation–π interactions has been investigated at the MP2/aug‐cc‐pVDZ level. Interesting cooperative and diminutive effects are observed when pnicogen bonds and cation–π interactions coexist in the same complex. These effects have been analyzed in terms of the structural, energetic, and charge‐transfer properties of the complexes. The variations in electron density at critical points of the intermolecular bond have been used to analyze bond strengthening or weakening. The nature of the interactions and the mechanisms of cooperative and diminutive effects have been studied by means of symmetry‐adapted perturbation theory and molecular electrostatic potentials.     

Interplay between Cation–π and Coinage‐Metal–Oxygen Interactions: An Ab Initio Study and Cambridge Structural Database Survey

The interplay between cation–π and coinage‐metal–oxygen interactions are investigated in the ternary systems N???PhCCM???O (N=Li⁺, Na⁺, Mg²⁺; M=Ag, Au; O=water, methanol, ethanol). A synergetic effect is observed when cation–π and coinage‐metal–oxygen interactions coexist in the same complex. The cation–π interaction in most triads has a greater enhancing effect on the coinage‐metal–oxygen interaction. This effect is analyzed in terms of the binding distance, interaction energy, and electrostatic potential in the complexes. Furthermore, the formation, strength, and nature of both the cation–π and coinage‐metal–oxygen interactions can be understood in terms of electrostatic potential and energy decomposition. In addition, experimental evidence for the coexistence of both interactions is obtained from the Cambridge Structural Database (CSD).     

The ability of Ex2Box4+ to interact with guests containing π‐electron‐rich and π‐electron‐poor moieties

The ability of Ex 2 Box⁴⁺ as a host, able to trap guests containing both π‐electron rich (polycyclic aromatic hydrocarbons‐PAHs) and π‐electron poor (quinoid‐ and nitro‐PAHs) moieties was investigated to shed light on the main factors that control the host–guest (HG) interaction. The nature of the HG interactions was elucidated by energy decomposition (EDA‐NOCV), noncovalent interaction (NCI), and magnetic response analyses. EDA‐NOCV reveals that dispersion contributions are the most significant to sustain the HG interaction, while electrostatic and orbital contributions are very tiny. In fact, no significant covalent character in the HG interactions was observed. The obtained results point strictly to NCIs, modulated by dispersion contributions. Regardless of whether the guests contain π‐electron‐rich or π‐electron‐poor moieties, and no significant charge‐transfer was observed. All in all, HG interactions between guests 3‐14 and host 2 are predominantly modulated by π‐π stacking.     

Surfing π Clouds for Noncovalent Interactions: Arenes versus Alkenes

A comparative study of molecular balances by NMR spectroscopy indicates that noncovalent functional‐group interactions with an arene dominate over those with an alkene, and that a π‐facial intramolecular hydrogen bond from a hydroxy group to an arene is favored by approximately 1.2 kJ mol^?1. The strongest interaction observed in this study was with the cyano group. Analysis of the series of groups CH₂CH₃, CH?CH₂, C?CH, and C?N shows a correlation between conformational free‐energy differences and the calculated charge on the C^α atom of these substituents, which is indicative of the electrostatic nature of their π interactions. Changes in the free‐energy differences of conformers show a linear dependence on the solvent hydrogen bond acceptor parameter β.