Combined Multiomics and In Silico Approach Uncovers PRKAR1A as a Putative Therapeutic Target in Multi-Organ Dysfunction Syndrome

[Image: see text] Despite all epidemiological, clinical, and experimental research efforts, therapeutic concepts in sepsis and sepsis-induced multi-organ dysfunction syndrome (MODS) remain limited and unsatisfactory. Currently, gene expression data sets are widely utilized to discover new biomarkers and therapeutic targets in diseases. In the present study, we analyzed MODS expression profiles (comprising 13 sepsis and 8 control samples) retrieved from NCBI-GEO and found 359 differentially expressed genes (DEGs), among which 170 were downregulated and 189 were upregulated. Next, we employed the weighted gene co-expression network analysis (WGCNA) to establish a MODS-associated gene co-expression network (weighted) and identified representative module genes having an elevated correlation with age. Based on the results, a turquoise module was picked as our hub module. Further, we constructed the PPI network comprising 35 hub module DEGs. The DEGs involved in the highest-confidence PPI network were utilized for collecting pathway and gene ontology (GO) terms using various libraries. Nucleotide di- and triphosphate biosynthesis and interconversion was the most significant pathway. Also, 3 DEGs within our PPI network were involved in the top 5 significantly enriched ontology terms, with hypercortisolism being the most significant term. PRKAR1A was the overlapping gene between top 5 significant pathways and GO terms, respectively. PRKAR1A was considered as a therapeutic target in MODS, and 2992 ligands were screened for binding with PRKAR1A. Among these ligands, 3 molecules based on CDOCKER score (molecular dynamics simulated-based score, which allows us to rank the binding poses according to their quality and to identify the best pose for each system) and crucial interaction with human PRKAR1A coding protein and protein kinase-cyclic nucleotide binding domains (PKA RI alpha CNB-B domain) via active site binding residues, viz. Val283, Val302, Gln304, Val315, Ile327, Ala336, Ala337, Val339, Tyr373, and Asn374, were considered as lead molecules.