Comprehensive computational modeling of coordination structures, thermodynamic stabilities, and luminescence spectra of uranyi-glycine-water complexes [UO2(Gly)(n)aq(m)](2+) (Gly = glycine, aq = H2O, n = 0-2, m = 0-5) in aqueous solution has been carried out using relativistic density functional approaches. The solvent is approximated by a dielectric continuum model and additional explicit water molecules. Detailed pictures are obtained by synergic combination of 4' experimental and theoretical data. The optimal equatorial coordination numbers E of uranyl are determined to be five. The energies of several complex conformations z are competitively close to each other. In non-basic solution the most probable complex forms are those with two water ligands replaced by the bidentate carboxyl groups of zwitterionic glycine. The N,O-chelation in non-basic solution is neither entropically nor enthalpically favored. The symmetric and antisymmetric stretch vibrations of the nearly linear O-U-O unit determine the luminescence features. The shapes of the vibrationally resolved experimental solution spectra are reproduced theoretically with an empirically fitted overall linewidth parameter. The calculated luminescence origins correspond to thermally populated, near-degenerate groups of the lowest electronically excited states of (3)Delta(g) and (3)Phi(g), character, originating from (U-O)sigma(u)->,(U-5f)delta(w)phi(u) configurations of the linear [OUO](2+) unit The intensity distributions of the vibrational progressions are consisten: with U-0 bond-length changes around 51/2 pm. The unusually high intensity of the short wavelength foot is explained by near-degeneracy of vibrationally and electronically excited states, and by intensity enhancement through the asymmetric O-U-O stretch mode. The combination of contemporary computational chemistry and experimental techniques leads to a detailed understanding of structures, thermodynamics, and luminescence of actinide compounds, including those with bioligands.