Inge
Bos
,
Eline
Brink
,
Lucile
Michels
and
Joris
Sprakel
*
Physical Chemistry and Soft Matter, Wageningen University & Research, Stippeneng 4, 6708 WE, Wageningen, The Netherlands. E-mail: joris.sprakel@wur.nl
First published on 17th February 2022
Single stranded DNA (ssDNA), or another polyanion, can be mixed with polycations to form liquid-like complex coacervates. When the polycations are replaced by cationic–neutral diblock copolymers, complex coacervate core micelles (C3Ms) can be formed instead. In both complex coacervates and C3Ms, dynamics plays an important role. Yet, to date, the effect of chain length on the dynamics effect is still not fully understood. The DNA complexes provide a versatile platform to further elucidate these chain length effects because the DNA is monodisperse and its length can be easily adapted. Therefore, we study in this paper the dynamics of fluorescently labelled ssDNA in both complex coacervate droplets and micelles. The DNA dynamics in the complex coacervate droplets is probed by fluorescence recovery after photobleaching (FRAP). We observe that the DNA diffusion coefficient depends more strongly on the DNA length than predicted by the sticky Rouse model and we show that this can be partly explained by changes in complex coacervate density, but that also other factors might play a role. We measure the molecular exchange of C3Ms by making use of Förster resonance energy transfer (FRET) and complement these measurements with Langevin dynamics simulations. We conclude that chain length polydispersity is the main cause of a broad distribution of exchange rates. We hypothesise that the different exchange rates that we observe for the monodisperse DNA are mainly caused by differences in dye interactions and show that the dye can indeed have a large effect on the C3M exchange. In addition, we show that a new description of the C3M molecular exchange is required that accounts among others for the effect of the length of the oppositely charged core species. Together our findings can help to better understand the dynamics in both specific DNA systems and in complex coacervate droplets and micelles in general.
When DNA is mixed with a cationic–neutral diblock copolymer, it can form complex coacervate core micelles (C3Ms) instead of complex coacervate droplets: the repulsive interactions between the neutral blocks prevent further growth of the complex coacervate and in this way thermodynamically stable nanostructures can be formed instead of the macroscopic complex coacervates that will eventually form when DNA is mixed with polycations without neutral blocks. The C3M core is formed by the DNA and the cationic blocks, while the neutral blocks form the micelle corona. The micelle corona forms a protective layer around the micelle core and prevents micelle coalescence.10 In this way, micelles with specific polymer aggregation numbers and sizes of typically 10 to 100 nm are formed11–13 that can protect their core components from external components. Their protective properties as well as their well-defined small size makes these DNA complex coacervate core micelles promising DNA-based medicine delivery tools.14–17
Studying the liquid-like DNA polyelectrolyte complexes is not only useful to better understand the formation of membraneless organelles and to improve the formation of DNA-based medicine delivery tools, but also to give fundamental insights into polyelectrolyte complexes in general. In particular, the effect of chain length and chain length polydispersity can be systemically studied with DNA because specific DNA sequences can be synthesised. This yields monodisperse DNA whose length can be systematically varied by changing the sequence. At the moment, the chain length effects on the complex coacervates and C3Ms are not completely understood. Especially the chain length effect on the coacervate dynamics is not well understood, as we will explain below.
The dynamics of complex coacervates determines the response of complex coacervate materials to deformation18,19 and might determine their response time to dissociation triggers.20 The chain length effect on the dynamics in complex coacervates is usually described by the sticky Rouse model21 where the ionic bonds act as sticky points that slow down the dynamics.18,19 In the sticky Rouse model, the overall polymer relaxation time scales with the squared number of intermolecular sticky bonds, which means for polyelectrolytes that the relaxation time should scale with the square of the polyelectrolyte length. For complex coacervates with matched chain length this description seems to work well,18 but for complexes where the lengths of the anionic and cationic chain differ, the situation is more complex. In some cases, the dynamics seems to be solely governed by one of the two polyelectrolytes, while the length of the other polyelectrolyte does not have any effect.18 Since only few studies have focused on polyion pairs with incommensurate lengths, an understanding of these asymmetry effects, and with that a complete description of coacervate dynamics, is lacking.
The molecular exchange dynamics of complex coacervate core micelles determines how often the core components are exposed to the surroundings and in this way the level of protection that the C3M offers to a cargo that it encapsulates. We have previously used Förster resonance energy transfer (FRET) to measure this molecular exchange of C3Ms and observed a broad range of exchange rates,22 similar to what has been observed earlier for other C3Ms.23 We hypothesised that this large difference in exchange rates is the result of chain polydispersity. Recently, the molecular exchange of C3Ms was also measured by small angle neutron scattering (SANS) and also here the presence of different exchange rates was explained by chain polydispersity.24 However, this polydispersity hypothesis has not been confirmed yet, since so far only the exchange of polydisperse polymers has been measured. Measuring the exchange of C3Ms with ssDNA allows to thoroughly test this hypothesis because the DNA is monodisperse and its length can be systematically varied.
Apart from the limited knowledge on the chain lengths effect on the dynamics in both complex coacervates and C3Ms, also little is known on how the dynamics inside bulk coacervates compares to that in the nanoscopically confined interior of C3Ms, and to what extent the dynamics in bulk coacervates is predictive for the exchange dynamics in micelles. Therefore, we study in this paper the dynamics of single stranded DNA in both complex coacervate droplets and micelles (Fig. 1) and we also compare these systems. For the complex coacervates droplets, we use fluorescence recovery after photobleaching (FRAP) to measure the diffusion coefficient of fluorescently labelled ssDNA. For the C3Ms, we use FRET-based exchange measurements to follow the exchange of both the ssDNA and the diblock copolymers and we complement these measurements with Langevin dynamics simulations. In this way, we show that both the DNA diffusion in complex coacervates and the DNA exchange of C3Ms depend on both the DNA chain length and the cationic chain length. We observe that the DNA diffusion coefficient shows a stronger dependence on the DNA length than predicted by the sticky Rouse model, which can only be partially explained by the chain length effect on coacervate density. In addition, we conclude from the comparison of the exchange of the monodisperse ssDNA and the polydisperse diblock copolymer that chain length polydispersity is indeed the main factor underlying the earlier observed broad distribution of exchange times. We hypothesise that the different exchange rates that we observe for the monodisperse DNA are mainly caused by differences in interactions of the donor and acceptor fluorophore label and show that the fluorophore can indeed have a large effect on the C3M exchange. Finally, we discuss a new description of the C3M exchange to account among others for the effect of the chain length of the oppositely charged core species on the exchange rate. Together our results can help to better understand both DNA specific systems and complex coacervates droplets and micelles in general.
The poly(L-lysine) homopolymers with number averaged lengths of 11, 55 and 234 amino acids (pLL11, pLL55 and pLL234 respectively) and the poly(ethyleneglycol)(5k)-b-poly(L-lysine)47 (PEG-5k-pLL47, where 5k indicates the 5300 g mol−1 molecular weight of the poly(ethyleneglycol) block and 47 indicates to the number averaged amino acid length of the poly(L-lysine) block) and poly(ethyleneglycol)(20k)-b-poly(L-lysine)212 (PEG-20k-pLL212) diblock copolymers were purchased from Alamanda Polymers as chloride salts. The overall polydispersities of the diblock copolymers as measured by the manufacturer were 1.06 for PEG-5k-pLL47 and 1.03 for PEG-20k-pLL212. Stock solutions of pLL homopolymers and diblock copolymers were prepared by dissolving in 1× phosphate-buffered saline (137 mM NaCl, 2.7 mM KCl, 8 mM Na2HPO4, 2 mM KH2PO4, pH 7.4) (PBS) and subsequent sonication for 10 minutes, as recommended by the manufacturer.
Sulfo-cyanine 3 NHS ester (sCy3) and sulfo-cyanine 5 NHS ester (sCy5) dyes for PEG-5k-pLL47 labelling were purchased from Lumiprobe.
Image analysis of the FRAP data was based on the method described by Taylor et al.25 The fluorescence intensity concentration C inside the bleach spot was normalised according to:
![]() | (1) |
To obtain the characteristic diffusion time τD, the bleaching corrected data was fitted with a 3D diffusion model in an infinite medium:25
![]() | (2) |
After image acquisition SPCImage 8.4 software (Becker & Hickl) was used to fit the fluorescence decay curves in each pixel. A two-component exponential decay was used to determine the average fluorescence lifetime in each pixel. A false colour scale was used to report lifetime values and yield the fluorescence lifetime images.
The majority of the sulfo-BODIPY partitioned inside the complex coacervate. Yet, a small fraction of the sulfo-BODIPY remained free in solution. To exclude the fluorescence lifetime of these free sulfo-BODIPY molecular rotors from the lifetime comparison, the lifetimes of pixels with a fluorescence intensity 2.5 times smaller than the average image intensity were excluded from the lifetime histograms. No atto488 fluorescence was detected outside the complex coacervate droplets and therefore the fluorescence lifetime histograms of the atto488 dye could be directly obtained from the FLIM images.
For the micelles with fluorescently labelled PEG-5k-pLL47 the same preparation protocol was used. The only difference is that labelled and unlabelled PEG-5k-pLL47 diblock copolymers were first mixed in the desired ratio and that afterwards the oligonucleotides were added. In this case, 10% of the PEG-5k-pLL47 polymers was coming from the fluorescently labelled PEG-5k-pLL47 stock solution and the other 90% was coming from the unlabelled PEG-5k-pLL47 stock solution with the same monomer concentration. The exact label fraction corresponding to this mix ratio is not known because we do not know the coupling efficiency of the PEG–pLL labelling reaction, but it is lower than 0.10 since the maximum average amount of fluorescent labels per polymer chain is 1 for this reaction.
The ssDNA/PEG–pLL micelles in 0.2× PBS solution were also prepared following the same thermal annealing protocol. For these micelles, only unlabelled ssDNA and PEG–pLL were added and no fluorescently labelled polymers were used.
![]() | (3) |
![]() | (4) |
![]() | (5) |
![]() | (6) |
An excitation wavelength of 490 nm was used to excite the (sulfo-)cyanine 3 donor dye. At this wavelength, the direct excitation of the (sulfo-)cyanine 5 acceptor dye is negligible. An excitation wavelength of 470 nm was used to excite the atto488 dye. To calculate the FRET efficiency, the emission of the atto532 acceptor dye was corrected for direct acceptor excitation at this wavelength. The degree of self-quenching of the atto532 acceptor dyes in the micelles was determined by using an excitation wavelength of 510 nm.
For the measurements of the equilibrated micelles, a single emission spectrum of the equilibrated sample was recorded. For the micelle exchange measurements, the equilibrated donor and acceptor micelles solutions were added to the cuvette in a 1 to 1 concentration ratio and mixed and placed in the spectrophotometer. Every minute an emission spectrum was recorded.
To determine the FRET efficiency, the spectrum (after correction for direct acceptor excitation for the atto488/atto532 FRET pair) was fitted with a linear combination of fixed log normal functions to determine the relative contribution of the donor and acceptor emission to the overall emission spectrum (ESI,† Section S3). Finally, the FRET efficiency E was calculated by:
![]() | (7) |
![]() | (8) |
Subsequently, the data of the C3M exchange experiments is fitted to the following function:
![]() | (9) |
![]() | (10) |
The average number of donors in acceptor micelles at a certain time t is given by:
![]() | (11) |
![]() | (12) |
![]() | (13) |
![]() | (14) |
![]() | (15) |
![]() | (16) |
![]() | (17) |
In this case, we have used two different exchange rates for both the donor and acceptor exchange. We imposed that the fraction of fast exchanging fluorophores xfast was the same for the donors and acceptors. In addition, we required that the fast exchange rate kfast was the same for the donor and the acceptor, while the slow exchange rate of the donor and acceptor could be different. This resulted in four different fit parameters: xfast, kfast, kD and kA, where kD and kA indicate the slow exchange rate of the donor and acceptor respectively. The fits were performed by using a nonlinear least-square fit algorithm, as implemented in the python software package SciPy, with the following boundaries for the fit parameters: xfast = [0,0.2], kfast = [10 h−1, 1000 h−1], kD = [0 h−1, 100 h−1] and kA = [0 h−1, 100 h−1]. For the initial guess of kD we have used a value twice as large as the initial guess of kA (in all cases, kD = 0.02 h−1 and kA = 0.01 h−1) since the use of kD = kA as initial guess decreased the quality of the fit with E(t) = E(∞) being reached too early.
For all 223 diblock copolymers in the simulation a neutral block length of 50 monomers was used. The cationic block was polydisperse and 48 diblock copolymers had a cationic block length of Npos = 14, for 47 diblock copolymers Npos = 16, for 40 diblock copolymers Npos = 18, for 40 diblock copolymers Npos = 20 and for 48 diblock copolymers Npos = 22. The polyanion was monodisperse. Two different simulations were performed: one with 400 polyanions with a length of Nneg = 10 and one with 200 polyanions with Nneg = 20. In both cases, the periodic box size was Ls = 163σs.
We started the simulations by placing the homopolymers, diblocks, and counterions randomly in the simulation box and used the GPU-optimized molecular dynamics software package HOOMD-blue32–35 to perform the Langevin dynamic simulations. The Coulomb interactions were calculated by using the particle–particle particle-mesh (PPPM) Ewald summation method34 with a real space cut off distance of 2.5σs. The neighbour lists were generated by using the linear bounding volumes hierarchies (LBVHs) method.35 We used γ = 1ms/τs as drag coefficient and Δt = 0.005τs as simulation time step where τs is the time unit in the system. We saved the configuration of the simulation every 2500 steps. The total simulation time was 2 × 107τs.
To determine the number of split events per polymer chain, the composition of the micelle clusters was followed over time. The micelle cluster identification was performed in the same way as described earlier.31 Every time a polymer was part of the smallest cluster formed after a splitting event of a micelle cluster, this was counted as a split event for this polymer. To reduce the effect of micelle growth on the split events analysis, the split events at a simulation time shorter than 2 × 106τs were excluded from the analysis.
Here we have used fluorescence recovery after photo-bleaching (FRAP) to systematically study the chain length effect on the diffusion of fluorescently labelled ssDNA in ssDNA/pLL coacervates (Fig. 2). We have varied both the pLL length and the ssDNA length independently from each other. The increase in FRAP recovery time for increasing chain lengths, shows that the DNA diffusion is not only affected by its own length (Fig. 2b), but also by the length of the oppositely charged poly(L-lysine) (Fig. 2c). This indicates that the relaxation time of a polyelectrolyte in a complex coacervate is not only determined by its own length and the ionic strength of the solution, as shown earlier,18,19,38 but that also other complex coacervate properties can affect the dynamics.
From the recovery curves the diffusion coefficient of the DNA in the complex coacervate droplets can be determined. Here, we have used the 3D diffusion model developed by Taylor et al.25 to obtain the apparent DNA diffusion coefficients, as explained in the materials and methods (Section 2.3). The DNA diffusion coefficient strongly decreases with increasing DNA length N (Fig. 2d). In fact, it decreases more strongly than the D ∝ N−1 dependence predicted by the sticky Rouse model, which is the model that is often used to describe polyelectrolyte dynamics in coacervates. Instead of an exponent of −1, we find an exponent of −2.6 ± 0.8, −2.9 ± 0.4 and −2.5 ± 0.1 for pLL11, pLL55 and pLL234 respectively. A possible explanation for this stronger decrease is that the density of the complex coacervate also increases with increasing DNA length and that this increase in density slows down the dynamics in the complex coacervate. This increase in density is predicted by multiple coacervation models39,40 and has been mentioned as explanation why the viscosity depends more strongly on the polyelectrolyte length than predicted by the sticky Rouse model.18 This density increase has also been mentioned to explain why the relaxation time in coarse-grained simulations of salt-free polyelectrolyte complexes of equal length showed a stronger dependence on the polyelectrolyte length than the sticky Rouse model.41 In the latter case, the increase in complex density was also confirmed by the same simulations. For the ssDNA/pLL complex coacervate droplets, an increase in density with increasing chain length might not only explain the stronger dependence of DNA diffusion coefficient on the DNA length, but it might also explain why the pLL chain length affects the DNA diffusion coefficient in the complex.
Increasing the DNA chain length results in an increase of fluorescence lifetime of the sulfo-BODIPY molecular rotor (Fig. 3), especially for the longest ssDNA length (nt88). This increase in lifetime indicates that the intramolecular rotation of the molecular rotor becomes more hindered and that the local complex coacervate density thus increases for longer ssDNA lengths. This shows that the deviations from the expected D ∝ N−1 dependence might indeed (partly) be explained by changes in the complex coacervate density.
In contrast, changes of the pLL chain length do not affect the lifetime of the rotor (Fig. 3m, n and o), even though the changes in the pLL lengths are larger than the changes in ssDNA length. This is unexpected because in theoretical descriptions of complex coacervation,39,40 longer chain lengths increase both the complex coacervate density and the critical salt concentration of the complex coacervates. This critical salt concentration is the salt concentration at which the complex coacervates fall apart and where the system transitions from a two phase system (below the critical salt concentration), consisting of a polymer rich phase (the coacervate) and a diluted phase, to a one phase system (above the critical salt concentration). Earlier measurements have shown that the critical salt concentration of ssDNA/pLL coacervates increases both by increasing the ssDNA length or by increasing the pLL length.42 Still, we do not observe any change in fluorescence lifetime of the molecular rotor in coacervate complexes with increasing pLL length. This might mean that the pLL length has no effect at all on the complex coacervate density or that the effect of the pLL length is too small to be observed by using the sulfo-BODIPY molecular rotor. Unfortunately, the precise sensitivity of the molecular rotor is difficult to determine from calibrations with solutions with a known viscosity because the fluorescence lifetime also depends on the solvent polarity, but to give a first indication of this sensitivity: in homogeneous solutions, the fluorescence lifetime of molecular rotors show a powerlaw dependence on the macroscopic viscosity, with an exponent that can differ per molecular rotor type.43
Irrespective of whether no density changes occur with changing the pLL length or that the density changes are too small to measure, the effect of pLL chain length is different from ssDNA length effect: while we do not observe any effect for the pLL length, we can observe an increase in the molecular rotor fluorescence lifetime increase with increasing DNA lengths. Similar asymmetric length effects on the complex coacervate density have been observed earlier in rheology measurements, where changing the length of one of the polyelectrolytes left the storage and loss modulus nearly unaffected, while increasing the length of the oppositely charged polyelectrolyte increased the both the storage and loss modulus.18 A higher modulus reflects a larger crosslink density and thus an increase in coacervate density.
To verify whether only minor changes in complex coacervate density occur when the pLL length is changed, we have performed FRET-FLIM measurements of nt44/pLL complex coacervates containing both nt44 labelled with a donor fluorophore and nt44 labelled with an acceptor fluorophore. In FRET-FLIM measurements, the lifetime of the donor fluorophore is measured, which will become shorter when the energy transfer from the donor to the acceptor is faster. This energy transfer becomes faster when the donor and acceptor are closer to each other. An increase in complex coacervate density will decrease the average distance between the donors and acceptors and therefore will result in a decrease of the average donor fluorescence lifetime. Comparison of the fluorescence lifetime histograms for the different pLL chain lengths (Fig. 4) shows that the pLL length only has little effect on the fluorescence lifetime of the donor, suggesting that the changes in complex coacervate density are indeed small. In fact, we can use the changes in donor fluorescence lifetime to roughly estimate the density changes: the average measured fluorescence lifetime of the donor 〈τF〉 is given by 〈τF〉 = 1/(kA + kD) where kA is the energy transfer rate to the acceptor and kD is the rate of photon emission from the donor in absence of the acceptor. The ratio between kD and kA depends on the distance r between the donor and acceptor fluorophore: kD/kA = (r/RF)6 where RF is the Förster radius. The ratio of transfer rate to the acceptor kA can thus be used to calculate the changes in average distance between the donor and acceptor fluorophore: r2/r1 = (kA,1/kA,2)1/6. Since the coacervate density ϕ ∝ r−3, we get ϕ2/ϕ1 = (kA,1/kA,2)−0.5. With an average donor fluorescence lifetime in complex coacervates without acceptor of 〈τF〉donor = 2.1 ns (and therefore kD = 1/〈τF〉donor ≈ 0.5 ns−1), an average donor fluorescence lifetime in a nt44/pLL11 complex coacervate of 〈τF〉pLL11 = 1.4 ns and an average donor fluorescence lifetime in a nt44/pLL234 complex coacervate of 〈τF〉pLL234 = 1.2 ns, we get ϕpLL234/ϕpLL11 ≈ 1.2. This factor is comparable to the density changes predicted in theoretical models for complex coacervation39,40 and might indicate that still minor changes in complex coacervate density occur with increasing pLL length, but that they are too small to be observed with the sulfo-BODIPY molecular rotor.
The effect of the pLL chain length on the complex coacervate density is relatively small compared to its effect on the DNA diffusion coefficient (Fig. 2): for nt44, the apparent diffusion coefficient decreased by a factor ∼3.4 when changing from pLL11 to pLL234, while the estimated density increase is only a factor ∼1.2. If changes in density are the only cause for the change in DNA diffusion coefficient with increasing pLL length, this would mean that the dynamics has to depend strongly on the complex coacervate density. In the classical model for the dynamics of associating polymers,21 there are two scenarios were the dynamics depends strongly on the concentration. The first scenario occurs in the sticky Rouse regime and corresponds to a transition from intramolecular to intermolecular bonds with increasing concentration. The second scenario occurs in the reptation regime, where an increase in concentration decreases the reptation tube diameter. Both scenarios do not apply to the ssDNA/pLL complex coacervates studied here. First, the reptation regime is only relevant when there are significant entanglements, which only occurs for long polyelectrolyte lengths.41 Here, we also see a strong length dependence for the short pLL and DNA lengths, which are too short for the reptation regime. In addition, the ionic bonds in (homo)polyelectrolyte complexes are always intermolecular since they are based on the attraction between oppositely charged species and therefore no intramolecular to intermolecular bond transition will occur. The probability of ionic bond formation can also increase with increasing concentration, but this probability p scales only weakly with the density (p ∝ ϕ for a theta solvent and p ∝ ϕ1.59 for a good solvent)21 and is therefore not enough to explain the relatively large change in DNA diffusion coefficient for small concentration changes. Our results thus seem to suggest that pLL length also affects the dynamics in a different way than by only changing the density. These other effects have not been discussed earlier for strongly charged polyelectrolyte complexes, but in a scaling model for weakly charged complexes a more complex length dependence has already been predicted for certain cases.44 Yet, we can also not exclude that the observed large difference between the density change and diffusion coefficient change is caused by experimental uncertainty. For example, there is a quite large variation in the donor fluorescence lifetime measured in the same complex coacervate (Fig. 4) and therefore the real density variation might be different from the estimated factor of 1.2. Therefore, it would be interesting to see in future research whether these small changes in density and relatively large changes in dynamics can also observed with other techniques, for example by using small angle X-ray scattering (SAXS) or small-angle neutron scattering (SANS) experiments to measure the correlation length45 or mesh size46 of the complex coacervate or rheology measurements to probe both the complex coacervate dynamics and its concentration.18,19
In conclusion, the chain length can have a complex effect on the density and dynamics of ssDNA/pLL complex coacervates. Changing the ssDNA chain length has a stronger effect on the coacervate density than changing the pLL length. The changes in complex coacervate density might (partly) explain the effect of the chain length on the dynamics in the complex coacervates. Yet, we cannot exclude that also other factors play a role, since the estimated changes in complex coacervate density seem relatively small compared to the observed changes in ssDNA diffusion coefficient.
Both the ssDNA length and the PEG–pLL diblock copolymer length seem to affect the molecular exchange rate of DNA between C3Ms (Fig. 5a and b): for shorter chain lengths, the observed increase in FRET efficiency is faster, which indicates a faster molecular exchange rate. To better observe the differences in exchange rate for the fast exchanging C3Ms, we have also performed the exchange measurements in a 0.2× PBS solution (Fig. 5c and d) instead of a 1× PBS solution. The lower ionic strength in the 0.2× PBS solution increases the strength of the electrostatic interactions in the C3M core and therefore decreases the exchange rate.
![]() | ||
Fig. 5 FRET-based exchange experiments of C3Ms with fluorescently labelled monodisperse DNA for different diblock copolymer lengths. In all cases, the donor fluorophore is atto488 and the acceptor fluorophore is atto532. The ssDNA length is 22 (a and c) or 44 (b and d) and the micelles are dissolved in a 1× PBS solution (a and b) or a 0.2× PBS solution (c and d). Shaded regions are an indication of the uncertainty in the exchange measurements with the region borders given by the minimum and maximum E(t)/E(∞) of three repetitions of the same exchange experiment. Symbols indicate an example of one of these three repetitions and solid black lines in (c) and (d) indicate fits of the data of this example exchange experiment to the analytical FRET model. The corresponding input parameters and the fit results are given in Table 1. |
We observe similar trends for the chain length effects on the exchange rate in both 1× PBS and 0.2× PBS solution, but we also observe very large variations between repetitions of the same experiment, as indicated by the shaded regions in Fig. 5. Because of this large uncertainty, we can conclude only that most likely both the ssDNA length and the PEG–pLL diblock copolymer length have a substantial effect on the DNA exchange rate, but we cannot exactly quantify this effect. The large variations between repetitions were unexpected: the formation of these ssDNA/PEG–pLL C3Ms is well-studied and the use of the thermal annealing protocol makes the formation of unlabelled ssDNA/PEG–pLL micelles well reproducible.13 Still, we expect that the origin of these large variations lies in differences in micelle preparation since the variations are much smaller for repetitions of the exchange experiments from the same micelle stock solutions (Fig. S5–S7, ESI†). These differences between micelle preparations might point to non-equilibrium micelle formation. The aim of the thermal annealing protocol was to prevent this non-equilibrium micelle formation, but maybe this protocol works insufficiently in the presence of dyes. Possibly, alternative annealing protocols, like salt annealing, could decrease the variation between the micelle preparations. Another explanation for the large variations between experiments is that differences in dye concentrations result in different exchange rates: when we compared the repetitions of the same mix experiment, we observed differences in dye concentration, in particular in the acceptor concentration (Fig. S5–S7, ESI†), which might be caused by concentration inhomogeneities in the fluorescently labelled DNA stock solutions. We will come back to this possible relation between uncertainty in the molecular exchange rate and dye concentration in the next section.
In some cases, the effect of the diblock copolymer length on the DNA exchange rate seems much larger than the effect of the polycation length on the DNA diffusion in complex coacervates: for ssDNA with a length of 44 nucleotides the measured exchange rates differ two orders of magnitude or more between PEG-5k-pLL47 and PEG-20k-pLL212 (Fig. 5b and d), while the average nt44 diffusion coefficient for pLL55 and pLL234 only differs a factor ∼1.4 (Fig. 2). This could suggest that the DNA exchange rate is not governed only by relaxation processes in the core since in that case the poly(L-lysine) length should have a similar effect on the DNA exchange rate and the DNA diffusion. Yet, for the ssDNA with a length of 22 nucleotides, the measured diblock copolymer effect is less pronounced and the variations of the exchange measurements for the different diblock copolymer are closer to each other. Therefore, to be entirely sure that the pLL chain length indeed has a much stronger effect on the molecular exchange rate than on the DNA diffusion, the uncertainty in the exchange rate determination first has to be decreased.
The measured DNA exchange rate is relatively fast compared to earlier measurements of the molecular exchange of C3Ms: the nt22/PEG–pLL47 micelles exchanged within 2 minutes, which is much faster than the exchange rate that we measured for C3Ms with poly(3-sulfopropyl methacrylate) (PSPMA) and poly(ethylene glycol)–poly(2-trimethylammonioethyl methacrylate) (PEG–PTMAEMA) at similar salt concentrations,22 where after 40 h the micelles were still not completely exchanged. The measured exchange rate of the nt22/PEG–pLL47 micelles is even faster than exchange rate of C3Ms with fluorescent proteins,47 while the protein exchange measurements were performed very close to the critical salt concentration at which the C3Ms disassemble and the nt22/PEG–pLL47 micelles are performed far away from the critical salt concentration (their critical salt concentration is ∼0.6 M NaCl48 and 1× PBS contains 0.137 M NaCl). In fact, the critical salt concentration is more comparable to the PSPMA/PEG–PTMAEMA micelles (∼0.79 M KCl), while the exchange rate is much faster than for the PSPMA/PEG–PTMAEMA micelles. Also for differences between the ssDNA/PEG–pLL micelles themselves, the exchange rate does not correlate with the critical salt concentration: an increase in diblock length seems to decrease the molecular exchange rates (Fig. 5), while it seems to leave the critical salt concentration unaffected.48 Together these observations show that the critical salt concentration is not a good predictor for the molecular exchange rate of the C3Ms.
Not only the exchange rates of the ssDNA/PEG–pLL micelles are different from the PSPMA/PEG–PTMAEMA micelles, but also the shape of exchange curve: the PSPMA/PEG–PTMAEMA showed a slow, gradual increase of the FRET efficiency as function of time plotted on logarithmic axes, while for the DNA exchange the increase is much steeper. A gradual increase on a double logarithmic scale corresponds to a broad distribution of exchange times. For the PSPMA/PEG–PTMAEMA measurements, we hypothesised that the broad range of exchange times was the result of chain polydispersity. The steeper increase for the monodisperse ssDNA exchange could support this hypothesis. Yet, we still observe different time scales for the C3Ms with slow monodisperse ssDNA exchange: there is an initial increase that is too fast to be measured and afterwards a second exchange time scale seems to occur. Subsequently, we observe a quasi-plateau and then the FRET efficiency increases again at longer logarithmic time scales. The initial fast exchange might correspond to a small population of unstable micelles, for example because these micelles are destabilised by the shear of mixing of donor and acceptor micelle solutions. A possible explanation for the two other time scales is that the ssDNA labelled with donor fluorophores have a different exchange rate than the ssDNA labelled with acceptor fluorophores.
As a first check of this explanation for the different time scales in the monodisperse ssDNA exchange, we have fitted the C3M exchange experiments with an analytical model for FRET-based micelle exchange experiments, which we developed recently.22 In this case, both the donor and acceptor fluorophores consisted of two different species, where the exchange rate of the fast exchanging species was the same for the donor and acceptor, while the exchange rate of the slow exchanging species could be different for the donor and acceptor. Further details of this fit procedure given in the materials and methods (Section 2.10). We note that the purpose of this fit is to explain the shape of the DNA exchange experiments curves and not to exactly determine the rate constants because the experimental uncertainty is too large for this. The precise values of the rate constants obtained from the fit (Table 1) are therefore not that relevant. The fits can describe the exchange experiments well (Fig. 5c and d), which suggests that the ssDNA exchange indeed occurs with only three different exchange rates instead of a broad distribution.
Sample | N | α | ν | x fast | k fast (h−1) | k D (h−1) | k A (h−1) |
---|---|---|---|---|---|---|---|
nt22/PEG-5k-pLL47 | 2454 | 0.03 | 0.01 | — | — | — | — |
nt22/PEG-20k-pLL212 | 3835 | 0.03 | 0.01 | 0.2 | 62 | 9 | 0.2 |
nt44/PEG-5k-pLL47 | 1079 | 0.06 | 0.01 | 0.2 | 133 | 41 | 1 |
nt44/PEG-20k-pLL212 | 1702 | 0.06 | 0.01 | 0.07 | 30 | 0.6 | 0.001 |
To further study the dye effect, we have also increased the dye concentration by a factor 2. A larger dye concentration results in a slower exchange (Fig. 6b). This shows that variations in dye concentration can indeed affect the exchange rate, as we suggested in the previous section to explain the large variations in exchange rate between repetitions of the same measurements. However, for the repetitions of the same exchange experiments the dye concentration effect seems to be reversed and the micelles with a lower acceptor concentration seem to exchange more slowly (Fig. S5–S7, ESI†). This could point to a complex dye effect where for example the acceptor concentration has a different effect on the exchange rate than the donor concentration. Alternatively, except from the dyes, also other effects in the micelle preparation might have played a role, like the formation of non-equilibrated micelles. To test whether other factors than the dyes play a role in the large variations in exchange experiments, the ssDNA/PEG–pLL exchange could be measured also by using small angle neutron scattering (SANS) since for these measurements no dyes are needed.
The large effects of the dyes on the molecular exchange rates of C3Ms show that one has to be careful with interpreting the results of FRET-based micelle exchange experiments and that in general SANS measurements are a better option to measure the exchange: SANS measurements are based on neutron scattering contrast and the introduction of deuterated protons will have a much smaller effect on the exchange rate than the introduction of dyes. Although the dye effect is a large disadvantage for fundamental studies to the molecular exchange of C3Ms, it shows that the molecular exchange can be tuned by just introducing an additional hydrophobic group. It was already shown earlier that additional hydrophobicity can affect the static properties of both complex coacervates50 and complex coacervate core micelles.51 In addition, by simulating the C3M exchange with Langevin dynamics simulations, we showed that the non-electrostatic attraction (like hydrophobic interactions) between the oppositely charged polyelectrolytes, can largely affect the exchange rate.31 These simulations also suggested that the non-electrostatic interactions between only one of the core species might not be enough to reduce the exchange rate. In that respect, it would be interesting to see whether the dyes have similar effects on the molecular exchange of other C3Ms systems or that also specific interactions with the poly(L-lysine) and/or the ssDNA have played a role in the large dye effect that we observed here.
We indeed observe a gradual increase in FRET efficiency for the C3Ms with fluorescently labelled PEG–pLL diblock copolymers, similar to what we observed earlier for the PSPMA exchange of PSPMA/PEG–PTMAEMA micelles.22 This further supports that chain polydispersity is the main cause for a broad distribution of exchange rates observed within a single exchange experiment. The difference between the ssDNA and PEG–pLL exchange also shows that the exchange of one of the core species cannot be used as a direct measure for the exchange of the oppositely charged core species.
The differences between the shape in the monodisperse ssDNA and polydisperse diblock copolymer exchange curves can be explained in two ways. In the first explanation, the ssDNA and the diblock copolymer exchange independently from each other. This explanation is unlikely because this does not agree with the combined exchange that we observed earlier in Langevin dynamics simulations.31 In addition, the diblock copolymer length seems to have a substantial effect on the ssDNA exchange rate (Fig. 5), which we do not expect for independent exchange. In the second explanation, the ssDNA and diblock copolymer exchange together, but the exchange of neutrally charged complexes is preferred and therefore the ssDNA exchanges mainly with polydiblock copolymers with a polycationic block length that is the same as (a multiple of) the ssDNA length. These diblock copolymers therefore have the largest exchange rate. Diblock copolymers with polycationic block lengths that result in the formation of complexes with a net charge will exchange less often and therefore have a lower exchange rate. In this way, the average exchange rate of all monodisperse DNA molecules can be the same, while the exchange rates of the polydisperse PEG–pLL diblock copolymers can differ depending on their polycationic length.
The fact that the charge ratio between the cationic block and the polyanion(s) has a larger impact on the diblock copolymer exchange rate than the cationic block length itself (Fig. 8a and b) indicates that the exchange rate is not only governed by the polyelectrolyte movements in the C3M core: inside a complex coacervate, the oppositely charged polyelectrolytes do not move simultaneously, as shown by different diffusion coefficients of oppositely charged polyelectrolytes in the same complex,52 and therefore for movements inside the C3M core, the absolute length of the polyelectrolytes should be more important than the charge ratio. Instead, the large dependence on the charge ratio suggests that also an activated process is involved in the exchange where the height of the energy barrier is mainly determined by the number of uncompensated charges in the expelled complex.
In contrast, the strong dependence of the exchange rate on the polyanion length, shows that the exchange rate is also not only governed by an activated process. For example, the same number of monomers is involved in the exchange of a polycationic block with Npos = 20 with the oppositely charged polyanion(s) with length Nneg = 10 or Nneg = 20. Therefore, the energy barrier is the same for both polyanion lengths. Still, we observe a much larger exchange rate for the shorter polyanion (Fig. 8a and b). Similar effects occur for other cationic block lengths and we have also observed the same effects in exchange experiments with fluorescently labelled ssDNA (Fig. 5). This indicates that relaxation processes in the C3M core also affect the exchange rate.
The chain length effect on the core relaxation part of the C3M exchange cannot be described by a simple (sticky) Rouse relaxation: for the polyanion, the average number of split events with a cationic block with length Npos = 20 decreases by a factor ∼8.8 when increasing the polyanion length from Nneg = 10 to Nneg = 20, which is larger than the τ ∝ N2 dependence of the sticky Rouse model. In contrast, increasing the cationic block length from Npos = 18 to Npos = 22, decreases the number of split events only by a factor of ∼1.03 (Nneg = 10) or ∼1.3 (Nneg = 20), while in the sticky Rouse model this should be a factor 1.5 and this factor can even become larger when also the effect of non-electrostatic interactions on the energy barrier is taken into account. At the moment, we cannot exactly explain this complex effect of the chain length on the core relaxation part of the exchange. Possibly, the average number of chains to which a polyelectrolyte binds partly determines the exchange rate. Shorter chains have a smaller radius of gyration and therefore bind to less chains on average, which might increase the expulsion rate. In that case, density effects might have further enhanced the effect of the polyanion length since a lower density could further decrease the average number of chains to which a polyelectrolyte binds. This density effect did not occur in the comparison of the cationic block length effect because here we compared cationic block lengths that are part of the same C3Ms and therefore the density is the same for both polycation lengths.
In summary, both an activated process and relaxation processes in the core determine the C3M exchange rate (Fig. 9), as was also suggested for the exchange of amphiphilic diblock copolymers micelles.53 However, the underlying mechanism of the C3M exchange are different than for the amphiphilic diblock copolymer micelles: for the C3M exchange, the exchange rate of a core species depends also on the length of the oppositely charged core species, while the exchange of an amphiphilic diblock copolymer only depends on its own length because there is only one core species. As a result, the energy barrier of the activated process in C3M exchange might not always increase with increasing polyelectrolyte length and the charge ratio can be more important than the absolute polyelectrolyte length. In addition, the contribution of the relaxation in the core to the exchange rate deviates from the Rouse relaxation suggested for amphiphilic diblock copolymer micelles. These differences between C3Ms and amphiphilic diblock copolymer micelles can explain why we could not use the amphiphilic diblock copolymer model to describe the exchange of PSPMA/PEG–PTMAEMA micelles22 and why in a recent SANS study of C3M exchange a lower polydispersity was needed than measured to obtain a good fit between the exchange data and the amphiphilic diblock copolymer model.24 Furthermore, these differences stress the importance of specifically studying C3Ms to further unravel the effect of chain lengths and other parameters on the C3M exchange rate.
For complex coacervate core micelles, the chain length effects were more difficult to quantify because we observed large variations in C3M exchange rates between repetitions of the same experiment. We hypothesise that these large variations are caused by dye concentration variations and we have shown that changing the dye concentration can indeed affect the C3M exchange rate. In addition, we have shown that the dye type can also affect the exchange rate and we have hypothesised that the different exchange rates observed for the exchange of monodisperse DNA are mainly caused by a differences in exchange rate for DNA labelled with a donor fluorophore and DNA labelled with an acceptor fluorophore. Measuring the exchange of ssDNA/PEG–pLL micelles with small angle neutron scattering (SANS) could be used to test our hypotheses on the dye effects.
Although the determination of the exact values of the C3M exchange rate constants was hindered by the large uncertainty of the measurements, the shape of the FRET-based C3M exchange curve still provided information on the exchange mechanisms: the broader distribution of exchange rates for the polydisperse PEG–pLL diblock copolymer than the monodisperse ssDNA suggests that chain length polydispersity is the main cause for a broad distribution of exchange times. The different exchange rates for ssDNA and PEG–pLL can be explained by a preferred exchange of ssDNA with PEG–pLL with a cationic block length that is the same as (a multiple of) the ssDNA length, which results in a lower exchange rate for PEG–pLL polymers with mismatched cationic block lengths. We have supported this explanation by performing Langevin dynamics simulations and we have shown that both relaxation processes in the core and an activated process determine the C3M exchange rate. Furthermore, we have concluded that the chain length effect on the C3M exchange is different from the chain length effect on the exchange of amphiphilic diblock copolymers. For C3Ms, the expulsion of neutral complexes is preferred and as a result the length ratio of the oppositely charged core species might play a larger role than the absolute chain length. Therefore, the length ratio of the oppositely charged polyelectrolytes might be used to tune exchange rate. However, the presence of polydisperse polymers make it more difficult to use this length ratio as tuning parameter since a mismatch in only the average chain length is not enough to prevent fast exchange, as was illustrated here by the fast ssDNA exchange in the presence of polydisperse PEG–pLL. Apart from the length ratio with the oppositely charged cationic block, the ssDNA length itself can also be used to adapt the exchange rate, with a lower exchange rate for longer ssDNA lengths.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1sm01787j |
This journal is © The Royal Society of Chemistry 2022 |