Position-specific Scoring Matrices Research Articles

DeepRank-GNN-esm: a graph neural network for scoring protein-protein models using protein language model.

Protein-Protein interactions (PPIs) play critical roles in numerous cellular processes. By modelling the 3D structures of the correspond protein complexes valuable insights can be obtained, providing, e.g. starting points for drug and protein design. One challenge in the modelling process is however the identification of near-native models from the large pool of generated models. To this end we have previously developed DeepRank-GNN, a graph neural network that integrates structural and sequence information to enable effective pattern learning at PPI interfaces. Its main features are related to the Position Specific Scoring Matrices (PSSMs), which are computationally expensive to generate, significantly limits the algorithm's usability. We introduce here DeepRank-GNN-esm that includes as additional features protein language model embeddings from the ESM-2 model. We show that the ESM-2 embeddings can actually replace the PSSM features at no cost in-, or even better performance on two PPI-related tasks: scoring docking poses and detecting crystal artifacts. This new DeepRank version bypasses thus the need of generating PSSM, greatly improving the usability of the software and opening new application opportunities for systems for which PSSM profiles cannot be obtained or are irrelevant (e.g. antibody-antigen complexes). DeepRank-GNN-esm is freely available from https://github.com/DeepRank/DeepRank-GNN-esm.

AAindex-PPII: Predicting polyproline type II helix structure based on amino acid indexes with an improved BiGRU-TextCNN model.

The polyproline-II (PPII) structure domain is crucial in organisms' signal transduction, transcription, cell metabolism, and immune response. It is also a critical structural domain for specific vital disease-associated proteins. Recognizing PPII is essential for understanding protein structure and function. To accurately predict PPII in proteins, we propose a novel method, AAindex-PPII, which only adopts amino acid index to characterize protein sequences and uses a Bidirectional Gated Recurrent Unit (BiGRU)-Improved TextCNN composite deep learning model to predict PPII in proteins. Experimental results show that, when tested on the same datasets, our method outperforms the state-of-the-art BERT-PPII method, achieving an AUC value of 0.845 on the strict data and an AUC value of 0.813 on the non-strict data, which is 0.024 and 0.03 higher than that of the BERT-PPII method. This study demonstrates that our proposed method is simple and efficient for PPII prediction without using pre-trained large models or complex features such as position-specific scoring matrices.

Viral community composition ofhypersaline lakes.

Despite their widespread distribution and remarkable antiquity no RNA viruses definitively associated with the domain Archaea have been identified. In contrast, 17 families of DNA viruses are known to infect archaea. In an attempt to uncover more of the elusive archaeal virosphere, we investigated the metatranscriptomes of hypersaline lakes that are a rich source of archaea. We sequenced RNA extracted from water filter samples of Lake Tyrrell (Victoria, Australia) and cultures seeded from four lakes in Antarctica. To identify highly divergent viruses in these data, we employed a variety of search tools, including Hidden Markov models (HMMs) and position-specific scoring matrices (PSSMs). From this, we identified 12 highly divergent, RNA virus-like candidate sequences from the virus phyla Artverviricota, Duplornaviricota, Kitrinoviricota, Negarnaviricota, and Pisuviricota, including those with similarity to the RNA-dependent RNA polymerase (RdRp). An additional analysis with an artificial intelligence (AI)-based approach that utilises both sequence and structural information identified seven putative and highly divergent RdRp sequences of uncertain phylogenetic position. A sequence matching the Pisuviricota from Deep Lake in Antarctica had the strongest RNA virus signal. Analyses of the dinucleotide representation of the virus-like candidates in comparison to that of potential host species were in some cases compatible with an association to archaeal or bacterial hosts. Notably, however, the use of archaeal CRISPR spacers as a BLAST database failed to detect any RNA viruses. We also described DNA viruses from the families Pleolipoviridae, Sphaerolipoviridae, Halspiviridae, and the class Caudoviricetes. Although we were unable to provide definitive evidence the existence of an RNA virus of archaea in these hypersaline lakes, this study lays the foundations for further investigations of highly divergent RNA viruses in natural environments.

Block Aligner: an adaptive SIMD-accelerated aligner for sequences and position-specific scoring matrices.

Efficiently aligning sequences is a fundamental problem in bioinformatics. Many recent algorithms for computing alignments through Smith-Waterman-Gotoh dynamic programming (DP) exploit Single Instruction Multiple Data (SIMD) operations on modern CPUs for speed. However, these advances have largely ignored difficulties associated with efficiently handling complex scoring matrices or large gaps (insertions or deletions). We propose a new SIMD-accelerated algorithm called Block Aligner for aligning nucleotide and protein sequences against other sequences or position-specific scoring matrices. We introduce a new paradigm that uses blocks in the DP matrix that greedily shift, grow, and shrink. This approach allows regions of the DP matrix to be adaptively computed. Our algorithm reaches over 5-10 times faster than some previous methods while incurring an error rate of less than 3% on protein and long read datasets, despite large gaps and low sequence identities. Our algorithm is implemented for global, local, and X-drop alignments. It is available as a Rust library (with C bindings) at https://github.com/Daniel-Liu-c0deb0t/block-aligner.

Abstract 2023: CoPheeMap: a co-regulation map of 30,000 phosphosite illuminates the dark cancer phosphoproteome

Abstract Mass spectrometry-based phosphoproteomics enables proteome-wide analysis of protein phosphorylation in biological samples. As a prime example, the Clinical Proteomic Tumor Analysis Consortium (CPTAC) has performed phosphoproteomic profiling for 1,191 tumors spanning 11 cancer types, generating quantitative data on 77,442 phosphosites. Despite the impressive data generation power, only less than 5% of these phosphosites have been associated with a regulatory kinase or biological function, and 90% of these annotated phosphosites were associated with 20% well-studied kinases. The “dark” phosphoproteome greatly limits our ability to gain functional insights into cancer signaling. Previous research has shown that co-regulation is a strong indicator of functional association. Here we leveraged machine learning and the vast amount of CPTAC data to build a co-regulation map of phosphosites to facilitate functional interpretation of phosphoproteomics findings. Based on a ground-truth dataset including 98,402 pairs of phosphosites known to be regulated by the same kinase (positives) and 1,317,273 pairs by kinases from distant kinase families (negatives), we developed an Extreme Gradient Boosting (XGBoost) classifier to distinguish the positive and negative pairs using CPTAC phosphoproteomic data, protein-protein interaction data, and phosphopeptide sequences. Applying the trained classifier to 3 billion phosphosite pairs identified 2,569,519 (0.08%) with high probability of co-regulation, i.e., 400 times more likely to connect positive pairs than negative pairs in an independent ground-truth dataset. These pairs constituted a co-regulation map of 30,499 unique phosphosites, called CoPheeMap.To demonstrate the utility of CoPheeMap, we integrated its network embedding features, embedding features from a kinase network, together with Position-Specific Scoring Matrices scores and site-kinase abundance associations to develop a XGBoost model to predict kinase substrate associations (KSAs). The resulted model CoPheeKSA showed superior performance with an AUROC of 0.97 and identified 12,000 high-quality novel KSAs involving 7,908 phosphosites. In another application, CoPheeMap based information propagation (CoPheeProp) assigned 5,000 phosphosites to different signaling pathways with high specificity, increasing existing knowledge by 9-fold. Applying CoPheeMap and the derived tools to CPTAC and other cancer phosphoproteomics datasets revealed regulatory and functional information for previously unannotated sites showing strong associations with somatic mutations and cancer phenotypes, leading to actionable mechanistic and therapeutic insights. Together, CoPheeMap, CoPheeKSA and CoPheeProp provide a systematic framework to illuminate the dark phosphoproteome and have broad applications ranging from basic cancer biology research to clinical investigations. Citation Format: Wen Jiang, Eric Jaehnig, Bing Zhang. CoPheeMap: a co-regulation map of 30,000 phosphosite illuminates the dark cancer phosphoproteome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2023.

Computational Prediction of Interactions Between SARS-CoV-2 and Human Protein Pairs by PSSM-Based Images

Identifying protein-protein interactions is essential to predict the behavior of the virus and to design antiviral drugs against an infection. Like other viruses, SARS-CoV-2 virus must interact with a host cell in order to survive. Such interaction results in an infection in the host organism. Knowing which human protein interacts with the SARS-CoV-2 protein is an essential step in preventing viral infection. In silico approaches provide a reference for in vitro validation to protein-protein interaction studies by finding interacting protein pair candidates. The representation of proteins is one of the key steps for protein interaction network prediction. In this study, we proposed an image representation of proteins based on position-specific scoring matrices (PSSM). PSSMs are matrices that are obtained from multiple sequence alignments. In each of its cells, there is information about the probability of the occurrence of amino acids or nucleotides. PSSM matrices were handled as gray-scale images and called PSSM images. The main motivation of the study is to investigate whether these PSSM images are a suitable protein representation method. To determine adequate image size, conversion to grayscale images was performed at different sizes. SARS-CoV-2-human protein interaction network prediction based on image classification with siamese neural network and Resnet50 was performed on PSSM image datasets of different sizes. The accuracy results obtained with 200x200 size images and siamese neural network as 0.915, and with 400x400 size images and Resnet50 as 0.922 showed that PSSM images can be used for protein representation.

Prediction of Coreceptor Tropism in HIV-1 Subtype C in Botswana.

It remains unknown whether the C-C motif chemokine receptor type 5 (CCR5) coreceptor is still the predominant coreceptor used by Human Immunodeficiency Virus-1 (HIV-1) in Botswana, where the HIV-1 subtype C predominates. We sought to determine HIV-1C tropism in Botswana using genotypic tools, taking into account the effect of antiretroviral treatment (ART) and virologic suppression. HIV-1 gp120 V3 loop sequences from 5602 participants were analyzed for viral tropism using three coreceptor use predicting algorithms/tools: Geno2pheno, HIV-1C Web Position-Specific Score Matrices (WebPSSM) and the 11/25 charge rule. We then compared the demographic and clinical characteristics of people living with HIV (PLWH) harboring R5- versus X4-tropic viruses using χ2 and Wilcoxon rank sum tests for categorical and continuous data analysis, respectively. The three tools congruently predicted 64% of viruses as either R5-tropic or X4-tropic. Geno2pheno and the 11/25 charge rule had the highest concordance at 89%. We observed a significant difference in ART status between participants harboring X4- versus R5-tropic viruses. X4-tropic viruses were more frequent among PLWH receiving ART (χ2 test, p = 0.03). CCR5 is the predominant coreceptor used by HIV-1C strains circulating in Botswana, underlining the strong potential for CCR5 inhibitor use, even in PLWH with drug resistance. We suggest that the tools for coreceptor prediction should be used in combination.

A library of sensitive position-specific scoring matrices for high-throughput identification of nuclear pore complex subunits.

The nuclear pore complex exhibits different manifestations across eukaryotes, with certain components being restricted to specific clades. Several studies have been conducted to delineate the nuclear pore complex composition in various model organisms. Due to its pivotal role in cell viability, traditional lab experiments, such as gene knockdowns, can prove inconclusive and need to be complemented by a high-quality computational process. Here, using an extensive data collection, we create a robust library of nucleoporin protein sequences and their respective family-specific position-specific scoring matrices. By extensively validating each profile in different settings, we propose that the created profiles can be used to detect nucleoporins in proteomes with high sensitivity and specificity compared to existing methods. This library of profiles and the underlying sequence data can be used for the detection of nucleoporins in target proteomes.

Optimized Data Set and Feature Construction for Substrate Prediction of Membrane Transporters.

α-Helical transmembrane proteins termed membrane transporters mediate the passage of small hydrophilic substrate molecules across biological lipid bilayer membranes. Annotating the specific substrates of the dozens to hundreds of individual transporters of an organism is an important task. In the past, machine learning classifiers have been successfully trained on pan-organism data sets to predict putative substrates of transporters. Here, we critically examine the selection of an optimal data set of protein sequence features for the classification task. We focus on membrane transporters of the three model organisms Escherichia coli, Arabidopsis thaliana, and Saccharomyces cerevisiae, as well as human. We show that organism-specific classifiers can be robustly trained if at least 20 samples are available for each substrate class. If information from position-specific scoring matrices is included, such classifiers have F1 scores between 0.85 and 1.00. For the largest data set (A. thaliana), a 4-class classifier yielded an F-score of 0.97. On a pan-organism data set composed of transporters of all four organisms, amino acid and sugar transporters were predicted with an F1 score of 0.91.

Identification of potential driver mutations in glioblastoma using machine learning.

Glioblastoma is a fast and aggressively growing tumor in the brain and spinal cord. Mutation of amino acid residues in targets proteins, which are involved in glioblastoma, alters the structure and function and may lead to disease. In this study, we collected a set of 9386 disease-causing (drivers) mutations based on the recurrence in patient samples and experimentally annotated as pathogenic and 8728 as neutral (passenger) mutations. We observed that Arg is highly preferred at the mutant sites of drivers, whereas Met and Ile showed preferences in passengers. Inspecting neighboring residues at the mutant sites revealed that the motifs YP, CP and GRH, are preferred in drivers, whereas SI, IQ and TVI are dominant in neutral. In addition, we have computed other sequence-based features such as conservation scores, Position Specific Scoring Matrices (PSSM) and physicochemical properties, and developed a machine learning-based method, GBMDriver (GlioBlastoma Multiforme Drivers), for distinguishing between driver and passenger mutations. Our method showed an accuracy and AUC of 73.59% and 0.82, respectively, on 10-fold cross-validation and 81.99% and 0.87 in a blind set of 1809 mutants. The tool is available at https://web.iitm.ac.in/bioinfo2/GBMDriver/index.html. We envisage that the present method is helpful to prioritize driver mutations in glioblastoma and assist in identifying therapeutic targets.

A Novel Ensemble Learning-Based Computational Method to Predict Protein-Protein Interactions from Protein Primary Sequences.

Simple SummaryProtein–protein interactions (PPIs) play a central role in the evolution and progression of various biological processes. In this article, we constructed a novel ensemble-learning-based model to predict potential PPIs, which only utilized the protein sequence information. The presented method used Discrete Hilbert transform to extract amino acid sequence information from position-specific scoring matrices. Then these extracted features were fed into rotation forest for training and predicting. When applying our method to the three datasets (Yeast, Human, and Oryza sativa) for detecting PPIs, we obtained excellent prediction performance. Furthermore, the comparison results indicated that our computational model is effective and robust in predicting potential PPI pairs.Protein–protein interactions (PPIs) are crucial for understanding the cellular processes, including signal cascade, DNA transcription, metabolic cycles, and repair. In the past decade, a multitude of high-throughput methods have been introduced to detect PPIs. However, these techniques are time-consuming, laborious, and always suffer from high false negative rates. Therefore, there is a great need of new computational methods as a supplemental tool for PPIs prediction. In this article, we present a novel sequence-based model to predict PPIs that combines Discrete Hilbert transform (DHT) and Rotation Forest (RoF). This method contains three stages: firstly, the Position-Specific Scoring Matrices (PSSM) was adopted to transform the amino acid sequence into a PSSM matrix, which can contain rich information about protein evolution. Then, the 400-dimensional DHT descriptor was constructed for each protein pair. Finally, these feature descriptors were fed to the RoF classifier for identifying the potential PPI class. When exploring the proposed model on the Yeast, Human, and Oryza sativa PPIs datasets, it yielded excellent prediction accuracies of 91.93, 96.35, and 94.24%, respectively. In addition, we also conducted numerous experiments on cross-species PPIs datasets, and the predictive capacity of our method is also very excellent. To further access the prediction ability of the proposed approach, we present the comparison of RoF with four powerful classifiers, including Support Vector Machine (SVM), Random Forest (RF), K-nearest Neighbor (KNN), and AdaBoost. We also compared it with some existing superiority works. These comprehensive experimental results further confirm the excellent and feasibility of the proposed approach. In future work, we hope it can be a supplemental tool for the proteomics analysis.

MFPS_CNN: Multi‐filter Pattern Scanning from Position‐specific Scoring Matrix with Convolutional Neural Network for Efficient Prediction of Ion Transporters

In cellular transportation mechanisms, the movement of ions across the cell membrane and its proper control are important for cells, especially for life processes. Ion transporters/pumps and ion channel proteins work as border guards controlling the incessant traffic of ions across cell membranes. We revisited the study of classification of transporters and ion channels from membrane proteins with a more efficient deep learning approach. Specifically, we applied multi-window scanning filters of convolutional neural networks on almost full-length position-specific scoring matrices for extracting useful information. In this way, we were able to retain important evolutionary information of the proteins. Our experiment results show that a convolutional neural network with a minimum number of convolutional layers can be enough to extract the conserved information of proteins which leads to higher performance. Our best prediction models were obtained after examining different data imbalanced handling techniques, and different protein encoding methods. We also showed that our models were superior to traditional deep learning approaches on the same datasets as well as other machine learning classification algorithms.

Using multiple convolutional window scanning of convolutional neural network for an efficient prediction of ATP‐binding sites in transport proteins

Protein multiple sequence alignment information has long been important features to know about functions of proteins inferred from related sequences with known functions. It is therefore one of the underlying ideas of Alpha fold 2, a breakthrough study and model for the prediction of three-dimensional structures of proteins from their primary sequence. Our study used protein multiple sequence alignment information in the form of position-specific scoring matrices as input. We also refined the use of a convolutional neural network, a well-known deep-learning architecture with impressive achievement on image and image-like data. Specifically, we revisited the study of prediction of adenosine triphosphate (ATP)-binding sites with more efficient convolutional neural networks. We applied multiple convolutional window scanning filters of a convolutional neural network on position-specific scoring matrices for as much as useful information as possible. Furthermore, only the most specific motifs are retained at each feature map output through the one-max pooling layer before going to the next layer. We assumed that this way could help us retain the most conserved motifs which are discriminative information for prediction. Our experiment results show that a convolutional neural network with not too many convolutional layers can be enough to extract the conserved information of proteins, which leads to higher performance. Our best prediction models were obtained after examining them with different hyper-parameters. Our experiment results showed that our models were superior to traditional use of convolutional neural networks on the same datasets as well as other machine-learning classification algorithms.

Predicting protein phosphorylation sites in soybean using interpretable deep tabular learning network.

Phosphorylation of proteins is one of the most significant post-translational modifications (PTMs) and plays a crucial role in plant functionality due to its impact on signaling, gene expression, enzyme kinetics, protein stability and interactions. Accurate prediction of plant phosphorylation sites (p-sites) is vital as abnormal regulation of phosphorylation usually leads to plant diseases. However, current experimental methods for PTM prediction suffers from high-computational cost and are error-prone. The present study develops machine learning-based prediction techniques, including a high-performance interpretable deep tabular learning network (TabNet) to improve the prediction of protein p-sites in soybean. Moreover, we use a hybrid feature set of sequential-based features, physicochemical properties and position-specific scoring matrices to predict serine (Ser/S), threonine (Thr/T) and tyrosine (Tyr/Y) p-sites in soybean for the first time. The experimentally verified p-sites data of soybean proteins are collected from the eukaryotic phosphorylation sites database and database post-translational modification. We then remove the redundant set of positive and negative samples by dropping protein sequences with >40% similarity. It is found that the developed techniques perform >70% in terms of accuracy. The results demonstrate that the TabNet model is the best performing classifier using hybrid features and with window size of 13, resulted in 78.96 and 77.24% sensitivity and specificity, respectively. The results indicate that the TabNet method has advantages in terms of high-performance and interpretability. The proposed technique can automatically analyze the data without any measurement errors and any human intervention. Furthermore, it can be used to predict putative protein p-sites in plants effectively. The collected dataset and source code are publicly deposited at https://github.com/Elham-khalili/Soybean-P-sites-Prediction.

Interpreting Neural Networks for Biological Sequences by Learning Stochastic Masks.

Sequence-based neural networks can learn to make accurate predictions from large biological datasets, but model interpretation remains challenging. Many existing feature attribution methods are optimized for continuous rather than discrete input patterns and assess individual feature importance in isolation, making them ill-suited for interpreting non-linear interactions in molecular sequences. Building on work in computer vision and natural language processing, we developed an approach based on deep learning - Scrambler networks - wherein the most salient sequence positions are identified with learned input masks. Scramblers learn to predict Position-Specific Scoring Matrices (PSSMs) where unimportant nucleotides or residues are scrambled by raising their entropy. We apply Scramblers to interpret the effects of genetic variants, uncover non-linear interactions between cis-regulatory elements, explain binding specificity for protein-protein interactions, and identify structural determinants of de novo designed proteins. We show that Scramblers enable efficient attribution across large datasets and result in high-quality explanations, often outperforming state-of-the-art methods.

Robust and accurate prediction of protein\u2013protein interactions by exploiting evolutionary information

Various biochemical functions of organisms are performed by protein–protein interactions (PPIs). Therefore, recognition of protein–protein interactions is very important for understanding most life activities, such as DNA replication and transcription, protein synthesis and secretion, signal transduction and metabolism. Although high-throughput technology makes it possible to generate large-scale PPIs data, it requires expensive cost of both time and labor, and leave a risk of high false positive rate. In order to formulate a more ingenious solution, biology community is looking for computational methods to quickly and efficiently discover massive protein interaction data. In this paper, we propose a computational method for predicting PPIs based on a fresh idea of combining orthogonal locality preserving projections (OLPP) and rotation forest (RoF) models, using protein sequence information. Specifically, the protein sequence is first converted into position-specific scoring matrices (PSSMs) containing protein evolutionary information by using the Position-Specific Iterated Basic Local Alignment Search Tool (PSI-BLAST). Then we characterize a protein as a fixed length feature vector by applying OLPP to PSSMs. Finally, we train an RoF classifier for the purpose of identifying non-interacting and interacting protein pairs. The proposed method yielded a significantly better results than existing methods, with 90.07% and 96.09% prediction accuracy on Yeast and Human datasets. Our experiment show the proposed method can serve as a useful tool to accelerate the process of solving key problems in proteomics.

Master Blaster: an approach to sensitive identification of remotely related proteins

Genome sequencing projects unearth sequences of all the protein sequences encoded in a genome. As the first step, homology detection is employed to obtain clues to structure and function of these proteins. However, high evolutionary divergence between homologous proteins challenges our ability to detect distant relationships. In the past, an approach involving multiple Position Specific Scoring Matrices (PSSMs) was found to be more effective than traditional single PSSMs. Cascaded search is another successful approach where hits of a search are queried to detect more homologues. We propose a protocol, ‘Master Blaster’, which combines the principles adopted in these two approaches to enhance our ability to detect remote homologues even further. Assessment of the approach was performed using known relationships available in the SCOP70 database, and the results were compared against that of PSI-BLAST and HHblits, a hidden Markov model-based method. Compared to PSI-BLAST, Master Blaster resulted in 10% improvement with respect to detection of cross superfamily connections, nearly 35% improvement in cross family and more than 80% improvement in intra family connections. From the results it was observed that HHblits is more sensitive in detecting remote homologues compared to Master Blaster. However, there are true hits from 46-folds for which Master Blaster reported homologs that are not reported by HHblits even using the optimal parameters indicating that for detecting remote homologues, use of multiple methods employing a combination of different approaches can be more effective in detecting remote homologs. Master Blaster stand-alone code is available for download in the supplementary archive.

RNANet: an automatically built dual-source dataset integrating homologous sequences and RNA structures.

MotivationApplied research in machine learning progresses faster when a clean dataset is available and ready to use. Several datasets have been proposed and released over the years for specific tasks such as image classification, speech-recognition and more recently for protein structure prediction. However, for the fundamental problem of RNA structure prediction, information is spread between several databases depending on the level we are interested in: sequence, secondary structure, 3D structure or interactions with other macromolecules. In order to speed-up advances in machine-learning based approaches for RNA secondary and/or 3D structure prediction, a dataset integrating all this information is required, to avoid spending time on data gathering and cleaning.ResultsHere, we propose the first attempt of a standardized and automatically generated dataset dedicated to RNA combining together: RNA sequences, homology information (under the form of position-specific scoring matrices) and information derived by annotation of available 3D structures (including secondary structure, canonical and non-canonical interactions and backbone torsion angles). The data are retrieved from public databases PDB, Rfam and SILVA. The paper describes the procedure to build such dataset and the RNA structure descriptors we provide. Some statistical descriptions of the resulting dataset are also provided.Availability and implementationThe dataset is updated every month and available online (in flat-text file format) on the EvryRNA software platform (https://evryrna.ibisc.univ-evry.fr/evryrna/rnanet). An efficient parallel pipeline to build the dataset is also provided for easy reproduction or modification.Supplementary information Supplementary data are available at Bioinformatics online.

Deep profiling of protease substrate specificity enabled by dual random and scanned human proteome substrate phage libraries

Proteolysis is a major posttranslational regulator of biology inside and outside of cells. Broad identification of optimal cleavage sites and natural substrates of proteases is critical for drug discovery and to understand protease biology. Here, we present a method that employs two genetically encoded substrate phage display libraries coupled with next generation sequencing (SPD-NGS) that allows up to 10,000-fold deeper sequence coverage of the typical six- to eight-residue protease cleavage sites compared to state-of-the-art synthetic peptide libraries or proteomics. We applied SPD-NGS to two classes of proteases, the intracellular caspases, and the ectodomains of the sheddases, ADAMs 10 and 17. The first library (Lib 10AA) allowed us to identify 104 to 105 unique cleavage sites over a 1,000-fold dynamic range of NGS counts and produced consensus and optimal cleavage motifs based position-specific scoring matrices. A second SPD-NGS library (Lib hP), which displayed virtually the entire human proteome tiled in contiguous 49 amino acid sequences with 25 amino acid overlaps, enabled us to identify candidate human proteome sequences. We identified up to 104 natural linear cut sites, depending on the protease, and captured most of the examples previously identified by proteomics and predicted 10- to 100-fold more. Structural bioinformatics was used to facilitate the identification of candidate natural protein substrates. SPD-NGS is rapid, reproducible, simple to perform and analyze, inexpensive, and renewable, with unprecedented depth of coverage for substrate sequences, and is an important tool for protease biologists interested in protease specificity for specific assays and inhibitors and to facilitate identification of natural protein substrates.

Identification of Cancer Biomarkers in Human Body Fluids by Using Enhanced Physicochemical-incorporated Evolutionary Conservation Scheme.

Cancer is one of the most serious diseases affecting human health. Among all current cancer treatments, early diagnosis and control significantly help increase the chances of cure. Detecting cancer biomarkers in body fluids now is attracting more attention within oncologists. In-silico predictions of body fluid-related proteins, which can be served as cancer biomarkers, open a door for labor-intensive and time-consuming biochemical experiments. In this work, we propose a novel method for high-throughput identification of cancer biomarkers in human body fluids. We incorporate physicochemical properties into the weighted observed percentages (WOP) and position-specific scoring matrices (PSSM) profiles to enhance their attributes that reflect the evolutionary conservation of the body fluid-related proteins. The least absolute selection and shrinkage operator (LASSO) feature selection strategy is introduced to generate the optimal feature subset. The ten-fold cross-validation results on training datasets demonstrate the accuracy of the proposed model. We also test our proposed method on independent testing datasets and apply it to the identification of potential cancer biomarkers in human body fluids. The testing results promise a good generalization capability of our approach.