We report a multisite cationic dummy atom model of the uranyl cation, suitable for atomistic simulation of the interaction of this ion with water and biomolecules. This reproduces geometry of model complexes with typical ligands, interactions and diffusion coefficient in water, and the structure of an adduct with a uranyl binding protein. A protocol for prediction of binding mode and strength of complexes formed with small peptides is proposed.