Enhanced stability of the model mini‐protein in amino acid ionic liquids and their aqueous solutions

Using molecular dynamics simulations, the structure of model mini‐protein was thoroughly characterized in the imidazolium‐based amino acid ionic liquids and their aqueous solutions. Complete substitution of water by organic cations and anions further results in hindered conformational flexibility of the mini‐protein. This observation suggests that amino acid‐based ionic liquids are able to defend proteins from thermally induced denaturation. We show by means of radial distributions that the mini‐protein is efficiently solvated by both solvents due to a good mutual miscibility. Amino acid‐based anions prevail in the first coordination sphere of positively charged sites of the mini‐protein whereas water molecules prevail in the first coordination sphere of negatively charged sites of the mini‐protein. © 2015 Wiley Periodicals, Inc.     

Effects of Cationic Structure on Cellulose Dissolution in Ionic Liquids: A Molecular Dynamics Study

In recent years, great progress has been made in the dissolution of cellulose with ionic liquids (ILs). However, the mechanism of cellulose dissolution, especially the role the IL cation played in the dissolution process, has not been clearly understood. Herein, the mixtures of cellulose with a series of imidazolium‐based chloride ionic liquids and 1‐butyl‐3‐methyl pyridinium chloride ([C₄mpy]Cl) were simulated to study the effect that varying the heterocyclic structure and alkyl chain length of the IL cation has on the dissolution of cellulose. It was shown that the dissolution of cellulose in [C₄mpy]Cl is better than that in [C₄mim]Cl. For imidazolium‐based ILs, the shorter the alkyl chain is, the higher the solubility will be. In addition, an all‐atom force field for 1‐allyl‐3‐methyl imidazolium cation ([Amim]⁺) was developed, for the first time, to investigate the effect the electron‐withdrawing group within the alkyl chain of the IL cation has on the dissolution of cellulose. It was found that the interaction energy between [Amim]⁺ and cellulose was greater than that between [C₃mim]⁺ and cellulose, indicating that the presence of electron‐withdrawing group in alkyl chain of the cation enhanced the interaction between the cation and cellulose due to the increase of electronegativity of the cations. These findings are used to assess the cationic effect on the dissolution of cellulose in ILs. They are also expected to be important for rational design of novel ILs for efficient dissolution of cellulose.     

Removal of heavy metal ions from aqueous solutions by alkaline-earth metal phosphates

The possibility of utilization of calcium or magnesium phosphates of various composition for heavy and non-ferrous metal extraction from aqueous solutions has been studied. The efficiency of the phosphates in removal of Pb(II), Cr(III) and Fe(III) ions has been shown to decrease in the following sequence: Mg₃(PO₄)₂>MgNH₄PO₄>Ca₃(PO₄)₂>CaHPO₄>Ca₁₀(PO₄)₆(OH)₂ which is inverse to their hydrolytic stability series. It has been established that phosphates of non-apatite structure are capable of binding up to 12 mmol g⁻¹ of the named heavy metals by a chemical interaction. Hydroxyapatite interacts with the polyvalent metal ions via either the above mentioned or ion-exchange mechanism, depending on preparation method used for the apatite and the nature of metal.     

Ab initio molecular dynamics studies of ionic dissolution and precipitation of sodium chloride and silver chloride in water clusters, NaCl(H2O)n and AgCl(H2O)n, n = 6, 10, and 14

An ab initio molecular dynamics method was used to compare the ionic dissolution of soluble sodium chloride (NaCl) in water clusters with the highly insoluble silver chloride (AgCl). The investigations focused on the solvation structures, dynamics, and energetics of the contact ion pair (CIP) and of the solvent-separated ion pair (SSIP) in NaCl(H(2)O)(n) and AgCl(H(2)O)(n) with cluster sizes of n = 6, 10 and 14. We found that the minimum cluster size required to stabilize the SSIP configuration in NaCl(H(2)O)(n) is temperature-dependent. For n = 6, both configurations are present as two distinct local minima on the free-energy profile at 100 K, whereas SSIP is unstable at 300 K. Both configurations, separated by a low barrier (<10 kJ mol(-1)), are identifiable on the free energy profiles of NaCl(H(2)O)(n) for n = 10 and 14 at 300 K, with the Na(+)/Cl(-) pairs being internally solvated in the water cluster and the SSIP configuration being slightly higher in energy (<5 kJ mol(-1)). In agreement with the low bulk solubility of AgCl, no SSIP minimum is observed on the free-energy profiles of finite AgCl(H(2)O)(n) clusters. The AgCl interaction is more covalent in nature, and is less affected by the water solvent. Unlike NaCl, AgCl is mainly solvated on the surface in finite water clusters, and ionic dissolution requires a significant reorganization of the solvent structure.     

Structure and cooperativity of the hydrogen bonds in sodium dihydrogen triacetate

Geometry and energetics of low energy conformers of sodium dihydrogen triacetate (SDHTA) and its anion are studied using density functional theory (DFT) at the Becke, Lee‐Yang‐Parr hybrid functional (BLYP) and Becke, three‐parameter, Lee‐Yang‐Parr hybrid functional (B3LYP) levels. For both cases, two structures of comparable energy are found, which have different symmetry with respect to the two hydrogen bonds (HBs). DFT‐based Born–Oppenheimer molecular dynamics simulations are performed for SDHTA, which show that both structures are visited at room temperature conditions. The trajectory analysis further reveals that the two HBs behave anticooperative, that is, on average elongation of one HB is accompanied by a compression of the other one. This is in accord with nuclear magnetic resonance (NMR) experimental studies for a similar counter ion–dihydrogen triacetate complex. © 2012 Wiley Periodicals, Inc.     

Acid-base chemistry of omeprazole in aqueous solutions

Omeprazole is a potent anti-acid drug. Its absorption and mode of action are closely related to its prototropic behavior. In the present study, omeprazole samples from different sources and in different forms were studied spectrophotometrically to obtain pK_a values. In the neutral to alkaline pH region, two consistent pK_a values of 7.1 and 14.7 were obtained from various samples. The assignment of these pK_a values was realized by comparison with the prototropic properties of N(1)-methylated omeprazole substituted on the nitrogen at the 1-position of the benzimidazole ring, which was found to have a pK_a of 7.5. The omeprazole pK_a of 14.7 is assigned to the dissociation of the hydrogen from the 1-position of the benzimidazole ring and the pK_a of 7.1 is assigned to the dissociation from the protonated pyridine nitrogen of omeprazole. The results presented are at variance with those of earlier work.     

稀醋酸/水溶液中单个醋酸与水分子的氢键结构

在醋酸/水体系的工业分离中,溶液中的氢键对分离效率有很大影响.本文采用两种第一性原理方法,即从头算分子动力学模拟(AIMD)和量子化学计算(QCC),对由单个醋酸和不同水分子所组成聚合体的氢键相互作用进行了研究,采用极化统一模型和自洽反应场模型计算得到了聚合体在水溶液中的热力学数据.从QCC计算的气相和水溶液中的聚合自由能表明六元环在两种状态下都为最优结构,热力学数据反映出的各种结构的相对稳定性与AIMD模拟的环分布符合得相当一致.研究表明,由于存在醋酸和水分子间的氢键作用,稀醋酸/水溶液中的醋酸分离要比在浓醋酸溶液中困难得多.     

Specific heats of aqueous solutions of NaCl,NaBr, and KCl: Comparisons with related thermal properties

Earlier unpublished measurements of the specific heat capacities of aqueous NaCl, KCl, and NaBr solutions from 5 to 85°C and from 0.05m to saturation are presented. A twin calorimeter was used. A precision of nearly 1 part in 10⁴ in the specific heat capacity is claimed. The results are compared with literature values (summaries or original data) for heat capacity, heat of dilution, and activity coefficient of these salts in solution by means of a polynomial in half-integer powers of molality and temperature. It is found that our values agree well with the more recent literature values of the heat capacities. Small systematic inconsistencies between the various types of data were found.The experimental results presented here are taken from the postdoctoral work of Dr. F. W. Lamb, 1946–47, and from the master's thesis of Dr. J. E. Tanner, Indiana University (1954), obtainable from University Microfilms, Inc., Ann Arbor, Michigan, order number M-763.     

On pseudo-clathrate hydrate structure in aqueous solutions of tetrahydrofuran

Dielectric relaxation in aqueous solutions of tetrahydrofuran is reported and the results interpreted in terms of the pseudo-clathrate hydrate model of the structure of the liquid state. The lifetime of any given pseudo-clathrate configuration is 10^–11 sec.Issued as NRCC Publication No. 17436     

Characterising the electronic structure of ionic liquids: an examination of the 1-butyl-3-methylimidazolium chloride ion pair

In this paper we analyse the electronic properties of gas-phase 1-butyl-3-methylimidazolium Cl ion pairs, [C(4)C(1)im]Cl, in order to deepen our understanding of ionic liquids in general. Examination of charge densities, natural bond orbitals (NBO), and delocalised molecular orbitals computed at the B3LYP and MP2/6-31(++)G(d,p) levels have enabled us to explain a number of experimental phenomena: the relative acidity of different sites on the imidazolium ring, variations in hydrogen-bond donor and acceptor abilities, the apparent contradiction of the hydrogen-bond-donor parameters for different types of solute, the low probability of finding a Cl(-) anion at the rear of the imidazolium ring and the expansion of the imidazolium ring in the presence of a strong hydrogen-bond acceptor. The unreactive but coordinating environment and large electrochemical window have also been accounted for, as has the strong electron-donating character of the carbon atoms to the rear of the ring in associated imidazolylidenes. The electronic structure of the [C(4)C(1)im](+) cation is best described by a C(4)==C(5) double bond at the rear, and a delocalised three-centre 4 e(-) component across the front (N(1)-C(2)-N(3)) of the imidazolium ring; delocalisation between these regions is also significant. Hydrogen-bond formation is driven by Coulombic stabilisation, which compensates for an associated destabilisation of the electronic part of the system. Interactions are dominated by a large positive charge at C(2) and the build up of pi-electron density above and below the ring, particularly that associated with the double bond between C(4) and C(5). The NBO partial charges have been computed and compared with those used in a number of classical simulations.     

Geometric H/D isotope effects and cooperativity of the hydrogen bonds in porphycene.

We investigate the primary, secondary, and vicinal hydrogen/deuterium (H/D) isotope effects on the geometry of the two intramolecular hydrogen bonds in porphycene. Multidimensional potential energy surfaces describing the anharmonic motion in the vicinity of the trans isomer are calculated for the different symmetric (HH/DD) and asymmetric (HD) isotopomers. From the solution of the nuclear Schr?dinger equation the ground-state wavefunction is obtained, which is further used to determine the quantum corrections to the classical equilibrium geometries of the hydrogen bonds and thus the geometric isotope effects. In particular, it is found that the hydrogen bonds are cooperative, that is, both expand simultaneously even in the case of an asymmetric isotopic substitution. The theoretical predictions compare favorably with NMR chemical-shift data.     

Hydration of the hydrogen ion in aqueous solutions from intensity of the OD band of HDO

The OD band of HDO has been studied in partly deuterated aqueous solutions of acids: hydrochloric, sulfuric, nitric, trifluoroacetic, and perchloric. Composition of the hydrogen ion is determined. Observations are made on hydration of the perchlorate and trifluoroacetate ions.     

What Far‐Infrared Spectra Can Contribute to the Development of Force Fields for Ionic Liquids Used in Molecular Dynamics Simulations

Symbiosis: Far‐infrared spectra can be used to check the quality of force fields for molecular dynamics simulations of ionic liquids. On the other hand, MD simulations can explain the molecular basis of measured properties for this new liquid material (see picture).

     

Influence of temperature on the structure of liquid water

Molecular dynamics simulations have been carried out for liquid water at 7 different temperatures to understand the nature of hydrogen bonding at molecular level through the investigation of the effects of temperature on the geometry of water molecules. The changes in bond length and bond angle of water molecules from gaseous state to liquid state have been observed, and the change in the bond angle of water molecules in liquid against temperature has been revealed, which has not been seen in literature so far. The analysis of the radial distribution functions and the coordinate numbers shows that, on an average, each water molecule in liquid acts as both receptor and donor, and forms at least two hydrogen bonds with its neigbors. The analysis of the results also indicates that the water molecules form clusters in liquid.     

Molecular Reorientational Dynamics of the Neat Ionic Liquid 1‐Butyl‐3‐methylimidazolium Hexafluorophosphate by Measurement of 13C Nuclear Magnetic Relaxation Data

The reorientational dynamics of the ionic liquid 1butyl‐3‐methylimidazolium hexafluorophosphate ([BMIM]PF₆) were studied over a wide range of temperatures by measurement of ¹³C spin–lattice relaxation rates and NOE factors. The reorientational dynamics were evaluated by performing fits to the experimental relaxation data. Thus, the overall reorientational motion was described by a Cole–Davidson spectral density with a Vogel–Fulcher–Tammann temperature dependence of the correlation times. The reorientational motion of the butyl chain was modelled by a combination of the latter model for the overall motion with a Bloembergen–Purcell–Pound spectral density and an Arrhenius temperature dependence for the internal motion. Except for C2 in the aromatic ring, an additional reduction of the spectral density by the Lipari–Szabo model had to be employed. This reduction is a consequence of fast molecular motions before the rotational diffusion process becomes effective. The C2 atom did not exhibit this reduction, because the librational motion of the corresponding C2? H vector is severely hindered due to hydrogen bonding with the hexafluorophosphate anion. The observed dynamic features of the [BMIM]⁺ cation confirm quantum‐chemical structures obtained in a former study.     

Enthalpic interactions of -alanine and -serine in aqueous urea solutions

The enthalpies of dilution of -alanine and -serine in various aqueous urea solutions have been determined by flow microcalorimetry at 298.15 K. The homogeneous enthalpic interaction coefficients over the whole range of aqueous urea solutions have been calculated according to the excess enthalpy concept. The results were interpreted from the point of view of solute–solute interactions moderated by solvent effects.     

Effective thickness of the diffusion layer during hydrogen ion reduction in aqueous hydrochloric acid solutions

The process of mass transport during hydrogen ion reduction in aqueous hydrochloric acid solutions is examined both with and without excess supporting electrolyte. The study of this process is based on a numerical solution to a system of equations of material balance and the movement of particles in solution under the influence of forces for diffusion, migration, and convection. The homogeneous chemical reaction of water dissociation is also taken into account. The results of calculations show that a diffusion layer forms near the electrode during the passage of current in these solutions and that the effective thickness of this layer is the same at any instant for all particles participating in mass transport in solution in spite of differences in their diffusion coefficients. The value of the diffusion coefficient measured in these multicomponent solutions by the methods of chronopotentiometry and rotating disk electrode should differ little from that of hydrogen ions in spite of the fact that other particles with different diffusion coefficients participate in the mass transport.