Yan_2021_Environ.Sci.Technol__

Quantitative structure-activity relationship (QSAR) modeling can be used to predict the toxicity of ionic liquids (ILs), but most QSAR models have been constructed by arbitrarily selecting one machine learning method and ignored the overall interactions between ILs and biological systems, such as proteins. In order to obtain more reliable and interpretable QSAR models and reveal the related molecular mechanism, we performed a systematic analysis of acetylcholinesterase (AChE) inhibition by 153 ILs using machine learning and molecular modeling. Our results showed that more reliable and stable QSAR models (R(2) > 0.85 for both cross-validation and external validation) were obtained by combining the results from multiple machine learning approaches. In addition, molecular docking results revealed that the cations and organic anions of ILs bound to specific amino acid residues of AChE through noncovalent interactions such as Pi interactions and hydrogen bonds. The calculation results of binding free energy showed that an electrostatic interaction (deltaE(ele) < -285 kJ/mol) was the main driving force for the binding of ILs to AChE. The overall findings from this investigation demonstrate that a systematic approach is much more convincing. Future research in this direction will help design the next generation of biosafe ILs.