The dehydration process of the widely studied drug theophylline is unravelled by combining molecular dynamics simulations with powder X-ray diffraction and solid-state NMR measurements and calculations. Based on the simulation, we determine a high-energy short-lived metastable crystal form, and additionally suggest the existence of a fast-converting form.