In silico assessment of human Calprotectin subunits (S100A8/A9) in presence of sodium and calcium ions using Molecular Dynamics simulation approach
Authors:
Nematollah Gheibi aff001; Mohammad Ghorbani aff002; Hanifeh Shariatifar aff003; Alireza Farasat aff001
Authors place of work:
Cellular and Molecular Research Center, Qazvin University of Medical Sciences, Qazvin, Iran
aff001; Department of Nanobiotechnology/ Biophysics, Faculty of Biological Science, Tarbiat Modares University, Tehran, Iran
aff002; Young Researchers and Elite Club, Tehran Medical Sciences, Islamic Azad University, Tehran, Iran
aff003
Published in the journal:
PLoS ONE 14(10)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0224095
Summary
Calprotectin is a heterodimeric protein complex which consists of two subunits including S100A8 and S100A9. This protein has a major role in different inflammatory disease and various types of cancers. In current study we aimed to evaluate the structural and thermodynamic changes of the subunits and the complex in presence of sodium and calcium ions using molecular dynamics (MD) simulation. Therefore, the residue interaction network (RIN) was visualized in Cytoscape program. In next step, to measure the binding free energy, the potential of mean force (PMF) method was performed. Finally, the molecular mechanics Poisson-Boltzmann surface area (MMPBSA) method was applied as an effective tool to calculate the molecular model affinities. The MD simulation results of the subunits represented their structural changes in presence of Ca2+. Moreover, the RIN and Hydrogen bond analysis demonstrated that cluster interactions between Calprotectin subunits in presence of Ca2+ were greater in comparison with Na+. Our findings indicated that the binding free energy of the subunits in presence of Ca2+ was significantly greater than Na+. The results revealed that Ca2+ has the ability to induce structural changes in subunits in comparison with Na+ which lead to create stronger interactions between. Hence, studying the physical characteristics of the human proteins could be considered as a powerful tool in theranostics and drug design purposes.
Keywords:
Simulation and modeling – Protein interactions – Biochemical simulations – Protein structure – Molecular dynamics – Free energy – Interaction networks – Sodium
Introduction
Calprotectin is a heterodimeric protein complex which consists of two subunits including S100A8 and S100A9. This protein has a great role in different inflammatory disease and various types of cancers. It is a heterodimeric complex of EF-hand Ca2+ binding protein. The protein has two subunits which called S100A8 (MRP8/ Calgranulin A) and S100A9 (MRP14/ Calgranulin B) and are fundamentally expressed by myeloid cells. Thus, they are called Myeloid Related Proteins (MRP) [1]. The S100A8 and S100A9 subunits consist of 93 and 113 amino acid residues with the molecular mass of 10.8 and 13.2 kDa respectively [2]. Both of the subunits have two EF-hands (EF-hand I and II) with the structural motifs of helix-loop-helix. The EF-hand I (N terminal motif) comprising helix I, the non-canonical Ca2+ binding loop and helix II while, EF-hand II (C terminal) comprises helix III, the canonical binding loop and helix IV. Thus, this structure provides two Ca2+ binding sites which are characterized by lower affinity of EF-hand domain I than EF-hand domain II. Despite the high structural similarity, discrepancies in residue compounds of their Ca2+ binding loop provide the EF-hand domains various affinities for Ca2+. While main chain carbonyl groups play the major role in binding EF-hand I, the EF-hand II binding carries out via acidic side chains with higher affinity for Ca2+. The linker regions of various lengths have the ability to separate two EF-hands in both proteins [3, 4]. In all S100 proteins family, the C-terminal tail of the S100A9 is the longest one [5]. Although the homology of the amino acid sequences between S100A8 and S100A9 is very low, but the tertiary structures of the different monomers are significantly similar which proves that all S100 proteins contain highly conserved structures [6]. In presence of Ca2+, the formation of heterodimer complex via non-covalent interactions is more probable. The process of Ca2+ binding in each subunits provides the structural conformation in each EF-hand helices. Thus, promoting the structural changes can lead to exposure of hydrophobic residues specially both alpha helices; H-I and H-IV in N and C terminal to create non-covalent interactions including hydrophobic bonding [7, 8]. The (S100A8/A9)2 heterotetramer is specifically created at increased level of Ca2+ [9]. Furthermore, some changes in different ionic concentration regulate several cellular procedures and also changes in cytosolic calcium concentration lead to various cellular processes modulation too. The Ca2+ binding proteins are considered as key molecules in cell cycle monitoring, differentiation and signal transduction. The aforementioned proteins have the ability to interact with target proteins in a calcium dependent manner which are generally a part of signaling pathways which may lead to create particular cellular responses [10]. The Calprotectin has also been mentioned as DAMPs (Damage Associated Molecular Patterns) with the ability of immune mechanism regulation including inflammatory processes, immune cell migration, immune cell adhesion and responses [11]. Previous studies have provided three specific cell surface proteins as functional receptors for S100A8 and S100A9 such as: Receptor for Advanced Glycation End Products (RAGE) [12], Toll-like receptor 4 (TRL4) [13] and Extracellular Matrix Metalloprotease Inducer (EMMPRIN) [14], even though data are conflicting. Moreover, the Calprotectin has been considered as an endogenous ligand when it is secreted by activated neutrophils to extracellular space which can lead to activate various signaling pathways including JAK/stat, PI3 kinase and MAP kinases (ERK, P38 and JNK). The expression of S100A8/A9 seems to be tissue and cell specific which is differentially increased in various abnormal conditions such as: Inflammatory Bowel Disease (IBD), Rheumatoid Arthritis (RA), Cystic Fibrosis (CF), different types of cancers and neurodegenerative disorders [15, 16]. The Molecular Dynamics (MD) simulation process has a major role in recognition of protein-ligand interactions and the protein conformational modifications at the atomic level. In this study, the MD simulation was applied to evaluate the structural and thermodynamic characteristics of S100A8, S100A9 monomers and S100A8/A9 dimer in physiologic concentration of calcium and sodium ions and also to explain the mechanism of S100A8/A9 dimerization. Although the structural changes were considered experimentally, nevertheless some parts of their mechanism of action were not examined thoroughly. In order to investigate thermodynamics and conformational changes of Calprotectin in different ionic conditions, the MD simulation was done. Additionally, the stability of Calprotectin complex was investigated by Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) and Potential of Mean Force (PMF) methods. Considering the importance of Calprotectin in inflammation process and various disease, a better understanding of the structure and its conformational changes could be useful in treatment purposes [15].
Methods
MD simulation
In current study, the MD simulation process was performed using GROMACS program version 5.1 and the CHARMM 36 force field was used for all simulations. The A9 (ID: 5I8N) [17] and A8 (ID: 5HLO) [18] proteins were obtained by RCSB Protein Data Bank (PDB). In this study, the above systems were solvated by transferable intermolecular potential with 3 points (TIP3P) water model in a cubic box with a distance of 10 Å from the furthest atom of the protein. After solvation, Na+ and Cl− ions were inserted to neutralize the system. Then, the concentration of 150 mM NaCl and CaCl2 were added to the systems [4, 19] and the energy minimization was carried out using the steepest descent method. Each system was equilibrated by 1 ns MD simulation in the canonical (NVT) ensemble and 1 ns MD simulation in the isothermal–isobaric (NPT) ensemble using position restraints on the heavy atoms of the protein to allow for the equilibration of the solvent. The Nose–Hoover thermostat constant was used for fixing the temperature of the system at 310 K. To maintain the pressure of the system at fixed 1 bar pressure, the Parrinello–Rahman pressure coupling method was used [20]. The electrostatic interactions were calculated using the Particle Mesh Ewald (PME) method with 1.0 nm short-range electrostatic and van der Waals cutoffs [21, 22]. To evaluate the effects of calcium ion on dimerization process, the A9 with and without calcium was used. Thus, the orientation of the complex (A8/A9, ID: 1XK4) [9] was applied to locate the above subunits (ID: 5I8N, ID: 5HLO) appropriately besides each other. Consequently, the process of 200 and 100 ns MD simulation was carried out for dimer (A8/A9) and monomers (each subunits) with time steps of 2 fs on equilibrated systems respectively. A summary of the MD simulation runs was shown in S1 Table.
RIN calculation
To calculate the residue interaction network (RIN), the average structure of last 10 ns of trajectory system was used. The ring web server [23] was applied to determine the bonds and non-covalent interactions including pi-pi stacking and electrostatic interactions. The Cytoscape program [24] was used to represent the RIN in which the protein amino acids and the bonds are exhibited as nodes and edges.
PMF
Umbrella sampling (US) is a method which could be applied to obtain the free energy pattern which is often referred as PMF along a special reaction coordinate including protein-protein separation distance. Applying physical reaction coordinate can cause further structural insights [25]. In this study the binding energy of A8 and A9 proteins was calculated from PMF via US method. In the first step, the MD simulation was done to drive A9 far away from A8 which was stable during the simulation process. In the next step about 50 configurations were created along the z axis coordinate. The z coordinates of COM (center of mass) interval among A8 and A9 differed by 0.5 Å in each configurations with the force constant of 10 kcal/mol.Å). The equilibration process for each window was done in a period of 10 ns. Moreover, a 10 ns production run was continued for sampling [26]. Consequently, the PMF pattern was obtained by using Weighted Histogram Analysis Method (WHAM) which was performed via GROMACS software as ‘gmx_wham’ command [27].
MMPBSA calculation
The MMPBSA is used as an effective tool to study the molecular models affinities including ligand-protein and protein-protein interactions [28, 29]. To evaluate the binding affinity of S100A8 to S100A9, the binding free energies were measured via MMPBSA to complete the structural analysis. The MMPBSA calculation were carried out on the last 10 ns trajectories via g-mmpbsa tool in GROMACS software [30]. The model of non-polar solvation is based on the solvent-accessible surface area (SASA) with the probe radius 1.4 Å [31]. To understand the MD process efficiently, the root mean square deviation (RMSD), root mean square fluctuation (RMSF), SASA and H-bond of each system were analyzed via GROMACS accessible tools during the simulation process. Consequently, all the images were visualized by Pymol software [32].
Results and discussion
Analysis the structural properties of A8, A9 and A8/A9 complex
The RMSD is an essential parameter in predicting the system equilibration during the simulation process [33–35]. To investigate the stability of the structure, the RMSD of Cα was calculated. The RMSD of A8 and A9 in presence of Na+ and Ca 2+ were analyzed during the simulation process. The RMSD of A8 protein in ionic condition was shown in Fig 1A. As depicted in the figure, there was no difference between the RMSD value of A8 protein in presence of Na+ and Ca2+. The RMSD value of the A8 protein in both ionic conditions was nearly similar. The A8 protein in both systems was equilibrated after 10 ns. The Fig 1B, illustrates the A9 protein in sodium and calcium ionic conditions respectively. The results demonstrate that the A9 protein was equilibrated in presence of Ca2+ and Na+ after 20 ns of MD simulation. The RMSD average for A9 protein in presence of Ca2+ and Na+ was almost 4Å and 2Å respectively. The results showed that the RMSD value of A9 in presence of Ca2+ was higher than Na+ which confirmed that the A9 protein acquired significant structural modifications in presence of Ca2+ than Na+. The RMSD value of A8/A9 complex in presence of Na+ and Ca2+ was illustrated in Fig 1C. As shown in the figure, the complex was equilibrated in presence of Na+ after 20 ns and Ca2+ after 50 ns. The RMSD average of the complex in presence of Na+ and Ca2+ were 12 Å and 3 Å respectively. These findings confirm that the RMSD value of the complex in presence of Ca2+ was lower than the time which the complex was in sodium condition. This was due to the strong interaction between two subunits in presence of Ca2+ which led to the complex stabilization.
The same as RMSD, the Radius of gyration (Rg) reveals the equilibrated manner of the protein and also the protein contraction during the simulation [21, 36, 37]. Thus, in current study, the Rg results confirm that in presence of Na+, the A9 protein is more compact than the time which the complex is in presence of Ca2+ (Fig 2A). These findings verify that broad structural changes are occurred in A9 protein in presence of Ca2+. The Rg results in presence of Na+ and Ca2+ in A8, were shown in Fig 2B. The Rg results of the aforementioned protein in presence of Na+ and Ca2+ were almost similar which revealed that the protein contraction in both conditions was similar. The Fig 2C shows the Rg results in presence of Na+ and Ca2+. As illustrated, the complex was more condensed in presence of Ca2+ than Na+.
The Root-mean-square fluctuations (RMSF) of Cα atoms can provide direct insights in to the protein’s flexibility and structural fluctuations [38, 39]. To study the flexibility of A9 protein in presence of Na+ and Ca2+, the RMSF was calculated. The RMSF of A9 protein in calcium and sodium ionic conditions was shown in Fig 3A. As depicted in the figure, the A9 protein has been reached the higher flexibility in 25–30 and 68–71 amino acid regions in presence of Ca2+ than Na+. Based on these results, in presence of Ca2+, the RMSF value of A9 has a little increase in above regions which are probably the calcium binding sites. It seems that presence of Ca2+ in system and interactions with these parts lead to structural changes of the A9 protein.
The RMSF value of the A8 protein in presence of Na+ and Ca2+ was almost similar (Fig 3B).
The RMSF value of the complex in presence of both ions was shown in Fig 4. The figure showed the comparison of each subunit of the complex in presence of Na+ and Ca2+. All the amino acids of the A8 subunit were less flexible in presence of Ca2+ than Na+ (Fig 4A). Furthermore, the main amino acid in A9 subunit were less flexible in presence of Ca2+ than Na+ (Fig 4B). These results verify that in presence of Ca2+, the subunits undergo a wide interaction.
The SASA method could be applied to evaluate and compare various molecules and conformations and determining the surface which concealed because of oligomerization [40]. Moreover, it indicates the accessible surface of the protein solvent or a part of a protein which exposes to solvent [41]. The SASA profile of A9 protein was illustrated in Fig 5A. Based on the figure, the amino acids of A9 protein (80–90 residues) were exposed in presence of Ca2+. In other words, the accessible surface of protein interaction region in presence of Ca2+ was higher than sodium ion. This parameter in A8 protein in presence of both ions was almost similar (Fig 5B). The accessible surface of the complex in presence of Ca2+ was significantly lower than the time which the complex was in Na+ environment (Fig 5C). Altogether, the RMSD, Rg, RMSF and SASA results confirm an extensive structural changes of the complex in presence of Ca2+ toward Na+.
The amino acids of A9 protein which were exposed in presence of calcium ion were shown in Fig 6.
RIN analysis
RIN, is a protein representation network which can provide the inter-residue and has a major role in proteins structural stability and functions. Hence, the RIN method could be used to evaluate the modifications in protein interaction network [42, 43]. In this study, the amino acids and their interactions were defined as nodes and edges. To provide the internal interaction network of the protein amino acids, the Ring program was used [23]. The resulted changes in internal interaction network of the complex in presence of Ca2+ and Na+ were shown using the program. In internal interaction network of the complex, the interactions of GLU4 and HIS91 were available in presence of Na+ while the broad interactions between the amino acids of the A9 and A8 were observable in presence of Ca2+ (S1 Fig). The analysis of the Calprotectin internal interaction network demonstrates an increased interaction between the A9 and A8 subunits in presence of Ca2+. It seems that the availability of the calcium leads to A9 structural modifications and the accessibility of the required amino acids for stable A9-A8 complex production.
Influence of ions on Calprotectin association
The US is a method which has the ability to calculate the free energy extensively and has long be applied to evaluate the ligand-receptor separation procedures and binding free energy determination [44, 45]. To recognize the interaction between two subunits in various ionic conditions, the PMF value was calculated. The PMF results of each ionic conditions were illustrated in Fig 7. As demonstrated, the binding energy of PMF indicates that the Calprotectin complex is more stable in presence of Ca2+ than Na+. Moreover, the binding energy of two subunits (A9 and A8) in presence of Ca2+ and Na+ was 23.4 and 13.4 kcal/mol respectively. The results verify that the binding energy of Calprotectin subunits in presence of Ca2+ is approximately 2 fold stronger than the time which the complex is in sodium solution. Hence, the binding energy of PMF is in consistent with MMPBSA results with similar report. Additionally, in presence of Ca2+, the Calprotectin subunits undergo structural changes which result in stronger reactions of the complex in compare with Na+. Furthermore, previous studies indicated that in presence of Ca2+, the hydrophobic amino acids of Calprotectin subunits became more accessible [9] which in current study, the structural modifications of Calprotectin subunits resulted in new amino acids accessibility. Thus, new interactions between the subunits were formed.
Binding energy of Calprotectin complex
In previous studies, the MM/PBSA method was considered to evaluate the ligand and receptor activity during the MD simulation [46, 47]. To evaluate the role of Ca2+ and Na+ in binding two subunits to each other, their binding energy was calculated using MM/PBSA method. The results demonstrated that in presence of Ca2+, the binding energy was significantly higher than the time which the complex was in presence of Na+. The binding energy of the Calprotectin subunits in presence of Ca2+ and Na+ were -287.292 kJ/mol and -194.388 kJ/mol respectively. In one study Korndörfer IP et al confirmed that in presence of Ca2+, the LEU9, ILE12, ILE13, PHE68, GLN69, and LEU72 amino acids of A8 and the THR87, TRP88, and HIS91 amino acids of A9 protein play an important role in binding the two subunits to each other [9]. These findings indicate the role of hydrophobic amino acids in binding two subunits to each other. Furthermore, in presence of Ca2+, the hydrophobic amino acids were more accessible and a network of hydrophobic interactions were created between the two subunits which led to a higher stability of the complex. Fig 8, illustrates that in presence of Na+, only the side chain of A8 protein (LEU9) locates in interaction region. However, in presence of Ca2+, the A9 moves toward A8, therefore the interaction region of heterodimer including THR87, TRP88, HIS91 (from A9 protein) and the GLN69 and LEU72 (from A8 protein) become closer to each other.
Influence of ionic conditions on Hydrogen bond interactions of the Calprotectin
The protein intermolecular H-bond formation plays a fundamental role in systems stability. The H-bond network is the exclusive solvent with the capability to monitor the dynamics, stability, structure and also the biomolecular functions [48]. The average number of H-bonds during the simulation when the heterodimer was present in presence of sodium and calcium ions were 3 and 11 respectively. The increase in the number of H-Bonds in presence of calcium indicates a stronger linkage and a higher stability of heterodimer in this manner (Fig 9).
Calcium entrance in to the A9 structure
The distance of calcium ion from A9 subunit was shown in Fig 10. As depicted in the figure, after 30 ns, the calcium ion was located at 0.6 nm distance of loop I (Fig 10A) and this interval was stable during the simulation. Another calcium ion was located at the interval of 0.6 nm from loop II after 15 ns and the distance was constant during the simulation. The 3D image of the calcium ions position in loop I and II was illustrated in Fig 10B. The calcium entrance in loop II in a short time may be due to the presence of negative charged amino acids which are abundant in loop II. (2 Asp in loop II in comparison with 1 Asp in loop I) [49].
Several studies proved that the calcium can induce conformational changes in protein subunits and facilitate the dimerization process [8, 49, 50]. The 3D structure of the subunits and the complex in presence of both ions was shown in Fig 11. In our study, it was demonstrated that the calcium ion can enter the structure (Fig 11A) which was confirmed by previous studies [4] but sodium possesses no binding site for the A9 protein (Fig 11B). The Fig 11C, illustrates the calcium availability in loop I and II of A8 subunit. Furthermore, the A8/A9 complex was shown in presence of Na+ and Ca2+ (Fig 11D and 11E). It should be noted that the calcium ion can induce the structural changes of the subunits which lead to a better biological activity of the complex. Moreover, less structural changes in NaCl in comparison to CaCl2 for A9 protein and the A8/A9 complex was observed.
Conclusion
The Calprotectin has a key role in various kinds of disease including cancers, autoimmune disorders and the involved processes such as inflammation. Hence, a better recognition of the protein and the mechanism of action is needed. In this study, we aimed to evaluate and characterize the protein structure in different ionic conditions. Thus, the MD simulation method was used to study the structural changes of the protein complex in presence of Na+ and Ca2+. The results represented a stronger binding and higher stability of the Calprotectin subunits in presence of Ca2+ which was the consequence of structural changes of the subunits. Moreover, such parameters including: RMSD, RMSF, SASA, RG and H-bonds confirmed the structural changes, while the PMF and MM-PBSA verified the stronger binding of the protein subunits. The results of this study confirm that presence of calcium leads to structural changes of the Calprotectin and eventually the accessibility of hydrophobic amino acids for stronger interactions and higher stability of the complex which declare the importance of ionic conditions in intracellular inflammation procedures. Hence, these findings could be useful in protein structure recognition which could be applied in drug design, diagnosis and targeted treatments of the disease in the future.
Supporting information
S1 Table [docx]
Summary of the MD runs.
Zdroje
1. Leukert N, Sorg C, Roth J. Molecular basis of the complex formation between the two calcium-binding proteins S100A8 (MRP8) and S100A9 (MRP14). Biological chemistry. 2005;386(5):429–34. doi: 10.1515/BC.2005.051 15927886
2. Shirley SH, von Maltzan K, Robbins PO, Kusewitt DF. Melanocyte and melanoma cell activation by calprotectin. Journal of skin cancer. 2014;2014.
3. Vogl T, Gharibyan AL, Morozova-Roche LA. Pro-inflammatory S100A8 and S100A9 proteins: self-assembly into multifunctional native and amyloid complexes. International journal of molecular sciences. 2012;13(3):2893–917. doi: 10.3390/ijms13032893 22489132
4. Reis RAG, Bortot LO, Caliri A. In silico assessment of S100A12 monomer and dimer structural dynamics: implications for the understanding of its metal-induced conformational changes. JBIC Journal of Biological Inorganic Chemistry. 2014;19(7):1113–20. doi: 10.1007/s00775-014-1149-y 24944024
5. Markowitz J, Carson WE III. Review of S100A9 biology and its role in cancer. Biochimica et Biophysica Acta (BBA)-Reviews on Cancer. 2013;1835(1):100–9.
6. Brunjes Brophy M, Nakashige TG, Gaillard A, Nolan EM. Contributions of the S100A9 C-terminal tail to high-affinity Mn (II) chelation by the host-defense protein human calprotectin. Journal of the American Chemical Society. 2013;135(47):17804–17. doi: 10.1021/ja407147d 24245608
7. Santamaria-Kisiel L, Rintala-Dempsey AC, Shaw GS. Calcium-dependent and-independent interactions of the S100 protein family. Biochemical Journal. 2006;396(2):201–14. doi: 10.1042/BJ20060195 16683912
8. Gheibi N, Asghari H, Chegini K, Sahmani M, Moghadasi M. The role of calcium in the conformational changes of the recombinant S100A8/S100A91. Molecular Biology. 2016;50(1):118–23.
9. Korndörfer IP, Brueckner F, Skerra A. The crystal structure of the human (S100A8/S100A9) 2 heterotetramer, calprotectin, illustrates how conformational changes of interacting α-helices can determine specific association of two EF-hand proteins. Journal of molecular biology. 2007;370(5):887–98. doi: 10.1016/j.jmb.2007.04.065 17553524
10. Ichikawa M, Williams R, Wang L, Vogl T, Srikrishna G. S100A8/A9 activate key genes and pathways in colon tumor progression. Molecular cancer research. 2011;9(2):133–48. doi: 10.1158/1541-7786.MCR-10-0394 21228116
11. Gomes LH, Raftery MJ, Yan WX, Goyette JD, Thomas PS, Geczy CL. S100A8 and S100A9—oxidant scavengers in inflammation. Free Radical Biology and Medicine. 2013;58:170–86. doi: 10.1016/j.freeradbiomed.2012.12.012 23277148
12. Leclerc E, Fritz G, Weibel M, Heizmann CW, Galichet A. S100B and S100A6 differentially modulate cell survival by interacting with distinct RAGE (receptor for advanced glycation end products) immunoglobulin domains. Journal of Biological Chemistry. 2007;282(43):31317–31. doi: 10.1074/jbc.M703951200 17726019
13. Ehrchen JM, Sunderkötter C, Foell D, Vogl T, Roth J. The endogenous Toll–like receptor 4 agonist S100A8/S100A9 (calprotectin) as innate amplifier of infection, autoimmunity, and cancer. Journal of leukocyte biology. 2009;86(3):557–66. doi: 10.1189/jlb.1008647 19451397
14. Hibino T, Sakaguchi M, Miyamoto S, Yamamoto M, Motoyama A, Hosoi J, et al. S100A9 is a novel ligand of EMMPRIN that promotes melanoma metastasis. Cancer research. 2013;73(1):172–83. doi: 10.1158/0008-5472.CAN-11-3843 23135911
15. Jung S-Y, Park Y-B, Ha Y-J, Lee K-H, Lee S-K. Serum calprotectin as a marker for disease activity and severity in adult-onset Still’s disease. The Journal of rheumatology. 2010;37(5):1029–34. doi: 10.3899/jrheum.091120 20231196
16. Shabani F, Farasat A, Mahdavi M, Gheibi N. Calprotectin (S100A8/S100A9): a key protein between inflammation and cancer. Inflammation Research. 2018;67(10):801–12. doi: 10.1007/s00011-018-1173-4 30083975
17. Chang C-C, Khan I, Tsai K-L, Li H, Yang L-W, Chou R-H, et al. Blocking the interaction between S100A9 and RAGE V domain using CHAPS molecule: A novel route to drug development against cell proliferation. Biochimica et Biophysica Acta (BBA)-Proteins and Proteomics. 2016;1864(11):1558–69.
18. Lin H, Andersen GR, Yatime L. Crystal structure of human S100A8 in complex with zinc and calcium. BMC structural biology. 2016;16(1):8. doi: 10.1186/s12900-016-0058-4 27251136
19. Batoulis H, Schmidt TH, Weber P, Schloetel J-G, Kandt C, Lang T. Concentration dependent ion-protein interaction patterns underlying protein oligomerization behaviours. Scientific reports. 2016;6:24131. doi: 10.1038/srep24131 27052788
20. Akya A, Farasat A, Ghadiri K, Rostamian M. Identification of HLA-I restricted epitopes in six vaccine candidates of Leishmania tropica using immunoinformatics and molecular dynamics simulation approaches. Infection, Genetics and Evolution. 2019;75:103953. doi: 10.1016/j.meegid.2019.103953 31284043
21. Farasat A, Rahbarizadeh F, Hosseinzadeh G, Sajjadi S, Kamali M, Keihan AH. Affinity enhancement of nanobody binding to EGFR: in silico site-directed mutagenesis and molecular dynamics simulation approaches. Journal of biomolecular structure and dynamics. 2017;35(8):1710–28. doi: 10.1080/07391102.2016.1192065 27691399
22. Ochoa R, Soler MA, Laio A, Cossio P. Assessing the capability of in silico mutation protocols for predicting the finite temperature conformation of amino acids. Physical Chemistry Chemical Physics. 2018;20(40):25901–9. doi: 10.1039/c8cp03826k 30289133
23. Piovesan D, Minervini G, Tosatto SC. The RING 2.0 web server for high quality residue interaction networks. Nucleic acids research. 2016;44(W1):W367–W74. doi: 10.1093/nar/gkw315 27198219
24. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome research. 2003;13(11):2498–504. doi: 10.1101/gr.1239303 14597658
25. Naughton FB, Kalli AC, Sansom MS. Modes of interaction of pleckstrin homology domains with membranes: toward a computational biochemistry of membrane recognition. Journal of molecular biology. 2018;430(3):372–88. doi: 10.1016/j.jmb.2017.12.011 29273202
26. Lemkul JA, Bevan DR. Assessing the stability of Alzheimer’s amyloid protofibrils using molecular dynamics. The Journal of Physical Chemistry B. 2010;114(4):1652–60. doi: 10.1021/jp9110794 20055378
27. Zeng S, Zhou G, Guo J, Zhou F, Chen J. Molecular simulations of conformation change and aggregation of HIV-1 Vpr13-33 on graphene oxide. Scientific reports. 2016;6:24906. doi: 10.1038/srep24906 27097898
28. Sulaiman KO, Kolapo TU, Onawole AT, Islam MA, Adegoke RO, Badmus SO. Molecular dynamics and combined docking studies for the identification of Zaire ebola virus inhibitors. Journal of Biomolecular Structure and Dynamics. 2019;37(12):3029–40. doi: 10.1080/07391102.2018.1506362 30058446
29. Wang C, Nguyen PH, Pham K, Huynh D, Le TBN, Wang H, et al. Calculating protein–ligand binding affinities with MMPBSA: Method and error analysis. Journal of computational chemistry. 2016;37(27):2436–46. doi: 10.1002/jcc.24467 27510546
30. Kumari R, Kumar R, Consortium OSDD, Lynn A. g_mmpbsa∙ A GROMACS tool for high-throughput MM-PBSA calculations. Journal of chemical information and modeling. 2014;54(7):1951–62. doi: 10.1021/ci500020m 24850022
31. Wang Z, Zhang Q, Meng F, Li S, Xu Q, Qi Z. Characterization of the ligand binding of PGRP-L in half-smooth tongue sole (Cynoglossus semilaevis) by molecular dynamics and free energy calculation. Electronic Journal of Biotechnology. 2018;31:93–9.
32. DeLano WL. The PyMOL molecular graphics system. http://wwwpymolorg. 2002.
33. Shen M, Guan J, Xu L, Yu Y, He J, Jones GW, et al. Steered molecular dynamics simulations on the binding of the appendant structure and helix-β2 in domain-swapped human cystatin C dimer. Journal of Biomolecular Structure and Dynamics. 2012;30(6):652–61. doi: 10.1080/07391102.2012.689698 22731964
34. Karain WI. Detecting transitions in protein dynamics using a recurrence quantification analysis based bootstrap method. BMC bioinformatics. 2017;18(1):525. doi: 10.1186/s12859-017-1943-y 29179670
35. Kaur G, Pandey B, Kumar A, Garewal N, Grover A, Kaur J. Drug targeted virtual screening and molecular dynamics of LipU protein of Mycobacterium tuberculosis and Mycobacterium leprae. Journal of Biomolecular Structure and Dynamics. 2019;37(5):1254–69. doi: 10.1080/07391102.2018.1454852 29557724
36. Lobanov MY, Bogatyreva N, Galzitskaya O. Radius of gyration as an indicator of protein structure compactness. Molecular Biology. 2008;42(4):623–8.
37. Anantram A, Kundaikar H, Degani M, Prabhu A. Molecular dynamic simulations on an inhibitor of anti-apoptotic Bcl-2 proteins for insights into its interaction mechanism for anti-cancer activity. Journal of Biomolecular Structure and Dynamics. 2019;37(12):3109–21. doi: 10.1080/07391102.2018.1508371 30526410
38. Chen J, Wang J, Zhu W. Molecular Mechanism and Energy Basis of Conformational Diversity of Antibody SPE7 Revealed by Molecular Dynamics Simulation and Principal Component Analysis. Scientific reports. 2016;6:36900. doi: 10.1038/srep36900 27830740
39. Mahapatra MK, Bera K, Singh DV, Kumar R, Kumar M. In silico modelling and molecular dynamics simulation studies of thiazolidine based PTP1B inhibitors. Journal of Biomolecular Structure and Dynamics. 2018;36(5):1195–211. doi: 10.1080/07391102.2017.1317026 28393626
40. Mitternacht S. FreeSASA: An open source C library for solvent accessible surface area calculations. F1000Research. 2016;5.
41. Chen H, Panagiotopoulos AZ. Molecular Modeling of Surfactant Micellization Using Solvent Accessible Surface Area. Langmuir. 2019.
42. Shukla H, Shukla R, Sonkar A, Tripathi T. Alterations in conformational topology and interaction dynamics caused by L418A mutation leads to activity loss of Mycobacterium tuberculosis isocitrate lyase. Biochemical and biophysical research communications. 2017;490(2):276–82. doi: 10.1016/j.bbrc.2017.06.036 28610921
43. Griffin JW, Bradshaw PC. In silico prediction of novel residues involved in amyloid primary nucleation of human I56T and D67H lysozyme. BMC structural biology. 2018;18(1):9. doi: 10.1186/s12900-018-0088-1 30029603
44. You W, Tang Z, Chang C-eA. Potential mean force from umbrella sampling simulations: what can we learn and what is missed? Journal of chemical theory and computation. 2019.
45. Bowman JD, Lindert S. Molecular Dynamics and Umbrella Sampling Simulations Elucidate Differences in Troponin C Isoform and Mutant Hydrophobic Patch Exposure. The Journal of Physical Chemistry B. 2018;122(32):7874–83. doi: 10.1021/acs.jpcb.8b05435 30070845
46. Liu Z, Zhang Y. Molecular dynamics simulations and MM–PBSA calculations of the lectin from snowdrop (Galanthus nivalis). Journal of molecular modeling. 2009;15(12):1501. doi: 10.1007/s00894-009-0502-5 19449157
47. Kumar A, Srivastava G, Negi AS, Sharma A. Docking, molecular dynamics, binding energy-MM-PBSA studies of naphthofuran derivatives to identify potential dual inhibitors against BACE-1 and GSK-3β. Journal of Biomolecular Structure and Dynamics. 2019;37(2):275–90. doi: 10.1080/07391102.2018.1426043 29310523
48. Mallamace D, Fazio E, Mallamace F, Corsaro C. The Role of Hydrogen Bonding in the Folding/Unfolding Process of Hydrated Lysozyme: A Review of Recent NMR and FTIR Results. International journal of molecular sciences. 2018;19(12):3825.
49. Leukert N, Vogl T, Strupat K, Reichelt R, Sorg C, Roth J. Calcium-dependent tetramer formation of S100A8 and S100A9 is essential for biological activity. Journal of molecular biology. 2006;359(4):961–72. doi: 10.1016/j.jmb.2006.04.009 16690079
50. Streicher WW, Lopez MM, Makhatadze GI. Modulation of quaternary structure of S100 proteins by calcium ions. Biophysical chemistry. 2010;151(3):181–6. doi: 10.1016/j.bpc.2010.06.003 20621410
Článek vyšel v časopise
PLOS One
2019 Číslo 10
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Correction: Low dose naltrexone: Effects on medication in rheumatoid and seropositive arthritis. A nationwide register-based controlled quasi-experimental before-after study
- Combining CDK4/6 inhibitors ribociclib and palbociclib with cytotoxic agents does not enhance cytotoxicity
- Risk factors associated with IgA vasculitis with nephritis (Henoch–Schönlein purpura nephritis) progressing to unfavorable outcomes: A meta-analysis
- Prevalence of pectus excavatum (PE), pectus carinatum (PC), tracheal hypoplasia, thoracic spine deformities and lateral heart displacement in thoracic radiographs of screw-tailed brachycephalic dogs