The Role of Lipid Coating in Transport of ... - ACS Publications

Subscriber access provided by TUFTS UNIV

Biological and Environmental Phenomena at the Interface

The Role of Lipid Coating in Transport of Nanodroplets across the Pulmonary Surfactant Layer Revealed by Molecular Dynamics Simulations Yan Xu, Shixin Li, Zhen Luo, Hao Ren, Xianren Zhang, Fang Huang, Yi Y. Zuo, and Tongtao Yue Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.8b01547 • Publication Date (Web): 09 Jul 2018 Downloaded from http://pubs.acs.org on July 10, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

The Role of Lipid Coating in Transport of Nanodroplets across the Pulmonary Surfactant Layer Revealed by Molecular Dynamics Simulations

Yan Xu,† Shixin Li,† Zhen Luo,† Hao Ren,† Xianren Zhang,‡ Fang Huang,† Yi Y. Zuo,§,|| and Tongtao Yue*,†

†

State Key Laboratory of Heavy Oil Processing, Center for Bioengineering and Biotechnology,

College of Chemical Engineering, China University of Petroleum (East China), Qingdao 266580, China ‡

State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical

Technology, Beijing 100029, China §

Department of Mechanical Engineering, University of Hawaii at Manoa, Honolulu, Hawaii

96822, United States ||

Department of Pediatrics, John A. Burns School of Medicine, University of Hawaii,

Honolulu, Hawaii 96826, United States

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT: Hydrophilic drugs can be delivered into lungs via nebulization for both local and systemic therapies. Once inhaled, ultrafine nanodroplets preferentially deposit in the alveolar region where they first interact with the pulmonary surfactant (PS) layer, with nature of the interaction determining both efficiency of the pulmonary drug delivery and extent of the PS perturbation. Here, we demonstrate by molecular dynamics simulations the transport of nanodroplets across the PS layer being improved by lipid coating. In the absence of lipids, bare nanodroplets deposit at the PS layer to release drugs that can be directly translocated across the PS layer. The translocation is quicker under higher surface tensions, but at the cost of opening pores that disrupt the ultrastructure of the PS layer. When the PS layer is compressed to lower surface tensions, the nanodroplet prompts collapse of the PS layer to induce severe PS perturbation. By coating the nanodroplet with lipids, the disturbance of the nanodroplet on the PS layer can be reduced. Moreover, the lipid-coated nanodroplet can be readily wrapped by the PS layer to form vesicular structures, which are expected to fuse with the cell membrane to release drugs into secondary organs. Properties of drug bioavailability, controlled drug release and enzymatic tolerance in real systems could be improved by the lipid coating on nanodroplets. Our results provide useful guidelines for molecular design of nanodroplets as carriers for pulmonary drug delivery.


Page 2 of 25

Page 3 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

INTRODUCTION Inhalation therapy has emerged for decades and already been used in clinical medicine.1 It appears to be a desirable route of administration for treatment of many respiratory disorders, such as asthma, pneumonia and cystic fibrosis.2 Such a ‘topical’ or targeted treatment with inhaled drugs allows the use of smaller doses that reduce the systemic side effect.3 Compared with other drug delivery routes, the human lung provides numerous advantages, including a large absorptive surface area, thin alveolar-capillary barrier and abundant underlying vasculature.4 Therefore, developing effective strategies to promote pulmonary drug delivery has attracted remarkable scientific interests in the past two decades.5 However, such a progress has been hindered by the lack of knowledge on how the inhlaed drugs or drug carriers interact with different biological entities in the pulmonary system. Hydrophilic drugs, such as insulin, can be delivered into lungs via nebulization to achieve efficient alveoli deposition.6 Nebulizers use compressed gas or ultrasound to break up drug solutions into ultrafine aerosol nanodroplets that can be directly inhaled from the mouthpiece of the device.7 The therapeutic efficacy is largely dependent on the drug stability, bioavailability, targeting and uptake efficiency, and biological activity.8 Once inhaled, nanodroplets need to cross the pulmonary surfactent (PS) layer to achieve pulmonary entry. The PS layer is a complex material lining at the air-water interface of lung alveoli. It is composed of approximately 90% lipids (phospholipids and cholesterol) and 10% proteins (SP-A, SP-B, SP-C and SP-D), and plays dual roles in tension reduction and innate host defense.9-14 After transport across the PS layer, nanodroplets enter the alveolar fluid and get ready to act on other biological entities. Obviously, the interaction of inhaled nanodroplets with the PS represents the initial nanobio interaction occurring in the lungs and critically influences the subsequent fate of drugs. Besides, both structures and biophysical functions of the endogenous PS could be perturbed by interacting with nanodroplets. There exist numerous experimental evidences convincingly showing that interactions with inhaled nanoparticles can adversely affect biophysical functions of the natural PS.15-17 Therefore, designing better nanodroplets as carriers for promoting the drug delivery efficiency and reducing the PS perturbation is in an urgent need. In spite of the importance of interactions between inhaled nanodroplets and the PS, relevant studies at the molecular level are still lacking. Here, we employ a model nanodroplet representing nebulized hydrophilic drugs to investigate interactions with the PS layer using molecular dynamics (MD) simulations. We also designed a nanodroplet coated with lipids that


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

make the nanodroplet surface more hydrophobic. Simulation results showed that bare nanodroplets rapidly adhered onto the PS layer to release drugs diffusing across the layer. The trans-layer diffusion was enhanced under higher surface tensions, but at the cost of opening pores that disrupt the ultrastrucutre of the PS layer. When the PS layer was under compression, the nanodroplet prompted collapse of the PS layer. By coating the nanodroplet with lipids, the PS perturbation by nanodroplet was reduced. Moreover, the lipid-coated nanodroplet can be readily wrapped by the PS layer to form vesicular structures, which may enhance fusion with cell membranes to improve drug release into secondary organs. Our simulations suggest that premixing drug solutions with lipids could let nebulizers generate lipid-coated nanodroplets which may improve the efficiency of pulmonary drug delivery and reduce the PS perturbation.

MODEL AND METHODOLOGY Coarse-grained Models and System Setup. The entire system was simulated using the coarse-grained (CG) models, which allow MD simulations to be conducted in adequate temporal and spatial scales. Among CG methods,18 the MARTINI force field,19 which maps four heavy atoms into one interactive bead, has been widely used to simulate biomolecules, such as lipids,20 proteins,21 and DNA,22 and especially the nano-bio interactions.23-24 The demonstrated ability to reproduce biological phenomena makes it an adequate tool to simulate the nanodroplet-PS interactions. The system setup includes a water slab in vacuum with two symmetric PS monolayers at the two air-water interfaces (Fig. 1a). Each layer contains 2,176 dipalmitoylphosphatidylcholine (DPPC), 704 palmitoyloleoylphosphatidylglycerol (POPG) and 704 cholesterol molecules, doped by 13 SP-B and 13 SP-C proteins (Fig. 1b). Both types of proteins were derived from their all-atom models in the protein data bank (PDB IDs: 2DWF (Mini-B); 1SPF (SP-C)). The aqueous phase contains 238,858 water beads and 1148 Na+ ions to neutralize the system. The bare nanodroplet was constructed by stacking 11,690 hydrophilic beads into a spherical shape of 14 nm, in which the relative position of each bead was derived from an equilibrated water box to ensure liquid phase of the droplet (Fig. 1c). To construct the lipid-coated nanodroplet, 1,150 DPPC molecules were physically adsorbed onto the nanodroplet (Fig. 1d). It accordingly has a larger diameter of 16 nm. The choice of DPPC for lipid coating is because it is the major component of the natrual PS, thus having higher biocompatibility. Note that the specific structure can be generated in reality via nebulizing a drug solution mixed with lipids.25 To simulate interactions between nanodroplets and the cell membrane, a lipid bilayer system containing 8,192 POPC molecules and 1,275,104 water beads was also prepared (Fig. 1e). This model has been widely used by other researchers to


Page 4 of 25

Page 5 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

represent the alveolar epithelial cell membrane.26 Simulation Details. Both bare and lipid-coated nanodroplets were equilibrated in air for 10 ns before positioned above the pre-equilibrated PS layer at the air-water interface. After energy minimization using the steepest descent algorithm, the CGMD simulations with constant particle number, surface tension and temperature were then carried out. For all simulations, a cutoff of 1.2 nm was used for van der Waals interactions. The Lennard-Jones potential was smoothly shifted to zero between 0.9 nm and 1.2 nm to reduce the cutoff noise. For electrostatic interactions, the coulombic potential, with a cutoff of 1.2 nm, was smoothly shifted to zero from 0 to 1.2 nm. Temperature was kept constant at 310 K using the Berendsen weak coupling algorithm with a time constant of 1 ps. The surface tension was varied between 0 and 40 mN m-1 via the Berendsen barostat to represent different respiration conditions.27 The compressibility was set to be 5 ×10−5 bar −1 in the lateral direction and 0 bar −1 in the normal direction. Periodic boundary conditions were implemented in all three directions. The time step of simulations was 20 fs, and the neighbor list for non-bonded interactions was updated every 10 steps. All simulations were performed using the open source code Gromacs 4.6.7.28 Snapshots were rendered by VMD.29

RESULTS AND DISCUSSION Bare Nanodroplets. We began CGMD simulations of the interaction between a bare nanodroplet of 14 nm and the PS layer under a fixed surface tension of 40 mN m-1, close to the alveolar surface tension during inspiration.30 As shown in Fig. 2a (the dynamic process can be found in Video S1), upon contacting the PS, the nanodroplet immediately adhered onto the layer. Unfavorable interactions between the hydrophilic nanoroplet and the hydrophobic PS tails were generated to disturb the ultrastructure of the PS layer. Shortly, a PS layer pore was opened, from which the nanodroplet was quickly transported (Fig. 2a), being reflected by a sudden decrease of the interaction energy between nanodroplet and water (Fig. 2b) and a simultaneous increase of the PS layer area (Fig. 2c). By contrast, no obvious change of the PS layer area was observed in the absence of nanodroplet (Fig. 2c). This is reminiscent of the instability of the natrual PS under highly humid conditions, as has been previously measured by Zuo et al. using the biophysical method.31 Recently, Hu et al. found that graphene oxide nanosheets have adverse effects on the ultrastructure of the PS by generating pores in the PS layer.32 Our simulations are consistent with those findings, suggesting that the unfavorable interactions between hydrophilic drugs or drug carriers and hydrophobic PS tails may induce the PS instability and the biophysical inhibition. ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

In our simulations, once a pore was generated, it rapidly expanded and developed into a more severe PS layer rupture (Fig. 2). Despite the fact that extra PS molecules exist in PS reserviors that can be transferred back to the interface to prevent the PS layer rupture,33 our results at least suggested that high dose of inhalation of nanodroplets may aggravate the PS perturbation.32 Even in the presence of PS reserviors, the opened PS layer pore may not close before completion of the nanodroplet translocation. To elaborate the PS perturbation by inhaled nanodroplets, we extracted the structural evolution of the PS layer before 30 ns and calculated the order parameters of all PS molecules (Fig. 3a, b). The order parameter has the form S =

1 (3cos 2 θ n − 1) , where θ n is the angle between bond connecting two adjacent tail 2

beads and the PS layer normal.32,

34

Before opening the PS layer pore, deposition of the

nanodroplet was found to disturb the local PS arrangement, as reflected in the transient order parameter diagram at 24 ns (Fig. 3b). The extent of the PS perturbation was quantified by counting the defected number of PS molecules with the order parameter smaller than 0.0. It appears that opening the PS layer pore was accompanied with a sudden increase of the defected number, while subsequent rupturing reduced the number by releasing surface energy of the PS layer (Fig. 3c). The larger defected number compared with that before opening the pore was due to the rearrangement of PS molecules around the nanodroplet, as revealed by radial distributions of both the density and the order parameter of PS molecules as a function of distance with center of the nanodroplet (Fig. 3d). Using a larger reference value of 0.25, middle between -0.5 and 1 from the formular, similar results were observed expect that a much larger number of PS molecules were identified to be defective even before rupturing of the PS layer (Fig. S1). Since the average order parameter of PS molecules under surface tension of 40 mN m-1 was 0.55, the reference value of 0.0 should be more appropriate to measure the PS perturbation by nanodroplets. During respiration, surface tension of the PS layer varies from tens mN m-1 to near zero.10, 30 Above simulations suggested that nanodroplets can be transported across the PS layer via opening destructive pores under deep inspiration (40 mN m-1). We next investigated the effect of surface tension on the nanodroplet interaction. As shown in Fig. 4a-c, the nanodroplet was found to adhere onto the PS layer under surface tensions of 30 mN m-1, 20 mN m-1 and 10 mN m-1. To decrease the interaction energy, the nanodroplet was partially wrapped by the PS layer, with the wrapping extent increasing with decrease of surface tension of the PS layer. A number of drug particles were found to diffuse across the PS layer to enter the water phase (Fig. 4a-c).


Page 6 of 25

Page 7 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

We calculated the number of drug particles diffusing across the PS layer, and found a roughly linear increase as a function of the simulation time (Fig. 4d). The transport rate being defined as r = ∂N / ∂t was found to increase with surface tension of the PS layer. In fact, the PS layer has been demonstrated to accelerate diffusion of oxygen through the interface, manifesting permeability of the PS layer for small molecules.35 At the molecular level, the

enhanced penetration of drug particles through the PS layer under higher surface tensions can be ascribed to the less compact stacking of PS molecules, similar with the enhanced permeability of the lipid bilayer membrane.36 Both experimental and simulation studies have demonstrated that nanoparticles can more easily penetrate through lipid bilayers under higher surface tensions.37-38 Once diffusing into the air or crossing the PS layer to enter the water phase, drug particles dispersed randomly and barely formed clusters. We thus used this property to count the number of drug particles in the nanodroplet. Specifically, one drug particle was regarded as part of the nanodroplet if the number of its neighbors was larger than 4. Two particles were defined as neighbors if their separation distance was smaller than 0.5 nm. As expected, we observed a corresponding decrease of the nanodroplet size (Fig. 4e), albeit with slight fluctuations due to the counting error. By summing the numbers of drug particles crossing the PS layer and staying in the nanodroplet, interestingly, they were found to keep nearly unchanged, regardless of different tensions (Fig. 4f). The fluctuations arose from both the diffusion of drug particles and the error of counting number of drug particles in the nanodroplet. These results suggested that inhaled nanodroplet deposit onto the PS layer to release drug particles that may directly diffuse across the PS layer. When the PS layer was compressed to a lower surface tension, the nanodroplet prompted collapse of the PS layer (Fig. 5 and Video S2), another indication of the PS perturbation. In the absence of nanodroplets, the PS layer was buckled under compression.39 Upon contacting the nanodroplet, the PS layer rapidly collapsed to wrap the nanodroplet to both release high surface energy of the PS layer under compression and decrease the interaction energy between nanodroplet and the PS (Fig. S2a, b). The wrapped nanodroplet was connected with the upper layer through a bilayer, which kept growing under constant compression of the layer to drive the nanodroplet continuously moving downward (Fig. S2c). Earlier simulations reported accumulation of water molecules in cave of the collapsed PS layer deposited with fullerenes,40 similar to the structure observed in our simulations. Note that rigid hydrophilic nanoparticles were found to directly penetrate through the PS layer regardless of different surface tensions.41 For hydrophilic nanoparticles entering the alveolar fluid, they were found to recruit PS molecules forming PS corona with a bilayer conformation, while a monolayer


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

conformation was formed on the hydrophobic NP surface.42 Those results suggest that nature of interactions between inhaled particles and the PS is determined by physicochemical properties of the particle, but also influenced by the local biological environment.

Lipid-coated Nanodroplets. To both improve efficiency of the pulmonary drug delivery and reduce the PS perturbation, we next designed a lipid-coated nanodroplet, with the inspiration coming from the promoted membrane transport of nanoparticles coated with polymer ligands.43-47 On the other hand, it has been evidenced that the PS layer is far from being a barrier against inhaled particles, but could be used as an efficient shuttle for pulmonary drug delivery.48 Therefore, rational design of surface coating on nanodroplets should better involve the PS layer facilitating the pulmonary drug delivery. Note that lipid nanoparticles have been synthesized by Liu et al. using the reverse micelle-double emulsion method,25 and showed performance for delivering insulin into the lungs.49 However, it still remains unclear about whether and how the lipid coating on nanodroplets promotes transport across the PS layer and reduce the PS perturbation. As shown in Fig. 6a-d, lipid-coated nanodroplets were partially or fully wrapped by the PS layer, depending on surface tension of the PS layer. Nine independent simulations were performed under different surface tensions ranging from 40 mN m-1 to 0 mN m-1. We found that the lipid-coated nanodroplet with a diameter of 16 nm and coated with 1,150 DPPC molecules can be fully wrapped by the PS layer if the surface tension was below 15 mN m-1 (Fig. 6e, f). Otherwise, only partial wrapping was achieved in the limited simulation time. Under a higher surface tension of 50 mN m-1, the PS layer immediately ruptured upon contacting the lipid-coated nanodroplet (Fig. S3a), owing to the PS depletion by interacting with the nanodroplet. By furthe increasing surface tension of the PS layer to 60 mN/m, it rapidly ruptured even in the absence of nanodroplets (Fig. S3b), being consistent with previous experimental and simulation results.50-51 We performed five independent simulations under a fixed surface tension of 40 mN m-1 to estimate properties of the interaction between lipid-coated nanodroplets and the PS layer under the inspiration condition.31 Repeated simulations showed that a number of drug particles slightly leaked from the nanodroplet during interaction with the PS layer (Fig. S4). Upon contacting the PS, the lipid-coated nanodroplet rapidly adhered onto and was partially wrapped by the PS layer (Fig. 7a). To further increase contact with the PS layer, it was deformed into an oblate shape, being reflected by an increase of the aspect ratio, defined as

AS = WXY / H Z , where WXY and H Z represent width and height of the nanodroplet, respectively (Fig. 7b). Notably, spreading deformation of the nanodroplet increased its surface ACS Paragon Plus Environment

Page 8 of 25

Page 9 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

tension,52 which induced a structural defect appearing at the wrapping front, as characterized by a sudden increase of the aspect ratio (Fig. 7b). We also observed a sudden striking decrease of drug particle number in the nanodroplet (Fig. 7b), suggesting that drug particles leaked from the nanodroplet and diffused in air or across the PS layer (Fig. 7c). Note that such behavior of slight drug leakage may not occur in real systems because extra PS molecules in reserviors can be transferred back to the layer to prevent drastic increase of surface tension of the PS layer. On the other hand, surface tension of the PS layer dynamically change during respiration. Thus the transient drug leakage may stop as further decrease of the surface tension. Nevertheless, considering that lipid-coated nanodroplets can be generated in experiments via nebulizing a drug solution admixed with lipids,25 our simulations suggest that a proper increase of the lipid concentration in drug solutions may increase density of lipids adsorbing onto nanodroplets, decreasing surface tension to prevent drug leakage during interactions with the PS layer. Similarly, it has been demonstrated that adding anionic surfactants into nebulisers decreases surface tension of the generated aerosols.53 Another feasible way to prevent the drug leakage is to decrease surface tension of the PS layer, because the lipid-coated nanodroplet can be wrapped by the PS layer more efficiently to reduce the nanodroplet deformation (Fig. 6). More importantly, when surface tension of the PS layer was decreased below 15 mN m-1, nanodroplets were fully wrapped by the PS layer (Fig. 8a, for the dynamic process, see Video S3), with nearly no nanodroplet deformation being accompanied (Fig. 8b). The more efficient wrapping compared with bare nanodroplets was attributed to the more favorable interactions between hydrophobic lipid tails coating on the nanodroplet and hydrophobic PS tails lining at the layer (Fig. S5). After completion of the wrapping, a bilayer vesicular structure was formed, in which the nanodroplet was encapsulated. The vesicle was stable and connected with the PS layer through a short neck, in which two SP-B proteins were anchored to help stablize the connection between the vesicle and the PS layer (Fig. 8c). We performed seven independent simulations, in five of which at least one SP-B protein was found to finally locate at the connection between the vesicle and the PS layer (Fig. S6). This observation suggests that SP-B proteins may preferentially locate at the connection between the PS layer and the formed vesicle, in agreement with one of the hypothesized functions of the SP-B protein in maintaining connection between the PS layer and the generated PS reserviors.54

Cell Membrane Interactions with Nanodroplets. Once transported across the PS layer, nanodroplets enter the alveolar fluid and get ready to interact with the epithelial cells or be cleared by the macrophages. Especially for the systemic therapy, they need to enter the


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

systemic circulation via cell uptake and translocation through epithelial cells or paracellularly via intercellular spaces particularly if tight junctions are disrupted.55 Both transport and diffusion of nanoparticles across the cell membrane can be modulated by changing the nanoparticle property, such as size, shape, elasticity and surface charge.56-58 Here, we consider the interaction between transported nanodroplets and the model epithelial cell membrane. Previous studies have revealed that vesicles can readily fuse with the plasma membrane of cells, requiring less time and energy than endocytosis,59-60 though the endocytic pathway can be preferred by increasing elasticity of the vesicles.61 Above simulations suggested that bare nanodroplets can be dispersively translocated across or wrapped by the PS layer (Fig. 4 and 5), whereas lipid-coated nanodroplets are readily wrapped by the PS layer to form bilayered vesicular structures (Fig. 6). The two wrapped structures were equilibrated in water for 10 ns before positioned above a preequilibrated lipid bilayer. It was found that the bare nanodroplet wrapped by the PS monolayer kept suspending above the lipid bilayer with no apparent interaction occuring in the limited simulation time (Fig. S7a). For the lipid-coated nanodroplet wrapped by the PS layer, it released a number of PS molecules into the upper leaflet of the membrane (Fig. S7b), indicating a vesicle fusion pathway. To facilitate fusion, an external force (1000 kJ/mol/nm) mimicking that generated by the trans-SNARE complex formation was exerted on centers of nanodroplets toward the membrane.62 Consequently, several PS molecules were transferred from the nanodroplet to the membrane (Fig. 9a). Such behavior increased surface tension of the nanodroplet to induce rupture occurring on top of the nanodroplet, releasing drug particles into the water phase (Fig. 9a). By contrast, the bilayered vesicular structure formed by the PS layer wrapping on the lipid-coated nanodroplet was more stable during interaction with the membrane. They are expected to fuse with the membrane to release drugs inside a cell.63-64 Otherwise, they might be internalized via endocytosis, penetrate through cell junctions to enter the systemic circulation,55, 65 or be cleared by the macrophages. Our result showed that only partial fusion was achieved in the limited simulation time (Fig. 9b),66 manifesting energy barriers exisiting in the fusion process.67 Note that the vesicle was formed by the PS layer wrapping on the lipid-coated nanodroplet. Such structure exhibits a lower surface tension to further increase the barrier for fusion with the plasma membrane.

CONCLUSIONS In summary, we have investigated transport of nanodroplets across the PS layer using CGMD


Page 10 of 25

Page 11 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

simulations. It was found that bare nanodroplets preferentially adhere onto the PS layer, releasing drugs to be directly translocated across the PS layer. The transport is quicker under higher surface tensions, but at the cost of opening pores that disrupt the ultrastructure of the PS layer. When the PS layer is compressed to lower surface tensions, the nanodroplet prompts collapse of the PS layer, another important indication of the PS perturbation. By coating the nanodroplet with lipids, it can be readily wrapped by the PS layer with slight PS perturbation. Moreover, the formed vesicular structures are expected to fuse with the cell membrane to release drugs into secondary organs. Properties of drug bioavailability, controlled drug release and enzymatic tolerance in real biological systems can be improved by the lipid coating strategy. Taken together, our simulations may assist in the design of more efficient and safe nanocarriers for pulmonary drug delivery.

ASSOCIATED CONTENT Supporting Information The Supporting information is available free of charge on the ACS Publication website at DOI: Comparison of the calcualted defected number using two different reference values; Detailed characterization of the PS layer collapse prompted by the bare nanodroplet; Interactions between the lipid-coated nanodroplet and the PS layer under higher surface tensions of 50 mN m-1 and 60 mN m-1; Repeated simulation results of the interaction between lipid-coated nanodroplets and the PS layer under a lower surface tension of 40 mN m-1; Comparison of the interaction energy between two different nanodroplets and the PS layer under different surface tensions; Repeated simulations results of preferential location of SP-B proteins in connection between the vesicle and the PS layer; Unbiased simulation results on interactions between nanodroplets and the plasma membrane (PDF). Rupturing of the PS layer induced by a bare nanodroplet (AVI) PS layer collapse prompted by a bare nanodroplet (AVI) PS layer wrapping on a lipid-coated nanodroplet (AVI)

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]

ORCID


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 25

Tongtao Yue: 0000-0002-8329-167X

Notes The authors declare no competing financial interest.

ACKNOWLEDGMENTS We are grateful for financial support from the National Natural Science foundation of China (no. 21303269), Science and Technology Major Project of Shandong Province (2016GSF117033),

the

Natural

Science

Foundation

of

Shandong

Province

(No.

ZR2018MC004), and Qingdao Science and Technology Project (no. 16-5-1-73-jch). This work was also partly supported by the Fundamental Research Funds for the Central Universities. MD simulations were performed at the National Super-computing Center in Shenzhen.

REFERENCES (1) Patton, J. S.; Byron, P. R. Inhaling medicines: delivering drugs to the body through the lungs. Nat. Rev. Drug Discov. 2007, 6 (1), 67-74. (2) Roy, I.; Vij, N. Nanodelivery in airway diseases: challenges and therapeutic applications. Nanomedicine 2010, 6 (2), 237-244. (3) Loira-Pastoriza, C.; Todoroff, J.; Vanbever, R. Delivery strategies for sustained drug release in the lungs. Adv. Drug Deliv. Rev. 2014, 75, 81-91. (4) Wearley, L. L. Recent progress in protein and peptide delivery by noninvasive routes. Crit. Rev. Ther. Drug Carrier Syst. 1991, 8 (4), 331-394. (5) Zhou, Q.; Tang, P.; Leung, S. S. Y.; Chan, J. G. Y.; Chan, H.-K. Emerging inhalation aerosol devices and strategies: Where are we headed? Adv. Drug Deliv. Rev. 2014, 75, 3-17. (6) Patton, J. S.; Bukar, J. G.; Eldon, M. A. Clinical Pharmacokinetics and Pharmacodynamics of Inhaled Insulin. Clin. Pharmacokinet. 2004, 43 (12), 781-801. (7) Dolovich, M. B.; Dhand, R. Aerosol drug delivery: developments in device design and clinical use. The Lancet 2011, 377 (9770), 1032-1045. (8) Tiwari, G.; Tiwari, R.; Sriwastawa, B.; Bhati, L.; Pandey, S.; Pandey, P.; Bannerjee, S. K. Drug delivery systems: An updated review. Int. J. Pharm. Investig. 2012, 2 (1), 2-11. (9) Bachofen, H.; Schürch, S. Alveolar surface forces and lung architecture. Comp. Biochem. Phys. A 2001, 129 (1), 183-193. (10) Zuo, Y. Y.; Veldhuizen, R. A. W.; Neumann, A. W.; Petersen, N. O.; Possmayer, F. Current perspectives in pulmonary surfactant — Inhibition, enhancement and evaluation. BBA - Biomembranes 2008, 1778 (10), 1947-1977. (11) Baoukina, S.; Tieleman, D. P. Computer simulations of lung surfactant. BBA Biomembranes 2016, 1858 (10), 2431-2440. (12) Wright, J. R. Pulmonary surfactant: a front line of lung host defense. J. Clin. Invest. 2003, 111 (10), 1453-1455. (13) Sano, H.; Kuroki, Y. The lung collectins, SP-A and SP-D, modulate pulmonary innate immunity. Mol. Immunol. 2005, 42 (3), 279-287. (14) Andersson, J. M.; Grey, C.; Larsson, M.; Ferreira, T. M.; Sparr, E. Effect of cholesterol on the molecular structure and transitions in a clinical-grade lung surfactant extract. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (18), E3592-E3601. ACS Paragon Plus Environment

Page 13 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(15) Fan, Q.; Wang, Y. E.; Zhao, X.; Loo, J. S. C.; Zuo, Y. Y. Adverse Biophysical Effects of Hydroxyapatite Nanoparticles on Natural Pulmonary Surfactant. ACS Nano 2011, 5 (8), 6410-6416. (16) Valle, R. P.; Huang, C. L.; Loo, J. S. C.; Zuo, Y. Y. Increasing Hydrophobicity of Nanoparticles Intensifies Lung Surfactant Film Inhibition and Particle Retention. ACS Sustain. Chem. Eng. 2014, 2 (7), 1574-1580. (17) Sachan, A. K.; Galla, H. J. Understanding the Mutual Impact of Interaction between Hydrophobic Nanoparticles and Pulmonary Surfactant Monolayer. Small 2014, 10 (6), 10691075. (18) Liu, P.; Izvekov, S.; Voth, G. A. Multiscale Coarse-Graining of Monosaccharides. J. Phys. Chem. B 2007, 111 (39), 11566-11575. (19) Marrink, S. J.; Risselada, H. J.; Yefimov, S.; Tieleman, D. P.; De Vries, A. H. The MARTINI force field: coarse grained model for biomolecular simulations. J. Phys. Chem. B 2007, 111 (27), 7812-7824. (20) Baoukina, S.; Monticelli, L.; Risselada, H. J.; Marrink, S. J.; Tieleman, D. P. The molecular mechanism of lipid monolayer collapse. Proc. Natl. Acad. Sci. U. S. A. 2008, 105 (31), 10803-10808. (21) Koldsø, H.; Sansom, M. S. Organization and dynamics of receptor proteins in a plasma membrane. J. Am. Chem. Soc. 2015, 137 (46), 14694-14704. (22) Uusitalo, J. J.; Ingólfsson, H. I.; Akhshi, P.; Tieleman, D. P.; Marrink, S. J. Martini coarse-grained force field: extension to DNA. J. Chem. Theory Comput. 2015, 11 (8), 39323945. (23) Wong-Ekkabut, J.; Baoukina, S.; Triampo, W.; Tang, I.-M.; Tieleman, D. P.; Monticelli, L. Computer simulation study of fullerene translocation through lipid membranes. Nat. Nanotechnol. 2008, 3 (6), 363-368. (24) Titov, A. V.; Král, P.; Pearson, R. Sandwiched graphene− membrane superstructures. ACS Nano 2009, 4 (1), 229-234. (25) Liu, J.; Gong, T.; Wang, C.; Zhong, Z.; Zhang, Z. Solid lipid nanoparticles loaded with insulin by sodium cholate-phosphatidylcholine-based mixed micelles: Preparation and characterization. Inter. J. Pharm. 2007, 340 (1), 153-162. (26) Li, Y.; Yuan, H.; von dem Bussche, A.; Creighton, M.; Hurt, R. H.; Kane, A. B.; Gao, H. Graphene microsheets enter cells through spontaneous membrane penetration at edge asperities and corner sites. Proc. Natl. Acad. Sci. U. S. A. 2013, 110 (30), 12295-12300. (27) Xu, Y.; Luo, Z.; Li, S.; Li, W.; Zhang, X.; Zuo, Y. Y.; Huang, F.; Yue, T. Perturbation of the pulmonary surfactant monolayer by single-walled carbon nanotubes: a molecular dynamics study. Nanoscale 2017, 9 (29), 10193-10204. (28) Hess, B.; Kutzner, C.; Van Der Spoel, D.; Lindahl, E. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 2008, 4 (3), 435-447. (29) Humphrey, W.; Dalke, A.; Schulten, K. VMD: visual molecular dynamics. J. Mol. Graph. 1996, 14 (1), 33-38. (30) Valle, R. P.; Wu, T.; Zuo, Y. Y. Biophysical influence of airborne carbon nanomaterials on natural pulmonary surfactant. ACS Nano 2015, 9 (5), 5413-5421. (31) Zuo, Y. Y.; Acosta, E.; Policova, Z.; Cox, P. N.; Hair, M. L.; Neumann, A. W. Effect of humidity on the stability of lung surfactant films adsorbed at air–water interfaces. BBA Biomembranes 2006, 1758 (10), 1609-1620. (32) Hu, Q.; Jiao, B.; Shi, X.; Valle, R. P.; Zuo, Y. Y.; Hu, G. Effects of graphene oxide nanosheets on the ultrastructure and biophysical properties of the pulmonary surfactant film. Nanoscale 2015, 7 (43), 18025-18029. (33) Schürch, S.; Qanbar, R.; Bachofen, H.; Possmayer, F. The surface-associated surfactant reservoir in the alveolar lining. Neonatology 1995, 67, 61-76. ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(34) Lin, X.; Bai, T.; Zuo, Y. Y.; Gu, N. Promote potential applications of nanoparticles as respiratory drug carrier: insights from molecular dynamics simulations. Nanoscale 2014, 6 (5), 2759-2767. (35) Olmeda, B.; Villén, L.; Cruz, A.; Orellana, G.; Perez-Gil, J. Pulmonary surfactant layers accelerate O2 diffusion through the air-water interface. BBA - Biomembranes 2010, 1798 (6), 1281-1284. (36) Olbrich, K.; Rawicz, W.; Needham, D.; Evans, E. Water permeability and mechanical strength of polyunsaturated lipid bilayers. Biophys. J. 2000, 79 (1), 321-327. (37) Hong, S.; Leroueil, P. R.; Janus, E. K.; Peters, J. L.; Kober, M.-M.; Islam, M. T.; Orr, B. G.; Baker Jr, J. R.; Banaszak Holl, M. M. Interaction of polycationic polymers with supported lipid bilayers and cells: nanoscale hole formation and enhanced membrane permeability. Bioconj. Chem. 2006, 17 (3), 728-734. (38) Yan, L.-T.; Yu, X. Enhanced permeability of charged dendrimers across tense lipid bilayer membranes. ACS Nano 2009, 3 (8), 2171-2176. (39) Baoukina, S.; Monticelli, L.; Risselada, H. J.; Marrink, S. J.; Tieleman, D. P. The molecular mechanism of lipid monolayer collapse. Proc. Natl. Acad. Sci. U. S. A. 2008, 105 (31), 10803-10808. (40) Barnoud, J.; Urbini, L.; Monticelli, L. C60 fullerene promotes lung monolayer collapse. J. R. Soc. Interface 2015, 12, 20140931 (41) Xu, Y.; Deng, L.; Ren, H.; Zhang, X.; Huang, F.; Yue, T. Transport of nanoparticles across pulmonary surfactant monolayer: a molecular dynamics study. Phys. Chem. Chem. Phys. 2017, 19 (27), 17568-17576. (42) Hu, Q.; Bai, X.; Hu, G.; Zuo, Y. Y. Unveiling the Molecular Structure of Pulmonary Surfactant Corona on Nanoparticles. ACS Nano 2017, 11 (7), 6832-6842. (43) Verma, A.; Uzun, O.; Hu, Y.; Hu, Y.; Han, H.-S.; Watson, N.; Chen, S.; Irvine, D. J.; Stellacci, F. Surface-structure-regulated cell-membrane penetration by monolayer-protected nanoparticles. Nat. Mater. 2008, 7 (7), 588-595. (44) Li, Y.; Li, X.; Li, Z.; Gao, H. Surface-structure-regulated penetration of nanoparticles across a cell membrane. Nanoscale 2012, 4 (12), 3768-3775. (45) Van Lehn, R. C.; Ricci, M.; Silva, P. H. J.; Andreozzi, P.; Reguera, J.; Voïtchovsky, K.; Stellacci, F.; Alexander-Katz, A. Lipid tail protrusions mediate the insertion of nanoparticles into model cell membranes. Nat. Commun. 2014, 5, 4482. (46) Li, Y.; Kröger, M.; Liu, W. K. Endocytosis of PEGylated nanoparticles accompanied by structural and free energy changes of the grafted polyethylene glycol. Biomaterials 2014, 35 (30), 8467-8478. (47) Xia, Q.-s.; Ding, H.-m.; Ma, Y.-q. Can dual-ligand targeting enhance cellular uptake of nanoparticles? Nanoscale 2017, 9 (26), 8982-8989. (48) Hidalgo, A.; Cruz, A.; Pérez-Gil, J. Barrier or carrier? Pulmonary surfactant and drug delivery. Eur. J. Pharm. and Biopharm. 2015, 95, 117-127. (49) Liu, J.; Gong, T.; Fu, H.; Wang, C.; Wang, X.; Chen, Q.; Zhang, Q.; He, Q.; Zhang, Z. Solid lipid nanoparticles for pulmonary delivery of insulin. Int. J. Pharm. 2008, 356, 333-344. (50) Zuo, Y. Y.; Keating, E.; Zhao, L.; Tadayyon, S. M.; Veldhuizen, R. A.; Petersen, N. O.; Possmayer, F. Atomic force microscopy studies of functional and dysfunctional pulmonary surfactant films. I. Micro-and nanostructures of functional pulmonary surfactant films and the effect of SP-A. Biophys. J. 2008, 94 (9), 3549-3564. (51) Baoukina, S.; Monticelli, L.; Marrink, S. J.; Tieleman, D. P. Pressure− area isotherm of a lipid monolayer from molecular dynamics simulations. Langmuir 2007, 23 (25), 1261712623. (52) Yue, T.; Zhang, X. Molecular modeling of the pathways of vesicle–membrane interaction. Soft Matter 2013, 9 (2), 559-569.


Page 14 of 25

Page 15 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(53) Mc Callion, O. N. M.; Taylor, K. M. G.; Thomas, M.; Taylor, A. J. The influence of surface tension on aerosols produced by medical nebulisers. Int. J. Pharm. 1996, 129 (1), 123136. (54) Baoukina, S.; Tieleman, D P. Lung Surfactant Protein SP-B Promotes Formation of Bilayer Reservoirs from Monolayer and Lipid Transfer between the Interface and Subphase. Biophys. J. 2011, 100 (7), 1678-1687. (55) Thorley, A. J.; Ruenraroengsak, P.; Potter, T. E.; Tetley, T. D. Critical Determinants of Uptake and Translocation of Nanoparticles by the Human Pulmonary Alveolar Epithelium. ACS Nano 2014, 8 (11), 11778-11789. (56) Mao, J.; Chen, P.; Liang, J.; Guo, R.; Yan, L.-T. Receptor-mediated endocytosis of two-dimensional nanomaterials undergoes flat vesiculation and occurs by revolution and selfrotation. ACS Nano 2016, 10 (1), 1493-1502. (57) Chen, P.; Huang, Z.; Liang, J.; Cui, T.; Zhang, X.; Miao, B.; Yan, L.-T. Diffusion and Directionality of Charged Nanoparticles on Lipid Bilayer Membrane. ACS Nano 2016, 10 (12), 11541-11547. (58) Liang, J.; Chen, P.; Dong, B.; Huang, Z.; Zhao, K.; Yan, L.-T. Ligand–Receptor Interaction-Mediated Transmembrane Transport of Dendrimer-like Soft Nanoparticles: Mechanisms and Complicated Diffusive Dynamics. Biomacromolecules 2016, 17 (5), 18341844. (59) Düzgüneş, N.; Nir, S. Mechanisms and kinetics of liposome–cell interactions. Adv. Drug Deliv. Rev. 1999, 40 (1-2), 3-18. (60) François-Martin, C.; Rothman, J. E.; Pincet, F. Low energy cost for optimal speed and control of membrane fusion. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (6), 1238-1241. (61) Guo, P.; Liu, D.; Subramanyam, K.; Wang, B.; Yang, J.; Huang, J.; Auguste, D. T.; Moses, M. A. Nanoparticle elasticity directs tumor uptake. Nat. Commun. 2018, 9 (1), 130. (62) Oelkers, M.; Witt, H.; Halder, P.; Jahn, R.; Janshoff, A. SNARE-mediated membrane fusion trajectories derived from force-clamp experiments. Proc. Natl. Acad. Sci. U. S. A. 2016, 113 (46), 13051-13056. (63) Shillcock, J. C.; Lipowsky, R. Tension-induced fusion of bilayer membranes and vesicles. Nat. Mater. 2005, 4 (3), 225-228. (64) Baoukina, S.; Tieleman, D. P. Direct Simulation of Protein-Mediated Vesicle Fusion: Lung Surfactant Protein B. Biophys. J. 2010, 99 (7), 2134-2142. (65) Takechi-Haraya, Y.; Goda, Y.; Sakai-Kato, K. Control of Liposomal Penetration into Three-Dimensional Multicellular Tumor Spheroids by Modulating Liposomal Membrane Rigidity. Mol. Pharm. 2017, 14, 2158-2165. (66) Kawamoto, S.; Klein, M. L.; Shinoda, W. Coarse-grained molecular dynamics study of membrane fusion: Curvature effects on free energy barriers along the stalk mechanism. J. Chem. Phys. 2015, 143 (24), 243112. (67) Grafmüller, A.; Shillcock, J.; Lipowsky, R. Pathway of membrane fusion with two tension-dependent energy barriers. Phys. Rev. Lett. 2007, 98 (21), 218101.


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. Coarse-grained models and simulation system setup. (a) Schematic of the PS layer system. (b) Coarse-grained lipid and protein models used in our simulations. (c) The bare nanodroplet model. (d) The lipid-coated nanodroplet model. (e) The lipid bilayer model.


Page 16 of 25

Page 17 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 2. The PS layer rupture and dispersive translocation of nanodroplet across the layer under tension of 40 mN m-1. (a) Time sequence of typical snapshots showing the process of PS layer rupture and nanodroplet translocation. (b) Time evolutions of interaction energies of nanodroplet with both water (black) and PS molecules (red). (c) Time evolutions of the monolayer area with (red) and without (black) deposition of nanodroplets.


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Characterization of the PS perturbation by deposited nanodroplets. (a) Time sequence of typical snapshots from the top view. The nanodroplet was not shown for clarity. (b) Time sequence of the PS order parameter diagram. (c) Time evolution of the defected number. (d) Radial distributions of both the PS density and PS order parameter as a function of distance from the nanodroplet at t = 30 ns. Surface tension of the PS layer is fixed at 40 mN m-1.


Page 18 of 25

Page 19 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 4. Effect of surface tension of the PS layer on the interaction with bare nanodroplets. (a-c) Final simulated snapshots from the side view under surface tensions of 30 mN m-1, 20 mN m-1 and 10 mN m-1, respectively. (d) Time evolutions of the number of drug particles diffusing across the PS layer under different surface tensions. (e) Time evolutions of the number of drug particles in the nanodroplets. (f) Time evolutions of the total number of drug particles both diffusing across the PS layer and staying in the nanodroplets.


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. Bare nanodroplet prompts collapse of the PS layer under zero surface tension.


Page 20 of 25

Page 21 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 6. Effect of surface tension of the PS layer on interactions with lipid-coated nanodroplets. (a-d) Final simulated snapshots under surface tensions of 40 mN m-1, 30 mN m1

, 20 mN m-1 and 10 mN m-1, respectively. (e) Time evolutions of the wrapping percentage of

lipid-coated nanodroplets by the PS layer under tensions ranging from 40 mN m-1 to 0 mN m-1 with an interval of 5 mN m-1. (f) Time evolutions of the distance between centers of nanodroplets and the PS layer.


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7. Partial PS layer wrapping on a lipid-coated nanodroplet under a higher surface tension of 40 mN m-1. (a) Time sequence of typical snapshots showing the partial wrapping, with the rupturing location marked with blue arrows. (b) Time evolutions of the nanodroplet drug number and the aspect ratio. The time of drug leakage reflected by a sudden decrease of the nanodroplet drug number and a striking increase of the aspect ratio is marked by a dashed blue line. (c) Enlarged local structure with the drug leakage labeled by a white dashed circle.


Page 22 of 25

Page 23 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Figure 8. Full PS layer wrapping on a lipid-coated nanodroplet under a lower surface tension of 10 mN m-1. (a) Time sequence of typical snapshots showing the wrapping process. (b) Time evolutions of the nanodroplet drug number and the aspect ratio. (c) Final structure showing locations of both SP-B (gray) and SP-C (orange) proteins in the formed vesicular structure.


Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 9. Cell membrane interactions with two different nanodroplets wrapped by the PS layer. (a) The bare nanodroplet undergoes rupture by interacting with the cell membrane. (b) The lipid-coated nanodroplet partially fuses with the cell membrane in the limited simulation time.


Page 24 of 25

Page 25 of 25 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

For Table of Contents Use Only.


The Role of Lipid Coating in Transport of ... - ACS Publications

Recommend Documents