Unlike most matrix metalloproteinase (MMP) inhibitors, which target the conserved catalytic zinc site, [email protected](OH)22#nbsp;indirectly inhibits MMP-9 activity by binding at the ligand specificity S1′ loop. The allosteric binding makes [email protected](OH)22#nbsp;a promising inhibitor selective for MMP-9. However, the hydrophobic nature of the aromatic carbon cage may cause [email protected](OH)22#nbsp;to self-aggregate in aqueous solutions, hence weakening the binding. In this study, we designed [email protected](OH)22#nbsp;derivatives aiming at improving the binding affinity for MMP-9. Upon a mutation that substitutes a new functional group (-PO42-, -CH2CO2-, -CO2-, -NH3+, or -CONH2) for a hydroxyl group on the fullerenol surface, we calculated the changes in the binding free energy to the catalytic domain of human MMP-9 using the free energy perturbation (FEP) method. We found that the higher the net charge of the functional group, the stronger the binding. Compared with [email protected](OH)22, [email protected](OH)21(PO4)2-#nbsp;binds at least 1.5 to 2.5 kcal mol-1#nbsp;more strongly to MMP-9. The binding is specifically controlled by electrostatic interactions between the phosphate group and the charged residues at the binding site. In addition to the net charge, the binding free energy can be delicately adjusted by other factors, such as the functionalization site on [email protected](OH)22, the local environment of the putative binding site of MMP-9, and the presence of ions near the charged functional group. The results of our study shed light on the potential of developing [email protected](OH)22#nbsp;derivatives as nanodrugs for treating the pathological diseases associated with unregulated MMP-9 activity.