In order to look for the origins of tautomeric preferences in neutral creatinine in vacuo, we examined prototropic conversions for model azoles, namely mono-hydroxy and mono-amino imidazoles, and also for their selected 1-methyl derivatives. All possible isomeric forms of creatinine and model compounds, resulting from intramolecular proton transfer (prototropy), conformational isomerism about –OH, and configurational isomerism about =NH, were studied in the gas phase (model of non-polar environment) by means of quantum-chemical methods. Because the bond-length alternation is a consequence of the resonance phenomenon, it was measured for all DFT-optimized structures by means of the harmonic oscillator model of electron delocalization (HOMED) index. Important HOMED analogies were discussed for investigated azoles and compared with those for previously studied cyclic azines, including pyrimidine nucleic acid bases. The internal effects were taken into account, and the stabilities of the investigated tautomers-rotamers were analyzed. Significant conclusions on the favored factors that can dictate the tautomeric preferences in creatinine were derived.