Search

Linear-B cell epitopes (LBCE) play a vital role in vaccine design; thus, efficiently detecting them from protein sequences is of primary importance. These epitopes consist of amino acids arranged in continuous or discontinuous patterns. Vaccines employ attenuated viruses and purified antigens. LBCE stimulate humoral immunity in the body, where B and T cells target circulating infections. To predict LBCE, the underlying protein sequences undergo a process of feature extraction, feature selection, and classification. Various system models have been proposed for this purpose, but their classification accuracy is only moderate. In order to enhance the accuracy of LBCE classification, this paper presents a novel 2-step metaheuristic variant-feature selection method that combines a linear support vector classifier (LSVC) with a Modified Genetic Algorithm (MGA). The feature selection model employs mono-peptide, dipeptide, and tripeptide features, focusing on the most diverse ones. These selected features are fed into a machine learning (ML)-based parallel ensemble classifier. The ensemble classifier combines correctly classified instances from various classifiers, including k-Nearest Neighbor (kNN), random forest (RF), logistic regression (LR), and support vector machine (SVM). The ensemble classifier came up with an impressively high accuracy of 99.3% as a result of its work. This accuracy is superior to the most recent models that are considered to be state-of-the-art for linear B-cell classification. As a direct consequence of this, the entire system model can now be utilised effectively in real-time clinical settings. [ABSTRACT FROM AUTHOR]

This paper outlines methods used to identify novel antigens for use in the development of serological assays. Specifically, we applied these methods to a neurogenic parasitic nematode of cervids called Parelaphostrongylus tenuis. This parasite is of particular concern in both wild and domestic ungulates as it causes significant neurological signs and definitive diagnosis is only possible post-mortem, necessitating the development of serologic assays for antemortem diagnosis. Proteins extracted from P. tenuis organisms were affinity isolated using antibodies enriched from seropositive moose (Alces alces). The proteins were analyzed using mass spectrometry and liquid chromatography to obtain amino acid sequences that were then cross-referenced to open reading frames predicted from an assembled transcriptome. An antigen of interest was assessed for immunogenic epitopes and subsequently synthesized into 10-mer synthetic overlapping peptides representing these regions. These synthetic peptides were then assessed for reactivity against positive and negative moose sera and demonstrated potential use as a serological assay in diagnostic laboratories. Known negative moose sera revealed significantly lower optical density when compared to the positive samples (p < 0.05). This method serves as a pipeline for the construction of diagnostic assays of pathogens in both human and veterinary medicine. [ABSTRACT FROM AUTHOR]

Protein structure prediction (PSP) has achieved significant progress lately via prediction of inter-residue distances using deep learning models and exploitation of the predictions during conformational search. In this context, prediction of large inter-residue distances and also prediction of distances between residues separated largely in the protein sequence remain challenging. To deal with these challenges, state-of-the-art inter-residue distance prediction algorithms have used large sets of coevolutionary and non-coevolutionary features. In this paper, we argue that the more the types of features used, the more the kinds of noises introduced and then the deep learning model has to overcome the noises to improve the accuracy of the predictions. Also, multiple features capturing similar underlying characteristics might not necessarily have significantly better cumulative effect. So we scrutinise the feature space to reduce the types of features to be used, but at the same time, we strive to improve the prediction accuracy. Consequently, for inter-residue real distance prediction, in this paper, we propose a deep learning model named scrutinised distance predictor (SDP), which uses only 2 coevolutionary and 3 non-coevolutionary features. On several sets of benchmark proteins, our proposed SDP method improves mean Local Distance Different Test (LDDT) scores at least by 10% over existing state-of-the-art methods. The SDP program along with its data is available from the website https://gitlab.com/mahnewton/sdp. [ABSTRACT FROM AUTHOR]

Jiangxi is one of the provinces in China most seriously affected by the haemorrhagic fever with renal syndrome (HFRS) epidemic. The aim of this paper was to systematically explore the HFRS epidemic in Jiangxi from the perspective of Hantavirus (HV) prevalence in rodents and humans and virus molecular characteristics. Individual information on all HFRS cases in Jiangxi from 2005 to 2021 was extracted from the China Information System for Disease Control and Prevention. All S and M fragment sequences of the Seoul virus and Hantan virus strains uploaded by Jiangxi and its neighbouring provinces and some representative sequences from provinces in China or some countries of Southeast Asia with the highest HV prevalence were retrieved and downloaded from NCBI GenBank. Periodogram and spatial autocorrelation were adopted for temporal periodicity and spatial clustering analysis of the HFRS epidemic. Joinpoint regression was utilized to explore the changing morbidity trend patterns of HFRS. Multiple sequence alignment and amino acid variation analysis were used to explore the homology and variation of strain prevalence in Jiangxi. Based on monthly morbidity time series, the periodogram analysis showed that the prevalence of HFRS had periodicities of 6 months and 12 months. Spatial autocorrelation analysis showed that HFRS distributed in Jiangxi was not random, with a "High-High" clustering area around Gaoan County. HFRS morbidity among the 0 ~ 15-year-old and ~ 61-year-old or older populations in Jiangxi increased significantly during the period of 2008–2015. Generally, HFRS morbidity was significantly positively correlated with the index of rat with virus (IRV) (r = 0.742) in the counties surrounding Gaoan from 2005 to 2019. HTNV strains in Jiangxi were in one independent branch, while the SEOV strains in Jiangxi were relatively more diverse. Both the YW89-15 and GAW30/2021 strains shared approximately 85% nucleotide homology and approximately 97% amino acid homology with their corresponding standard strains and vaccine strains. GAW30/2021 and YW89-15 had some amino acid site variations in nucleoprotein, glycoprotein precursor and RNA-dependent polymerase with their corresponding vaccine strains Z10 (HTNV) and Z37 (SEOV). The HFRS epidemic in Jiangxi has obvious temporal periodicity and spatial clustering, and the significant increase in the non-Immunization Expanded Program (EPI) targeted population (children and elderly) suggests that HFRS vaccination in this population needs to be considered. Although applying the EPI played a certain role in curbing the incidence of HFRS in Jiangxi from the perspective of ecological epidemiology, HTNV and SEOV strains prevalent in Jiangxi have some amino acid site variations compared to their corresponding vaccine strains, suggesting that HV variation needs to be continuously monitored in the future to observe vaccine protective efficiency. [ABSTRACT FROM AUTHOR]

Lysine malonylation is one of the most important post-translational modifications (PTMs). It affects the functionality of cells. Malonylation site prediction in proteins can unfold the mechanisms of cellular functionalities. Experimental methods are one of the due prediction approaches. But they are typically costly and time-consuming to implement. Recently, methods based on machine-learning solutions have been proposed to tackle this problem. Such practices have been shown to reduce costs and time complexities and increase accuracy. However, these approaches also have specific shortcomings, including inappropriate feature extraction out of protein sequences, high-dimensional features, and inefficient underlying classifiers. A machine learning-based method is proposed in this paper to cope with these problems. In the proposed approach, seven different features are extracted. Then, the extracted features are combined, ranked based on the Fisher's score (F-score), and the most efficient ones are selected. Afterward, malonylation sites are predicted using various classifiers. Simulation results show that the proposed method has acceptable performance compared with some state-of-the-art approaches. In addition, the XGBOOST classifier, founded on extracted features such as TFCRF, has a higher prediction rate than the other methods. The codes are publicly available at: https://github.com/jimy2020/Malonylation-site-prediction [ABSTRACT FROM AUTHOR]

Many pore-forming proteins originating from pathogenic bacteria are toxic against agricultural pests. They are the key ingredients in several pesticidal products for agricultural use, including transgenic crops. There is an urgent need to identify novel pore-forming proteins to combat development of resistance in pests to existing products, and to develop products that are effective against a broader range of pests. Existing computational methodologies to search for these proteins rely on sequence homology-based approaches. These approaches are based on similarities between protein sequences, and thus are limited in their usefulness for discovering novel proteins. In this paper, we outline a novel deep learning model trained on pore-forming proteins from the public domain. We compare different ways of encoding protein information during training, and contrast it with traditional approaches. We show that our model is capable of identifying known pore formers with no sequence similarity to the proteins used to train the model, and therefore holds promise for identifying novel pore formers. [ABSTRACT FROM AUTHOR]

The advent of the fourth industrial revolution, characterized by artificial intelligence (AI) as its central component, has resulted in the mechanization of numerous previously labor-intensive activities. The use of in silico tools has become prevalent in the design of biopharmaceuticals. Upon conducting a comprehensive analysis of the genomes of many organisms, it has been discovered that their tissues can generate specific peptides that confer protection against certain diseases. This study aims to identify a selected group of neuropeptides (NPs) possessing favorable characteristics that render them ideal for production as neurological biopharmaceuticals. Until now, the construction of NP classifiers has been the primary focus, neglecting to optimize these characteristics. Therefore, in this study, the task of creating ideal NPs has been formulated as a multi-objective optimization problem. The proposed framework, NPpred, comprises two distinct components: NSGA-NeuroPred and BERT-NeuroPred. The former employs the NSGA-II algorithm to explore and change a population of NPs, while the latter is an interpretable deep learning-based model. The utilization of explainable AI and motifs has led to the proposal of two novel operators, namely p-crossover and p-mutation. An online application has been deployed at https://neuropred.anvil.app for designing an ideal collection of synthesizable NPs from protein sequences. [ABSTRACT FROM AUTHOR]

Predicting the interaction affinity between drugs and target proteins is crucial for rapid and accurate drug discovery and repositioning. Therefore, more accurate prediction of DTA has become a key area of research in the field of drug discovery and drug repositioning. However, traditional experimental methods have disadvantages such as long operation cycles, high manpower requirements, and high economic costs, making it difficult to predict specific interactions between drugs and target proteins quickly and accurately. Some methods mainly use the SMILES sequence of drugs and the primary structure of proteins as inputs, ignoring the graph information such as bond encoding, degree centrality encoding, spatial encoding of drug molecule graphs, and the structural information of proteins such as secondary structure and accessible surface area. Moreover, previous methods were based on protein sequences to learn feature representations, neglecting the completeness of information. To address the completeness of drug and protein structure information, we propose a Transformer graph-based early fusion research approach for drug-target affinity prediction (GEFormerDTA). Our method reduces prediction errors caused by insufficient feature learning. Experimental results on Davis and KIBA datasets showed a better prediction of drugtarget affinity than existing affinity prediction methods. [ABSTRACT FROM AUTHOR]

Biosimilar insulin analogues are increasing market access for diabetic patients globally. Scientific establishment of biosimilarity is cornerstone of this key change in the medical landscape. BGL-ASP is a biosimilar insulin aspart developed by BioGenomics Limited, India. BioGenomics has considered a stepwise approach in generating the totality of evidence required to establish similarity with reference product. Insulin aspart is a recombinant rapid-acting human insulin analogue utilised in the treatment of type-1 and type-2 diabetes mellitus. The single amino acid substitution at position B28 where proline is replaced with aspartic acid results in a decreased propensity to form hexamers, thus increasing the absorption rate on subcutaneous administration compared to native insulin. In order to establish the safety and efficacy of BGL-ASP, the critical quality attributes (CQAs) of BGL-ASP are identified based on the impact created on biological activity, pharmacokinetic/pharmacodynamic (PK/PD), immunogenicity and safety. The CQAs of insulin aspart are related to product structure, purity and functionality and are characterised using a series of state-of-the-art orthogonal analytical tools. The primary protein sequence, the secondary, tertiary and quaternary structure are found to be highly similar for BGL-ASP and reference product. The product related impurities of insulin aspart and the assay content are determined using high performance liquid chromatography (HPLC) based analysis and is similar for BGL-ASP and reference insulin aspart sourced from United States of America (US), Europe Union (EU) and India. The safety, efficacy and immunogenicity of BGL-ASP is also found to be comparable with reference product and is confirmed through the clinical trials conducted as recommended by International Council for Harmonisation of Technical Requirements of Pharmaceuticals for Human Use (ICH) and European Medicines Agency (EMA) guidelines. The data encompassed in this study demonstrates that reference insulin aspart and BGL-ASP are highly similar in terms of structural, physicochemical, and biological properties, thus confirming its safety and efficacy for usage as potential alternative economical medicinal treatment for diabetes mellitus. [ABSTRACT FROM AUTHOR]

Protein–peptide interactions play a crucial role in various cellular processes and are implicated in abnormal cellular behaviors leading to diseases such as cancer. Therefore, understanding these interactions is vital for both functional genomics and drug discovery efforts. Despite a significant increase in the availability of protein–peptide complexes, experimental methods for studying these interactions remain laborious, time-consuming, and expensive. Computational methods offer a complementary approach but often fall short in terms of prediction accuracy. To address these challenges, we introduce PepCNN, a deep learning-based prediction model that incorporates structural and sequence-based information from primary protein sequences. By utilizing a combination of half-sphere exposure, position specific scoring matrices from multiple-sequence alignment tool, and embedding from a pre-trained protein language model, PepCNN outperforms state-of-the-art methods in terms of specificity, precision, and AUC. The PepCNN software and datasets are publicly available at https://github.com/abelavit/PepCNN.git. [ABSTRACT FROM AUTHOR]

Cerebrospinal fluid (CSF) leak can be diagnosed in clinical laboratories by detecting a diagnostic marker β2-transferrin (β2-Tf) in secretion samples. β2-Tf and the typical transferrin (Tf) proteoform in serum, β1-transferrin (β1-Tf), are Tf glycoforms. An innovative affinity capture technique for sample preparation, called microprobe-capture in-emitter elution (MPIE), was incorporated with high-resolution mass spectrometry (HR-MS) to study the Tf glycoforms and the primary structures of β1-Tf and β2-Tf. To implement MPIE, an analyte is first captured on the surface of a microprobe, and subsequently eluted from the microprobe inside an electrospray emitter. The capture process is monitored in real-time via next-generation biolayer interferometry (BLI). When electrospray is established from the emitter to a mass spectrometer, the analyte is immediately ionized via electrospray ionization (ESI) for HR-MS analysis. Serum, CSF, and secretion samples were analyzed using MPIE-ESI-MS. Based on the MPIE-ESI-MS results, the primary structures of β1-Tf and β2-Tf were elucidated. As Tf glycoforms, β1-Tf and β2-Tf share the amino acid sequence but contain varying N-glycans: (1) β1-Tf, the major serum-type Tf, has two G2S2 N-glycans on Asn413 and Asn611; and (2) β2-Tf, the major brain-type Tf, has an M5 N-glycan on Asn413 and a G0FB N-glycan on Asn611. The resolving power of the innovative MPIE-ESI-MS method was demonstrated in the study of β2-Tf as well as β1-Tf. Knowing the N-glycan structures on β2-Tf allows for the design of more novel test methods for β2-Tf in the future. [ABSTRACT FROM AUTHOR]

Protein therapeutics are susceptible to clipping via enzymatic and nonenzymatic mechanisms that create neo-N-termini. Typically, neo-N-termini are identified by chemical derivatization of the N-terminal amine with (N-Succinimidyloxycarbonylmethyl)tris(2,4,6-trimethoxyphenyl)phosphonium bromide (TMPP) followed by proteolysis and mass spectrometric analysis. Detection of the TMPP-labeled peptide is achieved by mapping the peptide sequence to the product ion spectrum derived from collisional activation. The site-specific localization of the TMPP tag enables unambiguous determination of the true N-terminus or neo-N-termini. In addition to backbone product ions, TMPP reporter ions at m/z 573, formed via collision-induced dissociation, can be diagnostic for the presence of a processed N-termini. However, reporter ions generated by collision-induced dissociation may be uninformative because of their low abundance. We demonstrate a novel high-throughput LC–MS method for the facile generation of the TMPP reporter ion at m/z 533 and, in some instances m/z 590, upon electron transfer dissociation. We further demonstrate the diagnostic utility of TMPP labeled peptides derived from a total cell lysate shows high degree of specificity towards selective N-terminal labeling over labeling of lysine and tyrosine and highly-diagnostic Receiver Operating Characteristic's (ROC) of TMPP reporter ions of m/z 533 and m/z 590. The abundant generation of these reporters enables subsequent MS/MS by intensity and m/z-dependent triggering of complementary ion activation modes such as collision-induced dissociation, high-energy collision dissociation, or ultraviolet photo dissociation for subsequent peptide sequencing. [ABSTRACT FROM AUTHOR]

High-throughput proteomic analysis of archaeological skeletal remains provides information about past fauna community compositions and species dispersals in time and space. Archaeological skeletal remains are a finite resource, however, and therefore it becomes relevant to optimize methods of skeletal proteome extraction. Ancient proteins in bone specimens can be highly degraded and consequently, extraction methods for well-preserved or modern bone might be unsuitable for the processing of highly degraded skeletal proteomes. In this study, we compared six proteomic extraction methods on Late Pleistocene remains with variable levels of proteome preservation. We tested the accuracy of species identification, protein sequence coverage, deamidation, and the number of post-translational modifications per method. We find striking differences in obtained proteome complexity and sequence coverage, highlighting that simple acid-insoluble proteome extraction methods perform better in highly degraded contexts. For well-preserved specimens, the approach using EDTA demineralization and protease-mix proteolysis yielded a higher number of identified peptides. The protocols presented here allowed protein extraction from ancient bone with a minimum number of working steps and equipment and yielded protein extracts within three working days. We expect further development along this route to benefit large-scale screening applications of relevance to archaeological and human evolution research. [ABSTRACT FROM AUTHOR]

Various approaches have used neural networks as probabilistic models for the design of protein sequences. These "inverse folding" models employ different objective functions, which come with trade-offs that have not been assessed in detail before. This study introduces probabilistic definitions of protein stability and conformational specificity and demonstrates the relationship between these chemical properties and the p (structure | seq) Boltzmann probability objective. This links the Boltzmann probability objective function to experimentally verifiable outcomes. We propose a novel sequence decoding algorithm, referred to as "BayesDesign", that leverages Bayes' Rule to maximize the p (structure | seq) objective instead of the p (seq | structure) objective common in inverse folding models. The efficacy of BayesDesign is evaluated in the context of two protein model systems, the NanoLuc enzyme and the WW structural motif. Both BayesDesign and the baseline ProteinMPNN algorithm increase the thermostability of NanoLuc and increase the conformational specificity of WW. The possible sources of error in the model are analyzed. [ABSTRACT FROM AUTHOR]

Phycobiliproteins (PBPs) are a class of water-soluble pigments with a variety of biological functions that are present in red macroalgae and cyanobacterial species. The crude forms of phycocyanin (C-PC) from the blue green alga Arthrospira platensis and allophycocyanin (APC) from the red macroalga Corallina officinalis were extracted and purified by ammonium sulphate precipitation, anion exchange chromatography, and size exclusion chromatography methods, respectively. The obtained C-PC and APC from A. platensis and C. officinalis were 0.31 mg/mL and 0.08 mg/mL, respectively, with molecular masses of "17.0 KDa and 19.0 KDa" and "15.0 KDa and 17.0 KDa" corresponding to α and β subunits, respectively. FT-IR was used to characterize the purified APC and C-PC in order to look into their structures. Highly purified extracts (A620/A280 > 4.0) were obtained from subtractions' PC3 and PC4 that were tested for their biological activities. APC and C-PC crude extracts plus their fractions exhibited potent anti-oxidant in different ratios by using three techniques. PC1 showed high anti-inflammatory (75.99 and 74.55%) and anti-arthritic (78.89 and 76.92%) activities for C. officinalis and A. platensis, respectively compared with standard drugs (72.02 and 71.5%). The methanolic and water extracts of both species showed greater antibacterial efficacy against Gram +ve than Gram −ve marine bacteria. Our study shed light on the potential medical uses of C-PC and APC extracted from the tested species as natural substances in a variety of foods and drugs. Further investigations are required to explore the diverse chemical natures of distinct PBPs from different cyanobacteria and red algae because their amino acid sequences vary among different algal species. [ABSTRACT FROM AUTHOR]

Antibodies with similar amino acid sequences, especially across their complementarity-determining regions, often share properties. Finding that an antibody of interest has a similar sequence to naturally expressed antibodies in healthy or diseased repertoires is a powerful approach for the prediction of antibody properties, such as immunogenicity or antigen specificity. However, as the number of available antibody sequences is now in the billions and continuing to grow, repertoire mining for similar sequences has become increasingly computationally expensive. Existing approaches are limited by either being low-throughput, non-exhaustive, not antibody specific, or only searching against entire chain sequences. Therefore, there is a need for a specialized tool, optimized for a rapid and exhaustive search of any antibody region against all known antibodies, to better utilize the full breadth of available repertoire sequences. We introduce Known Antibody Search (KA-Search), a tool that allows for the rapid search of billions of antibody variable domains by amino acid sequence identity across either the variable domain, the complementarity-determining regions, or a user defined antibody region. We show KA-Search in operation on the ∼ 2.4 billion antibody sequences available in the OAS database. KA-Search can be used to find the most similar sequences from OAS within 30 minutes and a representative subset of 10 million sequences in less than 9 seconds. We give examples of how KA-Search can be used to obtain new insights about an antibody of interest. KA-Search is freely available at https://github.com/oxpig/kasearch. [ABSTRACT FROM AUTHOR]

A novel betanucleorhabdovirus infecting Paris polyphylla var. yunnanensis, tentatively named Paris yunnanensis rhabdovirus 1 (PyRV1), was recently identified in Yunnan Province, China. The infected plants showed vein clearing and leaf crinkle at early stage of infection, followed by leaf yellowing and necrosis. Enveloped bacilliform particles were observed using electron microscopy. The virus was mechanically transmissible to Nicotiana bethamiana and N. glutinosa. The complete genome of PyRV1 consists of 13,509 nucleotides, the organization of which was typical of rhabdoviruses, containing six open reading frames encoding proteins N–P–P3–M–G–L on the anti-sense strand, separated by conserved intergenic regions and flanked by complementary 3′-leader and 5′-trailer sequences. The genome of PyRV1 shared highest nucleotide sequence identity (55.1%) with Sonchus yellow net virus (SYNV), and the N, P, P3, M, G, and L proteins showed 56.9%, 37.2%, 38.4%, 41.8%, 56.7%, and 49.4% amino acid sequence identities with respective proteins of SYNV, suggesting RyRV1 belongs to a new species of the genus Betanucleorhabdovirus. [ABSTRACT FROM AUTHOR]

Identification of protein–protein interactions (PPI) is among the critical problems in the domain of bioinformatics. Previous studies have utilized different AI-based models for PPI classification with advances in artificial intelligence (AI) techniques. The input to these models is the features extracted from different sources of protein information, mainly sequence-derived features. In this work, we present an AI-based PPI identification model utilizing a PPI network and protein sequences. The PPI network is represented as a graph where each node is a protein pair, and an edge is defined between two nodes if there exists a common protein between these nodes. Each node in a graph has a feature vector. In this work, we have used the language model to extract feature vectors directly from protein sequences. The feature vectors for protein in pairs are concatenated and used as a node feature vector of a PPI network graph. Finally, we have used the Graph-BERT model to encode the PPI network graph with sequence-based features and learn the hidden representation of the feature vector for each node. The next step involves feeding the learned representations of nodes to the fully connected layer, the output of which is fed into the softmax layer to classify the protein interactions. To assess the efficacy of the proposed PPI model, we have performed experiments on several PPI datasets. The experimental results demonstrate that the proposed approach surpasses the existing PPI works and designed baselines in classifying PPI. [ABSTRACT FROM AUTHOR]

The majority of the vitamin D that is present in the blood binds to vitamin D binding protein (VDBP) and circulates in the form of a complex (VDBP-Complex). Knowing the level of vitamin D in the body is crucial for vitamin D-related treatments so that the right dosage of vitamin D can be given. In other words, it is essential to distinguish between the protein VDBP and the complex form bound to vitamin D. As a novel way for the detection of VDBP-Complex, a more effective phage display methodology was applied in this study along with the addition of two approaches. In order to screen a sequence specific to the target only, the pre-binding method and after-binding method were performed. VDBP-Complex was directly coated on the petri dishes. In order to select phages that specifically bind to the VDBP-Complex, random phages were attached, and selected by 7 times of biopanning. Individual DNA sequences were analyzed for each biopanning to find specific peptide sequences for VDBP-Complex. The affinity of binding phages was verified by ELISA assay using an anti-M13 antibody. The phage having a sequence of SFTKTSTFTWRD (called as M3) has shown the highest binding affinity to VDBP-Complex. As a result of the removal test of VDBP-Complex using magnetic beads conjugated with M3 peptide, it was confirmed that significant decrease of VDBP-Complex. The unique characteristic of the M3 sequence was confirmed through a sequence-modified peptide (SFT motif). That is, it is expected that the M3 peptide may be used to determine the vitamin D levels in the blood. [ABSTRACT FROM AUTHOR]

Prokaryotic genomes evolve via horizontal gene transfer (HGT), mutations, and rearrangements. A noteworthy part of the HGT process is facilitated by genomic islands (GIs). While previous computational biology research has focused on developing tools to detect GIs in prokaryotic genomes, there has been little research investigating GI patterns and biological connections across species. We have pursued the novel idea of connecting GIs across prokaryotic and phage genomes via patterns of protein families. Such patterns are sequences of protein families frequently present in the genomes of multiple species. We combined the large data set from the IslandViewer4 database with protein families from Pfam while implementing a comprehensive strategy to identify patterns making use of HMMER, BLAST, and MUSCLE. we also implemented Python programs that link the analysis into a single pipeline. Research results demonstrated that related GIs often exist in species that are evolutionarily unrelated and in multiple bacterial phyla. Analysis of the discovered patterns led to the identification of biological connections among prokaryotes and phages. These connections suggest broad HGT connections across the bacterial kingdom and its associated phages. The discovered patterns and connections could provide the basis for additional analysis on HGT breadth and the patterns in pathogenic GIs. [ABSTRACT FROM AUTHOR]

Cancer has become a major factor threatening human life and health. Under the circumstance that traditional treatment methods such as chemotherapy and radiotherapy are not highly specific and often cause severe side effects and toxicity, new treatment methods are urgently needed. Anticancer peptide drugs have low toxicity, stronger efficacy and specificity, and have emerged as a new type of cancer treatment drugs. However, experimental identification of anticancer peptides is time-consuming and expensive, and difficult to perform in a high-throughput manner. Computational identification of anticancer peptides can make up for the shortcomings of experimental identification. In this study, a deep learning-based predictor named ACPred-BMF is proposed for the prediction of anticancer peptides. This method uses the quantitative and qualitative properties of amino acids, binary profile feature to numerical representation for the peptide sequences. The Bidirectional LSTM network architecture is used in the model, and the attention mechanism is also considered. To alleviate the black-box problem of deep learning model prediction, we visualized the automatically extracted features and used the Shapley additive explanations algorithm to determine the importance of features to further understand the anticancer peptide mechanism. The results show that our method is one of the state-of-the-art anticancer peptide predictors. A web server as the implementation of ACPred-BMF that can be accessed via: http://mialab.ruc.edu.cn/ACPredBMFServer/. [ABSTRACT FROM AUTHOR]

Mouse monoclonal antibody M4M was recently designed to block human TRPM4 channel. The polypeptide for generating M4M is composed of peptide A1 between the transmembrane segment 5 (S5) and the pore, and a second peptide A2 between the pore and the transmembrane segment 6 (S6). Using peptide microarray, a 4-amino acid sequence EPGF within the A2 was identified to be the binding epitope for M4M. Substitution of EPGF with other amino acids greatly reduced binding affinity. Structural analysis of human TRPM4 structure indicates that EPGF is located externally to the channel pore. A1 is close to the EPGF binding epitope in space, albeit separated by a 37-amino acid peptide. Electrophysiological study reveals that M4M could block human TRPM4, but with no effect on rodent TRPM4 which shares a different amino acid sequence ERGS for the binding motif. Our results demonstrate that M4M is a specific inhibitor for human TRPM4. [ABSTRACT FROM AUTHOR]

Alignment-free methods for sequence comparison and phylogeny inference have attracted a great deal of attention in recent years. Several algorithms have been implemented in diverse software packages. Despite the great number of existing methods, most of them are based on word statistics. Although they propose different filtering and weighting strategies and explore different metrics, their performance may be limited by the phylogenetic signal preserved in these words. Herein, we present a different approach based on the species-specific amino acid neighborhood preferences. These differential preferences can be assessed in the context of vector spaces. In this way, a distance-based method to build phylogenies has been developed and implemented into an easy-to-use R package. Tests run on real-world datasets show that this method can reconstruct phylogenetic relationships with high accuracy, and often outperforms other alignment-free approaches. Furthermore, we present evidence that the new method can perform reliably on datasets formed by non-orthologous protein sequences, that is, the method not only does not require the identification of orthologous proteins, but also does not require their presence in the analyzed dataset. These results suggest that the neighborhood preference of amino acids conveys a phylogenetic signal that may be of great utility in phylogenomics. [ABSTRACT FROM AUTHOR]

The World Health Organization categorized SARS-CoV-2 as a variant of concern, having numerous mutations in spike protein, which have been found to evade the effect of antibodies stimulated by the COVID-19 vaccine. The susceptibility to omicron variant by immunization-induced antibodies are direly required for risk evaluation. To avoid the risk of arising viral illness, the construction of a specific vaccine that triggers the production of targeted antibodies to combat infection remains highly imperative. The aim of the present study is to develop a particular vaccine exploiting bioinformatics approaches which can target B- and T-cells epitopes. Through this approach, novel epitopes of the S protein-SARS-CoV-2 were predicted for the development of a multiple epitope vaccine. Multiple epitopes were selected on the basis of toxicity, immunogenicity and antigenicity, and vaccine subunit was constructed having potential immunogenic properties. The epitopes were linked with 3 types of linker EAAAK, AAY and GPGPG for vaccine construction. Subsequently, vaccine structure was docked with the receptor and cloned in a pET-28a (+) vector. The constructed vaccine was ligated in pET-28a (+) vector in E. coli using the SnapGene tool for the expression study and a good immune response was observed. Several computational tools were used to predict and analyze the vaccine constructed by using spike protein sequence of omicrons. The current study identified a Multi-Epitope Vaccine (MEV) as a significant vaccine candidate that could potentially help the global world to combat SARS-CoV-2 infections. [ABSTRACT FROM AUTHOR]

Accurate identification of drug-targets in human body has great significance for designing novel drugs. Compared with traditional experimental methods, prediction of drug-targets via machine learning algorithms has enhanced the attention of many researchers due to fast and accurate prediction. In this study, we propose a machine learning-based method, namely XGB-DrugPred for accurate prediction of druggable proteins. The features from primary protein sequences are extracted by group dipeptide composition, reduced amino acid alphabet, and novel encoder pseudo amino acid composition segmentation. To select the best feature set, eXtreme Gradient Boosting-recursive feature elimination is implemented. The best feature set is provided to eXtreme Gradient Boosting (XGB), Random Forest, and Extremely Randomized Tree classifiers for model training and prediction. The performance of these classifiers is evaluated by tenfold cross-validation. The empirical results show that XGB-based predictor achieves the best results compared with other classifiers and existing methods in the literature. [ABSTRACT FROM AUTHOR]

The long-necked woodwasp superfamily Xiphydrioidea belongs to the suborder Symphyta (Hymenoptera). Here we newly characterize the complete mitochondrial genome of the South Korean Euxiphydria potanini (Xiphydriidae) using next-generation sequencing: 16,500 bp long with 84.27% A + T content and 37 typical mitochondrial genes including those encoding 13 PCGs, 2 rRNAs, 22 tRNAs, and one A + T rich region. We compare the patterns of symphytan mitochondrial gene arrangement with those of an ancestral insect form and found some synapomorphic rearrangements in phylogenetic context. We use a variety of nucleotide and amino acid sequence alignments (thirteen mtPCGs and/or eight nDNAs) alongside step-by-step exclusions of long-branched taxa to elucidate the phylogenetic position of Xiphydrioidea and phylogenetic relationships among the seven symphytan superfamilies, except for Anaxyeloidea of which no mtgenome was available. The monophyly of symphytan superfamilies (with weak support for Pamphilioidea), sister-group relationship of Xiphydrioidea and Cephoidea, and Symphyta being paraphyletic to Apocrita, etc. are consistently supported by maximum likelihood and Bayesian inference trees. We also discuss the problematic phylogenetic positions of Orussoidea and Siricoidea and propose a hypothetical scenario of morphological character transition during hymenopteran evolution based on morphological key characteristics, such as the cenchrus and the wasp-waist. [ABSTRACT FROM AUTHOR]

Self-supervised language modeling is a rapidly developing approach for the analysis of protein sequence data. However, work in this area is heterogeneous and diverse, making comparison of models and methods difficult. Moreover, models are often evaluated only on one or two downstream tasks, making it unclear whether the models capture generally useful properties. We introduce the ProteinGLUE benchmark for the evaluation of protein representations: a set of seven per-amino-acid tasks for evaluating learned protein representations. We also offer reference code, and we provide two baseline models with hyperparameters specifically trained for these benchmarks. Pre-training was done on two tasks, masked symbol prediction and next sentence prediction. We show that pre-training yields higher performance on a variety of downstream tasks such as secondary structure and protein interaction interface prediction, compared to no pre-training. However, the larger base model does not outperform the smaller medium model. We expect the ProteinGLUE benchmark dataset introduced here, together with the two baseline pre-trained models and their performance evaluations, to be of great value to the field of protein sequence-based property prediction. Availability: code and datasets from https://github.com/ibivu/protein-glue. [ABSTRACT FROM AUTHOR]

Predicting the local structural features of a protein from its amino acid sequence helps its function prediction to be revealed and assists in three-dimensional structural modeling. As the sequence-structure gap increases, prediction methods have been developed to bridge this gap. Additionally, as the size of the structural database and computing power increase, the performance of these methods have also significantly improved. Herein, we present a powerful new tool called S-Pred, which can predict eight-state secondary structures (SS8), accessible surface areas (ASAs), and intrinsically disordered regions (IDRs) from a given sequence. For feature prediction, S-Pred uses multiple sequence alignment (MSA) of a query sequence as an input. The MSA input is converted to features by the MSA Transformer, which is a protein language model that uses an attention mechanism. A long short-term memory (LSTM) was employed to produce the final prediction. The performance of S-Pred was evaluated on several test sets, and the program consistently provided accurate predictions. The accuracy of the SS8 prediction was approximately 76%, and the Pearson's correlation between the experimental and predicted ASAs was 0.84. Additionally, an IDR could be accurately predicted with an F1-score of 0.514. The program is freely available at https://github.com/arontier/S%5fPred%5fPaper and https://ad3.io as a code and a web server. [ABSTRACT FROM AUTHOR]

Glycosylation of proteins and lipids occurs in vertebrates, usually terminating with sialylation, which regulates the physicochemical and biological properties of these glycoconjugates. Although less commonly known, sialic acid residues also undergo various modifications, such as acetylation, methylation, and sulfation. However, except for acetylation, the enzymes or functions of the other modification processes are unknown. To the best of our knowledge, this study is the first to demonstrate the ubiquitous occurrence of sulfated sialic acids and two genes encoding the sialate: O-sulfotransferases 1 and 2 in vertebrates. These two enzymes showed about 50% amino acid sequence identity, and appeared to be complementary to each other in acceptor substrate preferences. Gene targeting experiments showed that the deficiency of these genes was lethal for medaka fish during young fry development and accompanied by different phenotypes. Thus, the sulfation of sialic acids is essential for the vertebrate development. [ABSTRACT FROM AUTHOR]

The Axin family of scaffolding proteins control diverse processes, such as facilitating the interactions between cellular components and providing specificity to signaling pathways. While several Axin family members have been discovered in metazoans and shown to play crucial roles, their mechanism of action are not well understood. The Caenorhabditis elegans Axin homolog, pry-1, is a powerful tool for identifying interacting genes and downstream effectors that function in a conserved manner to regulate Axin-mediated signaling. Our lab and others have established pry-1's essential role in developmental processes that affect the reproductive system, seam cells, and a posterior P lineage cell, P11.p. Additionally, pry-1 is crucial for lipid metabolism, stress responses, and aging. In this study, we expanded on our previous work on pry-1 by reporting a novel interacting gene named picd-1 (pry-1-interacting and Cabin1 domain-containing). PICD-1 protein shares sequence conservation with CABIN1, a component of the HUCA complex. Our findings have revealed that PICD-1 is involved in several pry-1-mediated processes, including stress response and lifespan maintenance. picd-1's expression overlapped with that of pry-1 in multiple tissues throughout the lifespan. Furthermore, PRY-1 and PICD-1 inhibited CREB-regulated transcriptional coactivator homolog CRTC-1, which promotes longevity in a calcineurin-dependent manner. Overall, our study has demonstrated that picd-1 is necessary for mediating pry-1 function and provides the basis to investigate whether Cabin-1 domain-containing protein plays a similar role in Axin signaling in other systems. [ABSTRACT FROM AUTHOR]

AMPylation is an emerging post-translational modification that occurs on the hydroxyl group of threonine, serine, or tyrosine via a phosphodiester bond. AMPylators catalyze this process as covalent attachment of adenosine monophosphate to the amino acid side chain of a peptide. Recent studies have shown that this post-translational modification is directly responsible for the regulation of neurodevelopment and neurodegeneration and is also involved in many physiological processes. Despite the importance of this post-translational modification, there is no peptide sequence dataset available for conducting computation analysis. Therefore, so far, no computational approach has been proposed for predicting AMPylation. In this study, we introduce a new dataset of this distinct post-translational modification and develop a new machine learning tool using a deep convolutional neural network called DeepAmp to predict AMPylation sites in proteins. DeepAmp achieves 77.7%, 79.1%, 76.8%, 0.55, and 0.85 in terms of Accuracy, Sensitivity, Specificity, Matthews Correlation Coefficient, and Area Under Curve for AMPylation site prediction task, respectively. As the first machine learning model, DeepAmp demonstrate promising results which highlight its potential to solve this problem. Our presented dataset and DeepAmp as a standalone predictor are publicly available at https://github.com/MehediAzim/DeepAmp. [ABSTRACT FROM AUTHOR]

Bio-sequence comparators are one of the most basic and significant methods for assessing biological data, and so, due to the importance of proteins, protein sequence comparators are particularly crucial. On the other hand, the complexity of the problem, the growing number of extracted protein sequences, and the growth of studies and data analysis applications addressing protein sequences have necessitated the development of a rapid and accurate approach to account for the complexities in this field. As a result, we propose a protein sequence comparison approach, called PCV, which improves comparison accuracy by producing vectors that encode sequence data as well as physicochemical properties of the amino acids. At the same time, by partitioning the long protein sequences into fix-length blocks and providing encoding vector for each block, this method allows for parallel and fast implementation. To evaluate the performance of PCV, like other alignment-free methods, we used 12 benchmark datasets including classes with homologous sequences which may require a simple preprocessing search tool to select the homologous data. And then, we compared the protein sequence comparison outcomes to those of alternative alignment-based and alignment-free methods, using various evaluation criteria. These results indicate that our method provides significant improvement in sequence classification accuracy, compared to the alternative alignment-free methods and has an average correlation of about 94% with the ClustalW method as our reference method, while considerably reduces the processing time. [ABSTRACT FROM AUTHOR]

Protein protein interactions (PPI) are crucial for protein functioning, nevertheless predicting residues in PPI interfaces from the protein sequence remains a challenging problem. In addition, structure-based functional annotations, such as the PPI interface annotations, are scarce: only for about one-third of all protein structures residue-based PPI interface annotations are available. If we want to use a deep learning strategy, we have to overcome the problem of limited data availability. Here we use a multi-task learning strategy that can handle missing data. We start with the multi-task model architecture, and adapted it to carefully handle missing data in the cost function. As related learning tasks we include prediction of secondary structure, solvent accessibility, and buried residue. Our results show that the multi-task learning strategy significantly outperforms single task approaches. Moreover, only the multi-task strategy is able to effectively learn over a dataset extended with structural feature data, without additional PPI annotations. The multi-task setup becomes even more important, if the fraction of PPI annotations becomes very small: the multi-task learner trained on only one-eighth of the PPI annotations—with data extension—reaches the same performances as the single-task learner on all PPI annotations. Thus, we show that the multi-task learning strategy can be beneficial for a small training dataset where the protein's functional properties of interest are only partially annotated. [ABSTRACT FROM AUTHOR]

Understanding the function of microbial proteins is essential to reveal the clinical potential of the microbiome. The application of high-throughput sequencing technologies allows for fast and increasingly cheaper acquisition of data from microbial communities. However, many of the inferred protein sequences are novel and not catalogued, hence the possibility of predicting their function through conventional homology-based approaches is limited, which indicates the need for further research on alignment-free methods. Here, we leverage a deep-learning-based representation of proteins to assess its utility in alignment-free analysis of microbial proteins. We trained a language model on the Unified Human Gastrointestinal Protein catalogue and validated the resulting protein representation on the bacterial part of the SwissProt database. Finally, we present a use case on proteins involved in SCFA metabolism. Results indicate that the deep learning model manages to accurately represent features related to protein structure and function, allowing for alignment-free protein analyses. Technologies that contextualize metagenomic data are a promising direction to deeply understand the microbiome. [ABSTRACT FROM AUTHOR]

A newly designed TEMPO-FRIPS reagent, 4-(2,2,6,6-tetramethylpiperidine-1-oxyl) methyl benzyl succinic acid N-hydroxysuccinimide ester or p-TEMPO–Bn–Sc–NHS, was synthesized to achieve single-step free radical-initiated peptide sequencing mass spectrometry (FRIPS MS) for a number of model peptides, including phosphopeptides. The p-TEMPO–Bn–Sc–NHS reagent was conjugated to target peptides, and the resulting peptides were subjected to collisional activation. The peptide backbone dissociation behaviors of the MS/MS and MS3 experiments were monitored in positive ion mode. Fragment ions were observed even at the single-step thermal activation of the p-TEMPO–Bn–Sc–peptides, showing mainly a-/x- and c-/z-type fragments and neutral loss ions. This confirms that radical-driven peptide backbone dissociations occurred with the p-TEMPO–Bn–Sc–peptides. Compared to the previous version of the TEMPO reagent, i.e., o-TEMPO–Bz–C(O)–NHS, the newly designed p-TEMPO–Bn–Sc–NHS has better conjugation efficiency for the target peptides owing to its improved structural flexibility and solubility in the experimental reagents. An energetic interpretation using the survival fraction as a function of applied normalized collision energy (NCE) ascertained the difference in the thermal activation between p-TEMPO–Bn–Sc– and o-TEMPO–Bz–C(O)– radical initiators. This study clearly demonstrates that the application of the p-TEMPO–Bn–Sc– radical initiator can improve the duty cycle, and this FRIPS MS approach has the potential to be implemented in proteomics studies, including phosphoproteomics. [ABSTRACT FROM AUTHOR]

Proteins are the essential biological macromolecules required to perform nearly all biological processes, and cellular functions. Proteins rarely carry out their tasks in isolation but interact with other proteins (known as protein–protein interaction) present in their surroundings to complete biological activities. The knowledge of protein–protein interactions (PPIs) unravels the cellular behavior and its functionality. The computational methods automate the prediction of PPI and are less expensive than experimental methods in terms of resources and time. So far, most of the works on PPI have mainly focused on sequence information. Here, we use graph convolutional network (GCN) and graph attention network (GAT) to predict the interaction between proteins by utilizing protein's structural information and sequence features. We build the graphs of proteins from their PDB files, which contain 3D coordinates of atoms. The protein graph represents the amino acid network, also known as residue contact network, where each node is a residue. Two nodes are connected if they have a pair of atoms (one from each node) within the threshold distance. To extract the node/residue features, we use the protein language model. The input to the language model is the protein sequence, and the output is the feature vector for each amino acid of the underlying sequence. We validate the predictive capability of the proposed graph-based approach on two PPI datasets: Human and S. cerevisiae. Obtained results demonstrate the effectiveness of the proposed approach as it outperforms the previous leading methods. The source code for training and data to train the model are available at https://github.com/JhaKanchan15/PPI_GNN.git. [ABSTRACT FROM AUTHOR]

Protein–ligand interactions (PLIs) are essential for biochemical functionality and their identification is crucial for estimating biophysical properties for rational therapeutic design. Currently, experimental characterization of these properties is the most accurate method, however, this is very time-consuming and labor-intensive. A number of computational methods have been developed in this context but most of the existing PLI prediction heavily depends on 2D protein sequence data. Here, we present a novel parallel graph neural network (GNN) to integrate knowledge representation and reasoning for PLI prediction to perform deep learning guided by expert knowledge and informed by 3D structural data. We develop two distinct GNN architectures: GNN F is the base implementation that employs distinct featurization to enhance domain-awareness, while GNN P is a novel implementation that can predict with no prior knowledge of the intermolecular interactions. The comprehensive evaluation demonstrated that GNN can successfully capture the binary interactions between ligand and protein's 3D structure with 0.979 test accuracy for GNN F and 0.958 for GNN P for predicting activity of a protein–ligand complex. These models are further adapted for regression tasks to predict experimental binding affinities and pIC 50 crucial for compound's potency and efficacy. We achieve a Pearson correlation coefficient of 0.66 and 0.65 on experimental affinity and 0.50 and 0.51 on pIC 50 with GNN F and GNN P , respectively, outperforming similar 2D sequence based models. Our method can serve as an interpretable and explainable artificial intelligence (AI) tool for predicted activity, potency, and biophysical properties of lead candidates. To this end, we show the utility of GNN P on SARS-Cov-2 protein targets by screening a large compound library and comparing the prediction with the experimentally measured data. [ABSTRACT FROM AUTHOR]

Protein language models have emerged as an alternative to multiple sequence alignment for enriching sequence information and improving downstream prediction tasks such as biophysical, structural, and functional properties. Here we show that a method called SPOT-1D-LM combines traditional one-hot encoding with the embeddings from two different language models (ProtTrans and ESM-1b) for the input and yields a leap in accuracy over single-sequence-based techniques in predicting protein 1D secondary and tertiary structural properties, including backbone torsion angles, solvent accessibility and contact numbers for all six test sets (TEST2018, TEST2020, Neff1-2020, CASP12-FM, CASP13-FM and CASP14-FM). More significantly, it has a performance comparable to profile-based methods for those proteins with homologous sequences. For example, the accuracy for three-state secondary structure (SS3) prediction for TEST2018 and TEST2020 proteins are 86.7% and 79.8% by SPOT-1D-LM, compared to 74.3% and 73.4% by the single-sequence-based method SPOT-1D-Single and 86.2% and 80.5% by the profile-based method SPOT-1D, respectively. For proteins without homologous sequences (Neff1-2020) SS3 is 80.41% by SPOT-1D-LM which is 3.8% and 8.3% higher than SPOT-1D-Single and SPOT-1D, respectively. SPOT-1D-LM is expected to be useful for genome-wide analysis given its fast performance. Moreover, high-accuracy prediction of both secondary and tertiary structural properties such as backbone angles and solvent accessibility without sequence alignment suggests that highly accurate prediction of protein structures may be made without homologous sequences, the remaining obstacle in the post AlphaFold2 era. [ABSTRACT FROM AUTHOR]

Neuropeptides, as pervasive intercellular signaling molecules in the CNS, modulate a variety of behavioral systems in both protostomes and deuterostomes. Allatostatins are neuropeptides in arthropods that inhibit the biosynthesis of juvenile hormones. Based on amino acid sequences, they are divided into three different types in arthropods: allatostatin A, allatostatin B, allatostatin C. Allatostatin C (AstC) was first isolated from Manduca sexta, and it has an important conserved feature of a disulfide bridge formed by two cysteine residues. Moreover, AstC appears to be the ortholog of mammalian somatostatin, and it has functions in common with somatostatin, such as modulating feeding behaviors. The AstC signaling system has been widely studied in arthropods, but minimally studied in molluscs. In this study, we seek to identify the AstC signaling system in the marine mollusc Aplysia californica. We cloned the AstC precursor from the cDNA of Aplysia. We predicted a 15-amino acid peptide with a disulfide bridge, i.e., AstC, using NeuroPred. We then cloned two putative allatostatin C-like receptors and through NCBI Conserved Domain Search we found that they belonged to the G protein-coupled receptor (GPCR) family. In addition, using an inositol monophosphate 1 (IP1) accumulation assay, we showed that Aplysia AstC could activate one of the putative receptors, i.e., the AstC-R, at the lowest EC50, and AstC without the disulfide bridge (AstC') activated AstC-R with the highest EC50. Moreover, four molluscan AstCs with variations of sequences from Aplysia AstC but with the disulfide bridge activated AstC-R at intermediate EC50. In summary, our successful identification of the Aplysia AstC precursor and its receptor (AstC-R) represents the first example in molluscs, and provides an important basis for further studies of the AstC signaling system in Aplysia and other molluscs. [ABSTRACT FROM AUTHOR]

Low complexity regions (LCRs) are protein sequences formed by a set of compositionally biased residues. LCRs are extremely abundant in cellular proteins and have also been reported in viruses, where they may partake in evasion of the host immune system. Analyses of 28,231 SARS-CoV-2 whole proteomes and of 261,051 spike protein sequences revealed the presence of four extremely conserved LCRs in the spike protein of several SARS-CoV-2 variants. With the exception of Iota, where it is absent, the Spike LCR-1 is present in the signal peptide of 80.57% of the Delta variant sequences, and in other variants of concern and interest. The Spike LCR-2 is highly prevalent (79.87%) in Iota. Two distinctive LCRs are present in the Delta spike protein. The Delta Spike LCR-3 is present in 99.19% of the analyzed sequences, and the Delta Spike LCR-4 in 98.3% of the same set of proteins. These two LCRs are located in the furin cleavage site and HR1 domain, respectively, and may be considered hallmark traits of the Delta variant. The presence of the medically-important point mutations P681R and D950N in these LCRs, combined with the ubiquity of these regions in the highly contagious Delta variant opens the possibility that they may play a role in its rapid spread. [ABSTRACT FROM AUTHOR]

Despite the structural and functional information contained in the statistical coupling between pairs of residues in a protein, coevolution associated with function is often obscured by artifactual signals such as genetic drift, which shapes a protein's phylogenetic history and gives rise to concurrent variation between protein sequences that is not driven by selection for function. Here, we introduce a background model for phylogenetic contributions of statistical coupling that separates the coevolution signal due to inter-clade and intra-clade sequence comparisons and demonstrate that coevolution can be measured on multiple phylogenetic timescales within a single protein. Our method, nested coevolution (NC), can be applied as an extension to any coevolution metric. We use NC to demonstrate that poorly conserved residues can nonetheless have important roles in protein function. Moreover, NC improved the structural-contact predictions of several coevolution-based methods, particularly in subsampled alignments with fewer sequences. NC also lowered the noise in detecting functional sectors of collectively coevolving residues. Sectors of coevolving residues identified after application of NC were more spatially compact and phylogenetically distinct from the rest of the protein, and strongly enriched for mutations that disrupt protein activity. Thus, our conceptualization of the phylogenetic separation of coevolution provides the potential to further elucidate relationships among protein evolution, function, and genetic diseases. [ABSTRACT FROM AUTHOR]

Neuromyelitis optica spectrum disorder (NMOSD) is a relapsing autoimmune disease characterized by the presence of pathogenic autoantibodies, anti-aquaporin 4 (AQP4) antibodies. Recently, HLA-DQA1*05:03 was shown to be significantly associated with NMOSD in a Japanese patient cohort. However, the specific mechanism by which HLA-DQA1*05:03 is associated with the development of NMOSD has yet to be elucidated. In the current study, we revealed that HLA-DQA1*05:03 exhibited significantly higher cell surface expression levels compared to other various DQA1 alleles, and that its expression strongly depended on the amino acid sequence of the α1 domain, with a preference for leucine at position 75. Moreover, in silico analysis indicated that the HLA-DQ encoded by HLA-DQA1*05:03 preferentially presents immunodominant AQP4 peptides, and that the peptide major histocompatibility complexes (pMHCs) are more energetically stable in the presence of HLA-DQA1*05:03 than other HLA-DQA1 alleles. In silico 3D structural models were also applied to investigate the validity of the energetic stability of pMHCs. Taken together, our findings indicate that HLA-DQA1*05:03 possesses a distinct property to play a pathogenic role in the development of NMOSD. [ABSTRACT FROM AUTHOR]

Owing to their ability to maintain a thermodynamically stable fold at extremely high temperatures, thermophilic proteins (TTPs) play a critical role in basic research and a variety of applications in the food industry. As a result, the development of computation models for rapidly and accurately identifying novel TTPs from a large number of uncharacterized protein sequences is desirable. In spite of existing computational models that have already been developed for characterizing thermophilic proteins, their performance and interpretability remain unsatisfactory. We present a novel sequence-based thermophilic protein predictor, termed SCMTPP, for improving model predictability and interpretability. First, an up-to-date and high-quality dataset consisting of 1853 TPPs and 3233 non-TPPs was compiled from published literature. Second, the SCMTPP predictor was created by combining the scoring card method (SCM) with estimated propensity scores of g-gap dipeptides. Benchmarking experiments revealed that SCMTPP had a cross-validation accuracy of 0.883, which was comparable to that of a support vector machine-based predictor (0.906–0.910) and 2–17% higher than that of commonly used machine learning models. Furthermore, SCMTPP outperformed the state-of-the-art approach (ThermoPred) on the independent test dataset, with accuracy and MCC of 0.865 and 0.731, respectively. Finally, the SCMTPP-derived propensity scores were used to elucidate the critical physicochemical properties for protein thermostability enhancement. In terms of interpretability and generalizability, comparative results showed that SCMTPP was effective for identifying and characterizing TPPs. We had implemented the proposed predictor as a user-friendly online web server at http://pmlabstack.pythonanywhere.com/SCMTPP in order to allow easy access to the model. SCMTPP is expected to be a powerful tool for facilitating community-wide efforts to identify TPPs on a large scale and guiding experimental characterization of TPPs. [ABSTRACT FROM AUTHOR]

Red palm weevil [Rhynchophorus ferrugineus (Olivier)], is native to South Asia and expanding its distribution range globally. Recent invasions of red palm weevil around the world, including Saudi Arabia, has become a global constraint for the production of palm species. Although, several control measures have been tested, none of them seemed successful against this invasive species. Therefore, we focused on silencing the reproduction control gene vitellogenin (Vg) based on RNA interference (RNAi) strategy for its possible application in the management of R. ferrugineus. The Vg is a major yolk protein precursor critical for oogenesis. To do this, fat body transcriptome of R. ferrugineus female adults was sequenced, which provided partial Vg gene transcript (FPKM 5731.60). A complete RfVg gene transcript of 5504 bp encoding 1787 amino acids was then sequenced using RCAE-PCR strategy and characterized. Phylogenetic analysis suggested that RfVg has closer ancestry to the coleopteran insects. The RfVg-based RNAi significantly suppressed the expressions of Vg gene. The 15, 20 and 25 days post-injection periods suppressed Vg expressions by 95, 96.6 and 99%, respectively. The suppressed Vg expressions resulted in the dramatic failure of Vg protein expression, which caused atrophied ovaries or no oogenesis and ultimately eggs were not hatched. These results suggest that knockdown of Vg gene involved in R. ferrugineus reproduction is a promising target for RNAi-based management of R. ferrugineus. [ABSTRACT FROM AUTHOR]

The amino acid sequence of a protein contains all the necessary information to specify its shape, which dictates its biological activities. However, it is challenging and expensive to experimentally determine the three-dimensional structure of proteins. The backbone torsion angles play a critical role in protein structure prediction, and accurately predicting the angles can considerably advance the tertiary structure prediction by accelerating efficient sampling of the large conformational space for low energy structures. Here we first time propose evolutionary signatures computed from protein sequence profiles, and a novel recurrent architecture, termed ESIDEN, that adopts a straightforward architecture of recurrent neural networks with a small number of learnable parameters. The ESIDEN can capture efficient information from both the classic and new features benefiting from different recurrent architectures in processing information. On the other hand, compared to widely used classic features, the new features, especially the Ramachandran basin potential, provide statistical and evolutionary information to improve prediction accuracy. On four widely used benchmark datasets, the ESIDEN significantly improves the accuracy in predicting the torsion angles by comparison to the best-so-far methods. As demonstrated in the present study, the predicted angles can be used as structural constraints to accurately infer protein tertiary structures. Moreover, the proposed features would pave the way to improve machine learning-based methods in protein folding and structure prediction, as well as function prediction. The source code and data are available at the website https://kornmann.bioch.ox.ac.uk/leri/resources/download.html. [ABSTRACT FROM AUTHOR]

The development of an accurate and reliable variant effect prediction tool is important for research in human genetic diseases. A large number of predictors have been developed towards this goal, yet many of these predictors suffer from the problem of data circularity. Here we present MTBAN (Mutation effect predictor using the Temporal convolutional network and the Born-Again Networks), a method for predicting the deleteriousness of variants. We apply a form of knowledge distillation technique known as the Born-Again Networks (BAN) to a previously developed deep autoregressive generative model, mutationTCN, to achieve an improved performance in variant effect prediction. As the model is fully unsupervised and trained only on the evolutionarily related sequences of a protein, it does not suffer from the problem of data circularity which is common across supervised predictors. When evaluated on a test dataset consisting of deleterious and benign human protein variants, MTBAN shows an outstanding predictive ability compared to other well-known variant effect predictors. We also offer a user-friendly web server to predict variant effects using MTBAN, freely accessible at http://mtban.kaist.ac.kr. To our knowledge, MTBAN is the first variant effect prediction tool based on a deep generative model that provides a user-friendly web server for the prediction of deleteriousness of variants. [ABSTRACT FROM AUTHOR]

Invertebrate antimicrobial peptides (AMPs) are at the forefront in the search for agents of therapeutic utility against multi-resistant microbial pathogens, and in recent years substantial advances took place in the in silico prediction of antimicrobial function of amino acid sequences. A yet neglected aspect is taxonomic bias in the performance of these tools. Owing to differences in the prediction algorithms and used training data sets between tools, and phylogenetic differences in sequence diversity, physicochemical properties and evolved biological functions of AMPs between taxa, notable discrepancies may exist in performance between the currently available prediction tools. Here we tested if there is taxonomic bias in the prediction power in 10 tools with a total of 20 prediction algorithms in 19 invertebrate taxa, using a data set containing 1525 AMP and 3050 non-AMP sequences. We found that most of the tools exhibited considerable variation in performance between tested invertebrate groups. Based on the per-taxa performances and on the variation in performances across taxa we provide guidance in choosing the best-performing prediction tool for all assessed taxa, by listing the highest scoring tool for each of them. [ABSTRACT FROM AUTHOR]

C4 photosynthesis is a remarkable complex trait, elucidations of the evolutionary trajectory of C4 photosynthesis from its ancestral C3 pathway can help us better understand the generic principles of the evolution of complex traits and guide the engineering of C3 crops for higher yields. Here, we used the genus Flaveria that contains C3, C3–C4, C4-like and C4 species as a system to study the evolution of C4 photosynthesis. We first mapped transcript abundance, protein sequence and morphological features onto the phylogenetic tree of the genus Flaveria, and calculated the evolutionary correlation of different features; we then predicted the relative changes of ancestral nodes of those features to illustrate the major events during the evolution of C4 photosynthesis. We found that gene expression and protein sequence showed consistent modification patterns in the phylogenetic tree. High correlation coefficients ranging from 0.46 to 0.9 among gene expression, protein sequence and morphology were observed. The greatest modification of those different features consistently occurred at the transition between C3-C4 species and C4-like species. Our results show highly coordinated changes in gene expression, protein sequence and morphological features, which support evolutionary major events during the evolution of C4 metabolism. [ABSTRACT FROM AUTHOR]