Systems biology-based analysis exploring shared biomarkers and pathogenesis of myocardial infarction combined with osteoarthritis.

Background: More and more evidence supports the association between myocardial infarction (MI) and osteoarthritis (OA). The purpose of this study is to explore the shared biomarkers and pathogenesis of MI complicated with OA by systems biology. Methods: Gene expression profiles of MI and OA were downloaded from the Gene Expression Omnibus (GEO) database. The Weighted Gene Co-Expression Network Analysis (WGCNA) and differentially expressed genes (DEGs) analysis were used to identify the common DEGs. The shared genes related to diseases were screened by three public databases, and the protein-protein interaction (PPI) network was built. GO and KEGG enrichment analyses were performed on the two parts of the genes respectively. The hub genes were intersected and verified by Least absolute shrinkage and selection operator (LASSO) analysis, receiver operating characteristic (ROC) curves, and single-cell RNA sequencing analysis. Finally, the hub genes differentially expressed in primary cardiomyocytes and chondrocytes were verified by RT-qPCR. The immune cell infiltration analysis, subtypes analysis, and transcription factors (TFs) prediction were carried out. Results: In this study, 23 common DEGs were obtained by WGCNA and DEGs analysis. In addition, 199 common genes were acquired from three public databases by PPI. Inflammation and immunity may be the common pathogenic mechanisms, and the MAPK signaling pathway may play a key role in both disorders. DUSP1, FOS, and THBS1 were identified as shared biomarkers, which is entirely consistent with the results of single-cell RNA sequencing analysis, and furher confirmed by RT-qPCR. Immune infiltration analysis illustrated that many types of immune cells were closely associated with MI and OA. Two potential subtypes were identified in both datasets. Furthermore, FOXC1 may be the crucial TF, and the relationship of TFs-hub genes-immune cells was visualized by the Sankey diagram, which could help discover the pathogenesis between MI and OA. Conclusion: In summary, this study first revealed 3 (DUSP1, FOS, and THBS1) novel shared biomarkers and signaling pathways underlying both MI and OA. Additionally, immune cells and key TFs related to 3 hub genes were examined to further clarify the regulation mechanism. Our study provides new insights into shared molecular mechanisms between MI and OA.

A practically efficient algorithm for identifying critical control proteins in directed probabilistic biological networks.

Network controllability is unifying the traditional control theory with the structural network information rooted in many large-scale biological systems of interest, from intracellular networks in molecular biology to brain neuronal networks. In controllability approaches, the set of minimum driver nodes is not unique, and critical nodes are the most important control elements because they appear in all possible solution sets. On the other hand, a common but largely unexplored feature in network control approaches is the probabilistic failure of edges or the uncertainty in the determination of interactions between molecules. This is particularly true when directed probabilistic interactions are considered. Until now, no efficient algorithm existed to determine critical nodes in probabilistic directed networks. Here we present a probabilistic control model based on a minimum dominating set framework that integrates the probabilistic nature of directed edges between molecules and determines the critical control nodes that drive the entire network functionality. The proposed algorithm, combined with the developed mathematical tools, offers practical efficiency in determining critical control nodes in large probabilistic networks. The method is then applied to the human intracellular signal transduction network revealing that critical control nodes are associated with important biological features and perturbed sets of genes in human diseases, including SARS-CoV-2 target proteins and rare disorders. We believe that the proposed methodology can be useful to investigate multiple biological systems in which directed edges are probabilistic in nature, both in natural systems or when determined with large uncertainties in-silico.

Non-dialyzable uremic toxins and renal tubular cell damage in CKD patients: a systems biology approach.

BACKGROUND: Chronic kidney disease presents global health challenges, with hemodialysis as a common treatment. However, non-dialyzable uremic toxins demand further investigation for new therapeutic approaches. Renal tubular cells require scrutiny due to their vulnerability to uremic toxins. METHODS: In this study, a systems biology approach utilized transcriptomics data from healthy renal tubular cells exposed to healthy and post-dialysis uremic plasma. RESULTS: Differential gene expression analysis identified 983 up-regulated genes, including 70 essential proteins in the protein-protein interaction network. Modularity-based clustering revealed six clusters of essential proteins associated with 11 pathological pathways activated in response to non-dialyzable uremic toxins. CONCLUSIONS: Notably, WNT1/11, AGT, FGF4/17/22, LMX1B, GATA4, and CXCL12 emerged as promising targets for further exploration in renal tubular pathology related to non-dialyzable uremic toxins. Understanding the molecular players and pathways linked to renal tubular dysfunction opens avenues for novel therapeutic interventions and improved clinical management of chronic kidney disease and its complications.

Transient frequency preference responses in cell signaling systems.

Ligand-receptor systems, covalent modification cycles, and transcriptional networks are the fundamental components of cell signaling and gene expression systems. While their behavior in reaching a steady-state regime under step-like stimulation is well understood, their response under repetitive stimulation, particularly at early time stages is poorly characterized. Yet, early-stage responses to external inputs are arguably as informative as late-stage ones. In simple systems, a periodic stimulation elicits an initial transient response, followed by periodic behavior. Transient responses are relevant when the stimulation has a limited time span, or when the stimulated component's timescale is slow as compared to the timescales of the downstream processes, in which case the latter processes may be capturing only those transients. In this study, we analyze the frequency response of simple motifs at different time stages. We use dose-conserved pulsatile input signals and consider different metrics versus frequency curves. We show that in ligand-receptor systems, there is a frequency preference response in some specific metrics during the transient stages, which is not present in the periodic regime. We suggest this is a general system-level mechanism that cells may use to filter input signals that have consequences for higher order circuits. In addition, we evaluate how the described behavior in isolated motifs is reflected in similar types of responses in cascades and pathways of which they are a part. Our studies suggest that transient frequency preferences are important dynamic features of cell signaling and gene expression systems, which have been overlooked.

Investigation of the relationship between COVID-19 and pancreatic cancer using bioinformatics and systems biology approaches.

BACKGROUND: The coronavirus disease 2019 (COVID-19) pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus, poses a huge threat to human health. Pancreatic cancer (PC) is a malignant tumor with high mortality. Research suggests that infection with SARS-CoV-2 may increase disease severity and risk of death in patients with pancreatic cancer, while pancreatic cancer may also increase the likelihood of contracting SARS-CoV-2, but the link is unclear. METHODS: This study investigated the transcriptional profiles of COVID-19 and PC patients, along with their respective healthy controls, using bioinformatics and systems biology approaches to uncover the molecular mechanisms linking the 2 diseases. Specifically, gene expression data for COVID-19 and PC patients were obtained from the Gene Expression Omnibus datasets, and common differentially expressed genes (DEGs) were identified. Gene ontology and pathway enrichment analyses were performed on the common DEGs to elucidate the regulatory relationships between the diseases. Additionally, hub genes were identified by constructing a protein-protein interaction network from the shared DEGs. Using these hub genes, we conducted regulatory network analyses of microRNA/transcription factors-genes relationships, and predicted potential drugs for treating COVID-19 and PC. RESULTS: A total of 1722 and 2979 DEGs were identified from the transcriptome data of PC (GSE119794) and COVID-19 (GSE196822), respectively. Among these, 236 common DEGs were found between COVID-19 and PC based on protein-protein interaction analysis. Functional enrichment analysis indicated that these shared DEGs were involved in pathways related to viral genome replication and tumorigenesis. Additionally, 10 hub genes, including extra spindle pole bodies like 1, holliday junction recognition protein, marker of proliferation Ki-67, kinesin family member 4A, cyclin-dependent kinase 1, topoisomerase II alpha, cyclin B2, ubiquitin-conjugating enzyme E2 C, aurora kinase B, and targeting protein for Xklp2, were identified. Regulatory network analysis revealed 42 transcription factors and 23 microRNAs as transcriptional regulatory signals. Importantly, lucanthone, etoposide, troglitazone, resveratrol, calcitriol, ciclopirox, dasatinib, enterolactone, methotrexate, and irinotecan emerged as potential therapeutic agents against both COVID-19 and PC. CONCLUSION: This study unveils potential shared pathogenic mechanisms between PC and COVID-19, offering novel insights for future research and therapeutic strategies for the treatment of PC and SARS-CoV-2 infection.

OmicScope unravels systems-level insights from quantitative proteomics data.

Shotgun proteomics analysis presents multifaceted challenges, demanding diverse tool integration for insights. Addressing this complexity, OmicScope emerges as an innovative solution for quantitative proteomics data analysis. Engineered to handle various data formats, it performs data pre-processing - including joining replicates, normalization, data imputation - and conducts differential proteomics analysis for both static and longitudinal experimental designs. Empowered by Enrichr with over 224 databases, OmicScope performs Over Representation Analysis (ORA) and Gene Set Enrichment Analysis (GSEA). Additionally, its Nebula module facilitates meta-analysis from independent datasets, providing a systems biology approach for enriched insights. Complete with a data visualization toolkit and accessible as Python package and a web application, OmicScope democratizes proteomics analysis, offering an efficient and high-quality pipeline for researchers.

Key gene modules and hub genes associated with pyrethroid and organophosphate resistance in Anopheles mosquitoes: a systems biology approach.

Indoor residual spraying (IRS) and insecticide-treated nets (ITNs) are the main methods used to control mosquito populations for malaria prevention. The efficacy of these strategies is threatened by the spread of insecticide resistance (IR), limiting the success of malaria control. Studies of the genetic evolution leading to insecticide resistance could enable the identification of molecular markers that can be used for IR surveillance and an improved understanding of the molecular mechanisms associated with IR. This study used a weighted gene co-expression network analysis (WGCNA) algorithm, a systems biology approach, to identify genes with similar co-expression patterns (modules) and hub genes that are potential molecular markers for insecticide resistance surveillance in Kenya and Benin. A total of 20 and 26 gene co-expression modules were identified via average linkage hierarchical clustering from Anopheles arabiensis and An. gambiae, respectively, and hub genes (highly connected genes) were identified within each module. Three specific genes stood out: serine protease, E3 ubiquitin-protein ligase, and cuticular proteins, which were top hub genes in both species and could serve as potential markers and targets for monitoring IR in these malaria vectors. In addition to the identified markers, we explored molecular mechanisms using enrichment maps that revealed a complex process involving multiple steps, from odorant binding and neuronal signaling to cellular responses, immune modulation, cellular metabolism, and gene regulation. Incorporation of these dynamics into the development of new insecticides and the tracking of insecticide resistance could improve the sustainable and cost-effective deployment of interventions.

Exploring the mechanism of luteolin improving immune and inflammatory responses in systemic sclerosis based on systems biology and cell experiments.

There is a growing trend of applying traditional Chinese medicine (TCM) to treat immune diseases. This study reveals the possible mechanism of luteolin, an active ingredient in the core prescription of TCM, in alleviating systemic sclerosis (SSc) inflammation. Bibliometrics was performed to retrieve the core keywords of SSc inflammation. The key inflammatory indicators in the serum samples of 50 SSc patients were detected by ELISA. Data mining was applied for correlation analysis, association rule analysis, and binary logistic regression analysis on the clinical indicators and medication of 50 SSc patients before and after treatment to determine the core prescription. Network pharmacology was used for identifying candidate genes and pathways; molecular docking was conducted to determine the core monomer components of the prescription, providing a basis for subsequent in vitro molecular mechanism research. The effect of luteolin on SSc-human dermal fibroblasts (HDF) viability and inflammatory factors was evaluated by means of ELISA, RT-PCR, and Western blot. The role of TNF in inflammation was explored by using a TNF overexpression vector, NF-κB inhibitor (PKM2), and SSc-HDF. The involvement of TNF/NF-κB pathway was validated by RT-PCR, Western blot, and immunofluorescence. TCM treatment partially corrected the inflammatory changes in SSc patients, indicating its anti-inflammatory effects in the body. Atractylodes, Yam, Astragalus root, Poria cocos, Pinellia ternata, Salvia miltiorrhiza, Safflower, Cassia twig, and Angelica were identified as the core prescriptions for improving inflammatory indicators. Luteolin was the main active ingredient in the prescription and showed a strong binding energy with TNF and NF-κB. Luteolin exerted anti-inflammatory effects in vitro by reducing inflammatory cytokines in SSc-HDF and inhibiting the activation of TNF/NF-κB. Mechanistically, luteolin inhibited the activation of the TNF/NF-κB pathway in SSc-HDF, as manifested by an increase in extranuclear p-P65 and TNF but a decrease in intranuclear p-P65. Interestingly, the addition of PKM2 augmented the therapeutic function of luteolin against inflammation in SSc-HDF. Our study showed the TCM alleviates the inflammatory response of SSc by inhibiting the activation of the TNF/NF-κB pathway and is an effective therapeutic agent for the treatment of SSc.

TopoDoE: a design of experiment strategy for selection and refinement in ensembles of executable gene regulatory networks.

BACKGROUND: Inference of Gene Regulatory Networks (GRNs) is a difficult and long-standing question in Systems Biology. Numerous approaches have been proposed with the latest methods exploring the richness of single-cell data. One of the current difficulties lies in the fact that many methods of GRN inference do not result in one proposed GRN but in a collection of plausible networks that need to be further refined. In this work, we present a Design of Experiment strategy to use as a second stage after the inference process. It is specifically fitted for identifying the next most informative experiment to perform for deciding between multiple network topologies, in the case where proposed GRNs are executable models. This strategy first performs a topological analysis to reduce the number of perturbations that need to be tested, then predicts the outcome of the retained perturbations by simulation of the GRNs and finally compares predictions with novel experimental data. RESULTS: We apply this method to the results of our divide-and-conquer algorithm called WASABI, adapt its gene expression model to produce perturbations and compare our predictions with experimental results. We show that our networks were able to produce in silico predictions on the outcome of a gene knock-out, which were qualitatively validated for 48 out of 49 genes. Finally, we eliminate as many as two thirds of the candidate networks for which we could identify an incorrect topology, thus greatly improving the accuracy of our predictions. CONCLUSION: These results both confirm the inference accuracy of WASABI and show how executable gene expression models can be leveraged to further refine the topology of inferred GRNs. We hope this strategy will help systems biologists further explore their data and encourage the development of more executable GRN models.

Biomarker Identification for Preterm Birth Susceptibility: Vaginal Microbiome Meta-Analysis Using Systems Biology and Machine Learning Approaches.

PROBLEM: The vaginal microbiome has a substantial role in the occurrence of preterm birth (PTB), which contributes substantially to neonatal mortality worldwide. However, current bioinformatics approaches mostly concentrate on the taxonomic classification and functional profiling of the microbiome, limiting their abilities to elucidate the complex factors that contribute to PTB. METHOD OF STUDY: A total of 3757 vaginal microbiome 16S rRNA samples were obtained from five publicly available datasets. The samples were divided into two categories based on pregnancy outcome: preterm birth (PTB) (N = 966) and term birth (N = 2791). Additionally, the samples were further categorized based on the participants' race and trimester. The 16S rRNA reads were subjected to taxonomic classification and functional profiling using the Parallel-META 3 software in Ubuntu environment. The obtained abundances were analyzed using an integrated systems biology and machine learning approach to determine the key microbes, pathways, and genes that contribute to PTB. The resulting features were further subjected to statistical analysis to identify the top nine features with the greatest effect sizes. RESULTS: We identified nine significant features, namely Shuttleworthia, Megasphaera, Sneathia, proximal tubule bicarbonate reclamation pathway, systemic lupus erythematosus pathway, transcription machinery pathway, lepA gene, pepX gene, and rpoD gene. Their abundance variations were observed through the trimesters. CONCLUSIONS: Vaginal infections caused by Shuttleworthia, Megasphaera, and Sneathia and altered small metabolite biosynthesis pathways such as lipopolysaccharide folate and retinal may increase the susceptibility to PTB. The identified organisms, genes, pathways, and their networks may be specifically targeted for the treatment of bacterial infections that increase PTB risk.

Discovery of a Therapeutic Agent for Glioblastoma Using a Systems Biology-Based Drug Repositioning Approach.

Glioblastoma (GBM), a highly malignant tumour of the central nervous system, presents with a dire prognosis and low survival rates. The heterogeneous and recurrent nature of GBM renders current treatments relatively ineffective. In our study, we utilized an integrative systems biology approach to uncover the molecular mechanisms driving GBM progression and identify viable therapeutic drug targets for developing more effective GBM treatment strategies. Our integrative analysis revealed an elevated expression of CHST2 in GBM tumours, designating it as an unfavourable prognostic gene in GBM, as supported by data from two independent GBM cohorts. Further, we pinpointed WZ-4002 as a potential drug candidate to modulate CHST2 through computational drug repositioning. WZ-4002 directly targeted EGFR (ERBB1) and ERBB2, affecting their dimerization and influencing the activity of adjacent genes, including CHST2. We validated our findings by treating U-138 MG cells with WZ-4002, observing a decrease in CHST2 protein levels and a reduction in cell viability. In summary, our research suggests that the WZ-4002 drug candidate may effectively modulate CHST2 and adjacent genes, offering a promising avenue for developing efficient treatment strategies for GBM patients.

Bioinformatics and systems biology approach to identify the pathogenetic link of neurological pain and major depressive disorder.

Neurological pain (NP) is always accompanied by symptoms of depression, which seriously affects physical and mental health. In this study, we identified the common hub genes (Co-hub genes) and related immune cells of NP and major depressive disorder (MDD) to determine whether they have common pathological and molecular mechanisms. NP and MDD expression data was downloaded from the Gene Expression Omnibus (GEO) database. Common differentially expressed genes (Co-DEGs) for NP and MDD were extracted and the hub genes and hub nodes were mined. Co-DEGs, hub genes, and hub nodes were analyzed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment. Finally, the hub nodes, and genes were analyzed to obtain Co-hub genes. We plotted Receiver operating characteristic (ROC) curves to evaluate the diagnostic impact of the Co-hub genes on MDD and NP. We also identified the immune-infiltrating cell component by ssGSEA and analyzed the relationship. For the GO and KEGG enrichment analyses, 93 Co-DEGs were associated with biological processes (BP), such as fibrinolysis, cell composition (CC), such as tertiary granules, and pathways, such as complement, and coagulation cascades. A differential gene expression analysis revealed significant differences between the Co-hub genes ANGPT2, MMP9, PLAU, and TIMP2. There was some accuracy in the diagnosis of NP based on the expression of ANGPT2 and MMP9. Analysis of differences in the immune cell components indicated an abundance of activated dendritic cells, effector memory CD8+ T cells, memory B cells, and regulatory T cells in both groups, which were statistically significant. In summary, we identified 6 Co-hub genes and 4 immune cell types related to NP and MDD. Further studies are needed to determine the role of these genes and immune cells as potential diagnostic markers or therapeutic targets in NP and MDD.

Design Principles for Perfect Adaptation in Biological Networks with Nonlinear Dynamics.

Establishing a mapping between the emergent biological properties and the repository of network structures has been of great relevance in systems and synthetic biology. Adaptation is one such biological property of paramount importance that promotes regulation in the presence of environmental disturbances. This paper presents a nonlinear systems theory-driven framework to identify the design principles for perfect adaptation with respect to external disturbances of arbitrary magnitude. Based on the prior information about the network, we frame precise mathematical conditions for adaptation using nonlinear systems theory. We first deduce the mathematical conditions for perfect adaptation for constant input disturbances. Subsequently, we translate these conditions to specific necessary structural requirements for adaptation in networks of small size and then extend to argue that there exist only two classes of architectures for a network of any size that can provide local adaptation in the entire state space, namely, incoherent feed-forward (IFF) structure and negative feedback loop with buffer node (NFB). The additional positiveness constraints further narrow the admissible set of network structures. This also aids in establishing the global asymptotic stability for the steady state given a constant input disturbance. The proposed method does not assume any explicit knowledge of the underlying rate kinetics, barring some minimal assumptions. Finally, we also discuss the infeasibility of certain IFF networks in providing adaptation in the presence of downstream connections. Moreover, we propose a generic and novel algorithm based on non-linear systems theory to unravel the design principles for global adaptation. Detailed and extensive simulation studies corroborate the theoretical findings.

Simulated complexes formed from a set of postsynaptic proteins suggest a localised effect of a hypomorphic Shank mutation.

BACKGROUND: The postsynaptic density is an elaborate protein network beneath the postsynaptic membrane involved in the molecular processes underlying learning and memory. The postsynaptic density is built up from the same major proteins but its exact composition and organization differs between synapses. Mutations perturbing protein: protein interactions generally occurring in this network might lead to effects specific for cell types or processes, the understanding of which can be especially challenging. RESULTS: In this work we use systems biology-based modeling of protein complex distributions in a simplified set of major postsynaptic proteins to investigate the effect of a hypomorphic Shank mutation perturbing a single well-defined interaction. We use data sets with widely variable abundances of the constituent proteins. Our results suggest that the effect of the mutation is heavily dependent on the overall availability of all the protein components of the whole network and no trivial correspondence between the expression level of the directly affected proteins and overall complex distribution can be observed. CONCLUSIONS: Our results stress the importance of context-dependent interpretation of mutations. Even the weakening of a generally occurring protein: protein interaction might have well-defined effects, and these can not easily be predicted based only on the abundance of the proteins directly affected. Our results provide insight on how cell-specific effects can be exerted by a mutation perturbing a generally occurring interaction even when the wider interaction network is largely similar.

Systems modeling of oncogenic G-protein and GPCR signaling reveals unexpected differences in downstream pathway activation.

Mathematical models of biochemical reaction networks are an important and emerging tool for the study of cell signaling networks involved in disease processes. One promising potential application of such mathematical models is the study of how disease-causing mutations promote the signaling phenotype that contributes to the disease. It is commonly assumed that one must have a thorough characterization of the network readily available for mathematical modeling to be useful, but we hypothesized that mathematical modeling could be useful when there is incomplete knowledge and that it could be a tool for discovery that opens new areas for further exploration. In the present study, we first develop a mechanistic mathematical model of a G-protein coupled receptor signaling network that is mutated in almost all cases of uveal melanoma and use model-driven explorations to uncover and explore multiple new areas for investigating this disease. Modeling the two major, mutually-exclusive, oncogenic mutations (Gαq/11 and CysLT2R) revealed the potential for previously unknown qualitative differences between seemingly interchangeable disease-promoting mutations, and our experiments confirmed oncogenic CysLT2R was impaired at activating the FAK/YAP/TAZ pathway relative to Gαq/11. This led us to hypothesize that CYSLTR2 mutations in UM must co-occur with other mutations to activate FAK/YAP/TAZ signaling, and our bioinformatic analysis uncovers a role for co-occurring mutations involving the plexin/semaphorin pathway, which has been shown capable of activating this pathway. Overall, this work highlights the power of mechanism-based computational systems biology as a discovery tool that can leverage available information to open new research areas.

Poincaré and SimBio: a versatile and extensible Python ecosystem for modeling systems.

MOTIVATION: Chemical reaction networks (CRNs) play a pivotal role in diverse fields such as systems biology, biochemistry, chemical engineering, and epidemiology. High-level definitions of CRNs enables to use various simulation approaches, including deterministic and stochastic methods, from the same model. However, existing Python tools for simulation of CRN typically wrap external C/C++ libraries for model definition, translation into equations and/or numerically solving them, limiting their extensibility and integration with the broader Python ecosystem. RESULTS: In response, we developed Poincaré and SimBio, two novel Python packages for simulation of dynamical systems and CRNs. Poincaré serves as a foundation for dynamical systems modeling, while SimBio extends this functionality to CRNs, including support for the Systems Biology Markup Language (SBML). Poincaré and SimBio are developed as pure Python packages enabling users to easily extend their simulation capabilities by writing new or leveraging other Python packages. Moreover, this does not compromise the performance, as code can be just-in-time compiled with Numba. Our benchmark tests using curated models from the BioModels repository demonstrate that these tools may provide a potentially superior performance advantage compared to other existing tools. In addition, to ensure a user-friendly experience, our packages use standard typed modern Python syntax that provides a seamless integration with integrated development environments. Our Python-centric approach significantly enhances code analysis, error detection, and refactoring capabilities, positioning Poincaré and SimBio as valuable tools for the modeling community. AVAILABILITY AND IMPLEMENTATION: Poincaré and SimBio are released under the MIT license. Their source code is available on GitHub (https://github.com/maurosilber/poincare and https://github.com/hgrecco/simbio) and can be installed from PyPI or conda-forge.

Advances in Parkinson's disease research - A computational network pharmacological approach.

Parkinson's disease (PD), the second most prevalent neurodegenerative disorder, is projected to see a significant rise in incidence over the next three decades. The precise treatment of PD remains a formidable challenge, prompting ongoing research into early diagnostic methodologies. Network pharmacology, a burgeoning field grounded in systems biology, examines the intricate networks of biological systems to identify critical signal nodes, facilitating the development of multi-target therapeutic molecules. This approach systematically maps the components of Parkinson's disease, thereby reducing its complexity. In this review, we explore the application of network pharmacology workflows in PD, discuss the techniques employed in this field, and evaluate the current advancements and status of network pharmacology in the context of Parkinson's disease. The comprehensive insights will pave newer paths to explore early disease biomarkers and to develop diagnosis with a holistic in silico, in vitro, in vivo and clinical studies.

Reprogramming of Energy Metabolism in Human PKD1 Polycystic Kidney Disease: A Systems Biology Analysis.

Multiple alterations of cellular metabolism have been documented in experimental studies of autosomal dominant polycystic kidney disease (ADPKD) and are thought to contribute to its pathogenesis. To elucidate the molecular pathways and transcriptional regulators associated with the metabolic changes of renal cysts in ADPKD, we compared global gene expression data from human PKD1 renal cysts, minimally cystic tissues (MCT) from the same patients, and healthy human kidney cortical tissue samples. We found gene expression profiles of PKD1 renal cysts were consistent with the Warburg effect with gene pathway changes favoring increased cellular glucose uptake and lactate production, instead of pyruvate oxidation. Additionally, mitochondrial energy metabolism was globally depressed, associated with downregulation of gene pathways related to fatty acid oxidation (FAO), branched-chain amino acid (BCAA) degradation, the Krebs cycle, and oxidative phosphorylation (OXPHOS) in renal cysts. Activation of mTORC1 and its two target proto-oncogenes, HIF-1α and MYC, was predicted to drive the expression of multiple genes involved in the observed metabolic reprogramming (e.g., GLUT3, HK1/HK2, ALDOA, ENO2, PKM, LDHA/LDHB, MCT4, PDHA1, PDK1/3, MPC1/2, CPT2, BCAT1, NAMPT); indeed, their predicted expression patterns were confirmed by our data. Conversely, we found AMPK inhibition was predicted in renal cysts. AMPK inhibition was associated with decreased expression of PGC-1α, a transcriptional coactivator for transcription factors PPARα, ERRα, and ERRγ, all of which play a critical role in regulating oxidative metabolism and mitochondrial biogenesis. These data provide a comprehensive map of metabolic pathway reprogramming in ADPKD and highlight nodes of regulation that may serve as targets for therapeutic intervention.