Several aspects of the molecular and electronic structure of biliverdin derivatives have been studied using density functional theory (DFT). The calculations have been performed for complexes of trianion (BvO2)3- and dianion [BvO(OH)]2-, derived from two tautomeric forms of biliverdin, BvO2H3 and [BvO(OH)]H2, with redox innocent metal ions: lithium(I), zinc(II), and gallium(III). One-electron-oxidized and reduced forms of each complex (cation and anion radicals) have been also considered. The molecular structures of all species investigated are characterized by a helical arrangement of tetrapyrrolic ligands with the metal ion lying in the plane formed by the two central pyrrole rings. The spin density distribution in four types of metallobiliverdin radicals--[(BvO2.)Mn+]n-2,[[BvO(OH).]Mn+]n-1 (cation radicals),[(BvO2.)Mn+]n-4,[[BvO(OH).]Mn+]n-3 (anion radicals)--has been investigated. In general, the absolute values of spin density on meso carbon atoms were larger than for the beta-carbon atoms. Sign alteration of spin density has been found for meso positions, and also for the beta-carbon atoms of at least two pyrrole rings. The calculated spin density maps accounted for the essential NMR spectroscopic features of iron biliverdin derivatives, including the considerable isotropic shifts detected for the meso resonances and shift alteration at the meso and beta-positions.