In Silico, In Vitro, and In Vivo Studies Indicate the Potential Use of Bolaamphiphiles for Therapeutic siRNAs Delivery.

Specific small interfering RNAs (siRNAs) designed to silence different oncogenic pathways can be used for cancer therapy. However, non-modified naked siRNAs have short half-lives in blood serum and encounter difficulties in crossing biological membranes due to their negative charge. These obstacles can be overcome by using siRNAs complexed with bolaamphiphiles, consisting of two positively charged head groups that flank an internal hydrophobic chain. Bolaamphiphiles have relatively low toxicities, long persistence in the blood stream, and most importantly, in aqueous conditions can form poly-cationic micelles thus, becoming amenable to association with siRNAs. Herein, two different bolaamphiphiles with acetylcholine head groups attached to an alkyl chain in two distinct configurations are compared for their abilities to complex with siRNAs and deliver them into cells inducing gene silencing. Our explicit solvent molecular dynamics (MD) simulations showed that bolaamphiphiles associate with siRNAs due to electrostatic, hydrogen bonding, and hydrophobic interactions. These in silico studies are supported by various in vitro and in cell culture experimental techniques as well as by some in vivo studies. Results demonstrate that depending on the application, the extent of siRNA chemical protection, delivery efficiency, and further intracellular release can be varied by simply changing the type of bolaamphiphile used.Molecular Therapy-Nucleic Acids (2013) 2, e80; doi:10.1038/mtna.2013.5; published online 19 March 2013.

Introduction

In the past decade, efforts to develop RNA-based therapeutic technologies have significantly intensified.[1,2,3,4,5] Triggering RNA interference (RNAi), in particular, has become one of the most widely used techniques for biomedical applications.[1,2,3,4,5,6,7,8,9,10,11,12,13] RNAi uses a mechanism of posttranscriptional sequence specific gene silencing by processing double-stranded RNAs into small-interfering RNAs (siRNAs) used as part of the RNA-induced silencing complex to selectively cleave target mRNA.[14] After the discovery that synthetic siRNAs can be exogenously introduced into cells to activate RNAi,[15,16] this approach has become a powerful method for selective suppression of specific genes of interest in different species, showing potential for use in cancer therapeutics.[3,5,6] However, the biomedical utility of the synthetic siRNAs is limited by several RNA structure-related factors such as the negative charge (uptake by cells that also have a negatively charged surface) and instability in the blood circulation (non-modified siRNAs have a very short half-life in the blood stream, mostly because of degradation by nucleases).[4] These impediments can be overcome by using polymeric or lipid-based carriers to shield the negative charge and provide protection against nuclease activity,[17,18,19,20,21,22] or by using a high-density lipoprotein.[23]

Novel synthetic bolaamphiphiles (bolas) have recently been shown to form complexes with a variety of small and large molecules including peptides,[24,25] proteins,[26] and plasmid DNAs[25,27] delivering them across biological membranes[28] including the blood brain barrier.[18] Bolas also have relatively lower toxicity than cationic liposomes.[29] Bolas have a hydrophobic alkyl chain flanked by positively charged head groups that can potentially interact with negatively charged nucleic acids and promote their delivery into cells. However, the nature of these interactions as well as the possibility to use bolas for delivery of therapeutic siRNAs remains largely unknown. To better understand the essence of bolas/siRNA interactions and to adjust the experimental conditions important for efficiencies of siRNA delivery and release inside cells, we used two structurally related bolas: GLH-19 and GLH-20 (). These bolas have identical hydrophobic domains with acetylcholine head groups attached to an alkyl chain in two different configurations; via either the acetyl group of the acetylcholine (GLH-19), or the nitrogen atom of the choline moiety (GLH-20). These bolas have been shown to deliver their cargo in vivo with different efficiencies[26,30] which can be due to their distinct interactions with the cargo molecules and the cell membranes. In this work, the interactions between these two structurally different bolas and siRNAs were theoretically studied by explicit solvent molecular dynamics (MD) simulations and then experimentally tested in vitro, in cell culture, and in vivo.

Results

MD simulation of bolaamphiphile micelles

As detailed in Methods section, 20 individual bola molecules were aggregated into a single micelle in 0.15 mol/l NaCl solution during 50 nanoseconds simulations. , show the radial distribution functions calculated for two nitrogen atoms (labeled as N1 and N2) representing the hydrophilic head groups of the micelle and two carbon atoms (labeled as C1 and C2) corresponding to the hydrophobic surface of the micelle. The radial distribution functions for these atoms were measured with respect to the center of the micelle. Thus, the different surface properties of the GLH-19 and GLH-20 micelles can be attributed to the different distances of their head groups from the center of each micelle (). For example, the radial distribution function of GLH-19 indicated that the nitrogen atoms were distributed at ~12 Å above the hydrophobic surface of the micelle (), whereas for the case of GLH-20, the distribution was ~4 Å above the surface (). Hence, the hydrophilic head groups of GLH-19 protruded farther from the surface of the micelle and were therefore expected to be more flexible than those of the GLH-20 (,). Consistent with this hypothesis, the mobilities of the head groups, presented as a mean square displacement (cm2) per unit of time (second), were calculated to be 1.8 × 10−7 cm2/s for GLH-19 and 1.0 × 10−7 cm2/s for GLH-20.

The differences between the properties of the GLH-19 and GLH-20 head groups affected the availability of the micelle surface for the chloride ions (Cl−) that are present in solution. The radial distribution function () indicated a higher density of Cl− ions on the GLH-20 surface, compared with GLH-19 for longer periods of time (see ). As a result, a more stable and solid cloud of Cl− ions was formed around the GLH-20 micelles (). Therefore, the electrostatic attractions between the GLH-20 micelle and the siRNA were expected to be more interfered with by the Cl− ion cloud than that of the GLH-19 micelle and the siRNA. Accordingly, a better association with siRNA was predicted for GLH-19.

MD simulation of bolaamphiphile/siRNA complexes

Each MD simulation was initiated by placing the center of mass of the individual siRNA duplex 40 Å away from the center of mass of the GLH-19 or GLH-20 micelle. This distance was sufficient to promote the electrostatic attraction between the positively charged hydrophilic head groups of the bolas and the negatively charged phosphates of the siRNA resulting in the formation of the bolas/siRNA complex. To track the behavior of bola/siRNA complexes during the simulation time, positions representing the 5′-side, middle, and 3′-side of the siRNA were selected. The centroids associated with the phosphate groups of C3, U14 and G25 (counted from the 5′-side of the sense strand) and their complements (G3, A14, and G24 of the antisense strand) were used. The relative distances between each centroid and the nearest bola head groups were continuously measured over the course of 50 nanoseconds (). Within the first ~13 nanoseconds of the simulation, the siRNA duplexes formed stable complexes with both the GLH-19 and GLH-20 micelles (Supplementary Videos S1 and S2). However, during that time, the head groups of GLH-19 were positioned significantly closer to the selected phosphates than those of GLH-20. This suggested a faster and stronger binding interaction with GLH-19 (,). Also, the extent of the interactions within the GLH-19/siRNA and GLH-20/siRNA complexes measured over the total simulation time was quite different. As can be seen from ,, all three phosphates of the siRNA were associated with GLH-19 during the entire simulation (Supplementary Video S1); whereas in the case of GLH-20, the 3′ side of the duplex (the right side of ) was completely released from the interaction during the time periods of 25–30 nanoseconds and 40–45 nanoseconds (Supplementary Video S2). Snapshots of GLH-19/siRNA and GLH-20/siRNA complexes taken at 43 nanoseconds and at 50 nanoseconds () clearly illustrate the difference in their binding affinities, with GLH-19 providing a higher association surface area and protection to the siRNA. These observations can be partly explained by the fact that the GLH-20 micelle was covered by a thicker cloud of Cl− ions (), which interferes with the interactions between the micelle and the siRNA consequently preventing complex formation. Due to the strong electrostatic interactions between the phosphates of the siRNAs, and the head groups of the bolas, the siRNA duplexes on the surface of each bola micelle were slightly bent () with a higher relative deformation for the siRNA occurring for the case of GLH-19/siRNA (,, Supplementary Video S1). The root mean square deviation of the siRNA duplex with respect to its initial structure reached 7.3 Å when it formed the GLH-19/siRNA complex, whereas the root mean square deviation of the siRNA duplex in the GLH-20/siRNA complex reached only 6.2 Å. Therefore, the siRNA that was associated with the GLH-19 micelle experienced a greater deformation due to the stronger interactions with GLH-19.

When the micelles formed stable complexes with the siRNAs (,), a large number of bola head groups were found near the siRNA phosphates due to the electrostatic interactions. We monitored the occupancies of the bola head groups within a 7 Å range of each RNA phosphate during the 30–50 nanoseconds time range of the simulation. Then, the sum of the occupancies for each base pair was analyzed (see Supplementary Figure S1). Both the GLH-19 and GLH-20 head groups populated the 5′-end (nucleotides 1–8 in the sense and the corresponding bases in the antisense strands including the two base overhang at the 3′ end), the center of the siRNA (nucleotides 12–19 in the sense and antisense strand) and its 3′-side (nucleotides 22–25 in the sense and its complements in the antisense strand). Due to the helical nature of RNA, only these particular regions were accessible for interactions with the bolas. However, the occupancies of the GLH-19 head groups near the siRNA phosphates were higher than GLH-20 in most regions. (Supplementary Figure S1, upper panel). Thus, the greater number of GLH-19 head groups that were associated with the phosphates (Supplementary Figure S1, lower panel) provided better protection and higher binding affinity to the siRNA.

The hydrogen bond formations between the siRNAs and the bolas were also monitored (see Supplementary Table S1 and Figure S3). Hydroxyls positioned next to the C1 and C2 atoms of the bolas can form stable hydrogen bonds with the backbone phosphates of the siRNAs. Hydroxyls in GLH-20 located close to the positively charged head groups were brought in close proximity to the negatively charged siRNA phosphates thus, forming a higher number of hydrogen bonds than the hydroxyls in GLH-19 (Supplementary Figure S3a,b). Head groups located close to the minor groove formed stable hydrogen bonds with the H22 and N2 atoms of guanosines (Supplementary Figure S3c,d). In addition, the terminal bases at the 5′ and 3′ ends also formed hydrogen bond interactions with the bola head groups. Interestingly, the sugar edge hydroxyl groups of the siRNA formed very few and relatively weak hydrogen bonds with the bolas. For a 20 nanoseconds time range, the total number of hydrogen bond interactions for GLH-19/siRNA and GLH-20/siRNA were 8 and 11, respectively (Supplementary Table S1). Besides the electrostatic and hydrogen bond interactions between the bola head groups and the siRNA, we found that the overhang and end bases of the siRNA duplex also interacted with hydrophobic chains of the bolas.

The free energies of binding for GLH-19/siRNA and GLH-20/siRNA were calculated using the molecular mechanics Poission Boltzmann/Generalized Born surface area (MM-PB/GBSA) methods (). The contribution of δEmech (which is (δEbond + δEangle + δEdihedral)) in the free energy of GLH-19/siRNA binding was higher due to the deformation of the siRNA. However, the contribution of the van der Waals interactions (δEvdW) in the free energy of binding was lower for the GLH-19/siRNA complexes. The GLH-20/siRNA complex had a lower contribution for the electrostatic interactions (δEele) in the binding energy than the GLH-19/siRNA. However, the solvation energy contributions in the binding free energies estimated by both PBSA and GBSA methods were more favorable in the case of GLH-19/siRNA. The difference in the entropy contribution for the GLH-19/siRNA complex and the GLH-20/siRNA complex was minimal. Therefore, the total binding free energies obtained by the PBSA (~11 kcal/mol) and GBSA (~23 kcal/mol) methods were stronger for the GLH-19/siRNA. Once more, these data indicated a greater binding affinity of the siRNA for GLH-19 than for GLH-20.

To demonstrate the generality of the approach, two additional 40 nanoseconds MD simulations were performed for each bola/siRNA complex. One simulation set used different initial velocities and the other used different micelle surfaces. In all arrangements, a higher number of siRNA phosphates populated the bola head groups corresponding to the GLH-19 micelles (Supplementary Figure S4). The binding affinities calculated by MM-PBSA for GLH-19/siRNA were lower by ~4 kcal/mol (for MD simulations with different initial velocities) and ~13 kcal/mol (for MD simulations with a different initial micelle surface) compared with GLH-20/siRNA. The value of δδG (δG (GLH-19/siRNA) − δG (GLH-20/siRNA)) derived from three independent trajectories was calculated to be −9.3 ± 2.8 (SEM) kcal/mol. This confirms a higher affinity of binding to siRNA for GLH-19 micelles.

As an alternative simulation to illustrate relative propensity of unconstrained bolas to interact with RNA, a shorter duplex (12-mer) was surrounded by individual bolas evenly distributed in the 3D simulation space (Supplementary Figure S5 and S6). To increase the statistical sampling, we placed 6, 15, and 24 single bola molecules around the RNA duplex and analyzed their interactions using 30 nanoseconds MD simulations. For all three cases, GLH-19 head groups showed higher population densities near the RNA phosphate groups. Interestingly, in the case of 24 individual bolas, the relatively short distances between them promoted their aggregation resulting in small micelle formations which in turn interacted with the RNA duplex (Supplementary Figure S6). All these MD simulations additionally prove a better binding affinity of GLH-19/RNA complexes.

Experimental concentration dependence and structural characterization of GLH-19/siRNA and GLH-20/siRNA complex formations

Binding affinity. To assess the abilities of GLH-19 and GLH-20 bolas to interact with nucleic acids, bolas at different concentrations (2, 2.5, 3, and 3.5 µg/ml) were introduced to fixed amounts of siRNAs (400 nmol/l final). Bolas at higher concentrations (10 µg/ml) were also tested and demonstrated consistent results (data not shown). After mixing and incubating, samples were analyzed by agarose gel electrophoresis as described in methods (). Upon staining with ethidium bromide, free siRNAs, not associated with bolas, appeared as single bands. However, the association of the positively charged bolas with the negatively charged siRNAs resulted in the formation of a complex too big to enter the gel thereby, preventing staining. To recover the siRNAs, a commonly used detergent Triton X100 was added to the bola/siRNA complexes. All the conditions tested show the formation of a complex between the bolas and the siRNAs because no free siRNA labeling was observed in the absence of detergent. Results indicated however that those complexes were very different in their sensitivity to detergent extraction of siRNA. For all the concentrations tested, the GLH-20/siRNA complexes were disrupted by the addition of Triton X100 leading to a complete recovery of the RNAs. Conversely, a very tight binding between GLH-19 and siRNAs was observed as the complex that could only be partially disrupted by the detergent at the lowest concentration (2 µg/ml), whereas at the higher concentrations (2.5, 3, and 3.5 µg/ml), they remained intact in the presence of Triton X100. These results were in a good agreement with the MD simulations that predicted that GLH-19/siRNA complexes have a stronger binding free energy than the GLH-20/siRNA complexes.

Digestion. To investigate the effect of complex formation on the nucleic acids' chemical stability, digestion assays with nucleases were carried out (). In these experiments, for reason of economy, we used a DNA duplex as a model system, with sequences that are identical to those used in the siRNA. The main difference between DNA and RNA is the presence of hydroxyl at the sugar edge. Our MD simulation showed that hydrogen bond interactions between the bolas and the hydroxyls in the siRNA were minor (see Supplementary Table S1). Therefore, the fundamental interactions between the DNA and bolas can be presumably very similar to that of the siRNAs and the corresponding bolas. The 3′-end of the DNA antisense and the 5′-end of the DNA sense strands were tagged with Alexa488 and IowaBlack fluorescence quencher respectively (see Supplementary Data). Therefore, the fluorescence of Alexa488 is quenched in a DNA duplex due to the proximity of the IowaBlack fluorescence quencher. However, upon digestion of the DNA duplex with DNase, IowaBlack gets separated from Alexa488 restoring its emission (, black line). Fluorescently quenched DNA duplexes (100 nmol/l final) were associated with GLH-19 and GLH-20 (10 µg/ml final) and DNase was added after 2 minutes of incubation at 37 °C. The results demonstrated some protection for the GLH-20/DNA (red line) complex with complete digestion of DNA only after 5 hours, whereas for GLH-19/DNA (blue line), no significant degradation of the DNA duplexes was observed during the same time period. It should be noticed that a significant part of the siRNA helix was exposed to water in the GLH-20/siRNA complex during the MD simulation, whereas the GLH-19 head groups maintained stable interactions with the entire siRNA under the same conditions. Thus, the experimental results are in a good agreement with the computationally predicted considerably higher accessibility of the duplex in the case of the GLH-20 complexes.

Cryogenic electron microscopy imaging and dynamic light scattering. We further studied the bola/siRNA complexes using cryogenic electron microscopy (cryo-EM) and to the best of our knowledge, this was one of the first investigations of siRNA/lipid micelle complexes by this technique. The cryo-EM images () showed that a minority of GLH-19/siRNA complexes (upper image) was buried in thick ice and a majority was located on the carbon grid, whereas the GLH-20/siRNA complexes (lower image) were mostly found on the carbon grid. Interactions of bolas with the carbon surface thus emphasize their hydrophobic properties () with GLH-20/siRNA complexes being more hydrophobic. The size of the GLH-19/siRNA and GLH-20/siRNA complexes varied from 100 nm to 600 nm with the majority of the complexes being ~250 nm in diameter.

Using dynamic light scattering, the diameters of individual micelles were measured to be 135 nm for GLH-19 and 250 nm for GLH-20 (Supplementary Figure S7). The sizes of bola/siRNA complexes were measured to be slightly higher: 224 and 305 nm for GLH-19/siRNA and GLH-20/siRNA, respectively () and are in a good agreement with the cryo-EM observations. The slight increases in size upon siRNA addition can be attributed to the partial growth of the poly-cationic micelles in the presence of multiple negatively charged siRNAs which is consistent with computational predictions (Supplementary Figure S5 and S6).

Transfection and gene silencing efficiencies for bola/siRNA complexes. To study the potential of bolas for siRNA delivery, GLH-19 and GLH-20 were associated with fluorescently labeled (RNA/DNA) duplexes. As fluorophore modified RNA strands are considerably more expensive than DNAs, RNA/DNA hybrids with labeled DNAs were used. The 3′-end of the DNA antisense strands were fluorescently labeled with Alexa488 and duplexed with the sense siRNA strand with a sequence capable of silencing enhanced green fluorescent protein (eGFP). These complexes were then transfected into human breast cancer cells MDA-MB-231 as described in methods. On the next day, transfection efficiencies were visualized by confocal fluorescence microscopy () and statistically analyzed by fluorescence-activated cell sorting, revealing significant intracellular uptake with slightly higher efficiency for the GLH-19 complexes (). The lesser uptake for the GLH-20/siRNA complexes can be explained by the relatively lower binding affinity of GLH-20 which may promote partial dissociation of the complexes in the transfection media thus, preventing some fractions of the siRNAs from entering the cells.

To assess the ability of bola/siRNA complexes to release siRNA inside the cells, experiments with human breast cancer cells stably expressing enhanced green fluorescent protein (MDA-MB-231/eGFP) were carried out (,). Cells were transfected with siRNAs against eGFP[27] complexed with GLH-19 and GLH-20. After 3 days, the levels of eGFP expression were analyzed with fluorescence microscopy () and fluorescence-activated cell sorting (). The results demonstrated significant and comparable levels of silencing in eGFP production for both types of bolas/siRNA complexes, which taken together with the higher levels of intracellular uptake for the GLH-19 complexes, confirmed the MD predictions of lower binding affinity for GLH-20/siRNA, which thus allows for a better release of the siRNA inside the cell.

In vivo experiments. To assess the ability of bolas to be used for siRNA delivery in vivo, experiments with tumor bearing athymic nude mice were carried out (). siRNAs fluorescently labeled with IRDye700 and associated with GLH-19 were systemically delivered to the mouse by tail-vein injections. Next, in vivo siRNA tumor uptake was evaluated by whole body fluorescence imaging. Fluorescent signals of major organs (spleen, liver, kidney, intestines, and gall bladder) were measured and normalized to their weights. Fluorescent signals in the spleen, lungs, heart, and brain were not detected (data not shown). The results demonstrated a relatively high uptake of the siRNA within the tumor compared with the other major organs within the time course of 3 hours thus, suggesting the potential for using bolas for in vivo siRNA delivery.

Discussion

The purpose of this study was to investigate the potential use of bolaamphiphile micelles as a carrier for siRNA delivery. We characterized the interactions between siRNA and two types of bolaamphiphiles, GLH-19 and GLH-20 using MD simulations and various experimental techniques. With the MD simulations, we predicted better protection against nuclease degradation for siRNAs associated with the GLH-19 micelles. The MM-PBSA and MM-GBSA methods predicted a higher binding affinity for the GLH-19/siRNA complex. Consistent with the computational results, gel experiments indicated stronger binding and more stable interactions for the GLH-19/siRNA complexes, which in addition showed almost no degradation in the presence of nucleases. Our cryo-EM studies helped to characterize the bolaamphiphile/siRNA complexes indicating that GLH-19/siRNA and GLH-20/siRNA have different relative hydrophobicities with GLH-19/siRNA being less hydrophobic. On the basis of the cell culture transfections and the in vivo live fluorescence imaging results, we can conclude that both GLH-19 and GLH-20 bolaamphiphiles have a great potential to be used as carriers for siRNA delivery with GLH-19 being a better candidate. Bola/siRNA complexes significantly increase the chemical stability of siRNAs and provide excellent intracellular uptake followed by specific gene silencing. Moreover, depending on the application, the extent of chemical protection of the siRNA can be altered by simply changing the carrier.

Materials and Methods

MD simulations. All MD simulations and experimental methods are detailed in the Supplementary Data. Explicit solvent MD simulations were used to study the interactions between a bola micelle and an siRNA duplex in the presence of 0.15 mol/l NaCl. Sequences for the siRNA duplex in our MD and experimental studies were derived from previous studies.[1,7,31,32,33] The initial structure of the siRNA molecule was equilibrated by 30 nanoseconds MD simulations before it was used to examine the interactions with the bola micelle. Bola micelle structures were prepared as follows. MD simulations of a single GLH-19 and a single GLH-20 bola were performed for 40 nanoseconds, and 20 different snapshots of the structure were taken from the MD trajectories of each bola. The 20 structures of the same bola were placed 5 Å from each other, and used for an additional 50 nanoseconds of MD simulations. During these MD simulations, the 20 bolas gradually formed a micelle structure in which the hydrophobic chains aggregated, whereas the positive head groups were distributed on the micelle surface. These equilibrated GLH-19 and GLH-20 micelles were used to examine the interactions with the siRNA duplex (Supplementary Figure S2). The initial separation distance between the center of mass of the equilibrated bola and the center of mass of the equilibrated siRNA duplex was set to 40 Å. Water and ions were introduced and allowed to be equilibrated with a gradual release of harmonic constraints on the siRNAs and bolas followed by a heating step bringing the system to 300 K (see Supplementary Data for more details). Two additional configurations, one with different initial velocities and another with different initial bola micelle surfaces (achieved by rotating the micelles around the siRNA), were tested for 40 nanoseconds of MD simulations. Also, we examined the interactions between individual bola molecules and RNAs by placing 6, 15, and 24 individual bolas around a 12-mer duplex, mimicking the process of RNA coating with bolaamphiphiles.

More details of the MD simulations are discussed in the Supplementary Data. For all MD simulations, the sizes of the bola micelles (~10 nm) were smaller compared with the experimentally found sizes (~200 nm). This discrepancy is a result of the limitations in computational power and computing time, but it is believed that the simulations are reasonable approximations of larger systems.

Hydrogen bond interactions between the bolas and the siRNA were monitored using 3.5 Å and 120° cutoffs. The population density of the bola head groups and the Cl− ions on the micelle surfaces were calculated by a radial distribution function. The binding free energies of the siRNA-bola complexes were calculated by molecular mechanics MM-PB/GBSA methods.

Synthesis of bolaamphiphiles. The bolas, GLH-19 and GLH-20, were synthesized in Ben Gurion University as previously described[25] and detailed in the Supplementary Data.

Duplex assemblies, bolaamphiphiles/siRNA complex formations, and qualitative agarose gel electrophoresis. Sequences for the siRNA duplex were used from previous studies[1,7,31,32,33] and were designed to silence eGFP. Prepared siRNA-bolas complexes were analyzed at room temperature on 2% agarose gels in the presence of 89 mmol/l Tris-borate, pH 8.3. A Hitachi FMBIO II Multi-View Imager was used to visualize ethidium bromide stained oligos.

Dynamic light scattering experiments. For dynamic light scattering experiments,[34] sample solutions containing either bola micelles prepared RNA/bola complexes were used. The samples were measured at 25 °C with a Zetasizernano (Malvern Instruments, Worcestershire, UK) equipped with a 633 nm laser.

Stability assays with RQ1 DNase fluorescent studies. To study the stability of the preformed duplex/bola complexes in the presence of nucleases, DNA duplexes containing one 3′ antisense strand modified with Alexa488 and a 5′ sense strand modified with IowaBlack FQ (see Supplementary Data) were preincubated with bolas (GLH-19 or GLH-20) and the dequenching of Alexa488 upon digestion with RQ1 DNase (Promega, Madison, WI) was followed by fluorescence resonance energy transfer experiments.

Cryo-EM experiments. Quantifoil Copper 200 mesh R 3.5/1 grids were washed overnight with acetone. To prepare a frozen, hydrated grid, 2.5 µl of sample was applied to the grid, blotted, and plunged into liquid ethane using Vitrobot III (FEI, Hillsboro, OR). Images were collected at liquid nitrogen temperature (~100 K) on a JEM-2010F (JEOL, Tokyo, Japan) transmission electron cryo-microscope equipped with a field emission gun. JEM-2010F was operating at 200 kV and was equipped with a Gatancryo-holder (model 626) (Gatan, Pleasanton, CA). Images were recorded on DE-12, a 12.6 megapixel (3,072 ×4,096) Direct Detection Device sensor (Direct Electron LP, San Diego, CA). Samples were imaged at 13,900× effective magnification targeted at 3–4 µm under focus. We used a total specimen exposure for each image of 30 e−/Å2 second.

Transfection of human breast cancer cells with bolaamphiphile/siRNA complexes. All siRNA transfections in this project were performed with the human breast cancer cell line MDA-MB-231 (with or without eGFP) using bolas (GLH-19 or GLH-20).

Fluorescent light microscopy. To assess the silencing efficiency, cells were imaged after the transfection with a LSM 710 confocal microscope (Carl Zeiss) and a Nikon 200 TE inverted microscope (Nikon, Melville, NJ) respectively.

Flow cytometry experiments. For statistical analysis of the flow cytometry experiments, at least 30,000 MDA-MB-231 cells (with or without eGFP) were analyzed by fluorescence-activated cell sorting analysis.

In vivo experiments. Animal studies were performed according to the Frederick National Laboratory for Cancer Research (Frederick, MD) Animal Care and Use Committee guidelines. Imaging studies were performed on MDA-MB-231 tumor bearing athymic nude mice (Charles River Laboratories, Frederick, MD).

Figure S1. Total occupancy of bola head groups near the siRNA phosphate groups in terms of base pairs for the 30–50 nanoseconds range of the simulation. Figure S2. Structures and nomenclature of siRNA duplex and bolas used. Figure S3. Examples of typical hydrogen bond interactions between hydroxyls of one bola (other bola molecules in each micelle are not shown) and the phosphates of the siRNA backbone for (a) GLH-19 and (b) GLH-20. Figure S4. The number of bola head groups that are associated with the siRNA phosphate groups in the additional MD simulations, which used different velocities (upper panel) and different micelle surfaces (lower panel). Figure S5. MD simulations of individual bolas interaction with RNA duplex. Figure S6. MD simulations of individual bolas interaction with RNA duplex. Figure S7. Size histograms from dynamic light scattering experiments indicate average diameters for GLH-19 (blue) and GLH-20 (red) micelles to be 135 and 246 nm respectively. Table S1. Hydrogen bond interactions between the siRNAs and bolaamphiphiles (bolas). Video S1. Video S2. Data.

Table 1

Results from molecular mechanics Poission-Boltzmann/Generalized Born surface area methods