Nanoparticle-Engendered Rupture of Lipid Membranes - The Journal

Aug 7, 2018 - Tension-induced rupture of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid membranes with encapsulated hydrophobic ...
0 downloads 0 Views 5MB Size
Letter pubs.acs.org/JPCL

Cite This: J. Phys. Chem. Lett. 2018, 9, 4872−4877

Nanoparticle-Engendered Rupture of Lipid Membranes Sean Burgess,† Aleksey Vishnyakov,‡ Christopher Tsovko,† and Alexander V. Neimark*,† †

Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854, United States ‡ Skolkovo Institute of Science and Technology, Nobel st. 1, 121205 Moscow, Russia

J. Phys. Chem. Lett. 2018.9:4872-4877. Downloaded from pubs.acs.org by UNIV OF SUNDERLAND on 10/09/18. For personal use only.

S Supporting Information *

ABSTRACT: Tension-induced rupture of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid membranes with encapsulated hydrophobic nanoparticles is elucidated using dissipative particle dynamics simulations. The dynamics of hole formation is studied, and a nanoparticle size-dependent relationship is established for the probability of membrane rupture within a given time as a function of the membrane tension. Two mechanisms of hole formation are explored: homogeneous nucleation and heterogeneous nucleation at the nanoparticle surface. While the kinetics of homogeneous nucleation in unloaded membranes complies with the predictions of the classical Deryagin−Gutop (DG) theory, the heterogeneous nucleation causes progressively lower lysis tensions as the particle size increases. The thermodynamics of heterogeneous nucleation is treated by introducing an effective contact angle at the three-phase, solvent−membrane−solid boundary into the DG theory. The proposed approach helps quantitatively interpret the simulation results and predict the membrane stability in real experiments with significantly larger (by many orders of magnitude) observation times and spatial scales.

I

assessment of cell membrane damage by commercial nanoparticles. Lipid-based membranes are prone to rupture under stress or being affected by external stimuli. The basic thermodynamic theory of membrane rupture was suggested by Deryagin and Gutop (DG)6 based on the classical nucleation theory (CNT) applied to a homogeneous two-dimensional film. The DG model implies that film rupture is triggered by nucleation, due to thermal fluctuations, of an unstable hole (pore), which grows spontaneously upon achieving a certain critical size. The hole nucleation is controlled by the membrane surface tension, σ, and the line tension, Γ, which represents the excess free energy of the membrane edge per unit length. The work of formation of a round hole of radius r equals the difference between the work of formation of the hole edge of length 2πr and the work of extension of the membrane area by πr2

nteractions of nanoparticles (NPs) with lipid bilayers (LBs), which constitute the foundation of cell membranes, play an important role in emerging nanotechnologies, such as drug delivery, biomedical imaging, and design of bioinspired membranes for water purification.1−3 At the same time, an explosion in production of nanomaterials and their commercial applications raises multiple health-related concerns regarding nanoparticle toxicity.4 Adhesion of inhaled or digested nanoparticles may lead to membrane rupture and cell apoptosis.5 A better understanding of the nanoparticleengendered mechanisms of membrane instabilities is required for developing new nanoparticle biotechnologies as well as for evaluating health threats related to nanoparticle manufacturing. In this Letter, we explore tension-induced rupture of lipid membranes with encapsulated nanoparticles using dissipative particle dynamics (DPD) simulations. To compensate for the lack of relevant experimental data, we produced a set of data on the dynamics of hole nucleation in 1,2-dimyristoyl-snglycero-3-phosphocholine (DMPC) membranes with encapsulated NPs of different sizes by using an original simulation setup, which mimics experiments performed with unloaded membranes. To explain the found correlation between the probability of rupture and the particle size, we developed a thermodynamic theory of the heterogeneous nucleation by introducing a notion of the contact angle at the three-phase, solvent−membrane−solid, boundary controlled by the balance of respective line tensions. These findings shed light on the mechanism of heterogeneous nucleation of holes in membranes triggered by adhered or encapsulated nanoparticles that may have important implications both for the therapeutical use of nanoparticles as drug and imaging carriers and for the © 2018 American Chemical Society

E(r , σ ) = 2πr Γ − πr 2σ

(1)

E(r, σ) achieves a maximum at rc = Γ/σ, which represents the size of critical hole nucleus: holes of r < rc may collapse, while holes of r > rc expand irreversibly. The work of formation of the critical nucleus Ec = E(rc , σ ) = πrc 2σ = π Γ 2/σ

(2)

represents the free energy, or nucleation (activation) barrier, of expanding hole formation, which is the main quantitative parameter for predicting the probability and kinetics of film Received: May 31, 2018 Accepted: August 7, 2018 Published: August 7, 2018 4872

DOI: 10.1021/acs.jpclett.8b01696 J. Phys. Chem. Lett. 2018, 9, 4872−4877

Letter

The Journal of Physical Chemistry Letters

silica. Ganzenmuller et al.39 explored shock-wave rupture of lipid bilayers. Several papers considered relevant phenomena of NP translocation 34,38,42,43 and wrapping by lipid bilayers;32,33,40,44 the latter problem is important for analyses of stability of NP dispersions.45 Yue et al.32 and Zhang et al.33 employed DPD simulations to study membrane rupture triggered by rotating NPs brought into contact with LBs. They determined a critical surface tension, above which the membrane ruptures. Although these studies shed light on the mechanisms of NP−membrane interactions, the problem of how particle encapsulation affects stability of lipid membranes under tension remained unexplored, and a robust theoretical model capable of predicting the dependence of the lysis tension on the NP size and chemistry is lacking. In this work, we analyze the stability of NP-loaded membranes in silico and treat the heterogeneous nucleation in terms of an original theoretical model by introducing an effective contact angle on the three-phase, hole−membrane− solid boundary into the DG theory. To elucidate the specifics of tension-induced rupture of lipid membranes with encapsulated nanoparticles, we implement standard DPD simulation method in the form46 suggested by Groot and Warren.47 As an instructive example, we study the dynamics of hole formation in DMPC lipid bilayers with loaded hydrophobic nanoparticles ranging from 1 to 4 nm in radius upon incremental membrane stretching; the simulation details are presented in the Supporting Information, section I. In DPD simulations, we mimic the classical experiments with a film placed of a frame with a movable plank, to which a constant force is applied to control the membrane tension (Figure 1); the details of the simulation setup are given in

rupture. The nucleation barrier determines the rate constant of film rupture k f = A f,0 exp( −Ec /kBT ) = A f,0 exp(−π Γ 2/σkBT )

(3)

Af,0 is the kinetic rate prefactor that depends on the film area and lipid mobility.7 It is worth noting that the DG model is purely macroscopic and does not account for the molecular features of the lipid interfaces that are important at the nanoscale level of critical nuclei. Nevertheless, the DG model in its original and modified forms has been widely used for interpretation of the results of experimental8−14 as well as molecular simulation7,15−20 studies of the lysis tension defined as characteristic tension of membrane rupture at given conditions. In particular, Levadny et al.11 studied the tension-driven rupture of giant unilamellar vesicles (GUVs) and identified the lysis tension by the micropipette suction method.21 By repeating experiments, the authors determined the probability Pf(σ, t) of the membrane rupture within given time t at fixed tension σ. The rupture probability Pf(σ, t) as a function of σ was correlated with the predictions of the classical Kramers dynamics theory22 that implies the Poissonian stochastic process of hole formation with the rate constant kf determined by the energy barrier via eq 3 Pf (σ , t ) = 1 − exp( −k f t )

(4)

From this correlation, the authors calculated the line tension, Γ. This approach with various modifications was applied for interpretation of experimental studies of homogeneous lipid membranes of different type. Rupture of unloaded homogeneous membranes has been studied by various theoretical and molecular simulation methods15−20,23,24 attempting to identify the critical nucleation conditions and to calculate the membrane line tension. The problem of transformation of the results of molecular simulations for predicting the experimental conditions of membrane rupture remains open. The reported values of the lysis tension tend to be larger in simulated systems (σ = 20− 40 pN/nm)16,20,24 when compared to experimental systems (σ = 5−10 pN/nm)11,12,25−27 because of significant, by many orders of magnitude, difference in observation time and system size. Because of a strong correlation between the line and surface tensions,12,19 the values of the line tension determined in simulations (Γ = 12−40 pN) are larger than those in experiments (Γ = 4−28 pN). While the homogeneous hole nucleation in unloaded lipid membranes was studied in great detail,9−12,16,24−28 we are not aware of any experimental studies of heterogeneous hole nucleation in lipid membranes engendered by adhered or encapsulated NPs, which is the focus of our work. Various simulation studies monitored the process of particle-assisted membrane rupture,29−31 yet very few consider the NP effect on the nucleation barriers of hole formation.32−34 By using SCFT, Ting and Wang34 calculated free energy barriers for translocation of hydrophobic NPs across LBs. They showed different mechanisms of NP translocation: NP encapsulation followed by NP release and membrane relaxation and NP encapsulation followed by membrane rupture. It is worth noting that DPD simulations have been extensively used and found efficient for modeling various processes involving lipid membranes and nanoparticles.35−41 Doi et al.37 investigated the specifics of water interlayers between lipid membrane and

Figure 1. Snapshots (side and top views) illustrating the simulation setup for monitoring the dynamics of the tension-driven membrane rupture. Lipid membrane is placed on a frame with a movable plank to which a constant force is applied to control the membrane tension. The force, and respectively the membrane tension, is incrementally increased until the hole nucleation event occurs and the membrane ruptures at a certain lysis tension. Left: membrane with a spherical NP of 2.15 nm in radius. Right: membrane loaded with cylindrical tabletlike NP of 3.15 nm in radius and 2.80 nm in height. Bead colors: brown, lipid tails; pink, lipid junction; white, lipid head; red, NP; dark green, inner static bar; cyan, outer static bar; purple, inner pulling bar; light green, outer pulling bar; water not shown.

Supporting Information, section I. Upon creation of the equilibrated initial configuration, the applied force, Fp, is increased in small increments causing the respective increase of the membrane tension, σ. The increments of the applied force are chosen to be small enough to secure the independence of the time distribution of the lysis tension of the rate of tension increase. At each incremental step, statistics are collected over 4873

DOI: 10.1021/acs.jpclett.8b01696 J. Phys. Chem. Lett. 2018, 9, 4872−4877

Letter

The Journal of Physical Chemistry Letters

Figure 2. Example of heterogeneous nucleation at a 3.15 nm spherical NP: (a) subcritical embryonic hole, (b) hole of nearly critical size of rc = 0.77 nm, (c) expanding overcritical hole, (d) large hole at the verge of membrane rupture. Lipid beads above NP are transparent.

Figure 3. Simulation data (symbols) and theoretical dependence of the rupture probability Pr(σ, t)|t=176ns as a function of the applied tension σ for DMPC membrane loaded with spherical (left) and cylindrical (right) particles of different radii. The error bars are calculated using the normal approximation interval with a 95% confidence interval. For a 1 nm NP, no effects of heterogeneous nucleation are found and the theoretical line (magenta) overlaps with the theoretical line for unloaded membrane (green).

Figure 4. Snapshot (left) of a hole formed at the NP surface. The hole filled by water has a quasi-circular shape similar to a 2D bubble residing on a curved solid surface. 2D schematics of this configuration (right) illustrates the definition of the contact angle θ at the three-phase solvent−lipid− solid points of contact (points A and A*). The value of cos θ is determined from the balance of separate line tensions on lipid−solvent, Γls; particle−solvent, Γps; and particle−lipid, Γpl, boundaries according to eq 5. Bead colors: brown, lipid tails; pink, lipid junction; white, lipid head; red, NP; dark green, inner static bar; cyan, outer static bar; purple, inner pulling bar; light green, outer pulling bar; water not shown.

100 independent repeats with 200 000 steps at each σ to determine the probability, Pr(t, σ), of membrane rupture within given time t at fixed tension σ (membrane rupture probability). A typical example of the evolution of heterogeneous hole nucleation is shown in Figure 2. The results of simulations with spherical and cylindrical particles of different size are shown in Figure 3. The right-most green data in the left image of Figure 3 corresponds to the homogeneous nucleation in an unloaded membrane. The green solid line represents the predictions of the DG model according to eqs 3 and 4 with fitted

membrane−solvent line tension Γfs= 38.7 pN and dynamic prefactor Af,0 = 9.4 × 1018 /s. The lysis tension σ* (defined as the tension at which the rupture probability Pr(σ, t) = 0.5) within given observation time of 176 ns (which is equivalent to 200 000 DPD time steps), is σ* = 41 pN/nm. The obtained value of the line tension is within the range of reported simulation data.15,16,24 As seen from Figure 3, encapsulated NPs affect the membrane rupture and the lysis tension decreases as the NP size increases. The smallest NPs of 1 nm radius, which are completely immersed within the hydrophobic interior of the 4874

DOI: 10.1021/acs.jpclett.8b01696 J. Phys. Chem. Lett. 2018, 9, 4872−4877

Letter

The Journal of Physical Chemistry Letters membrane, do not affect the dynamics of hole nucleation significantly. The events of heterogeneous nucleation near the NP were not observed. In contrast, in the case of NPs with radius 4 nm (diameter twice as large as the membrane thickness), all nucleation events occurred near the NP surface. In the case of 2 and 3 nm nanoparticles, both homogeneous and heterogeneous nucleation events were observed, but with larger probability for heterogeneous nucleation. The statistics of homogeneous and heterogeneous nucleation events are presented in the Supporting Information, section II. A typical example of heterogeneous nucleation is shown in the left-hand panel of Figure 4. The hole filled by solvent (not shown) is of a quasi-circular shape resembling a 2D bubble residing on the curved NP surface with a certain contact angle at the solid−solvent−lipid points of contact. To consider the heterogeneous nucleation in the spirit of the DG theory, we present the hole as a circle with excluded lens due to the intersection with the NP, as sketched in the right-hand panel of Figure 4. The mechanical equilibrium in this three-phase, twodimensional model is determined by the balance of the line tensions Γls, Γps, and Γpl of lipid−solvent (blue arc AB*A*), particle−solvent (red arc AA*), and particle−lipid (black arc) boundaries, respectively. The balance of line tensions defines the effective contact angle, θ, at the three-phase solvent−lipid− solid points of contact (points A and A*) via the Young equation: cos(θ) = (Γps − Γpl)/Γls

k p = A p exp( −Ep,c /kBT )

(7)

The proposed extension of the DG model to heterogeneous nucleation of holes in LB at the surface of encapsulated NPs is purely macroscopic and can be justified only for large particles of size significantly exceeding the LB thickness. However, as shown below it provides a useful insight on the heterogeneous mechanism of membrane rupture on the same level of accuracy as the original DG model for the homogeneous nucleation. In the case of NP-loaded membranes, the rupture may be initiated by either heterogeneous or homogeneous nucleation, and the rupture probability Pr(t,σ) is determined with the sum of the respective rates kp (eq 3) and kf (eq 7) Pr(σ , t ) = 1 − exp( −(k f + k p)t )

(8)

with the prefactors depending on the system size and NP size. To determine the contact angle θ, we applied the proposed model of heterogeneous nucleation to reproduce the simulation data for a 4 nm NP because in this case no events of homogeneous nucleation were observed. The blue lines in Figure 3 represent the predictions of the expanded DG model for the 4 nm NP systems, calculated according to eqs 6 and 7 with fitted contact angle of θ = 61° and dynamic prefactors Ap,0 = 9.43 × 1012 /s for spherical and Ap,0 = 1.24 × 1012 /s for cylindrical NPs. The calculated nucleation barrier for spherical NP equals Ep,c = 14.1 kT with respective critical hole radius rc = 0.87 nm at the lysis tension of 33.8 pN/nm. The calculated energy barrier for cylindrical NP is somewhat smaller, Ep,c = 12.7 kT, with a critical hole radius of 0.79 nm at the lysis tension of 36.7 pN/nm. Note that we assume the membrane− solvent line tension Γfs = 38.7 pN determined for the homogeneous nucleation and ignore its potential dependence on the membrane tension σ. In the case of 2 and 3 nm NPs, we have to account for both mechanisms of homogeneous and heterogeneous nucleation according to eq 8. The respective theoretical predictions are shown by solid lines in Figure 3. The calculations are performed without using any adjustable parameters using the values of line tension and contact angle determined above. However, the dynamic prefactors are adjusted to account for the system geometry. The prefactor Ap of the rate of the heterogeneous nucleation is scaled with respect to the NP radius Rp as

(5)

This 2D geometrical model allows one to calculate the change of the system free energy, Ep(r, Rp, Γls, θ, σ), due to the formation of a hole of radius r at the NP of radius Rp at the constant tension conditions, which depends on the lipid− solvent line tension, Γls; contact angle, θ; and given tension, σ: Ep(r , R p , Γls, θ , σ ) = [2πr − lABA*]Γls + lAB*A*(Γps − Γpl) − [πr 2 − AAB*A*B]σ = [2πr − lABA* + lAB*A* cos(θ )]Γls − [πr 2 − AAB*A*B]σ (6)

Here, lABA* and lAB*A* are the lengths of arcs ABA* and AB*A*, respectively, and AAB*A*B is the area of lens AB*A*B, all of which depend on r,Rp, and θ (detailed equations are given in the Supporting Information, section III). The first term in the left-hand side (LHS) of the first line of eq 6 corresponds to the work of formation of the hole−lipid boundary of length 2πr − lABA*; the second term corresponds to the difference of line energy of solid−lipid and solid− solvent boundaries of length ABA; the third term corresponds to the surface energy gained because of formation of the hole of area πr2 − AAB*A*B. The second line of eq 6 utilizes the definition of the contact angle (eq 5). Equation 6 presents the work of formation of the hole of radius r at the surface of NP of radius R in LB with given Γls, θ, and σ. Similar to the DG relationship (eq 1), this dependence is nonmonotonic with a maximum equaled to the energy barrier Ep,c of the heterogeneous nucleation achieved at the critical hole radius rc, which depends on Rp, Γls, θ, and σ. An example of such dependence is given in the Supporting Information (Figure S4). The energy barrier (eq 6) determines the rate of heterogeneous nucleation in the same fashion as the DG model does in the case of homogeneous nucleation.

A p = A p,0R p/R p,0

(9)

where Ap,0 is the prefactor for the NP of radius Rp,0 that is taken as the reference (here, Rp,0 = 4 nm). In the case of spherical NPs, the NP radius is additionally adjusted to account for the NP curvature, as described in the Supporting Information, section V. The rate of homogeneous nucleation prefactor, Af,0, is adjusted by subtracting from the membrane area the area occupied by NP where the homogeneous nucleation cannot be initiated A f = A f,0(S − πR p2)/S

(10)

where S is the total surface area of the membrane and Af,0 is the prefactor for the unloaded membrane (eq 3). The calculated energy barriers, lysis tension, and critical radii for spherical and cylindrical NPs are presented in the Supporting Information, section IV. Overall, the agreement with the simulation data shown in Figure 3 confirms that the proposed model, despite its 4875

DOI: 10.1021/acs.jpclett.8b01696 J. Phys. Chem. Lett. 2018, 9, 4872−4877

The Journal of Physical Chemistry Letters



simplistic 2D representation of the NP−LB interface, captures the specifics of the heterogeneous hole nucleation and the mechanisms of rupture of NP-loaded membranes. The introduced contact angle provides a rationale for quantitative assessments of NP−LB interactions. It can be estimated from the data of in silico or in vitro experiments on tension-driven membrane rupture based on the extended DG model presented here. It is worth noting that in order to use the simulation results to predict experiments, it is necessary to account for the difference by many decades in the system size and the observation time. Because of this difference, as well as an elevated level of environmental fluctuations, the lysis tension observed in experiments is significantly smaller than in simulations. The reported data for the lysis tension in different types of unloaded pure experimental lipid membranes varies in the range of 5−10 pN/nm, and the respective range of calculated line tension lies within 4−28 pN.11,12,25−27 A simple estimate can be done based on eq 4 by adjusting the parameters determined in simulation taking into account the experimental membrane size and observation time. Assuming that the prefactor in eq 3 is proportional to the membrane area, Aexp = AsimSexp /Ssim and using the characteristic membrane area Sexp = 380 μm2 (compared with simulated Ssim = 500 nm2) and the experimental observation time of texp = 360 sec (compared with simulated tsim = 176 ns), as reported in ref 11, the predicted experimental lysis tension is σ* = 17.7 pN/nm, which is significantly lower than the value σ* = 41 pN/nm determined in simulations yet still higher than the experimental result. Another factor affecting the reduction of the lysis tension is the decrease of the line tension with the decrease of the surface tension,12,19 and further improvement of the proposed model requires a relationship between the line and surface tensions. For example, if we assume a smaller value of the line tension of 24.5 pN (instead of 38.7 pN determined in simulation), the lysis tension reduces to σ* = 7.0 pN/nm, which is within the experimentally observed range. Dissipative particle dynamics is used to simulate the lipid bilayer system under tension. All systems are run in a box size of 30 × 60 × 25 Rc (20 × 39 × 16 nm), bead density of 3 beads/Rc3. A time step of 0.01τ is chosen with τ = 88 ps. The velocity Verlet algorithm is implemented in an NVT-ensemble, setting T = 298 K. A 3-water-per-bead coarse graining with the bead size Rc = 0.65 nm is implemented. Calculations are performed using DL-MESO software48 run on the Extreme Science and Engineering Discovery Environment (XSEDE).49 To create the initial configuration, 1170 lipid molecules are placed between the frame bars in the bath of 120 000 water molecules, and the system is equilibrated under zero surface tension for one million steps (880 ns) to form a stable homogeneous LB. The LB surface tension is increased in increments of 5 pN/nm at a rate of 0.0228 pN/nm/ns, to a point where preliminary scanning suggests ripping may occur. At this point, the increments are decreased to 1 pN/nm at a rate of 0.0057 pN/nm/ns. All data-gathering points are run for 200 000 steps (176 ns). The DPD model and further simulation details are presented in the Supporting Information, section 1. System images are visualized using Visual Molecular Dynamics.50

Letter

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.8b01696. Dissipative particle dynamics description, system setup, DMPC model justification, interaction parameters, nanoparticle radii, heterogeneous rupture preference, derivation for expanded DG theory, energy barrier of hole formation, calculated lysis tensions and energy barriers for all systems studied, modifications for spherical NPs, tension loading rate justification, and list of abbreviations (PDF) Homogeneous rupture simulation (MPG) Heterogeneous (3 nm) rupture simulation (MPG)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Sean Burgess: 0000-0001-8209-9662 Aleksey Vishnyakov: 0000-0002-1621-3858 Alexander V. Neimark: 0000-0002-3443-0389 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors thank the National Science Foundation for providing support through the CBET Grant No. 1264702 “Adhesion and Translocation of Nanoparticles through Lipid Membranes” and computing resources through the Extreme Science and Engineering Discovery Environment (XSEDE).



REFERENCES

(1) Kumar Teli, M.; Mutalik, S.; Rajanikant, G. K. Nanotechnology and Nanomedicine: Going Small Means Aiming Big. Curr. Pharm. Des. 2010, 16, 1882−1892. (2) Rosenholm, J.; Sahlgren, C.; Linden, M. Cancer-Cell Targeting and Cell-Specific Delivery by Mesoporous Silica Nanoparticles. J. Mater. Chem. 2010, 20, 2707−2713. (3) Gonzalez-Perez, A.; Persson, K. M. Bioinspired Materials for Water Purification. Materials 2016, 9, 447−458. (4) Nel, A.; Xia, T.; Madler, L.; Li, N. Toxic Potential of Materials at the Nanolevel. Science 2006, 311, 622−627. (5) Nel, A. E.; Madler, L.; Velegol, D.; Xia, T.; Hoek, E. M. V.; Somasundaran, P.; Klaessig, F.; Castranova, V.; Thompson, M. Understanding Biophysicochemical Interactions at the Nano-Bio Interface. Nat. Mater. 2009, 8, 543−557. (6) Deryagin, B. V.; Gutop, Y. V. Theory of the Breakdown (Rupture) of Free Films. Kolloidn Zh. 1962, 24, 370−374. (7) Akimov, S. A.; Volynsky, P. E.; Galimzyanov, T. R.; Kuzmin, P. I.; Pavlov, K. V.; Batishchev, O. V. Pore Formation in Lipid Membrane I: Continuous Reversible Trajectory from Intact Bilayer Through Hydrophobic Defect to Transversal Pore. Sci. Rep. 2017, 7, 12152. (8) Karal, M. A. S.; Yamazaki, M. Communication: Activation Energy of Tension-Induced Pore Formation in Lipid Membranes. J. Chem. Phys. 2015, 143, 081103. (9) Karal, M. A. S.; Levadnyy, V.; Tsuboi, T. A.; Belaya, M.; Yamazaki, M. Electrostatic Interaction Effects on Tension-Induced Pore Formation in Lipid Membranes. Phys. Rev. E 2015, 92, 012708. (10) Evans, E.; Smith, B. A. Kinetics of Hole Nucleation in Biomembrane Rupture. New J. Phys. 2011, 13, 095010.

4876

DOI: 10.1021/acs.jpclett.8b01696 J. Phys. Chem. Lett. 2018, 9, 4872−4877

Letter

The Journal of Physical Chemistry Letters

Molecular Dynamics Study. ACS Appl. Mater. Interfaces 2017, 9, 26665−26673. (34) Ting, C. L.; Wang, Z. G. Minimum Free Energy Paths for a Nanoparticle Crossing the Lipid Membrane. Soft Matter 2012, 8, 12066−12071. (35) Espanol, P.; Warren, P. B. Perspective: Dissipative Particle Dynamics. J. Chem. Phys. 2017, 146, 150901. (36) Guigas, G.; Morozova, D.; Weiss, M. Exploring Membrane and Protein Dynamics with Dissipative Particle Dynamics. Adv. Protein Chem. Struct. Biol. 2011, 85, 143−182. (37) Doi, H.; Okuwaki, K.; Mochizuki, Y.; Ozawa, T. A New Treatment for Water Near the Interface Between Lipid Membrane and Silica in Dissipative Particle Dynamics Simulation. J. Comput. Chem., Jpn. 2017, 16, 28−31. (38) Zhang, L. Y.; Wang, X. Q. Coarse-Grained Modeling of Vesicle Responses to Active Rotational Nanoparticles. Nanoscale 2015, 7, 13458−13467. (39) Ganzenmuller, G. C.; Hiermaier, S.; Steinhauser, M. O. ShockWave Induced Damage in Lipid Bilayers: A Dissipative Particle Dynamics Simulation Study. Soft Matter 2011, 7, 4307−4317. (40) Yue, T. T.; Li, S. X.; Xu, Y.; Zhang, X. R.; Huang, F. Interplay Between Nanoparticle Wrapping and Clustering of Inner Anchored Membrane Proteins. J. Phys. Chem. B 2016, 120, 11000−11009. (41) Zhang, L. Y.; Wang, X. Q. Nanotube-Enabled Vesicle-Vesicle Communication: A Computational Model. J. Phys. Chem. Lett. 2015, 6, 2530−2537. (42) Su, C. F.; Merlitz, H.; Rabbel, H.; Sommer, J. U. Nanoparticles of Various Degrees of Hydrophobicity Interacting with Lipid Membranes. J. Phys. Chem. Lett. 2017, 8, 4069−4076. (43) Gkeka, P.; Sarkisov, L.; Angelikopoulos, P. Homogeneous Hydrophobic-Hydrophilic Surface Patterns Enhance Permeation of Nanoparticles Through Lipid Membranes. J. Phys. Chem. Lett. 2013, 4, 1907−1912. (44) Rossi, G.; Barnoud, J.; Monticelli, L. Polystyrene Nanoparticles Perturb Lipid Membranes. J. Phys. Chem. Lett. 2014, 5, 241−246. (45) Alipour, E.; Halverson, D.; McWhirter, S.; Walker, G. C. Phospholipid Bilayers: Stability and Encapsulation of Nanoparticles. Annu. Rev. Phys. Chem. 2017, 68, 261−283. (46) Hoogerbrugge, P. J.; Koelman, J. Simulating Microscopic Hydrodynamic Phenomena with Dissipative Particle Dynamics. EPL 1992, 19, 155−160. (47) Groot, R. D.; Warren, P. B. Dissipative Particle Dynamics: Bridging the Gap Between Atomistic and Mesoscopic Simulation. J. Chem. Phys. 1997, 107, 4423−4435. (48) Seaton, M. A.; Anderson, R. L.; Metz, S.; Smith, W. DL_MESO: Highly Scalable Mesoscale Simulations. Mol. Simul. 2013, 39, 796−821. (49) Towns, J.; Cockerill, T.; Dahan, M.; Foster, I.; Gaither, K.; Grimshaw, A.; Hazlewood, V.; Lathrop, S.; Lifka, D.; Peterson, G. D.; et al. XSEDE: Accelerating Scientific Discovery. Comput. Sci. Eng. 2014, 16, 62−74. (50) Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33−38.

(11) Levadny, V.; Tsuboi, T.; Belaya, M.; Yamazaki, M. Rate Constant of Tension-Induced Pore Formation in Lipid Membranes. Langmuir 2013, 29, 3848−3852. (12) Ly, H. V.; Longo, M. L. The Influence of Short-Chain Alcohols on Interfacial Tension, Mechanical Properties, Area/Molecule, and Permeability of Fluid Lipid Bilayers. Biophys. J. 2004, 87, 1013−1033. (13) Alam Shibly, S. U.; Ghatak, C.; Karal, M. A. S.; Moniruzzaman, M.; Yamazaki, M. Experimental Estimation of Membrane Tension Induced by Osmotic Pressure. Biophys. J. 2016, 111, 2190−2201. (14) Karal, M. A. S.; Levadnyy, V.; Yamazaki, M. Analysis of Constant Tension-Induced Rupture of Lipid Membranes Using Activation Energy. Phys. Chem. Chem. Phys. 2016, 18, 13487−13495. (15) Tolpekina, T. V.; den Otter, W. K.; Briels, W. J. Nucleation Free Energy of Pore Formation in an Amphiphilic Bilayer Studied by Molecular Dynamics Simulations. J. Chem. Phys. 2004, 121, 12060− 12066. (16) Tolpekina, T. V.; den Otter, W. K.; Briels, W. J. Simulations of Stable Pores in Membranes: System Size Dependence and Line Tension. J. Chem. Phys. 2004, 121, 8014−8020. (17) Leontiadou, H.; Mark, A. E.; Marrink, S. J. Molecular Dynamics Simulations of Hydrophilic Pores in Lipid Bilayers. Biophys. J. 2004, 86, 2156−2164. (18) Wohlert, J.; den Otter, W. K.; Edholm, O.; Briels, W. J. Free Energy of a Trans-Membrane Pore Calculated from Atomistic Molecular Dynamics Simulations. J. Chem. Phys. 2006, 124, 154905. (19) Akimov, S. A.; Volynsky, P. E.; Galimzyanov, T. R.; Kuzmin, P. I.; Pavlov, K. V.; Batishchev, O. V. Pore Formation in Lipid Membrane II: Energy Landscape Under External Stress. Sci. Rep. 2017, 7, 12509. (20) Ting, C. L.; Appelo, D.; Wang, Z. G. Minimum Energy Path to Membrane Pore Formation and Rupture. Phys. Rev. Lett. 2011, 106, 168101. (21) Hochmuth, R. M. Micropipette Aspiration of Living Cells. J. Biomech. 2000, 33, 15−22. (22) Hanggi, P.; Talkner, P.; Borkovec, M. Reaction-Rate Theory Fifty Years After Kramers. Rev. Mod. Phys. 1990, 62, 251−341. (23) Zhdanov, V. P. Mechanism of Rupture of Single Adsorbed Vesicles. Chem. Phys. Lett. 2015, 641, 20−22. (24) Jiang, F. Y.; Bouret, Y.; Kindt, J. T. Molecular Dynamics Simulations of the Lipid Bilayer Edge. Biophys. J. 2004, 87, 182−192. (25) Portet, T.; Dimova, R. A New Method for Measuring Edge Tensions and Stability of Lipid Bilayers: Effect of Membrane Composition. Biophys. J. 2010, 99, 3264−3273. (26) Zhelev, D. V.; Needham, D. Tension-Stabilized Pores in Giant Vesicles - Determination of Pore-Size and Pore Line Tension. Biochim. Biophys. Acta, Biomembr. 1993, 1147, 89−104. (27) Kim, D. H.; Frangos, J. A. Effects of Amyloid Beta-Peptides on the Lysis Tension of Lipid Bilayer Vesicles Containing Oxysterols. Biophys. J. 2008, 95, 620−628. (28) Guigas, G.; Morozova, D.; Weiss, M. Exploring Membrane and Protein Dynamics with Dissipative Particle Dynamics. In Advances in Protein Chemistry and Structural Biology, Vol. 85: Computational Chemistry Methods in Structural Biology; Christov, C., Ed.; Elsevier Academic Press Inc: San Diego, CA, 2011; pp 143−182. (29) Alexeev, A.; Uspal, W. E.; Balazs, A. C. Harnessing Janus Nanoparticles to Create Controllable Pores in Membranes. ACS Nano 2008, 2, 1117−1122. (30) Ginzburg, V. V.; Balijepalli, S. Modeling the Thermodynamics of the Interaction of Nanoparticles with Cell Membranes. Nano Lett. 2007, 7, 3716−3722. (31) Vacha, R.; Martinez-Veracoechea, F. J.; Frenkel, D. ReceptorMediated Endocytosis of Nanoparticles of Various Shapes. Nano Lett. 2011, 11, 5391−5395. (32) Yue, T. T.; Zhang, X. R.; Huang, F. Molecular Modeling of Membrane Responses to the Adsorption of Rotating Nanoparticles: Promoted Cell Uptake and Mechanical Membrane Rupture. Soft Matter 2015, 11, 456−465. (33) Zhang, L. Y.; Zhao, Y. P.; Wang, X. Q. Nanoparticle-Mediated Mechanical Destruction of Cell Membranes: A Coarse-Grained 4877

DOI: 10.1021/acs.jpclett.8b01696 J. Phys. Chem. Lett. 2018, 9, 4872−4877