Ram Rup Sarkar

Comparison of human cell signaling pathway databases—evolution, drawbacks and challenges

Elucidating the complexities of cell signaling pathways is of immense importance to gain understanding about various biological phenomenon, such as dynamics of gene/protein expression regulation, cell fate determination, embryogenesis and disease progression. The successful completion of human genome project has also helped experimental and theoretical biologists to analyze various important pathways. To advance this study, during the past two decades, systematic collections of pathway data from experimental studies have been compiled and distributed freely by several databases, which also integrate various computational tools for further analysis. Despite significant advancements, there exist several drawbacks and challenges, such as pathway data heterogeneity, annotation, regular update and automated image reconstructions, which motivated us to perform a thorough review on popular and actively functioning 24 cell signaling databases. Based on two major characteristics, pathway information and technical details, freely accessible data from commercial and academic databases are examined to understand their evolution and enrichment. This review not only helps to identify some novel and useful features, which are not yet included in any of the databases but also highlights their current limitations and subsequently propose the reasonable solutions for future database development, which could be useful to the whole scientific community.

Structural and Logical Analysis of a Comprehensive Hedgehog Signaling Pathway to Identify Alternative Drug Targets for Glioma, Colon and Pancreatic Cancer

Hedgehog is an evolutionarily conserved developmental pathway, widely implicated in controlling various cellular responses such as cellular proliferation and stem cell renewal in human and other organisms, through external stimuli. Aberrant activation of this pathway in human adult stem cell line may cause different types of cancers. Hence, targeting this pathway in cancer therapy has become indispensable, but the non availability of detailed molecular interactions, complex regulations by extra- and intra-cellular proteins and cross talks with other pathways pose a serious challenge to get a coherent understanding of this signaling pathway for making therapeutic strategy. This motivated us to perform a computational study of the pathway and to identify probable drug targets. In this work, from available databases and literature, we reconstructed a complete hedgehog pathway which reports the largest number of molecules and interactions to date. Using recently developed computational techniques, we further performed structural and logical analysis of this pathway. In structural analysis, the connectivity and centrality parameters were calculated to identify the important proteins from the network. To capture the regulations of the molecules, we developed a master Boolean model of all the interactions between the proteins and created different cancer scenarios, such as Glioma, Colon and Pancreatic. We performed perturbation analysis on these cancer conditions to identify the important and minimal combinations of proteins that can be used as drug targets. From our study we observed the under expressions of various oncoproteins in Hedgehog pathway while perturbing at a time the combinations of the proteins GLI1, GLI2 and SMO in Glioma; SMO, HFU, ULK3 and RAS in Colon cancer; SMO, HFU, ULK3, RAS and ERK12 in Pancreatic cancer. This reconstructed Hedgehog signaling pathway and the computational analysis for identifying new combinatory drug targets will be useful for future in-vitro and in-vivo analysis to control different cancers.

Comparison of codon usage bias across Leishmania and Trypanosomatids to understand mRNA secondary structure, relative protein abundance and pathway functions.

Understanding the variations in gene organization and its effect on the phenotype across different Leishmania species, and to study differential clinical manifestations of parasite within the host, we performed large scale analysis of codon usage patterns between Leishmania and other known Trypanosomatid species. We present the causes and consequences of codon usage bias in Leishmania genomes with respect to mutational pressure, translational selection and amino acid composition bias. We establish GC bias at wobble position that governs codon usage bias across Leishmania species, rather than amino acid composition bias. We found that, within Leishmania, homogenous codon context coding for less frequent amino acid pairs and codons avoiding formation of folding structures in mRNA are essentially chosen. We predicted putative differences in global expression between genes belonging to specific pathways across Leishmania. This explains the role of evolution in shaping the otherwise conserved genome to demonstrate species-specific function-level differences for efficient survival.

Optimal combinations of control strategies and cost-effective analysis for visceral leishmaniasis disease transmission

Visceral leishmaniasis (VL) is a deadly neglected tropical disease that poses a serious problem in various countries all over the world. Implementation of various intervention strategies fail in controlling the spread of this disease due to issues of parasite drug resistance and resistance of sandfly vectors to insecticide sprays. Due to this, policy makers need to develop novel strategies or resort to a combination of multiple intervention strategies to control the spread of the disease. To address this issue, we propose an extensive SIR-type model for anthroponotic visceral leishmaniasis transmission with seasonal fluctuations modeled in the form of periodic sandfly biting rate. Fitting the model for real data reported in South Sudan, we estimate the model parameters and compare the model predictions with known VL cases. Using optimal control theory, we study the effects of popular control strategies namely, drug-based treatment of symptomatic and PKDL-infected individuals, insecticide treated bednets and spray of insecticides on the dynamics of infected human and vector populations. We propose that the strategies remain ineffective in curbing the disease individually, as opposed to the use of optimal combinations of the mentioned strategies. Testing the model for different optimal combinations while considering periodic seasonal fluctuations, we find that the optimal combination of treatment of individuals and insecticide sprays perform well in controlling the disease for the time period of intervention introduced. Performing a cost-effective analysis we identify that the same strategy also proves to be efficacious and cost-effective. Finally, we suggest that our model would be helpful for policy makers to predict the best intervention strategies for specific time periods and their appropriate implementation for elimination of visceral leishmaniasis.

Revealing the mystery of metabolic adaptations using a genome scale model of Leishmania infantum

Human macrophage phagolysosome and sandfly midgut provide antagonistic ecological niches for Leishmania parasites to survive and proliferate. Parasites optimize their metabolism to utilize the available inadequate resources by adapting to those environments. Lately, a number of metabolomics studies have revived the interest to understand metabolic strategies utilized by the Leishmania parasite for optimal survival within its hosts. For the first time, we propose a reconstructed genome-scale metabolic model for Leishmania infantum JPCM5, the analyses of which not only captures observations reported by metabolomics studies in other Leishmania species but also divulges novel features of the L. infantum metabolome. Our results indicate that Leishmania metabolism is organized in such a way that the parasite can select appropriate alternatives to compensate for limited external substrates. A dynamic non-essential amino acid motif exists within the network that promotes a restricted redistribution of resources to yield required essential metabolites. Further, subcellular compartments regulate this metabolic re-routing by reinforcing the physiological coupling of specific reactions. This unique metabolic organization is robust against accidental errors and provides a wide array of choices for the parasite to achieve optimal survival.

Drug Targets and Biomarker Identification from Computational Study of Human Notch Signaling Pathway

Notch signaling pathway is widely implicated in controlling various cellular functions, cell fate determination, and stem cell renewal in human but aberrant activity in cancer stem cells may cause different types of cancers. Understanding the complexity of this pathway to identify important targets for cancer therapy and to suppress the pathway activity without affecting the normal functions is of utmost importance to clinical and experimental pharmacologists. For developing therapeutic strategy, non availability of detailed molecular interactions, complex regulations and cross talks with other pathways pose a serious challenge to get a coherent understanding of this pathway. This motivated us to reconstruct the largest human cell specific Notch pathway with more number of molecules and interactions available from literatures and databases. To identify probable drug targets and biomarkers for cancer prognosis, we also performed computational study of the pathway using structural and logical analysis and identified important hub proteins, cross talks and feedback mechanisms. The model simulation is validated using reported mRNA expression profile in Glioblastoma cell line and the predictions not only show significant accuracy but also able to identify the undetermined expressions. From our simulation, to identify novel combinations of drug targetable proteins and better substitute for GAMMA SECRETASE inhibition, we proposed two alternative scenarios: partial suppression of Notch target proteins by NICD1 & HIF1A; and complete suppression by NICD1 & MAML, in Glioblastoma cell line. This reconstructed Notch signaling pathway and the computational analysis for identifying new biomarkers and combinatory drug targets will be useful for future in-vitro and in-vivo analysis to control different cancers.

Identification of Th1/Th2 regulatory switch to promote healing response during leishmaniasis: a computational approach

Leishmania devices its survival strategy by suppressing the host’s immune functions. The antigen molecules produced by Leishmania interferes with the host’s cell signaling cascades and consequently changes the protein expression pattern of the antigen-presenting cell (APC). This creates an environment suitable for the switching of the T-cell responses from a healing Th1 response to a non-healing Th2 response that is favorable for the continued survival of the parasite inside the host APC. Using a reconstructed signaling network of the intracellular and intercellular reactions between a Leishmania infected APC and T-cell, we propose a computational model to predict the inhibitory effect of the Leishmania infected APC on the T-cell and to identify the regulators of this Th1-/Th2-switching behavior as observed during Leishmania infection. In this work, we hypothesize that a complete removal of the parasite could only be achieved with a simultaneous up-regulation of the healing Th1 response and stimulation of nitric oxide (NO) production from the APCs, and downregulation of the non-healing Th2 response and thereby propose several unique combinations of protein molecules that could elicit this anti-Leishmania immune response. Our results indicate that TLR3 may play a positive role in eliciting NO synthesis, while TLR2 may be responsible for inhibiting an anti-Leishmania immune response. Also, TLR3 overexpression (in the APC), when combined with SHP2 inhibition (in the T cell), produces an anti-Leishmania response that is better than the conventional IFN-gamma or IL12 treatment. A similar anti-Leishmania response is also obtained in another combination where TLR3 (in APC) is overexpressed, and SHC and MKP (of T cell) are inhibited and activated, respectively. Through our study, we also observe that Leishmania infection may induce an upregulation of IFN-beta production from the APC that may lead to an upregulation of the RAP1 and SOCS3 proteins inside the T cell, the potential inhibitors of MAPK and JAK-STAT signaling pathways, respectively, via the TYK2-mediated pathway. This study not only enhances our knowledge in understanding the Th1/Th2 regulatory switch to promote healing response during leishmaniasis but also helps to identify novel combinations of proteins as potential immunomodulators.

Integrative Approaches to Understand the Mastery in Manipulation of Host Cytokine Networks by Protozoan Parasites with Emphasis on Plasmodium and Leishmania Species

Diseases by protozoan pathogens pose a significant public health concern, particularly in tropical and subtropical countries, where these are responsible for significant morbidity and mortality. Protozoan pathogens tend to establish chronic infections underscoring their competence at subversion of host immune processes, an important component of disease pathogenesis and of their virulence. Modulation of cytokine and chemokine levels, their crosstalks and downstream signaling pathways, and thereby influencing recruitment and activation of immune cells is crucial to immune evasion and subversion. Many protozoans are now known to secrete effector molecules that actively modulate host immune transcriptome and bring about alterations in host epigenome to alter cytokine levels and signaling. The complexity of multi-dimensional events during interaction of hosts and protozoan parasites ranges from microscopic molecular levels to macroscopic ecological and epidemiological levels that includes disrupting metabolic pathways, cell cycle (Toxoplasma and Theileria sp.), respiratory burst, and antigen presentation (Leishmania spp.) to manipulation of signaling hubs. This requires an integrative systems biology approach to combine the knowledge from all these levels to identify the complex mechanisms of protozoan evolution via immune escape during host–parasite coevolution. Considering the diversity of protozoan parasites, in this review, we have focused on Leishmania and Plasmodium infections. Along with the biological understanding, we further elucidate the current efforts in generating, integrating, and modeling of multi-dimensional data to explain the modulation of cytokine networks by these two protozoan parasites to achieve their persistence in host via immune escape during host–parasite coevolution.

Exploring the differences in metabolic behavior of astrocyte and glioblastoma: a flux balance analysis approach

Brain cancers demonstrate a complex metabolic behavior so as to adapt the external hypoxic environment and internal stress generated by reactive oxygen species. To survive in these stringent conditions, glioblastoma cells develop an antagonistic metabolic phenotype as compared to their predecessors, the astrocytes, thereby quenching the resources expected for nourishing the neurons. The complexity and cumulative effect of the large scale metabolic functioning of glioblastoma is mostly unexplored. In this study, we reconstruct a metabolic network comprising of pathways that are known to be deregulated in glioblastoma cells as compared to the astrocytes. The network, consisted of 147 genes encoding for enzymes performing 247 reactions distributed across five distinct model compartments, was then studied using constrained-based modeling approach by recreating the scenarios for astrocytes and glioblastoma, and validated with available experimental evidences. From our analysis, we predict that glycine requirement of the astrocytes are mostly fulfilled by the internal glycine–serine metabolism, whereas glioblastoma cells demand an external uptake of glycine to utilize it for glutathione production. Also, cystine and glucose were identified to be the major contributors to glioblastoma growth. We also proposed an extensive set of single and double lethal reaction knockouts, which were further perturbed to ascertain their role as probable chemotherapeutic targets. These simulation results suggested that, apart from targeting the reactions of central carbon metabolism, knockout of reactions belonging to the glycine–serine metabolism effectively reduce glioblastoma growth. The combinatorial targeting of glycine transporter with any other reaction belonging to glycine–serine metabolism proved lethal to glioblastoma growth.

Temporal protein expression pattern in intracellular signalling cascade during T-cell activation: A computational study

Various T-cell co-receptor molecules and calcium channel CRAC play a pivotal role in the maintenance of cell’s functional responses by regulating the production of effector molecules (mostly cytokines) that aids in immune clearance and also maintaining the cell in a functionally active state. Any defect in these co-receptor signalling pathways may lead to an altered expression pattern of the effector molecules. To study the propagation of such defects with time and their effect on the intracellular protein expression patterns, a comprehensive and largest pathway map of T-cell activation network is reconstructed manually. The entire pathway reactions are then translated using logical equations and simulated using the published time series microarray expression data as inputs. After validating the model, the effect of in silico knock down of co-receptor molecules on the expression patterns of their downstream proteins is studied and simultaneously the changes in the phenotypic behaviours of the T-cell population are predicted, which shows significant variations among the proteins expression and the signalling routes through which the response is propagated in the cytoplasm. This integrative computational approach serves as a valuable technique to study the changes in protein expression patterns and helps to predict variations in the cellular behaviour.