Actinide endohedral inter-metalloid clusters of the group 15 elements

Inter-metalloid clusters in Zintl chemistry have been extensively studied due to their unique electronic structures and potential applications. In this work, we explored a series of actinide endohedral inter-metalloid clusters of the group 15 elements [An@Bi12]4- (An = Th-U) and [An@Sb12]4- using density functional theory (DFT). [Th@Bi12]4- and [U@Bi12]4- exhibit Cs symmetry, while [Pa@Bi12]4- and [An@Sb12]4- (An = Th-U) have C1 structures. Bonding analyses such as bond order, molecular orbitals (MO) and quantum theory of atoms in molecules (QTAIM) show covalent An-Bi/An-Sb bonding in the clusters. All these clusters are highly stable according to the studied formation reactions and may be accessible experimentally. Compared with [An@Bi12]4-, [An@Sb12]4- possesses stronger bonding interactions, mainly arising from the higher electrostatic interaction energy. For clusters with the same group 15 elements, the bonding interactions increase gradually from Th to U, which is mainly determined by the covalent interactions of An-Bi/An-Sb bonding. This work is expected to provide potential avenues for the construction of robust inter-metalloid clusters of the group 15 elements.

Density functional theory study of crown ether-magnesium complexes: from a solvated ion to an ion trap

Metal ion detection rests on host-guest recognition. We propose a theoretical protocol for designing an optimal trap for a desired metal cation. A host for magnesium ions was sought for among derivatives of crown ethers 12-crown-4, 15-crown-5, and 18-crown-6. Mg-crown complexes and their hydrated counterparts with water molecules bound to the cation were optimized using density functional theory. Based on specific geometric criteria, Interacting quantum atoms analysis and density functional theory-based molecular dynamics of Mg-crown complexes immersed in water, crown ethers for optimal accommodation of Mg2+ in aqueous solution were identified. Selectivity of the chosen crowns towards Na+, K+, and Ca2+ ions is addressed.

Electron density-based protocol to recover the interacting quantum atoms components of intermolecular binding energy

A new approach for obtaining interacting quantum atoms-defined components of binding energy of intermolecular interactions, which bypasses the use of standard six-dimensional integrals and two-particle reduced density matrix (2-RDM) reconstruction, is proposed. To examine this approach, three datasets calculated within the density functional theory framework using the def2-TZVP basis have been explored. The first two, containing 53 weakly bound bimolecular associates and 13 molecular clusters taken from the crystal, were used in protocol refinement, and the third one containing other 20 bimolecular and three cluster systems served as a validation reference. In addition, to verify the performance of the proposed approach on an exact 2-RDM, calculations within the coupled cluster formalism were performed for part of the first set systems using the cc-pVTZ basis set. The process of optimization of the proposed parametric model is considered, and the role of various energy contributions in the formation of non-covalent interactions is discussed with regard to the obtained trends.

The Ultrashort Spike-Ring Interaction in Substituted Iron Maiden Molecules

The in forms of molecular iron maidens are known for their unique ultrashort interaction between the apical hydrogen atom or its small substituent and the surface of the benzene ring. It is generally believed that this forced ultrashort X⋯π contact is associated with high steric hindrance, which is responsible for specific properties of iron maiden molecules. The main aim of this article is to investigate the influence of significant charge enrichment or depletion of the benzene ring on the characteristics of the ultrashort C-X⋯π contact in iron maiden molecules. For this purpose, three strongly electron-donating (-NH2) or strongly electron-withdrawing (-CN) groups were inserted into the benzene ring of in-[34,10][7]metacyclophane and its halogenated (X = F, Cl, Br) derivatives. It is shown that, despite such extremely electron-donating or electron-accepting properties, the considered iron maiden molecules surprisingly reveal quite high resistance to changes in electronic properties.

Actinide-doped boron clusters: from borophenes to borospherenes

Similar to graphene and fullerene, metal-doping has been considered to be an effective approach to the construction of highly stable boron clusters. In this work, a series of actinide metal-doped boron clusters AnB₃₆ (An = Pa, Np, Pu, Am, Cm, Bk, and Cf) have been explored using extensive first-principles calculations. We found that the quasi-planar structure of B₃₆ transforms to an endohedral borospherene An@B₃₆ after actinide metal doping. Actinoborospherenes exhibit C_2h symmetry with Pa, Np, and Pu dopants and for Am, Cm, Bk and Cf dopants with larger atomic radii, the symmetry of An@B₃₆ is reduced to C_i. Bonding property analyses such as bond order, molecular orbital (MO) and quantum theory of atoms in molecules (QTAIM) analysis show that the covalency of the An-B bonds in C_2h An@B₃₆ (An = Pa, Np, and Pu) is higher than that in C_i An@B₃₆ (An = Am, Cm, Bk, and Cf). These endohedral borospherenes are robust according to thermodynamic and dynamic analyses. As expected, the C_i An@B₃₆ clusters are less stable compared to C_2h An@B₃₆, which is consistent with the stronger covalent bonds of the latter. These results indicate that the existence of the actinide-boron bonding is essential for the high stability of the An@B₃₆ clusters, confirming that the fullerene-like boron cages can be stabilized by actinide encapsulation. This work is expected to provide potential routes for the construction of robust borospherenes.

The IQA Energy Partition in a Drug Design Setting: A Hepatitis C Virus RNA-Dependent RNA Polymerase (NS5B) Case Study

The interaction of the thumb site II of the NS5B protein of hepatitis C virus and a pair of drug candidates was studied using a topological energy decomposition method called interacting quantum atoms (IQA). The atomic energies were then processed by the relative energy gradient (REG) method, which extracts chemical insight by computation based on minimal assumptions. REG reveals the most important IQA energy contributions, by atom and energy type (electrostatics, sterics, and exchange–correlation), that are responsible for the behaviour of the whole system, systematically from a short-range ligand–pocket interaction until a distance of approximately 22 Å. The degree of covalency in various key interatomic interactions can be quantified. No exchange–correlation contribution is responsible for the changes in the energy profile of both pocket–ligand systems investigated in the ligand–pocket distances equal to or greater than that of the global minimum. Regarding the hydrogen bonds in the system, a “neighbour effect” was observed thanks to the REG method, which states that a carbon atom would rather not have its covalent neighbour oxygen form a hydrogen bond. The combination of IQA and REG enables the automatic identification of the pharmacophore in the ligands. The coarser Interacting Quantum Fragments (IQF) enables the determination of which amino acids of the pocket contribute most to the binding and the type of energy of said binding. This work is an example of the contribution topological energy decomposition methods can make to fragment-based drug design.

Binding Energy Partition of Promising IRAK-4 Inhibitor (Zimlovisertib) for the Treatment of COVID-19 Pneumonia

The technique of Fragment-Based Drug Design (FBDD) considers the interactions of different moieties of molecules with biological targets for the rational construction of potential drugs. One basic assumption of FBDD is that the different functional groups of a ligand interact with a biological target in an approximately additive, that is, independent manner. We investigated the interactions of different fragments of ligands and Interleukin-1 Receptor-Associated Kinase 4 (IRAK-4) throughout the FBDD design of Zimlovisertib, a promising anti-inflammatory, currently in trials to be used for the treatment of COVID-19 pneumonia. We utilised state-of-the-art methods of wave function analyses mainly the Interacting Quantum Atoms (IQA) energy partition for this purpose. By means of IQA, we assessed the suitability of every change to the ligand in the five stages of FBDD which led to Zimlovisertib on a quantitative basis. We determined the energetics of the interaction of different functional groups in the ligands with the IRAK-4 protein target and thereby demonstrated the adequacy (or lack thereof) of the changes made across the design of this drug. This analysis permits to verify whether a given alteration of a prospective drug leads to the intended tuning of non-covalent interactions with its protein objective. Overall, we expect that the methods exploited in this paper will prove valuable in the understanding and control of chemical modifications across FBDD processes.

The physical nature of the ultrashort spike-ring interaction in iron maiden molecules

The so-called 'iron maiden' molecules belong to one of the most interesting subgroups of cyclophanes due to the presence of the ultrashort interaction between the CX apical bond and the benzene ring. This article presents an in-depth theoretical study of 16 'iron maiden' molecules, in which X = H, F, Cl or Br and the side chains are of various lengths and types: CSC, CSCC, CCC, and CCCC. It is shown that the H → F → Cl → Br substitution leads to a significant expansion of the 'iron maiden' molecule. Shorter chains lead to more pronounced effects, while insertion of sulfur atoms into the side chains lowers them. Structural changes are associated with an increase in energetic destabilization of X. Moreover, unlike for H, in the case of X = halogen, the out → in isomerization is energetically disadvantageous. The 'iron maiden' molecules are characterized by the presence of only three X⋯C_Ar bond paths. Particularly noteworthy are unusually large (even up to 32) values of the X⋯C_Ar bond ellipticity, which results from flat electron density distribution. The X⋯π interaction in each of the investigated 'iron maiden' molecule turned out to be multi-center, stabilizing and almost purely covalent in nature as indicated by the definitely dominant percentage (94.8%-101.6%) of the exchange-correlation energy. The spatial hindrance within the 'iron maiden' molecules appears to be not so much due to the X⋯π repulsion, but due to unfavorable steric interactions between X and the CC side bonds. It is also confirmed that some CH⋯HC interactions in aliphatic chains can be very weakly stabilizing.

Proton leap: shuttling of protons onto benzonitrile

The detailed description of chemical transformations in the interstellar medium allows deciphering the origin of a number of small and medium - sized organic molecules. We present density functional theory analysis of proton transfer from the trihydrogen cation and the ethenium cation to benzonitrile, a recently discovered species in the Taurus Molecular Cloud 1. Detailed energy transformations along the reaction paths were analyzed using the interacting quantum atoms methodology, which elucidated how the proton carrier influences the lightness to deliver the proton to benzonitrile's nitrogen atom. The proton carriers' deformation energy represents the largest destabilizing effect, whereas a proton's promotion energy, the benzonitrile-proton Coulomb attraction, as well as non-classical benzonitrile-proton and carrier-proton interaction are the dominant stabilizing energy components. As two ion-molecule reactions proceed without energy barriers, rate constants were estimated using the classical capture theory and were found to be an order of magnitude larger for the reaction with the trihydrogen cation compared to that with the ethenium cation (∼10^-8 and 10^-9 cm³ s^-1, respectively). These results were obtained both with quantum chemical and ab initio molecular dynamics simulations (the latter performed at 10 K and 100 K), confirming that up to 100 K both systems choose energetically undemanding routes by tracking the corresponding minimum energy paths. A concept of a turning point is introduced, which is an equivalent to the transition state in barrierless reactions.

Highly stable actinide(III) complexes supported by doubly aromatic ligands

Owing to the electron-deficient nature of boron atom, the structures and properties of boron clusters can be enriched by doping various metal atoms, including lanthanide metal atoms. Nevertheless, the viability...

Yet another perspective on hole interactions

Hole interactions are known by different names depending on the key atom of the bond (halogen bond, chalcogen bond, hydrogen bond, etc.), and the geometry of the interaction (σ if in line, π if perpendicular to the Lewis acid plane). However, its origin starts with the creation of a Lewis acid by an underlying covalent bond, which forms an electrostatic depletion and a virtual antibonding orbital, which can create non-covalent interactions with Lewis bases. In this (maybe subjective) perspective, we will claim that hole interactions must be defined via the molecular orbital origin of the molecule. Under this premise we can better explore the richness of such bonding patterns. For that, we will study old, recent and new systems, trying to pinpoint some misinterpretations that are often associated with them. We will use as exemplars the triel bonds, a couple of metal complexes, a discussion on convergent σ-holes, and many cases of anti-electrostatic hole interactions.

Macrolactonization Reactions Driven by a Pentafluorobenzoyl Group*

Macrolactones constitute a privileged class of natural and synthetic products with a broad range of applications in the fine chemicals and pharmaceutical industry. Despite all the progress made towards their synthesis, notably from seco-acids, a macrolactonization promoter system that is effective, selective, flexible, readily available, and, insofar as possible, compatible with manifold functional groups is still lacking. Herein, we describe a strategy that relies on the formation of a mixed anhydride incorporating a pentafluorophenyl group which, due to its high electronic activation enables a convenient access to macrolactones, macrodiolides and esters with a broad versatility. Kinetic studies and DFT computations were performed to rationalize the reactivity of the pentafluorophenyl group in macrolactonization reactions.

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.

Orbital-Free Quantum Crystallographic View on Noncovalent Bonding: Insights into Hydrogen Bonds, π⋅⋅⋅π and Reverse Electron Lone Pair⋅⋅⋅π Interactions

A detailed analysis of a complete set of the local potentials that appear in the Euler equation for electron density is carried out for noncovalent interactions in the crystal of a uracil derivative using experimental X-ray charge density. The interplay between the quantum theory of atoms in molecules and crystals and the local potentials and corresponding inner-crystal electronic forces of electrostatic and kinetic origin is explored. Partitioning of crystal space into atomic basins and atomic-like potential basins led us to the definite description of interatomic interaction and charge transfer. Novel physically grounded bonding descriptors derived within the orbital-free quantum crystallography provided the detailed examination of π-stacking and intricate C=O⋅⋅⋅π interactions and nonclassical hydrogen bonds present in the crystal. The donor-acceptor character of these interactions is revealed by analysis of Pauli and von Weizsäcker potentials together with well-known functions, e. g., deformation electron density and electron localization function. In this way, our analysis throws light on aspects of these closed-shell interactions hitherto hidden from the description.

The significance of the metal cation in guanine-quartet – metalloporphyrin complexes

The distinct positions of the divalent metal ions with respect to the porphyrin ring are responsible for different interaction energies between metalloporphyrins and the guanine quartet.

Theoretical prediction of chiral actinide endohedral borospherenes

Actinide encapsulation can form chiral borospherenes, and the covalent character of An–B bonds dominates the formation of these actinoborospherenes.

Anatomy of Base Pairing in DNA by Interacting Quantum Atoms

The formation of purine and pyrimidine base pairs (BPs), which contributes to shaping of the canonical and noncanonical 3D structures of nucleic acids, is one the most investigated phenomena in chemistry and life sciences. In this contribution, the anatomy of the bond energy (BDE) of the base-pairing interaction in 39 different arrangements found experimentally or predicted for DNA structures containing the four common nucleobases (A, C, G, T) in their neutral or protonated forms is described in light of the theory of interacting quantum atoms within the context of the quantum theory of atoms in molecules. The interplay of individual energy components involved in the three stages of the bond formation process (structural deformation, electron-density promotion, and intermolecular interaction) is studied. We recognized that for the neutral BPs, variations in the kinetic and electrostatic contributions to the BDE are rather negligible, leaving the exchange-correlation energy as the main stabilizing component. It is shown that the contribution of the exchange-correlation term can be recovered by including atoms that are formally assumed to be hydrogen bonded (primary interaction). In contrast, to recover the electrostatic component of interaction, one must consider both the primary and secondary (formally nonbonded atoms) interatomic interactions. The results of our study were employed to design new types of BPs with altered bonding anatomy. We demonstrate that improving the electrostatic characteristics of the BPs does not necessarily result in greater interaction energies if weak secondary hydrogen bonding is destroyed. However, the main tuning factor for systems with conserved interacting faces (primary interactions) is the electrostatic component of the interaction energy resulting from the secondary atom-atom electrostatics.

Interacting Quantum Atoms-A Review

The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.

Why Can Cationic Halogen Bond Donors Activate the Ritter-Type Solvolysis of Benzhydryl Bromide but Cationic Hydrogen Bond Donors Can Not?

It is found by experiment that the cationic halogen bond donors (cationic iodoimidazolium compounds) can activate the Ritter-type solvolysis of benzhydryl bromide, while the cationic hydrogen bond donors (cationic imidazolium compounds) could not. To understand the activation mechanism, various noncovalent interactions between benzhydryl bromide and a series of activators in solution, including halogen bond, hydrogen bond, lone pair···π/π+, and C-H···π/π+, were explored theoretically. Our study revealed that the activation difference can be contributed by the variation of the noncovalent interactions. For halogen bond donors, the successful activation is attributed by halogen bond and lone pair···π. The halogen bonds mainly provide the stabilization energy of the ion-pair complex with the help of lone pair···π. For hydrogen-bond donors, the contribution of the hydrogen bond is unable to compensate the like-charge repulsion arising from the generation of the carbocation, leading to the unsuccessful activation. In general, lone pair···π makes a difference.

Unsaturated binuclear homoleptic nickel carbonyl anions Ni2(CO)n− (n = 4–6) featuring double three-center two-electron Ni–C–Ni bonds

The two nickel atoms in the Ni₂(CO)_n⁻ (n = 4–6) complexes are joined by two bridging carbonyl ligands via the sharing three-center two-electron Ni–C–Ni bond in turn to achieve the (16,16), (16,18), and eventually the favored (18,18) configurations.

Multicenter electron-sharing σ-bonding in the AgFe(CO)4− complex

For the AgFe(CO)₄⁻ anion, the silver atom is covalently bonded to the anionic tetracarbonyl-iron, an isolobal analogue of the methyl radical, via a peculiar decentralized electron-sharing σ bond.

Solid State and Solution Study on the Formation of Inorganic Anion Complexes with a Series of Tetrazine-Based Ligands

Four molecules (L1-L4) constituted by an s-tetrazine ring appended with two identical aliphatic chains of increasing length bearing terminal morpholine groups were studied as anion receptors in water. The basicity properties of these molecules were also investigated. Speciation of the anion complexes formed in solution and determination of their stability constants were performed by means of potentiometric (pH-metric) titrations, while further information was obtained by NMR and isothermal titration calorimetry (ITC) measurements. The crystal structures of two neutral ligands (L3, L4) and of their H2L3(ClO4)2∙2H2O, H2L4(ClO4)2∙2H2O, H2L3(PF6)2, and H2L3(PF6)2∙2H2O anion complexes were determined by single crystal X-ray diffraction. The formation of anion-π interactions is the leitmotiv of these complexes, both in solution and in the solid state, although hydrogen bonding and/or formation of salt-bridges can contribute to their stability. Evidence of the ability of these ligands to form anion-π interactions is given by the observation that even the neutral (not-protonated) molecules bind anions in water to form complexes of significant stability, including elusive OH- anions.

Atypical Lone Pair-π Interaction with Quinone Methides in a Series of Imido-Ferrociphenol Anticancer Drug Candidates

Ferrociphenols, especially those possessing a heterocycle at the terminus of an aliphatic chain, display strong anticancer activity through a novel redox mechanism that generates active metabolites such as quinone methides (QMs). X-ray crystallography and UV/Vis spectroscopy reveal that the specific lone pair (lp)-π interaction between a carbonyl group of the imide and the quinone motif of the QM plays an important role in the exceptional cytotoxic behaviour of their imido-ferrociphenol precursors. This intramolecular lp-π interaction markedly enhanced the stability of the QMs and lowered the pK_a values of the corresponding phenol/phenolate couples. As the first example of such a non-covalent interaction that stabilizes QMs remotely, it not only expands the scope of the lp-π interaction in supramolecular chemistry, but also represents a new mode of stabilization of a QM. This unprecedented application of lp-π interactions in imido-ferrociphenol anticancer drug candidates may also have great potential in drug discovery and organocatalyst design.

Co-doped triel-pnicogen graphene as metal-free catalyst for CO oxidation: Role of multi-center covalency

Third and fifth group atoms, named triel and pnicogen, respectively, were used for graphene doping. AlP and AlN structures were selected as co-doped graphene-based catalysts. The electronic structure and catalytic properties of binary AlN, AlP co-doped graphene were investigated through density functional theory (DFT). Results show that the AlP co-doped graphene strictly enhances the oxygen reactivity compared to AlN one. The CO oxidation on AlP and AlN co-doped graphene sheets is mainly done through Eley-Rideal mechanism as follows: CO + O₂➔CO₂ + O_ads and CO + O_ads➔ CO₂. The CO oxidation reaction paths over the AlN and AlP-co-doped graphene have revealed that they can be regarded as the competent catalyst for CO oxidation in the presence of O₂. Mechanistically, both AlP and AlN co-doped graphene catalysts are appropriately active in the first step while the second step is too hard to do regarding the multi-center covalency character between O₂ and AlP co-doped graphene and hence its catalytic efficiency is significantly lower compared to AlN co-doped graphene sheet. Thus, the substitution of a C-C bond with Al-N is an effective way to design the graphene-based catalysts for CO oxidation.

When hydroquinone meets methoxy radical: Hydrogen abstraction reaction from the viewpoint of interacting quantum atoms

Interacting Quantum Atoms methodology is used for a detailed analysis of hydrogen abstraction reaction from hydroquinone by methoxy radical. Two pathways are analyzed, which differ in the orientation of the reactants at the corresponding transition states. Although the discrepancy between the two barriers amounts to only 2 kJ/mol, which implies that the two pathways are of comparable probability, the extent of intra-atomic and inter-atomic energy changes differs considerably. We thus demonstrated that Interacting Quantum Atoms procedure can be applied to unravel distinct energy transfer routes in seemingly similar mechanisms. Identification of energy components with the greatest contribution to the variation of the overall energy (intra-atomic and inter-atomic terms that involve hydroquinone's oxygen and the carbon atom covalently bound to it, the transferring hydrogen and methoxy radical's oxygen), is performed using the Relative energy gradient method. Additionally, the Interacting Quantum Fragments approach shed light on the nature of dominant interactions among selected fragments: both Coulomb and exchange-correlation contributions are of comparable importance when considering interactions of the transferring hydrogen atom with all other atoms, whereas the exchange-correlation term dominates interaction between methoxy radical's methyl group and hydroquinone's aromatic ring. This study represents one of the first applications of Interacting Quantum Fragments approach on first order saddle points. © 2018 Wiley Periodicals, Inc.

Interplay between hydrogen bonding and n→π* interaction in an analgesic drug salicin

The competition and cooperation between weak intermolecular interactions are important in determining the conformational preferences of molecules. Understanding the relative strengths of these effects in the context of potential drug candidates is therefore essential. We use a combination of gas-phase spectroscopy and quantum-chemical calculations to elucidate the nature of such interactions for the analgesic salicin [2-(hydroxymethyl)phenyl β-d-glucopyranoside], an analog of aspirin found in willow bark. Of several possible conformers, only three are observed experimentally, and these are found to correspond with the three lowest energy conformers obtained from density functional theory calculations and simulated Franck-Condon spectra. Natural bond orbital analyses show that these are characterized by a subtle interplay between weak n→π* interaction and conventional strong hydrogen bond, with additional insights into this interaction provided by analysis of quantum theory of atoms in molecules and symmetry-adapted perturbation theory calculations. In contrast, the higher energy conformers, which are not observed experimentally, are mostly stabilized by the hydrogen bond with negligible contribution of n→π* interaction. The n→π* interaction results in a preference for the benzyl alcohol group of salicin to adopt a gauche conformation, a characteristic also found when salicin is bound to the β-glucosidase enzyme. As such, understanding the interplay between these weak interactions has significance in the rationalization of protein structures.

Electron Density Analysis of Metal Clusters with Semi-Interstitial Main Group Atoms. Chemical Bonding in [Co6X(CO)16]- Species

In this work, we propose a careful and thorough analysis of the chemical bond nature in high nuclearity metal carbonyl clusters having semi-interstitial main group atoms. We investigated the species [Co₆X(CO)₁₆]^- (X = As, P), known for a rather interesting conformational flexibility of the cluster (leading to open or closed cages) and a corresponding polymorphism in the solid state (observed at least for X = As). The factors that trigger the molecular isomerism and the nature of X-Co and Co-Co interactions emerge from theoretical calculations and high resolution X-ray diffraction. Both energy and charge density atomic partitioning (QTAIM, EDA, IQA) are employed for this analysis, with the aim of revealing the stabilizing/destabilizing factors of the interaction between the cage and the semi-interstitial atoms in the various conformations.

Protein environment affects the water–tryptophan binding mode. MD, QM/MM, and NMR studies of engrailed homeodomain mutants

Water molecules can interact with the π-face of tryptophan either forming an O–H⋯π hydrogen bond or by a lone-pair⋯π interaction. Surrounding amino acids can favor the one or the other interaction type.

New electron delocalization tools to describe the aromaticity in porphyrinoids

There are several possible pathways in the macrocycle of large porphyrinoids and, among aromaticity indices, only AV_minis capable of recognizing the most aromatic one.

Lone pair-π interactions in biological systems: occurrence, function, and physical origin

Lone pair-π interactions are now recognized as a supramolecular bond whose existence in biological systems is documented by a growing number of examples. They are commonly attributed to electrostatic forces. This review attempts to highlight some recent discoveries evidencing the important role which lone pair-π interactions, and anion-π interactions in particular, play in stabilizing the structure and affecting the function of biomolecules. Special attention is paid to studies exploring the physical origin of these at first glance counterintuitive interactions between a lone pair of electrons of one residue and the π-cloud of another. Recent theoretical work went beyond the popular electrostatic model and inquired the extent to which orbital interactions have to be taken into account. In at least one biologically relevant case-that of anion-flavin interactions-a substantial charge-transfer component has been shown to operate.

Exploiting non-covalent π interactions for catalyst design

Molecular recognition, binding and catalysis are often mediated by non-covalent interactions involving aromatic functional groups. Although the relative complexity of these so-called π interactions has made them challenging to study, theory and modelling have now reached the stage at which we can explain their physical origins and obtain reliable insight into their effects on molecular binding and chemical transformations. This offers opportunities for the rational manipulation of these complex non-covalent interactions and their direct incorporation into the design of small-molecule catalysts and enzymes.