It is established that density functional theory (DFT) + U is a better choice compared to DFT for describing the correlated electron metal center in organometallics. The value of the Hubbard U parameter may be determined from linear response, either by considering the response of the metal site alone or by additionally considering the response of other sites in the compound. We analyze here in detail the influence of ligand shells of increasing size on the U parameter calculated from the linear response for five transition metal phthalocyanines. We show that the calculated multiple-site U ligand atoms that are mainly responsible for this difference are is larger than the single-site U by as much as 1 eV and the ligand atoms that are mainly responsible for this difference are the isoindole nitrogen atoms directly bonded to the central metal atom. This suggests that a different U value may be required for computations of chemisorbed molecules compared to physisorbed and gas-phase cases.