We performed extensive density functional calculations on various metallofullerene complexes and their polyanions to gain insight into novel η1 and η2[6:5] metal (M)–C60 bonding modes. For LnMC60 (L = ligand), the η1 mode is calculated to be the most stable, followed by η2[6:5] and η2[6:6] for –3 anions, in contrast to η2[6:6] >> η2[6:5] ≈ η1 for neutral cases. This observation is responsible for the transformation from η2[6:6] to η1 for LnM3C60, such as [Os3(CO)9C60], upon successive electron reductions. Our energy partitioning analysis (EPA) indicates that the π‐type character of η2[6:6] is much larger than that of η2[6:5]. An electron addition decreases the π‐type interaction of both the η2[6:6] and η2[6:5] modes by about 35 %, whereas it has little effect on σ‐type interactions. Because of the large proportion of π‐character in η2[6:6] coordination, the stability of η2[6:6] coordination decreases steeply as electron reductions continue. On the basis of the EPA results, we could explain why the reaction of [Os3(CO)8(CNR)(μ3‐η2[6:6],η2[6:6],η2[6:6]‐C60)] (R = CH2Ph) with CNR (4e donor) produces [Os3(CO)8(CNR)(μ3‐CNR)(μ3‐η1,η2[6:5],η1‐C60)]. The η1 and η2[6:5] bonding modes of M–C60 are crucial to fully understand the bonding nature of M–C60 bonds in exohedral metallofullerene complexes.