On the hierarchical classification of G protein-coupled receptors

The role of leucine and isoleucine in tuning the hydropathy of class A GPCRs

Leucine and Isoleucine are two amino acids that differ only by the positioning of one methyl group. This small difference can have important consequences in α-helices, as the β-branching of Ile results in helix destabilization. We set out to investigate whether there are general trends for the occurrences of Leu and Ile residues in the structures and sequences of class A GPCRs (G protein-coupled receptors). GPCRs are integral membrane proteins in which α-helices span the plasma membrane seven times and which play a crucial role in signal transmission. We found that Leu side chains are generally more exposed at the protein surface than Ile side chains. We explored whether this difference might be attributed to different functions of the two amino acids and tested if Leu tunes the hydrophobicity of the transmembrane domain based on the Wimley-White whole-residue hydrophobicity scales. Leu content decreases the variation in hydropathy between receptors and correlates with the non-Leu receptor hydropathy. Both measures indicate that hydropathy is tuned by Leu. To test this idea further, we generated protein sequences with random amino acid compositions using a simple numerical model, in which hydropathy was tuned by adjusting the number of Leu residues. The model was able to replicate the observations made with class A GPCR sequences. We speculate that the hydropathy of transmembrane domains of class A GPCRs is tuned by Leu (and to some lesser degree by Lys and Val) to facilitate correct insertion into membranes and/or to stably anchor the receptors within membranes.

Mechanisms and molecular targets surrounding the potential therapeutic effects of psychedelics

Psychedelics, also known as classical hallucinogens, have been investigated for decades due to their potential therapeutic effects in the treatment of neuropsychiatric and substance use disorders. The results from clinical trials have shown promise for the use of psychedelics to alleviate symptoms of depression and anxiety, as well as to promote substantial decreases in the use of nicotine and alcohol. While these studies provide compelling evidence for the powerful subjective experience and prolonged therapeutic adaptations, the underlying molecular reasons for these robust and clinically meaningful improvements are still poorly understood. Preclinical studies assessing the targets and circuitry of the post-acute effects of classical psychedelics are ongoing. Current literature is split between a serotonin 5-HT2A receptor (5-HT2AR)-dependent or -independent signaling pathway, as researchers are attempting to harness the mechanisms behind the sustained post-acute therapeutically relevant effects. A combination of molecular, behavioral, and genetic techniques in neuropharmacology has begun to show promise for elucidating these mechanisms. As the field progresses, increasing evidence points towards the importance of the subjective experience induced by psychedelic-assisted therapy, but without further cross validation between clinical and preclinical research, the why behind the experience and its translational validity may be lost.

To Affinity and Beyond: A Personal Reflection on the Design and Discovery of Drugs

Faced with new and as yet unmet medical need, the stark underperformance of the pharmaceutical discovery process is well described if not perfectly understood. Driven primarily by profit rather than societal need, the search for new pharmaceutical products—small molecule drugs, biologicals, and vaccines—is neither properly funded nor sufficiently systematic. Many innovative approaches remain significantly underused and severely underappreciated, while dominant methodologies are replete with problems and limitations. Design is a component of drug discovery that is much discussed but seldom realised. In and of itself, technical innovation alone is unlikely to fulfil all the possibilities of drug discovery if the necessary underlying infrastructure remains unaltered. A fundamental revision in attitudes, with greater reliance on design powered by computational approaches, as well as a move away from the commercial imperative, is thus essential to capitalise fully on the potential of pharmaceutical intervention in healthcare.

EnsembleFam: towards more accurate protein family prediction in the twilight zone

Background

Current protein family modeling methods like profile Hidden Markov Model (pHMM), k-mer based methods, and deep learning-based methods do not provide very accurate protein function prediction for proteins in the twilight zone, due to low sequence similarity to reference proteins with known functions.

Results

We present a novel method EnsembleFam, aiming at better function prediction for proteins in the twilight zone. EnsembleFam extracts the core characteristics of a protein family using similarity and dissimilarity features calculated from sequence homology relations. EnsembleFam trains three separate Support Vector Machine (SVM) classifiers for each family using these features, and an ensemble prediction is made to classify novel proteins into these families. Extensive experiments are conducted using the Clusters of Orthologous Groups (COG) dataset and G Protein-Coupled Receptor (GPCR) dataset. EnsembleFam not only outperforms state-of-the-art methods on the overall dataset but also provides a much more accurate prediction for twilight zone proteins.

Conclusions

EnsembleFam, a machine learning method to model protein families, can be used to better identify members with very low sequence homology. Using EnsembleFam protein functions can be predicted  using just sequence information with better accuracy than state-of-the-art methods.

pdCSM-GPCR: predicting potent GPCR ligands with graph-based signatures

MOTIVATION

G protein-coupled receptors (GPCRs) can selectively bind to many types of ligands, ranging from light-sensitive compounds, ions, hormones, pheromones and neurotransmitters, modulating cell physiology. Considering their role in many essential cellular processes, they are one of the most targeted protein families, with over a third of all approved drugs modulating GPCR signalling. Despite this, the large diversity of receptors and their multipass transmembrane architectures make the identification and development of novel specific, and safe GPCR ligands a challenge. While computational approaches have the potential to assist GPCR drug development, they have presented limited performance and generalization capabilities. Here, we explored the use of graph-based signatures to develop pdCSM-GPCR, a method capable of rapidly and accurately screening potential GPCR ligands.

RESULTS

Bioactivity data (IC50, EC50, Ki and Kd) for individual GPCRs were curated. After curation, we used the data for developing predictive models for 36 major GPCR targets, across 4 classes (A, B, C and F). Our models compose the most comprehensive computational resource for GPCR bioactivity prediction to date. Across stratified 10-fold cross-validation and blind tests, our approach achieved Pearson's correlations of up to 0.89, significantly outperforming previous methods. Interpreting our results, we identified common important features of potent GPCRs ligands, which tend to have bicyclic rings, leading to higher levels of aromaticity. We believe pdCSM-GPCR will be an invaluable tool to assist screening efforts, enriching compound libraries and ranking candidates for further experimental validation.

AVAILABILITY AND IMPLEMENTATION

pdCSM-GPCR predictive models and datasets used have been made available via a freely accessible and easy-to-use web server at http://biosig.unimelb.edu.au/pdcsm_gpcr/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics Advances online.

An improved deep learning model for hierarchical classification of protein families

Although genes carry information, proteins are the main role player in providing all the functionalities of a living organism. Massive amounts of different proteins involve in every function that occurs in a cell. These amino acid sequences can be hierarchically classified into a set of families and subfamilies depending on their evolutionary relatedness and similarities in their structure or function. Protein characterization to identify protein structure and function is done accurately using laboratory experiments. With the rapidly increasing huge amount of novel protein sequences, these experiments have become difficult to carry out since they are expensive, time-consuming, and laborious. Therefore, many computational classification methods are introduced to classify proteins and predict their functional properties. With the progress of the performance of the computational techniques, deep learning plays a key role in many areas. Novel deep learning models such as DeepFam, ProtCNN have been presented to classify proteins into their families recently. However, these deep learning models have been used to carry out the non-hierarchical classification of proteins. In this research, we propose a deep learning neural network model named DeepHiFam with high accuracy to classify proteins hierarchically into different levels simultaneously. The model achieved an accuracy of 98.38% for protein family classification and more than 80% accuracy for the classification of protein subfamilies and sub-subfamilies. Further, DeepHiFam performed well in the non-hierarchical classification of protein families and achieved an accuracy of 98.62% and 96.14% for the popular Pfam dataset and COG dataset respectively.

Development and validation of multiple machine learning algorithms for the classification of G-protein-coupled receptors using molecular evolution model-based feature extraction strategy

Machine learning is one of the most potential ways to realize the function prediction of the incremental large-scale G-protein-coupled receptors (GPCR). Prior research reveals that the key to determining the overall classification accuracy of GPCR is extracting valuable features and filtering out redundancy. To achieve a more efficient classification model, we put the feature synonym problem into consideration and create a new method based on functional word clustering and integration. Through evaluating the evolution correlation between features using the transition scores in mature molecular substitution matrices, candidate features are clustered into synonym groups. Each group of the clustered features is then integrated and represented by a unique key functional word. These retained key functional words are used to form a feature knowledge base. The original GPCR sequences are then transferred into feature vectors based on a feature re-extraction strategy according to the features in the knowledge base before the training and testing stage. We create multiple machine learning models based on Naïve Bayesian (NB), random forest (RF), support vector machine (SVM), and multi-layer perceptron (MLP) algorithms. The established model is applied to classify two public data sets containing 8354 and 12,731 GPCRs, respectively. These models achieve significant performance in almost all evaluation criteria in comparison with state-of-the art. This work demonstrated the potential of the novel feature extraction strategy and provided an effective theoretical design for the hierarchical classification of GPCRs.

ADACT: a tool for analysing (dis)similarity among nucleotide and protein sequences using minimal and relative absent words

MOTIVATION

Researchers and practitioners use a number of popular sequence comparison tools that use many alignment-based techniques. Due to high time and space complexity and length-related restrictions, researchers often seek alignment-free tools. Recently, some interesting ideas, namely, Minimal Absent Words (MAW) and Relative Absent Words (RAW), have received much interest among the scientific community as distance measures that can give us alignment-free alternatives. This drives us to structure a framework for analysing biological sequences in an alignment-free manner.

RESULTS

In this application note, we present Alignment-free Dissimilarity Analysis & Comparison Tool (ADACT), a simple web-based tool that computes the analogy among sequences using a varied number of indexes through the distance matrix, species relation list and phylogenetic tree. This tool basically combines absent word (MAW or RAW) computation, dissimilarity measures, species relationship and thus brings all required software in one platform for the ease of researchers and practitioners alike in the field of bioinformatics. We have also developed a restful API.

AVAILABILITY AND IMPLEMENTATION

ADACT has been hosted at http://research.buet.ac.bd/ADACT/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The LPA3 Receptor: Regulation and Activation of Signaling Pathways

The lysophosphatidic acid 3 receptor (LPA₃) participates in different physiological actions and in the pathogenesis of many diseases through the activation of different signal pathways. Knowledge of the regulation of the function of the LPA₃ receptor is a crucial element for defining its roles in health and disease. This review describes what is known about the signaling pathways activated in terms of its various actions. Next, we review knowledge on the structure of the LPA₃ receptor, the domains found, and the roles that the latter might play in ligand recognition, signaling, and cellular localization. Currently, there is some information on the action of LPA₃ in different cells and whole organisms, but very little is known about the regulation of its function. Areas in which there is a gap in our knowledge are indicated in order to further stimulate experimental work on this receptor and on other members of the LPA receptor family. We are convinced that knowledge on how this receptor is activated, the signaling pathways employed and how the receptor internalization and desensitization are controlled will help design new therapeutic interventions for treating diseases in which the LPA₃ receptor is implicated.

Deep hierarchical embedding for simultaneous modeling of GPCR proteins in a unified metric space

GPCR proteins belong to diverse families of proteins that are defined at multiple hierarchical levels. Inspecting relationships between GPCR proteins on the hierarchical structure is important, since characteristics of the protein can be inferred from proteins in similar hierarchical information. However, modeling of GPCR families has been performed separately for each of the family, subfamily, and sub-subfamily level. Relationships between GPCR proteins are ignored in these approaches as they process the information in the proteins with several disconnected models. In this study, we propose DeepHier, a deep learning model to simultaneously learn representations of GPCR family hierarchy from the protein sequences with a unified single model. Novel loss term based on metric learning is introduced to incorporate hierarchical relations between proteins. We tested our approach using a public GPCR sequence dataset. Metric distances in the deep feature space corresponded to the hierarchical family relation between GPCR proteins. Furthermore, we demonstrated that further downstream tasks, like phylogenetic reconstruction and motif discovery, are feasible in the constructed embedding space. These results show that hierarchical relations between sequences were successfully captured in both of technical and biological aspects.

Engineering G protein-coupled receptor signalling in yeast for biotechnological and medical purposes

G protein-coupled receptors (GPCRs) comprise the largest class of membrane proteins in the human genome, with a common denominator of seven-transmembrane domains largely conserved among eukaryotes. Yeast is naturally armoured with three different GPCRs for pheromone and sugar sensing, with the pheromone pathway being extensively hijacked for characterising heterologous GPCR signalling in a model eukaryote. This review focusses on functional GPCR studies performed in yeast and on the elucidated hotspots for engineering, and discusses both endogenous and heterologous GPCR signalling. Key emphasis will be devoted to studies describing important engineering parameters to consider for successful coupling of GPCRs to the yeast mating pathway. We also review the various means of applying yeast for studying GPCRs, including the use of yeast armed with heterologous GPCRs as a platform for (i) deorphanisation of orphan receptors, (ii) metabolic engineering of yeast for production of bioactive products and (iii) medical applications related to pathogen detection and drug discovery. Finally, this review summarises the current challenges related to expression of functional membrane-bound GPCRs in yeast and discusses the opportunities to continue capitalising on yeast as a model chassis for functional GPCR signalling studies.

Improving protein domain classification for third-generation sequencing reads using deep learning

BACKGROUND

With the development of third-generation sequencing (TGS) technologies, people are able to obtain DNA sequences with lengths from 10s to 100s of kb. These long reads allow protein domain annotation without assembly, thus can produce important insights into the biological functions of the underlying data. However, the high error rate in TGS data raises a new challenge to established domain analysis pipelines. The state-of-the-art methods are not optimized for noisy reads and have shown unsatisfactory accuracy of domain classification in TGS data. New computational methods are still needed to improve the performance of domain prediction in long noisy reads.

RESULTS

In this work, we introduce ProDOMA, a deep learning model that conducts domain classification for TGS reads. It uses deep neural networks with 3-frame translation encoding to learn conserved features from partially correct translations. In addition, we formulate our problem as an open-set problem and thus our model can reject reads not containing the targeted domains. In the experiments on simulated long reads of protein coding sequences and real TGS reads from the human genome, our model outperforms HMMER and DeepFam on protein domain classification.

CONCLUSIONS

In summary, ProDOMA is a useful end-to-end protein domain analysis tool for long noisy reads without relying on error correction.

Biased Ghrelin Receptor Signaling and the Dopaminergic System as Potential Targets for Metabolic and Psychological Symptoms of Anorexia Nervosa

Anorexia Nervosa (AN) is a complex disease that impairs the metabolic, mental and physiological health of affected individuals in a severe and sometimes lethal way. Many of the common symptoms in AN patients, such as reduced food intake, anxiety, impaired gut motility or overexercising are connected to both the orexigenic gut hormone ghrelin and the dopaminergic system. Targeting the ghrelin receptor (GhrR) to treat AN seems a promising possibility in current research. However, GhrR signaling is highly complex. First, the GhrR can activate four known intracellular pathways Gαq, Gαi/o, Gα12/13 and the recruitment of β-arrestin. Biased signaling provides the possibility to activate or inhibit only one or a subset of the intracellular pathways of a pleiotropic receptor. This allows specific targeting of physiological functions without adverse effects. Currently little is known on how biased signaling could specifically modulate GhrR effects. Second, GhrR signaling has been shown to be interconnected with the dopaminergic system, particularly in the context of AN symptoms. This review highlights that a biased agonist for the GhrR may be a promising target for the treatment of AN, however extensive and systematic translational studies are still needed and the connection to the dopaminergic system has to be taken into account.

G-Protein coupled receptors: structure and function in drug discovery

The G-protein coupled receptors (GPCRs) superfamily comprise similar proteins arranged into families or classes thus making it one of the largest in the mammalian genome. GPCRs take part in many vital physiological functions making them targets for numerous novel drugs. GPCRs share some distinctive features, such as the seven transmembrane domains, they also differ in the number of conserved residues in their transmembrane domain. Here we provide an introductory and accessible review detailing the computational advances in GPCR pharmacology and drug discovery. An overview is provided on family A-C GPCRs; their structural differences, GPCR signalling, allosteric binding and cooperativity. The dielectric constant (relative permittivity) of proteins is also discussed in the context of site-specific environmental effects.

DeepFam: deep learning based alignment-free method for protein family modeling and prediction

Motivation

A large number of newly sequenced proteins are generated by the next-generation sequencing technologies and the biochemical function assignment of the proteins is an important task. However, biological experiments are too expensive to characterize such a large number of protein sequences, thus protein function prediction is primarily done by computational modeling methods, such as profile Hidden Markov Model (pHMM) and k-mer based methods. Nevertheless, existing methods have some limitations; k-mer based methods are not accurate enough to assign protein functions and pHMM is not fast enough to handle large number of protein sequences from numerous genome projects. Therefore, a more accurate and faster protein function prediction method is needed.

Results

In this paper, we introduce DeepFam, an alignment-free method that can extract functional information directly from sequences without the need of multiple sequence alignments. In extensive experiments using the Clusters of Orthologous Groups (COGs) and G protein-coupled receptor (GPCR) dataset, DeepFam achieved better performance in terms of accuracy and runtime for predicting functions of proteins compared to the state-of-the-art methods, both alignment-free and alignment-based methods. Additionally, we showed that DeepFam has a power of capturing conserved regions to model protein families. In fact, DeepFam was able to detect conserved regions documented in the Prosite database while predicting functions of proteins. Our deep learning method will be useful in characterizing functions of the ever increasing protein sequences.

Availability and implementation

Codes are available at https://bhi-kimlab.github.io/DeepFam.

DEEPre: sequence-based enzyme EC number prediction by deep learning

Motivation

Annotation of enzyme function has a broad range of applications, such as metagenomics, industrial biotechnology, and diagnosis of enzyme deficiency-caused diseases. However, the time and resource required make it prohibitively expensive to experimentally determine the function of every enzyme. Therefore, computational enzyme function prediction has become increasingly important. In this paper, we develop such an approach, determining the enzyme function by predicting the Enzyme Commission number.

Results

We propose an end-to-end feature selection and classification model training approach, as well as an automatic and robust feature dimensionality uniformization method, DEEPre, in the field of enzyme function prediction. Instead of extracting manually crafted features from enzyme sequences, our model takes the raw sequence encoding as inputs, extracting convolutional and sequential features from the raw encoding based on the classification result to directly improve the prediction performance. The thorough cross-fold validation experiments conducted on two large-scale datasets show that DEEPre improves the prediction performance over the previous state-of-the-art methods. In addition, our server outperforms five other servers in determining the main class of enzymes on a separate low-homology dataset. Two case studies demonstrate DEEPre’s ability to capture the functional difference of enzyme isoforms.

Availability and implementation

The server could be accessed freely at http://www.cbrc.kaust.edu.sa/DEEPre.

Supplementary information

Supplementary data are available at Bioinformatics online.

GPCR-SAS: A web application for statistical analyses on G protein-coupled receptors sequences

G protein-coupled receptors (GPCRs) are one of the largest protein families in mammals. They mediate signal transduction across cell membranes and are important targets for the pharmaceutical industry. The G Protein-Coupled Receptors-Sequence Analysis and Statistics (GPCR-SAS) web application provides a set of tools to perform comparative analysis of sequence positions between receptors, based on a curated structural-informed multiple sequence alignment. The analysis tools include: (i) percentage of occurrence of an amino acid or motif and entropy at a position or range of positions, (ii) covariance of two positions, (iii) correlation between two amino acids in two positions (or two sequence motifs in two ranges of positions), and (iv) snake-plot representation for a specific receptor or for the consensus sequence of a group of selected receptors. The analysis of conservation of residues and motifs across transmembrane (TM) segments may guide the design of more selective ligands or help to rationalize activation mechanisms, among others. As an example, here we analyze the amino acids of the "transmission switch", that initiates receptor activation following ligand binding. The tool is freely accessible at http://lmc.uab.cat/gpcrsas/.

Hierarchical classification strategy for Phenotype extraction from epidermal growth factor receptor endocytosis screening

BACKGROUND

Endocytosis is regarded as a mechanism of attenuating the epidermal growth factor receptor (EGFR) signaling and of receptor degradation. There is increasing evidence becoming available showing that breast cancer progression is associated with a defect in EGFR endocytosis. In order to find related Ribonucleic acid (RNA) regulators in this process, high-throughput imaging with fluorescent markers is used to visualize the complex EGFR endocytosis process. Subsequently a dedicated automatic image and data analysis system is developed and applied to extract the phenotype measurement and distinguish different developmental episodes from a huge amount of images acquired through high-throughput imaging. For the image analysis, a phenotype measurement quantifies the important image information into distinct features or measurements. Therefore, the manner in which prominent measurements are chosen to represent the dynamics of the EGFR process becomes a crucial step for the identification of the phenotype. In the subsequent data analysis, classification is used to categorize each observation by making use of all prominent measurements obtained from image analysis. Therefore, a better construction for a classification strategy will support to raise the performance level in our image and data analysis system.

RESULTS

In this paper, we illustrate an integrated analysis method for EGFR signalling through image analysis of microscopy images. Sophisticated wavelet-based texture measurements are used to obtain a good description of the characteristic stages in the EGFR signalling. A hierarchical classification strategy is designed to improve the recognition of phenotypic episodes of EGFR during endocytosis. Different strategies for normalization, feature selection and classification are evaluated.

CONCLUSIONS

The results of performance assessment clearly demonstrate that our hierarchical classification scheme combined with a selected set of features provides a notable improvement in the temporal analysis of EGFR endocytosis. Moreover, it is shown that the addition of the wavelet-based texture features contributes to this improvement. Our workflow can be applied to drug discovery to analyze defected EGFR endocytosis processes.

Eukaryotic G protein-coupled receptors as descendants of prokaryotic sodium-translocating rhodopsins

Abstract

Microbial rhodopsins and G-protein coupled receptors (GPCRs, which include animal rhodopsins) are two distinct (super) families of heptahelical (7TM) membrane proteins that share obvious structural similarities but no significant sequence similarity. Comparison of the recently solved high-resolution structures of the sodium-translocating bacterial rhodopsin and various Na⁺-binding GPCRs revealed striking similarity of their sodium-binding sites. This similarity allowed us to construct a structure-guided sequence alignment for the two (super)families, which highlighted their evolutionary relatedness. Our analysis supports a common underlying molecular mechanism for both families that involves a highly conserved aromatic residue playing a pivotal role in rotation of the 6th transmembrane helix.

Reviewers

This article was reviewed by Oded Beja, G. P. S. Raghava and L. Aravind.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0091-4) contains supplementary material, which is available to authorized users.

Helix 8 of the angiotensin- II type 1A receptor interacts with phosphatidylinositol phosphates and modulates membrane insertion

The carboxyl-terminus of the type 1 angiotensin II receptor (AT_1A) regulates receptor activation/deactivation and the amphipathic Helix 8 within the carboxyl-terminus is a high affinity interaction motif for plasma membrane lipids. We have used dual polarisation interferometry (DPI) to examine the role of phosphatidylinositdes in the specific recognition of Helix 8 in the AT_1A receptor. A synthetic peptide corresponding to Leu³⁰⁵ to Lys³²⁵ (Helix 8 AT_1A) discriminated between PIPs and different charges on lipid membranes. Peptide binding to PtdIns(4)P-containing bilayers caused a dramatic change in the birefringence (a measure of membrane order) of the bilayer. Kinetic modelling showed that PtdIns(4)P is held above the bilayer until the mass of bound peptide reaches a threshold, after which the peptides insert further into the bilayer. This suggests that Helix 8 can respond to the presence of PI(4)P by withdrawing from the bilayer, resulting in a functional conformational change in the receptor.

Data mining strategies to improve multiplex microbead immunoassay tolerance in a mouse model of infectious diseases

Multiplex methodologies, especially those with high-throughput capabilities generate large volumes of data. Accumulation of such data (e.g., genomics, proteomics, metabolomics etc.) is fast becoming more common and thus requires the development and implementation of effective data mining strategies designed for biological and clinical applications. Multiplex microbead immunoassay (MMIA), on xMAP or MagPix platform (Luminex), which is amenable to automation, offers a major advantage over conventional methods such as Western blot or ELISA, for increasing the efficiencies in serodiagnosis of infectious diseases. MMIA allows detection of antibodies and/or antigens efficiently for a wide range of infectious agents simultaneously in host blood samples, in one reaction vessel. In the process, MMIA generates large volumes of data. In this report we demonstrate the application of data mining tools on how the inherent large volume data can improve the assay tolerance (measured in terms of sensitivity and specificity) by analysis of experimental data accumulated over a span of two years. The combination of prior knowledge with machine learning tools provides an efficient approach to improve the diagnostic power of the assay in a continuous basis. Furthermore, this study provides an in-depth knowledge base to study pathological trends of infectious agents in mouse colonies on a multivariate scale. Data mining techniques using serodetection of infections in mice, developed in this study, can be used as a general model for more complex applications in epidemiology and clinical translational research.

The influence of alignment-free sequence representations on the semi-supervised classification of class C G protein-coupled receptors: semi-supervised classification of class C GPCRs

G protein-coupled receptors (GPCRs) are integral cell membrane proteins of relevance for pharmacology. The tertiary structure of the transmembrane domain, a gate to the study of protein functionality, is unknown for almost all members of class C GPCRs, which are the target of the current study. As a result, their investigation must often rely on alignments of their amino acid sequences. Sequence alignment entails the risk of missing relevant information. Various approaches have attempted to circumvent this risk through alignment-free transformations of the sequences on the basis of different amino acid physicochemical properties. In this paper, we use several of these alignment-free methods, as well as a basic amino acid composition representation, to transform the available sequences. Novel semi-supervised statistical machine learning methods are then used to discriminate the different class C GPCRs types from the transformed data. This approach is relevant due to the existence of orphan proteins to which type labels should be assigned in a process of deorphanization or reverse pharmacology. The reported experiments show that the proposed techniques provide accurate classification even in settings of extreme class-label scarcity and that fair accuracy can be achieved even with very simple transformation strategies that ignore the sequence ordering.

GPCRsort-responding to the next generation sequencing data challenge: prediction of G protein-coupled receptor classes using only structural region lengths

Next generation sequencing (NGS) and the attendant data deluge are increasingly impacting molecular life sciences research. Chief among the challenges and opportunities is to enhance our ability to classify molecular target data into meaningful and cohesive systematic nomenclature. In this vein, the G protein-coupled receptors (GPCRs) are the largest and most divergent receptor family that plays a crucial role in a host of pathophysiological pathways. For the pharmaceutical industry, GPCRs are a major drug target and it is estimated that 60%-70% of all medicines in development today target GPCRs. Hence, they require an efficient and rapid classification to group the members according to their functions. In addition to NGS and the Big Data challenge we currently face, an emerging number of orphan GPCRs further demand for novel, rapid, and accurate classification of the receptors since the current classification tools are inadequate and slow. This study presents the development of a new classification tool for GPCRs using the structural features derived from their primary sequences: GPCRsort. Comparison experiments with the current known GPCR classification techniques showed that GPCRsort is able to rapidly (in the order of minutes) classify uncharacterized GPCRs with 97.3% accuracy, whereas the best available technique's accuracy is 90.7%. GPCRsort is available in the public domain for postgenomics life scientists engaged in GPCR research with NGS: http://bioserver.ceng.metu.edu.tr/GPCRSort .

Bioinformatics tools for predicting GPCR gene functions

The automatic classification of GPCRs by bioinformatics methodology can provide functional information for new GPCRs in the whole 'GPCR proteome' and this information is important for the development of novel drugs. Since GPCR proteome is classified hierarchically, general ways for GPCR function prediction are based on hierarchical classification. Various computational tools have been developed to predict GPCR functions; those tools use not simple sequence searches but more powerful methods, such as alignment-free methods, statistical model methods, and machine learning methods used in protein sequence analysis, based on learning datasets. The first stage of hierarchical function prediction involves the discrimination of GPCRs from non-GPCRs and the second stage involves the classification of the predicted GPCR candidates into family, subfamily, and sub-subfamily levels. Then, further classification is performed according to their protein-protein interaction type: binding G-protein type, oligomerized partner type, etc. Those methods have achieved predictive accuracies of around 90 %. Finally, I described the future subject of research of the bioinformatics technique about functional prediction of GPCR.

Targeting G protein coupled receptor-related pathways as emerging molecular therapies

G protein coupled receptors (GPCRs) represent the most important targets in modern pharmacology because of the different functions they mediate, especially within brain and peripheral nervous system, and also because of their functional and stereochemical properties. In this paper, we illustrate, via a variety of examples, novel advances about the GPCR-related molecules that have been shown to play diverse roles in GPCR pathways and in pathophysiological phenomena. We have exemplified how those GPCRs’ pathways are, or might constitute, potential targets for different drugs either to stimulate, modify, regulate or inhibit the cellular mechanisms that are hypothesized to govern some pathologic, physiologic, biologic and cellular or molecular aspects both in vivo and in vitro. Therefore, influencing such pathways will, undoubtedly, lead to different therapeutical applications based on the related pharmacological implications. Furthermore, such new properties can be applied in different fields. In addition to offering fruitful directions for future researches, we hope the reviewed data, together with the elements found within the cited references, will inspire clinicians and researchers devoted to the studies on GPCR’s properties.

Identification of Candidate Vaccine Antigens In Silico

The identification of immunogenic whole-protein antigens is fundamental to the successful discovery of candidate subunit vaccines and their rapid, effective, and efficient transformation into clinically useful, commercially successful vaccine formulations. In the wider context of the experimental discovery of vaccine antigens, with particular reference to reverse vaccinology, this chapter adumbrates the principal computational approaches currently deployed in the hunt for novel antigens: genome-level prediction of antigens, antigen identification through the use of protein sequence alignment-based approaches, antigen detection through the use of subcellular location prediction, and the use of alignment-independent approaches to antigen discovery. Reference is also made to the recent emergence of various expert systems for protein antigen identification.

Classification of GPCRs using family specific motifs

The classification of G-Protein Coupled Receptor (GPCR) sequences is an important problem that arises from the need to close the gap between the large number of orphan receptors and the relatively small number of annotated receptors. Equally important is the characterization of GPCR Class A subfamilies and gaining insight into the ligand interaction since GPCR Class A encompasses a very large number of drug-targeted receptors. In this work, we propose a method for Class A subfamily classification using sequence-derived motifs which characterizes the subfamilies by discovering receptor-ligand interaction sites. The motifs that best characterize a subfamily are selected by the Distinguishing Power Evaluation (DPE) technique we propose. The experiments performed on GPCR sequence databases show that our method outperforms state-of-the-art classification techniques for GPCR Class A subfamily prediction. An important contribution of our work is to discover key receptor-ligand interaction sites which is very important for drug design.

The G protein-coupled receptors in the pufferfish Takifugu rubripes

BACKGROUND

Guanine protein-coupled receptors (GPCRs) constitute a eukaryotic transmembrane protein family and function as "molecular switches" in the second messenger cascades and are found in all organisms between yeast and humans. They form the single, biggest drug-target family due to their versatility of action and their role in several physiological functions, being active players in detecting the presence of light, a variety of smells and tastes, amino acids, nucleotides, lipids, chemicals etc. in the environment of the cell. Comparative genomic studies on model organisms provide information on target receptors in humans and their function. The Japanese teleost Fugu has been identified as one of the smallest vertebrate genomes and a compact model to study the human genome, owing to the great similarity in its gene repertoire with that of human and other vertebrates. Thus the characterization of the GPCRs of Fugu would provide insights to the evolution of the vertebrate genome.

RESULTS

We classified the GPCRs in the Fugu genome and our analysis of its 316 membrane-bound receptors, available on the public databases as well as from literature, detected 298 GPCRs that were grouped into five main families according to the GRAFS classification system (namely, Glutamate, Rhodopsin, Adhesion, Frizzled and Secretin). We also identified 18 other GPCRs that could not be grouped under the GRAFS family and hence were classified as 'Other 7TM' receptors. On comparison of the GPCR information from the Fugu genome with those in the human and chicken genomes, we detected 96.83% (306/316) and 96.51% (305/316) orthology in GPCRs among the Fugu-human genomes and Fugu-chicken genomes, respectively.

CONCLUSIONS

This study reveals the position of pisces in vertebrate evolution from the GPCR perspective. Fugu can act as a reference model for the human genome for other protein families as well, going by the high orthology observed for GPCRs between Fugu and human. The evolutionary comparison of GPCR sequences between key vertebrate classes of mammals, birds and fish will help in identifying key functional residues and motifs so as to fill in the blanks in the evolution of GPCRs in vertebrates.

Identification of a novel endoplasmic reticulum export motif within the eighth α-helical domain (α-H8) of the human prostacyclin receptor

The human prostacyclin receptor (hIP) undergoes agonist-dependent trafficking involving a direct interaction with Rab11a GTPase. The region of interaction was localised to a 14 residue Rab11a binding domain (RBD) within the proximal carboxyl-terminal (C)-tail domain of the hIP, consisting of Val(299)-Val(307) within the eighth helical domain (α-H8) adjacent to the palmitoylated residues at Cys(308)-Cys(311). However, the factors determining the anterograde transport of the newly synthesised hIP from the endoplasmic reticulum (ER) to the plasma membrane (PM) have not been identified. The aim of the current study was to identify the major ER export motif(s) within the hIP initially by investigating the role of Lys residues in its maturation and processing. Through site-directed and Ala-scanning mutational studies in combination with analyses of protein expression and maturation, functional analyses of ligand binding, agonist-induced intracellular signalling and confocal image analyses, it was determined that Lys(297), Arg(302) and Lys(304) located within α-H8 represent the critical determinants of a novel ER export motif of the hIP. Furthermore, while substitution of those critical residues significantly impaired maturation and processing of the hIP, replacement of the positively charged Lys with Arg residues, and vice versa, was functionally permissible. Hence, this study has identified a novel 8 residue ER export motif within the functionally important α-H8 of the hIP. This ER export motif, defined by "K/R(X)(4)K/R(X)K/R," has a strict requirement for positively charged, basic Lys/Arg residues at the 1st, 6th and 8th positions and appears to be evolutionarily conserved within IP sequences from mouse to man.

Expression of G protein-coupled receptors and related proteins in HEK293, AtT20, BV2, and N18 cell lines as revealed by microarray analysis

BACKGROUND

G protein coupled receptors (GPCRs) are one of the most widely studied gene superfamilies. Thousands of GPCR research studies have utilized heterologous expression systems such as human embryonic kidney cells (HEK293). Though often treated as 'blank slates', these cell lines nevertheless endogenously express GPCRs and related signaling proteins. The outcome of a given GPCR study can be profoundly influenced by this largely unknown complement of receptors and/or signaling proteins. Little easily accessible information exists that describes the expression profiles of the GPCRs in cell lines. What is accessible is often limited in scope - of the hundreds of GPCRs and related proteins, one is unlikely to find information on expression of more than a dozen proteins in a given cell line. Microarray technology has allowed rapid analysis of mRNA levels of thousands of candidate genes, but though often publicly available, the results can be difficult to efficiently access or even to interpret.

RESULTS

To bridge this gap, we have used microarrays to measure the mRNA levels of a comprehensive profile of non-chemosensory GPCRs and over a hundred GPCR signaling related gene products in four cell lines frequently used for GPCR research: HEK293, AtT20, BV2, and N18.

CONCLUSIONS

This study provides researchers an easily accessible mRNA profile of the endogenous signaling repertoire that these four cell lines possess. This will assist in choosing the most appropriate cell line for studying GPCRs and related signaling proteins. It also provides a better understanding of the potential interactions between GPCRs and those signaling proteins.

An improved classification of G-protein-coupled receptors using sequence-derived features

BACKGROUND

G-protein-coupled receptors (GPCRs) play a key role in diverse physiological processes and are the targets of almost two-thirds of the marketed drugs. The 3 D structures of GPCRs are largely unavailable; however, a large number of GPCR primary sequences are known. To facilitate the identification and characterization of novel receptors, it is therefore very valuable to develop a computational method to accurately predict GPCRs from the protein primary sequences.

RESULTS

We propose a new method called PCA-GPCR, to predict GPCRs using a comprehensive set of 1497 sequence-derived features. The principal component analysis is first employed to reduce the dimension of the feature space to 32. Then, the resulting 32-dimensional feature vectors are fed into a simple yet powerful classification algorithm, called intimate sorting, to predict GPCRs at five levels. The prediction at the first level determines whether a protein is a GPCR or a non-GPCR. If it is predicted to be a GPCR, then it will be further predicted into certain family, subfamily, sub-subfamily and subtype by the classifiers at the second, third, fourth, and fifth levels, respectively. To train the classifiers applied at five levels, a non-redundant dataset is carefully constructed, which contains 3178, 1589, 4772, 4924, and 2741 protein sequences at the respective levels. Jackknife tests on this training dataset show that the overall accuracies of PCA-GPCR at five levels (from the first to the fifth) can achieve up to 99.5%, 88.8%, 80.47%, 80.3%, and 92.34%, respectively. We further perform predictions on a dataset of 1238 GPCRs at the second level, and on another two datasets of 167 and 566 GPCRs respectively at the fourth level. The overall prediction accuracies of our method are consistently higher than those of the existing methods to be compared.

CONCLUSIONS

The comprehensive set of 1497 features is believed to be capable of capturing information about amino acid composition, sequence order as well as various physicochemical properties of proteins. Therefore, high accuracies are achieved when predicting GPCRs at all the five levels with our proposed method.

Classification of G-protein coupled receptors based on support vector machine with maximum relevance minimum redundancy and genetic algorithm

Background

Because a priori knowledge about function of G protein-coupled receptors (GPCRs) can provide useful information to pharmaceutical research, the determination of their function is a quite meaningful topic in protein science. However, with the rapid increase of GPCRs sequences entering into databanks, the gap between the number of known sequence and the number of known function is widening rapidly, and it is both time-consuming and expensive to determine their function based only on experimental techniques. Therefore, it is vitally significant to develop a computational method for quick and accurate classification of GPCRs.

Results

In this study, a novel three-layer predictor based on support vector machine (SVM) and feature selection is developed for predicting and classifying GPCRs directly from amino acid sequence data. The maximum relevance minimum redundancy (mRMR) is applied to pre-evaluate features with discriminative information while genetic algorithm (GA) is utilized to find the optimized feature subsets. SVM is used for the construction of classification models. The overall accuracy with three-layer predictor at levels of superfamily, family and subfamily are obtained by cross-validation test on two non-redundant dataset. The results are about 0.5% to 16% higher than those of GPCR-CA and GPCRPred.

Conclusion

The results with high success rates indicate that the proposed predictor is a useful automated tool in predicting GPCRs. GPCR-SVMFS, a corresponding executable program for GPCRs prediction and classification, can be acquired freely on request from the authors.

Identification of protein functions using a machine-learning approach based on sequence-derived properties

Background

Predicting the function of an unknown protein is an essential goal in bioinformatics. Sequence similarity-based approaches are widely used for function prediction; however, they are often inadequate in the absence of similar sequences or when the sequence similarity among known protein sequences is statistically weak. This study aimed to develop an accurate prediction method for identifying protein function, irrespective of sequence and structural similarities.

Results

A highly accurate prediction method capable of identifying protein function, based solely on protein sequence properties, is described. This method analyses and identifies specific features of the protein sequence that are highly correlated with certain protein functions and determines the combination of protein sequence features that best characterises protein function. Thirty-three features that represent subtle differences in local regions and full regions of the protein sequences were introduced. On the basis of 484 features extracted solely from the protein sequence, models were built to predict the functions of 11 different proteins from a broad range of cellular components, molecular functions, and biological processes. The accuracy of protein function prediction using random forests with feature selection ranged from 94.23% to 100%. The local sequence information was found to have a broad range of applicability in predicting protein function.

Conclusion

We present an accurate prediction method using a machine-learning approach based solely on protein sequence properties. The primary contribution of this paper is to propose new PNPRD features representing global and/or local differences in sequences, based on positively and/or negatively charged residues, to assist in predicting protein function. In addition, we identified a compact and useful feature subset for predicting the function of various proteins. Our results indicate that sequence-based classifiers can provide good results among a broad range of proteins, that the proposed features are useful in predicting several functions, and that the combination of our and traditional features may support the creation of a discriminative feature set for specific protein functions.

Role of helix 8 in G protein-coupled receptors based on structure-function studies on the type 1 angiotensin receptor

G protein-coupled receptors (GPCRs) are transmembrane receptors that convert extracellular stimuli to intracellular signals. The type 1 angiotensin II receptor is a widely studied GPCR with roles in blood pressure regulation,water and salt balance and cell growth. The complex molecular and structural changes that underpin receptor activation and signaling are the focus of intense research. Increasingly, there is an appreciation that the plasma membrane participates in receptor function via direct, physical interactions that reciprocally modulate both lipid and receptor and provide microdomains for specialized activities. Reversible protein:lipid interactions are commonly mediated by amphipathic -helices in proteins and one such motif - a short helix, referred to as helix VIII/8 (H8), located at the start of the carboxyl (C)-terminus of GPCRs - is gaining recognition for its importance to GPCR function. Here, we review the identification of H8 in GPCRs and examine its capacity to sense and interact with diverse proteins and lipid environment, most notably with acidic lipids that include phosphatidylinositol phosphates.

GPCRTree: online hierarchical classification of GPCR function

Background

G protein-coupled receptors (GPCRs) play important physiological roles transducing extracellular signals into intracellular responses. Approximately 50% of all marketed drugs target a GPCR. There remains considerable interest in effectively predicting the function of a GPCR from its primary sequence.

Findings

Using techniques drawn from data mining and proteochemometrics, an alignment-free approach to GPCR classification has been devised. It uses a simple representation of a protein's physical properties. GPCRTree, a publicly-available internet server, implements an algorithm that classifies GPCRs at the class, sub-family and sub-subfamily level.

Conclusion

A selective top-down classifier was developed which assigns sequences within a GPCR hierarchy. Compared to other publicly available GPCR prediction servers, GPCRTree is considerably more accurate at every level of classification. The server has been available online since March 2008 at URL: .