Potential Energy Curves for Cation−π Interactions: Off-Axis Configurations Are Also Attractive

Substituent effects on non-covalent interactions with aromatic rings: insights from computational chemistry

Non-covalent interactions with aromatic rings pervade modern chemical research. The strength and orientation of these interactions can be tuned and controlled through substituent effects. Computational studies of model complexes have provided a detailed understanding of the origin and nature of these substituent effects, and pinpointed flaws in entrenched models of these interactions in the literature. Here, we provide a brief review of efforts over the last decade to unravel the origin of substituent effects in π-stacking, XH/π, and ion/π interactions through detailed computational studies. We highlight recent progress that has been made, while also uncovering areas where future studies are warranted.

Probing eudesmane cation-π interactions in catalysis by aristolochene synthase with non-canonical amino acids

Stabilization of the reaction intermediate eudesmane cation (3) through interaction with Trp 334 during catalysis by aristolochene synthase from Penicillium roqueforti was investigated by site-directed incorporation of proteinogenic and non-canonical aromatic amino acids. The amount of germacrene A (2) generated by the mutant enzymes served as a measure of the stabilization of 3. 2 is a neutral intermediate, from which 3 is formed during PR-AS catalysis by protonation of the C6,C7 double bond. The replacement of Trp 334 with para-substituted phenylalanines of increasing electron-withdrawing properties led to a progressive accumulation of 2 that showed a good correlation with the interaction energies of simple cations such as Na(+) with substituted benzenes. These results provide compelling evidence for the stabilizing role played by Trp 334 in aristolochene synthase catalysis for the energetically demanding transformation of 2 to 3.

Quantitative assessment of substituent effects on cation-π interactions using molecular electrostatic potential topography

A molecular electrostatic potential (MESP) topography based approach has been proposed to quantify the substituent effects on cation-π interactions in complexes of mono-, di-, tri-, and hexasubstituted benzenes with Li(+), Na(+), K(+), and NH(4)(+). The MESP minimum (V(min)) on the π-region of C(6)H(5)X showed strong linear dependency to the cation-π interaction energy, E(M(+)). Further, cation-π distance correlated well with V(min)-π distance. The difference between V(min) of C(6)H(5)X and C(6)H(6) (ΔV(min)) is proposed as a good parameter to quantify the substituent effect on cation-π interaction. Compared to benzene, electron-donating groups stabilize the di-, tri-, and hexasubstituted cation-π complexes while electron-withdrawing groups destabilize them. In multiple substituted complexes, E(M(+)) is almost equal (∼95%) to the sum of the individual substituent contributions (E(M(+)) ≈ Σ(ΔE(M(+)))), suggesting that substituent effect on cation-π interactions is largely additive. The ΔV(min) of C(6)H(5)X systems and additivity feature have been used to make predictions on the interaction energies of 80 multiple substituted cation-π complexes with above 97% accuracy. The average mean absolute deviation of the V(min)-predicted interaction energy, E(M(+))(V) from the calculated E(M(+)) is -0.18 kcal/mol for Li(+), -0.09 kcal/mol for Na(+), -0.43 kcal/mol for K(+), and -0.67 kcal/mol for NH(4)(+), which emphasize the predictive power of V(min) as well as the additive feature of the substituent effect.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu

The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.