The local electronic structure of glycine in neutral, basic, and acidic aqueous solution is studied experimentally by X-ray photoelectron spectroscopy and theoretically by molecular dynamics simulations accompanied by first-principle electronic structure and spectrum calculations. Measured and computed nitrogen and carbon is binding energies are assigned to different local atomic environments, which are shown to be sensitive to the protonation/deprotonation of the amino and carboxyl functional groups at different pH values. We report the first accurate computation of core-level chemical shifts of an aqueous solute in various protonation states and explicitly show how the distributions of photoelectron binding energies (core-level peak widths) are related to the details of the hydrogen bond configurations, i.e. the geometries of the water solvation shell and the associated electronic screening. The comparison between the experiments and calculations further enables the separation of protonation-induced (covalent) and solvent-induced (electrostatic) screening contributions to the chemical shifts in the aqueous phase. The present core-level line shape analysis facilitates an accurate interpretation of photoelectron spectra from larger biomolecular solutes than glycine.