Fluctuations in physics observables in high energy ion collisions have been a topic of particular interest in recent years as they may provide important signals regarding the formation of quark-gluon plasma and the existence of a critical point. We provide predictions for basic initial state characteristics of proton-nucleus collisions from Glauber Monte Carlo models. The following systems were simulated and analysed: p + 12C, p + 14N, p + 63Cu, and p + 208Pb at wide energy range. We apply GLISSANDO accordingly fitted to tasks defined in this paper.