Simulation of the osmosis-based drug encapsulation in erythrocytes

Drug-loaded erythrocytes have been proposed for the treatment of disease. A common way to load drugs into erythrocytes is to apply osmotic shock. Currently, osmosis-based drug encapsulation is studied mainly experimentally, whereas a related theoretical model is still incomplete. In this study, a set of equations is developed to simulate the osmosis-based drug-encapsulation process. First, the modeling is validated with hemolysis rates and the drug-loaded quantities to be found in the literature. Then, the variation of the erythrocyte volume, formation of the pore on the erythrocyte membrane, and quantities of drug loaded into and hemoglobin released from erythrocytes are studied. Finally, an optimized operating condition for encapsulating drugs is proposed. The results show that the volume of erythrocytes exposed to hypotonic NaCl solution increases first and then abruptly decreases because of the pore formation; afterwards, it again increases and then decreases slowly. In the presence of the pore, the drug is loaded by diffusion, whereas the leak-induced convection goes against the loading. For an allowed 45% hemolysis rate, with a 10% hematocrit, the optimized NaCl concentration is 0.44%, the optimized time for sealing the loaded erythrocytes with hypertonic NaCl solution is at 6.5 s, and the quantity of albumin (drug) loaded is 4.5 mg/ml cells.