[Identification the key factor of pulmonary fibrosis following silica nanoparticles exposure based on bioinformatics analysis].

Objective: To investigate the main mechanisms of pulmonary fibrosis following silica nanoparticles (SiNPs) exposure through constructing the macrophage-fibroblast model in vitro, which simulated the process of pulmonary fibrosis. Methods: In January 2021, human mononuclear leukemia cells (THP-1) were treated with 0, 25, 50, 100 μg/ml SiNPs for 24 h. The supernatant of THP-1 cells was collected and applied to human embryonic lung fibroblast cells (MRC-5) which divided into control and low, medium and high dose groups at the logarithmic growth stage for 24 h. MRC-5 cell viability was detected by CCK8. The hydroxyproline (Hyp), interleukin 6 (IL-6), interleukin 1 beta (IL-1β) and tumor necrosis factor-alpha (TNF-α) expression were detected in the supernatants of MRC-5. The changed proteins were detected by liquid-phase mass spectrometry in high dose group. GeneCard database were applied to identity the differential pulmonary fibrosis proteins in high dose group. Gene Ontology (GO) was performed to identity the key biological process in differential pulmonary fibrosis proteins of high dose group. The String database was used to construct the protein-protein interactions (PPI) network of differential pulmonary fibrosis proteins. The APP of CytoHubba was applied to calculate the key protein of differential pulmonary fibrosis proteins in PPI network. Correlation coefficients between key differential pulmonary fibrosis proteins were calculated using Pearson correlation analysis. Western blotting was applied to detect the expression of key proteins of differential pulmonary fibrosis proteins in different groups. Results: CCK8 results showed that MRC-5 cell viability was increasing in low, medium and high dose groups compared with control group (P<0.05). The expression levels of Hyp and IL-1β in different group were increased compared with control group, the expression levels of IL-6 and TNF-α were increased in high dose group compared with control group (P<0.05). GeneCard database identified 26 differential pulmonary fibrosis proteins, which were mainly involved in extracellular matrix hydrolysis, cell inflammatory response, tissue repair, cell proliferation, inflammation response by GO analysis. The APP of CytoHubba was calculated that matrix metalloproteinase 9 (MMP9) and tissue inhibitor metalloproteinase 1 (TIMP1) played an important role in PPI network. The results of correlation analysis showed that MMP9 was correlated with the expression of matrix metalloproteinase 1 (MMP1), matrix metalloproteinase 3 (MMP3), TIMP1 and epidermal growth factor receptor (EGFR) (r=0.97, 0.98, 0.94, 0.93, P<0.05). Western blotting results showed that TIMP1 protein expression was increased in low, medium and high dose groups, while MMP9 protein expression was increased only in high dose group (P<0.05) . Conclusion: Differential expression proteins related with pulmonary fibrosis in MRC-5 cells mainly regulate biological processes of extracellular matrix hydrolysis, tissue repair, and cellular inflammation response following SiNPs exposure. MMP9 and TIMP1 may be the key proteins, which affected the fibrosis process in vitro pulmonary fibrosis model.