Search

Lysine crotonylation is a newly discovered and reversible posttranslational modification involved in various biological processes, especially metabolism regulation. A total of 5159 lysine crotonylation sites in 2272 protein groups were identified. Twenty-seven motifs were found to be the preferred amino acid sequences for crotonylation sites. Functional annotation analyses revealed that most crotonylated proteins play important roles in metabolic processes and photosynthesis. Bioinformatics analysis suggested that lysine crotonylation preferentially targets a variety of important biological processes, including ribosome, glyoxylate and dicarboxylate metabolism, carbon fixation in photosynthetic organisms, proteasome and the TCA cycle, indicating lysine crotonylation is involved in the common mechanism of metabolic regulation. A protein interaction network analysis revealed that diverse interactions are modulated by protein crotonylation. These results suggest that lysine crotonylation is involved in a variety of biological processes. HSP70 is a crucial protein involved in protecting plant cells and tissues from thermal or abiotic stress responses, and HSP70 protein was found to be crotonylated in paper mulberry. This systematic analysis provides the first comprehensive analysis of lysine crotonylation in paper mulberry and provides important resources for further study on the regulatory mechanism and function of the lysine crotonylated proteome. [ABSTRACT FROM AUTHOR]

Background: The initial version of SEDA assists life science researchers without programming skills with the preparation of DNA and protein sequence FASTA files for multiple bioinformatics applications. However, the initial version of SEDA lacks a command-line interface for more advanced users and does not allow the creation of automated analysis pipelines. Results: The present paper discusses the updates of the new SEDA release, including the addition of a complete command-line interface, new functionalities like gene annotation, a framework for automated pipelines, and improved integration in Linux environments. Conclusion: SEDA is an open-source Java application and can be installed using the different distributions available (https://www.sing-group.org/seda/download.html) as well as through a Docker image (https://hub.docker.com/r/pegi3s/seda). It is released under a GPL-3.0 license, and its source code is publicly accessible on GitHub (https://github.com/sing-group/seda). The software version at the time of submission is archived at Zenodo (version v1.6.0, http://doi.org/10.5281/zenodo.10201605). [ABSTRACT FROM AUTHOR]

Bicyclic peptides have attracted the interest of pharmaceutical companies because of their remarkable properties, putting them on a new path in medicine. Their conformational rigidity improves proteolytic stability and leads to rapid penetration into tissues via any possible route of administration. Moreover, elimination of renal metabolism is of great importance, for example, for people with a history of liver diseases. In addition, each ring can function independently, making bicyclic peptides extremely versatile molecules for further optimization. In this paper, we compared the potentiometric and spectroscopic properties studied by UV–vis, MCD, and EPR of four synthetic analogues of the bi-cyclic peptide c(PKKHP-c(CFWKTC)-PKKH) (BCL). In particular, we correlated the structural and spectral properties of complexes with coordinating abilities toward Cu(II) ions of MCL1 (Ac-PKKHPc(CFWKTC)PKKH-NH2) that contains the unbinding cycle and N- and C-terminal linear parts with two histidine residues, one per part; two monocyclic ligands containing one histidine residue, both in the N-terminal position, i.e., MCL2 (Ac-PKKHPc(CFWKTC)PKKS-NH2) and in the C-terminal position, i.e., MCL3 (Ac-PKKSPc(CFWKTC)PKKH-NH2), respectively; and the linear structure LNL (Ac-PKKHPSFWKTSPKKH-NH2). Potentiometric results have shown that the bicyclic structure promotes the involvement of the side chain imidazole donors in Cu(II) binding. On the other hand, the results obtained for the mono-cyclic analogues lead to the conclusion that the coordination of the histidine moiety as an anchoring group is promoted by its location in the peptide sequence further from the nonbinding cycle, strongly influencing the involvement of the amide donors in Cu(II) coordination. [ABSTRACT FROM AUTHOR]

Protein-ligand binding affinity plays a pivotal role in drug development, particularly in identifying potential ligands for target disease-related proteins. Accurate affinity predictions can significantly reduce both the time and cost involved in drug development. However, highly precise affinity prediction remains a research challenge. A key to improve affinity prediction is to capture interactions between proteins and ligands effectively. Existing deep-learning-based computational approaches use 3D grids, 4D tensors, molecular graphs, or proximity-based adjacency matrices, which are either resource-intensive or do not directly represent potential interactions. In this paper, we propose atomic-level distance features and attention mechanisms to capture better specific protein-ligand interactions based on donor-acceptor relations, hydrophobicity, and π -stacking atoms. We argue that distances encompass both short-range direct and long-range indirect interaction effects while attention mechanisms capture levels of interaction effects. On the very well-known CASF-2016 dataset, our proposed method, named Distance plus Attention for Affinity Prediction (DAAP), significantly outperforms existing methods by achieving Correlation Coefficient (R) 0.909, Root Mean Squared Error (RMSE) 0.987, Mean Absolute Error (MAE) 0.745, Standard Deviation (SD) 0.988, and Concordance Index (CI) 0.876. The proposed method also shows substantial improvement, around 2% to 37%, on five other benchmark datasets. The program and data are publicly available on the website https://gitlab.com/mahnewton/daap. Scientific Contribution Statement This study innovatively introduces distance-based features to predict protein-ligand binding affinity, capitalizing on unique molecular interactions. Furthermore, the incorporation of protein sequence features of specific residues enhances the model's proficiency in capturing intricate binding patterns. The predictive capabilities are further strengthened through the use of a deep learning architecture with attention mechanisms, and an ensemble approach, averaging the outputs of five models, is implemented to ensure robust and reliable predictions. [ABSTRACT FROM AUTHOR]

RAMOSA1 (RA1) is a Cys2-His2-type (C2H2) zinc finger transcription factor that controls plant meristem fate and identity and has played an important role in maize domestication. Despite its importance, the origin of RA1 is unknown, and the evolution in plants is only partially understood. In this paper, we present a well-resolved phylogeny based on 73 amino acid sequences from 48 embryophyte species. The recovered tree topology indicates that, during grass evolution, RA1 arose from two consecutive SUPERMAN duplications, resulting in three distinct grass sequence lineages: RA1-like A, RA1-like B, and RA1; however, most of these copies have unknown functions. Our findings indicate that RA1 and RA1-like play roles in the nucleus despite lacking a traditional nuclear localization signal. Here, we report that copies diversified their coding region and, with it, their protein structure, suggesting different patterns of DNA binding and protein–protein interaction. In addition, each of the retained copies diversified regulatory elements along their promoter regions, indicating differences in their upstream regulation. Taken together, the evidence indicates that the RA1 and RA1-like gene families in grasses underwent subfunctionalization and neofunctionalization enabled by gene duplication. [ABSTRACT FROM AUTHOR]

Purpose: Sequenced Protein–Protein Interaction (PPI) prediction represents a pivotal area of study in biology, playing a crucial role in elucidating the mechanistic underpinnings of diseases and facilitating the design of novel therapeutic interventions. Conventional methods for extracting features through experimental processes have proven to be both costly and exceedingly complex. In light of these challenges, the scientific community has turned to computational approaches, particularly those grounded in deep learning methodologies. Despite the progress achieved by current deep learning technologies, their effectiveness diminishes when applied to larger, unfamiliar datasets. Results: In this study, the paper introduces a novel deep learning framework, termed DL-PPI, for predicting PPIs based on sequence data. The proposed framework comprises two key components aimed at improving the accuracy of feature extraction from individual protein sequences and capturing relationships between proteins in unfamiliar datasets. 1. Protein Node Feature Extraction Module: To enhance the accuracy of feature extraction from individual protein sequences and facilitate the understanding of relationships between proteins in unknown datasets, the paper devised a novel protein node feature extraction module utilizing the Inception method. This module efficiently captures relevant patterns and representations within protein sequences, enabling more informative feature extraction. 2. Feature-Relational Reasoning Network (FRN): In the Global Feature Extraction module of our model, the paper developed a novel FRN that leveraged Graph Neural Networks to determine interactions between pairs of input proteins. The FRN effectively captures the underlying relational information between proteins, contributing to improved PPI predictions. DL-PPI framework demonstrates state-of-the-art performance in the realm of sequence-based PPI prediction. [ABSTRACT FROM AUTHOR]

Next-generation single-molecule protein sequencing technologies have the potential to significantly accelerate biomedical research. These technologies offer sensitivity and scalability for proteomic analysis. One auspicious method is fluorosequencing, which involves: cutting naturalized proteins into peptides, attaching fluorophores to specific amino acids, and observing variations in light intensity as one amino acid is removed at a time. The original peptide is classified from the sequence of light-intensity reads, and proteins can subsequently be recognized with this information. The amino acid step removal is achieved by attaching the peptides to a wall on the C-terminal and using a process called Edman Degradation to remove an amino acid from the N-Terminal. Even though a framework (Whatprot) has been proposed for the peptide classification task, processing times remain restrictive due to the massively parallel data acquisicion system. In this paper, we propose a new beam search decoder with a novel state formulation that obtains considerably lower processing times at the expense of only a slight accuracy drop compared to Whatprot. Furthermore, we explore how our novel state formulation may lead to even faster decoders in the future. Author summary: Proteomic analyses frequently rely on mass spectrometry, a method characterized by its limited dynamic range, potentially overlooking low-abundant proteins. To address this limitation, single-molecule protein sequencing methods offer a solution. Fluorosequencing is a cutting-edge single-molecule protein sequencing method, which can distinguish peptides or protein molecules massively parallelly. This method has attracted interest from investors, as evidenced by the recent funding of Erisyon, a company developing this technology. This technique contains a challenging classification task: determining the original peptide sequence from light-intensity observations obtained after several Edman cycles. A classifier based on a combination of k Nearest Neighbors (kNN) with Hidden Markov Models (HMM) had been shown to have close-to-optimal accuracy with tractable complexity. We propose in this paper a new algorithm that reduces computation time significantly at the expense of a slight reduction in accuracy compared to state-of-the-art method. [ABSTRACT FROM AUTHOR]

CRNDE is considered an oncogene expressed as long non-coding RNA. Our previous paper is the only one reporting CRNDE as a micropeptide-coding gene. The amino acid sequence of this micropeptide (CRNDEP) has recently been confirmed by other researchers. This study aimed at providing a mass spectrometry (MS)-based validation of the CRNDEP sequence and an investigation of how the differential expression of CRNDE(P) influences the metabolism and chemoresistance of ovarian cancer (OvCa) cells. We also assessed cellular localization changes of CRNDEP, looked for its protein partners, and bioinformatically evaluated its RNA-binding capacities. Herein, we detected most of the CRNDEP sequence by MS. Moreover, our results corroborated the oncogenic role of CRNDE, portraying it as the gene impacting carcinogenesis at the stages of DNA transcription and replication, affecting the RNA metabolism, and stimulating the cell cycle progression and proliferation, with CRNDEP being detected in the centrosomes of dividing cells. We also showed that CRNDEP is located in nucleoli and revealed interactions of this micropeptide with p54, an RNA helicase. Additionally, we proved that high CRNDE(P) expression increases the resistance of OvCa cells to treatment with microtubule-targeted cytostatics. Furthermore, altered CRNDE(P) expression affected the activity of the microtubular cytoskeleton and the formation of focal adhesion plaques. Finally, according to our in silico analyses, CRNDEP is likely capable of RNA binding. All these results contribute to a better understanding of the CRNDE(P) role in OvCa biology, which may potentially improve the screening, diagnosis, and treatment of this disease. [ABSTRACT FROM AUTHOR]

Motivation Eukaryotic linear motifs (ELMs), or Short Linear Motifs, are protein interaction modules that play an essential role in cellular processes and signaling networks and are often involved in diseases like cancer. The ELM database is a collection of manually curated motif knowledge from scientific papers. It has become a crucial resource for investigating motif biology and recognizing candidate ELMs in novel amino acid sequences. Users can search amino acid sequences or UniProt Accessions on the ELM resource web interface. However, as with many web services, there are limitations in the swift processing of large-scale queries through the ELM web interface or API calls, and, therefore, integration into protein function analysis pipelines is limited. Results To allow swift, large-scale motif analyses on protein sequences using ELMs curated in the ELM database, we have extended the gget suite of Python and command line tools with a new module, gget elm , which does not rely on the ELM server for efficiently finding candidate ELMs in user-submitted amino acid sequences and UniProt Accessions. gget elm increases accessibility to the information stored in the ELM database and allows scalable searches for motif-mediated interaction sites in the amino acid sequences. Availability and implementation The manual and source code are available at https://github.com/pachterlab/gget. [ABSTRACT FROM AUTHOR]

In proteomics, a crucial aspect is to identify peptide sequences. De novo sequencing methods have been widely employed to identify peptide sequences, and numerous tools have been proposed over the past two decades. Recently, deep learning approaches have been introduced for de novo sequencing. Previous methods focused on encoding tandem mass spectra and predicting peptide sequences from the first amino acid onwards. However, when predicting peptides using tandem mass spectra, the peptide sequence can be predicted not only from the first amino acid but also from the last amino acid due to the coexistence of b-ion (or a- or c-ion) and y-ion (or x- or z-ion) fragments in the tandem mass spectra. Therefore, it is essential to predict peptide sequences bidirectionally. Our approach, called NovoB, utilizes a Transformer model to predict peptide sequences bidirectionally, starting with both the first and last amino acids. In comparison to Casanovo, our method achieved an improvement of the average peptide-level accuracy rate of approximately 9.8% across all species. Author summary: Understanding the characteristics of data is very important in deep learning methods. When predicting sentences, the transformer model naturally predicts from the first word. For this reason, previous methods predicted peptide sequences from the first amino acid. However, in tandem mass spectra, it is possible to predict peptide sequences bidirectionally. This method shows better results than previous approaches because it can better encode tandem mass spectra. We have demonstrated that good results can be achieved simply by understanding the characteristics of such data and using the model appropriately. We hope that this paper will help various readers improve the performance capabilities of their models. Furthermore, given that bidirectional peptide sequence prediction is crucial in de novo peptide sequence analysis, we hope that this approach will be applied to both existing and future methods utilizing deep learning techniques. [ABSTRACT FROM AUTHOR]

Interpreting biological phenomena at the molecular and cellular levels reveals the ways in which information that is specific to living organisms is processed: from the genetic record contained in a strand of DNA, to the translation process, and then to the construction of proteins that carry the flow and processing of information as well as reveal evolutionary mechanisms. The processing of a surprisingly small amount of information, i.e., in the range of 1 GB, contains the record of human DNA that is used in the construction of the highly complex system that is the human body. This shows that what is important is not the quantity of information but rather its skillful use—in other words, this facilitates proper processing. This paper describes the quantitative relations that characterize information during the successive steps of the "biological dogma", illustrating a transition from the recording of information in a DNA strand to the production of proteins exhibiting a defined specificity. It is this that is encoded in the form of information and that determines the unique activity, i.e., the measure of a protein's "intelligence". In a situation of information deficit at the transformation stage of a primary protein structure to a tertiary or quaternary structure, a particular role is served by the environment as a supplier of complementary information, thus leading to the achievement of a structure that guarantees the fulfillment of a specified function. Its quantitative evaluation is possible via using a "fuzzy oil drop" (FOD), particularly with respect to its modified version. This can be achieved when taking into account the participation of an environment other than water in the construction of a specific 3D structure (FOD-M). The next step of information processing on the higher organizational level is the construction of the proteome, where the interrelationship between different functional tasks and organism requirements can be generally characterized by homeostasis. An open system that maintains the stability of all components can be achieved exclusively in a condition of automatic control that is realized by negative feedback loops. This suggests a hypothesis of proteome construction that is based on the system of negative feedback loops. The purpose of this paper is the analysis of information flow in organisms with a particular emphasis on the role of proteins in this process. This paper also presents a model introducing the component of changed conditions and its influence on the protein folding process—since the specificity of proteins is coded in their structure. [ABSTRACT FROM AUTHOR]

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]

Background: The Drug–Target Interaction (DTI) prediction uses a drug molecule and a protein sequence as inputs to predict the binding affinity value. In recent years, deep learning-based models have gotten more attention. These methods have two modules: the feature extraction module and the task prediction module. In most deep learning-based approaches, a simple task prediction loss (i.e., categorical cross entropy for the classification task and mean squared error for the regression task) is used to learn the model. In machine learning, contrastive-based loss functions are developed to learn more discriminative feature space. In a deep learning-based model, extracting more discriminative feature space leads to performance improvement for the task prediction module. Results: In this paper, we have used multimodal knowledge as input and proposed an attention-based fusion technique to combine this knowledge. Also, we investigate how utilizing contrastive loss function along the task prediction loss could help the approach to learn a more powerful model. Four contrastive loss functions are considered: (1) max-margin contrastive loss function, (2) triplet loss function, (3) Multi-class N-pair Loss Objective, and (4) NT-Xent loss function. The proposed model is evaluated using four well-known datasets: Wang et al. dataset, Luo's dataset, Davis, and KIBA datasets. Conclusions: Accordingly, after reviewing the state-of-the-art methods, we developed a multimodal feature extraction network by combining protein sequences and drug molecules, along with protein–protein interaction networks and drug–drug interaction networks. The results show it performs significantly better than the comparable state-of-the-art approaches. [ABSTRACT FROM AUTHOR]

Background: Lysine glutarylation (Kglu) is one of the most important Post-translational modifications (PTMs), which plays significant roles in various cellular functions, including metabolism, mitochondrial processes, and translation. Therefore, accurate identification of the Kglu site is important for elucidating protein molecular function. Due to the time-consuming and expensive limitations of traditional biological experiments, computational-based Kglu site prediction research is gaining more and more attention. Results: In this paper, we proposed GBDT_KgluSite, a novel Kglu site prediction model based on GBDT and appropriate feature combinations, which achieved satisfactory performance. Specifically, seven features including sequence-based features, physicochemical property-based features, structural-based features, and evolutionary-derived features were used to characterize proteins. NearMiss-3 and Elastic Net were applied to address data imbalance and feature redundancy issues, respectively. The experimental results show that GBDT_KgluSite has good robustness and generalization ability, with accuracy and AUC values of 93.73%, and 98.14% on five-fold cross-validation as well as 90.11%, and 96.75% on the independent test dataset, respectively. Conclusion: GBDT_KgluSite is an effective computational method for identifying Kglu sites in protein sequences. It has good stability and generalization ability and could be useful for the identification of new Kglu sites in the future. The relevant code and dataset are available at https://github.com/flyinsky6/GBDT_KgluSite. [ABSTRACT FROM AUTHOR]

LncRNA–protein interactions are ubiquitous in organisms and play a crucial role in a variety of biological processes and complex diseases. Many computational methods have been reported for lncRNA–protein interaction prediction. However, the experimental techniques to detect lncRNA–protein interactions are laborious and time-consuming. Therefore, to address this challenge, this paper proposes a reweighting boosting feature selection (RBFS) method model to select key features. Specially, a reweighted apporach can adjust the contribution of each observational samples to learning model fitting; let higher weights are given more influence samples than those with lower weights. Feature selection with boosting can efficiently rank to iterate over important features to obtain the optimal feature subset. Besides, in the experiments, the RBFS method is applied to the prediction of lncRNA–protein interactions. The experimental results demonstrate that our method achieves higher accuracy and less redundancy with fewer features. [ABSTRACT FROM AUTHOR]

Protein structure prediction represents a significant challenge in the field of bioinformatics, with the prediction of protein structures using backbone dihedral angles recently achieving significant progress due to the rise of deep neural network research. However, there is a trend in protein structure prediction research to employ increasingly complex neural networks and contributions from multiple models. This study, on the other hand, explores how a single model transparently behaves using sequence data only and what can be expected from the predicted angles. To this end, the current paper presents data acquisition, deep learning model definition, and training toward the final protein backbone angle prediction. The method applies a simple fully connected neural network (FCNN) model that takes only the primary structure of the protein with a sliding window of size 21 as input to predict protein backbone ϕ and ψ dihedral angles. Despite its simplicity, the model shows surprising accuracy for the ϕ angle prediction and somewhat lower accuracy for the ψ angle prediction. Moreover, this study demonstrates that protein secondary structure prediction is also possible with simple neural networks that take in only the protein amino-acid residue sequence, but more complex models are required for higher accuracies. [ABSTRACT FROM AUTHOR]

Simple Summary: It is widely accepted that climate change is a rapidly increasing threat to the aquaculture industry, with the biggest issue being a slow but continuous rise in water temperature and an associated decrease in dissolved oxygen in some areas. In New Zealand fish farms, during the hottest period of the year (January/February), increased water temperatures have been shown to be detrimental to Chinook salmon performance and are associated with mortalities. We report on the responses of Chinook salmon to a chronic, long-term increase in temperature and decrease in oxygen over a three-month period. We also identified, characterised, and developed relevant biomarkers for Chinook salmon and applied them to monitor responses to changes in the fish's environment. This paper provides evidence on how the use of molecular tools, alongside general health and haematological parameters, have the potential to better inform the aquaculture industry of stock health status, thereby enhancing resilience, efficiency, and productivity. In New Zealand, during the hottest periods of the year, some salmon farms in the Marlborough Sounds reach water temperatures above the optimal range for Chinook salmon. High levels of mortality are recorded during these periods, emphasising the importance of understanding thermal stress in this species. In this study, the responses of Chinook salmon (Oncorhynchus tshawytscha) to chronic, long-term changes in temperature and dissolved oxygen were investigated. This is a unique investigation due to the duration of the stress events the fish were exposed to. Health and haematological parameters were analysed alongside gene expression results to determine the effects of thermal stress on Chinook salmon. Six copies of heat shock protein 90 (HSP90) were discovered and characterised: HSP90AA1.1a, HSP90AA1.2a, HSP90AA1.1b, HSP90AA1.2b, HSP90AB1a and HSP90AB1b, as well as two copies of SOD1, named SOD1a and SOD1b. The amino acid sequences contained features similar to those found in other vertebrate HSP90 and SOD1 sequences, and the phylogenetic tree and synteny analysis provided conclusive evidence of their relationship to other vertebrate HSP90 and SOD1 genes. Primers were designed for qPCR to enable the expression of all copies of HSP90 and SOD1 to be analysed. The expression studies showed that HSP90 and SOD1 were downregulated in the liver and spleen in response to longer term exposure to high temperatures and lower dissolved oxygen. HSP90 was also downregulated in the gill; however, the results for SOD1 expression in the gill were not conclusive. This study provides important insights into the physiological and genetic responses of Chinook salmon to temperature and oxygen stress, which are critical for developing sustainable fish aquaculture in an era of changing global climates. [ABSTRACT FROM AUTHOR]

Background: Obtaining accurate drug–target binding affinity (DTA) information is significant for drug discovery and drug repositioning. Although some methods have been proposed for predicting DTA, the features of proteins and drugs still need to be further analyzed. Recently, deep learning has been successfully used in many fields. Hence, designing a more effective deep learning method for predicting DTA remains attractive. Results: Dynamic graph DTA (DGDTA), which uses a dynamic graph attention network combined with a bidirectional long short-term memory (Bi-LSTM) network to predict DTA is proposed in this paper. DGDTA adopts drug compound as input according to its corresponding simplified molecular input line entry system (SMILES) and protein amino acid sequence. First, each drug is considered a graph of interactions between atoms and edges, and dynamic attention scores are used to consider which atoms and edges in the drug are most important for predicting DTA. Then, Bi-LSTM is used to better extract the contextual information features of protein amino acid sequences. Finally, after combining the obtained drug and protein feature vectors, the DTA is predicted by a fully connected layer. The source code is available from GitHub at https://github.com/luojunwei/DGDTA. Conclusions: The experimental results show that DGDTA can predict DTA more accurately than some other methods. [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]

Background: Drug-target binding affinity (DTA) prediction is important for the rapid development of drug discovery. Compared to traditional methods, deep learning methods provide a new way for DTA prediction to achieve good performance without much knowledge of the biochemical background. However, there are still room for improvement in DTA prediction: (1) only focusing on the information of the atom leads to an incomplete representation of the molecular graph; (2) the self-supervised learning method could be introduced for protein representation. Results: In this paper, a DTA prediction model using the deep learning method is proposed, which uses an undirected-CMPNN for molecular embedding and combines CPCProt and MLM models for protein embedding. An attention mechanism is introduced to discover the important part of the protein sequence. The proposed method is evaluated on the datasets Ki and Davis, and the model outperformed other deep learning methods. Conclusions: The proposed model improves the performance of the DTA prediction, which provides a novel strategy for deep learning-based virtual screening methods. [ABSTRACT FROM AUTHOR]

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]

Background: Drug–target affinity (DTA) prediction is a critical step in the field of drug discovery. In recent years, deep learning-based methods have emerged for DTA prediction. In order to solve the problem of fusion of substructure information of drug molecular graphs and utilize multi-scale information of protein, a self-supervised pre-training model based on substructure extraction and multi-scale features is proposed in this paper. Results: For drug molecules, the model obtains substructure information through the method of probability matrix, and the contrastive learning method is implemented on the graph-level representation and subgraph-level representation to pre-train the graph encoder for downstream tasks. For targets, a BiLSTM method that integrates multi-scale features is used to capture long-distance relationships in the amino acid sequence. The experimental results showed that our model achieved better performance for DTA prediction. Conclusions: The proposed model improves the performance of the DTA prediction, which provides a novel strategy based on substructure extraction and multi-scale features. [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]

Simple Summary: As COVID-19 shifts from pandemic to endemic, emerging variants may be more or less virulent. Predicting whether an emerging COVID-19 variant has of high risk of causing severe disease is needed to plan for potential burdens on hospital capacity and protecting vulnerable populations. However, it takes time to do laboratory and animal experiments to determine whether a new genetic variant might be more severe, and the results may not be representative of when the virus infects humans. By the time there is epidemiological data on the severity of disease associated with a new variant, it can be too late for designing an optimal public health response. There is a critical need for computer models that can predict severe disease risk from genetic sequence data, which can be obtained from just the first few infections in a potential incoming wave. Two key challenges make computer modeling difficult: (1) sequence changes are complex, and (2) using historical data to predict future disease requires accounting for the confounding effects of changing patient demographics, improving therapeutics, and increased vaccination. In this paper, we introduce a novel interpretable deep learning architecture to solve this problem, demonstrating that it can make robust predictions for emerging variants. Through the COVID-19 pandemic, SARS-CoV-2 has gained and lost multiple mutations in novel or unexpected combinations. Predicting how complex mutations affect COVID-19 disease severity is critical in planning public health responses as the virus continues to evolve. This paper presents a novel computational framework to complement conventional lineage classification and applies it to predict the severe disease potential of viral genetic variation. The transformer-based neural network model architecture has additional layers that provide sample embeddings and sequence-wide attention for interpretation and visualization. First, training a model to predict SARS-CoV-2 taxonomy validates the architecture's interpretability. Second, an interpretable predictive model of disease severity is trained on spike protein sequence and patient metadata from GISAID. Confounding effects of changing patient demographics, increasing vaccination rates, and improving treatment over time are addressed by including demographics and case date as independent input to the neural network model. The resulting model can be interpreted to identify potentially significant virus mutations and proves to be a robust predctive tool. Although trained on sequence data obtained entirely before the availability of empirical data for Omicron, the model can predict the Omicron's reduced risk of severe disease, in accord with epidemiological and experimental data. [ABSTRACT FROM AUTHOR]

Background: The stability of protein sequence structure plays an important role in the prevention and treatment of diseases. Results: In this paper, particle swarm optimization and tabu search are combined to propose a new method for protein structure prediction. The experimental results show that: for four groups of artificial protein sequences with different lengths, this method obtains the lowest potential energy value and stable structure prediction results, and the effect is obviously better than the other two comparison methods. Taking the first group of protein sequences as an example, our method improves the prediction of minimum potential energy by 127% and 7% respectively. Conclusions: Therefore, the method proposed in this paper is more suitable for the prediction of protein structural stability. [ABSTRACT FROM AUTHOR]

Protein sequence similarity analysis (PSSA) is a significant task in bioinformatics, which can obtain information about unknown sequences such as protein structures and homology relationships. Protein sequence refers to the series of amino acids with rich physical and chemical properties, namely the basic structure of proteins. However, sequence similarity analysis and phylogenetic analysis between different species which have complex amino acid sequences is a challenging problem. In this paper, nine properties of amino acids were considered and the sequence was converted into numerical values by principal component analysis (PCA); with Haar Wavelet Transform, and Higuchi fractal dimension (HFD), a new feature vector is constructed to represent the sequence; Spearman distance was selected to calculate the distance matrix and the phylogenetic tree was constructed. In this paper, two representative protein sequences (9 ND5 (NADH dehydrogenase 5) and 8 ND6 (NADH dehydrogenase 6)) were selected for similarity analysis and phylogenetic analysis, and compared with MEGA software and other existing methods. The extensive results show that our method is outperforming and results consistent with the known facts. [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]

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]

In recent years, protein structure problems have become a hotspot for understanding protein folding and function mechanisms. It has been observed that most of the protein structure works rely on and benefit from co-evolutionary information obtained by multiple sequence alignment (MSA). As an example, AlphaFold2 (AF2) is a typical MSA-based protein structure tool which is famous for its high accuracy. As a consequence, these MSA-based methods are limited by the quality of the MSAs. Especially for orphan proteins that have no homologous sequence, AlphaFold2 performs unsatisfactorily as MSA depth decreases, which may pose a barrier to its widespread application in protein mutation and design problems in which there are no rich homologous sequences and rapid prediction is needed. In this paper, we constructed two standard datasets for orphan and de novo proteins which have insufficient/none homology information, called Orphan62 and Design204, respectively, to fairly evaluate the performance of the various methods in this case. Then, depending on whether or not utilizing scarce MSA information, we summarized two approaches, MSA-enhanced and MSA-free methods, to effectively solve the issue without sufficient MSAs. MSA-enhanced model aims to improve poor MSA quality from the data source by knowledge distillation and generation models. MSA-free model directly learns the relationship between residues on enormous protein sequences from pre-trained models, bypassing the step of extracting the residue pair representation from MSA. Next, we evaluated the performance of four MSA-free methods (trRosettaX-Single, TRFold, ESMFold and ProtT5) and MSA-enhanced (Bagging MSA) method compared with a traditional MSA-based method AlphaFold2, in two protein structure-related prediction tasks, respectively. Comparison analyses show that trRosettaX-Single and ESMFold which belong to MSA-free method can achieve fast prediction (⁠|$\sim\! 40$| s) and comparable performance compared with AF2 in tertiary structure prediction, especially for short peptides, |$\alpha $| -helical segments and targets with few homologous sequences. Bagging MSA utilizing MSA enhancement improves the accuracy of our trained base model which is an MSA-based method when poor homology information exists in secondary structure prediction. Our study provides biologists an insight of how to select rapid and appropriate prediction tools for enzyme engineering and peptide drug development. Contact guofei@csu.edu.cn , jj.tang@siat.ac.cn [ABSTRACT FROM AUTHOR]

Classifying epitopes is essential since they can be applied in various fields, including therapeutics, diagnostics and peptide-based vaccines. To determine the epitope or peptide against an antibody, epitope mapping with peptides is the most extensively used method. However, this method is more time-consuming and inefficient than using present methods. The ability to retrieve data on protein sequences through laboratory procedures has led to the development of computational models that predict epitope binding based on machine learning and deep learning (DL). It has also evolved to become a crucial part of developing effective cancer immunotherapies. This paper proposes an architecture to generalize this case since various research strives to solve a low-performance classification problem. A proposed DL model is the fusion architecture, which combines two architectures: Transformer architecture and convolutional neural network (CNN), called MITNet and MITNet-Fusion. Combining these two architectures enriches feature space to correlate epitope labels with the binary classification method. The selected epitope–T-cell receptor (TCR) interactions are GILG, GLCT and NLVP, acquired from three databases: IEDB, VDJdb and McPAS-TCR. The previous input data was extracted using amino acid composition, dipeptide composition, spectrum descriptor and the combination of all those features called AADIP composition to encode the input data to DL architecture. For ensuring consistency, fivefold cross-validations were performed using the area under curve metric. Results showed that GILG, GLCT and NLVP received scores of 0.85, 0.87 and 0.86, respectively. Those results were compared to prior architecture and outperformed other similar deep learning models. [ABSTRACT FROM AUTHOR]

Spider venoms are composed, among other substances, of peptide toxins whose selectivity for certain physiological targets has made them powerful tools for applications such as bioinsecticides, analgesics, antiarrhythmics, antibacterials, antifungals and antimalarials, among others. Bioinsecticides are an environmentally friendly alternative to conventional agrochemicals. In this paper, the primary structure of an insecticidal peptide was obtained from the venom gland transcriptome of the ctenid spider Phoneutria depilata (Transcript ID PhdNtxNav24). The peptide contains 53 amino acids, including 10 Cys residues that form 5 disulfide bonds. Using the amino acid sequence of such peptide, a synthetic gene was constructed de novo by overlapping PCRs and cloned into an expression vector. A recombinant peptide, named delta-ctenitoxin (rCtx-4), was obtained. It was expressed, folded, purified and validated using mass spectrometry (7994.61 Da). The insecticidal activity of rCtx-4 was demonstrated through intrathoracic injection in crickets (LD50 1.2 μg/g insect) and it was not toxic to mice. rCtx-4 is a potential bioinsecticide that could have a broad spectrum of applications in agriculture. [ABSTRACT FROM AUTHOR]

Owing to the rapid changes in the antigenicity of influenza viruses, it is difficult for humans to obtain lasting immunity through antiviral therapy. Hence, tracking the dynamic changes in the antigenicity of influenza viruses can provide a basis for vaccines and drug treatments to cope with the spread of influenza viruses. In this paper, we developed a novel quantitative prediction method to predict the antigenic distance between virus strains using attribute network embedding techniques. An antigenic network is built to model and combine the genetic and antigenic characteristics of the influenza A virus H3N2, using the continuous distributed representation of the virus strain protein sequence (ProtVec) as a node attribute and the antigenic distance between virus strains as an edge weight. The results show a strong positive correlation between supplementing genetic features and antigenic distance prediction accuracy. Further analysis indicates that our prediction model can comprehensively and accurately track the differences in antigenic distances between vaccines and influenza virus strains, and it outperforms existing methods in predicting antigenic distances between strains. [ABSTRACT FROM AUTHOR]

In this paper, the influence of external terahertz electromagnetic fields with different frequencies of 4 THz, 10 THz, 15 THz, and 20 THz on the permeability of the Kv1.2 voltage-gated potassium ion channel on the nerve cell membrane was studied using the combined model of the "Constant Electric Field-Ion Imbalance" method by molecular dynamics. We found that although the applied terahertz electric field does not produce strong resonance with the –C=O groups of the conservative sequence T-V-G-Y-G amino acid residue of the selective filter (SF) of the channel, it would affect the stability of the electrostatic bond between potassium ions and the carbonyl group of T-V-G-Y-G of SF, and it would affect the stability of the hydrogen bond between water molecules and oxygen atoms of the hydroxyl group of the 374THR side chain at the SF entrance, changing the potential and occupied states of ions in the SF and the occurrence probability of the permeation mode of ions and resulting in the change in the permeability of the channel. Compared with no external electric field, when the external electric field with 15 THz frequency is applied, the lifetime of the hydrogen bond is reduced by 29%, the probability of the "soft knock on" mode is decreased by 46.9%, and the ion flux of the channel is activated by 67.7%. Our research results support the view that compared to "direct knock-on", "soft knock-on" is a slower permeation mode. [ABSTRACT FROM AUTHOR]

The Major Histocompatibility Complex (MHC) binds to the derived peptides from pathogens to present them to killer T cells on the cell surface. Developing computational methods for accurate, fast, and explainable peptide-MHC binding prediction can facilitate immunotherapies and vaccine development. Various deep learning-based methods rely on separate feature extraction from the peptide and MHC sequences and ignore their pairwise binding information. This paper develops a capsule neural network-based method to efficiently capture the peptide-MHC complex features to predict the peptide-MHC class I binding. Various evaluations confirmed our method outperformance over the alternative methods, while it can provide accurate prediction over less available data. Moreover, for providing precise insights into the results, we explored the essential features that contributed to the prediction. Since the simulation results demonstrated consistency with the experimental studies, we concluded that our method can be utilized for the accurate, rapid, and interpretable peptide-MHC binding prediction to assist biological therapies. A capsule neural network-based method captures pairwise information to accurately predict MHC-I and peptide binding. The comprehensive test of the model indicated both its superiority in performance and interpretability. [ABSTRACT FROM AUTHOR]

Simple Summary: This article is mainly about that molecular characterization and expression pattern of leptin in yellow cheek carp (Elopichthys bambusa) and its transcriptional changes in response to fasting and refeeding. In this study, the authors used PCR to clone the CDS of leptin in yellow cheek carp, and analyzed the sequence differences of the gene with other species, constructed the phylogenetic tree, used real-time PCR for analyzing the expression of leptin in different tissues, including the expression of leptin in the brain and liver after fasting–refeeding of yellow cheek carp. This paper found that the full-length cDNA sequence of Eblep was 1140 bp and the length of the open reading frame (ORF), which can encode a protein of 174 amino acids, was 525 bp. The Eblep mRNA transcript was detected in all tested tissues, with the highest expression in the liver and lowest expression in the spleen. It was found that the change in the mRNA expression of EbLep may be an adaptive strategy for different energy levels by studying the expression of EbLep mRNA in the brain and liver under fasting and refeeding. Leptin, a secretory protein encoded by obese genes, plays an important role in regulating feeding and energy metabolism in fish. To study the structure and function of the Leptin gene in yellow cheek carp (Elopichthys bambusa), the full-length cDNA sequence of leptin was cloned, named EbLep. The full-length cDNA of Eblep was 1140 bp, and the length of the open reading frame (ORF), which can encode a protein of 174 amino acids, was 525 bp. The signal peptide was predicted to contain 33 amino acids. Sequence alignment showed that the amino acid sequence of Leptin was conserved in cyprinid fish. Despite large differences between primary structures, the tertiary structure of the EbLep protein was similar to that of the human protein and had four α-helices. The EbLep mRNA transcript was detected in all tested tissues, with the highest expression in the liver and lowest expression in the spleen. In this study, short-term fasting significantly increased the mRNA expression of EbLep in the liver, which returned to a normal level after 6 days of refeeding and was significantly lower than the normal level after 28 days of refeeding. In the brain, the mRNA expression of EbLep significantly decreased during short-term fasting and significantly increased to a higher value than the control group after 1 h of refeeding. It then rapidly decreased to a lower value than the control group after 6 h of refeeding, returning to the normal level after 1 day of refeeding, and significantly decreasing to a lower value than the control group after 28 days of refeeding. To sum up, the change in the mRNA expression of EbLep in the brain and liver may be an adaptive strategy for different energy levels. [ABSTRACT FROM AUTHOR]

As biological macromolecules, proteins are involved in important cellular functions ranging from DNA replication and biosynthesis to metabolic signalling and environmental sensing. Protein sequencing can help understand the relationship between protein function and structure, and provide key information for disease diagnosis and new drug design. Nanopore sensors are a novel technology to achieve the goal of label‐free and high‐throughput protein sequencing. In recent years, nanopore‐based biosensors have been widely used in the detection and analysis of biomolecules such as DNA, RNA, and proteins. At the same time, computer simulations can describe the transport of proteins through nanopores at the atomic level. This paper reviews the applications of nanopore sensors in protein sequencing over the past decade and the solutions to key problems from a computer simulation perspective, with the aim of pointing the way to the future of nanopore protein sequencing. [ABSTRACT FROM AUTHOR]

Simple Summary: Accurate identification of potential targets for drugs to interact with can accelerate drug development. The identification of drug–target interactions can provide insights into hidden drug efficacy. This paper presents a prediction model based on feature similarity fusion that can identify crucial features of drugs and targets to help predict drug–target interactions. The key to drug discovery is the identification of a target and a corresponding drug compound. Effective identification of drug–target interactions facilitates the development of drug discovery. In this paper, drug similarity and target similarity are considered, and graphical representations are used to extract internal structural information and intermolecular interaction information about drugs and targets. First, drug similarity and target similarity are fused using the similarity network fusion (SNF) method. Then, the graph isomorphic network (GIN) is used to extract the features with information about the internal structure of drug molecules. For target proteins, feature extraction is carried out using TextCNN to efficiently capture the features of target protein sequences. Three different divisions (CVD, CVP, CVT) are used on the standard dataset, and experiments are carried out separately to validate the performance of the model for drug–target interaction prediction. The experimental results show that our method achieves better results on AUC and AUPR. The docking results also show the superiority of the proposed model in predicting drug–target interactions. [ABSTRACT FROM AUTHOR]

The paper contains a discussion on mathematical modifying and redesigning DNA with the use of Markov chains. We give a simple mathematical technique for overwriting missing parts of DNA. With a certain probability (without even knowing the function of the missing codon) we can find a synonymous codon, so that there is no frequency change in amino acid sequences of proteins. We use Markov Chain to analyze the dependencies in DNA sequence of the human gene Alpha 1,3-Galactosyltransfe rase 2. We include a theoretical introduction which facilitates the understanding of the paper for non-mathematicians, especially for biologists not familiar with the theory of Markov chains. [ABSTRACT FROM AUTHOR]

Background: In protein sequences—as there are 61 sense codons but only 20 standard amino acids—most amino acids are encoded by more than one codon. Although such synonymous codons do not alter the encoded amino acid sequence, their selection can dramatically affect the expression of the resulting protein. Codon optimization of synthetic DNA sequences is important for heterologous expression. However, existing solutions are primarily based on choosing high-frequency codons only, neglecting the important effects of rare codons. In this paper, we propose a novel recurrent-neural-network based codon optimization tool, ICOR, that aims to learn codon usage bias on a genomic dataset of Escherichia coli. We compile a dataset of over 7,000 non-redundant, high-expression, robust genes which are used for deep learning. The model uses a bidirectional long short-term memory-based architecture, allowing for the sequential context of codon usage in genes to be learned. Our tool can predict synonymous codons for synthetic genes toward optimal expression in Escherichia coli. Results: We demonstrate that sequential context achieved via RNN may yield codon selection that is more similar to the host genome. Based on computational metrics that predict protein expression, ICOR theoretically optimizes protein expression more than frequency-based approaches. ICOR is evaluated on 1,481 Escherichia coli genes as well as a benchmark set of 40 select DNA sequences whose heterologous expression has been previously characterized. ICOR's performance is measured across five metrics: the Codon Adaptation Index, GC-content, negative repeat elements, negative cis-regulatory elements, and codon frequency distribution. Conclusions: The results, based on in silico metrics, indicate that ICOR codon optimization is theoretically more effective in enhancing recombinant expression of proteins over other established codon optimization techniques. Our tool is provided as an open-source software package that includes the benchmark set of sequences used in this study. [ABSTRACT FROM AUTHOR]

Drug-target binding affinity prediction plays a key role in the early stage of drug discovery. Numerous experimental and data-driven approaches have been developed for predicting drug-target binding affinity. However, experimental methods highly rely on the limited structural-related information from drug-target pairs, domain knowledge, and time-consuming assays. On the other hand, learning-based methods have shown an acceptable prediction performance. However, most of them utilize several simple and complex types of proteins and drug compounds data, ranging from the protein sequences to the topology of a graph representation of drug compounds, employing multiple deep neural networks for encoding and feature extraction, and so, leads to the computational overheads. In this study, we propose a unified measure for protein sequence encoding, named BiComp, which provides compression-based and evolutionary-related features from the protein sequences. Specifically, we employ Normalized Compression Distance and Smith-Waterman measures for capturing complementary information from the algorithmic information theory and biological domains, respectively. We utilize the proposed measure to encode the input proteins feeding a new deep neural network-based method for drug-target binding affinity prediction, named BiComp-DTA. BiComp-DTA is evaluated utilizing four benchmark datasets for drug-target binding affinity prediction. Compared to the state-of-the-art methods, which employ complex models for protein encoding and feature extraction, BiComp-DTA provides superior efficiency in terms of accuracy, runtime, and the number of trainable parameters. The latter achievement facilitates execution of BiComp-DTA on a normal desktop computer in a fast fashion. As a comparative study, we evaluate BiComp's efficiency against its components for drug-target binding affinity prediction. The results have shown superior accuracy of BiComp due to the orthogonality and complementary nature of Smith-Waterman and Normalized Compression Distance measures for protein sequences. Such a protein sequence encoding provides efficient representation with no need for multiple sources of information, deep domain knowledge, and complex neural networks. Author summary: Drugs alter the protein conformation and functionality through binding to the target proteins. Measuring the binding tendency between drugs and proteins, known as drug-target binding affinity prediction, is an important task in developing new drug candidates and novel therapeutics in the early stage of drug discovery. Measuring the affinity strength through biological assays is expensive and time-consuming. On the other hand, considerable data-driven and computational methods are proposed to predict binding affinity prediction through employing complicated models, and multiple types of proteins and drug compounds data. In this paper, we present a deep learning-based method, named BiComp-DTA, which utilizes both biological and sequence content features extracted from the protein sequences using a new unified measure, named BiComp. Encoded protein sequences and drug's information are passed to the fully-connected and hybrid neural networks, including CNN and separable CNN layers, respectively. Simulation results using benchmark datasets indicate that BiComp-DTA provides an accurate and fast prediction for drug-target binding affinity prediction without employing multiple complicated networks, several types of proteins and drug compounds data, and external tools for input data modeling. [ABSTRACT FROM AUTHOR]

Drug discovery for a protein target is a laborious and costly process. Deep learning (DL) methods have been applied to drug discovery and successfully generated novel molecular structures, and they can substantially reduce development time and costs. However, most of them rely on prior knowledge, either by drawing on the structure and properties of known molecules to generate similar candidate molecules or extracting information on the binding sites of protein pockets to obtain molecules that can bind to them. In this paper, DeepTarget, an end-to-end DL model, was proposed to generate novel molecules solely relying on the amino acid sequence of the target protein to reduce the heavy reliance on prior knowledge. DeepTarget includes three modules: Amino Acid Sequence Embedding (AASE), Structural Feature Inference (SFI), and Molecule Generation (MG). AASE generates embeddings from the amino acid sequence of the target protein. SFI inferences the potential structural features of the synthesized molecule, and MG seeks to construct the eventual molecule. The validity of the generated molecules was demonstrated by a benchmark platform of molecular generation models. The interaction between the generated molecules and the target proteins was also verified on the basis of two metrics, drug–target affinity and molecular docking. The results of the experiments indicated the efficacy of the model for direct molecule generation solely conditioned on amino acid sequence. [ABSTRACT FROM AUTHOR]

The formation of mature mRNA requires cutting introns and splicing exons. The occurrence of splicing involves the participation of the spliceosome. Common spliceosomes mainly include five snRNPs: U1, U2, U4/U6, and U5. SF3a2, an essential component of spliceosome U2 snRNP, participates in splicing a series of genes. There is no definition of SF3a2 in plants. The paper elaborated on SF3a2s from a series of plants through protein sequence similarity. We constructed the evolutionary relationship of SF3a2s in plants. Moreover, we analyzed the similarities and differences in gene structure, protein structure, the cis-element of the promoter, and expression pattern; we predicted their interacting proteins and constructed their collinearity. We have preliminarily analyzed SF3a2s in plants and clarified the evolutionary relationship between different species; these studies can better serve for in-depth research on the members of the spliceosome in plants. [ABSTRACT FROM AUTHOR]

Directed laboratory evolution applies iterative rounds of mutation and selection to explore the protein fitness landscape and provides rich information regarding the underlying relationships between protein sequence, structure, and function. Laboratory evolution data consist of protein sequences sampled from evolving populations over multiple generations and this data type does not fit into established supervised and unsupervised machine learning approaches. We develop a statistical learning framework that models the evolutionary process and can infer the protein fitness landscape from multiple snapshots along an evolutionary trajectory. We apply our modeling approach to dihydrofolate reductase (DHFR) laboratory evolution data and the resulting landscape parameters capture important aspects of DHFR structure and function. We use the resulting model to understand the structure of the fitness landscape and find numerous examples of epistasis but an overall global peak that is evolutionarily accessible from most starting sequences. Finally, we use the model to perform an in silico extrapolation of the DHFR laboratory evolution trajectory and computationally design proteins from future evolutionary rounds. Author summary: Laboratory evolution has revolutionized our understanding of protein structure, function, and evolution, and has generated countless useful proteins broad applications in medicine, biocatalysis, and biotechnology. These experiments explore protein sequence space through iterative rounds of mutation and selection and can provide rich data of populations traversing the fitness landscape. In this paper, we present a statistical learning framework that models the evolutionary process and can infer the structure of the underlying protein fitness landscape from multiple snapshots along a laboratory evolution trajectory. We generate a dihydrofolate reductase (DHFR) laboratory evolution data set and apply our modeling approach to infer the landscape parameters. The estimated parameters pinpoint key residues that dictate DHFR structure and function. We use the resulting model to understand the local and global structure of the fitness landscape and to perform in silico directed evolution for protein engineering. [ABSTRACT FROM AUTHOR]

While autoregulative adaptation is a common feature of living tissues, only a few feedback-controlled adaptive biomaterials are available so far. This paper herein reports a new polymer hydrogel platform designed to release anti-inflammatory molecules in response to the inflammatory activation of human blood. In this system, anti-inflammatory peptide drugs, targeting either the complement cascade, a complement receptor, or cyclophilin A, are conjugated to the hydrogel by a peptide sequence that is cleaved by elastase released from activated granulocytes. As a proof of concept, the adaptive drug delivery from the gel triggered by activated granulocytes and the effect of the released drug on the respective inflammatory pathways are demonstrated. Adjusting the gel functionalization degree is shown to allow for tuning the drug release profiles to effective doses within a micromolar range. Feedback-controlled delivery of covalently conjugated drugs from a hydrogel matrix is concluded to provide valuable safety features suitable to equip medical devices with highly active anti-inflammatory agents without suppressing the general immunosurveillance. [ABSTRACT FROM AUTHOR]

Genetic factors affect an individual’s risk of developing obesity, but in most cases each genetic variant has a small effect. Discovery of genes that regulate obesity may provide clues about its underlying biological processes and point to new ways the disease can be treated. Preclinical animal models facilitate genetic discovery in obesity because environmental factors can be better controlled compared with the human population. We studied inbred mouse strains to identify novel genes affecting obesity and glucose metabolism. BTBR T+ Itpr3tf/J (BTBR) mice are fatter and more glucose intolerant than C57BL/6J (B6) mice. Prior genetic studies of these strains identified an obesity locus on chromosome 2. Using congenic mice, we found that obesity was affected by a ∼316 kb region, with only two known genes, pyruvate dehydrogenase kinase 1 (Pdk1) and integrin α 6 (Itga6). Both genes had mutations affecting their amino acid sequence and reducing mRNA levels. Both genes have known functions that could modulate obesity, lipid metabolism, insulin secretion, and/or glucose homeostasis. We hypothesized that genetic variation in or near Pdk1 or Itga6 causing reduced Pdk1 and Itga6 expression would promote obesity and impaired glucose tolerance. We used knockout mice lacking Pdk1 or Itga6 fed an obesigenic diet to test this hypothesis. Under the conditions we studied, we were unable to detect an individual contribution of either Pdk1 or Itga6 to body weight. During our studies, with conditions outside our control, we were unable to reproduce some of our previous body weight data. However, we identified a previously unknown role for Pdk1 in cardiac cholesterol metabolism providing the basis for future investigations. The studies described in this paper highlight the importance and the challenge using physiological outcomes to study obesity genes in mice. [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]

Amyloid aggregates play a major role in diseases as well as in normal physiological function. Receptor-interacting protein kinases 1 and 3 (RIP1/RIP3) aggregates complexes in cellular necroptosis is one example of protein aggregation in normal cellular function. Although recently there have been several studies on full kinase proteins aggregation, the aggregation potential of small peptide sequences of RIP1/RIP3, the physicochemical properties, and the potential in biomedical applications have not been explored. Hence, in this paper, we study the aggregation propensity of peptides consisting of four and twelve amino acid sequences in the RHIM region of RIP1/RIP3 proteins that are known to drive the beta-sheet formation and the subsequent aggregation. The aggregation kinetics, physicochemical characterization, mechanosensitive properties, cellular effects, and potential as a cancer drug depot have been investigated. The results show that the number and concentration of amino acids play a role in amyloid-like aggregates' properties. Further, the aggregates when formulated with cisplatin-induced significant lung cancer cell toxicity compared to an equal amount of cisplatin with and without ultrasound. The study would serve as a platform for further investigation on RIP1/RIP3 peptide and protein aggregates, their role in multiple cellular functions and diseases, and their potential as drug depots. [ABSTRACT FROM AUTHOR]

Multiple sequence alignment approaches refer to algorithmic solutions for the alignment of biological sequences. Since multiple sequence alignment has exponential time complexity when a dynamic programming approach is applied, a substantial number of parallel computing approaches have been implemented in the last two decades to improve their performance. In this paper, we present a systematic literature review of parallel computing approaches applied to multiple sequence alignment algorithms for proteins, published in the open literature from 1988 to 2022; we extracted articles from four scientific databases: ACM Digital Library, IEEE Xplore, Science Direct and SpringerLink, and four journals: Bioinformatics, PLOS Computational Biology, PLOS ONE, and Scientific Reports. Additionally, in order to cover other potential databases and journals, we performed a transversal search through Google Scholar. We conducted a selection process that yielded 106 research articles; then, we analyzed these articles and defined a classification framework. Additionally, we point out some directions and trends for parallel computing approaches for multiple sequence alignment, as well as some unsolved problems. [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]

Helicases are nucleic acid-unwinding enzymes involved in the maintenance of genome integrity. Helicases share several "helicase motifs" that are highly conserved amino acid sequences and are classified into six superfamilies (SFs). The helicase SFs are further grouped into two classes based on their functional units. One class that includes SFs 3-6 functions as a hexamer that can form a ring around DNA. Another class that includes SFs 1 and 2 functions in a non-hexameric form. The high homology in the primary and tertiary structures among SF1 helicases suggests that SF1 helicases have a common underlying mechanism. However, two opposing models for the functional unit, monomer and dimer models, have been proposed to explain DNA unwinding by SF1 helicases. This paper briefly describes the classification of helicase SFs and discusses the structural homology and the two opposing nonhexameric helicase models of SF1 helicases by focusing on Escherichia coli SF1 helicase UvrD, which plays a significant role in both nucleotide-excision repair and methyl-directed mismatch repair. This paper reviews past and recent studies on UvrD, including the author's single-molecule direct visualization of wild-type UvrD and a UvrD mutant lacking the C-terminal 40 amino acids (UvrD-40C), the latter of which was used in genetic and biochemical assays that supported the monomer model. The visualization revealed that multiple UvrD-40C molecules jointly unwind DNA, presumably in an oligomeric form, similar to wild-type UvrD. Therefore, singlemolecule direct visualization of nucleic acid-binding proteins can provide quantitative and kinetic information to reveal their fundamental mechanisms. [ABSTRACT FROM AUTHOR]