Silk fibroin hydrogel not only has biocompatibility, but also has environmental response ability. It plays an important role in the development of material. The gelation mechanism of silk fibroin hydrogel is very important to textile and medicine fields. The molecular dynamics simulation was used to discuss the structure and non-bond interaction of silk fibroin hydrogel. The results show that the non-bond interactions between silk fibroin molecules and water molecules have certain influence on the formation of silk fibroin hydrogel. According to the hydrogen bond analysis, the hydrogen bonds are mainly formed between random coil peptide fragments at the two ends of silk fibroin molecules and residues 252-254 are the key residues. The electrostatic and polar solvation interactions between silk fibroin molecules plays a major role in cross-linking of the coil segments of two silk fibroin molecules.
Citation: Xuewei Jiang, Lu Zheng, Huhe Wu, Jun Zhang. Investigations into the role of non-bond interaction on gelation mechanism of silk fibroin hydrogel[J]. Mathematical Biosciences and Engineering, 2021, 18(4): 4071-4083. doi: 10.3934/mbe.2021204
| [1] | Li Li, Xueqing Wang, Xiaoting Liu, Rui Guo, Ruidong Zhang . Integrative analysis reveals key mRNAs and lncRNAs in monocytes of osteoporotic patients. Mathematical Biosciences and Engineering, 2019, 16(5): 5947-5971. doi: 10.3934/mbe.2019298 |
| [2] | Jiu-Xin Tan, Shi-Hao Li, Zi-Mei Zhang, Cui-Xia Chen, Wei Chen, Hua Tang, Hao Lin . Identification of hormone binding proteins based on machine learning methods. Mathematical Biosciences and Engineering, 2019, 16(4): 2466-2480. doi: 10.3934/mbe.2019123 |
| [3] | Jinmiao Song, Shengwei Tian, Long Yu, Qimeng Yang, Qiguo Dai, Yuanxu Wang, Weidong Wu, Xiaodong Duan . RLF-LPI: An ensemble learning framework using sequence information for predicting lncRNA-protein interaction based on AE-ResLSTM and fuzzy decision. Mathematical Biosciences and Engineering, 2022, 19(5): 4749-4764. doi: 10.3934/mbe.2022222 |
| [4] | Hong Yuan, Jing Huang, Jin Li . Protein-ligand binding affinity prediction model based on graph attention network. Mathematical Biosciences and Engineering, 2021, 18(6): 9148-9162. doi: 10.3934/mbe.2021451 |
| [5] | Cristina De Ambrosi, Annalisa Barla, Lorenzo Tortolina, Nicoletta Castagnino, Raffaele Pesenti, Alessandro Verri, Alberto Ballestrero, Franco Patrone, Silvio Parodi . Parameter space exploration within dynamic simulations of signaling networks. Mathematical Biosciences and Engineering, 2013, 10(1): 103-120. doi: 10.3934/mbe.2013.10.103 |
| [6] | Xiao Tu, Qinran Zhang, Wei Zhang, Xiufen Zou . Single-cell data-driven mathematical model reveals possible molecular mechanisms of embryonic stem-cell differentiation. Mathematical Biosciences and Engineering, 2019, 16(5): 5877-5896. doi: 10.3934/mbe.2019294 |
| [7] | Bochuan Du, Pu Tian . Factorization in molecular modeling and belief propagation algorithms. Mathematical Biosciences and Engineering, 2023, 20(12): 21147-21162. doi: 10.3934/mbe.2023935 |
| [8] | Abdulhamed Alsisi, Raluca Eftimie, Dumitru Trucu . Non-local multiscale approach for the impact of go or grow hypothesis on tumour-viruses interactions. Mathematical Biosciences and Engineering, 2021, 18(5): 5252-5284. doi: 10.3934/mbe.2021267 |
| [9] | Linard Hoessly, Carsten Wiuf . Fast reactions with non-interacting species in stochastic reaction networks. Mathematical Biosciences and Engineering, 2022, 19(3): 2720-2749. doi: 10.3934/mbe.2022124 |
| [10] | O. E. Adebayo, S. Urcun, G. Rolin, S. P. A. Bordas, D. Trucu, R. Eftimie . Mathematical investigation of normal and abnormal wound healing dynamics: local and non-local models. Mathematical Biosciences and Engineering, 2023, 20(9): 17446-17498. doi: 10.3934/mbe.2023776 |
Silk fibroin hydrogel not only has biocompatibility, but also has environmental response ability. It plays an important role in the development of material. The gelation mechanism of silk fibroin hydrogel is very important to textile and medicine fields. The molecular dynamics simulation was used to discuss the structure and non-bond interaction of silk fibroin hydrogel. The results show that the non-bond interactions between silk fibroin molecules and water molecules have certain influence on the formation of silk fibroin hydrogel. According to the hydrogen bond analysis, the hydrogen bonds are mainly formed between random coil peptide fragments at the two ends of silk fibroin molecules and residues 252-254 are the key residues. The electrostatic and polar solvation interactions between silk fibroin molecules plays a major role in cross-linking of the coil segments of two silk fibroin molecules.
Silk is a kind of natural protein fiber with high strength. There are two silk fibroins and sericin around them [1,2]. Silk fibroin has excellent biocompatibility and degradability. It has been widely used in textile and biomaterial fields [3,4,5]. It has great potential in developing sustainable and environment-interactive materials [6].
Silk fibroin can be easily processed into various materials, such as nanoparticles and microcapsules [1,7,8,9,10,11]. Hydrogel is composed of polymer with three-dimensional network structure and medium filled in the gaps of polymer segments. The medium can be gas or liquid, and most of them are liquid in general. Therefore, hydrogels can be regarded as swelling bodies with three-dimensional polymer networks containing liquid. Almost all-natural gels and most synthetic gels are hydrogels with water as liquid medium. The silk fibroin hydrogel is a swelling body formed by water contained in the three-dimensional network. At present, silk fibroin hydrogel has been developed to construct functional tissues and has been used to solve many problems in drug release [12,13,14,15,16,17].
Silk fibroin hydrogel not only has biocompatibility, but also has environmental response ability, which plays an important role in the development of smart material. However, the gel time limit its application. It takes several days or even several weeks for the aqueous solution of silk fibroin with the concentration of 0.6-15% to transform into gel at room temperature [1,18]. Therefore, it is necessary to explore the gelation mechanism and promotes the formation of silk fibroin hydrogel.
Due to extensive hydrogen bonding, silk fibroin is a hydrophobic protein, 70% of the total components are alanine and glycine, with high crystallinity [19,20,21]. The silk fibroin molecule consists of light chain, heavy chain and glycoprotein. The ratio of them is 6:6:1 [22,23,24]. It is difficult to detect and quantify the interaction between these molecules by experimental methods. Molecular dynamics can provide atomic-scale details of atomic interaction and structure property, supplementing the details of the relation between them.
A four-stage physical gelation process was proposed by Liu et al [25]. Silk fibroin hydrogels have a multi-domain network with weak domain-domain interactions. By studying the process of silk fibroin gelation, they believed that silk fibroin gelation was divided into four steps. First, silk fibroin molecules aggregate into a nucleus. Then the crystal nucleus grows into β-crystallites. Subsequently, the nanofibers further grow and branch into a single-domain. Finally, these silk fiber networks penetrate each other and form a stable multi-domain. This process shows that the weak interaction plays a very important role in the formation of silk fibroin hydrogel. Therefore, an atomistic silk fibroin hydrogel model comprising of two anti-parallel β-sheet has been built in this work. The radial distribution, secondary structure, hydrogen bond and non-bond interaction have been used to describe the weak interaction and structure of the silk fibroin hydrogel based on the molecular dynamics [26].
The initial model was taken from the protein database (ID: 3UA0, Figure 1a). The molecular model comprises of two crystal structure of the N-terminal domain of silk fibroin with double anti-parallel β-sheet [2]. The missing part of this model was supplemented by the homologous modeling method with Chimera (Figure 1b) [27,28]. Two silk fibroin molecules were stacked in the minimum space (Figure 1c). One trajectory at temperature 300K had been performed for 0-200ns. The equilibrium of the system was analyzed by energy and RMSD. The rationality of the model was proved, and the gel model was discussed by the hydrogen bonds between silk fibroin molecules. The secondary structure, radial distribution function, hydrogen bond and other parameters were used to characterize the structure and interactions of the gel.
The model was simulated by molecular dynamics at 300 K with Amber software [29]. The process of energy minimization and heating had been performed firstly, and then performed the equilibrium and product simulation in this work (Figure 1d). In the minimization stage, constrained and non-constrained energy minimization were performed. The steepest descent method [30] was used in the first 2000 steps, and then converted to the 2000-step conjugate gradient method [31,32]. In the minimization of constraining energy, a force of 10.0 kcal/mol was used to constrain the position of the residues. The temperature of system was heated to 300 K from 0 K within 300 ps, and the same force was used for restraint. To prevent the temperature rising too fast, the process was divided into two steps. All bonds involving hydrogen atom were restrained by SHAKE algorithm [33]. The NPT ensemble was used in 5 ns equilibrium simulate. The simulation time of the product simulation is 200 ns with NPT ensemble. The integral time was set to 2 fs and the cutoff distance of non-bonding was set to 10 Å. The sampling interval was 10 ps. Berendsen algorithm [34] was used to control the pressure at 1 bar, and the temperature was controlled with the Langevin method [35] at 300K. Collision frequency for Langevin thermostat was 5.0 ps-1. Pressure relaxation time was set to 10.0 ps. The force field ff14SB [36] and tip3p water cube box of 10.0 Å were used.
Root Mean Square Deviation (RMSD) has been designed to measure the coordinate root mean square difference of a specific atom relative to the reference structure, and to characterize the difference between the structure and the initial structure of protein. RMSD value of α-carbon atom of protein molecular chain is an important basis to measure the conformational change of protein in the simulation process [37]. The smaller the value, the smaller the deviation.
RMSD and potential energy of the system were adopted to verify whether the system has reached equilibrium (Figure 2). The RMSD and potential energy tends to the constant values 12.5 Å and -181000 kcal/mol respectively during 150-200 ns. This indicates that the system reached equilibrium after 150 ns.
The label Total is the amount of every intermolecular hydrogen bonds (Figure 3). The label M2-M1 indicates the hydrogen bonds which donor are provided by molecule 2 and acceptor are provided by molecule 1. The label M1-M2 is the opposite of the label M2-M1. The intermolecular hydrogen bonds are formed between silk fibroin molecules and this suggests that cross-linking is formed [38]. This model can be used as a silk fibroin gel model.
The radial distribution function was estimated by cpptraj of Amber. The high probability indicates that the particle is likely to agglomerate around the given particle [39,40]. The label M1-M2 represents the radial distribution function between O or N atoms of two silk fibroin molecules (Figure 4). The label M-W represents the radial distribution function between O or N atoms of silk fibroin molecules and O atoms of water molecules.
The distance 2.6-3 Å represents hydrogen bonds and the distance 3-5 Å represents van der Waals forces. The Figure 5 shown that hydrogen bonds are formed between two silk fibroin molecules and between silk fibroin molecules and water molecules. Weak non-bond interactions such as van der Waals forces also exist. The radial distribution probability between two silk fibroin molecules is lower than that between silk fibroin molecules and water molecules. This indicates that the weak interaction of former is weaker than that of the later.
There is a gap between two silk fibroin molecules at the beginning and the random coil gradually fold together and fill the gap (Figure 5a). Some α-helix structures are around random coil during the folding process and β-sheets of the two molecules are relatively stable (Figure 5b). These results denote that cross-linking of silk fibroin molecules mainly occurred in random coil segments.
The distribution probability of 2.6-3 Å increases gradually during simulation (Figure 6a). This indicates that the non-bond interaction between silk fibroin molecules increases gradually with the formation and crosslinking of silk fibroin, while the non-bond interaction between silk fibroin and water molecules doesn't change significantly (Figure 6b). These results suggest that cross-linking occurs between silk fibroins. Because non-bond interaction between silk fibroin increases gradually and is weak in the early stage, the non-bond interaction between silk fibroin is weaker than that between silk fibroin and water molecules (Figure 4).
It can be seen from the secondary structure probability diagram (Figure 7) that most of them are anti-parallel β-sheet and random coil, with a small amount of helix structure. It is consistent with the conclusion reported in the literature [41]. The β-sheet is a kind of stable secondary structure, which has a positive effect on the stability of silk fibroin structure [42,43]. Most of the interactions between molecules occur in random coils, thus they play a positive role in the gelation formation.
The structure of each residue can be divided into β-sheet and random coil at the stage of 100-200 ns. For the sake of statistics, more than three consecutive residues with anti-parallel β-sheet are defined as β-sheet segments, and others are defined as random coil segments (Figure 8). The types and number of residues on β-sheet and random coil were counted [44]. Overall, the residues number of random coils is more than that of the β-sheets. There are more nonpolar residues in β-sheet peptide and more polar residues in random coil peptide (Figure 8). Compared with hydrophobic residues, both β-sheet and random coil have more hydrophilic residues. According to the analysis of structure and hydrogen bond, the polar residues of coil and hydrogen bonds of segments promote the aggregation of two silk fibroin molecules and form a cross-linked structure in the gap region.
The total number of intermolecular hydrogen bonds is about 3-5 (Figure 3). The detail of these hydrogen bonds is necessary to discuss (Table 1 and Figure 9). The residues ALA 50 and MET 253 are hydrophobic [44,45,46,47,48], and the others are hydrophilic. The occupancy of hydrogen bond between residue ASN 272 and TYR 241 is more than 50% in the process of 150-200 ns [49]. This result denotes that it is important for the formation of gel. The hydrogen bond between residue GLN 491 and ALA 50 appears after 170 ns and it is consistent with the literature report [50]. Because of the occupancy of intermolecular hydrogen bonds between the three pairs of residues in M2-M1 are more than 76% in the process of 150-200 ns, they are more important for the formation of gel. With the simulation, the number of intermolecular hydrogen bonds increases gradually. This indicates that silk fibroin molecules aggregate gradually, which promotes further formation of osmotic network.
| Position | Acceptor | Donor | Frames during 150-200 ns | Occupancy in 150-200 ns |
| M1-M2 | ASN_272@O | TYR_241@N | 2531 | 50.62% |
| GLN_491@O | ALA_50@N | 653 | 13.06% | |
| M2-M1 | GLN_239@O | ASN_272@N | 3846 | 76.92% |
| SER_252@O | ARG_254@N | 4438 | 88.76% | |
| SER_252@O | MET_253@N | 4290 | 85.80% |
DownLoad:
CSV
In the initial phase of the simulation, a large number of hydrogen bonds formed between the end of silk fibroin molecules and waters (Figure 10). With the number of hydrogen bonds between them decreasing gradually, the number of hydrogen bonds between random coils of silk fibroin molecules increased gradually (Figure 3). There is no interaction among coil segments and they have great freedom to move freely in the early stage. The hydrogen bonds between the coil segments and water can pull silk fibroin molecules to aggregate, shorten the average distance of molecules, and create conditions for intermolecular hydrogen bonds forming. The intermolecular hydrogen bonds between silk fibroin molecules increase gradually and form cross-linked structures.
The probability of hydrophobic residues was 38% at the tail-end of molecule 1 and 42.1% at the head-end of molecule 2. It is easier for the tail-end of molecule 1 to form hydrogen bonds with water molecules. Because of the random coil at the head-end of molecule 2 (First2) is longer than that at the tail-end of molecule 1 (End1), there are more hydrophilic residues at End2, which can form hydrogen bonds easily with water molecules. Therefore, the ability to form intermolecular hydrogen bonds with water molecules has little difference (Figure 10).
The non-bond interaction of the system was calculated by MMPBSA of AmberTools [51,52]. The contribution of each residue belonging to above hydrogen bonds was analyzed according interaction energies. The negative energy indicates attractive force, while positive indicates repulsive. Because the hydrogen bonding is closely related to electrostatic interaction, the secondary electrostatic interaction model was established by Van der Lubbe [53]. It has been used to predict and explain the hydrogen bond strength of self-assembled systems, and the experimental results are often consistent with the prediction of this model. The interaction energies of these residues were sorted according to electrostatic interactions and shown in Table 2.
| Residue | Number | Electrostatic(kcal/mol) | Polar Solvation(kcal/mol) | Non-Polar Solvation(kcal/mol) |
| MET | 253 | -182.147 | 180.327 | -0.869 |
| ARG | 254 | -140.236 | 134.356 | -0.431 |
| SER | 252 | -60.254 | 56.319 | -0.604 |
| GLN | 491 | -11.906 | 10.831 | -0.559 |
| ASN | 272 | -7.943 | 7.722 | -0.424 |
| TYR | 241 | -7.171 | 8.013 | -0.828 |
| GLN | 239 | -4.295 | 5.718 | -0.521 |
| ALA | 50 | -2.895 | 2.844 | -0.461 |
DownLoad:
CSV
The contribution of residues to electrostatic interaction is closely related to interaction of polar solvation. The residues with high electrostatic interaction have high polar solvation energy. The residues SER 252 MET 253 and ARG 254 have the highest electrostatic interaction among these 8 residues. The stronger the electrostatic interaction between residue pairs, the easier to form hydrogen bonds. The residue ALA 50 has the lowest electrostatic interaction, and it is difficult to form hydrogen bonds. Therefore, residues 252-254 can be considered as key residues.
Most of coil are non-sheet peptides, and they are located on gap region in conformation (Figure 5a). The interaction of electrostatic electricity and polar solvation plays an important role in the aggregation of these residues. With the aggregation of silk fibroin molecules, hydrogen bonds cross-and linking among coils are formed gradually and leads to the formation of gel structures. The contribution of electrostatic explains the reason that the silk fibroin solution can be induced to gel rapidly by surfactants [54,55,56].
Diffusion coefficient is an important physical parameter of biomolecules, which defines the diffusivity and permeability of the solvent [57]. The larger the molecular diffusion coefficient, the stronger the molecular migration ability [58]. Diffusion coefficient can not only describe the solution viscosity, but also relate to the intermolecular interaction [59]. Diffusion coefficient D can be estimated by Einstein relation:
| (1) |
where, the unit of D is cm2/s. MSD is mean square displacement from initial position. The dimension degree n is 3 for all MSDs. The slope of MSD is obtained by linear regression. The diffusion coefficient of silk fibroin molecules is about 2.445 × 10-5 cm2/s and that of water is about 2.6624 × 10-5 cm2/s. There is little difference between the diffusion coefficient of water molecules and that of pure water solution (2.275 × 10-5 cm2/s) [60]. The diffusion coefficient of water molecules is slightly higher than that of silk fibroin molecules, but it is still the same order of magnitude, which reflects the good molecular migration ability of silk fibroin hydrogel. The diffusion coefficient of silk fibroin hydrogel in different solvent in is around 0.2-8.1 × 10-7 cm2/s [61,62,63]. The interactions between small molecules and silk fibroin lead to these diffusion coefficients are much lower than that of this work. These results indicate that the diffusivity of silk fibroin hydrogel formation depends on the interactions between molecules and solvents.
In this paper, the gelation process of silk fibroin and the role of non-bond interaction in the gelation process were analyzed by molecular dynamics. The formation of silk fibroin hydrogel depends on non-bond interactions between silk fibroin molecules and solvent. The hydrogen bonds were mainly formed between the random coil peptide fragments at the ends of different molecules. In the early stage of simulation, hydrogen bonds with water were formed, which leads to silk fibroin molecules assemble and hydrogen bonds form. The residues 252-254 of the random coil were identified as key residues according hydrogen bond and non-bond interactions. In the gelation process of silk fibroin, the electrostatic interaction and polar solvation interaction among the random coil segments play an important role in the cross-linking and effect on diffusion coefficient. This is important for further understanding of the gelation mechanism.
We would like to thank the support of the Scientific Research Project of Hubei key laboratory of biomass fiber and ecological dyeing and finishing. This research was funded by the Scientific Research Project of Hubei key laboratory of biomass fiber and ecological dyeing and finishing under, grant number STRZ2018004.
| [1] | X. Li, S. Yan, J. Qu, M. Li, D. Ye, R. You, et al., Soft freezing-induced self-assembly of silk fibroin for tunable gelation, Int. J. Biol. Macromol., 117 (2018), 691-695. |
| [2] |
Y. X. He, N. N. Zhang, W. F. Li, N. Jia, B. Y. Chen, K. Zhou, et al., N-terminal domain of bombyx mori fibroin mediates the assembly of silk in response to pH decrease, J. Mol. Biol., 418 (2012), 197-207. doi: 10.1016/j.jmb.2012.02.040
|
| [3] |
A. E. Thurber, F. G. Omenetto, D. L. Kaplan, In vivo bioresponses to silk proteins, Biomaterials, 71 (2015), 145-157. doi: 10.1016/j.biomaterials.2015.08.039
|
| [4] |
Y. Feng, X. Li, M. Li, D Ye, Q. Zhang, R You, et al., Facile preparation of biocompatible silk fibroin/cellulose nanocomposite films with high mechanical performance, ACS Sustainable Chem. Eng., 5 (2017), 6227-6236. doi: 10.1021/acssuschemeng.7b01161
|
| [5] |
X. Li, J. Zhang, Y. Feng, S. Yan, Q. Zhang, R. You, Tuning the structure and performance of silk biomaterials by combining mulberry and non-mulberry silk fibroin, Polym. Degrad. Stab., 147 (2018), 57-63. doi: 10.1016/j.polymdegradstab.2017.11.013
|
| [6] |
H. Tao, D. L. Kaplan, F. G. Omenetto, Silk materials—a road to sustainable high technology, Adv. Mater., 24 (2012), 2824-2837. doi: 10.1002/adma.201104477
|
| [7] |
A. Motta, L. Fambri, C. Migliaresi, Regenerated silk fibroin films: Thermal and dynamic mechanical analysis, Macromol. Chem. Phys., 203 (2002), 1658-1665. doi: 10.1002/1521-3935(200207)203:10/11<1658::AID-MACP1658>3.0.CO;2-3
|
| [8] |
Z. Shao, F. Vollrath. Surprising strength of silkworm silk, Nature, 418 (2002), 741. doi: 10.1038/418741a
|
| [9] |
H. J. Kim, U. J. Kim, G. Vunjak-Novakovic, B. H. Min, D. L. Kaplan, Influence of macroporous protein scaffolds on bone tissue engineering from bone marrow stem cells, Biomaterials, 26 (2005), 4442-4452. doi: 10.1016/j.biomaterials.2004.11.013
|
| [10] |
S. Cao, R. Tang, G. Sudlow, Z. Wang, K. K. Liu, J. Luan, et al., Shape-dependent biodistribution of biocompatible silk microcapsules, ACS Appl. Mater. Interfaces, 11 (2019), 5499-5508. doi: 10.1021/acsami.8b17809
|
| [11] | J. H. Choi, D. K. Kim, J. E. Song, J. M. Oliveira, R. L. Reis, G. Khang, Silk fibroin-based scaffold for bone tissue engineering, in Novel Biomaterials for Regenerative Medicine, Springer, (2018), 371-387. |
| [12] |
N. H. Romano, D Sengupta, C. Chung, S. C. Heilshorn, Protein-engineered biomaterials: Nanoscale mimics of the extracellular matrix, Biochim. Biophys. Acta, 1810 (2011), 339-349. doi: 10.1016/j.bbagen.2010.07.005
|
| [13] | D. D. McKinnon, D. W. Domaille, T. E. Brown, K. A. Kyburz, E. Kiyotake, J. N. Cha, et al., Measuring cellular forces using bis-aliphatic hydrazone crosslinked stress-relaxing hydrogels, Soft Matter, 10 (2014), 9230-9236. |
| [14] |
R. V. Rughani, M. C. Branco, D. J. Pochan, J. P. Schneider, De novo design of a shear-thin recoverable peptide-based hydrogel capable of intrafibrillar photopolymerization, Macromolecules, 43 (2010), 7924-7930. doi: 10.1021/ma1014808
|
| [15] |
L. Yu, J. Ding, Injectable hydrogels as unique biomedical materials, Chem. Soc. Rev., 37 (2008), 1473-1481. doi: 10.1039/b713009k
|
| [16] |
S. Bai, X. Zhang, Q. Lu, W. Sheng, L. Liu, B. Dong, et al., Reversible hydrogel-solution system of silk with high beta-sheet content, Biomacromolecules, 15 (2014), 3044-3051. doi: 10.1021/bm500662z
|
| [17] |
H. Meng, J. Zheng, X. Wen, Z. Cai, J. Zhang, T. Chen, pH- and sugar-induced shape memory hydrogel based on reversible phenylboronic acid-diol ester bonds, Macromol. Rapid Commun., 36 (2015), 533-537. doi: 10.1002/marc.201400648
|
| [18] |
X. Wang, J. Kluge, G. G. Leisk, D. L. Kaplan, Sonication-induced gelation of silk fibroin for cell encapsulation, Biomaterials, 29 (2008), 1054-1064. doi: 10.1016/j.biomaterials.2007.11.003
|
| [19] | K. Mita, S. Ichimura, T. C. James, Highly repetitive structure and its organization of the silk fibroin gene, J. Mol. Evol., 38 (1994), 583-592. |
| [20] |
G. H. Altman, F. Diaz, C. Jakuba, T. Calabro, R. L. Horan, J. Chen, et al., Silk-based biomaterials, Biomaterials, 24 (2003), 401-416. doi: 10.1016/S0142-9612(02)00353-8
|
| [21] |
B. Liu, Y. Song, L. Jin, Z. Wang, D. Pu, S. Lin, et al., Silk structure and degradation, Colloids Surf. B, 131 (2015), 122-128. doi: 10.1016/j.colsurfb.2015.04.040
|
| [22] |
V. Pandey, T. Haider, P. Jain, P. N. Gupta, V. Soni, Silk as a leading-edge biological macromolecule for improved drug delivery, J. Drug Delivery Sci. Technol., 55 (2020), 101294. doi: 10.1016/j.jddst.2019.101294
|
| [23] |
C. W. P. Foo, D. L. Kaplan, Genetic engineering of fibrous proteins: spider dragline silk and collagen, Adv. Drug Delivery Rev., 54 (2002), 1131-1143. doi: 10.1016/S0169-409X(02)00061-3
|
| [24] | S. Inoue, K. Tanaka, F. Arisaka, S. Kimura, K. Ohtomo, S. Mizuno, Silk fibroin of bombyx mori is secreted, assembling a high molecular mass elementary unit consisting of H-chain, L-chain, and P25, with a 6: 6: 1 molar ratio, J. Biol. Chem., 275 (2000), 40517-40528. |
| [25] |
Z. Chen, H. Zhang, Z. Lin, Y. Lin, J. H. Esch, X. Y. Liu, Programing performance of silk fibroin materials by controlled nucleation, Adv. Funct. Mater., 26 (2016), 8978-8990. doi: 10.1002/adfm.201602908
|
| [26] |
Y. Yang, W. Wu, H. Liu, H. Xu, Y. Zhong, L. Zhang, et al. Aggregation behaviors of thermo-responsive methylcellulose in water: A molecular dynamics simulation study, J. Mol. Graphics Modell., 97 (2020), 107554. doi: 10.1016/j.jmgm.2020.107554
|
| [27] |
E. F. Pettersen, T. D. Goddard, C. C. Huang, G. S. Couch, D. M. Greenblatt, E. C. Meng, et al., UCSF Chimera—A visualization system for exploratory research and analysis, J. Comput. Chem., 25 (2004), 1605-1612. doi: 10.1002/jcc.20084
|
| [28] |
Y. Wang, Q. Wei, S. Wang, W. Chai, Y. Zhang, Structural and water diffusion of Poly(acryl amide)/Poly(vinyl alcohol) blend films: Experiment and molecular dynamics simulations, J. Mol. Graphics Modell., 71 (2017), 40-49. doi: 10.1016/j.jmgm.2016.11.001
|
| [29] |
S. Plimpton, Fast parallel algorithms for short-range molecular dynamics, J. Comput. Phys., 117 (1995), 1-19. doi: 10.1006/jcph.1995.1039
|
| [30] |
H. Attouch, R. Cominetti, A dynamical approach to convex minimization coupling approximation with the steepest descent method, J. Diff. Equations, 128 (1996), 519-540. doi: 10.1006/jdeq.1996.0104
|
| [31] | W. H. Press, B. P. Flannery, S. A. Teukolsky, W. T. Vetterling, P. B. Kramer, Numerical recipes, The art of scientific computing, Phys. Today, 40 (1987), 120-122. |
| [32] |
B. T. Polyak, The conjugate gradient method in extremal problems, Ussr Comput. Math. Math. Phys., 9 (1969), 94-112. doi: 10.1016/0041-5553(69)90035-4
|
| [33] |
V. Krä utler, W. F. Van Gunsteren, P. H. Hünenberger, A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations, J. Comput. Chem., 22 (2001), 501-508. doi: 10.1002/1096-987X(20010415)22:5<501::AID-JCC1021>3.0.CO;2-V
|
| [34] |
H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. DiNola, J. R. Haak, Molecular-dynamics with coupling to an external bath, J. Chem. Phys., 81 (1984), 3684-3690. doi: 10.1063/1.448118
|
| [35] |
S. E. Feller, Y. Zhang, R. W. Pastor, Constant pressure molecular dynamics simulation: The Langevin piston method, J. Chem. Phys., 103 (1995), 4613-4621. doi: 10.1063/1.470648
|
| [36] |
J. A. Maier, C. Martinez, K. Kasavajhala, L. Wickstrom, K. E. Hauser, C. Simmerling, ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB, J. Chem. Theory Comput., 11 (2015), 3696-3713. doi: 10.1021/acs.jctc.5b00255
|
| [37] |
M. P. Eastwood, C. Hardin, Z. Luthey-Schulten, P. G. Wolynes, Evaluating protein structure-prediction schemes using energy landscape theory, IBM J. Res. Dev., 45 (2001), 475-497. doi: 10.1147/rd.453.0475
|
| [38] |
H. Zhao, Y. Li, A novel pH/temperature-responsive hydrogel based on tremella polysaccharide and poly(N-isopropylacrylamide), Colloids Surf. A, 586 (2020), 124270. doi: 10.1016/j.colsurfa.2019.124270
|
| [39] |
S. Mathesan, A. Rath, P. Ghosh, Molecular mechanisms in deformation of cross-linked hydrogel nanocomposite, Mater. Sci. Eng. C, 59 (2016), 157-167. doi: 10.1016/j.msec.2015.09.087
|
| [40] |
S. G. Lee, G. F. Brunello, S. S. Jang, D. G. Bucknall, Molecular dynamics simulation study of P (VP-co-HEMA) hydrogels: effect of water content on equilibrium structures and mechanical properties, Biomaterials, 30 (2009), 6130-6141. doi: 10.1016/j.biomaterials.2009.07.035
|
| [41] |
R. E. Marsh, R. B. Corey, L. Pauling, An investigation of the structure of silk fibroin, Biochim. Biophys. Acta, 16 (1955), 1-34. doi: 10.1016/0006-3002(55)90178-5
|
| [42] |
G. Freddi, P. Monti, M. Nagura, Y. Gotoh, M. Tsukada, Structure and molecular conformation of tussah silk fibroin films: Effect of heat treatment, J. Polym. Sci. Part B, 35 (1997), 841-847. doi: 10.1002/(SICI)1099-0488(19970415)35:5<841::AID-POLB13>3.0.CO;2-A
|
| [43] |
M. Tsukada, G. Freddi, N. Kasai, P. Monti, Structure and molecular conformation of tussah silk fibroin films treated with water-methanol solutions: Dynamic mechanical and thermomechanical behavior, J. Polym. Sci. Part B, 36 (1998), 2717-2724. doi: 10.1002/(SICI)1099-0488(19981115)36:15<2717::AID-POLB5>3.0.CO;2-O
|
| [44] |
W. R. Taylor, The classification of amino acid conservation, J. Theor. Biol., 119 (1986), 205-218. doi: 10.1016/S0022-5193(86)80075-3
|
| [45] |
J. W. West, D. E. Patton, T. Scheuer, Y. Wang, A. L. Goldin, W. A. Catterall, A cluster of hydrophobic amino acid residues required for fast Na(+)-channel inactivation, Proc. Natl. Acad. Sci., 89 (1992), 10910-10914. doi: 10.1073/pnas.89.22.10910
|
| [46] |
J. L. Regier, F. Shen, S. J. Triezenberg, Pattern of aromatic and hydrophobic amino acids critical for one of two subdomains of the VP16 transcriptional activator, Proc. Natl. Acad. Sci., 90 (1993), 883-887. doi: 10.1073/pnas.90.3.883
|
| [47] | J. Yu, H. K. Chan, T. Gengenbach, J. A Denman, Protection of hydrophobic amino acids against moisture-induced deterioration in the aerosolization performance of highly hygroscopic spray-dried powders, Eur. J. Pharm. Biopharm., 119 (2017), 224-234. |
| [48] |
K. Chen, Y. Jiang, L. Du, L. Kurgan, Prediction of integral membrane protein type by collocated hydrophobic amino acid pairs, J. Comput. Chem., 30 (2009), 163-172. doi: 10.1002/jcc.21053
|
| [49] | A. T. Zikri, H. D. Pranowo, W. Haryadi, Stability, hydrogen bond occupancy analysis and binding free energy calculation from flavonol docked in DAPK1 active site using molecular dynamic simulation approaches, Indones. J. Chem., 21 (2021). |
| [50] |
C. H. Görbitz, Hydrogen-bond distances and angles in the structures of amino acids and peptides, Acta Crystallogr. Sec. B, 45 (1989), 390-395. doi: 10.1107/S0108768189003666
|
| [51] |
B. R. Miller, T. D. McGee, J. M. Swails, N. Homeyer, H. Gohlke, A. E. Roitberg, MMPBSA.py: An efficient program for end-state free energy calculations, J. Chem. Theory Comput., 8 (2012), 3314-3321. doi: 10.1021/ct300418h
|
| [52] |
M. G. Montalbán, S. Chakraborty, J. Peñ a-García, H. Verli, G. Villora, H. Pérez-Sánchez, et al., Molecular insight into silk fibroin based delivery vehicle for amphiphilic drugs: Synthesis, characterization and molecular dynamics studies, J. Mol. Liq., 299 (2020), 112156. doi: 10.1016/j.molliq.2019.112156
|
| [53] |
S. C. C. van der Lubbe, F. Zaccaria, X. Sun, C. F. Guerra, Secondary electrostatic interaction model revised: Prediction comes mainly from measuring charge accumulation in hydrogen-bonded monomers, J. Am. Chem. Soc., 141 (2019), 4878-4885. doi: 10.1021/jacs.8b13358
|
| [54] |
X. Wu, J. Hou, M. Li, J. Wang, D. L. Kaplan, S. Lu, Sodium dodecyl sulfate-induced rapid gelation of silk fibroin, Acta Biomater., 8 (2012), 2185-2192. doi: 10.1016/j.actbio.2012.03.007
|
| [55] |
J. H. Park, M. H. Kim, L. Jeong, D. Cho, O. H. Kwon, W. H. Park, Effect of surfactants on sol-gel transition of silk fibroin, J. Sol. Gel. Sci. Technol., 71 (2014), 364-371. doi: 10.1007/s10971-014-3379-4
|
| [56] |
Z. Li, Z. Zheng, Y. Yang, G. Fang, J. Yao, Z. Shao, et al., Robust protein hydrogels from silkworm silk, ACS Sustainable Chem. Eng., 4 (2016), 1500-1506. doi: 10.1021/acssuschemeng.5b01463
|
| [57] |
H. Mrad, S. Alix, S. Migneault, A. Koubaa, P. Perré, Numerical and experimental assessment of water absorption of wood-polymer composites, Measurement, 115 (2018), 197-203. doi: 10.1016/j.measurement.2017.10.011
|
| [58] | W. Fan, P. Zhao, J. Kang, G. Nan, Application of molecular simulation technology to emulsified asphalt study, J. China Univ. Pet., 38 (2014), 179-185. |
| [59] |
L. He, B. Niemeyer, A novel correlation for protein diffusion coefficients based on molecular weight and radius of gyration, Biotechnol. Prog., 19 (2003), 544-548. doi: 10.1021/bp0256059
|
| [60] |
Z. J. Derlacki, A. J. Easteal, A.V. J. Edge, L. A. Woolf, Z. Roksandic, Diffusion coefficients of methanol and water and the mutual diffusion coefficient in methanol-water solutions at 278 and 298 K, J. Phys. Chem. B, 89 (1985), 5318-5322. doi: 10.1021/j100270a039
|
| [61] |
Y. Huang, K. Bailey, S. Wang, X. Feng, Silk fibroin films for potential applications in controlled release, React. Funct. Polym., 116 (2017), 57-68. doi: 10.1016/j.reactfunctpolym.2017.05.007
|
| [62] |
A. Rattanamanee, H. Niamsup, L. Srisombat. W. Punyodom, R. Watanesk, S. Watanesk, Role of chitosan on some physical properties and the urea controlled release of the silk fibroin/gelatin hydrogel, J. Polym. Environ., 23 (2015), 334-340. doi: 10.1007/s10924-014-0703-6
|
| [63] |
T. Asakura, M. Demura, H. Ogawa, K. Matsushita, M. Imanari, NMR imaging of diffusion of small organic molecules in silk fibroin gel, Macromolecules, 24 (1991), 620-622. doi: 10.1021/ma00002a045
|
| 1. | Saphia A. L. Matthew, Refaya Rezwan, Jirada Kaewchuchuen, Yvonne Perrie, F. Philipp Seib, Mixing and flow-induced nanoprecipitation for morphology control of silk fibroin self-assembly, 2022, 12, 2046-2069, 7357, 10.1039/D1RA07764C | |
| 2. | San Seint Seint Aye, Zhi-Hong Zhang, Xin Yu, Heidi Yu, Wen-Dong Ma, Kai Yang, Xin Liu, Jian Li, Jing-Liang Li, Silk Hydrogel Electrostatically Functionalized with a Polycationic Antimicrobial Peptide: Molecular Interactions, Gel Properties, and Antimicrobial Activity, 2022, 38, 0743-7463, 50, 10.1021/acs.langmuir.1c01312 | |
| 3. | Saphia A. L. Matthew, Refaya Rezwan, Yvonne Perrie, F. Philipp Seib, Volumetric Scalability of Microfluidic and Semi-Batch Silk Nanoprecipitation Methods, 2022, 27, 1420-3049, 2368, 10.3390/molecules27072368 | |
| 4. | Liangyu Li, Zhen Liu, Ronghui Qi, Molecular dynamics simulations in hydrogel research and its applications in energy utilization: A review, 2024, 3, 27729702, 100072, 10.1016/j.enrev.2024.100072 | |
| 5. | Zhuangwei Zhang, Yucheng You, Shuoqi Jiang, Yuncheng Tang, Zhongwei Jia, Jian-Xin Li, In silico screening of dyes for enhanced silk coloring using the feeding approach, 2024, 42, 23525541, 101804, 10.1016/j.scp.2024.101804 |
Xuewei Jiang, Lu Zheng, Huhe Wu, Jun Zhang. Investigations into the role of non-bond interaction on gelation mechanism of silk fibroin hydrogel[J]. Mathematical Biosciences and Engineering, 2021, 18(4): 4071-4083. doi: 10.3934/mbe.2021204
| Position | Acceptor | Donor | Frames during 150-200 ns | Occupancy in 150-200 ns |
| M1-M2 | ASN_272@O | TYR_241@N | 2531 | 50.62% |
| GLN_491@O | ALA_50@N | 653 | 13.06% | |
| M2-M1 | GLN_239@O | ASN_272@N | 3846 | 76.92% |
| SER_252@O | ARG_254@N | 4438 | 88.76% | |
| SER_252@O | MET_253@N | 4290 | 85.80% |
DownLoad:
CSV
| Residue | Number | Electrostatic(kcal/mol) | Polar Solvation(kcal/mol) | Non-Polar Solvation(kcal/mol) |
| MET | 253 | -182.147 | 180.327 | -0.869 |
| ARG | 254 | -140.236 | 134.356 | -0.431 |
| SER | 252 | -60.254 | 56.319 | -0.604 |
| GLN | 491 | -11.906 | 10.831 | -0.559 |
| ASN | 272 | -7.943 | 7.722 | -0.424 |
| TYR | 241 | -7.171 | 8.013 | -0.828 |
| GLN | 239 | -4.295 | 5.718 | -0.521 |
| ALA | 50 | -2.895 | 2.844 | -0.461 |
DownLoad:
CSV
| Position | Acceptor | Donor | Frames during 150-200 ns | Occupancy in 150-200 ns |
| M1-M2 | ASN_272@O | TYR_241@N | 2531 | 50.62% |
| GLN_491@O | ALA_50@N | 653 | 13.06% | |
| M2-M1 | GLN_239@O | ASN_272@N | 3846 | 76.92% |
| SER_252@O | ARG_254@N | 4438 | 88.76% | |
| SER_252@O | MET_253@N | 4290 | 85.80% |
| Residue | Number | Electrostatic(kcal/mol) | Polar Solvation(kcal/mol) | Non-Polar Solvation(kcal/mol) |
| MET | 253 | -182.147 | 180.327 | -0.869 |
| ARG | 254 | -140.236 | 134.356 | -0.431 |
| SER | 252 | -60.254 | 56.319 | -0.604 |
| GLN | 491 | -11.906 | 10.831 | -0.559 |
| ASN | 272 | -7.943 | 7.722 | -0.424 |
| TYR | 241 | -7.171 | 8.013 | -0.828 |
| GLN | 239 | -4.295 | 5.718 | -0.521 |
| ALA | 50 | -2.895 | 2.844 | -0.461 |
