Paola
Posocco
*a,
Sabrina
Pricl
*a,
Simon
Jones
b,
Anna
Barnard
b and
David K.
Smith
*b
aMolecular Simulations Engineering (MOSE) Laboratory, Department of Materials and Natural Resources (DMNR), University of Trieste, Piazzale Europa 1, 34127, Trieste, Italy. E-mail: POSOCCO@dicamp.units.it; SABRINA.PRICL@dicamp.units.it
bDepartment of Chemistry, University of York, Heslington, York, UK YO10 5DD. E-mail: dks3@york.ac.uk; Fax: +44 1904 432516
First published on 9th July 2010
This edge article reports the multiscale modelling of spermine-functionalised dendrons designed to bind DNA and deliver it into cells. The modelling provides an insight into the mode of self-assembly of the dendrons, in particular aggregate charge density and shape, and hence suggests explanations for some of the unexpected experimental observations. In particular, the self-assembly model helps explain why the first generation dendron is more effective in binding DNA than the second generation analogue, even though the second generation system has a greater number of surface spermine ligands. Effective self assembly of the first generation dendron leads to a high charge density assembled structure – more effective than the larger number of ligands on the second generation dendron – i.e., less is more. The modelling also suggests an alternative shape of self-assembly for the system with two hydrophobic cholesterol units (rather than one) at the dendron focal point – this may help explain why this system shows much enhanced gene delivery in vitro.
Recently, we have developed an interest in the ability of relatively small, multivalent dendrons, with spermine surface groups to bind DNA.6 Spermine is a tetraamine present in all eukaryotic cells at millimolar levels, and plays a key role in the natural processes of DNA folding.7 In previous work we demonstrated that higher generation dendrons bound DNA more strongly, due to their ability to establish a larger number of contacts with the target – i.e. multivalency in action. Furthermore, we have used a combined experimental and theoretical approach to demonstrate that for higher generation systems, the binding was independent of salt concentration, and explained this observation by developing a new paradigm in multivalency, in which some of the surface ligands ‘sacrificed’ their interactions with DNA in order to shield the remaining spermine–DNA interactions from the competitive influence of sodium cations.8
In general terms, however, the synthesis of high generation, well-defined dendrimers with the required levels of purity, is demanding, and we have therefore become interested in developing alternative approaches to enhancing multivalency, whilst limiting the synthetic input. With this goal in mind, we reasoned that rather than using covalent synthesis to put the multivalent array in place, simple dendron self-assembly processes could achieve the same goal, with the self-assembled dendrimer behaving somewhat like a higher generation covalently constructed system (Fig. 1).9 We reasoned that by combining simple dendron chemistry with self-assembly, we could achieve a cost-effective approach to tunable nano-assemblies with high-affinity for biological targets. Recently, there has been increasing interest in the development of self-assembling multivalent arrays to bind diverse biological targets such as glycoproteins,10 integrins11 and collagen.12 In general, to achieve this goal, target-specific ligands are appended with a hydrophobic chain and then inserted into liposome assemblies, with the overall assemblies being multivalent. It should also be pointed out that cationic lipids are well-known to bind and deliver DNA.13 We therefore decided to design surfactant-like dendrons capable of self-assembly into higher order aggregates. There has been considerable interest in the development of self-assembling dendrons.14 Indeed, the groups of Diederich, Florence and others have previously used this type of strategy to generate dendrons with enhanced gene delivery potential.15 We therefore decided to modify our dendrons and provide them with some inherent capacity to self-assemble. In this edge article, we use multiscale modelling to gain a unique understanding of the impacts of self-assembly on biological activity.
![]() | ||
Fig. 1 Schematic figure showing self-assembling dendrons for DNA binding and gene delivery. |
![]() | ||
Fig. 2 Compounds used for DNA binding and delivery studies: Z-G1, Z-G2, Chol-G1, Chol-G2 and Chol2G1. |
DNA binding ability was assessed using ethidium bromide (EthBr) exclusion assays (using calf thymus DNA) and gel electrophoresis studies (using pGL3 plasmid DNA). A summary of the key data from these assays is presented in Table 1. The N:
P ratios represent the ratio between protonatable cationic amines on the dendron and anionic phosphate groups on the DNA required to achieve 50% EthBr displacement (in the EthBr exclusion assay) or to retain the plasmid DNA in the loading well (in the gel electrophoresis assay).
Data from EthBr exclusion assay | Data from gel electrophoresis | Gene delivery | |
---|---|---|---|
N![]() ![]() |
N![]() ![]() |
Max. % of PEI control | |
a The N![]() ![]() |
|||
Z-G1 | 2.70 | 2.21 | <2b |
Chol-G1 | 0.52 | 0.34 | ca. 25b |
Chol2-G1 | 0.49 | 0.85 | 60 |
Z-G2 | 0.76 | 1.07 | ca. 5b |
Chol-G2 | 1.35 | 0.51 | 10 |
It is evident that when there is a simple Z-protecting group at the focal point, the N:
P ratios for DNA binding decrease on going from first to second generation dendron. This indicates that for Z-G2, a smaller amount of amine is required to condense the DNA, and that it is therefore a significantly more effective DNA binder than Z-G1 (Fig. 2). This is in agreement with the hypothesis that the surface spermine groups act as a multivalent array and collaborate in binding the DNA, with the presence of more surface ligands enhancing binding. We have previously used molecular dynamics methods to model the 1
:
1 interaction between the dendron and DNA double helix and hence gain a detailed understanding of the multivalency effects in this system.8 However, when a cholesterol unit is at the focal point of the dendron, surprisingly, Chol-G1 (and Chol2-G1) have significantly lower N
:
P ratios than Chol-G2. This was a surprising result, indicating that the first generation systems were much more effective DNA binders than their second generation analogues – i.e., less is more – even though the cholesterol unit cannot be directly involved in forming interactions with the DNA. It is also worth noting that whilst changing the Z group for a cholesterol unit significantly improved the DNA binding of the G1 dendrons, the same structural modification had a lesser effect on the binding ability of the G2 dendron. The nature of this binding enhancement for Chol-G1 could not be meaningfully modelled using simple molecular dynamics methods on the 1
:
1 dendron
:
DNA complex, because the spermine-functionalised surfaces of Z-G1 and Chol-G1 are, in their own right, identical, yet the degrees of binding are very different.
Gene delivery ability of the dendrons was assessed in (e.g.) HEK293, C2C12 and MDA-MB-231 cells, and compared against branched PEI as a control.6b,16 As is evident, many of the dendrons, although effective DNA binders, were ineffective in gene delivery. Compounds Z-G1 and Z-G2 both required the addition of chloroquine, a known endosomal escape agent,17 in order for transfection to be observed. However, the cholesterol-modified derivatives were better in terms of gene delivery. Compound Chol-G1 still required the addition of chloroquine, but demonstrated quite effective gene delivery, while Chol-G2 and Chol2G1 showed very good levels of transfection even in the absence of chloroquine – notably Chol2G1 was particularly effective and furthermore, a combination of the two vectors improved transfection levels up to ca. 200% of the positive control.16
The cytotoxicity of these dendrons was assessed using a cell titre blue assay.16 All dendrons were relatively non-toxic at transfection concentrations. It was noteworthy that Chol-G1 and Chol-G2 were slightly more toxic than PEI giving 50% cell viability at 10 μg mL−1, whilst Chol2-G1 was somewhat less toxic (50% cell viability at 25 μg mL−1).
In summary, therefore, we found that the addition of a hydrophobic unit at the focal point of the dendron had a significant impact on the ability of the compound to bind and deliver DNA. In our original paper,16 we proposed that the hydrophobic units were able to self-assemble, and hence generate multivalent arrays of spermine ligands in situ, and that therefore Chol-G1 was more effective than Z-G1. However, there were a number of interesting experimental observations which we were unable to fully explain. In particular:
(i) Why are lower generation Chol-G1 compounds more effective DNA binders than their higher generation Chol-G2 analogues – even though they have fewer spermine ligands on their surfaces?
(ii) Why was Chol2G1 significantly more effective in terms of gene delivery into cells than the Chol-G1 – how does the extra cholesterol unit assist gene delivery?
We therefore hoped to probe the nature of the DNA binding and propose answers to these questions using multiscale modelling methods.
Initially, we modelled Chol-G1 and Chol-G2 using this approach, in the absence of DNA, in order to gain an insight into why Chol-G1 unexpectedly showed greater affinity for DNA, even though it has fewer spermine surface ligands. On refining Chol-G1 it was evident that the self-assembly of these ligands into aggregates was a favoured process. In Fig. 3a, the green beads represent the cholesterol units and the yellow beads represent the spermine dendron part. Clearly, the cholesterol units self-assemble within the water. In Fig. 3b these cholesterol aggregates are represented volumetrically as green spheres. On increasing the volume fraction (concentration) of dendron within the box (Fig. 4), the aggregates increase in number, but do not significantly change in terms of morphology or composition. Comparing the modelling of Chol-G2 gives results which appear, at first sight, similar to those obtained with Chol-G1 (Fig. 4). Spherical aggregates are formed, driven by the hydrophobicity of the (green) cholesterol units, which want to separate themselves from the aqueous phase. Analogous to Chol-G1, increasing the volume fraction of the dendron does not significantly change the proposed aggregate morphology. Stable supramolecular dendritic micelles related to these have recently been used to explore biotin–avidin binding.18
![]() | ||
Fig. 3 Mesoscale modelling of Chol-G1 showing aggregation into micellar objects. (a) Yellow beads represent dendrons and green beads represent cholesterol units; (b) green spheres represent hydrophobic regions and yellow beads represent dendron units. |
![]() | ||
Fig. 4 Mesoscale modelling of Chol-G1 and Chol-G2 with increasing volume fraction. Chol-G1 and Chol-G2 form stable micellar aggregates across all volume fractions. |
Following the expectations stemming from geometrical considerations, both Chol-G1 and Chol-G2 predominantly form spherical monodisperse micelles with average diameters of 3.4 and 3.8 nm, respectively. This is a net consequence of their conical shape, with a large cationic headgroup connected to a comparatively small lipophilic tail (vide infra). However, when the aggregates of Chol-G1 and Chol-G2 are considered in more detail, there is a significant difference between them. Fig. 5 represents the effective positive charge density on the surfaces of the aggregates. Clearly Chol-G1 gives rise to a more highly charged aggregate (represented in red) than Chol-G2 (represented in blue). This may seem counterintuitive, because Chol-G2 has twenty seven protonatable amine groups on the surface, whereas Chol-G1 has only nine. However, the self-assembly process is significantly more effective for Chol-G1 than the Chol-G2 analogue, and this means that Chol-G1 gives rise to a much more effectively packed aggregate with a tightly packed positively charged surface. Conversely, for the aggregates of Chol-G2, the surface charge is spread more diffusely. This argument is supported by evaluation of the surface charge densities σm for both micellar systems. In the case of Chol-G1, the estimated aggregation number Nagg = 21 leads to a σm value of 5.3 e nm−2, whilst for Chol-G2Nagg = 7 and, hence, σm = 4.2. In this way, we propose that self-assembly controls the surface charge of the aggregate, and hence the relative affinities of these systems for DNA. In this sense, less is more, and the smaller dendron is actually better able to bind DNA once the multivalency becomes expressed on the nanoscale through dendron self-assembly
![]() | ||
Fig. 5 Mesoscale modelling of Chol-G1 and Chol-G2 representing surface positive charge of the aggregates. Chol-G1 forms aggregates with a significantly higher charge density. |
We then performed modelling studies on Chol2-G1, which is actually the most effective dendron for the purposes of cellular gene delivery. In Fig. 6, the cholesterol units are represented in green and the dendron units in purple. At low volume fractions, this dendron once again formed roughly spherical or short elongated micellar aggregates. However, on increasing the volume fraction, the mode of assembly changed dramatically, to give a cylindrical micellar hexagonal phase morphology, in which the hydrophobic cholesterol tails are located inside the cylinders with the hydrophilic dendritic heads being on the surface. Neither Chol-G1 or Chol-G2 exhibited this type of behaviour, with spherical micelles (supramolecular dendrimers) being the dominant and stable aggregate form. We propose that the combined action of a greater degree of hydrophobicity of Chol2-G1 and of a reduced like-charge repulsion in the dendritic cationic headgroup exerted by the ionic strength encourages this phase change in the nature of the aggregate. Notably, with respect to lipid-based molecules used as nanovectors, members of this Chol-based dendron series have a hydrophobic tail the characteristics of which are (i) a high number of carbon atoms (28 for 1 Chol-moiety vs. 10–14 for, e.g., aliphatic fatty alcohols), (ii) higher steric hindrance, and (iii) higher rigidity of the hydrophobic part, as the four condensed rings constitute a fairly stiff skeleton.
![]() | ||
Fig. 6 Mesoscale modelling of Chol2-G1 with increasing volume fraction. The self-assembled morphology evolves from a roughly spherical micellar form to hexagonally packed cylindrical micelles. |
Israelachvili's rules can be used to predict micellar phases based on molecular-level structural aspects of lipids.19 As the geometry and the degree of order of the eventual aggregate architectures are determined by the enforced curvature in the assembly from the relative sizes of the polar and apolar domains of the amphiphiles, the dimensionless packing parameter P, originally developed for small amphiphiles in water,19 can be generalized and used to define the relative size of the hydrophobic region of the Chol-modified dendrons. Accordingly, the balance between the solvent-phobic and solvent-philic interactions gives rise to an optimal surface area of the solvent-phobic block at the interface between the polar and apolar parts of the molecule (α0). This, together with the length lc and the volume of the insoluble domain νh, contributes to the packing parameter defined as . As a general guideline, spherical micelles are formed when P ≤ 1/3; cylindrical micelles are formed at 1/3 < P ≤ 1/2; if 1/2 < P < 1, bilayers with a spontaneous curvature (vesicles) are produced; for P = 1, planar bilayers will be formed; and, lastly, if P > 1, micellar aggregates with a reverse curvature will be generated.
Simple molecular modeling considerations allowed us to calculate the packing parameter for all Chol-based dendrons in hydrated conditions: the volume νh (391 Å3), and the length lc (16.4 Å) of the hydrophobic tail are close to the data reported in the literature.20 The estimated values of the surface areas at the interface α0 are 96 Å2 for Chol-G1 and Chol2-G1, and 194 Å2 for Chol-G2. By using these values, packing parameters P of 0.24 and 0.12 could be calculated for Chol-G1 and Chol-G2, respectively, thus predicting spherical aggregate formation. For Chol2-G1, however, the estimated value of P = 0.47 theoretically supports the formation of cylindrical micelles for this molecular series. Cylindrical dendritic micellar assemblies have been widely explored by Percec and co-workers.14d
As the concentration of the solution is increased, Chol-dendron–solvent interactions become more significant. In order to minimize the free energy, the local packing of the aggregate may change, leading to a decrease in the curvature of the local assembly. In other terms, derivatized-dendrons, which at low concentration form spherical aggregates, assemble into cylindrical and eventually membrane-like aggregates as the concentration increases. Hence, again as a general rule, on moving from low to high concentrations, spherical micelles pack into cubic phases, followed by hexagonally packed cylinders and, finally, lamellae. This is the case for Chol2-G1, for which the balance between hydrophobic and hydrophilic is such that, in order to avoid contact between the solvent and the cholesterol moiety, a higher level of molecular frustration, and hence a higher curvature is attained on increasing dendron concentration.
These morphological proposals also provide a rationale for the results of the EthBr exclusion assay shown in Table 1. In the case of the assembled systems characterized by smaller micelles (i.e., Chol-G1 and Chol-G2), the overall DNA surface will be less covered (patch-like coverage) which allows for better residual binding of small molecules like EthBr. Therefore, in order to “fill the gaps”, higher concentrations of dendrons are necessary to attain the same result displayed by the cylindrical micelles formed by Chol2-G1. However, Chol-G1 has a significantly higher surface charge density than Chol2-G1 (for which Nagg = 34 and σm = 1.7 e nm−2) and this effect partly compensates for the ‘patchy’ DNA coverage of Chol-G1, making it comparable as a DNA binder to its more hydrophobic counterpart, Chol2-G1.
Importantly, it has previously been demonstrated for cationic lipid gene delivery agents that the mode of aggregation plays a key role during endosomal escape, with hexagonal phases being known to destabilise endosomal membranes most effectively.21 We therefore propose that the multiscale modelling of the Chol2-G1 may offer an explanation as to why this dendron shows the most effective gene delivery and lower cytotoxicity (vide infra).
![]() | ||
Fig. 7 Mesoscale modelling of the interaction between DNA (red) and the aggregates of Chol-G1 (green) – direct contact is highlighted by the yellow circle. |
![]() | ||
Fig. 8 Mesoscale modelling of the interaction between DNA (red) and the aggregates of (a) Chol-G1, (b) Chol-G2 and (c) Chol2-G1. |
To a first approximation, the average number of DNA base pairs covered by the Chol-G1 and Chol-G2 micelles of Fig. 8a,b can be estimated by the simple relationship:
![]() | (1) |
Dendron Chol2-G1 can form cylindrical dendritic micelles; the DNA chains inserted into the regular interspace of the cylinders. We note that such a system would be expected to (i) bind and compact DNA efficiently, (ii) protect DNA from degradation, (iii) achieve efficient DNA release.21 It is generally known that the DNA double helix behaves as a semi-flexible coil; thus, the cylindrical shape of the micelles formed by Chol2-G1 is better able to comply with the wrapping pattern of the semi-rigid DNA helix than the small spherical micelles generated by the other dendrons. As can be easily understood from Fig. 8, this mesoscopic organization of Chol2-G1 with DNA results in more efficient and well-organized DNA compaction and, hence, potential protection from degradation during DNA delivery.
Applying eqn (1) to this system yields a value of nbp = 44.7. In line with the underlying discussion, the cylindrical shape of the resulting Chol2-G1 dendron aggregates, although having overall lower surface charge density (σm = 1.7 e nm−2), can not only exploit the surface charge distribution more efficiently by following the seemingly cylindrical shape of the semi-rigid DNA duplex, but can also complex a considerably larger number of DNA bases per aggregate with respect to its less hydrophobic counterparts.
It is also not difficult to infer from the mesoscale simulation images, that the spherical micellar systems actually contain more water space than the cylindrical ones. Again at first approximation, given the average center-to-center distance of the cylinders dcyl/cyl = 5.9 nm, and some simple geometrical considerations, the volume fraction of water ϕcylw in the hexagonal lattice formed by Chol2-G1 can be calculated as 0.42. Although the same approach cannot be applied to the alternative micellar systems due to the absence of a long-range regular structure, the comparison of the density maps of water for the Chol2-G1/- and Chol-G2/DNA systems support this concept (Fig. 9). The higher water content in the Chol-G2 spherical micellar system may imply that a significant fraction of dendron and DNA charges are wasted in terms of mutual interactions. This is in line with the illustrations in Fig. 7 and 8, as a notable fraction of the periphery of the spherical dendron micelles fail to contact DNA and vice versa.
![]() | ||
Fig. 9 Three-dimensional density distribution of water molecules in the Chol2-G1/- (top, left) and Chol-G2/DNA (bottom, right) complexes. The higher level of intermicellar hydration is apparent in the Chol-G2-based system. |
Very importantly, the coupling of simulation/experimental evidence concurs to show that these Chol2-G1/DNA complexes, capable of forming a cylindrical phase, efficiently transfect cells. It is tempting to speculate that the existence of a regular and continuous DNA sub-organisation may facilitate the release of the DNA cargo, as in principle, all DNA is accessible once a part of it is exposed to the inside of the cell. This is in contrast to both the spherical micellar systems, where DNA in less efficiently wrapped and, hence, less protected and organised, and the widely employed lipid-based lamellae-forming nanovectors, in which the enveloping lipid bilayers have to be removed layer-by-layer to release all of the encapsulated DNA.
The main findings from the present work allow us to speculate on the possible structure–activity relationship according to the following line of thought. Endocytosis is the dominant mechanism of cellular entry by dendron-based nanocarriers.3 After cellular uptake, the nanovectors and their cargoes must escape from endosomes, in order for DNA to progress towards the cell nucleus and avoid the rather rapid degradation of the endosome content. Endosomal escape is usually initiated by destabilisation of the endosome membrane that results in flip-flop of the anionic lipids predominantly located on the cytoplasmic face of the membrane.13 These anionic lipids can then laterally diffuse into the complex and form charge-neutralised ion pairs with the cationic part of the nanocarriers. This displaces the cargo DNA from the complex and generates a mechanical and osmotic stress that further contributes to the rupture of the endosomal bilayer and release of DNA into the cytoplasm. The importance of endosomal escape as a barrier to transfection can be assessed by performing transfection experiments in the presence of chloroquine.17 Several experiments performed with cationic lipids as DNA carriers supported the idea that nanovector/DNA complexes with lower surface charge densities σm may remain trapped within the endosome,15d as they are not able to induce either flip-flop of the membrane lipids or osmotic pressure increase for membrane bursting. At the same time, however, σm is not the only transfection modulating physicochemical parameter – inherent biomembrane fusogenicity and complex rigidity play major roles in transfection efficiency.
With these ideas in mind, we can postulate that, in the case of Chol-G1, for which an enhancement of the transfection efficiency was achieved only in the presence of chloroquine, the high charge density which will favor endosomal escape is counteracted by the more compact packing of the amphiphiles in the spherical micellar aggregates (Nagg = 21). This, in turn, imparts an overall higher rigidity to the corresponding nanovectors which are thus less likely to fuse efficaciously with both the cellular and endosomal membranes. In the case of Chol-G2, the larger cationic headgroup and consequent higher Coulombic repulsion among like-charges leads to the formation of micelles with a smaller aggregation number (Nagg = 7) and a smaller surface charge than Chol-G1. This results in a lower DNA affinity but, at the same time, in a less tightly packed and hence more flexible micelle with respect to its lower generation counterpart. Therefore, the enhanced transfection efficiency of Chol-G2 over Chol-G1 (see Table 1) could be ascribed to a better capacity of this nanovector to fuse with the anionic lipid components of the endosomal membrane and to partially effective complex dissociation in the endosome/cytoplasm due to the lower DNA–micelle attractive interactions. Lastly, in the case of the most effective vector Chol2-G1, we propose that the coupling of a different mesoscopic assembly, aggregate dimensions, surface charge value and, last but not least, DNA encasing and wrapping patterns probably represents the best compromise between binding and transfection efficiency.21 The result obtained through the use of a relatively small and rigid amphiphile of unusual shape (a large hydrophobic portion coupled with a very large hydrophilic head) is a cylindrical nanostructure showing an enhancement of binding for the target DNA and improved transfection efficiency. This enhancement may be attributed to three effects: (1) the reduced competition of water molecules for the binding of DNA within the partially desolvated interface, (2) the presence of a flexible linker moiety which endows the resulting cylindrical micelle with the proper degree of flexibility for effective fusion with the cellular/endosomal membrane, and (3) the enhanced multivalency effect derived from the presence of multiple DNA anchoring points in close proximity and structurally well matched to the DNA.
We propose that these factors, as observed from the modelling, help explain the more effective gene delivery observed with this vector.
Initially, the ability of Chol-G1 to self-assemble in the absence of DNA was tested using a Nile red solubilisation fluorescence assay, in which hydrophobic Nile red dye is used to report on the formation of a hydrophobic domain within an assembled nanostructure.23 Nile red is only solubilised when the dendron itself self-assembles into a nanostructure with a hydrophobic domain. Fig. 10 demonstrates that self-assembling behaviour was observed for Chol-G1, with a discontinuity in the fluorescence intensity of Nile red, plotted against increasing dendron concentration, occurring at the critical aggregation concentration (CAC), above which assembly occurs. Using this approach, the CAC could be quantified as 15 μM. This demonstrates that these dendrons are capable of self-assembly processes.
![]() | ||
Fig. 10 Fluorescence intensity data from Nile red solubilisation assay to determine the critical aggregation concentration (CAC) of Chol-G1. |
We used dynamic light scattering (DLS) methods to monitor the aggregation of all five dendrons in the presence of pGL3 plasmid DNA. The dimensions of the uncomplexed plasmid DNA were 66 nm. As is clear from Fig. 11, there are significant differences between the dendrons in terms of their ability to aggregate with plasmid DNA. In the presence of plasmid DNA, all of the dendrons modified with cholesterol units gave rise to larger aggregates than the previously investigated Z-G1 and Z-G2.6b Dendrons Chol-G1 and Chol-G2 both formed assemblies which increased in size as the loading of dendron went up (and the N:
P ratio consequently increased). Interestingly, Chol-G1 formed significantly larger aggregates with DNA than Chol-G2, even though it is a smaller dendron. This is the reverse of what was observed for Z-G1 and Z-G2 where the latter compound gave rise to larger aggregates with plasmid DNA on account of its multivalency. These observations may be indicative of the greater surface charge of the self-assembled Chol-G1 supramolecular dendrimers compared with those of Chol-G2, as was proposed from the multiscale modeling – hence leading to more effective DNA binding and the formation of larger combined aggregates. This is in agreement with simple self-assembly arguments, in that the small but highly charged aggregates formed from Chol-G1, encourage the formation of larger overall aggregate structures, whereas the less highly charged Chol-G2 assemblies give rise to less aggregation of the plasmid DNA. As the loadings of both dendrons increase, and with it the overall positive charge, the observed diameter, and hence degree of dendron aggregation also increase. Intriguingly, the interaction of Chol2-G1 with DNA as monitored by DLS was relatively independent of N
:
P ratio, and was markedly different in nature to the aggregation of either Chol-G1 or Chol-G2, giving aggregates with an apparently smaller diameter. However, the aggregates were still significantly larger than those observed for Z-G1. This would be in agreement with the results of multiscale modeling in which Chol2-G1 was observed to assemble into a different nanoscale morphology, with much more effective packing of the DNA, as a consequence of the greater degree of hydrophobicity at the dendron focal point. These light scattering studies are, therefore, in qualitative agreement with the multiscale modelling study.
![]() | ||
Fig. 11 Dynamic light scattering results indicating the polyamine![]() ![]() ![]() ![]() |
Ethidium bromide (EthBr) exclusion assays were performed in PIPES buffer at pH 7.4, in the presence of 150 mM NaCl. The concentration of calf thymus DNA (per nucleotide) was 2.0 μM, and the concentration of EthBr was 2.54 μM. The concentration of ligand required to displace 50% of the EthBr was determined by fluorescence spectroscopy. Full details can be found in ref. 6a, 16.
Gel electrophoresis assays were performed with pGL3 plasmid DNA and increasing amounts of dendron. Full details can be found in ref. 6b, 16.
Nile red solubilisation fluorescence assays were performed in aqueous solution with increasing concentrations of dendron, and excitation at 550 nm. Nile red emission was monitored at 635 nm.
Dynamic light scattering studies were performed in PBS buffer using dendrons and pGL3 plasmid DNA. The solutions were vortexed and then left to stand for 10 min before the sizes were measured on a Brookhaven instruments corporation 90Plus particle size analyzer at 25 °C, at 1 min per run and 5 runs per sample.
Gene delivery studies were performed using a variety of cell lines, including HeLa, HEK293, C2C12 and MDA-MB-231 cells using a luciferase assay with pGL3. Data quoted in Table 1 are the maximal observed transfection levels at any N:
P loading of the dendron. Once again, full details of degrees of transfection and cell lines at different loadings can be found in the original experimental literature – ref. 6b and 16.
Similar to MD, in a DPD run, the time evolution of each DPD particle can be calculated by solving Newton's second law:
![]() | (2) |
All of the above forces are acting within a sphere of interaction or cut-off radiusrc, which is the length scale parameter of the system. γ and σ are the coefficients of the dissipative and random forces, respectively. In analogy with the fluctuation-dissipation theorem, Español and Warren26 obtained the detailed balanced condition for the DPD as:
![]() | (3) |
This fundamental DPD parameter aij can be determined as a function of the DPD number density ρ and system temperature by matching the compressibility condition:27,28
![]() | (4) |
The dissipative and random weight functions take the general form:29
![]() | (5) |
In this edge article, we used standard DPD particles to represent the chain beads, and the DNA molecule was modeled as a string of DPD particles which are sequentially connected through the WLC spring with the pairwise interaction force given by:
![]() | (6) |
However, the DPD particles can pass through each other without restriction. In order to enforce self-avoidance of the chain, we therefore applied the Lennard-Jones potential ULJ between beads as in Symeonidis et al.31 It is noted that the Lennard-Jones potential applied here is defined at a mesoscopic level and not at the microscopic level as in molecular dynamics. When the beads get close, the repulsion force derived from the Lennard-Jones potential will dominate and push them apart. The repulsion force is calculated as follows:
![]() | (7) |
Moreover, in order to correctly derive the electrostatic interactions between charged beads (i.e. present in dendron and DNA molecules), the electrostatic force FEij between charged beads i and j was analyzed following the approach reported in Groot's work.32 According to this study, the electrostatic field is solved by smearing the charges over a lattice grid, the size of which is determined by a balance between the fast implementation and the correct representation of the electrostatic field.
Concerning the modeling of the mesoscale architecture, we modeled the dendron molecules at a coarse-grained level using a branched and flexible amphiphilic chain made up of three bead types: a neutral bead type, Ch, as a cholesterol building block, a positively charged bead type Sp, as spermine repeating unit, and a neutral bead type L, linking the two parts together (when present). According to the chemical structures reported in Fig. 2, and satisfying the theoretical criteria for which each bead type should be of similar volume, the coarse-grained models of each dendritic system were obtained, as shown in Fig. 12. Solvent molecules were simulated by single bead types W, and an appropriate number of counterions of a charge of ±1 were added to preserve charge neutrality and to account for the experimental solution ionic strength. The inclusion of explicit counterions was necessary because counterion condensation and the interactions between the counterions and the charged groups may affect the complex structure to a certain extent. Each DNA chain was represented as a series of 500 negatively charged particles D linked linearly through springs.
![]() | ||
Fig. 12 Schematic representation of the DPD models of the Chol-modified dendrons considered in this work. The colours of the beads stand for hydrophilic spermine beads Sp (yellow), hydrophobic cholesterol bead Ch (green), and linker bead L (purple). |
Having set the reference volume of one bead and the density of the system to a value of ρ = 3, then the cut-off radius rc, which sets the length scale of the system, corresponds approximately to 0.76 nm. All simulations were performed in a 3D-periodic cubic box of 20 × 20 × 20 r3c, equivalent to a bulk system of 3510 nm3. The appropriate number of dendrons and DNA molecules were added to the simulation box in order to fit experimental concentrations.
A characteristic time scale was then defined as . Accordingly, the physical time scale unit in the present simulations was τ ≈ 0.009 ns. Since 5 × 106 time steps were performed in each DPD run, the total physical time of each calculation was about 2.3 μs.
The intra- and intermolecular interactions between DPD particles are expressed by the conservative parameter aij, defined by eqn (4). This quantity accounts for the underlying chemistry of the system considered. In this work, we employed a well-validated strategy that correlates the interaction energies estimated from atomistic molecular dynamics (MD) simulations to the mesoscale aij parameter values.33 Following this computational recipe, for each dendron type a 30 base-pair double-stranded DNA with a randomly generated sequence was generated with the Nucgen module of AMBER 9.34 The molecular models of the three dendron structures (Chol-G1, Chol2-G1, and Chol-G2) were then built and geometry optimized using the Antichamber module of AMBER 9 and the GAFF force field.35
The corresponding bimolecular complexes were generated by placing the dendron molecule close to the DNA major groove, solvating the resultant assembly with a box of TIP3P water molecules36 extending at least 10 Å from the solute in each direction. The appropriate number of counterions (Na+ and Cl−) were added to neutralize the system and to generate the 0.15 M ionic strength. Each system was then relaxed and equilibrated by MD in the microcanonical ensemble (constant volume–constant temperature, or NVT) for 200 ps at 300 K, using the Pmemd module of AMBER 9. Equilibrium density was then reached by performing 5 ns isobaric–isothermal (NPT) MD simulations. These were followed by 10 ns production runs under NVT conditions (T = 300 K), using an integration time step for the equations of motion of 2 fs, the Langevin thermostat, a cut-off for long-range interaction of 12 Å, and the Particle Mesh Ewald technique37 to treat the system electrostatics.
The interaction energies between the system's molecular constituents were estimated using the Molecular Mechanics/Poisson–Boltzmann Surface Area (MM/PBSA) methodology,38 and then rescaled onto the corresponding mesoscale segments adapting the procedure described in detail in ref. 30. The bead–bead interaction parameter for water was set equal to aww = 25, in agreement with the correct value of DPD density ρ = 3 (see eqn (4)).24 The maximum level of hydrophobic/hydrophilic repulsion was captured by setting the interaction parameter aij between the positively charged spermine bead Sp and the lipophilic cholesterol bead Ch as 80. The counterions were set to have the interaction parameters of water.29 Once these parameters were assigned, all the remaining bead–bead interaction parameters for the DPD simulations were easily obtained, starting from the atomistic interaction energy values, as described in the previous work.30 The entire set of DPD interaction parameters employed in this work are summarized in Table 2.
a ij | Sp | Ch | L | D | W |
---|---|---|---|---|---|
Sp | 45 | 80 | 30 | 10 | 15 |
Ch | 80 | 25 | 40 | 80 | 80 |
L | 30 | 40 | 25 | 30 | 35 |
D | 10 | 80 | 30 | 45 | 15 |
W | 15 | 80 | 35 | 15 | 25 |
This journal is © The Royal Society of Chemistry 2010 |