Like-Charge Attraction of Molecular Cations in Water: Subtle Balance between Interionic Interactions and Ionic Solvation Effect

Molecular-scale understanding of diluent effects on ligand assembly for metal ion separations

Use of metal-selective ligands in solvent extraction is instrumental in extraction of critical materials and recycling, yet, diluent effects on extraction performance are not well understood. Experimental and empirical solvent parameters have been proposed to correlate with extraction performance, but are often inadequate predictors. We follow the hypothesis that the diluents' primary influence on extraction efficiency is whether or not it hinders assembly of the bulky extracting ligands into a geometry necessary for metal complexation. This behavior is readily accessible with molecular dynamics (MD), where the atomistic description of molecules can be applied to arbitrary extractant-solvent molecules and their mixtures. Several simulated quantities are considered, from both pairwise and graph theoretical analyses, and compared to experimental distribution ratio data for americium extraction by TODGA in a series of inert, non-interacting diluents. These simple properties, especially the formation of closed triplets corresponding to the 3 : 1 ligand : metal stoichiometric solvate, suggest a potential predictive power of this approach. This methodology provides a path forward to comprehensively understand and predict diluent effects in more complex systems involving different extracting ligands and multi-component diluent mixtures.

Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids

The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.

Specific Proton-Donor Properties of Glycine Betaine. Metric Parameters and Enthalpy of Noncovalent Interactions in its Dimer, Water Complexes and Crystalline Hydrate

Trimethylglycine (glycine betaine, GB) is an important organic osmolyte that accumulates in various plant species in response to environmental stresses and has significant potential as a bioactive agent with low environmental impact. It is assumed that the hydration of GB is playing an important role in the protective mechanism. The hydration and aggregation properties of GB have not yet been studied in detail at the atomistic level. In this work, noncovalent interactions in the GB dimer and its complexes with water and crystalline monohydrate are studied. Depending on the object, periodic and non-periodic DFT calculations are used. Particular attention is paid to the metric parameters and enthalpies of intermolecular hydrogen bonds. The identification of noncovalent interactions is carried out by means of the Bader analysis of periodic or non-periodic electron density. The enthalpy of hydrogen bonds is estimated using the Rosenberg formula (PCCP 2 (2000) 2699). The specific proton donor properties of glycine betaine are due to its ability to form intermolecular C-H∙∙∙O bonds with the oxygen atom of a water molecule or the carboxylate group of a neighboring GB. The enthalpy of these bonds can be significantly greater than 10 kJ/mol. The water molecule that forms a hydrogen bond with the carboxylate group of GB also interacts with its CH groups through lone pairs of electrons. The C-H∙∙∙O bonds contribute up to 40% of the total entropy of the GB-water interaction, which is about 45 kJ/mol. The possibility of identifying C-H∙∙∙O bonds by the proton nuclear magnetic resonance method is discussed.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Diclofenac Ion Hydration: Experimental and Theoretical Search for Anion Pairs

Self-assembly of organic ions in aqueous solutions is a hot topic at the present time, and substances that are well-soluble in water are usually studied. In this work, aqueous solutions of sodium diclofenac are investigated, which, like most medicinal compounds, is poorly soluble in water. Classical MD modeling of an aqueous solution of diclofenac sodium showed equilibrium between the hydrated anion and the hydrated dimer of the diclofenac anion. The assignment and interpretation of the bands in the UV, NIR, and IR spectra are based on DFT calculations in the discrete-continuum approximation. It has been shown that the combined use of spectroscopic methods in various frequency ranges with classical MD simulations and DFT calculations provides valuable information on the association processes of medical compounds in aqueous solutions. Additionally, such a combined application of experimental and calculation methods allowed us to put forward a hypothesis about the mechanism of the effect of diclofenac sodium in high dilutions on a solution of diclofenac sodium.

Controlling "like-likes-like" charge attraction in hydroxy-functionalized ionic liquids by polarizability of the cations, interaction strength of the anions and varying alkyl chain length

We provide comprehensive understanding of "like-likes-like" charge attraction in hydroxy-functionalized ionic liquids (ILs) by means of infrared spectroscopy (IR), quantum chemistry and differential scanning calorimetry (DSC). We show that hydrogen bonding between cation and cation (c-c) is possible despite the repulsive forces between ions of like charge. Already at room temperature, the (c-c) hydrogen bonds can compete with the regular Coulomb-enhanced hydrogen bonds between cation and anion (c-a). For a large set of well-selected ILs, we show that "like-charge attraction" between the OH-functionalized cations is controllable by the polarizability of the cation, the interaction strength of the anion and the length of the hydroxyalkyl chain. In particular, we clarify whether tethering the OH group away from the positive charge center of the cationic ring with longer hydroxyalkyl chains compensates for unfavourable cation/anion combinations with respect to (c-c) cluster formation. For that purpose, we synthesized and characterized twelve ionic liquids including the differently polarizable cations, 1-(n-hydroxyalkyl)-1-methylpiperidinium [HOC_nMPip]⁺ and 1-(n-hydroxyalkyl)-pyridinium [HOC_nPy]⁺, as well as the weakly and strongly interacting anions, bis(trifluoromethanesulfonyl)imide [NTf₂]^- and methanesulfonate [OMs]^-, respectively. On top, we varied the hydroxyalkyl chain length (HOC_n) (n = 2-5). We systematically show how these three molecular ion parameters affect like-charge attraction. The use of polarizable cations, weakly interacting anions, and long alkyl chain tethers results in (c-c) clustering already at room temperature. Kinetic trapping is not a prerequisite for the existence of (c-c) cluster species in ILs. Moreover, we demonstrate that micro structuring affects macroscopic behavior of this type of ILs. We observed that substantial (c-c) interaction prevents ILs from crystallizing. Instead, these ILs supercool and finally form a glass.

When hydrogen bonding overcomes Coulomb repulsion: from kinetic to thermodynamic stability of cationic dimers

Quantum chemical calculations have been employed to study the kinetic and thermodynamic stability of hydroxy-functionalized 1-(3-hydroxyalkyl)pyridinium cationic dimers. For [Py-(CH2)n-OH+]2 structures with n = 2-17 we have calculated the robust local minima with clear dissociation barriers preventing their "Coulomb explosion" into separated cations. For n = 15 hydrogen bonding and dispersion forces fully compensate for the repulsive Coulomb forces between the cations allowing for the quantification of the pure hydrogen bond in the order of 20 kJ mol-1. The increasing kinetic stability even turns to thermodynamic stability with further elongated hydroxyalkyl chains. Now, quantum-type short-range attraction wins over classical long-range electrostatic repulsion resulting in negative binding energies and providing the first thermodynamically stable cationic dimers. The electronic, structural and spectroscopic signatures of the cationic dimers could be correlated to NBO parameters, supporting the existence of anti-electrostatic hydrogen bonds (AEHB) as recently suggested by Weinhold. In principle, these pure cationic dimers should be detectable in gas-phase experiments at low temperatures without the need of mediating molecules or counteranions.

Like-likes-Like: Cooperative Hydrogen Bonding Overcomes Coulomb Repulsion in Cationic Clusters with Net Charges up to Q=+6e

Quantum chemical calculations have been employed to study kinetically stable cationic clusters, wherein the monovalent cations are trapped by hydrogen bonding despite strongly repulsive electrostatic forces. We calculated linear and cyclic clusters of the hydroxy-functionalized cation N-(3-hydroxypropyl) pyridinium, commonly used as cation in ionic liquids. The largest kinetically stable cluster was a cyclic hexamer that very much resembles the structural motifs of molecular clusters, as known for water and alcohols. Surprisingly, strong cooperative hydrogen bonds overcome electrostatic repulsion and result in cationic clusters with a high net charge up to Q=+6e. The structural, spectroscopic, and electronic signatures of the cationic and related molecular clusters of 3-phenyl-1-propanol could be correlated to NBO parameters, supporting the existence of "anti-electrostatic" hydrogen bonds (AEHB), as recently suggested by Weinhold. We also showed that dispersion forces enhance the cationic cluster formation and compensate the electrostatic repulsion of one additional positive charge.

Counting charges on membrane-bound peptides

Quantifying the number of charges on peptides bound to interfaces requires reliable estimates of (i) surface coverage and (ii) surface charge, both of which are notoriously difficult parameters to obtain, especially at solid/water interfaces. Here, we report the thermodynamics and electrostatics governing the interactions of l-lysine and l-arginine octamers (Lys₈ and Arg₈) with supported lipid bilayers prepared.

Quantifying the number of charges on peptides bound to interfaces requires reliable estimates of (i) surface coverage and (ii) surface charge, both of which are notoriously difficult parameters to obtain, especially at solid/water interfaces. Here, we report the thermodynamics and electrostatics governing the interactions of l-lysine and l-arginine octamers (Lys₈ and Arg₈) with supported lipid bilayers prepared from a 9 : 1 mixture of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (sodium salt) (DMPG) from second harmonic generation (SHG) spectroscopy, quartz crystal microbalance with dissipation monitoring (QCM-D) and nanoplasmonic sensing (NPS) mass measurements, and atomistic simulations. The combined SHG/QCM-D/NPS approach provides interfacial charge density estimates from mean field theory for the attached peptides that are smaller by a factor of approximately two (0.12 ± 0.03 C m^–2 for Lys₈ and 0.10 ± 0.02 C m^–2 for Arg₈) relative to poly-l-lysine and poly-l-arginine. These results, along with atomistic simulations, indicate that the surface charge density of the supported lipid bilayer is neutralized by the attached cationic peptides. Moreover, the number of charges associated with each attached peptide is commensurate with those found in solution; that is, Lys₈ and Arg₈ are fully ionized when attached to the bilayer. Computer simulations indicate Lys₈ is more likely than Arg₈ to “stand-up” on the surface, interacting with lipid headgroups through one or two sidechains while Arg₈ is more likely to assume a “buried” conformation, interacting with the bilayer through up to six sidechains. Analysis of electrostatic potential and charge distribution from atomistic simulations suggests that the Gouy–Chapman model, which is widely used for mapping surface potential to surface charge, is semi-quantitatively valid; despite considerable orientational preference of interfacial water, the apparent dielectric constant for the interfacial solvent is about 30, due to the thermal fluctuation of the lipid–water interface.

Unveiling the Aggregation Behavior of Doxorubicin Hydrochloride in Aqueous Solution of 1-Octyl-3-methylimidazolium Chloride and the Effect of Bile Salt on These Aggregates: A Microscopic Study

In this article, we have unveiled the aggregation behavior of a potent chemotherapeutic drug, doxorubicin hydrochloride (Dox) in a well-known imidazolium based surface active ionic liquid (SAIL), 1-octyl-3-methylimidazolium chloride (C₈mimCl). The aggregates formed by Dox in C₈mimCl have been characterized using dynamic light scattering (DLS), fluorescence lifetime imaging microscopy (FLIM), high-resolution transmission electron microscopy (HR-TEM), analytical transmission electron microscopy (analytical TEM), field emission scanning electron microscopy (FESEM), atomic force microscopy (AFM), and Fourier-transform infrared spectroscopy (FTIR) measurements. It is found that Dox forms large spherical aggregates in the presence of C₈mimCl SAIL. We have also explored the driving force behind this aggregation behavior of Dox in C₈mimCl. Furthermore, it is observed that in the presence of a common bile salt, sodium cholate (NaCh), Dox/C₈mimCl spherical aggregates disrupt to form rodlike fibrillar aggregates. Therefore, formation of spherical aggregates and also its disruption into rodlike fibrillar aggregates have been performed, and this is expected to open a new scope for the design of a new generation smart drug delivery system where the drug itself aggregates to form the delivery system.

The influence of like-charge attraction on the structure and dynamics of ionic liquids: NMR chemical shifts, quadrupole coupling constants, rotational correlation times and failure of Stokes–Einstein–Debye

The formation of clusters of like-charge influences the structure and dynamics of ionic liquids.

Orientation of Methylguanidinium Ions at the Water-Air Interface

We use heterodyne-detected vibrational sum-frequency generation (HD-VSFG) to determine the orientation of the molecular plane of methylguanidinium ions at the surface of aqueous solutions. We measure the VSFG response of the symmetric and antisymmetric methyl stretch vibrations of the methylguanidinium ion with different polarization combinations. We find that for at least 50% of the methylguanidinium ions the molecular plane is at an angle >20° with respect to the surface plane. Hence, for only a minor fraction of the ions does the molecular plane have an orientation (near-)parallel to the surface plane, in contrast to the predictions of recent molecular dynamics simulation studies.

Controlling the kinetic and thermodynamic stability of cationic clusters by the addition of molecules or counterions

We discuss the stability of cationic clusters when adding molecules or counterions, and predict their occurrence in gas phase experiments.

Formation of Long, Multicenter π-[TCNE]22- Dimers in Solution: Solvation and Stability Assessed through Molecular Dynamics Simulations

Purely organic radical ions dimerize in solution at low temperature, forming long, multicenter bonds, despite the metastability of the isolated dimers. Here, we present the first computational study of these π-dimers in solution, with explicit consideration of solvent molecules and finite temperature effects. By means of force-field and ab initio molecular dynamics and free energy simulations, the structure and stability of π-[TCNE]₂^2- (TCNE=tetracyanoethylene) dimers in dichloromethane have been evaluated. Although the dimers dissociate at room temperature, they are stable at 175 K and their structure is similar to the one in the solid state, with a cofacial arrangement of the radicals at an interplanar separation of approximately 3.0 Å. The π-[TCNE]₂^2- dimers form dissociated ion pairs with the NBu₄⁺ counterions, and their first solvation shell comprises approximately 20 CH₂ Cl₂ molecules. Among them, the eight molecules distributed along the equatorial plane of the dimer play a key role in stabilizing the dimer through bridging C-H⋅⋅⋅N contacts. The calculated free energy of dimerization of TCNE^.- in solution at 175 K is -5.5 kcal mol^-1 . These results provide the first quantitative model describing the pairing of radical ions in solution, and demonstrate the key role of solvation forces on the dimerization process.

3D-RISM-MP2 Approach to Hydration Structure of Pt(II) and Pd(II) Complexes: Unusual H-Ahead Mode vs Usual O-Ahead One

Solvation of transition metal complexes with water has been one of the fundamental topics in physical and coordination chemistry. In particular, Pt(II) complexes have recently attracted considerable interest for their relation to anticancer activity in cisplatin and its analogues, yet the interaction of the water molecule and the metal center has been obscured. The challenge from a theoretical perspective remains that both the microscopic solvation effect and the dynamical electron correlation (DEC) effect have to be treated simultaneously in a reasonable manner. In this work we derive the analytical gradient for the three-dimensional reference interaction site model Møller-Plesset second order (3D-RISM-MP2) free energy. On the basis of the three-regions 3D-RISM self-consistent field (SCF) method recently proposed by us, we apply a new layer of the Z-vector method to the CP-RISM equation as well as point-charge approximation to the derivatives with respect to the density matrix elements in the RISM-CPHF equation to remarkably reduce the computational cost. This method is applied to study the interaction of H2O with the d(8) square planar transition metal complexes in aqueous solution, trans-[Pt(II)Cl2(NH3)(glycine)] (1a), [Pt(II)(NH3)4](2+) (1b), [Pt(II)(CN)4](2-) (1c), and their Pd(II) analogues 2a, 2b, and 2c, respectively, to elucidate whether the usual H2O interaction through O atom (O-ahead mode) or unusual one through H atom (H-ahead mode) is stable in these complexes. We find that the interaction energy of the coordinating water and the transition metal complex changes little when switching from gas to aqueous phase, but the solvation free energy differs remarkably between the two interaction modes, thereby affecting the relative stability of the H-ahead and O-ahead modes. Particularly, in contrast to the expectation that the O-ahead mode is preferred due to the presence of positive charges in 1b, the H-ahead mode is also found to be more stable. The O-ahead mode is found to be more stable than the H-ahead one only in 2b. The energy decomposition analysis (EDA) at the 3D-RISM-MP2 level revealed that the O-ahead mode is stabilized by the electrostatic (ES) interaction, whereas the H-ahead one is mainly stabilized by the DEC effect. The ES interaction is also responsible for the difference between the Pd(II) and Pt(II) complexes; because the electrostatic potential is more negative along the z-axis in the Pt(II) complex than in the Pd(II) one, the O-ahead mode prefers the Pd(II) complexes, whereas the H-ahead becomes predominant in the Pt(II) complexes.

Structure–mechanical property correlations of hydrogel forming β-sheet peptides

This review discusses about β-sheet peptide structure at the molecular level and the bulk mechanical properties of the corresponding hydrogels.

Like-charge ion pairs of hydronium and hydroxide in aqueous solution?

Hydroxide can form stronger like-ion pairs than hydronium in aqueous solution mostly due to its versatile coordination ability with solvents.