The highly crystalline nature of cellulose results in poor processability and solubility, necessitating the search for solvents that can efficiently dissolve this material. Thus, ionic liquids (ILs) have recently been shown to be well suited for this purpose, although the corresponding dissolution mechanism has not been studied in detail. Herein, we adopt a molecular dynamics (MD) approach to study the dissolution of model cellulose crystal structures in imidazolium-based ILs and gain deep mechanistic insights, demonstrating that dissolution involves IL penetration-induced cleavage of hydrogen bonds between cellulose molecular chains. Moreover, we reveal that in ILs with high cellulose dissolving power (powerful solvents, such as 1-allyl-3-methylimidazolium chloride and 1-ethyl-3-methylimidazolium chloride), the above molecular chains are peeled from the crystal phase and subsequently dispersed in the solvent, whereas no significant structural changes are observed in poor-dissolving-power solvents. Finally, we utilize MD trajectory analysis to show that the solubility of microcrystalline cellulose is well correlated with the number of intermolecular hydrogen bonds in cellulose crystals. The obtained results allow us to conclude that both anions and cations of high-dissolving-power ILs contribute to the stepwise breakage of hydrogen bonds between cellulose chains, whereas this breakage does not occur to a sufficient extent in poorly solubilizing ILs.