Much work has been devoted to analysing thermodynamic models for solid dispersions with a view to identifying regions in the phase diagram where amorphous phase separation or drug recrystallization can occur. However, detailed partial differential equation non-equilibrium models that track the evolution of solid dispersions in time and space are lacking. Hence theoretical predictions for the timescale over which phase separation occurs in a solid dispersion are not available. In this paper, we address some of these deficiencies by (i) constructing a general multicomponent diffusion model for a dissolving solid dispersion; (ii) specializing the model to a binary drug/polymer system in storage; (iii) deriving an effective concentration dependent drug diffusion coefficient for the binary system, thereby obtaining a theoretical prediction for the timescale over which phase separation occurs; (iv) calculating the phase diagram for the Felodipine/HPMCAS system; and (iv) presenting a detailed numerical investigation of the Felodipine/HPMCAS system assuming a Flory-Huggins activity coefficient. The numerical simulations exhibit numerous interesting phenomena, such as the formation of polymer droplets and strings, Ostwald ripening/coarsening, phase inversion, and droplet-to-string transitions. A numerical simulation of the fabrication process for a solid dispersion in a hot melt extruder was also presented. Statement of Significance Solid dispersions are products that contain mixtures of drug and other materials e.g. polymer. These are liable to separate-out over time- a phenomenon known as phase separation. This means that it is possible the product differs both compositionally and structurally between the time of manufacture and the time it is taken by the patient, leading to poor bioavailability and so ultimately the shelf-life of the product has to be reduced. Theoretical predictions for the timescale over which phase separation occurs are not currently available. Also lacking are detailed partial differential equation non-equilibrium models that track the evolution of solid dispersions in time and space. This study addresses these issues, before presenting a detailed investigation of a particular drug-polymer system. (C) 2019 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

We propose a discrete in continuous mathematical model describing the in vitro growth process of biophsy-derived mammalian cardiac progenitor cells growing as clusters in the form of spheres (Cardiospheres). The approach is hybrid: discrete at cellular scale and continuous at molecular level. In the present model cells are subject to the self-organizing collective dynamics mechanism and, additionally, they can proliferate and differentiate, also depending on stochastic processes. The two latter processes are triggered and regulated by chemical signals present in the environment. Numerical simulations show the structure and the development of the clustered progenitors and are in a good agreement with the results obtained from in vitro experiments.

We propose a novel in silico model for computing drug release from multi-layer capsules. The diffusion problem in such heterogeneous layer-by-layer composite medium is described by a system of coupled partial differential equations, which we solve analytically using separation of variables. In addition to the conventional partitioning and mass transfer interlayer conditions, we consider a surface finite mass transfer resistance, which corresponds to the case of a coated capsule. The drug concentration in the core and through all the layers, as well as in the external release medium, is given in terms of a Fourier series that we compute numerically to describe and characterize the drug release mechanism.

The nucleotides ATP and ADP regulate many aspects of endothelial cell (EC) biology, including intracellular calcium concentrations, focal adhesion activation, cytoskeletal organization, and cellular motility. In vivo, ECs are constantly under flow, and the concentration of ATP/ADP on the EC surface is determined by the combined effects of nucleotide convective and diffusive transport as well as hydrolysis by ectonucleotidases on the EC surface. In addition, experiments have demonstrated that flow induces ATP release from the cells. Previous computational models have incorporated the above effects and thus described nucleotide concentration at the EC surface. However, it remains unclear what physical processes are responsible for nucleotide regulation. While some EC responses to flow have been shown to be directly driven by shear stress, others appear to also involve a non-negligible contribution of transport. In the present work, we develop a mathematical model and perform numerical simulations to investigate the relative contributions of shear stress and transport to nucleotide concentration at the EC surface. Because in vitro experiments are performed by using confluent cells in some cases and subconfluent cells in other cases, we also investigate the effect of cell density on the results. The outcomes of the simulations demonstrate a complex interplay between shear stress and transport such that transport has a significant contribution at certain shear stress values but not at others. The effect of transport on nucleotide concentration increases with cell density. The present findings enhance our understanding of the mechanisms that govern the regulation at the EC surface under flow. The implications of these findings for downstream responses such as cellular motility merit future investigation.

We present a general mechanistic model of mass diffusion for a composite sphere placed in a large ambient medium. The multi-layer problem is described by a system of diffusion equations coupled via interlayer boundary conditions such as those imposing a finite mass resistance at the extemal surface of the sphere. While the work is applicable to the generic problem of heat or mass transfer in a multi-layer sphere, the analysis and results are presented in the context of drug kinetics for desorbing and absorbing spherical microcapsules. We derive an analytical solution for the concentration in the sphere and in the surrounding medium that avoids any artificial truncation at a finite distance. The closed-form solution in each concentric layer is expressed in terms of a suitably-defined inverse Laplace transform that can be evaluated numerically. Concentration profiles and drug mass curves in the spherical layers and in the external environment are presented and the dependency of the solution on the mass transfer coefficient at the surface of the sphere analyzed.

In this study, we designed a novel drug-eluting coating for vascular implants consisting of a core coating of the anti-proliferative drug docetaxel ( DTX) and a shell coating of the platelet glycoprotein IIb/IIIa receptor monoclonal antibody SZ-21. The core/shell structure was sprayed onto the surface of 316L stainless steel stents using a coaxial electrospray process with the aim of creating a coating that exhibited a differential release of the two drugs. The prepared stents displayed a uniform coating consisting of nano/micro particles. In vitro drug release experiments were performed, and we demonstrated that a biphasic mathematical model was capable of capturing the data, indicating that the release of the two drugs conformed to a diffusion-controlled release system. We demonstrated that our coating was capable of inhibiting the adhesion and activation of platelets, as well as the proliferation and migration of smooth muscle cells ( SMCs), indicating its good biocompatibility and anti-proliferation qualities. In an in vivo porcine coronary artery model, the SZ-21/DTX drug-loaded hydrophobic core/hydrophilic shell particle coating stents were observed to promote re-endothelialization and inhibit neointimal hyperplasia. This core/shell particle-coated stent may serve as part of a new strategy for the differential release of different functional drugs to sequentially target thrombosis and in-stent restenosis during the vascular repair process and ensure rapid re-endothelialization in the field of cardiovascular disease.

The topic of controlled drug release has received much attention in recent years, for example in the design of tablets and in local drug delivery devices such as stents, transdermal patches and orthopaedic implants. In recent years, we have developed a series of models for such devices to describe drug release from a polymeric platform and transport in surrounding biological tissues. These works have culminated in the development of a mathematical model that demonstrates agreement with in-vivo drug release and tissue uptake data, for the case of a drug-eluting stent. If, on the one hand, these fully coupled models are indeed necessary to understand the spatio-temporal drug concentration in the surrounding environment, on the other hand it is clear that device manufacturers cannot intervene on the underlying biology. What they can control, however, are the properties of the polymeric platform to ensure the desired drug release profile is achieved. Indeed, the release profile is known to be a key predictor of device performance. Therefore, in the present work we focus instead primarily on the properties of the drug-containing coating. We consider two particular aspects of the drug coating design. Firstly, the delivery of two therapeutic agents, what we refer to as dual drug delivery. Depending on the particular application in question, it may be desirable for the drugs to be released at similar rates, or perhaps one of the drugs released rapidly with the other being eluted over a longer period of time. In the case of drug-eluting stents, for example, devices which release an anti-proliferative and a 'pro-healing' drug have been proposed, whilst a combination of the two has also been suggested. Secondly, motivated by today's advances in micro and nanotechnology, we propose variable porosity multi-layer coatings as an additional means of controlling the dual drug delivery. In this talk we present our mathematical model of dual drug delivery from a durable polymer coated device. We demonstrate how the release rate of each drug may in principle be controlled by varying the underlying microstructure of polymer coating or by changing the initial loading configuration of the two drugs . Our results show the role of the relevant material parameters used to tailor the release curves to a given application.

Many drugs currently on the market or in development are poorly water-soluble. This presents a serious challenge to the pharmaceutical industry because orally delivered drugs that are poorly soluble tend to pass through the gastrointestinal tract before they can fully dissolve, leading to poor drug bioavailability. One strategy to improve drug solubility is to use a solid dispersion. A solid dispersion typically consists of a hydrophobic drug embedded in a hydrophilic polymer matrix. When the dispersion dissolves in the stomach, drug-polymer interactions maintain the drug at supersaturated levels, thereby accelerating drug dissolution. Unfortunately, despite extensive research, the dissolution behaviour of solid dispersions is only partially understood. This makes the design of successful solid dispersions a somewhat hit and miss affair. Clearly, the construction of reliable mathematical models that capture the key interactions between the drug, polymer and solvent molecules in a dissolving solid dispersion would greatly assist with their rational design. In this presentation, we develop mathematical models describing the storage and dissolution of solid dispersions. The models consist of coupled systems of nonlinear partial differential equations. We then analyze in detail a particular problem describing a solid dispersion in storage. The drug-polymer interaction in the dispersion is modelled using Flory-Huggins theory , and we use the model to identify regimes in the model parameter space that lead to stable, metastable and unstable storage behaviour (phase separation).We illustrate the various phenomena arising using numerical simulations.

We propose a novel in-silico model for computing drug release from multi-layer capsules. The diffusion problem in such heterogeneous layer-by-layer composite medium is described by a system of coupled partial differential equations, which we solve analytically using separation of variables. In addition to the conventional partitioning and mass transfer interlayer conditions, we consider also the case of finite mass transfer resistance, which corresponds to the case of a coated capsule. The drug concentration in the core and through all layers, as well as in the external release medium, is given in terms of a Fourier series that we compute numerically to describe pharmacokinetics and to characterize the drug release mechanism.

We present a lattice Boltzmann model for charged leaky dielectric multiphase fluids in the context of electrified jet simulations, which are of interest for a number of production technologies including electrospinning. The role of nonlinear rheology on the dynamics of electrified jets is considered by exploiting the Carreau model for pseudoplastic fluids. We report exploratory simulations of charged droplets at rest and under a constant electric field, and we provide results for charged jet formation under electrospinning conditions.

The dynamics of deformable liquid-filled bodies (e.g., droplets, capsules, lipid vesicles) suspended in a fluid flow is a fascinating fundamental problem with increasing relevance for technological applications, for example, in drug delivery, or designing lab-on-chip devices. In the present chapter we review two main families of lattice Boltzmann models for multicomponent flows, their mechanical properties, and transport phenomena, with special focus on their application to biofluidic problems, such as the dynamics, merging, and breakup of microfluidic droplets and the motion of deformable membranes and vesicles under geometrical confinement.

Electrospinning is a nanotechnology process whereby an external electric field is used to accelerate and stretch a charged polymer jet, so as to produce fibers with nanoscale diameters. In quest of a further reduction in the cross section of electrified jets hence of a better control on the morphology of the resulting electrospun fibers, we explore the effects of an external rotating electric field orthogonal to the jet direction. Through intensive particle simulations, it is shown that by a proper tuning of the electric field amplitude and frequency, a reduction of up to a 30% in the aforementioned radius can be obtained, thereby opening new perspectives in the design of future ultra-thin electrospun fibers. Applications can be envisaged in the fields of nanophotonic components as well as for designing new and improved filtration materials.

We consider the mathematical modeling of the coupled mechanical-transport problem found in arterial physiopathology. The chapter provides an extensive description of the physiology and the microstructure of the arterial wall and explains how these are coupled to the molecular transport problem. The details of the most recent developments of the mechanistic descriptions of the multi-scale multi-physics are presented.

In this paper we investigate the extent to which variable porosity drug-eluting coatings can provide better control over drug release than coatings where the porosity is constant throughout. In particular, we aim to establish the potential benefits of replacing a single-layer with a two-layer coating of identical total thickness and initial drug mass. In our study, what distinguishes the layers (other than their individual thickness and initial drug loading) is the underlying microstructure, and in particular the effective porosity and the tortuosity of the material. We consider the effect on the drug release profile of varying the initial distribution of drug, the relative thickness of the layers and the relative resistance to diffusion offered by each layer's composition. Our results indicate that the contrast in properties of the two layers can be used as a means of better controlling the release, and that the quantity of drug delivered in the early stages can be modulated by varying the distribution of drug across the layers. We conclude that microstructural and loading differences between multi-layer variable porosity coatings can be used to tune the properties of the coating materials to obtain the desired drug release profile for a given application.

The topic of controlled drug release has received much attention in recent years, for example in the design of tablets and in local drug delivery devices such as stents, transdermal patches, therapeutic contact lenses and orthopaedic implants. In recent years, we have developed a series of models for such devices to describe drug release from a polymeric platform, drug transport in surrounding biological tissues and fully coupled models of them. These works have culminated in the development of the first mathematical model to demonstrate agreement with in-vivo drug release and tissue uptake data, for the case of a drug-eluting stent . If, on the one hand, these fully coupled models are indeed necessary to understand the spatio-temporal drug concentration in the surrounding environment, on the other hand it is clear that device manufacturers cannot intervene on the underlying biology. What they can control, however, are the properties of the polymeric platform to ensure the desired drug release profile is achieved. Indeed, the release profile is known to be a key predictor of device performance. Therefore, in the present work we take a step back from the fully coupled computational models and focus instead solely on the properties of the drug-containing coating. We consider two particular aspects of the drug coating design. Firstly, the delivery of two therapeutic agents, what we refer to as dual drug delivery. Depending on the particular application in question, it may be desirable for the drugs to be released at similar rates, or perhaps one of the drugs released rapidly with the other being eluted over a longer period of time. In the case of drug-eluting stents, for example, devices which release an anti-proliferative and a `pro-healing' drug have been proposed, whilst a combination of two of the early drug-eluting stent drugs - paclitaxel and sirolimus - has also been suggested. Secondly, motivated by today's advances in micro and nanotechnology, we propose variable porosity multi-layer coatings as an additional means of controlling the dual drug delivery.

In this paper we describe a theoretical mathematical model of dual drug delivery from a durable polymer coated medical device. We demonstrate how the release rate of each drug may in principle be controlled by altering the initial loading configuration of the two drugs. By varying the underlying microstructure of polymer coating, further control may be obtained, providing the opportunity to tailor the release profile of each drug for the given application.

In this paper we describe a combined in-vitro experimental and mathematical modelling approach to predict in-vivo drug-eluting stent kinetics. We coated stents with a mixture of sirolimus and a novel acrylic-based polymer in two different ratios. Our results indicate differential release kinetics between low and high drug dose formulations. Furthermore, mathematical model simulations of target receptor saturation suggest potential differences in efficacy.

We present a multi-layer mathematical model to describe the transdermal drug release from an iontophoretic system. The Nernst-Planck equation describes the basic convection-diffusion process, with the electric potential obtained by solving the Laplace's equation. These equations are complemented with suitable interface and boundary conditions in a multi-domain. The stability of the mathematical problem is discussed in different scenarios and a finite-difference method is used to solve the coupled system. Numerical experiments are included to illustrate the drug dynamics under different conditions.