Aqueous ionic liquids and their influence on peptide conformations: denaturation and dehydration mechanisms.

Low concentrated aqueous ionic liquids (ILs) and their influence on protein structures have attracted a lot of interest over the last few years. This can be mostly attributed to the fact that aqueous ILs, depending on the ion species involved, can be used as protein protectants or protein denaturants. Atomistic molecular dynamics (MD) simulations are performed in order to study the influence of different aprotic ILs on the properties of a short hairpin peptide. Our results reveal distinct binding and denaturation effects for 1-ethyl-3-methylimidazolium (EMIM) in combination with different anions, namely, chloride (CL), tetrafluoroborate (BF4) and acetate (ACE). The simulation outcomes demonstrate that the studied ILs with larger anions reveal a more pronounced accumulation behavior of the individual ion species around the peptide, which is accomplished by a stronger dehydration effect. We can relate these findings to the implications of the Kirkwood-Buff theory, which provides a thermodynamic explanation for the denaturation strength in terms of the IL accumulation behavior. The results for the spatial distribution functions, the binding energies and the local/bulk partition coefficients are in good agreement with metadynamics simulations in order to determine the energetically most stable peptide conformations. The free energy landscapes indicate a decrease of the denaturation strength in the order EMIM/ACE, EMIM/BF4 and EMIM/CL, which coincides with a decreasing size of the anion species. An analysis of the potential binding energies reveals that this effect is mainly of enthalpic nature.

ForConX: A forcefield conversion tool based on XML.

The force field conversion from one MD program to another one is exhausting and error-prone. Although single conversion tools from one MD program to another exist not every combination and both directions of conversion are available for the favorite MD programs Amber, Charmm, Dl-Poly, Gromacs, and Lammps. We present here a general tool for the force field conversion on the basis of an XML document. The force field is converted to and from this XML structure facilitating the implementation of new MD programs for the conversion. Furthermore, the XML structure is human readable and can be manipulated before continuing the conversion. We report, as testcases, the conversions of topologies for acetonitrile, dimethylformamide, and 1-ethyl-3-methylimidazolium trifluoromethanesulfonate comprising also Urey-Bradley and Ryckaert-Bellemans potentials. © 2017 Wiley Periodicals, Inc.

Molecular dynamics analysis of the effect of electronic polarization on the structure and single-particle dynamics of mixtures of ionic liquids and lithium salts.

We report a molecular dynamics study on the effect of electronic polarization on the structure and single-particle dynamics of mixtures of the aprotic ionic liquid 1-ethyl-3-methylimidazolium bis-(trifluoromethylsulfonyl)-imide ([EMIM][TFSI]) doped with a lithium salt with the same anion at 298 K and 1 bar. In particular, we analyze the effect of electron density fluctuations on radial distribution functions, velocity autocorrelation functions, cage correlation functions, mean-squared displacements, and vibrational densities of states, comparing the predictions of the quantum-chemistry-based Atomistic Polarizable Potential for Liquids, Electrolytes, & Polymers (APPLE&P) with those of its nonpolarizable version and those of the standard non-polarizable Optimized Potentials for Liquid Simulations-All Atom (OPLS-AA). We found that the structure of the mixture is scarcely modified by the fluctuations in electron charge of their constituents, but their transport properties are indeed quite drastically changed, with larger mobilities being predicted for the different species in the bulk mixtures with the polarizable force field. Specifically, the mean-squared displacements are larger for the polarizable potentials at identical time intervals and the intermediate subdiffusive plateaus are greatly reduced, so the transition to the diffusive regime takes place much earlier than in the non-polarizable media. Moreover, the correlations of the added cations inside their cages are weakened out earlier and their vibrational densities of states are slightly red-shifted, reflecting the weakening effect of the electronic polarization on the Coulomb coupling in these dense ionic media. The comparison of OPLS-AA with non-polarizable APPLE&P indicates that adding polarization to OPLS-AA is not sufficient to achieve results close to experiments.

Imidazolium-Based Lipid Analogues and Their Interaction with Phosphatidylcholine Membranes.

4,5-Dialkylated imidazolium lipid salts are a new class of lipid analogues showing distinct biological activities. The potential effects of the imidazolium lipids on artificial lipid membranes and the corresponding membrane interactions was analyzed. Therefore, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) was employed to create an established lipid monolayer model and a bilayer membrane. Mixed monolayers of DPPC and 4,5-dialkylimidazolium lipids differing by their alkyl chain length (C7, C11, and C15) were characterized by surface pressure-area (π-A) isotherms using a Wilhelmy film balance in combination with epifluorescence microscopy. Monolayer hysteresis for binary mixtures was examined by recording triplicate consecutive compression-expansion cycles. The lipid miscibility and membrane stability of DPPC/imidazolium lipids were subsequently evaluated by the excess mean molecular area (ΔAex) and the excess Gibbs free energy (ΔGex) of mixing. Furthermore, the thermotropic behavior of mixed liposomes of DPPC/imidazolium lipids was investigated by differential scanning calorimetry (DSC). The C15-imidazolium lipid (C15-IMe·HI) forms a thermodynamically favored and kinetically reversible Langmuir monolayer with DPPC and exhibits a rigidification effect on both DPPC monolayer and bilayer structures at low molar fractions (X ≤ 0.3). However, the incorporation of the C11-imidazolium lipid (C11-IMe·HI) causes the formation of an unstable and irreversible Langmuir-Gibbs monolayer with DPPC and disordered DPPC liposomes. The C7-imidazolium lipid (C7-IMe·HI) displays negligible membrane activity. To better understand these results on a molecular level, all-atom molecular dynamics (MD) simulations were performed. The simulations yield two opposing molecular mechanisms governing the different behavior of the three imidazolium lipids: a lateral ordering effect and a free volume/stretching effect. Overall, our study provides the first evidence that the membrane interaction of the C15 and C11 derivatives modulates the structural organization of lipid membranes. On the contrary, for the C7 derivative its membrane activity is too low to contribute to its earlier reported potent cytotoxicity.

Atomistic insights into deep eutectic electrolytes: the influence of urea on the electrolyte salt LiTFSI in view of electrochemical applications.

The influence of urea on the conducting salt lithium bis-(trifluoromethanesulfonyl)-imide (LiTFSI) in terms of lithium ion coordination numbers and lithium ion transport properties is studied via atomistic molecular dynamics simulations. Our results indicate that the presence of urea favors the formation of a deep eutectic electrolyte with pronounced ion conductivities which can be explained by a competition between urea and TFSI in occupying the first coordination shell around lithium ions. All simulation findings verify that high urea concentrations lead to a significant increase of ionic diffusivities and an occurrence of relatively high lithium transference numbers in good agreement with experimental results. The outcomes of our study point at the possible application of deep eutectic electrolytes as ion conducting materials in lithium ion batteries.

Properties of Apolar Solutes in Alkyl Imidazolium-Based Ionic Liquids: The Importance of Local Interactions.

The solvation and the dynamic properties of apolar model solutes in alkyl imidazolium-based ionic liquids (IL) are studied by using all-atom molecular dynamics simulations. In regards to specific IL effects, we focused on the often used 1-ethyl-3-methyl imidazolium cation in combination with the anions tetrafluoroborate, acetate, and bis(trifluoromethanesulfonyl)imide. Our findings reveal that the size of the anion crucially influences the accumulation behavior of the cations, which results in modified IL solvation properties. Deviations between the different alkyl imidazolium-based IL combinations can be also observed with regard to the results for the radial distribution functions, the number of surrounding molecules, and the molecular orientation. The analysis of the van Hove function further shows pronounced differences in the dynamic behavior of the solutes. The simulations verify that the solute mobilities are mainly influenced by the composition of the local solvent shell and the properties of the underlying Lennard-Jones interactions. Additional simulations with regard to modified short-range dispersion energies for alkyl imidazolium-based ILs validate our conclusions.

The influence of cations on lithium ion coordination and transport in ionic liquid electrolytes: a MD simulation study.

The dynamical and structural properties in two ionic liquid electrolytes (ILEs) based on 1-ethyl-3-methylimidazolium bis-(trifluoromethanesulfonyl)-imide ([emim][TFSI]) and N-methyl-N-propylpyrrolidinium bis-(trifluoromethanesulfonyl)imide([pyr13][TFSI]) were compared as a function of lithium bis-(trifluoromethanesulfonyl)-imide (LiTFSI) salt concentrations using atomistic molecular dynamics (MD) simulations. The many-body polarizable APPLE&P force field has been utilized. The influence of anion polarization on the structure of the first coordination shell of Li(+) was examined. In particular, the reduction of the oxygen of the TFSI anion (OTFSI) polarizability from 1.36 Å(3) to 1.00 Å(3) resulted in an increased fraction of the TFSI anion bidentate coordination to the Li(+). While the overall dynamics in [pyr13][TFSI]-based ILEs was slower than in [emim][TFSI]-based ILEs, the exchange of TFSI anions in and out of the first coordination shell of Li(+) was found to be faster in pyr13-based systems. The Li(+) ion transference number is higher for these systems as well. These trends can be related to the difference in interaction of TFSI with the IL cation which is stronger for pyr13 than for emim.

Peptides in the presence of aqueous ionic liquids: tunable co-solutes as denaturants or protectants?

We studied the stability of a small ß-hairpin peptide under the influence of an aqueous 1-ethyl-3-methylimidazolium acetate ([EMIM](+)[ACE](-)) solution via all-atom molecular dynamics simulations in combination with metadynamics. Our free energy results indicate a denaturation of the peptide structure in the presence of the ionic liquid which is validated by a significant broadening of the end-to-end distance. The radial distribution functions between the ions and the peptide were used for the calculation of the preferential binding coefficients in terms of the Kirkwood-Buff theory. A significant structure dependent binding behavior of acetate to the peptide was found which can be interpreted as the main reason for the denaturation of the native conformation. The outcomes of our simulations allow us to propose a simple mechanism to explain the unfolding of the peptide with regard to the specific properties of ionic liquids. Our results are in good agreement with experimental findings and demonstrate the benefits of ionic liquids as tunable co-solutes with regard to their influence on protein structural properties.

Solvent effects of 1-ethyl-3-methylimidazolium acetate: solvation and dynamic behavior of polar and apolar solutes.

We study the solvation properties of the ionic liquid 1-ethyl-3-methylimidazolium acetate ([EMIM](+)[ACE](-)) and the resulting dynamic behavior for differently charged model solutes at room temperature via atomistic molecular dynamics (MD) simulations of 300 ns length and 200 ns equilibration time. The solutes are simple model spheres which are either positively or negatively charged with a valency of one, or uncharged. The numerical findings indicate a distinct solvation behavior with the occurrence of well-pronounced solvation shells whose composition significantly depends on the charge of the solute. All the results of our simulations evidence the existence of a long-range perturbation effect in presence of the solutes. Our findings validate the dominance of electrostatic interactions with regard to unfavorable entropic ordering effects which elucidates the enthalpic character of the solvation process in ionic liquids for charged solutes.

Comparing induced point-dipoles and Drude oscillators.

Classical Molecular Dynamics simulations describing electrostatic interactions only by point charges can be augmented by the inclusion of atomic polarisabilities modelling charge flexibility. Two widely used models, Drude oscillators and induced point-dipoles, are compared in a systematic study using their respective implementations in CHARMM and AMBER. The question of necessity and importance of polarisable hydrogen atoms is raised and two implementations, in an implicit or explicit manner, are compared to the case of non-polarisable hydrogen atoms. For all these polarisability models, the strength of the respective atomic polarisabilities was incremented in steps of ten percent up to their full values. The influence of polarisability on the structure and dynamics of the ionic liquid EMIM(⊕)CF3SO, which is chosen as a test case, is studied thoroughly. Using appropriate model functions, the respective dynamical and structural data are fitted. Thus, a small set of parameters is deduced, which highlights the effect of polarisability. Generally, flexibility of the charge distribution leads to enhanced fluidity and less pronounced structure. As this usually occurs when adding a co-solvent to an ionic liquid, the inclusion of polarisability can be seen in much the same way in that it acts like an inner solvent.

A combined theoretical and experimental study of the influence of different anion ratios on lithium ion dynamics in ionic liquids.

In this paper, we investigate via experimental and simulation techniques the transport properties, in terms of total ionic conductivity and ion diffusion coefficients, of ionic liquids doped with lithium salts. They are composed of two anions, bis(fluorosulfonyl)imide (FSI) and bis(trifluoromethanesulfonyl)imide (TFSI), and two cations, N-ethyl-N-methylimidazolium (emim) and lithium ions. The comparison of the experimental results with the simulations shows very good agreement over a wide temperature range and a broad range of compositions. The addition of TFSI gives rise to the formation of lithium dimers (Li(+)-TFSI(-)-Li(+)). A closer analysis of such dimers shows that involved lithium ions move nearly as fast as single lithium ions, although they have a different coordination and much slower TFSI exchange rates.