Synergy between experimental and computational studies of aromatic stacking interactions

Direct Through-Space Substituent-π Interactions in Noncovalent Arene-Fullerene Assemblies

The arene-arene interactions between electron-rich and deficient aromatics have been less understood. Herein, we focus on a [60]fullerene π-surface as an electron-deficient aromatics. Using a 1H signal of H2O@C60 as a magnetic probe, the presence of benzene-fullerene interactions was confirmed. To investigate substituent effects on the noncovalent arene-fullerene interactions, NMR titration experiments were carried out using an open-[60]fullerene and a series of substituted benzenes, i. e., PhX (X=NO2, CN, Cl, OMe, H, CH3, and NH2), demonstrating a 1 : 2 stoichiometry with a positive correlation between stabilization energies upon the first association (ΔG1) and Hammet constants (σm). The destabilization of the self-assembled structure for X=OMe with a σ-withdrawing nature clearly showed direct through-space substituent-π interactions describable by the Wheeler-Houk model while the second association was suggested to be considerably perturbed by the secondary effects.

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.

Host-guest systems for the SAMPL9 blinded prediction challenge: phenothiazine as a privileged scaffold for binding to cyclodextrins

Model systems are widely used in biology and chemistry to gain insight into more complex systems. In the field of computational chemistry, researchers use host-guest systems, relatively simple exemplars of noncovalent binding, to train and test the computational methods used in drug discovery. Indeed, host-guest systems have been developed to support the community-wide blinded SAMPL prediction challenges for over a decade. While seeking new host-guest systems for the recent SAMPL9 binding prediction challenge, which is the focus of the present PCCP Themed Collection, we identified phenothiazine as a privileged scaffold for guests of β cyclodextrin (βCD) and its derivatives. Building on this observation, we used calorimetry and NMR spectroscopy to characterize the noncovalent association of native βCD and three methylated derivatives of βCD with five phenothiazine drugs. The strongest association observed, that of thioridazine and one of the methyl derivatives, exceeds the well-known high affinity of rimantidine with βCD. Intriguingly, however, methylation of βCD at the 3 position abolished detectible binding for all of the drugs studied. The dataset has a clear pattern of entropy-enthalpy compensation. The NMR data show that all of the drugs position at least one aromatic proton at the secondary face of the CD, and most also show evidence of deep penetration of the binding site. The results of this study were used in the SAMPL9 blinded binding affinity-prediction challenge, which are detailed in accompanying papers of the present Themed Collection. These data also open the phenothiazines and, potentially, chemically similar drugs, such as the tricyclic antidepressants, as relatively potent binders of βCD, setting the stage for future SAMPL challenge datasets and for possible applications as drug reversal agents.

The effect of weak π-π interactions on single-molecule electron transport properties of the tetraphenylethene molecule and its derivatives: a first-principles study

Intramolecular π-π interactions are a significant research focus in fields such as chemistry, biology, and materials science. Different configurations of benzene-benzene moieties within a molecule can affect the magnitude of their π-π interactions, consequently influencing the electronic transport capabilities of the molecule. In this study, we designed three π-conjugated molecules, TPEM, TPEEM, and TEEPM, based on tetraphenylethene (TPE). These three molecules exhibit three distinct π-conjugated structures: linear cis-π-conjugation, linear trans-π-conjugation, and cross-π-conjugation. Thereinto, TPEM and TPEEM molecules share the same TPE core, with identical π-π interaction distances, while the TEEPM molecule has acetylene groups between the TPE units, thereby increasing the π-π interaction distances between the benzene moieties. Using density functional theory calculations combined with non-equilibrium Green's function (DFT+NEGF), our results reveal that the conductance order of different π-conjugated structures in TPEM and TPEEM molecules is as follows: cis > cross ≈ trans. Through analysis of transmission spectra, transmission pathways, and the innermost π orbitals, we find that in TPEM and TPEEM molecules, the cis- and cross-π-conjugated structures exhibit π-π interactions between benzene moieties and provide special through-space electron transport pathways, enhancing their electronic transport capabilities in coordination with the bonded molecular framework, whereas their trans-conjugated structures only allow electron transport along the molecular backbone. In contrast, in TEEPM molecule, due to the absence of π-π interactions, the conductance of different π-conjugated structures is primarily determined by the molecular backbone and follows the order: trans > cis > cross. These findings provide a theoretical basis for designing single-molecule electronic devices with multiple electron channels based on intramolecular π-π interactions.

Aromatic interactions with heterocycles in water

Conformationally well-defined supramolecular complexes that can be studied in different solvents provide a platform for separating and quantifying free energy contributions due to functional group interactions and desolvation. Here 1:1 complexes formed between four different calix[4]pyrrole receptors and eleven different pyridine N-oxide guests have been used to dissect the factors that govern aromatic interactions with heterocycles in water and in chloroform solution. 1H NMR spectroscopy shows that the three-dimensional structures of the complexes are fixed by four H-bonding interactions between the pyrrole donors at the bottom of the receptor and the N-oxide acceptor on the guest, locking the geometrical arrangement of interacting functional groups in the binding pocket at the other end of the receptor. An aromatic heterocycle on the guest makes two stacking interactions and two edge-to-face interactions with the side walls of the receptor. Chemical double mutant cycles were used to measure the free energy contribution of these four aromatic interactions to the overall stability of the complex. In chloroform, the aromatic interactions measured with pyridine, pyrimidine, furan, thiophene and thiazole are similar to the interactions with a phenyl group, but the effect of introducing a heteroatom depends on where it sits with respect to the aromatic side-walls of the cavity. A nitrogen lone pair directed into a π-face of the side-walls of the binding site leads to repulsive interactions of up to 8 kJ mol-1. In water, the heterocycle aromatic interactions are all significantly more favourable (by up to 12 kJ mol-1). For the non-polar heterocycles, furan and thiophene, the increase in interaction energy correlates directly with hydrophobicity, as measured by the free energy of transfer of the heterocycle from n-hexadecane into water (ΔG°(water-hex)). For the heterocycles with polar nitrogen H-bond acceptors, water can access cracks in the walls of the receptor binding site to solvate the edges of the heterocycles without significantly affecting the geometry of the aromatic interactions, and these nitrogen-water H-bonds stabilise the complexes by about 15 kJ mol-1. The results highlight the complexity of the solvation processes that govern molecular recognition in water.

Two are not enough: synthetic strategies and applications of unnatural base pairs

Nucleic acid chemistry is a rapidly evolving field, and the need for novel nucleotide modifications and artificial nucleotide building blocks for diagnostic and therapeutic use, material science or for studying cellular processes continues unabated. This review focusses on the development and application of unnatural base pairs as part of an expanded genetic alphabet. Not only recent developments in "nature-like" artificial base pairs are presented, but also current synthetic methods to get access to C-glycosidic nucleotides. Wide-ranging viability in synthesis is a prerequisite for the successful use of unnatural base pairs in a broader spectrum and will be discussed.

Systematic Study of Heteroarene Stacking Using a Congeneric Set of Molecular Glues for Procaspase-6

Heteroaromatic stacking interactions are important in drug binding, supramolecular chemistry, and materials science, making protein-ligand model systems of these interactions of considerable interest. Here we studied 30 congeneric ligands that each present a distinct heteroarene for stacking between tyrosine residues at the dimer interface of procaspase-6. Complex X-ray crystal structures of 10 analogs showed that stacking geometries were well conserved, while high-accuracy computations showed that heteroarene stacking energy was well correlated with predicted overall ligand binding energies. Empirically determined KD values in this system thus provide a useful measure of heteroarene stacking with tyrosine. Stacking energies are discussed in the context of torsional strain, the number and positioning of heteroatoms, tautomeric state, and coaxial orientation of heteroarene in the stack. Overall, this study provides an extensive data set of empirical and high-level computed binding energies in a versatile new protein-ligand system amenable to studies of other intermolecular interactions.

Substituent effects on aromatic interactions in water

Molecular recognition in water involves contributions due to polar functional group interactions, partial desolvation of polar and non-polar surfaces and changes in conformational flexibility, presenting a challenge for rational design and interpretation of supramolecular behaviour. Conformationally well-defined supramolecular complexes that can be studied in both water and non-polar solvents provide a platform for disentangling these contributions. Here 1 : 1 complexes formed between four different calix[4]pyrrole receptors and thirteen different pyridine N-oxide guests have been used to dissect the factors that govern substituent effects on aromatic interactions in water. H-bonding interactions between the receptor pyrrole donors and the guest N-oxide acceptor at one end of the complex lock the geometrical arrangement of a cluster of aromatic interactions at the other end of the complex, so that a phenyl group on the guest makes two edge-to-face and two stacking interactions with the four aromatic side-walls of the receptor. The thermodynamic contribution of these aromatic interactions to the overall stability of the complex was quantified by chemical double mutant cycles using isothermal titration calorimetry and ¹H NMR competition experiments. Aromatic interactions between the receptor and a phenyl group on the guest stabilise the complex by a factor of 1000, and addition of substituents to the guest phenyl group further stabilises the complex by an additional factor of up to 1000. When a nitro substituent is present on the guest phenyl group, the complex has a sub-picomolar dissociation constant (370 fM). The remarkable substituent effects observed in water for these complexes can be rationalised by comparison with the magnitude of the corresponding substituent effects measured in chloroform. In chloroform, the double mutant cycle free energy measurements of the aromatic interactions correlate well with the substituent Hammett parameters. Electron-withdrawing substituents increase the strength of the interactions by a factor of up to 20, highlighting the role of electrostatics in stabilising both the edge-to-face and stacking interactions. The enhanced substituent effects observed in water are due to entropic contributions associated with the desolvation of hydrophobic surfaces on the substituents. The flexible alkyl chains that line the open end of the binding site assist the desolvation of the non-polar π-surfaces of polar substituents, like nitro, but at the same time allow water to interact with the polar H-bond acceptor sites on the substituent. This flexibility allows polar substituents to maximise non-polar interactions with the receptor and polar interactions with the solvent, leading to remarkably high binding affinities.

A molecular descriptor of intramolecular noncovalent interaction for regulating optoelectronic properties of organic semiconductors

In recent years, intramolecular noncovalent interaction has become an important means to modulate the optoelectronic performances of organic/polymeric semiconductors. However, it lacks a deep understanding and a direct quantitative relationship among the molecular geometric structure, strength of noncovalent interaction, and optoelectronic properties in organic/polymeric semiconductors. Herein, upon systematical theoretical calculations on 56 molecules with and without noncovalent interactions (X···Y, X = O, S, Se, Te; Y = C, F, O, S, Cl), we reveal the essence of the interactions and the dependence of its strength on the molecular geometry. Importantly, a descriptor S is established as a function of several basic geometric parameters to well characterize the noncovalent interaction energy, which exhibits a good inverse correlation with the reorganization energies of the photo-excited states or electron-pumped charged states in organic/polymeric semiconductors. In particular, the experimental ¹H, ⁷⁷Se, and ¹²⁵Te NMR, the optical absorption and emission spectra, and single crystal structures of eight compounds fully confirm the theoretical predictions. This work provides a simple descriptor to characterize the strength of noncovalent intramolecular interactions, which is significant for molecular design and property prediction.

Development and Implementation of Atomically Anisotropic First-Principles Force Fields: A Benzene Case Study

π-interactions are an important motif in chemical and biochemical systems. However, due to their anisotropic electron densities and complex balance of intermolecular interactions, aromatic molecules represent an ongoing challenge for accurate and transferable force field development. Historically, ab initio force fields for aromatics have not exhibited good accuracy with respect to bulk properties or have only been used to study gas-phase dimers. Using benzene as a proof of concept, herein we show how our own ab initio MASTIFF force field incorporates an atomically anisotropic description of intermolecular interactions to yield an accurate and robust model for aromatic interactions irrespective of phase. Compared to existing models, the MASTIFF benzene force field not only is accurate for liquid phase properties but also offers transferability to the gas and solid phases. Additionally, we introduce a computationally efficient OpenMM plugin which enables customizable anisotropic intermolecular functional forms and which can be generically used in any MD simulation where a model for nonspherical atomic features is required. Overall, our results demonstrate the importance of atomic-level anisotropy in enabling next-generation ab initio force field development.

Synthesis of a Pyridone-Based Phthalimide Fleximer and Its Characterization and Supramolecular Property Evaluation

In this study, a novel pyridone-based phthalimide fleximer, that is, ethyl 5-cyano-6-(3-(1,3-dioxoisoindolin-2-yl)propoxy)-4-(3-methoxyphenyl)-2-methylnicotinate, was synthesized, and its structure was established by the single-crystal X-ray diffraction method. The supramolecular self-assembly of the titled compound through noncovalent interactions was then investigated thoroughly. The titled compound crystallized with two symmetry-independent molecules (A and B, Z' = 2). In agreement with experimental observations, our density functional theory calculations also showed that the titled compound has a flexible motif and can occur in various conformations, including molecules A and B. The investigation of the supramolecular framework revealed that the molecules are notably bound by the nonclassical C-H···O and C-H···N hydrogen bonds and C-H···π interactions. Hirshfeld surface analysis was carried out to quantify the various intermolecular interactions. The dual anti-inflammatory activity of the tilted compound was also explored by molecular docking in the active sites of 5-LOX and COX-2 receptors, which revealed good binding affinities of -9.0 and -8.6 kcal/mol, respectively.

Codelivery of π-π Stacked Dual Anticancer Drugs Based on Aloe-Derived Nanovesicles for Breast Cancer Therapy

To overcome the low efficacy of conventional monotherapeutic approaches that use a single drug, functional nanocarriers loaded with an amalgamation of anticancer drugs have been promising in cancer therapy. Herein, aloe-derived nanovesicles (gADNVs) are modified with an active integrin-targeted peptide (Arg-Gly-Asp, RGD) by the postinsertion technique to deliver indocyanine green (ICG) and doxorubicin (DOX) for efficient breast cancer therapy. We presented for the first time that the π-π stacking interaction can turn the "competitive" relationship of ICG and DOX inside gADNVs into a "cooperative" relationship and enhance their loading efficiency. The dual-drug codelivery nanosystem, denoted as DIARs, was well stable and leakproof, exhibiting high tumor-targeting capability both in vitro and in vivo. Meanwhile, this nanosystem showed significant inhibition of cell growth and migration and induced cell apoptosis with the combination of phototherapy and chemotherapy. Intravenous administration of DIARs exhibited high therapeutic efficacy in a 4T1 tumor-bearing mouse model and exhibited no obvious damage to other organs. Overall, our DIAR nanosystem constitutively integrated the natural and economical gADNVs, π-π stacking interaction based on efficient drug loading, and tumor-targeted RGD modification to achieve an effective combination therapy for breast cancer.

Complexes of Alkaline and Ammonium Cations with Dopamine and Eumelanin Precursors: Dissecting the Role of Noncovalent Cation-π and Cation-Lone Pair (σ-Type) Interactions

Cation-π interactions and their possible competition with other noncovalent interactions (NCI) might play a key role in both dopamine- and eumelanin-based bioinspired materials. In this contribution, to unravel the delicate interplay between cation-π interactions and other possible competing forces, the configurational space of noncovalent complexes formed by dopamine or eumelanin precursors (o-benzoquinone, DHI and a semiquinone dimer) and three different cations (Na⁺, K⁺, and NH₄⁺) is sampled by means of accurate ab initio calculations. To this end, we resort to the mp2^mod method, recently validated by us for benzene-, phenol-, and catechol-cation complexes, whose computational convenience allows for an extensive exploration of the cation-molecule interaction energy surface, by sampling a total of more than 10⁴ arrangements. The mp2^mod interaction energy landscapes reveal that, besides the expected cation-π driven arrangements, for all considered molecule-cation pairs the most stable complexes are found when the cation lies within the plane containing the six-membered ring, thus maximizing the σ-type interaction with the oxygen's lone pairs. Due to the loss of aromaticity, the σ-type/cation-π strength ratio is remarkably large in o-benzoquinone, where cation-π complexes seem unlikely to be formed. The above features are shared among all considered cations but are significantly larger when considering the smaller Na⁺. Besides delivering a deeper insight onto the NCI network established by the considered precursors in the presence of ions, the present results can serve as a reference database to validate or refine lower level methods, as, for instance, the force fields employed in classical simulations.

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.

Augmented Self-Association by Electrostatic Forces in Thienopyrrole-Fused Thiadiazoles that Contain an Ester instead of an Ether Linker

During the self-assembly of π-conjugated molecules, linkers and substituents can potentially add supportive noncovalent intermolecular interactions to π-stacking interactions. Here, we report the self-assembly behavior of thienopyrrole-fused thiadiazole (TPT) fluorescent dyes that possess ester or ether linkers and dodecyloxy side chains in solution and the condensed phase. A comparison of the self-association behavior of the ester- and ether-bridged compounds in solution using detailed UV-vis, fluorescence, and NMR spectroscopic studies revealed that the subtle replacement of the ether linkers by ester linkers leads to a distinct increase in the association constant (ca. 3-4 fold) and the enthalpic contribution (ca. 3 kcal mol-1). Theoretical calculations suggest that the ester linkers, which are in close proximity to one another due to the π-stacking interactions, induce attractive electrostatic forces and augment self-association. The self-assembly of TPT dyes into well-defined 1D clusters with high aspect ratios was observed, and their morphologies and crystallinity were investigated using SEM and X-ray diffraction analyses. TPTs with ester linkers exhibit a columnar liquid crystalline mesophase in the condensed phase.

Dissection of the Polar and Non-Polar Contributions to Aromatic Stacking Interactions in Solution

Aromatic stacking interactions have been a matter of study and debate due to their crucial role in chemical and biological systems. The strong dependence on orientation and solvent together with the relatively small interaction energies have made evaluation and rationalization a challenge for experimental and theoretical chemists. We have used a supramolecular cage formed by two tris(pyridylmethyl)amines units to build chemical Double Mutant Cycles (DMC) for the experimental measurement of the free energies of π-stacking interactions. Extrapolating the substituent effects to remove the contribution due to electrostatic interactions reveals that there is a substantial contribution to the measured stacking interaction energies which is due to non-polar interactions (-3 to -6 kJ mol-1 ). The perfectly flat nature of the surface of an aromatic ring gives π-stacking an inherent advantage over non-polar interactions with alkyl groups and accounts for the wide-spread prevalence of stacking interactions in Nature.

Arene-Perfluoroarene Interactions in Solution

A systematic study of arene-perfluoroarene interactions in solution is presented. Using a combination of NMR titration experiments, X-ray crystallography, and computational analysis, we analyze the effects of fluorination, substituents, ring size, and solvation on the arene-perfluoroarene interaction. We find that fluorination, extension of the π systems, and enhancement of solvent polarity greatly stabilize the stacking energy up to 3 orders of magnitude (K_a = <1 to 6000 M^-1), with the highest K_a achieved for the interaction of water-soluble variants of perfluoronaphthalene and anthracene in buffered D₂O (pD = 12). Combining computational and experimental results, we conclude that this impressive binding energy is a result of enthalpically favorable electrostatic and dispersion interactions as well as the entropically driven hydrophobic effect. The enhanced understanding of arene-perfluoroarene interactions in aqueous solution sets the stage for the implementation of this abiotic intermolecular interaction in biology and medicine.

Anions effect construction of 1D naphthalene diimide supramolecular chains by π interactions and fluorescence detect iodide anion

Five diverse 1D supramolecular chains, {[4-pmntd]}_n(1), {H₂[4-pmntd]•2Br^-}_n(2), {H₂[4-pmntd]•2NO₃^-}_n(3), {H₂[4-pmntd]•2ClO₄^-}_n(4), {H₂[4-pmntd]•2BF₄^-}_n(5), (where 4-pmntd was N,N'-bis (4-pyridylmethyl)naphthalene diimide) were synthesized and characterized by X-ray single-crystal structure analysis, IR spectroscopy, elemental analyses, thermogravimetric analyses, fluorescence detection. The anions effect construction of their 1D chain structural diversity through different π interactions. Compound 1 through the adjacent pyridine rings parallel π∙∙∙π interactions formed 1D linear chain structure. Compound 2 through Br^- anion∙∙∙π interactions and halogenbond interactions formed 1D zigzag chain structure. Compound 3 through lone pair∙∙∙π interactions of naphthalene diimide and the adjacent carboxyl group formed 1D stairway chain structure. Compound 4 through ClO₄^- anion∙∙∙π interactions formed 1D ribbon chain structure. Compound 5 through parallel π∙∙∙π interactions of the adjacent naphthalene diimide planes and pyridine rings formed 1D ladder chain structure. The five compounds' fluorescence properties and thermal stabilities were investigated. The compound 2 solution could fluorescence detection for iodide anion via fluorescence quenching.

N-Arylimide Molecular Balances: A Comprehensive Platform for Studying Aromatic Interactions in Solution

Noncovalent interactions of aromatic surfaces play a key role in many biological processes and in determining the properties and utility of synthetic materials, sensors, and catalysts. However, the study of aromatic interactions has been challenging because these interactions are usually very weak and their trends are modulated by many factors such as structural, electronic, steric, and solvent effects. Recently, N-arylimide molecular balances have emerged as highly versatile and effective platforms for studying aromatic interactions in solution. These molecular balances can accurately measure weak noncovalent interactions in solution via their influence on the folded-unfolded conformational equilibrium. The structure (i.e., size, shape, π-conjugation, and substitution) and nature (i.e., element, charge, and polarity) of the π-surfaces and interacting groups can be readily varied, enabling the study of a wide range of aromatic interactions. These include aromatic stacking, heterocyclic aromatic stacking, and alkyl-π, chalcogen-π, silver-π, halogen-π, substituent-π, and solvent-π interactions. The ability to measure a diverse array of aromatic interactions within a single model system provides a unique perspective and insights as the interaction energies, stability trends, and solvent effects for different types of interactions can be directly compared. Some broad conclusions that have emerged from this comprehensive analysis include: (1) The strongest aromatic interactions involve groups with positive charges such as pyridinium and metal ions which interact with the electrostatically negative π-face of the aromatic surface via cation-π or metal-π interactions. Attractive electrostatic interactions can also form between aromatic surfaces and groups with partial positive charges. (2) Electrostatic interactions involving aromatic surfaces can be switched from repulsive to attractive using electron-withdrawing substituents or heterocycles. These electrostatic trends appear to span many types of aromatic interactions involving a polar group interacting with a π-surface such as halogen-π, chalcogen-π, and carbonyl-π. (3) Nonpolar groups form weak but measurable stabilizing interactions with aromatic surfaces in organic solvents due to favorable dispersion and/or solvophobic effects. A good predictor of the interaction strength is provided by the change in solvent-accessible surface area. (4) Solvent effects modulate the aromatic interactions in the forms of solvophobic effects and competitive solvation, which can be modeled using solvent cohesion density and specific solvent-solute interactions.

Discrete π Stack of a Tweezer-Shaped Naphthalenediimide-Anthracene Conjugate

The design and synthesis of a tweezer-shaped naphthalenediimide (NDI)-anthracene conjugate (2NDI) are reported. In the structure of the closed form (π_NDI ⋅⋅⋅π_NDI stack) of 2NDI, which was elucidated by single-crystal XRD, the existence of C-H⋅⋅⋅O hydrogen bonding involving the nearest carbonyl oxygen atom of an NDI unit was suggested. The tunability of π_NDI ⋅⋅⋅π_NDI interactions was studied by means of UV/Vis absorption, fluorescence and NMR spectroscopy and molecular modelling. This revealed that the π_NDI ⋅⋅⋅π_NDI interactions in 2NDI affect the absorption and emission properties depending on the temperature. Furthermore, in polar solvents, 2NDI prefers the stronger π_NDI ⋅⋅⋅π_NDI stack, whereas the π_NDI ⋅⋅⋅π_NDI interaction is diminished in nonpolar solvents. Importantly, the conformational variations of 2NDI can be reversibly switched by variation in temperature, and this suggests potential application for fluorogenic molecular switches upon temperature changes.

Hydrogen-Mediated Noncovalent Interactions in Solids: What Can NMR Crystallography Tell About?

Hydrogen atoms play a crucial role in the aggregation of organic (bio)molecules through diverse number of noncovalent interactions that they mediate, such as electrostatic in proton transfer systems, hydrogen bonding, and CH-π interactions, to mention only the most prominent. To identify and adequately describe such low-energy interactions, increasingly sensitive methods have been developed over time, among which quantum chemical computations have witnessed impressive advances in recent years. For reaching the present state-of-the-art, computations had to rely on a pool of relevant experimental data, needed at least for validation, if not also for other purposes. In the case of molecular crystals, the best illustration for the synergy between computations and experiment is given by the so-called NMR crystallography approach. Originally designed to increase the confidence level in crystal structure determination of organic compounds from powders, NMR crystallography is able now to offer also a wealth of information regarding the noncovalent interactions that drive molecules to pack in a given crystalline pattern or another. This is particularly true for the noncovalent interactions which depend on the exact location of labile hydrogen atoms in the system: in such cases, NMR crystallography represents a valuable characterization tool, in some cases complementing even the standard single-crystal X-ray diffraction technique. A concise introduction in the field is made in this mini-review, which is aimed at providing a comprehensive picture with respect to the current accuracy level reached by NMR crystallography in the characterization of hydrogen-mediated noncovalent interactions in organic solids. Different types of practical applications are illustrated with the example of molecular crystals studied by our research group, but references to other representative developments reported in the literature are also made. By summarizing the major concepts and methodological progresses, the present work is also intended to be a guide to the practical potential of this relatively recent analytical tool for the scientists working in areas where crystal engineering represents the main approach for rational design of novel materials.

Physical Mechanisms Governing Substituent Effects on Arene-Arene Interactions in a Protein Milieu

Arene-arene interactions play important roles in protein-ligand complex formation. Here, we investigate the characteristics of arene-arene interactions between small organic molecules and aromatic amino acids in protein interiors. The study is based on X-ray crystallographic data and quantum mechanical calculations using the enzyme acetylcholinesterase and selected inhibitory ligands as a model system. It is shown that the arene substituents of the inhibitors dictate the strength of the interaction and the geometry of the resulting complexes. Importantly, the calculated interaction energies correlate well with the measured inhibitor potency. Non-hydrogen substituents strengthened all interaction types in the protein milieu, in keeping with results for benzene dimer model systems. The interaction energies were dispersion-dominated, but substituents that induced local dipole moments increased the electrostatic contribution and thus yielded more strongly bound complexes. These findings provide fundamental insights into the physical mechanisms governing arene-arene interactions in the protein milieu and thus into molecular recognition between proteins and small molecules.

The Solvent Effect on Weak Interactions in Supramolecular Polymers: Differences between Small Molecular Probes and Supramolecular Polymers

In this minireview, weak interactions that occur in supramolecular polymers are discussed. Combination of weak and strong interactions plays an important role in the construction of supramolecular polymers. It is beneficial to separate the contributions of the weak interactions, as well as each solvent effect on the weak interactions. However, it is generally difficult to observe each solvent effect separately at work in each interaction. Small molecular probes are useful to estimate the contributions of the weak interaction. But, the results should be treated with caution when applied to supramolecular polymer systems. To overcome the problems, a new solvent parameter, solvation ability (SA), is introduced, which was determined on the balance point of extended and stacked forms of porphyrin-based interconvertible supramolecular polymers.

Mechanistic Study of the Solvent-Dependent Formation of Extended and Stacked Supramolecular Polymers Composed of Bis(imidazolylporphyrinatozinc) Molecules

Bis(imidazolylporphyrinatozinc) molecules linked through a 1,3-butadiynylene moiety respond to the solvents they are dissolved in to afford exclusively extended (E) or stacked (S) supramolecular polymers. This system is expected to be a solvation/desolvation indicator. However, the principles underlying the solvent-dependent formation of the two types of polymers and the mechanism of the transformation between them are unclear. The formation of the polymers is considered to depend on the two types of complementary coordination bonds that can be formed and the π-π interactions between the porphyrins. In this study, the contributions and solvent dependence of both the coordination bonds and the π-π interactions have been investigated. The results clearly indicate that the coordination bonds are weakly or little solvent-dependent, and that the π-π interactions function effectively only in the inner porphyrins of the S-polymer and are strongly solvent-dependent. Thermodynamic analysis revealed that the formation of the E- or S-polymer in solution is determined by the total energies and the type of solvent used. The transformation of the E- to S-polymer was investigated by gel permeation chromatography. The kinetics of the transformation were also determined. The role of the terminal imidazolylporphyrinatozinc moieties was also investigated: The results indicate that the transformation from the E- to S-polymer occurs by an exchange mechanism between the polymers, induced by attack of terminal free imidazolyl groups on a polymer to zinc porphyrins on other polymers.

Photooxidative Generation of Dodecaborate-Based Weakly Coordinating Anions

Redox-active proanions of the type B₁₂(OCH₂Ar)₁₂ [Ar = C₆F₅ (1), 4-CF₃C₆H₄ (2), 3,5-(CF₃)₂C₆H₃ (3)] are introduced in the context of an experimental and computational study of the visible-light-initiated polymerization of a family of styrenes. Neutral, air-stable proanions 1-3 were found to initiate styrene polymerization through single-electron oxidation under blue-light irradiation, resulting in polymers with number-average molecular weights (M_n) ranging from ∼6 to 100 kDa. Shorter polymer products were observed in the majority of experiments, except in the case of monomers containing 4-X (X = F, Cl, Br) substituents on the styrene monomer when polymerized in the presence of 1 in CH₂Cl₂. Only under these specific conditions are longer polymers (>100 kDa) observed, strongly supporting the formulation that reaction conditions significantly modulate the degree of ion pairing between the dodecaborate anion and cationic chain end. This also suggests that 1-3 behave as weakly coordinating anions (WCA) upon one-electron reduction because no incorporation of the cluster-based photoinitiators is observed in the polymeric products analyzed. Overall, this work is a conceptual realization of a single reagent that can serve as a strong photooxidant, subsequently forming a WCA.

Supramolecular cage encapsulation as a versatile tool for the experimental quantification of aromatic stacking interactions

The widespread presence of aromatic stacking interactions in chemical and biological systems, combined with their relatively small energetic contribution, have led to a plethora of theoretical and experimental studies for their quantification and rationalization. Typically, π-π aromatic interactions are studied as a function of substituents to gather information about the interaction mechanism. While experiments suggest that aromatic interactions are dominated by local electrostatic contacts between π-electron density and CH groups, theoretical work has raised the possibility that direct electrostatic interactions between local dipoles of the substituents may play a role. We describe a supramolecular cage that binds two aromatic carboxylates in a stacked geometry such that the aromatic substituents are remote in space. Chemical Double Mutant Cycles (DMCs) were used to measure fifteen different aromatic stacking interactions as a function of substituent (NMe2, OMe, Me, Cl and NO2). When both aromatic rings have electron-withdrawing nitro substituents, the interaction is attractive (-2.8 kJ mol-1) due to reduced π-electron repulsion. When both aromatic rings have electron-donating di-methylamino substituents, the interaction is repulsive (+2.0 kJ mol-1) due to increased π-electron repulsion. The results show that aromatic stacking interactions are dominated by short range electrostatic contacts rather than substituent dipole interactions.

Study of through-space substituent–π interactions using N-phenylimide molecular balances

Substituent–π interactions associated with aromatic stacking interactions were experimentally measured using a small N-phenylimide molecular balance model system.

Density Functional Study on the Adsorption of 5-Membered N-Heterocycles on B/N/BN-Doped Graphene: Coronene as a Model System

Adsorption of seven 5-membered N-heterocycles on B/N/BN-doped graphene (with coronene as a model system) has been studied using density functional theory (DFT). The geometry of the complexes validated the involvement of both π···π stacking and N-H···π interaction in the adsorption process. The stability of the complexes is measured in terms of stabilization energy, and the results suggested that the complexes are stable enough (stabilization energies are in the range of 7.61-14.77 kcal mol^-1). Studies confirmed the stability of complexes in the solvent phase too irrespective of the dielectric of the solvent. Dispersive force is the major mode of interaction in stabilizing the complexes. Natural bond orbital analysis indicated a small contribution from electrostatic and covalent interactions. Thermochemical analysis revealed that the complexation is exothermic in nature and favorable at a lower temperature. Adsorption of N-heterocycles exerts a nominal impact on the electronic properties of the undoped/doped graphene. The study presents a simple approach to introduce an arbitrary functionality to undoped/doped graphene by preserving its electronic properties.

Size-dependent rate acceleration in the silylation of secondary alcohols: the bigger the faster

Relative rates for the reaction of secondary alcohols carrying large aromatic moieties with silyl chlorides carrying equally large substituents have been determined in organic solvents. Introducing thoroughly matching pairs of big dispersion energy donor (DED) groups enhanced rate constants up to four times, notably depending on the hydrogen bond donor ability of the solvent. A linear correlation between computed dispersion energy contributions to the stability of the silyl ether products and experimental relative rate constants was found. These results indicate a cooperation between solvophobic effects and DED-groups in the kinetic control of silylation reactions.

Synthesis of New 5-Aryl-benzo[f][1,7]naphthyridines via a Cascade Process (Ugi-3CR/Intramolecular Aza-Diels-Alder Cycloaddition)/Aromatization

A series of eight new 5-aryl-benzo[f][1,7]naphthyridines were synthesized in 17 to 64% overall yields via an improved MW-assisted cascade-like one pot process (Ugi–three component reaction/intramolecular aza-Diels-Alder cycloaddition) coupled to an aromatization process from tri-functional dienophile-containing ester-anilines, substituted benzaldehydes and the chain-ring tautomerizable 2-isocyano-1-morpholino-3-phenylpropan-1-one as starting reagents, under mild conditions. The doubly activated dienophile and the aza-diene functionalities of the eight new Ugi-adducts were exploited to perform an in situ aza-Diels-Alder cycloaddition/aromatization (dehydration/oxidation) process, toward the complex polysubstituted 5-aryl-polyheterocycles, which could be taken as starting point for further SAR studies because the benzo[f][1,7]naphthyridine is the core of various bioactive products. It is relevant to emphasize that the synthesis or isolation of benzo[f][1,7]naphthyridines containing a substituted aromatic ring in the C-5 position, has not been published before.

Reversible mechanofluorochromism of aniline-terminated phenylene ethynylenes

Seven three-ring phenylene-ethynylene (PE) structural analogs, differing only in the lengths of alkyl chains on terminal aniline substituents, show 50-62 nm bathochromic shifts in emission maxima in response to mechanical force (mechanofluorochromism, MC). These shifts are fully reversible with heat or solvent fuming. Shearing of these solids yields a transition from green-emitting crystalline phases to orange-emitting amorphous phases as established by differential scanning calorimetry and X-ray diffraction. Molecules with shorter alkyl chain lengths required higher temperatures to recover the hypsochromically shifted crystalline phases after grinding, while the recovery with chain lengths longer than butyl occurred at room temperature. In addition to this structure-dependent thermochromism, these compounds retain their MC properties in polymer hosts to various extents. The crystalline phases of these materials have PE chromophores that are twisted due to non-covalent perfluoroarene-arene (ArF-ArH) interactions involving perfluorophenyl pendants and the terminal rings of the PE chromophore, resulting in interrupted conjugation and an absence of chromophore aggregation. The MC behavior of an analog without the perfluoroarene rings is severely attenuated. This work demonstrates the general utility of twisted PEs as stimuli-responsive moieties and reveals clear structure-property relationships regarding the effects of alkyl chain length on these materials.

Guest Back-Folding: A Molecular Design Strategy That Produces a Deep-Red Fluorescent Host/Guest Pair with Picomolar Affinity in Water

One of the major goals of modern supramolecular chemistry, with important practical relevance in many technical fields, is the development of synthetic host/guest partners with ultrahigh affinity and selectivity in water. Currently, most association pairs exhibit micromolar affinity or weaker, and there are very few host/guest systems with K_a > 10⁹ M^-1, apparently due to a barrier imposed by enthalpy/entropy compensation. This present study investigated the threading of a water-soluble tetralactam cyclophane by a deep-red fluorescent squaraine guest with flanking polyethylene glycol chains, an association process that is dominated by a highly favorable enthalpic driving force. A squaraine structure was rationally designed to permit guest back-folding as a strategy to greatly expand the hydrophobic surface area that could be buried upon complexation. Guided by computational modeling, an increasing number of N-benzyl groups were appended to the squaraine core, so that, after threading, the aromatic rings could fold back and stack against the cyclophane periphery. The final design iteration exhibited an impressive combination of fluorescence and supramolecular properties, including ratiometric change in deep-red emission, picomolar affinity ( K_a = 5.1 × 10¹⁰ M^-1), and very rapid threading ( k_on = 7.9 × 10⁷ M^-1 s^-1) in water at 25 °C. Similar excellent behavior was observed in serum solution. A tangible outcome of this study is a new cyclophane/squaraine association pair that will be a versatile platform for many different types of fluorescence-based imaging and diagnostics applications. From a broader perspective, guest back-folding of aromatic groups is a promising new supramolecular stabilization strategy to overcome enthalpy/entropy compensation and produce ultrahigh affinity [2]pseudorotaxane complexes in water and biological media.

Noncovalent Interactions in the Catechol Dimer

Noncovalent interactions play a significant role in a wide variety of biological processes and bio-inspired species. It is, therefore, important to have at hand suitable computational methods for their investigation. In this paper, we report on the contribution of dispersion and hydrogen bonds in both stacked and T-shaped catechol dimers, with the aim of delineating the respective role of these classes of interactions in determining the most stable structure. By using second-order Møller⁻Plesset (MP2) calculations with a small basis set, specifically optimized for these species, we have explored a number of significant sections of the interaction potential energy surface and found the most stable structures for the dimer, in good agreement with the highly accurate, but computationally more expensive coupled cluster single and double excitation and the perturbative triples (CCSD(T))/CBS) method.