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.

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 [912]. 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 volume ) that contains a number of donor and 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 of 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, , and by the acceptor liposomes, . It is then initially and , as well as in equilibrium and , where and 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., for and ) 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.

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 as spatially 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 number of donor liposomes that carry drug molecules. An analogous definition is used for the number of acceptor liposomes that carry drug molecules. The index is confined to the region where is the maximal number of drug molecules that a liposome can carry. The time-dependent distribution functions and represent a full microscopic knowledge of the kinetics of drug transfer. The total numbers of donor liposomes , acceptor liposomes , drug molecules residing in donor liposomes , and drug molecules residing in acceptor liposomes , can be calculated on the basis of the distribution functions and according to Mathematically, and are the zeroth-moments of the distributions functions and whereas and appear as the corresponding first moments. We assume that and are constant (i.e., independent of time), and so then is the total number of liposomes . 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, , 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.

Collisions require two liposomes to come to close proximity. The magnitude of drug transport between, say, donor liposomes and is thus where 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 drug molecules to a donor that carries initially drug molecules, the distribution function changes according to , , , and . Hence, the numbers and decrease whereas and increase. Figure 3 shows an illustration of this scheme for and . The transfer rate between the populations and will also depend on the corresponding numbers of drug molecules and . We assume the drug molecules within each liposome form an ideal mixture so that the transfer rate is directly proportional to . In writing a rate equation for donor population , we have to account for all possible ways of collisions between donor liposomes of index with other liposomes (donors and acceptors) of index . These considerations lead us to where we have defined the function In (2), is the unit rate of drug transfer through collisions between two chemically equivalent liposomes, and denotes the time derivative of a physical quantity . 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.

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 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 each donor and acceptor liposome will contain the same number of drug molecules in thermal equilibrium). We do not attempt to calculate from a microscopic model; yet below we show how 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 obtain from by replacing , , and , 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 implying and thus ensuring conservation of the number of donor and acceptor liposomes ( and ) as well as of the total number of drug molecules (). To characterize the total numbers and of drug molecules that reside in donor and acceptor liposomes, respectively, we carry out the summations and using (2) and (4). The result are the two first-order differential equations where we have introduced the definition of the apparent rate constant Initially, all drug molecules are incorporated in the donor liposomes, implying and . The solution of (6) is then Hence, indeed appears as the inverse characteristic time for the transfer process. In contrast to previous models [14], depends only on the total concentration of liposomes 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 and ). This includes but is not restricted to sink conditions (where ).

Thermodynamic equilibrium corresponds to the long-time limit, , at which we have and with From (9), we obtain the difference between the numbers of drug molecules per donor and acceptor liposome, . This agrees with our interpretation of in (2) and (4). We note that for chemically identical donor and acceptor liposomes, it is and all liposomes carry the same number of drug molecules in equilibrium, implying . The largest possible value of is for which we obtain and . The smallest possible value of is implying and . Hence, .

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 with rate constants and . The corresponding kinetic behavior is then governed by the equations and where and are the numbers of drug molecules carried by donor and acceptor liposomes, respectively. With and we obtain which has indeed the same structure as (8). Comparison of (8) with (11) reveals and . The equilibrium constant of the reaction in (10) is thus Comparing this with (where is Boltzmann's constant and is the absolute temperature) allows us to compute the change in standard Gibbs free energy for the transfer of a single drug molecule from a donor to an acceptor liposome. The enthalpic and entropic contributions to will be influenced by , which is, generally, temperature dependent. Let us briefly discuss two cases. First, if donor and acceptor liposomes are chemically identical, then and has only an entropic contribution. Specifically, for , we find because a given drug molecule has more donor liposomes to reside in than acceptor liposomes. Second, the limiting cases for , namely, and , yield (thus, with all drugs migrating to the acceptor liposomes) and (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-tail in their concentration dependence. For example, the kinetic behavior of the dimerization reaction 2 monomer → dimer follows the equation where is the concentration of the reactant (i.e., the monomers) and the rate constant. With an initial concentration the time behavior becomes , implying 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 and 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] with rate constants , , , and 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 . We assume the probability of a drug molecule to leave donor liposomes of index to be proportional to the total number of drug molecules in that liposome population. Similarly, the probability of a drug molecule to enter donor liposomes of index is assumed to be proportional to the total number 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 in the aqueous phase (here, is the total number of drug molecules residing in the aqueous phase). This leads to the following rate equations for (with for or ). A similar equation can be written for the acceptor liposomes. Based on (15), it can be verified that , thus ensuring conservation of (and similarly for ). Carrying out the summation using (15) leads to 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) 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, . 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 , , and . We finally assume the same rate for the uptake of drug molecules from the aqueous phase into donor and acceptor liposomes, implying . 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 Equation (18) are now identical to (6) if we identify and where appears as the rate constant. Here again, as for (6), the validity of (18) is not subject to a restriction with respect to and .

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 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 . 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, , , and . Without that approximation, we obtain instead of (18) a nonlinear set of differential equations For the special case that donor and acceptor liposomes are chemically similar, , we obtain a simple exponential behavior Here, high drug loading simply increases the rate constant for the diffusion mechanism by the factor . In the general case , 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 and . For (weak loading regime; broken lines in Figure 4) we observe the simple exponential behavior according to (18) with equilibrium values and . For the initial loading of the donor liposomes is maximal. This leads to both a faster decay and a shift in the equilibrium distribution, reaching and . 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 and . The equilibrium is shifted toward a more uniform distribution of drug molecules between donor and acceptor liposomes (in agreement with Figure 4).

3.2. Sigmoidal Behavior

Our model presented so far is unable to predict sigmoidal behavior. That is, no inflection point can be observed in and . 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 in (3). However, if drug molecules within a given liposome interact with each other, the simple relation 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 Clearly, for weak loading ( and ) 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 and . 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 given in (22).

Figure 5 shows the behavior of and as function of (with ), derived for . For simplicity, we have set which results in an equipartitioning of drug molecules between donor and acceptor liposomes (). We start with liposomes. The acceptor liposomes are initially empty whereas each donor liposome contains initially drug molecules (out of a maximal number ). Different curves in Figure 5 correspond to (a), (b), (c), (d), and (e). As long as the drug loading is weak (curves (a) and (b)), the solution is simply exponential, characterized by (see (8) with ). Here, the kinetics is independent of the total number of drug molecules (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 , 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.

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 by and the number of drug molecules residing in the inner () and outer () leaflets of donor liposomes, respectively. Similarly, and refer to the number of drug molecules residing in the inner () and outer () 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 Here, and are the two rate constants corresponding to the transfer of drugs between the two leaflets of the donor liposomes (and similarly for and referring to the acceptor liposomes). The rate constants and are identical to those for the single-state model, where is given in (19). Based on (23), the rate equations can be written as 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 as well as . In this case, the rate equations depend on only two parameters, the rate constants and . If we assume all drug molecules initially reside in the donor liposomes, the initial conditions are , and , where is the total number of drug molecules in the system. The solution of (25) can be expressed as The solution is thus a combination of exponential decays with corresponding effective rate constants and . 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 and can be calculated from and through Hence, a measurement of and could be used to obtain the two model parameters ( and ). Figure 6 displays a plot of , , , , , , calculated for and . 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 (), the slow process is the flip-flop of drug molecules between the two leaflets of the liposomes. Hence, at intermediate times, say at , 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 (), half of all drug molecules have migrated to the acceptor liposomes. Clearly, the presence of the two different rate constants ( and ) leads to the biexponential behavior of and in Figure 6.

We briefly discuss two limiting cases for (26). First, for the flip-flop of drug molecules between the inner and outer leaves is infinitely slow, implying , , . In this case, we recover the kinetics of (8), yet with only (instead of ) drug molecules participating in the transfer and identical donor and acceptor liposomes (). Second, for flip-flop becomes infinitely fast and (26) read . 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 to .

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.

Acknowledgments

The authers thank Drs. Alexander Wagner, Martin Holzer, and Rolf Schubert for illuminating discussions. S. May acknowledges support from NIH through Grant GM077184.