QuanTP: A Software Resource for Quantitative Proteo-Transcriptomic Comparative Data Analysis and Informatics.

Next-generation sequencing technologies, coupled to advances in mass-spectrometry-based proteomics, have facilitated system-wide quantitative profiling of expressed mRNA transcripts and proteins. Proteo-transcriptomic analysis compares the relative abundance levels of transcripts and their corresponding proteins, illuminating discordant gene product responses to perturbations. These results reveal potential post-transcriptional regulation, providing researchers with important new insights into underlying biological and pathological disease mechanisms. To carry out proteo-transcriptomic analysis, researchers require software that statistically determines transcript-protein abundance correlation levels and provides results visualization and interpretation functionality, ideally within a flexible, user-friendly platform. As a solution, we have developed the QuanTP software within the Galaxy platform. The software offers a suite of tools and functionalities critical for proteo-transcriptomics, including statistical algorithms for assessing the correlation between single transcript-protein pairs as well as across two cohorts, outlier identification and clustering, along with a diverse set of results visualizations. It is compatible with analyses of results from single experiment data or from a two-cohort comparison of aggregated replicate experiments. The tool is available in the Galaxy Tool Shed through a cloud-based instance and a Docker container. In all, QuanTP provides an accessible and effective software resource, which should enable new multiomic discoveries from quantitative proteo-transcriptomic data sets.

FusionHub: A unified web platform for annotation and visualization of gene fusion events in human cancer.

Gene fusion is a chromosomal rearrangement event which plays a significant role in cancer due to the oncogenic potential of the chimeric protein generated through fusions. At present many databases are available in public domain which provides detailed information about known gene fusion events and their functional role. Existing gene fusion detection tools, based on analysis of transcriptomics data usually report a large number of fusion genes as potential candidates, which could be either known or novel or false positives. Manual annotation of these putative genes is indeed time-consuming. We have developed a web platform FusionHub, which acts as integrated search engine interfacing various fusion gene databases and simplifies large scale annotation of fusion genes in a seamless way. In addition, FusionHub provides three ways of visualizing fusion events: circular view, domain architecture view and network view. Design of potential siRNA molecules through ensemble method is another utility integrated in FusionHub that could aid in siRNA-based targeted therapy. FusionHub is freely available at https://fusionhub.persistent.co.in.

Structure and function of a highly active Bile Salt Hydrolase (BSH) from Enterococcus faecalis and post-translational processing of BSH enzymes.

Bile Salt Hydrolase (BSH), a member of Cholylglycine hydrolase family, catalyzes the de-conjugation of bile acids and is evolutionarily related to penicillin V acylase (PVA) that hydrolyses a different substrate such as penicillin V. We report the three-dimensional structure of a BSH enzyme from the Gram-positive bacteria Enterococcus faecalis (EfBSH) which has manifold higher hydrolase activity compared to other known BSHs and displays unique allosteric catalytic property. The structural analysis revealed reduced secondary structure content compared to other known BSH structures, particularly devoid of an anti-parallel ß-sheet in the assembly loop and part of a ß-strand is converted to increase the length of a substrate binding loop 2. The analysis of the substrate binding pocket showed reduced volume owing to altered loop conformations and increased hydrophobicity contributed by a higher ratio of hydrophobic to hydrophilic groups present. The aromatic residues F18, Y20 and F65 participate in substrate binding. Thus, their mutation affects enzyme activity. Docking and Molecular Dynamics simulation studies showed effective polar complementarity present for the three hydroxyl (-OH) groups of GCA substrate in the binding site contributing to higher substrate specificity and efficient catalysis. These are unique features characteristics of this BSH enzyme and thought to contribute to its higher activity and specificity towards bile salts as well as allosteric effects. Further, mechanism of autocatalytic processing of Cholylglycine Hydrolases by the excision of an N-terminal Pre-peptide was examined by inserting different N-terminal pre-peptides in EfBSH sequence. The results suggest that two serine residues next to nucleophile cysteine are essential for autocalytic processing to remove precursor peptide. Since pre-peptide is absent in EfBSH the mutation of these serines is tolerated. This suggests that an evolution-mediated subordination of the pre-peptide excision site resulted in loss of pre-peptide in EfBSH and other related Cholylglycine hydrolases.

Elucidation of molecular mechanism of stability of the heme-regulated eIF2α kinase upon binding of its ligand, hemin in its catalytic kinase domain.

The eIF2α kinase activity of the heme-regulated inhibitor (HRI) is regulated by heme which makes it a unique member of the family of eIF2α kinases. Since heme concentrations create an equilibrium for the kinase to be active/inactive, it becomes important to study the heme binding effects upon the kinase and understanding its mechanism of functionality. In the present study, we report the thermostability achieved by the catalytic kinase domain of HRI (HRI.CKD) upon ligand (heme) binding. Our CD data demonstrates that the HRI.CKD retains its secondary structure at higher temperatures when it is in ligand bound state. HRI.CKD when incubated with hemin loses its monomeric state and attains a higher order oligomeric form resulting in its stability. The HRI.CKD fails to refold into its native conformation upon mutation of H377A/H381A, thereby confirming the necessity of these His residues for correct folding, stability, and activity of the kinase. Though our in silico study demonstrated these His being the ligand binding sites in the kinase insert region, the spectra-based study did not show significant difference in heme affinity for the wild type and His mutant HRI.CKD.

Fusarium oxysporum mediates systems metabolic reprogramming of chickpea roots as revealed by a combination of proteomics and metabolomics.

Molecular changes elicited by plants in response to fungal attack and how this affects plant-pathogen interaction, including susceptibility or resistance, remain elusive. We studied the dynamics in root metabolism during compatible and incompatible interactions between chickpea and Fusarium oxysporum f. sp. ciceri (Foc), using quantitative label-free proteomics and NMR-based metabolomics. Results demonstrated differential expression of proteins and metabolites upon Foc inoculations in the resistant plants compared with the susceptible ones. Additionally, expression analysis of candidate genes supported the proteomic and metabolic variations in the chickpea roots upon Foc inoculation. In particular, we found that the resistant plants revealed significant increase in the carbon and nitrogen metabolism; generation of reactive oxygen species (ROS), lignification and phytoalexins. The levels of some of the pathogenesis-related proteins were significantly higher upon Foc inoculation in the resistant plant. Interestingly, results also exhibited the crucial role of altered Yang cycle, which contributed in different methylation reactions and unfolded protein response in the chickpea roots against Foc. Overall, the observed modulations in the metabolic flux as outcome of several orchestrated molecular events are determinant of plant's role in chickpea-Foc interactions.

Structural analysis of a penicillin V acylase from Pectobacterium atrosepticum confirms the importance of two Trp residues for activity and specificity.

Penicillin V acylases (PVA) catalyze the deacylation of the beta-lactam antibiotic phenoxymethylpenicillin (Pen V). They are members of the Ntn hydrolase family and possess an N-terminal cysteine as the main catalytic nucleophile residue. They form the evolutionarily related cholylglycine hydrolase (CGH) group which includes bile salt hydrolases (BSH) responsible for bile deconjugation. Even though a few PVA and BSH structures have been reported, no structure of a functional PVA from Gram-negative bacteria is available. Here, we report the crystal structure of a highly active PVA from Gram-negative Pectobacterium atrosepticum (PaPVA) at 2.5Å resolution. Structural comparison with PVAs from Gram-positive bacteria revealed that PaPVA had a distinctive tetrameric structure and active site organization. In addition, mutagenesis of key active site residues and biochemical characterization of the resultant variants elucidated the role of these residues in substrate binding and catalysis. The importance of residue Trp23 and Trp87 side chains in binding and correct positioning of Pen V by PVAs was confirmed using mutagenesis and substrate docking with a 15ns molecular dynamics simulation. These results establish the unique nature of Gram-negative CGHs and necessitate further research about their substrate spectrum.

Sequence and structure-based comparative analysis to assess, identify and improve the thermostability of penicillin G acylases.

Penicillin acylases are enzymes employed by the pharmaceutical industry for the manufacture of semi-synthetic penicillins. There is a continuous demand for thermostable and alkalophilic enzymes in such applications. We have carried out a computational analysis of known penicillin G acylases (PGAs) in terms of their thermostable nature using various protein-stabilizing factors. While the presence of disulfide bridges was considered initially to screen putative thermostable PGAs from the database, various other factors such as high arginine to lysine ratio, less content of thermolabile amino acids, presence of proline in ß-turns, more number of ion-pair and other non-bonded interactions were also considered for comparison. A modified consensus approach designed could further identify stabilizing residue positions by site-specific comparison between mesostable and thermostable PGAs. A most likely thermostable enzyme identified from the analysis was PGA from Paracoccus denitrificans (PdPGA). This was cloned, expressed and tested for its thermostable nature using biochemical and biophysical experiments. The consensus site-specific sequence-based approach predicted PdPGA to be more thermostable than Escherichia coli PGA, but not as thermostable as the PGA from Achromobacter xylosoxidans. Experimental data showed that PdPGA was comparatively less thermostable than Achromobacter xylosoxidans PGA, although thermostability factors favored a much higher stability. Despite being mesostable, PdPGA being active and stable at alkaline pH is an advantage. Finally, several residue positions could be identified in PdPGA, which upon mutation selectively could improve the thermostability of the enzyme.

Engineering Proteins for Thermostability with iRDP Web Server.

Engineering protein molecules with desired structure and biological functions has been an elusive goal. Development of industrially viable proteins with improved properties such as stability, catalytic activity and altered specificity by modifying the structure of an existing protein has widely been targeted through rational protein engineering. Although a range of factors contributing to thermal stability have been identified and widely researched, the in silico implementation of these as strategies directed towards enhancement of protein stability has not yet been explored extensively. A wide range of structural analysis tools is currently available for in silico protein engineering. However these tools concentrate on only a limited number of factors or individual protein structures, resulting in cumbersome and time-consuming analysis. The iRDP web server presented here provides a unified platform comprising of iCAPS, iStability and iMutants modules. Each module addresses different facets of effective rational engineering of proteins aiming towards enhanced stability. While iCAPS aids in selection of target protein based on factors contributing to structural stability, iStability uniquely offers in silico implementation of known thermostabilization strategies in proteins for identification and stability prediction of potential stabilizing mutation sites. iMutants aims to assess mutants based on changes in local interaction network and degree of residue conservation at the mutation sites. Each module was validated using an extensively diverse dataset. The server is freely accessible at http://irdp.ncl.res.in and has no login requirements.

Structure mediation in substrate binding and post-translational processing of penicillin acylases: Information from mutant structures of Kluyvera citrophila penicillin G acylase.

Penicillin acylases are industrially important enzymes for the production of 6-APA, which is used extensively in the synthesis of secondary antibiotics. The enzyme translates into an inactive single chain precursor that subsequently gets processed by the removal of a spacer peptide connecting the chains of the mature active heterodimer. We have cloned the penicillin G acylase from Kluyvera citrophila (KcPGA) and prepared two mutants by site-directed mutagenesis. Replacement of N-terminal serine of the ß-subunit with cysteine (Serß1Cys) resulted in a fully processed but inactive enzyme. The second mutant in which this serine is replaced by glycine (Serß1Gly) remained in the unprocessed and inactive form. The crystals of both mutants belonged to space group P1 with four molecules in the asymmetric unit. The three-dimensional structures of these mutants were refined at resolutions 2.8 and 2.5 Å, respectively. Comparison of these structures with similar structures of Escherichia coli PGA (EcPGA) revealed various conformational changes that lead to autocatalytic processing and consequent removal of the spacer peptide. The large displacements of residues such as Arg168 and Arg477 toward the N-terminal cleavage site of the spacer peptide or the conformational changes of Arg145 and Phe146 near the active site in these structures suggested probable steps in the processing dynamics. A comparison between the structures of the processed Serß1Cys mutant and that of the processed form of EcPGA showed conformational differences in residues Argα145, Pheα146, and Pheß24 at the substrate binding pocket. Three conformational transitions of Argα145 and Pheα146 residues were seen when processed and unprocessed forms of KcPGA were compared with the substrate bound structure of EcPGA. Structure mediation in activity difference between KcPGA and EcPGA toward acyl homoserine lactone (AHL) is elucidated.

Metabolic profiling of chickpea-Fusarium interaction identifies differential modulation of disease resistance pathways.

Chickpea is the third most widely grown legume in the world and mainly used as a vegetarian source of human dietary protein. Fusarium wilt, caused by Fusarium oxysporum f. sp. ciceri (Foc), is one of the major threats to global chickpea production. Host resistance is the best way to protect crops from diseases; however, in spite of using various approaches, the mechanism of Foc resistance in chickpea remains largely obscure. In the present study, non-targeted metabolic profiling at several time points of resistant and susceptible chickpea cultivars using high-resolution liquid chromatography-mass spectrometry was applied to better understand the mechanistic basis of wilt resistance or susceptibility. Multivariate analysis of the data (OPLS-DA) revealed discriminating metabolites in chickpea root tissue after Foc inoculation such as flavonoids, isoflavonoids, alkaloids, amino acids and sugars. Foc inoculated resistant plants had more flavonoids and isoflavonoids along with their malonyl conjugates. Many antifungal metabolites that were induced after Foc infection viz., aurantion-obstine ß-glucosides and querecitin were elevated in resistant cultivar. Overall, diverse genetic and biochemical mechanisms were operational in the resistant cultivar for Foc defense as compared to the susceptible plant. The resistant chickpea plants employed the above-mentioned metabolic pathways as potential defense strategy against Foc.

An improved method for specificity annotation shows a distinct evolutionary divergence among the microbial enzymes of the cholylglycine hydrolase family.

Bile salt hydrolases (BSHs) are gut microbial enzymes that play a significant role in the bile acid modification pathway. Penicillin V acylases (PVAs) are enzymes produced by environmental microbes, having a possible role in pathogenesis or scavenging of phenolic compounds in their microbial habitats. The correct annotation of such physiologically and industrially important enzymes is thus vital. The current methods relying solely on sequence homology do not always provide accurate annotations for these two members of the cholylglycine hydrolase (CGH) family as BSH/PVA enzymes. Here, we present an improved method [binding site similarity (BSS)-based scoring system] for the correct annotation of the CGH family members as BSH/PVA enzymes, which along with the phylogenetic information incorporates the substrate specificity as well as the binding site information. The BSS scoring system was developed through the analysis of the binding sites and binding modes of the available BSH/PVA structures with substrates glycocholic acid and penicillin V. The 198 sequences in the dataset were then annotated accurately using BSS scores as BSH/PVA enzymes. The dataset presented contained sequences from Gram-positive bacteria, Gram-negative bacteria and archaea. The clustering obtained for the dataset using the method described above showed a clear distinction in annotation of Gram-positive bacteria and Gram-negative bacteria. Based on this clustering and a detailed analysis of the sequences of the CGH family in the dataset, we could infer that the CGH genes might have evolved in accordance with the hypothesis stating the evolution of diderms and archaea from the monoderms.

In-silico analysis of binding site features and substrate selectivity in plant flavonoid-3-O glycosyltransferases (F3GT) through molecular modeling, docking and dynamics simulation studies.

Flavonoids are a class of plant secondary metabolites that act as storage molecules, chemical messengers, as well as participate in homeostasis and defense processes. They possess pharmaceutical properties important for cancer treatment such as antioxidant and anti-tumor activities. The drug-related properties of flavonoids can be improved by glycosylation. The enzymes glycosyltransferases (GTs) glycosylate acceptor molecules in a regiospecific manner with the help of nucleotide sugar donor molecules. Several plant GTs have been characterized and their amino acid sequences determined. However, three-dimensional structures of only a few are reported. Here, phylogenetic analysis using amino acid sequences have identified a group of GTs with the same regiospecific activity. The structures of these closely related GTs were modeled using homologous GT structures. Their substrate binding sites were elaborated by docking flavonoid acceptor and UDP-sugar donor molecules in the modeled structures. Eight regions near the acceptor binding site in the N- and C- terminal domain of GTs have been identified that bind and specifically glycosylate the 3-OH group of acceptor flavonoids. Similarly, a conserved motif in the C-terminal domain is known to bind a sugar donor substrate. In certain GTs, the substitution of a specific glutamine by histidine in this domain changes the preference of sugar from glucose to galactose as a result of changed pattern of interactions. The molecular modeling, docking, and molecular dynamics simulation studies have revealed the chemical and topological features of the binding site and thus provided insights into the basis of acceptor and donor recognition by GTs.

A new role for penicillin acylases: degradation of acyl homoserine lactone quorum sensing signals by Kluyvera citrophila penicillin G acylase.

Use of penicillin acylases for the production of semi-synthetic penicillins is well-known. Escherichia coli penicillin G acylase (EcPGA) has been extensively used for this purpose; however, Kluyvera citrophila penicillin G acylase (KcPGA) is assumed to be a better substitute, owing to its increased resilience to extreme pH conditions and ease of immobilization. In the present article we report a new dimension for the amidase activity of KcPGA by demonstrating its ability to cleave bacterial quorum sensing signal molecules, acyl homoserine lactones (AHL) with acyl chain length of 6-8 with or without oxo-substitution at third carbon position. Initial evidence of AHL degrading capability of KcPGA was obtained using CV026 based bioassay method. Kinetic studies performed at pH 8.0 and 50 °C revealed 3-oxo-C6 HSL to be the best substrate for the enzyme with V(max) and K(m) values of 21.37+0.85 mM/h/mg of protein and 0.1+0.01 mM, respectively. C6 HSL was found to be the second best substrate with V(max) and K(m) value of 10.06+0.27 mM/h/mg of protein and 0.28+0.02 mM, respectively. Molecular modeling and docking studies performed on the active site of the enzyme support these findings by showing the fitting of AHLs perfectly within the hydrophobic pocket of the enzyme active site.

Structural modelling of substrate binding and inhibition in penicillin V acylase from Pectobacterium atrosepticum.

Penicillin V acylases (PVAs) and bile salt hydrolases (BSHs) have considerable sequence and structural similarity; however, they vary significantly in their substrate specificity. We have identified a PVA from a Gram-negative organism, Pectobacterium atrosepticum (PaPVA) that turned out to be a remote homolog of the PVAs and BSHs reported earlier. Even though the active site residues were conserved in PaPVA it showed high specificity towards penV and interestingly the penV acylase activity was inhibited by bile salts. Comparative modelling and docking studies were carried out to understand the structural differences of the binding site that confer this characteristic property. We show that PaPVA exhibits significant differences in structure, which are in contrast to those of known PVAs and such enzymes from Gram-negative bacteria require further investigation.

Dimensionality reduction in computational demarcation of protein tertiary structures.

Predictive classification of major structural families and fold types of proteins is investigated deploying logistic regression. Only five to seven dimensional quantitative feature vector representations of tertiary structures are found adequate. Results for benchmark sample of non-homologous proteins from SCOP database are presented. Importance of this work as compared to homology modeling and best-known quantitative approaches is highlighted.