{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,26]],"date-time":"2026-03-26T08:02:15Z","timestamp":1774512135354,"version":"3.50.1"},"reference-count":48,"publisher":"Oxford University Press (OUP)","issue":"2","license":[{"start":{"date-parts":[[2026,3,25]],"date-time":"2026-03-25T00:00:00Z","timestamp":1774396800000},"content-version":"vor","delay-in-days":24,"URL":"https:\/\/creativecommons.org\/licenses\/by-nc\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,3,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>Computational methods for therapeutic target discovery face challenges in integrating multi-scale biological data and predicting system-wide therapeutic effects. We present a computational framework that integrates single-cell transcriptomics, weighted gene co-expression network analysis (WGCNA), and computational perturbation simulation to systematically discover novel therapeutic targets. The framework constructs a knowledge graph comprising 29 896 nodes (23 530 genes and 6366 pathways) with 322 136 edges, integrating gene\u2013gene, gene\u2013pathway, and module relationships. Using perturbation simulation algorithms, we systematically explored 265 candidate targets, scoring each based on perturbation magnitude, module response, therapeutic effect, and statistical significance. Applied to single-cell RNA sequencing data from 38 patients (51 415 cells), the framework identified 66 novel therapeutic targets, including nine very high novelty targets. Computational validation demonstrates efficient scalability (knowledge graph construction: &amp;lt;5 min; 265 perturbation simulations: &amp;lt;2 min) and robust performance across different module sizes. This approach represents a novel computational method for systems-level therapeutic target discovery, with generalizable applications to other complex diseases.<\/jats:p>","DOI":"10.1093\/bib\/bbag129","type":"journal-article","created":{"date-parts":[[2026,3,26]],"date-time":"2026-03-26T06:18:05Z","timestamp":1774505885000},"source":"Crossref","is-referenced-by-count":0,"title":["Computational framework for therapeutic target discovery via perturbation simulation: application to cystic fibrosis airway disease"],"prefix":"10.1093","volume":"27","author":[{"given":"George","family":"Sun","sequence":"first","affiliation":[{"name":"High Dimensional Predictive Biological Lab is located at 1 Lampe Drive , North Carolina State University, Raleigh, NC, 27695 ,","place":["United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yi-Hui","family":"Zhou","sequence":"additional","affiliation":[{"name":"High Dimensional Predictive Biological Lab is located at 1 Lampe Drive , North Carolina State University, Raleigh, NC, 27695 ,","place":["United States"]},{"name":"Departments of Biological Sciences and Statistics , North Carolina State University, 1 Lampe Drive, North Carolina State University, Raleigh, NC, 27695 ,","place":["United States"]}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2026,3,25]]},"reference":[{"key":"2026032602180214500_ref1","doi-asserted-by":"publisher","first-page":"682","DOI":"10.1038\/nchembio.118","article-title":"Network pharmacology: the next paradigm in drug discovery","volume":"4","author":"Hopkins","year":"2008","journal-title":"Nat Chem Biol"},{"key":"2026032602180214500_ref2","doi-asserted-by":"publisher","first-page":"56","DOI":"10.1038\/nrg2918","article-title":"Network medicine: a network-based approach to human disease","volume":"12","author":"Barab\u00e1si","year":"2011","journal-title":"Nat Rev Genet"},{"key":"2026032602180214500_ref3","doi-asserted-by":"publisher","first-page":"1119","DOI":"10.1038\/nbt1338","article-title":"Drug\u2013target network","volume":"25","author":"Yildirim","year":"2007","journal-title":"Nat Biotechnol"},{"key":"2026032602180214500_ref4","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1038\/s41467-018-05116-5","article-title":"Network-based approach to prediction and population-based validation of in silico drug repurposing","volume":"9","author":"Cheng","year":"2012","journal-title":"Nat Commun"},{"key":"2026032602180214500_ref5","doi-asserted-by":"publisher","first-page":"1257601","DOI":"10.1126\/science.1257601","article-title":"Uncovering disease-disease relationships through the incomplete interactome","volume":"347","author":"Menche","year":"2015","journal-title":"Science"},{"key":"2026032602180214500_ref6","doi-asserted-by":"publisher","first-page":"710","DOI":"10.1038\/ng1589","article-title":"An integrative genomics approach to infer causal associations between gene expression and disease","volume":"37","author":"Schadt","year":"2005","journal-title":"Nat Genet"},{"key":"2026032602180214500_ref7","doi-asserted-by":"publisher","first-page":"371","DOI":"10.1038\/s41551-024-01312-5","article-title":"Interpretable identification of cancer genes across biological networks via transformer-powered graph representation learning","volume":"9","author":"Xiao-Rui","year":"2025","journal-title":"Nat Biomed Eng"},{"key":"2026032602180214500_ref8","doi-asserted-by":"publisher","first-page":"121360","DOI":"10.1016\/j.ins.2024.121360","article-title":"Regulation-aware graph learning for drug repositioning over heterogeneous biological network","volume":"686","author":"Zhao","year":"2025","journal-title":"Inf Sci"},{"key":"2026032602180214500_ref9","doi-asserted-by":"publisher","first-page":"4281","DOI":"10.1109\/JBHI.2024.3383591","article-title":"Motif-aware miRNA-disease association prediction via hierarchical attention network","volume":"28","author":"Zhao","year":"2024","journal-title":"IEEE J Biomed Health Inform"},{"key":"2026032602180214500_ref10","doi-asserted-by":"publisher","first-page":"1888","DOI":"10.1016\/j.cell.2019.05.031","article-title":"Comprehensive integration of single-cell data.","volume-title":"Cell","author":"Stuart","year":"2019"},{"key":"2026032602180214500_ref11","doi-asserted-by":"publisher","first-page":"3573","DOI":"10.1016\/j.cell.2021.04.048","article-title":"Integrated analysis of multimodal single-cell data","volume":"184","author":"Hao","year":"2021","journal-title":"Cell"},{"key":"2026032602180214500_ref12","first-page":"190","article-title":"Single-cell reconstruction of human basal cell diversity in normal and idiopathic pulmonary fibrosis lungs","volume":"63","author":"Carraro","year":"2020","journal-title":"Am J Respir Cell Mol Biol"},{"key":"2026032602180214500_ref13","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1186\/1471-2105-9-559","article-title":"WGCNA: an r package for weighted correlation network analysis","volume":"9","author":"Langfelder","year":"2008","journal-title":"BMC Bioinformatics"},{"key":"2026032602180214500_ref14","doi-asserted-by":"publisher","first-page":"1","DOI":"10.2202\/1544-6115.1128","article-title":"A general framework for weighted gene co-expression network analysis","volume":"4","author":"Zhang","year":"2005","journal-title":"Stat Appl Genet Mol Biol"},{"key":"2026032602180214500_ref15","doi-asserted-by":"crossref","DOI":"10.1093\/acprof:oso\/9780199206650.001.0001","volume-title":"Networks: An Introduction","author":"Newman","year":"2010"},{"key":"2026032602180214500_ref16","doi-asserted-by":"publisher","first-page":"45","DOI":"10.1038\/nrg3849","article-title":"Cystic fibrosis genetics: from molecular understanding to clinical application","volume":"16","author":"Cutting","year":"2015","journal-title":"Nat Rev Genet"},{"key":"2026032602180214500_ref17","doi-asserted-by":"publisher","first-page":"1066","DOI":"10.1126\/science.2475911","article-title":"Identification of the cystic fibrosis gene: cloning and characterization of complementary DNA","volume":"245","author":"Riordan","year":"1989","journal-title":"Science"},{"key":"2026032602180214500_ref18","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1038\/nrdp.2015.49","article-title":"Cystic fibrosis","volume":"1","author":"Cohen","year":"2015","journal-title":"Nat Rev Dis Primers"},{"key":"2026032602180214500_ref19","doi-asserted-by":"crossref","first-page":"65","DOI":"10.1016\/S2213-2600(19)30337-6","article-title":"The future of cystic fibrosis care: a global perspective","volume":"8","author":"Bell","year":"2020","journal-title":"Lancet Respir Med"},{"key":"2026032602180214500_ref20","doi-asserted-by":"publisher","first-page":"1324","DOI":"10.1164\/rccm.202209-1653OC","article-title":"Genetic modifiers of cystic fibrosis lung disease severity: whole-genome analysis of 7,840 patients","volume":"207","author":"Zhou","year":"2023","journal-title":"Am J Respir Crit Care Med"},{"key":"2026032602180214500_ref21","doi-asserted-by":"publisher","first-page":"1663","DOI":"10.1056\/NEJMoa1105185","article-title":"A CFTR potentiator in patients with cystic fibrosis and the G551D mutation","volume":"365","author":"Ramsey","year":"2011","journal-title":"N Engl J Med"},{"key":"2026032602180214500_ref22","doi-asserted-by":"publisher","first-page":"220","DOI":"10.1056\/NEJMoa1409547","article-title":"Lumacaftor\u2013ivacaftor in patients with cystic fibrosis homozygous for Phe508del CFTR","volume":"373","author":"Wainwright","year":"2015","journal-title":"N Engl J Med"},{"key":"2026032602180214500_ref23","doi-asserted-by":"publisher","first-page":"1809","DOI":"10.1056\/NEJMoa1908639","article-title":"Elexacaftor\u2013tezacaftor\u2013ivacaftor for cystic fibrosis with a single Phe508del allele","volume":"381","author":"Taylor-Cousar","year":"2019","journal-title":"N Engl J Med"},{"key":"2026032602180214500_ref24","doi-asserted-by":"publisher","first-page":"109748","DOI":"10.1016\/j.compbiomed.2025.109748","article-title":"Identifying novel therapeutic targets in cystic fibrosis through advanced single-cell transcriptomics analysis","volume":"187","author":"Sun","year":"2025","journal-title":"Comput Biol Med"},{"key":"2026032602180214500_ref25","doi-asserted-by":"publisher","first-page":"25","DOI":"10.1038\/75556","article-title":"Gene ontology: tool for the unification of biology","volume":"25","author":"Ashburner","year":"2000","journal-title":"Nat Genet"},{"key":"2026032602180214500_ref26","doi-asserted-by":"publisher","first-page":"27","DOI":"10.1093\/nar\/28.1.27","article-title":"KEGG: Kyoto encyclopedia of genes and genomes","volume":"28","author":"Kanehisa","year":"2000","journal-title":"Nucleic Acids Res"},{"key":"2026032602180214500_ref27","doi-asserted-by":"publisher","first-page":"183","DOI":"10.1385\/1-59259-861-7:183","article-title":"Well-differentiated human airway epithelial cell cultures","volume":"107","author":"Fulcher","year":"2005","journal-title":"Methods Mol Med"},{"key":"2026032602180214500_ref28","first-page":"1253","article-title":"Airway epithelial cell differentiation and human disease","volume":"76","author":"John Ross","year":"2019","journal-title":"Cell Mol Life Sci"},{"key":"2026032602180214500_ref29","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1186\/s13059-019-1874-1","article-title":"Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression","volume":"20","author":"Hafemeister","year":"2019","journal-title":"Genome Biol"},{"key":"2026032602180214500_ref30","doi-asserted-by":"publisher","first-page":"1289","DOI":"10.1038\/s41592-019-0619-0","article-title":"Fast, sensitive and accurate integration of single-cell data with harmony","volume":"16","author":"Korsunsky","year":"2019","journal-title":"Nat Methods"},{"key":"2026032602180214500_ref31","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biol"},{"key":"2026032602180214500_ref32","doi-asserted-by":"publisher","first-page":"289","DOI":"10.1111\/j.2517-6161.1995.tb02031.x","article-title":"Controlling the false discovery rate: a practical and powerful approach to multiple testing","volume":"57","author":"Benjamini","year":"1995","journal-title":"J R Stat Soc B Methodol"},{"key":"2026032602180214500_ref33","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1089\/omi.2011.0118","article-title":"clusterProfiler: an R package for comparing biological themes among gene clusters","volume":"16","author":"Guangchuang","year":"2012","journal-title":"OMICS"},{"key":"2026032602180214500_ref34","doi-asserted-by":"crossref","article-title":"Exploring network structure, dynamics, and function using networkX","author":"Hagberg","DOI":"10.25080\/TCWV9851"},{"key":"2026032602180214500_ref35","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1038\/s41467-019-09186-x","article-title":"Network-based prediction of drug combinations","volume":"10","author":"Chen","year":"2019","journal-title":"Nat Commun"},{"key":"2026032602180214500_ref36","doi-asserted-by":"publisher","first-page":"861","DOI":"10.1038\/nrg3074","article-title":"Non-coding RNAs in human disease","volume":"12","author":"Esteller","year":"2011","journal-title":"Nat Rev Genet"},{"key":"2026032602180214500_ref37","doi-asserted-by":"publisher","first-page":"129","DOI":"10.1080\/00031305.2016.1154108","article-title":"The ASA statement on p-values: context, process, and purpose","volume":"70","author":"Wasserstein","year":"2016","journal-title":"Am Stat"},{"key":"2026032602180214500_ref38","doi-asserted-by":"publisher","first-page":"773","DOI":"10.1109\/JBHI.2025.3585290","article-title":"LLM-DDI: leveraging large language models for drug-drug interaction prediction on biomedical knowledge graph","volume":"30","author":"Li","year":"2025","journal-title":"IEEE J Biomed Health Inform"},{"key":"2026032602180214500_ref39","doi-asserted-by":"publisher","first-page":"5730","DOI":"10.1109\/TSMC.2025.3572738","article-title":"Link-based attributed graph clustering via approximate generative Bayesian learning","volume":"55","author":"Yang","year":"2025","journal-title":"IEEE Trans Syst Man Cybern Syst"},{"key":"2026032602180214500_ref40","first-page":"991","article-title":"Long noncoding RNA in genome regulation: prospects and mechanisms","volume":"12","author":"Liang","year":"2015","journal-title":"RNA Biol"},{"key":"2026032602180214500_ref41","doi-asserted-by":"publisher","first-page":"863","DOI":"10.1038\/ncb2090","article-title":"Defective CFTR induces aggressive formation and airway inflammation in cystic fibrosis through ROS-mediated autophagy inhibition","volume":"12","author":"Luciani","year":"2012","journal-title":"Nat Cell Biol"},{"key":"2026032602180214500_ref42","first-page":"1816","article-title":"Mitochondrial dysfunction in cystic fibrosis","volume":"1863","author":"Valdivieso","year":"2017","journal-title":"Biochim Biophys Acta Mol Basis Dis"},{"key":"2026032602180214500_ref43","first-page":"4730","article-title":"Mitochondrial dysfunction in cystic fibrosis airway epithelium","volume":"20","author":"Vergoossen","year":"2019","journal-title":"Int J Mol Sci"},{"key":"2026032602180214500_ref44","first-page":"586","article-title":"Nf-$\\kappa $b signaling in cystic fibrosis","volume":"24","author":"Cohen","year":"2018","journal-title":"Curr Opin Pulm Med"},{"key":"2026032602180214500_ref45","doi-asserted-by":"publisher","first-page":"S28","DOI":"10.1164\/ajrccm.164.supplement_2.2106061","article-title":"Remodeling in asthma and chronic obstructive lung disease","volume":"164","author":"Jeffery","year":"2001","journal-title":"Am J Respir Crit Care Med"},{"key":"2026032602180214500_ref46","first-page":"21","article-title":"Airway remodelling in asthma: from benchside to clinical practice","volume":"10","author":"Bergeron","year":"2003","journal-title":"Can Respir J"},{"key":"2026032602180214500_ref47","doi-asserted-by":"publisher","first-page":"167653","DOI":"10.1109\/ACCESS.2019.2953499","article-title":"A survey on digital twin: definitions, characteristics, applications, and design implications","volume":"7","author":"Barricelli","year":"2019","journal-title":"IEEE Access"},{"key":"2026032602180214500_ref48","doi-asserted-by":"publisher","first-page":"609","DOI":"10.1038\/520609a","article-title":"Personalized medicine: time for one-person trials","volume":"520","author":"Schork","year":"2015","journal-title":"Nature"}],"container-title":["Briefings in Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bib\/article-pdf\/27\/2\/bbag129\/67508928\/bbag129.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bib\/article-pdf\/27\/2\/bbag129\/67508928\/bbag129.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,3,26]],"date-time":"2026-03-26T06:18:08Z","timestamp":1774505888000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bib\/article\/doi\/10.1093\/bib\/bbag129\/8540362"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,3,1]]},"references-count":48,"journal-issue":{"issue":"2","published-print":{"date-parts":[[2026,3,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bib\/bbag129","relation":{},"ISSN":["1467-5463","1477-4054"],"issn-type":[{"value":"1467-5463","type":"print"},{"value":"1477-4054","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,3]]},"published":{"date-parts":[[2026,3,1]]},"article-number":"bbag129"}}