Exploiting Cation Structure and Water Content in Modulating the Acidity of Ammonium Hydrogen Sulfate Protic Ionic Liquids.

In this paper, we investigated the effect of cation structure and water content on proton dissociation in alkylammonium [HSO4]- protic ionic liquids (ILs) doped with 20 wt % water and correlated this with experimental Hammett acidities. For pure systems, increased cation substitution resulted in a reduction in the number of direct anion-anion neighbors leading to larger numbers of small aggregates, which is further enhanced with addition of water. We also observed spontaneous proton dissociation from [HSO4]- to water only for primary amine-based protic ILs, preceded by the formation of an anion trimer motif. Investigation using DFT calculations revealed spontaneous proton dissociation from [HSO4]- to water can occur for each of the protic ILs investigated; however, this is dependent on the size of the anion aggregates. These findings are important in the fields of catalysis and lignocellulosic biomass, where solvent acidity is a crucial parameter in biomass fractionation and lignin chemistry.

Unveiling the Rational Development of Stimuli-Responsive Silk Fibroin-Based Ionogel Formulations.

We present an approach for the rational development of stimuli-responsive ionogels which can be formulated for precise control of multiple unique ionogel features and fill niche pharmaceutical applications. Ionogels are captivating materials, exhibiting self-healing characteristics, tunable mechanical and structural properties, high thermal stability, and electroconductivity. However, the majority of ionogels developed require complex chemistry, exhibit high viscosity, poor biocompatibility, and low biodegradability. In our work, we overcome these limitations. We employ a facile production process and strategically integrate silk fibroin, the biocompatible ionic liquids (ILs) choline acetate ([Cho][OAc]), choline dihydrogen phosphate ([Cho][DHP]), and choline chloride ([Cho][Cl]), traditional pharmaceutical excipients, and the model antiepileptic drug phenobarbital. In the absence of ILs, we failed to observe gel formation; yet in the presence of ILs, thermoresponsive ionogels formed. Systems were assessed via visual tests, transmission electron microscopy, confocal reflection microscopy, dynamic light scattering, zeta potential and rheology measurements. We formed diverse ionogels of strengths ranging between 18 and 642 Pa. Under 25 °C storage, formulations containing polyvinylpyrrolidone (PVP) showed an ionogel formation period ranging over 14 days, increasing in the order of [Cho][DHP], [Cho][OAc], and [Cho][Cl]. Formulations lacking PVP showed an ionogel formation period ranging over 32 days, increasing in the order of [Cho][OAc], [Cho][DHP] and [Cho][Cl]. By heating from 25 to 60 °C, immediately following preparation, thermoresponsive ionogels formed below 41 °C in the absence of PVP. Based on our experimental results and density functional theory calculations, we attribute ionogel formation to macromolecular crowding and confinement effects, further enhanced upon PVP inclusion. Holistically, applying our rational development strategy enables the production of ionogels of tunable physicochemical and rheological properties, enhanced drug solubility, and structural and energetic stability. We believe our rational development approach will advance the design of biomaterials and smart platforms for diverse drug delivery applications.

Ionic Liquid-Based Strategy for Predicting Protein Aggregation Propensity and Thermodynamic Stability.

Novel drug candidates are continuously being developed to combat the most life-threatening diseases; however, many promising protein therapeutics are dropped from the pipeline. During biological and industrial processes, protein therapeutics are exposed to various stresses such as fluctuations in temperature, solvent pH, and ionic strength. These can lead to enhanced protein aggregation propensity, one of the greatest challenges in drug development. Recently, ionic liquids (ILs), in particular, biocompatible choline chloride ([Cho]Cl)-based ILs, have been used to hinder stress-induced protein conformational changes. Herein, we develop an IL-based strategy to predict protein aggregation propensity and thermodynamic stability. We examine three key variables influencing protein misfolding: pH, ionic strength, and temperature. Using dynamic light scattering, zeta potential, and variable temperature circular dichroism measurements, we systematically evaluate the structural, thermal, and thermodynamic stability of fresh immunoglobin G4 (IgG4) antibody in water and 10, 30, and 50 wt % [Cho]Cl. Additionally, we conduct molecular dynamics simulations to examine IgG4 aggregation propensity in each system and the relative favorability of different [Cho]Cl-IgG4 packing interactions. We re-evaluate each system following 365 days of storage at 4 °C and demonstrate how to predict the thermodynamic properties and protein aggregation propensity over extended storage, even under stress conditions. We find that increasing [Cho]Cl concentration reduced IgG4 aggregation propensity both fresh and following 365 days of storage and demonstrate the potential of using our predictive IL-based strategy and formulations to radically increase protein stability and storage.

An experimental approach probing the conformational transitions and energy landscape of antibodies: a glimmer of hope for reviving lost therapeutic candidates using ionic liquid.

Understanding protein folding in different environmental conditions is fundamentally important for predicting protein structures and developing innovative antibody formulations. While the thermodynamics and kinetics of folding and unfolding have been extensively studied by computational methods, experimental methods for determining antibody conformational transition pathways are lacking. Motivated to fill this gap, we prepared a series of unique formulations containing a high concentration of a chimeric immunoglobin G4 (IgG4) antibody with different excipients in the presence and absence of the ionic liquid (IL) choline dihydrogen phosphate. We determined the effects of different excipients and IL on protein thermal and structural stability by performing variable temperature circular dichroism and bio-layer interferometry analyses. To further rationalise the observations of conformational changes with temperature, we carried out molecular dynamics simulations on a single antibody binding fragment from IgG4 in the different formulations, at low and high temperatures. We developed a methodology to study the conformational transitions and associated thermodynamics of biomolecules, and we showed IL-induced conformational transitions. We showed that the increased propensity for conformational change was driven by preferential binding of the dihydrogen phosphate anion to the antibody fragment. Finally, we found that a formulation containing IL with sugar, amino acids and surfactant is a promising candidate for stabilising proteins against conformational destabilisation and aggregation. We hope that ultimately, we can help in the quest to understand the molecular basis of the stability of antibodies and protein misfolding phenomena and offer new candidate formulations with the potential to revive lost therapeutic candidates.

Experimental validation of calculated atomic charges in ionic liquids.

A combination of X-ray photoelectron spectroscopy and near edge X-ray absorption fine structure spectroscopy has been used to provide an experimental measure of nitrogen atomic charges in nine ionic liquids (ILs). These experimental results are used to validate charges calculated with three computational methods: charges from electrostatic potentials using a grid-based method (ChelpG), natural bond orbital population analysis, and the atoms in molecules approach. By combining these results with those from a previous study on sulfur, we find that ChelpG charges provide the best description of the charge distribution in ILs. However, we find that ChelpG charges can lead to significant conformational dependence and therefore advise that small differences in ChelpG charges (<0.3 e) should be interpreted with care. We use these validated charges to provide physical insight into nitrogen atomic charges for the ILs probed.

NEXAFS spectroscopy of ionic liquids: experiments versus calculations.

Experimental near edge X-ray absorption fine structure (NEXAFS) spectra are reported for 12 ionic liquids (ILs) encompassing a range of chemical structures for both the sulfur 1s and nitrogen 1s edges and compared with time-dependent density functional theory (TD-DFT) calculations. The energy scales for the experimental data were carefully calibrated against literature data. Gas phase calculations were performed on lone ions, ion pairs and ion pair dimers, with a wide range of ion pair conformers considered. For the first time, it is demonstrated that TD-DFT is a suitable method for simulating NEXAFS spectra of ILs, although the number of ions included in the calculations and their conformations are important considerations. For most of the ILs studied, calculations on lone ions in the gas phase were sufficient to successfully reproduce the experimental NEXAFS spectra. However, for certain ILs - for example, those containing a protic ammonium cation - calculations on ion pairs were required to obtain a good agreement with experimental spectra. Furthermore, significant conformational dependence was observed for the protic ammonium ILs, providing insight into the predominant liquid phase cation-anion interactions. Among the 12 investigated ILs, we find that four have an excited state that is delocalised across both the cation and the anion, which has implications for any process that depends on the excited state, for example, radiolysis. Considering the collective experimental and theoretical data, we recommend that ion pairs should be the minimum number of ions used for the calculation of NEXAFS spectra of ILs.

Atomic charges of sulfur in ionic liquids: experiments and calculations.

Experimental near edge X-ray absorption fine structure (NEXAFS) spectra, X-ray photoelectron (XP) spectra and Auger electron spectra are reported for sulfur in ionic liquids (ILs) with a range of chemical structures. These values provide experimental measures of the atomic charge in each IL and enable the evaluation of the suitability of NEXAFS spectroscopy and XPS for probing the relative atomic charge of sulfur. In addition, we use Auger electron spectroscopy to show that when XPS binding energies differ by less than 0.5 eV, conclusions on atomic charge should be treated with caution. Our experimental data provides a benchmark for calculations of the atomic charge of sulfur obtained using different methods. Atomic charges were computed for lone ions and ion pairs, both in the gas phase (GP) and in a solvation model (SMD), with a wide range of ion pair conformers considered. Three methods were used to compute the atomic charges: charges from the electrostatic potential using a grid based method (ChelpG), natural bond orbital (NBO) population analysis and Bader's atoms in molecules (AIM) approach. By comparing the experimental and calculated measures of the atomic charge of sulfur, we provide an order for the sulfur atoms, ranging from the most negative to the most positive atomic charge. Furthermore, we show that both ChelpG and NBO are reasonable methods for calculating the atomic charge of sulfur in ILs, based on the agreement with both the XPS and NEXAFS spectroscopy results. However, the atomic charges of sulfur derived from ChelpG are found to display significant, non-physical conformational dependence. Only small differences in individual atomic charge of sulfur were observed between lone ion (GP) and ion pair IL(SMD) model systems, indicating that ion-ion interactions do not strongly influence individual atomic charges.

Linking the structures, free volumes, and properties of ionic liquid mixtures.

The formation of ionic liquid (IL) mixtures has been proposed as an approach to rationally fine-tune the physicochemical properties of ILs for a variety of applications. However, the effects of forming such mixtures on the resultant properties of the liquids are only beginning to be understood. Towards a more complete understanding of both the thermodynamics of mixing ILs and the effect of mixing these liquids on their structures and physicochemical properties, the spatial arrangement and free volume of IL mixtures containing the common [C4C1im]+ cation and different anions have been systematically explored using small angle X-ray scattering (SAXS), positron annihilation lifetime spectroscopy (PALS) and 129Xe NMR techniques. Anion size has the greatest effect on the spatial arrangement of the ILs and their mixtures in terms of the size of the non-polar domains and inter-ion distances. It was found that differences in coulombic attraction between oppositely charged ions arising from the distribution of charge density amongst the atoms of the anion also significantly influences these inter-ion distances. PALS and 129Xe NMR results pertaining to the free volume of these mixtures were found to strongly correlate with each other despite the vastly different timescales of these techniques. Furthermore, the excess free volumes calculated from each of these measurements were in excellent agreement with the excess volumes of mixing measured for the IL mixtures investigated. The correspondence of these techniques indicates that the static and dynamic free volume of these liquid mixtures are strongly linked. Consequently, fluxional processes such as hydrogen bonding do not significantly contribute to the free volumes of these liquids compared to the spatial arrangement of ions arising from their size, shape and coulombic attraction. Given the relationship between free volume and transport properties such as viscosity and conductivity, these results provide a link between the structures of IL mixtures, the thermodynamics of mixing and their physicochemical properties.

Doubly ionic hydrogen bond interactions within the choline chloride-urea deep eutectic solvent.

Deep eutectic solvents (DESs) are exemplars of systems with the ability to form neutral, ionic and doubly ionic H-bonds. Herein, the pairwise interactions of the constituent components of the choline chloride-urea DES are examined. Evidence is found for a tripodal CHCl doubly ionic H-bond motif. Moreover it is found that the covalency of doubly ionic H-bonds can be greater than, or comparable with, neutral and ionic examples. In contrast to many traditional solvents, an "alphabet soup" of many different types of H-bond (OHO[double bond, length as m-dash]C, NHO[double bond, length as m-dash]C, OHCl, NHCl, OHNH, CHCl, CHO[double bond, length as m-dash]C, NHOH and NHNH) can form. These H-bonds exhibit substantial flexibility in terms of number and strength. It is anticipated that H-bonding will have a significant impact on the entropy of the system and thus could play an important role in the formation of the eutectic. The 2 : 1 urea : choline-chloride eutectic point of this DES is often associated with the formation of a [Cl(urea)2](-) complexed anion. However, urea is found to form a H-bonded urea[choline](+) complexed cation that is energetically competitive with [Cl(urea)2](-). The negative charge on [Cl(urea)2](-) is found to remain localised on the chloride, moreover, the urea[choline](+) complexed cation forms the strongest H-bond studied here. Thus, there is potential to consider a urea[choline](+)·urea[Cl](-) interaction.

A structural investigation of ionic liquid mixtures.

The structures of mixtures of ionic liquids (ILs) featuring a common 1-butyl-3-methylimidazolium ([C4C1im](+)) cation but different anions have been investigated both experimentally and computationally. (1)H and (13)C NMR of the ILs and their mixtures has been performed both on the undiluted liquids and those diluted by CD2Cl2. These experiments have been complemented by quantum chemical density functional theory calculations and molecular dynamics simulations. These techniques have identified the formation of preferential interactions between H(2) of the imidazolium cation and the most strongly hydrogen bond (H-bond) accepting anion. In addition, a preference for the more weakly H-bond accepting anion to interact above the imidazolium ring through anion-π(+) interactions has been identified. The modelling of these data has identified that the magnitude of these preferences are small, of the order of only a few kJ mol(-1), for all IL mixtures. No clustering of the anions around a specific cation could be observed, indicating that these interactions arise from the reorientation of the cation within a randomly assigned network of anions. π(+)-π(+) stacking of the imidazolium cations was also studied and found to be promoted by ILs with a strong H-bond accepting anion. Stacking interactions are easily disrupted by the introduction of small proportions (<50 mol%) of a weakly coordinating anion due to their propensity to form anion-π(+) interactions. These results suggest that the formation of IL mixtures with different anions leads to subtle structural changes of much lower energy than the Coulombic ordering of ions, accounting for why most IL mixtures exhibit ideal, or nearly ideal, behaviour.

Hydrogen bonding and π-π interactions in imidazolium-chloride ionic liquid clusters.

A systematic electronic structure analysis of hydrogen bonding (H-bonding), anion-π(+) and π(+)-π(+) interactions present in [C1C1im]Cl ion-pairs (IPs) and selected [C1C1im]2Cl2 IP-dimers has been carried out. Interactions have been characterised using a combination of QTAIM, NCIPLOT, NBO and qualitative MO theory. IP-dimers form non-directional charge quadrupolar arrangements due to Coulombic interactions. These are found to associate either as clusters or as loosely associated IP-IP structures. Large conformational changes are found to occur for very little cost in energy, indicating that charge screening is essentially independent of the cation ring orientation. H-bond formation is accompanied by charge transfer and polarisation of the entire [C1C1im](+) ring. Charge transfer does not follow the same trend for the CHelpG, QTAIM and NBO methods. Weak "stacked" π(+)-π(+) interactions are stabilised in the presence of anions, which locate between and at the periphery of the rings, novel strongly bent H-bonds are also present. Primary (ring; C-H···Cl(-)) H-bonds and anion-π(+) (C(2)···Cl(-)) interactions are found to decay more rapidly with distance than secondary (aliphatic; C(M)-H···Cl(-)) H-bonds. This leads to an increase in the relative importance of secondary H-bond interactions in the IP-dimers. Moreover, rotation of the methyl groups within the "stacked" π(+)-π(+) IP-dimers facilitates the formation of (stronger) linear secondary H-bonds. Thus, compared to isolated IPs, secondary H-bonds may play an increased role within the condensed phase. Overall we find that structural fluidity is facilitated by fluctuating hydrogen bond, π(+)-π(+) and anion-π(+) interactions.

An OHD-RIKES and simulation study comparing a benzylmethylimidazolium ionic liquid with an equimolar mixture of dimethylimidazolium and benzene.

The principal difference between 1-benzyl-3-methyl-imidazolium triflimide [BzC1im][NTf2] and an equimolar mixture of benzene and dimethylimidazolium triflimide [C1C1im][NTf2] is that in the former the benzene moieties are tied to the imidazolium ring, while in the latter they move independently. We use femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES) and molecular simulations to explore some properties of these two systems. The Kerr spectra show small differences in the spectral densities; the simulations also show very similar environments for both the imidazolium rings and the phenyl or benzene parts of the molecules. The low frequency vibrational densities of states are also similar in the model systems. In order to perform the simulations we developed a model for the [BzC1im](+) cation and found that the barriers to rotation of the two parts of the molecule are low.

Hydrogen bonding in ionic liquids.

Ionic liquids (IL) and hydrogen bonding (H-bonding) are two diverse fields for which there is a developing recognition of significant overlap. Doubly ionic H-bonds occur when a H-bond forms between a cation and anion, and are a key feature of ILs. Doubly ionic H-bonds represent a wide area of H-bonding which has yet to be fully recognised, characterised or explored. H-bonds in ILs (both protic and aprotic) are bifurcated and chelating, and unlike many molecular liquids a significant variety of distinct H-bonds are formed between different types and numbers of donor and acceptor sites within a given IL. Traditional more neutral H-bonds can also be formed in functionalised ILs, adding a further level of complexity. Ab initio computed parameters; association energies, partial charges, density descriptors as encompassed by the QTAIM methodology (ρBCP), qualitative molecular orbital theory and NBO analysis provide established and robust mechanisms for understanding and interpreting traditional neutral and ionic H-bonds. In this review the applicability and extension of these parameters to describe and quantify the doubly ionic H-bond has been explored. Estimating the H-bonding energy is difficult because at a fundamental level the H-bond and ionic interaction are coupled. The NBO and QTAIM methodologies, unlike the total energy, are local descriptors and therefore can be used to directly compare neutral, ionic and doubly ionic H-bonds. The charged nature of the ions influences the ionic characteristics of the H-bond and vice versa, in addition the close association of the ions leads to enhanced orbital overlap and covalent contributions. The charge on the ions raises the energy of the Ylp and lowers the energy of the X-H σ* NBOs resulting in greater charge transfer, strengthening the H-bond. Using this range of parameters and comparing doubly ionic H-bonds to more traditional neutral and ionic H-bonds it is clear that doubly ionic H-bonds cover the full range of weak through to very strong H-bonds.

The impact of anion electronic structure: similarities and differences in imidazolium based ionic liquids.

In this paper the structural and energetic landscapes of ion-pair dimer conformers of 1,3-dimethylimidazolium based ionic liquids have been explored ([C1C1im][A])2, A = Cl(-), [NO3](-), [MeSO4](-), [OTf](-) and [BF4](-)). A common low-energy conformer has been selected for full electronic structure analysis. We have compared and contrasted each cluster based on the relative hydrogen bonding ability (ß-value) of the anion, which varies experimentally as Cl(-) > [NO3](-) ≈ [MeSO4](-) > [OTf](-) ≈ [BF4](-). Correlations between experimental ß-values, computed binding energies, charge transfer and various hydrogen bonding data have been made and outliers have been explained in terms of environmental effects present in the liquid phase. This is most evident in the structurally similar [MeSO4](-) and [OTf](-) anions that have very similar hydrogen bonding motifs, but significantly different ß-values. Moreover, detailed analysis of the cluster molecular orbitals, for each anion, reveals a subtle interplay between two modes of interaction, an in-plane traditional H-bonding and inter-planar anion-π interaction. Inter-planar anion-π interactions are particularly prominent for the [NO3](-) cluster. We have rationalized how the full range of interactions could impact on the structuring of ILs at surfaces and the effect these may have on viscosity.

Competitive pi interactions and hydrogen bonding within imidazolium ionic liquids.

In this paper we have explored the structural and energetic landscape of potential π(+)-π(+) stacked motifs, hydrogen-bonding arrangements and anion-π(+) interactions for gas-phase ion pair (IP) conformers and IP-dimers of 1,3-dimethylimidazolium chloride, [C1C1im]Cl. We classify cation-cation ring stacking as an electron deficient π(+)-π(+) interaction, and a competitive anion on-top IP motif as an anion-donor π(+)-acceptor interaction. 21 stable IP-dimers have been obtained within an energy range of 0-126 kJ mol(-1). The structures have been found to exhibit a complex interplay of structural features. We have found that low energy IP-dimers are not necessarily formed from the lowest energy IP conformers. The sampled range of IP-dimers exhibits new structural forms that cannot be recovered by examining the ion-pairs alone, moreover the IP-dimers are recovering additional key features of the local liquid structure. Including dispersion is shown to impact both the relative energy ordering and the geometry of the IPs and IP-dimers, however the impact is found to be subtle and dependent on the underlying functional.

Hydrogen bonding in 1-butyl- and 1-ethyl-3-methylimidazolium chloride ionic liquids.

A detailed investigation of hydrogen bonding in the pure ionic liquids [C4C1im]Cl and [C2C1im]Cl has been carried out using primarily molecular dynamics techniques. Analyses of the individual atom-atom pair radial distribution functions, and in particular those for C···Cl(-), have revealed that hydrogen bonding to the first methylene or methyl units of the substituent groups is important. Multiple geometric criteria for defining a hydrogen bond have been applied, and in particular the choice of the cutoff angle has been carefully examined. The interpretation of hydrogen bonding within these ionic liquids is highly angle dependent, and justification is provided for why it may be appropriate to employ a wider angle criteria than the 30° used for water or alcohol systems. The different types of hydrogen bond formed are characterized, and "top" conformations where the Cl anion resides above (or below) the imidazolium ring are investigated. The number of hydrogen bonds undertaken by each hydrogen atom (and the chloride anion) is quantified, and the propensity to form zero, one, or two hydrogen bonds is established. The effects of an increase in temperature on the static hydrogen bonding are also briefly examined.

Conformational preferences in diglycosyl disulfides: NMR and molecular modeling studies.

The conformations of several ß1→ß1' diglycosyl disulfides were investigated by NMR and computational methods. Experimental data, such as NOEs, proton-proton and proton-carbon-13 coupling constants, measured for solutions in DMSO, are in good agreement with values obtained by MD simulations in explicit DMSO. The disulfide torsion angles (C1-S-S-C1') preferentially sample values close to either +90° or -90° (+g or -g) and appear as the main metric that determines the conformational behavior of these glycomimetics. There is more conformational freedom around the C1-S and C1'-S' bonds (Φ and Ω torsions, respectively) and population cluster analysis allowed to identify up to four allowed conformational regions for each of the +g or -g forms. Population analysis of the hydroxylic group rotamers, based on proton-proton and proton-carbon-13 couplings as well as on calculated hydrogen bonding statistics, did not reveal any significant intramolecular hydrogen bonds in DMSO solution.

Using solvent binding and dielectric friction to interpret the hydration behavior of complex anions.

We investigate the hydration structure and water/ion dynamics about complex anions using a revised platinum group metal chloro-anion force field. Nanosecond atomistic molecular dynamics simulations were performed for the platinum group metal chloro-anion complexes. This investigation makes the first attempt at describing diffusion trends of polyatomic complex anions with counterions such as these using both hydrodynamic and dielectric friction properties of the anion solution. The transition metal anion complex diffusion rates are shown to be correlated to their first solvent shell radial distribution function peaks, their mean water residence times, and their solvation volumes as calculated by Voronoi tessellation of the simulation cell. The general trend is for slower diffusion rates to result from larger hydration shell volumes. This diffusion rate trend calculated from Stokes' law is best described using the solventberg approach with well-chosen effective solvated radii. However, to improve the diffusion constant estimates when they are compared with those calculated from computer simulations, the dielectric friction is required.

Experimentally consistent ion association predicted for metal solutions from free energy simulations.

The calculation of association constants from computer simulations has historically been complicated because of difficulties in validating metal ion force fields for solution simulations. Here we develop a method that produces a force field for divalent metal ions in metal sulfate solutions (i.e., Mg(2+)SO(4)(2-), Ca(2+)SO(4)(2-), Mn(2+)SO(4)(2-), Fe(2+)SO(4)(2-), Co(2+)SO(4)(2-), Ni(2+)SO(4)(2-), Cu(2+)SO(4)(2-), and Zn(2+)SO(4)(2-)). Using free energy of perturbation calculations, we are able to calibrate the potential of mean force W(r) for these metal sulfate solutions. The calibrated free energy profiles then allow us to produce association constants for contact, solvent-shared, and solvent-separated ion pairs that are in excellent agreement with available ultrasonic and dielectric spectroscopic data. This metal solution force field is accurate for the calculation of relative free energies from physical and biophysical computer simulations.