To examine the interaction of uranyl with nitrogen containing groups of humic substances, the model complexes, L N ¼ NH 2 CH 3 , N(CH 3 ) 3 , and NC 5 H 5 in aqueous solution were studied computationally with an all-electron relativistic density functional method. Results are compared with the corresponding penta-aqua complex of uranyl. Although pyridine coordinates with about the same strength as L ¼ H 2 O, methylamine binds $10 kJ mol À1 stronger and trimethylamine $40 kJ mol À1 weaker than a fifth aqua ligand. Yet, each of these ligands L N donates about the same amount of charge to uranyl as L ¼ H 2 O. UAN bonds are $10 pm longer than the UAO bonds of the aqua ligands. From the present model results, one does not expect that, when compared with carboxyl groups, monodentate N-containing functional groups contribute significantly to uranyl complexation by humic substances.