An automated computational framework (MoDoop) was developed to predict the biopolymer solubilities in ionic liquids (ILs) on the basis of conductor-like screening model for real solvents calculations of two thermodynamic properties: logarithmic activity coefficient (ln γ) at infinite dilution and excess enthalpy (HE) of mixture. The calculation was based on the optimized two-dimensional structures of biopolymer models and ILs by searching the lowest-energy conformer and optimizing molecular geometry. Three lignin models together with one IL dataset were used to evaluate the prediction ability of the developed method. The evaluation results show that ln γ is a more reliable property to predict lignin solubilities in ILs and the p-coumaryl alcohol model is considered as the best model to represent lignin molecules. The developed MoDoop approach is efficient for rapid in silico screening of suitable ionic liquids to dissolve biopolymers.