Recombinant human keratinocyte growth factor (rhKGF) is a highly aggregation-prone therapeutic protein. The present study aimed to reduce aggregation propensity of rhKGF by engineering the aggregation hotspots. Initially, 21 mutants were designed based on the previously-identified aggregation-prone regions (APRs) and then four of them including mutants No. 4 (L91K, I119K), 7 (V13S, L91K), 14 (L91D, I119D), and 21 (A51E) were selected based on molecular dynamics (MD) simulations for further experimental studies. The recombinantly produced rhKGF and mutants were analyzed regarding secondary structure, thermal stability, aggregation propensity, and biological activity. Far-UV CD spectroscopy showed that the mutants have similar secondary structure with rhKGF. A51E mutant showed enhanced stability and decreased monomer loss under heat stress suggesting its reduced aggregation propensity compared to rhKGF. Mutant No. 14 showed higher stability and less aggregation tendency than mutant No. 4 indicating that only mutations decreasing pI of rhKGF are effective in reducing its aggregation tendency. All of the mutants were at least as potent as rhKGF in stimulating proliferation of MCF-7 epithelial cells. Our results identified A51E as an equally potent, more stable, and less aggregation-prone analog of rhKGF which could be a promising alternative drug candidate for the commercially available rhKGF (Palifermin).