Identification of potential microRNAs regulating metabolic plasticity and cellular phenotypes in glioblastoma.

MicroRNAs (miRNAs) play important role in regulating cellular metabolism, and are currently being explored in cancer. As metabolic reprogramming in cancer is a major mediator of phenotypic plasticity, understanding miRNA-regulated metabolism will provide opportunities to identify miRNA targets that can regulate oncogenic phenotypes by taking control of cellular metabolism. In the present work, we studied the effect of differentially expressed miRNAs on metabolism, and associated oncogenic phenotypes in glioblastoma (GBM) using patient-derived data. Networks of differentially expressed miRNAs and metabolic genes were created and analyzed to identify important miRNAs that regulate major metabolism in GBM. Graph network-based approaches like network diffusion, backbone extraction, and different centrality measures were used to analyze these networks for identification of potential miRNA targets. Important metabolic processes and cellular phenotypes were annotated to trace the functional responses associated with these miRNA-regulated metabolic genes and associated phenotype networks. miRNA-regulated metabolic gene subnetworks of cellular phenotypes were extracted, and important miRNAs regulating these phenotypes were identified. The most important outcome of the study is the target miRNA combinations predicted for five different oncogenic phenotypes that can be tested experimentally for miRNA-based therapeutic design in GBM. Strategies implemented in the study can be used to generate testable hypotheses in other cancer types as well, and design context-specific miRNA-based therapy for individual patient. Their usability can be further extended to other gene regulatory networks in cancer and other genetic diseases.

Molecular basis of drug resistance in smoothened receptor: An in silico study of protein resistivity and specificity.

Smoothened (SMO) antagonist Vismodegib effectively inhibits the Hedgehog pathway in proliferating cancer cells. In early stage of treatment, Vismodegib exhibited promising outcomes to regress the tumors cells, but ultimately relapsed due to the drug resistive mutations in SMO mostly occurring before (primary mutations G497W) or after (acquired mutations D473H/Y) anti-SMO therapy. This study investigates the unprecedented insights of structural and functional mechanism hindering the binding of Vismodegib with sensitive and resistant mutant variants of SMO (SMOMut ). Along with the basic dynamic understanding of Vismodegib-SMO complexes, network propagation theory based on heat diffusion principles is first time applied here to identify the modules of residues influenced by the individual mutations. The allosteric modulation by GLY497 residue in Vismodegib bound SMO wild-type (SMOWT ) conformation depicts the interconnections of intermediate residues of SMO with the atom of Vismodegib and identify two important motifs (E-X-P-L) and (Q-A-N-V-T-I-G) mediating this allosteric regulation. In this study a novel computational framework based on the heat diffusion principle is also developed, which identify significant residues of allosteric site causing drug resistivity in SMOMut . This framework could also be useful for assessing the potential allosteric sites of different other proteins. Moreover, previously reported novel inhibitor "ZINC12368305," which is proven to make an energetically favorable complex with SMOWT is chosen as a control sample to assess the impact of receptor mutation on its binding and subsequently identify the important factors that govern binding disparity between Vismodegib and ZINC12368305 bound SMOWT/Mut conformations.

Effect of copper contamination on zooplankton epidemics.

Infectious disease and chemical contamination are increasingly becoming vital issues in many ecosystems. However, studies integrating the two are surprisingly rare. Contamination not only affects the inherent host-resource interaction which influences the epidemic process but may also directly affect epidemiological traits via changes in host's behaviour. The fact that heavy metal such as copper is also an essential trace element for organisms, further increase complexity which make predicting the resultant effect of contamination and disease spread difficult. Motivated by this, we model the effect of copper enrichment on a phytoplankton-zooplankton-fungus system. We show that extremely deficient or toxic copper may have a destabilizing effect on the underlying host-resource dynamics due to increased relative energy fluxes as a result of low host mortality due to fish predation. Further, on incorporating disease into the system, we find that the system can become disease-free for an intermediate range of copper concentration whereas it may persist for very less copper enrichment. Also, we predict that there may exist vulnerable regions of copper concentration near the toxic and deficient levels, where the parasite can invade the system for a comparatively lower spore yield. Overall, our results demonstrate that, the effect of contamination may be fundamental to understanding disease progression in community ecology.

Evolutionary Perspectives of Genotype-Phenotype Factors in Leishmania Metabolism.

The sandfly midgut and the human macrophage phagolysosome provide antagonistic metabolic niches for the endoparasite Leishmania to survive and populate. Although these environments fluctuate across developmental stages, the relative changes in both these environments across parasite generations might remain gradual. Such environmental restrictions might endow parasite metabolism with a choice of specific genotypic and phenotypic factors that can constrain enzyme evolution for successful adaptation to the host. With respect to the available cellular information for Leishmania species, for the first time, we measure the relative contribution of eight inter-correlated predictors related to codon usage, GC content, gene expression, gene length, multi-functionality, and flux-coupling potential of an enzyme on the evolutionary rates of singleton metabolic genes and further compare their effects across three Leishmania species. Our analysis reveals that codon adaptation, multi-functionality, and flux-coupling potential of an enzyme are independent contributors of enzyme evolutionary rates, which can together explain a large variation in enzyme evolutionary rates across species. We also hypothesize that a species-specific occurrence of duplicated genes in novel subcellular locations can create new flux routes through certain singleton flux-coupled enzymes, thereby constraining their evolution. A cross-species comparison revealed both common and species-specific genes whose evolutionary divergence was constrained by multiple independent factors. Out of these, previously known pharmacological targets and virulence factors in Leishmania were identified, suggesting their evolutionary reasons for being important survival factors to the parasite. All these results provide a fundamental understanding of the factors underlying adaptive strategies of the parasite, which can be further targeted.

Machine learning prediction of interaction energies in rigid water clusters.

Classical force fields form a computationally efficient avenue for calculating the energetics of large systems. However, due to the constraints of the underlying analytical form, it is sometimes not accurate enough. Quantum mechanical (QM) methods, although accurate, are computationally prohibitive for large systems. In order to circumvent the bottle-neck of interaction energy estimation of large systems, data driven approaches based on machine learning (ML) have been employed in recent years. In most of these studies, the method of choice is artificial neural networks (ANN). In this work, we have shown an alternative ML method, support vector regression (SVR), that provides comparable accuracy with better computational efficiency. We have further used many body expansion (MBE) along with SVR to predict interaction energies in water clusters (decamers). In the case of dimer and trimer interaction energies, the root mean square errors (RMSEs) of the SVR based scheme are 0.12 kcal mol-1 and 0.34 kcal mol-1, respectively. We show that the SVR and MBE based scheme has a RMSE of 2.78% in the estimation of decamer interaction energy against the parent QM method in a computationally efficient way.

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.

Metabolic reprogramming and signalling cross-talks in tumour-immune interaction: a system-level exploration.

Tumour-immune microenvironment (TIME) is pivotal in tumour progression and immunoediting. Within TIME, immune cells undergo metabolic adjustments impacting nutrient supply and the anti-tumour immune response. Metabolic reprogramming emerges as a promising approach to revert the immune response towards a pro-inflammatory state and conquer tumour dominance. This study proposes immunomodulatory mechanisms based on metabolic reprogramming and employs the regulatory flux balance analysis modelling approach, which integrates signalling, metabolism and regulatory processes. For the first time, a comprehensive system-level model is constructed to capture signalling and metabolic cross-talks during tumour-immune interaction and regulatory constraints are incorporated by considering the time lag between them. The model analysis identifies novel features to enhance the immune response while suppressing tumour activity. Particularly, altering the exchange of succinate and oxaloacetate between glioma and macrophage enhances the pro-inflammatory response of immune cells. Inhibition of glutamate uptake in T-cells disrupts the antioxidant mechanism of glioma and reprograms metabolism. Metabolic reprogramming through adenosine monophosphate (AMP)-activated protein kinase (AMPK), coupled with glutamate uptake inhibition, was identified as the most impactful combination to restore T-cell function. A comprehensive understanding of metabolism and gene regulation represents a favourable approach to promote immune cell recovery from tumour dominance.

A realistic host-vector transmission model for describing malaria prevalence pattern.

Malaria continues to be a major public health concern all over the world even after effective control policies have been employed, and considerable understanding of the disease biology have been attained, from both the experimental and modelling perspective. Interactions between different general and local processes, such as dependence on age and immunity of the human host, variations of temperature and rainfall in tropical and sub-tropical areas, and continued presence of asymptomatic infections, regulate the host-vector interactions, and are responsible for the continuing disease prevalence pattern.In this paper, a general mathematical model of malaria transmission is developed considering short and long-term age-dependent immunity of human host and its interaction with pathogen-infected mosquito vector. The model is studied analytically and numerically to understand the role of different parameters related to mosquitoes and humans. To validate the model with a disease prevalence pattern in a particular region, real epidemiological data from the north-eastern part of India was used, and the effect of seasonal variation in mosquito density was modelled based on local climactic data. The model developed based on general features of host-vector interactions, and modified simply incorporating local environmental factors with minimal changes, can successfully explain the disease transmission process in the region. This provides a general approach toward modelling malaria that can be adapted to control future outbreaks of malaria.

Mechano-immunology in microgravity.

Life on Earth has evolved to thrive in the Earth's natural gravitational field; however, as space technology advances, we must revisit and investigate the effects of unnatural conditions on human health, such as gravitational change. Studies have shown that microgravity has a negative impact on various systemic parts of humans, with the effects being more severe in the human immune system. Increasing costs, limited experimental time, and sample handling issues hampered our understanding of this field. To address the existing knowledge gap and provide confidence in modelling the phenomena, in this review, we highlight experimental works in mechano-immunology under microgravity and different computational modelling approaches that can be used to address the existing problems.

Emerging landscape of molecular interaction networks:Opportunities, challenges and prospects.

Network biology finds application in interpreting molecular interaction networks and providing insightful inferences using graph theoretical analysis of biological systems. The integration of computational biomodelling approaches with different hybrid network-based techniques provides additional information about the behaviour of complex systems. With increasing advances in high-throughput technologies in biological research, attempts have been made to incorporate this information into network structures, which has led to a continuous update of network biology approaches over time. The newly minted centrality measures accommodate the details of omics data and regulatory network structure information. The unification of graph network properties with classical mathematical and computational modelling approaches and technologically advanced approaches like machine-learning- and artificial intelligence-based algorithms leverages the potential application of these techniques. These computational advances prove beneficial and serve various applications such as essential gene prediction, identification of drug-disease interaction and gene prioritization. Hence, in this review, we have provided a comprehensive overview of the emerging landscape of molecular interaction networks using graph theoretical approaches. With the aim to provide information on the wide range of applications of network biology approaches in understanding the interaction and regulation of genes, proteins, enzymes and metabolites at different molecular levels, we have reviewed the methods that utilize network topological properties, emerging hybrid network-based approaches and applications that integrate machine learning techniques to analyse molecular interaction networks. Further, we have discussed the applications of these approaches in biomedical research with a note on future prospects.

Theoretical investigation of malaria prevalence in two Indian cities using the response surface method.

BACKGROUND: Elucidation of the relationships between malaria incidence and climatic and non-climatic factors in a region is of utmost importance in understanding the causative factors of disease spread and design of control strategies. Very often malaria prevalence data is restricted to short time scales (months to few years). This demands application of rigorous statistical modelling techniques for analysis and prediction. The monthly malaria prevalence data for three to five years from two cities in southern India, situated in two different climatic zones, are studied to capture their dependence on climatic factors. METHODS: The statistical technique of response surface method (RSM) is applied for the first time to study any epidemiological data. A new step-by-step model reduction technique is proposed to refine the initial model obtained from RSM. This provides a simpler structure and gives better fit. This combined approach is applied to two types of epidemiological data (Slide Positivity Rates values and Total Malaria cases), for two cities in India with varying strengths of disease prevalence and environmental conditions. RESULTS: The study on these data sets reveals that RSM can be used successfully to elucidate the important environmental factors influencing the transmission of the disease by analysing short epidemiological time series. The proposed approach has high predictive ability over relatively long time horizons. CONCLUSIONS: This method promises to provide reliable forecast of malaria incidence across varying environmental conditions, which may help in designing useful control programmes for malaria.

Mathematical models of malaria--a review.

Mathematical models have been used to provide an explicit framework for understanding malaria transmission dynamics in human population for over 100 years. With the disease still thriving and threatening to be a major source of death and disability due to changed environmental and socio-economic conditions, it is necessary to make a critical assessment of the existing models, and study their evolution and efficacy in describing the host-parasite biology. In this article, starting from the basic Ross model, the key mathematical models and their underlying features, based on their specific contributions in the understanding of spread and transmission of malaria have been discussed. The first aim of this article is to develop, starting from the basic models, a hierarchical structure of a range of deterministic models of different levels of complexity. The second objective is to elaborate, using some of the representative mathematical models, the evolution of modelling strategies to describe malaria incidence by including the critical features of host-vector-parasite interactions. Emphasis is more on the evolution of the deterministic differential equation based epidemiological compartment models with a brief discussion on data based statistical models. In this comprehensive survey, the approach has been to summarize the modelling activity in this area so that it helps reach a wider range of researchers working on epidemiology, transmission, and other aspects of malaria. This may facilitate the mathematicians to further develop suitable models in this direction relevant to the present scenario, and help the biologists and public health personnel to adopt better understanding of the modelling strategies to control the disease.

Design of regulation and dynamics in simple biochemical pathways.

Complex regulation of biochemical pathways in a cell is brought about by the interaction of simpler regulatory structures. Among the basic regulatory designs, feedback inhibition of gene expression is the most common motif in gene regulation and a ubiquitous control structure found in nature. In this work, we have studied a common structural feature (delayed feedback) in gene organisation and shown, both theoretically and experimentally, its subtle but important functional role in gene expression kinetics in a negatively auto-regulated system. Using simple deterministic and stochastic models with varying levels of realism, we present detailed theoretical representations of negatively auto-regulated transcriptional circuits with increasing delays in the establishment of feedback of repression. The models of the circuits with and without delay are studied analytically as well as numerically for variation of parameters and delay lengths. The positive invariance, boundedness of the solutions, local and global asymptotic stability of both the systems around the unique positive steady state are studied analytically. Existence of transient temporal dynamics is shown mathematically. Comparison of the two types of model circuits shows that even though the long-term dynamics is stable and not affected by delays in repression, there is interesting variation in the transient dynamical features with increasing delays. Theoretical predictions are validated through experimentally constructed gene circuits of similar designs. This combined theoretical and experimental study helps delineate the opposing effects of delay-induced instability, and the stability-enhancing property of negative feedback in the pathway behaviour, and gives rationale for the abundance of similar designs in real biochemical pathways.

Delineating infection strategies of Leishmania donovani secretory proteins in Human through host-pathogen protein Interactome prediction.

Interactions of Leishmania donovani secretory virulence factors with the host proteins and their interplay during the infection process in humans is poorly studied in Visceral Leishmaniasis. Lack of a holistic study of pathway level de-regulations caused due to these virulence factors leads to a poor understanding of the parasite strategies to subvert the host immune responses, secure its survival inside the host and further the spread of infection to the visceral organs. In this study, we propose a computational workflow to predict host-pathogen protein interactome of L.donovani secretory virulence factors with human proteins combining sequence-based Interolog mapping and structure-based Domain Interaction mapping techniques. We further employ graph theoretical approaches and shortest path methods to analyze the interactome. Our study deciphers the infection paths involving some unique and understudied disease-associated signaling pathways influencing the cellular phenotypic responses in the host. Our statistical analysis based in silico knockout study unveils for the first time UBC, 1433Z and HS90A mediator proteins as potential immunomodulatory candidates through which the virulence factors employ the infection paths. These identified pathways and novel mediator proteins can be effectively used as possible targets to control and modulate the infection process further aiding in the treatment of Visceral Leishmaniasis.

Study of Combinatorial Drug Synergy of Novel Acridone Derivatives With Temozolomide Using in-silico and in-vitro Methods in the Treatment of Drug-Resistant Glioma.

Drug resistance is one of the critical challenges faced in the treatment of Glioma. There are only limited drugs available in the treatment of Glioma and among them Temozolomide (TMZ) has shown some effectiveness in treating Glioma patients, however, the rate of recovery remains poor due to the inability of this drug to act on the drug resistant tumor sub-populations. Hence, in this study three novel Acridone derivative drugs AC2, AC7, and AC26 have been proposed. These molecules when combined with TMZ show major tumor cytotoxicity that is effective in suppressing growth of cancer cells in both drug sensitive and resistant sub-populations of a tumor. In this study a novel mathematical model has been developed to explore the various drug combinations that may be useful for the treatment of resistant Glioma and show that the combinations of TMZ and Acridone derivatives have a synergistic effect. Also, acute toxicity studies of all three acridone derivatives were carried out for 14 days and were found safe for oral administration of 400 mg/kg body weight on albino Wistar rats. Molecular Docking studies of acridone derivatives with P-glycoprotein (P-gp), multiple resistant protein (MRP), and O6-methylguanine-DNA methyltransferase (MGMT) revealed different binding affinities to the transporters contributing to drug resistance. It is observed that while the Acridone derivatives bind with these drug resistance causing proteins, the TMZ can produce its cytotoxicity at a much lower concentration leading to the synergistic effect. The in silico analysis corroborate well with our experimental findings using TMZ resistant (T-98) and drug sensitive (U-87) Glioma cell lines and we propose three novel drug combinations (TMZ with AC2, AC7, and AC26) and dosages that show high synergy, high selectivity and low collateral toxicity for the use in the treatment of drug resistant Glioma, which could be future drugs in the treatment of Glioblastoma.

Differential suitability of reactive oxygen species and the role of glutathione in regulating paradoxical behavior in gliomas: A mathematical perspective.

Manipulative strategies of ROS in cancer are often exhibited as changes in the redox and thiol ratio of the cells. Cellular responses to oxidative insults are generated in response to these changes which are triggered due to the rerouting of the metabolic framework to maintain survival under stress. However, mechanisms of these metabolic re-routing are not clearly understood and remained debatable. In the present work, we have designed a context-based dynamic metabolic model to establish that the coordinated functioning of glutathione peroxidase (GTHP), glutathione oxidoreductase (GTHO) and NADPH oxidase (NOX) is crucial in determining cancerous transformation, specifically in gliomas. Further, we propose that the puzzling duality of ROS (represented by changes in h2o2 in the present model) in exhibiting varying cellular fates can be determined by considering simultaneous changes in nadph/nadp+ and gsh/gssg that occur during the reprogramming of metabolic reactions. This will be helpful in determining the pro-apoptotic or anti-apoptotic fate of gliomas and can be useful in designing effective pro-oxidant and/or anti-oxidant therapeutic approaches against gliomas.

Essential gene prediction using limited gene essentiality information-An integrative semi-supervised machine learning strategy.

Essential gene prediction helps to find minimal genes indispensable for the survival of any organism. Machine learning (ML) algorithms have been useful for the prediction of gene essentiality. However, currently available ML pipelines perform poorly for organisms with limited experimental data. The objective is the development of a new ML pipeline to help in the annotation of essential genes of less explored disease-causing organisms for which minimal experimental data is available. The proposed strategy combines unsupervised feature selection technique, dimension reduction using the Kamada-Kawai algorithm, and semi-supervised ML algorithm employing Laplacian Support Vector Machine (LapSVM) for prediction of essential and non-essential genes from genome-scale metabolic networks using very limited labeled dataset. A novel scoring technique, Semi-Supervised Model Selection Score, equivalent to area under the ROC curve (auROC), has been proposed for the selection of the best model when supervised performance metrics calculation is difficult due to lack of data. The unsupervised feature selection followed by dimension reduction helped to observe a distinct circular pattern in the clustering of essential and non-essential genes. LapSVM then created a curve that dissected this circle for the classification and prediction of essential genes with high accuracy (auROC > 0.85) even with 1% labeled data for model training. After successful validation of this ML pipeline on both Eukaryotes and Prokaryotes that show high accuracy even when the labeled dataset is very limited, this strategy is used for the prediction of essential genes of organisms with inadequate experimentally known data, such as Leishmania sp. Using a graph-based semi-supervised machine learning scheme, a novel integrative approach has been proposed for essential gene prediction that shows universality in application to both Prokaryotes and Eukaryotes with limited labeled data. The essential genes predicted using the pipeline provide an important lead for the prediction of gene essentiality and identification of novel therapeutic targets for antibiotic and vaccine development against disease-causing parasites.

Prediction of survival rate and effect of drugs on cancer patients with somatic mutations of genes: An AI-based approach.

The causal role of somatic mutation and its interrelationship with gene expression profile during tumor development has already been observed, which plays a major role to decide the cancer grades and overall survival. Accurate and robust prediction of tumor grades and patients' overall survival are important for prognosis, risk factors identification and betterment of the treatment strategy, especially for highly lethal tumors, like gliomas. Here, with the help of more accurate and widely used machine learning-based approaches, we propose an integrative computational pipeline that incorporates somatic mutations and gene expression profile for survival and grade prediction of glioma patients and simultaneously relates it to the drugs to be administered. This study gives us a clear understanding that the same drug is not effective for the treatment of same grade of cancer if the gene mutations are different. The alteration in a specific gene plays a very important role in tumor progression and should also be considered for the selection of appropriate drugs. This proposed framework includes all the necessary factors required for enhancement of therapeutic designs and could be useful for clinicians in determining an accurate and personalized treatment strategy for individual patients suffering from different life threatening diseases.

Exploring Notch Pathway to Elucidate Phenotypic Plasticity and Intra-tumor Heterogeneity in Gliomas.

The phenotypic plasticity and self-renewal of adult neural (aNSCs) and glioblastoma stem cells (GSCs) are both known to be governed by active Notch pathway. During development, GSCs can establish differential hierarchy to produce heterogeneous groups of tumor cells belong to different grades, which makes the tumor ecosystem more complex. However, the molecular events regulating these entire processes are unknown hitherto. In this work, based on the mechanistic regulations of Notch pathway activities, a novel computational framework is introduced to inspect the intra-cellular reactions behind the development of normal and tumorigenic cells from aNSCs and GSCs, respectively. The developmental dynamics of aNSCs/GSCs are successfully simulated and molecular activities regulating the phenotypic plasticity and self-renewal processes in normal and tumor cells are identified. A novel scoring parameter "Activity Ratio" score is introduced to find out driver molecules responsible for the phenotypic plasticity and development of different grades of tumor. A new quantitative method is also developed to predict the future risk of Glioblastoma tumor of an individual with appropriate grade by using the transcriptomics profile of that individual as input. Also, a novel technique is introduced to screen and rank the potential drug-targets for suppressing the growth and differentiation of tumor cells.

Delay-induced model for tumor-immune interaction and control of malignant tumor growth.

The paper deals with the qualitative analysis of the solutions of a system of delay differential equations describing the interaction between tumor and immune cells. The asymptotic stability of the possible steady states is showed and the occurrence of limit cycle of the system around the interior equilibrium is proved by the application of Hopf bifurcation theorem by using the delay as a bifurcation parameter. The length of the delay parameter for preserving stability of the system is also estimated, which gives the idea about the mode of action for controlling oscillations in malignant tumor cell growth. The theoretical and numerical outcomes have been supported through experimental results from literatures. This approach gives new insight of modeling tumor-immune interactions and provides significant control strategies to overcome the large oscillations in tumor cells.