Abhinav Kaushik

Altered Pathway Analyzer: A gene expression dataset analysis tool for identification and prioritization of differentially regulated and network rewired pathways.

Gene connection rewiring is an essential feature of gene network dynamics. Apart from its normal functional role, it may also lead to dysregulated functional states by disturbing pathway homeostasis. Very few computational tools measure rewiring within gene co-expression and its corresponding regulatory networks in order to identify and prioritize altered pathways which may or may not be differentially regulated. We have developed Altered Pathway Analyzer (APA), a microarray dataset analysis tool for identification and prioritization of altered pathways, including those which are differentially regulated by TFs, by quantifying rewired sub-network topology. Moreover, APA also helps in re-prioritization of APA shortlisted altered pathways enriched with context-specific genes. We performed APA analysis of simulated datasets and p53 status NCI-60 cell line microarray data to demonstrate potential of APA for identification of several case-specific altered pathways. APA analysis reveals several altered pathways not detected by other tools evaluated by us. APA analysis of unrelated prostate cancer datasets identifies sample-specific as well as conserved altered biological processes, mainly associated with lipid metabolism, cellular differentiation and proliferation. APA is designed as a cross platform tool which may be transparently customized to perform pathway analysis in different gene expression datasets. APA is freely available at http://bioinfo.icgeb.res.in/APA.

Proteomic Identification and Analysis of Arginine-Methylated Proteins of Plasmodium falciparum at Asexual Blood Stages

Plasmodium falciparum undergoes a tightly regulated developmental process in human erythrocytes, and recent studies suggest an important regulatory role of post-translational modifications (PTMs). As compared with Plasmodium phosphoproteome, little is known about other PTMs in the parasite. In the present study, we performed a global analysis of asexual blood stages of Plasmodium falciparum to identify arginine-methylated proteins. Using two different methyl arginine-specific antibodies, we immunoprecipitated the arginine-methylated proteins from the stage-specific parasite lysates and identified 843 putative arginine-methylated proteins by LC-MS/MS. Motif analysis of the protein sequences unveiled that the methylation sites are associated with the previously known methylation motifs such as GRx/RGx, RxG, GxxR, or WxxxR. We identified Plasmodium homologues of known arginine-methylated proteins in trypanosomes, yeast, and human. Hydrophilic interaction liquid chromatography (HILIC) was performed on the immunoprecipitates from the trophozoite stage to enrich arginine-methylated peptides. Mass spectrometry analysis of immunoprecipitated and HILIC fractions identified 55 arginine-methylated peptides having 62 methylated arginine sites. Functional classification revealed that the arginine-methylated proteins are involved in RNA metabolism, protein synthesis, intracellular protein trafficking, proteolysis, protein folding, chromatin organization, hemoglobin metabolic process, and several other functions. Summarily, the findings suggest that protein methylation of arginine residues is a widespread phenomenon in Plasmodium, and the PTM may play an important regulatory role in a diverse set of biological pathways, including host-pathogen interactions.

miRMOD: a tool for identification and analysis of 5′ and 3′ miRNA modifications in Next Generation Sequencing small RNA data

In the past decade, the microRNAs (miRNAs) have emerged to be important regulators of gene expression across various species. Several studies have confirmed different types of post-transcriptional modifications at terminal ends of miRNAs. The reports indicate that miRNA modifications are conserved and functionally significant as it may affect miRNA stability and ability to bind mRNA targets, hence affecting target gene repression. Next Generation Sequencing (NGS) of the small RNA (sRNA) provides an efficient and reliable method to explore miRNA modifications. The need for dedicated software, especially for users with little knowledge of computers, to determine and analyze miRNA modifications in sRNA NGS data, motivated us to develop miRMOD. miRMOD is a user-friendly, Microsoft Windows and Graphical User Interface (GUI) based tool for identification and analysis of 5′ and 3′ miRNA modifications (non-templated nucleotide additions and trimming) in sRNA NGS data. In addition to identification of miRNA modifications, the tool also predicts and compares the targets of query and modified miRNAs. In order to compare binding affinities for the same target, miRMOD utilizes minimum free energies of the miRNA:target and modified-miRNA:target interactions. Comparisons of the binding energies may guide experimental exploration of miRNA post-transcriptional modifications. The tool is available as a stand-alone package to overcome large data transfer problems commonly faced in web-based high-throughput (HT) sequencing data analysis tools. miRMOD package is freely available at http://bioinfo.icgeb.res.in/miRMOD.

Widespread occurrence of lysine methylation in Plasmodium falciparum proteins at asexual blood stages

Post-transcriptional and post-translational modifications play a major role in Plasmodium life cycle regulation. Lysine methylation of histone proteins is well documented in several organisms, however in recent years lysine methylation of proteins outside histone code is emerging out as an important post-translational modification (PTM). In the present study we have performed global analysis of lysine methylation of proteins in asexual blood stages of Plasmodium falciparum development. We immunoprecipitated stage specific Plasmodium lysates using anti-methyl lysine specific antibodies that immunostained the asexual blood stage parasites. Using liquid chromatography and tandem mass spectrometry analysis, 570 lysine methylated proteins at three different blood stages were identified. Analysis of the peptide sequences identified 605 methylated sites within 422 proteins. Functional classification of the methylated proteins revealed that the proteins are mainly involved in nucleotide metabolic processes, chromatin organization, transport, homeostatic processes and protein folding. The motif analysis of the methylated lysine peptides reveals novel motifs. Many of the identified lysine methylated proteins are also interacting partners/substrates of PfSET domain proteins as revealed by STRING database analysis. Our findings suggest that the protein methylation at lysine residues is widespread in Plasmodium and plays an important regulatory role in diverse set of the parasite pathways.

Gene Network Rewiring to Study Melanoma Stage Progression and Elements Essential for Driving Melanoma.

Metastatic melanoma patients have a poor prognosis, mainly attributable to the underlying heterogeneity in melanoma driver genes and altered gene expression profiles. These characteristics of melanoma also make the development of drugs and identification of novel drug targets for metastatic melanoma a daunting task. Systems biology offers an alternative approach to re-explore the genes or gene sets that display dysregulated behaviour without being differentially expressed. In this study, we have performed systems biology studies to enhance our knowledge about the conserved property of disease genes or gene sets among mutually exclusive datasets representing melanoma progression. We meta-analysed 642 microarray samples to generate melanoma reconstructed networks representing four different stages of melanoma progression to extract genes with altered molecular circuitry wiring as compared to a normal cellular state. Intriguingly, a majority of the melanoma network-rewired genes are not differentially expressed and the disease genes involved in melanoma progression consistently modulate its activity by rewiring network connections. We found that the shortlisted disease genes in the study show strong and abnormal network connectivity, which enhances with the disease progression. Moreover, the deviated network properties of the disease gene sets allow ranking/prioritization of different enriched, dysregulated and conserved pathway terms in metastatic melanoma, in agreement with previous findings. Our analysis also reveals presence of distinct network hubs in different stages of metastasizing tumor for the same set of pathways in the statistically conserved gene sets. The study results are also presented as a freely available database at http://bioinfo.icgeb.res.in/m3db/. The web-based database resource consists of results from the analysis presented here, integrated with cytoscape web and user-friendly tools for visualization, retrieval and further analysis.

An RNAi-based high-throughput screening assay to identify small molecule inhibitors of hepatitis B virus replication

Persistent or chronic infection with the hepatitis B virus (HBV) represents one of the most common viral diseases in humans. The hepatitis B virus deploys the hepatitis B virus X protein (HBx) as a suppressor of host defenses consisting of RNAi-based silencing of viral genes. Because of its critical role in countering host defenses, HBx represents an attractive target for antiviral drugs. Here, we developed and optimized a loss-of-function screening procedure, which identified a potential pharmacophore that abrogated HBx RNAi suppression activity. In a survey of 14,400 compounds in the Maybridge Screening Collection, we prioritized candidate compounds via high-throughput screening based on reversal of green fluorescent protein (GFP)–reported, RNAi-mediated silencing in a HepG2/GFP-shRNA RNAi sensor line. The screening yielded a pharmacologically active compound, N-(2,4-difluorophenyl)-N′-[3-(1H-imidazol-1-yl) propyl] thiourea (IR415), which blocked HBx-mediated RNAi suppression indicated by the GFP reporter assay. We also found that IR415 reversed the inhibitory effect of HBx protein on activity of the Dicer endoribonuclease. We further confirmed the results of the primary screen in IR415-treated, HBV-infected HepG2 cells, which exhibited a marked depletion of HBV core protein synthesis and down-regulation of pre-genomic HBV RNA. Using a molecular interaction analysis system, we confirmed that IR415 selectively targets HBx in a concentration-dependent manner. The screening assay presented here allows rapid and improved detection of small-molecule inhibitors of HBx and related viral proteins. The assay may therefore potentiate the development of next-generation RNAi pathway-based therapeutics and promises to accelerate our search for novel and effective drugs in antiviral research.

Role of CYP1B1, p.E229K and p.R368H mutations among 120 families with sporadic juvenile onset open-angle glaucoma

BackgroundTo determine the frequency of CYP1B1 p.E229K and p.R368H, gene mutations in a cohort of sporadic juvenile onset open-angle glaucoma (JOAG) patients and to evaluate their genotype/phenotype correlation.MethodsUnrelated JOAG patients whose first-degree relatives had been examined and found to be unaffected were included in the study. The patients and their parents were screened for p.E229K and p.R368H mutations. The phenotypic characteristics were compared between probands carrying the mutations and those who did not carry these mutations.ResultsOut of 120 JOAG patients included in the study, the p.E229K mutation was seen in 9 probands (7.5%) and p.R368H in 7 (5.8%). The average age of onset of the disease (p = 0.3) and the highest untreated IOP (p = 0.4) among those carrying mutations was not significantly different from those who did not have these mutations. The proportion of probands with angle dysgenesis among those with p.E229K and p.R368H mutations was 70% (11 out of 16) in comparison to 65% (67 out of 104) of those who did not harbour these mutations (p = 0.56). Similarly, the probands with moderate to high myopia among those with p.E229K and p.R368H mutations was 20% (3 out of 16) in comparison to 18% (18 out of 104) of those who did not harbour these mutations (p = 0.59).ConclusionThe frequency of p.E229K and p.R368H mutations of the CYP1B1 gene is low even among sporadic JOAG patients. Moreover, there is no clinical correlation between the presence of these mutations and disease severity.

Exploring Heme and Hemoglobin Binding Regions of Plasmodium Heme Detoxification Protein for New Antimalarial Discovery

Hemoglobin degradation/hemozoin formation, essential steps in the Plasmodium life cycle, are targets of existing antimalarials. The pathway still offers vast possibilities to be explored for new antimalarial discoveries. Here, we characterize heme detoxification protein, PfHDP, a major protein involved in hemozoin formation, as a novel drug target. Using in silico and biochemical approaches, we identified two heme binding sites and a hemoglobin binding site in PfHDP. Treatment of Plasmodium falciparum 3D7 parasites with peptide corresponding to the hemoglobin binding domain in PfHDP resulted in food vacuole abnormalities similar to that seen with a cysteine protease inhibitor, E-64 (I-1). Screening of compounds that bound the modeled PfHDP structure in the heme/hemoglobin-binding pockets from Maybridge Screening Collection identified a compound, ML-2, that inhibited parasite growth in a dose-dependent manner, thus paving the way for testing its potential as a new drug candidate. These results provide functional insights into the role of PfHDP in Hz formation and further suggest that PfHDP could be an important drug target to combat malaria.

Protein folding grand challenge: hydrophobic vs. hydrophilic forces

In a recent publication (Ben-Naim, 2012), Arieh Ben-Naim has revisited Protein Folding Problem (PFP) – one of the grand challenges in computational biology. The author supports Levinthal’s conclusion from his work that protein folding doesn’t follow a random path and it is speeded and guided by certain local interactions (Levinthal, 1968). It may be recalled that in his famous work, Levinthal did not mention: what is the nature of these local interactions and whether these local interactions are target or cause-guided. According to the author, these local interactions are the hydrophilic forces originating from solvent molecules. In fact, the author has been advocating the role of hydrophilic forces for a long time, as evident by his several publications on the subject. In this publication, he discusses and dispels the previous landmark theories that “derailed” numerous attempts to solve the PFP, including theories that support target-biased folding or significant role of hydrophobic interactions, before reinforcing his own hypothesis. The author begins the first section of the paper by rejecting the target-based theory put forward by Zwanzig, Szabo, and Bagchi (1992), based on Dawkins’ metaphors, illustrating random events biased towards specific goals. He rejected the theory by stating that there are no protein folding goals marked by correct or incorrect bonds and evolution does not favor correct bonds, and therefore target-based approach cannot answer PFP. His argument is true as there could not be any predefined and specific 3D targets for protein folding, which has to be achieved during protein folding process or even evolution for that matter. He also mentions that nature never had the problem of folding proteins into native structures, rejecting other similar theories, including the one about possible universal code for folding (Kolata, 1986). This is correct as the process of evolution favors the conservation of indispensable proteins, and moreover essentially in conformations (not necessarily corresponding to folded states) that make the proteins functionally active – something that cannot be predefined but happens as a result of evolutionary selection pressures. Next, the author explained how Anfinsen’s thermodynamic hypothesis (Anfinsen, 1973) was misinterpreted by stating that native structure of protein can be obtained by searching for global minimum energy value in GEL or EL, towards which a protein is attracted. He argues that for a stable conformation of a protein, the second law of thermodynamics requires that Gibbs energy reach a single absolute minimum, and GEL in a minima but not essentially a global minimum. Hence, the second law does not support the funnel shape of the GEL and it is misleading to assume that folding process must proceed toward the target of the global GEL minimum. Abandoning the target-based approach, the author advocates cause-based protein folding in the remaining sections of the paper. In fact, the author presents a novel metaphor to build up his theory – the probability of a drunken person reaching a particular house in a big city. By the metaphor, he tries to convince that the probability of the person to reach the same house could be enhanced only if guided by external forces. In following sections, he debates that the protein folding causes to be hydrophilic as opposed to the classical belief that hydrophobic forces are more important in protein folding. He argues that there is no strong evidence in support of the claims that hydrophobic forces play an important role in protein folding. For example, the only evidence reported by Kauzmann (1959) was challenged by him long back (Ben Naim, 1992). In fact, the author out-rightly rejects Dill’s (Dill, Ozkan, Shell, & Weikl, 2008) hypothesis point-by-point. We believe that the author is right in pointing out that there were lacunae in the support of role of hydrophobic forces in protein folding, as Dill described only the role of

A hydrophobic patch surrounding Trp154 in human neuroserpin controls the helix F dynamics with implications in inhibition and aggregation

Neuroserpin (NS) mediated inhibition of tissue-type plasminogen activator (tPA) is important for brain development, synapse formation and memory. Aberrations in helix F and β-sheet A movement during inhibition can directly lead to epilepsy or dementia. Conserved W154 residue in a hydrophobic patch between helix F and β-sheet A is ideally placed to control their movement during inhibition. Molecular Dynamics (MD) simulation on wild type (WT) NS and its two variants (W154A and W154P) demonstrated partial deformation in helix F and conformational differences in strands 1A and 2A only in W154P. A fluorescence and Circular Dichroism (CD) analysis with purified W154 variants revealed a significant red-shift and an increase in α-helical content in W154P as compared to W154A and WT NS. Kinetics of tPA inhibition showed a decline in association rates (ka) for W154A as compared to WT NS with indication of complex formation. Appearance of cleaved without complex formation in W154P indicates that the variant acts as substrate due to conformational misfolding around helix F. Both the variants however showed increased rate of aggregation as compared to WT NS. The hydrophobic patch identified in this study may have importance in helix F dynamics of NS.