BDB: biopanning data bank

Sequence homology score-based deep fuzzy network for identifying therapeutic peptides

The detection of therapeutic peptides is a topic of immense interest in the biomedical field. Conventional biochemical experiment-based detection techniques are tedious and time-consuming. Computational biology has become a useful tool for improving the detection efficiency of therapeutic peptides. Most computational methods do not consider the deviation caused by noise. To improve the generalization performance of therapeutic peptide prediction methods, this work presents a sequence homology score-based deep fuzzy echo-state network with maximizing mixture correntropy (SHS-DFESN-MMC) model. Our method is compared with the existing methods on eight types of therapeutic peptide datasets. The model parameters are determined by 10 fold cross-validation on their training sets and verified by independent test sets. Across the 8 datasets, the average area under the receiver operating characteristic curve (AUC) values of SHS-DFESN-MMC are the highest on both the training (0.926) and independent sets (0.923).

Phage Display Panning on Silica: Optimization of Elution Conditions for Selection of Strong Binders

Phage display panning is a powerful tool to select strong peptide binders to a given target, and when applied to inorganic materials (e.g., silica) as a target, it provides information on binding events and molecular recognition at the peptide-mineral interface. The panning process has limitations with the phage chemical elution being affected by bias toward positively charged binders, resulting in the potential loss of information on binder diversity; the presence of fast growing phages with an intrinsic growth advantage; and the presence of false positives from target unrelated peptides. To overcome some of these limitations, we developed a panning approach based on the sequential use of different eluents (Gly-HCl, pH-2.2; MgCl2, pH-6.1; and TEA, pH-11.0), or pH conditions (Gly-HCl 2.2 < pH < 11.0) that allows the identification of a diverse and comprehensive pool of strong binders. We have assessed and tested the authenticity of the identified silica binders via a complementary experimental (in vivo phage recovery rates and TEM imaging) and bioinformatics approach. We provide experimental evidence of the nonspecificity of the Gly-HCl eluent as typically used. Using a fluorimetric assay, we investigate in vitro binding of two peptides that differ by pI-S4 (HYIDFRW, pI 7.80) and S5 (YSLKQYQ, pI 9.44)─modified at the C terminal with an amide group to simulate net charges in the phage display system, confirming the vital role of electrostatic interactions as driving binding forces in the phage panning process. The presented optimized phage panning approach provides an opportunity to match known surface interactions at play with suitable elution conditions; to select only sequences relevant to a particular interfacial system. The approach has the potential to open up avenues to design interfacial systems to advance our understanding of peptide-assisted mineral growth, among other possibilities.

Material-specific binding peptides empower sustainable innovations in plant health, biocatalysis, medicine and microplastic quantification

Material-binding peptides (MBPs) have emerged as a diverse and innovation-enabling class of peptides in applications such as plant-/human health, immobilization of catalysts, bioactive coatings, accelerated polymer degradation and analytics for micro-/nanoplastics quantification. Progress has been fuelled by recent advancements in protein engineering methodologies and advances in computational and analytical methodologies, which allow the design of, for instance, material-specific MBPs with fine-tuned binding strength for numerous demands in material science applications. A genetic or chemical conjugation of second (biological, chemical or physical property-changing) functionality to MBPs empowers the design of advanced (hybrid) materials, bioactive coatings and analytical tools. In this review, we provide a comprehensive overview comprising naturally occurring MBPs and their function in nature, binding properties of short man-made MBPs (<20 amino acids) mainly obtained from phage-display libraries, and medium-sized binding peptides (20-100 amino acids) that have been reported to bind to metals, polymers or other industrially produced materials. The goal of this review is to provide an in-depth understanding of molecular interactions between materials and material-specific binding peptides, and thereby empower the use of MBPs in material science applications. Protein engineering methodologies and selected examples to tailor MBPs toward applications in agriculture with a focus on plant health, biocatalysis, medicine and environmental monitoring serve as examples of the transformative power of MBPs for various industrial applications. An emphasis will be given to MBPs' role in detecting and quantifying microplastics in high throughput, distinguishing microplastics from other environmental particles, and thereby assisting to close an analytical gap in food safety and monitoring of environmental plastic pollution. In essence, this review aims to provide an overview among researchers from diverse disciplines in respect to material-(specific) binding of MBPs, protein engineering methodologies to tailor their properties to application demands, re-engineering for material science applications using MBPs, and thereby inspire researchers to employ MBPs in their research.

CD47Binder: Identify CD47 Binding Peptides by Combining Next-Generation Phage Display Data and Multiple Peptide Descriptors

CD47/SIRPα pathway is a new breakthrough in the field of tumor immunity after PD-1/PD-L1. While current monoclonal antibody therapies targeting CD47/SIRPα have demonstrated some anti-tumor effectiveness, there are several inherent limitations associated with these formulations. In the paper, we developed a predictive model that combines next-generation phage display (NGPD) and traditional machine learning methods to distinguish CD47 binding peptides. First, we utilized NGPD biopanning technology to screen CD47 binding peptides. Second, ten traditional machine learning methods based on multiple peptide descriptors and three deep learning methods were used to build computational models for identifying CD47 binding peptides. Finally, we proposed an integrated model based on support vector machine. During the five-fold cross-validation, the integrated predictor demonstrated specificity, accuracy, and sensitivity of 0.755, 0.764, and 0.772, respectively. Furthermore, an online bioinformatics tool called CD47Binder has been developed for the integrated predictor. This tool is readily accessible on http://i.uestc.edu.cn/CD47Binder/cgi-bin/CD47Binder.pl .

Systemic application of bone-targeting peptidoglycan hydrolases as a novel treatment approach for staphylococcal bone infection

IMPORTANCE

The rising prevalence of antimicrobial resistance in S. aureus has rendered treatment of staphylococcal infections increasingly difficult, making the discovery of alternative treatment options a high priority. Peptidoglycan hydrolases, a diverse group of bacteriolytic enzymes, show high promise as such alternatives due to their rapid and specific lysis of bacterial cells, independent of antibiotic resistance profiles. However, using these enzymes for the systemic treatment of local infections, such as osteomyelitis foci, needs improvement, as the therapeutic distributes throughout the whole host, resulting in low concentrations at the actual infection site. In addition, the occurrence of intracellularly persisting bacteria can lead to relapsing infections. Here, we describe an approach using tissue-targeting to increase the local concentration of therapeutic enzymes in the infected bone. The enzymes were modified with a short targeting moiety that mediated accumulation of the therapeutic in osteoblasts and additionally enables targeting of intracellularly surviving bacteria.

Identification and validation of Sertoli cell homing peptides as molecular steering for testis targeted drug delivery

The testicle, an organ privileged with immunity because of Blood-Testis Barrier (BTB), poses a major impediment to developing and delivering drugs to the testes. These problems can be prevented by targeting testicular cells using specific ligands, such as homing peptides. This is the first study to demonstrate the successful selection of Sertoli cell homing peptides using a phage display peptide library. The identification of peptides is performed with Sanger sequencing and high-throughput NGS. The Sertoli cell and testis targeting potential of the SCHP1 and SCHP2 was confirmed using confocal microscopy and flow cytometry of the FITC-labelled peptides and in vivo bio-distribution of the corresponding Cy5.5-tagged peptides. Secondary structures were predicted in the setting of different polarity by circular dichroism. The results suggest that SCHP1 and SCHP2 can effectively target Sertoli cells. In vivo bio-distribution in mouse models indicated significantly higher uptake of SCHP1 and SCHP2 by testes compared with the heart, brain, and spleen. SCHP1 and SCHP2 can be adopted as molecular steering for targeted male contraceptive delivery, treatment of testicular cancer, and male infertility. Further development of the peptides into peptidomimetics may increase their stability, and information on the molecular targets of these peptides may reveal their therapeutic potential.

PepQSAR: a comprehensive data source and information platform for peptide quantitative structure-activity relationships

Peptide quantitative structure-activity relationships (pQSARs) have been widely applied to the statistical modeling and empirical prediction of peptide activity, property and feature. In the procedure, the peptide structure is characterized at sequence level using amino acid descriptors (AADs) and then correlated with observations by machine learning methods (MLMs), consequently resulting in a variety of quantitative regression models used to explain the structural factors that govern peptide activities, to generalize peptide properties of unknown from known samples, and to design new peptides with desired features. In this study, we developed a comprehensive platform, termed PepQSAR database, which is a systematic collection and decomposition of various data sources and abundant information regarding the pQSARs, including AADs, MLMs, data sets, peptide sequences, measured activities, model statistics, and literatures. The database also provides a comparison function for the various previously built pQSAR models reported by different groups via distinct approaches. The structured and searchable PepQSAR database is expected to provide a useful resource and powerful tool for the computational peptidology community, which is freely available at http://i.uestc.edu.cn/PQsarDB .

Identification of a novel colon adenocarcinoma cell targeting peptide using phage display library biopanning

Phage display is well recognized as a promising high-throughput screening tool for the discovery of novel cancer-targeting peptides. Here, we screened a phage display library of 7-mer random peptides through in vitro biopanning to isolate peptide ligands binding to SW480 human colon adenocarcinoma cells. Three rounds of negative and positive selection caused a remarkable enrichment of colon cancer cell-binding phage clones with a significant enhancement of phage recovery efficiency (about 157-fold). A number of phage clones were picked out from the eluted phages of last selection round and sequenced. According to the results of cell binding assay and phage cell-based ELISA, one of the isolated peptides denoted as CCBP1 (with the sequence HAMRAQP) was indicated to have the highest binding efficiency, selectivity, and specificity toward colon cancer cells with no significant binding to control cells. Peptide competitive inhibition assay revealed that binding of the phage-displayed CCBP1 is competitively inhibited by the same free peptide, suggesting that CCBP1 specific binding to the target cell is independent of the phage context. Taken together, our findings provide support for the notion that CCBP1 binds specifically to colon cancer cells and might be a potential lead candidate for targeted delivery of imaging agents or therapeutic genes/drugs to colon tumors.

Analysis of Compositional Bias in a Commercial Phage Display Peptide Library by Next-Generation Sequencing

The principal presumption of phage display biopanning is that the naïve library contains an unbiased repertoire of peptides, and thus, the enriched variants derive from the affinity selection of an entirely random peptide pool. In the current study, we utilized deep sequencing to characterize the widely used Ph.DTM-12 phage display peptide library (New England Biolabs). The next-generation sequencing (NGS) data indicated the presence of stop codons and a high abundance of wild-type clones in the naïve library, which collectively result in a reduced effective size of the library. The analysis of the DNA sequence logo and global and position-specific frequency of amino acids demonstrated significant bias in the nucleotide and amino acid composition of the library inserts. Principal component analysis (PCA) uncovered the existence of four distinct clusters in the naïve library and the investigation of peptide frequency distribution revealed a broad range of unequal abundances for peptides. Taken together, our data provide strong evidence for the notion that the naïve library represents substantial departures from randomness at the nucleotide, amino acid, and peptide levels, though not undergoing any selective pressure for target binding. This non-uniform sequence representation arises from both the M13 phage biology and technical errors of the library construction. Our findings highlight the paramount importance of the qualitative assessment of the naïve phage display libraries prior to biopanning.

Peptide Shuttles for Blood–Brain Barrier Drug Delivery

The blood–brain barrier (BBB) limits the delivery of therapeutics to the brain but also represents the main gate for nutrient entrance. Targeting the natural transport mechanisms of the BBB offers an attractive route for brain drug delivery. Peptide shuttles are able to use these mechanisms to increase the transport of compounds that cannot cross the BBB unaided. As peptides are a group of biomolecules with unique physicochemical and structural properties, the field of peptide shuttles has substantially evolved in the last few years. In this review, we analyze the main classifications of BBB–peptide shuttles and the leading sources used to discover them.

PDL1Binder: Identifying programmed cell death ligand 1 binding peptides by incorporating next-generation phage display data and different peptide descriptors

Monoclonal antibody drugs targeting the PD-1/PD-L1 pathway have showed efficacy in the treatment of cancer patients, however, they have many intrinsic limitations and inevitable drawbacks. Peptide inhibitors as alternatives might compensate for the drawbacks of current PD-1/PD-L1 interaction blockers. Identifying PD-L1 binding peptides by random peptide library screening is a time-consuming and labor-intensive process. Machine learning-based computational models enable rapid discovery of peptide candidates targeting the PD-1/PD-L1 pathway. In this study, we first employed next-generation phage display (NGPD) biopanning to isolate PD-L1 binding peptides. Different peptide descriptors and feature selection methods as well as diverse machine learning methods were then incorporated to implement predictive models of PD-L1 binding. Finally, we proposed PDL1Binder, an ensemble computational model for efficiently obtaining PD-L1 binding peptides. Our results suggest that predictive models of PD-L1 binding can be learned from deep sequencing data and provide a new path to discover PD-L1 binding peptides. A web server was implemented for PDL1Binder, which is freely available at http://i.uestc.edu.cn/pdl1binder/cgi-bin/PDL1Binder.pl.

Restriction of the Global IgM Repertoire in Antiphospholipid Syndrome

The typical anti-phospholipid antibodies (APLA) in the anti-phospholipid syndrome (APS) are reactive with the phospholipid-binding protein β2GPI as well as a growing list of other protein targets. The relation of APLA to natural antibodies and the fuzzy set of autoantigens involved provoked us to study the changes in the IgM repertoire in APS. To this end, peptides selected by serum IgM from a 7-residue linear peptide phage display library (PDL) were deep sequenced. The analysis was aided by a novel formal representation of the Igome (the mimotope set reflecting the IgM specificities) in the form of a sequence graph. The study involved women with APLA and habitual abortions (n=24) compared to age-matched clinically healthy pregnant women (n=20). Their pooled Igomes (297 028 mimotope sequences) were compared also to the global public repertoire Igome of pooled donor plasma IgM (n=2 796 484) and a set of 7-mer sequences found in the J regions of human immunoglobulins (n=4 433 252). The pooled Igome was represented as a graph connecting the sequences as similar as the mimotopes of the same monoclonal antibody. The criterion was based on previously published data. In the resulting graph, identifiable clusters of vertices were considered related to the footprints of overlapping antibody cross-reactivities. A subgraph based on the clusters with a significant differential expression of APS patients’ mimotopes contained predominantly specificities underrepresented in APS. The differentially expressed IgM footprints showed also an increased cross-reactivity with immunoglobulin J regions. The specificities underexpressed in APS had a higher correlation with public specificities than those overexpressed. The APS associated specificities were strongly related also to the human peptidome with 1 072 mimotope sequences found in 7 519 human proteins. These regions were characterized by low complexity. Thus, the IgM repertoire of the APS patients was found to be characterized by a significant reduction of certain public specificities found in the healthy controls with targets representing low complexity linear self-epitopes homologous to human antibody J regions.

High-Throughput Monoclonal Antibody Discovery from Phage Libraries: Challenging the Current Preclinical Pipeline to Keep the Pace with the Increasing mAb Demand

Simple Summary

Monoclonal antibodies are increasingly used for a broad range of diseases. Rising demand must face with time time-consuming and laborious processes to isolate novel monoclonal antibodies. Next-generation sequencing coupled to phage display provides timely and sustainable high throughput selection strategy to rapidly access novel target. Here, we describe the current NGS-guided strategies to identify potential binders from enriched sub-libraires by applying a user-friendly informatic pipeline to identify and discard false positive clones. Rescue step and strategies to boost mAb yield are also discussed to improve the limiting selection and screening steps.

Abstract

Monoclonal antibodies are among the most powerful therapeutics in modern medicine. Since the approval of the first therapeutic antibody in 1986, monoclonal antibodies keep holding great expectations for application in a range of clinical indications, highlighting the need to provide timely and sustainable access to powerful screening options. However, their application in the past has been limited by time-consuming and expensive steps of discovery and production. The screening of antibody repertoires is a laborious step; however, the implementation of next-generation sequencing-guided screening of single-chain antibody fragments has now largely overcome this issue. This review provides a detailed overview of the current strategies for the identification of monoclonal antibodies from phage display-based libraries. We also discuss the challenges and the possible solutions to improve the limiting selection and screening steps, in order to keep pace with the increasing demand for monoclonal antibodies.

A Survey for Predicting ATP Binding Residues of Proteins Using Machine Learning Methods

Protein-ligand interactions are necessary for majority protein functions. Adenosine-5'-triphosphate (ATP) is one such ligand that plays vital role as a coenzyme in providing energy for cellular activities, catalyzing biological reaction and signaling. Knowing ATP binding residues of proteins is helpful for annotation of protein function and drug design. However, due to the huge amounts of protein sequences influx into databases in the post-genome era, experimentally identifying ATP binding residues is cost-ineffective and time-consuming. To address this problem, computational methods have been developed to predict ATP binding residues. In this review, we briefly summarized the application of machine learning methods in detecting ATP binding residues of proteins. We expect this review will be helpful for further research.

TUPDB: Target-Unrelated Peptide Data Bank

The isolation of target-unrelated peptides (TUPs) through biopanning remains as a major problem of phage display selection experiments. These TUPs do not have any actual affinity toward targets of interest, which tend to be mistakenly identified as target-binding peptides. Therefore, an information portal for storing TUP data is urgently needed. Here, we present a TUP data bank (TUPDB), which is a comprehensive, manually curated database of approximately 73 experimentally verified TUPs and 1963 potential TUPs collected from TUPScan, the BDB database, and public research articles. The TUPScan tool has been integrated in TUPDB to facilitate TUP analysis. We believe that TUPDB can help identify and remove TUPs in future reports in the biopanning community. The database is of great importance to improving the quality of phage display-based epitope mapping and promoting the development of vaccines, diagnostics, and therapeutics. The TUPDB database is available at http://i.uestc.edu.cn/tupdb .

Highly paralleled emulsion droplets for efficient isolation, amplification, and screening of cancer biomarker binding phages

Based on the linkage of genotype and phenotype, display technology has been widely used to generate specific ligands for profiling, imaging, diagnosis and therapy applications. However, due to the lack of effective monoclonal manipulation and affinity evaluation methods, traditional display technology has to undergo tedious steps of selection, clone isolation, amplification, sequencing, synthesis and characterization to obtain the binding sequences. To directly acquire high-affinity clones, we propose a double monoclonal display approach (dm-Display) for peptide screening based on highly paralleled monoclonal manipulation in emulsion droplets. dm-Display can monoclonally link the genotype, phenotype and affinity to realize integrated monoclonal separation, amplification, recognition and staining in one droplet so that discrete high-affinity clones can be quickly extracted. Monoclonal manipulations highly-parallelly occur in millions of droplets so that molecular screening of a highly diverse phage library is achieved. We have screened specific peptide ligands against CD71 and GPC1, proving the feasibility and generality of dm-Display. As a highly efficient ligand screening platform, dm-Display will promote the further development of molecular screening.

Streamlined phage display library protocols for identification of insect gut binding peptides highlight peptide specificity

Phage display libraries have been used to isolate insect gut binding peptides for use as pathogen transmission blocking agents, and to provide artificial anchors for increased toxicity of bacteria-derived pesticidal proteins. Previously, phage clones displaying enriched peptides were sequenced by Sanger sequencing. Here we present a streamlined protocol for identification of insect gut binding peptides, using insect-appropriate feeding strategies, with next generation sequencing and tailored bioinformatics analyses. The bioinformatics pipeline is designed to eliminate poorly enriched and false positive peptides, and to identify peptides predicted to be stable and hydrophilic. In addition to developing streamlined protocols, we also sought to address whether candidate gut binding peptides can bind to insects from more than one order, which is an important consideration for safe, practical use of peptide-modified pesticidal proteins. To this end, we screened phage display libraries for peptides that bind to the gut epithelia of two pest insects, the Asian citrus psyllid, Diaphorina citri (Hemiptera) and beet armyworm, Spodoptera exigua (Lepidoptera), and one beneficial insect, the western honey bee, Apis mellifera (Hymenoptera). While unique peptide sequences totaling 13,427 for D. citri, 89,561 for S. exigua and 69,053 for A. mellifera were identified from phage eluted from the surface of the insect guts, final candidate pools were comprised of 53, 107 and 1423 peptides respectively. The benefits of multiple rounds of biopanning, along with peptide binding properties in relation to practical use of peptide-modified pesticidal proteins for insect pest control are discussed.

Application of Next Generation Sequencing (NGS) in Phage Displayed Peptide Selection to Support the Identification of Arsenic-Binding Motifs

Next generation sequencing (NGS) in combination with phage surface display (PSD) are powerful tools in the newly equipped molecular biology toolbox for the identification of specific target binding biomolecules. Application of PSD led to the discovery of manifold ligands in clinical and material research. However, limitations of traditional phage display hinder the identification process. Growth-based library biases and target-unrelated peptides often result in the dominance of parasitic sequences and the collapse of library diversity. This study describes the effective enrichment of specific peptide motifs potentially binding to arsenic as proof-of-concept using the combination of PSD and NGS. Arsenic is an environmental toxin, which is applied in various semiconductors as gallium arsenide and selective recovery of this element is crucial for recycling and remediation. The development of biomolecules as specific arsenic-binding sorbents is a new approach for its recovery. Usage of NGS for all biopanning fractions allowed for evaluation of motif enrichment, in-depth insight into the selection process and the discrimination of biopanning artefacts, e.g., the amplification-induced library-wide reduction in hydrophobic amino acid proportion. Application of bioinformatics tools led to the identification of an SxHS and a carboxy-terminal QxQ motif, which are potentially involved in the binding of arsenic. To the best of our knowledge, this is the first report of PSD combined with NGS of all relevant biopanning fractions.

NeuroCS: A Tool to Predict Cleavage Sites of Neuropeptide Precursors

BACKGROUND

Neuropeptides are a class of bioactive peptides produced from neuropeptide precursors through a series of extremely complex processes, mediating neuronal regulations in many aspects. Accurate identification of cleavage sites of neuropeptide precursors is of great significance for the development of neuroscience and brain science.

OBJECTIVE

With the explosive growth of neuropeptide precursor data, it is pretty much needed to develop bioinformatics methods for predicting neuropeptide precursors' cleavage sites quickly and efficiently.

METHODS

We started with processing the neuropeptide precursor data from SwissProt and NueoPedia into two sets of data, training dataset and testing dataset. Subsequently, six feature extraction schemes were applied to generate different feature sets and then feature selection methods were used to find the optimal feature subset of each. Thereafter the support vector machine was utilized to build models for different feature types. Finally, the performance of models were evaluated with the independent testing dataset.

RESULTS

Six models are built through support vector machine. Among them the enhanced amino acid composition-based model reaches the highest accuracy of 91.60% in the 5-fold cross validation. When evaluated with independent testing dataset, it also showed an excellent performance with a high accuracy of 90.37% and Area under Receiver Operating Characteristic curve up to 0.9576.

CONCLUSION

The performance of the developed model was decent. Moreover, for users' convenience, an online web server called NeuroCS is built, which is freely available at http://i.uestc.edu.cn/NeuroCS/dist/index.html#/. NeuroCS can be used to predict neuropeptide precursors' cleavage sites effectively.

Novel anti-cancer candidates from a combinatorial peptide library

STAT3 is attractive target for development of anti-cancer therapeutics as it is implicated in nearly all forms of human tumors. To identify novel leads, we screened a combinatorial peptide library displayed on the surface of M13 bacteriophage. After three rounds of biopanning, a dodecapeptide with the YYVSWPPDMMHY sequence was found to be enriched by 36% while another with a short consensus motif was displayed in 20% of the phages. Binding analysis by isothermal titration calorimetry shows the most displayed peptide interacted with a K_d of 1.79 μM, which on modification of its structure to mimic the natural binding partners of STAT3 brought the affinity to high nanomolar range (K_d  = 500 nM). Using a panel of tumor cell lines, we show that the peptides prevented the proliferation of triple-negative breast cancer cells with a moderate activity (GI₅₀  = 50 μM). Furthermore, gene expression analysis shows the peptide reduced the expression of oncoproteins critical for tumor cell proliferation, angiogenesis, and metastasis. To find novel STAT3-interacting proteins, we searched the non-redundant sequences of the National Center for Biotechnology Information database which allowed us to identify potential binding partners of the protein. In sum, our data show the identified agents could serve as useful therapeutics candidates for further development.

Recent Advances in Computational Methods for Identifying Anticancer Peptides

Anticancer peptide (ACP) is a kind of small peptides that can kill cancer cells without damaging normal cells. In recent years, ACP has been pre-clinically used for cancer treatment. Therefore, accurate identification of ACPs will promote their clinical applications. In contrast to labor-intensive experimental techniques, a series of computational methods have been proposed for identifying ACPs. In this review, we briefly summarized the current progress in computational identification of ACPs. The challenges and future perspectives in developing reliable methods for identification of ACPs were also discussed. We anticipate that this review could provide novel insights into future researches on anticancer peptides.

A Survey for Predicting Enzyme Family Classes Using Machine Learning Methods

Enzymes are proteins that act as biological catalysts to speed up cellular biochemical processes. According to their main Enzyme Commission (EC) numbers, enzymes are divided into six categories: EC-1: oxidoreductase; EC-2: transferase; EC-3: hydrolase; EC-4: lyase; EC-5: isomerase and EC-6: synthetase. Different enzymes have different biological functions and acting objects. Therefore, knowing which family an enzyme belongs to can help infer its catalytic mechanism and provide information about the relevant biological function. With the large amount of protein sequences influxing into databanks in the post-genomics age, the annotation of the family for an enzyme is very important. Since the experimental methods are cost ineffective, bioinformatics tool will be a great help for accurately classifying the family of the enzymes. In this review, we summarized the application of machine learning methods in the prediction of enzyme family from different aspects. We hope that this review will provide insights and inspirations for the researches on enzyme family classification.

SSH: A Tool for Predicting Hydrophobic Interaction of Monoclonal Antibodies Using Sequences

Therapeutic antibodies are one of the most important parts of the pharmaceutical industry. They are widely used in treating various diseases such as autoimmune diseases, cancer, inflammation, and infectious diseases. Their development process however is often brought to a standstill or takes a longer time and is then more expensive due to their hydrophobicity problems. Hydrophobic interactions can cause problems on half-life, drug administration, and immunogenicity at all stages of antibody drug development. Some of the most widely accepted and used technologies for determining the hydrophobic interactions of antibodies include standup monolayer adsorption chromatography (SMAC), salt-gradient affinity-capture self-interaction nanoparticle spectroscopy (SGAC-SINS), and hydrophobic interaction chromatography (HIC). However, to measure SMAC, SGAC-SINS, and HIC for hundreds of antibody drug candidates is time-consuming and costly. To save time and money, a predictor called SSH is developed. Based on the antibody's sequence only, it can predict the hydrophobic interactions of monoclonal antibodies (mAbs). Using the leave-one-out crossvalidation, SSH achieved 91.226% accuracy, 96.396% sensitivity or recall, 84.196% specificity, 87.754% precision, 0.828 Mathew correlation coefficient (MCC), 0.919 f-score, and 0.961 area under the receiver operating characteristic (ROC) curve (AUC).

Development and Application of Computational Methods in Phage Display Technology

BACKGROUND

Phage display is a powerful and versatile technology for the identification of peptide ligands binding to multiple targets, which has been successfully employed in various fields, such as diagnostics and therapeutics, drug-delivery and material science. The integration of next generation sequencing technology with phage display makes this methodology more productive. With the widespread use of this technique and the fast accumulation of phage display data, databases for these data and computational methods have become an indispensable part in this community. This review aims to summarize and discuss recent progress in the development and application of computational methods in the field of phage display.

METHODS

We undertook a comprehensive search of bioinformatics resources and computational methods for phage display data via Google Scholar and PubMed. The methods and tools were further divided into different categories according to their uses.

RESULTS

We described seven special or relevant databases for phage display data, which provided an evidence-based source for phage display researchers to clean their biopanning results. These databases can identify and report possible target-unrelated peptides (TUPs), thereby excluding false-positive data from peptides obtained from phage display screening experiments. More than 20 computational methods for analyzing biopanning data were also reviewed. These methods were classified into computational methods for reporting TUPs, for predicting epitopes and for analyzing next generation phage display data.

CONCLUSION

The current bioinformatics archives, methods and tools reviewed here have benefitted the biopanning community. To develop better or new computational tools, some promising directions are also discussed.

AGONOTES: A Robot Annotator for Argonaute Proteins

The argonaute protein (Ago) exists in almost all organisms. In eukaryotes, it functions as a regulatory system for gene expression. In prokaryotes, it is a type of defense system against foreign invasive genomes. The Ago system has been engineered for gene silencing and genome editing and plays an important role in biological studies. With an increasing number of genomes and proteomes of various microbes becoming available, computational tools for identifying and annotating argonaute proteins are urgently needed. We introduce AGONOTES (Argonaute Notes). It is a web service especially designed for identifying and annotating Ago. AGONOTES uses the BLASTP similarity search algorithm to categorize all submitted proteins into three groups: prokaryotic argonaute protein (pAgo), eukaryotic argonaute protein (eAgo), and non-argonaute protein (non-Ago). Argonaute proteins can then be aligned to the corresponding standard set of Ago sequences using the multiple sequence alignment program MUSCLE. All functional domains of Ago can further be curated from the alignment results and visualized easily through Bio::Graphic modules in the BioPerl bundle. Compared with existing tools such as CD-Search and available databases such as UniProt and AGONOTES showed a much better performance on domain annotations, which is fundamental in studying the new Ago. AGONOTES can be freely accessed at http://i.uestc.edu.cn/agonotes/. AGONOTES is a friendly tool for annotating Ago domains from a proteome or a series of protein sequences.

From biomedicinal to in silico models and back to therapeutics: a review on the advancement of peptidic modeling

Bioactive peptides participate in numerous metabolic functions of living organisms and have emerged as potential therapeutics on a diverse range of diseases. Albeit peptide design does not go without challenges, overwhelming advancements on in silico methodologies have increased the scope of peptide-based drug design and discovery to an unprecedented amount. Within an in silico model versus an experimental validation scenario, this review aims to summarize and discuss how different in silico techniques contribute at present to the design of peptide-based molecules. Published in silico results from 2014 to 2018 were selected and discriminated in major methodological groups, allowing a transversal analysis, promoting a landscape vision and asserting its increasing value in drug design.

Combinatorial Avidity Selection of Mosaic Landscape Phages Targeted at Breast Cancer Cells-An Alternative Mechanism of Directed Molecular Evolution

Low performance of actively targeted nanomedicines required revision of the traditional drug targeting paradigm and stimulated the development of novel phage-programmed, self-navigating drug delivery vehicles. In the proposed smart vehicles, targeting peptides, selected from phage libraries using traditional principles of affinity selection, are substituted for phage proteins discovered through combinatorial avidity selection. Here, we substantiate the potential of combinatorial avidity selection using landscape phage in the discovery of Short Linear Motifs (SLiMs) and their partner domains. We proved an algorithm for analysis of phage populations evolved through multistage screening of landscape phage libraries against the MDA-MB-231 breast cancer cell line. The suggested combinatorial avidity selection model proposes a multistage accumulation of Elementary Binding Units (EBU), or Core Motifs (CorMs), in landscape phage fusion peptides, serving as evolutionary initiators for formation of SLiMs. Combinatorial selection has the potential to harness directed molecular evolution to create novel smart materials with diverse novel, emergent properties.

HMMCAS: A Web Tool for the Identification and Domain Annotations of CAS Proteins

The CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated proteins) adaptive immune systems are discovered in many bacteria and most archaea. These systems are encoded by cas (CRISPR-associated) operons that have an extremely diverse architecture. The most crucial step in the depiction of cas operons composition is the identification of cas genes or Cas proteins. With the continuous increase of the newly sequenced archaeal and bacterial genomes, the recognition of new Cas proteins is becoming possible, which not only provides candidates for novel genome editing tools but also helps to understand the prokaryotic immune system better. Here, we describe HMMCAS, a web service for the detection of CRISPR-associated structural and functional domains in protein sequences. HMMCAS uses hmmscan similarity search algorithm in HMMER3.1 to provide a fast, interactive service based on a comprehensive collection of hidden Markov models of Cas protein family. It can accurately identify the Cas proteins including those fusion proteins, for example the Cas1-Cas4 fusion protein in Candidatus Chloracidobacterium thermophilum B (Cab. thermophilum B). HMMCAS can also find putative cas operon and determine which type it belongs to. HMMCAS is freely available at http://i.uestc.edu.cn/hmmcas.

PhD7Faster 2.0: predicting clones propagating faster from the Ph.D.-7 phage display library by coupling PseAAC and tripeptide composition

Selection from phage display libraries empowers isolation of high-affinity ligands for various targets. However, this method also identifies propagation-related target-unrelated peptides (PrTUPs). These false positive hits appear because of their amplification advantages. In this report, we present PhD7Faster 2.0 for predicting fast-propagating clones from the Ph.D.-7 phage display library, which was developed based on the support vector machine. Feature selection was performed against PseAAC and tripeptide composition using the incremental feature selection method. Ten-fold cross-validation results show that PhD7Faster 2.0 succeeds a decent performance with the accuracy of 81.84%, the Matthews correlation coefficient of 0.64 and the area under the ROC curve of 0.90. The permutation test with 1,000 shuffles resulted in p < 0.001. We implemented PhD7Faster 2.0 into a publicly accessible web tool (http://i.uestc.edu.cn/sarotup3/cgi-bin/PhD7Faster.pl) and constructed standalone graphical user interface and command-line versions for different systems. The standalone PhD7Faster 2.0 is able to detect PrTUPs within small datasets as well as large-scale datasets. This makes PhD7Faster 2.0 an enhanced and powerful tool for scanning and reporting faster-growing clones from the Ph.D.-7 phage display library.

SAROTUP: a suite of tools for finding potential target-unrelated peptides from phage display data

SAROTUP (Scanner And Reporter Of Target-Unrelated Peptides) 3.1 is a significant upgrade to the widely used SAROTUP web server for the rapid identification of target-unrelated peptides (TUPs) in phage display data. At present, SAROTUP has gathered a suite of tools for finding potential TUPs and other purposes. Besides the TUPScan, the motif-based tool, and three tools based on the BDB database, i.e., MimoScan, MimoSearch, and MimoBlast, three predictors based on support vector machine, i.e., PhD7Faster, SABinder and PSBinder, are integrated into SAROTUP. The current version of SAROTUP contains 27 TUP motifs and 823 TUP sequences. We also developed the standalone SAROTUP application with graphical user interface (GUI) and command line versions for processing deep sequencing phage display data and distributed it as an open source package, which can perform perfectly locally on almost all systems that support C++ with little or no modification. The web interfaces of SAROTUP have also been redesigned to be more self-evident and user-friendly. The latest version of SAROTUP is freely available at http://i.uestc.edu.cn/sarotup3.

Molecular Design of Peptide-Fc Fusion Drugs

Background:

Peptide-Fc fusion drugs, also known as peptibodies, are a category of biological therapeutics in which the Fc region of an antibody is genetically fused to a peptide of interest. However, to develop such kind of drugs is laborious and expensive. Rational design is urgently needed.

Methods:

We summarized the key steps in peptide-Fc fusion technology and stressed the main computational resources, tools, and methods that had been used in the rational design of peptide-Fc fusion drugs. We also raised open questions about the computer-aided molecular design of peptide-Fc.

Results:

The design of peptibody consists of four steps. First, identify peptide leads from native ligands, biopanning, and computational design or prediction. Second, select the proper Fc region from different classes or subclasses of immunoglobulin. Third, fuse the peptide leads and Fc together properly. At last, evaluate the immunogenicity of the constructs. At each step, there are quite a few useful resources and computational tools.

Conclusion:

Reviewing the molecular design of peptibody will certainly help make the transition from peptide leads to drugs on the market quicker and cheaper.

Recent Advances on the Machine Learning Methods in Identifying DNA Replication Origins in Eukaryotic Genomics

The initiate site of DNA replication is called origins of replication (ORI) which is regulated by a set of regulatory proteins and plays important roles in the basic biochemical process during cell growth and division in all living organisms. Therefore, the study of ORIs is essential for understanding the cell-division cycle and gene expression regulation so that scholars can develop a new strategy against genetic diseases by using the knowledge of DNA replication. Thus, the accurate identification of ORIs will provide key clues for DNA replication research and clinical medicine. Although, the conventional experiments could provide accurate results, they are time-consuming and cost ineffective. On the contrary, bioinformatics-based methods can overcome these shortcomings. Especially, with the emergence of DNA sequences in the post-genomic era, it is highly expected to develop high throughput tools to identify ORIs based on sequence information. In this review, we will summarize the current progress in computational prediction of eukaryotic ORIs including the collection of benchmark dataset, the application of machine learning-based techniques, the results obtained by these methods, and the construction of web servers. Finally, we gave the future perspectives on ORIs prediction. The review provided readers with a whole background of ORIs prediction based on machine learning methods, which will be helpful for researchers to study DNA replication in-depth and drug therapy of genetic defect.

Identification of Antioxidant Proteins With Deep Learning From Sequence Information

Antioxidant proteins have been found closely linked to disease control for its ability to eliminate excess free radicals. Because of its medicinal value, the study of identifying antioxidant proteins is on the upsurge. Many machine-learning classifiers have performed poorly owing to the nonlinear and unbalanced nature of biological data. Recently, deep learning techniques showed advantages over many state-of-the-art machine learning methods in various fields. In this study, a deep learning based classifier was proposed to identify antioxidant proteins based on mixed g-gap dipeptide composition feature vector. The classifier employed deep autoencoder to extract nonlinear representation from raw input. The t-Distributed Stochastic Neighbor Embedding (t-SNE) was used for dimensionality reduction. Support vector machine was finally performed for classification. The classifier achieved F 1 score of 0.8842 and MCC of 0.7409 in 10-fold cross validation. Experimental results show that our proposed method outperformed the traditional machine learning methods and could be a promising tool for antioxidant protein identification. For the convenience of others' scientific research, we have developed a user-friendly web server called IDAod for antioxidant protein identification, which can be accessed freely at http://bigroup.uestc.edu.cn/IDAod/.

Special issue on Computational Resources and Methods in Biological Sciences

This special issue covers a wide range of topics in computational biology, such as database construction, sequence analysis and function prediction with machine learning methods, disease-related diagnosis, drug-target and drug discovery, and electronic health record system construction.

Comprehensive in-silico prediction of damage associated SNPs in Human Prolidase gene

Prolidase is cytosolic manganese dependent exopeptidase responsible for the catabolism of imido di and tripeptides. Prolidase levels have been associated with a number of diseases such as bipolar disorder, erectile dysfunction and varied cancers. Single nucleotide polymorphism present in coding region of proteins (nsSNPs) has the potential to alter the primary structure as well as function of the protein. Hence, it becomes necessary to differentiate the potential harmful nsSNPs from the neutral ones. 19 nsSNPs were predicted as damaging by in-silico analysis of 298 nsSNPs retrieved from dbSNP database. Consurf analysis showed 18 out of 19 substitutions were present in the conserved regions. 4 substitutions (D276N, D287N, E412K, and G448R) that observed to have damaging effect are present in catalytic pocket. Four SNPs listed in splice site region were found to affect splicing of mRNA by altering acceptor site. On 3′UTR scan of 77 SNPs listed in SNP database, 9 SNPs were lead to alter miRNA target sites. These results provide a filtered data to explore the effect of uncharacterized nsSNP and SNP related to UTRs and splice site of prolidase to find their association with the disease susceptibility and to design the target dependent drugs for therapeutics.

A Fast Projection-Based Algorithm for Clustering Big Data

With the fast development of various techniques, more and more data have been accumulated with the unique properties of large size (tall) and high dimension (wide). The era of big data is coming. How to understand and discover new knowledge from these data has attracted more and more scholars' attention and has become the most important task in data mining. As one of the most important techniques in data mining, clustering analysis, a kind of unsupervised learning, could group a set data into objectives(clusters) that are meaningful, useful, or both. Thus, the technique has played very important role in knowledge discovery in big data. However, when facing the large-sized and high-dimensional data, most of the current clustering methods exhibited poor computational efficiency and high requirement of computational source, which will prevent us from clarifying the intrinsic properties and discovering the new knowledge behind the data. Based on this consideration, we developed a powerful clustering method, called MUFOLD-CL. The principle of the method is to project the data points to the centroid, and then to measure the similarity between any two points by calculating their projections on the centroid. The proposed method could achieve linear time complexity with respect to the sample size. Comparison with K-Means method on very large data showed that our method could produce better accuracy and require less computational time, demonstrating that the MUFOLD-CL can serve as a valuable tool, at least may play a complementary role to other existing methods, for big data clustering. Further comparisons with state-of-the-art clustering methods on smaller datasets showed that our method was fastest and achieved comparable accuracy. For the convenience of most scholars, a free soft package was constructed.

RicyerDB: A Database For Collecting Rice Yield-related Genes with Biological Analysis

The Rice Yield-related Database (RicyerDB) was created to complement with related research of influence rice (Oryza sativa L.) yield in multiple traits by manually curating the related databases and literature, and genomics and proteomics information that could be useful for comprehensive understanding of the rice biology. RicyerDB provides a more valuable resource in which to efficiently investigate, browse and analyze yield-related genes. The whole data set can be easily queried and downloaded through the webpage. In addition, RicyerDB also constructed a protein-protein interaction network with biological analysis. The combined rice database opens a new path to facilitate researchers achieving information on rice gene in terms of their effects on traits important for rice breeding. The web server is freely available at: http://server.malab.cn/Ricyer/index.html.

iRSpot-Pse6NC: Identifying recombination spots in Saccharomyces cerevisiae by incorporating hexamer composition into general PseKNC

Meiotic recombination caused by meiotic double-strand DNA breaks. In some regions the frequency of DNA recombination is relatively higher, while in other regions the frequency is lower: the former is usually called "recombination hotspot", while the latter the "recombination coldspot". Information of the hot and cold spots may provide important clues for understanding the mechanism of genome revolution. Therefore, it is important to accurately predict these spots. In this study, we rebuilt the benchmark dataset by unifying its samples with a same length (131 bp). Based on such a foundation and using SVM (Support Vector Machine) classifier, a new predictor called "iRSpot-Pse6NC" was developed by incorporating the key hexamer features into the general PseKNC (Pseudo K-tuple Nucleotide Composition) via the binomial distribution approach. It has been observed via rigorous cross-validations that the proposed predictor is superior to its counterparts in overall accuracy, stability, sensitivity and specificity. For the convenience of most experimental scientists, the web-server for iRSpot-Pse6NC has been established at http://lin-group.cn/server/iRSpot-Pse6NC, by which users can easily obtain their desired result without the need to go through the detailed mathematical equations involved.

HBPred: a tool to identify growth hormone-binding proteins

Hormone-binding protein (HBP) is a kind of soluble carrier protein and can selectively and non-covalently interact with hormone. HBP plays an important role in life growth, but its function is still unclear. Correct recognition of HBPs is the first step to further study their function and understand their biological process. However, it is difficult to correctly recognize HBPs from more and more proteins through traditional biochemical experiments because of high experimental cost and long experimental period. To overcome these disadvantages, we designed a computational method for identifying HBPs accurately in the study. At first, we collected HBP data from UniProt to establish a high-quality benchmark dataset. Based on the dataset, the dipeptide composition was extracted from HBP residue sequences. In order to find out the optimal features to provide key clues for HBP identification, the analysis of various (ANOVA) was performed for feature ranking. The optimal features were selected through the incremental feature selection strategy. Subsequently, the features were inputted into support vector machine (SVM) for prediction model construction. Jackknife cross-validation results showed that 88.6% HBPs and 81.3% non-HBPs were correctly recognized, suggesting that our proposed model was powerful. This study provides a new strategy to identify HBPs. Moreover, based on the proposed model, we established a webserver called HBPred, which could be freely accessed at http://lin-group.cn/server/HBPred.

Research on folding diversity in statistical learning methods for RNA secondary structure prediction

How to improve the prediction accuracy of RNA secondary structure is currently a hot topic. The existing prediction methods for a single sequence do not fully consider the folding diversity which may occur among RNAs with different functions or sources. This paper explores the relationship between folding diversity and prediction accuracy, and puts forward a new method to improve the prediction accuracy of RNA secondary structure. Our research investigates the following: 1. The folding feature based on stochastic context-free grammar is proposed. By using dimension reduction and clustering techniques, some public data sets are analyzed. The results show that there is significant folding diversity among different RNA families. 2. To assign folding rules to RNAs without structural information, a classification method based on production probability is proposed. The experimental results show that the classification method proposed in this paper can effectively classify the RNAs of unknown structure. 3. Based on the existing prediction methods of statistical learning models, an RNA secondary structure prediction framework is proposed, namely "Cluster - Training - Parameter Selection - Prediction". The results show that, with information on folding diversity, prediction accuracy can be significantly improved.

Compositional Bias in Naïve and Chemically-modified Phage-Displayed Libraries uncovered by Paired-end Deep Sequencing

Understanding the composition of a genetically-encoded (GE) library is instrumental to the success of ligand discovery. In this manuscript, we investigate the bias in GE-libraries of linear, macrocyclic and chemically post-translationally modified (cPTM) tetrapeptides displayed on the M13KE platform, which are produced via trinucleotide cassette synthesis (19 codons) and NNK-randomized codon. Differential enrichment of synthetic DNA {S}, ligated vector {L} (extension and ligation of synthetic DNA into the vector), naïve libraries {N} (transformation of the ligated vector into the bacteria followed by expression of the library for 4.5 hours to yield a "naïve" library), and libraries chemically modified by aldehyde ligation and cysteine macrocyclization {M} characterized by paired-end deep sequencing, detected a significant drop in diversity in {L} → {N}, but only a minor compositional difference in {S} → {L} and {N} → {M}. Libraries expressed at the N-terminus of phage protein pIII censored positively charged amino acids Arg and Lys; libraries expressed between pIII domains N1 and N2 overcame Arg/Lys-censorship but introduced new bias towards Gly and Ser. Interrogation of biases arising from cPTM by aldehyde ligation and cysteine macrocyclization unveiled censorship of sequences with Ser/Phe. Analogous analysis can be used to explore library diversity in new display platforms and optimize cPTM of these libraries.

Biopanning data bank 2018: hugging next generation phage display

Database URL

The BDB database is available at http://immunet.cn/bdb.

PSBinder: A Web Service for Predicting Polystyrene Surface-Binding Peptides

Polystyrene surface-binding peptides (PSBPs) are useful as affinity tags to build a highly effective ELISA system. However, they are also a quite common type of target-unrelated peptides (TUPs) in the panning of phage-displayed random peptide library. As TUP, PSBP will mislead the analysis of panning results if not identified. Therefore, it is necessary to find a way to quickly and easily foretell if a peptide is likely to be a PSBP or not. In this paper, we describe PSBinder, a predictor based on SVM. To our knowledge, it is the first web server for predicting PSBP. The SVM model was built with the feature of optimized dipeptide composition and 87.02% (MCC = 0.74; AUC = 0.91) of peptides were correctly classified by fivefold cross-validation. PSBinder can be used to exclude highly possible PSBP from biopanning results or to find novel candidates for polystyrene affinity tags. Either way, it is valuable for biotechnology community.