- Email: [email protected]
Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles
CSBJ-00275; No of Pages 8 Computational and Structural Biotechnology Journal xxx (xxxx) xxx
Contents lists available at ScienceDirect
journal homepage: www.elsevier.com/locate/csbj
4 5 6
a
7
a r t i c l e
8 9 10 11 12 17 16 15 14 13 37 38 39 40 41 42 43
Article history: Received 29 September 2018 Received in revised form 20 December 2018 Accepted 25 December 2018 Available online xxxx
Department of Civil Engineering, National Taiwan University, Taipei, Taiwan Inistitute of Polymer Science and Engineering, National Taiwan University, Taipei, Taiwan
i n f o
a b s t r a c t
Biodegradable hydrogels have become promising materials for many biological applications in the past years. Recently, novel waterborne biodegradable polyurethane (WDPU) nanoparticles have been synthesized by a green water-based process, and serve as fundamental building blocks to form materials with great biocompatibility, biodegradability, and mechanical properties. However, the molecular structures and mechanisms of the WDPU nanoparticles and the relationship between the chemical compositions of the polymer segments and the material properties of the biodegradable hydrogels at macro-scale are still not well understood. In this study, we explore the fundamental mechanisms of WDPU nanoparticles through a full atomistic simulation approach to understand how the chemical compositions at the molecular level affect the molecular structures and material properties of WDPU nanoparticles. Specifically, we compare two WDPUs, i.e. PCL75LL25 and PCL75DL25, of the same hard segment composition and very similar soft segment composition [75% poly(e-caprolatone) and 25% polylactide], except the lactide in the former is L-form and in the latter is D,L-form. Our results show that the material properties of the biodegradable hydrogel can be designed by tuning the chemical compositions of the polymer segments. We find that the PCL75DL25 and PCL75LL25 have distinct molecular structures and physical crosslinks within the nanoparticles. The molecular structure of WDPU with PDLLA as soft segments is more extended, leading to more physical crosslinks between PCL segments. This study provide fundamental insights into the molecular structures and mechanisms of WDPU nanoparticles and help enabling the design of material properties of biocompatible hydrogel. © 2018 Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
D
b
F
Chien-Hui Wen a, Shun-Chieh Hsu a,b, Shan-hui Hsu b,⁎, Shu-Wei Chang a,⁎
O
Q23Q1
R O
2
Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles
P
1
E
C
T
E
Keywords: Molecular dynamics Polymer Hydrogel Nanoparticle biomaterial green material
R
36
52 53 54 55 56 57 58 59 60 61 62 63 64
In the past few decades, the development and applications of hydrogel have attracted increasing interest, specifically as a new class of biomaterials [1–8]. Hydrogels, namely hydrophilic gels, are networks of polymer chains that are extensively swollen with water [1,9–14]. Among hydrogels, biodegradable hydrogels are of utmost interest because they exhibit similar mechanical properties as natural soft tissues and can be degraded in an aqueous environment after their useful lifetime. Biodegradable hydrogels include natural hydrogels and synthetic type. Natural hydrogels have been gradually replaced by synthetic biodegradable hydrogels because of their strong water absorption capacity, long service life time, and a wide range of available chemical materials [10]. Recently, a novel waterborne biodegradable polyurethane (WDPU) has been synthesized and shown to have great potential in biomedical applications [15–20]. It is synthesized by a green water-based process, and has great biocompatibility, biodegradability, and mechanical
N C
50 51
1. Introduction
U
49
O R
47 45 44 46 48
18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
⁎ Corresponding authors. E-mail addresses: [email protected] (S. Hsu), [email protected] (S.-W. Chang).
properties. The WDPU has been demonstrated that it can be synthesized using biodegradable polyesters as the soft segment without chemical crosslinking [18]. The properties of the WDPU nanoparticle could be adjusted by the chemical composition and ratio of the polyesters in the soft segment [15], such as PCL diol, ad PLA diol. The behavior of thermally induced self-assembly and gel formation and the degradation rate of the hydrogel can be tuned by changing the soft segments of WDPU, indicating that WDPU hydrogels with a various of material properties can be designed by changing the composition of soft segments for different biomedical applications. It has also been shown that the WDPU hydrogels are able to be applied to 3D printing for neutral stem cells (NSCs) carrier, and shown great potential in central nervous system repair [16]. The integration of biodegradable hydrogel and 3D printing technology has opened great opportunities for the design of smart biocompatible scaffolds for many applications due to the ability to access complex internal structures [15–20]. WDPU nanoparticles are the fundamental building blocks of the WDPU hydrogels. The molecular structures of WDPU nanoparticles play an important role in the molecular mechanisms of sol-gel transition and the material properties of the WDPU hydrogel. However, the molecular structures and mechanisms of the WDPU and the
https://doi.org/10.1016/j.csbj.2018.12.007 2001-0370/© 2018 Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85
94
2. Materials and Methods
95
2.1. Simulation Details
96
2.1.1. Model Construction Fig. 1(a) illustrates the polymer chains constructed in this work. Three different soft segments are used in this study: PCL diol, poly(D,Llactide) diol (PDLLA diol, Mn ~ 2000 Da) and poly(L-lactide) diol (PLLA diol, Mn ~ 2000). The WDPU with PCL diol as both soft segments is
P D E T C E R R
99 100
O
97 98
C
92
N
90 91
U
88 89
named WDPUPCL. The WDPU with PDLLA diol as both soft segments is named WDPUPDLA. The WDPU with PLLAdiol as both soft segments is named WDPUPLLA. Each WDPU polymer chain is consisted of 4 repeating units with two soft segments. PCL and PDLLA are mixed in 3:1 M ratio to construct the PCL75DL25 nanoparticle while PCL and PLLA are mixed in 3:1 M ratio to construct the PCL75LL25 nanoparticle (Fig. 1(a)). We use Materials Studio to construct WDPU monomers according to the previously published synthesis procedure [15]. Each chain contains 3070 atoms and the molecular weight is about 39,018 Da. We use Amorphous Cell to construct the initial amorphous cells for PCL75DL25 and PCL75LL25 as shown in Fig. 1(b). Each initial unit cell consists 10 chains of polymers. To equilibrate the amorphous structure, we firstly run NPT simulation at 500 K for 0.6 ns to speed up the process, then cool down the temperature to 300 K for 0.2 ns to reach equilibrium at the room temperature. The equilibrated unit cell is then replicated in all three directions to construct a larger amorphous packing nanoparticle (Fig. 1(c))
R O
93
relationship between the chemical compositions of the polymer segments and the material properties of the biodegradable hydrogels at macro-scale are still not clear. In this paper, we employed full atomistic simulations to explore the molecular structure and fundamental mechanisms of WDPU nanoparticles. WDPU nanoparticles with 75% PCL diol and 25% PLA diol soft segments (L-form or D,L-form) are simulated to investigate the relationships between the chemical composition of soft segments and the physical properties of the WDPU nanoparticles.
F
86 87
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
O
2
Fig. 1. (a) Chemical composition of the polymer chains. (b) WDPU polymer chains with 75% PCL and 25% PLA are constructed. (c) Illustration of the amorphous unit cells.
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143
h
i
1−e−αðb−b0 Þ þ
2.2.3. Surface Area Solvent-accessible surface area of atoms in the nanoparticle is calculated by using the assigned radius for each atom with an extended radius for 1.4 Å to find the points on a sphere that are exposed to solvent. To normalize the surface area, we calculate the surface area of WDPU nanoparticles (Am), and divided by surface area when the nanoparticle is spherical (As). Values close to one imply a higher spherical symmetry while a larger value indicates that the nanoparticle has unsmooth surface or irregular shape. We compute the average of the surface area of the last 0.1 ns after the simulations of 2 ns. 2.2.4. End to End Distance End to end distance is a parameter to describe the conformation of polymer chains. To calculate the end to end distance of polymer chains, we define the atomic positions of the Nitrogen atoms on N-C-C-N as the ends of a segment, and the distance between these two atomic positions are considered as the end to end distance of segments. 3. Results
X
H θ ðθ−θ0 Þ2 X XX þ H ϕ ½1 þ s cosðnϕÞ þ Hχ χ 2 þ F bb0 ðb−b0 Þ χ b X b0 X X X ϕ0 0 F θθ0 ðθ−θ0 Þ θ0 −θ00 þ F bθ ðb−b0 Þ b −b0 þ 0 θ bXθX θ X 0 ðθ−θ0 Þ þ F ϕθθ0 cosϕðθ−θ0 Þ θ −θ00 þ F χχ0 χχ 0 Db bX
166
θ
151 152
2.2.1. Radius of Gyration Radius of gyration, Rg, is a useful representation to measure and characterize the size of complex polymer configurations. We compute the average of the radius of gyration of the last 0.1 ns after the simulations of 2 ns. We use Eq. (2) to calculate the radius of gyration of ! WDPU nanoparticles. In Eq. (2), r i is the position of the i-th atom and ! r cm is the center of mass of the WDPU nanoparticle.
N C
149 150
ð1Þ
O R
2.2. Analysis of Molecular Structures
147 148
χ0
R
145
146
χ
" ϕ 6 # X qi q j r 12 r þ ε −2 r r εr ij
E
þ
X
C
T
Epot ¼
X
U
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi u N u 1 X ! 2 r i −! r cm Rg ¼ t N i¼1
ð2Þ
154
155 156 157 158 159 160 161 162 163
ð3Þ
F
124
Imin Iaverage
2.2.2. Eccentricity To understand if the differences on the radius of gyration are resulted from the shape of nanoparticles, we calculate the eccentricity (ρeccen) of the nanoparticles to investigate the structures of nanoparticles. The eccentricity measures the symmetry of WDPU nanoparticles. If the value is close to zero, it implies a high symmetry, while if the value is larger, it implies less symmetry. We use Eq. (3) to calculate the eccentricity of WDPU nanoparticles [26]. In Eq. (3), Imin is the minimum value of the tri-axial moment of inertia and Iavg is the average value of the tri-axial moment of inertia. We compute the average of
167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182
3.1. Nanoparticle Characterization
183
To characterize the structures of WDPU nanoparticles, we analyze the radius of gyrations, eccentricity, and surface area of PCL75DL25 and PCL75LL25 nanoparticles. Fig. 2 (a) shows the results. We find that the PCL75LL25 nanoparticle has a radius of gyration with a value of 70.01 ± 0.03 Å, while the PCL75DL25 nanoparticle has a slightly smaller radius of gyration of 68.90 ± 0.03 Å. The results of the eccentricity show that the PCL75DL25 nanoparticle has a value of ρeccen = 0.039 ± 0.004 and the value of the PCL75LL25 nanoparticle is 0.079 ± 0.004, suggesting that the shape of the PCL75LL25 nanoparticle is more irregular. For the surface area, there is no significant differences on the surface area between the PCL75DL25 and the PCL75LL25 nanoparticles. It is worth noting that these results show that both PCL75DL25 and PCL75LL25 have unsmoothed surface and the PCL75LL25 nanoparticle has a more irregular and unsymmetrical shape. Since both polymers have the same molecular weights, the differences on the radius of gyration might be resulted from either the change on the packing of molecules within the nanoparticles or the change on the surface morphology. We analyze the density distributions of both nanoparticles as shown in Fig. 2 (b). In Fig. 2 (b), we find that both nanoparticles have the same core density when the radius is smaller than 50 Å. These results indicate that there is no difference on the density in the center between different WDPU nanoparticles, while there is significant differences on the surface morphology between different WDPU nanoparticles. From the distance of 55 Å, the density of both PCL75DL25 and PCL75LL25 nanoparticles start to decrease indicating that the packing of polymer chains on the surfaces of both nanoparticles are much looser than the molecular packing in the center of the nanoparticles. Although both nanoparticles have a radius of gyration ~70 Å, the density of PCL75DL25 and PCL75LL25 decrease slowly and goes to 0 when the radius exceeds 110 Å. It is worth noting that from 50 to 80 Å, the PCL75DL25 has a larger density while a reverse trend is found for the distance from 80 to 120 Å. The density of PCL75DL25 reaches 0 at a distance smaller than the PCL75LL25. These results suggest that the selfassembly of PCL75LL25 nanoparticle is more loosely on the surface when compare with the PCL75DL25 nanoparticle. That is, the slightly larger radius gyration of the PCL75LL25 is not resulted from different packing density in the center of the nanoparticle but resulted from the differences on the surface morphology of the WDPU nanoparticles. We
184 185
E
122 123
P eccen ¼ 1−
O
2.1.2. Force Field We use consistent valence force field (CVFF) [21] to simulate the molecular structure of the WDPU nanoparticles. The CVFF originally applied to biological system [22] and is a prototype of the consistent force field (CFF) and its later derivatives (the polymer consistent force field PCFF and COMPASS). Terms in CVFF included a bond stretching or compression, an angle bending, torsion angle twisting and out-of-plane deformation of a planar system, five cross-coupling terms between bond deformations, angle bending, a bond deformation and an angle bending, a torsion angle twisting and the two associated angle bending, and coupling between out-of-plane deformations, and repulsive, dispersive, and Coulombic interactions between nonbonded atoms [23]. The energy description is shown in Eq. (1). The first four terms in Eq. (1) included a Morse potential for bond stretching and compression (b), a harmonic term for angle bending (θ), cosine torsional (ϕ) and out-ofplane torsional terms (χ). Db, Hθ, Hϕ, and Hχ are the force constants for the corresponding intramolecular deformations. b0 and θ0 are the equilibrium bond length and angle. The next five terms are cross terms that account for interactions between the four types of internal coordinates. Fij are the force constants for the cross terms. The last two terms represents Lennard-Jones and Coulombic terms for nonbonded interaction. ε and r∗ are the parameters for the nonbonded repulsive and dispersive interactions, and qi are the partial charges on each atom [23–25].
R O
121
164
the eccentricity of the last 0.1 ns after the simulations of 2 ns.
P
to simulate the self-assembly of the nanoparticle in a larger box in vacuum. The final model consists of 245,600 atoms.
D
119 120
3
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221
4
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
(a) Eccentricity
70.5
0.1
70
0.08
69.5 69 68.5
Surface area 2.54 2.52 2.5
0.06
A m/As
Eccentricity
0.04
2.42
0
2.4
PCL75DL25 PCL75LL25
PCL75DL25
PCL75DL25
PCL75LL25
R O
2 1.8 1.6
1.4 1.2
P
Density
PCL75LL25
F
PCL75DL25 PCL75LL25
(b)
2.46 2.44
0.02
68
2.48
O
Radius of gyration (Å)
Radius of gyration
1
0.8
D
0.6
10
30
50
70
T
0.2 0 -10
E
0.4
90
110
130
150
C
distance from the center of particle d (Å)
222 223
R R
E
Fig. 2. Nanoparticle characterizations. (a) Radius of gyration, eccentricity and surface area of WDPN nanoparticles. The error bars are the standard deviations from the ensemble average of the last 0.1 ns simulations. (b) Analysis of the density distribution of WDPU nanoparticles. The results show that the PCL75LL25 nanoparticle has a more irregular shape. There is no differences on the core density between the PCL75DL25 and PCL75LL25 nanoparticles but there is significant differences on the surface morphology when the PLLA soft segment is replaced with PDLLA soft segments. The molecular packing on the surface of nanoparticle is less dense for the PCL75DL25 nanoparticle.
231
3.2. Molecular Structure of WDPU Nanoparticles
232
The characterizations of WDPU nanoparticles show that the PCL75LL25 nanoparticle has smaller density in the surface layer and higher value of eccentricity. The packing and self-assembly of the nanoparticles are determined by the differences on the properties of each chain and the interactions between polymer chains. To understand the molecular structural differences between the PCL75DL25 and PCL75LL25 nanoparticles, we further analyze the end to end distances of polymer chains as shown in Fig. 3 (a). We find that although there is not much difference on the characterization of the shape between the PCL75DL25 and PCL75LL25 nanoparticles, there are significant differences at the molecular level between the PCL75DL25 and PCL75LL25 nanoparticles. The end to end distances of the polymer chains in the PCL75DL25 nanoparticles (99.976 ± 0.134 Å) is almost two times larger than the end to end distance of the polymer chains in
233 234 235 236 237 238 239 240 241 242 243 244 245
C
228 229
N
226 227
U
224 225
O
230
anticipate that the loose surface packing on both nanoparticles play important roles in the sol-gel transition. The loose surface packing allows polymer chains on the surface to move more freely and interact with other nanoparticles in solutions at higher temperature to initiate the physical crosslinks between nanoparticles and eventually form networks and further initiate the hydrogel formation. The larger radius of gyration and the looser packing on the surface of the PCL75LL25 nanoparticle provide molecular insights into its faster sol-gel transition time observed in experimental results.
the PCL75LL25 nanoparticles (54.990 ± 0.228 Å). It is worth noting that all the polymer chains in both WDPU nanoparticles have similar contour length, indicating that the polymer chains in the PCL75DL25 nanoparticle is much stiffer than the polymer chains in the PCL75LL25 nanoparticle. These results suggest that PLLA soft segments significantly cause shorter end to end distances of polymer chains while the PDLLA soft segments do not have the same effects. Since all the polymer chains in both nanoparticles have 75% PCL segments and 25% PLA segments (PDLLA soft segments in PCL75DL25 and PLLA soft segments in PCL75LL25), we further analyze the segment to segment distances of PCL segments and PLA segments in both nanoparticles as shown in Fig. 3 (b). The distributions of the segment to segment distances are shown in Fig. 3 (c) and (d). The results reveal that there are significant differences on the segment to segment distances of the PCL soft segments. The PLLA soft segments significantly cause shorter segment to segment distances of PCL soft segments in the PCL75LL25 nanoparticle while the PCL soft segments have much larger segment to segment distances in the PCL75DL25 nanoparticle. These results suggest that PDLLA increases the stiffness of the PCL soft segments in the PCL75DL25 nanoparticle. The shorter end to end distances, i.e. lower stiffness, for the PCL75LL25 nanoparticle suggest that when solvating in the solution, the polymer chains are easier to reach out and form physical crosslinks with other nanoparticles. We further analyze the radial distribution function (RDF) between different polymer chains to study the packing behaviors of the PCL75DL25 and PCL75LL25 nanoparticles. Fig. 4 (a) shows the RDF
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
(b) 90
120
Segment to Segment distance (Å)
100 80 60 40 20
80 70 50 40 30 20 10 0
0 PCL75DL25
(c)
PCL
PCL75LL25
O
DL LL
30
R O
25 Numbers
PLA
PLA Segment to Segment distance distribution (dPLA) 35
20 15 10
P
5 0 0
10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170
E
DL LL
D
PCL Segment to Segment distance distribution (dPCL)
(d) 25
T
20 15 10
C
Numbers
DL LL
60
F
End to end distance (Å)
(a)
5
E
5 0
10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170
R
0
276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293
N C
274 275
between the polymer chains. The peak of the RDF occurs at 17.75 Å for the PCL75DL25 nanoparticle while at 16.75 Å for the PCL75LL25 nanoparticle, suggesting the packing of both nanoparticles are similar from the point of view of the polymer chains. It is worth noting that the RDF of the PCL75DL25 nanoparticle drops to zero at a shorter distance when compared with the RDF of the PCL75LL25 nanoparticle. These results are consistent with the analysis of the density distribution of the nanoparticles as shown in Fig. 2 (b). Although there is no significant difference on the RDF of the polymer chains, we find that there is significant difference if we analyze the RDF of the PCL and PLA soft segments separately for the nanoparticles. Fig. 4 (b) and (c) show the RDF of the PCL soft segments and PLA soft segments. The RDF of the PCL soft segments are similar for both nanoparticles while there is a significant difference on the RDF of the PLA soft segments as shown in Fig. 4 (c). The first peak of the RDF of the PLA soft segments occurs at 8.25 Å for the PCL75DL25 nanoparticle while at 23.75 Å for the PCL75LL25 nanoparticle. Together with the results from the analysis of the end to end distances and segment to segment distances as shown in Fig. 3, our results reveal that the PCL and PLA soft segments are more extended in the PCL75DL25 nanoparticle resulted in a packing of polymer chains with shorter distances between PLA soft segments in the PCL75DL25 nanoparticle.
U
272 273
O R
Fig. 3. Analysis of the end to end distance of polymer chains in WDPU nanoparticles. (a) The end to end distances of polymer chains in WDPU nanoparticles. (b) PCL and PLA segment to segment distances in WDPU nanoparticles. The error bars are the standard deviations from the ensemble average of the last 0.1 ns simulations. (c) the distribution of the end to end distances of PLA segments in WDPU nanoparticles. (d) the distribution of the end to end distances of PCL segments in WDPU nanoparticles.
3.3. Intermolecular and Intramolecular Hydrogen Bonds in WDPU 294 Nanoparticles 295 The analysis of hydrogen bonds in the WDPU nanoparticles are shown in Fig. 5 to investigate the differences on the molecular interactions between the PCL75DL25 and PCL75LL25 nanoparticles. We compute the average of the hydrogen bonds of the last 0.1 ns after the simulations of 2 ns, with cutoff 30 Å and angle 30. The number of hydrogen bonds in the PCL75DL25 nanoparticle is 108 ± 6.48 while the number of hydrogen bonds in the PCL75LL25 nanoparticle is 67.909 ± 6.685. That is, there are about 50% more physical crosslinks in the PCL75DL25 nanoparticle. To reveal the molecular origin of the significant difference on the number of hydrogen bonds in PCL75DL25 and PCL75LL25 nanoparticles, we further analyze the hydrogen bonds of PCL75DL25 and PCL75LL25 by categorizing the type of hydrogen bonds into three groups. We named the segment with PCL as soft segments WDPUPCL, and the segment with PLA as soft segment WDPUPLA. Thus the three categories are: (1) hydrogen bonds between WDPUPCL and WDPUPCL; (2) hydrogen bonds between WDPUPLA and WDPUPLA; (3) hydrogen bonds between WDPUPCL and WDPUPLA. The results are shown in Fig. 5 (c).
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314
6
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
(a) 0.025 PCL75DL25
PCL75LL25
0.02 0.015 0.01 0.005
0
(b)
40
60
PCL75DL25
80
100 distance (Å)
120
140
(c)
PCL75LL25
160
180
PCL75DL25
0.001
0.0006
0.0008
0.0005
PCL75LL25
R O
0.0004
0.0006
0.0003 0.0004 0.0002 0.0002
0.0001
0 50
100
150
200
0
distance (Å)
P
0 0
200
F
20
O
0
50
100 distance (Å)
150
200
(b)
T
140
C
100
E
80 60 40
R R
Number of hbonds
120
DL
20 1.6
1.7
O
0 1.5
1.8
140 120
Number of hbonds
(a)
E
D
Fig. 4. Analysis of the radial distribution function between different polymer chains. (a) Radial distribution of the polymer chains. (b) Radial distribution of the PCL segments. (c) Radial distribution of the PLA segments.
LL
100 80 60 40 20 0
1.9
PCL75DL25
2
PCL75LL25
Time(ns)
80 60 40 20
DL
LL
0
20
Number of hbonds
Number of hbonds
U
100
Number of hbonds
N
C
(c)
15 10 5 DL
LL
Time(ns)
25 20 15 10 5
DL
LL
0
0
1.5 1.6 1.7 1.8 1.9 2
30
1.5 1.6 1.7 1.8 1.9 Time(ns)
2
1.5 1.6 1.7 1.8 1.9
2
Time(ns)
Fig. 5. Analysis of the number of hydrogen bonds in the PCL75DL25 and PCL75LL25 nanoparticles. The hydrogen bond is defined when the distance between the donor and acceptor is less than 3.0 Å and the angle is less than 30 degrees. (a) The number of total hydrogen bonds versus simulation time. (b) The average number of total hydrogen bonds. The error bars are the standard deviations from the ensemble average of the last 0.1 ns simulations. (c) Analysis of the number of three different kinds of hydrogen bonds.
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
336 337 338 339 340 341 342 343
F
O
R O
334 335
P
332 333
D
330 331
E
328 329
T
326 327
C
324 325
345 346
E
322 323
Full atomistic simulations of the PCL75DL25 and PCL75LL25 nanoparticles reveal that both nanoparticles do not have perfect spherical shapes. The PCL75LL25 nanoparticle is more irregular and asymmetry. These results are consistent with the experiment results that the size of PCL80LL20 nanoparticle is larger than the size of PCL80DL20 nanoparticle, and the structure of PCL80LL20 particle is relatively loose [16]. The analysis of the density distribution of the nanoparticle show that there is no difference on the center of the nanoparticles. Both nanoparticles have loose packing on the surface allowing them to interact with other nanoparticles in the solution to form hydrogels at high solid content and high temperature. We find that the PCL75LL25 has a larger surface layer and lower packing density in the surface layer. This may be the molecular origin of the faster sol-gel transition for the PCL75LL25. In order to understand the molecular structures of the WDPU nanoparticles, end to end distance of each polymer chains and hydrogen bonds are analyzed. We find that there are significant differences on the end to end distance and number of hydrogen bonds between the PCL75DL25 and PCL75LL25 nanoparticles. The polymer chains in the PCL75DL25 has larger end to end distances. The structure of WDPU with PDLLA as soft segments is more extended, while the distance between PLLA and PCL is close, which may reduce the occurrence of hydrogen bonding in PCL. There are more hydrogen bonds between PCL segments in the PCL75DL25 nanoparticle. Our results reveal that the optical rotation of PLA segments would significantly affect the molecular structures of PCL segments and the number of physical crosslinks between PCL segments. These results provide a possible molecular mechanism on how the WDPU nanoparticles form hydrogels and reveal the
R
321
344
O R
319 320
4. Discussion and Conclusions
N C
317 318
Surprisingly, our results reveal that the major difference does not occur between WDPUPLA and WDPUPLA. The major difference occurs between WDPUPCL instead, although both nanoparticles have the same WDPUPCL soft segments. These results suggest that the PDLLA soft segment increases the number of hydrogen bonds between the PCL soft segments in the PCL75DL25 nanoparticle. Our findings provide a molecular evidence that the PLLA soft segments induces the crystalline of PCL soft segments as observed in previous experimental studies [27]. We anticipate that the difference in the location and number of hydrogen bond affects the configurations of the WDPU nanoparticles. From previous analysis, we know that the major difference of hydrogen bond number is between WDPUPCL segments, suggesting that the PLA soft segment enhances the interactions between PCL segments in the nanoparticles. Further analysis of the simulation results shows that the conformation of PCL75DL25 and PCL75LL25 are different. Fig. 6 illustrates the interactions between polymer chains of PCL75DL25 and PCL75LL25, the color pink and red represent WDPUPLA, and WDPUPCL is in blue and cyan. For PCL75DL25, WDPUPLA have larger end to end distances and tend to tangle together resulting a separation of the PCL soft segments, while for PCL75LL25, PLA segments and PCL segments prefer to assemble together resulting shorter end to end distances. Fig. 6 also illustrates the intramolecular hydrogen bonds of PCL75DL25 and PCL75LL25 nanoparticles. These results show that the optical rotation of PLA would significantly affect the conformation of the whole WDPU and the number of physical crosslinks within the nanoparticle. PDLLA would result a configuration in which molecules are more extended, while the configuration of WDPU with PLLA as soft segment is more staggered.
U
315 316
7
Fig. 6. Illustration of the intermolecular and intramolecular hydrogen bond conformations in WDPU nanoparticles. (a) PCL75DL25 and (b) PCL75LL25.
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371
381
The authors declare no conflict of interest.
382
Acknowledgments
383 384 385 Q3
This research was supported by Ministry of Science and Technology (MOST 106-2221-E-002-136), Taiwan, R.O.C. and by National Taiwan University NTU-CC-107L891102 and NTU-CC-108L891102.
386
References
N
C
O
R R
E
C
T
[1] Thoniyot P, et al. Nanoparticle–hydrogel composites: concept, design, and applications of these promising, multi-functional materials. Adv Sci 2015;2:1–2. [2] Nair LS, Laurencin CT. Biodegradable polymers as biomaterials. Prog Polym Sci 2007; 32(8):762–98. [3] Ulery BD, Nair LS, Laurencin CT. Biomedical applications of biodegradable polymers. J Polym Sci B 2011;49(12):832–64. [4] Zohuriaan-Mehr MJ, Kabiri K. Superabsorbent polymer materials: a review. Iran Polym J 2008;17(6):451–77. [5] Hoffman AS. Hydrogels for biomedical applications. Adv Drug Deliv Rev 2002;54(1): 3–12. [6] Schmaljohann D. Thermo- and pH-responsive polymers in drug delivery. Adv Drug Deliv Rev 2006;58(15):1655–70. [7] Kohjiya S, et al. Drug release behavior from polyurethane gel. React Polym 1991;15: 165–75. [8] Park D, Wu W, Wang YD. A functionalizable reverse thermal gel based on a polyurethane/PEG block copolymer. Biomaterials 2011;32(3):777–86. [9] Ahmed EM, et al. An innovative method for preparation of nanometal hydroxide superabsorbent hydrogel. Carbohydr Polym 2013;91(2):693–8.
U
387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404
F
378
O
Conflicts of Interest
376 377
R O
380
374 375
[10] Ahmed EM. Hydrogel: preparation, characterization, and applications: a review. J Adv Res 2015;6(2):105–21. [11] Beebe DJ, et al. Functional hydrogel structures for autonomous flow control inside microfluidic channels. Nature 2000;404(6778):588–90. [12] Patil J, et al. Ionotropic gelation and polyelectrolyte complexation: the novel techniques to design hydrogel particulate sustained, modulated drug delivery system: a review. Dig J Nanomater Biostruct 2010;5(1):241–8. [13] Peppas NA. Biomedical applications of hydrogels handbook. Springer Science & Business Media; 2010. [14] Satarkar NS, Biswal D, Hilt JZ. Hydrogel nanocomposites: a review of applications as remote controlled biomaterials. Soft Matter 2010;6(11):2364–71. [15] Hsu S-h, et al. Water-based synthesis and processing of novel biodegradable elastomers for medical applications. J Mater Chem B 2014;2(31):5083–92. [16] Hsieh F-Y, Lin H-H, Hsu S-h. 3D bioprinting of neural stem cell-laden thermoresponsive biodegradable polyurethane hydrogel and potential in central nervous system repair. Biomaterials 2015;71:48–57. [17] Hsieh F-Y, et al. A novel biodegradable self-healing hydrogel to induce blood capillary formation. NPG Asia Mater 2017;9(3):e363. [18] Hsu S-h, Chang W-C, Yen C-T. Novel flexible nerve conduits made of water-based biodegradable polyurethane for peripheral nerve regeneration. J Biomed Mater Res A 2017;105(5):1383–92. [19] Ou C-W, et al. Characterization of biodegradable polyurethane nanoparticles and thermally induced self-assembly in water dispersion. ACS Appl Mater Interfaces 2014;6(8):5685–94. [20] Hung K-C, Tseng C-S, Hsu S-h. Synthesis and 3D printing of biodegradable polyurethane elastomer by a water-based process for cartilage tissue engineering applications. Adv Healthc Mater 2014;3(10):1578–87. [21] Hobza P, et al. Performance of empirical potentials (AMBER, CFF95, CVFF, CHARMM, OPLS, POLTEV), semiempirical quantum chemical methods (AM1, MNDO/M, PM3), and ab initio Hartree–Fock method for interaction of DNA bases: Comparison with nonempirical beyond Hartree–Fock results. J Comput Chem 1997;18(9):1136–50. [22] Dauber-Osguthorpe P, et al. Structure and energetics of ligand binding to proteins: Escherichia coli dihydrofolate reductase-trimethoprim, a drug-receptor system. Proteins 1988;4(1):31–47. [23] Mark JE. Physical properties of polymers handbook. , vol. 1076Springer; 2007. [24] Manual C. Force field based simulations. San Diego, CA: MSI; 1998. [25] Ma QT, et al. Guest dependent pressure behavior of the flexible MIL-53(Cr): a computational exploration. Dalton Trans 2012;41(14):3915–9. [26] Bruce CD, et al. Molecular dynamics simulation of sodium dodecyl sulfate micelle in water: micellar structural characteristics and counterion distribution. J Phys Chem B 2002;106(15):3788–93. [27] Lin HH, et al. Preparation and characterization of a biodegradable polyurethane hydrogel and the hybrid gel with soy protein for 3D cell-laden bioprinting. J Mater Chem B 2016;4(41):6694–705.
P
379
molecular origin of the faster sol-gel transition for the PCL75LL25. Our findings suggest that adding a small amount of PDLLA soft segment could significantly increase the number of hydrogen bonds within the nanoparticles and increases the sol-gel transition behavior. The fundamental insights into the molecular structure and mechanism of WDPU nanoparticles could help enabling the design of material properties of biocompatible hydrogel by precisely tuning the chemical compositions in the soft segments of the WDPU.
D
372 373
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
E
8
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446
Contents lists available at ScienceDirect
journal homepage: www.elsevier.com/locate/csbj
4 5 6
a
7
a r t i c l e
8 9 10 11 12 17 16 15 14 13 37 38 39 40 41 42 43
Article history: Received 29 September 2018 Received in revised form 20 December 2018 Accepted 25 December 2018 Available online xxxx
Department of Civil Engineering, National Taiwan University, Taipei, Taiwan Inistitute of Polymer Science and Engineering, National Taiwan University, Taipei, Taiwan
i n f o
a b s t r a c t
Biodegradable hydrogels have become promising materials for many biological applications in the past years. Recently, novel waterborne biodegradable polyurethane (WDPU) nanoparticles have been synthesized by a green water-based process, and serve as fundamental building blocks to form materials with great biocompatibility, biodegradability, and mechanical properties. However, the molecular structures and mechanisms of the WDPU nanoparticles and the relationship between the chemical compositions of the polymer segments and the material properties of the biodegradable hydrogels at macro-scale are still not well understood. In this study, we explore the fundamental mechanisms of WDPU nanoparticles through a full atomistic simulation approach to understand how the chemical compositions at the molecular level affect the molecular structures and material properties of WDPU nanoparticles. Specifically, we compare two WDPUs, i.e. PCL75LL25 and PCL75DL25, of the same hard segment composition and very similar soft segment composition [75% poly(e-caprolatone) and 25% polylactide], except the lactide in the former is L-form and in the latter is D,L-form. Our results show that the material properties of the biodegradable hydrogel can be designed by tuning the chemical compositions of the polymer segments. We find that the PCL75DL25 and PCL75LL25 have distinct molecular structures and physical crosslinks within the nanoparticles. The molecular structure of WDPU with PDLLA as soft segments is more extended, leading to more physical crosslinks between PCL segments. This study provide fundamental insights into the molecular structures and mechanisms of WDPU nanoparticles and help enabling the design of material properties of biocompatible hydrogel. © 2018 Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
D
b
F
Chien-Hui Wen a, Shun-Chieh Hsu a,b, Shan-hui Hsu b,⁎, Shu-Wei Chang a,⁎
O
Q23Q1
R O
2
Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles
P
1
E
C
T
E
Keywords: Molecular dynamics Polymer Hydrogel Nanoparticle biomaterial green material
R
36
52 53 54 55 56 57 58 59 60 61 62 63 64
In the past few decades, the development and applications of hydrogel have attracted increasing interest, specifically as a new class of biomaterials [1–8]. Hydrogels, namely hydrophilic gels, are networks of polymer chains that are extensively swollen with water [1,9–14]. Among hydrogels, biodegradable hydrogels are of utmost interest because they exhibit similar mechanical properties as natural soft tissues and can be degraded in an aqueous environment after their useful lifetime. Biodegradable hydrogels include natural hydrogels and synthetic type. Natural hydrogels have been gradually replaced by synthetic biodegradable hydrogels because of their strong water absorption capacity, long service life time, and a wide range of available chemical materials [10]. Recently, a novel waterborne biodegradable polyurethane (WDPU) has been synthesized and shown to have great potential in biomedical applications [15–20]. It is synthesized by a green water-based process, and has great biocompatibility, biodegradability, and mechanical
N C
50 51
1. Introduction
U
49
O R
47 45 44 46 48
18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35
⁎ Corresponding authors. E-mail addresses: [email protected] (S. Hsu), [email protected] (S.-W. Chang).
properties. The WDPU has been demonstrated that it can be synthesized using biodegradable polyesters as the soft segment without chemical crosslinking [18]. The properties of the WDPU nanoparticle could be adjusted by the chemical composition and ratio of the polyesters in the soft segment [15], such as PCL diol, ad PLA diol. The behavior of thermally induced self-assembly and gel formation and the degradation rate of the hydrogel can be tuned by changing the soft segments of WDPU, indicating that WDPU hydrogels with a various of material properties can be designed by changing the composition of soft segments for different biomedical applications. It has also been shown that the WDPU hydrogels are able to be applied to 3D printing for neutral stem cells (NSCs) carrier, and shown great potential in central nervous system repair [16]. The integration of biodegradable hydrogel and 3D printing technology has opened great opportunities for the design of smart biocompatible scaffolds for many applications due to the ability to access complex internal structures [15–20]. WDPU nanoparticles are the fundamental building blocks of the WDPU hydrogels. The molecular structures of WDPU nanoparticles play an important role in the molecular mechanisms of sol-gel transition and the material properties of the WDPU hydrogel. However, the molecular structures and mechanisms of the WDPU and the
https://doi.org/10.1016/j.csbj.2018.12.007 2001-0370/© 2018 Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85
94
2. Materials and Methods
95
2.1. Simulation Details
96
2.1.1. Model Construction Fig. 1(a) illustrates the polymer chains constructed in this work. Three different soft segments are used in this study: PCL diol, poly(D,Llactide) diol (PDLLA diol, Mn ~ 2000 Da) and poly(L-lactide) diol (PLLA diol, Mn ~ 2000). The WDPU with PCL diol as both soft segments is
P D E T C E R R
99 100
O
97 98
C
92
N
90 91
U
88 89
named WDPUPCL. The WDPU with PDLLA diol as both soft segments is named WDPUPDLA. The WDPU with PLLAdiol as both soft segments is named WDPUPLLA. Each WDPU polymer chain is consisted of 4 repeating units with two soft segments. PCL and PDLLA are mixed in 3:1 M ratio to construct the PCL75DL25 nanoparticle while PCL and PLLA are mixed in 3:1 M ratio to construct the PCL75LL25 nanoparticle (Fig. 1(a)). We use Materials Studio to construct WDPU monomers according to the previously published synthesis procedure [15]. Each chain contains 3070 atoms and the molecular weight is about 39,018 Da. We use Amorphous Cell to construct the initial amorphous cells for PCL75DL25 and PCL75LL25 as shown in Fig. 1(b). Each initial unit cell consists 10 chains of polymers. To equilibrate the amorphous structure, we firstly run NPT simulation at 500 K for 0.6 ns to speed up the process, then cool down the temperature to 300 K for 0.2 ns to reach equilibrium at the room temperature. The equilibrated unit cell is then replicated in all three directions to construct a larger amorphous packing nanoparticle (Fig. 1(c))
R O
93
relationship between the chemical compositions of the polymer segments and the material properties of the biodegradable hydrogels at macro-scale are still not clear. In this paper, we employed full atomistic simulations to explore the molecular structure and fundamental mechanisms of WDPU nanoparticles. WDPU nanoparticles with 75% PCL diol and 25% PLA diol soft segments (L-form or D,L-form) are simulated to investigate the relationships between the chemical composition of soft segments and the physical properties of the WDPU nanoparticles.
F
86 87
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
O
2
Fig. 1. (a) Chemical composition of the polymer chains. (b) WDPU polymer chains with 75% PCL and 25% PLA are constructed. (c) Illustration of the amorphous unit cells.
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143
h
i
1−e−αðb−b0 Þ þ
2.2.3. Surface Area Solvent-accessible surface area of atoms in the nanoparticle is calculated by using the assigned radius for each atom with an extended radius for 1.4 Å to find the points on a sphere that are exposed to solvent. To normalize the surface area, we calculate the surface area of WDPU nanoparticles (Am), and divided by surface area when the nanoparticle is spherical (As). Values close to one imply a higher spherical symmetry while a larger value indicates that the nanoparticle has unsmooth surface or irregular shape. We compute the average of the surface area of the last 0.1 ns after the simulations of 2 ns. 2.2.4. End to End Distance End to end distance is a parameter to describe the conformation of polymer chains. To calculate the end to end distance of polymer chains, we define the atomic positions of the Nitrogen atoms on N-C-C-N as the ends of a segment, and the distance between these two atomic positions are considered as the end to end distance of segments. 3. Results
X
H θ ðθ−θ0 Þ2 X XX þ H ϕ ½1 þ s cosðnϕÞ þ Hχ χ 2 þ F bb0 ðb−b0 Þ χ b X b0 X X X ϕ0 0 F θθ0 ðθ−θ0 Þ θ0 −θ00 þ F bθ ðb−b0 Þ b −b0 þ 0 θ bXθX θ X 0 ðθ−θ0 Þ þ F ϕθθ0 cosϕðθ−θ0 Þ θ −θ00 þ F χχ0 χχ 0 Db bX
166
θ
151 152
2.2.1. Radius of Gyration Radius of gyration, Rg, is a useful representation to measure and characterize the size of complex polymer configurations. We compute the average of the radius of gyration of the last 0.1 ns after the simulations of 2 ns. We use Eq. (2) to calculate the radius of gyration of ! WDPU nanoparticles. In Eq. (2), r i is the position of the i-th atom and ! r cm is the center of mass of the WDPU nanoparticle.
N C
149 150
ð1Þ
O R
2.2. Analysis of Molecular Structures
147 148
χ0
R
145
146
χ
" ϕ 6 # X qi q j r 12 r þ ε −2 r r εr ij
E
þ
X
C
T
Epot ¼
X
U
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi u N u 1 X ! 2 r i −! r cm Rg ¼ t N i¼1
ð2Þ
154
155 156 157 158 159 160 161 162 163
ð3Þ
F
124
Imin Iaverage
2.2.2. Eccentricity To understand if the differences on the radius of gyration are resulted from the shape of nanoparticles, we calculate the eccentricity (ρeccen) of the nanoparticles to investigate the structures of nanoparticles. The eccentricity measures the symmetry of WDPU nanoparticles. If the value is close to zero, it implies a high symmetry, while if the value is larger, it implies less symmetry. We use Eq. (3) to calculate the eccentricity of WDPU nanoparticles [26]. In Eq. (3), Imin is the minimum value of the tri-axial moment of inertia and Iavg is the average value of the tri-axial moment of inertia. We compute the average of
167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182
3.1. Nanoparticle Characterization
183
To characterize the structures of WDPU nanoparticles, we analyze the radius of gyrations, eccentricity, and surface area of PCL75DL25 and PCL75LL25 nanoparticles. Fig. 2 (a) shows the results. We find that the PCL75LL25 nanoparticle has a radius of gyration with a value of 70.01 ± 0.03 Å, while the PCL75DL25 nanoparticle has a slightly smaller radius of gyration of 68.90 ± 0.03 Å. The results of the eccentricity show that the PCL75DL25 nanoparticle has a value of ρeccen = 0.039 ± 0.004 and the value of the PCL75LL25 nanoparticle is 0.079 ± 0.004, suggesting that the shape of the PCL75LL25 nanoparticle is more irregular. For the surface area, there is no significant differences on the surface area between the PCL75DL25 and the PCL75LL25 nanoparticles. It is worth noting that these results show that both PCL75DL25 and PCL75LL25 have unsmoothed surface and the PCL75LL25 nanoparticle has a more irregular and unsymmetrical shape. Since both polymers have the same molecular weights, the differences on the radius of gyration might be resulted from either the change on the packing of molecules within the nanoparticles or the change on the surface morphology. We analyze the density distributions of both nanoparticles as shown in Fig. 2 (b). In Fig. 2 (b), we find that both nanoparticles have the same core density when the radius is smaller than 50 Å. These results indicate that there is no difference on the density in the center between different WDPU nanoparticles, while there is significant differences on the surface morphology between different WDPU nanoparticles. From the distance of 55 Å, the density of both PCL75DL25 and PCL75LL25 nanoparticles start to decrease indicating that the packing of polymer chains on the surfaces of both nanoparticles are much looser than the molecular packing in the center of the nanoparticles. Although both nanoparticles have a radius of gyration ~70 Å, the density of PCL75DL25 and PCL75LL25 decrease slowly and goes to 0 when the radius exceeds 110 Å. It is worth noting that from 50 to 80 Å, the PCL75DL25 has a larger density while a reverse trend is found for the distance from 80 to 120 Å. The density of PCL75DL25 reaches 0 at a distance smaller than the PCL75LL25. These results suggest that the selfassembly of PCL75LL25 nanoparticle is more loosely on the surface when compare with the PCL75DL25 nanoparticle. That is, the slightly larger radius gyration of the PCL75LL25 is not resulted from different packing density in the center of the nanoparticle but resulted from the differences on the surface morphology of the WDPU nanoparticles. We
184 185
E
122 123
P eccen ¼ 1−
O
2.1.2. Force Field We use consistent valence force field (CVFF) [21] to simulate the molecular structure of the WDPU nanoparticles. The CVFF originally applied to biological system [22] and is a prototype of the consistent force field (CFF) and its later derivatives (the polymer consistent force field PCFF and COMPASS). Terms in CVFF included a bond stretching or compression, an angle bending, torsion angle twisting and out-of-plane deformation of a planar system, five cross-coupling terms between bond deformations, angle bending, a bond deformation and an angle bending, a torsion angle twisting and the two associated angle bending, and coupling between out-of-plane deformations, and repulsive, dispersive, and Coulombic interactions between nonbonded atoms [23]. The energy description is shown in Eq. (1). The first four terms in Eq. (1) included a Morse potential for bond stretching and compression (b), a harmonic term for angle bending (θ), cosine torsional (ϕ) and out-ofplane torsional terms (χ). Db, Hθ, Hϕ, and Hχ are the force constants for the corresponding intramolecular deformations. b0 and θ0 are the equilibrium bond length and angle. The next five terms are cross terms that account for interactions between the four types of internal coordinates. Fij are the force constants for the cross terms. The last two terms represents Lennard-Jones and Coulombic terms for nonbonded interaction. ε and r∗ are the parameters for the nonbonded repulsive and dispersive interactions, and qi are the partial charges on each atom [23–25].
R O
121
164
the eccentricity of the last 0.1 ns after the simulations of 2 ns.
P
to simulate the self-assembly of the nanoparticle in a larger box in vacuum. The final model consists of 245,600 atoms.
D
119 120
3
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221
4
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
(a) Eccentricity
70.5
0.1
70
0.08
69.5 69 68.5
Surface area 2.54 2.52 2.5
0.06
A m/As
Eccentricity
0.04
2.42
0
2.4
PCL75DL25 PCL75LL25
PCL75DL25
PCL75DL25
PCL75LL25
R O
2 1.8 1.6
1.4 1.2
P
Density
PCL75LL25
F
PCL75DL25 PCL75LL25
(b)
2.46 2.44
0.02
68
2.48
O
Radius of gyration (Å)
Radius of gyration
1
0.8
D
0.6
10
30
50
70
T
0.2 0 -10
E
0.4
90
110
130
150
C
distance from the center of particle d (Å)
222 223
R R
E
Fig. 2. Nanoparticle characterizations. (a) Radius of gyration, eccentricity and surface area of WDPN nanoparticles. The error bars are the standard deviations from the ensemble average of the last 0.1 ns simulations. (b) Analysis of the density distribution of WDPU nanoparticles. The results show that the PCL75LL25 nanoparticle has a more irregular shape. There is no differences on the core density between the PCL75DL25 and PCL75LL25 nanoparticles but there is significant differences on the surface morphology when the PLLA soft segment is replaced with PDLLA soft segments. The molecular packing on the surface of nanoparticle is less dense for the PCL75DL25 nanoparticle.
231
3.2. Molecular Structure of WDPU Nanoparticles
232
The characterizations of WDPU nanoparticles show that the PCL75LL25 nanoparticle has smaller density in the surface layer and higher value of eccentricity. The packing and self-assembly of the nanoparticles are determined by the differences on the properties of each chain and the interactions between polymer chains. To understand the molecular structural differences between the PCL75DL25 and PCL75LL25 nanoparticles, we further analyze the end to end distances of polymer chains as shown in Fig. 3 (a). We find that although there is not much difference on the characterization of the shape between the PCL75DL25 and PCL75LL25 nanoparticles, there are significant differences at the molecular level between the PCL75DL25 and PCL75LL25 nanoparticles. The end to end distances of the polymer chains in the PCL75DL25 nanoparticles (99.976 ± 0.134 Å) is almost two times larger than the end to end distance of the polymer chains in
233 234 235 236 237 238 239 240 241 242 243 244 245
C
228 229
N
226 227
U
224 225
O
230
anticipate that the loose surface packing on both nanoparticles play important roles in the sol-gel transition. The loose surface packing allows polymer chains on the surface to move more freely and interact with other nanoparticles in solutions at higher temperature to initiate the physical crosslinks between nanoparticles and eventually form networks and further initiate the hydrogel formation. The larger radius of gyration and the looser packing on the surface of the PCL75LL25 nanoparticle provide molecular insights into its faster sol-gel transition time observed in experimental results.
the PCL75LL25 nanoparticles (54.990 ± 0.228 Å). It is worth noting that all the polymer chains in both WDPU nanoparticles have similar contour length, indicating that the polymer chains in the PCL75DL25 nanoparticle is much stiffer than the polymer chains in the PCL75LL25 nanoparticle. These results suggest that PLLA soft segments significantly cause shorter end to end distances of polymer chains while the PDLLA soft segments do not have the same effects. Since all the polymer chains in both nanoparticles have 75% PCL segments and 25% PLA segments (PDLLA soft segments in PCL75DL25 and PLLA soft segments in PCL75LL25), we further analyze the segment to segment distances of PCL segments and PLA segments in both nanoparticles as shown in Fig. 3 (b). The distributions of the segment to segment distances are shown in Fig. 3 (c) and (d). The results reveal that there are significant differences on the segment to segment distances of the PCL soft segments. The PLLA soft segments significantly cause shorter segment to segment distances of PCL soft segments in the PCL75LL25 nanoparticle while the PCL soft segments have much larger segment to segment distances in the PCL75DL25 nanoparticle. These results suggest that PDLLA increases the stiffness of the PCL soft segments in the PCL75DL25 nanoparticle. The shorter end to end distances, i.e. lower stiffness, for the PCL75LL25 nanoparticle suggest that when solvating in the solution, the polymer chains are easier to reach out and form physical crosslinks with other nanoparticles. We further analyze the radial distribution function (RDF) between different polymer chains to study the packing behaviors of the PCL75DL25 and PCL75LL25 nanoparticles. Fig. 4 (a) shows the RDF
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
(b) 90
120
Segment to Segment distance (Å)
100 80 60 40 20
80 70 50 40 30 20 10 0
0 PCL75DL25
(c)
PCL
PCL75LL25
O
DL LL
30
R O
25 Numbers
PLA
PLA Segment to Segment distance distribution (dPLA) 35
20 15 10
P
5 0 0
10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170
E
DL LL
D
PCL Segment to Segment distance distribution (dPCL)
(d) 25
T
20 15 10
C
Numbers
DL LL
60
F
End to end distance (Å)
(a)
5
E
5 0
10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170
R
0
276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293
N C
274 275
between the polymer chains. The peak of the RDF occurs at 17.75 Å for the PCL75DL25 nanoparticle while at 16.75 Å for the PCL75LL25 nanoparticle, suggesting the packing of both nanoparticles are similar from the point of view of the polymer chains. It is worth noting that the RDF of the PCL75DL25 nanoparticle drops to zero at a shorter distance when compared with the RDF of the PCL75LL25 nanoparticle. These results are consistent with the analysis of the density distribution of the nanoparticles as shown in Fig. 2 (b). Although there is no significant difference on the RDF of the polymer chains, we find that there is significant difference if we analyze the RDF of the PCL and PLA soft segments separately for the nanoparticles. Fig. 4 (b) and (c) show the RDF of the PCL soft segments and PLA soft segments. The RDF of the PCL soft segments are similar for both nanoparticles while there is a significant difference on the RDF of the PLA soft segments as shown in Fig. 4 (c). The first peak of the RDF of the PLA soft segments occurs at 8.25 Å for the PCL75DL25 nanoparticle while at 23.75 Å for the PCL75LL25 nanoparticle. Together with the results from the analysis of the end to end distances and segment to segment distances as shown in Fig. 3, our results reveal that the PCL and PLA soft segments are more extended in the PCL75DL25 nanoparticle resulted in a packing of polymer chains with shorter distances between PLA soft segments in the PCL75DL25 nanoparticle.
U
272 273
O R
Fig. 3. Analysis of the end to end distance of polymer chains in WDPU nanoparticles. (a) The end to end distances of polymer chains in WDPU nanoparticles. (b) PCL and PLA segment to segment distances in WDPU nanoparticles. The error bars are the standard deviations from the ensemble average of the last 0.1 ns simulations. (c) the distribution of the end to end distances of PLA segments in WDPU nanoparticles. (d) the distribution of the end to end distances of PCL segments in WDPU nanoparticles.
3.3. Intermolecular and Intramolecular Hydrogen Bonds in WDPU 294 Nanoparticles 295 The analysis of hydrogen bonds in the WDPU nanoparticles are shown in Fig. 5 to investigate the differences on the molecular interactions between the PCL75DL25 and PCL75LL25 nanoparticles. We compute the average of the hydrogen bonds of the last 0.1 ns after the simulations of 2 ns, with cutoff 30 Å and angle 30. The number of hydrogen bonds in the PCL75DL25 nanoparticle is 108 ± 6.48 while the number of hydrogen bonds in the PCL75LL25 nanoparticle is 67.909 ± 6.685. That is, there are about 50% more physical crosslinks in the PCL75DL25 nanoparticle. To reveal the molecular origin of the significant difference on the number of hydrogen bonds in PCL75DL25 and PCL75LL25 nanoparticles, we further analyze the hydrogen bonds of PCL75DL25 and PCL75LL25 by categorizing the type of hydrogen bonds into three groups. We named the segment with PCL as soft segments WDPUPCL, and the segment with PLA as soft segment WDPUPLA. Thus the three categories are: (1) hydrogen bonds between WDPUPCL and WDPUPCL; (2) hydrogen bonds between WDPUPLA and WDPUPLA; (3) hydrogen bonds between WDPUPCL and WDPUPLA. The results are shown in Fig. 5 (c).
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314
6
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
(a) 0.025 PCL75DL25
PCL75LL25
0.02 0.015 0.01 0.005
0
(b)
40
60
PCL75DL25
80
100 distance (Å)
120
140
(c)
PCL75LL25
160
180
PCL75DL25
0.001
0.0006
0.0008
0.0005
PCL75LL25
R O
0.0004
0.0006
0.0003 0.0004 0.0002 0.0002
0.0001
0 50
100
150
200
0
distance (Å)
P
0 0
200
F
20
O
0
50
100 distance (Å)
150
200
(b)
T
140
C
100
E
80 60 40
R R
Number of hbonds
120
DL
20 1.6
1.7
O
0 1.5
1.8
140 120
Number of hbonds
(a)
E
D
Fig. 4. Analysis of the radial distribution function between different polymer chains. (a) Radial distribution of the polymer chains. (b) Radial distribution of the PCL segments. (c) Radial distribution of the PLA segments.
LL
100 80 60 40 20 0
1.9
PCL75DL25
2
PCL75LL25
Time(ns)
80 60 40 20
DL
LL
0
20
Number of hbonds
Number of hbonds
U
100
Number of hbonds
N
C
(c)
15 10 5 DL
LL
Time(ns)
25 20 15 10 5
DL
LL
0
0
1.5 1.6 1.7 1.8 1.9 2
30
1.5 1.6 1.7 1.8 1.9 Time(ns)
2
1.5 1.6 1.7 1.8 1.9
2
Time(ns)
Fig. 5. Analysis of the number of hydrogen bonds in the PCL75DL25 and PCL75LL25 nanoparticles. The hydrogen bond is defined when the distance between the donor and acceptor is less than 3.0 Å and the angle is less than 30 degrees. (a) The number of total hydrogen bonds versus simulation time. (b) The average number of total hydrogen bonds. The error bars are the standard deviations from the ensemble average of the last 0.1 ns simulations. (c) Analysis of the number of three different kinds of hydrogen bonds.
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
336 337 338 339 340 341 342 343
F
O
R O
334 335
P
332 333
D
330 331
E
328 329
T
326 327
C
324 325
345 346
E
322 323
Full atomistic simulations of the PCL75DL25 and PCL75LL25 nanoparticles reveal that both nanoparticles do not have perfect spherical shapes. The PCL75LL25 nanoparticle is more irregular and asymmetry. These results are consistent with the experiment results that the size of PCL80LL20 nanoparticle is larger than the size of PCL80DL20 nanoparticle, and the structure of PCL80LL20 particle is relatively loose [16]. The analysis of the density distribution of the nanoparticle show that there is no difference on the center of the nanoparticles. Both nanoparticles have loose packing on the surface allowing them to interact with other nanoparticles in the solution to form hydrogels at high solid content and high temperature. We find that the PCL75LL25 has a larger surface layer and lower packing density in the surface layer. This may be the molecular origin of the faster sol-gel transition for the PCL75LL25. In order to understand the molecular structures of the WDPU nanoparticles, end to end distance of each polymer chains and hydrogen bonds are analyzed. We find that there are significant differences on the end to end distance and number of hydrogen bonds between the PCL75DL25 and PCL75LL25 nanoparticles. The polymer chains in the PCL75DL25 has larger end to end distances. The structure of WDPU with PDLLA as soft segments is more extended, while the distance between PLLA and PCL is close, which may reduce the occurrence of hydrogen bonding in PCL. There are more hydrogen bonds between PCL segments in the PCL75DL25 nanoparticle. Our results reveal that the optical rotation of PLA segments would significantly affect the molecular structures of PCL segments and the number of physical crosslinks between PCL segments. These results provide a possible molecular mechanism on how the WDPU nanoparticles form hydrogels and reveal the
R
321
344
O R
319 320
4. Discussion and Conclusions
N C
317 318
Surprisingly, our results reveal that the major difference does not occur between WDPUPLA and WDPUPLA. The major difference occurs between WDPUPCL instead, although both nanoparticles have the same WDPUPCL soft segments. These results suggest that the PDLLA soft segment increases the number of hydrogen bonds between the PCL soft segments in the PCL75DL25 nanoparticle. Our findings provide a molecular evidence that the PLLA soft segments induces the crystalline of PCL soft segments as observed in previous experimental studies [27]. We anticipate that the difference in the location and number of hydrogen bond affects the configurations of the WDPU nanoparticles. From previous analysis, we know that the major difference of hydrogen bond number is between WDPUPCL segments, suggesting that the PLA soft segment enhances the interactions between PCL segments in the nanoparticles. Further analysis of the simulation results shows that the conformation of PCL75DL25 and PCL75LL25 are different. Fig. 6 illustrates the interactions between polymer chains of PCL75DL25 and PCL75LL25, the color pink and red represent WDPUPLA, and WDPUPCL is in blue and cyan. For PCL75DL25, WDPUPLA have larger end to end distances and tend to tangle together resulting a separation of the PCL soft segments, while for PCL75LL25, PLA segments and PCL segments prefer to assemble together resulting shorter end to end distances. Fig. 6 also illustrates the intramolecular hydrogen bonds of PCL75DL25 and PCL75LL25 nanoparticles. These results show that the optical rotation of PLA would significantly affect the conformation of the whole WDPU and the number of physical crosslinks within the nanoparticle. PDLLA would result a configuration in which molecules are more extended, while the configuration of WDPU with PLLA as soft segment is more staggered.
U
315 316
7
Fig. 6. Illustration of the intermolecular and intramolecular hydrogen bond conformations in WDPU nanoparticles. (a) PCL75DL25 and (b) PCL75LL25.
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371
381
The authors declare no conflict of interest.
382
Acknowledgments
383 384 385 Q3
This research was supported by Ministry of Science and Technology (MOST 106-2221-E-002-136), Taiwan, R.O.C. and by National Taiwan University NTU-CC-107L891102 and NTU-CC-108L891102.
386
References
N
C
O
R R
E
C
T
[1] Thoniyot P, et al. Nanoparticle–hydrogel composites: concept, design, and applications of these promising, multi-functional materials. Adv Sci 2015;2:1–2. [2] Nair LS, Laurencin CT. Biodegradable polymers as biomaterials. Prog Polym Sci 2007; 32(8):762–98. [3] Ulery BD, Nair LS, Laurencin CT. Biomedical applications of biodegradable polymers. J Polym Sci B 2011;49(12):832–64. [4] Zohuriaan-Mehr MJ, Kabiri K. Superabsorbent polymer materials: a review. Iran Polym J 2008;17(6):451–77. [5] Hoffman AS. Hydrogels for biomedical applications. Adv Drug Deliv Rev 2002;54(1): 3–12. [6] Schmaljohann D. Thermo- and pH-responsive polymers in drug delivery. Adv Drug Deliv Rev 2006;58(15):1655–70. [7] Kohjiya S, et al. Drug release behavior from polyurethane gel. React Polym 1991;15: 165–75. [8] Park D, Wu W, Wang YD. A functionalizable reverse thermal gel based on a polyurethane/PEG block copolymer. Biomaterials 2011;32(3):777–86. [9] Ahmed EM, et al. An innovative method for preparation of nanometal hydroxide superabsorbent hydrogel. Carbohydr Polym 2013;91(2):693–8.
U
387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404
F
378
O
Conflicts of Interest
376 377
R O
380
374 375
[10] Ahmed EM. Hydrogel: preparation, characterization, and applications: a review. J Adv Res 2015;6(2):105–21. [11] Beebe DJ, et al. Functional hydrogel structures for autonomous flow control inside microfluidic channels. Nature 2000;404(6778):588–90. [12] Patil J, et al. Ionotropic gelation and polyelectrolyte complexation: the novel techniques to design hydrogel particulate sustained, modulated drug delivery system: a review. Dig J Nanomater Biostruct 2010;5(1):241–8. [13] Peppas NA. Biomedical applications of hydrogels handbook. Springer Science & Business Media; 2010. [14] Satarkar NS, Biswal D, Hilt JZ. Hydrogel nanocomposites: a review of applications as remote controlled biomaterials. Soft Matter 2010;6(11):2364–71. [15] Hsu S-h, et al. Water-based synthesis and processing of novel biodegradable elastomers for medical applications. J Mater Chem B 2014;2(31):5083–92. [16] Hsieh F-Y, Lin H-H, Hsu S-h. 3D bioprinting of neural stem cell-laden thermoresponsive biodegradable polyurethane hydrogel and potential in central nervous system repair. Biomaterials 2015;71:48–57. [17] Hsieh F-Y, et al. A novel biodegradable self-healing hydrogel to induce blood capillary formation. NPG Asia Mater 2017;9(3):e363. [18] Hsu S-h, Chang W-C, Yen C-T. Novel flexible nerve conduits made of water-based biodegradable polyurethane for peripheral nerve regeneration. J Biomed Mater Res A 2017;105(5):1383–92. [19] Ou C-W, et al. Characterization of biodegradable polyurethane nanoparticles and thermally induced self-assembly in water dispersion. ACS Appl Mater Interfaces 2014;6(8):5685–94. [20] Hung K-C, Tseng C-S, Hsu S-h. Synthesis and 3D printing of biodegradable polyurethane elastomer by a water-based process for cartilage tissue engineering applications. Adv Healthc Mater 2014;3(10):1578–87. [21] Hobza P, et al. Performance of empirical potentials (AMBER, CFF95, CVFF, CHARMM, OPLS, POLTEV), semiempirical quantum chemical methods (AM1, MNDO/M, PM3), and ab initio Hartree–Fock method for interaction of DNA bases: Comparison with nonempirical beyond Hartree–Fock results. J Comput Chem 1997;18(9):1136–50. [22] Dauber-Osguthorpe P, et al. Structure and energetics of ligand binding to proteins: Escherichia coli dihydrofolate reductase-trimethoprim, a drug-receptor system. Proteins 1988;4(1):31–47. [23] Mark JE. Physical properties of polymers handbook. , vol. 1076Springer; 2007. [24] Manual C. Force field based simulations. San Diego, CA: MSI; 1998. [25] Ma QT, et al. Guest dependent pressure behavior of the flexible MIL-53(Cr): a computational exploration. Dalton Trans 2012;41(14):3915–9. [26] Bruce CD, et al. Molecular dynamics simulation of sodium dodecyl sulfate micelle in water: micellar structural characteristics and counterion distribution. J Phys Chem B 2002;106(15):3788–93. [27] Lin HH, et al. Preparation and characterization of a biodegradable polyurethane hydrogel and the hybrid gel with soy protein for 3D cell-laden bioprinting. J Mater Chem B 2016;4(41):6694–705.
P
379
molecular origin of the faster sol-gel transition for the PCL75LL25. Our findings suggest that adding a small amount of PDLLA soft segment could significantly increase the number of hydrogen bonds within the nanoparticles and increases the sol-gel transition behavior. The fundamental insights into the molecular structure and mechanism of WDPU nanoparticles could help enabling the design of material properties of biocompatible hydrogel by precisely tuning the chemical compositions in the soft segments of the WDPU.
D
372 373
C.-H. Wen et al. / Computational and Structural Biotechnology Journal xxx (xxxx) xxx
E
8
Please cite this article as: C.-H. Wen, S.-C. Hsu, S. Hsu, et al., Molecular Structures and Mechanisms of Waterborne Biodegradable Polyurethane Nanoparticles, Comput Struct Biotechnol J, https://doi.org/10.1016/j.csbj.2018.12.007
405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446