Guy
Jacoby
a,
Irina
Portnaya
b,
Dganit
Danino
b,
Haim
Diamant
c and
Roy
Beck
*a
aThe Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, Ramat Aviv, Tel Aviv 6997801, Israel. E-mail: roy@tauex.tau.ac.il
bDepartment of Chemical Engineering, Technion-Israel Institute of Technology, Haifa 3200003, Israel
cThe Raymond and Beverly School of Chemistry, Tel Aviv University, Tel Aviv 6997801, Israel
First published on 19th November 2019
Metastable states in first-order phase-transitions have been traditionally described by classical nucleation theory (CNT). However, recently an increasing number of systems displaying such a transition have not been successfully modelled by CNT. The delayed crystallization of phospholipids upon super-cooling is an interesting case, since the extended timescales allow access into the dynamics. Herein, we demonstrate the controllable behavior of the long-lived metastable liquid-crystalline phase of dilauroyl-phosphatidylethanolamine (DLPE), arranged in multi-lamellar vesicles, and the ensuing cooperative transition to the crystalline state. Experimentally, we find that the delay in crystallization is a bulk phenomenon, which is tunable and can be manipulated to span two orders of magnitude in time by changing the quenching temperature, solution salinity, or adding a secondary phospholipid. Our results reveal the robust persistence of the metastability, and showcase the apparent deviation from CNT. This distinctive suppression of the transition may be explained by the resistance of the multi-lamellar vesicle to deformations caused by nucleated crystalline domains. Since phospholipids are used as a platform for drug-delivery, a programmable design of cargo hold and release can be of great benefit.
However, there is an increasing number of systems displaying nucleation processes that do not conform to this classical picture. Complex dynamics can arise due to intermediate states leading to multi-step nucleation,8–10 or long-range interactions that can result in macroscopic nucleation11 or a cooperative delayed transition.12 Such dynamics may require an extension of the classical theory or in some cases a comprehensive revision.13 Examples of such complex dynamics can be found in self-assembled amphiphilic systems, which display long-lived metastable phases upon temperature change.
Amphiphiles are molecules that contain hydrophilic and hydrophobic chemical groups. The key characteristic associated with amphiphilic molecules is their ability to spontaneously self-assemble into macro-molecular structures.14 Biological amphiphiles, such as lipids, self-assemble into a wide variety of mesophases. Most notably, lipids constitute the membranes of cells and organelles, and are involved in many important biological functions. Alongside basic research into their physical and biochemical properties, self-assembly is also utilized for designing modern biomedical applications such as drug delivery.15
In particular, the phospholipid amphiphiles have several predominant lamellar phases, such as the disordered liquid-crystalline phase (Lα), gel phase (Lβ) and ordered crystalline phase (Lc), which differ by their degree of spatial symmetry. Transitions between these phases can be induced by changing the temperature, but the pathways depend on the physicochemical properties of the molecules and their thermal history. Previously, the phospholipid DLPE was shown to have long-lived metastable phases with lifetimes on the order of hours or even days.16–18 These time-scales are orders of magnitude longer than the rapid transitions as in the case of melting. The reports were based mostly on X-ray scattering and differential scanning calorimetry (DSC), both very useful techniques for phase detection and characterization. However, they were performed as static measurements at different points in time, separated by long periods of unrecorded incubation. These limited observations only allowed for qualitative descriptions of the dynamics.
Herein, we present our experimental investigation of the metastable Lα to Lc phase-transition using time-resolved solution X-ray scattering (SXS) and DSC measurements. We demonstrate the cooperative and controllable behavior of the transition dynamics. We highlight and discuss the deviations from CNT, and rationalize them based on the free-energy cost of deforming the Lα vesicles by crystal nucleation.
Here, the metastable Lα to Lc phase-transition was recorded by time-resolved SXS. In the experiments, the X-ray beam illuminating the sample has a cross sectional area of approximately 0.64 mm2, which produces a bulk-averaged scattering signal. The scattering from a sample that has not been pre-heated is first recorded at TQ (Tβ→α < TQ < Tc→α) as a reference point to the initial low-energy state of the system. The scattering pattern of this initial state pertains to a 3D crystal with an orthorhombic unit cell, with the largest dimension corresponding to the lamellar repeating distance, and the shorter ones to the in-plane ordering of the lipids. The scattering from the sample is then recorded after heating to 60 °C (above Tc→α), where the lipid membranes have transitioned to the smectic Lα phase, i.e. a lamellar structure with liquid-like disorder in-plane. The phase change is accompanied by a change in morphology of the lipid particles, from faceted crystals to curved multi-lamellar vesicles (MLVs).18 The in-plane order-to-disorder transition is marked by the disappearance of the correlation peaks at wide scattering angles (Fig. 1).
After three hours at 60 °C, the sample is quenched (at a rate of 2.5 °C min−1) back to TQ and measured every hour until it has finished transitioning back to Lc. The time-resolved scattering spectra show the evolution of the correlation peaks, which reports on structural changes and the phase-transition (Fig. 1a). This heating–cooling procedure is the typical experiment conducted to record the metastability. Surprisingly, regardless of the parameters changed between samples, such as buffer salinity or the inclusion of a secondary lipid into the system, there are several prominent features in the evolution of the scattering spectrum after cooling back to TQ. The first is the extended period of time in which the Lα phase remains metastable, which we denote as the delay time τ (Fig. 1b). The second is the time the system spends transitioning between phases, denoted as the transition period τ*, which is by and large an order of magnitude shorter than τ. The third notable feature occurring in most of the experiments accounted for in this work is a substantial change in the lamellar scattering intensity during the transition period. The intensity remains constant during the delay time, however it decreases, up to an order of magnitude, during the transition period. This indicates a decrease in the average number of lamellae per particle during the transition. This implies that the phase-transition is accompanied by a macro-scale structural change of the particles.
The delay time until crystallization (τ) and the duration of the ensuing transition (τ*) are extracted from the time-dependent correlation peak intensities, which are a direct measurement of the order parameters. Specifically, we fit a sigmoid function to the time-dependent intensity of the lamellar correlation peak (001) and to a mixed correlation peak in the wide angles (206) (Fig. 1b), to measure the out- and in-plane order parameters respectively.
In the classical nucleation theory, a single timescale should be observed, namely, that associated with the rate of nucleation. However, given that SXS produces bulk-averaged signals, one can immediately notice that the dynamics presented in the X-ray spectra do not readily conform to the classical picture of a single stochastic process culminating in a phase-transition. Instead of a gradual increase in crystalline scattering we observe two distinct timescales, τ and τ*, which are orders of magnitude larger than the melting transition times, with no detectable scattering at wide angles prior to the transition. This implies a collective bulk transition rather than stochastic events of crystallization within the macroscopic illuminated area.
We set out to explore the different system parameters that can affect the dynamics of the transition. By changing system parameters such as the lipid stoichiometry and chemical structure, salinity and the quenching temperature (TQ) we found that we were able to manipulate the metastability in a pre-determined and controllable fashion.
The lifetime of the metastable phase depends on the strength of the thermodynamic force driving the transition. Close to the transition temperature the energy barrier that the system must overcome is high, which results in a slow rate of nucleation. As the temperature is lowered, the barrier becomes smaller and the rate increases. This is demonstrated here by the increasing persistence of the metastable Lα phase, closer to Tc→α. Fig. 2 shows an exponential increase in the average delay time as a function of the quenching temperature, TQ, for samples of pure DLPE. Remarkably, for most quenching temperatures the spread of experimental results is very small, further supporting our claim that τ is an intrinsic property of the system's dynamics set by its macroscopic parameters. However, due to larger fluctuations closer to the critical temperature (Tc→α) we notice a large spread of τ at TQ = 41 °C. There, two samples did not transition within the duration of the experiment (800 hours). In addition, at TQ = 31 °C the Lα phase rapidly transitioned to Lβ, which then became metastable. Evidently, the delay time for the Lα phase at TQ = 32 °C is shorter than at TQ = 31 °C for the Lβ gel phase, as the liquid phase is expected to be more labile.20
Fig. 2 Delay time as a function of quenching temperature for pure DLPE. Empty squares represent measurements of individual capillaries at each quenching temperature. Solid red curve is the fit using eqn (2), excluding the measurements at TQ = 31 °C (empty stars represent the Lα → Lβ transition) and TQ = 41 °C due to the high variance caused by fluctuations. Red vertical dashed-line indicates the estimated Tc→α transition temperature. |
In previous studies, the metastability was examined in samples of pure DLPE.16,17 However, we found that the addition of a secondary phospholipid not only preserves the metastable phase, but also extends its lifetime (Fig. 3). Moreover, the delay time is sensitive to changes in the hydrocarbon chain length and headgroup. Phosphatidylglycerols (PGs) were chosen as a secondary charged lipid due to the stabilizing effect they have on PE bilayers,21 and specifically, DLPG was previously used along with DLPE as the building blocks for a drug delivery system.22–24 PGs with 12 (DLPG), 14 (DMPG), 16 (DPPG) or 18 (DSPG) carbons in their saturated hydrocarbon chains were chosen as chain length variants. In addition, the zwitterionic dilauroyl-phosphatidylcholine (DLPC) was chosen as a headgroup variant. The delay time seems to increase as a function of chain length for the PGs, and it is greatly increased when the PG headgroup is swapped with a PC.
Since electrostatics are known to have a central role in stabilizing lipid lamellar systems, and the delay time is observed to increase with the fraction of charged PG lipids in our system, we tested the effect of the solution salt concentration on the delay time of samples with different DLPE:DLPG ratios. As shown in Fig. 4, changing the average membrane charge density produces two features in the salt dependence: (a) there seems to be a minimum of the delay time at approximately 150 mM, splitting the dependency into two regimes, and (b) the delay time increases with the fraction of DLPG at a given salt concentration.
We propose that these two regimes originate from two different phenomena. At low salt concentrations (<150 mM), the decrease in delay time towards the minimum can be attributed to the decrease in the electrostatic screening length. At 150 mM the electrostatic screening length is comparable to the DLPG headgroup diameter (≈8 Å).25 Segregation of non-DLPE lipids, which is essential for recovering the homogeneous DLPE crystals, is facilitated by the screening of interacting charged PG headgroups. On the contrary, high salt concentrations (>300 mM) can lead to adsorption of ions on the charged membrane. This can lead to an increase in τ, since ions must evacuate from between the lamellae, yet ion transport across membranes is unfavored. The adsorption of charges can result in an increase in the membrane's bending rigidity, which in turn can strengthen the metastability (see Section 4.2). The results show that the delay time of samples containing charged headgroups responded to changes in salt concentration, yet no significant dependence was observed in samples containing only DLPE. We would like to accentuate the extended lifetime of the metastable phase at the highest concentration measured (500 mM), which exceeded 500 hours in the case of 90:10 DLPE:DLPG (mole%).
The structural study of the delayed nucleation phenomenon, using time-resolved SXS, does not directly report on the thermodynamic processes. To address this, calorimetric measurements are commonly used to investigate the thermodynamics of lipid systems, mostly in the form of differential scanning calorimetry (DSC). However, since we are investigating a time-delayed transition at a fixed temperature, we employed DSC in a non-trivial quasi-isothermal manner. After samples were incubated at 60 °C for 3 hours, the temperature was lowered to TQ = 37 °C, and the samples were scanned back-and-forth between 36 and 37 °C at a very slow rate (0.1 °C hour−1 on heating, 0.43 °C hour−1 on cooling).
In Fig. 5 we compare the quasi-isothermal energy flux measurements performed on a sample of pure DLPE and a sample containing 90:10 DLPE:DLPG (mole%). The results in both cases show an exothermic signature, as expected for the transition back to the low-energy crystalline phase, with similar timescales to those in our X-ray scattering measurements (Fig. 3). The sample with pure DLPE shows a broad and slow change in the excess heat capacity, peaking at 28 hours, while the mixed sample shows a much narrower peak, centered at 55 hours. The corresponding average delay times in the X-ray scattering experiments are 20 and 45 hours, respectively. In addition, a DSC measurement was performed to determine the enthalpy of transition from the Lc to the Lα phase for a pure sample of DLPE, which yielded hc = 11.1 kcal mole−1.
Fig. 5 Time-resolved quasi-isothermal DSC measurements of two samples: (a) pure DLPE, (b) 90:10 DLPE:DLPG (mole%). The exothermic peak seen in both occurs at similar times as in the scattering experiments. The red curves (shown in both) are heating scans and the blue curves shown only in (a) are cooling scans. There is a mismatch between cooling and heating scan data possibly due to the different scanning rates (see Experimental section), thus the cooling scans were shifted in (a) and omitted from (b), for clarity. Raw data are shown in ESI† (Fig. S1). |
In the process of sample measurement and analysis of the phase-transition dynamics, the extracted temporal parameters represent the transformation occurring in the illuminated volume. However, this volume includes only a portion of the lipid pellet at the bottom of the capillary. If one follows a nucleation and growth framework, it is important to assess whether the transformation initiates concurrently throughout the sample or propagates successively from a starting point. To test this, we prepared a sample with a large pellet of several millimeters at the bottom of the capillary and measured the delay time at different locations along the vertically held capillary. The capillary size and experimental protocol (SXS) remained as previously described. Fig. 6a shows the delay time as a function of the spatial coordinate along the capillary. The transition seems to propagate outward from a certain location, with neighboring locations transitioning at later times.
A similar experiment was performed on a horizontally held capillary. There, the transition began at the water–pellet interface and propagated at a steady velocity of approximately 100 μm hour−1 towards the end of the pellet (Fig. 6b). In both experiments a new feature in the scattering spectra could be observed: a coordinated increase in lamellar scattering over a period of time prior to the transition (Fig. 6c and d). We denote the beginning of this period by τB, the point in time from which a slow increase in scattering culminates in a sharp drop of intensity. During this period, the amplitude of the peak increases significantly, which implies an increase in the average number of scattering lamellae. In addition, the width of the peak slightly increases as well (Fig. S2, ESI†). The width of the scattering peak is inversely proportional to the number of lamellae in correlation per membrane stack.26 Therefore, if stacks of correlated membranes would grow during the period before the transition (starting at τB), the peak amplitude would increase, but the width would decrease. However, the data shows an increase of both. This implies an increase in the number of different, independent stacks of correlated membranes, prior to crystallization. Only after the drop in the lamellar scattering intensity is there a detectable change in wide-angle scattering. Therefore, the metastable state remains during the build-up period. Surprisingly, this structural reorganization is coordinated over several millimeters in the sample (Fig. 6c and d).
Metastable phases are often very sensitive to energy fluctuations, as even minute inputs of energy can result in a transition to the stable phase. Since the lipid metastable phase is stable against various changes in system parameters, we tested its stability against external inputs of energy by subjecting lipid dispersions to mechanical agitation in the form of rigorous pipetting. A lipid dispersion of approx. 1.5 ml, at 30 mg ml−1, was prepared as a bulk dispersion from which samples would be pipetted out, and measured intermittently. It was incubated at 37 °C for one hour, followed by 3 hours at 60 °C, as performed regularly with the SXS samples. The incubator was then set to TQ = 37 °C and a sample was drawn from the bulk dispersion after t = 1, 2, 3, 5.5 and 19.5 hours by pipetting out approximately 100 μl and placing into a capillary. The capillary was then placed in the SXS temperature chamber, pre-heated to 37 °C, and measured after Δt minutes (Fig. 7). The control sample, taken from the bulk dispersion before it was placed in the incubator, underwent the regular temperature procedure in the SXS temperature chamber, and transitioned after τ = 19.5 hours. The samples from the bulk dispersion, taken during first few hours, transitioned approximately an hour after being pipetted out. The sample taken after 5.5 hours was in the middle of transitioning when measured initially (Δt = 0 min), and the 19.5 hour sample had already transitioned (Fig. 7). This experiment demonstrates that the lifetime of the metastable phase is significantly shortened by mechanical agitation applied after thermal incubation.
Lastly, phospholipids are utilized in bio-medicine as building blocks of vesicles designed for specific targeting and controlled release. When designing such drug-delivery systems, it is crucial to assert the stability of the carrier with its cargo. Since DLPE and DLPG have been used as the lipid components of such systems,22–24 we tested the stability of the metastable phase in the presence of cargo. Fig. 8 shows the delay time of samples containing 90:10 DLPE:DLPG (mole%) and the hydrophobic drug Prednisolone; an established and commercially available steroid used to treat a wide range of conditions and illnesses. The results show that the addition of the drug had a large impact on the delay time, shortening it by almost an order-of-magnitude. Since lipid systems continue to serve as appealing ingredients in drug delivery systems, controlled delayed nucleation may serve as a novel designing factor to deposit cargo in a predetermined timing. Nonetheless, the effect of cargo on transition dynamics should not be overlooked when designing such systems.
Fig. 8 The delay time decreases in the presence of prednisolone, a commercial hydrophobic drug, but the delay in transition still persists. |
We present additional findings that we would like to discuss in the context of deviations from CNT. Firstly, it is important to mention that at no point in the preparation of the samples was there any effort to homogenise the particles' sizes. And yet, despite this heterogeneity, the delay time was shown to be reproducible with a peak in the probability at a non-zero value.18 Secondly, our results shown in Fig. 6 demonstrate that the transition occurs at different times at different locations, yet the onset of the structural change in samples with large pellets, marked by τB, is macroscopically coordinated over millimeters. This length-scale is orders of magnitude larger than any microscopic length-scale associated with lipid self-assembly. Lastly, in heterogeneous nucleation impurities are considered preferential nucleation sites due to a lower surface energy penalty compared to the homogeneous case. In our system, the inclusion of a secondary lipid only served to hinder crystallization.
During the incubation time at the high temperature (≥60 °C) water molecules and ions enter in between the membranes.18 Concurrently, secondary lipid molecules enter the liquid membrane and disrupt its homogeneity. Upon cooling, the water molecules and ions must evacuate, and lipid segregation must occur to re-form the network of connections as in the initial homogeneous crystal. Here, we demonstrated that the persistence of the metastable phase is sensitive to the properties of the secondary lipid. Not only does the metastability with its features persist, but the chemical structure of the lipid has a large impact on the change of the delay time (Fig. 3).
A recent study experimentally showed long-range interlayer alignment of phase-separated intralayer domains, across hundreds of lamellae in multi-component supported lipid membranes.27 A follow-up study proposed a theoretical explanation to the interlayer correlation between phase-separated domains.28 Using a model of stacked 2D Ising spins to represent the stacked lipid membranes, they showed that the system forms a continuous columnar structure in equilibrium, for any finite interaction across adjacent layers. Such an interlayer interaction should be a key component in cooperative nucleation in MLVs, for which a mechanism is proposed below.
Because of the rigidity of the crystalline phase, the formation of a crystalline domain in a membrane flattens that region, which affects adjacent membranes and thus locally deforms the MLV. The free-energy penalty per unit area due to the deformation is proportional to the effective surface tension of the MLV, arising from its compression (B) and bending (K) moduli.29 This penalty makes the free-energy of the crystalline phase effectively higher, shifting the transition from T(0)c for an isolated membrane to a lower temperature Tc for a membrane in a curved MLV. The change in the free-energy per unit area between the two phases is Δg(Tc) ≃ 2γ/3. Using the relation Δg ≃ (hc/a)(T(0)c − T)/T(0)c, the temperature shift can be estimated from the measured enthalpy of transition (hc = 11.1 kcal mole−1), the area per lipid (a ≃ 0.5 nm2), and assuming γ larger than 0.1 mN m−1. We get a decrease in the transition temperature, proportional to γ, of more than 10 degrees.
Thus, under conditions where an isolated membrane would crystallize, a single membrane in the MLV would not. On the other hand, if all the membranes in the MLV were to crystallize, the total free-energy would inevitably decrease. Hence, the MLV must ultimately crystallize, but it can do so only through a multi-membrane cooperative process. This cooperativity is essential for departing from a single Poisson process, typical to CNT.
The compression and/or bending moduli of the Lα phase, and therefore also γ, increase with increasing lipid chain length and increasing membrane charge. Such modifications, by increasing γ, should deepen the metastability. This is consistent with the experimental observations reported here, which show the increase of τ with an increase in the fraction of charged headgroups or chain length (Fig. 3).
Let us consider the collective metastability described above in slightly more detail. Since the deformation of the membrane stack is localized,29 its free-energy penalty is intensive in the number of membranes, whereas the free-energy gain due to crystallization is extensive. Hence, there is a critical number of membranes, nc ∼ γ/(Δg), beyond which the suppression of the transition is overcome. The multi-membrane critical nucleus has the size Rc ∼ λnc, where λ is the MLV's periodicity. It corresponds to a multi-membrane nucleation barrier, Fc ∼ σRc2 ∼ σγ2λ2/(Δg)2, where σ is an effective surface tension; a combination of γ and the line tension of intra-membrane crystalline domains.
To sum up this analysis, we obtain the multi-membrane nucleation barrier as
(1) |
τ = τ0eFc/(kBT) = τ0eb/(Tc−T)2, | (2) |
Fig. 2 shows the fit of the function in eqn (2) to the data, where Tc was fixed to 43 °C. The fitting parameter b yields a value of ≈22 K2, consistent with the qualitative estimate shown above. Changing Tc to be 42 or 44 °C changes b to be approx. 8 or 44, respectively. A detailed theory along the lines presented in this section will be presented elsewhere.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9sm01834d |
This journal is © The Royal Society of Chemistry 2020 |