Differential gene expression analysis combined with molecular dynamics simulation study to elucidate the novel potential biomarker involved in pulmonary TB

Tuberculosis (TB) is a lethal multisystem disease that attacks the lungs' first line of defense. A substantial threat to public health and a primary cause of death is pulmonary TB. This study aimed to identify and investigate the probable differentially expressed genes (DEGs) primarily involved in Pulmonary TB. Accordingly, three inde-pendent gene expression data sets, numbered GSE139825, GSE139871, and GSE54992, were utilized for this purpose. The identified DEGs were used for bioinformatics-based analysis, including physical gene interaction, Gene Ontology (GO), network analysis and pathway studies using the Kyoto Encyclopedia of Genes and Genomes pathway (KEGG). The computational analysis predicted that TNFAIP6 is the significant DEG in the gene expression profiling of TB datasets. According to gene ontology analysis, TNFAIP6 is also essential in injury and inflammation. Further, TNFA1P6 is strongly linked to arsenic poisoning, evident from the results of Networ-kAnalyst, a comprehensive and interactive platform for gene expression profiling via network visual analytics. As a result, the TNFAIP6 gene was ultimately chosen as a candidate DEG and subsequently employed for in silico structural characterization studies. The tertiary structure of TNFAIP6 was modelled using the ROBETTA server, followed by validation with SAVES and ProSA webserver. Additionally, structural dynamic studies, including molecular dynamics simulation (MDS) and essential dynamics analysis, including principal component (PC) based free energy landscape (FEL) analysis, was used for checking the stability of TNFAIP6 models. The dynamics result established the structural rigidity of modelled TNFAIP6 through RMSD, RMSF and RoG results. The FEL analysis revealed the restricted conformational flexibility of TNFAIP6 by displaying a single minimum energy basin in the contour plot. The comprehensive computational analysis established that TNFAIP6 could serve as a viable biomarker to assess the severity of pulmonary TB.

Automated summary

This study aimed to identify and investigate the differentially expressed genes (DEGs) primarily involved in Pulmonary TB. TNFAIP6 was identified as a significant DEG in TB datasets, and it was found to be closely associated with injury, inflammation, and arsenic poisoning. Further structural characterization studies confirmed the stability of TNFAIP6 models and revealed its potential as a viable biomarker for assessing the severity of pulmonary TB.