Mathematical modelling of controlled drug release from polymer micro-spheres : incorporating the effects of swelling, diffusion and dissolution via moving boundary problems

Controlled drug delivery is a key topic in modern pharmacotherapy, where controlled drug delivery devices are required to prolong the period of release, maintain a constant release rate, or release the drug with a predetermined release profile. In the pharmaceutical industry, the development process of a controlled drug delivery device may be facilitated enormously by the mathematical modelling of drug release mechanisms, directly decreasing the number of necessary experiments. Such mathematical modelling is difficult because several mechanisms are involved during the drug release process. The main drug release mechanisms of a controlled release device are based on the device’s physiochemical properties, and include diffusion, swelling and erosion. In this thesis, four controlled drug delivery models are investigated. These four models selectively involve the solvent penetration into the polymeric device, the swelling of the polymer, the polymer erosion and the drug diffusion out of the device but all share two common key features. The first is that the solvent penetration into the polymer causes the transition of the polymer from a glassy state into a rubbery state. The interface between the two states of the polymer is modelled as a moving boundary and the speed of this interface is governed by a kinetic law. The second feature is that drug diffusion only happens in the rubbery region of the polymer, with a nonlinear diffusion coefficient which is dependent on the concentration of solvent. These models are analysed by using both formal asymptotics and numerical computation, where front-fixing methods and the method of lines with finite difference approximations are used to solve these models numerically. This numerical scheme is conservative, accurate and easily implemented to the moving boundary problems and is thoroughly explained in Section 3.2. From the small time asymptotic analysis in Sections 5.3.1, 6.3.1 and 7.2.1, these models exhibit the non-Fickian behaviour referred to as Case II diffusion, and an initial constant rate of drug release which is appealing to the pharmaceutical industry because this indicates zeroorder release. The numerical results of the models qualitatively confirms the experimental behaviour identified in the literature. The knowledge obtained from investigating these models can help to develop more complex multi-layered drug delivery devices in order to achieve sophisticated drug release profiles. A multi-layer matrix tablet, which consists of a number of polymer layers designed to provide sustainable and constant drug release or bimodal drug release, is also discussed in this research. The moving boundary problem describing the solvent penetration into the polymer also arises in melting and freezing problems which have been modelled as the classical onephase Stefan problem. The classical one-phase Stefan problem has unrealistic singularities existed in the problem at the complete melting time. Hence we investigate the effect of including the kinetic undercooling to the melting problem and this problem is called the one-phase Stefan problem with kinetic undercooling. Interestingly we discover the unrealistic singularities existed in the classical one-phase Stefan problem at the complete melting time are regularised and also find out the small time behaviour of the one-phase Stefan problem with kinetic undercooling is different to the classical one-phase Stefan problem from the small time asymptotic analysis in Section 3.3. In the case of melting very small particles, it is known that surface tension effects are important. The effect of including the surface tension to the melting problem for nanoparticles (no kinetic undercooling) has been investigated in the past, however the one-phase Stefan problem with surface tension exhibits finite-time blow-up. Therefore we investigate the effect of including both the surface tension and kinetic undercooling to the melting problem for nanoparticles and find out the the solution continues to exist until complete melting. The investigation of including kinetic undercooling and surface tension to the melting problems reveals more insight into the regularisations of unphysical singularities in the classical one-phase Stefan problem. This investigation gives a better understanding of melting a particle, and contributes to the current body of knowledge related to melting and freezing due to heat conduction.

Negative nuclear Overhauser effects as probes of macromolecular structure

Nuclear Overhauser effects (NOE) and circular dichroism (CD) techniques have been used to probe @-turn conformations in acyclic and cyclic peptides containingPro-Xsequences. The model peptides studied are of the type Piv-Pro-X-NHMe (X = Aib, D-Ala, Gly, Val, and Leu) and Boc-Cys-Pro-X-C s NHMe (X = Aib, L-Ala, D-Ala, Gly, and Leu). In the acyclic series, observation of NOES between Pro C"H and X-NH, together with solvent and temperature dependence of NH chemical shifts, establishes a 4 - 1 hydrogen bond stabilized type I1 @-turn in the Gly, D-Ala, and Aib peptides, in CDC13 and (CD3)2S0. A positive n-r* CD band at -225-230 nm appears to be characteristic of this structure. For the acyclic Pro-Leu peptide the observation of NOE's for both Pro and Leu C"H resonances on saturation of Leu NH is compatible with a type V bend or consecutive y-turn conformation. In the cyclic disulfide series the Pro-Aib and Pro-D-Ala peptides favor type I1 @-turns, whereas all other peptides adopt type I (111) conformations. All the cyclic disulfides exhibit an intense negative CD band at -228-230 nm. The results suggest thatgeneralcorrelations between CD spectral type and specific 0-turn conformations may not be obtained. Evidence for solvent-dependent structural changes in the Pro-Aib sequence in both cyclic and acyclic peptides is presented.

Simultaneous observation of positive and negative nuclear Overhauser effects in oligopeptides due to segmental motio

Nuclear Overhauser effect (NOE) studies of the symmetrical cystine peptides (Formula: see text) (n = 1-3) in dimethylsulfoxide, have resulted in the simultaneous observation of both positive and negative NOEs. Positive NOEs are observed on the Trp C2H and C4H protons of the indole ring upon irradiation of Trp C alpha H and C beta H2 resonances in the peptides where n = 1 and 2. Negative NOEs are observed between backbone NH and C alpha H protons. The magnitudes of the observed NOEs are sensitive to changes in molecular size and solvent viscosity. The results demonstrate that NOEs may be a useful probe of sidechain segmental motion in oligopeptides.

Evaluation of solute-solvent interactions from solvent blue-shifts of n → π* transitions of C=O, C=S, NO2 and N=N groups: hydrogen bond energies of various donor-acceptor systems

Effects of non-polar, polar and proton-donating solvents on the n → π* transitions of C=O, C=S, NO2 and N=N groups have been investigated. The shifts of the absorption maxima in non-polar and polar solvents have been related to the electrostatic interactions between solute and solvent molecules, by employing the theory of McRAE. In solvents which can donate protons the solvent shifts are mainly determined by solute-solvent hydrogen bonding. Isobestic points have been found in the n → π* bonds of ethylenetrithio-carbonate in heptane-alcohol and heptane-chloroform solvent systems, indicating the existence of equilibria between the hydrogen bonded and the free species of the solute. Among the different proton-donating solvents studied water produces the largest blue-shifts. The blue-shifts in alcohols decrease in the order 2,2,2-trifluoroethanol, methanol, ethanol, isopropanol and t-butanol, the blue-shift in trifluoroethanol being nearly equal to that in water. This trend is exactly opposite to that for the self-association of alcohols. It is suggested that electron-withdrawing groups not merely decrease the extent of self-association of alcohols, but also increase the ability to donate hydrogen bonds. The approximate hydrogen-bond energies for several donor-acceptor systems have been estimated. In a series of aliphatio ketones and nitro compounds studied, the blue-shifts and consequently the hydrogen bond energies decrease with the decrease in the electron-withdrawing power of the alkyl groups. It is felt that electron-withdrawing groups render the chromophores better proton acceptors, and the alcohols better donors. A linear relationship between n → π* transition frequency and the infrared frequency of ethylenetrithiocarbonate has been found. It is concluded that stabilization of the electronic ground states of solute molecules by electrostatic and/or hydrogen-bond interactions determines the solvent shifts.

Inertial effects in solvation dynamics

The dynamics of solvation of newly created charged species in dense dipolar liquids can proceed at a high speed with time constants often in the subpicosecond domain. The motion of the solvent molecules can be in the inertial limit at such short times. In this paper we present a microscopic study of the effects of inertial motion of solvent molecules on the solvation dynamics of a newly created ion in a model dipolar liquid. Interesting dynamical behavior emerges when the relative contribution of the translational modes in the wave-vector-dependent longitudinal relaxation time is significant. Especially, the theory predicts that the time correlation function of the solvation energy can become oscillatory in some limiting situations. In general, the dynamics becomes faster in the presence of the inertial contribution. We discuss the experimental situations where the inertial effects can be noticeable.

Effects of molecular size in solvation dynamics

The effects of molecular size on the dynamics of polar solvation are studied by using a microscopic theory which includes the translational relaxation modes of the solvent consistently. It is shown that while in the absence of the translational contribution the solvation rate increases with the size of the solute (in agreement with the conclusions of the nonequilibrium MSA theory),a complete reversal of the solute size dependence occurs when translational modes make a significant contribution to the solvent polarization relaxation.

Dielectric Relaxation Spectroscopy for Evaluation of the Influence of Solvent Dynamics on Ion Transport in Succinonitrile-Salt Plastic Crystalline Electrolytes

Influence of succinonitrile (SN) dynamics on ion transport in SN-lithium perchlorate (LiClO4) electrolytes is discussed here via dielectric relaxation spectroscopy. Dielectric relaxation spectroscopy (similar to 2 x 10(-3) Hz to 3 MHz) of SN and SN-LiClO4 was studied as a function of salt content (up to 7 mol % or 1 M) and temperature (-20 to +60 degrees C). Analyses of real and imaginary parts of permittivity convincingly reveal the influence Of trans gauche isomerism and solvent-salt association (solvation) effects on ion transport. The relaxation processes are highly dependent on the salt concentration and temperature. While pristine SN display only intrinsic dynamics (i.e., trans-gauche isomerism) which enhances with an increase in temperature, SN-LiClO4 electrolytes especially at high salt concentrations (similar to 0.04-1 M) show salt-induced relaxation processes. In the concentrated electrolytes, the intrinsic dynamics was observed to be a function of salt content, becoming faster with an increase in salt concentration. Deconvolution of the imaginary part of the permittivity spectra using Havriliak-Negami (HN) function show a relaxation process corresponding to the above phenomena. The permittivity data analyzed using HN and Kohlrausch-Williams-Watta (KWW) functions show non-Debye relaxation processes and enhancement in the trans phase (enhanced solvent dynamics) as a function of salt concentration and temperature.

Effects of Ultrafast Solvation on the Rate of Adiabatic Outer-Sphere Electron Transfer Reactions

Recent studies have demonstrated that solvation dynamics in many common dipolar liquids contain an initial, ultrafast Gaussian component which may contribute even more than 60% to the total solvation energy. It is also known that adiabatic electron transfer reactions often probe the high-frequency components of the relevant solvent friction (Hynes, J. T. J. Phys. Chem. 1986, 90, 3701). In this paper, we present a theoretical study of the effects of the ultrafast solvent polar modes on the adiabatic electron transfer reactions by using the formalism of Hynes. Calculations have been carried out for a model system and also for water and acetonitrile. It is found that, in general, the ultrafast modes can greatly enhance the rate of electron transfer, even by more than an order of magnitude, over the rate obtained by using only the slow overdamped modes usually considered. For water, this acceleration of the rate can be attributed to the high-frequency intermolecular vibrational and librational modes. For a weakly adiabatic reaction, the rate is virtually indistinguishable from the rate predicted by the Marcus transition state theory. Another important result is that even in this case of ultrafast underdamped solvation, energy diffusion appears to be efficient so that electron transfer reaction in water is controlled essentially by the barrier crossing dynamics. This is because the reactant well frequency is-directly proportional to the rate of the initial Gaussian decay of the solvation time correlation function. As a result, the value of the friction at the reactant well frequency rarely falls below the value required for the Kramers turnover except when the polarizability of the water molecules may be neglected. On the other hand, in acetonitrile, the rate of electron transfer reaction is found to be controlled by the energy diffusion dynamics, although a significant contribution to the rate comes also from the barrier crossing rate. Therefore, the present study calls for a need to understand the relaxation of the high-frequency modes in dipolar liquids.

Self-consistent microscopic treatment of the effects of self-motion of the probe on ionic and dipolar solvation dynamics

A simple but self-consistent microscopic theory for the time dependent solvation energy of both ions and dipoles is presented which includes, for the first time, the details of the self-motion of the probe on its own solvation dynamics. The theory leads to several interesting predictions. The most important of them is that, for dipolar solvation, both the rotational and the translational motions of the dipolar solute probe can significantly accelerate the rate of solvation. In addition, the rotational self-motion of the solute can also give rise to an additional mechanism of nonexponentiality in solvation time correlation functions in otherwise slow liquids. A comparison between the present theoretical predictions and the recent experimental studies of Maroncelli et al. on solvation dynamics of aniline in l-propanol seems to indicate that the said experiments have missed the initial solvent response up to about 45 ps. After mapping the experimental results on the redefined time scale, the theoretical results can explain the experimental results for solvation of aniline in 1-propanol very well. For ionic solvation, the translational motion is significant for light solutes only. For example, for Li+ in water, translational motion speeds up the solvation by about 20%. The present theory demonstrates that in dipolar solvation the partial quenching of the self-motion due to the presence of specific solute-solvent interactions (such as H-bonding) may lead to a much slower solvation than that when the self-motion is present. This point has been discussed. In addition, we present the theoretical results for solvation of aniline in propylene carbonate, Here, the solvation is predicted to be complete within 15-20 ps.

Ultrafast solvation dynamics in water: Isotope effects and comparison with experimental results

A detailed theoretical study of solvation dynamics in water is presented. The motivation of the present study comes from the recent experimental observation that the dynamics of solvation of an ion in water is ultrafast and the solvation time correlation function decays with a time constant of about 55 fs. The slower decay in the long time can be described by a sum of two exponentials with time constants equal to 126 and 880 fs. The molecular theory (developed earlier) predicts a time constant equal to 52 fs for the initial Gaussian decay and time constants equal to 134 and 886 fs for the two exponential components at the long time. This nearly perfect agreement is obtained by using the most detailed dynamical information available in the literature. The present study emphasizes the importance of the intermolecular vibrational band originating from the O...O stretching mode of the O�H...O units in the initial dynamics and raises several interesting questions regarding the nature of the decay of this mode. We have also studied the effects of isotope substitution on solvation dynamics. It is predicted that a significant isotope effect may be observed in the long time. The experimental results have also been compared with the prediction of the dynamic mean spherical approximation (DMSA); the agreement is not satisfactory at the long time. It is further found that the molecular theory and the DMSA lead to virtually identical results if the translational modes of the solvent molecules are neglected in the former. DMSA has also been used to investigate the dynamics of solvation of a dipolar solute in water. It is found that the dynamics of dipolar solvation exhibit features rather different from those of ion solvation. © 1995 American Institute of Physics.

Substituent effect on fluorescence signaling of the cell permeable HSO4- receptors through single point to ratiometric response in green solvent

Two new 2-(2-aminophenyl)benzimidazole-based HSO4- ion selective receptors, 6-(4-nitrophenyl)-5,6-dihydrobenzo4,5]imidazo1,2-c]quinazoline (L1H) and 6-(4-methoxyphenyl)-5,6-dihydrobenzo4,5]imidazo1,2-c] quinazoline (L2H), and their 1 : 1 molecular complexes with HSO4- were prepared in a facile synthetic method and characterized by physicochemical and spectroscopic techniques along with the detailed structural analysis of L1H by single crystal X-ray crystallography. Both receptors (L1H and L2H) behave as highly selective chemosensor for HSO4- ions at biological pH in ethanol-water HEPES buffer (1/5) (v/v) medium over other anions such as F-, Cl-, Br-, I-, AcO-, H2PO4-, N-3(-) and ClO4-. Theoretical and experimental studies showed that the emission efficiency of the receptors (L1H and L2H) was tuned successfully through single point to ratiometric detection by employing the substituent effects. Using 3 sigma method the LOD for HSO4- ions were found to be 18.08 nM and 14.11 nM for L1H and L2H, respectively, within a very short responsive time (15-20 s) in 100 mM HEPES buffer (ethanol-water: 1/5, v/v). Comparison of the utility of the probes (L1H and L2H) as biomarkers for the detection of intracellular HSO4- ions concentrations under a fluorescence microscope has also been included and both probes showed no cytotoxic effect.

Periodically Grafted Amphiphilic Copolymers: Effects of Steric Crowding and Reversal of Amphiphilicity

Two series of periodically clickable polyesters were prepared; one of them carries alkylene segments along its backbone, whereas the other carries poly(ethylene glycol) (PEG) segments. These polyesters were clicked with either MPEG-350 azide or docosyl (C22) azide to yield periodically grafted amphiphilic copolymers (PGACs) carrying either flexible hydrophilic or crystallizable hydrophobic backbone segments. The immiscibility between hydrocarbon and PEG segments causes both of these systems to fold in either a zigzag or hairpin-like conformation; the hairpin-like conformation appears to be preferred when flexible PEG segments are present in the backbone. The folded chains further reorganize in the solid state to develop a lamellar morphology that permits the collocation of the PEG and hydrocarbon (HC) segments within alternate domains; evidence for the self-segregation was gained from DSC, SAXS, and AFM studies. SAXS studies revealed the formation of an extended lamellar structure, whereas AFM images showed uniform layered morphology with layer heights that matched reasonably well with the interlamellar spacing obtained from the SAXS study. Labeling One representative PGAC, carrying crystallizable long alkylene segments in the backbone and pendant PEG-350 side chains, with a small mole fraction of pyrene fluorophore permitted the examination of the conformational transition that occurs upon going from a good to a poor solvent; this single-chain folded conformation, we postulate, is the intermediate that organizes into the lamellar morphology.

Nanoparticle induced miscibility in LCST polymer blends: critically assessing the enthalpic and entropic effects

The use of copolymer and polymer blends widened the possibility of creating materials with multilayered architectures. Hierarchical polymer systems with a wide array of micro and nanostructures are generated by thermally induced phase separation (TIPS) in partially miscible polymer blends. Various parameters like the interaction between the polymers, concentration, solvent/non-solvent ratio, and quenching temperature have to be optimized to obtain these micro/nanophase structures. Alternatively, the addition of nanoparticles is another strategy to design materials with desired hetero-phase structures. The dynamics of the polymer nanocomposite depends on the statistical ordering of polymers around the nanoparticle, which is dependent on the shape of the nanoparticle. The entropic loss due to deformation of polymer chains, like the repulsive interactions due to coiling and the attractive interactions in the case of swelling has been highlighted in this perspective article. The dissipative particle dynamics has been discussed and is correlated with the molecular dynamics simulation in the case of polymer blends. The Cahn Hillard Cook model on variedly shaped immobile fillers has shown difference in the propagation of the composition wave. The nanoparticle shape has a contributing effect on the polymer particle interaction, which can change the miscibility window in the case of these phase separating polymer blends. Quantitative information on the effect of spherical particles on the demixing temperature is well established and further modified to explain the percolation of rod shaped particles in the polymer blends. These models correlate well with the experimental observations in context to the dynamics induced by the nanoparticle in the demixing behavior of the polymer blend. The miscibility of the LCST polymer blend depends on the enthalpic factors like the specific interaction between the components, and the solubility product and the entropic losses occurring due to the formation of any favorable interactions. Hence, it is essential to assess the entropic and enthalpic interactions induced by the nanoparticles independently. The addition of nanoparticles creates heterogeneity in the polymer phase it is localized. This can be observed as an alteration in the relaxation behavior of the polymer. This changes the demixing behavior and the interaction parameter between the polymers. The compositional changes induced due to the incorporation of nanoparticles are also attributed as a reason for the altered demixing temperature. The particle shape anisotropy causes a direction dependent depletion, which changes the phase behavior of the blend. The polymer-grafted nanoparticles with varying grafting density show tremendous variation in the miscibility of the blend. The stretching of the polymer chains grafted on the nanoparticles causes an entropy penalty in the polymer blend. A comparative study on the different shaped particles is not available up to date for understanding these aspects. Hence, we have juxtaposed the various computational studies on nanoparticle dynamics, the shape effect of NPs on homopolymers and also the cases of various polymer blends without nanoparticles to sketch a complete picture on the effect of various particles on the miscibility of LCST blends.

Non-Gravitational Effects with Density-Matching in Evaluating the Influence of Sedimentation on Colloidal Coagulation

The method of density matching between the solid and liquid phases is often adopted to effectively eliminate the effect of sedimentation of suspensions in studies on dynamic behaviour of a colloidal system. However, the associated changes in the solvent composition may bring side effects to the properties investigated and therefore might lead to a faulty conclusion if the relevant correction is not made. To illustrate the importance of this side effect, we present an example of the sedimentation influence on the coagulation rate of suspensions of 2 μm (diameter) polystyrene. The liquid mixtures, in the proper proportions of water (H2O), deuterium oxide (D2O) and methanol (MeOH) as the liquid phase, density-matched and unmatched experiments are performed. Besides the influence of viscosity, the presence of methanol in solvent media, used to enhance the sedimentation effect, causes significant changes (reduction) in rapid coagulation rates compared to that in pure water. Without the relevant corrections for those non-gravitational factors it seems that gravitational sedimentation would retard the coagulation. The magnitude of the contribution from the non-gravitational factor is quantitatively determined, making the relevant correction possible. After necessary corrections for all factors, our experiments show that the influence of the sedimentation on coagulation rates at the initial stage of the coagulation is not observable.

Effects of particle size ratio on single particle motion in colloidal dispersions

The motion of a single Brownian particle of arbitrary size through a dilute colloidal dispersion of neutrally buoyant bath spheres of another characteristic size in a Newtonian solvent is examined in two contexts. First, the particle in question, the probe particle, is subject to a constant applied external force drawing it through the suspension as a simple model for active and nonlinear microrheology. The strength of the applied external force, normalized by the restoring forces of Brownian motion, is the Péclet number, Pe. This dimensionless quantity describes how strongly the probe is upsetting the equilibrium distribution of the bath particles. The mean motion and fluctuations in the probe position are related to interpreted quantities of an effective viscosity of the suspension. These interpreted quantities are calculated to first order in the volume fraction of bath particles and are intimately tied to the spatial distribution, or microstructure, of bath particles relative to the probe. For weak Pe, the disturbance to the equilibrium microstructure is dipolar in nature, with accumulation and depletion regions on the front and rear faces of the probe, respectively. With increasing applied force, the accumulation region compresses to form a thin boundary layer whose thickness scales with the inverse of Pe. The depletion region lengthens to form a trailing wake. The magnitude of the microstructural disturbance is found to grow with increasing bath particle size -- small bath particles in the solvent resemble a continuum with effective microviscosity given by Einstein's viscosity correction for a dilute dispersion of spheres. Large bath particles readily advect toward the minimum approach distance possible between the probe and bath particle, and the probe and bath particle pair rotating as a doublet is the primary mechanism by which the probe particle is able to move past; this is a process that slows the motion of the probe by a factor of the size ratio. The intrinsic microviscosity is found to force thin at low Péclet number due to decreasing contributions from Brownian motion, and force thicken at high Péclet number due to the increasing influence of the configuration-averaged reduction in the probe's hydrodynamic self mobility. Nonmonotonicity at finite sizes is evident in the limiting high-Pe intrinsic microviscosity plateau as a function of bath-to-probe particle size ratio. The intrinsic microviscosity is found to grow with the size ratio for very small probes even at large-but-finite Péclet numbers. However, even a small repulsive interparticle potential, that excludes lubrication interactions, can reduce this intrinsic microviscosity back to an order one quantity. The results of this active microrheology study are compared to previous theoretical studies of falling-ball and towed-ball rheometry and sedimentation and diffusion in polydisperse suspensions, and the singular limit of full hydrodynamic interactions is noted.

Second, the probe particle in question is no longer subject to a constant applied external force. Rather, the particle is considered to be a catalytically-active motor, consuming the bath reactant particles on its reactive face while passively colliding with reactant particles on its inert face. By creating an asymmetric distribution of reactant about its surface, the motor is able to diffusiophoretically propel itself with some mean velocity. The effects of finite size of the solute are examined on the leading order diffusive microstructure of reactant about the motor. Brownian and interparticle contributions to the motor velocity are computed for several interparticle interaction potential lengths and finite reactant-to-motor particle size ratios, with the dimensionless motor velocity increasing with decreasing motor size. A discussion on Brownian rotation frames the context in which these results could be applicable, and future directions are proposed which properly incorporate reactant advection at high motor velocities.