Artificial pancreas system (APS) is a viable option to treat diabetic patients. Researchers, however, have not conclusively determined the best control method for APS. Due to intra-/inter-variability of insulin absorption and action, an individualized algorithm is required to control blood glucose level (BGL) for each patient. To this end, we developed model predictive control (MPC) based on artificial neural networks (ANNs), which combines ANN for BGL prediction based on inputs and MPC for BGL control based on the ANN (NN-MPC). First, we developed a mathematical model for diabetic rats, which was used to identify individual virtual subjects by fitting to empirical data collected through an APS, including BGL data, insulin injection, and food intake. Then, the virtual subjects were used to generate datasets for training ANNs. The NN-MPC determines control actions (insulin injection) based on BGL predicted by the ANN. To evaluate the NN-MPC, we conducted experiments using four virtual subjects under three different scenarios. Overall, the NN-MPC maintained BGL within the normal range about 90% of the time with a mean absolute deviation of 4.7 mg/dl from a desired BGL. Our findings suggest that the NN-MPC can provide subject-specific BGL control in conjunction with a closed-loop APS. Graphical abstract ᅟ.