P.
Parvizi-Bahktar
a,
J.
Mendez-Campos
a,
L.
Raju
b,
N. A.
Khalique
b,
E.
Jubeli
bc,
H.
Larsen
d,
D.
Nicholson
e,
M. D.
Pungente
f and
T. M.
Fyles
*a
aDepartment of Chemistry, University of Victoria, P.O. Box 3065, Stn CSC, Victoria, BC V8W 3V6, Canada. E-mail: tmf@uvic.ca
bResearch Division, Weill Cornell Medicine in Qatar, Education City, P.O. Box 24144, Doha, Qatar
cUniversité Paris-Sud, EA 401, IFR 141, Faculté de Pharmacie, F-92296 Châtenay Malabry, France
dDepartment of Physics, University of Stavanger, 4036 Stavanger, Norway
eDepartment of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
fPremedical Unit, Weill Cornell Medicine in Qatar, Education City, P.O. Box 24144, Doha, Qatar
First published on 15th February 2016
The efficiency of the transfection of a plasmid DNA encoding a galactosidase promoted by a series of pyridinium lipids in mixtures with other cationic lipids and neutral lipids was assessed in CHO-K1 cells. We identify key molecular parameters of the lipids in the mixture – clogP, lipid length, partial molar volume – to predict the morphology of the lipid-DNA lipoplex and then correlate these same parameters with transfection efficiency in an in vitro assay. We define a Transfection Index that provides a linear correlation with normalized transfection efficiency over a series of 90 different lipoplex compositions. We also explore the influence of the same set of molecular parameters on the cytotoxicity of the formulations.
Non-viral nucleic acid delivery via lipid–nucleic acid complexes (lipoplexes) is significantly less effective compared to viral delivery but the low immunogenic potential and generally moderate toxicity of the active agents make it a very appealing alternative.5,9,10 Since the first report by Felgner et al.11 in 1987 describing the use of the cationic lipid N-[1-(2,3-dioleyloxy)propyl]-N,N,N-trimethylammonium (DOTMA) as a vehicle for nucleic acid delivery, the following decades brought modifications and refinements mainly to lipid architectures and particle formulations with the hope of vastly improving gene transfer efficiencies.9,10,12,13 The majority of structural modifications to the lipid architecture focused on features having a significant impact on particle integrity and packing morphology, specifically the lipid hydrophobic domain,14,15 the functional nature of the charged head group,16–18 and the type of chemical bond or linkage through which these groups are connected.19–22 Although a few general trends emerged such as the advantage of using a double-chain hydrophobic domain23 and the replacement of saturated with monounsaturated chains,13 it is difficult to correlate cationic lipid structural details with transfection efficacy using simple “golden rules” and the strategies to improve performance are still mainly empirical.
Typically the anionic nucleic acid is rendered hydrophobic via electrostatic complex formation with cationic amphiphiles or polymers which, in the presence of neutral co-lipids, form into particles, liposomes (vesicles) or other aggregates. The packaged nucleic acid is protected within a hydrophobic environment, the particles themselves are not initially recognized by the immune system, and the individual components are potentially of low toxicity24 and can be manufactured to high purity with conventional GMP.9 The functional supramolecular sophistication of these systems is not very evolved and there is still considerable uncertainty over which supramolecular organizations or properties hold the keys to highly effective transfection.9 Supramolecular structure and assembly follows from the molecular properties of the components so is potentially amenable to optimization through the structures of the components. Our focus here is on low molecular weight cationic lipids and neutral co-lipids. Our premise is that the molecular properties of these small molecules will control supramolecular organization and function and thus will control transfection efficiency.
Which molecular and supramolecular properties are known to play a role in transfection by cationic lipids? Extra-cellular stability is obviously important; the particle must resist both enzymatic and thermal degradation (aggregation, decomposition), and in vivo must resist opsonation (labelling for destruction). This is usually accomplished through a surface coating of lipid-anchored polyethylene glycol (“pegylation”) that provides a steric cushion that limits the interactions with proteins in the extra-cellular medium,25,26 although cationic lipid head group optimization can also produce serum-resistant formulations.27 Physical factors such as particle size and surface charge density are clearly important at the level of lipoplex transport and uptake in vivo.7,8,28–32 Cellular uptake of larger particles in vitro occurs readily.33–35 Cells are diverse and utilize parallel uptake pathways each having particular triggering events,28 but the basic events from the particle perspective are similar between pathways; the particle is recognized for uptake, is pulled into the cell, and is initially held within an endosome or similar compartment within the cell. These stages can handle particles up to 400–500 nm in diameter,28 although in vivo applications likely require particles smaller than 125 nm (ref. 36) to avoid filtering by the spleen. The initial recognition of the particle may depend on the particle lipid phase morphology as well.8,37
Lipid phase behaviour plays a key role in the release of the nucleic acid cargo from within the endosome.8,37–42 Release is clearly correlated with the propensity to form inverted lipid phases, either hexagonal or cubic. This propensity is proposed to facilitate migration of lipids to the endosomal membrane, to reorganize the lipids locally to an inverted phase leading to rupture of the lamellar bilayer that contains the particle and its cargo.40,42,43 This key phase behaviour is temperature-dependent, and transfection is optimal when the inverted phase transition is close to the temperature of transfection (typically 37 °C).41,42 Lipid phase morphology is directly related to lipid structure through the shape or packing parameter S = V/(a0 × lc); where V is the equilibrium hydrocarbon chain volume, a0 is the equilibrium head group area at the polar-non-polar interface, and lc is the equilibrium hydrocarbon chain length.44,45 Pure lipids with ∼0.7 < S < 1.0 form lamellar phases; S > 1.0 leads to inverted phases. Lipid mixtures behave as pure lipids and have a shape parameter Smix given as the molar weighted average of the S values of the components.46
A common structural variable in synthesis-based structure–activity studies of cationic lipids is the hydrocarbon chain (length, degree of unsaturation).10,13 This parameter will influence S to a limited extent but plays a more significant role in varying the hydrophobicity and partitioning behaviour of the lipid. Quantitative structure–activity explorations47,48 of transfection efficiency include partitioning using a calculated logP term (clog
P) based on additive structural increments. The rationale is that a more hydrophobic lipid will increase the structural stability of the initial particle, but at some point will inhibit cargo release leading to a bell-shaped dependence48 usually described as hydrophilic–lipophilic balance.49
We previously reported the development of a Transfection Index (TI) that was based upon S, clogP, and the geometric constraints of packing lipids and DNA into a lipoplex having overall hexagonal (HII) morphology.50 The goal was to develop a predictive tool based upon additive molecular structural parameters, ideally without recourse to any adjustable parameters. To the extent that we uncovered a linear correlation between TI and transfection efficiency by a limited number of previously reported pyridinium diesters51 in a range of formulations, we were encouraged that the structural parameters selected and the overall approach had some potential. However, we need to evaluate a broader set of compounds and formulations to establish if the previous correlation was adventitious.
We report here transfection, cell viability, and structural experiments for a larger set of pyridinium diesters, including binary and ternary mixtures of hydrocarbon chains. This dataset can be correlated using a modified TI that takes into account the potential role lipid phase behaviour plays in transfection. We also explore the correlation of formulation and structural parameters on cytotoxicity.
The yields of the intermediate diester pyridines are variable, primarily due to chromatographic losses. The diester pyridines are very weakly basic54 and difficult to directly N-alkylate. In our hands, alkylating agents that produce nucleophilic counterions, e.g. methyl iodide, also produce a small amount of addition product (N-alkyl-2-X-dihydropyridines). To circumvent this problem we use the forcing trimethyloxonium tetrafluoroborate for the methylation step which occurs quantitatively.
Below we explore binary blends of symmetrical diesters as a means to vary S and clogP. We were interested as well in mixed chain species as a comparison. A statistical synthesis with mixed starting alcohols e.g. oleyl alcohol and HO-brC20:0, produces a ternary mixture (diC18:1, dibrC20:0, and (C18:1)(brC20:0)). Electrospray MS, either of the protonated pyridines or of the methyl pyridiniums, shows only single molecular ions whose intensity ratios relative to a fixed internal standard (tetrapentylammonium) are independent of total concentration over a ten-fold concentration range. We take this as evidence that matrix suppression is negligible in this system and that the intensity ratios directly reflect the composition of the samples. In the example above an initial 1
:
1 ratio of alcohols produced a 16
:
28
:
56 mixture of diC18:1, dibrC20:0, and (C18:1)(brC20:0). Details of the synthetic and analytical procedures and compound characterization data are available in the ESI.†
In some experiments, it was obvious that the more hydrophobic pyridinium lipids (diC20:0 and dibrC20:0) were not entirely incorporated into the lipid dispersion. Electrospray-mass spectrometry (ESI-MS) analysis established that the ratio of dibrC20:0:
EPC in an aqueous dispersion is significantly decreased from the formulated 1
:
1 ratio. Data are available in the ESI.†
Aliquots of the liposome dispersion were added to plasmid DNA at molar charge ratios (CR; pyridinium N: phosphate P) of 0.5, 1.5, 3, 5 and 10. Dynamic light scattering analysis showed general increase in average diameter from liposome to lipoplex (0.3–10 μm) that is dependent on the charge ratio. A similar trend of increased lipoplex size compared to liposome is reported as due to fusion in the buffer solution, electrostatic adhesion of neighbouring liposomes and/or from lipoplex aggregation.55,56 In the highly diluted samples used in this study, DLS measurements are influenced by the presence of a small number of large aggregates. The largest particles are associated with intermediate values of CR (1.5–5) in most cases. The polydispersity of all formulations is high (0.2 < PDI < 0.9); the polydispersity of the majority of formulations were lower than those of the corresponding liposomes. Previous studies of in vitro transfection revealed that large lipoplexes of the sizes observed here are effective in transfection.33–35 The hydration method of lipoplex formation results in high variability in size and dispersity. No relation could be found between lipoplex size and either transfection efficiency or cell viability, but it is likely that some of the observed variation arises from this source. Particle sizing data are available in Table S1 in the ESI.†
Pyridinium lipida | Co-lipid | Charge ratio (N![]() ![]() |
Lipid phase | Lattice parameter (δ, Å) | S mix |
---|---|---|---|---|---|
a Lipoplexes formulated with a 3![]() ![]() ![]() ![]() ![]() ![]() |
|||||
diC12:0 | DOPE | 1.5, 3.0 | Hexagonal | 63 | 1.09 |
diC14:0 | DOPE | 1.5, 3.0 | Hexagonal | 69 | 1.09 |
diC16:0 | DOPE | 1.5 | Hexagonal | 68 | 1.09 |
diC18:0 | DOPE | 1.5 | Hexagonal | 76–84 | 1.09 |
diC18:1 | DOPE | 1.5 | Hexagonal | 68 | 1.10 |
diC12:0 | Cholesterol | 1.5, 3.0 | Hexagonal | 80, 84 | 1.08 |
diC14:0 | Cholesterol | 3.0 | Hexagonal | 84 | 1.08 |
diC16:0 | Cholesterol | 1.5 | Hexagonal | 82 | 1.09 |
diC18:0 | Cholesterol | 3.0 | Lamellar | 70 | 1.09 |
diC18:1 | Cholesterol | 1.5 | Hexagonal | 82 | 1.10 |
diC9:0 | DOPE | 1.5 | Lamellar | 63 | 1.08 |
diisoC9:0 | DOPE | 1.5 | Lamellar | 68 | 1.23 |
diC20:0 | DOPE | 3.0 | Hexagonal | 64 | 1.10 |
dibrC20:0 | DOPE | 3.0 | Hexagonal | 68 | 1.38 |
diC9:0 | Cholesterol | 1.5 | Lamellar | 68 | 1.08 |
diisoC9:0 | Cholesterol | 1.5 | Lamellar | 68–85 | 1.22 |
diC20:0 | Cholesterol | 3.0 | Lamellar | 71 | 1.09 |
dibrC20:0 | Cholesterol | 3.0 | Lamellar | 70 | 1.37 |
Direct protection of encapsulated pDNA from hydrolysis by DNase I was also assessed. Fig. 1B gives a representative example; full data are available on pages 18–23 in the ESI.† EPC:
cholesterol and EPC
:
DOPE afford some level of protection at CR 1.5 and above and again, the majority of pyridinium-containing formulations behave similarly. Notable exceptions include diC9:0 and diisoC9:0 which show only low levels of DNA protection at CR 3 and above, and diC20:0 and dibrC20:0 which are protective only at CR 5 and above. Some formulations based on mixtures and mixed chain species diC18:1, dibrC20:0, and (C18:1)(brC20:0) show a maximal protection at intermediate 1.5 < CR < 5 with CR 10 producing a low level of protection. This issue is discussed in a subsequent section. Data for these experiments are available in the Biological data – Figures section of the ESI.†
Transfection experiments with controls and replicates limit the number of independent formulations that can be assessed in parallel within a single experiment to about thirty, or six compounds at five levels of CR. Independent day replicates gave the same relative transfection efficiencies, but direct merging of data from different days required the use of a normalized transfection efficiency (TEnorm) defined as the average transfection efficiency (n = 6) relative to the maximum and minimum levels of transfection observed in the experiment. The maximum was usually the Lipofectamine2000™ result and the minimum was the cells alone without added DNA or lipids. The values of TEnorm from comparable experiments conducted months apart differed by less than 15%. This is comparable to the overall relative standard error for different conditions within a single experiment which indicates that the normalization process used does not introduce significant bias at the same time as it allows comparisons between experiments conducted over many months. The numerical transfection data (raw and normalized) are available in the ESI.†
Our earlier approach results in a mixture of meanings for the terms involved. As noted above, V, a0, and lc are the equilibrium tail volume, head group area, and critical chain length of the lipid. These are experimental parameters and do not directly reflect a physical volume, area, or length.60 By use of an additive model, we implicitly assume they do take a physical meaning. The success of the previous predictions depends in part on the cancellation of errors between the three terms in the derived parameter S.50 We wish to continue to treat the additive approach as if it generated physical parameters, so we require estimates for cholesterol based upon the same approach, recognizing that we freeze what we know to be a variable quantity – the cholesterol head group area. The length of cholesterol is certainly 17.5 Å either directly from experiment or by bond length additive methods. We note that the idea of “head-group” either in its equilibrium sense or a physical sense relates to an incompressible area below which steric repulsions rapidly increase the potential energy.44 For cholesterol this is the span C-2 to C-6. Using this region as head group we derived a value of a0 = 27.5 Å2 and V = 545 Å3. Together these give a value of S = 1.14, in reasonable agreement with the previous value of 1.20.
Our previous approach ignored the counterion of the cationic lipid in the determination of the head group area. A closer reading of the background to additive partial molar volumes suggested the parameters applied only to complete salts.61 We have incorporated a chloride ion in the current model. A table of the molecular parameters utilized is given in Table S4 in the ESI.†
partition term = clog![]() ![]() ![]() | (1) |
shape term = |S+ − 1| | (2) |
volume term = (nexp/nlat) − (nexp/ncyl); nexp < nlat | (3a) |
volume term = 1 − (nexp/ncyl); nexp ≥ nlat | (3b) |
The clogP terms of eqn (1) are computed for the mixture of lipids, the mixture of cationic lipids, and the mixture of neutral lipids with subscripts mix, +, and 0 respectively. We note a problem with this term that allows it to approach ∞ as the difference in the denominator approaches 0; we resolve this issue below. Eqn (2) focussed solely on the computed shape parameter of the cationic lipid. The volume filling term (eqn (3a) and (3b)) is based on the number of lipids per DNA base pair (n) required to fulfil the volumetric requirements of a hexagonal lattice of lipid-coated DNA and of a cylinder of bilayer-coated DNA (subscripts lat and cyl respectively) relative to the experimentally set number of lipids per base pair (subscript exp).
Despite the adequate correlation that TIPSV provided for our initial dataset, that result was clearly adventitious as illustrated in Fig. 2A for a “training” subset of pyridinium lipids of varying linear hydrocarbon chain length diC12:0 to diC20:0 at a CR = 3 with two different neutral lipids. A closer look at the data suggested that there was an optimum chain length at di14:0 (Fig. 2B). Addition of data for diC18:1 and dibrC20:0 shows that the variation in transfection cannot be related to partition; diC20:0 and dibrC20:0 have closely similar clogP but differ significantly in activity. The parameter of merit is length of the longest carbon chain, although in itself this is incomplete as diC18:0 is much less active that di18:1.
The problem with the partition term was noted above so we explored alternative expressions. Eqn (4) is one transformation of the same terms used previously and can reach a maximum value of 1, similar to the volume term of eqn (3a) and (3b).
partition term = 1/(clog![]() ![]() ![]() | (4) |
This change made no change to the fortunate correlations observed previously, or to the lack of correlation in Fig. 2A. We note as well that the shape term (eqn (2)) is not an independent term as it contains the same geometric parameters of volume, length and area as are used in the volume filling term (eqn (3a) and (3b)). However, these are technical refinements that do not address the core issues.
Hydrocarbon chain length directly influences two physical properties of lipid mixtures. One is the melting behaviour of the pure lipid components which scales with carbon number for a range of head groups and lipid structures.62 Another parameter is the miscibility of lipid mixtures which is directly related to chain-length mismatch in simple cases.63 Previous studies that have concurrently investigated the relationship between lipid phase behaviour and transfection have shown that the melting temperature of a binary mixture of EPC with saturated C14 and C16 chains plays a critical role in transfection efficiency.13,42,64 The most effective compositions are those where the melting temperature is close to the temperature of the transfection experiment (37 °C). Lipid mobility is enhanced at the melting temperature via the inevitable formation of defects as the structured phase melts.62
Our ternary lipid mixtures are unlikely to have simple phase behaviours and it is impossible to predict their true nature. However, it is likely that as the chain length of the pyridinium component increases, there will be an increase in the “melting” temperatures of the mixture or a lipid redistribution between coexisting phases that may be present. If there are fluid phases present, then variation in chain length will have little effect on transfection. However, above some threshold value of pyridinium chain length the melting temperature may shift to higher than the experimental temperature with an increasingly negative effect on transfection efficiency due to decreased lipid mobility.
We attempt to capture this in a melting term (eqn (5a)–(5c)):
melting term = 1/(Cn+ − Cnmix); Cn+ ≥ Cnmix | (5a) |
melting term = 1; Cn+ < Cnmix | (5b) |
melting term = 1; |1/(Cn+ − Cnmix| > 1 | (5c) |
The focus is on Cn, the number of independent rotors in the longest continuous chain of either the pyridinium lipid (subscript +) or the lipid mixture (subscript mix; determined as the molar weighted average of the Cn of the components present). This approach follows directly from the correlations derived from both melting and mixing of pure lipids and binary mixtures.62,63 This term will vary only under the conditions of eqn (5a) to diminish the overall TI under conditions where the pyridinium lipid significantly pushes the mixture towards higher melting temperature. If the pyridinium chain length is shorter than the molar weighted average of the mixture then the pyridinium lipid will play a melting temperature reducing role in the overall phase behaviour and the term will take its maximum value (eqn (5b)). Similarly, if the pyridinium is roughly comparable to the mixture and the denominator of (5a) becomes large, eqn (5c) ensures the term does not exceed the maximum of 1.
For linear hydrocarbon chains Cn is the total carbon length less the two end groups as this defines the number of independent rotors. For an oleyl chain (18:
1) this value is 8 as previously reported.62,63 Cholesterol behaves less as a chain and more as a diluent to reduce the cooperative length of the phase transition;62,63 we assume it behaves as if it were a diluent of the same overall length as a C16 hydrocarbon chain as these two species have very similar value of lc (17.4 Å) and thus the same value of Cn ( = 14).
We can therefore write an expression for TIPVM as a product of there terms:
TIPVM = partition term × shape term × melting term |
The required terms are eqn (1), (3a)–(3b) and (5a)–(5c) respectively. This produces values between 0 and 0.1 for the range of clogP values in the experiments so we multiply by ten to get a convenient scaling. Does this approach work any better than in Fig. 2A? The results for the same dataset are given in Fig. 3. Visually the correlation is more persuasive. Statistically we fit a line to the data using the known standard errors in the normalized transfection efficiencies to give the results shown with an r2 = 0.89. There is obvious scatter but there is a basis for deeper exploration of TIPVM beyond the initial 14 formulations of the “training” subset in Fig. 3.
![]() | ||
Fig. 3 Normalized transfection efficiency as a function of TIPVM. The experimental data is the same as Fig. 2. |
Although TIPVM has removed the lipid shape parameter as an explicit term, pyridinium lipid shape is presumed to be an important factor at the point of lipoplex release of cargo within the cell.13,40 The range of S for the series of compounds correlated in Fig. 3 is quite limited (S = 1.21 ± 0.02) resulting in a narrow range of values in the transfection experiments (Smix = 1.09 ± 0.01); we require compounds with high values of S as a means to probe TIPVM more extensively. The bis-isononyl ester diisoC9:0 (S = 1.67) proved to be quite cytotoxic (see below), but dibrC20:0 (S = 2.16) had acceptably low toxicity and at high CR was an active transfection agent. As noted above, the lipid films and initial liposome dispersions of the dibrC20:0 resulted in lipid fractionation. The lamellar phases observed by SAXD for samples with high Smix (Table 1) are consistent with this type of fractionation.
Fortunately, binary blends of diC18:1 and dibrC20:0 were apparently jointly taken into the liposomes so it was possible to explore the effect of S using mixtures of these two components. For example, a 1:
1 blend of these components produces an effective S = 1.70 which in turn produces Smix = 1.24 when mixed with EPC and a co-lipid. Fig. 4A shows normalized TE as a function of TIPVM for binary blends with values of Smix ranging from 1.24–1.33. The correlation is similar (r2 = 0.92) to that of Fig. 3.
As a further check we directly prepared ternary mixtures from a 1:
1 mixture of the starting alcohols. ESI MS analysis of the product mixture gave the composition 16
:
28
:
56 mol% diC18:1, dibrC20:0 and the mixed chain species (C18:1)(brC20:0) with an effective S = 1.70, which in turn produced lipoplexes with Smix = 1.28. A binary blend of diC18:1 and dibrC20:0 in a ratio of 32
:
68 produces the same final Smix. Fig. 4B shows normalized TE as a function of TIPVM for these ternary mixtures and binary blends, all having the same value of Smix. The correlation in this case is poorer than previously (r2 = 0.70) and appears to be dependent upon the neutral co-lipid. For a given co-lipid the correlations are comparable to one another and the previous datasets (DOPE r2 = 0.86; cholesterol r2 = 0.77). It is possible that some differential lipid fractionation occurred as these systems do involve dibrC20:0.
The full dataset consisted of 123 independent conditions conducted in independently duplicated triplicates (see ESI†). Some conditions were quite toxic to the cells leaving 90 conditions where cell survival exceeded 20% (50 with DOPE; 40 with cholesterol as co-lipid). The selection of 20% survivorship is both arbitrary and quite low. It is largely guided by two factors: (1) Lipofectamine2000™ is quite toxic in this cell line and in some experiments cell viability in the Lipofectamine2000™ treatment is about this level; (2) We are interested in capturing the full effect of structural variables on transfection so we need to consider as wide a range of compounds as possible, including those with toxicity that precludes any subsequent practical use.
The correlation of the full dataset with TIPVM weighted by the known standard errors in normalized TE is given in Fig. 5. The overall r2 for this weighted data is 0.7; there is no difference between the cholesterol and DOPE subsets considered separately. The scatter is obvious in the dataset but does indicate an underlying correlation with TIPVM which in turn is based upon a handful of parameters which are not adjusted.
![]() | ||
Fig. 5 Normalized transfection efficiency as a function of TIPVM for all available data on pyridinium lipids with EPC and co-lipid. |
A second factor controlling the toxicity of a formulation is the lipophilic character of lipid mixture (Fig. 6B). Cholesterol has a substantially lower clogP (9.9) than DOPE (14.8) so the clog
Pmix of a given pyridinium lipid in a mixture has an offset that is dependent upon the co-lipid; the offset is 1.9 at the ratio of cationic to neutral lipids used. Although there appears to be a general increase in cell survival as clog
Pmix increases, the correlation is poor and other factors must be important. One of those factors appears in the comparison of the toxicity of binary blends and ternary mixtures of diC18:1, dibrC20:0, and (C18:1)(brC20:0) as given in Fig. 6C; clog
Pmix alone is insufficient to rationalize the observed toxicity. Fig. 6C and 4B are derived from the same experiments, the Smix is fixed at 1.28, and the clog
Pmix values for a given co-lipid are the same (DOPE 13.3; cholesterol 11.4). Even with close similarities the two binary blend sample series are much less toxic than the two ternary mixture series. Note as well that this is independent of the relative transfection efficiencies as given in Fig. 4B; transfection and toxicity are uncorrelated. The key difference between the binary blended samples and the ternary mixture is the presence of (18
:
1)(br20:0) in the latter. The toxicity evident in Fig. 6C must be due to this lipophilic component.
The two most toxic pyridinium compounds uncovered are diisoC9:0 and (18:
1)(brC20
:
0). They lie at the extremes of the range of clog
P for the set of compounds investigated but have high and similar S values (1.67 and 1.60 respectively). They are also the least “lipid-like” compounds in this study (branching in the chains; chain-length mis-match). It is possible that neither of these compounds is well-tolerated in cell membranes and as a result these compounds are more mobile in the cytoplasm and can reach cellular targets more easily than less mobile components.
The experimental work in this study is based on a 1:
1 ratio of a pyridinium compound and EPC. To further extract the toxic effect of the pyridinum compound from the toxic effect of EPC this ratio would have to vary. More detailed conclusions on the structural factors that control the toxicity of a given pyridinum compound await that dataset.
It is doubtful that we have identified all the key components even within this relatively limited system as we have restricted ourselves to pyridinium-based cationic lipids. There is ample evidence that the overall charge on the lipoplex plays a significant role in transfection.27,29,32 We have been able to ignore this factor as the pyridinium lipids have fixed charges and are mixed in fixed proportions to other fixed-charge species (EPC, DNA) to produce – we assume – lipoplexes with approximately constant surface charge. The same type of study could examine the issue of charge as a contributing term using ionisable lipids (polyamines) in conjunction with related quaternized ammonium compounds.
Another as yet unexplored dimension is the influence of cell type on transfection with defined lipoplexes. Commercial transfection kits promise a broad spectrum of activity. Is this a consequence of a general structure-derived source of activity as our correlation here suggests? The same question arises with respect to the nature of the nucleic acids.
In the end, the robust packaging of the cargo is the key issue. We are a long way from packaging that is too robust based on supramolecular assembly from lipidic components. As noted above, we have relatively poor control over the particle size and dispersity with the preparative method utilized. It is possible that some of the scatter between replicates and some of the deviations from the overall trend are simply due to variations in the physical parameters of the assemblies. Microfluidic mixing is an emerging technology for the formation of lipid-nucleic acid nanoparticles for siRNA gene silencing.65,66 Such techniques offer considerable control over particle size and are an important next step in establishing that the correlations reported here are themselves not simply adventitious.
Footnote |
† Electronic supplementary information (ESI) available: Synthetic methods and compound characterization; biological methods and data. See DOI: 10.1039/c6ob00041j |
This journal is © The Royal Society of Chemistry 2016 |