While R has an established package (caret) to run simulated annealing, Python seems to be lacking such a solution. If you know of any Python package that does feature selection with simulated annealing, do let me know!