In this study, we systematically investigated the interaction of fullerenes/fullerenols with model proteins using a widely used computational docking program Autodock 4.0. We found that pi-pi interaction existed in all the proteins-fullerene/fullerenol systems investigated here, and contributed greatly to the overall interaction energy. We also found that with the increase of the carbon cage size, the binding strength between proteins and fullerenes/fullerenols increased constantly. In addition, our results show that functionalization of fullerenes with polar groups, such as hydroxyl groups, decreases the binding between proteins and fullerene derivatives. In other words, the more hydroxyl groups on fullerenols, the weaker binding between proteins and fullerenols.