Solvation free energy of H2S in five cholinium-based amino acid ILs was computed using Bennet Acceptance Ratio (BAR) method under molecular dynamics (MD) simulation. From the predicted free energy, degree of removal of acid gases in liquids system can be measured and quantified by looking at excess chemical potential (μex) and Henry's law constant (kH). In order to obtain accurate value of predicted free energy, the optimised potential for liquid simulation (OPLS) force field used was validated against experimental density. The predicted density of five cholinium-based amino acid ILs showed a great agreement with the experimental results showing percentage error of <3.0. The lowest value of Henry's law constant (kH) of H2S was obtained by cholinium phenylalanate ChlPHY (11.20 atm) while the highest value was 44.12 atm gained by cholinium serinate ChlSRI. To support our understanding, the radial distribution function (RDF) of pure ILs and RDF for a specific atom of H2S in these ILs was evaluated. As observed, there were other factors that can affect the solubility of H2S in these ILs such as the length of alkyl chain of anion, interaction energy between cation and anion as well as the hydrogen bonding interaction between polar sites of ILs with H2S molecules. Significantly, the longer hydrocarbon chain of anion, the weaker the interaction energy between the two ionic species. This in turn increases the free volume between the two ions and consequently promotes an excellent solubility of H2S in ILs. Furthermore, strong hydrogen bond interaction between the H2S and ILs molecules also contributes to high solubility of H2S in this ionic environment. © 2019 Elsevier B.V.