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.

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.

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(H₂ O)₂₀ ]^q+ clusters starting with the well known 5¹² isomer of (water)₂₀ , placing a single mono, di, or trivalent X^q+ 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.

Microsolvation of electrons by a handful of ammonia molecules

Microsolvation of electrons in ammonia is studied here via anionic NH_{3 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.

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.

Solvation Free Energies in Subsystem Density Functional Theory

For many chemical processes the accurate description of solvent effects are vitally important. Here, we describe a hybrid ansatz for the explicit quantum mechanical description of solute-solvent and solvent-solvent interactions based on subsystem density functional theory and continuum solvation schemes. Since explicit solvent molecules may compromise the scalability of the model and transferability of the predicted solvent effect, we aim to retain both, for different solutes as well as for different solvents. The key for the transferability is the consistent subsystem decomposition of solute and solvent. The key for the scalability is the performance of subsystem DFT for increasing numbers of subsystems. We investigate molecular dynamics and stationary point sampling of solvent configurations and compare the resulting (Gibbs) free energies to experiment and theoretical methods. We can show that with our hybrid model reaction barriers and reaction energies are accurately reproduced compared to experimental data.

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 20H2O → (H2O)20 and CH4 + 20H2O → 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.

An efficient and robust procedure to calculate absorption spectra of aqueous charged species applied to NO2. Phys Chem Chem Phys 2021;23:14857-14872. [PMID: 34223573 DOI: 10.1039/d1cp00652e]

Accurate calculation of absorption spectra of aqueous NO2- requires rigorously sampling the quantum potential energy surfaces for microsolvation of NO2- with at least five explicit water molecules and embedding the resulting clusters in a continuum solvent accounting for the statistical weighted contributions of individual isomers. This method, which we address as ASCEC + PCM, introduces several desired features when compared against MD simulations derived QM/MM spectra: comparatively fewer explicit solvent molecules to be treated with expensive QM methods, the identification of equilibrium structures in the quantum PES to be used in further vibrational spectroscopy, and the unequivocal identification of cluster orbitals undergoing electronic transitions and charge transfer that originate the spectral bands.

A molecular twist on hydrophobicity

A thorough exploration of the molecular basis for hydrophobicity with a comprehensive set of theoretical tools and an extensive set of organic solvent S/water binary systems is discussed in this work. Without a single exception, regardless of the nature or structure of S, all quantum descriptors of bonding yield stabilizing S⋯water interactions, therefore, there is no evidence of repulsion and thus no reason for etymological hydrophobicity at the molecular level. Our results provide molecular insight behind the exclusion of S molecules by water, customarily invoked to explain phase separation and the formation of interfaces, in favor of a complex interplay between entropic, enthalpic, and dynamic factors. S⋯water interfaces are not just thin films separating the two phases; instead, they are non-isotropic regions with density gradients for each component whose macroscopic stability is provided by a large number of very weak dihydrogen contacts. We offer a definition of interface as the region in which the density of the components in the A/B binary system is not constant. At a fundamental level, our results contribute to better current understanding of hydrophobicity.

Notwithstanding the very weak nature of individual contacts, it is the cumulative effect of a large number of interactions (green NCI surfaces) which provides macroscopic stability to the interfaces.

Towards a converged strategy for including microsolvation in reaction mechanism calculations

A major part of chemical conversions is carried out in the fluid phase, where an accurate modeling of the involved reactions requires to also take into account solvation effects. Implicit solvation models often cover these effects with sufficient accuracy but can fail drastically when specific solvent-solute interactions are important. In those cases, microsolvation, i.e., the explicit inclusion of one or more solvent molecules, is a commonly used strategy. Nevertheless, microsolvation also introduces new challenges-a consistent workflow as well as strategies how to systematically improve prediction performance are not evident. For the COSMO and COSMO-RS solvation models, this work proposes a simple protocol to decide if microsolvation is needed and how the corresponding molecular model has to look like. To demonstrate the improved accuracy of the approach, specific application examples are presented and discussed, i.e., the computation of aqueous pK_a values and a mechanistic study of the methanol mediated Morita-Baylis-Hillman reaction.

A detailed look at the bonding interactions in the microsolvation of monoatomic cations

Global and local descriptors of the properties of intermolecular bonding, formally derived from independent methodologies (QTAIM, NCI, NBO, density differences) afford a highly complex picture of the bonding interactions responsible for microsolvation of monoatomic cations. In all cases, the dominant factor dictating geometries and interaction strengths is the electrophilic power of the metal cation. The formal charge disrupts the hydrogen bonding network otherwise present in pristine water clusters, making the hydrogen bonds considerably stronger, even inducing some degree of covalency. All MO interactions are highly ionic, with strengths than in some cases approach that of the reference LiCl bond. Accumulation of electron density in the region connecting MO is observed, thus, ionic bonding in the microsolvation of monoatomic cations is not as simple as an electrostatic interaction between opposing charges.

Dissection of the Mechanism of the Wittig Reaction

A wide variety of descriptors of the evolution of bonding, rooted in the formalism of quantum mechanics, but otherwise conceptually and methodologically independent of each other (based on the quantum theory of atoms in molecules and natural bond orbitals), consistently indicate that in the mechanism of the salt-free Wittig reaction, regardless of the nature of the ylide, regardless of the nature of the transition state, and regardless of the positioning of the substituents around the reactive center, the degree of advance in the formation of the emerging C-C bond as early as at the transition state for the oxaphosphetane formation step is firmly tied to the stereochemistry of the final alkene. In addition to the fast evolution of the emerging C-C bond, very early in the reaction, a long range, weak interaction between a lone pair in the oxygen atom of the carbonyl group and an empty p orbital in the phosphorous atom, resulting from the polarization of the P═C bond in the ylide (n_O → π_P═C^*), clamps the P═C and C═O bonds to the positions required for the subsequent formation of oxaphosphetanes, thus explaining the formation of cyclic intermediates rather than betaines. Each step of the Wittig reaction is a highly asynchronous process.

Antioxidant activity of curcuminoids

An exploration of the antioxidant power of curcumin, demethoxycurcumin, and bisdemethoxycurcumin, three natural antioxidants found in Curcuma longa, is reported in this work. We exhaust all structural possibilities leading to intramolecular hydrogen bonding and evaluate 15 isomers in total. Calculations were carried out in the gas phase and in the presence of solvents (water, to mimic biological media, and ethanol, to reproduce experimental assays) following the hydrogen atom transfer (HAT) and single electron transfer (SET) mechanisms. CH3OH-O hydrogen bonds are directly related to the antioxidant power via both mechanisms. We provide evidence to explain the experimental observations and to understand the fundamental factors driving antioxidant activity from a molecular perspective. Noticeably, the solvent enhances the antioxidant power in every case. All structures considered here are predicted to have better antioxidant abilities than phenol, and come very close to or surpass vitamin E.

Fe3+chelating quinoline–hydrazone hybrids with proven cytotoxicity, leishmanicidal, and trypanocidal activities

Neglected tropical diseases cause great concern in developing countries where there are millions of reported infected humans. Our calculations support a direct relationship between biological activity and the Fe³⁺chelating ability of the shown set of quinoline–hydrazone hybrids.