I am attempting to predict the conformations of an enzyme after a single mutation. I have been using pymol to place the mutation into the protein and then I use the rosetta relax application to optimize the structure. Is this a valid way of performing this technique? Or should I be using rosetta scripts and repack the sidechains instead. I am worried that the relax protocol may not be able to optimize the residue if pymol puts it in a conformation that is far away from the optimal conformation.