Our aim is to predict computationally the thermostability of enzymes in which mutations, proposed by machine learning algorithms, have been introduced. We are using Free Energy Perturbation (FEP) calculations to estimate the DDG of folding of mutated systems, relating the DG of folding to the enzyme melting temperature.