Search

Convolutional neural networks (CNNs) are currently among the most widely-used deep neural network (DNN) architectures available and achieve state-of-the-art performance for many problems. Originally applied to computer vision tasks, CNNs work well with any data with a spatial relationship, besides images, and have been applied to different fields. However, recent works have highlighted numerical stability challenges in DNNs, which also relates to their known sensitivity to noise injection. These challenges can jeopardise their performance and reliability. This paper investigates DeepGOPlus, a CNN that predicts protein function. DeepGOPlus has achieved state-of-the-art performance and can successfully take advantage and annotate the abounding protein sequences emerging in proteomics. We determine the numerical stability of the model's inference stage by quantifying the numerical uncertainty resulting from perturbations of the underlying floating-point data. In addition, we explore the opportunity to use reduced-precision floating point formats for DeepGOPlus inference, to reduce memory consumption and latency. This is achieved by instrumenting DeepGOPlus' execution using Monte Carlo Arithmetic, a technique that experimentally quantifies floating point operation errors and VPREC, a tool that emulates results with customizable floating point precision formats. Focus is placed on the inference stage as it is the primary deliverable of the DeepGOPlus model, widely applicable across different environments. All in all, our results show that although the DeepGOPlus CNN is very stable numerically, it can only be selectively implemented with lower-precision floating-point formats. We conclude that predictions obtained from the pre-trained DeepGOPlus model are very reliable numerically, and use existing floating-point formats efficiently. [ABSTRACT FROM AUTHOR]

Protein hotspot residues are key sites that mediate protein-protein interactions. Accurate identification of these residues is essential for understanding the mechanism from protein to function and for designing drug targets. Current research has mostly focused on using machine learning methods to predict hot spots from known interface residues, which artificially extract the corresponding features of amino acid residues from sequence, structure, evolution, energy, and other information to train and test machine learning models. The process is cumbersome, time-consuming and laborious to some extent. This paper proposes a novel idea that develops a pre-trained protein sequence embedding model combined with a one-dimensional convolutional neural network, called Embed-1dCNN, to predict protein hotspot residues. In order to obtain large data samples, this work integrates and extracts data from the datasets of ASEdb, BID, SKEMPI and dbMPIKT to generate a new dataset, and adopts the SMOTE algorithm to expand positive samples to form the training set. The experimental results show that the method achieves an F1 score of 0.82 on the test set. Compared with other hot spot prediction methods, our model achieved better prediction performance. [ABSTRACT FROM AUTHOR]

Introduction: Emerging tick-borne viruses of medical and veterinary importance are increasingly being reported globally. This resurgence emphasizes the need for sustained surveillance to provide insights into tick-borne viral diversity and associated potential public health risks. We report on a virus tentatively designated Kinna virus (KIV) in the family Phenuiviridae and genus Bandavirus. The virus was isolated from a pool of Amblyomma gemma ticks from Kinna in Isiolo County, Kenya. High throughput sequencing of the virus isolate revealed close relatedness to the Guertu virus. The virus genome is consistent with the described genomes of other members of the genus Bandavirus, with nucleotides lengths of 6403, 3332 and 1752 in the Large (L), Medium (M) and Small (S) segments respectively. Phylogenetic analysis showed that the virus clustered with Guertu virus although it formed a distinct and well supported branch. The RdRp amino acid sequence had a 93.3% identity to that of Guertu virus, an indication that the virus is possibly novel. Neutralizing antibodies were detected in 125 (38.6%, 95% CI 33.3–44.1%) of the human sera from the communities in this region. In vivo experiments showed that the virus was lethal to mice with death occurring 6–9 days post-infection. The virus infected mammalian cells (Vero cells) but had reduced infectivity in the mosquito cell line (C636) tested. Conclusion: Isolation of this novel virus with the potential to cause disease in human and animal populations necessitates the need to evaluate its public health significance and contribution to disease burden in the affected regions. This also points to the need for continuous monitoring of vector and human populations in high-risk ecosystems to update pathogen diversity. [ABSTRACT FROM AUTHOR]

Drug discovery relies on predicting drug-target interaction (DTI), which is an important challenging task. The purpose of DTI is to identify the interaction between drug chemical compounds and protein targets. Traditional wet lab experiments are time-consuming and expensive, that's why in recent years, the use of computational methods based on machine learning has attracted the attention of many researchers. Actually, a dry lab environment focusing more on computational methods of interaction prediction can be helpful in limiting search space for wet lab experiments. In this paper, a novel multi-stage approach for DTI is proposed that called SRX-DTI. In the first stage, combination of various descriptors from protein sequences, and a FP2 fingerprint that is encoded from drug are extracted as feature vectors. A major challenge in this application is the imbalanced data due to the lack of known interactions, in this regard, in the second stage, the One-SVM-US technique is proposed to deal with this problem. Next, the FFS-RF algorithm, a forward feature selection algorithm, coupled with a random forest (RF) classifier is developed to maximize the predictive performance. This feature selection algorithm removes irrelevant features to obtain optimal features. Finally, balanced dataset with optimal features is given to the XGBoost classifier to identify DTIs. The experimental results demonstrate that our proposed approach SRX-DTI achieves higher performance than other existing methods in predicting DTIs. The datasets and source code are available at: https://github.com/Khojasteh-hb/SRX-DTI. [ABSTRACT FROM AUTHOR]

With the growing interest in using phages to combat antimicrobial resistance, computational methods for predicting phage-host interactions have been explored to help shortlist candidate phages. Most existing models consider entire proteomes and rely on manual feature engineering, which poses difficulty in selecting the most informative sequence properties to serve as input to the model. In this paper, we framed phage-host interaction prediction as a multiclass classification problem that takes as input the embeddings of a phage's receptor-binding proteins, which are known to be the key machinery for host recognition, and predicts the host genus. We explored different protein language models to automatically encode these protein sequences into dense embeddings without the need for additional alignment or structural information. We show that the use of embeddings of receptor-binding proteins presents improvements over handcrafted genomic and protein sequence features. The highest performance was obtained using the transformer-based protein language model ProtT5, resulting in a 3% to 4% increase in weighted F1 and recall scores across different prediction confidence thresholds, compared to using selected handcrafted sequence features. [ABSTRACT FROM AUTHOR]

Protein secondary structure prediction is extremely important for determining the spatial structure and function of proteins. In this paper, we apply an optimized convolutional neural network and long short-term memory neural network models to protein secondary structure prediction, which is called OCLSTM. We use an optimized convolutional neural network to extract local features between amino acid residues. Then use the bidirectional long short-term memory neural network to extract the remote interactions between the internal residues of the protein sequence to predict the protein structure. Experiments are performed on CASP10, CASP11, CASP12, CB513, and 25PDB datasets, and the good performance of 84.68%, 82.36%, 82.91%, 84.21% and 85.08% is achieved respectively. Experimental results show that the model can achieve better results. [ABSTRACT FROM AUTHOR]

Viral-host protein protein interaction (PPI) analysis is essential to decode the molecular mechanism of viral pathogen and host immunity processes which eventually help to control viral diseases and optimize therapeutics. The state-of-the-art viral-host PPI predictor leverages unsupervised embedding learning technique (doc2vec) to generate statistical representations of viral-host protein sequences and a Random Forest classifier for interaction prediction. However, doc2vec approach generates the statistical representations of viral-host protein sequences by merely modelling the local context of residues which only partially captures residue semantics. The paper in hand proposes a novel technique for generating better statistical representations of viral and host protein sequences based on the infusion of comprehensive local and global contextual information of the residues. While local residue context aware encoding captures semantic relatedness and short range dependencies of residues. Global residue context aware encoding captures comprehensive long-range residues dependencies, positional invariance of residues, and unique residue combination distribution important for interaction prediction. Using concatenated rich statistical representations of viral and host protein sequences, a robust machine learning framework "LGCA-VHPPI" is developed which makes use of a deep forest model to effectively model complex non-linearity of viral-host PPI sequences. An in-depth performance comparison of the proposed LGCA-VHPPI framework with existing diverse sequence encoding schemes based viral-host PPI predictors reveals that LGCA-VHPPI outperforms state-of-the-art predictor by 6%, 2%, and 2% in terms of matthews correlation coefficient over 3 different benchmark viral-host PPI prediction datasets. [ABSTRACT FROM AUTHOR]

Mammalian protein arginine methyltransferase 7 (PRMT7) has been shown to target substrates with motifs containing two arginine residues separated by one other residue (RXR motifs). In particular, the repression domain of human histone H2B (29-RKRSR-33) has been a key substrate in determining PRMT7 activity. We show that incubating human PRMT7 and [3H]-AdoMet with full-length Xenopus laevis histone H2B, containing the substitutions K30R and R31K (RKRSR to RRKSR), results in greatly reduced methylation activity. Using synthetic peptides, we have now focused on the enzymology behind this specificity. We show for the human and Xenopus peptide sequences 23–37 the difference in activity results from changes in the Vmax rather than the apparent binding affinity of the enzyme for the substrates. We then characterized six additional peptides containing a single arginine or a pair of arginine residues flanked by glycine and lysine residues. We have corroborated previous findings that peptides with an RXR motif have much higher activity than peptides that contain only one Arg residue. We show that these peptides have similar apparent km values but significant differences in their Vmax values. Finally, we have examined the effect of ionic strength on these peptides. We found the inclusion of salt had little effect on the Vmax value but a considerable increase in the apparent km value, suggesting that the inhibitory effect of ionic strength on PRMT7 activity occurs largely by decreasing apparent substrate-enzyme binding affinity. In summary, we find that even subtle substitutions in the RXR recognition motif can dramatically affect PRMT7 catalysis. [ABSTRACT FROM AUTHOR]

Background: Identification of DNA/Protein motifs is a crucial problem for biologists. Computational techniques could be of great help in this identification. In this direction, many computational models for motifs have been proposed in the literature. Methods: One such important model is the motif model. In this paper we describe a motif search web tool that predominantly employs this motif model. This web tool exploits the state-of-the art algorithms for solving the motif search problem. Results: The online tool has been helping scientists identify many unknown motifs. Many of our predictions have been successfully verified as well. We hope that this paper will expose this crucial tool to many more scientists. Availability and requirements: Project name: PMS - Panoptic Motif Search Tool. Project home page: http://pms.engr.uconn.edu or http://motifsearch.com. Licence: PMS tools will be readily available to any scientist wishing to use it for non-commercial purposes, without restrictions. The online tool is freely available without login. [ABSTRACT FROM AUTHOR]

In this paper, we use network approaches to analyze the relations between protein sequence features for the top hierarchical classes of CATH and SCOP. We use fundamental connectivity measures such as correlation (CR), normalized mutual information rate (nMIR), and transfer entropy (TE) to analyze the pairwise-relationships between the protein sequence features, and use centrality measures to analyze weighted networks constructed from the relationship matrices. In the centrality analysis, we find both commonalities and differences between the different protein 3D structural classes. Results show that all top hierarchical classes of CATH and SCOP present strong non-deterministic interactions for the composition and arrangement features of Cystine (C), Methionine (M), Tryptophan (W), and also for the arrangement features of Histidine (H). The different protein 3D structural classes present different preferences in terms of their centrality distributions and significant features. [ABSTRACT FROM AUTHOR]

Approximate string matching with k-differences has a number of practical applications, ranging from pattern recognition to computational biology. This paper proposes an efficient memory-access algorithm for parallel approximate string matching with k-differences on Graphics Processing Units (GPUs). In the proposed algorithm, all threads in the same GPUs warp share data using warp-shuffle operation instead of accessing the shared memory. Moreover, we implement the proposed algorithm by exploiting the memory structure of GPUs to optimize its performance. Experiment results for real DNA packages revealed that the performance of the proposed algorithm and its implementation archived up to 122.64 and 1.53 times compared to that of sequential algorithm on CPU and previous parallel approximate string matching algorithm on GPUs, respectively. [ABSTRACT FROM AUTHOR]

Background: Identifying conserved regions in protein sequences is a fundamental operation, occurring in numerous sequence-driven analysis pipelines. It is used as a way to decode domain-rich regions within proteins, to compute protein clusters, to annotate sequence function, and to compute evolutionary relationships among protein sequences. A number of approaches exist for identifying and characterizing protein families based on their domains, and because domains represent conserved portions of a protein sequence, the primary computation involved in protein family characterization is identification of such conserved regions. However, identifying conserved regions from large collections (millions) of protein sequences presents significant challenges. Methods: In this paper we present a new, alignment-free method for detecting conserved regions in protein sequences called NADDA (No-Alignment Domain Detection Algorithm). Our method exploits the abundance of exact matching short subsequences (k-mers) to quickly detect conserved regions, and the power of machine learning is used to improve the prediction accuracy of detection. We present a parallel implementation of NADDA using the MapReduce framework and show that our method is highly scalable. Results: We have compared NADDA with Pfam and InterPro databases. For known domains annotated by Pfam, accuracy is 83%, sensitivity 96%, and specificity 44%. For sequences with new domains not present in the training set an average accuracy of 63% is achieved when compared to Pfam. A boost in results in comparison with InterPro demonstrates the ability of NADDA to capture conserved regions beyond those present in Pfam. We have also compared NADDA with ADDA and MKDOM2, assuming Pfam as ground-truth. On average NADDA shows comparable accuracy, more balanced sensitivity and specificity, and being alignment-free, is significantly faster. Excluding the one-time cost of training, runtimes on a single processor were 49s, 10,566s, and 456s for NADDA, ADDA, and MKDOM2, respectively, for a data set comprised of approximately 2500 sequences. [ABSTRACT FROM AUTHOR]

Onchocerciasis (river blindness) is a debilitating tropical disease that causes significant eye and skin damage, afflicting millions worldwide. As global efforts shift from disease management to elimination, vaccines have become crucial supplementary tools. The Onchocerciasis Vaccine for Africa (TOVA) Initiative was established in 2015, to advance at least one vaccine candidate initially targeting onchocerciasis in infants and children below 5 years of age, through Phase I human trials by 2025. Notably, Ov-RAL-2 and Ov-103 antigens have shown great promise during pre-clinical development, however, the overall success rate of vaccine candidates during clinical development remains relatively low due to certain adverse effects and immunogenic limitations. This study, thus, aimed at predicting the safety and immunogenicity of Ov-RAL-2 and Ov-103 potential onchocerciasis vaccine candidates prior to clinical trials. Advanced molecular simulation models and analytical immunoinformatics algorithms were applied to predict potential adverse side effects and efficacy of these antigens in humans. The analyses revealed that both Ov-RAL-2 and Ov-103 demonstrate favourable safety profiles as toxicogenic and allergenic epitopes were found to be absent within each antigen. Also, both antigens were predicted to harbour substantial numbers of a wide range of distinct epitopes (antibodies, cytokines, and T- Cell epitopes) associated with protective immunity against onchocerciasis. In agreement, virtual vaccination simulation forecasted heightened, but sustained levels of primary and secondary protective immune responses to both vaccine candidates over time. Ov-103 was predicted to be non-camouflageable, as it lacked epitopes identical to protein sequences in the human proteome. Indeed, both antigens were able to bind with high affinity and activate the innate immune TLR4 receptor, implying efficient immune recognition. These findings suggest that Ov-RAL-2 and Ov-103 can induce sufficient protective responses through diverse humoral and cellular mechanisms. Overall, our study provides additional layer of evidence for advancing the clinical development of both vaccine candidates against onchocerciasis. [ABSTRACT FROM AUTHOR]

The American manatee (Trichechus manatus), experiencing population declines due to various threats, is the focus of conservation efforts that include the capture, rehabilitation, and release of orphaned calves when their mothers are unable to care for them. These efforts are compromised by the use of commercially available milk substitutes that lack essential components found in natural manatee breast milk, particularly immunoglobulin A (IgA). IgA plays a crucial role in nurturing the immune mucosal system and fostering a healthy microbiota. However, research on IgA in non-maternally fed manatees is limited due to the lack of species-specific reagents. To address this gap, our study employs immuno-informatics analysis to compare IgA sequences from manatees with those from other species, aiming to explore epitope similarity and sharing. We compared the protein sequence of manatee IgA with available IgA sequences, assessing similarity at the sequence, 3D structures, and epitope levels. Our findings reveal that human IgA exhibits the highest similarity in terms of sequence and 3D structure. Additionally, epitope analysis shows high conservation, identity, and similarity of predicted epitopes compared to human IgA. Future studies should focus on functional analysis using human IgA polyclonal reagents to detect manatee IgA in breast milk. Our findings highlight the potential of comparative analysis in advancing the understanding of immunology in non-human animals and overcoming challenges associated with the scarcity of species-specific reagents. [ABSTRACT FROM AUTHOR]

Glutathione transferase P1 (GSTP1) is a multi-functional protein that protects cells from electrophiles by catalyzing their conjugation with glutathione, and contributes to the regulation of cell proliferation, apoptosis, and signalling. GSTP1, usually described as a cytosolic enzyme, can localize to other cell compartments and we have reported its strong association with the plasma membrane. In the current study, the hypothesis that GSTP1 is palmitoylated and this modification facilitates its dynamic localization and function was investigated. Palmitoylation is the reversible post-translational addition of a 16-C saturated fatty acid to proteins, most commonly on Cys residues through a thioester bond. GSTP1 in MCF7 cells was modified by palmitate, however, GSTP1 Cys to Ser mutants (individual and Cys-less) retained palmitoylation. Treatment of palmitoylated GSTP1 with 0.1 N NaOH, which cleaves ester bonds, did not remove palmitate. Purified GSTP1 was spontaneously palmitoylated in vitro and peptide sequencing revealed that Cys48 and Cys102 undergo S-palmitoylation, while Lys103 undergoes the rare N-palmitoylation. N-palmitoylation occurs via a stable NaOH-resistant amide bond. Analysis of subcellular fractions of MCF7-GSTP1 cells and a modified proximity ligation assay revealed that palmitoylated GSTP1 was present not only in the membrane fraction but also in the cytosol. GSTP1 isolated from E. coli, and MCF7 cells (grown under fatty acid free or regular conditions), associated with plasma membrane-enriched fractions and this association was not altered by palmitoyl CoA. Overall, GSTP1 is modified by palmitate, at multiple sites, including at least one non-Cys residue. These modifications could contribute to regulating the diverse functions of GSTP1. [ABSTRACT FROM AUTHOR]

The IUCN Red List of Threatened Species contains 175 Brazilian bat species that are threatened by extinction in some degree. From this perspective, it is essential to expand the knowledge about the genetic diversity of vulnerable bats. Genomic sequencing can be useful to generate robust and informative genetic references, increasing resolution when analyzing relationships among populations, species, or higher taxonomic levels. In this study, we sequenced and characterized in detail the first complete mitochondrial genomes of Furipterus horrens, Lonchorhina aurita, and Natalus macrourus, and investigated their phylogenetic position based on amino acid sequences of protein-coding genes (PCGs). The mitogenomes of these species are 16,516, 16,697, and 16,668 bp in length, respectively, and each comprises 13 PCGs, 22 tRNA genes, two rRNA genes, and a putative control region (CR). In the three species, genes were arranged similarly to all other previously described bat mitogenomes, and nucleotide composition was also consistent with the reported range. The length and arrangement of rrnS and rrnL were also consistent with those of other bat species, showing a positive AT-skew and a negative GC-skew. Except for trnS1, for which we did not observe the DHU arm, all other tRNAs showed the cloverleaf secondary structure in the three species. In addition, the mitogenomes showed minor differences in start and stop codons, and in all PCGs, codons ending in adenine were more common compared to those ending in guanine. We found that PCGs of the three species use multiple codons to encode each amino acid, following the previously documented pattern. Furthermore, all PCGs are under purifying selection, with atp8 experiencing the most relaxed purifying selection. Considering the phylogenetic reconstruction, F. horrens was recovered as sister to Noctilio leporinus, L. aurita and Tonatia bidens shared a node within Phyllostomidae, and N. macrourus appeared as sister to Molossidae and Vespertilionidae. [ABSTRACT FROM AUTHOR]

In this paper, we describe SAFlex (Structural Alphabet Flexibility), an extension of an existing structural alphabet (HMM-SA), to better explore increasing protein three dimensional structure information by encoding conformations of proteins in case of missing residues or uncertainties. An SA aims to reduce three dimensional conformations of proteins as well as their analysis and comparison complexity by simplifying any conformation in a series of structural letters. Our methodology presents several novelties. Firstly, it can account for the encoding uncertainty by providing a wide range of encoding options: the maximum a posteriori, the marginal posterior distribution, and the effective number of letters at each given position. Secondly, our new algorithm deals with the missing data in the protein structure files (concerning more than 75% of the proteins from the Protein Data Bank) in a rigorous probabilistic framework. Thirdly, SAFlex is able to encode and to build a consensus encoding from different replicates of a single protein such as several homomer chains. This allows localizing structural differences between different chains and detecting structural variability, which is essential for protein flexibility identification. These improvements are illustrated on different proteins, such as the crystal structure of an eukaryotic small heat shock protein. They are promising to explore increasing protein redundancy data and obtain useful quantification of their flexibility., Competing Interests: The authors have declared that no competing interests exist.

DNA-binding proteins are fundamentally important in cellular processes. Several computational-based methods have been developed to improve the prediction of DNA-binding proteins in previous years. However, insufficient work has been done on the prediction of DNA-binding proteins from protein sequence information. In this paper, a novel predictor, DNABP (DNA-binding proteins), was designed to predict DNA-binding proteins using the random forest (RF) classifier with a hybrid feature. The hybrid feature contains two types of novel sequence features, which reflect information about the conservation of physicochemical properties of the amino acids, and the binding propensity of DNA-binding residues and non-binding propensities of non-binding residues. The comparisons with each feature demonstrated that these two novel features contributed most to the improvement in predictive ability. Furthermore, to improve the prediction performance of the DNABP model, feature selection using the minimum redundancy maximum relevance (mRMR) method combined with incremental feature selection (IFS) was carried out during the model construction. The results showed that the DNABP model could achieve 86.90% accuracy, 83.76% sensitivity, 90.03% specificity and a Matthews correlation coefficient of 0.727. High prediction accuracy and performance comparisons with previous research suggested that DNABP could be a useful approach to identify DNA-binding proteins from sequence information. The DNABP web server system is freely available at http://www.cbi.seu.edu.cn/DNABP/., Competing Interests: The authors have declared that no competing interests exist.

Celiac disease (CD) is an inflammatory disorder of the upper small intestine caused by the ingestion of storage proteins (prolamins and glutelins) from wheat, barley, rye, and, in rare cases, oats. CD patients need to follow a gluten-free diet by consuming gluten-free products with gluten contents of less than 20 mg/kg. Currently, the recommended method for the quantitative determination of gluten is an enzyme-linked immunosorbent assay (ELISA) based on the R5 monoclonal antibody. Because the R5 ELISA mostly detects the prolamin fraction of gluten, a new independent method is required to detect prolamins as well as glutelins. This paper presents the development of a method to quantitate 16 wheat marker peptides derived from all wheat gluten protein types by liquid chromatography tandem mass spectrometry (LC-MS/MS) in the multiple reaction monitoring mode. The quantitation of each marker peptide in the chymotryptic digest of a defined amount of the respective reference wheat protein type resulted in peptide-specific yields. This enabled the conversion of peptide into protein type concentrations. Gluten contents were expressed as sum of all determined protein type concentrations. This new method was applied to quantitate gluten in wheat starches and compared to R5 ELISA and gel-permeation high-performance liquid chromatography with fluorescence detection (GP-HPLC-FLD), which resulted in a strong correlation between LC-MS/MS and the other two methods. [ABSTRACT FROM AUTHOR]

DNA-binding proteins play pivotal roles in alternative splicing, RNA editing, methylating and many other biological functions for both eukaryotic and prokaryotic proteomes. Predicting the functions of these proteins from primary amino acids sequences is becoming one of the major challenges in functional annotations of genomes. Traditional prediction methods often devote themselves to extracting physiochemical features from sequences but ignoring motif information and location information between motifs. Meanwhile, the small scale of data volumes and large noises in training data result in lower accuracy and reliability of predictions. In this paper, we propose a deep learning based method to identify DNA-binding proteins from primary sequences alone. It utilizes two stages of convolutional neutral network to detect the function domains of protein sequences, and the long short-term memory neural network to identify their long term dependencies, an binary cross entropy to evaluate the quality of the neural networks. When the proposed method is tested with a realistic DNA binding protein dataset, it achieves a prediction accuracy of 94.2% at the Matthew’s correlation coefficient of 0.961. Compared with the LibSVM on the arabidopsis and yeast datasets via independent tests, the accuracy raises by 9% and 4% respectively. Comparative experiments using different feature extraction methods show that our model performs similar accuracy with the best of others, but its values of sensitivity, specificity and AUC increase by 27.83%, 1.31% and 16.21% respectively. Those results suggest that our method is a promising tool for identifying DNA-binding proteins. [ABSTRACT FROM AUTHOR]

This study presents the complete genome sequence of a novel nege-like virus identified in whiteflies (Bemisia tabaci MEAM1), provisionally designated as whitefly negevirus 1 (WfNgV1). The virus possesses a single-stranded RNA genome comprising 11,848 nucleotides, organized into four open reading frames (ORFs). These ORFs encode the putative RNA-dependent-RNA-polymerase (RdRp, ORF 1), a glycoprotein (ORF 2), a structural protein with homology to those in the SP24 family, (ORF 3), and a protein of unknown function (ORF 4). Phylogenetic analysis focusing on RdRp and SP24 amino acid sequences revealed a close relationship between WfNgV1 and Bemisia tabaci negevirus 1, a negevirus sequence recently discovered in whiteflies from Israel. Both viruses form a clade sharing a most recent common ancestor with the proposed nelorpivirus and centivirus taxa. The putative glycoprotein from ORF 2 and SP24 (ORF 3) of WfNgV1 exhibit the characteristic topologies previously reported for negevirus counterparts. This marks the first reported negevirus-like sequence from whiteflies in the Americas. [ABSTRACT FROM AUTHOR]

S-glutathionylation, the reversible formation of mixed disulfides between glutathione(GSH) and cysteine residues in proteins, is a specific form of post-translational modification that plays important roles in various biological processes, including signal transduction, redox homeostasis, and metabolism inside cells. Experimentally identifying S-glutathionylation sites is labor-intensive and time consuming, whereas bioinformatics methods provide an alternative way to this problem by predicting S-glutathionylation sites in silico. The bioinformatics approaches give not only candidate sites for further experimental verification but also bio-chemical insights into the mechanism of S-glutathionylation. In this paper, we firstly collect experimentally determined S-glutathionylated proteins and their corresponding modification sites from the literature, and then propose a new method for predicting S-glutathionylation sites by employing machine learning methods based on protein sequence data. Promising results are obtained by our method with an AUC (area under ROC curve) score of 0.879 in 5-fold cross-validation, which demonstrates the predictive power of our proposed method. The datasets used in this work are available at http://csb.shu.edu.cn/SGDB. [ABSTRACT FROM AUTHOR]

Background: Phylogenetic study of protein sequences provides unique and valuable insights into the molecular and genetic basis of important medical and epidemiological problems as well as insights about the origins and development of physiological features in present day organisms. Consensus phylogenies based on the bootstrap and other resampling methods play a crucial part in analyzing the robustness of the trees produced for these analyses. Methodology: Our focus was to increase the number of bootstrap replications that can be performed on large protein datasets using the maximum parsimony, distance matrix, and maximum likelihood methods. We have modified the PHYLIP package using MPI to enable large-scale phylogenetic study of protein sequences, using a statistically robust number of bootstrapped datasets, to be performed in a moderate amount of time. This paper discusses the methodology used to parallelize the PHYLIP programs and reports the performance of the parallel PHYLIP programs that are relevant to the study of protein evolution on several protein datasets. Conclusions: Calculations that currently take a few days on a state of the art desktop workstation are reduced to calculations that can be performed over lunchtime on a modern parallel computer. Of the three protein methods tested, the maximum likelihood method scales the best, followed by the distance method, and then the maximum parsimony method. However, the maximum likelihood method requires significant memory resources, which limits its application to more moderately sized protein datasets. [ABSTRACT FROM AUTHOR]

Since the importance of DNA-binding proteins in multiple biomolecular functions has been recognized, an increasing number of researchers are attempting to identify DNA-binding proteins. In recent years, the machine learning methods have become more and more compelling in the case of protein sequence data soaring, because of their favorable speed and accuracy. In this paper, we extract three features from the protein sequence, namely NMBAC (Normalized Moreau-Broto Autocorrelation), PSSM-DWT (Position-specific scoring matrix—Discrete Wavelet Transform), and PSSM-DCT (Position-specific scoring matrix—Discrete Cosine Transform). We also employ feature selection algorithm on these feature vectors. Then, these features are fed into the training SVM (support vector machine) model as classifier to predict DNA-binding proteins. Our method applys three datasets, namely PDB1075, PDB594 and PDB186, to evaluate the performance of our approach. The PDB1075 and PDB594 datasets are employed for Jackknife test and the PDB186 dataset is used for the independent test. Our method achieves the best accuracy in the Jacknife test, from 79.20% to 86.23% and 80.5% to 86.20% on PDB1075 and PDB594 datasets, respectively. In the independent test, the accuracy of our method comes to 76.3%. The performance of independent test also shows that our method has a certain ability to be effectively used for DNA-binding protein prediction. The data and source code are at . [ABSTRACT FROM AUTHOR]

Protein classification is one of the critical problems in bioinformatics. Early studies used geometric distances and polygenetic-tree to classify proteins. These methods use binary trees to present protein classification. In this paper, we propose a new protein classification method, whereby theories of information and networks are used to classify the multivariate relationships of proteins. In this study, protein universe is modeled as an undirected network, where proteins are classified according to their connections. Our method is unsupervised, multivariate, and alignment-free. It can be applied to the classification of both protein sequences and structures. Nine examples are used to demonstrate the efficiency of our new method. [ABSTRACT FROM AUTHOR]

In this paper, we have proposed a novel alignment-free method for comparing the similarity of protein sequences. We first encode a protein sequence into a 440 dimensional feature vector consisting of a 400 dimensional Pseudo-Markov transition probability vector among the 20 amino acids, a 20 dimensional content ratio vector, and a 20 dimensional position ratio vector of the amino acids in the sequence. By evaluating the Euclidean distances among the representing vectors, we compare the similarity of protein sequences. We then apply this method into the ND5 dataset consisting of the ND5 protein sequences of 9 species, and the F10 and G11 datasets representing two of the xylanases containing glycoside hydrolase families, i.e., families 10 and 11. As a result, our method achieves a correlation coefficient of 0.962 with the canonical protein sequence aligner ClustalW in the ND5 dataset, much higher than those of other 5 popular alignment-free methods. In addition, we successfully separate the xylanases sequences in the F10 family and the G11 family and illustrate that the F10 family is more heat stable than the G11 family, consistent with a few previous studies. Moreover, we prove mathematically an identity equation involving the Pseudo-Markov transition probability vector and the amino acids content ratio vector. [ABSTRACT FROM AUTHOR]

The red abalone (Haliotis rufescens) represents North America's most important aquaculture species. Its hepatopancreas is rich in cellulases and other polysaccharide-degrading enzymes, which provide it the remarkable ability to digest cellulose-rich macroalgae; nevertheless, its cellulolytic systems are poorly explored. This manuscript describes some functional and structural properties of an endogenous trimeric glycosylated endoglucanase from H. rufescens. The purified enzyme showed a molecular mass of 23.4 kDa determined by MALDI-TOF mass spectrometry, which behaved as a homotrimer in gel filtration chromatography and zymograms. According to the periodic acid-Schiff reagent staining, detecting sugar moieties in SDS-PAGE gel confirmed that abalone cellulase is a glycoprotein. Hydrolysis of cello-oligosaccharides and p-nitrophenyl-β-D-glucopyranosides confirmed its endo/exoactivity. A maximum enzyme activity toward 0.5% (w/v) carboxymethylcellulose of 53.9 ± 1.0 U/mg was achieved at 45°C and pH 6.0. We elucidated the abalone cellulase primary structure using proteases and mass spectrometry methods. Based on these results and using a bioinformatic approach, we identified the gene encoding this enzyme and deduced its full-length amino acid sequence; the mature protein comprised 177 residues with a calculated molecular mass of 19.1 kDa and, according to sequence similarity, it was classified into the glycosyl-hydrolase family 45 subfamily B. An AlphaFold theoretical model and docking simulations with cellopentaose confirmed that abalone cellulase is a β-sheet rich protein, as also observed by circular dichroism experiments, with conserved catalytic residues: Asp26, Asn109, and Asp134. Interestingly, the AlphaFold-Multimer analysis indicated a trimeric assembly for abalone cellulase, which supported our experimental findings. The discovery and characterization of these enzymes may contribute to developing efficient cellulose bioconversion processes for biofuels and sustainable bioproducts. [ABSTRACT FROM AUTHOR]

Monogenic diabetes is characterized as a group of diseases caused by rare variants in single genes. Like for other rare diseases, multiple genes have been linked to monogenic diabetes with different measures of pathogenicity, but the information on the genes and variants is not unified among different resources, making it challenging to process them informatically. We have developed an automated pipeline for collecting and harmonizing data on genetic variants linked to monogenic diabetes. Furthermore, we have translated variant genetic sequences into protein sequences accounting for all protein isoforms and their variants. This allows researchers to consolidate information on variant genes and proteins linked to monogenic diabetes and facilitates their study using proteomics or structural biology. Our open and flexible implementation using Jupyter notebooks enables tailoring and modifying the pipeline and its application to other rare diseases. [ABSTRACT FROM AUTHOR]

The black soldier fly (BSF), Hermetia illucens, has the potential to serve as a valuable resource for waste bioconversion due to the ability of the larvae to thrive in a microbial-rich environment. Being an ecological decomposer, the survival of BSF larvae (BSFL) relies on developing an efficient defense system. Cathepsin L (CTSL) is a cysteine protease that plays roles in physiological and pathological processes. In this study, the full-length of CTSL was obtained from BSF. The 1,020-bp open reading frame encoded a preprotein of 339 amino acids with a predicted molecular weight of 32 kDa. The pro-domain contained the conserved ERFNIN, GNYD, and GCNGG motifs, which are all characteristic of CTSL. Homology revealed that the deduced amino acid sequence of BSF CTSL shared 74.22–72.99% identity with Diptera flies. Immunohistochemical (IHC) analysis showed the CTSL was predominantly localized in the gut, especially in the midgut. The mRNA expression of CTSL in different larval stages was analyzed by quantitative real-time PCR (RT-qPCR), which revealed that CTSL was expressed in the second to sixth instar, with the highest expression in the fifth instar. Following an immune challenge in vivo using Escherichia coli (E. coli), CTSL mRNA was significantly up-regulated at 6 h post-stimulation. The Z-Phe-Arg-AMC was gradually cleaved by the BSFL extract after 3 h post-stimulation. These results shed light on the potential role of CTSL in the defense mechanism that helps BSFL to survive against pathogens in a microbial-rich environment. [ABSTRACT FROM AUTHOR]

Analbuminemia is characterized by the near absence of albumin in the plasma. Different methods are available for measuring albumin levels, but they do not necessarily agree with one another. It is a concern that analbuminemic samples could be falsely characterized due to the incorrect estimation of albumin. The objective of the work was to evaluate the performance of different assays in detecting analbuminemia. Albumin knockout (Alb-/-) mouse plasma was used to test the suitability of different albumin assays for their ability to properly characterize extreme albumin deficiency. Bromocresol green (BCG), bromocresol purple (BCP), enzyme-linked immunosorbent assay (ELISA), liquid chromatography-tandem mass spectrometry (LC-MS/MS), and gel electrophoresis were tested. The LC-MS/MS assay exhibited broad coverage of the amino acid sequence of albumin and indicated 8,400-fold lower (P<0.0001) albumin expression in Alb-/- than wildtype (WT), demonstrating its suitability for identifying extreme albumin deficiency. ELISA estimated albumin at 1.5±0.1 g/dL in WT and was below the detection limit in all Alb-/- samples. Gel electrophoresis yielded consistent results with LC-MS/MS and ELISA. The BCG assay overestimated albumin with apparently appreciable albumin concentrations in Alb-/- mice, yet the assay still indicated a significant difference between genotypes (Alb-/-, 1.2±0.05 g/dL, WT, 3.7±0.1 g/dL, P<0.0001). BCP drastically overestimated albumin and could not successfully identify the known analbuminemic phenotype of Alb-/- mice. By using Alb-/- plasma as a reference material and LC-MS/MS as a reference method, ELISA and gel electrophoresis appear appropriate for identifying analbuminemia, while BCG and BCP are not suitable. It is concluded that dye-binding assays should be avoided when extreme hypoalbuminemia or analbuminemia is suspected. [ABSTRACT FROM AUTHOR]

FASTA and FASTQ are basic and ubiquitous formats for storing nucleotide and protein sequences. Common manipulations of FASTA/Q file include converting, searching, filtering, deduplication, splitting, shuffling, and sampling. Existing tools only implement some of these manipulations, and not particularly efficiently, and some are only available for certain operating systems. Furthermore, the complicated installation process of required packages and running environments can render these programs less user friendly. This paper describes a cross-platform ultrafast comprehensive toolkit for FASTA/Q processing. SeqKit provides executable binary files for all major operating systems, including Windows, Linux, and Mac OSX, and can be directly used without any dependencies or pre-configurations. SeqKit demonstrates competitive performance in execution time and memory usage compared to similar tools. The efficiency and usability of SeqKit enable researchers to rapidly accomplish common FASTA/Q file manipulations. SeqKit is open source and available on Github at . [ABSTRACT FROM AUTHOR]

Pectate lyases are known to play a key role in pectin degradation by catalyzing the random cleavage of internal polymer linkages (endo-pectinases). In this paper, four novel cDNAs, designated Hg-pel-3, Hg-pel-4, Hg-pel-6 and Hg-pel-7, that encode pectate lyases were cloned and characterized from the soybean cyst nematode, Heterodera glycines. The predicted protein sequences of HG-PEL-3, HG-PEL-4 and HG-PEL-6 differed significantly in both their amino acid sequences and their genomic structures from other pectate lyases of H. glycines (HG-PEL-1, HG-PEL-2 and HG-PEL-7). A phylogenetic study revealed that the pectate lyase proteins of H. glycines are clustered into distinct clades and have distinct numbers and positioning of introns, which suggests that the pectate lyase genes of H. glycines may have evolved from at least two ancestral genes. A Southern blot analysis revealed that multiple Hg-pel-6-like genes were present in the H. glycines genome. In situ hybridization showed that four novel pectate lyases (Hg-pel-3, Hg-pel-4, Hg-pel-6 and Hg-pel-7) were actively transcribed in the subventral esophageal gland cells. A semi-quantitative RT-PCR assay supported the finding that the expression of these genes was strong in the egg, pre-parasitic second-stage juvenile (J2) and early parasitic J2 stages and that it declined in further developmental stages of the nematode. This expression pattern suggests that these proteins play a role in the migratory phase of the nematode life cycle. Knocking down Hg-pel-6 using in vitro RNA interference resulted in a 46.9% reduction of the number of nematodes that invaded the plants and a 61.5% suppression of the development of H. glycines females within roots compared to the GFP-dsRNA control. Plant host-derived RNAi induced the silencing of the Hg-pel-6gene, which significantly reduced the nematode infection levels at 7 Days post inoculation (dpi). Similarly, this procedure reduced the number of female adults at 40 dpi, which suggests the important roles of this gene in the early stages of parasitism. Our combined data suggest that two types of pectate lyases are present in the H. glycines genome and may have different roles during infection. [ABSTRACT FROM AUTHOR]

We introduce a new representation and feature extraction method for biological sequences. Named bio-vectors (BioVec) to refer to biological sequences in general with protein-vectors (ProtVec) for proteins (amino-acid sequences) and gene-vectors (GeneVec) for gene sequences, this representation can be widely used in applications of deep learning in proteomics and genomics. In the present paper, we focus on protein-vectors that can be utilized in a wide array of bioinformatics investigations such as family classification, protein visualization, structure prediction, disordered protein identification, and protein-protein interaction prediction. In this method, we adopt artificial neural network approaches and represent a protein sequence with a single dense n-dimensional vector. To evaluate this method, we apply it in classification of 324,018 protein sequences obtained from Swiss-Prot belonging to 7,027 protein families, where an average family classification accuracy of 93%±0.06% is obtained, outperforming existing family classification methods. In addition, we use ProtVec representation to predict disordered proteins from structured proteins. Two databases of disordered sequences are used: the DisProt database as well as a database featuring the disordered regions of nucleoporins rich with phenylalanine-glycine repeats (FG-Nups). Using support vector machine classifiers, FG-Nup sequences are distinguished from structured protein sequences found in Protein Data Bank (PDB) with a 99.8% accuracy, and unstructured DisProt sequences are differentiated from structured DisProt sequences with 100.0% accuracy. These results indicate that by only providing sequence data for various proteins into this model, accurate information about protein structure can be determined. Importantly, this model needs to be trained only once and can then be applied to extract a comprehensive set of information regarding proteins of interest. Moreover, this representation can be considered as pre-training for various applications of deep learning in bioinformatics. The related data is available at Life Language Processing Website: and Harvard Dataverse: . [ABSTRACT FROM AUTHOR]

Playing crucial roles in various cellular processes, such as recognition of specific nucleotide sequences, regulation of transcription, and regulation of gene expression, DNA-binding proteins are essential ingredients for both eukaryotic and prokaryotic proteomes. With the avalanche of protein sequences generated in the postgenomic age, it is a critical challenge to develop automated methods for accurate and rapidly identifying DNA-binding proteins based on their sequence information alone. Here, a novel predictor, called "iDNA-Prot|dis", was established by incorporating the amino acid distance-pair coupling information and the amino acid reduced alphabet profile into the general pseudo amino acid composition (PseAAC) vector. The former can capture the characteristics of DNA-binding proteins so as to enhance its prediction quality, while the latter can reduce the dimension of PseAAC vector so as to speed up its prediction process. It was observed by the rigorous jackknife and independent dataset tests that the new predictor outperformed the existing predictors for the same purpose. As a user-friendly web-server, iDNA-Prot|dis is accessible to the public at http://bioinformatics.hitsz.edu.cn/iDNA-Prot&#95;dis/. Moreover, for the convenience of the vast majority of experimental scientists, a step-by-step protocol guide is provided on how to use the web-server to get their desired results without the need to follow the complicated mathematic equations that are presented in this paper just for the integrity of its developing process. It is anticipated that the iDNA-Prot|dis predictor may become a useful high throughput tool for large-scale analysis of DNA-binding proteins, or at the very least, play a complementary role to the existing predictors in this regard.

The most informative probability distribution functions (PDFs) describing the Ramachandran phi-psi dihedral angle pair, a fundamental descriptor of backbone conformation of protein molecules, are derived from high-resolution X-ray crystal structures using an information-theoretic approach. The Information Maximization Device (IMD) is established, based on fundamental information-theoretic concepts, and then applied specifically to derive highly resolved phi-psi maps for all 20 single amino acid and all 8000 triplet sequences at an optimal resolution determined by the volume of current data. The paper shows that utilizing the latent information contained in all viable high-resolution crystal structures found in the Protein Data Bank (PDB), totaling more than 77,000 chains, permits the derivation of a large number of optimized sequence-dependent PDFs. This work demonstrates the effectiveness of the IMD and the superiority of the resulting PDFs by extensive fold recognition experiments and rigorous comparisons with previously published triplet PDFs. Because it automatically optimizes PDFs, IMD results in improved performance of knowledge-based potentials, which rely on such PDFs. Furthermore, it provides an easy computational recipe for empirically deriving other kinds of sequence-dependent structural PDFs with greater detail and precision. The high-resolution phi-psi maps derived in this work are available for download.

Knowledge-based potentials are energy functions derived from the analysis of databases of protein structures and sequences. They can be divided into two classes. Potentials from the first class are based on a direct conversion of the distributions of some geometric properties observed in native protein structures into energy values, while potentials from the second class are trained to mimic quantitatively the geometric differences between incorrectly folded models and native structures. In this paper, we focus on the relationship between energy and geometry when training the second class of knowledge-based potentials. We assume that the difference in energy between a decoy structure and the corresponding native structure is linearly related to the distance between the two structures. We trained two distance-based knowledge-based potentials accordingly, one based on all inter-residue distances (PPD), while the other had the set of all distances filtered to reflect consistency in an ensemble of decoys (PPE). We tested four types of metric to characterize the distance between the decoy and the native structure, two based on extrinsic geometry (RMSD and GTD-TS*), and two based on intrinsic geometry (Q* and MT). The corresponding eight potentials were tested on a large collection of decoy sets. We found that it is usually better to train a potential using an intrinsic distance measure. We also found that PPE outperforms PPD, emphasizing the benefits of capturing consistent information in an ensemble. The relevance of these results for the design of knowledge-based potentials is discussed. [ABSTRACT FROM AUTHOR]

Phylogenetic analysis of protein sequences provides a powerful means of identifying novel protein functions and subfamilies, and for identifying and resolving annotation errors. However, automation of functional clustering based on phylogenetic trees has been challenging and most of it is done manually. Clustering phylogenetic trees usually requires the delineation of tree-based thresholds (e.g., distances), leading to an ad hoc problem. We propose a new phylogenetic clustering approach that identifies clusters without using ad hoc distances or other pre-defined values. Our workflow combines uniform manifold approximation and projection (UMAP) with Gaussian mixture models as a k-means like procedure to automatically group sequences into clusters. We then apply a "second pass" clade identification algorithm to resolve non-monophyletic groups. We tested our approach with several well-curated protein families (outer membrane porins, acyltransferase, and nuclear receptors) and showed our automated methods recapitulated known subfamilies. We also applied our methods to a broad range of different protein families from multiple databases, including Pfam, PANTHER, and UniProt, and to alignments of RNA viral genomes. Our results showed that AutoPhy rapidly generated monophyletic clusters (subfamilies) within phylogenetic trees evolving at very different rates both within and among phylogenies. The phylogenetic clusters generated by AutoPhy resolved misannotations and identified new protein functional groups and novel viral strains. [ABSTRACT FROM AUTHOR]

Spinocerebellar ataxia 2 (SCA2) is a neurodegenerative disorder caused by the expansion of the poly-glutamine (polyQ) tract of Ataxin-2 (ATXN2). Other polyQ-containing proteins such as ATXN7 and huntingtin are associated with the development of neurodegenerative diseases when their N-terminal polyQ domains are expanded. Furthermore, they undergo proteolytic processing events that produce N-terminal fragments that include the polyQ stretch, which are implicated in pathogenesis. Interestingly, N-terminal ATXN2 fragments were reported in a brain extract from a SCA2 patient, but it is currently unknown whether an expanded polyQ domain contributes to ATXN2 proteolytic susceptibility. Here, we used transient expression in HEK293 cells to determine whether ATXN2 is a target for specific N-terminal proteolysis. We found that ATXN2 proteins with either normal or expanded polyQ stretches undergo proteolytic cleavage releasing an N-terminal polyQ-containing fragment. We identified a short amino acid sequence downstream of the polyQ domain that is necessary for N-terminal cleavage of full-length ATXN2 and sufficient to induce proteolysis of a heterologous protein. However, this sequence is not required for cleavage of a short ATXN2 isoform produced from an alternative start codon located just upstream of the CAG repeats encoding the polyQ domain. Our study extends our understanding of ATXN2 posttranslational regulation by revealing that this protein can be the target of specific proteolytic cleavage events releasing polyQ-containing products that are modulated by the N-terminal domain of ATXN2. N-terminal ATXN2 proteolysis of expanded polyQ domains might contribute to SCA2 pathology, as observed in other neurodegenerative disorders caused by polyQ domain expansion. [ABSTRACT FROM AUTHOR]

Proteins are fundamental components of diverse cellular systems and play crucial roles in a variety of disease processes. Consequently, it is crucial to comprehend their structure, function, and intricate interconnections. Classifying proteins into families or groups with comparable structural and functional characteristics is a crucial aspect of this comprehension. This classification is crucial for evolutionary research, predicting protein function, and identifying potential therapeutic targets. Sequence alignment and structure-based alignment are frequently ineffective techniques for identifying protein families.This study addresses the need for a more efficient and accurate technique for feature extraction and protein classification. The research proposes a novel method that integrates bispectrum characteristics, deep learning techniques, and machine learning algorithms to overcome the limitations of conventional methods. The proposed method uses numbers to represent protein sequences, utilizes bispectrum analysis, uses different topologies for convolutional neural networks to pull out features, and chooses robust features to classify protein families. The goal is to outperform existing methods for identifying protein families, thereby enhancing classification metrics. The materials consist of numerous protein datasets, whereas the methods incorporate bispectrum characteristics and deep learning strategies. The results of this study demonstrate that the proposed method for identifying protein families is superior to conventional approaches. Significantly enhanced quality metrics demonstrated the efficacy of the combined bispectrum and deep learning approaches. These findings have the potential to advance the field of protein biology and facilitate pharmaceutical innovation. In conclusion, this study presents a novel method that employs bispectrum characteristics and deep learning techniques to improve the precision and efficiency of protein family identification. The demonstrated advancements in classification metrics demonstrate this method's applicability to numerous scientific disciplines. This furthers our understanding of protein function and its implications for disease and treatment. [ABSTRACT FROM AUTHOR]

Microtubules are polymeric filaments, constructed of α-β tubulin heterodimers that underlie critical subcellular structures in eukaryotic organisms. Four homologous proteins (γ-, δ-, ε- and ζ-tubulin) additionally contribute to specialized microtubule functions. Although there is an immense volume of publicly available data pertaining to tubulins, it is difficult to assimilate all potentially relevant information across diverse organisms, isotypes, and categories of data. We previously assembled an extensive web-based catalogue of published missense mutations to tubulins with >1,500 entries that each document a specific substitution to a discrete tubulin, the species where the mutation was described and the associated phenotype with hyperlinks to the amino acid sequence and citation(s) for research. This report describes a significant update and expansion of our online resource (TubulinDB.bio.uci.edu) to nearly 18,000 entries. It now encompasses a cross-referenced catalog of post-translational modifications (PTMs) to tubulin drawn from public datasets, primary literature, and predictive algorithms. In addition, tubulin protein structures were used to define local interactions with bound ligands (GTP, GDP and diverse microtubule-targeting agents) and amino acids at the intradimer interface, within the microtubule lattice and with associated proteins. To effectively cross-reference these datasets, we established a universal tubulin numbering system to map entries into a common framework that accommodates specific insertions and deletions to tubulins. Indexing and cross-referencing permitted us to discern previously unappreciated patterns. We describe previously unlinked observations of loss of PTM sites in the context of cancer cells and tubulinopathies. Similarly, we expanded the set of clinical substitutions that may compromise MAP or microtubule-motor interactions by collecting tubulin missense mutations that alter amino acids at the interface with dynein and doublecortin. By expanding the database as a curated resource, we hope to relate model organism data to clinical findings of pathogenic tubulin variants. Ultimately, we aim to aid researchers in hypothesis generation and design of studies to dissect tubulin function. [ABSTRACT FROM AUTHOR]

Nitrotyrosine is one of the post-translational modifications (PTMs) in proteins that occurs when their tyrosine residue is nitrated. Compared with healthy people, a remarkably increased level of nitrotyrosine is detected in those suffering from rheumatoid arthritis, septic shock, and coeliac disease. Given an uncharacterized protein sequence that contains many tyrosine residues, which one of them can be nitrated and which one cannot? This is a challenging problem, not only directly related to in-depth understanding the PTM’s mechanism but also to the nitrotyrosine-based drug development. Particularly, with the avalanche of protein sequences generated in the postgenomic age, it is highly desired to develop a high throughput tool in this regard. Here, a new predictor called “iNitro-Tyr” was developed by incorporating the position-specific dipeptide propensity into the general pseudo amino acid composition for discriminating the nitrotyrosine sites from non-nitrotyrosine sites in proteins. It was demonstrated via the rigorous jackknife tests that the new predictor not only can yield higher success rate but also is much more stable and less noisy. A web-server for iNitro-Tyr is accessible to the public at http://app.aporc.org/iNitro-Tyr/. For the convenience of most experimental scientists, we have further provided a protocol of step-by-step guide, by which users can easily get their desired results without the need to follow the complicated mathematics that were presented in this paper just for the integrity of its development process. It has not escaped our notice that the approach presented here can be also used to deal with the other PTM sites in proteins. [ABSTRACT FROM AUTHOR]

One of the fundamental tasks in biology is to identify the functions of all proteins to reveal the primary machinery of a cell. Knowledge of the subcellular locations of proteins will provide key hints to reveal their functions and to understand the intricate pathways that regulate biological processes at the cellular level. Protein subcellular location prediction has been extensively studied in the past two decades. A lot of methods have been developed based on protein primary sequences as well as protein-protein interaction network. In this paper, we propose to use the protein-protein interaction network as an infrastructure to integrate existing sequence based predictors. When predicting the subcellular locations of a given protein, not only the protein itself, but also all its interacting partners were considered. Unlike existing methods, our method requires neither the comprehensive knowledge of the protein-protein interaction network nor the experimentally annotated subcellular locations of most proteins in the protein-protein interaction network. Besides, our method can be used as a framework to integrate multiple predictors. Our method achieved 56% on human proteome in absolute-true rate, which is higher than the state-of-the-art methods. [ABSTRACT FROM AUTHOR]

The aggregation of proteins or peptides in amyloid fibrils is associated with a number of clinical disorders, including Alzheimer's, Huntington's and prion diseases, medullary thyroid cancer, renal and cardiac amyloidosis. Despite extensive studies, the molecular mechanisms underlying the initiation of fibril formation remain largely unknown. Several lines of evidence revealed that short amino-acid segments (hot spots), located in amyloid precursor proteins act as seeds for fibril elongation. Therefore, hot spots are potential targets for diagnostic/therapeutic applications, and a current challenge in bioinformatics is the development of methods to accurately predict hot spots from protein sequences. In this paper, we combined existing methods into a meta-predictor for hot spots prediction, called MetAmyl for METapredictor for AMYLoid proteins. MetAmyl is based on a logistic regression model that aims at weighting predictions from a set of popular algorithms, statistically selected as being the most informative and complementary predictors. We evaluated the performances of MetAmyl through a large scale comparative study based on three independent datasets and thus demonstrated its ability to differentiate between amyloidogenic and non-amyloidogenic polypeptides. Compared to 9 other methods, MetAmyl provides significant improvement in prediction on studied datasets. We further show that MetAmyl is efficient to highlight the effect of point mutations involved in human amyloidosis, so we suggest this program should be a useful complementary tool for the diagnosis of these diseases. [ABSTRACT FROM AUTHOR]

Identifying shared sequence segments along amino acid sequences generally requires a collection of closely related proteins, most often curated manually from the sequence datasets to suit the purpose at hand. Currently developed statistical methods are strained, however, when the collection contains remote sequences with poor alignment to the rest, or sequences containing multiple domains. In this paper, we propose a completely unsupervised and automated method to identify the shared sequence segments observed in a diverse collection of protein sequences including those present in a smaller fraction of the sequences in the collection, using a combination of sequence alignment, residue conservation scoring and graph-theoretical approaches. Since shared sequence fragments often imply conserved functional or structural attributes, the method produces a table of associations between the sequences and the identified conserved regions that can reveal previously unknown protein families as well as new members to existing ones. We evaluated the biological relevance of the method by clustering the proteins in gold standard datasets and assessing the clustering performance in comparison with previous methods from the literature. We have then applied the proposed method to a genome wide dataset of 17793 human proteins and generated a global association map to each of the 4753 identified conserved regions. Investigations on the major conserved regions revealed that they corresponded strongly to annotated structural domains. This suggests that the method can be useful in predicting novel domains on protein sequences. [ABSTRACT FROM AUTHOR]

Cell survival in complex, vascularized tissues, has been implicated as a major bottleneck in advancement of therapies based on cardiac tissue engineering. This limitation motivates the search for small, inexpensive molecules that would simultaneously be cardio-protective and vasculogenic. Here, we present peptide sequence QHREDGS, based upon the fibrinogen-like domain of angiopoietin-1, as a prime candidate molecule. We demonstrated previously that QHREDGS improved cardiomyocyte metabolism and mitigated serum starved apoptosis. In this paper we further demonstrate the potency of QHREDGS in its ability to enhance endothelial cell survival, metabolism and tube formation. When endothelial cells were exposed to the soluble form of QHREDGS, improvements in endothelial cell barrier functionality, nitric oxide production and cell metabolism (ATP levels) in serum starved conditions were found. The functionality of the peptide was then examined when conjugated to collagen-chitosan hydrogel, a potential carrier for in vivo application. The presence of the peptide in the hydrogel mitigated paclitaxel induced apoptosis of endothelial cells in a dose dependent manner. Furthermore, the peptide modified hydrogels stimulated tube-like structure formation of encapsulated endothelial cells. When integrin αvβ3 or α5β1were antibody blocked during cell encapsulation in peptide modified hydrogels, tube formation was abolished. Therefore, the dual protective nature of the novel peptide QHREDGS may position this peptide as an appealing augmentation for collagen-chitosan hydrogels that could be used for biomaterial delivered cell therapies in the settings of myocardial infarction. [ABSTRACT FROM AUTHOR]

S100A4 and S100A9 proteins have been described as playing roles in the control of tumor growth and metastasis. We show here that a chemical probe, oxyclozanide (OX), selected for inhibiting the interaction between S100A9 and the receptor for advanced glycation end-products (RAGE) interacts with both S100A9 and S100A4. Furthermore, we show that S100A9 and S100A4 interact with RAGE and TLR4; interactions that can be inhibited by OX. Hence, S100A4 and S100A9 display similar functional elements despite their primary sequence diversity. This was further confirmed by showing that S100A4 and S100A9 dimerize both in vitro and in vivo. All of these interactions required levels of Zn++ that are found in the extracellular space but not intracellularly. Interestingly, S100A4 and S100A9 are expressed by distinct CD11b+ subpopulations both in healthy animals and in animals with either inflammatory disease or tumor burden. The functions of S100A9 and S100A4 described in this paper, including heterodimerization, may therefore reflect S100A9 and S100A4 that are released into the extra-cellular milieu. [ABSTRACT FROM AUTHOR]

Protein sequence databases are the pillar upon which modern proteomics is supported, representing a stable reference space of predicted and validated proteins. One example of such resources is UniProt, enriched with both expertly curated and automatic annotations. Taken largely for granted, similar mature resources such as UniProt are not available yet in some other “omics” fields, lipidomics being one of them. While having a seasoned community of wet lab scientists, lipidomics lies significantly behind proteomics in the adoption of data standards and other core bioinformatics concepts. This work aims to reduce the gap by developing an equivalent resource to UniProt called ‘LipidHome’, providing theoretically generated lipid molecules and useful metadata. Using the ‘FASTLipid’ Java library, a database was populated with theoretical lipids, generated from a set of community agreed upon chemical bounds. In parallel, a web application was developed to present the information and provide computational access via a web service. Designed specifically to accommodate high throughput mass spectrometry based approaches, lipids are organised into a hierarchy that reflects the variety in the structural resolution of lipid identifications. Additionally, cross-references to other lipid related resources and papers that cite specific lipids were used to annotate lipid records. The web application encompasses a browser for viewing lipid records and a ‘tools’ section where an MS1 search engine is currently implemented. LipidHome can be accessed at http://www.ebi.ac.uk/apweiler-srv/lipidhome. [ABSTRACT FROM AUTHOR]

In functional metagenomics, BLAST homology search is a common method to classify metagenomic reads into protein/domain sequence families such as Clusters of Orthologous Groups of proteins (COGs) in order to quantify the abundance of each COG in the community. The resulting functional profile of the community is then used in downstream analysis to correlate the change in abundance to environmental perturbation, clinical variation, and so on. However, the short read length coupled with next-generation sequencing technologies poses a barrier in this approach, essentially because similarity significance cannot be discerned by searching with short reads. Consequently, artificial functional families are produced, in which those with a large number of reads assigned decreases the accuracy of functional profile dramatically. There is no method available to address this problem. We intended to fill this gap in this paper. We revealed that BLAST similarity scores of homologues for short reads from COG protein members coding sequences are distributed differently from the scores of those derived elsewhere. We showed that, by choosing an appropriate score cut-off, we are able to filter out most artificial families and simultaneously to preserve sufficient information in order to build the functional profile. We also showed that, by incorporated application of BLAST and RPS-BLAST, some artificial families with large read counts can be further identified after the score cutoff filtration. Evaluated on three experimental metagenomic datasets with different coverages, we found that the proposed method is robust against read coverage and consistently outperforms the other E-value cutoff methods currently used in literatures. [ABSTRACT FROM AUTHOR]

The analysis of biological information from protein sequences is important for the study of cellular functions and interactions, and protein fold recognition plays a key role in the prediction of protein structures. Unfortunately, the prediction of protein fold patterns is challenging due to the existence of compound protein structures. Here, we processed the latest release of the Structural Classification of Proteins (SCOP, version 1.75) database and exploited novel techniques to impressively increase the accuracy of protein fold classification. The techniques proposed in this paper include ensemble classifying and a hierarchical framework, in the first layer of which similar or redundant sequences were deleted in two manners; a set of base classifiers, fused by various selection strategies, divides the input into seven classes; in the second layer of which, an analogous ensemble method is adopted to predict all protein folds. To our knowledge, it is the first time all protein folds can be intelligently detected hierarchically. Compared with prior studies, our experimental results demonstrated the efficiency and effectiveness of our proposed method, which achieved a success rate of 74.21%, which is much higher than results obtained with previous methods (ranging from 45.6% to 70.5%). When applied to the second layer of classification, the prediction accuracy was in the range between 23.13% and 46.05%. This value, which may not be remarkably high, is scientifically admirable and encouraging as compared to the relatively low counts of proteins from most fold recognition programs. The web server Hierarchical Protein Fold Prediction (HPFP) is available at http://datamining.xmu.edu.cn/software/hpfp. [ABSTRACT FROM AUTHOR]

Background: Shotgun proteomics data analysis usually relies on database search. Because commonly employed protein sequence databases of most species do not contain sufficient protein information, the application of shotgun proteomics to the research of protein sequence profile remains a big challenge, especially to the species whose genome has not been sequenced yet. Methodology/Principal Findings: In this paper, we present a workflow with integrated database to partly address this problem. First, we downloaded the homologous species database. Next, we identified the transcriptome of the sample, created a protein sequence database based on the transcriptome data, and integtrated it with homologous species database. Lastly, we developed a workflow for identifying peptides simultaneously from shotgun proteomics data. Conclusions/Significance: We used datasets from orange leaves samples to demonstrate our workflow. The results showed that the integrated database had great advantage on orange shotgun proteomics data analysis compared to the homologous species database, an 18.5% increase in number of proteins identification. [ABSTRACT FROM AUTHOR]