Membrane Wrapping Pathway of Injectable Hydrogels: From Vertical

Jul 11, 2019 - This website uses cookies to improve your user experience. By continuing to use the site, you are accepting our use of cookies. Read th...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF SOUTHERN INDIANA

Biological and Environmental Phenomena at the Interface

Membrane Endocytosis Pathway of Injectable Hydrogels: from Vertically Capillary Adhesion to Laterally Compressed Wrapping Xianyu Song, Chongzhi QIAO, Teng Zhao, Bo Bao, Shuangliang Zhao, Jing Xu, and Honglai Liu Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.9b01395 • Publication Date (Web): 11 Jul 2019 Downloaded from pubs.acs.org on July 17, 2019

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 32 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

Graphic abstract

By means of large-scale computational simulations, the membrane internalization of injectable hydrogels was studied, and this complicated process transforms from vertically capillary adhesion to laterally compressed wrapping.

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

1

Membrane Wrapping Pathway of Injectable Hydrogels: from Vertical Capillary

2

Adhesion to Lateral Compressed Wrapping

3

Xianyu Song1, Chongzhi Qiao1, Teng Zhao1, Bo Bao1, Shuangliang Zhao1,*, Jing Xu1, and

4

Honglai Liu2,*

5

1

State Key Laboratory of Chemical Engineering and School of Chemical Engineering, East

6

China University of Science and Technology, Shanghai, 200237, China

7

2

8

Engineering, East China University of Science and Technology, Shanghai, 200237, China

State Key Laboratory of Chemical Engineering and School of Chemistry and Molecular

9 10 To whom correspondence should be addressed. Email: [email protected] (SZ) or

11

*

12

[email protected] (HL)

13 14 15 16 17 18 19 20 21

1 / 31

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32 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

22

Abstract

23

Membrane wrapping pathway of injectable hydrogels (IHs) plays a vital role in the

24

nanocarrier effectiveness and biomedical safety. Whereas considerable progress in

25

understanding this complicated process has been made, the mechanism behind this process

26

has remained elusive. Herein, with the help of large-scale dissipative particle dynamics

27

simulations, we explore the molecular mechanism of membrane wrapping by systematically

28

examining the IH architectures and hydrogel-lipid binding strengths. To the best of our

29

knowledge, this is the first report on the membrane wrapping pathway on which IHs

30

transform from vertical capillary adhesion to lateral compressed wrapping. This

31

transformation results from the elastocapillary deformation of networked gels and nanoscale

32

confinement of bilayer membrane, and it takes long time for the IHs to be fully wrapped

33

owing to the high energy barriers and wrapping-induced shape deformation. A collapsed

34

morphologies and small compressed angles are identified in the IHs with thick shell or strong

35

binding strength to lipids. In addition, the IHs binding intensively to membrane exhibit

36

special nanoscale mixing and favorable deformability during the wrapping process. Our study

37

provides a detailed mechanistic understanding of the influence of architecture and binding

38

strength on the IHs’ membrane wrapping efficiency. This work may serve as rational

39

guidance for the design and fabrication of IH-based drug carriers and tissue engineering.

40

Keywords: membrane wrapping; injectable hydrogels; lipid membrane; computational

41

simulations; elastocapillary deformation

42 2 / 31

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

43

1. Introduction

44

Cell-based therapy is considered to be one of the most booming strategies for treating

45

various diseases.1-2 The major challenge in cell-based therapy lies in designing ideal

46

nanocarriers to minimize the cytotoxic side effects and deliver the therapeutic cells to target

47

sites efficiently and precisely.3 The diverse chemical compositions and elastic network

48

structure endow injectable hydrogels (IHs) with unique properties, such as deformability,

49

reversibility, self-healing, and stimuli responsiveness, when compared to ordinary rigid

50

nanomaterials.

51

capable of gelling in situ, and thus they readily integrate with the cell membrane and provide

52

a local therapeutic opportunity at the targeted sites. Due to their great value and huge

53

potential in application, such as biocompatibility, ease of systemic administration, and

54

minimal invasion,6-7 previous studies have demonstrated that IHs could provide a versatile

55

platform for biomedical applications, especially for the cell-based therapy.8

4-5

Once penetrating successfully into the cell membrane, the hydrogels are

56

By introducing the chemophotothermal agents9-10 and advanced nanomaterials (i.e.,

57

nanoparticles,11-12 carbon nanotubes,13 graphene oxide nanosheets14), various synergistic

58

therapy strategies were proposed. Notably, there has been inspiring achievements in

59

developing IHs for facile cellular encapsulation and controlled delivery toward smart control

60

and scalable engineering.5, 15-16 Moreover, these cell-based therapy systems exhibit excellent

61

multi-function and high drug loading as well as prolonged drug retention. Nevertheless, there

62

is still a lack of molecular understanding of the elastocapillary deformation of the IH systems.

63

It has been reported that polymer capsules with weak stiffness were efficiently wrapped by 3 / 31

ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32 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

64

the cell membrane,17 and during this process, the membrane perturbation or fluctuation may

65

be significantly affected.18 This special membrane wrapping pathway impedes possible

66

transmembrane translocation because the energy barriers to be overcome during the

67

membrane wrappings is smaller.19

68

Besides, the adhesive stability and structure-property relationship of the IHs-membrane

69

remain unclear in their apllications. For example, the inability to retain the structural integrity

70

after injection limits their applications in the cell-based therapy.20 Recently, Yang et al.4

71

found that the morphologies of IHs play a significant role in the invasive process.

72

Specifically, the thinning IHs with two-dimensional geometry facilitates its encapsulation by

73

a membrane, thereby greatly improving the efficacy of cell-based therapy.21-23 In addition, it

74

is helpful to monitor the noninvasive pathway of IHs over time for optimizing the therapeutic

75

efficacy and safety,21 and it’s of particular interest to understand the interaction between the

76

IHs and membranes, e.g., how the IHs adapt to the surrounding heterogeneous matrix of a

77

bilayer membrane, and whether this environment affects the IHs’ targeting and function.

78

Though the incorporation of fluorescence dyes into IHs21-22 and in-situ imaging

79

technologies23 have been developed to monitor the in-vivo membrane wrapping of IHs, the

80

underlying mechanism still remains elusive, which is of great importance towards the

81

highly-efficient design of IHs-based drug carriers and tissue engineering.24

82

To tackle these problems, large-scale dissipative particle dynamics (DPD) simulations

83

are performed to investigate the membrane wrapping pathway of the IHs. With the aid of

84

coarse-graining, the specific information regarding the membrane wrapping dynamics of the 4 / 31

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

85

IHs can be obtained over large time scale.25-26 Furthermore, by varying the morphological

86

characteristics and hydrogel-lipid binding strength, the membrane wrapping pathway and

87

corresponding membrane perturbation are examined.

88

2. Simulation method and computational details

89

2.1 Coarse-grained models

90

DPD method represents a classical coarse-grained simulation, in which a cluster of

91

atoms is grouped as one single DPD bead. Specifically, the lipid molecule is coarse-grained

92

by simplifying the lipid heads (H bead) and tails (T bead) into hydrophilic and hydrophobic

93

beads, respectively. There are different ways of coarse-graining in the existing literature to

94

construct lipid membranes including the double-tail model and linear chain model.27-30 Each

95

model has its own merits. For example, the linear lipid model was employed to study the

96

bilayer properties, vesicle structure, and the drug delivery processes, which showed good

97

agreements with experimental observations.31-33 However, Shillcock et al.31 argued that the

98

lamellar phase composed of linear chain lipids is highly dependent on the tail stiffness of

99

lipid, which may affect the equilibrium structure. Nevertheless, here we adopt the linear chain

100

model HT3 in our simulations for reducing the computational cost. Furthermore, the

101

employment of this model allows us to describe the experimental process,34 in which the size

102

of a single injectable hydrogel ranges typically between 50 nm and 5 m.35 The

103

coarse-grained model of lipid and the construction of bilayer membrane are presented in

104

Figure 1A. One of the major advantages of the present simulations lies in that a large-scale

105

membrane with 114.0×114.0 nm2 should be constructed with sufficient lipid molecules, 5 / 31

ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32 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

106

Langmuir

which allows for mimicking the dynamical process of membrane wrapping.

107 108

Figure 1. Schematic representation of computational models for bilayer membrane and IH:

109

(A) coarse-grained models of bilayer membrane and single lipid molecule, (B) lateral view

110

and (C) half-sectional view of the coarse-grained models of spherical hydrogel capsule. The

111

bilayer membrane consists of 28800 lipids. Red beads signify the lipid heads, yellow beads

112

represent the lipid tails, and the cyan beads indicate the hydrogel components.

113

The hydrogel is coarse-grained by simplifying each polymer monomer into one single

114

DPD bead (HY bead) with uniform size  , and these beads are internally connected by

115

tetra-functional cross-linkers.36-37 Herein, the bead size  equals to the characteristic length

116

scale rc as discussed below. In order to examine the influence of hydrogel architecture on

117

the membrane wrapping, two types of hydrogels are designed, namely, the spherical stuffed

118

hydrogel and spherical hydrogel capsules with different shell thicknesses.

119 120

By tailoring an ideal diamond cubic crystal structure, a spherical networked hydrogel with the radius Rgel  20 is firstly obtained. Then, the inner core with radius Rs  5 is 6 / 31

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

Page 8 of 32

121

eliminated. The obtained hydrogel capsule is shown in Figure 1B and 1C. The detailed

122

fabrication method of hydrogel models can be found in previous studies.36-37 We note that

123

this model presents a comparable size with practical hydrogel as stated below. In our

124

simulations, the moderate crosslinked hydrogel with the crosslink density  =0.0217 is used,

125

where 

126

polymer particle and crosslinker particle, respectively.

127

2.2 Simulation methodology

is defined as  = N p (N p +N c ) 34 with N p and N c being the numbers of

128

Based on the dynamics of soft particles interacting with each other through conservative

129

force FijC , dissipative force FijD , and random forces FijR , the time evolution of each DPD

130

bead i subject to the total forces Fi follows the standard Newtonian equation of motion:38

Fi  mi

131

dvi dr   FijC  FijD FijR , vi  i , dt i  j dt

(1)

132

where vi , ri , and mi are the velocity, position, and mass of the 𝑖-th bead, respectively. The

133

conservative force is soft-repulsive, and it reads:

FijC  aij (1  rij )rˆij ,

134

(2)

135

where aij is the maximum repulsion force between the bead i and j, rij  ( ri  rj ) rc ,

136

rij  rij , and rˆij  rij rij . The dissipative force FijD is given as:

FijD  D ( rij )( rˆij  vij )rˆij ,

137

(3)

138

where  corresponds to the strength of friction, D ( rij )  (1  rij )2 , and vij = vi - v j . The

139

random force is taken as:

FijR  R ( rij )ij rˆij ,

140 141

where ij

(4)

represents a zero-mean Gaussian random variable of unit variance, 7 / 31

ACS Paragon Plus Environment

Page 9 of 32 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

2k BT  , where  is the friction

142

R ( rij )  1  rij , and  is a noise amplitude and equals to

143

factor.39 In addition, the corresponding harmonic force U ijS on the bead pair i and j is given

144

as:

U ijS  Cb ( r0  rij )rij ,

145

(5)

146

where Cb and r0 are the force constant and equilibrium distance, respectively. Herein, Cb =

147

50 k BT rc 2 and r0 = 0.7 rc .27

148

Meanwhile, the equation of the bond angle potential is given by

149

U ijA =Ca (ij  0 )2

150

where Ca and 0 are the bending stiffness and equilibrium angle, respectively. Here Ca =

151

3.0 k BT and 0 = 180°.27

152

2.3 Simulation details

(6)

153

Thermoresponsive microgels are among the most investigated stimuli-responsive

154

systems in drug delivery applications.40 Poly(N-isopropylacrylamide) (PNIPAM) represents

155

one of the thermally sensitive polymers and displays a significant volumetric change at an

156

LCST (about 32 ℃ ). With the aids of thermoresponsive hydrogels, drug agents can be

157

effectively shielded and protected until reaching a treatment site due to the microgels’

158

swelling properties below the LCST, and then release a predefined dosage of the drug within

159

a locally heated tumour (about 40~42 ℃ ) because of the collapsed behavior of PNIPAM

160

hydrogels above the LCST.41 We have investigated the swelling/deswelling behavior of

161

PNIPAM microgels.36-37 Previous studies have demonstrated that PNIPAM chains become

162

insoluble in aqueous solutions and tend to aggregate above LCST.42 Thus, the 8 / 31

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

163

PNIPAM-based hydrogels are the hydrophobic state when they locate at the heated tumor and

164

penetrate into the membrane. For any two beads of the same type, we adopt the conservative

165

parameter aii =25 k BT rc .43 For any two beads of different types in lipid models, the

166

conservative parameter aW -H is 25 k BT rc , and aW -T  aT -H is 100 k BT rc ,44-45 where W, H,

167

T signify the water bead, the lipid head bead, and lipid tail bead, respectively. The DPD

168

interaction parameter aW -H = 25 k BT rc indicates that water bead and lipid head bead are

169

miscible and belongs to the same group of DPD component. To represent the hydrophobic

170

state of PNIPAM-based hydrogels, we take a conservative parameter aT -HY =25 k BT rc .

171

Meanwhile, binding = aH  HY  aW  HY vary from 30 to 85 k BT rc to describe the different

172

hydrogel-lipid binding strengths, which is aroused from the certain surface or composition

173

modifications of PNIPAM-based hydrogels, such as block polymer grafting with different

174

types of polymer monomers or grafting ratios.46-48 The conservative parameters are listed in

175

Table 1.

176

In the DPD simulations, we take rc as the characteristic length scale and k BT as the

177

characteristic energy scale, and the characteristic time scale is defined as  = mrc 2 k BT .49

178

The velocity-Verlet integration algorithm is adopted with t  0.04 ,50 and a total of 6 × 105

179

DPD simulation steps are performed to obtain the equilibrium state. On comparing a

180

reasonable estimation for lipid bilayer thickness of 4 nm in experimental results, the thickness of

181

the linear lipid molecule is around 3.5 rc .27 Thereby, the basic unit length in DPD simulations is

182

about rc =1.14 nm. Mapping the diffusion coefficient of around 5 μm2s−1 of the POPC

183

bilayer, the simulation time-scale is about 0.032 ns/step.51 Therefore the IH with a radius of 9 / 31

ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32 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

184

20  corresponds a real size of 45.6 nm, comparable with practical hydrogels. For instance,

185

the PNIPAM-based hydrogels, displaying a significant volumetric change at a lower critical

186

solution temperature (LCST, 32℃), is of about 55nm above LCST.52-54 In our study, all

187

simulations are accomplished in a cubic box with a size of 100 × 100 × 80 rc with periodic

188

boundary conditions applied in all three directions, where the reduced particle density of

189

simulation system is 3.0. All simulations are conducted in the NVT ensembles using

190

LAMMPS package.55

3

191

Table 1. Conservative parameters expressed in k BT rc units W (water) H (lipid head) T (lipid tail) HY (hydrogels) W

25

H

25

25

T

100

100

25

HY

30~85

30~85

25

192

3. Results and discussion

193

3.1 Membrane wrapping mechanism of IHs

25

194

Here, we first investigate the membrane wrapping behavior of spherical stuffed IH. The

195

resultant morphological evolution of membrane wrapping is shown in Figure 2A. Figure

196

2A1-A3 show that the spherical IH is partially engulfed by the membrane and forms a

197

gyro-shaped morphology. Specifically, we observe a sharpened adhesion area at the wrapping

198

position due to the vertical structural elongation of the polymeric matrix. In this case, the

199

generated geometry allows a comprehensive wrapping with a lower energy barrier. The

200

subtle observation in Figure 2A1-A3 demonstrates that the membrane forms small

201

undulations and then buckles. This result manifests a typical capillary phenomenon. Herein, 10 / 31

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

202

we define the initial wrapping process (t≤0.32µs) as vertical capillary adhesion stage.

203

Likewise, experimental and theoretical observations have reported the capillary adhesion

204

phenomenon in the rigid nanoparticle-membrane systems.56-57 Originating from the van der

205

Waals attractions between the IH and membrane, the formation of vertical capillary adhesion

206

provides a driving force for the membrane to wrap the corner of IHs. As increasing the

207

simulation time (t>0.32µs), the membrane progressively engulfs the IH, and a rugby-shaped

208

morphology is observed, as shown in Figure 2A4-A6. In this stage, the sharpened adhesion

209

area disappears, and the lateral compressed state is generated in the middle of the bilayer

210

membrane. The comprehensive balance between membrane tension and hydrogel

211

deformability plays a significant role in this stage. We identify the following process as

212

lateral compressed wrapping stage. In generally, all these stages are controlled by the

213

elastocapillary deformation of networked gels and nanoscale confinement of bilayer

214

membrane.

215

The density profiles over simulation time are computed. Figure 2B1-B6 show that the

216

density distribution of the IH at the Z direction becomes more concentrated as increasing

217

the simulation time. Ultimately, the density distributions of the IH and membrane coincide at

218

the latter simulation time. In summary, all these results endorse a unique membrane wrapping

219

pathway that the IHs transform from vertical capillary adhesion to lateral compressed

220

wrapping during the wrapping process. Similar observations can be found in the IHs with

221

different architecture characteristics, as included in the Supplementary Information (SI).

222

Importantly, this exclusive membrane wrapping pathway is just identified in the soft IHs 11 / 31

ACS Paragon Plus Environment

Page 12 of 32

Page 13 of 32 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

223

rather than in ordinary rigid nanomaterials.11-12 13 14 Besides, the designed IHs finally embed

224

into the lipid membrane, rather than penetrating through the membrane, which is also

225

different from the membrane endocytosis of rigid nanoparticles. Though affected by the

226

bending modulus of the membrane,58 the membrane translocation induced by the IHs is

227

majorly determined by the shape deformability of the IHs, which is not favorable for

228

membrane rupture. Conversely, the intrinsic rotation of rigid nanoparticles offers a strong

229

driving force for membrane rupture.59-60 The dynamical properties of this structural transition

230

would be further discussed below.

12 / 31

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

231 232

Figure 2. Membrane wrapping pathway of IHs: (A1-A6) representative snapshots of

233

membrane wrapping of spherical hydrogel with a radius Rgel  20 at different simulation

234

time, (B1-B6) density profiles of the spherical hydrogel at the direction of simulation box-Z

235

and (C1-C6) 2D density mappings of IHs corresponding to snapshot figures in (A1-A6). The

236

color bar represents hydrogel density in the units of  2 . The solvent beads are not shown

237

for clarity and the same as below. 13 / 31

ACS Paragon Plus Environment

Page 14 of 32

Page 15 of 32 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

238

The two-dimensional (2D) density maps of the IHs are given in Figure 2C1-C6. The

239

ring-like density distributions are obtained in the spherical stuffed IHs. This demonstrates a

240

hollow structure is formed in the IH matrixes, which may result from strong hydrogel-lipid

241

binding interactions and good deformability of the IHs. The chain configuration can be

242

characterized by the radius-of-gyration paralleled to individual axes, which is defined as 61

243

Rg 

1 n  ( Ri  Rc )2 n i 1

(7)

244

where  corresponds to each component (  being x, y, z), Ri is the position of the i-th

245

polymer bead, and Rc is the position of the mass center of the polymer at the  direction,

246

and

247

radius-of-gyration for the IHs is observed at the directions of axis-X and axis-Y, while a

248

decreased radius-of-gyration is obtained at the direction of axis-Z. This significant hydrogel

249

deformation originates from the strong van der Waals attraction between the IH and

250

membrane along the Z direction. Namely, once the hydrogels contact with the membrane,

251

the vertical elongation of the IH is quickly identified. Consequently, at the initial wrapping

252

process (t≤0.32µs), the radius-of-gyration for the IHs at the direction of axis-Z is larger than

253

those at the direction of axis-X and axis-Y. The inverse variation trend is observed at the

254

following wrapping process (t>0.32µs). Moreover, the total radius-of-gyration for the IHs

255

also indicates a structural transition at 0.32µs, as shown in Figure 3B. The evolution of

256

average height for the IHs slowly attenuates when t≤0.32µs, then rapidly decreases, and

257

eventually reaches a platform. All these results confirm a structural transition in the wrapping

258

dynamics as mentioned above.

denotes the ensemble average. As shown in Figure 3A, an increased

14 / 31

ACS Paragon Plus Environment

Langmuir

(A)

R g /

2.8

2.6

2.4

Axis-X Axis-Y Axis-Z

0.1 259

1 Time/s

10

70

(B) 4.68

65 60

Total Rg/

4.64

55

4.6

50 45

4.56

40

Height/

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 16 of 32

35

4.52 0.1 260

1 Time/s

10

30

261

Figure 3. Membrane wrapping dynamics of IHs: (A) radius-of-gyration of the spherical

262

hydrogel at different directions of simulation box and (B) average height and total

263

radius-of-gyration of the spherical hydrogel.

264

3.2 Effects of hydrogel shell thickness on wrappings

265

Microgel capsules have been utilized in various encapsulation and controlled release

266

applications.36, 62-63 By varying the hydrogel shell thickness ranging from Rs =8.5 to 20  ,

267

the effect of hydrogel architectures on the membrane wrapping is studied. The membrane

268

wrapping of the hydrogel capsule with different thickness is presented in Figure 4A1-A4. 15 / 31

ACS Paragon Plus Environment

Page 17 of 32 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

269

Similarly, the equilibrated morphologies of the IHs also demonstrate the lateral compression

270

with different degrees. Compared with the spherical stuffed IHs, the hollow structure isn’t

271

observed in the hydrogel capsules, as confirmed in Figure 4B1-B4. The high affinity

272

between the IH and membrane induces the formation of hollow structure in the stuffed IHs,

273

while the favorable lateral deformation makes the hollow IH tend to be collapsed state.

274

Moreover, it seems that the hydrogel nanosheet is generated at the shell thickness Ds 

275

8.5  , as illustrated in Figure 4A4 and 4B4. This is because the hydrogel capsules are

276

favorable to be laterally compressed, as confirmed by the total radius-of-gyration for the IHs.

277 278

Figure 4. Membrane wrapping pathway of the IHs with different shell thicknesses. (A1-A4)

279

Representative snapshots of spherical hydrogel with different shell thicknesses Ds  20.0,

280

15.0, 11.25, 8.5  , respectively. The hydrogel with a shell thickness of

281

corresponds to the spherical stuffed hydrogel. (B1-B4) 2D density maps of the IHs

282

correspond to the snapshot Figures A1-A4. The color bar indicates the hydrogel density in the 16 / 31

ACS Paragon Plus Environment

Ds  20.0 

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 18 of 32

283

units of  2 . (C) The evolution of total radius-of-gyration for the IHs. (D) evolution of

284

contact angle for the IHs embedded in the bilayer membrane. The equilibrium simulation

285

time is 19.2 s.

286 287

Figure 4C displays that the thicker capsule generates larger total radius-of-gyration for

288

the IHs. Meanwhile, the structural transition can be further confirmed by the evolution of

289

total radius-of-gyration. Furthermore, we introduce the compressed angle conception to

290

quantitatively interpret the observed relation between structural compression and IHs

291

architecture. The definition of compressed angle is illustrated in Figure 4D. The decrease of

292

the shell thickness results in a decreased compression angle. It is because it’s easier for a

293

thinner IH to generate shape deformation. The compressed angle of the spherical stuffed IHs

294

is smaller than that of the hydrogel capsule with a shell thickness Rs =15  . This result can

295

be interpreted by the fact that the spherical stuffed IHs form a hollow structure as confined at

296

the membrane (see Figure 4B1). In order to quantify the lateral compression strength, the

297

geometrical shape parameter  is introduced, which is defined as 64

298 299

1 2

 = g( Rgx2  Rgy2 ) Rgz2 .

(8)

Note that   1 indicates a vertical elongation at the axis-Z direction and   1

300

represents a spherical shape or isotropic object, while   1 suggests a lateral compression.

301

Figure 5A shows the dynamical shape variation of the IHs, which indicates the IHs undergo

302

structure transition over simulation time. Especially, the structure transition time depends on

303

the hydrogel shell thickness: ts  20.0 < ts 15.0 < ts 11.25 < ts 8.5 . The evolution of the 17 / 31

ACS Paragon Plus Environment

Page 19 of 32 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

304

wrapping ratio for each IH is recorded in Figure 5B. Here the wrapping ratio  H is defined

305

as the ratio of the instantaneous height variation of IHs to the maximal height variation of

306

IHs. It should be noted that the wrapping ratio defined here is different from the conventional

307

one,  S , defined by contacting surface conception.65-66 The relation and difference between

308

them are  S =f S ( H ,  )g H , where H is height variation of IHs, and  represents the

309

shape deformation parameter of IHs.

310

coefficient. We take f S ( H ,  ) as a constant for simplicity. Nevertheless, the wrapping ratio

311

introduced in this work can quantify the wrapping process, as confirmed by the scaling

312

factors. (see Figures 5B and 8B) As increasing the simulation time, the wrapping ratio for

313

each IH gradually increases, and then reaches a platform. At the monotonous increasing

314

region, the wrapping ratio shows approximate scaling from  ~ t 0.928 to  ~ t1.908 . It seems

315

that the hydrogel shell thickness has a weak influence on the wrapping time tw (defined as

316

the time required for a successful IH internalization). The lateral compression inhibits or

317

frustrates the wrapping process due to the wrapping-induced shape deformation, thus, the

318

wrapping time tw is dominated by the lateral compressed wrapping stage. The theoretical

319

computation of wrapping time is included in the SI.

f S ( H ,  ) is the H-dependent and  -dependent

18 / 31

ACS Paragon Plus Environment

Langmuir

Structural compression

2 (A)1.8 1.6

Rs=20.0 Rs=15.0 Rs=11.25

1.4

Rs=8.5 laterally compressed wrapping

1.2 1

vertically capillary adhesion

0.8 0.1

320

1 Time/s

10 Rs=20.0

(B)

Rs=15.0 Rs=11.25

1

Rs=8.5

tio Wrapping ra

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 20 of 32

0.1 Kslope2=0.928

0.01

321

Rs=8.5 Rs=11.25 Rs=15.0 Rs=20.0

Kslope1=1.908

10

T

1 /s ime

0.1

322

Figure 5. Membrane wrapping dynamics of different hydrogel shell thickness: (A) structural

323

compression of IHs and (B) wrapping ratio of IHs.

324

3.3 Effects of hydrogel-lipid binding strength on wrappings

325

The conservative parameter aij between both beads of different types corresponds to

326

the mutual solubility, and herein the conservative parameter binding = aH  HY = aW  HY is used

327

to describe the hydrogel-lipid binding strength. The conservative parameter binding varies

328

from 30 to 85 k BT rc . A small conservative parameter suggests a strong hydrogel-lipid 19 / 31

ACS Paragon Plus Environment

Page 21 of 32 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

329

binding strength and vice versa. The bilayer membrane wrapping with different

330

hydrogel-lipid binding strengths is simulated, as depicted in Figure 6. The lateral

331

compression with different degrees of deformation is observed in the equilibrated

332

morphologies of the IHs. Intriguingly, Figure 6A1 reveals a collapsed structure under high

333

hydrogel-lipid binding strength (i.e., binding = 30k BT rc ). Comparatively, the 2D density map

334

shows an obvious corona at the edge of the spherical structure, as shown in Figure 6B1.

335

When at a moderated hydrogel-lipid binding strength (i.e., 45 k BT rc  binding  65 k BT rc ),

336

the hollow morphology is formed in the bilayer membrane (illustrated in Figure 6B2 and

337

B3). However, at a low hydrogel-lipid binding strength (i.e., binding =85k BT rc ), a spherical

338

encapsulated morphology is formed (illustrated in Figure 6A4), which can be confirmed by

339

the corresponding 2D density map of the IHs (illustrated in Figure 6B4). Compared with a

340

low hydrogel-lipid binding strength, the nanoscale mixing is identified under high

341

hydrogel-lipid binding strength. The sectional-view and top-view of the IHs are given in

342

Figure 7. From Figure 7A1-A2, we examine that the IH dissolves into the membrane

343

accompanying with special co-dissolution phenomena. However, in the case of medium

344

hydrogel-lipid binding strength, the deformable wrapping occurs, as depicted in Figure

345

7B1-B2.

20 / 31

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

346 347

Figure 6. Membrane wrapping under different hydrogel-lipid binding strengths. (A1-A4)

348

binding =30, 45, 65, 85 kBT rc , respectively. (B1-B4) 2D density map of IHs corresponding

349

to snapshot figures in A1-A4. The color bar represents the hydrogel density in the units of

350

 2 . (C) The evolution of total radius-of-gyration for the IHs. (D) evolution of contact angle

351

for the IHs embedded in the bilayer membrane. The equilibrium simulation time is 19.2 s.

352

21 / 31

ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32 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

353

Figure 7. Comparison of membrane wrapping morphologies under different hydrogel-lipid

354

binding strengths: (A1-A2) binding = 30 k BT rc and (B1-B2) the binding = 65 k BT rc .

355

Though different hydrogel-lipid binding strengths are simulated, an obvious structural

356

transition is also identified in Figure 6C. Moreover, high hydrogel-lipid binding strength

357

leads to an increased radius-of-gyration owing to the confined compressibility of the IHs. The

358

computed compression angle is displayed in Figure 6D, which indicates the IH morphology

359

transforms from the wetting compression state to adhesive encapsulation, as reflected in

360

Figure 6A1-A4. The phenomenon of the wetting compression state was also reported by

361

previous studies when the hydrogels are placed at the confinement space.67-69 Figure 8A

362

presents the evolution of structural compression of the IHs. We observe the large

363

hydrogel-lipid binding strength gives rise to short transition time. The evolution of wrapping

364

ratio for each IH is given in Figure 8B. It is demonstrated that the hydrogel-lipid binding

365

strength has a significant role in the wrapping time. However, the approximative scaling from

366

 ~ t1.294 to  ~ t1.819 is obtained at the monotonous increase region. Our simulation results

367

indicate that the hydrogel-lipid binding strength of binding  85k BT / rc could optimize the

368

wrapping time. With high hydrogel-lipid binding strength, a longer wrapping time is needed

369

because of the nanoscale mixing and favorable deformability during the wrapping process.

370

On comparison of different hydrogel-lipid binding strengths, severe structure compression

371

and deformation of IHs can be found at high hydrogel-lipid binding strength of

372

binding  30k BT / rc , as shown in Figure 8A. In fact, stronger deformability results in longer

373

wrapping time owing to higher energy barriers.18, 70 22 / 31

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

Page 24 of 32

374

Satarifard and Zhang et al.18, 59 have studied the effect of surface tension on membrane

375

wrapping. Herein, the stress distribution for the IHs is computed, as shown in Figure S1. The

376

stress distribution with a surface normal in the Y-direction inside the membrane is

377

characterized via the equation of ( X , Z )  PZZ  0.5g( PXX  PYY ) , where the PXX , PYY , PZZ ,

378

are the pressure tensor at the X, Y, and Z directions.18,

379

tension at the contact site implies that the bilayer membrane experiences a severe

380

heterogeneous deformation and favorable mechanical rupture.55 As shown in Figure S1, the

381

scattered distribution of stress demonstrates the interaction between the membrane and IHs

382

doesn’t readily make membrane rupture, and thus the membrane endocytosis induced by the

383

IHs has not been found.71 In addition, the affluence of surface tension on the membrane

384

wrapping is quantified through theoretical calculation, in which the relationship between the

385

surface tension and lateral compressed wrapping time tlcw is identified. The computation

386

details are presented in the SI. As a result, the lateral compressed wrapping time tlcw is

387

mainly dependent on the surface tension: tlcw ∼ 2 IH  IH , where  are the characteristic

388

radii of IHs,  IH is the viscosity of the IHs, and  IH is the surface tension of IHs. Moreover,

389

the surface tension also determines the adhesive morphology of polymer matrix. Generally, a

390

large membrane tension triggers an axisymmetric morphology in a polymer matrix, while a

391

small membrane tension induces a non-axisymmetric shape.18

59

392

23 / 31

ACS Paragon Plus Environment

A dramatic change of surface

Page 25 of 32

2 (A) 1.8

binding=30kBT/rc binding=45kBT/rc

1.6

binding=65kBT/rc binding=85kBT/rc

1.4 1.2

laterally compressed wrapping

1

0.1

393

1 Time/s

10 (B)

binding=30kBT/rc binding=45kBT/rc binding=65kBT/rc

1

Kslope1=1.294

binding=85kBT/rc

tio

0.8

vertically capillary adhesion

Wrapping ra

Structural compression

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

0.1 Kslope2=1.819

0.01

394

0.1 Time 1 /s

30kBT/rc 45kBT/rc 65kBT/rc 10

85kBT/rc

395

Figure 8. Membrane wrapping dynamics of the IHs with different hydrogel-lipid binding

396

strengths: (A) structural compression of IHs and (B) wrapping ratio of IHs.

397 398

Figure 9. Membrane wrapping pathway of IHs.  HS ,  MS ,  HM correspond to the surface 24 / 31

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

399

stress tensors of the biphasic interface, VHS , VMS , VHM indicate the tangent vectors of the

400

biphasic interface, and the subscript H , S , M represent hydrogel, solvent, membrane,

401

respectively.

402

4. Conclusions

403

With the help of large-scale DPD simulations, we have systematically investigated the

404

membrane wrapping process of IHs by considering different IH architectures and

405

hydrogel-lipid binding strengths. A theatrical structure transition from vertical capillary

406

adhesion to lateral compressed wrapping is discovered, which is significantly different from

407

the situation with non-elastic nanomaterials. Originating from the van der Waals attraction of

408

hydrogel-lipid, the formation of vertical capillary adhesion offers a driving force for the

409

membrane to wrap the corners of IHs (see Figure 9A). The conformational variation of the

410

IHs at the boundary of the membrane is mainly determined by a vector balance of the surface

411

stresses,72 as illustrated in Figure 9B. Correspondingly, the IHs gradually embed into the

412

membrane (see Figure 9C). However, the lateral compression inhibits the membrane

413

wrapping process due to the wrapping-induced shape deformation, and thus, the wrapping

414

time for vertical capillary adhesion stage is much shorter than that for lateral compressed

415

wrapping stage. All these results can be explained by the unique collaboration of

416

lastocapillary deformation of networked gels and nanoscale confinement of bilayer

417

membrane. Moreover, our results demonstrate the stuffed IH forms a hollow structure during

418

the wrapping process. Conversely, a collapsed architecture is discovered in the IH capsules

419

because of their strong compressibility and deformability compared with stuffed IH. In 25 / 31

ACS Paragon Plus Environment

Page 26 of 32

Page 27 of 32 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

420

addition, when at large hydrogel-lipid binding strength, the special nano-mixing occurs, and

421

this would result in a long wrapping time. We are currently extending our computational

422

approach to address these mechanistic membrane wrappings of IHs in a systematic manner as

423

will be discussed in subsequent work.

424

Conflicts of interest

425

The authors have no conflicts of interest to declare.

426

Acknowledgments

427

This work is supported by National Natural Science Foundation of China (Nos. 21878078,

428

and 21808056), by National Natural Science Foundation of China for Innovative Research

429

Groups (No. 51621002), and by the Shanghai Science and Technology Innovation Action

430

Plan (18160743700).

431 432

References

433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448

1.Segers, V. F. M.; Lee, R. T., Stem-cell therapy for cardiac disease. Nature 2008, 451, 937. 2.Daley, G. Q.; Scadden, D. T., Prospects for Stem Cell-Based Therapy. Cell 2008, 132 (4), 544-548. 3.Lü, S.; Gao, C.; Xu, X.; Bai, X.; Duan, H.; Gao, N.; Feng, C.; Xiong, Y.; Liu, M., Injectable and Self-Healing Carbohydrate-Based Hydrogel for Cell Encapsulation. Acs Applied Materials & Interfaces 2015, 7 (23), 13029-13037. 4.Yang, L.; Liu, A.; de Ruiter, M. V.; Hommersom, C. A.; Katsonis, N.; Jonkheijm, P.; Cornelissen, J. J. L. M., Compartmentalized supramolecular hydrogels based on viral nanocages towards sophisticated cargo administration. Nanoscale 2018, 10 (8), 4123-4129. 5.Zhang, Y. S.; Khademhosseini, A., Advances in engineering hydrogels. Science 2017, 356 (6337), eaaf3627. 6.Li, L.; Yan, B.; Yang, J.; Chen, L.; Zeng, H., Novel Mussel-Inspired Injectable Self-Healing Hydrogel with Anti-Biofouling Property. Advanced Materials 2015, 27 (7), 1294-1299. 7.Norouzi, M.; Nazari, B.; Miller, D. W., Injectable hydrogel-based drug delivery systems 26 / 31

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

449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490

for local cancer therapy. Drug Discovery Today 2016, 21 (11), 1835-1849. 8.Seliktar, D., Designing Cell-Compatible Hydrogels for Biomedical Applications. Science 2012, 336 (6085), 1124-1128. 9.Maiti, D.; Chao, Y.; Dong, Z.; Yi, X.; He, J.; Liu, Z.; Yang, K., Development of a thermosensitive protein conjugated nanogel for enhanced radio-chemotherapy of cancer. Nanoscale 2018, 10 (29), 13976-13985. 10.Liu, Y.; Ai, K.; Liu, J.; Deng, M.; He, Y.; Lu, L., Dopamine-Melanin Colloidal Nanospheres: An Efficient Near-Infrared Photothermal Therapeutic Agent for In Vivo Cancer Therapy. Advanced Materials 2013, 25 (9), 1353-1359. 11.Jin, H.; Liu, X.; Gui, R.; Wang, Z., Facile synthesis of gold nanorods/hydrogels core/shell nanospheres for pH and near-infrared-light induced release of 5-fluorouracil and chemo-photothermal therapy. Colloids and Surfaces B: Biointerfaces 2015, 128, 498-505. 12.Yue, T.; Zhou, H.; Sun, H.; Li, S.; Zhang, X.; Cao, D.; Yi, X.; Yan, B., Why are nanoparticles trapped at cell junctions when the cell density is high? Nanoscale 2019, 11 (14), 6602-6609. 13.Zhou, M.; Liu, S.; Jiang, Y.; Ma, H.; Shi, M.; Wang, Q.; Zhong, W.; Liao, W.; Xing, M. M. Q., Doxorubicin-Loaded Single Wall Nanotube Thermo-Sensitive Hydrogel for Gastric Cancer Chemo-Photothermal Therapy. Advanced Functional Materials 2015, 25 (29), 4730-4739. 14.Zhu, X.; Zhang, Y.; Huang, H.; Zhang, H.; Hou, L.; Zhang, Z., Functionalized graphene oxide-based thermosensitive hydrogel for near-infrared chemo-photothermal therapy on tumor. Journal of Biomaterials Applications 2016, 30 (8), 1230-1241. 15.Yue, T.; Zhang, X., Cooperative Effect in Receptor-Mediated Endocytosis of Multiple Nanoparticles. Acs Nano 2012, 6 (4), 3196-3205. 16.Li, S.; Luo, Z.; Xu, Y.; Ren, H.; Deng, L.; Zhang, X.; Huang, F.; Yue, T., Interaction pathways between soft lipid nanodiscs and plasma membranes: A molecular modeling study. Biochimica et Biophysica Acta (BBA) - Biomembranes 2017, 1859 (10), 2096-2105. 17.Hartmann, R.; Weidenbach, M.; Neubauer, M.; Fery, A.; Parak, W. J., Stiffness-Dependent In Vitro Uptake and Lysosomal Acidification of Colloidal Particles. Angewandte Chemie International Edition 2015, 54 (4), 1365-1368. 18.Satarifard, V.; Grafmüller, A.; Lipowsky, R., Nanodroplets at Membranes Create Tight-Lipped Membrane Necks via Negative Line Tension. Acs Nano 2018, 12 (12), 12424-12435. 19.Shen, Z.; Ye, H.; Yi, X.; Li, Y., Membrane Wrapping Efficiency of Elastic Nanoparticles during Endocytosis: Size and Shape Matter. Acs Nano 2019, 13 (1), 215-228. 20.Thakur, A.; Jaiswal, M. K.; Peak, C. W.; Carrow, J. K.; Gentry, J.; Dolatshahi-Pirouz, A.; Gaharwar, A. K., Injectable shear-thinning nanoengineered hydrogels for stem cell delivery. Nanoscale 2016, 8 (24), 12362-12372. 21.Liu, J.; Qi, C.; Tao, K.; Zhang, J.; Zhang, J.; Xu, L.; Jiang, X.; Zhang, Y.; Huang, L.; Li, Q.; Xie, H.; Gao, J.; Shuai, X.; Wang, G.; Wang, Z.; Wang, L., Sericin/Dextran Injectable 27 / 31

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32 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

491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532

Langmuir

Hydrogel as an Optically Trackable Drug Delivery System for Malignant Melanoma Treatment. Acs Applied Materials & Interfaces 2016, 8 (10), 6411-6422. 22.Qi, C.; Liu, J.; Jin, Y.; Xu, L.; Wang, G.; Wang, Z.; Wang, L., Photo-crosslinkable, injectable sericin hydrogel as 3D biomimetic extracellular matrix for minimally invasive repairing cartilage. Biomaterials 2018, 163, 89-104. 23.Lei, K.; Ma, Q.; Yu, L.; Ding, J., Functional biomedical hydrogels for in vivo imaging. Journal of Materials Chemistry B 2016, 4 (48), 7793-7812. 24.Li, J.; Mooney, D. J., Designing hydrogels for controlled drug delivery. Nature Reviews Materials 2016, 1, 16071. 25.Song, X.; Zhao, S.; Fang, S.; Ma, Y.; Duan, M., Mesoscopic Simulations of Adsorption and Association of PEO-PPO-PEO Triblock Copolymers on a Hydrophobic Surface: From Mushroom Hemisphere to Rectangle Brush. Langmuir 2016, 32 (44), 11375-11385. 26.Gartner, T. E.; Jayaraman, A., Modeling and Simulations of Polymers: A Roadmap. Macromolecules 2019, 52 (3), 755-786. 27.Shen, Z.; Loe, D. T.; Awino, J. K.; Kröger, M.; Rouge, J. L.; Li, Y., Self-assembly of core-polyethylene glycol-lipid shell (CPLS) nanoparticles and their potential as drug delivery vehicles. Nanoscale 2016, 8 (31), 14821-14835. 28.Li, D.-W.; Liu, X. Y.; Feng, Y. P., Bond-Angle-Potential-Dependent Dissipative Particle Dynamics Simulation and Lipid Inverted Phase. The Journal of Physical Chemistry B 2004, 108 (30), 11206-11213. 29.Yan, Z.; Li, S.; Luo, Z.; Xu, Y.; Yue, T., Membrane nanotube pearling restricted by confined polymers. Soft Matter 2018, 14 (46), 9383-9392. 30.Tian, F.; Yue, T.; Dong, W.; Yi, X.; Zhang, X., Size-dependent formation of membrane nanotubes: continuum modeling and molecular dynamics simulations. Physical Chemistry Chemical Physics 2018, 20 (5), 3474-3483. 31.Shillcock, J. C.; Lipowsky, R., Equilibrium structure and lateral stress distribution of amphiphilic bilayers from dissipative particle dynamics simulations. The Journal of Chemical Physics 2002, 117 (10), 5048-5061. 32.Laradji, M.; Sunil Kumar, P. B., Dynamics of Domain Growth in Self-Assembled Fluid Vesicles. Physical Review Letters 2004, 93 (19), 198105. 33.Guigas, G.; Morozova, D.; Weiss, M., Chapter 4 - Exploring Membrane and Protein Dynamics with Dissipative Particle Dynamics. In Advances in Protein Chemistry and Structural Biology, Christov, C., Ed. Academic Press: 2011; Vol. 85, pp 143-182. 34.Li, Z.; Tang, Y.-H.; Li, X.; Karniadakis, G. E., Mesoscale modeling of phase transition dynamics of thermoresponsive polymers. Chemical Communications 2015, 51 (55), 11038-11040. 35.Pelton, R., Temperature-sensitive aqueous microgels. Advances in Colloid and Interface Science 2000, 85 (1), 1-33. 36.Song, X.; Bao, B.; Tao, J.; Zhao, S.; Han, X.; Liu, H., Deswelling Dynamics of Thermoresponsive Microgel Capsules and Their Ultrasensitive Sensing Applications: A Mesoscopic Simulation Study. The Journal of Physical Chemistry C 2019, 123 (3), 1828-1838. 28 / 31

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

533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574

37.Song, X.; Qiao, C.; Tao, J.; Bao, B.; Han, X.; Zhao, S., Interfacial Engineering for Thermoresponsive Microgel Capsules: Polymeric Wetting versus Colloidal Adhesion. Macromolecules 2018, 52 (10), 3869-3880. 38.Li, S.; Yan, Z.; Luo, Z.; Xu, Y.; Huang, F.; Hu, G.; Zhang, X.; Yue, T., Directional and Rotational Motions of Nanoparticles on Plasma Membranes as Local Probes of Surface Tension Propagation. Langmuir 2019, 35 (15), 5333-5341. 39.Xiang, W. J.; Zhao, S. L.; Song, X. Y.; Fang, S. W.; Wang, F.; Zhong, C.; Luo, Z. Y., Amphiphilic nanosheet self-assembly at the water/oil interface: computer simulations. Physical Chemistry Chemical Physics 2017, 19 (11), 7576-7586. 40.Masoud, H.; Alexeev, A., Controlled release of nanoparticles and macromolecules from responsive microgel capsules. Acs Nano 2011, 6 (1), 212-219. 41.Mura, S.; Nicolas, J.; Couvreur, P., Stimuli-responsive nanocarriers for drug delivery. Nature materials 2013, 12 (11), 991. 42.Mei, A.; Guo, X.; Ding, Y.; Zhang, X.; Xu, J.; Fan, Z.; Du, B., PNIPAm-PEO-PPO-PEO-PNIPAm Pentablock Terpolymer: Synthesis and Chain Behavior in Aqueous Solution. Macromolecules 2010, 43 (17), 7312-7320. 43.Huo, J.; Jiang, H.; Chen, Z.; Zhou, J., Homoporous polymer membrane via forced surface segregation: A computer simulation study. Chemical Engineering Science 2018, 191, 490-499. 44.Ding, H.-m.; Ma, Y.-q., Interactions between Janus particles and membranes. Nanoscale 2012, 4 (4), 1116-1122. 45.Lv, K.; Li, Y., Indentation of Graphene-Covered Atomic Force Microscopy Probe Across a Lipid Bilayer Membrane: Effect of Tip Shape, Size, and Surface Hydrophobicity. Langmuir 2018, 34 (26), 7681-7689. 46.Schmidt, S.; Hellweg, T.; von Klitzing, R., Packing Density Control in P(NIPAM-co-AAc) Microgel Monolayers: Effect of Surface Charge, pH, and Preparation Technique. Langmuir 2008, 24 (21), 12595-12602. 47.Feng, C.; Shen, Z.; Li, Y.; Gu, L.; Zhang, Y.; Lu, G.; Huang, X., PNIPAM-b-(PEA-g-PDMAEA) double-hydrophilic graft copolymer: Synthesis and its application for preparation of gold nanoparticles in aqueous media. Journal of Polymer Science Part A: Polymer Chemistry 2009, 47 (7), 1811-1824. 48.Scherzinger, C.; Lindner, P.; Keerl, M.; Richtering, W., Cononsolvency of Poly(N,N-diethylacrylamide) (PDEAAM) and Poly(N-isopropylacrylamide) (PNIPAM) Based Microgels in Water/Methanol Mixtures: Copolymer vs Core−Shell Microgel. Macromolecules 2010, 43 (16), 6829-6833. 49.Bai, X.; Xu, M.; Liu, S.; Hu, G., Computational Investigations of the Interaction between the Cell Membrane and Nanoparticles Coated with a Pulmonary Surfactant. Acs Applied Materials & Interfaces 2018, 10 (24), 20368-20376. 50.Ayton, G.; Voth, G. A., Bridging Microscopic and Mesoscopic Simulations of Lipid Bilayers. Biophysical Journal 2002, 83 (6), 3357-3370. 51.Shillcock, J. C.; Lipowsky, R., Tension-induced fusion of bilayer membranes and vesicles. Nature materials 2005, 4, 225. 29 / 31

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32 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

575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616

Langmuir

52.Khanal, A.; Bui, M.-P. N.; Seo, S. S., Microgel-encapsulated methylene blue for the treatment of breast cancer cells by photodynamic therapy. Journal of breast cancer 2014, 17 (1), 18-24. 53.Senff, H.; Richtering, W., Temperature sensitive microgel suspensions: Colloidal phase behavior and rheology of soft spheres. The Journal of Chemical Physics 1999, 111 (4), 1705-1711. 54.Wong, J. E.; Gaharwar, A. K.; Müller-Schulte, D.; Bahadur, D.; Richtering, W., Dual-stimuli responsive PNiPAM microgel achieved via layer-by-layer assembly: Magnetic and thermoresponsive. Journal of Colloid and Interface Science 2008, 324 (1-2), 47-54. 55.Plimpton, S.; Crozier, P.; Thompson, A., LAMMPS-large-scale atomic/molecular massively parallel simulator. Sandia National Laboratories 2007, 18, 43. 56.Champion, J. A.; Mitragotri, S., Role of target geometry in phagocytosis. Proceedings of the National Academy of Sciences of the United States of America 2006, 103 (13), 4930-4934. 57.Yi, X.; Shi, X.; Gao, H., Cellular Uptake of Elastic Nanoparticles. Physical Review Letters 2011, 107 (9), 098101. 58.Takechi-Haraya, Y.; Goda, Y.; Sakai-Kato, K., Control of Liposomal Penetration into Three-Dimensional Multicellular Tumor Spheroids by Modulating Liposomal Membrane Rigidity. Molecular pharmaceutics 2017, 14 (6), 2158-2165. 59.Zhang, L.; Zhao, Y.; Wang, X., Nanoparticle-Mediated Mechanical Destruction of Cell Membranes: A Coarse-Grained Molecular Dynamics Study. Acs Applied Materials & Interfaces 2017, 9 (32), 26665-26673. 60.Zhang, L.; Becton, M.; Wang, X., Designing Nanoparticle Translocation through Cell Membranes by Varying Amphiphilic Polymer Coatings. The Journal of Physical Chemistry B 2015, 119 (9), 3786-3794. 61.Aoyagi, T.; Takimoto, J.-i.; Doi, M., Molecular dynamics study of polymer melt confined between walls. The Journal of Chemical Physics 2001, 115 (1), 552-559. 62.Sagbas, S.; Sahiner, N., Modifiable natural gum based microgel capsules as sustainable drug delivery systems. Carbohydrate Polymers 2018, 200, 128-136. 63.Seiffert, S.; Thiele, J.; Abate, A. R.; Weitz, D. A., Smart Microgel Capsules from Macromolecular Precursors. Journal of the American Chemical Society 2010, 132 (18), 6606-6609. 64.Cheng, S. Y.; Chou, G.; Buie, C.; Vaughn, M. W.; Compton, C.; Cheng, K. H., Maximally asymmetric transbilayer distribution of anionic lipids alters the structure and interaction with lipids of an amyloidogenic protein dimer bound to the membrane surface. Chemistry and Physics of Lipids 2016, 196, 33-51. 65.Yue, T.; Zhang, X., Molecular understanding of receptor-mediated membrane responses to ligand-coated nanoparticles. Soft Matter 2011, 7 (19), 9104-9112. 66.Yue, T.; Wang, X.; Huang, F.; Zhang, X., An unusual pathway for the membrane wrapping of rodlike nanoparticles and the orientation-and membrane wrapping-dependent nanoparticle interaction. Nanoscale 2013, 5 (20), 9888-9896. 30 / 31

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

617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636

67.Gumerov, R. A.; Rumyantsev, A. M.; Rudov, A. A.; Pich, A.; Richtering, W.; Möller, M.; Potemkin, I. I., Mixing of Two Immiscible Liquids within the Polymer Microgel Adsorbed at Their Interface. Acs Macro Letters 2016, 5 (5), 612-616. 68.Lee, M. W.; Han, M.; Bossa, G. V.; Snell, C.; Song, Z.; Tang, H.; Yin, L.; Cheng, J.; May, S.; Luijten, E.; Wong, G. C. L., Interactions between Membranes and “Metaphilic” Polypeptide Architectures with Diverse Side-Chain Populations. Acs Nano 2017, 11 (3), 2858-2871. 69.Gumerov, R. A.; Rudov, A. A.; Richtering, W.; Möller, M.; Potemkin, I. I., Amphiphilic Arborescent Copolymers and Microgels: From Unimolecular Micelles in a Selective Solvent to the Stable Monolayers of Variable Density and Nanostructure at a Liquid Interface. Acs Applied Materials & Interfaces 2017, 9 (37), 31302-31316. 70.Li, Y.; Zhang, X.; Cao, D., Nanoparticle hardness controls the internalization pathway for drug delivery. Nanoscale 2015, 7 (6), 2758-2769. 71.Mihut, A. M.; Dabkowska, A. P.; Crassous, J. J.; Schurtenberger, P.; Nylander, T., Tunable Adsorption of Soft Colloids on Model Biomembranes. Acs Nano 2013, 7 (12), 10752-10763. 72.Style, R. W.; Jagota, A.; Hui, C.-Y.; Dufresne, E. R., Elastocapillarity: Surface Tension and the Mechanics of Soft Solids. Annual Review of Condensed Matter Physics 2017, 8 (1), 99-118.

31 / 31

ACS Paragon Plus Environment

Page 32 of 32