A combined experimental and computational study of the molecular interactions between anionic ibuprofen and water

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.

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.

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.

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.

Absorption spectra of xanthines in aqueous solution: a computational study

We present a detailed computational analysis of the UV/Vis spectra of caffeine, paraxanthine and theophylline in aqueous solution. A hierarchy of solvation approaches for modeling the aqueous environment have been tested, ranging from the continuum model to the non-polarizable and polarizable quantum mechanical (QM)/molecular mechanics (MM) models, with and without the explicit inclusion of water molecules in the QM portion. The computed results are directly compared with the experimental data, thus highlighting the role of electrostatic, polarization and hydrogen boding solute-solvent interactions.

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.

Entropy drives the insertion of ibuprofen into model membranes

Entropy drives the insertion of ibuprofen into cell membranes.

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.

Microsolvation of NO3−: structural exploration and bonding analysis

A rich and complex structural diversity is uncovered in the microsolvation of the nitrate anion.

Hydrogen bonds in methane–water clusters

Water–methane clusters are stable at low temperatures as those found in Mars. Water cages enveloping methane are stable, although they present small probability to occur.

Microsolvation of methylmercury: structures, energies, bonding and NMR constants (199Hg,13C and17O)

Hartree–Fock (HF) and second order perturbation theory (MP2) calculations within the scalar and full relativistic frames were carried out in order to determine the equilibrium geometries and interaction energies between cationic methylmercury (CH₃Hg⁺) and up to three water molecules.

Understanding the azeotropic diethyl carbonate-water mixture by structural and energetic characterization of DEC(H2O)(n) heteroclusters

Diethyl carbonate (DEC) is an oxygenated fuel additive. During its synthesis through a promising green process, a DEC-water azeotrope is formed, which decreases DEC production efficiency in the gas phase. Molecular information about this system is scarce but could be of benefit in understanding (and potentially improving) the synthetic process. Therefore, we report a detailed computational study of the conformers of DEC, and their microsolvation with up to four water molecules, with the goal of understanding the observed 1:3 DEC:H2O molar ratio. The most stable DEC conformers (with mutual energy differences < 1.5 kcal mol(-1)) contribute to the energetic and structural properties of the complexes. An exhaustive stochastic exploration of each potential energy surface of DEC-(H2O)n, (where n = 1, 2, 3, 4) heteroclusters discovered 3, 8, 7, and 4 heterodimers, heterotrimers, heterotetramers, and heteropentamers, respectively, at the MP2/6-311++G(d,p) level of theory. DEC conformers and energies of the most stable structures at each heterocluster size were refined using CCSD(T)/6-311++G(d,p). Energy decomposition, electron density topology, and cooperative effects analyses were carried out to determine the relationship between the geometrical features of the heteroclusters and the non-covalent interaction types responsible for their stabilization. Our findings show that electrostatic and exchange energies are responsible for heterocluster stabilization, and also suggest a mutual weakening among hydrogen bonds when more than three water molecules are present. All described results are complementary and suggest a structural and energetic explanation at the molecular level for the experimental molar ratio of 1:3 (DEC:H2O) for the DEC-water azeotrope.

Theoretical tools to distinguish O-ylides from O-ylidic complexes in carbene–solvent interactions

Orbital interactions and bond indices are among the theoretical tools suitable to distinguishO-ylides fromO-ylidic carbene–solvent complexes.

Composition-dependent association behavior in the mixture of isopropanol and trichloromethane: a volumetric, vibration spectroscopic and quantum chemical study

Herein the intermolecular associative behaviors in the binary mixture of isopropanol and trichloromethane have been studied via a combined excess volumetric, vibration spectroscopic and quantum chemical approach.