Protein structure prediction in the deep learning era

SERT-StructNet: Protein secondary structure prediction method based on multi-factor hybrid deep model

Protein secondary structure prediction (PSSP) is a pivotal research endeavour that plays a crucial role in the comprehensive elucidation of protein functions and properties. Current prediction methodologies are focused on deep-learning techniques, particularly focusing on multi-factor features. Diverging from existing approaches, in this study, we placed special emphasis on the effects of amino acid properties and protein secondary structure propensity scores (SSPs) on secondary structure during the meticulous selection of multi-factor features. This differential feature-selection strategy results in a distinctive and effective amalgamation of the sequence and property features. To harness these multi-factor features optimally, we introduced a hybrid deep feature extraction model. The model initially employs mechanisms such as dilated convolution (D-Conv) and a channel attention network (SENet) for local feature extraction and targeted channel enhancement. Subsequently, a combination of recurrent neural network variants (BiGRU and BiLSTM), along with a transformer module, was employed to achieve global bidirectional information consideration and feature enhancement. This approach to multi-factor feature input and multi-level feature processing enabled a comprehensive exploration of intricate associations among amino acid residues in protein sequences, yielding a Q 3 accuracy of 84.9% and an Sov score of 85.1%. The overall performance surpasses that of the comparable methods. This study introduces a novel and efficient method for determining the PSSP domain, which is poised to deepen our understanding of the practical applications of protein molecular structures.

Engineering the thermal stability of a polyphosphate kinase by ancestral sequence reconstruction to expand the temperature boundary for an industrially applicable ATP regeneration system

ATP-dependent energy-consuming enzymatic reactions are widely used in cell-free biocatalysis. However, the direct addition of large amounts of expensive ATP can greatly increase cost, and enzymatic production is often difficult to achieve as a result. Although a polyphosphate kinase (PPK)-polyphosphate-based ATP regeneration system has the potential to solve this challenge, the generally poor thermal stability of PPKs limits the widespread use of this method. In this paper, we evaluated the thermal stability of a PPK from Sulfurovum lithotrophicum (SlPPK2). After directed evolution and computation-supported design, we found that SlPPK2 is very recalcitrant and cannot acquire beneficial mutations. Inspired by the usually outstanding stability of ancestral enzymes, we reconstructed the ancestral sequence of the PPK family and used it as a guide to construct three heat-stable variants of SlPPK2, of which the L35F/T144S variant has a half-life of more than 14 h at 60°C. Molecular dynamics simulations were performed on all enzymes to analyze the reasons for the increased thermal stability. The results showed that mutations at these two positions act synergistically from the interior and surface of the protein, leading to a more compact structure. Finally, the robustness of the L35F/T144S variant was verified in the synthesis of nucleotides at high temperature. In practice, the use of this high-temperature ATP regeneration system can effectively avoid byproduct accumulation. Our work extends the temperature boundary of ATP regeneration and has great potential for industrial applications.IMPORTANCEATP regeneration is an important basic applied study in the field of cell-free biocatalysis. Polyphosphate kinase (PPK) is an enzyme tool widely used for energy regeneration during enzymatic reactions. However, the thermal stability of the PPKs reported to date that can efficiently regenerate ATP is usually poor, which greatly limits their application. In this study, the thermal stability of a difficult-to-engineer PPK from Sulfurovum lithotrophicum was improved, guided by an ancestral sequence reconstruction strategy. The optimal variant has a 4.5-fold longer half-life at 60°C than the wild-type enzyme, thus enabling the extension of the temperature boundary for ATP regeneration. The ability of this variant to regenerate ATP was well demonstrated during high-temperature enzymatic production of nucleotides.

RNA threading with secondary structure and sequence profile

MOTIVATION

RNA threading aims to identify remote homologies for template-based modeling of RNA 3D structure. Existing RNA alignment methods primarily rely on secondary structure alignment. They are often time- and memory-consuming, limiting large-scale applications. In addition, the accuracy is far from satisfactory.

RESULTS

Using RNA secondary structure and sequence profile, we developed a novel RNA threading algorithm, named RNAthreader. To enhance the alignment process and minimize memory usage, a novel approach has been introduced to simplify RNA secondary structures into compact diagrams. RNAthreader employs a two-step methodology. Initially, integer programming and dynamic programming are combined to create an initial alignment for the simplified diagram. Subsequently, the final alignment is obtained using dynamic programming, taking into account the initial alignment derived from the previous step. The benchmark test on 80 RNAs illustrates that RNAthreader generates more accurate alignments than other methods, especially for RNAs with pseudoknots. Another benchmark, involving 30 RNAs from the RNA-Puzzles experiments, exhibits that the models constructed using RNAthreader templates have a lower average RMSD than those created by alternative methods. Remarkably, RNAthreader takes less than two hours to complete alignments with ∼5000 RNAs, which is 3-40 times faster than other methods. These compelling results suggest that RNAthreader is a promising algorithm for RNA template detection.

AVAILABILITY AND IMPLEMENTATION

https://yanglab.qd.sdu.edu.cn/RNAthreader.

Decrypting orphan GPCR drug discovery via multitask learning

The drug discovery of G protein-coupled receptors (GPCRs) superfamily using computational models is often limited by the availability of protein three-dimensional (3D) structures and chemicals with experimentally measured bioactivities. Orphan GPCRs without known ligands further complicate the process. To enable drug discovery for human orphan GPCRs, multitask models were proposed for predicting half maximal effective concentrations (EC50) of the pairs of chemicals and GPCRs. Protein multiple sequence alignment features, and physicochemical properties and fingerprints of chemicals were utilized to encode the protein and chemical information, respectively. The protein features enabled the transfer of data-rich GPCRs to orphan receptors and the transferability based on the similarity of protein features. The final model was trained using both agonist and antagonist data from 200 GPCRs and showed an excellent mean squared error (MSE) of 0.24 in the validation dataset. An independent test using the orphan dataset consisting of 16 receptors associated with less than 8 bioactivities showed a reasonably good MSE of 1.51 that can be further improved to 0.53 by considering the transferability based on protein features. The informative features were identified and mapped to corresponding 3D structures to gain insights into the mechanism of GPCR-ligand interactions across the GPCR family. The proposed method provides a novel perspective on learning ligand bioactivity within the diverse human GPCR superfamily and can potentially accelerate the discovery of therapeutic agents for orphan GPCRs.

Improved protein structure prediction with trRosettaX2, AlphaFold2, and optimized MSAs in CASP15

We present the monomer and multimer structure prediction results of our methods in CASP15. We first designed an elaborate pipeline that leverages complementary sequence databases and advanced database searching algorithms to generate high-quality multiple sequence alignments (MSAs). Top MSAs were then selected for the subsequent step of structure prediction. We utilized trRosettaX2 and AlphaFold2 for monomer structure prediction (group name Yang-Server), and AlphaFold-Multimer for multimer structure prediction (group name Yang-Multimer). Yang-Server and Yang-Multimer are ranked at the top and the fourth, respectively, for monomer and multimer structure prediction. For 94 monomers, the average TM-score of the predicted structure models by Yang-Server is 0.876, compared to 0.798 by the default AlphaFold2 (i.e., the group NBIS-AF2-standard). For 42 multimers, the average DockQ score of the predicted structure models by Yang-Multimer is 0.464, compared to 0.389 by the default AlphaFold-Multimer (i.e., the group NBIS-AF2-multimer). Detailed analysis of the results shows that several factors contribute to the improvement, including improved MSAs, iterated modeling for large targets, interplay between monomer and multimer structure prediction for intertwined structures, etc. However, the structure predictions for orphan proteins and multimers remain challenging, and breakthroughs in this area are anticipated in the future.

Deep transfer learning for inter-chain contact predictions of transmembrane protein complexes

Membrane proteins are encoded by approximately a quarter of human genes. Inter-chain residue-residue contact information is important for structure prediction of membrane protein complexes and valuable for understanding their molecular mechanism. Although many deep learning methods have been proposed to predict the intra-protein contacts or helix-helix interactions in membrane proteins, it is still challenging to accurately predict their inter-chain contacts due to the limited number of transmembrane proteins. Addressing the challenge, here we develop a deep transfer learning method for predicting inter-chain contacts of transmembrane protein complexes, named DeepTMP, by taking advantage of the knowledge pre-trained from a large data set of non-transmembrane proteins. DeepTMP utilizes a geometric triangle-aware module to capture the correct inter-chain interaction from the coevolution information generated by protein language models. DeepTMP is extensively evaluated on a test set of 52 self-associated transmembrane protein complexes, and compared with state-of-the-art methods including DeepHomo2.0, CDPred, GLINTER, DeepHomo, and DNCON2_Inter. It is shown that DeepTMP considerably improves the precision of inter-chain contact prediction and outperforms the existing approaches in both accuracy and robustness.

Modeling conformational states of proteins with AlphaFold

Many proteins exert their function by switching among different structures. Knowing the conformational ensembles affiliated with these states is critical to elucidate key mechanistic aspects that govern protein function. While experimental determination efforts are still bottlenecked by cost, time, and technical challenges, the machine-learning technology AlphaFold showed near experimental accuracy in predicting the three-dimensional structure of monomeric proteins. However, an AlphaFold ensemble of models usually represents a single conformational state with minimal structural heterogeneity. Consequently, several pipelines have been proposed to either expand the structural breadth of an ensemble or bias the prediction toward a desired conformational state. Here, we analyze how those pipelines work, what they can and cannot predict, and future directions.

Protein structure prediction with energy minimization and deep learning approaches

In this paper we discuss the advantages and problems of two alternatives for ab initio protein structure prediction. On one hand, recent approaches based on deep learning, which have significantly improved prediction results for a wide variety of proteins, are discussed. On the other hand, methods based on protein conformational energy minimization and with different search strategies are analyzed. In this latter case, our methods based on a memetic combination between differential evolution and the fragment replacement technique are included, incorporating also the possibility of niching in the evolutionary search. Different proteins have been used to analyze the pros and cons in both approaches, proposing possibilities of integration of both alternatives.

The nucleotide excision repair proteins through the lens of molecular dynamics simulations

Mutations that affect the proteins responsible for the nucleotide excision repair (NER) pathway can lead to diseases such as xeroderma pigmentosum, trichothiodystrophy, Cockayne syndrome, and Cerebro-oculo-facio-skeletal syndrome. Hence, understanding their molecular behavior is needed to elucidate these diseases' phenotypes and how the NER pathway is organized and coordinated. Molecular dynamics techniques enable the study of different protein conformations, adaptable to any research question, shedding light on the dynamics of biomolecules. However, as important as they are, molecular dynamics studies focused on DNA repair pathways are still becoming more widespread. Currently, there are no review articles compiling the advancements made in molecular dynamics approaches applied to NER and discussing: (i) how this technique is currently employed in the field of DNA repair, focusing on NER proteins; (ii) which technical setups are being employed, their strengths and limitations; (iii) which insights or information are they providing to understand the NER pathway or NER-associated proteins; (iv) which open questions would be suited for this technique to answer; and (v) where can we go from here. These questions become even more crucial considering the numerous 3D structures published regarding the NER pathway's proteins in recent years. In this work, we tackle each one of these questions, revising and critically discussing the results published in the context of the NER pathway.

Unstructured Polypeptides as a Versatile Drug Delivery Technology

Although polyethylene glycol (PEG), or "PEGylation" has become a widely applied approach for improving the efficiency of drug delivery, the immunogenicity and non-biodegradability of this synthetic polymer have prompted an evident need for alternatives. To overcome these caveats and to mimic PEG -or other natural or synthetic polymers- for the purpose of drug half-life extension, unstructured polypeptides are designed. Due to their tunable length, biodegradability, low immunogenicity and easy production, unstructured polypeptides have the potential to replace PEG as the preferred technology for therapeutic protein/peptide delivery. This review provides an overview of the evolution of unstructured polypeptides, starting from natural polypeptides to engineered polypeptides and discusses their characteristics. Then, it is described that unstructured polypeptides have been successfully applied to numerous drugs, including peptides, proteins, antibody fragments, and nanocarriers, for half-life extension. Innovative applications of unstructured peptides as releasable masks, multimolecular adaptors and intracellular delivery carriers are also discussed. Finally, challenges and future perspectives of this promising field are briefly presented. STATEMENT OF SIGNIFICANCE: : Polypeptide fusion technology simulating PEGylation has become an important topic for the development of long-circulating peptide or protein drugs without reduced activity, complex processes, and kidney injury caused by PEG modification. Here we provide a detailed and in-depth review of the recent advances in unstructured polypeptides. In addition to the application of enhanced pharmacokinetic performance, emphasis is placed on polypeptides as scaffolders for the delivery of multiple drugs, and on the preparation of reasonably designed polypeptides to manipulate the performance of proteins and peptides. This review will provide insight into future application of polypeptides in peptide or protein drug development and the design of novel functional polypeptides.