Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Soft Matter. Author manuscript; available in PMC 2013 January 1.
Published in final edited form as:
PMCID: PMC3375911

Simulation of fusion-mediated nanoemulsion interactions with model lipid bilayers


Perfluorocarbon-based nanoemulsion particles have become promising platforms for the delivery of therapeutic and diagnostic agents to specific target cells in a non-invasive manner. A “contact-facilitated” delivery mechanism has been proposed wherein the emulsifying phospholipid monolayer on the nanoemulsion surface contacts and forms a lipid complex with the outer monolayer of target cell plasma membrane, allowing cargo to diffuse to the surface of target cell. While this mechanism is supported by experimental evidence, its molecular details are unknown. The present study develops a coarse-grained model of nanoemulsion particles that are compatible with the MARTINI force field. Simulations using this coarse-grained model have demonstrated multiple fusion events between the particles and a model vesicular lipid bilayer. The fusion proceeds in the following sequence: dehydration at the interface, close apposition of the particles, protrusion of hydrophobic molecules to the particle surface, transient lipid complex formation, absorption of nanoemulsion into the liposome. The initial monolayer disruption acts as a rate-limiting step and is strongly influenced by particle size as well as by the presence of phospholipids supporting negative spontaneous curvature. The core-forming perfluorocarbons play critical roles in initiating the fusion process by facilitating protrusion of hydrophobic moieties into the interface between the two particles. This study directly supports the hypothesized nanoemulsion delivery mechanism and provides the underlying molecular details that enable engineering of nanoemulsions for a variety of medical applications.


Nanoscale particles, of 200 to 300 nm in diameter, have been developed to sequester loaded cargo until a specific target, recognized by the targeting surface ligands, is reached. The delivery of therapeutic agents including drugs, imaging agents, and macromolecules using nanoscale particles has advantages over conventional small molecule treatment. Delivery using nanoscale particles prevents premature degradation of the therapeutic agents, concentrates the agents at a specific target tissue or cell type, and aids agents in crossing through biological barriers such as the epithelium, endothelium, and plasma membrane [1]. Nanoparticle design must be precisely tailored to deliver each particular agent to its appropriate destination. To achieve efficient delivery to subcellular organelles, the scientific community must understand the molecular mechanisms of how nanoparticles interact with plasma membranes.

Phagocytosis and endocytosis are the most common internalization mechanisms for nanoparticles such as liposomes and polymer-based nanoparticles [1]. However, these internalization mechanisms involve lysosomal degradation that attenuates therapeutic efficiency and as such are not suited for the delivery of therapeutic biomolecules. Alternative strategies have been developed where cargo molecules are directly inserted into the cytoplasmic space through transient pores in the plasma membrane. Direct insertion methods include the conjugation of therapeutic cargo to cell-penetrating peptides [2, 3], electroporation [4, 5], and therapeutic ultrasound with microbubbles [68]. However, all of these methods have the potential to cause cell damage by disrupting the plasma membrane [3, 5, 911].

Another delivery mechanism is provided by perfluorocarbon-based nanoemulsion particles (PFC-NEPs) which are stabilized by an emulsifying phospholipid monolayer. The so called “contact-facilitated” delivery mechanism of PFC-NEPs involves neither lysosomal pathways nor substantial perturbations of the membrane. This delivery mechanism is hypothesized to start with the formation of a hemifusion complex between the monolayer of PFC-NEP and the outer monolayer of target cell plasma membrane. Cargo molecules then diffuse to the plasma membrane through the hemifusion complex and are finally internalized by lipid raft mediated endocytosis [1114]. The contact-facilitated delivery mechanism is particularly useful for the delivery of biomolecules that are highly susceptible to enzymatic reactions. However, the molecular details of this mechanism are as yet undetermined due to experimental difficulties of structure determination of small, fluid, and highly heterogeneous systems.

Molecular dynamics simulations have been widely used to determine membrane structures at both atomistic [1518] and coarse-grained levels [1922]. Atomistic simulations are useful for collecting accurate structural details but are often too costly for examining larger-scale membrane behavior such as self-assembly of lipids into bilayers or vesicles, bilayer phase changes, domain formation, pore formation, and membrane fusion. Therefore, coarse-grained models have been extensively used to simulate these mesoscopic phenomena [1922]. Proposed membrane fusion mechanisms involving a hemifusion stalk [23] have been demonstrated by the appearance of proposed intermediate structures in coarse-grained simulations [24, 25]. Furthermore, such simulations have identified important molecular-scale structural motifs, such as splayed lipids in inducing membrane fusion [26, 27]. Over the last few years, coarse-grained simulations have been adopted to model atomistic structures of biological emulsions such as low density lipoproteins and high density lipoproteins that support these simulations to model nanoemulsion particles and their interaction with bilayer [2830].

This research focuses on an initial step of the contact-facilitated delivery mechanism [14], where two monolayers form a hemifusion complex which resembles the hemifusion stalk [31]. Coarse-grained molecular dynamics simulations have been used to directly observe hemifusion complex formation between a PFC based nanoemulsion particle and a liposome that models a target cell. These simulations were also used to test the dependence of nanoparticle-liposome fusion on particle size and lipid composition in order to identify important structural (physico-chemical) features of the particles required for optimization of fusion.


Development of a coarse-grained PFOB model

To simplify a model system, the MARTINI force field usually maps four atoms into a coarse-grained bead, or interaction site. The MARTINI force field contains four main types of interaction sites: polar (P), non-polar (N), apolar (C), and charged (Q). These main types are subdivided into 18 different subtypes based on their hydrogen bonding capability and level of polarity. Within the 18 different subtypes, 10 different levels of van der Waals interaction potentials are defined: super-attractive, attractive, almost-attractive, semi-attractive, intermediate, almost-intermediate, semi-repulsive, almost-repulsive, repulsive, and super-repulsive [32]. These potentials were obtained by varying the Lennard-Jones energetic well-depth parameter ε from 2.0 to 5.6 kJ/mol with a size parameter σ of 0.47 nm. The only exception is the super-repulsive potential where the size parameter is increased (σ = 0.62 nm) to preserve the strongly unfavorable interaction between charged (Q) and apolar (C) interaction sites. Despite its simplicity, the MARTINI force field has been widely and successfully used to simulate mesoscopic phenomena including membrane fusion, pore formation, self-assembly, and domain formation. Furthermore, the MARTINI force field is easy to implement and is transferable for other small molecules. These issues motivated the development of a coarse-grained model of a perfluorocarbon molecule, perfluorooctylbromide (PFOB), in accordance with the MARTINI force field (see Table 1).

Table 1
PFOB coarse-grained force field parameters.

PFOB was modeled by three interaction sites, named as follows: BRC (BrCF2CF2CF2-), CBM (-CF2-CF2-CF2-), and CBE (-CF2-CF2-CF3) such that the coarse-grained model was a linear 3-particle system, BRC-CBM-CBE. Bond length was determined based on the Lennard-Jones (LJ) size parameters and bending angle was determined based on the existing MARTINI representation of the CG lipid aliphatic chain. To accelerate computation, the MARTINI force field normally uses the same mass for all interaction sites and ignores partial charges [32]. However, to preserve the asymmetric nature of PFOB, true mass and partial charges were used for each site in the force field development and simulations, as discussed in the Discussion section. These three sites were modeled by a new interaction site type, CB, which was developed in this work to preserve the stronger hydrophobic character of perfluorocarbons as compared with hydrocarbons. The initial guess for the LJ parameters (σ = 0.47 nm and ε = 3.5 kJ mol−1) was taken from the polar interaction site type C1 of the MARTINI force field [32]. These parameters were then systematically varied starting from the initial guess to reproduce the density and heat of vaporization of liquid PFOB [3336]. The final optimized CB parameters were σ= 0.5 nm and ε = 3.5 kJ mol−1.

To describe the interaction of PFOB with other molecules such as lipids and water, the LJ parameters between CB and other interaction sites were prepared using the constant size parameter of σ= 0.5 nm and different well depth parameters (ε) based on the parameterized level of attractiveness between the two sites. The slightly greater than standard σ value imposes less favorable interaction level for all pairs than the corresponding pairs with the apolar (C) site. The super-repulsive potential between a charged (Q) and CB was treated as in the MARTINI force field, using a σ value of 0.62 nm.

Force field descriptions for other molecular components

PFOB was modeled by the new coarse-grained parameters described above. Two different phospholipids were used in the simulations: 1-palmitoyl-2-oleolyl-phosphatidylcholine (POPC), and 1-palmitoyl-2-oleolyl-phosphatidylethanolamine (POPE). The parameters for these lipids and water molecules were taken from the standard MARTINI force field [37]. Recently, the MARTINI force field changed the interaction between the charged sites (Q) and apolar sites (C) from repulsive in previous versions to super-repulsive. This change was made to prevent the penetration of ionic sites (Q) through the bilayers [37]. However, this modification also strongly inhibits membrane fusion events that require the protrusion of lipid tails out of the membrane. These events are inhibited in the new force field because of the stronger super-repulsive interaction between apolar lipid tails and the charged head group sites (the choline and the phosphate group) [26]. Therefore, for the current study, interaction values between the charged sites of lipid head group and the lipid tails were reduced from the super-repulsive value to their former repulsive parameters to enable the protrusion of hydrophobic tails in the simulated fusion process. Without this reduction, fusion is inhibited for these lipidic systems, preventing the observation of fusion events on the timescale of available simulation durations [26].

Generation of PFOB-NEP particles and liposomes

To generate a PFOB-NEP droplet, varying numbers of phospholipids, PFOB, and water molecules were randomly mixed in an initial box with dimensions 25.3 nm × 25.3 nm × 25.3 nm. A water shell of 5 nm in thickness surrounded this mixture in all directions so that the hydrophobic molecules were isolated inside the simulation box and would not form a lamellar phase with their periodic images. The system sizes and compositions that were successful in forming an intact PFOB-NEP are listed in Table 2. This initial simulation was extended, following the procedure below, until a PFOB-NEP droplet was formed.

Table 2
Properties of particles tested for fusion. See Figure 1 for geometry definitions.

The protocol of Marrink et al. [21] was used to form the liposome: a small patch was taken from a large bilayer and put in a simulation box. Water molecules were then added into the box so that the small bilayer would be solvated in all directions. The composition of each of these pre-liposome formulations is listed in Table 2. The mean dimensions of the resulting simulation box were 26.5 nm × 26.5 nm × 26.5 nm along the x, y and z directions. This initial simulation was extended, following the procedure below, until a unilamellar liposome was formed.

Molecular dynamics simulations for particle formation and analyses of trajectories were performed with GROMACS version 4.0 [38]. As the first step in the simulations, an energy minimization was performed using a steepest descent method and the system was then gradually heated from 50 K to 303 K through a series of short molecular dynamics (MD) simulations. After these temperature equilibration steps, production simulations were run until the spherical PFOB-NEP and liposomes were formed. A distance cutoff of 1.2 nm was used for non-bonded interactions in combination with a standard shift function in which both the energy and force go to zero at the cutoff distance. The LJ potential is shifted smoothly from 0.9 nm to the cutoff distance of 1.2 nm and the electrostatic potential is shifted from 0.0 nm to the cutoff distance of 1.2 nm. The simulation was performed in an isobaric-isothermal ensemble (NpT) with the Berendsen pressure coupling scheme with a coupling time of 4 ps. The simulation temperature was 303 K via the Berendsen coupling scheme with the coupling time of 0.4 ps [39]. Time steps of 20 fs were used in production simulations.

Liposome:PFOB-NEP fusion simulations

The geometry of each particle tested for fusion is provided in Table 2. To start the fusion simulation, two spherical particles were placed together with less than 1 nm distance between their closest lipids and surrounded by water, as shown in Figure 1. Following the bilayer fusion procedure of Kasson et al.,[24, 40] two lipids, one from each particle, were selected to be connected by an intermediate water site which formed an artificial chemical linker between the water and the phosphate sites of each lipid with an equilibrium length of 0.5 nm. A series of short simulations were performed to gradually reduce the initial linker bond length to the equilibrium length as described below. A list of liposome:PFOB-NEP fusion pairs is provided in Table 4. Five independent replica simulations were run for each fusion pair to enhance sampling efficiency and improve simulation statistics.

Figure 1
The structure of a linked NEP and liposome and the structures of molecules contained in the system. A) The structures of phospholipid (POPC), perfluorooctylbromide (PFOB) and water molecules, from top to bottom, presented at an atomistic (left) and coarse-grained ...
Table 4
Summary of fusion simulations. The fusion initiation modes are denoted by letter with descriptions in Figure 7 and the text.

The same energy minimization and temperature equilibration steps were repeated by applying the simulation parameters described above for particle generation. After temperature equilibration, a series of short simulations of 1 ns length were performed while decreasing the artificial linker length by 0.2 nm after each simulation until the equilibrium length of 0.5 nm was achieved. These short simulations were followed by production simulations, each of which was performed for at most 4 µs. However, if a fusion complex was formed between the two particles before the 4 µs time limit, the simulation was stopped.


Evaluation of the new coarse-grained PFOB model

To determine the bulk thermodynamic parameters including density and heat of vaporization, a simulation containing 559 PFOB molecules and a simulation of a single PFOB (at the equilibrium box size) was carried out for 100 ns at 298 K. These two simulations were used to compute the PFOB heat of vaporization as described previously [41]. As summarized in Table 3, the resulting PFOB density was 1.954 kg L−1 and the heat of vaporization was 44.66 kJ mol−1, values that are in reasonable agreement with the experimental results of 1.925 kg L−1 for density and 42.67 kJ mol−1 for heat of vaporization [3336].

Table 3
Bulk properties of liquid PFOB. The references for each experimental bulk property are Riess et al. [33] and Andre et al. [34] for density (ρ); Riess et al. [33], Gregor et al. [35], Faithful et al. [36] for heat of vaporization (Δvap ...

Additionally, the phase separation of PFOB in water was tested in the absence and in the presence of an emulsifying POPC monolayer. The results were compared to the structures from previous atomistic simulations [41]. The density profiles showed good agreement between the coarse-grained and atomistic simulations, further supporting the utility and applicability of the new coarse-grained model (Figure 2 and Supporting Movie 3).

Figure 2
Interface structure comparison between atomistic and coarse-grained simulations A) in the absence and B) in the presence of an emulsifying phospholipid monolayer. A) The normalized number density profiles of water are shown in dotted red (coarse-grained) ...

Spontaneous formation of liposomes and PFOB-NEP particles

Systems of randomly dispersed PFOB-NEP components were prepared as described above. From these random starting structures, the emulsion formed rapidly: even after accounting for the approximately 4-fold faster time scale of the CG model [37], formation occurs in less than a few hundred nanoseconds. As shown in Figure 3 and Supporting Movie 1, PFOB molecules aggregated while amphipathic phospholipids were expelled to the surface of the growing PFOB particles. This process progressed in several phases and eventually resulted in small monolayer PFOB particles which merged and became a single continuous entity in a toroidal shape. The toroidal hole gradually decreased and PFOB molecules took a spherical shape with phospholipids covering the surface. The structure of the complete PFOB-NEP was defined by its outer radius ro, which is the distance between the center of mass and the peak of phosphate group density, and rc, which is the distance between the center of mass and the peak of terminal methyl carbon sites (see Figure 1). Incomplete PFOB-NEPs containing a water droplet inside the emulsion were occasionally observed. The water droplet was also enclosed by a phospholipid monolayer and, even after significant further simulation (close to 1 µs), it did not disappear. These structures would form when the system contained surplus lipids compared to the amount of PFOB. Therefore, when this was observed, new simulations were run with reduced numbers of phospholipids to obtain complete PFOB-NEPs.

Figure 3
Spontaneous emulsion formation. Snapshots were taken at every 1.5 ns along the trajectory for the first 6 images (A–F). G was taken at 13.5 ns and H was taken at 22 ns, by which time the PFOB-NEP formation was complete. PFOB is shown in green, ...

In addition to the emulsions, liposomes were spontaneously generated from the structures described above [21]. The L2 liposome (see Table 2) used in the simulation was the smallest that could be formed spontaneously from the starting planar bilayer structures: smaller initial structures failed to form liposomes, presumably due to the high curvature stress (data not shown). In addition to the parameters ro and rc, defined above, the structure of a liposome also characterized by 2pz, which is the distance between two phosphate density peaks (see Figure 1).

Liposome:PFOB-NEP fusion

To characterize the fusion between PFOB-NEPs and liposomes prepared by the new coarse-grained simulations, particle size and lipid composition were varied. The tested liposome sizes were 8.4 nm and 7.4 nm in radius (ro). The lipid composition of PFOB-NEP varied from 20 to 50 mol % of POPE with the remainder POPC. The tested combinations of size and composition with relevant characteristics of each fusion process are summarized in Table 4. An example fusion event is shown as an animation in Supporting Movie 2.

Fusion between the liposome and the PFOB-NEP particle was observed in all simulations except one (L2-N4-5) within the simulation time limit of 4 µs. An abnormal fusion event was observed in the L2-N4-2 simulation where a PFOB-NEP spontaneously fused with the periodic image of the linked liposome and subsequent fusion processes were hindered because of the periodic constraints. For this reason, L2-N4-2 fusion simulation was not further analyzed. All other fusion events followed the process presented in Figure 4. These fusion processes started with the merging of the two monolayers followed by PFOB molecules quickly moving into the inter-monolayer space of the fused liposome. The lipids from the PFOB-NEP were mixed into the outer monolayer of the liposome while the PFOB became evenly distributed in the inter-monolayer space. As a result of the fusion, the two particles were integrated into a single liposome with its 2pz distance increased by ~2 nm (see Supporting Figure S2).

Figure 4
Structural changes over the course of fusion. The snapshots are chosen at different time points to show the structural changes over the course of fusion event. Corresponding time points are defined in the following Figure 5A. PFOB is shown in green, the ...

The distance between the center of mass of each particle is denoted by RCM and defined in Figure 1. Figure 5A shows example RCM profiles during a fusion event: RCM fluctuates at a mean value that is equivalent to the sum of the two particle radii for a certain lag time then it drastically decreases. A fusion initiation was determined by the onset of drastic decay in RCM – this period was designated as A. While the lag time for fusion initiation varied as summarized in Table 4, the subsequent changes in RCM were fast and indistinguishable among different fusion events (see Figure 5B). The factors that lead to different fusion initiation are discussed below.

Figure 5
A reaction coordinate used to analyze the fusion process. A) The evolution of RCM, the distance between the center of mass of each particle, is plotted as a function of time. Once the fusion is initiated, RCM rapidly decreases to zero. The inset figure ...

To understand the role of the liposome and nanoemulsion components in initiating fusion, the molecular composition of the interface between the two particles was assessed. To determine the composition, the whole system was transformed by a translation and rotation so that the molecular axis connecting the center of mass of each particle was on the z-axis. Next, the system was divided into thin slabs of 0.1 nm in thickness along the z-axis. To locate the slab at the interface, the number of lipid sites (CG beads) was counted for each slab. The contact interface was identified as the position with the lowest lipid density within the slab. [27]

To determine whether initiation of fusion was accompanied by dehydration at the contact interface, the number of water sites at the interface was counted. For this counting, a distance cutoff of 3 nm from the z-axis was applied so that fluctuations in bulk water would not obscure behavior at the fusion site. Figure 6 shows the changes in the number of waters between the monolayers at the fusion site. The number of waters between the monolayers dropped from 19 ± 4 to < 1 when fusion started. This drastic decrease in the number of waters suggested that the hydrophobic exposure described in the next section is dependent upon dehydration at the interface.

Figure 6
Changes in molecular composition at the interface. A) The evolution of the number of water (red) and each hydrophobic moiety (blue: lipid tails of liposome, cyan: lipid tails of PFOB-NEP, green:PFOB) at the interface is plotted with the fusion initiation ...

As shown in Figure 6A, the appearance of the hydrophobic moieties followed dehydration at the interface. Exposure of hydrophobic moieties to the membrane surface is known to be a rate limiting step of membrane fusion [26]. To explore the participation of hydrophobic molecules in initiating fusion, the appearance of lipid tails and PFOB was assessed by counting the CG beads of each molecule at the interface. The lipids of the PFOB-NEP and the liposome were evaluated separately. The sequence of hydrophobic moiety appearance differed between fusion events. Figure 6B shows the frequency with which each hydrophobic moiety appeared first, second, and third at the surface. PFOB was most frequently observed as the first hydrophobic moiety at the interface followed by the lipid tails of the PFOB-NEP. Additional details of fusion are presented in Supporting Information Figure S1.

Using this analysis, 23 fusion events were categorized into four distinct modes (labeled A–D in subsequent discussion) by the order of exposure of hydrophobic groups. For this classification, hydrophobic exposure events were considered to be simultaneous if they occurred within 0.3 ns of each other, thus setting a lower time limit for the mechanistic analysis. The modes assigned to each fusion simulation are presented in Table 4 and example structures of each mode are shown in Figure 7. In mode A, which was observed only in one fusion event (L1-N1-3), PFOBs first appeared at the interface with the exposure of lipid tails following after 0.4 ns. In mode B, observed in 5 fusion events (L1-N1-2, L1-N1-4, L2-N3-3, L2-N3-5, L2-N4-1, and L2-N4-4), PFOB and lipid tails from the PFOB-NEP first appeared at the interface, followed by the appearance of liposome lipid tails after 0.3 to 0.5 ns. In mode C, observed in 4 fusion events (L1-N1-1, L2-N2-1, L2-N3-1, and L2-N3-2), lipid tails from both PFOB-NEP and liposome were first observed at the interface while the appearance of PFOB was delayed by 0.5 to 0.9 ns. Finally, in mode D, which included the remaining 13 fusion events, all three hydrophobic moieties appeared at the interface simultaneously.

Figure 7
Different fusion initiation modes. A) Primary action of PFOB prior to the protrusion of lipids. B) Protrusion of PFOB-NEP lipid in conjunction with PFOB molecules (solid red) C) Protrusion of lipids in the absence of PFOB molecules. D) Interaction among ...

Dependence of fusion lag-time on lipid composition

The fusion lag time dependency on lipid composition was examined while the particle size was kept constant. To test the influence of negative intrinsic curvature, the fraction of POPE was decreased from 50 to 20 mol %. Figure 8B indicates that the fraction of POPE was negatively correlated with the lag time. The lag time differences of L2-N4 from L2-N1 and L2-N2 were statistically significant at a 95% confidence level. A similar analysis of fusion rate on particle size (Figure 8A) was performed and but the results were not statistically significant at a 95% confidence level.

Figure 8
Fusion dependence on particle size and lipid composition. A) The lag times for the fusion initiation are compared among fusion simulations that differ in comprising particle size at a constant lipid composition of POPC:POPE=1:1. B) The lag times are compared ...

Order parameter dependence on particle size and lipid composition

Order parameters of lipid tails were examined to determine if particle size and different lipid composition affected the flexibility of lipid tail conformation in the NEP as addressed for liposomes [42]. The deuterium order parameters SCD were computed with the following relation [43]:


where Sij is an element of an order parameter tensor S and is defined by


for i, j = x, y, z, where θi is the angle between the ith molecular axis and the sphere normal that is defined by a vector from a center of mass of the particle to the phosphate group of each lipid. The angular bracket implies averaging over time and molecules. The molecular axes are defined per CG bead along the tail. For the nth bead (Cn), the z axis was defined by a vector from Cn−1 to Cn+1, the y axis by a vector perpendicular to z and in the plane formed by Cn−1, Cn, and Cn+1, and the x axis by a vector perpendicular to y and z.

The deuterium order parameters were computed for each CG bead along the lipid tails of the liposome (L1, L2) and the PFOB-NEP (N1, N2, N3, N4) and are shown in Figure 9. The order parameters of palmitoyl (left) and oleoyl chain (right) are plotted separately. It was observed that greater particle sizes were positively correlated with greater order parameters, which indicated that conformational freedom was reduced as the particle size increased. Even though decreasing the POPE fraction did not affect the tail order substantially, the order parameters of oleoyl chain (C2 and C3), were significantly greater in the particles containing reduced amount of POPE as shown in Figure 10.

Figure 9
Lipid tail order parameters of liposomes and PFOB-NEPs. The lipid tail orders of sn-1 palmitoyl chain (A) and sn-2 oleoyl chain (B) of liposomes are plotted. The lipid tail orders of sn-1 palmitoyl chain (C) and sn-2 oleoyl chain (D) of PFOB-NEPs are ...
Figure 10
Statistically significant differences in tail orders among PFOB-NEPs. A) The tail orders at the C2 position (A) and C3 position of the oleoyl chain are plotted. The differences of the designated pairs are statistically significant at a 95% confidence ...


To better understand the “contact-facilitated” delivery mechanism proposed for NEP-membrane interactions, force field parameters were developed at a coarse-grained level, generated spherical PFOB nanoemulsions and liposomes using equilibrium simulations, and ran fusion simulations after physically linking a PFOB-NEP and a liposome. 24 fusion events were observed across different combinations of PFOB-NEPs and liposomes, which strongly supports the contact-facilitated delivery mechanism. To the best of our knowledge, this is the first demonstration of fusion between liposomes and nanoemulsion particles enclosed by a phospholipid monolayer. This process shared some features with the well-established fusion process between two bilayers but also exhibited distinctive features that are discussed below.

Force field parameters

The newly developed coarse-grained PFOB model was able to reproduce bulk thermodynamic properties of PFOB in reasonable agreement with experimental observations. The emulsion-solvent interfacial structures determined with the coarse-grained model showed good agreement with the structure determined by atomistic simulations both in the absence and presence of an emulsifying phospholipid monolayer. Additionally, the spontaneous formation of PFOB emulsion droplets enclosed by a phospholipid monolayer starting from random mixtures of PFOB, phospholipids, and water further supported the validity of this new model. Even though non-zero partial charges were assigned to the three interaction sites to preserve the asymmetric nature of PFOB, PFOB behaved symmetrically without showing preferential orientation near the PFC-NEP surface (data not shown) as observed in atomistic simulations [44]. The small partial charges were not influential under these simplified simulation conditions. Hence, in the future, simulation performance of PFOB-NEPs could be improved by excluding these charges with no penalty to the simulation accuracy.

System configurations and approximations

To observe fusion events within the limit of available computing resources, a few constraints were imposed on the system configuration.

First, small particles were used, which has been common for vesicle fusion simulations [21, 24, 27, 40, 45]. Small size can both speed up the simulation and facilitate fusion between membranes due to pronounced curvature. While the smallest liposome size of ~15 nm in diameter corresponded to the minimum liposome size feasible in experiments [46], the smallest PFOB-NEP size ~13 nm in diameter was far below the experimentally-observed minimum PFC-NEP size of ~50 nm in diameter [47]. These small sized particles emphasized the effect of local curvature on the surface PFC-NEPs and target cell membranes to initiate fusion. Indeed, it was recently reported that synaptotagmin fusion protein could induce high local positive curvature and reduce the energy required to initiate membrane fusion [48, 49].

Second, the lipid:PFOB molecule ratios for PFOB-NEP formation were determined empirically under the condition that initial random mixtures of certain ratio would not form water droplets inside the emulsions. This condition could have resulted in dilute lipid packing on the surface that may affect the fusion event. Therefore, future studies should introduce lipids in water so that they spontaneously inserted into the preformed NEP surface to overcome this limitation.

Third, the NEP and liposome particles were placed close together and were connected by an artificial linker that maintained the initial close proximity between the two particles. Such linkage was necessary because spontaneous apposition of the two particles is a rare event beyond the current timescales of molecular simulation. Such encounter is entropically disfavored and energetically hindered due to electrostatic repulsion and dehydration penalties but can be facilitated by additional molecules such as ligand-receptor interactions [50]; e.g., to target cell-specific surface molecules [51]. In this work, these ligand-receptor interactions were mimicked by an artificial linker, as used in many recent fusion simulations [24, 27, 40].

Mechanisms for NEP-liposome complex formation

The 23 successful fusion events revealed four different modes for the disruption of the monolayers, as shown in Figure 7 and listed in Table 4. Three of the modes (A, B, and D) showed that PFOB played a primary role in initiating membrane disruption and lipid complex formation. In mode A, PFOB leaked prior to the reorientation of lipids into the interface and subsequently triggered membrane disruption. In modes B and D, PFOB also facilitated the reorientation and protrusion of lipid tails by solvating the hydrophobic tails exposed to surface. However, the boundaries between these modes were somewhat indistinct and, in any mode, the interplay between phospholipids and PFOBs was prevalent. This again emphasizes the critical role of PFOB in the initiation of fusion.

The “contact-facilitated” delivery mechanism originally hypothesized that the NEP and liposome particles would stay intact with a hemifusion complex spanning the gap and providing a continuous hydrophobic surface between the two particles [14]. This hypothesis was based on experimental observations [9, 12, 50, 52, 53] as well as on the lipophobic properties of PFOB that suggested it would maintain the emulsion core even after a hemifusion stalk formed. In contrast to this hypothesis, these simulations showed that hemifusion complexes are very short-lived, quickly expanding radially to allow the massive flow of PFOB into the intermonolayer space, finally leading to the complete absorption of PFOB-NEP into the apposed liposome.

There are a few possible reasons for this difference between simulation and experiment. First, the small particle sizes used in the simulations generate high curvature strain and surface tension that strongly drives the system to minimize the surface area for a given volume. In these simulations, disappearance of the hemifusion stalk through complete absorption is a likely mechanism to maximize such reduction in curvature and surface area. Such high curvature stress is much less likely when particles of larger, and more realistic, sizes are used. Second, the current systems are pure and lack proteins. The flow of PFOB from the NEP into the liposome separates the two monolayers and thickens the membrane. In plasma membranes containing high fractions of transmembrane protein, thickening of membrane would not be as energetically favorable because of increased hydrophobic mismatches with membrane proteins [54]. Finally, these simulations used a single flexible linker to keep the two particles in close proximity while the close apposition of PFC-NEPs to their target cells are mediated by multivalent ligand-receptor interactions with multiple copies of homing liands on PFC-NEP surface [55]. Such multivalent ligand-receptor interactions should keep PFC-NEPs from being absorbed into the target cells and thus further support the stability of the hemifusion stalk.

Fusion dependence on lipid composition

Lipids of negative spontaneous curvature are known to facilitate membrane fusion between the two outer monolayers because they can accommodate the highly negative curvature that occurs during hemifusion stalk formation [56]. This dependency was observed in the fusion simulations: lower fractions of POPE (~20 mol %) prolonged the lag time for the initiation of fusion compared to the PFOB-NEPs containing higher POPE fractions (50 and 40 mol %). The lag time trend with POPE fraction was statistically significant. In particular fusion was not observed between L2 and N4 in one simulation (L2-N4-5) with a low fraction of POPE molecules, until the simulation reached the 4 µs maximum time limit of available computing resources. Correlated increases in order parameters with increasing liposome size and decreasing POPE proportion help to explain the cause of the prolonged lag time, which is further discussed below.

Correlation of lipid tail order with fusion mechanism and rate

As expected, lipids were more ordered in the larger liposome (L1) than in the smaller one (L2). On the other hand, even though the area per head group was the same among PFOB-NEPs (N1 through N4), the higher fraction of POPE caused the lipid tails to become more disordered.

The lag time prolongation was more strongly correlated with increased tail order of PFOB-NEPs (4.1-fold longer) rather than to tail order of the liposome (1.7-fold longer). This observation suggests that membrane disruption in the PFOB-NEP monolayer may act as a rate-limiting step which, in turn, triggers liposome outer monolayer disruption. This speculation is supported by the fact that the hydrophobic moieties of the PFOB-NEP, including lipid tails and PFOB, were more frequently found as the first molecule at the NEP-liposome interface – rather than the lipid tails of liposome.


This paper presents a new coarse-grained PFOB model that preserves the thermodynamic properties of pure PFOBs as well as the phase separation in water due to its strong hydrophobicity. The validity of the model was also confirmed by the spontaneous emulsion formations starting from randomly mixed initial structures.

The reported simulations using this new CG model directly showed that fusion could occur between the monolayer of PFOB-NEP and the outer monolayer of liposome. This observation strongly supports the “contact-facilitated” delivery mechanism, which was hypothesized to proceed by forming an intermediate structure similar to the hemifusion stalk. Indeed, this fusion event between a PFOB-NEP and a liposome shared many aspects in common with the fusion mechanisms for two bilayers. The NEP-liposome fusion showed a similar dependence on the particle size and lipid composition: smaller sizes of higher negative-curvature lipid content resulted in faster fusion. Fusion was initiated by forming a hemifusion complex made of hydrophobic moieties protruded to the surface. Initial membrane disruption and hemifusion complex formation were the rate-limiting steps while subsequent steps occurred very quickly and barely differed among the observed fusion events.

This fusion process between a PFOB-NEP and a liposome has several distinctive features. It appeared that freely mobile PFOB molecules played critical roles in facilitating the fusion process. PFOB enhanced the protrusion of hydrophobic moieties from the PFOB-NEP to the interface region by solvating lipid tails. Finally it was observed that the flow of PFOB into the inter-monolayer space of the apposed liposome led to the complete absorption of a PFOB-NEP into the connected liposome.

This work has laid a foundation for new studies to understand the delivery mechanism of cargo molecules to target cell through membrane fusion. Systematic analysis to test the fusion dependency on particle size and lipid dependency as well as core-forming perfluorocarbon molecules would help to better characterize the mechanisms and to aid in the rational design of particles for optimal delivery.

Supplementary Material

Supporting Information

Supporting Movie 1

Supporting Movie 2

Supporting Movie 3


The authors would like to thank Michael Daily for his help with manuscript preparation. This work was sponsored by NIH grant U54 CA11934205 to SAW, NAB, and GML, and NIH grant U01 NS073457-01 to NAB and GML. Computational resources were provided by the Texas Advanced Computing Center through Teragrid Grants TG-MCB060053 and TG-MCA08X003 as well as the National Biomedical Computation Resource (NIH P41 RR0860516).


1. Hillaireau H, Couvreur P. Nanocarriers’ entry into the cell: relevance to drug delivery. Cellular and Molecular Life Sciences. 2009;66(17):2873–2896. [PubMed]
2. Rinne J, et al. Internalization of novel non-viral vector TAT-streptavidin into human cells. BMC Biotechnology. 2007;7(1):1. [PMC free article] [PubMed]
3. Kerkis A, et al. Properties of cell penetrating peptides (CPPs) IUBMB life. 2006;58(1):7–13. [PubMed]
4. Fei Z, et al. Gene Transfection of Mammalian Cells Using Membrane Sandwich Electroporation. Analytical Chemistry. 2007;79(15):5719–5722. [PubMed]
5. Ionescu-Zanetti C, Blatz A, Khine M. Electrophoresis-assisted single-cell electroporation for efficient intracellular delivery. Biomedical Microdevices. 2008;10(1):113–116. [PubMed]
6. Guzman HR, et al. Ultrasound-mediated disruption of cell membranes. I. Quantification of molecular uptake and cell viability. The Journal of the Acoustical Society of America. 2001;110(1):588–596. [PubMed]
7. Taniyama Y, et al. Local Delivery of Plasmid DNA Into Rat Carotid Artery Using Ultrasound. Circulation. 2002;105(10):1233–1239. [PubMed]
8. Feinstein SB, et al. The powerful microbubble: from bench to bedside, from intravascular indicator to therapeutic delivery system, and beyond. American Journal of Physiology - Heart and Circulatory Physiology. 2004;287(2):H450–H457. [PubMed]
9. Soman NR, et al. Acoustic activation of targeted liquid perfluorocarbon nanoparticles does not compromise endothelial integrity. IEEE Transactions on Nanobioscience. 2006;5(2):69–75. [PubMed]
10. Miller DL, Quddus J. Diagnostic ultrasound activation of contrast agent gas bodies induces capillary rupture in mice. Proceedings of the National Academy of Sciences. 2000;97(18):10179–10184. [PubMed]
11. Wible JH, Jr, et al. Microbubbles induce renal hemorrhage when exposed to diagnostic ultrasound in anesthetized rats. Ultrasound in Medicine &amp; Biology. 2002;28(11–12):1535–1546. [PubMed]
12. Soman NR, et al. New mechanisms for non-porative ultrasound stimulation of cargo delivery to cell cytosol with targeted perfluorocarbon nanoparticles. Nanotechnology. 2008;19(18):185102–185102. [PMC free article] [PubMed]
13. Kaneda MM, et al. Mechanisms of nucleotide trafficking during siRNA delivery to endothelial cells using perfluorocarbon nanoemulsions. Biomaterials. 2010;31(11):3079–3086. [PMC free article] [PubMed]
14. Partlow KC, Lanza GM, Wickline SA. Exploiting lipid raft transport with membrane targeted nanoparticles: A strategy for cytosolic drug delivery. Biomaterials. 2008;29(23) p.-. [PMC free article] [PubMed]
15. Lee S-J, Song Y, Baker NA. Molecular dynamics simulations of asymmetric NaCl and KCl solutions separated by phosphatidylcholine bilayers: potential drops and structural changes induced by strong Na+-lipid interactions and finite size effects. Biophysical Journal. 2008;94(9):3565–3576. [PubMed]
16. Marrink SJ, et al. Simulation of the Spontaneous Aggregation of Phospholipids into Bilayers. Journal of the American Chemical Society. 2001;123(35):8638–8639. [PubMed]
17. Olsen BN, Schlesinger PH, Baker NA. Perturbations of membrane structure by cholesterol and cholesterol derivatives are determined by sterol orientation. Journal of the American Chemical Society. 2009;131(13):4854–4865. [PMC free article] [PubMed]
18. Song Y, Guallar V, Baker NA. Molecular dynamics simulations of salicylate effects on the micro- and mesoscopic properties of a dipalmitoylphosphatidylcholine bilayer. Biochemistry. 2005;44(41):13425–13438. [PMC free article] [PubMed]
19. Marrink SJ, de Vries AH, Mark AE. Coarse Grained Model for Semiquantitative Lipid Simulations. Journal of Physical Chemistry B. 2004;108(2):750–760.
20. Shelley JC, et al. A Coarse Grain Model for Phospholipid Simulations. The Journal of Physical Chemistry B. 2001;105(19):4464–4470.
21. Marrink SJ, Mark AE. Molecular Dynamics Simulation of the Formation, Structure, and Dynamics of Small Phospholipid Vesicles. Journal of the American Chemical Society. 2003;125(49):15233–15242. [PubMed]
22. Izvekov S, Voth GA. A Multiscale Coarse-Graining Method for Biomolecular Systems. Journal of Physical Chemistry B. 2005;109(7):2469–2473. [PubMed]
23. Kozlovsky Y, Kozlov MM. Stalk Model of Membrane Fusion: Solution of Energy Crisis. Biophysical Journal. 2002;82(2):882–895. [PubMed]
24. Kasson PM, et al. Ensemble molecular dynamics yields submillisecond kinetics and intermediates of membrane fusion. Proceedings of the National Academy of Sciences of the United States of America. 2006;103(32):11916–11921. [PubMed]
25. Smeijers AF, et al. A detailed look at vesicle fusion. The journal of physical chemistry. B. 2006;110(26):13212–13219. [PubMed]
26. Mirjanian D, et al. Splaying of aliphatic tails plays a central role in barrier crossing during liposome fusion. J Phys Chem B. 2010;114(34):11061–11068. [PMC free article] [PubMed]
27. Kasson PM, Lindahl E, Pande VS. Atomic-Resolution Simulations Predict a Transition State for Vesicle Fusion Defined by Contact of a Few Lipid Tails. PLoS Comput Biol. 2010;6(6):e1000829. [PMC free article] [PubMed]
28. Murtola T, et al. Low density lipoprotein: structure, dynamics, and interactions of apoB-100 with lipids. Soft Matter. 7(18)
29. Vuorela T, et al. Role of Lipids in Spheroidal High Density Lipoproteins. PLoS Comput Biol. 6(10):e1000964. [PMC free article] [PubMed]
30. Henner√© Gl, et al. Molecular dynamics study of a phospholipid monolayer at a water/triglyceride interface: towards lipid emulsion modelling. Chem. Phys. Lipids. 2009;157(2):86–93. [PubMed]
31. Chernomordik LV, Kozlov MM. Membrane Hemifusion: Crossing a Chasm in Two Leaps. Cell. 2005;123(3):375–382. [PubMed]
32. Marrink SJ, et al. The MARTINI Force Field:[APL funtional symbol quad question] Coarse Grained Model for Biomolecular Simulations. The Journal of Physical Chemistry B. 2007;111(27):7812–7824. [PubMed]
33. Riess JG. Oxygen Carriers (“Blood Substitutes”)Raison d'Etre, Chemistry, and Some Physiology Blut ist ein ganz besondrer Saft 1. Chemical Reviews. 2001;101(9):2797–2920. [PubMed]
34. Andre M, Nelson T, Mattrey R. Physical and Acoustical Properties of Perfluorooctylbromide, An Ultrasound Contrast Agent. Investigative Radiology. 1990;25(9):983–987. [PubMed]
35. Gregor T, et al. Aerosolization of perfluorocarbons during mechanical ventilation: an in vitro study. Intensive Care Medicine. 2003;29(8):1354–1360. [PubMed]
36. Faithfull NS, Weers JG. USPTO, Editor. United States; 1997. Partial liquid breathing of fluorocarbons.
37. Marrink S, et al. The MARTINI Force Field:[APL funtional symbol quad question] Coarse Grained Model for Biomolecular Simulations. The Journal of Physical Chemistry B. 2007;111(27):7812–7824. [PubMed]
38. Hess B, et al. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. Journal of Chemical Theory and Computation. 2008;4(3):435–447. [PubMed]
39. Berendsen HJC, et al. Molecular dynamics with coupling to an external bath. The Journal of Chemical Physics. 1984;81(8):3684–3690.
40. Kasson P, Pande VS. Control of membrane fusion mechanism by lipid composition: predictions from ensemble molecular dynamics. PLoS computational biology. 2005;(2007):e220–e220. preprint. [PMC free article] [PubMed]
41. Lee S-J, et al. Characterization of Perfluorooctylbromide-Based Nanoemulsion Particles Using Atomistic Molecular Dynamics Simulations. The Journal of Physical Chemistry B. 2010;114(31):10086–10096. [PMC free article] [PubMed]
42. Risselada J, Marrink S. Curvature effects on lipid packing and dynamics in liposomes revealed by coarse grained molecular dynamics simulations. Phys. Chem. Chem. Phys. 2009;11(12):2056–2067. [PubMed]
43. Egberts E, Berendsen HJC. Molecular dynamics simulation of a smectic liquid crystal with atomic detail. Journal of Chemical Physics. 1988;89(6):3718–3732.
44. Lee S-J, et al. Characterization of Perfluorooctylbromide-Based Nanoemulsion Particles Using Atomistic Molecular Dynamics Simulations. The Journal of Physical Chemistry B. 2010;114(31):10086–10096. [PMC free article] [PubMed]
45. Knecht V, Marrink S-J. Molecular Dynamics Simulations of Lipid Vesicle Fusion in Atomic Detail. Biophysical Journal. 2007;92(12):4254–4261. [PubMed]
46. Brouillette CG, et al. Minimal size phosphatidylcholine vesicles: effects of radius of curvature on head group packing and conformation. Biochemistry. 1982;21(19):4569–4575. [PubMed]
47. Keipert PE, et al. Influence of Perflubron Emulsion Particle Size on Blood Half-Life and Febrile Response in Rats. Artificial Cells, Blood Substitutes and Biotechnology. 1994;22(4):1169–1174. [PubMed]
48. Martens S, Kozlov MM, McMahon HT. How Synaptotagmin Promotes Membrane Fusion. Science. 2007;316(5828):1205–1208. [PubMed]
49. Hui E, et al. Synaptotagmin-Mediated Bending of the Target Membrane Is a Critical Step in Ca2+-Regulated Fusion. Cell. 2009;138(4):709–721. [PMC free article] [PubMed]
50. Soman NR, et al. Synthesis and Characterization of Stable Fluorocarbon Nanostructures as Drug Delivery Vehicles for Cytolytic Peptides. Nano Lett. 2008 p.-. [PMC free article] [PubMed]
51. Crowder K, et al. Sonic activation of molecularly-targeted nanoparticles accelerates transmembrane lipid delivery to cancer cells through contact-mediated mechanisms: implications for enhanced local drug delivery. Ultrasound in medicine & biology. 2005;31(12):1693–1700. [PubMed]
52. Soman NR, et al. Molecularly targeted nanocarriers deliver the cytolytic peptide melittin specifically to tumor cells in mice, reducing tumor growth. The Journal of clinical investigation. 2009;119(9):2830–2842. [PMC free article] [PubMed]
53. Soman NR, et al. Noncavitational mechanisms of interaction of ultrasound with targeted perfluorocarbon nanoparticles: implications for drug delivery. Vol. 3. IEEE; 2005. pp. 1712–1715.
54. Antoinette KJ. Hydrophobic mismatch between proteins and lipids in membranes. Biochimica et Biophysica Acta (BBA) - Reviews on Biomembranes. 1998;137(6–3):401–416. [PubMed]
55. Kaneda M, et al. Perfluorocarbon Nanoemulsions for Quantitative Molecular Imaging and Targeted Therapeutics. Annals of Biomedical Engineering. 2009;37(10):1922–1933. [PMC free article] [PubMed]
56. Chernomordik LV, et al. The shape of lipid molecules and monolayer membrane fusion. Biochimica et Biophysica Acta (BBA) - Biomembranes. 1985;812(3):643–655.
57. Song W, Rossky PJ, Maroncelli M. Modeling alkane + perfluoroalkane interactions using all-atom potentials: Failure of the usual combining rules. Journal of Chemical Physics. 2003;119(17):9145–9162.