For protein phylogenies, the model should be PROTGAMMA[JTT|WAG]F . There appears to be no off the shelf option for a plain ML calculation.