We investigated the solvation of Ubiquitin (PDB id. 1UBQ Equilibrium molecular dynamics simulations were performed using the GROMACS 2023.3 software . The initial configurations of the systems were generated using Packmol . The ionic liquids (ILs) parameters were described using the virtual-site OPLS force fields , while the protein was modeled with the OPLS-AA force fields . The TIP3P model was used for water molecules .