In this work, the Bingel-Hirschaddition of diethylbromomalonate to all non-equivalent bonds of Sc 3 N@D 3h -C 78 was studied using density functional theory calculations.T he regioselectivities observed computationally allowed the proposal of aset of rules,the predictive aromaticity criteria (PAC), to identify the most reactive bonds of ag iven endohedral metallofullerene based on as imple evaluation of the cage structure.T he predictions based on the PACa re fully confirmed by both the computational and experimental exploration of the Bingel-Hirschr eaction of Sc 3 N@D 5h -C 80, thus indicating that these rules are rather general and applicable to other isolated pentagon rule endohedral metallofullerenes.Many endohedral metallofullerenes (EMFs) have been reported in the literature and they range from encapsulated single atoms to clusters of up to seven or eight atoms. [1] Computations have been crucial for the correct assignment and characterization of these compounds as they are usually obtained in low yields.[2] In EMFs,aformal charge transfer of up to six electrons from the metal cluster to the fullerene cage takes place,a nd is mainly responsible for their special properties and reactivity. [3] In the last years,m any EMFs with adjacent five-membered rings (5-MRs) have been synthesized. These EMFs do not obey the so-called isolated pentagon rule (IPR).[4] It has been recently demonstrated that the main reason behind the noncompliance of the IPR for EMFs is the aromaticity of the system:t he more aromatic the negatively charged fullerene cages are,t he more stable.[5] In addition, it has been proven that aromaticity can play ak ey role in the EMFs chemical reactivity.[5b]Thee xohedral functionalization of EMFs has been extensively studied both experimentally and computationally. [2,6] Thepreferred addition sites are usually characterized by:a )short C À Cb ond lengths,b )relatively high pyramidalization angles,a nd c) appropriately shaped LUMOs. [6c,e] Depending on the fullerene cage and metal cluster encapsulated, the addition is preferred at sites that are far from (e.g. Sc 3 N@I h -C 80 )o rc lose to (e.g. Y 3 N@I h -C 80 )t he metal cluster. [6a,b] Since the initial Bingel-Hirsch (BH) addition to Gd@C 60 , [7] many EMFs have been functionalized using the same protocol.[8] However, there are some recent experimental and computational examples of BH additions to EMFs wherein the usually employed reactivity parameters are not able to explain the experimental observations.F or instance, the BH addition to the non-IPR Sc 3 N@D 3 (6140)-C 68 [8b,9] was found to occur at a[6,6] bond close to the Sc center.In ar ecent report, Garcia-Borràs et al. [10] demonstrated that the thermodynamics of the BH addition to Gd 3 N@C s -(51365)-C 84 ,h aving one adjacent pentagon pair (APP), [11] Y 3 N@C 2 (22010)-C 78 with two APPs, [12] and Sc 3 N@D 3 (6140)-C 68 with three APPs [13] is governed by the additive local aromaticity of the rings of the final adducts.[10] However, Alegret et al. found that BH additi...