Exploring the hub gene CERS6 as a therapeutic target in type 1 diabetes through a bioinformatics and network analyst approach.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:41326554", "kind": "science-record", "revision": 3, "content_hash": "sha256:8ba8bdff43add055bd34a2bff488a5cd12900a2ec31cd22ffe07d87b31c05011", "created_at": "2026-11-01T20:31:21Z", "updated_at": "2026-11-01T22:15:06Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Exploring the hub gene CERS6 as a therapeutic target in type 1 diabetes through a bioinformatics and network analyst approach.", "doi": "10.1038/s41598-025-30045-x", "pmid": "41326554", "authors": ["Hassan Ayaz", "Tajamul Hussain", "Asia Nawaz", "Muhammad Suleman", "Salman Alrokayan", "Dilber Uzun Ozsahin", "Khalid Muhammad", "Yasir Waheed"], "journal": "Scientific reports", "year": 2025, "abstract": "Insulin-producing β-cells are destroyed in type 1 diabetes mellitus (T1DM), a chronic autoimmune disease that results in complete insulin insufficiency and metabolic dysfunction. According to a survival study that used p values, some hub genes are important for predicting and diagnosing illness. Scientists have inferred medicines to identify possible therapies that interact with the identified hub genes. The GSE10586 gene expression dataset from the Gene Expression Omnibus (GEO) was used for this investigation, which included 27 samples from 15 healthy controls and 12 diabetic patients. Normalization methods such as variance stabilization normalization (VSN) were used as part of the data pretreatment. A protein‒protein interaction (PPI) network was constructed, principal component analysis (PCA) was performed, heatmaps were created, and the Limma algorithm was used to analyze differential gene expression. Using DAVID v6.8 and KEGG pathway annotations, the functional enrichment of differentially expressed genes (DEGs) was evaluated. Furthermore, a computational study revealed CERS6 to be one of the potential hub genes. Four drugs, methotrexate, eliglustat, myriocin and statin, were the focus of further studies on the basis of predictions made via ChemSpider and PubChem database analysis. To determine the optimal binding positions of these drugs with CERS6, we used molecular docking techniques. The binding affinity of methotrexate was 8.48 kcal/mol, that of myriocin was 7.85 kcal/mol, that of eliglustat was - 6.62 kcal/mol, and that of serine was - 4.90 kcal/mol against the binding pocket's active residues. To determine how consistently each drug interacted with the CERS6 protein over time, molecular dynamics (MD) simulations were run. Throughout the simulation intervals, both medications were confirmed to be stable, with minor alterations in the CERS6 protein loop region. Therefore, the investigation of structure-based drug design has potential for identifying specific therapeutic targets. Ten hub genes were identified via network analysis of differentially expressed genes. These hub genes could serve as novel targets for T1DM detection, prognosis, and targeting. CERS6 exhibited the highest degree of interaction. Methotrexate, eliglustat, myriocin and statins were identified as potential drugs for CERS6. Overall, these findings provide valuable insights that could pave the way for new experimental strategies in T1DM therapy.", "claims": [{"id": "claim-data-access", "statement": "The data can be provided on suitable request to the corresponding author.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:41326554", "locator": "sections[32] \"Data availability\" paragraph 2", "quote": "The data can be provided on suitable request to the corresponding author."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Ten hub genes were identified via network analysis of differentially expressed genes.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:41326554", "locator": "sections[0] \"Abstract\" paragraph 1", "quote": "Ten hub genes were identified via network analysis of differentially expressed genes."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "Ten hub genes were identified via network analysis of differentially expressed genes.", "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:41326554"], "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": "complete", "blocked_on": null, "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": ["data:geo:GSE10586"], "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": ["data:geo:GSE10586"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:41326554", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1038/s41598-025-30045-x", "sha256": "sha256:dbc682dbcdd63b0bcccaebd9f92a84ebd4d34753126cadc88d69d4096f84fd54", "bytes": 74505, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1038/s41598-025-30045-x", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC12775468", "source": "europe_pmc_fulltext", "coverage": "whole_body", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE10586", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE10586", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE10586", "version": null, "note": "The data can be provided on suitable request to the corresponding author.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The data can be provided on suitable request to the corresponding author.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM266963", "GSM266964", "GSM266965", "GSM266966", "GSM266967", "GSM266968", "GSM266969", "GSM266970", "GSM266971", "GSM266972", "GSM266973", "GSM266974", "GSM266975", "GSM266976", "GSM266977", "GSM266978", "GSM266979", "GSM266980", "GSM266981", "GSM266982", "GSM266983", "GSM266984", "GSM266985", "GSM266986", "GSM266987", "GSM266988", "GSM266989"], "sample_count": 27, "group_sizes": null, "group_sizes_basis": null, "platform": "570", "taxon": "Homo sapiens", "title": "Expression data from recent-onset T1D patients and controls", "resolution_checked_at": "2026-11-01T20:16:08Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE10586 not downloaded", "asset:data:geo:GSE10586 licence not read"]}, "analysis": {"state": "complete", "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": "differential_expression", "specific": "limma_voom", "label": "Differential expression analysis", "dsl_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["count_matrix_with_condition"], "evidence": {"score": 0.7, "band": "medium", "quote": "Limma"}}, {"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": "pathway enrichment"}}], "targets": [], "gaps": []}, "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": "complete", "tier": "R1", "blocking_axis": "none", "blocked_on": "nothing:replicable_now", "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": "R1", "reachable": true, "accessions": ["GSE10586"], "resolved_accessions": ["GSE10586"], "statements": ["The data can be provided on suitable request to the corresponding author."], "basis": "The paper names GSE10586 in Gene Expression Omnibus (GEO) database of NCBI and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["pca_fit"], "assets": ["source:pmid:41326554", "data:geo:GSE10586"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE10586"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T20:07:35Z", "independence": "single_reader"}, "evidence_check": {"method": null, "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 18657, "completion_tokens": 1358, "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:41326554", "type": "paper", "label": "pmid:41326554"}, {"id": "source:pmid:41326554", "type": "manuscript", "label": "PMID 41326554 full text"}, {"id": "data:geo:GSE10586", "type": "dataset", "label": "GEO GSE10586"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:41326554", "type": "paper", "label": "pmid:41326554"}, {"id": "source:pmid:41326554", "type": "manuscript", "label": "PMID 41326554 full text"}, {"id": "data:geo:GSE10586", "type": "dataset", "label": "GEO GSE10586"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-41326554", "type": "extraction", "used": ["source:pmid:41326554"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE10586", "type": "accession_resolution", "used": ["data:geo:GSE10586"], "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:41326554", "parent_hash": "sha256:6b5c033a2815554f6f7976a567f641a1482666b6f60655f2050116b283157e78", "fork_of": null, "generation": 2}, "gaps": ["asset:data:geo:GSE10586 licence not read", "asset:data:geo:GSE10586 not downloaded", "execution.not_run", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "data_generation", "definition_id": "chip_seq", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["nscription factors were first identified via ChIP-seq (ChEA) studies, and the main TFs were subseq"], "order_basis": "paper_text_order:6724", "id": "assay:1"}, {"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": "consumes_earlier_output", "id": "assay:2"}]}