Application of machine learning on understanding biomolecule interactions in cellular machinery

Machine learning (ML) applications have become ubiquitous in all fields of research including protein science and engineering. Apart from protein structure and mutation prediction, scientists are focusing on knowledge gaps with respect to the molecular mechanisms involved in protein binding and interactions with other components in the experimental setups or the human body. Researchers are working on several wet-lab techniques and generating data for a better understanding of concepts and mechanics involved. The information like biomolecular structure, binding affinities, structure fluctuations and movements are enormous which can be handled and analyzed by ML. Therefore, this review highlights the significance of ML in understanding the biomolecular interactions while assisting in various fields of research such as drug discovery, nanomedicine, nanotoxicity and material science. Hence, the way ahead would be to force hand-in hand of laboratory work and computational techniques.

In-silico Studies Calculated a New Chitin Oligomer Binding Site Inside Vicilin: A Potent Antifungal and Insecticidal Agent

Vicilins are major seed storage proteins and show differential binding affinities toward sugar moieties of fungal cell wall and insect gut epithelium. Hence, purpose of study is the thorough in-silico characterization of interactions between vicilin and chitin oligomer followed by fungal and insecticidal bioassays. This work covers the molecular simulation studies explaining the interactions between Pisum sativum vicilin (PsV) and chitin oligomer followed by protein bioassay against different pathogens. LC-MS/MS of purified PsV (∼50 kDa) generated residual data along high pea vicilin homology (UniProtKB ID; P13918). Predicted model (PsV) indicated the characteristic homotrimer joined through head-to-tail association and each monomer is containing a bicupin domain. PsV site map analysis showed a new site (Site 4) into which molecular docking confirmed the strong binding of chitin oligomer (GlcNAc)₄. Molecular dynamics simulation data (50 ns) indicated that chitin-binding site was comprised of 8 residues (DKEDRNEN). However, aspartate and glutamate significantly contributed in the stability of ligand binding. Computational findings were further verified via significant growth inhibition of Aspergillus flavus, A. niger, and Fusarium oxysporum against PsV. Additionally, the substantial adult population of Brevicoryne brassicae was reduced and different life stages of Tribolium castaneum also showed significant mortality.

Development of nano affinity columns for the study of ligand (including SARS-CoV-2 related proteins) binding to heparan sulfate proteoglycans

The interactions of heparan sulfate proteoglycans (HSPGs) present on the cell surface with target proteins lead to cell signaling and they are considered as viral receptors. The analysis of the recognition mechanism between HSPG and its potential ligands and high-throughput screening in drug discovery thus remain important challenges. Glycidyl methacrylate-based monoliths were thus prepared in situ in miniaturized capillary columns (internal diameter 75 μm) and HSPG was grafted onto them by the use of the Schiff base method. The quantity of grafted HSPG was in the nanogram range (11 nanograms per cm of capillary length). This is of significant importance when working with less available or expensive biological material. Other advantages of our miniaturized capillary column are as follows: (i) the immobilization process of HSPG onto the organic monolithic support was reliable and reproducible. (ii) The resultant affinity capillary column showed a strong resistance to changes in temperature and pH and a negligible non-specific interaction. So as to confirm the proper functioning of our miniaturized capillary column, the molecular recognition by HSPG of five selected compounds including three ligands of interest related to SARS-CoV-2 was studied.

Molecular representations in AI-driven drug discovery: a review and practical guide

The technological advances of the past century, marked by the computer revolution and the advent of high-throughput screening technologies in drug discovery, opened the path to the computational analysis and visualization of bioactive molecules. For this purpose, it became necessary to represent molecules in a syntax that would be readable by computers and understandable by scientists of various fields. A large number of chemical representations have been developed over the years, their numerosity being due to the fast development of computers and the complexity of producing a representation that encompasses all structural and chemical characteristics. We present here some of the most popular electronic molecular and macromolecular representations used in drug discovery, many of which are based on graph representations. Furthermore, we describe applications of these representations in AI-driven drug discovery. Our aim is to provide a brief guide on structural representations that are essential to the practice of AI in drug discovery. This review serves as a guide for researchers who have little experience with the handling of chemical representations and plan to work on applications at the interface of these fields.

Structural Analysis of Binding Determinants of Salmonella typhimurium Trehalose-6-phosphate Phosphatase Using Ground-State Complexes

Trehalose-6-phosphate phosphatase (T6PP) catalyzes the dephosphorylation of trehalose 6-phosphate (T6P) to the disaccharide trehalose. The enzyme is not present in mammals but is essential to the viability of multiple lower organisms as trehalose is a critical metabolite, and T6P accumulation is toxic. Hence, T6PP is a target for therapeutics of human pathologies caused by bacteria, fungi, and parasitic nematodes. Here, we report the X-ray crystal structures of Salmonella typhimurium T6PP (StT6PP) in its apo form and in complex with the cofactor Mg²⁺ and the substrate analogue trehalose 6-sulfate (T6S), the product trehalose, or the competitive inhibitor 4-n-octylphenyl α-d-glucopyranoside 6-sulfate (OGS). OGS replaces the substrate phosphoryl group with a sulfate group and the glucosyl ring distal to the sulfate group with an octylphenyl moiety. The structures of these substrate-analogue and product complexes with T6PP show that specificity is conferred via hydrogen bonds to the glucosyl group proximal to the phosphoryl moiety through Glu123, Lys125, and Glu167, conserved in T6PPs from multiple species. The structure of the first-generation inhibitor OGS shows that it retains the substrate-binding interactions observed for the sulfate group and the proximal glucosyl ring. The OGS octylphenyl moiety binds in a unique manner, indicating that this subsite can tolerate various chemotypes. Together, these findings show that these conserved interactions at the proximal glucosyl ring binding site could provide the basis for the development of broad-spectrum therapeutics, whereas variable interactions at the divergent distal subsite could present an opportunity for the design of potent organism-specific therapeutics.

Chemical Modifications of Vicilins Interfere with Chitin-Binding Affinity and Toxicity to Callosobruchus maculatus (Coleoptera: Chrysomelidae) Insect: A Combined In Vitro and In Silico Analysis

Vicilins are related to cowpea seed resistance toward Callosobruchus maculatus due to their ability to bind to chitinous structures lining larval midgut. However, this binding mechanism is not fully understood. Here, we identified chitin binding sites and investigated how in vitro and in silico chemical modifications interfere with vicilin chitin binding and insect toxicity. In vitro assays showed that unmodified vicilin strongly binds to chitin matrices, mainly with acetylated chitin. Chemical modifications of specific amino acids (tryptophan, lysine, tyrosine), as well as glutaraldehyde cross-linking, decreased the evaluated parameters. In silico analyses identified at least one chitin binding site in vicilin monomer, the region between Arg208 and Lys216, which bears the sequence REGIRELMK and forms an α helix, exposed in the 3D structure. In silico modifications of Lys223 (acetylated at its terminal nitrogen) and Trp316 (iodinated to 7-iodine-L-tryptophan or oxidized to β-oxy-indolylalanine) decreased vicilin chitin binding affinity. Glucose, sucrose, and N-acetylglucosamine also interfered with vicilin chitin binding affinity.

The current structural glycome landscape and emerging technologies

Carbohydrates represent one of the building blocks of life, along with nucleic acids, proteins and lipids. Although glycans are involved in a wide range of processes from embryogenesis to protein trafficking and pathogen infection, we are still a long way from deciphering the glycocode. In this review, we aim to present a few of the challenges that researchers working in the area of glycobiology can encounter and what strategies can be utilised to overcome them. Our goal is to paint a comprehensive picture of the current saccharide landscape available in the Protein Data Bank (PDB). We also review recently updated repositories relevant to the topic proposed, the impact of software development on strategies to structurally solve carbohydrate moieties, and state-of-the-art molecular and cellular biology methods that can shed some light on the function and structure of glycans.

Glucose-induced structural changes and anomalous diffusion of elastin

Elastin is the principal protein component of elastic fiber, which renders essential elasticity to connective tissues and organs. Here, we adopted a multi-technique approach to study the transport, viscoelastic, and structural properties of elastin exposed to various glucose concentrations (X=[gluc]/[elastin]). Laser light scattering experiments revealed an anomalous behavior (anomaly exponent, β <0.6) of elastin. In this regime (β <0.6), the diffusion constant decreases by 40% in the presence of glucose (X> 10), which suggests the structural change in elastin. We have observed a peculiar inverse temperature transition of elastin protein, which is a measure of structural change, at 40 °C through rheology experiments. Moreover, we observe its shift towards lower temperature with a higher X. FTIR revealed that the presence of glucose (X < 10) favors the formation of β-sheet structure in elastin. However, for X > 10, dominative crowding effect reduces the mobility of protein and favors the increase in β-turns and γ-turns by 25 ± 1% over the β-sheet (β-sheet decreases by 12 ± 0.8%) and α-helix (α-helix decreases by 13 ± 0.8%). The stiffness of protein is estimated through Flory characteristic ratio, C_∞ and found to be increasing with X. These glucose-based structural changes in the elastin may explain the role of glucose in age-related issues of the skin.

Computational Prediction of Carbohydrate-Binding Proteins and Binding Sites

Protein-carbohydrate interaction is essential for biological systems, and carbohydrate-binding proteins (CBPs) are important targets when designing antiviral and anticancer drugs. Due to the high cost and difficulty associated with experimental approaches, many computational methods have been developed as complementary approaches to predict CBPs or carbohydrate-binding sites. However, most of these computational methods are not publicly available. Here, we provide a comprehensive review of related studies and demonstrate our two recently developed bioinformatics methods. The method SPOT-CBP is a template-based method for detecting CBPs based on structure through structural homology search combined with a knowledge-based scoring function. This method can yield model complex structure in addition to accurate prediction of CBPs. Furthermore, it has been observed that similarly accurate predictions can be made using structures from homology modeling, which has significantly expanded its applicability. The other method, SPRINT-CBH, is a de novo approach that predicts binding residues directly from protein sequences by using sequence information and predicted structural properties. This approach does not need structurally similar templates and thus is not limited by the current database of known protein-carbohydrate complex structures. These two complementary methods are available at https://sparks-lab.org. © 2018 by John Wiley & Sons, Inc.

Development of a sugar-binding residue prediction system from protein sequences using support vector machine

Several methods have been proposed for protein-sugar binding site prediction using machine learning algorithms. However, they are not effective to learn various properties of binding site residues caused by various interactions between proteins and sugars. In this study, we classified sugars into acidic and nonacidic sugars and showed that their binding sites have different amino acid occurrence frequencies. By using this result, we developed sugar-binding residue predictors dedicated to the two classes of sugars: an acid sugar binding predictor and a nonacidic sugar binding predictor. We also developed a combination predictor which combines the results of the two predictors. We showed that when a sugar is known to be an acidic sugar, the acidic sugar binding predictor achieves the best performance, and showed that when a sugar is known to be a nonacidic sugar or is not known to be either of the two classes, the combination predictor achieves the best performance. Our method uses only amino acid sequences for prediction. Support vector machine was used as a machine learning algorithm and the position-specific scoring matrix created by the position-specific iterative basic local alignment search tool was used as the feature vector. We evaluated the performance of the predictors using five-fold cross-validation. We have launched our system, as an open source freeware tool on the GitHub repository (https://doi.org/10.5281/zenodo.61513).

Sequence-Based Prediction of Protein-Carbohydrate Binding Sites Using Support Vector Machines

Carbohydrate-binding proteins play significant roles in many diseases including cancer. Here, we established a machine-learning-based method (called sequence-based prediction of residue-level interaction sites of carbohydrates, SPRINT-CBH) to predict carbohydrate-binding sites in proteins using support vector machines (SVMs). We found that integrating evolution-derived sequence profiles with additional information on sequence and predicted solvent accessible surface area leads to a reasonably accurate, robust, and predictive method, with area under receiver operating characteristic curve (AUC) of 0.78 and 0.77 and Matthew's correlation coefficient of 0.34 and 0.29, respectively for 10-fold cross validation and independent test without balancing binding and nonbinding residues. The quality of the method is further demonstrated by having statistically significantly more binding residues predicted for carbohydrate-binding proteins than presumptive nonbinding proteins in the human proteome, and by the bias of rare alleles toward predicted carbohydrate-binding sites for nonsynonymous mutations from the 1000 genome project. SPRINT-CBH is available as an online server at http://sparks-lab.org/server/SPRINT-CBH .

Identification of Glucose-Binding Pockets in Human Serum Albumin Using Support Vector Machine and Molecular Dynamics Simulations

Human Serum Albumin (HSA) has been suggested to be an alternate biomarker to the existing Hemoglobin-A1c (HbA1c) marker for glycemic monitoring. Development and usage of HSA as an alternate biomarker requires the identification of glycation sites, or equivalently, glucose-binding pockets. In this work, we combine molecular dynamics simulations of HSA and the state-of-art machine learning method Support Vector Machine (SVM) to predict glucose-binding pockets in HSA. SVM uses the three dimensional arrangement of atoms and their chemical properties to predict glucose-binding ability of a pocket. Feature selection reveals that the arrangement of atoms and their chemical properties within the first 4Å from the centroid of the pocket play an important role in the binding of glucose. With a 10-fold cross validation accuracy of 84 percent, our SVM model reveals seven new potential glucose-binding sites in HSA of which two are exposed only during the dynamics of HSA. The predictions are further corroborated using docking studies. These findings can complement studies directed towards the development of HSA as an alternate biomarker for glycemic monitoring.

The evolution of drug discovery: from phenotypes to targets, and back

Cumulative scientific and technological advances over the past two centuries have transformed drug discovery from a largely serendipitous process into the high tech pipelines of today.

MOWGLI: prediction of protein-MannOse interacting residues With ensemble classifiers usinG evoLutionary Information

Proteins interact with carbohydrates to perform various cellular interactions. Of the many carbohydrate ligands that proteins bind with, mannose constitute an important class, playing important roles in host defense mechanisms. Accurate identification of mannose-interacting residues (MIR) may provide important clues to decipher the underlying mechanisms of protein-mannose interactions during infections. This study proposes an approach using an ensemble of base classifiers for prediction of MIR using their evolutionary information in the form of position-specific scoring matrix. The base classifiers are random forests trained by different subsets of training data set Dset128 using 10-fold cross-validation. The optimized ensemble of base classifiers, MOWGLI, is then used to predict MIR on protein chains of the test data set Dtestset29 which showed a promising performance with 92.0% accurate prediction. An overall improvement of 26.6% in precision was observed upon comparison with the state-of-art. It is hoped that this approach, yielding enhanced predictions, could be eventually used for applications in drug design and vaccine development.

Carbohydrate-binding protein identification by coupling structural similarity searching with binding affinity prediction

Carbohydrate-binding proteins (CBPs) are potential biomarkers and drug targets. However, the interactions between carbohydrates and proteins are challenging to study experimentally and computationally because of their low binding affinity, high flexibility, and the lack of a linear sequence in carbohydrates as exists in RNA, DNA, and proteins. Here, we describe a structure-based function-prediction technique called SPOT-Struc that identifies carbohydrate-recognizing proteins and their binding amino acid residues by structural alignment program SPalign and binding affinity scoring according to a knowledge-based statistical potential based on the distance-scaled finite-ideal gas reference state (DFIRE). The leave-one-out cross-validation of the method on 113 carbohydrate-binding domains and 3442 noncarbohydrate binding proteins yields a Matthews correlation coefficient of 0.56 for SPalign alone and 0.63 for SPOT-Struc (SPalign + binding affinity scoring) for CBP prediction. SPOT-Struc is a technique with high positive predictive value (79% correct predictions in all positive CBP predictions) with a reasonable sensitivity (52% positive predictions in all CBPs). The sensitivity of the method was changed slightly when applied to 31 APO (unbound) structures found in the protein databank (14/31 for APO versus 15/31 for HOLO). The result of SPOT-Struc will not change significantly if highly homologous templates were used. SPOT-Struc predicted 19 out of 2076 structural genome targets as CBPs. In particular, one uncharacterized protein in Bacillus subtilis (1oq1A) was matched to galectin-9 from Mus musculus. Thus, SPOT-Struc is useful for uncovering novel carbohydrate-binding proteins. SPOT-Struc is available at http://sparks-lab.org.

Community-based network study of protein-carbohydrate interactions in plant lectins using glycan array data

Lectins play major roles in biological processes such as immune recognition and regulation, inflammatory responses, cytokine signaling, and cell adhesion. Recently, glycan microarrays have shown to play key roles in understanding glycobiology, allowing us to study the relationship between the specificities of glycan binding proteins and their natural ligands at the omics scale. However, one of the drawbacks in utilizing glycan microarray data is the lack of systematic analysis tools to extract information. In this work, we attempt to group various lectins and their interacting carbohydrates by using community-based analysis of a lectin-carbohydrate network. The network consists of 1119 nodes and 16769 edges and we have identified 3 lectins having large degrees of connectivity playing the roles of hubs. The community based network analysis provides an easy way to obtain a general picture of the lectin-glycan interaction and many statistically significant functional groups.

Prediction of vitamin interacting residues in a vitamin binding protein using evolutionary information

BACKGROUND

The vitamins are important cofactors in various enzymatic-reactions. In past, many inhibitors have been designed against vitamin binding pockets in order to inhibit vitamin-protein interactions. Thus, it is important to identify vitamin interacting residues in a protein. It is possible to detect vitamin-binding pockets on a protein, if its tertiary structure is known. Unfortunately tertiary structures of limited proteins are available. Therefore, it is important to develop in-silico models for predicting vitamin interacting residues in protein from its primary structure.

RESULTS

In this study, first we compared protein-interacting residues of vitamins with other ligands using Two Sample Logo (TSL). It was observed that ATP, GTP, NAD, FAD and mannose preferred {G,R,K,S,H}, {G,K,T,S,D,N}, {T,G,Y}, {G,Y,W} and {Y,D,W,N,E} residues respectively, whereas vitamins preferred {Y,F,S,W,T,G,H} residues for the interaction with proteins. Furthermore, compositional information of preferred and non-preferred residues along with patterns-specificity was also observed within different vitamin-classes. Vitamins A, B and B6 preferred {F,I,W,Y,L,V}, {S,Y,G,T,H,W,N,E} and {S,T,G,H,Y,N} interacting residues respectively. It suggested that protein-binding patterns of vitamins are different from other ligands, and motivated us to develop separate predictor for vitamins and their sub-classes. The four different prediction modules, (i) vitamin interacting residues (VIRs), (ii) vitamin-A interacting residues (VAIRs), (iii) vitamin-B interacting residues (VBIRs) and (iv) pyridoxal-5-phosphate (vitamin B6) interacting residues (PLPIRs) have been developed. We applied various classifiers of SVM, BayesNet, NaiveBayes, ComplementNaiveBayes, NaiveBayesMultinomial, RandomForest and IBk etc., as machine learning techniques, using binary and Position-Specific Scoring Matrix (PSSM) features of protein sequences. Finally, we selected best performing SVM modules and obtained highest MCC of 0.53, 0.48, 0.61, 0.81 for VIRs, VAIRs, VBIRs, PLPIRs respectively, using PSSM-based evolutionary information. All the modules developed in this study have been trained and tested on non-redundant datasets and evaluated using five-fold cross-validation technique. The performances were also evaluated on the balanced and different independent datasets.

CONCLUSIONS

This study demonstrates that it is possible to predict VIRs, VAIRs, VBIRs and PLPIRs from evolutionary information of protein sequence. In order to provide service to the scientific community, we have developed web-server and standalone software VitaPred (http://crdd.osdd.net/raghava/vitapred/).

Prediction of protein-mannose binding sites using random forest

Mannose is an abundant cell surface monosaccharide and has an important role in many biochemical processes. It binds to a great diversity of receptor proteins. In this study we have employed Random Forest for prediction of mannose binding sites. Mannosebinding site is taken to be a sphere around the centroid of the ligand and the sphere is subdivided into different layers and atom wise and residue wise features were extracted for each layer. The method achieves 95.59 % of accuracy using Random Forest with 10 fold cross validation. Prediction of mannose binding site analysis will be quite useful in drug design.

Local functional descriptors for surface comparison based binding prediction

Background

Molecular recognition in proteins occurs due to appropriate arrangements of physical, chemical, and geometric properties of an atomic surface. Similar surface regions should create similar binding interfaces. Effective methods for comparing surface regions can be used in identifying similar regions, and to predict interactions without regard to the underlying structural scaffold that creates the surface.

Results

We present a new descriptor for protein functional surfaces and algorithms for using these descriptors to compare protein surface regions to identify ligand binding interfaces. Our approach uses descriptors of local regions of the surface, and assembles collections of matches to compare larger regions. Our approach uses a variety of physical, chemical, and geometric properties, adaptively weighting these properties as appropriate for different regions of the interface. Our approach builds a classifier based on a training corpus of examples of binding sites of the target ligand. The constructed classifiers can be applied to a query protein providing a probability for each position on the protein that the position is part of a binding interface. We demonstrate the effectiveness of the approach on a number of benchmarks, demonstrating performance that is comparable to the state-of-the-art, with an approach with more generality than these prior methods.

Conclusions

Local functional descriptors offer a new method for protein surface comparison that is sufficiently flexible to serve in a variety of applications.

Prediction of carbohydrate binding sites on protein surfaces with 3-dimensional probability density distributions of interacting atoms

Non-covalent protein-carbohydrate interactions mediate molecular targeting in many biological processes. Prediction of non-covalent carbohydrate binding sites on protein surfaces not only provides insights into the functions of the query proteins; information on key carbohydrate-binding residues could suggest site-directed mutagenesis experiments, design therapeutics targeting carbohydrate-binding proteins, and provide guidance in engineering protein-carbohydrate interactions. In this work, we show that non-covalent carbohydrate binding sites on protein surfaces can be predicted with relatively high accuracy when the query protein structures are known. The prediction capabilities were based on a novel encoding scheme of the three-dimensional probability density maps describing the distributions of 36 non-covalent interacting atom types around protein surfaces. One machine learning model was trained for each of the 30 protein atom types. The machine learning algorithms predicted tentative carbohydrate binding sites on query proteins by recognizing the characteristic interacting atom distribution patterns specific for carbohydrate binding sites from known protein structures. The prediction results for all protein atom types were integrated into surface patches as tentative carbohydrate binding sites based on normalized prediction confidence level. The prediction capabilities of the predictors were benchmarked by a 10-fold cross validation on 497 non-redundant proteins with known carbohydrate binding sites. The predictors were further tested on an independent test set with 108 proteins. The residue-based Matthews correlation coefficient (MCC) for the independent test was 0.45, with prediction precision and sensitivity (or recall) of 0.45 and 0.49 respectively. In addition, 111 unbound carbohydrate-binding protein structures for which the structures were determined in the absence of the carbohydrate ligands were predicted with the trained predictors. The overall prediction MCC was 0.49. Independent tests on anti-carbohydrate antibodies showed that the carbohydrate antigen binding sites were predicted with comparable accuracy. These results demonstrate that the predictors are among the best in carbohydrate binding site predictions to date.

Automated identification of protein-ligand interaction features using Inductive Logic Programming: a hexose binding case study

BACKGROUND

There is a need for automated methods to learn general features of the interactions of a ligand class with its diverse set of protein receptors. An appropriate machine learning approach is Inductive Logic Programming (ILP), which automatically generates comprehensible rules in addition to prediction. The development of ILP systems which can learn rules of the complexity required for studies on protein structure remains a challenge. In this work we use a new ILP system, ProGolem, and demonstrate its performance on learning features of hexose-protein interactions.

RESULTS

The rules induced by ProGolem detect interactions mediated by aromatics and by planar-polar residues, in addition to less common features such as the aromatic sandwich. The rules also reveal a previously unreported dependency for residues cys and leu. They also specify interactions involving aromatic and hydrogen bonding residues. This paper shows that Inductive Logic Programming implemented in ProGolem can derive rules giving structural features of protein/ligand interactions. Several of these rules are consistent with descriptions in the literature.

CONCLUSIONS

In addition to confirming literature results, ProGolem's model has a 10-fold cross-validated predictive accuracy that is superior, at the 95% confidence level, to another ILP system previously used to study protein/hexose interactions and is comparable with state-of-the-art statistical learners.

Identification of mannose interacting residues using local composition

Background

Mannose binding proteins (MBPs) play a vital role in several biological functions such as defense mechanisms. These proteins bind to mannose on the surface of a wide range of pathogens and help in eliminating these pathogens from our body. Thus, it is important to identify mannose interacting residues (MIRs) in order to understand mechanism of recognition of pathogens by MBPs.

Results

This paper describes modules developed for predicting MIRs in a protein. Support vector machine (SVM) based models have been developed on 120 mannose binding protein chains, where no two chains have more than 25% sequence similarity. SVM models were developed on two types of datasets: 1) main dataset consists of 1029 mannose interacting and 1029 non-interacting residues, 2) realistic dataset consists of 1029 mannose interacting and 10320 non-interacting residues. In this study, firstly, we developed standard modules using binary and PSSM profile of patterns and got maximum MCC around 0.32. Secondly, we developed SVM modules using composition profile of patterns and achieved maximum MCC around 0.74 with accuracy 86.64% on main dataset. Thirdly, we developed a model on a realistic dataset and achieved maximum MCC of 0.62 with accuracy 93.08%. Based on this study, a standalone program and web server have been developed for predicting mannose interacting residues in proteins (http://www.imtech.res.in/raghava/premier/).

Conclusions

Compositional analysis of mannose interacting and non-interacting residues shows that certain types of residues are preferred in mannose interaction. It was also observed that residues around mannose interacting residues have a preference for certain types of residues. Composition of patterns/peptide/segment has been used for predicting MIRs and achieved reasonable high accuracy. It is possible that this novel strategy may be effective to predict other types of interacting residues. This study will be useful in annotating the function of protein as well as in understanding the role of mannose in the immune system.

Structural motif screening reveals a novel, conserved carbohydrate-binding surface in the pathogenesis-related protein PR-5d

BACKGROUND

Aromatic amino acids play a critical role in protein-glycan interactions. Clusters of surface aromatic residues and their features may therefore be useful in distinguishing glycan-binding sites as well as predicting novel glycan-binding proteins. In this work, a structural bioinformatics approach was used to screen the Protein Data Bank (PDB) for coplanar aromatic motifs similar to those found in known glycan-binding proteins.

RESULTS

The proteins identified in the screen were significantly associated with carbohydrate-related functions according to gene ontology (GO) enrichment analysis, and predicted motifs were found frequently within novel folds and glycan-binding sites not included in the training set. In addition to numerous binding sites predicted in structural genomics proteins of unknown function, one novel prediction was a surface motif (W34/W36/W192) in the tobacco pathogenesis-related protein, PR-5d. Phylogenetic analysis revealed that the surface motif is exclusive to a subfamily of PR-5 proteins from the Solanaceae family of plants, and is absent completely in more distant homologs. To confirm PR-5d's insoluble-polysaccharide binding activity, a cellulose-pulldown assay of tobacco proteins was performed and PR-5d was identified in the cellulose-binding fraction by mass spectrometry.

CONCLUSIONS

Based on the combined results, we propose that the putative binding site in PR-5d may be an evolutionary adaptation of Solanaceae plants including potato, tomato, and tobacco, towards defense against cellulose-containing pathogens such as species of the deadly oomycete genus, Phytophthora. More generally, the results demonstrate that coplanar aromatic clusters on protein surfaces are a structural signature of glycan-binding proteins, and can be used to computationally predict novel glycan-binding proteins from 3 D structure.

An Inductive Logic Programming Approach to Validate Hexose Binding Biochemical Knowledge

Hexoses are simple sugars that play a key role in many cellular pathways, and in the regulation of development and disease mechanisms. Current protein-sugar computational models are based, at least partially, on prior biochemical findings and knowledge. They incorporate different parts of these findings in predictive black-box models. We investigate the empirical support for biochemical findings by comparing Inductive Logic Programming (ILP) induced rules to actual biochemical results. We mine the Protein Data Bank for a representative data set of hexose binding sites, non-hexose binding sites and surface grooves. We build an ILP model of hexose-binding sites and evaluate our results against several baseline machine learning classifiers. Our method achieves an accuracy similar to that of other black-box classifiers while providing insight into the discriminating process. In addition, it confirms wet-lab findings and reveals a previously unreported Trp-Glu amino acids dependency.

Transcriptome analysis of agmatine and putrescine catabolism in Pseudomonas aeruginosa PAO1

Polyamines (putrescine, spermidine, and spermine) are major organic polycations essential for a wide spectrum of cellular processes. The cells require mechanisms to maintain homeostasis of intracellular polyamines to prevent otherwise severe adverse effects. We performed a detailed transcriptome profile analysis of Pseudomonas aeruginosa in response to agmatine and putrescine with an emphasis in polyamine catabolism. Agmatine serves as the precursor compound for putrescine (and hence spermidine and spermine), which was proposed to convert into 4-aminobutyrate (GABA) and succinate before entering the tricarboxylic acid cycle in support of cell growth, as the sole source of carbon and nitrogen. Two acetylpolyamine amidohydrolases, AphA and AphB, were found to be involved in the conversion of agmatine into putrescine. Enzymatic products of AphA were confirmed by mass spectrometry analysis. Interestingly, the alanine-pyruvate cycle was shown to be indispensable for polyamine utilization. The newly identified dadRAX locus encoding the regulator alanine transaminase and racemase coupled with SpuC, the major putrescine-pyruvate transaminase, were key components to maintaining alanine homeostasis. Corresponding mutant strains were severely hampered in polyamine utilization. On the other hand, an alternative gamma-glutamylation pathway for the conversion of putrescine into GABA is present in some organisms. Subsequently, GabD, GabT, and PA5313 were identified for GABA utilization. The growth defect of the PA5313 gabT double mutant in GABA suggested the importance of these two transaminases. The succinic-semialdehyde dehydrogenase activity of GabD and its induction by GABA were also demonstrated in vitro. Polyamine utilization in general was proven to be independent of the PhoPQ two-component system, even though a modest induction of this operon was induced by polyamines. Multiple potent catabolic pathways, as depicted in this study, could serve pivotal roles in the control of intracellular polyamine levels.