Exploring hub pyroptosis-related genes, molecular subtypes, and potential drugs in ankylosing spondylitis by comprehensive bioinformatics analysis and molecular docking.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:37386410", "kind": "science-record", "revision": 3, "content_hash": "sha256:5999870a2c5c07254e4ea5f6e70aecf39f660d983423566431fda6ceec25545c", "created_at": "2026-11-01T19:02:08Z", "updated_at": "2026-11-01T22:09:47Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Exploring hub pyroptosis-related genes, molecular subtypes, and potential drugs in ankylosing spondylitis by comprehensive bioinformatics analysis and molecular docking.", "doi": "10.1186/s12891-023-06664-8", "pmid": "37386410", "authors": ["Xin Li", "Xiangying Li", "Hongqiang Wang", "Xiang Zhao"], "journal": "BMC musculoskeletal disorders", "year": 2023, "abstract": "BACKGROUND: Ankylosing spondylitis (AS) is a chronic inflammatory autoimmune disease, and the diagnosis and treatment of AS have been limited because its pathogenesis is still unclear. Pyroptosis is a proinflammatory type of cell death that plays an important role in the immune system. However, the relationship between pyroptosis genes and AS has never been elucidated.\nMETHODS: GSE73754, GSE25101, and GSE221786 datasets were collected from the Gene Expression Omnibus (GEO) database. Differentially expressed pyroptosis-related genes (DE-PRGs) were identified by R software. Machine learning and PPI networks were used to screen key genes to construct a diagnostic model of AS. AS patients were clustered into different pyroptosis subtypes according to DE-PRGs using consensus cluster analysis and validated using principal component analysis (PCA). WGCNA was used for screening hub gene modules between two subtypes. Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were used for enrichment analysis to elucidate underlying mechanisms. The ESTIMATE and CIBERSORT algorithms were used to reveal immune signatures. The connectivity map (CMAP) database was used to predict potential drugs for the treatment of AS. Molecular docking was used to calculate the binding affinity between potential drugs and the hub gene.\nRESULTS: Sixteen DE-PRGs were detected in AS compared to healthy controls, and some of these genes showed a significant correlation with immune cells such as neutrophils, CD8 + T cells, and resting NK cells. Enrichment analysis showed that DE-PRGs were mainly related to pyroptosis, IL-1β, and TNF signaling pathways. The key genes (TNF, NLRC4, and GZMB) screened by machine learning and the protein-protein interaction (PPI) network were used to establish the diagnostic model of AS. ROC analysis showed that the diagnostic model had good diagnostic properties in GSE73754 (AUC: 0.881), GSE25101 (AUC: 0.797), and GSE221786 (AUC: 0.713). Using 16 DE-PRGs, AS patients were divided into C1 and C2 subtypes, and these two subtypes showed significant differences in immune infiltration. A key gene module was identified from the two subtypes using WGCNA, and enrichment analysis suggested that the module was mainly related to immune function. Three potential drugs, including ascorbic acid, RO 90-7501, and celastrol, were selected based on CMAP analysis. Cytoscape showed GZMB as the highest-scoring hub gene. Finally, molecular docking results showed that GZMB and ascorbic acid formed three hydrogen bonds, including ARG-41, LYS-40, and HIS-57 (affinity: -5.3 kcal/mol). GZMB and RO-90-7501 formed one hydrogen bond, including CYS-136 (affinity: -8.8 kcal/mol). GZMB and celastrol formed three hydrogen bonds, including TYR-94, HIS-57, and LYS-40 (affinity: -9.4 kcal/mol).\nCONCLUSIONS: Our research systematically analyzed the relationship between pyroptosis and AS. Pyroptosis may play an essential role in the immune microenvironment of AS. Our findings will contribute to a further understanding of the pathogenesis of AS.", "claims": [{"id": "claim-method", "statement": "Finally, receiver operating characteristic (ROC) curve analysis was performed on the GSE73754, GSE25101, and GSE221786 datasets through the “pROC” package [27].", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:37386410", "locator": "sections[9] \"Construction and evaluation of the diagnostic model and nomogram\" paragraph 2", "quote": "Finally, receiver operating characteristic (ROC) curve analysis was performed on the GSE73754, GSE25101, and GSE221786 datasets through the “pROC” package [27]."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The datasets analysed during the present study are available in the GEO repository [GSE73754, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE73754], [GSE25101, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25101], [GSE221786, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE221786].", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:37386410", "locator": "sections[35] \"notes\" paragraph 1", "quote": "The datasets analysed during the present study are available in the GEO repository [GSE73754, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE73754], [GSE25101, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25101], [GSE221786, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE221786]."}], "state": "complete"}, {"id": "claim-target-result", "statement": "ROC analysis showed that the diagnostic model had good diagnostic properties in GSE73754 (AUC: 0.881), GSE25101 (AUC: 0.797), and GSE221786 (AUC: 0.713).", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:37386410", "locator": "sections[0] \"Abstract\" paragraph 3", "quote": "ROC analysis showed that the diagnostic model had good diagnostic properties in GSE73754 (AUC: 0.881), GSE25101 (AUC: 0.797), and GSE221786 (AUC: 0.713)."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "objective": "ROC analysis showed that the diagnostic model had good diagnostic properties in GSE73754 (AUC: 0.881), GSE25101 (AUC: 0.797), and GSE221786 (AUC: 0.713).", "steps": [{"id": "read_published_estimate", "order": 1, "name": "Read the published estimate from the source", "kind": "other", "state": "complete", "blocked_on": null, "description": "The claim, method and evidence span published with the source article.", "duration_hours": null, "parameters": {}, "inputs": ["source:pmid:37386410"], "outputs": ["the published estimate and its evidence"], "instruments": [], "reagents": []}, {"id": "obtain_source_data", "order": 2, "name": "Obtain the source data from the named repository", "kind": "data-preparation", "state": "blocked", "blocked_on": "data:accession_absent", "description": "The GEO series accession resolves and its sample accessions are listed. This lane did not download the dataset or read its licence.", "duration_hours": null, "parameters": {}, "inputs": [], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "complete", "blocked_on": null, "description": "Apply the paper's statistical method to the resolved series and compare the result with the published target.", "duration_hours": null, "parameters": {}, "inputs": [], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": ["step:obtain_source_data blocked on data:accession_absent"]}, "data": {"state": "partial", "assets": [{"id": "source:pmid:37386410", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1186/s12891-023-06664-8", "sha256": "sha256:227d15274f1bc52cb42211f3bf619fbe8bd1de649cc60d989ca15bd4c4d8aebf", "bytes": 48483, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1186/s12891-023-06664-8", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC10308648", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}], "gaps": []}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "principal_component_analysis", "specific": "pca", "label": "Principal component analysis", "dsl_operations": ["pca_fit"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["wide_numeric_matrix"], "evidence": {"score": 0.7, "band": "medium", "quote": "principal component analysis"}}, {"id": "method-2", "family": "clustering", "specific": "kmeans", "label": "Clustering", "dsl_operations": ["kmeans_fit"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["feature_matrix"], "evidence": {"score": 0.8, "band": "high", "quote": "K-means"}}, {"id": "method-3", "family": "network_analysis", "specific": "ppi_network", "label": "Network / pathway analysis", "dsl_operations": ["build_graph", "centrality_scores", "pathway_enrich"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["edge_list"], "evidence": {"score": 0.7, "band": "medium", "quote": "PPI"}}, {"id": "method-4", "family": "network_analysis", "specific": "pathway_enrichment", "label": "Network / pathway analysis", "dsl_operations": ["pathway_enrich"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["edge_list"], "evidence": {"score": 0.8, "band": "high", "quote": "GO enrichment"}}, {"id": "method-5", "family": "diagnostic_evaluation", "specific": "roc", "label": "Diagnostic test evaluation", "dsl_operations": ["confusion_counts", "roc_curve", "sensitivity_specificity", "solve_confusion_matrix", "classification_metrics", "published_precision_envelope"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["reference_standard_and_test"], "evidence": {"score": 0.7, "band": "medium", "quote": "ROC analysis"}}, {"id": "method-6", "family": null, "specific": null, "label": "receiver operating characteristic (ROC) curve analysis", "dsl_operations": [], "missing_operations": [], "expressible": null, "parameters": {}, "data_shapes": [], "evidence": {"score": 1.0, "band": "verbatim", "quote": "Finally, receiver operating characteristic (ROC) curve analysis was performed on the GSE73754, GSE25101, and GSE221786 datasets through the “pROC” package [27]."}}], "targets": [], "gaps": ["method.reader_named_unmapped"]}, "capability": {"state": "complete", "requirements": [{"capability_id": "capability:computational_analysis", "label": "Computational analysis platform", "assay_families": ["meta_analysis", "machine_learning", "survival_analysis", "statistical_regression", "pca", "molecular_docking", "network_analysis", "differential_expression", "gwas"], "availability_state": "owned", "tier_unlocked": "R1", "candidate": "Automated Assays Laravel application and simulated adapters"}], "gaps": []}, "readiness": {"state": "blocked", "tier": "R4", "blocking_axis": "data", "blocked_on": "data:accession_absent", "axes": {"method": {"tier": "R1", "expressible": true, "required_operations": ["pca_fit"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: pca."}, "capability": {"tier": "R1", "capability_id": "capability:computational_analysis", "availability_state": "owned", "basis": "The computational analysis is owned. The physical bench the paper used is not, so a physical validation cannot be run here."}, "data": {"tier": "R4", "reachable": false, "accessions": [], "resolved_accessions": [], "statements": ["The datasets analysed during the present study are available in the GEO repository [GSE73754, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE73754], [GSE25101, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25101], [GSE221786, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE221786]."], "basis": "The paper names no accession to its source data, so the data has to come from the authors."}}, "can_run": {"dsl_operations": ["pca_fit"], "assets": ["source:pmid:37386410"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:None"]}, "gaps": ["data:accession_absent"]}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T18:45:25Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 12899, "completion_tokens": 1447, "cached_tokens": 256}}, "human_critique": null, "stages": [{"id": "extraction", "state": "complete", "reader": "MiniMax-M3.1-Flash", "provider": "minimax"}], "gaps": ["review.independent_reader", "review.machine_critique", "review.human_critique"]}, "execution": {"state": "absent", "status": "not_run", "mode": null, "runs": [], "gaps": ["execution.not_run"]}, "provenance": {"entities": [{"id": "pmid:37386410", "type": "paper", "label": "pmid:37386410"}, {"id": "source:pmid:37386410", "type": "manuscript", "label": "PMID 37386410 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:37386410", "type": "paper", "label": "pmid:37386410"}, {"id": "source:pmid:37386410", "type": "manuscript", "label": "PMID 37386410 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-37386410", "type": "extraction", "used": ["source:pmid:37386410"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-None", "type": "accession_resolution", "used": ["data:geo:None"], "agent": "agent:corpus-assay-engine"}], "agents": [{"id": "agent:corpus-assay-engine", "type": "software", "label": "Corpus assay engine"}, {"id": "agent:minimax-reader", "type": "software", "label": "MiniMax literature reader"}], "source_run_id": "reading-sci-2026-11-01", "dsl_catalogue_fingerprint": "343fb30756ba6ec3", "capability_table_fingerprint": "2cb3e42f047d2ecf", "snapshot_id": "ba34b49490d9bfaf"}, "lineage": {"parent_id": "pmid:37386410", "parent_hash": "sha256:b86708f3d98dc2b91e1c4b394c19e9f06486db2ce9b0211050ab22dda675cc78", "fork_of": null, "generation": 2}, "gaps": ["execution.not_run", "method.reader_named_unmapped", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:obtain_source_data blocked on data:accession_absent"], "assays": [{"role": "analysis", "definition_id": "principal_component_analysis", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["pca_fit"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 0.7, "evidence": ["principal component analysis"], "order_basis": "analysis_only", "id": "assay:1"}]}