Understanding the fundamental role of π/π, σ/σ, and σ/π dispersion interactions in shaping carbon-based materials

Understanding Emergent Complexity from a Single-Molecule Perspective

Molecules, with structural, scaling, and interaction diversities, are crucial for the emergence of complex behaviors. Interactions are essential prerequisites for complex systems to exhibit emergent properties that surpass the sum of individual component characteristics. Tracing the origin of complex molecular behaviors from interactions is critical to understanding ensemble emergence, and requires insights at the single-molecule level. Electrical signals from single-molecule junctions enable the observation of individual molecular behaviors, as well as intramolecular and intermolecular interactions. This technique provides a foundation for bottom-up explorations of emergent complexity. This Perspective highlights investigations of various interactions via single-molecule junctions, including intramolecular orbital and weak intermolecular interactions and interactions in chemical reactions. It also provides potential directions for future single-molecule junctions in complex system research.

Crystal Growth Regulation of Ruddlesden-Popper Perovskites via Self-Assembly of Semiconductor Spacers for Efficient Solar Cells

The crystal growth and orientation of two-dimensional (2D) perovskite films significantly impact solar cell performance. Here, we incorporated robust quadrupole-quadrupole interactions to govern the crystal growth of 2D Ruddlesden-Popper (RP) perovskites. This was achieved through the development of two unique semiconductor spacers, namely PTMA and 5FPTMA, with different dipole moments. The ((5FPTMA)0.1 (PTMA)0.9 )2 MAn-1 Pbn I3n+1 (nominal n=5, 5F/PTMA-Pb) film shows a preferred vertical orientation, reduced grain boundaries, and released residual strain compared to (PTMA)2 MAn-1 Pbn I3n+1 (nominal n=5, PTMA-Pb), resulting in a decreased exciton binding energy and reduced electron-phonon coupling coefficients. In contrast to PTMA-Pb device with an efficiency of 15.66 %, the 5F/PTMA-Pb device achieved a champion efficiency of 18.56 %, making it among the best efficiency for 2D RP perovskite solar cells employing an MA-based semiconductor spacer. This work offers significant insights into comprehending the crystal growth process of 2D RP perovskite films through the utilization of quadrupole-quadrupole interactions between semiconductor spacers.

Ligand 1 H NMR Chemical Shifts as Accurate Reporters for Protein-Ligand Binding Interfaces in Solution

The availability of high-resolution 3D structural information is crucial for investigating guest-host systems across a wide range of fields. In the context of drug discovery, the information is routinely used to establish and validate structure-activity relationships, grow initial hits from screening campaigns, and to guide molecular docking. For the generation of protein-ligand complex structural information, X-ray crystallography is the experimental method of choice, however, with limited information on protein flexibility. An experimentally verified structural model of the binding interface in the native solution-state would support medicinal chemists in their molecular design decisions. Here we demonstrate that protein-bound ligand 1 H NMR chemical shifts are highly sensitive and accurate probes for the immediate chemical environment of protein-ligand interfaces. By comparing the experimental ligand 1 H chemical shift values with those computed from the X-ray structure using quantum mechanics methodology, we identify significant disagreements for parts of the ligand between the two experimental techniques. We show that quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) ensembles can be used to refine initial X-ray co-crystal structures resulting in a better agreement with experimental 1 H ligand chemical shift values. Overall, our findings highlight the usefulness of ligand 1 H NMR chemical shift information in combination with a QM/MM MD workflow for generating protein-ligand ensembles that accurately reproduce solution structural data.

Structures and Supramolecular Properties of Inclusion Complexes of Anthracene-Triptycene Nanocages with Fullerene Guests and Their Dynamic Motion as Molecular Gyroscopes

Three derivatives of macrocyclic cage compounds consisting of diarylanthracene and triptycene units were synthesized. These nanocages formed host-guest complexes with C₆₀ and other fullerene guests as confirmed by ¹ H NMR and fluorescence spectroscopy. The association constant of the mesityl and 2,4,6-tributoxyphenyl derivatives with C₆₀ was determined to be 2.2 × 10⁴  L mol^-1 , which was larger than that of the pentafluorophenyl derivative. Direct experimental evidence of the complexation was obtained by X-ray diffraction analysis: the guest C₆₀ molecule was included in the cavity via multipoint CH⋅⋅⋅π interactions. Dynamic disorders of the included C₆₀ molecule in variable-temperature X-ray analysis indicated uniaxial motion, such as gyroscopic motion. The unique dynamic behavior of the spherical C₆₀ rotor anchored by the cage stator via CH⋅⋅⋅π interactions in the crystal, as well as substituent effects on the association properties, are discussed with the aid of DFT calculations.

σ-σ Stacked supramolecular junctions

Intermolecular charge transport plays an essential role in organic electronic materials and biological systems. To date, experimental investigations of intermolecular charge transport in molecular materials and electronic devices have been restricted to conjugated systems in which π-π stacking interactions are involved. Herein we demonstrate that the σ-σ stacking interactions between neighbouring non-conjugated molecules offer an efficient pathway for charge transport through supramolecular junctions. The conductance of σ-σ stacked molecular junctions formed between two non-conjugated cyclohexanethiol or single-anchored adamantane molecules is comparable to that of π-π stacked molecular junctions formed between π-conjugated benzene rings. The current-voltage characteristics and flicker noise analysis demonstrate the existence of stacked molecular junctions formed between the electrode pairs and exhibit the characteristics of through-space charge transport. Density functional theory calculations combined with the non-equilibrium Green's function method reveal that efficient charge transport occurs between two molecules configured with σ-σ stacking interactions.

London Dispersion Favors Sterically Hindered Diarylthiourea Conformers in Solution

We present an experimental and computational study on the conformers of N,N'-diphenylthiourea substituted with different dispersion energy donor (DED) groups. While the unfolded anti-anti conformer is the most relevant for thiourea catalysis, intramolecular noncovalent interactions counterintuitively favor the folded syn-syn conformer, as evident from a combination of low-temperature nuclear magnetic resonance measurements and computations. In order to quantify the noncovalent interactions, we utilized local energy decomposition analysis and symmetry-adapted perturbation theory at the DLPNO-CCSD(T)/def2-TZVPP and sSAPT0/6-311G(d,p) levels of theory. Additionally, we applied a double-mutant cycle to experimentally study the effects of bulky substituents on the equilibria. We determined London dispersion as the key interaction that shifts the equilibria towards the syn-syn conformers. This preference is likely a factor why such thiourea derivatives can be poor catalysts.

Electrostatic penetration effects stand at the heart of aromatic π interactions

The nature of the interaction in benzene-containing dimers has been analysed by means of Symmetry Adapted Perturbation Theory (SAPT). The total interaction energy and the preference for the dimers to adopt slipped structures are, apparently, consequence of the balance between repulsion and dispersion. However, our results indicate that this only holds when trends are analysed using fixed intermolecular distances. Employing the most favourable separations between rings it turns out that the changes on the total interaction energy are mostly controlled by electrostatics, while repulsion and dispersion cancel each other to a great extent. Most of the electrostatic contribution is accounted for by electrostatic penetration, so a description based on multipoles should not be employed to rationalise the interaction in benzene-containing dimers. The changes on the interaction energy in benzene-containing dimers are steered by electrostatic penetration which, though often overlooked, plays an essential role for the description of aromatic π interactions.

Towards a predictive model for polymer solubility using the noncovalent interaction index: polyethylene as a case study

In this work we present the development of a novel, quantitative solubility descriptor based on the non-covalent interaction index. It is presented as a more insightful alternative to Hansen's solubility parameters and the COSMO model to assess and predict polymer solubility in different solvents. To this end, we studied the solvation behaviour as a function of the chain length of a single chain of arguably the most simple polymer, polyethylene, in anisole (solvent) and methanol (poor solvent) via molecular dynamics simulations. It was found that in anisole the solute maximized its interface with the solvent, whereas in methanol the macromolecule formed rod-like structures by folding on itself once the chain length surpassed a certain barrier. We assessed this behaviour - which can be related to solubility - quantitatively and qualitatively via well-known descriptors, namely the solvation free energy, and the solvent accessible surface area. In addition, we propose the non-covalent interaction (NCI) index as a versatile descriptor, providing information on the strength, as well as the nature, of the solute-solvent interactions, the solute's intramolecular interactions and on the solute's conformation, both qualitatively and quantitatively. Finally, as a quantitative measure for solubility, defined in this context as the solute's tendency to maximize its interactions with the solvent, we propose two new NCI-based descriptors: the relative integrated NCI density and the integrated NCI difference. The former represents the quantitative difference in solute-solvent interactions between a fully extended coil and the actual conformation during simulation and the latter the quantitative difference between the intermolecular (solute-solvent) and the intramolecular (in the solute) non-covalent interactions. The easy interpretation and calculation of these novel quantities open up the possibility of fast, reliable and insightful high-throughput screening of different (anti)solvent and solute combinations.

A Combined Experimental/Quantum-Chemical Study of Tetrel, Pnictogen, and Chalcogen Bonds of Linear Triatomic Molecules

Linear triatomic molecules (CO₂, N₂O, and OCS) are scrutinized for their propensity to form perpendicular tetrel (CO₂ and OCS) or pnictogen (N₂O) bonds with Lewis bases (dimethyl ether and trimethyl amine) as compared with their tendency to form end-on chalcogen bonds. Comparison of the IR spectra of the complexes with the corresponding monomers in cryogenic solutions in liquid argon enables to determine the stoichiometry and the nature of the complexes. In the present cases, perpendicular tetrel and pnictogen 1:1 complexes are identified mainly on the basis of the lifting of the degenerate ν 2 bending mode with the appearance of both a blue and a red shift. Van ′t Hoff plots of equilibrium constants as a function of temperature lead to complexation enthalpies that, when converted to complexation energies, form the first series of experimental complexation energies on sp¹ tetrel bonds in the literature, directly comparable to quantum-chemically obtained values. Their order of magnitude corresponds with what can be expected on the basis of experimental work on halogen and chalcogen bonds and previous computational work on tetrel bonds. Both the order of magnitude and sequence are in fair agreement with both CCSD(T) and DFA calculations, certainly when taking into account the small differences in complexation energies of the different complexes (often not more than a few kJ mol⁻¹) and the experimental error. It should, however, be noted that the OCS chalcogen complexes are not identified experimentally, most probably owing to entropic effects. For a given Lewis base, the stability sequence of the complexes is first successfully interpreted via a classical electrostatic quadrupole–dipole moment model, highlighting the importance of the magnitude and sign of the quadrupole moment of the Lewis acid. This approach is validated by a subsequent analysis of the molecular electrostatic potential, scrutinizing the σ and π holes, as well as the evolution in preference for chalcogen versus tetrel bonds when passing to “higher” chalcogens in agreement with the evolution of the quadrupole moment. The energy decomposition analysis gives further support to the importance/dominance of electrostatic effects, as it turns out to be the largest attractive term in all cases considered, followed by the orbital interaction and the dispersion term. The natural orbitals for chemical valence highlight the sequence of charge transfer in the orbital interaction term, which is dominated by an electron-donating effect of the N or O lone-pair(s) of the base to the central atom of the triatomics, with its value being lower than in the case of comparable halogen bonding situations. The effect is appreciably larger for TMA, in line with its much higher basicity than DME, explaining the comparable complexation energies for DME and TMA despite the much larger dipole moment for DME.

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.

Novel Coordination Mode in the Potassium Mefenamate Trihydrate Polymeric Structure

As a result of the synthesis of mefenamic acid with potassium hydroxide, a salt with a polymeric structure is formed. The one-dimensional polymeric structure was studied by single crystal X-ray diffraction. The potassium cation is coordinated to one oxygen atom of the carboxylate group and six water oxygen atoms. Potassium ions are bridged by oxygen atoms of water molecules. The crystal structure was used as an input to QTAIM and NCI approaches to investigate the K-O interactions linking the cation with the water oxygen and carboxylate groups. The weak K-O interactions of the potassium cation and water oxygen atoms were strong enough to form a polymeric structure. The flexibility of the weak interactions is responsible for a novel coordination mode in the potassium mefenamate trihydrate.

The Crystal Structure and Intermolecular Interactions in Fenamic Acids-Acridine Complexes

In order to improve pharmaceutical properties of drugs, complexes are synthesized as combinations with other chemical substances. The complexes of fenamic acid and its derivatives, such as mefenamic-, tolfenamic- and flufenamic acid, with acridine were obtained and the X-ray structures were discussed. Formation of the crystals is determined by the presence of the intermolecular O–H^…N hydrogen bond that occur between fenamic acids and acridine. Intermolecular interactions stabilizing the crystals such as π^…π stacking, C–H^…X (X = O, Cl) intermolecular hydrogen bonds as well as C–H^…π and other dispersive interactions were analyzed by theoretical methods: the quantum theory of atoms in molecules (QTAIM) and noncovalent interaction (NCI) approaches.

The nature of the intermolecular interaction in (H2X)2 (X = O, S, Se)

Hydrogen bonds (HBs) are crucial non-covalent interactions in chemistry. Recently, the occurrence of an HB in (H2S)2 has been reported (Arunan et al., Angew. Chem., Int. Ed., 2018, 57, 15199), challenging the textbook view of H2S dimers as mere van der Waals clusters. We herein try to shed light on the nature of the intermolecular interactions in the H2O, H2S, and H2Se dimers via correlated electronic structure calculations, Symmetry Adapted Perturbation Theory (SAPT) and Quantum Chemical Topology (QCT). Although (H2S)2 and (H2Se)2 meet some of the criteria for the occurrence of an HB, potential energy curves as well as SAPT and QCT analyses indicate that the nature of the interaction in (H2O)2 is substantially different (e.g. more anisotropic) from that in (H2S)2 and (H2Se)2. QCT reveals that the HB in (H2O)2 includes substantial covalent, dispersion and electrostatic contributions, while the last-mentioned component plays only a minor role in (H2S)2 and (H2Se)2. The major contributions to the interactions of the dimers of H2S and H2Se are covalency and dispersion as revealed by the exchange-correlation components of QCT energy partitions. The picture yielded by SAPT is somewhat different but compatible with that offered by QCT. Overall, our results indicate that neither (H2S)2 nor (H2Se)2 are hydrogen-bonded systems, showing how the nature of intermolecular contacts involving hydrogen atoms evolves in a group down the periodic table.

Interplay between non-covalent interactions in 1D supramolecular polymers based on 1,4-bis(iodoethynyl)benzene

A study of the solution-phase, solid-state structures of halogen-bonded co-crystals of 1,4-bis(iodoethynyl)benzene (p-BIB) with three salts, namely, decyltrimethylammonium bromide (DTMABr), tetrapropylammonium bromide (TPABr), and tetrabutylammonium bromide (TBABr), has been carried out, along with theoretical calculations. Isothermal titration calorimetry (ITC) showed that the binding constant of bromide with p-BIB in THF is not strongly dependent on the cation, and that the entropic term clearly dominates the enthalpic one in the free energy of binding. In the three crystal structures, the bromide anion acts as a doubly connected node for halogen bonding interactions, which results in linear or angular open chains. The intrachain angles (IBr^-I) of the 1D supramolecular polymers based on p-BIB depend on the geometry and size of the cation and vary from 180° for DMTABr to 75° for TBABr. Non-covalent interaction (NCI) analysis of selected motifs and optimized crystals demonstrates that the balance between halogen bonds, hydrogen bonds, and van der Waals interactions, especially type-I halogenhalogen contacts, determines the crystal structures.

New insight into the role of glutathione reductase in glutathione peroxidase-like activity determination by coupled reductase assay: Molecular Docking Study

Previously we have shown that among 15 substituted salicyloyl (2-hydroxybenzoyl) 5-seleninic acids (SSAs) 4 compounds with longer side chains or a cyclohexyl group exhibit no glutathione peroxidase (GPx)-like activity in the coupled reductase assay. Experimental inhibition of glutathione reductase (GR) by the selenenylsulfide (a main intermediate in the catalytic cycle for GPx-like activity determination) of one of the inactive compounds led us to assess the interactions between 15 selenenylsulfide compounds and the active site of GR by molecular docking. Docking results showed that S and Se atoms in selenenylsulfides of the compounds with no GPx-like activity were beyond 5 Å from S atom of Cys-58 or N atom of imidazole ring of His-467 (Root Mean Square Distances for general assessment of 3 major distances were over 4.8 Å) in the active site, so that they could not be catalyzed to be reduced by GR. Furthermore, their docking scores over 89 Kcal/mol meant that the selenenylsulfides were bound too strongly to the active site to leave it, leading eventually to inhibition of GR. We also applied the molecular docking to other GPx mimics such as ebselen, cyclic seleninate esters and di(propylaminomethylphenyl) diselenides to explain the differences in their GPx-like activity depending to the assays used. Our results suggest that the reduction of a selenenylsulfide by GR plays a positive role in GPx-like activity of GPx mimics in the coupled assay and recommended the prediction of possibility and strength of GPx-like activity by molecular docking before entering experimental research.

Computational screen-out strategy for electrically pumped organic laser materials

Electrically pumped organic lasing is one of the most challenging issues in organic optoelectronics. We present a systematic theoretical investigation to screen out electrical pumping lasing molecules over a wide range of organic materials. With the electronic structure information obtained from time-dependent density functional theory, we calculate multiple photophysical parameters of a set of optical pumping organic laser molecules in our self-developed molecular material property prediction package (MOMAP) to judge whether the electrically pumped lasing conditions can be satisfied, namely, to avoid reabsorption from excitons and/or polarons, and the accumulation of triplet excitons. In addition, a large oscillator strength of S₁ and weak intermolecular π–π interaction are preferred. With these criteria, we are able to conclude that BP3T, BSBCz, and CzPVSBF compounds are promising candidates for electrically pumped lasing, and the proposed computational strategy could serve as a general protocol for molecular design of organic lasing materials.

Though the goal of current organic solid-state laser research remains the realization of electrically pumped lasing, identifying organic semiconductors with ideal properties remains a challenge. Here, the authors report a computational strategy for screening electrical pumping lasing molecules.

What Is the Main Feature Distinguishing the Through-Space Interactions in Cyclophanes from Their Aliphatic Analogues?

Classical cyclophanes with two benzene rings have been compared with cyclophanes with one benzene ring replaced with an aliphatic part and aliphatic compounds, which are cyclophane analogues. Analysis of geometry, atomic charges, and aromatic and steric energy and investigation of intramolecular noncovalent interactions and charge mobility show that there is no special feature that distinguishes the classical cyclophanes from aliphatic analogues, so the definition of cyclophanes can be extended to other compounds.

On the strength of hydrogen bonding within water clusters on the coordination limit

Hydrogen bonds (HB) are arguably the most important noncovalent interactions in chemistry. We study herein how differences in connectivity alter the strength of HBs within water clusters of different sizes. We used for this purpose the interacting quantum atoms energy partition, which allows for the quantification of HB formation energies within a molecular cluster. We could expand our previously reported hierarchy of HB strength in these systems (Phys. Chem. Chem. Phys., 2016, 18, 19557) to include tetracoordinated monomers. Surprisingly, the HBs between tetracoordinated water molecules are not the strongest HBs despite the widespread occurrence of these motifs (e.g., in ice I_h ). The strongest HBs within H₂ O clusters involve tricoordinated monomers. Nonetheless, HB tetracoordination is preferred in large water clusters because (a) it reduces HB anticooperativity associated with double HB donors and acceptors and (b) it results in a larger number of favorable interactions in the system. Finally, we also discuss (a) the importance of exchange-correlation to discriminate among the different examined types of HBs within H₂ O clusters, (b) the use of the above-mentioned scale to quickly assess the relative stability of different isomers of a given water cluster, and (c) how the findings of this research can be exploited to indagate about the formation of polymorphs in crystallography. Overall, we expect that this investigation will provide valuable insights into the subtle interplay of tri- and tetracoordination in HB donors and acceptors as well as the ensuing interaction energies within H₂ O clusters.

Stacking interactions of resonance-assisted hydrogen-bridged rings and C6-aromatic rings

Stacking interactions between six-membered resonance-assisted hydrogen-bridged (RAHB) rings and C6-aromatic rings were systematically studied by analyzing crystal structures in the Cambridge Structural Database (CSD). The interaction energies were calculated by quantum-chemical methods. Although the interactions are stronger than benzene/benzene stacking interactions (-2.7 kcal mol-1), the strongest calculated RAHB/benzene stacking interaction (-3.7 kcal mol-1) is significantly weaker than the strongest calculated RAHB/RAHB stacking interaction (-4.7 kcal mol-1), but for a particular composition of RAHB rings, RAHB/benzene stacking interactions can be weaker or stronger than the corresponding RAHB/RAHB stacking interactions. They are also weaker than the strongest calculated stacking interaction between five-membered saturated hydrogen-bridged rings and benzene (-4.4 kcal mol-1) and between two five-membered saturated hydrogen-bridged rings (-4.9 kcal mol-1). SAPT energy decomposition analyses show that the strongest attractive term in RAHB/benzene stacking interactions is dispersion, however, it is mostly canceled by a repulsive exchange term; hence the geometries of the most stable structures are determined by an electrostatic term.

NCIPLOT4: Fast, Robust, and Quantitative Analysis of Noncovalent Interactions

The NonCovalent Interaction index (NCI) enables identification of attractive and repulsive noncovalent interactions from promolecular densities in a fast manner. However, the approach remained up to now qualitative, only providing visual information. We present a new version of NCIPLOT, NCIPLOT4, which allows quantifying the properties of the NCI regions (volume, charge) in small and big systems in a fast manner. Examples are provided of how this new twist enables characterization and retrieval of local information in supramolecular chemistry and biosystems at the static and dynamic levels.

The relative position of π-π interacting rings notably changes the nature of the substituent effect

The substituent effect in monosubstituted benzene dimers mostly follows changes on electrostatics mainly controlled by the direct interaction of the substituent and the other phenyl ring, whereas the contribution from the interacting rings is smaller. As the substituent is located further away the two contributions become of similar magnitude, so the global result is a combination of both effects. These trends are confirmed in larger systems containing a contact between phenyl rings; at closer distances the interaction of the substituent and the other ring clearly dominates over changes associated with the substituted ring, but as the substituent is located further away its contribution decreases and the contribution from the ring becomes more relevant. Care should be taken in larger systems because the observed energy change can also be affected by interactions with other regions of the molecule not directly involved in the π-π interaction.

Unveiling the Impact of Aggregation on Optical Anisotropy of Triazaacephenanthrylene Single Crystals. A Combined Quantum Crystallography and Conceptual Density Functional Theory Approach

Triazaacephenanthrylene (TAAP) triclinic single crystals show substantial optical anisotropy of absorption and fluorescence. The maximum effect can be correlated with the direction perpendicular to the plane of chromophores connected in a head-to-tail manner via weak dispersive interactions. This phenomenon is uncommon as usually the existence of postulated π···π interactions between the molecules forming dimers or stacks cause quenching of fluorescence. Herein we present a comprehensive study of inter- and intramolecular interactions in the crystal of TAAP enriched with the investigation of aromaticity. Our results show that intramolecular interactions stabilize the overall conformation of the molecule whereas dispersive forces determine the aggregation between TAAP molecules. In fact, there is no conventional π···π interaction between the molecules in the dimer. Instead, we observed a close contact between the lone pair of the bridgehead N10B atom and π-deficient pyrazine ring from an adjacent molecule. Optical anisotropy in TAAP crystals was directly correlated with the alignment of the molecular transition dipole moments caused by specific molecular self-assembly.

Substituent effect on inter-ring interaction in paracyclophanes

The theoretical calculations, namely multipole-derived charge analysis, quantum theory of atom in molecules, and non-bonding interaction (NCI), were performed for [2.2]paracyclophanes, [2.2]paracyclophane-7,9-dienes, and [3.3]paracyclophanes optimized at B3LYP/6-311++G** level, including dispersion correction. The substituent effect of the electron donor N(Me)2 and electron acceptor NO2 group and the influence of the length of bridges joining the aromatic ring on aromatic ring interaction energy (AIE) and strain energy were discussed. The local and electrostatic character of the substituent effect in paracyclophanes was shown. The presence of the weak orbital through-space C···C interaction between the [3.3]paracyclophane ring and weak CH···O hydrogen bonds between the substituents in the different rings was shown.

Prototypical cyclohexane dimers: spectroscopic evidence for σ stacking at low temperatures

In this work, the first unambiguous spectroscopic evidence for the existence of σ stacking interactions in cyclohexane dimers has been provided using matrix isolation infrared spectroscopy.

Curvature and size effects hinder halogen bonds with extended π systems

The curvature of aromatic systems strengthens the interaction by the concave face while it weakens by the convex one. Parallel structures are favoured over halogen bonded ones.

Classifying the chemical bonds involving the noble-gas atoms

The Ng–X bonds are classified into covalent (Cov), and different types of non-covalent (nCov), or partially-covalent (pCov) interactions.

The Na-O bond in sodium fenamate

The one-dimensional polymeric structure of sodium diaquafenamate-water (1/1) was studied by X-ray diffraction. The sodium cation is coordinated to one oxygen atom of the carboxylate group and to four water oxygen atoms. To characterize the Na-O bonds, the quantum theory of atoms in molecules (QTAIM) and noncovalent interaction (NCI) approaches have been used. Both methods confirmed that the Na-O bonds are very weak, comparable with the weak N-H...O intramolecular hydrogen bond. The polymeric structure is stabilized by the interaction of the sodium cation with the surrounding water molecules.

Uncovering the Molecular Interactions in the Catalytic Loop That Modulate the Conformational Dynamics in Protein Tyrosine Phosphatase 1B

Active-site loops are integral to the function of numerous enzymes. They enable substrate and product binding and release, sequester reaction intermediates, and recruit catalytic groups. Here, we examine the catalytic loop in the enzyme protein tyrosine phosphatase 1B (PTP1B). PTP1B has a mobile so-called WPD loop (named for its three N-terminal residues) that initiates the dephosphorylation of phosphor-tyrosine substrates upon loop closure. We have combined X-ray crystallography, solution NMR, and pre-steady-state kinetics experiments on wild-type and five WPD loop mutants to identify the relationships between the loop structure, dynamics, and function. The motions of the WPD loop are modulated by the formation of weak molecular interactions, where perturbations of these interactions modulate the conformational equilibrium landscape. The point mutants in the WPD loop alter the loop equilibrium position from a predominantly open state (P185A) to 50:50 (F182A), 35:65 (P188A), and predominantly closed states (T177A and P188A). Surprisingly, there is no correlation between the observed catalytic rates in the loop mutants and changes to the WPD loop equilibrium position. Rather, we observe a strong correlation between the rate of dephosphorylation of the phosphocysteine enzyme intermediate and uniform millisecond motions, not only within the loop but also in the adjacent α-helical domain of PTP1B. Thus, the control of loop motion and thereby catalytic activity is dispersed and resides within not only the loop sequence but also the surrounding protein architecture. This has broad implications for the general mechanistic understanding of enzyme reactions and the role that flexible loops play in the catalytic cycle.

A Benchmark of Density Functional Approximations For Thermochemistry and Kinetics of Hydride Reductions of Cyclohexanones

The performance of density functionals and wavefunction methods for describing the thermodynamics and kinetics of hydride reductions of 2-substituted cyclohexanones has been evaluated for the first time. A variety of exchange correlation functionals ranging from generalized gradient approximations to double hybrids have been tested and their performance to describe the facial selectivity of hydride reductions of cyclohexanones has been carefully assessed relative to the CCSD(T) method. Among the tested methods, an approach in which single-point energy calculations using the double hybrid B2PLYP-D3 functional on ωB97X-D optimized geometries provides the most accurate transition state energies for these kinetically-controlled reactions. Moreover, the role of torsional strain, temperature, solvation, noncovalent interactions on the stereoselectivity of these reductions was elucidated. Our results indicate a prominent role of the substituent on the cis/trans ratios driven by the delicate interplay between torsional strain and dispersion interactions.

Weak interactions in furan dimers

Dimers of furan, 2,3-dihydrofuran, 2,5-dihydrofuran and tetrahydrofuran were investigated with the use of theoretical methods to determine the interactions that keep the molecules together. The QTAIM and NCI methods confirmed that for furan dimers the C-H⋯O hydrogen bond and stacking interactions can form the dimers with similar energy. For 2,3-dihydrofuran, 2,5-dihydrofuran and tetrahydrofuran, the decisive mechanism of dimer formation is the stacking interaction between the furan rings.

On the Nature of σ-σ, σ-π, and π-π Stacking in Extended Systems

Stacking interactions have been evaluated, employing computational methods, in dimers formed by analogous aliphatic and aromatic species of increasing size. Changes in stability as the systems become larger are mostly controlled by the balance of increasing repulsion and dispersion contributions, while electrostatics plays a secondary but relevant role. The interaction energy increases as the size of the system grows, but it does much faster in π-π dimers than in σ-π complexes and more remarkably than in σ-σ dimers. The main factor behind the larger stability of aromatic dimers compared to complexes containing aliphatic molecules is related to changes in the properties of the aromatic systems due to electron delocalization leading to larger dispersion contributions. Besides, an extra stabilization in π-π complexes is due to the softening of the repulsive wall in aromatic species that allows the molecules to come closer.

Steric "attraction": not by dispersion alone

Non-covalent interactions between neutral, sterically hindered organic molecules generally involve a strong stabilizing contribution from dispersion forces that in many systems turns the 'steric repulsion' into a 'steric attraction'. In addition to London dispersion, such systems benefit from electrostatic stabilization, which arises from a short-range effect of charge penetration and gets bigger with increasing steric bulk. In the present work, we quantify this contribution for a diverse set of molecular cores, ranging from unsubstituted benzene and cyclohexane to their derivatives carrying tert-butyl, phenyl, cyclohexyl and adamantyl substituents. While the importance of electrostatic interactions in the dimers of sp²-rich (e.g., π-conjugated) cores is well appreciated, less polarizable assemblies of sp³-rich systems with multiple short-range CH···HC contacts between the bulky cyclohexyl and adamantyl moieties are also significantly influenced by electrostatics. Charge penetration is drastically larger in absolute terms for the sp²-rich cores, but still has a non-negligible effect on the sp³-rich dimers, investigated herein, both in terms of their energetics and equilibrium interaction distances. These results emphasize the importance of this electrostatic effect, which has so far been less recognized in aliphatic systems compared to London dispersion, and are therefore likely to have implications for the development of force fields and methods for crystal structure prediction.

Topology and Electronic Density Driven Generation of Alkali Cation Complexes

The formation and characterization of K⁺ and Cs⁺ complexes originating from the cooperativity of three non-covalent interactions is explored. The tridimensional preorganization of the naphthothiophene platform displays a favorable well-defined bay region combining a π fragment and a carbonyl moiety flanking a central sulfur atom. A joint theoretical and experimental infrared multiple photon dissociation (IRMPD) study allowed deciphering the key contribution of the orthogonal phenyl fragment to the elaboration of alkali metal complexes. In combination with S and CO interactions, the π-cation interaction significantly enhances the binding energies of naphthothiophene derivatives.

Bending Carbon Nanoforms for Supramolecular Recognition: A Topological Study on Hemifullerene-Based Aggregates

Buckybowls have risen as appealing fullerene fragment derivatives. Their intrinsic curvature has been exploited in the generation of host-guest supramolecular assemblies, not only through concave-convex complementarity but also through less-known concave-concave staggered arrangements. Whereas the stabilization of bowl-in-bowl dispositions has been ascribed to efficient π-π forces together with favorable dipole-dipole interactions, a detailed analysis on the forces guiding the formation of the staggered arrangements is missing so far. Herein, we present a thorough theoretical characterization of bowl-in-bowl vs staggered hemifullerene-based homodimers and heterodimers with the electron-donor truxTTF molecule, as test cases, under the density functional theory and by means of chemical bonding techniques. Our results clearly reveal strong and localized noncovalent signatures, together with an enhanced orbital interaction, associated with CH-π and sulfur-mediated interactions governing the staggered formation. Bending the fullerene fragment is demonstrated to favor the stabilization in both homo- and heterodimers, in good accord with the depletion in the π-electron density calculated upon increasing the buckybowl curvature. The optimal buckybowl curvature for the highest interaction energy is, however, dependent on the type of supramolecular assembly (bowl-in-bowl vs staggered) and the concave region to which hemifullerene approaches truxTTF. Interestingly, two regimes are found as a function of buckybowl curvature for hemifullerene homodimers: bowl-in-bowl dispositions are calculated more stable at low curvatures whereas staggered dimers prevail for highly curved buckybowls. Our results highlight the potential of discrete CH-π and sulfur-mediated interactions to generate unconventional staggered supramolecular arrangements toward the development of a new and unexplored host-guest chemistry.

Exploiting the interactions of aromatic units for folding and assembly in aqueous environments

A variety of non-covalent interactions (including hydrogen bonding, ionic interactions, metal coordination and desolvation/solvation) have been utilized to organize oligomers into well-defined structures. Herein is described a survey of aromatic foldamers that capitalize on electrostatic complementarity of substituted aromatic units to drive folding and assembly in aqueous environments. A brief description of recent advances in the understanding of aromatic interactions is provided, followed by examples of foldamers that exploit interactions between aromatic units to drive their assembly in predictable fashion. The history of our aromatic foldamers is traced from the first structure designed to fold into a pleated structure in an aqueous environment to a heteroduplex system more related to nucleic acids. Taken together, the results demonstrate that electrostatic complementarity of aromatic units provides a versatile framework for driving predictable folding and assembly in aqueous environments.

Ultra-low friction mechanism of highly sp3-hybridized amorphous carbon controlled by interfacial molecule adsorption

The friction behaviors of highly sp³-hybridized carbon films, including ultra-nanocrystalline diamond and diamond-like carbon materials, strongly depend on atmospheres.