Medical implant-related infections remain notoriously difficult to treat due to the formation of bacterial biofilms. Systemic antibiotic delivery is often ineffective and antibiotic-eluting technologies remain immature in this field, at least in part due to limitations in adequately controlling the antibiotic release rate. A confounding factor is the lack of understanding of the most efficacious antibiotic release profile. In this paper, we introduce a novel theoretical framework that leverages functionally graded materials to achieve tunable, spatially controlled antibiotic delivery – addressing both of these key challenges. Specifically, we develop a new coupled nonlinear partial differential equation model that simultaneously captures antibiotic release from a functionally graded material coating and its transport dynamics within an evolving biofilm. Our results reveal that functionally graded material coatings can outperform homogeneous coatings in sustaining local antibiotic concentrations and suppressing biofilm growth. This study thus establishes functionally graded materials as a promising, previously underexplored design paradigm for infection-resistant medical implants and provides a quantitative basis for optimizing antibiotic release profiles in biofilm-prone environments.

Releasing capsules are widely employed in biomedical applications as smart carriers of therapeutic agents, including drugs and bioactive compounds. Such delivery vehicles typically consist of a loaded core, enclosed by one or multiple concentric coating strata. In this work, we extended existing mechanistic models to account for such multi-strata structures, including possible concurrent erosion of the capsule itself, and we characterized the release kinetics of the active substance into the surrounding medium. We presented a computational study of drug release from a spherical microcapsule, modeled through a non-linear diffusion equation incorporating radial asymmetric diffusion and space- and time-discontinuous coefficients, as suggested by the experimental data specifically collected for this study. The problem was solved numerically using a finite volume scheme on a grid with adaptive spatial and temporal resolution. Analytical expressions for concentration and cumulative release were derived for all strata, enabling the exploration of parameter sensitivity—such as coating permeability and internal diffusivity—on the overall release profile. The resulting release curves provide mechanistic insight into the transport processes and offer design criteria for achieving controlled release. Model predictions were benchmarked against in vitro experimental data obtained under physiologically relevant conditions, showing good agreement and validating the key features of the model. The proposed model thus serves as a practical tool for predicting the behavior of composite coated particles, supporting performance evaluation and the rational design of next-generation drug delivery systems with reduced experimental effort.

In this work, we deal with a mathematical model describing the dissolution process of irregularly shaped particles. In particular, we consider a complete dissolution model accounting for surface kinetics, convective diffusion, and relative velocity between fluid and dissolving particles, for three drugs with different solubility and wettability: theophylline, griseofulvin, and nimesulide. The possible subsequent recrystallization process in the bulk fluid is also considered. The governing differential equations are revisited in the context of the level-set method and Hamilton-Jacobi equations, then they are solved numerically. This choice allows us to deal with the simultaneous dissolution of hundreds of different polydisperse particles. We show the results of many computer simulations which investigate the impact of the particle size, shape, area/volume ratio, and the dependence of the interfacial mass transport coefficient on the surface curvature.

Advances in material design have led to the rapid development of novel materials with increasing complexity and functions in bioengineering. In particular, functionally graded materials (FGMs) offer important advantages in various fields of application. In this work, we consider a heterogeneous reaction-diffusion model for an FGM spherical drug release system that generalizes the multi-layer configuration to arbitrary spatially-variable coefficients. Our model proposes a possible form for the drug diffusivity and reaction rate functions exhibiting fixed average material properties and a drug release profile that can be continuously varied between the limiting cases of a homogeneous system (constant coefficients) and two-layer system (stepwise coefficients). A semi-analytical solution is then used to solve the model, which provides closed-form expressions for the drug concentration and drug release profiles in terms of generalized Fourier series. Our results show how the release rate of the proposed FGM drug release system can be controlled and continuously varied between a fast (homogeneous) and slow (two-layer) release while maintaining the same averaged values for the diffusivity and reaction rate.

Objective: Drug delivery from a drug-loaded device into an adjacent tissue is a complicated process involving drug transport through diffusion and advection, coupled with drug binding kinetics responsible for drug uptake in the tissue. This work presents a theoretical model to predict drug delivery from a device into a multilayer tissue, assuming linear reversible drug binding in the tissue layers. Methods: The governing mass conservation equations based on diffusion, advection and drug binding in a multilayer cylindrical geometry are written, and solved using Laplace transformation. The model is used to understand the impact of various non-dimensional parameters on the amounts of bound and unbound drug concentrations as functions of time. Results: Good agreement for special cases considered in past work is demonstrated. The effect of forward and reverse binding reaction rates on the multilayer drug binding process is studied in detail. The effect of the nature of the external boundary condition on drug binding and drug loss is also studied. For typical parameter values, results indicate that only a small fraction of drug delivered binds in the tissue. Additionally, the amount of bound drug rises rapidly with time due to early dominance of the forward reaction, reaches a maxima and then decays due to the reverse reaction. Conclusions: The general model presented here can account for other possible effects such as drug absorption within the device. Besides generalizing past work on drug delivery modeling, this work also offers analytical tools to understand and optimize practical drug delivery devices.

Osteosarcoma (OS) is a rare form of primary bone cancer, impacting approximately 3.4 × 106 individuals worldwide each year, primarily afflicting children. Given the limitations of existing cancer therapies, the emergence of nanotheranostic platforms has generated considerable research interest in recent decades. These platforms seamlessly integrate therapeutic potential of drug compounds with the diagnostic capabilities of imaging probes within a single construct. This innovation has opened avenues for enhanced drug delivery to targeted sites while concurrently enabling real-time monitoring of the vehicle's trajectory. In this study, we developed a nanotheranostic system employing the layer-by-layer (LbL) technique on a core containing doxorubicin (DOXO) and in-house synthesized carbon quantum dots. By utilizing chitosan and chondroitin sulfate as polyelectrolytes, we constructed a multilayered coating to encapsulate DOXO and docetaxel, achieving a coordinated co-delivery of both drugs. The LbL-functionalized nanoparticles exhibited an approximate size of 150 nm, manifesting a predominantly uniform and spherical morphology, with an encapsulation efficiency of 48% for both drugs. The presence of seven layers in these systems facilitated controlled drug release over time, as evidenced by in vitro release tests. Finally, the impact of the LbL-functionalized nanoparticles was evaluated on U2OS and Saos-2 osteosarcoma cells. The synergistic effect of the two drugs was found to be crucial in inducing cell death, particularly in Saos-2 cells treated with nanoparticles at concentrations higher than 10 μg/ml. Transmission electron microscopy analysis confirmed the internalization of the nanoparticles into both cell types through endocytic mechanisms, revealing an underlying mechanism of necrosis-induced cell death.

Objective There is increasing interest in simultaneous endovascular delivery of more than one drug from a drug-loaded stent into a diseased artery. There may be an opportunity to obtain a therapeutically desirable uptake profile of the two drugs over time by appropriate design of the initial drug distribution in the stent. Due to the non-linear, coupled nature of diffusion and reversible specific/non-specific binding of both drugs as well as competition between the drugs for a fixed binding site density, a comprehensive numerical investigation of this problem is critically needed. Methods This paper presents numerical computation of dual drug delivery in a stent-artery system, accounting for diffusion as well as specific and non-specific reversible binding. The governing differential equations are discretized in space, followed by integration over time using a stiff numerical solver. Three different cases of initial dual drug distribution are considered. Results For the particular case of sirolimus and paclitaxel, results show that competition for a limited non-specific binding site density and the significant difference in the forward/backward reaction coefficients play a key role in determining the nature of drug uptake. The nature of initial distribution of the two drugs in the stent is also found to influence the binding process, which can potentially be used to engineer a desirable dual drug uptake profile. Conclusions These results help improve the fundamental understanding of endovascular dual drug delivery. In addition, the numerical technique and results presented here may be helpful for designing and optimizing other drug delivery problems as well.

Dissolution of drug from its solid form to a dissolved form is an important consideration in the design and optimization of drug delivery devices, particularly owing to the abundance of emerging compounds that are extremely poorly soluble. When the solid dosage form is encapsulated, for example by the porous walls of an implant, the impact of the encapsulant drug transport properties is a further confounding issue. In such a case, dissolution and diffusion work in tandem to control the release of drug. However, the interplay between these two competing processes in the context of drug delivery is not as well understood as it is for other mass transfer problems, particularly for practical controlled-release considerations such as an encapsulant layer around the drug delivery device. To address this gap, this work presents a mathematical model that describes controlled release from a drug-loaded device surrounded by a passive porous layer. A solution for the drug concentration distribution is derived using the method of eigenfunction expansion. The model is able to track the dissolution front propagation, and predict the drug release curve during the dissolution process. The utility of the model is demonstrated through comparison against experimental data representing drug release from a cylindrical drugloaded orthopedic fixation pin, where the model is shown to capture the data very well. Analysis presented here reveals how the various geometrical and physicochemical parameters influence drug dissolution and, ultimately, the drug release profile. It is found that the non-dimensional initial concentration plays a key role in determining whether the problem is diffusion-limited or dissolution-limited, whereas the nature of the problem is largely independent of other parameters including diffusion coefficient and encapsulant thickness. We expect the model will prove to be a useful tool for those designing encapsulated drug delivery devices, in terms of optimizing the design of the device to achieve a desired drug release profile.

Functionally graded materials (FGMs), possessing properties that vary smoothly from one region to another,have been receiving increasing attention in recent years, particularly in the aerospace, automotive andbiomedical sectors. However, they have yet to reach their full potential. In this paper, we explore the potentialof FGMs in the context of drug delivery, where the unique material characteristics offer the potential of finetuningdrug-release for the desired application. Specifically, we develop a mathematical model of drug releasefrom a thin film FGM, based upon a spatially-varying drug diffusivity. We demonstrate that, depending on thefunctional form of the diffusivity (related to the material properties) a wide range of drug release profilesmay be obtained. Interestingly, the shape of these release profiles are not, in general, achievable from ahomogeneous medium with a constant diffusivity.

Endothelial cell (EC) migration is crucial for a wide range of processes including vascular wound healing, tumor angiogenesis, and the development of viable endovascular implants. We have previously demonstrated that ECs cultured on 15-?m wide adhesive line patterns exhibit three distinct migration phenotypes: (a) "running" cells that are polarized and migrate continuously and persistently on the adhesive lines with possible spontaneous directional changes, (b) "undecided" cells that are highly elongated and exhibit periodic changes in the direction of their polarization while maintaining minimal net migration, and (c) "tumbling-like" cells that migrate persistently for a certain amount of time but then stop and round up for a few hours before spreading again and resuming migration. Importantly, the three migration patterns are associated with distinct profles of cell length. Because of the impact of adenosine triphosphate (ATP) on cytoskeletal organization and cell polarization, we hypothesize that the observed diferences in EC length among the three diferent migration phenotypes are driven by diferences in intracellular ATP levels. In the present work, we develop a mathematical model that incorporates the interactions between cell length, cytoskeletal (F-actin) organization, and intracellular ATP concentration. An optimization procedure is used to obtain the model parameter values that best ft the experimental data on EC lengths. The results indicate that a minimalist model based on diferences in intracellular ATP levels is capable of capturing the diferent cell length profiles observed experimentally.

The rate of drug delivery to cells and the subsequent rate of drug metabolism are dependent on the cell membrane permeability to the drug. In some cases, tissue may be composed of different types of cells that exhibit order of magnitude differences in their membrane permeabilities. This paper presents a brief review of the components of the tissue scale three-compartment pharmacokinetic model of drug delivery to single-cell-type populations. The existing model is extended to consider tissue composed of two different cell types. A case study is presented of infusion mediated delivery of doxorubicin to a tumor that is composed of a drug reactive cell type and of a drug resistive cell type. The membrane permeabilities of the two cell types differ by an order of magnitude. A parametric investigation of the population composition is conducted and it is shown that the drug metabolism of the low permeability cells are negatively influenced by the fraction of the tissue composed of the permeable drug reactive cells. This is because when the population is composed mostly of drug permeable cells, the extracellular space is rapidly depleted of the drug. This has two compounding effects: (i) locally there is simply less drug available to the neighboring drug resistant cells, and (ii) the depletion of the drug from the extracellular space near the vessel-tissue interface leaves less drug to be transported to both cell types farther away from the vessel.

Controlled release of a drug contained in a spherical polymer capsule is of significant interest in many fields of medicine. There is growing interest in tailoring the erosion properties of the drug to help control and optimize the drug release process. Theoretical understanding of the nature of drug release from a bioerodible capsule is, therefore, important for designing effective drug delivery systems. While drug release from a fixed-radius capsule is relatively easier to model, the shrinking nature of a bioerodible capsule due to surface erosion presents several difficulties in theoretical modeling. This work presents a closed-form solution for the drug concentration distribution and drug delivery characteristics from a spherical capsule undergoing linear surface erosion. This problem is solved by a transformation that converts the moving boundary problem into a fixed boundary problem. For uniform initial drug distribution, the solution is shown to depend on a single non-dimensional parameter. The theoretical model is used to develop an understanding of the impact of varying the drug diffusion coefficient and rate of erosion on drug delivery characteristics. It is found that, in general, the nature of drug release in a bioerodible sphere is determined by a delicate balance between two simultaneously occurring processes - erosion and diffusion. This work improves the theoretical understanding of diffusion in drug delivery systems by accounting for the practical erosion phenomena, and may contribute towards the design and optimization of drug delivery systems.

Objective: Customization of the rate of drug delivered based on individual patient requirements is of paramount importance in the design of drug delivery devices. Advances in manufacturing may enable multilayer drug delivery devices with different initial drug distributions in each layer. However, a robust mathematical understanding of how to optimize such capabilities is critically needed. The objective of this work is to determine the initial drug distribution needed in a spherical drug delivery device such as a capsule in order to obtain a desired drug release profile. Methods: This optimization problem is posed as an inverse mass transfer problem, and optimization is carried out using the solution of the forward problem. Both non-erodible and erodible multilayer spheres are analyzed. Cases with polynomial forms of initial drug distribution are also analyzed. Optimization is also carried out for a case where an initial burst in drug release rate is desired, followed by a constant drug release rate. Results: More than 60% reduction in root-mean-square deviation of the actual drug release rate from the ideal constant drug release rate is reported. Typically, the optimized initial drug distribution in these cases prevents or minimizes large drug release rate at early times, leading to a much more uniform drug release overall. Conclusions: Results demonstrate potential for obtaining a desired drug delivery profile over time by carefully engineering the drug distribution in the drug delivery device. These results may help engineer devices that offer customized drug delivery by combining advanced manufacturing with mathematical optimization.

The rate of drug delivery to cells and the subsequent rate of drug metabolism are dependent on the cell membrane permeability to the drug. In some cases, tissue may be composed of different types of cells that exhibit order of magnitude differences in their membrane permeabilities. This paper presents a brief review of the components of the tissue scale three-compartment pharmacokinetic model of drug delivery to single-cell-type populations. The existing model is extended to consider tissue composed of two different cell types. A case study is presented of infusion mediated delivery of doxorubicin to a tumor that is composed of a drug reactive cell type and of a drug resistive cell type. The membrane permeabilities of the two cell types differ by an order of magnitude. A parametric investigation of the population composition is conducted and it is shown that the drug metabolism of the low permeability cells are negatively influenced by the fraction of the tissue composed of the permeable drug reactive cells. This is because when the population is composed mostly of drug permeable cells, the extracellular space is rapidly depleted of the drug. This has two compounding effects: (i) locally there is simply less drug available to the neighboring drug resistant cells, and (ii) the depletion of the drug from the extracellular space near the vessel-tissue interface leaves less drug to be transported to both cell types farther away from the vessel.

We develop a lumped parameter model to describe and predict the mass release of (absorption from) an arbitrary shaped body of any dimension in a large environment. Through the one-to-one analogy between diffusion-dominated mass transfer systems and electrical circuits we provide exact solutions in terms of averaged concentrations and mass released. An estimate of the equivalent resistance and of the release time is also given, and shown to be inversely proportional to the diffusivity. The proposed electrical analogue approach allows a time constant to be defined and provides an easy extension to a multi-layer and multi-phase cases in planar and spherical geometries. The simulation results are compared with those obtained from the solution of the corresponding analytical, numerical and experimental solutions, showing a satisfactory accuracy and a good agreement.

The prediction of drug dissolution profiles is crucial for elucidating the pharmacokinetic behaviour of drugs and the bioavailability of dosage forms. In this work, we develop a mathematical model to describe the dissolution process of irregularly shaped particles. We use a complete dissolution model that accounts for both surface kinetics and convective diffusion. The mechanistic relationship between the mass transfer coefficient and the local curvature is derived from the fundamental physical laws governing these processes. Our model theoretically shows that the dissolution rate depends nonlinearly on the surface curvature. The subsequent recrystallization process in the bulk fluid is also considered. The main result of this work is its simplicity, since only two coupled nonlinear ordinary differential equations are needed to describe the dissolution process. Another remarkable advantage is the possibility to determine the model parameters using common independent techniques, so that the importance of the wettability of solids on the dissolution process can be evaluated. Finally, the proposed model demonstrated the importance of particle shape in describing the experimental dissolution data of theophylline monohydrate.

Drug delivery carriers are considered an encouraging approach for the localized treatment of disease with minimum effect on the surrounding tissue. Particularly, layer-by-layer releasing particles have gained increasing interest for their ability to develop multifunctional systems able to control the release of one or more therapeutical drugs and biomolecules. Although experimental methods can offer the opportunity to establish cause and effect relationships, the data collection can be excessively expensive or/and time-consuming. For a better understanding of the impact of different design conditions on the drug-kinetics and release profile, properly designed mathematical models can be greatly beneficial. In this work, we develop a continuum-scale mathematical model to evaluate the transport and release of a drug from a microparticle based on an inner core covered by a polymeric shell. The present mathematical model includes the dissolution and diffusion of the drug and accounts for a mechanism that takes into consideration the drug biomolecules entrapped into the polymeric shell. We test a sensitivity analysis to evaluate the influence of changing the model conditions on the total system behavior. To prove the effectiveness of this proposed model, we consider the specific application of antibacterial treatment and calibrate the model against the data of the release profile for an antibiotic drug, metronidazole. The results of the numerical simulation show that ~85% of the drug is released in 230 h, and its release is characterized by two regimes where the drug dissolves, diffuses, and travels the external shell layer at a shorter time, while the drug is released from the shell to the surrounding medium at a longer time. Within the sensitivity analysis, the outer layer diffusivity is more significant than the value of diffusivity in the core, and the increase of the dissolution parameters causes an initial burst release of the drug. Finally, changing the shape of the particle to an ellipse produces an increased percentage of drugs released with an unchanged release time.