# Modeling the Release Kinetics of Poorly Water-Soluble Drug Molecules from Liposomal Nanocarriers

Authors: Stephan Loew; Alfred Fahr; Sylvio May
Journal: Journal of Drug Delivery
 (2011)

---

## Abstract

Liposomes are frequently used as pharmaceutical nanocarriers to deliver poorly water-soluble drugs such as temoporfin, cyclosporine A, amphotericin B, and paclitaxel to their target site. Optimal drug delivery depends on understanding the release kinetics of the drug molecules from the host liposomes during the journey to the target site and at the target site. Transfer of drugs in model systems consisting of donor liposomes and acceptor liposomes is known from experimental work to typically exhibit a first-order kinetics with a simple exponential behavior. In some cases, a fast component in the initial transfer is present, in other cases the transfer is sigmoidal. We present and analyze a theoretical model for the transfer that accounts for two physical mechanisms, collisions between liposomes and diffusion of the drug molecules through the aqueous phase. Starting with the detailed distribution of drug molecules among the individual liposomes, we specify the conditions that lead to an apparent first-order kinetic behavior. We also discuss possible implications on the transfer kinetics of (1) high drug loading of donor liposomes, (2) attractive interactions between drug molecules within the liposomes, and (3) slow transfer of drugs between the inner and outer leaflets of the liposomes.


---

## Body

## 1. Introduction

Poor solubility in water is a well-recognized obstacle for efficient oral or parenteral drug administration [1, 2]. Liposomes are among the most widely used type of pharmaceutical nanocarriers for small and poorly water-soluble drug molecules [3]. These drugs preferentially partition into the hydrophobic compartment that is formed by the hydrocarbon tails of the liposomal lipids. Liposomes have been used in their first generation (conventional liposomes) predominantly as long-circulating transport vehicles [4, 5], followed by a second generation that improved the circulation time further by decorating the surface with PEG-chains (stealth liposomes [6]). Third-generation liposomes are now being engineered to contain targeting ligands [7] and to carry out stimuli-sensitive triggering of the drug release [8].An important property of liposome-based drug delivery is the release kinetics of the drug from the host, which has been investigated for a number of model systems [9–12]. Experimental investigations of the transfer of temoporfin between two different types of liposomes (i.e., from donor liposomes to acceptor liposomes) have recently been carried out using a mini ion exchange column technique [13]. The column separates donor from acceptor liposomes and thus allows to monitor the time dependence of the drug transfer. It is observed that, typically, the transfer follows an apparent first-order behavior, characterized by a single exponential function. This is remarkable given the complexity of the system, with the drug molecules being able to migrate from the donor to the acceptor liposomes via different physical mechanisms. In fact, there are two mechanisms that, in general, act simultaneously. The first mechanism is the transfer of drugs upon collisions between two liposomes. In this case, the drug molecules directly migrate from one liposome to another with minimal exposure to the aqueous phase. The second mechanism refers to the transfer of drugs via diffusion through the aqueous phase. We note that the collision mechanism has been invoked, for example, to explain the transfer of lipids [14] and cholesterol [15] between vesicles, and the transfer of fatty acids between vesicles and fatty acid binding proteins [16]. Also the diffusion mechanism was found to be consistent with the transport of lipids [17]. In some cases, both mechanisms were suggested to contribute to the transport of lipids between vesicles [18] and to the transport of lipophilic drugs from oil-in-water emulsions to cells [19] and from plasma proteins to lipid vesicles [20]. In our preceding experimental work, where we have investigated the kinetics of temoporfin transport from donor to acceptor liposomes [13], we found that above a certain concentration (corresponding to a liposome-to-liposome distance of about 200 nm for our specific system) the transfer was dominated by collisions; for smaller concentrations transport through diffusion was prevalent.The objective of the present work is to introduce and discuss a detailed kinetic model for the release properties of poorly water-soluble drug molecules from liposomal nanocarriers. Despite a large number of experimental studies about the kinetics of lipid and drug transfer between liposomes and other nanocarriers, there is little theoretical work available that addresses the nature of the transfer kinetics. Our theoretical formalism is based on a detailed distribution function of drug molecules among the individual liposomes. Kinetic rate equations for that distribution function account for two transport mechanisms: collisions between liposomes and drug diffusion through the aqueous phase. We specify a set of conditions at which our microscopic model produces an apparent first-order kinetics with simple exponential behavior, as used in previous work [14, 19]. We point out that our kinetic model can be applied to any kind of small mobile pharmaceutical nanocarrier, including liposomes, micelles [21], colloids [22], and nanoparticles [23].In the second part of our work, we discuss conditions that lead to deviations from simple exponential behavior. First, for the diffusion mechanism, high drug loading tends to increase the transfer rate. The kinetics remains exponential only if donor and acceptor liposomes are chemically similar. Second, the presence of attractive interactions between drug molecules within the liposomes (which can result in the formation of aggregates [24]) is expected to slow down the transfer kinetics. We note that not much molecular detail is presently known about how poorly water-soluble drug molecules inside a lipid bilayer interact. However, modeling studies of rigid membrane-embedded inclusions such as transmembrane proteins or peptides suggest a general tendency of the host membrane to mediate attractive interactions between inclusions that may lead to the formation of aggregates [25]. These attractive interactions may be driven by elastic deformations of the host membrane [26], by depletion of the flexible lipid chains from the region in between rigid inclusions [27], and by fluctuations via the Casimir effect [28]. Our analysis for the collision mechanism suggests that aggregate formation can give rise to sigmoidal behavior of the transfer kinetics. Third, drug molecules (even if they are poorly water soluble) do not necessarily reside predominantly in the innermost region of the membrane's hydrocarbon region. For example, some aromatic compounds such as indole are well known for their preference of the membrane's interfacial region between the headgroups and the hydrocarbon chains [29, 30]. Other aromatic compounds such as benzene are distributed throughout the hydrocarbon chain region without a bias toward the polar/apolar interface [30]. Among the reasons for the preferential partitioning of indole are electrostatic interactions, hydrogen bond formation, and the steric shape of the molecule. For lipid monolayers, there is evidence that drug partitioning also depends on the lateral pressure [31]. Generally, whenever a drug molecule interacts more favorably with the interfacial or headgroup region than with the hydrocarbon tail region, the corresponding partitioning preference can be lumped into at least two energetically preferred states that correspond to the inner and outer leaflet of the membrane. Transfer between the two states (i.e., flip-flop) then introduces an additional characteristic time [32]. We note that two- or multiple-state modeling has been invoked previously to model the partitioning of amino-acid analogues in membranes [33] and the permeation of drug molecules through membranes [34, 35].


## 2. Model

Consider an aqueous solution (of fixed volumeV) that contains a number of Nd donor and Na acceptor liposomes. Donor and acceptor liposomes may be either two chemically different types of liposomes (i.e., composed of different lipids) or equivalent liposomes (i.e., containing the same lipid composition). In the latter case, the distinction between donor and acceptor liposomes refers only to their initial loading; at the end of the transport process (i.e., at thermal equilibrium), both types would be indistinguishable.The donor liposomes initially carry a total number ofM poorly water-soluble drug molecules. Transfer of drug molecules from donor to acceptor liposomes will take place with time until, eventually, an equilibrium partitioning is reached. We describe the time dependence of this transfer by the number of drug molecules carried by the donor liposomes, Md(t), and by the acceptor liposomes, Ma(t). It is then initially Md(t=0)=M and Ma(t=0)=0, as well as in equilibrium Md(t→∞)=Mdeq and Ma(t→∞)=Maeq, where Mdeq and Maeq denote the equilibrium number of drug molecules carried by donor and acceptor liposomes, respectively. We point out that, although we refer to the drug carriers as liposomes, our model is more general. That is, it can be applied to different types of mobile carriers such as micelles, colloids, nanoparticles, or polymeric aggregates, given the carrier possesses some capacity to host poorly water-soluble drug molecules.In the following, we suggest a model for the time dependence of the transfer process (i.e., forMd(t) and Ma(t)) that leads to a first-order kinetics, characterized by a simple exponential function. We consider a “single-state model” where there is only a single energetic state available for each drug molecule in a given liposome. The single-state model excludes the presence of intraliposomal kinetics (the extension to a two-state model will be discussed below). We account for two different transport mechanisms: (i) transport through collisions between liposomes and (ii) transport via diffusion of drug molecules through the aqueous phase. Both mechanisms are schematically illustrated in Figure 1.Figure 1
Transfer of a drug molecule (black bullets) from donor liposome (dark-shaded) to acceptor liposome (light-shaded) upon the collision of the two liposomes or upon diffusion of the drug molecule through the aqueous phase. The displayed scheme refers to the special situation of initially empty acceptor liposomes but analogous schemes apply to any other initial situation.Our transport model of drugs from donor to acceptor liposomes employs the framework of chemical reaction kinetics. We note that due to the generally slow release kinetics of poorly water-soluble drugs, we can treat the aqueous solution asspatially uniform at all times. Hence, no combined diffusion-reaction kinetics [36] needs to be included in our model.


### 2.1. Transfer through Collisions Only

Our model for the collision-mediated drug transfer between liposomes starts with the detailed distribution of drug molecules among all liposomes. We introduce the numberdj of donor liposomes that carry j drug molecules. An analogous definition is used for the number aj of acceptor liposomes that carry j drug molecules. The index j is confined to the region 0≤j≤m where m is the maximal number of drug molecules that a liposome can carry. The time-dependent distribution functions dj=dj(t) and aj=aj(t) represent a full microscopic knowledge of the kinetics of drug transfer. The total numbers of donor liposomes Nd, acceptor liposomes Na, drug molecules residing in donor liposomes Md, and drug molecules residing in acceptor liposomes Ma, can be calculated on the basis of the distribution functions dj=dj(t) and aj=aj(t) according to(1)Nd=∑j=0mdj,Na=∑j=0maj,Md=∑j=0mjdj,Ma=∑j=0mjaj.
Mathematically, Nd and Na are the zeroth-moments of the distributions functions dj=dj(t) and aj=aj(t) whereas Md and Ma appear as the corresponding first moments. We assume that Nd and Na are constant (i.e., independent of time), and so then is the total number of liposomes N=Nd+Na. This is appropriate if fusion and fission between liposomes can be ignored. Due to our focus on poorly water-soluble drug molecules, it is also justified to assume that the total number of drug molecules carried by all liposomes, M=Md+Ma, is constant. That is, we neglect the small fraction of drug molecules that reside in the aqueous phase without being bound to a liposome. Figure 2  schematically illustrates a specific exemplification of the system.Figure 2
Exemplification of our system:Nd=6 donor liposomes (dark shaded) and Na=9 acceptor liposomes (light shaded) reside in an aqueous space of volume V; each liposome can carry at most m=3 drug molecules (black bullets). For the displayed hypothetical snapshot (taken at a certain time t), the distribution function of drug molecules among the donor liposomes is d0=1, d1=2, d2=2, d3=1, leading to a total of Md=∑j=0mjdj=9 drug molecules residing in donor liposomes. Analogously for the acceptor liposomes, the distribution function is a0=5, a1=3, a2=1, a3=0, implying Ma=∑j=0mjaj=5.Collisions require two liposomes to come to close proximity. The magnitude of drug transport between, say, donor liposomesdi and dj is thus ~di×dj/V where V is the volume of the aqueous solution. The underlying transfer process is thus second order. If a single drug molecule is transferred from a donor that carries initially i drug molecules to a donor that carries initially j drug molecules, the distribution function changes according to di→di-1, di-1→di-1+1, dj→dj-1, and dj+1→dj+1+1. Hence, the numbers di and dj decrease whereas di-1 and dj+1 increase. Figure 3  shows an illustration of this scheme for i=5 and j=1. The transfer rate between the populations di and dj will also depend on the corresponding numbers of drug molecules i and j. We assume the drug molecules within each liposome form an ideal mixture so that the transfer rate is directly proportional to |i-j|. In writing a rate equation for donor population dj, we have to account for all possible ways of collisions between donor liposomes of index j with other liposomes (donors and acceptors) of index i. These considerations lead us to(2)VKcollḋj=∑i=0jdi[dj+1g(j+1,i)-djg(j,i)]+∑i=jmdi[dj-1g(i,j-1)-djg(i,j)]+∑i=kjai-k[dj+1g(j+1,i)-djg(j,i)]+∑i=jm+kai-k[dj-1g(i,j-1)-djg(i,j)],
where we have defined the function(3)g(i,j)=i-j.
In (2), Kcoll is the unit rate of drug transfer through collisions between two chemically equivalent liposomes, and ẋ=dx/dt denotes the time derivative of a physical quantity x(t). The first two lines in (2) account for collisions of donor liposomes with other donor liposomes. The last two lines in (2) account for collisions of donor liposomes with acceptor liposomes.Figure 3
Transfer of a drug molecule (black bullets) upon the collision of two liposomes (here assumed to be two donor liposomes). The drug distribution function changes from initiallyd1=1, d2=0, d3=0, d4=0, d5=1 to d1=0, d2=1, d3=0, d4=1, d5=0. This represents an example (for i=5 and j=1) of the general scheme di→di-1, di-1→di-1+1, dj→dj-1, and dj+1→dj+1+1.Note that  (2) allows for a difference in the chemical nature of donor and acceptor liposomes. This chemical mismatch is accounted for by the integer k in the last two lines of (2), which expresses the difference in the number of drug molecules between a donor and acceptor liposomes in thermal equilibrium, (i.e., for k=0 each donor and acceptor liposome will contain the same number of drug molecules in thermal equilibrium). We do not attempt to calculate k from a microscopic model; yet below we show how k is related to the change in standard Gibbs free energy for the process of transferring drug molecules from donor to acceptor liposomes.Due to symmetry, we obtainȧj from ḋj by replacing di→ai, ai→di, and k→-k,(4)VKcollȧj=∑i=0jai[aj+1g(j+1,i)-ajg(j,i)]+∑i=jmai[aj-1g(i,j-1)-ajg(i,j)]+∑i=-kjdi+k[aj+1g(j+1,i)-ajg(j,i)]+∑i=jm-kdi+k[aj-1g(i,j-1)-ajg(i,j)].
Equations (2) and (4) constitute a microscopic model for the kinetic behavior of drug transport from donor to acceptor liposomes through the collision mechanism; it can be verified that(5)∑j=0mḋj=∑j=0mȧj=∑j=0mj(ȧj+ḋj)=0,
implying Ṅd=Ṅa=Ṁ=0 and thus ensuring conservation of the number of donor and acceptor liposomes (Nd and Na) as well as of the total number of drug molecules (M=Md+Ma). To characterize the total numbers Md and Ma of drug molecules that reside in donor and acceptor liposomes, respectively, we carry out the summations ∑j=0mjḋj and ∑j=0mjȧj using  (2) and (4). The result are the two first-order differential equations(6)Ṁd=KN(MaNd-MdNa+kNaNd),Ṁa=KN(MdNa-MaNd-kNaNd),
where we have introduced the definition of the apparent rate constant(7)K=KcollNV.
Initially, all drug molecules are incorporated in the donor liposomes, implying Md(t=0)=M and Ma(t=0)=0. The solution of  (6) is then(8)Ma(t)=M-Md(t)=(1-e-Kt)NaN(M-kNd).
Hence, K indeed appears as the inverse characteristic time for the transfer process. In contrast to previous models [14], K depends only on the total concentration of liposomes N/V but not on the concentrations of donor or acceptor liposomes individually. We also mention that  (6) and the solution in (8) are valid for any number of donor and acceptor liposomes (i.e, any choice of Nd and Na). This includes but is not restricted to sink conditions (where Na≫Nd).Thermodynamic equilibrium corresponds to the long-time limit,t→∞, at which we have Md=Mdeq and Ma=Maeq with(9)MdeqM=NdN(1+kNaM),MaeqM=NaN(1-kNdM).
From  (9), we obtain the difference between the numbers of drug molecules per donor and acceptor liposome, (Mdeq/Nd)-(Maeq/Na)=k. This agrees with our interpretation of k in  (2) and (4). We note that for chemically identical donor and acceptor liposomes, it is k=0 and all liposomes carry the same number of drug molecules in equilibrium, implying Mdeq/Nd=Maeq/Na=M/N. The largest possible value of k is k=M/Nd for which we obtain Maeq=0 and Mdeq=M. The smallest possible value of k is k=-M/Na implying Maeq=M and Mdeq=0. Hence, -M/Na≤k≤M/Nd.The solution in (8) corresponds to a simple exponential decay of the number of drug molecules in the donor liposomes. This suggests that we can express the transfer kinetics of drug molecules from donor (D) to acceptor (A) liposomes as the chemical reaction scheme(10)D⇌K2K1A,
with rate constants K1 and K2. The corresponding kinetic behavior is then governed by the equations Ṁd=-K1Md+K2Ma and Ṁa=K1Md-K2Ma where Md=Md(t) and Ma=Ma(t) are the numbers of drug molecules carried by donor and acceptor liposomes, respectively. With Md(t=0)=M and Ma(t=0)=0 we obtain(11)Ma(t)=M-Md(t)=(1-e-(K1+K2)t)(K1K1+K2)NaNM,
which has indeed the same structure as (8). Comparison of (8) with (11) reveals K1=(1-kNd/M)KNa/N and K2=(1+kNa/M)KNd/N. The equilibrium constant Keq=K1/K2 of the reaction in (10) is thus(12)Keq=NaNdM-kNdM+kNa.
Comparing this with Keq=exp(-Δg0/kBT) (where kB is Boltzmann's constant and T is the absolute temperature) allows us to compute the change in standard Gibbs free energy(13)Δg0=kBTln⁡(M/Na+kM/Nd-k),
for the transfer of a single drug molecule from a donor to an acceptor liposome. The enthalpic and entropic contributions to Δg0 will be influenced by k, which is, generally, temperature dependent(k=k(t)). Let us briefly discuss two cases. First, if donor and acceptor liposomes are chemically identical, then k=0 and Δg0=kBTln(Nd/Na) has only an entropic contribution. Specifically, for Nd>Na, we find Δg0>0 because a given drug molecule has more donor liposomes to reside in than acceptor liposomes. Second, the limiting cases for k, namely, k=-M/Na and k=M/Nd, yield Δg0→-∞ (thus, with all drugs migrating to the acceptor liposomes) and Δg0→∞ (thus with all drugs remaining in the donor liposomes), respectively.We point out that our model predicts a simple exponential time behavior despite the presence of drug transfer through a second-order two-body collision process (i.e., collisions between two liposomes). Chemical reactions that deplete the reactants through binary collisions generally display a long time-tailc(t)~1/t in their concentration dependence. For example, the kinetic behavior of the dimerization reaction 2 monomer → dimer follows the equation ċ=k~c2   where c(t) is the concentration of the reactant (i.e., the monomers) and k~ the rate constant. With an initial concentration c(t=0)=c0 the time behavior becomes c(t)=c0/(1+k~t), implying c(t)~1/t for long times. For our system, however, the numbers of donor and acceptor liposomes remain unchanged. Thus, collisions do not deplete the reactants, and the concentration dependencies of Md(t) and Ma(t) become exponential in time.


### 2.2. Transfer through Diffusion Only

Diffusion allows for transfer of drug molecules directly through the aqueous phase, without the need of collisions between liposomes. Denoting the additional state in the aqueous phase by W (in addition to donor (D) and acceptor (A)) the corresponding transport scheme (again, as in (10), formally expressed as a chemical reaction) can be written as [14, 37](14)D⇌KduptKdrel⁡W⇌Karel⁡KauptA,
with rate constants Kdrel, Karel, Kdupt, and Kaupt for the drug release (“rel”) and uptake (“upt”) in donor (“d”) and acceptor (“a”) liposomes. To formulate the rate equations, it is useful to first consider the drug distribution function dj(t). We assume the probability of a drug molecule to leave donor liposomes of index j to be proportional to the total number jdj of drug molecules in that liposome population. Similarly, the probability of a drug molecule to enter donor liposomes of index j is assumed to be proportional to the total number (m-j)dj of empty binding sites in that liposome population. Because the uptake is based on collisions of liposomes with drug molecules in the aqueous solution, the rate should also be proportional to the drug concentration Mw/V in the aqueous phase (here, Mw is the total number of drug molecules residing in the aqueous phase). This leads to the following rate equations(15)ḋj=Kdrel⁡[(j+1)dj+1-jdj]+KduptMwV[(m-(j-1))dj-1-(m-j)dj],
for 0≤j≤m (with dj=aj=0 for j<0 or j>m). A similar equation can be written for the acceptor liposomes. Based on (15), it can be verified that ∑j=0mḋj=0, thus ensuring conservation of Nd (and similarly for Na). Carrying out the summation Ṁd=∑j=0mjḋj using (15) leads to(16)Ṁd=-Kdrel⁡Md+KduptMwV(mNd-Md).
This equation simply expresses the proportionality of the release to the total number of bound drug molecules and the proportionality of the uptake to the total number of free binding sites. Consistent with (16) we complete the set of rate equations corresponding to the scheme in (14) (17)Ṁw=Kdrel⁡Md-KduptMwV(mNd-Md)+Karel⁡Ma-KauptMwV(mNa-Ma),Ṁa=-Karel⁡Ma+KauptMwV(mNa-Ma).
To obtain first-order behavior, we make three assumptions. The first is a steady-state approximation for the number of drug molecules in the aqueous phase, Ṁw=0. The solubility limit of poorly water-soluble drugs is small so that, effectively, any release of drugs from one liposome is accompanied by an immediate uptake by another (or the same [38]) liposome. The second assumption is weak drug loading of all liposomes; this amounts to Md≪mNd, Ma≪mNa, and M≪mN. We finally assume the same rate for the uptake of drug molecules from the aqueous phase into donor and acceptor liposomes, implying Kdupt=Kaupt. This is strictly valid only for chemically equivalent donor and acceptor liposomes but should generally be a reasonable approximation. That is, we expect the energy barrier for entering a liposome from the aqueous phase to be small (as compared to the energy barrier for the release from a liposome), irrespective of the liposome's chemical structure. Subject to our three assumptions (16) and (17) become equivalent to(18)Ṁd=-Kdrel⁡NaNMd+Karel⁡NdNMa,Ṁa=Kdrel⁡NaNMd-Karel⁡NdNMa.
Equation (18) are now identical to (6) if we identify Kdrel=Kdiff(1-(kNd)/M) and Karel=Kdiff(1+kNa/M) where Kdiff=K appears as the rate constant. Here again, as for (6), the validity of (18) is not subject to a restriction with respect to Nd and Na.


## 2.1. Transfer through Collisions Only

Our model for the collision-mediated drug transfer between liposomes starts with the detailed distribution of drug molecules among all liposomes. We introduce the numberdj of donor liposomes that carry j drug molecules. An analogous definition is used for the number aj of acceptor liposomes that carry j drug molecules. The index j is confined to the region 0≤j≤m where m is the maximal number of drug molecules that a liposome can carry. The time-dependent distribution functions dj=dj(t) and aj=aj(t) represent a full microscopic knowledge of the kinetics of drug transfer. The total numbers of donor liposomes Nd, acceptor liposomes Na, drug molecules residing in donor liposomes Md, and drug molecules residing in acceptor liposomes Ma, can be calculated on the basis of the distribution functions dj=dj(t) and aj=aj(t) according to(1)Nd=∑j=0mdj,Na=∑j=0maj,Md=∑j=0mjdj,Ma=∑j=0mjaj.
Mathematically, Nd and Na are the zeroth-moments of the distributions functions dj=dj(t) and aj=aj(t) whereas Md and Ma appear as the corresponding first moments. We assume that Nd and Na are constant (i.e., independent of time), and so then is the total number of liposomes N=Nd+Na. This is appropriate if fusion and fission between liposomes can be ignored. Due to our focus on poorly water-soluble drug molecules, it is also justified to assume that the total number of drug molecules carried by all liposomes, M=Md+Ma, is constant. That is, we neglect the small fraction of drug molecules that reside in the aqueous phase without being bound to a liposome. Figure 2  schematically illustrates a specific exemplification of the system.Figure 2
Exemplification of our system:Nd=6 donor liposomes (dark shaded) and Na=9 acceptor liposomes (light shaded) reside in an aqueous space of volume V; each liposome can carry at most m=3 drug molecules (black bullets). For the displayed hypothetical snapshot (taken at a certain time t), the distribution function of drug molecules among the donor liposomes is d0=1, d1=2, d2=2, d3=1, leading to a total of Md=∑j=0mjdj=9 drug molecules residing in donor liposomes. Analogously for the acceptor liposomes, the distribution function is a0=5, a1=3, a2=1, a3=0, implying Ma=∑j=0mjaj=5.Collisions require two liposomes to come to close proximity. The magnitude of drug transport between, say, donor liposomesdi and dj is thus ~di×dj/V where V is the volume of the aqueous solution. The underlying transfer process is thus second order. If a single drug molecule is transferred from a donor that carries initially i drug molecules to a donor that carries initially j drug molecules, the distribution function changes according to di→di-1, di-1→di-1+1, dj→dj-1, and dj+1→dj+1+1. Hence, the numbers di and dj decrease whereas di-1 and dj+1 increase. Figure 3  shows an illustration of this scheme for i=5 and j=1. The transfer rate between the populations di and dj will also depend on the corresponding numbers of drug molecules i and j. We assume the drug molecules within each liposome form an ideal mixture so that the transfer rate is directly proportional to |i-j|. In writing a rate equation for donor population dj, we have to account for all possible ways of collisions between donor liposomes of index j with other liposomes (donors and acceptors) of index i. These considerations lead us to(2)VKcollḋj=∑i=0jdi[dj+1g(j+1,i)-djg(j,i)]+∑i=jmdi[dj-1g(i,j-1)-djg(i,j)]+∑i=kjai-k[dj+1g(j+1,i)-djg(j,i)]+∑i=jm+kai-k[dj-1g(i,j-1)-djg(i,j)],
where we have defined the function(3)g(i,j)=i-j.
In (2), Kcoll is the unit rate of drug transfer through collisions between two chemically equivalent liposomes, and ẋ=dx/dt denotes the time derivative of a physical quantity x(t). The first two lines in (2) account for collisions of donor liposomes with other donor liposomes. The last two lines in (2) account for collisions of donor liposomes with acceptor liposomes.Figure 3
Transfer of a drug molecule (black bullets) upon the collision of two liposomes (here assumed to be two donor liposomes). The drug distribution function changes from initiallyd1=1, d2=0, d3=0, d4=0, d5=1 to d1=0, d2=1, d3=0, d4=1, d5=0. This represents an example (for i=5 and j=1) of the general scheme di→di-1, di-1→di-1+1, dj→dj-1, and dj+1→dj+1+1.Note that  (2) allows for a difference in the chemical nature of donor and acceptor liposomes. This chemical mismatch is accounted for by the integer k in the last two lines of (2), which expresses the difference in the number of drug molecules between a donor and acceptor liposomes in thermal equilibrium, (i.e., for k=0 each donor and acceptor liposome will contain the same number of drug molecules in thermal equilibrium). We do not attempt to calculate k from a microscopic model; yet below we show how k is related to the change in standard Gibbs free energy for the process of transferring drug molecules from donor to acceptor liposomes.Due to symmetry, we obtainȧj from ḋj by replacing di→ai, ai→di, and k→-k,(4)VKcollȧj=∑i=0jai[aj+1g(j+1,i)-ajg(j,i)]+∑i=jmai[aj-1g(i,j-1)-ajg(i,j)]+∑i=-kjdi+k[aj+1g(j+1,i)-ajg(j,i)]+∑i=jm-kdi+k[aj-1g(i,j-1)-ajg(i,j)].
Equations (2) and (4) constitute a microscopic model for the kinetic behavior of drug transport from donor to acceptor liposomes through the collision mechanism; it can be verified that(5)∑j=0mḋj=∑j=0mȧj=∑j=0mj(ȧj+ḋj)=0,
implying Ṅd=Ṅa=Ṁ=0 and thus ensuring conservation of the number of donor and acceptor liposomes (Nd and Na) as well as of the total number of drug molecules (M=Md+Ma). To characterize the total numbers Md and Ma of drug molecules that reside in donor and acceptor liposomes, respectively, we carry out the summations ∑j=0mjḋj and ∑j=0mjȧj using  (2) and (4). The result are the two first-order differential equations(6)Ṁd=KN(MaNd-MdNa+kNaNd),Ṁa=KN(MdNa-MaNd-kNaNd),
where we have introduced the definition of the apparent rate constant(7)K=KcollNV.
Initially, all drug molecules are incorporated in the donor liposomes, implying Md(t=0)=M and Ma(t=0)=0. The solution of  (6) is then(8)Ma(t)=M-Md(t)=(1-e-Kt)NaN(M-kNd).
Hence, K indeed appears as the inverse characteristic time for the transfer process. In contrast to previous models [14], K depends only on the total concentration of liposomes N/V but not on the concentrations of donor or acceptor liposomes individually. We also mention that  (6) and the solution in (8) are valid for any number of donor and acceptor liposomes (i.e, any choice of Nd and Na). This includes but is not restricted to sink conditions (where Na≫Nd).Thermodynamic equilibrium corresponds to the long-time limit,t→∞, at which we have Md=Mdeq and Ma=Maeq with(9)MdeqM=NdN(1+kNaM),MaeqM=NaN(1-kNdM).
From  (9), we obtain the difference between the numbers of drug molecules per donor and acceptor liposome, (Mdeq/Nd)-(Maeq/Na)=k. This agrees with our interpretation of k in  (2) and (4). We note that for chemically identical donor and acceptor liposomes, it is k=0 and all liposomes carry the same number of drug molecules in equilibrium, implying Mdeq/Nd=Maeq/Na=M/N. The largest possible value of k is k=M/Nd for which we obtain Maeq=0 and Mdeq=M. The smallest possible value of k is k=-M/Na implying Maeq=M and Mdeq=0. Hence, -M/Na≤k≤M/Nd.The solution in (8) corresponds to a simple exponential decay of the number of drug molecules in the donor liposomes. This suggests that we can express the transfer kinetics of drug molecules from donor (D) to acceptor (A) liposomes as the chemical reaction scheme(10)D⇌K2K1A,
with rate constants K1 and K2. The corresponding kinetic behavior is then governed by the equations Ṁd=-K1Md+K2Ma and Ṁa=K1Md-K2Ma where Md=Md(t) and Ma=Ma(t) are the numbers of drug molecules carried by donor and acceptor liposomes, respectively. With Md(t=0)=M and Ma(t=0)=0 we obtain(11)Ma(t)=M-Md(t)=(1-e-(K1+K2)t)(K1K1+K2)NaNM,
which has indeed the same structure as (8). Comparison of (8) with (11) reveals K1=(1-kNd/M)KNa/N and K2=(1+kNa/M)KNd/N. The equilibrium constant Keq=K1/K2 of the reaction in (10) is thus(12)Keq=NaNdM-kNdM+kNa.
Comparing this with Keq=exp(-Δg0/kBT) (where kB is Boltzmann's constant and T is the absolute temperature) allows us to compute the change in standard Gibbs free energy(13)Δg0=kBTln⁡(M/Na+kM/Nd-k),
for the transfer of a single drug molecule from a donor to an acceptor liposome. The enthalpic and entropic contributions to Δg0 will be influenced by k, which is, generally, temperature dependent(k=k(t)). Let us briefly discuss two cases. First, if donor and acceptor liposomes are chemically identical, then k=0 and Δg0=kBTln(Nd/Na) has only an entropic contribution. Specifically, for Nd>Na, we find Δg0>0 because a given drug molecule has more donor liposomes to reside in than acceptor liposomes. Second, the limiting cases for k, namely, k=-M/Na and k=M/Nd, yield Δg0→-∞ (thus, with all drugs migrating to the acceptor liposomes) and Δg0→∞ (thus with all drugs remaining in the donor liposomes), respectively.We point out that our model predicts a simple exponential time behavior despite the presence of drug transfer through a second-order two-body collision process (i.e., collisions between two liposomes). Chemical reactions that deplete the reactants through binary collisions generally display a long time-tailc(t)~1/t in their concentration dependence. For example, the kinetic behavior of the dimerization reaction 2 monomer → dimer follows the equation ċ=k~c2   where c(t) is the concentration of the reactant (i.e., the monomers) and k~ the rate constant. With an initial concentration c(t=0)=c0 the time behavior becomes c(t)=c0/(1+k~t), implying c(t)~1/t for long times. For our system, however, the numbers of donor and acceptor liposomes remain unchanged. Thus, collisions do not deplete the reactants, and the concentration dependencies of Md(t) and Ma(t) become exponential in time.


## 2.2. Transfer through Diffusion Only

Diffusion allows for transfer of drug molecules directly through the aqueous phase, without the need of collisions between liposomes. Denoting the additional state in the aqueous phase by W (in addition to donor (D) and acceptor (A)) the corresponding transport scheme (again, as in (10), formally expressed as a chemical reaction) can be written as [14, 37](14)D⇌KduptKdrel⁡W⇌Karel⁡KauptA,
with rate constants Kdrel, Karel, Kdupt, and Kaupt for the drug release (“rel”) and uptake (“upt”) in donor (“d”) and acceptor (“a”) liposomes. To formulate the rate equations, it is useful to first consider the drug distribution function dj(t). We assume the probability of a drug molecule to leave donor liposomes of index j to be proportional to the total number jdj of drug molecules in that liposome population. Similarly, the probability of a drug molecule to enter donor liposomes of index j is assumed to be proportional to the total number (m-j)dj of empty binding sites in that liposome population. Because the uptake is based on collisions of liposomes with drug molecules in the aqueous solution, the rate should also be proportional to the drug concentration Mw/V in the aqueous phase (here, Mw is the total number of drug molecules residing in the aqueous phase). This leads to the following rate equations(15)ḋj=Kdrel⁡[(j+1)dj+1-jdj]+KduptMwV[(m-(j-1))dj-1-(m-j)dj],
for 0≤j≤m (with dj=aj=0 for j<0 or j>m). A similar equation can be written for the acceptor liposomes. Based on (15), it can be verified that ∑j=0mḋj=0, thus ensuring conservation of Nd (and similarly for Na). Carrying out the summation Ṁd=∑j=0mjḋj using (15) leads to(16)Ṁd=-Kdrel⁡Md+KduptMwV(mNd-Md).
This equation simply expresses the proportionality of the release to the total number of bound drug molecules and the proportionality of the uptake to the total number of free binding sites. Consistent with (16) we complete the set of rate equations corresponding to the scheme in (14) (17)Ṁw=Kdrel⁡Md-KduptMwV(mNd-Md)+Karel⁡Ma-KauptMwV(mNa-Ma),Ṁa=-Karel⁡Ma+KauptMwV(mNa-Ma).
To obtain first-order behavior, we make three assumptions. The first is a steady-state approximation for the number of drug molecules in the aqueous phase, Ṁw=0. The solubility limit of poorly water-soluble drugs is small so that, effectively, any release of drugs from one liposome is accompanied by an immediate uptake by another (or the same [38]) liposome. The second assumption is weak drug loading of all liposomes; this amounts to Md≪mNd, Ma≪mNa, and M≪mN. We finally assume the same rate for the uptake of drug molecules from the aqueous phase into donor and acceptor liposomes, implying Kdupt=Kaupt. This is strictly valid only for chemically equivalent donor and acceptor liposomes but should generally be a reasonable approximation. That is, we expect the energy barrier for entering a liposome from the aqueous phase to be small (as compared to the energy barrier for the release from a liposome), irrespective of the liposome's chemical structure. Subject to our three assumptions (16) and (17) become equivalent to(18)Ṁd=-Kdrel⁡NaNMd+Karel⁡NdNMa,Ṁa=Kdrel⁡NaNMd-Karel⁡NdNMa.
Equation (18) are now identical to (6) if we identify Kdrel=Kdiff(1-(kNd)/M) and Karel=Kdiff(1+kNa/M) where Kdiff=K appears as the rate constant. Here again, as for (6), the validity of (18) is not subject to a restriction with respect to Nd and Na.


## 3. Discussion

Both transfer mechanisms, through liposome collisions and via diffusion through the aqueous phase, lead to the same first-order kinetic behavior; see (6) and (18). The rate constant of the combined process is(19)K=KcollNV+Kdiff.
Its dependence on the total liposome concentration allows the experimental determination of the transfer mechanism [13]. We note that the first-order behavior predicted by (6) and (18) requires several assumptions to be fulfilled: low liposome loading with drug molecules, rate constants that are strictly proportional to concentrations of drug molecules, and no intraliposomal kinetics with a rate similar to K. In the following, we discuss how the kinetic behavior is predicted to change if any of these assumptions is not fulfilled.


### 3.1. Extension to High Drug Loading

While high drug loading obviously increases the number of available drug molecules (and thus increases the efficiency of liposomal carriers [39]) it also affects the kinetics of the drug release. Our present model predicts such a dependence for the diffusion mechanism whereas the kinetics for the collision mechanism is not affected. Recall that the transition from  (16) and (17) to (18) was based on the approximation of weak drug loading, Md≪mNd, Ma≪mNa, and M≪mN. Without that approximation, we obtain instead of (18) a nonlinear set of differential equations(20)Ṁd=-Kdrel⁡Na/N-(Ma/M)(M/mN)1-M/mNMd+Karel⁡Nd/N-(Md/M)(M/mN)1-M/mNMa,Ṁa=Kdrel⁡Na/N-(Ma/M)(M/mN)1-M/mNMd-Karel⁡Nd/N-(Md/M)(M/mN)1-M/mNMa.
For the special case that donor and acceptor liposomes are chemically similar, Kdrel=Karel=Kdiff, we obtain a simple exponential behavior(21)Ma(t)=M-Md(t)=(1-e-Kdifft/(1-M/mN))NaNM.
Here, high drug loading simply increases the rate constant for the diffusion mechanism by the factor 1/(1-M/(mN)). In the general case Kdrel≠Karel, and no simple exponential decay is predicted for high loading of the liposomes with drug molecules. Figure 4 shows a numerical example, based on (20) with Kdrel/Karel=3 and Nd/N=Na/N=0.5. For M≪mN (weak loading regime; broken lines in Figure 4) we observe the simple exponential behavior according to (18) with equilibrium values Mdeq/M=1/4 and Maeq/M=3/4. For M/(mN)=0.5 the initial loading of the donor liposomes is maximal. This leads to both a faster decay and a shift in the equilibrium distribution, reaching Mdeq/M=(3-1)/2=0.366 and Maeq/M=(3-3)/2=0.634. The reason for the increased rate constant is the reduced ability of highly loaded liposomes to take up drug molecules. Hence, if drug molecules are released from initially highly loaded donor liposomes they will be taken up exclusively by acceptor liposomes. The increase in the transfer rate at high loading also affects the equilibrium values Mdeq/M and Maeq/M. The equilibrium is shifted toward a more uniform distribution of drug molecules between donor and acceptor liposomes (in agreement with Figure 4).Figure 4
Numerical solutions of  (20), derived for M/(Nm)=0 (broken lines) and M/(Nm)=0.5 (solid lines). The remaining parameters are Kdrel/Karel=3, Nd/N=Na/N=0.5. The time t is plotted in units of 1/Karel.


### 3.2. Sigmoidal Behavior

Our model presented so far is unable to predict sigmoidal behavior. That is, no inflection point can be observed inMd(t) and Ma(t). Behind this prediction is our assumption that the transfer rates are strictly proportional to the concentration difference of the drug molecules. For the collision mechanism, this is expressed by our definition of the function g(i,j) in (3). However, if drug molecules within a given liposome interact with each other, the simple relation g(i,j)=i-j will no longer be valid. More specifically, attractive interactions between drug molecules within liposomes will increase the energy barrier to remove a drug molecule. This becomes relevant at high drug loading. Hence, in the presence of attractive interactions, it will be more unlikely that a drug molecule is transferred from a highly loaded donor liposome to an empty acceptor liposome.We discuss the consequences of attractive interactions for the collision mechanisms, which is described by (2) and (4). To account for the decrease in the rate constant at high loading we replace (3) by(22)g(i,j)=(i-j)(1-im)(1-jm).
Clearly, for weak loading (i≪m and j≪m) the original first-order model leading to the exponential behavior in (8) is recovered. For large loading of either donor or acceptor liposomes, the transfer rate becomes small. We note that using (22) does not lead to a set of differential equations in terms of only Md(t) and Ma(t). Here, we do not attempt to provide an analytical solution to the problem. Instead, we illustrate its predictions by numerically solving (2) and (4) with g(i,j) given in (22).Figure5  shows the behavior of Md(t) and Ma(t) as function of tK (with K=KcollN/V), derived for m=100. For simplicity, we have set k=0 which results in an equipartitioning of drug molecules between donor and acceptor liposomes (Md/Nd=Ma/Na=M/N). We start with Nd=Na=100 liposomes. The acceptor liposomes are initially empty whereas each donor liposome contains initially l drug molecules (out of a maximal number m=100). Different curves in Figure 5  correspond to l=2 (a), l=10 (b), l=50 (c), l=90 (d), and l=98 (e). As long as the drug loading is weak (curves (a) and (b)), the solution is simply exponential, characterized by Ma/M=1-Md/M=(1-e-Kt)Na/N (see (8) with k=0). Here, the kinetics is independent of the total number of drug molecules M=lNd (which is why curves (a) and (b) virtually overlap). If the initial loading of the donor liposomes becomes larger (curve (c)) the kinetics slows down. Eventually, once the initial loading approaches its maximal value mNd, the behavior slows down even more and, in addition, becomes sigmoidal. Attractive drug-drug interactions slow down the release from initially highly loaded donor liposomes; at later times (when the donor liposomes are no longer highly loaded), the release becomes faster. This leads to sigmoidal behavior.Figure 5
Fraction of drug molecules contained in donor liposomes (Md(t)/M; upper set of curves) and acceptor liposomes (Ma(t)/M; lower set of curves) as function of the scaled time Kt. The curves represent numerical solutions of  (2) and (4) with (22), derived for k=0 and m=100 with the initial conditions dj(0)=0 for j≠l, dl(0)=100, aj(0)=0 for j>0, a0(0)=100. Different curves correspond to l=2 (a), l=10 (b), l=50 (c), l=90 (d), and l=98 (e).


### 3.3. Extension to a Two-State Model

In the final part of this work, we briefly discuss an extension of our model to account for two distinct states of the drug molecule inside each liposome. A simple rationale for the presence of two distinct states is provided by the bilayer structure of the liposomes. That is, a drug molecule may preferentially be bound to either the inner or outer monolayer, having to flip-flop in order to change the host monolayer. The typical flip-flop time can be large if the drug has some amphiphilicity or surface activity instead of being strongly lipophilic [40]. Drug molecules residing in the inner monolayer cannot be transported directly to another liposome; they first have to migrate to the outer monolayer.We denote byMdI and MdO the number of drug molecules residing in the inner (DI) and outer (DO) leaflets of donor liposomes, respectively. Similarly, MaI and MaO refer to the number of drug molecules residing in the inner (AI) and outer (AO) leaflets of acceptor liposomes. The reaction scheme in (10) can then be generalized to account for the inter leaflet transport in donor and acceptor liposomes(23)DI⇌K2dK1dDO⇌K2K1AO⇌K1aK2aAI.
Here, K1d and K2d are the two rate constants corresponding to the transfer of drugs between the two leaflets of the donor liposomes (and similarly for K1a and K2a referring to the acceptor liposomes). The rate constants K1=(1-kNd/M)KNa/N and K2=(1+kNa/M)KNd/N are identical to those for the single-state model, where K is given in (19). Based on (23), the rate equations can be written as(24)ṀdO=KN(MaONd-MdONa+kNaNd)-K2dMdO+K1dMdI,ṀdI=K2dMdO-K1dMdI,ṀaO=KN(MdONa-MaONd-kNaNd)-K2aMaO+K1aMaI,ṀaI=K2aMaO-K1aMaI.
In the limit of a symmetric lipid bilayer, the two rate constants for flip-flop of a drug molecule from the inner to the outer leaf and from the outer to the inner leaf are identical (we note that the two leaflets of a liposomal bilayer are not strictly equivalent which, in a more refined model, would entail two different rate constants for flip-flop; this dependence on liposome curvature is neglected here). If we assume furthermore that donor and acceptor liposomes are chemically similar, we may write K1d=K2d=K1a=K2a=G as well as k=0. In this case, the rate equations(25)ṀdO=KN(MaONd-MdONa)-G(MdO-MdI),ṀdI=G(MdO-MdI),ṀaO=KN(MdONa-MaONd)-G(MaO-MaI),ṀaI=G(MaO-MaI)
depend on only two parameters, the rate constants K and G. If we assume all drug molecules initially reside in the donor liposomes, the initial conditions are MdO(t=0)=MdI(t=0)=M/2, and MaO(t=0)=MaI(t=0)=0, where M is the total number of drug molecules in the system. The solution of  (25) can be expressed as(26)MdI(t)=M2[NdN+NaNω2e-ω1t-ω1e-ω2tω2-ω1],MdO(t)-MdI(t)=M2KNaNe-ω2t-e-ω1tω2-ω1,MaI(t)=MNa2N[1-ω2e-ω1t-ω1e-ω2tω2-ω1],MaO(t)-MaI(t)=M2KNaNe-ω1t-e-ω2tω2-ω1.
The solution is thus a combination of exponential decays with corresponding effective rate constants ω1 and ω2. Such biexponential behavior has been observed for the spontaneous transfer of certain lipids between phosphatidylcholine vesicles [41] and also for the release behavior of an imidazole derivate from liposomes [42]. The effective rate constants ω1 and ω2 can be calculated from G and K through(27)2G+K=ω1+ω2,4G2+K2=(ω2-ω1)2.
Hence, a measurement of ω1 and ω2 could be used to obtain the two model parameters (K and G). Figure 6  displays a plot of MdO(t), MdI(t), MaO(t), MaI(t), Ma(t), Md(t), calculated for G/K=1/10 and Na/N=Nd/N=0.5. All drug molecules are initially distributed equally among the two leaflets of the donor liposomes. Release of drug molecules from the outer leaf of the donor liposomes is fast (K=10G), the slow process is the flip-flop of drug molecules between the two leaflets of the liposomes. Hence, at intermediate times, say at t=3/K, the outer leaflets have almost reached their equilibrium values whereas the inner layers remain still fairly close to their initial values. After reaching thermal equilibrium (t→∞), half of all drug molecules have migrated to the acceptor liposomes. Clearly, the presence of the two different rate constants (K and G) leads to the biexponential behavior of Md and Ma in Figure 6.Figure 6
Fractions of drug molecules in inner and outer leaflets of donor and acceptor liposomes. The quantitiesMdO(t), MdI(t), MaO(t), and MaI(t) are plotted according to  (26) for G/K=1/10 and Na/N=Nd/N=0.5. The broken lines show the biexponential behaviors of the sums Md=MdO+MdI and Ma=MaO+MaI. The time is plotted in units of the inverse rate constant K. Note also ω1=1.11K and ω2=0.09K are the effective rate constants for the decay.We briefly discuss two limiting cases for (26). First, for G=0 the flip-flop of drug molecules between the inner and outer leaves is infinitely slow, implying MdI(t)=M/2, MaI(t)=0, MaO(t)=M/2-MdO(t)=(1-e-Kt)(MNa)/(2N). In this case, we recover the kinetics of (8), yet with only M/2 (instead of M) drug molecules participating in the transfer and identical donor and acceptor liposomes (k=0). Second, for G→∞ flip-flop becomes infinitely fast and (26) read MaI(t)=MaO(t)=M/2-MdI(t)=M/2-MdO(t)=(1-e-Kt/2)(MNa)/(2N). Because 50% of the drug molecules reside in the inner leaflets, they do not contribute to the outer-leaflet-concentration-differences that drive the transfer kinetics. Hence, the apparent rate constant is reduced from K to K/2.


## 3.1. Extension to High Drug Loading

While high drug loading obviously increases the number of available drug molecules (and thus increases the efficiency of liposomal carriers [39]) it also affects the kinetics of the drug release. Our present model predicts such a dependence for the diffusion mechanism whereas the kinetics for the collision mechanism is not affected. Recall that the transition from  (16) and (17) to (18) was based on the approximation of weak drug loading, Md≪mNd, Ma≪mNa, and M≪mN. Without that approximation, we obtain instead of (18) a nonlinear set of differential equations(20)Ṁd=-Kdrel⁡Na/N-(Ma/M)(M/mN)1-M/mNMd+Karel⁡Nd/N-(Md/M)(M/mN)1-M/mNMa,Ṁa=Kdrel⁡Na/N-(Ma/M)(M/mN)1-M/mNMd-Karel⁡Nd/N-(Md/M)(M/mN)1-M/mNMa.
For the special case that donor and acceptor liposomes are chemically similar, Kdrel=Karel=Kdiff, we obtain a simple exponential behavior(21)Ma(t)=M-Md(t)=(1-e-Kdifft/(1-M/mN))NaNM.
Here, high drug loading simply increases the rate constant for the diffusion mechanism by the factor 1/(1-M/(mN)). In the general case Kdrel≠Karel, and no simple exponential decay is predicted for high loading of the liposomes with drug molecules. Figure 4 shows a numerical example, based on (20) with Kdrel/Karel=3 and Nd/N=Na/N=0.5. For M≪mN (weak loading regime; broken lines in Figure 4) we observe the simple exponential behavior according to (18) with equilibrium values Mdeq/M=1/4 and Maeq/M=3/4. For M/(mN)=0.5 the initial loading of the donor liposomes is maximal. This leads to both a faster decay and a shift in the equilibrium distribution, reaching Mdeq/M=(3-1)/2=0.366 and Maeq/M=(3-3)/2=0.634. The reason for the increased rate constant is the reduced ability of highly loaded liposomes to take up drug molecules. Hence, if drug molecules are released from initially highly loaded donor liposomes they will be taken up exclusively by acceptor liposomes. The increase in the transfer rate at high loading also affects the equilibrium values Mdeq/M and Maeq/M. The equilibrium is shifted toward a more uniform distribution of drug molecules between donor and acceptor liposomes (in agreement with Figure 4).Figure 4
Numerical solutions of  (20), derived for M/(Nm)=0 (broken lines) and M/(Nm)=0.5 (solid lines). The remaining parameters are Kdrel/Karel=3, Nd/N=Na/N=0.5. The time t is plotted in units of 1/Karel.


## 3.2. Sigmoidal Behavior

Our model presented so far is unable to predict sigmoidal behavior. That is, no inflection point can be observed inMd(t) and Ma(t). Behind this prediction is our assumption that the transfer rates are strictly proportional to the concentration difference of the drug molecules. For the collision mechanism, this is expressed by our definition of the function g(i,j) in (3). However, if drug molecules within a given liposome interact with each other, the simple relation g(i,j)=i-j will no longer be valid. More specifically, attractive interactions between drug molecules within liposomes will increase the energy barrier to remove a drug molecule. This becomes relevant at high drug loading. Hence, in the presence of attractive interactions, it will be more unlikely that a drug molecule is transferred from a highly loaded donor liposome to an empty acceptor liposome.We discuss the consequences of attractive interactions for the collision mechanisms, which is described by (2) and (4). To account for the decrease in the rate constant at high loading we replace (3) by(22)g(i,j)=(i-j)(1-im)(1-jm).
Clearly, for weak loading (i≪m and j≪m) the original first-order model leading to the exponential behavior in (8) is recovered. For large loading of either donor or acceptor liposomes, the transfer rate becomes small. We note that using (22) does not lead to a set of differential equations in terms of only Md(t) and Ma(t). Here, we do not attempt to provide an analytical solution to the problem. Instead, we illustrate its predictions by numerically solving (2) and (4) with g(i,j) given in (22).Figure5  shows the behavior of Md(t) and Ma(t) as function of tK (with K=KcollN/V), derived for m=100. For simplicity, we have set k=0 which results in an equipartitioning of drug molecules between donor and acceptor liposomes (Md/Nd=Ma/Na=M/N). We start with Nd=Na=100 liposomes. The acceptor liposomes are initially empty whereas each donor liposome contains initially l drug molecules (out of a maximal number m=100). Different curves in Figure 5  correspond to l=2 (a), l=10 (b), l=50 (c), l=90 (d), and l=98 (e). As long as the drug loading is weak (curves (a) and (b)), the solution is simply exponential, characterized by Ma/M=1-Md/M=(1-e-Kt)Na/N (see (8) with k=0). Here, the kinetics is independent of the total number of drug molecules M=lNd (which is why curves (a) and (b) virtually overlap). If the initial loading of the donor liposomes becomes larger (curve (c)) the kinetics slows down. Eventually, once the initial loading approaches its maximal value mNd, the behavior slows down even more and, in addition, becomes sigmoidal. Attractive drug-drug interactions slow down the release from initially highly loaded donor liposomes; at later times (when the donor liposomes are no longer highly loaded), the release becomes faster. This leads to sigmoidal behavior.Figure 5
Fraction of drug molecules contained in donor liposomes (Md(t)/M; upper set of curves) and acceptor liposomes (Ma(t)/M; lower set of curves) as function of the scaled time Kt. The curves represent numerical solutions of  (2) and (4) with (22), derived for k=0 and m=100 with the initial conditions dj(0)=0 for j≠l, dl(0)=100, aj(0)=0 for j>0, a0(0)=100. Different curves correspond to l=2 (a), l=10 (b), l=50 (c), l=90 (d), and l=98 (e).


## 3.3. Extension to a Two-State Model

In the final part of this work, we briefly discuss an extension of our model to account for two distinct states of the drug molecule inside each liposome. A simple rationale for the presence of two distinct states is provided by the bilayer structure of the liposomes. That is, a drug molecule may preferentially be bound to either the inner or outer monolayer, having to flip-flop in order to change the host monolayer. The typical flip-flop time can be large if the drug has some amphiphilicity or surface activity instead of being strongly lipophilic [40]. Drug molecules residing in the inner monolayer cannot be transported directly to another liposome; they first have to migrate to the outer monolayer.We denote byMdI and MdO the number of drug molecules residing in the inner (DI) and outer (DO) leaflets of donor liposomes, respectively. Similarly, MaI and MaO refer to the number of drug molecules residing in the inner (AI) and outer (AO) leaflets of acceptor liposomes. The reaction scheme in (10) can then be generalized to account for the inter leaflet transport in donor and acceptor liposomes(23)DI⇌K2dK1dDO⇌K2K1AO⇌K1aK2aAI.
Here, K1d and K2d are the two rate constants corresponding to the transfer of drugs between the two leaflets of the donor liposomes (and similarly for K1a and K2a referring to the acceptor liposomes). The rate constants K1=(1-kNd/M)KNa/N and K2=(1+kNa/M)KNd/N are identical to those for the single-state model, where K is given in (19). Based on (23), the rate equations can be written as(24)ṀdO=KN(MaONd-MdONa+kNaNd)-K2dMdO+K1dMdI,ṀdI=K2dMdO-K1dMdI,ṀaO=KN(MdONa-MaONd-kNaNd)-K2aMaO+K1aMaI,ṀaI=K2aMaO-K1aMaI.
In the limit of a symmetric lipid bilayer, the two rate constants for flip-flop of a drug molecule from the inner to the outer leaf and from the outer to the inner leaf are identical (we note that the two leaflets of a liposomal bilayer are not strictly equivalent which, in a more refined model, would entail two different rate constants for flip-flop; this dependence on liposome curvature is neglected here). If we assume furthermore that donor and acceptor liposomes are chemically similar, we may write K1d=K2d=K1a=K2a=G as well as k=0. In this case, the rate equations(25)ṀdO=KN(MaONd-MdONa)-G(MdO-MdI),ṀdI=G(MdO-MdI),ṀaO=KN(MdONa-MaONd)-G(MaO-MaI),ṀaI=G(MaO-MaI)
depend on only two parameters, the rate constants K and G. If we assume all drug molecules initially reside in the donor liposomes, the initial conditions are MdO(t=0)=MdI(t=0)=M/2, and MaO(t=0)=MaI(t=0)=0, where M is the total number of drug molecules in the system. The solution of  (25) can be expressed as(26)MdI(t)=M2[NdN+NaNω2e-ω1t-ω1e-ω2tω2-ω1],MdO(t)-MdI(t)=M2KNaNe-ω2t-e-ω1tω2-ω1,MaI(t)=MNa2N[1-ω2e-ω1t-ω1e-ω2tω2-ω1],MaO(t)-MaI(t)=M2KNaNe-ω1t-e-ω2tω2-ω1.
The solution is thus a combination of exponential decays with corresponding effective rate constants ω1 and ω2. Such biexponential behavior has been observed for the spontaneous transfer of certain lipids between phosphatidylcholine vesicles [41] and also for the release behavior of an imidazole derivate from liposomes [42]. The effective rate constants ω1 and ω2 can be calculated from G and K through(27)2G+K=ω1+ω2,4G2+K2=(ω2-ω1)2.
Hence, a measurement of ω1 and ω2 could be used to obtain the two model parameters (K and G). Figure 6  displays a plot of MdO(t), MdI(t), MaO(t), MaI(t), Ma(t), Md(t), calculated for G/K=1/10 and Na/N=Nd/N=0.5. All drug molecules are initially distributed equally among the two leaflets of the donor liposomes. Release of drug molecules from the outer leaf of the donor liposomes is fast (K=10G), the slow process is the flip-flop of drug molecules between the two leaflets of the liposomes. Hence, at intermediate times, say at t=3/K, the outer leaflets have almost reached their equilibrium values whereas the inner layers remain still fairly close to their initial values. After reaching thermal equilibrium (t→∞), half of all drug molecules have migrated to the acceptor liposomes. Clearly, the presence of the two different rate constants (K and G) leads to the biexponential behavior of Md and Ma in Figure 6.Figure 6
Fractions of drug molecules in inner and outer leaflets of donor and acceptor liposomes. The quantitiesMdO(t), MdI(t), MaO(t), and MaI(t) are plotted according to  (26) for G/K=1/10 and Na/N=Nd/N=0.5. The broken lines show the biexponential behaviors of the sums Md=MdO+MdI and Ma=MaO+MaI. The time is plotted in units of the inverse rate constant K. Note also ω1=1.11K and ω2=0.09K are the effective rate constants for the decay.We briefly discuss two limiting cases for (26). First, for G=0 the flip-flop of drug molecules between the inner and outer leaves is infinitely slow, implying MdI(t)=M/2, MaI(t)=0, MaO(t)=M/2-MdO(t)=(1-e-Kt)(MNa)/(2N). In this case, we recover the kinetics of (8), yet with only M/2 (instead of M) drug molecules participating in the transfer and identical donor and acceptor liposomes (k=0). Second, for G→∞ flip-flop becomes infinitely fast and (26) read MaI(t)=MaO(t)=M/2-MdI(t)=M/2-MdO(t)=(1-e-Kt/2)(MNa)/(2N). Because 50% of the drug molecules reside in the inner leaflets, they do not contribute to the outer-leaflet-concentration-differences that drive the transfer kinetics. Hence, the apparent rate constant is reduced from K to K/2.


## 4. Conclusions

In this work, we have presented a detailed model for the transfer kinetics of poorly water-soluble drug molecules between liposomal carrier systems. Apart from liposomes, the scope of the model includes other types of small and mobile pharmaceutical nanocarriers, such as micelles, colloids, and nanoparticles. Starting from a microscopic distribution function of drug molecules among donor and acceptor liposomes, we have specified the conditions that lead to an apparent first-order kinetic behavior. These include low drug loading of the liposomes, strict proportionality of all rate constants to drug concentrations, no aggregation phenomena of drugs within liposomes, and no overlap of the intraliposomal flip-flop kinetics. Systems that do not fulfill these conditions do not, generally, exhibit an apparent first-order kinetics. Instead the behavior may become biexponential or sigmoidal. High drug loading may preserve the first order kinetics but with increased apparent rate constant.An optimal drug delivery system should keep the drug load on the way to the target and release it only after arrival at the target. Understanding the kinetics and mechanisms of drug release from liposomal (and other) nanocarriers is thus a prerequisite to systematically improving drug delivery systems.



---

*Source: 376548-2011-06-07.xml*