A series of systematic DFT calculations were conducted on Keggin [SiW(9)M(3)O(40)](n-), M = Mo, V, and Nb; and Wells-Dawson anions [P(2)M(18)O(62)],(6-) M = W and Mo; [P(2)M(15)M(3)'O(62)](m-), M = W and Mo, M' = W, Mo, and V to analyze the redox properties and the basicity of the external oxygen sites in polyoxometalates with nonequivalent addenda metals. The energy and composition of the lowest unoccupied orbitals, formally delocalized over the addenda atoms, determine the redox properties of a polyoxometalate. When a Mo(6+) substitutes one W(6+) in the 1:12 tungstate, the energy of the LUMO decreases and the cluster is more easily reduced. The tungstoniobates behave differently because the niobium orbitals insert into the tungsten band and the reduction of [SiW(9)Nb(3)O(40)](7-) yields the blue species SiW(9)Nb(3) 1e and not the cluster SiW(9)Nb(2)Nb(IV). In Wells-Dawson structures, the polar and equatorial sites have different electron affinities and the reduction preferentially occurs in the equatorial sites. Inserting ions with larger electron affinities into the polar sites can modify this traditional conduct. Hence, the trisubstituted [P(2)W(15)V(3)O(62)](9-) anion is reduced in the vanadium polar sites. By means of molecular electrostatic potential maps and the relative energy of the various protonated forms of [SiW(9)V(3)O(40)](7-) and [SiW(9)Mo(3)O(40)](4-), we established the basicity scale: OV(2) > OMo(2) > OW(2) > OV > OW > OMo. Finally, a continuum model for the solvent enabled us to compare anions with different total charges.