The bifunctional catalytic role of water clusters in the formation of acid rain

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.

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.

On the rearrangements of biologically-relevant vinyl allene oxides to cis-cyclopentenones, ketols, and Favorskii-type carboxylic acids

In addition to stereodefined cis-cyclopentenones, the rearrangement of naturally-occurring vinyl allene oxides can provide ketols, cyclopropylcarbinols, and Favorskii-type bis-(Z)-but-2-en-1-yl acetic acids. These processes have been studied by DFT computations using (Z)-but-1-en-1-yl allene oxides as model systems. Prior studies on the stepwise cascade process starting from (Z)-but-1-en-1-yl allene oxides established as key steps the ring opening of the oxirane to give oxidopentadienyl biradicals, and their isomerization through formation of alkenylcyclopropanone intermediates prior to the conrotatory electrocyclic ring closure to cis-configured cyclopentenones. Under neutral or under acidic conditions, the corresponding ketols and cyclopropylcarbinols have been computationally characterized as resulting from S_N2, S_N1 and S_N1'-type processes, showing that the rearrangement of vinyl allene oxides is pH-dependent. Moreover, stereoconvergent base-induced Favorskii-type rearrangements to provide bis-(Z)-but-1-en-1-yl substituted acetic acids have also been justified. Since the model system captures the structural features of the vinyl allene oxides of biological relevance, our computations provide the most comprehensive overview of the complex reactivity of these natural species.

Bonding analysis of water clusters using quasi-atomic orbitals

The quasi-atomic orbital (QUAO) bonding analysis introduced by Ruedenberg and co-workers is used to develop an understanding of the hydrogen bonds in small water clusters, from the dimer through the hexamer (bag, boat, book, cyclic, prism and cage conformers). Using kinetic bond orders as a metric, it is demonstrated that as the number of waters in simple cyclic clusters increases, the hydrogen bonds strengthen, from the dimer through the cyclic hexamer. However, for the more complex hexamer isomers, the strength of the hydrogen bonds varies, depending on whether the cluster contains double acceptors and/or double donors. The QUAO analysis also reveals the three-center bonding nature of hydrogen bonds in water clusters.

Water clusters as bifunctional catalysts in organic chemistry: the hydrolysis of oxirane and its methyl derivatives

This contribution explores the bifunctional catalytic activity of water clusters ((H₂O)_n with n = 1-5) in organic chemistry similar to that observed in the formation of H₂SO₄ in acid rain (Chem. Commun., 53, 3516, (2017)). We considered for this purpose the Hydrolysis of Epoxides (HE), in particular, that of oxirane and its methyl derivatives. Surrounding water molecules with H-bond cooperative effects decrease the activation energy of the rate-limiting step of HE in condensed phase, especially when they lead to an anti-periplanar attack on the alkoxide leaving group. Furthermore, the water molecules have a bifunctional catalytic role in HE by (i) increasing the nucleophilic and electrophilic character of the attacking oxygen atom and the leaving group of the reaction, respectively, and (ii) placing the reactants in a suitable disposition for the substitution to occur. Overall, this investigation provides relevant insights into the collective action of water molecules on organic reactions in neutral, basic and acid media.

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.

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.

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.

Thermally Induced Hydrogen-Bond Rearrangements in Small Water Clusters and the Persistent Water Tetramer

Small water clusters absorb heat and catalyze pivotal atmospheric reactions. Yet, experiments produced conflicting results on water cluster distribution under atmospheric conditions. Additionally, it is unclear which "phase transitions" such clusters exhibit, at what temperatures, and what are their underlying molecular mechanisms. We find that logarithmically small tails in the radial probability densities of (H₂O) _n clusters (n = 2 - 6) provide direct testimony for such transitions. Using the best available water potential (MB-pol), an advanced thermostating algorithm (g-BAOAB), and sufficiently long trajectories, we map the "bifurcation", "melting", and (hitherto unexplored) "vaporization" transitions, finding that both melting and vaporization proceed via a "monomer on a ring" conformer, exhibiting huge distance fluctuations at the vaporization temperatures (T _v). T _v may play a role in determining the atmospheric cluster size distribution such that the dimer and tetramer, with their exceptionally low/high T _v values, are under/over-represented in these distributions, as indeed observed in nondestructive mass spectrometric measurements.

Understanding the nature of bonding interactions in the carbonic acid dimers

Carbonic acid dimer, (CA)₂, (H₂CO₃)₂, helps to explain the existence of this acid as a stable species, different to a simple sum between carbon dioxide and water. Five distinct, well characterized types of intermolecular interactions contribute to the stabilization of the dimers, namely, C=O⋯H-O, H-O⋯H-O, C=O⋯C=O, C=O⋯O-H, and C-O⋯O-H. In many cases, the stabilizing hydrogen bonds are of at least the same strength as in the water dimer. We dissect the nature of intermolecular interactions and assess their influence on stability. For a set of 40 (H₂CO₃)₂ isomers, C=O⋯H-O hydrogen bonds between the carbonyl oxygen in one CA molecule and the acidic hydrogen in the hydroxyl group at a second CA molecule are the major stabilizing factors because they exhibit the shortest interaction distances, the largest orbital interaction energies, and the largest accumulation of electron densities around the corresponding bond critical points. In most cases, these are closed-shell hydrogen bonds, however, in a few instances, some covalent character is induced. Bifurcated hydrogen bonds are a common occurrence in the dimers of carbonic acid, resulting in a complex picture with multiple orbital interactions of various strengths. Two anti-anti monomers interacting via the strongest C=O⋯H-O hydrogen bonds are the ingredients for the formation of the lowest energy dimers. Graphical Abstract Carbonic acid dimer, (CA)₂, (H₂CO₃)₂, helps explaining the existence of this acid as a stable species, different to a simple sum between carbon dioxide and water. Five distinct, well-characterized types of intermolecular interactions contribute to the stabilization of the dimers, namely, C=O⋯H-O, H-O⋯O-H, C=O⋯C=O, C=O⋯O-C, and C-O⋯O-C. In many cases, the stabilizing hydrogen bonds are of at least the same strength as in the water dimer.

Partition of electronic excitation energies: the IQA/EOM-CCSD method

Different developments in chemistry and emerging technologies have generated a renewed interest in the properties of molecular excited states. We present herein the partition of black-box, size-consistent equation-of-motion coupled cluster singles and doubles (EOM-CCSD) excitation energies within the framework of the interacting quantum atoms (IQA) formalism. We denote this method as IQA/EOM-CCSD. We illustrate this approach by considering small molecules used often in the study of excited states. This investigation shows how the combination of IQA and EOM-CCSD may provide valuable insights into the molecular changes induced by electron excitation via the real space distribution of the energy of an absorbed photon in a molecular system. Our results reveal (i) the most energetically deformed atomic basins and (ii) the most affected covalent and non-covalent interactions within a molecule due to a given photoexcitation. In other words, this kind of analysis provides insights into the spatial energetic redistribution accompanying an electronic excitation, with interesting foreseeable applications in the rational design of photoexcitations with tailored chemical effects. Altogether, we expect that the IQA/EOM-CCSD excitation energy partition will prove useful in the understanding of systems and processes of interest in photophysics and photochemistry.

An interacting quantum atom study of model SN 2 reactions (X- ···CH3 X, X = F, Cl, Br, and I)

The quantum chemical topology method has been used to analyze the energetic profiles in the X-  + CH3 X → XCH3  + X- SN 2 reactions, with X = F, Cl, Br, and I. The evolution of the electron density properties at the BCPs along the reaction coordinate has been analysed. The interacting quantum atoms (IQA) method has been used to evaluate the intra-atomic and interatomic energy variations along the reaction path. The different energetic terms have been examined by the relative energy gradient method and the ANANKE program, which enables automatic and unbiased IQA analysis. Four of the six most important IQA energy contributions were needed to reproduce the reaction barrier common to all reactions. The four reactions considered share many common characteristics but when X = F a number of particularities occur. © 2017 Wiley Periodicals, Inc.

Neutral and anionic phosphate-diesters as molecular templates for the encapsulation of a water dimer

Neutral and anionic phosphate diesters act as molecular templates for the encapsulation of water dimers.

Hydrogen-Bond Weakening through π Systems: Resonance-Impaired Hydrogen Bonds (RIHB)

Redefining interactions: The concept of the resonance-impaired hydrogen bond (RIHB) as an interaction in which a conjugated π system strongly impairs the formation of a hydrogen bond (HB) is introduced. A typical HB involving charged species can have a formation energy of tens of kcal mol^-1 , whereas the corresponding value for the examined RIHB is only 2.6 kcal mol^-1 . Quantum chemical topology tools are used to analyse the low formation energy of the studied RIHBs.