Current bottlenecks in the large-scale commercial use of many ionic liquids (ILs) include their high costs, low biodegradability, and often unknown toxicities. As a proactive effort to better understand the molecular mechanisms of ionic liquid toxicities, the work herein presents a comprehensive molecular simulation study on the interactions of 1-n-alkyl-3-methylimidazolium-based ILs with a phosphatidylcholine (PC) lipid bilayer. We explore the effects of increasing alkyl chain length (n = 4, 8, and 12) in the cation and anion hydrophobicity on the interactions with the lipid bilayer. Bulk atomistic molecular dynamics (MD) simulations performed at millimolar (mM) IL concentrations show spontaneous insertion of cations into the lipid bilayer regardless of the alkyl chain length and a favorable orientational preference once a cation is inserted. Cations also exhibit the ability to "flip" inside the lipid bilayer (as is common for amphiphiles) if partially inserted with an unfavorable orientation. Moreover, structural analysis of the lipid bilayer show that cationic insertion induces roughening of the bilayer surface, which may be a precursor to bilayer disruption. To overcome the limitation in the timescale of our simulations, free energies for a single IL cation and anion insertion have been determined based on potential of mean force calculations. These results show a decrease in free energy in response to both short and long alkyl chain IL cation insertion, and likewise for a single hydrophobic anion insertion, but an increase in free energy for the insertion of a hydrophilic chloride anion. Both bulk MD simulations and free energy calculations suggest that toxicity mechanisms toward biological systems are likely caused by ILs behaving as ionic surfactants. [Yoo et al., Soft Matter, 2014].

Ionic liquids (ILs) are salts that remain liquid down to low temperatures, and sometimes well below room temperature. ILs have been called “green solvents” because of their extraordinarily low vapor pressure and excellent solvation power, but ecotoxicology studies have shown that some ILs exhibit greater toxicity than traditional solvents. A fundamental understanding of the molecular mechanisms responsible for IL toxicity remains elusive. Here we show that one mode of IL toxicity on unicellular organisms is driven by swelling of the cell membrane. Cytotoxicity assays, confocal laser scanning microscopy, and molecular simulations reveal that IL cations nucleate morphological defects in the microbial cell membrane at concentrations near the half maximal effective concentration (EC50) of several microorganisms. Cytotoxicity increases with increasing alkyl chain length of the cation due to the ability of the longer alkyl chain to more easily embed in, and ultimately disrupt, the cell membrane.

Cassandra is an open source atomistic Monte Carlo software package that is effective in simulating the thermodynamic properties of fluids and solids. The different features and algorithms used in Cassandra are described, along with implementation details and theoretical underpinnings to various methods used. Benchmark and example calculations are shown, and information on how users can obtain the package and contribute to it are provided. © 2017 Wiley Periodicals, Inc.