At present, polyhydroxyalkanoates (PHAs) have been considered as a promising alternative to conventional plastics due to their diverse variability in structure and rapid biodegradation. To ensure cost competitiveness in the market, thermoseparating aqueous two-phase extraction (ATPE) with the advantages of being mild and environmental-friendly was suggested as the primary isolation and purification tool for PHAs. Utilizing two-level full factorial design, this work studied the influence and interaction between four independent variables on the partitioning behavior of PHAs. Based on the experimental results, feed forward neural network (FFNN) was used to develop an empirical model of PHAs based on the ATPE thermoseparating input-output parameter. In this case, bootstrap resampling technique was used to generate more data. At the conditions of 15 wt % phosphate salt, 18 wt % ethylene oxide–propylene oxide (EOPO), and pH 10 without the addition of NaCl, the purification and recovery of PHAs achieved a highest yield of 93.9%. Overall, the statistical analysis demonstrated that the phosphate concentration and thermoseparating polymer concentration were the most significant parameters due to their individual influence and synergistic interaction between them on all the response variables. The final results of the FFNN model showed the ability of the model to seamlessly generalize the relationship between the input–output of the process.