Encapsulation of charged halogens by the 512 water cage.

This study focuses on the encapsulation of the entire series of halides by the 512 cage of twenty water molecules and on the characterization of water to water and water to anion interactions. State-of-the-art computations are used to determine equilibrium geometries, energy related quantities, and thermal stability towards dissociation and to dissect the nature and strength of intermolecular interactions holding the clusters as stable units. Two types of structures are revealed: heavily deformed cages for F- indicating a preference for microsolvation, and slightly deformed cages for the remaining anions indicating a preference for encapsulation. The primary variable dictating the properties of the clusters is the charge density of the central halide, with the most severe effects observed for the F- case. For the remaining halides, the anion may be safely viewed as a sort of "big electron" with little local disruptive power, enough to affect the network of non-covalent hydrogen bonds in the cage, but not enough to break it. Gibbs energies for dissociation either into cavity and halide or into water molecules and halide suggest that, in a similar way as to methane clathrate, a more weakly bonded complex that has been detected in the gas phase, all halide containing clathrate-like structures should be amenable to experimental detection in the gas phase at moderate temperature and pressure conditions.

Interactions and reactivity in crystalline intermediates of mechanochemical cyclorhodation reactions.

There is experimental evidence that solid mixtures of the rhodium dimer [Cp*RhCl2]2 and benzo[h] quinoline (BHQ) produce two different polymorphic molecular cocrystals called 4α and 4ß under ball milling conditions. The addition of NaOAc to the mixture leads to the formation of the rhodacycle [Cp*Rh-(BHQ)Cl], where the central Rh atom retains its tetracoordinate character. Isolate 4ß reacts with NaOAc leading to the same rhodacycle while isolate 4α does not under the same conditions. We show that the puzzling difference in reactivity between the two cocrystals can be traced back to fundamental aspects of the intermolecular interactions between the BHQ and [Cp*RhCl2]2 fragments in the crystalline environment. To support this view, we report a number of descriptors of the nature and strength of chemical bonds and intermolecular interactions in the extended solids and in a cluster model. We calculate formal quantum mechanical descriptors based on electronic structure, electron density, and binding and interaction energies including an energy decomposition analysis. Without exception, all descriptors point to 4ß being a transient structure higher in energy than 4α with larger local and global electrophilic and nucleophilic powers, a more favorable spatial and energetic distribution of the frontier orbitals, and a more fragile crystal structure.

The Role of Crystalline Intermediates in Mechanochemical Cyclorhodation Reactions Elucidated by in-Situ X-ray Powder Diffraction and Computation.

The occurrence of crystalline intermediates in mechanochemical reactions might be more widespread than previously assumed. For example, a recent study involving the acetate-assisted C-H activation of N-Heterocycles with [Cp*RhCl2 ]2 by ball milling revealed the formation of transient cocrystals between the reagents prior to the C-H activation step. However, such crystalline intermediates were only observed through stepwise intervallic ex-situ analysis, and their exact role in the C-H activation process remained unclear. In this study, we monitored the formation of discrete, stoichiometric cocrystals between benzo[h]quinoline and [Cp*RhCl2 ]2 by ball milling using in-situ synchrotron X-ray powder diffraction. This continuous analysis revealed an initial cocrystal that transformed into a second crystalline form. Computational studies showed that differences in noncovalent interactions made the [Cp*RhCl2 ]2 unit in the later-appearing cocrystal more reactive towards NaOAc. This demonstrated the advantage of cocrystal formation before the acetate-assisted metalation-deprotonation step, and how the net cooperative action of weak interactions between the reagents in mechanochemical experiments can lead to stable supramolecular assemblies, which can enhance substrate activation under ball-milling conditions. This could explain the superiority of some mechanochemical reactions, such as acetate-assisted C-H activation, compared to their solution-based counterparts.

Aromaticity: Quo Vadis.

Aromaticity is one of the most deeply rooted concepts in chemistry. But why, if two-thirds of existing compounds can be classified as aromatic, is there no consensus on what aromaticity is? σ-, π-, δ-, spherical, Möbius, or all-metal aromaticity… why are so many attributes needed to specify a property? Is aromaticity a dubious concept? This perspective aims to reflect where the aromaticity community is and where it is going.

Encapsulation of Astatide by a water cage.

A 512 cage of (H2O)20 consisting of 30 hydrogen bonds encapsulates Astatide with little geometrical distortion. The cage is marginally destabilized but the non-covalent interactions are actually strengthened. Host⋯cage interactions in the [At@(H2O)20]- cluster are anti-electrostatic, placing both negatively charged atoms in direct contact as in Atδ-⋯δ-O-Hδ+. An orbital interaction analysis reveals that explicit host⋯cage contacts are "inverted" hydrogen bonds. That is, the same type of donor→acceptor charge transfer as in hydrogen bonding, with no proton bridging the two negative charges.

The transition state region in nonsynchronous concerted reactions.

The critical and vanishing points of the reaction force F(ξ) = -dV(ξ)/dξ yield five important coordinates (ξR, ξR* , ξTS, ξP* , ξP) along the intrinsic reaction coordinate (IRC) for a given concerted reaction or reaction step. These points partition the IRC into three well-defined regions, reactants (ξR→ξR* ), transition state (ξR* →ξP* ), and products (ξP* →ξP), with traditional roles of mostly structural changes associated with the reactants and products regions and mostly electronic activity associated with the transition state (TS) region. Following the evolution of chemical bonding along the IRC using formal descriptors of synchronicity, reaction electron flux, Wiberg bond orders, and their derivatives (or, more precisely, the intensity of the electron activity) unambiguously indicates that for nonsynchronous reactions, electron activity transcends the TS region and takes place well into the reactants and products regions. Under these circumstances, an extension of the TS region toward the reactants and products regions may occur.

Water Maintains the UV-Vis Spectral Features During the Insertion of Anionic Naproxen and Ibuprofen into Model Cell Membranes.

UV-vis spectra of anionic ibuprofen and naproxen in a model lipid bilayer of the cell membrane are investigated using computational techniques in combination with a comparative analysis of drug spectra in purely aqueous environments. The simulations aim at elucidating the intricacies behind the negligible changes in the maximum absorption wavelength in the experimental spectra. A set of configurations of the systems constituted by lipid, water, and drugs or just water and drugs are obtained from classical Molecular Dynamics simulations. UV-vis spectra are computed in the framework of atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) approaches together with Time-Dependent Density Functional Theory (TD-DFT). Our results suggest that the molecular orbitals involved in the electronic transitions are the same, regardless of the chemical environment. A thorough analysis of the contacts between the drug and water molecules reveals that no significant changes in UV-vis spectra are a consequence of ibuprofen and naproxen molecules being permanently microsolvated by water molecules, despite the presence of lipid molecules. Water molecules microsolvate the charged carboxylate group as expected but also microsolvate the aromatic regions of the drugs.

Acid Dissociation in (HX)n (H2 O)n Clusters (X=F, Cl, Br, I; n=2, 3).

In this work, we analyze the interactions between two or three hydrogen halide molecules and the same number of water moieties through a systematic exploration of their potential energy surfaces. Our results indicate that the most stable HF and HCl aggregates do not experience dissociation of any of the acid fragments, even with three water molecules. In contrast, in the HBr and HI clusters, one of the acid fragments does dissociate. While the global minimum of (HBr)3 (H2 O)3 is a hydrogen-bridged bihalide anion (BrHBr- ), which is persistent at temperatures up to 203 K, the lowest energy structure of (HI)3 (H2 O)3 has a separated ion pair, but the motif with a bihalide anion (IHI- ) is only 0.2 kcal mol-1 above the global minimum. Among the more stable structures is a broad spectrum of contacts, including water⋯water, HX⋯water, and HX⋯HX hydrogen bonds, halogen bonds, ionic and long-range X⋯H contacts.

Dissecting Bonding Interactions in Cysteine Dimers.

Neutral (n) and zwitterionic (z) forms of cysteine monomers are combined in this work to extensively explore the potential energy surfaces for the formation of cysteine dimers in aqueous environments represented by a continuum. A simulated annealing search followed by optimization and characterization of the candidate structures afforded a total of 746 structurally different dimers held together via 80 different types of intermolecular contacts in 2894 individual non-covalent interactions as concluded from Natural Bond Orbitals (NBO), Quantum Theory of Atoms in Molecules (QTAIM) and Non-Covalent Interactions (NCI) analyses. This large pool of interaction possibilities includes the traditional primary hydrogen bonds and salt bridges which actually dictate the structures of the dimers, as well as the less common secondary hydrogen bonds, exotic X⋯Y (X = C, N, O, S) contacts, and H⋯H dihydrogen bonds. These interactions are not homogeneous but have rather complex distributions of strengths, interfragment distances and overall stabilities. Judging by their Gibbs bonding energies, most of the structures located here are suitable for experimental detection at room conditions.

Formation and evolution of C-C, C-O, C[double bond, length as m-dash]O and C-N bonds in chemical reactions of prebiotic interest.

A series of prebiotic chemical reactions yielding the precursor building blocks of amino acids, proteins and carbohydrates, starting solely from HCN and water is studied here. We closely follow the formation and evolution of the pivotal C-C, C-O, C[double bond, length as m-dash]O, and C-N bonds, which dictate the chemistry of the molecules of life. In many cases, formation of these bonds is set in motion by proton transfers in which individual water molecules act as catalysts so that water atoms end up in the products. Our results indicate that the prebiotic formation of carbon dioxide, formaldehyde, formic acid, formaldimine, glycolaldehyde, glycine, glycolonitrile, and oxazole derivatives, among others, are best described as highly nonsynchronous concerted single step processes. Nonetheless, for all reactions involving double proton transfer, the formation and breaking of O-H bonds around a particular O atom occur in a synchronous fashion, apparently independently from other primitive processes. For the most part, the first process to initiate seems to be the double proton transfer in the reactions where they are present, then bond breaking/formation around the reactive carbon in the carbonyl group and finally rupture of the C-N bonds in the appropriate cases, which are the most reluctant to break. Remarkably, within the limitations of our non-dynamical computational model, the wide ranges of temperature and pressure in which these reactions occur, downplay the problematic determination of the exact constraints on the early Earth.

Microsolvation of electrons by a handful of ammonia molecules.

Microsolvation of electrons in ammonia is studied here via anionic NH3 n - clusters with n = 2-6. Intensive samplings of the corresponding configurational spaces using second-order perturbation theory with extended basis sets uncover rich and complex energy landscapes, heavily populated by many local minima in tight energy windows as calculated from highly correlated coupled cluster methods. There is a marked energetical preference for structures that place the excess electron external to the molecular frame, effectively coordinating it with the three protons from a single ammonia molecule. Overall, as the clusters grow in size, the lowest energy dimer serves as the basic motif over which additional ammonia molecules are attached via unusually strong charge-assisted hydrogen bonds. This is a priori quite unexpected because, on electrostatic grounds, the excess electron would be expected to be in contact with as many protons as possible. Accordingly, a full quantum mechanical treatment of the bonding interactions under the tools provided by the quantum theory of atoms in molecules is carried out in order to dissect and understand the nature of intermolecular contacts. Vertical detachment energies reveal bound electrons even for n = 2.

Initial Recognition and Attachment of the Zika Virus to Host Cells: A Molecular Dynamics and Quantum Interaction Approach.

The zika virus (ZIKV), transmitted to humans from the bites of Aedes Aegypti and Aedes Albopictus mosquitoes produces Zika fever and neurodegenerative disorders that despite affecting millions of people, most recently in Africa and the Americas, has been declared a neglected tropical disease by the World Health Organization. In this work, atomistic molecular dynamics simulations followed by rigorous analysis of the intermolecular interactions reveal crucial aspects of the initial virus⋯cell molecular recognition and attachment, events that trigger the infectious cycle. Previous experimental studies have shown that Dermatan Sulfate (DS) and Chondroitin Sulfate A (CSA), two glycosaminoglycans which are actually epimers to each other and that are structural constituents of receptors expressed in cell membranes, are the preferred anchorage sites, with a marked preference for DS. Our calculations rationalize this preference from a molecular perspective as follows: when free of the virus, DS has one sulfate group that does not participate in intramolecular strong hydrogen bonds, thus, it is readily available to interact with the envelope protein of the virus (Zika-E), then, after formation of the complexes, Zika-E⋯DS exhibits ten strong salt brides connecting the two fragments against only six salt bridges and two hydrogen bonds in Zika-E⋯CSA. Our results complement the current view of the interaction between the virus and the receptor glycosoaminoglycans revealing that the negatively charged carboxylate groups in CSA and DS are just as important as the sulfates because of the formation of equally strong salt bridges with the positively charged Arginine and Lysine aminoacids in the envelope protein of the virus.

Microsolvation versus Encapsulation in Mono, Di, and Trivalent Cations.

The effects of the formal charge in the stability and bonding of water cavities when solvating a cation are studied here using [X(H2 O)20 ]q+ clusters starting with the well known 512 isomer of (water)20 , placing a single mono, di, or trivalent Xq+ cation at the interior, and then optimizing and characterizing the resulting clusters. Highly correlated interaction and deformation energies are calculated using the CCSD(T)-DLPNO formalism. Bonding interactions are characterized using the tools provided by the quantum theory of atoms in molecules, natural bond orbitals, and non-covalent surfaces. Our results indicate that water to water hydrogen bonds are sensibly strengthened resulting in strong cooperative effects, which amount to ≈ 2 ${ \approx 2}$  kcal/mol per hydrogen bond in the bare cavity and to larger values for the systems including the cations. Approximate encapsulation, that is, surrounding the cation by a network of hydrogen bonds akin to the well known methane clathrate seems to be preferred by cations with smaller charge densities while microsolvation, that is, cluster structures having explicit X⋯O contacts seem to be preferred by cations with larger charge densities which severely deform the cavity.

Alkane C-H activation and ligand exchange on silica supported d0 metal alkylidenes: relevance to alkane metathesis.

In this work, we study the ligand exchange process between an alkane and a series of silica supported metal alkylidenes, which may occur by different pathways: C-H addition, σ-bond metathesis, and α-H abstraction. The results indicate that the α-H abstraction pathway is the preferred one, regardless of the catalyst and ligands. This is in contrast to the expected preference for the C-H addition route. When looking for the origin of this preference, our calculations revealed that the α-H abstraction pathway is driven by entropy, which favors the initial dissociation of the alkyl ligand from the catalyst.

Analysis of Conformational Preferences in Caffeine.

High level DLPNO−CCSD(T) electronic structure calculations with extended basis sets over B3LYP−D3 optimized geometries indicate that the three methyl groups in caffeine overcome steric hindrance to adopt uncommon conformations, each one placing a C−H bond on the same plane of the aromatic system, leading to the C−H bonds eclipsing one carbonyl group, one heavily delocalized C−N bond constituent of the fused double ring aromatic system, and one C−H bond from the imidazole ring. Deletion of indiscriminate and selective non-Lewis orbitals unequivocally show that hyperconjugation in the form of a bidirectional −CH3 ⇆ aromatic system charge transfer is responsible for these puzzling conformations. The structural preferences in caffeine are exclusively determined by orbital interactions, ruling out electrostatics, induction, bond critical points, and density redistribution because the steric effect, the allylic effect, the Quantum Theory of Atoms in Molecules (QTAIM), and the non-covalent interactions (NCI), all predict wrong energetic orderings. Tiny rotational barriers, not exceeding 1.3 kcal/mol suggest that at room conditions, each methyl group either acts as a free rotor or adopts fluxional behavior, thus preventing accurate determination of their conformations. In this context, our results supersede current experimental ambiguity in the assignation of methyl conformation in caffeine and, more generally, in methylated xanthines and their derivatives.

Ring Vibrations to Sense Anionic Ibuprofen in Aqueous Solution as Revealed by Resonance Raman.

We unravel the potentialities of resonance Raman spectroscopy to detect ibuprofen in diluted aqueous solutions. In particular, we exploit a fully polarizable quantum mechanics/molecular mechanics (QM/MM) methodology based on fluctuating charges coupled to molecular dynamics (MD) in order to take into account the dynamical aspects of the solvation phenomenon. Our findings, which are discussed in light of a natural bond orbital (NBO) analysis, reveal that a selective enhancement of the Raman signal due to the normal mode associated with the C-C stretching in the ring, νC=C, can be achieved by properly tuning the incident wavelength, thus facilitating the recognition of ibuprofen in water samples.

The Role of Spike Protein Mutations in the Infectious Power of SARS-COV-2 Variants: A Molecular Interaction Perspective.

Specific S477N, N501Y, K417N, K417T, E484K mutations in the receptor binding domain (RBD) of the spike protein in the wild type SARS-COV-2 virus have resulted, among others, in the following variants: B.1.160 (20A or EU2, first reported in continental Europe), B1.1.7 (α or 20I501Y.V1, first reported in the United Kingdom), B.1.351 (ß or 20H/501Y.V2, first reported in South Africa), B.1.1.28.1 (γ or P.1 or 20J/501Y.V3, first reported in Brazil), and B.1.1.28.2 (ζ, or P.2 or 20B/S484K, also first reported in Brazil). From the analysis of a set of bonding descriptors firmly rooted in the formalism of quantum mechanics, including Natural Bond Orbitals (NBO), Quantum Theory of Atoms In Molecules (QTAIM) and highly correlated energies within the Domain Based Local Pair Natural Orbital Coupled Cluster Method (DLPNO-CCSD(T)), and from a set of computed electronic spectral patterns with environmental effects, we show that the new variants improve their ability to recognize available sites to either hydrogen bond or to form salt bridges with residues in the ACE2 receptor of the host cells. This results in significantly improved initial virus⋅⋅⋅cell molecular recognition and attachment at the microscopic level, which trigger the infectious cycle.

To Be or Not To Be? that is the Entropic, Enthalpic, and Molecular Interaction Dilemma in the Formation of (Water)20 Clusters and Methane Clathrate.

A detailed analysis under a comprehensive set of theoretical and computational tools of the thermodynamical factors and of the intermolecular interactions behind the stabilization of a well known set of (water)20 cavities and of the methane clathrate is offered in this work. Beyond the available reports of experimental characterization at extreme conditions of most of the systems studied here, all clusters should be amenable to experimental detection at 1 atm and moderate temperatures since 280 K marks the boundary at which, ignoring reaction paths, formation of all clusters is no longer spontaneous from the 20H2 O→(H2 O)20 and CH4 +20H2 O→CH4 @512 processes. As a function of temperature, a complex interplay leading to the free energy of formation occurs between the destabilizing entropic contributions, mostly due to cluster vibrations, and the stabilizing enthalpic contributions, due to intermolecular interactions and the PV term, is best illustrated by the highly symmetric 512 cage consistently showing signs of stronger intermolecular bonding despite having smaller binding energy than the other clusters. A fluxional wall of attractive non-covalent interactions, arising because of the cumulative effect of a large number of tiny individual charge transfers to the interstitial region, plays a pivotal role stabilizing the CH4 @512 clathrate.

Thermodynamics and Intermolecular Interactions during the Insertion of Anionic Naproxen into Model Cell Membranes.

The insertion process of Naproxen into model dimyristoylphosphatidylcholine (DMPC) membranes is studied by resorting to state-of-the-art classical and quantum mechanical atomistic computational approaches. Molecular dynamics simulations indicate that anionic Naproxen finds an equilibrium position right at the polar/nonpolar interphase when the process takes place in aqueous environments. With respect to the reference aqueous phase, the insertion process faces a small energy barrier of ≈5 kJ mol-1 and yields a net stabilization of also ≈5 kJ mol-1. Entropy changes along the insertion path, mainly due to a growing number of realizable microstates because of structural reorganization, are the main factors driving the insertion. An attractive fluxional wall of noncovalent interactions is characterized by all-quantum descriptors of chemical bonding (natural bond orbitals, quantum theory of atoms in molecules, noncovalent interaction, density differences, and natural charges). This attractive wall originates in the accumulation of tiny transfers of electron densities to the interstitial region between the fragments from a multitude of individual intermolecular contacts stabilizing the tertiary drug/water/membrane system.

A Comprehensive Picture of the Structures, Energies, and Bonding in the Alanine Dimers.

High level quantum mechanical computations and extensive stochastic searches of the potential energy surfaces of the Alanine dimers uncover rich and complex structural and interaction landscapes. A total of 416 strongly bound (up 13.4 kcal mol-1 binding energies at the DLPNO-CCSD(T)/6-311++G(d,p) level corrected by the basis set superposition error and by the zero point vibrational energies over B3LYP-D3 geometries), close energy equilibrium structures were located, bonded via 32 specific types of intermolecular contacts including Y⋅⋅⋅H-X primary and Y⋅⋅⋅H-C secondary hydrogen bonds, H⋅⋅⋅H dihydrogen contacts, and non conventional anti-electrostatic Y δ-⋯ X δ- interactions. The putative global minimum is triply degenerate, corresponding to the structure of the common dimer of a carboxylic acid. All quantum descriptors of chemical bonding point to a multitude of weak individual interactions within each dimer, whose cumulative effect results in large binding energies and in an attractive fluxional wall of non-covalent interactions in the interstitial region between the monomers.