Human milk fat substitutes (HMFSs) were synthesized by lipozyme RM IM-catalyzed acidolysis of chemically interesterified palm stearin (mp = 58 °C) with mixed FAs from rapeseed oil, sunflower oil, palm kernel oil, stearic acid, and myristic acid in a solvent-free system. Response surface methodology (RSM) was used to model and optimize the reactions, and the factors chosen were reaction time, temperature, substrate molar ratio, and enzyme load. The optimal conditions generated from the models were as follows: reaction time, 3.4 h; temperature, 57 °C; substrate molar ratio, 14.6 mol/mol; and enzyme load, 10.7 wt % (by the weight of total substrates). Under these conditions, the contents of palmitic acid (PA) and PA at sn-2 position (sn-2 PA) were 29.7 and 62.8%, respectively, and other observed FAs were all within the range of FAs of HMF. The product was evaluated by the cited model, and a high score (85.8) was obtained, which indicated a high degree of similarity of the product to HMF.