BioLiP: a semi-manually curated database for biologically relevant ligand-protein interactions

Protein-peptide binding residue prediction based on protein language models and cross-attention mechanism

Accurate identifications of protein-peptide binding residues are essential for protein-peptide interactions and advancing drug discovery. To address this problem, extensive research efforts have been made to design more discriminative feature representations. However, extracting these explicit features usually depend on third-party tools, resulting in low computational efficacy and suffering from low predictive performance. In this study, we design an end-to-end deep learning-based method, E2EPep, for protein-peptide binding residue prediction using protein sequence only. E2EPep first employs and fine-tunes two state-of-the-art pre-trained protein language models that can extract two different high-latent feature representations from protein sequences relevant for protein structures and functions. A novel feature fusion module is then designed in E2EPep to fuse and optimize the above two feature representations of binding residues. In addition, we have also design E2EPep+, which integrates E2EPep and PepBCL models, to improve the prediction performance. Experimental results on two independent testing data sets demonstrate that E2EPep and E2EPep + could achieve the average AUC values of 0.846 and 0.842 while achieving an average Matthew's correlation coefficient value that is significantly higher than that of existing most of sequence-based methods and comparable to that of the state-of-the-art structure-based predictors. Detailed data analysis shows that the primary strength of E2EPep lies in the effectiveness of feature representation using cross-attention mechanism to fuse the embeddings generated by two fine-tuned protein language models. The standalone package of E2EPep and E2EPep + can be obtained at https://github.com/ckx259/E2EPep.git for academic use only.

Characterization and kinetics of a cathepsin B-inhibiting protein from Musa acuminata Colla peel

Hyperexpression of cathepsin B caused by an imbalance of endogenous inhibitors is involved in multiple pathologies, hence making it a key therapeutic target. Protease inhibitors are effective biomolecules that regulate protease activities and are considered potential therapeutic agents in various diseases. Plant protease inhibitors have been reported as an effective complementary alternative drug. A proteinaceous cathepsin B inhibitor (CBI-BP) has been isolated from Musa acuminata Colla (banana) peel with a molecular weight of 27.9 kDa on SDS-PAGE. The purity of the CBI-BP was confirmed on the native- PAGE. The isolated CBI-BP showed an IC50 value of 8.14 μg and a Ki value of 10.59 μg (0.19 μM). Cathepsin B inhibition kinetics indicated that CBI-BP follows a mixed-type of cathepsin B inhibition. Its inhibition activity was also confirmed by reverse zymography. The inhibitor was stable from pH 2.6-10.0 with maximum activity at pH 7.2, temperature 25-100 °C and exhibited thermostability for 60 min at 70 °C. MALDI/TOF/MS analysis of CBI-BP showed 40 % similarity to the GH18 domain-containing protein (A0A4S8JRM9) from Musa balbisiana. Although in-silico docking studies showed binding of A0A4S8JRM9 to cathepsin B affects the binding energy of the substrate to cathepsin B but is not reported for any anti-cathepsin B activity. This suggests that isolated CBI-BP might be a novel protein with anti-cathepsin B activity. Thus the isolated CBI-BP may be further explored as possible anti-cathepsin B drug.

Predicting Affinity Through Homology (PATH): Interpretable Binding Affinity Prediction with Persistent Homology

Accurate binding affinity prediction is crucial to structure-based drug design. Recent work used computational topology to obtain an effective representation of protein-ligand interactions. While algorithms using algebraic topology have proven useful in predicting properties of biomolecules, previous algorithms employed uninterpretable machine learning models which failed to explain the underlying geometric and topological features that drive accurate binding affinity prediction. Moreover, they had high computational complexity which made them intractable for large proteins. We present the fastest known algorithm to compute persistent homology features for protein-ligand complexes using opposition distance, with a runtime that is independent of the protein size. Then, we exploit these features in a novel, interpretable algorithm to predict protein-ligand binding affinity. Our algorithm achieves interpretability through an effective embedding of distances across bipartite matchings of the protein and ligand atoms into real-valued functions by summing Gaussians centered at features constructed by persistent homology. We name these functions internuclear persistent contours (IPCs) . Next, we introduce persistence fingerprints , a vector with 10 components that sketches the distances of different bipartite matching between protein and ligand atoms, refined from IPCs. Let the number of protein atoms in the protein-ligand complex be n , number of ligand atoms be m , and ω ≈ 2.4 be the matrix multiplication exponent. We show that for any 0 < ε < 1, after an 𝒪 ( mn log( mn )) preprocessing procedure, we can compute an ε -accurate approximation to the persistence fingerprint in 𝒪 ( m log 6 ω ( m/ε )) time, independent of protein size. This is an improvement in time complexity by a factor of 𝒪 (( m + n ) 3 ) over any previous binding affinity prediction that uses persistent homology. We show that the representational power of persistence fingerprint generalizes to protein-ligand binding datasets beyond the training dataset. Then, we introduce PATH , Predicting Affinity Through Homology, a two-part algorithm consisting of PATH + and PATH - . PATH + is an interpretable, small ensemble of shallow regression trees for binding affinity prediction from persistence fingerprints. We show that despite using 1,400-fold fewer features, PATH + has comparable performance to a previous state-of-the-art binding affinity prediction algorithm that uses persistent homology. Moreover, PATH + has the advantage of being interpretable. We visualize the features captured by persistence fingerprint for variant HIV-1 protease complexes and show that persistence fingerprint captures binding-relevant structural mutations. PATH - , in turn, uses regression trees over IPCs to differentiate between binding and decoy complexes. Finally, we benchmarked PATH versus established binding affinity prediction algorithms spanning physics-based, knowledge-based, and deep learning methods, revealing that PATH has comparable or better performance with less overfitting, compared to these state-of-the-art methods. The source code for PATH is released open-source as part of the osprey protein design software package.

DIONYSUS: a database of protein-carbohydrate interfaces

Protein-carbohydrate interactions govern a wide variety of biological processes and play an essential role in the development of different diseases. Here, we present DIONYSUS, the first database of protein-carbohydrate interfaces annotated according to structural, chemical and functional properties of both proteins and carbohydrates. We provide exhaustive information on the nature of interactions, binding site composition, biological function and specific additional information retrieved from existing databases. The user can easily search the database using protein sequence and structure information or by carbohydrate binding site properties. Moreover, for a given interaction site, the user can perform its comparison with a representative subset of non-covalent protein-carbohydrate interactions to retrieve information on its potential function or specificity. Therefore, DIONYSUS is a source of valuable information both for a deeper understanding of general protein-carbohydrate interaction patterns, for annotation of the previously unannotated proteins and for such applications as carbohydrate-based drug design. DIONYSUS is freely available at www.dsimb.inserm.fr/DIONYSUS/.

Enterobacter cloacae-mediated polymer biodegradation: in-silico analysis predicts broad spectrum degradation potential by Alkane monooxygenase

Plastic pollution poses a significant environmental challenge. In this study, the strain Enterobacter cloacae O5-E, a bacterium displaying polyethylene-degrading capabilities was isolated. Over a span of 30 days, analytical techniques including x-ray diffractometry, scanning electron microscopy, optical profilometry, hardness testing and mass spectrometric analysis were employed to examine alterations in the polymer. Results revealed an 11.48% reduction in crystallinity, a 50% decrease in hardness, and a substantial 25-fold increase in surface roughness resulting from the pits and cracks introduced in the polymer by the isolate. Additionally, the presence of degradational by-products revealed via gas chromatography ascertains the steady progression of degradation. Further, recognizing the pivotal role of alkane monooxygenase in plastic degradation, the study expanded to detect this enzyme in the isolate molecularly. Molecular docking studies were conducted to assess the enzyme's affinity with various polymers, demonstrating notable binding capability with most polymers, especially with polyurethane (- 5.47 kcal/mol). These findings highlight the biodegradation potential of Enterobacter cloacae O5-E and the crucial involvement of alkane monooxygenase in the initial steps of the degradation process, offering a promising avenue to address the global plastic pollution crisis.

PMSFF: Improved Protein Binding Residues Prediction through Multi-Scale Sequence-Based Feature Fusion Strategy

Proteins perform different biological functions through binding with various molecules which are mediated by a few key residues and accurate prediction of such protein binding residues (PBRs) is crucial for understanding cellular processes and for designing new drugs. Many computational prediction approaches have been proposed to identify PBRs with sequence-based features. However, these approaches face two main challenges: (1) these methods only concatenate residue feature vectors with a simple sliding window strategy, and (2) it is challenging to find a uniform sliding window size suitable for learning embeddings across different types of PBRs. In this study, we propose one novel framework that could apply multiple types of PBRs Prediciton task through Multi-scale Sequence-based Feature Fusion (PMSFF) strategy. Firstly, PMSFF employs a pre-trained language model named ProtT5, to encode amino acid residues in protein sequences. Then, it generates multi-scale residue embeddings by applying multi-size windows to capture effective neighboring residues and multi-size kernels to learn information across different scales. Additionally, the proposed model treats protein sequences as sentences, employing a bidirectional GRU to learn global context. We also collect benchmark datasets encompassing various PBRs types and evaluate our PMSFF approach to these datasets. Compared with state-of-the-art methods, PMSFF demonstrates superior performance on most PBRs prediction tasks.

Accurately predicting enzyme functions through geometric graph learning on ESMFold-predicted structures

Enzymes are crucial in numerous biological processes, with the Enzyme Commission (EC) number being a commonly used method for defining enzyme function. However, current EC number prediction technologies have not fully recognized the importance of enzyme active sites and structural characteristics. Here, we propose GraphEC, a geometric graph learning-based EC number predictor using the ESMFold-predicted structures and a pre-trained protein language model. Specifically, we first construct a model to predict the enzyme active sites, which is utilized to predict the EC number. The prediction is further improved through a label diffusion algorithm by incorporating homology information. In parallel, the optimum pH of enzymes is predicted to reflect the enzyme-catalyzed reactions. Experiments demonstrate the superior performance of our model in predicting active sites, EC numbers, and optimum pH compared to other state-of-the-art methods. Additional analysis reveals that GraphEC is capable of extracting functional information from protein structures, emphasizing the effectiveness of geometric graph learning. This technology can be used to identify unannotated enzyme functions, as well as to predict their active sites and optimum pH, with the potential to advance research in synthetic biology, genomics, and other fields.

DP-site: A dual deep learning-based method for protein-peptide interaction site prediction

BACKGROUND

Protein-peptide interaction prediction is an important topic for several applications including various biological processes, understanding drug discovery, protein function abnormal cellular behaviors, and treating diseases. Over the years, studies have shown that experimental methods have improved the identification of this bio-molecular interaction. However, predicting protein-peptide interactions using these methods is laborious, time-consuming, dependent on third-party tools, and costly.

METHOD

To address these previous drawbacks, this study introduces a computational framework called DP-Site. The proposed framework concentrates on using a compound of a dual pipeline along with a combination predictor. A deep convolutional neural network for feature extraction and classification is embedded in pipeline 1. In addition, pipeline 2 includes a deep long-short-term memory-based and a random forest classifier for feature extraction and classification. In this investigation, the evolutionary, structure-based, sequence-based, and physicochemical information of proteins is utilized for identifying protein-peptide interaction at the residue level.

RESULTS

The proposed method is evaluated on both the ten-fold cross-validation and independent test sets. The robust and consistent results between cross-validation and independent test sets confirm the ability of the proposed method to predict peptide binding residues in proteins. Moreover, experimental findings demonstrate that DP-Site has significantly outperformed other state-of-the-art sequence-based and structure-based methods. The proposed method achieves a remarkable balance between a specificity of 0.799 and a sensitivity of 0.770, along with the best f-measure of 0.661 and the highest precision of 0.580 using an independent test set.

CONCLUSIONS

The outcome of various experiments confirms the proficiency of the proposed method and outperforms state-of-the-art sequence-based and structure-based methods in terms of the mentioned criteria. DP-Site can be accessed at https://github.com/shafiee 95/shima.shafiee.DP-Site.

Deciphering the Language of Protein-DNA Interactions: A Deep Learning Approach Combining Contextual Embeddings and Multi-Scale Sequence Modeling

Deciphering the mechanisms governing protein-DNA interactions is crucial for understanding key cellular processes and disease pathways. In this work, we present a powerful deep learning approach that significantly advances the computational prediction of DNA-interacting residues from protein sequences. Our method leverages the rich contextual representations learned by pre-trained protein language models, such as ProtTrans, to capture intrinsic biochemical properties and sequence motifs indicative of DNA binding sites. We then integrate these contextual embeddings with a multi-window convolutional neural network architecture, which scans across the sequence at varying window sizes to effectively identify both local and global binding patterns. Comprehensive evaluation on curated benchmark datasets demonstrates the remarkable performance of our approach, achieving an area under the ROC curve (AUC) of 0.89 - a substantial improvement over previous state-of-the-art sequence-based predictors. This showcases the immense potential of pairing advanced representation learning and deep neural network designs for uncovering the complex syntax governing protein-DNA interactions directly from primary sequences. Our work not only provides a robust computational tool for characterizing DNA-binding mechanisms, but also highlights the transformative opportunities at the intersection of language modeling, deep learning, and protein sequence analysis. The publicly available code and data further facilitate broader adoption and continued development of these techniques for accelerating mechanistic insights into vital biological processes and disease pathways. In addition, the code and data for this work are available at https://github.com/B1607/DIRP.

Advances in the Application of Protein Language Modeling for Nucleic Acid Protein Binding Site Prediction

Protein and nucleic acid binding site prediction is a critical computational task that benefits a wide range of biological processes. Previous studies have shown that feature selection holds particular significance for this prediction task, making the generation of more discriminative features a key area of interest for many researchers. Recent progress has shown the power of protein language models in handling protein sequences, in leveraging the strengths of attention networks, and in successful applications to tasks such as protein structure prediction. This naturally raises the question of the applicability of protein language models in predicting protein and nucleic acid binding sites. Various approaches have explored this potential. This paper first describes the development of protein language models. Then, a systematic review of the latest methods for predicting protein and nucleic acid binding sites is conducted by covering benchmark sets, feature generation methods, performance comparisons, and feature ablation studies. These comparisons demonstrate the importance of protein language models for the prediction task. Finally, the paper discusses the challenges of protein and nucleic acid binding site prediction and proposes possible research directions and future trends. The purpose of this survey is to furnish researchers with actionable suggestions for comprehending the methodologies used in predicting protein-nucleic acid binding sites, fostering the creation of protein-centric language models, and tackling real-world obstacles encountered in this field.

Accurate prediction of protein function using statistics-informed graph networks

Understanding protein function is pivotal in comprehending the intricate mechanisms that underlie many crucial biological activities, with far-reaching implications in the fields of medicine, biotechnology, and drug development. However, more than 200 million proteins remain uncharacterized, and computational efforts heavily rely on protein structural information to predict annotations of varying quality. Here, we present a method that utilizes statistics-informed graph networks to predict protein functions solely from its sequence. Our method inherently characterizes evolutionary signatures, allowing for a quantitative assessment of the significance of residues that carry out specific functions. PhiGnet not only demonstrates superior performance compared to alternative approaches but also narrows the sequence-function gap, even in the absence of structural information. Our findings indicate that applying deep learning to evolutionary data can highlight functional sites at the residue level, providing valuable support for interpreting both existing properties and new functionalities of proteins in research and biomedicine.

AI-Driven Deep Learning Techniques in Protein Structure Prediction

Protein structure prediction is important for understanding their function and behavior. This review study presents a comprehensive review of the computational models used in predicting protein structure. It covers the progression from established protein modeling to state-of-the-art artificial intelligence (AI) frameworks. The paper will start with a brief introduction to protein structures, protein modeling, and AI. The section on established protein modeling will discuss homology modeling, ab initio modeling, and threading. The next section is deep learning-based models. It introduces some state-of-the-art AI models, such as AlphaFold (AlphaFold, AlphaFold2, AlphaFold3), RoseTTAFold, ProteinBERT, etc. This section also discusses how AI techniques have been integrated into established frameworks like Swiss-Model, Rosetta, and I-TASSER. The model performance is compared using the rankings of CASP14 (Critical Assessment of Structure Prediction) and CASP15. CASP16 is ongoing, and its results are not included in this review. Continuous Automated Model EvaluatiOn (CAMEO) complements the biennial CASP experiment. Template modeling score (TM-score), global distance test total score (GDT_TS), and Local Distance Difference Test (lDDT) score are discussed too. This paper then acknowledges the ongoing difficulties in predicting protein structure and emphasizes the necessity of additional searches like dynamic protein behavior, conformational changes, and protein-protein interactions. In the application section, this paper introduces some applications in various fields like drug design, industry, education, and novel protein development. In summary, this paper provides a comprehensive overview of the latest advancements in established protein modeling and deep learning-based models for protein structure predictions. It emphasizes the significant advancements achieved by AI and identifies potential areas for further investigation.

MicroRNA-Dependent Mechanisms Underlying the Function of a β-Amino Carbonyl Compound in Glioblastoma Cells

Glioblastoma (GB) is an aggressive brain malignancy characterized by its invasive nature. Current treatment has limited effectiveness, resulting in poor patients' prognoses. β-Amino carbonyl (β-AC) compounds have gained attention due to their potential anticancerous properties. In vitro assays were performed to evaluate the effects of an in-house synthesized β-AC compound, named SHG-8, upon GB cells. Small RNA sequencing (sRNA-seq) and biocomputational analyses investigated the effects of SHG-8 upon the miRNome and its bioavailability within the human body. SHG-8 exhibited significant cytotoxicity and inhibition of cell migration and proliferation in U87MG and U251MG GB cells. GB cells treated with the compound released significant amounts of reactive oxygen species (ROS). Annexin V and acridine orange/ethidium bromide staining also demonstrated that the compound led to apoptosis. sRNA-seq revealed a shift in microRNA (miRNA) expression profiles upon SHG-8 treatment and significant upregulation of miR-3648 and downregulation of miR-7973. Real-time polymerase chain reaction (RT-qPCR) demonstrated a significant downregulation of CORO1C, an oncogene and a player in the Wnt/β-catenin pathway. In silico analysis indicated SHG-8's potential to cross the blood-brain barrier. We concluded that SHG-8's inhibitory effects on GB cells may involve the deregulation of various miRNAs and the inhibition of CORO1C.

PUResNetV2.0: a deep learning model leveraging sparse representation for improved ligand binding site prediction

Accurate ligand binding site prediction (LBSP) within proteins is essential for drug discovery. We developed ProteinUNetResNetV2.0 (PUResNetV2.0), leveraging sparse representation of protein structures to improve LBSP accuracy. Our training dataset included protein complexes from 4729 protein families. Evaluations on benchmark datasets showed that PUResNetV2.0 achieved an 85.4% Distance Center Atom (DCA) success rate and a 74.7% F1 Score on the Holo801 dataset, outperforming existing methods. However, its performance in specific cases, such as RNA, DNA, peptide-like ligand, and ion binding site prediction, was limited due to constraints in our training data. Our findings underscore the potential of sparse representation in LBSP, especially for oligomeric structures, suggesting PUResNetV2.0 as a promising tool for computational drug discovery.

SOFB is a comprehensive ensemble deep learning approach for elucidating and characterizing protein-nucleic-acid-binding residues

Proteins and nucleic-acids are essential components of living organisms that interact in critical cellular processes. Accurate prediction of nucleic acid-binding residues in proteins can contribute to a better understanding of protein function. However, the discrepancy between protein sequence information and obtained structural and functional data renders most current computational models ineffective. Therefore, it is vital to design computational models based on protein sequence information to identify nucleic acid binding sites in proteins. Here, we implement an ensemble deep learning model-based nucleic-acid-binding residues on proteins identification method, called SOFB, which characterizes protein sequences by learning the semantics of biological dynamics contexts, and then develop an ensemble deep learning-based sequence network to learn feature representation and classification by explicitly modeling dynamic semantic information. Among them, the language learning model, which is constructed from natural language to biological language, captures the underlying relationships of protein sequences, and the ensemble deep learning-based sequence network consisting of different convolutional layers together with Bi-LSTM refines various features for optimal performance. Meanwhile, to address the imbalanced issue, we adopt ensemble learning to train multiple models and then incorporate them. Our experimental results on several DNA/RNA nucleic-acid-binding residue datasets demonstrate that our proposed model outperforms other state-of-the-art methods. In addition, we conduct an interpretability analysis of the identified nucleic acid binding residue sequences based on the attention weights of the language learning model, revealing novel insights into the dynamic semantic information that supports the identified nucleic acid binding residues. SOFB is available at https://github.com/Encryptional/SOFB and https://figshare.com/articles/online_resource/SOFB_figshare_rar/25499452 .

A comprehensive survey on protein-ligand binding site prediction

Protein-ligand binding site prediction is critical for protein function annotation and drug discovery. Biological experiments are time-consuming and require significant equipment, materials, and labor resources. Developing accurate and efficient computational methods for protein-ligand interaction prediction is essential. Here, we summarize the key challenges associated with ligand binding site (LBS) prediction and introduce recently published methods from their input features, computational algorithms, and ligand types. Furthermore, we investigate the specificity of allosteric site identification as a particular LBS type. Finally, we discuss the prospective directions for machine learning-based LBS prediction in the near future.

Revolutionizing GPCR-ligand predictions: DeepGPCR with experimental validation for high-precision drug discovery

G-protein coupled receptors (GPCRs), crucial in various diseases, are targeted of over 40% of approved drugs. However, the reliable acquisition of experimental GPCRs structures is hindered by their lipid-embedded conformations. Traditional protein-ligand interaction models falter in GPCR-drug interactions, caused by limited and low-quality structures. Generalized models, trained on soluble protein-ligand pairs, are also inadequate. To address these issues, we developed two models, DeepGPCR_BC for binary classification and DeepGPCR_RG for affinity prediction. These models use non-structural GPCR-ligand interaction data, leveraging graph convolutional networks and mol2vec techniques to represent binding pockets and ligands as graphs. This approach significantly speeds up predictions while preserving critical physical-chemical and spatial information. In independent tests, DeepGPCR_BC surpassed Autodock Vina and Schrödinger Dock with an area under the curve of 0.72, accuracy of 0.68 and true positive rate of 0.73, whereas DeepGPCR_RG demonstrated a Pearson correlation of 0.39 and root mean squared error of 1.34. We applied these models to screen drug candidates for GPR35 (Q9HC97), yielding promising results with three (F545-1970, K297-0698, S948-0241) out of eight candidates. Furthermore, we also successfully obtained six active inhibitors for GLP-1R. Our GPCR-specific models pave the way for efficient and accurate large-scale virtual screening, potentially revolutionizing drug discovery in the GPCR field.

Q-BioLiP: A Comprehensive Resource for Quaternary Structure-based Protein-ligand Interactions

Since its establishment in 2013, BioLiP has become one of the widely used resources for protein-ligand interactions. Nevertheless, several known issues occurred with it over the past decade. For example, the protein-ligand interactions are represented in the form of single chain-based tertiary structures, which may be inappropriate as many interactions involve multiple protein chains (known as quaternary structures). We sought to address these issues, resulting in Q-BioLiP, a comprehensive resource for quaternary structure-based protein-ligand interactions. The major features of Q-BioLiP include: (1) representing protein structures in the form of quaternary structures rather than single chain-based tertiary structures; (2) pairing DNA/RNA chains properly rather than separation; (3) providing both experimental and predicted binding affinities; (4) retaining both biologically relevant and irrelevant interactions to alleviate the wrong justification of ligands' biological relevance; and (5) developing a new quaternary structure-based algorithm for the modelling of protein-ligand complex structure. With these new features, Q-BioLiP is expected to be a valuable resource for studying biomolecule interactions, including protein-small molecule interaction, protein-metal ion interaction, protein-peptide interaction, protein-protein interaction, protein-DNA/RNA interaction, and RNA-small molecule interaction. Q-BioLiP is freely available at https://yanglab.qd.sdu.edu.cn/Q-BioLiP/.

Leveraging machine learning models for peptide-protein interaction prediction

Peptides play a pivotal role in a wide range of biological activities through participating in up to 40% protein-protein interactions in cellular processes. They also demonstrate remarkable specificity and efficacy, making them promising candidates for drug development. However, predicting peptide-protein complexes by traditional computational approaches, such as docking and molecular dynamics simulations, still remains a challenge due to high computational cost, flexible nature of peptides, and limited structural information of peptide-protein complexes. In recent years, the surge of available biological data has given rise to the development of an increasing number of machine learning models for predicting peptide-protein interactions. These models offer efficient solutions to address the challenges associated with traditional computational approaches. Furthermore, they offer enhanced accuracy, robustness, and interpretability in their predictive outcomes. This review presents a comprehensive overview of machine learning and deep learning models that have emerged in recent years for the prediction of peptide-protein interactions.

Computational exploration of novel antimicrobial modalities targeting fucose-binding lectins and ribosomes in Mycobacterium smegmatis using tRNA-encoded peptides

tRNA-Encoded Peptides (tREPs), encoded by small open reading frames (smORFs) within tRNA genes, have recently emerged as a new class of functional peptides exhibiting antiparasitic activity. The discovery of tREPs has led to a re-evaluation of the role of tRNAs in biology and has expanded our understanding of the genetic code. This presents an immense, unexplored potential in the realm of tRNA-peptide interactions, paving the way for groundbreaking discoveries and innovative applications in various biological functions. This study explores the antimicrobial potential of tREPs against protein targets by employing a computational method that uses verified data sources and highly recognized predictive algorithms to provide a sorted list of likely antimicrobial peptides, which were then filtered for toxicity, cell permeability, allergenicity and half-life. These peptides were then docked with screened protein targets and computationally validated using molecular dynamics (MD) simulations for 150 ns and the binding free energy was estimated. The peptides Pep2 (VVLWRKPRVRKTG) and Pep6 (HRLRLRRRKPWW) exhibited good binding affinities of -110.5 +/- 2.5 and -129.0 +/- 3.9, respectively, with RMSD values of 0.4 and 0.25 nm against the fucose-binding lectin (7NEF) and the 30S ribosome of Mycobacterium smegmatis (5O5J) protein targets. The 7NEF-Pep2 and 5O5J-Pep6 complexes indicated higher negative binding free energies of -52.55 kcal/mol and -55.52 kcal/mol respectively, as calculated by Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA). Thus, the tREPs derived peptides designed as a part of this study, provide novel approaches for potential anti-bacterial therapeutic modalities.Communicated by Ramaswamy H. Sarma.

Generalized biomolecular modeling and design with RoseTTAFold All-Atom

Deep-learning methods have revolutionized protein structure prediction and design but are presently limited to protein-only systems. We describe RoseTTAFold All-Atom (RFAA), which combines a residue-based representation of amino acids and DNA bases with an atomic representation of all other groups to model assemblies that contain proteins, nucleic acids, small molecules, metals, and covalent modifications, given their sequences and chemical structures. By fine-tuning on denoising tasks, we developed RFdiffusion All-Atom (RFdiffusionAA), which builds protein structures around small molecules. Starting from random distributions of amino acid residues surrounding target small molecules, we designed and experimentally validated, through crystallography and binding measurements, proteins that bind the cardiac disease therapeutic digoxigenin, the enzymatic cofactor heme, and the light-harvesting molecule bilin.

A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites

Protein evolution is constrained by structure and function, creating patterns in residue conservation that are routinely exploited to predict structure and other features. Similar constraints should affect variation across individuals, but it is only with the growth of human population sequencing that this has been tested at scale. Now, human population constraint has established applications in pathogenicity prediction, but it has not yet been explored for structural inference. Here, we map 2.4 million population variants to 5885 protein families and quantify residue-level constraint with a new Missense Enrichment Score (MES). Analysis of 61,214 structures from the PDB spanning 3661 families shows that missense depleted sites are enriched in buried residues or those involved in small-molecule or protein binding. MES is complementary to evolutionary conservation and a combined analysis allows a new classification of residues according to a conservation plane. This approach finds functional residues that are evolutionarily diverse, which can be related to specificity, as well as family-wide conserved sites that are critical for folding or function. We also find a possible contrast between lethal and non-lethal pathogenic sites, and a surprising clinical variant hot spot at a subset of missense enriched positions.

A High-Quality Data Set of Protein-Ligand Binding Interactions Via Comparative Complex Structure Modeling

High-quality protein-ligand complex structures provide the basis for understanding the nature of noncovalent binding interactions at the atomic level and enable structure-based drug design. However, experimentally determined complex structures are scarce compared with the vast chemical space. In this study, we addressed this issue by constructing the BindingNet data set via comparative complex structure modeling, which contains 69,816 modeled high-quality protein-ligand complex structures with experimental binding affinity data. BindingNet provides valuable insights into investigating protein-ligand interactions, allowing visual inspection and interpretation of structural analogues' structure-activity relationships. It can also be used for evaluating machine-learning-based scoring functions. Our results indicate that machine learning models trained on BindingNet could reduce the bias caused by buried solvent-accessible surface area, as we previously found for models trained on the PDBbind data set. We also discussed strategies to improve BindingNet and its potential utilization for benchmarking the molecular docking methods and ligand binding free energy calculation approaches. The BindingNet complements PDBbind in constructing a sufficient and unbiased protein-ligand binding data set and is freely available at http://bindingnet.huanglab.org.cn.

Predicting Protein Functions Based on Heterogeneous Graph Attention Technique

In bioinformatics, protein function prediction stands as a fundamental area of research and plays a crucial role in addressing various biological challenges, such as the identification of potential targets for drug discovery and the elucidation of disease mechanisms. However, known functional annotation databases usually provide positive experimental annotations that proteins carry out a given function, and rarely record negative experimental annotations that proteins do not carry out a given function. Therefore, existing computational methods based on deep learning models focus on these positive annotations for prediction and ignore these scarce but informative negative annotations, leading to an underestimation of precision. To address this issue, we introduce a deep learning method that utilizes a heterogeneous graph attention technique. The method first constructs a heterogeneous graph that covers the protein-protein interaction network, ontology structure, and positive and negative annotation information. Then, it learns embedding representations of proteins and ontology terms by using the heterogeneous graph attention technique. Finally, it leverages these learned representations to reconstruct the positive protein-term associations and score unobserved functional annotations. It can enhance the predictive performance by incorporating these known limited negative annotations into the constructed heterogeneous graph. Experimental results on three species (i.e., Human, Mouse, and Arabidopsis) demonstrate that our method can achieve better performance in predicting new protein annotations than state-of-the-art methods.

Ab initio modeling of human IRS1 protein to find novel target to dock with drug MH to mitigate T2DM diabetes by insulin signaling

IRS1 is a cytoplasmic adaptor protein that helps in cellular growth, glucose metabolism, proliferation, and differentiation. Highly disordered (insulin receptor substrate 1) IRS1 protein sequence (mol.wt- 131,590.97 da) has been used to develop model using ab initio modeling technique by I-Tassar tool and Discovery Studio/ DogSite Server to decipher a novel active site. The constructed protein model has been submitted with PMDB Id- PM0082210. GRAVY index of IRS1 model ( - 0.675) indicated surface protein-water interaction. Protparam tool instability index (75.22) demonstrated disorderedness combined with loops owing to prolines/glycines. After refinement, the Ramachandran plot showed that 88 percent of AAs were present in the allowed region and only 0.5% in the disallowed region. Novel IRS1 model protein has 10 α-helices, 22 β-sheets, 20 β-hairpins, 5 β-bulges, 47 strands, 105 β-turns, and 8 γ-turns. Docking of IRS1 with drug MH demonstrated interaction of Ser-70, Thr-18, and Pro-69 with C-H bonds; Gln-71, and Glu-113 with hydrogen bonds; while both Glu-114 and Glu-113 with salt-bridge connection. Permissible 1.0-1.5 Å range of RMSD fluctuation between 20 and 45 ns was obtained in simulation of IRS1 and IRS1-met complex confirmed that both complexes were stable during whole simulation process. RMSF result showed that except positions 57AA and 114AA, the binding of drug had no severe effects on the flexibility of the IRS1 and IRS1-met complex. The RoG value of compactness and rigidity showed little change in IRS1 protein. SASA value of IRS1 indicated non-significant fluctuation between IRS1 and drug MH means ligand (drug) and IRS1 receptor form stable structure. Hydrogen bond strength of IRS1 and IRS1-met was 81.2 and 76.4, respectively, which suggested stable interaction.

EquiPNAS: improved protein-nucleic acid binding site prediction using protein-language-model-informed equivariant deep graph neural networks

Protein language models (pLMs) trained on a large corpus of protein sequences have shown unprecedented scalability and broad generalizability in a wide range of predictive modeling tasks, but their power has not yet been harnessed for predicting protein-nucleic acid binding sites, critical for characterizing the interactions between proteins and nucleic acids. Here, we present EquiPNAS, a new pLM-informed E(3) equivariant deep graph neural network framework for improved protein-nucleic acid binding site prediction. By combining the strengths of pLM and symmetry-aware deep graph learning, EquiPNAS consistently outperforms the state-of-the-art methods for both protein-DNA and protein-RNA binding site prediction on multiple datasets across a diverse set of predictive modeling scenarios ranging from using experimental input to AlphaFold2 predictions. Our ablation study reveals that the pLM embeddings used in EquiPNAS are sufficiently powerful to dramatically reduce the dependence on the availability of evolutionary information without compromising on accuracy, and that the symmetry-aware nature of the E(3) equivariant graph-based neural architecture offers remarkable robustness and performance resilience. EquiPNAS is freely available at https://github.com/Bhattacharya-Lab/EquiPNAS.

Novel endolithic bacteria of phylum Chloroflexota reveal a myriad of potential survival strategies in the Antarctic desert

The ice-free McMurdo Dry Valleys of Antarctica are dominated by nutrient-poor mineral soil and rocky outcrops. The principal habitat for microorganisms is within rocks (endolithic). In this environment, microorganisms are provided with protection against sub-zero temperatures, rapid thermal fluctuations, extreme dryness, and ultraviolet and solar radiation. Endolithic communities include lichen, algae, fungi, and a diverse array of bacteria. Chloroflexota is among the most abundant bacterial phyla present in these communities. Among the Chloroflexota are four novel classes of bacteria, here named Candidatus Spiritibacteria class. nov. (=UBA5177), Candidatus Martimicrobia class. nov. (=UBA4733), Candidatus Tarhunnaeia class. nov. (=UBA6077), and Candidatus Uliximicrobia class. nov. (=UBA2235). We retrieved 17 high-quality metagenome-assembled genomes (MAGs) that represent these four classes. Based on genome predictions, all these bacteria are inferred to be aerobic heterotrophs that encode enzymes for the catabolism of diverse sugars. These and other organic substrates are likely derived from lichen, algae, and fungi, as metabolites (including photosynthate), cell wall components, and extracellular matrix components. The majority of MAGs encode the capacity for trace gas oxidation using high-affinity uptake hydrogenases, which could provide energy and metabolic water required for survival and persistence. Furthermore, some MAGs encode the capacity to couple the energy generated from H2 and CO oxidation to support carbon fixation (atmospheric chemosynthesis). All encode mechanisms for the detoxification and efflux of heavy metals. Certain MAGs encode features that indicate possible interactions with other organisms, such as Tc-type toxin complexes, hemolysins, and macroglobulins.IMPORTANCEThe ice-free McMurdo Dry Valleys of Antarctica are the coldest and most hyperarid desert on Earth. It is, therefore, the closest analog to the surface of the planet Mars. Bacteria and other microorganisms survive by inhabiting airspaces within rocks (endolithic). We identify four novel classes of phylum Chloroflexota, and, based on interrogation of 17 metagenome-assembled genomes, we predict specific metabolic and physiological adaptations that facilitate the survival of these bacteria in this harsh environment-including oxidation of trace gases and the utilization of nutrients (including sugars) derived from lichen, algae, and fungi. We propose that such adaptations allow these endolithic bacteria to eke out an existence in this cold and extremely dry habitat.

Machine Learning for Sequence and Structure-Based Protein-Ligand Interaction Prediction

Developing new drugs is too expensive and time -consuming. Accurately predicting the interaction between drugs and targets will likely change how the drug is discovered. Machine learning-based protein-ligand interaction prediction has demonstrated significant potential. In this paper, computational methods, focusing on sequence and structure to study protein-ligand interactions, are examined. Therefore, this paper starts by presenting an overview of the data sets applied in this area, as well as the various approaches applied for representing proteins and ligands. Then, sequence-based and structure-based classification criteria are subsequently utilized to categorize and summarize both the classical machine learning models and deep learning models employed in protein-ligand interaction studies. Moreover, the evaluation methods and interpretability of these models are proposed. Furthermore, delving into the diverse applications of protein-ligand interaction models in drug research is presented. Lastly, the current challenges and future directions in this field are addressed.

Leveraging Machine Learning Models for Peptide-Protein Interaction Prediction

Peptides play a pivotal role in a wide range of biological activities through participating in up to 40% protein-protein interactions in cellular processes. They also demonstrate remarkable specificity and efficacy, making them promising candidates for drug development. However, predicting peptide-protein complexes by traditional computational approaches, such as Docking and Molecular Dynamics simulations, still remains a challenge due to high computational cost, flexible nature of peptides, and limited structural information of peptide-protein complexes. In recent years, the surge of available biological data has given rise to the development of an increasing number of machine learning models for predicting peptide-protein interactions. These models offer efficient solutions to address the challenges associated with traditional computational approaches. Furthermore, they offer enhanced accuracy, robustness, and interpretability in their predictive outcomes. This review presents a comprehensive overview of machine learning and deep learning models that have emerged in recent years for the prediction of peptide-protein interactions.

Expression and Purification of Cp3GT: Structural Analysis and Modeling of a Key Plant Flavonol-3-O Glucosyltransferase from Citrus paradisi

Glycosyltransferases (GTs) are pivotal enzymes in the biosynthesis of various biological molecules. This study focuses on the scale-up, expression, and purification of a plant flavonol-specific 3-O glucosyltransferase (Cp3GT), a key enzyme from Citrus paradisi, for structural analysis and modeling. The challenges associated with recombinant protein production in Pichia pastoris, such as proteolytic degradation, were addressed through the optimization of culture conditions and purification processes. The purification strategy employed affinity, anion exchange, and size exclusion chromatography, leading to greater than 95% homogeneity for Cp3GT. In silico modeling, using D-I-TASSER and COFACTOR integrated with the AlphaFold2 pipeline, provided insights into the structural dynamics of Cp3GT and its ligand binding sites, offering predictions for enzyme-substrate interactions. These models were compared to experimentally derived structures, enhancing understanding of the enzyme's functional mechanisms. The findings present a comprehensive approach to produce a highly purified Cp3GT which is suitable for crystallographic studies and to shed light on the structural basis of flavonol specificity in plant GTs. The significant implications of these results for synthetic biology and enzyme engineering in pharmaceutical applications are also considered.

The evolution and functional divergence of 10 Apolipoprotein D-like genes in Nilaparvata lugens

Apolipoprotein D (ApoD), a member of the lipocalin superfamily of proteins, is involved in lipid transport and stress resistance. Whereas only a single copy of the ApoD gene is found in humans and some other vertebrates, there are typically several ApoD-like genes in insects. To date, there have been relatively few studies that have examined the evolution and functional differentiation of ApoD-like genes in insects, particularly hemi-metabolous insects. In this study, we identified 10 ApoD-like genes (NlApoD1-10) with distinct spatiotemporal expression patterns in Nilaparvata lugens (BPH), which is an important pest of rice. NlApoD1-10 were found to be distributed on 3 chromosomes in a tandem array of NlApoD1/2, NlApoD3-5, and NlApoD7/8, and show sequence and gene structural divergence in the coding regions, indicating that multiple gene duplication events occurred during evolution. Phylogenetic analysis revealed that NlApoD1-10 can be clustered into 5 clades, with NlApoD3-5 and NlApoD7/8 potentially evolving exclusively in the Delphacidae family. Functional screening using an RNA interference approach revealed that only NlApoD2 was essential for BPH development and survival, whereas NlApoD4/5 are highly expressed in testes, and might play roles in reproduction. Moreover, stress response analysis revealed that NlApoD3-5/9, NlApoD3-5, and NlApoD9 were up-regulated after treatment with lipopolysaccharide, H2 O2 , and ultraviolet-C, respectively, indicating their potential roles in stress resistance.

HybridDBRpred: improved sequence-based prediction of DNA-binding amino acids using annotations from structured complexes and disordered proteins

Current predictors of DNA-binding residues (DBRs) from protein sequences belong to two distinct groups, those trained on binding annotations extracted from structured protein-DNA complexes (structure-trained) vs. intrinsically disordered proteins (disorder-trained). We complete the first empirical analysis of predictive performance across the structure- and disorder-annotated proteins for a representative collection of ten predictors. Majority of the structure-trained tools perform well on the structure-annotated proteins while doing relatively poorly on the disorder-annotated proteins, and vice versa. Several methods make accurate predictions for the structure-annotated proteins or the disorder-annotated proteins, but none performs highly accurately for both annotation types. Moreover, most predictors make excessive cross-predictions for the disorder-annotated proteins, where residues that interact with non-DNA ligand types are predicted as DBRs. Motivated by these results, we design, validate and deploy an innovative meta-model, hybridDBRpred, that uses deep transformer network to combine predictions generated by three best current predictors. HybridDBRpred provides accurate predictions and low levels of cross-predictions across the two annotation types, and is statistically more accurate than each of the ten tools and baseline meta-predictors that rely on averaging and logistic regression. We deploy hybridDBRpred as a convenient web server at http://biomine.cs.vcu.edu/servers/hybridDBRpred/ and provide the corresponding source code at https://github.com/jianzhang-xynu/hybridDBRpred.

BioLiP2: an updated structure database for biologically relevant ligand-protein interactions

With the progress of structural biology, the Protein Data Bank (PDB) has witnessed rapid accumulation of experimentally solved protein structures. Since many structures are determined with purification and crystallization additives that are unrelated to a protein's in vivo function, it is nontrivial to identify the subset of protein-ligand interactions that are biologically relevant. We developed the BioLiP2 database (https://zhanggroup.org/BioLiP) to extract biologically relevant protein-ligand interactions from the PDB database. BioLiP2 assesses the functional relevance of the ligands by geometric rules and experimental literature validations. The ligand binding information is further enriched with other function annotations, including Enzyme Commission numbers, Gene Ontology terms, catalytic sites, and binding affinities collected from other databases and a manual literature survey. Compared to its predecessor BioLiP, BioLiP2 offers significantly greater coverage of nucleic acid-protein interactions, and interactions involving large complexes that are unavailable in PDB format. BioLiP2 also integrates cutting-edge structural alignment algorithms with state-of-the-art structure prediction techniques, which for the first time enables composite protein structure and sequence-based searching and significantly enhances the usefulness of the database in structure-based function annotations. With these new developments, BioLiP2 will continue to be an important and comprehensive database for docking, virtual screening, and structure-based protein function analyses.

Predictive Modeling of Proteins Encoded by a Plant Virus Sheds a New Light on Their Structure and Inherent Multifunctionality

Plant virus genomes encode proteins that are involved in replication, encapsidation, cell-to-cell, and long-distance movement, avoidance of host detection, counter-defense, and transmission from host to host, among other functions. Even though the multifunctionality of plant viral proteins is well documented, contemporary functional repertoires of individual proteins are incomplete. However, these can be enhanced by modeling tools. Here, predictive modeling of proteins encoded by the two genomic RNAs, i.e., RNA1 and RNA2, of grapevine fanleaf virus (GFLV) and their satellite RNAs by a suite of protein prediction software confirmed not only previously validated functions (suppressor of RNA silencing [VSR], viral genome-linked protein [VPg], protease [Pro], symptom determinant [Sd], homing protein [HP], movement protein [MP], coat protein [CP], and transmission determinant [Td]) and previously identified putative functions (helicase [Hel] and RNA-dependent RNA polymerase [Pol]), but also predicted novel functions with varying levels of confidence. These include a T3/T7-like RNA polymerase domain for protein 1AVSR, a short-chain reductase for protein 1BHel/VSR, a parathyroid hormone family domain for protein 1EPol/Sd, overlapping domains of unknown function and an ABC transporter domain for protein 2BMP, and DNA topoisomerase domains, transcription factor FBXO25 domain, or DNA Pol subunit cdc27 domain for the satellite RNA protein. Structural predictions for proteins 2AHP/Sd, 2BMP, and 3A? had low confidence, while predictions for proteins 1AVSR, 1BHel*/VSR, 1CVPg, 1DPro, 1EPol*/Sd, and 2CCP/Td retained higher confidence in at least one prediction. This research provided new insights into the structure and functions of GFLV proteins and their satellite protein. Future work is needed to validate these findings.

M-Ionic: prediction of metal-ion-binding sites from sequence using residue embeddings

MOTIVATION

Understanding metal-protein interaction can provide structural and functional insights into cellular processes. As the number of protein sequences increases, developing fast yet precise computational approaches to predict and annotate metal-binding sites becomes imperative. Quick and resource-efficient pre-trained protein language model (pLM) embeddings have successfully predicted binding sites from protein sequences despite not using structural or evolutionary features (multiple sequence alignments). Using residue-level embeddings from the pLMs, we have developed a sequence-based method (M-Ionic) to identify metal-binding proteins and predict residues involved in metal binding.

RESULTS

On independent validation of recent proteins, M-Ionic reports an area under the curve (AUROC) of 0.83 (recall = 84.6%) in distinguishing metal binding from non-binding proteins compared to AUROC of 0.74 (recall = 61.8%) of the next best method. In addition to comparable performance to the state-of-the-art method for identifying metal-binding residues (Ca2+, Mg2+, Mn2+, Zn2+), M-Ionic provides binding probabilities for six additional ions (i.e. Cu2+, Po43-, So42-, Fe2+, Fe3+, Co2+). We show that the pLM embedding of a single residue contains sufficient information about its neighbours to predict its binding properties.

AVAILABILITY AND IMPLEMENTATION

M-Ionic can be used on your protein of interest using a Google Colab Notebook (https://bit.ly/40FrRbK). The GitHub repository (https://github.com/TeamSundar/m-ionic) contains all code and data.

A Novel COCH p.D544Vfs*3 Variant Associated with DFNA9 Sensorineural Hearing Loss Causes Pathological Multimeric Cochlin Formation

COCH (coagulation factor C homology) is one of the most frequently mutated genes of autosomal dominant non-syndromic hearing loss. Variants in COCH could cause DFNA9, which is characterized by late-onset hearing loss with variable degrees of vestibular dysfunction. In this study, we report a Chinese family with a novel COCH variant (c.1687delA) causing p.D544Vfs*3 in the cochlin. Comprehensive audiometric tests and vestibular function assessments were taken to acquire the phenotypic profile of the subjects. Next-generation sequencing was conducted and segregation analysis was carried out using Sanger sequencing. The proband presented mild vestibular symptoms and normal functional assessment results in almost every test, while the variant co-segregated with hearing impairment in the pedigree. The variant was located beyond the vWFA2 domain, which was predicted to affect the post-translational cleavage of the cochlin via molecular modeling analysis. Notably, in the overexpressing study, by transient transfecting the HEK 293T cells, we found that the p.D544Vfs*3 variant increased the formation of multimeric cochlin. Our result enriched the spectrum of DFNA9-linked pathological COCH variants and suggested that variants, causative of cochlin multimerization, could be related to DFNA9 with sensorineural hearing loss rather than serious vestibular symptoms.

Rieske Oxygenase-Catalyzed Oxidative Late-Stage Functionalization during Complex Antifungal Polyketide Biosynthesis

Rieske oxygenases (ROs) from natural product biosynthetic pathways are a poorly studied group of enzymes with significant potential as oxidative functionalization biocatalysts. A study on the ROs JerL, JerP, and AmbP from the biosynthetic pathways of jerangolid A and ambruticin VS-3 is described. Their activity was successfully reconstituted using whole-cell bioconversion systems coexpressing the ROs and their respective natural flavin-dependent reductase (FDR) partners. Feeding authentic biosynthetic intermediates and synthetic surrogates to these strains confirmed the involvement of the ROs in hydroxymethylpyrone and dihydropyran formation and revealed crucial information about the RO's substrate specificity. The pronounced dependence of JerL and JerP on the presence of a methylenolether allowed the precise temporal assignment of RO catalysis to the ultimate steps of jerangolid biosynthesis. JerP and AmbP stand out among the biosynthetic ROs studied so far for their ability to catalyze clean tetrahydropyran desaturation without further functionalizing the formed electron-rich double bonds. This work highlights the remarkable ability of ROs to highly selectively oxidize complex molecular scaffolds.

Protein-protein interaction site prediction by model ensembling with hybrid feature and self-attention

BACKGROUND

Protein-protein interactions (PPIs) are crucial in various biological functions and cellular processes. Thus, many computational approaches have been proposed to predict PPI sites. Although significant progress has been made, these methods still have limitations in encoding the characteristics of each amino acid in sequences. Many feature extraction methods rely on the sliding window technique, which simply merges all the features of residues into a vector. The importance of some key residues may be weakened in the feature vector, leading to poor performance.

RESULTS

We propose a novel sequence-based method for PPI sites prediction. The new network model, PPINet, contains multiple feature processing paths. For a residue, the PPINet extracts the features of the targeted residue and its context separately. These two types of features are processed by two paths in the network and combined to form a protein representation, where the two types of features are of relatively equal importance. The model ensembling technique is applied to make use of more features. The base models are trained with different features and then ensembled via stacking. In addition, a data balancing strategy is presented, by which our model can get significant improvement on highly unbalanced data.

CONCLUSION

The proposed method is evaluated on a fused dataset constructed from Dset186, Dset_72, and PDBset_164, as well as the public Dset_448 dataset. Compared with current state-of-the-art methods, the performance of our method is better than the others. In the most important metrics, such as AUPRC and recall, it surpasses the second-best programmer on the latter dataset by 6.9% and 4.7%, respectively. We also demonstrated that the improvement is essentially due to using the ensemble model, especially, the hybrid feature. We share our code for reproducibility and future research at https://github.com/CandiceCong/StackingPPINet .

Al-Dosari M, Sinha N, Parvez MK, Sehgal D. Methotrexate, an anti-inflammatory drug, inhibits Hepatitis E viral replication

Hepatitis E Virus (HEV) is a positively oriented RNA virus having a 7.2 kb genome. HEV consists of three open reading frames (ORF1-3). Of these, ORF1 codes for the enzymes Methyltransferase (Mtase), Papain-like cysteine protease (PCP), RNA helicase, and RNA-dependent RNA polymerase (RdRp). Unavailability of a vaccine or effective drug against HEV and considering the side effects associated with the off-label use of ribavirin (RBV) and pegylated interferons, an alternative approach is required by the modulation of specific enzymes to prevent the infection. HEV helicase is involved in unwinding the double-stranded RNA, RNA processing, transcriptional regulation, and pre-mRNA processing. Therefore, we screened FDA-approved compounds from the ZINC15 database against the modelled 3D structure of HEV helicase and found that methotrexate and compound A (Pubchem ID BTB07890) inhibit the NTPase and dsRNA unwinding activity leading to inhibition of HEV RNA replication. This may be further authenticated by in vivo study.

DeepProSite: structure-aware protein binding site prediction using ESMFold and pretrained language model

MOTIVATION

Identifying the functional sites of a protein, such as the binding sites of proteins, peptides, or other biological components, is crucial for understanding related biological processes and drug design. However, existing sequence-based methods have limited predictive accuracy, as they only consider sequence-adjacent contextual features and lack structural information.

RESULTS

In this study, DeepProSite is presented as a new framework for identifying protein binding site that utilizes protein structure and sequence information. DeepProSite first generates protein structures from ESMFold and sequence representations from pretrained language models. It then uses Graph Transformer and formulates binding site predictions as graph node classifications. In predicting protein-protein/peptide binding sites, DeepProSite outperforms state-of-the-art sequence- and structure-based methods on most metrics. Moreover, DeepProSite maintains its performance when predicting unbound structures, in contrast to competing structure-based prediction methods. DeepProSite is also extended to the prediction of binding sites for nucleic acids and other ligands, verifying its generalization capability. Finally, an online server for predicting multiple types of residue is established as the implementation of the proposed DeepProSite.

AVAILABILITY AND IMPLEMENTATION

The datasets and source codes can be accessed at https://github.com/WeiLab-Biology/DeepProSite. The proposed DeepProSite can be accessed at https://inner.wei-group.net/DeepProSite/.

Automated benchmarking of combined protein structure and ligand conformation prediction

The prediction of protein-ligand complexes (PLC), using both experimental and predicted structures, is an active and important area of research, underscored by the inclusion of the Protein-Ligand Interaction category in the latest round of the Critical Assessment of Protein Structure Prediction experiment CASP15. The prediction task in CASP15 consisted of predicting both the three-dimensional structure of the receptor protein as well as the position and conformation of the ligand. This paper addresses the challenges and proposed solutions for devising automated benchmarking techniques for PLC prediction. The reliability of experimentally solved PLC as ground truth reference structures is assessed using various validation criteria. Similarity of PLC to previously released complexes are employed to judge PLC diversity and the difficulty of a PLC as a prediction target. We show that the commonly used PDBBind time-split test-set is inappropriate for comprehensive PLC evaluation, with state-of-the-art tools showing conflicting results on a more representative and high quality dataset constructed for benchmarking purposes. We also show that redocking on crystal structures is a much simpler task than docking into predicted protein models, demonstrated by the two PLC-prediction-specific scoring metrics created. Finally, we introduce a fully automated pipeline that predicts PLC and evaluates the accuracy of the protein structure, ligand pose, and protein-ligand interactions.

Assessing the fast prediction of peptide conformers and the impact of non-natural modifications

We present an assessment of different approaches to predict peptide structures using modeling tools. Several small molecule, protein, and peptide-focused methodologies were used for the fast prediction of conformers for peptides shorter than 30 amino acids. We assessed the effect of including restraints based on annotated or predicted secondary structure motifs. A number of peptides in bound conformations and in solution were collected to compare the tools. In addition, we studied the impact of changing single amino acids to non-natural residues using molecular dynamics simulations. Deep learning methods such as AlphaFold2, or the combination of physics-based approaches with secondary structure information, produce the most accurate results for natural sequences. In the case of peptides with non-natural modifications, modeling the peptide containing natural amino acids first and then modifying and simulating the peptide using benchmarked force fields is a recommended pipeline. The results can guide the modeling of oligopeptides for drug discovery projects.

Enzyme function and evolution through the lens of bioinformatics

Enzymes have been shaped by evolution over billions of years to catalyse the chemical reactions that support life on earth. Dispersed in the literature, or organised in online databases, knowledge about enzymes can be structured in distinct dimensions, either related to their quality as biological macromolecules, such as their sequence and structure, or related to their chemical functions, such as the catalytic site, kinetics, mechanism, and overall reaction. The evolution of enzymes can only be understood when each of these dimensions is considered. In addition, many of the properties of enzymes only make sense in the light of evolution. We start this review by outlining the main paradigms of enzyme evolution, including gene duplication and divergence, convergent evolution, and evolution by recombination of domains. In the second part, we overview the current collective knowledge about enzymes, as organised by different types of data and collected in several databases. We also highlight some increasingly powerful computational tools that can be used to close gaps in understanding, in particular for types of data that require laborious experimental protocols. We believe that recent advances in protein structure prediction will be a powerful catalyst for the prediction of binding, mechanism, and ultimately, chemical reactions. A comprehensive mapping of enzyme function and evolution may be attainable in the near future.

ZeroBind: a protein-specific zero-shot predictor with subgraph matching for drug-target interactions

Existing drug-target interaction (DTI) prediction methods generally fail to generalize well to novel (unseen) proteins and drugs. In this study, we propose a protein-specific meta-learning framework ZeroBind with subgraph matching for predicting protein-drug interactions from their structures. During the meta-training process, ZeroBind formulates training a protein-specific model, which is also considered a learning task, and each task uses graph neural networks (GNNs) to learn the protein graph embedding and the molecular graph embedding. Inspired by the fact that molecules bind to a binding pocket in proteins instead of the whole protein, ZeroBind introduces a weakly supervised subgraph information bottleneck (SIB) module to recognize the maximally informative and compressive subgraphs in protein graphs as potential binding pockets. In addition, ZeroBind trains the models of individual proteins as multiple tasks, whose importance is automatically learned with a task adaptive self-attention module to make final predictions. The results show that ZeroBind achieves superior performance on DTI prediction over existing methods, especially for those unseen proteins and drugs, and performs well after fine-tuning for those proteins or drugs with a few known binding partners.

PepCNN deep learning tool for predicting peptide binding residues in proteins using sequence, structural, and language model features

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

CoDock-Ligand: combined template-based docking and CNN-based scoring in ligand binding prediction

For ligand binding prediction, it is crucial for molecular docking programs to integrate template-based modeling with a precise scoring function. Here, we proposed the CoDock-Ligand docking method that combines template-based modeling and the GNINA scoring function, a Convolutional Neural Network-based scoring function, for the ligand binding prediction in CASP15. Among the 21 targets, we obtained successful predictions in top 5 submissions for 14 targets and partially successful predictions for 4 targets. In particular, for the most complicated target, H1114, which contains 56 metal cofactors and small molecules, our docking method successfully predicted the binding of most ligands. Analysis of the failed systems showed that the predicted receptor protein presented conformational changes in the backbone and side chains of the binding site residues, which may cause large structural deviations in the ligand binding prediction. In summary, our hybrid docking scheme was efficiently adapted to the ligand binding prediction challenges in CASP15.

Protein-DNA binding sites prediction based on pre-trained protein language model and contrastive learning

Protein-DNA interaction is critical for life activities such as replication, transcription and splicing. Identifying protein-DNA binding residues is essential for modeling their interaction and downstream studies. However, developing accurate and efficient computational methods for this task remains challenging. Improvements in this area have the potential to drive novel applications in biotechnology and drug design. In this study, we propose a novel approach called Contrastive Learning And Pre-trained Encoder (CLAPE), which combines a pre-trained protein language model and the contrastive learning method to predict DNA binding residues. We trained the CLAPE-DB model on the protein-DNA binding sites dataset and evaluated the model performance and generalization ability through various experiments. The results showed that the area under ROC curve values of the CLAPE-DB model on the two benchmark datasets reached 0.871 and 0.881, respectively, indicating superior performance compared to other existing models. CLAPE-DB showed better generalization ability and was specific to DNA-binding sites. In addition, we trained CLAPE on different protein-ligand binding sites datasets, demonstrating that CLAPE is a general framework for binding sites prediction. To facilitate the scientific community, the benchmark datasets and codes are freely available at https://github.com/YAndrewL/clape.

Targeting the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) with synthetic/designer unnatural nucleoside analogs: an in silico study

CONTEXT

Since the outbreak of COVID-19 in December 2019, it developed into a pandemic affecting all the countries and millions of people around the globe. Until now, there is no medicine available to contain the spread of the virus. As an aid to drug discovery, the molecular docking and molecular dynamic tools were applied extensively. In silico studies made it possible for rapid screening of potential molecules as possible inhibitors/drugs against the targeted proteins. As a continuation of our drug discovery research, we have carried out molecular docking studies of our 12 reported unnatural nucleosides and 14 designer Avigan analogs with SARS-CoV-2, RNA-dependent RNA polymerase (RdRp), which we want to report herein. The same calculation was also carried out, taking 11 known/under trail/commercial nucleoside drug molecules for a comparison of the binding interactions in the catalytic site of RdRp. The docking results and binding efficiencies of our reported nucleosides and designer nucleosidic were compared with the binding energy of commercially available drugs such as remdesevir and favipiravir. Furthermore, we evaluated the protein-drug binding efficiency and stability of the best docked molecules by molecular dynamic studies (MD). From our study, we have found that few of our proposed drugs show promising binding efficiency at the catalytic pocket of SARS-CoV-2 RdRp and can be a promising RdRp inhibitor drug candidate. Hence, this study will be of importance to make progress toward developing successful nucleoside-based drugs and conduct the antiviral test in the wet lab to understand their efficacy against COVID-19.

METHOD

All the docking studies were carried out with AutoDock 4.2, AutoDock Vina and Molegro Virtual Docker. Following the docking studies, the MD simulations were carried out following the standard protocol with the GROMACS ver. 2019.6. by applying the CHARMM36 all-atom biomolecular force field. The drug-protein interaction was studied using the Biovia Discovery Studio suite, Ligplot software, and Protein-Ligand Interaction Profiler (PLIP).

How Ligands Interact with the Kinase Hinge

ATP-competitive kinase inhibitors form hydrogen bond interactions with the kinase hinge region at the adenine binding site. Thus, it is crucial to explore hinge-ligand recognition as part of a rational drug design strategy. Here, harnessing known ligand-bound kinase structures and experimental assay resources, we first created a kinase structure-assay database (KSAD) containing 2705 nM ligand-bound kinase complexes. Then, using KSAD, we systematically investigate hinge-ligand binding patterns using interaction fingerprints, thereby delineating 15 different hydrogen-bond interaction modes. We believe these results will be valuable for de novo drug design and/or scaffold hopping of kinase-targeted drugs.

CSM-Potential2: A comprehensive deep learning platform for the analysis of protein interacting interfaces

Proteins are molecular machinery that participate in virtually all essential biological functions within the cell, which are tightly related to their 3D structure. The importance of understanding protein structure-function relationship is highlighted by the exponential growth of experimental structures, which has been greatly expanded by recent breakthroughs in protein structure prediction, most notably RosettaFold, and AlphaFold2. These advances have prompted the development of several computational approaches that leverage these data sources to explore potential biological interactions. However, most methods are generally limited to analysis of single types of interactions, such as protein-protein or protein-ligand interactions, and their complexity limits the usability to expert users. Here we report CSM-Potential2, a deep learning platform for the analysis of binding interfaces on protein structures. In addition to prediction of protein-protein interactions binding sites and classification of biological ligands, our new platform incorporates prediction of interactions with nucleic acids at the residue level and allows for ligand transplantation based on sequence and structure similarity to experimentally determined structures. We anticipate our platform to be a valuable resource that provides easy access to a range of state-of-the-art methods to expert and non-expert users for the study of biological interactions. Our tool is freely available as an easy-to-use web server and API available at https://biosig.lab.uq.edu.au/csm_potential.