System level modeling and analysis of TNF-α mediated sphingolipid signaling pathway in neurological disorders for the prediction of therapeutic targets

Sphingomyelin (SM) belongs to a class of lipids termed sphingolipids. The disruption in the sphingomyelin signaling pathway is associated with various neurodegenerative disorders. TNF-α, a potent pro-inflammatory cytokine generated in response to various neurological disorders like Alzheimer’s disease (AD), Parkinson’s disease (PD), and Multiple Sclerosis (MS), is an eminent regulator of the sphingomyelin metabolic pathway. The immune-triggered regulation of the sphingomyelin metabolic pathway via TNF-α constitutes the sphingomyelin signaling pathway. In this pathway, sphingomyelin and its downstream sphingolipids activate various signaling cascades like PI3K/AKT and MAPK/ERK pathways, thus, controlling diverse processes coupled with neuronal viability, survival, and death. The holistic analysis of the immune-triggered sphingomyelin signaling pathway is imperative to make necessary predictions about its pivotal components and for the formulation of disease-related therapeutics. The current work offers a comprehensive in silico systems analysis of TNF-α mediated sphingomyelin and downstream signaling cascades via a model-based quantitative approach. We incorporated the intensity values of genes from the microarray data of control individuals from the AD study in the input entities of the pathway model. Computational modeling and simulation of the inflammatory pathway enabled the comprehensive study of the system dynamics. Network and sensitivity analysis of the model unveiled essential interaction parameters and entities during neuroinflammation. Scanning of the key entities and parameters allowed us to determine their ultimate impact on neuronal apoptosis and survival. Moreover, the efficacy and potency of the FDA-approved drugs, namely Etanercept, Nivocasan, and Scyphostatin allowed us to study the model’s response towards inhibition of the respective proteins/enzymes. The network analysis revealed the pivotal model entities with high betweenness and closeness centrality values including recruit FADD, TNFR_TRADD, act CASP2, actCASP8, actCASP3 and 9, cytochrome C, and RIP_RAIDD which profoundly impacted the neuronal apoptosis. Whereas some of the entities with high betweenness and closeness centrality values like Gi-coupled receptor, actS1PR, Sphingosine, S1P, actAKT, and actERK produced a high influence on neuronal survival. However, the current study inferred the dual role of ceramide, both on neuronal survival and apoptosis. Moreover, the drug Nivocasan effectively reduces neuronal apoptosis via its inhibitory mechanism on the caspases.