calculate binding energy of doped system

If I have a pure fullerene and I want to replace a metal instead of its carbon, what steps should I do to evaluate stability of doping system?  The next question is how to get the binding energy?

Is it right to do the following steps?

Do I first optimize the pure material, then replace the given metal with carbon in the optimized material?