Theoretical modeling and mechanism of drug release from controlled release implants by microstructural image characterization
Embodiments disclosed herein relate to a model for predicting the release profile of a controlled release device. The implant modeling system and models disclosed herein allow the accurate prediction of a release profile for a controlled release device based on features extracted from micro-resolution imagery. The models combine microstructural features that can be extracted at the XRCT resolution, including pore volume and connectivity, using erosion-dilation image analysis. This strategy allows prediction of release curves of the controlled release device using XRCT despite its resolution limitations.
1. A computer-implemented method for estimating a release profile of a device for releasing a compound, the method comprising:
receiving an x-ray computed tomography (XRCT) image of the device, the device comprising particles of the compound embedded in a plurality of pores in a matrix;
generating an unlabeled mask from the XRCT image, the unlabeled mask identifying pores and necks connecting a set of the pores;
performing an erosion process on the unlabeled mask to generate an eroded mask, the erosion process eroding the pores and removing the necks connecting the set of the pores resulting in disconnected pores;
labeling the disconnected pores in the eroded mask;
performing a dilation process on the eroded mask to generate a dilated mask, the dilation process enlarging the labeled pores that were previously eroded and unlabeled;
extracting, from the dilated mask, a set of features including pore volumes and connectivities for the labeled pores;
retrieving a minimalistic model for estimating the release profile of the device based on the pore volume and connectivity, wherein the minimalistic model determines a correlation between a release coefficient for a Higuchi-type release and an XRCT parameter, wherein the XRCT parameter is based on the pore volumes and connectivities;
and estimating the release profile of the device based on the correlation determined using the retrieved minimalistic model and the pore volumes and connectivities.
2. The computer-implemented method of claim 1 , wherein the XRCT image is a voxel-based virtual 3D model of the device with a maximum resolution of one micron.
3. The computer-implemented method of claim 1 , wherein the connectivity for a given pore is a connectivity number representing a number of other pores the given pore is connected to.
4. The computer-implemented method of claim 1 , wherein the connectivities are a connectivity matrix identifying other pores that each pore of the labeled pores is connected to.
5. The computer-implemented method of claim 1 , further comprising:
constructing, based on the pore volumes and connectivities, an initial minimalistic model of the device and calibrating the initial minimalistic model,
wherein the retrieving of the minimalistic model for estimating the release profile of the device comprises retrieving the calibrated minimalistic model.
6. The computer-implemented method of claim 1 , further comprising adjusting one or more parameters of a process for manufacturing the device based on the estimated release profile.
7. The computer-implemented method of claim 1 , wherein the device is a controlled release implant and the compound is an active pharmaceutical ingredient.
8. A system for estimating a release profile of a device for releasing a compound, the system comprising:
an x-ray computed tomography (XRCT) imaging system configured to capture an XRCT image of the device, the device comprising particles of the compound embedded in a plurality of pores in a matrix;
at least one processor;
and a non-transitory computer readable medium comprising stored instructions, the instructions when executed by the at least one processor cause the at least one processor to:
generate an unlabeled mask from the XRCT image, the unlabeled mask identifying pores and necks connecting a set of the pores;
perform an erosion process on the unlabeled mask to generate an eroded mask, the erosion process eroding the pores and removing the necks connecting the set of the pores resulting in disconnected pores;
label the disconnected pores in the eroded mask;
perform a dilation process on the eroded mask to generate a dilated mask, the dilation process enlarging the labeled pores that were previously eroded and unlabeled;
extract, from the dilated mask, a set of features including pore volumes and connectivities for the labeled pores;
retrieve a minimalistic model for estimating the release profile of the device based on the pore volume and connectivity, wherein the minimalistic model determines a correlation between a release coefficient for a Higuchi-type release and an XRCT parameter, wherein the XRCT parameter is based on the pore volumes and connectivities;
and estimate the release profile of the device based on the correlation determined using the retrieved minimalistic model and the pore volumes and connectivities.
9. The system of claim 8 , wherein the XRCT image is a voxel-based virtual 3D model of the device with a maximum resolution of one micron.
10. The system of claim 8 , wherein the connectivity for a given pore is a connectivity number representing a number of other pores the given pore is connected to.
11. The system of claim 8 , wherein the connectivities are a connectivity matrix identifying the other pores each pore of the labeled pores is connected to.
12. The system of claim 8 , wherein the instructions when executed bythe at least one processor further cause the at least one processor to:
construct, based on the pore volumes and connectivities, an initial minimalistic model of the device and calibrate the initial minimalistic model,
wherein the retrieval of the minimalistic model for estimating the release profile of the device comprises retrieving the calibrated minimalistic model.
13. The system of claim 8 , wherein the instructions when executed by the at least one processor further cause the at least one processor to adjust one or more parameters of a process for manufacturing the device based on the estimated release profile.
14. The system of claim 8 , wherein the device is a controlled release implant and the compound is an active pharmaceutical ingredient.
15. A non-transitory computer-readable medium comprising instructions for estimating a release profile of a device for releasing a compound, the instructions, when executed by a computing system, causing the computing system to perform operations including:
receiving an x-ray computed tomography (XRCT) image of the device, the device comprising particles of the compound embedded in a plurality of pores in a matrix;
generating an unlabeled mask from the XRCT image, the unlabeled mask identifying pores and necks connecting a set of the pores;
performing an erosion process on the unlabeled mask to generate an eroded mask, the erosion process eroding the pores and removing the necks connecting the set of the pores resulting in disconnected pores;
labeling the disconnected pores in the eroded mask;
performing a dilation process on the eroded mask to generate a dilated mask, the dilation process enlarging the labeled pores that were previously eroded and unlabeled;
extracting, from the dilated mask, a set of features including pore volumes and connectivities for the labeled pores;
retrieving a minimalistic model for estimating the release profile of the device based on the pore volume and connectivity, wherein the minimalistic model determines a correlation between a release coefficient for a Higuchi-type release and an XRCT parameter, wherein the XRCT parameter is based on the pore volumes and connectivities;
and estimating the release profile of the device based on the correlation determined using the retrieved minimalistic model and the pore volumes and connectivities.
16. The non-transitory computer-readable medium of claim 15 , wherein the XRCT image is a voxel-based virtual 3D model of the device with a maximum resolution of one micron.
17. The non-transitory computer-readable medium of claim 15 , wherein the connectivity for a given pore is a connectivity number representing a number of other pores the given pore is connected to.
18. The non-transitory computer-readable medium of claim 15 , wherein the connectivities are a connectivity matrix identifying other pores that each pore of the labeled pores is connected to.
19. The non-transitory computer-readable medium of claim 15 , wherein the operations further include:
constructing, based on the pore volumes and connectivities, an initial minimalistic model of the device;
and calibrating the initial minimalistic model,
wherein the retrieving of the minimalistic model for estimating the release profile of the device comprises retrieving the calibrated minimalistic model.
20. The non-transitory computer-readable medium of claim 15 , wherein the operations further include adjusting one or more parameters of a process for manufacturing the device based on the estimated release profile.