In chemistry, pi stacking (also called π–π stacking) refers to the presumptive attractive, noncovalent pi interactions (orbital overlap) between the pi bonds of aromatic rings. According to some authors direct stacking of aromatic rings (the "sandwich interaction") is electrostatically repulsive.
What is more commonly observed (see figure to the right) is either a staggered stacking (parallel displaced) or pi-teeing (perpendicular T-shaped) interaction both of which are electrostatic attractive[1] For example, the most commonly observed interactions between aromatic rings of amino acid residues in proteins is a staggered stacked followed by a perpendicular orientation. Sandwiched orientations are relatively rare.[2]
Pi stacking is repulsive as it places carbon atoms with partial negative charges from one ring on top of other partial negatively charged carbon atoms from the second ring and hydrogen atoms with partial positive charges on top of other hydrogen atoms that likewise carry partial positive charges. In staggered stacking, one of the two aromatic rings is offset sideways so that the carbon atoms with partial negative charge in the first ring are placed above hydrogen atoms with partial positive charge in the second ring so that the electrostatic interactions become attractive. Likewise, pi-teeing interactions in which the two rings are oriented perpendicular to either other is electrostatically attractive as it places partial positively charged hydrogen atoms in close proximity to partially negatively charged carbon atoms. An alternative explanation for the preference for staggered stacking is due to the balance between van der Waals interactions (attractive dispersion plus Pauli repulsion).[3]
These staggered stacking and π-teeing interactions between aromatic rings are important in nucleobase stacking within DNA and RNA molecules, protein folding, template-directed synthesis, materials science, and molecular recognition. Despite the wide use of term pi stacking in the scientific literature, there is no theoretical justification for its use.[4]
The benzene dimer is the prototypical system for the study of pi stacking, and is experimentally bound by 8–12 kJ/mol (2–3 kcal/mol) in the gas phase with a separation of 4.96 Å between the centers of mass for the T-shaped dimer. The small binding energy makes the benzene dimer difficult to study experimentally, and the dimer itself is only stable at low temperatures and is prone to cluster.[5]
Other evidence against pi stacking comes from X-ray crystallography. Perpendicular and offset parallel configurations can be observed in the crystal structures of many simple aromatic compounds. Similar offset parallel or perpendicular geometries were observed in a survey of high-resolution x-ray protein crystal structures in the Protein Data Bank.[6] Analysis of the aromatic amino acids phenylalanine, tyrosine, histidine, and tryptophan indicates that dimers of these side chains have many possible stabilizing interactions at distances larger than the average van der Waals radii.
The preferred geometries of the benzene dimer have been modeled at a high level of theory with MP2-R12/A computations and very large counterpoise-corrected aug-cc-PVTZ basis sets. The two most stable conformations are the parallel displaced and T-shaped, which are essentially isoenergetic. In contrast, the sandwich configuration maximizes overlap of the pi system, which destabilizes the interaction. The sandwich configuration represents an energetic saddle point, which is consistent with the relative rarity of this configuration in x-ray crystal data.
The relative binding energies of these three geometric configurations of the benzene dimer can be explained by a balance of quadrupole/quadrupole and London dispersion forces. While benzene does not have a dipole moment, it has a strong quadrupole moment.[7] The local C–H dipole means that there is positive charge on the atoms in the ring and a correspondingly negative charge representing an electron cloud above and below the ring. The quadrupole moment is reversed for hexafluorobenzene due to the electronegativity of fluorine. The benzene dimer in the sandwich configuration is stabilized by London dispersion forces but destabilized by repulsive quadrupole/quadrupole interactions. By offsetting one of the benzene rings, the parallel displaced configuration reduces these repulsive interactions and is stabilized. The large polarizability of aromatic rings lead to dispersive interactions as major contribution to stacking effects. These play a major role for interactions of nucleobases e.g. in DNA.[8] The T-shaped configuration enjoys favorable quadrupole/quadrupole interactions, as the positive quadrupole of one benzene ring interacts with the negative quadrupole of the other. The benzene rings are furthest apart in this configuration, so the favorable quadrupole/quadrupole interactions evidently compensate for diminished dispersion forces.
The ability to fine-tune pi stacking interactions would be useful in numerous synthetic efforts. One example would be to increase the binding affinity of a small-molecule inhibitor to an enzyme pocket containing aromatic residues. The effects of heteroatoms and substituents on pi stacking interactions is difficult to model and a matter of debate.
An early model for the role of substituents in pi stacking interactions was proposed by Hunter and Sanders.[9] They used a simple mathematical model based on sigma and pi atomic charges, relative orientations, and van der Waals interactions to qualitatively determine that electrostatics are dominant in substituent effects. According to their model, electron-withdrawing groups reduce the negative quadrupole of the aromatic ring and thereby favor parallel displaced and sandwich conformations. Contrastingly, electron donating groups increase the negative quadrupole, which may increase the interaction strength in a T-shaped configuration with the proper geometry. Based on this model, the authors proposed a set of rules governing pi stacking interactions which prevailed until more sophisticated computations were applied.
Experimental evidence for the Hunter–Sanders model was provided by Siegel et al. using a series of substituted syn- and .[10] In these compounds the aryl groups "face-off" in a stacked geometry due to steric crowding, and the barrier to epimerization was measured by nuclear magnetic resonance spectroscopy. The authors reported that aryl rings with electron-withdrawing substituents had higher barriers to rotation. The interpretation of this result was that these groups reduced the electron density of the aromatic rings, allowing more favorable sandwich pi stacking interactions and thus a higher barrier. In other words, the electron-withdrawing groups resulted in "less unfavorable" electrostatic interactions in the ground state.
Hunter et al. applied a more sophisticated chemical double mutant cycle with a hydrogen-bonded "zipper" to the issue of substituent effects in pi stacking interactions.[11] This technique has been used to study a multitude of noncovalent interactions. The single mutation, in this case changing a substituent on an aromatic ring, results in secondary effects such as a change in hydrogen bond strength. The double mutation quantifies these secondary interactions, such that even a weak interaction of interest can be dissected from the array. Their results indicate that more electron-withdrawing substituents have less repulsive pi stacking interactions. Correspondingly, this trend was exactly inverted for interactions with pentafluorophenylbenzene, which has a quadrupole moment equal in magnitude but opposite in sign as that of benzene. The findings provide direct evidence for the Hunter–Sanders model. However, the stacking interactions measured using the double mutant method were surprisingly small, and the authors note that the values may not be transferable to other systems.
In a follow-up study, Hunter et al. verified to a first approximation that the interaction energies of the interacting aromatic rings in a double mutant cycle are dominated by electrostatic effects.[12] However, the authors note that direct interactions with the ring substituents, discussed below, also make important contributions. Indeed, the interplay of these two factors may result in the complicated substituent- and geometry-dependent behavior of pi stacking interactions.
The Hunter–Sanders model has been criticized by numerous research groups offering contradictory experimental and computational evidence of pi stacking interactions that are not governed primarily by electrostatic effects.[13]
The clearest experimental evidence against electrostatic substituent effects was reported by Rashkin and Waters.[14] They used meta- and para-substituted N-benzyl-2-(2-fluorophenyl)-pyridinium bromides, which stack in a parallel displaced conformation, as a model system for pi stacking interactions. In their system, a methylene linker prohibits favorable T-shaped interactions. As in previous models, the relative strength of pi stacking interactions was measured by NMR as the rate of rotation about the biaryl bond, as pi stacking interactions are disrupted in the transition state. Para-substituted rings had small rotational barriers which increased with increasingly electron-withdrawing groups, consistent with prior findings. However, meta-substituted rings had much larger barriers of rotation despite having nearly identical electron densities in the aromatic ring. The authors explain this discrepancy as direct interaction of the edge of hydrogen atoms of one ring with the electronegative substituents on the other ring. This claim is supported by chemical shift data of the proton in question.
Much of the detailed analyses of the relative contributions of factors in pi stacking have been borne out by computation. Sherill and Sinnokrot reported a surprising finding using high-level theory that all substituted benzene dimers have more favorable binding interactions than a benzene dimer in the sandwich configuration.[15] Later computational work from the Sherill group revealed that the substituent effects for the sandwich configuration are additive, which points to a strong influence of dispersion forces and direct interactions between substituents.[16] It was noted that interactions between substituted benzenes in the T-shaped configuration were more complex. Finally, Sherill and Sinnokrot argue in their review article that any semblance of a trend based on electron donating or withdrawing substituents can be explained by exchange-repulsion and dispersion terms.[17]
Houk and Wheeler also provide compelling computational evidence for the importance of direct interaction in pi stacking.[18] In their analysis of substituted benzene dimers in a sandwich conformation, they were able to recapitulate their findings using an exceedingly simple model where the substituted benzene, Ph–X, was replaced by H–X. Remarkably, this crude model resulted in the same trend in relative interaction energies, and correlated strongly with the values calculated for Ph–X. This finding suggests that substituent effects in the benzene dimer are due to direct interaction of the substituent with the aromatic ring, and that the pi system of the substituted benzene is not involved. This latter point is expanded upon below.
In summary, it would seem that the relative contributions of electrostatics, dispersion, and direct interactions to the substituent effects seen in pi stacking interactions are highly dependent on geometry and experimental design. The lack of consensus on the matter may simply reflect the complexity of the issue.
The conventional understanding of pi stacking involves quadrupole interactions between delocalized electrons in p-orbitals. In other words, aromaticity should be required for this interaction to occur. However, several groups have provided contrary evidence, calling into question whether pi stacking is a unique phenomenon or whether it extends to other neutral, closed-shell molecules.
In an experiment not dissimilar from others mentioned above, Paliwal and coauthors constructed a molecular torsion balance from an aryl ester with two conformational states.[19] The folded state had a well-defined pi stacking interaction with a T-shaped geometry, whereas the unfolded state had no aryl–aryl interactions. The NMR chemical shifts of the two conformations were distinct and could be used to determine the ratio of the two states, which was interpreted as a measure of intramolecular forces. The authors report that a preference for the folded state is not unique to aryl esters. For example, the cyclohexyl ester favored the folded state more so than the phenyl ester, and the tert-butyl ester favored the folded state by a preference greater than that shown by any aryl ester. This suggests that aromaticity is not a strict requirement for favorable interaction with an aromatic ring.
Other evidence for non-aromatic pi stacking interactions results include critical studies in theoretical chemistry, explaining the underlying mechanisms of empirical observations. Grimme reported that the interaction energies of smaller dimers consisting of one or two rings are very similar for both aromatic and saturated compounds.[20] This finding is of particular relevance to biology, and suggests that the contribution of pi systems to phenomena such as stacked nucleobases may be overestimated. However, it was shown that an increased stabilizing interaction is seen for large aromatic dimers. As previously noted, this interaction energy is highly dependent on geometry. Indeed, large aromatic dimers are only stabilized relative to their saturated counterparts in a sandwich geometry, while their energies are similar in a T-shaped interaction.
A more direct approach to modeling the role of aromaticity was taken by Bloom and Wheeler.[21] The authors compared the interactions between benzene and either 2-methylnaphthalene or its non-aromatic isomer, 2-methylene-2,3-dihydronaphthalene. The latter compound provides a means of conserving the number of p-electrons while, however, removing the effects of delocalization. Surprisingly, the interaction energies with benzene are higher for the non-aromatic compound, suggesting that pi-bond localization is favorable in pi stacking interactions. The authors also considered a homodesmotic dissection of benzene into ethylene and 1,3-butadiene and compared these interactions in a sandwich with benzene. Their calculation indicates that the interaction energy between benzene and homodesmotic benzene is higher than that of a benzene dimer in both sandwich and parallel displaced conformations, again highlighting the favorability of localized pi-bond interactions. These results strongly suggest that aromaticity is not required for pi stacking interactions in this model.
Even in light of this evidence, Grimme concludes that pi stacking does indeed exist.[20] However, he cautions that smaller rings, particularly those in T-shaped conformations, do not behave significantly differently from their saturated counterparts, and that the term should be specified for larger rings in stacked conformations which do seem to exhibit a cooperative pi electron effect.
One demonstration of stacking is found in the buckycatcher. This molecular tweezer is based on two concave buckybowls with a perfect fit for one convex fullerene molecule. Complexation takes place simply by evaporating a toluene solution containing both compounds. In solution an association constant of 8600 M−1 is measured based on changes in NMR chemical shifts.
Pi stacking is prevalent in protein crystal structures, and also contributes to the interactions between small molecules and proteins. As a result, pi–pi and cation–pi interactions are important factors in rational drug design.[22] One example is the FDA-approved acetylcholinesterase (AChE) inhibitor tacrine which is used in the treatment of Alzheimer's disease. Tacrine is proposed to have a pi stacking interaction with the indolic ring of Trp84, and this interaction has been exploited in the rational design of novel AChE inhibitors.[23]
Several variants of pi coordinated phenyls have even been tested using transition metals for stacking η6-phenyltropanes, using cyclopentadienyl and tricarbonyl in place of a benzene. Which in the case of the tricarbonyl doubled the compound's affinity for its intended ligand site (posited as due to resultant electrostatic influences being more conducive to the target).[24]
π systems are building blocks in supramolecular assembly because they often engage in noncovalent interactions. An example of π–π interactions in supramolecular assembly is the synthesis of catenane. The major challenge for the synthesis of catenane is to interlock molecules in a controlled fashion. Stoddart and co-workers developed a series of systems utilizing the strong π–π interactions between electron-rich benzene derivatives and electron-poor pyridinium rings.[25] [2]Catanene was synthesized by reacting bis(pyridinium) (A), bisparaphenylene-34-crown-10 (B), and 1, 4-bis(bromomethyl)benzene (C) (Fig. 2). The π–π interaction between A and B directed the formation of an interlocked template intermediate that was further cyclized by substitution reaction with compound C to generate the [2]catenane product.