In this work, we introduce a simulation-based method for predicting the melting point of ionic liquids without prior knowledge of their crystal structure. We run molecular dynamics simulations of biofriendly, choline cation-based ionic liquids and apply the method to predict their melting point. The root-mean-square error of the predicted values is below 24 K. We advocate that such precision is sufficient for designing ionic liquids with relatively low melting points. The workflow for simulations is available for everyone and can be adopted for any species from the wide chemical space of ionic liquids.