Japanese encephalitis virus hijacks ER-associated degradation regulators for its replication.

Flaviviruses target their replication on membranous structures derived from the ER, where both viral and host proteins play crucial structural and functional roles. Here, we have characterized the involvement of the ER-associated degradation (ERAD) pathway core E3 ligase complex (SEL1L-HRD1) regulator proteins in the replication of Japanese encephalitis virus (JEV). Through high-resolution immunofluorescence imaging of JEV-infected HeLa cells, we observe that the virus replication complexes marked by NS1 strongly colocalize with the ERAD adapter SEL1L, lectin OS9, ER-membrane shuttle factor HERPUD1, E3 ubiquitin ligase HRD1 and rhomboid superfamily member DERLIN1. NS5 positive structures also show strong overlap with SEL1L. While these effectors show significant transcriptional upregulation, their protein levels remain largely stable in infected cells. siRNA mediated depletion of OS9, SEL1L, HERPUD1 and HRD1 significantly inhibit viral RNA replication and titres, with SEL1L depletion showing the maximum attenuation of replication. By performing protein translation arrest experiments, we show that SEL1L, and OS9 are stabilised upon JEV infection. Overall results from this study suggest that these ERAD effector proteins are crucial host-factors for JEV replication.

Genetic Encoding of Phosphorylated Amino Acids into Proteins.

Reversible phosphorylation is a fundamental mechanism for controlling protein function. Despite the critical roles phosphorylated proteins play in physiology and disease, our ability to study individual phospho-proteoforms has been hindered by a lack of versatile methods to efficiently generate homogeneous proteins with site-specific phosphoamino acids or with functional mimics that are resistant to phosphatases. Genetic code expansion (GCE) is emerging as a transformative approach to tackle this challenge, allowing direct incorporation of phosphoamino acids into proteins during translation in response to amber stop codons. This genetic programming of phospho-protein synthesis eliminates the reliance on kinase-based or chemical semisynthesis approaches, making it broadly applicable to diverse phospho-proteoforms. In this comprehensive review, we provide a brief introduction to GCE and trace the development of existing GCE technologies for installing phosphoserine, phosphothreonine, phosphotyrosine, and their mimics, discussing both their advantages as well as their limitations. While some of the technologies are still early in their development, others are already robust enough to greatly expand the range of biologically relevant questions that can be addressed. We highlight new discoveries enabled by these GCE approaches, provide practical considerations for the application of technologies by non-GCE experts, and also identify avenues ripe for further development.

High-throughput tagging of endogenous loci for rapid characterization of protein function.

To facilitate the interrogation of protein function at scale, we have developed high-throughput insertion of tags across the genome (HITAG). HITAG enables users to rapidly produce libraries of cells, each with a different protein of interest C-terminally tagged. HITAG is based on a modified strategy for performing Cas9-based targeted insertions, coupled with an improved approach for selecting properly tagged lines. Analysis of the resulting clones generated by HITAG reveals high tagging specificity, with most successful tagging events being indel free. Using HITAG, we fuse mCherry to a set of 167 stress granule-associated proteins and elucidate the features that drive a subset of proteins to strongly accumulate within these transient RNA-protein granules.

Human epididymis protein 4, a novel potential biomarker for diagnostic and prognosis monitoring of lung cancer.

OBJECTIVE: This study aimed to explore the application value of human epididymis protein 4 (HE4) in diagnosing and monitoring the prognosis of lung cancer. METHODS: First, TCGA (The Cancer Genome Atlas) databases were used to analyze whey-acidic-protein 4-disulfide bond core domain 2 (WFDC2) gene expression levels in lung cancer tissues. Then, a total of 160 individuals were enrolled, categorized into three groups: the lung cancer group (n = 80), the benign lesions group (n = 40), and the healthy controls group (n = 40). Serum HE4 levels and other biomarkers were quantified using an electro-chemiluminescent immunoassay. Additionally, the expression of HE4 in tissues was analyzed through immunohistochemistry (IHC). In vitro cultures of human airway epithelial (human bronchial epithelial [HBE]) cells and various lung cancer cell lines (SPC/PC9/A594/H520) were utilized to detect HE4 levels via western blot (WB). RESULTS: Analysis of the TCGA and UALCAN (The University of Alabama at Birmingham Cancer Data Analysis Portal) databases showed that WFDC2 gene expression levels were upregulated in lung cancer tissues (p < 0.01). Compared with the control group and the benign group, HE4 was significantly higher in the serum of patients with lung cancer (p < 0.001). Receiver operating characteristic (ROC) analysis confirmed that HE4 had better diagnostic efficacy than classical markers in the differential diagnosis of lung cancer and benign lesions and had the highest diagnostic value in lung adenocarcinoma (area under the ROC curve [AUC] = 0.826). HE4 increased in early lung cancer and positively correlated with poor prognosis (p < 0.001). Moreover, the results of WB and IHC revealed that the expression of HE4 was increased in lung cancer cells (SPC/A549/H520) and lung cancer tissues but decreased in PC9 cells with a lack of exon EGFR19 (p < 0.05). CONCLUSION: Serum HE4 emerges as a promising novel biomarker for the diagnosis and prognosis assessment of lung cancer.

Robust genetic codes enhance protein evolvability.

The standard genetic code defines the rules of translation for nearly every life form on Earth. It also determines the amino acid changes accessible via single-nucleotide mutations, thus influencing protein evolvability-the ability of mutation to bring forth adaptive variation in protein function. One of the most striking features of the standard genetic code is its robustness to mutation, yet it remains an open question whether such robustness facilitates or frustrates protein evolvability. To answer this question, we use data from massively parallel sequence-to-function assays to construct and analyze 6 empirical adaptive landscapes under hundreds of thousands of rewired genetic codes, including those of codon compression schemes relevant to protein engineering and synthetic biology. We find that robust genetic codes tend to enhance protein evolvability by rendering smooth adaptive landscapes with few peaks, which are readily accessible from throughout sequence space. However, the standard genetic code is rarely exceptional in this regard, because many alternative codes render smoother landscapes than the standard code. By constructing low-dimensional visualizations of these landscapes, which each comprise more than 16 million mRNA sequences, we show that such alternative codes radically alter the topological features of the network of high-fitness genotypes. Whereas the genetic codes that optimize evolvability depend to some extent on the detailed relationship between amino acid sequence and protein function, we also uncover general design principles for engineering nonstandard genetic codes for enhanced and diminished evolvability, which may facilitate directed protein evolution experiments and the bio-containment of synthetic organisms, respectively.

Mutual annotation-based prediction of protein domain functions with Domain2GO.

Identifying unknown functional properties of proteins is essential for understanding their roles in both health and disease states. The domain composition of a protein can reveal critical information in this context, as domains are structural and functional units that dictate how the protein should act at the molecular level. The expensive and time-consuming nature of wet-lab experimental approaches prompted researchers to develop computational strategies for predicting the functions of proteins. In this study, we proposed a new method called Domain2GO that infers associations between protein domains and function-defining gene ontology (GO) terms, thus redefining the problem as domain function prediction. Domain2GO uses documented protein-level GO annotations together with proteins' domain annotations. Co-annotation patterns of domains and GO terms in the same proteins are examined using statistical resampling to obtain reliable associations. As a use-case study, we evaluated the biological relevance of examples selected from the Domain2GO-generated domain-GO term mappings via literature review. Then, we applied Domain2GO to predict unknown protein functions by propagating domain-associated GO terms to proteins annotated with these domains. For function prediction performance evaluation and comparison against other methods, we employed Critical Assessment of Function Annotation 3 (CAFA3) challenge datasets. The results demonstrated the high potential of Domain2GO, particularly for predicting molecular function and biological process terms, along with advantages such as producing interpretable results and having an exceptionally low computational cost. The approach presented here can be extended to other ontologies and biological entities to investigate unknown relationships in complex and large-scale biological data. The source code, datasets, results, and user instructions for Domain2GO are available at https://github.com/HUBioDataLab/Domain2GO. Additionally, we offer a user-friendly online tool at https://huggingface.co/spaces/HUBioDataLab/Domain2GO, which simplifies the prediction of functions of previously unannotated proteins solely using amino acid sequences.

Genetic code robustness and protein evolvability are correlated and protein-specific.

The relationship between genetic code robustness and protein evolvability is unknown. A new study in PLOS Biology using in silico rewiring of genetic codes and functional protein data identified a positive correlation between code robustness and protein evolvability that is protein-specific.

Predicting hotspots for disease-causing single nucleotide variants using sequences-based coevolution, network analysis, and machine learning.

To enable personalized medicine, it is important yet highly challenging to accurately predict disease-causing mutations in target proteins at high throughput. Previous computational methods have been developed using evolutionary information in combination with various biochemical and structural features of protein residues to discriminate neutral vs. deleterious mutations. However, the power of these methods is often limited because they either assume known protein structures or treat residues independently without fully considering their interactions. To address the above limitations, we build upon recent progress in machine learning, network analysis, and protein language models, and develop a sequences-based variant site prediction workflow based on the protein residue contact networks: 1. We employ and integrate various methods of building protein residue networks using state-of-the-art coevolution analysis tools (RaptorX, DeepMetaPSICOV, and SPOT-Contact) powered by deep learning. 2. We use machine learning algorithms (Random Forest, Gradient Boosting, and Extreme Gradient Boosting) to optimally combine 20 network centrality scores to jointly predict key residues as hot spots for disease mutations. 3. Using a dataset of 107 proteins rich in disease mutations, we rigorously evaluate the network scores individually and collectively (via machine learning). This work supports a promising strategy of combining an ensemble of network scores based on different coevolution analysis methods (and optionally predictive scores from other methods) via machine learning to predict hotspot sites of disease mutations, which will inform downstream applications of disease diagnosis and targeted drug design.

Insertion-Deletion Events Are Depleted in Protein Regions with Predicted Secondary Structure.

A fundamental goal in evolutionary biology and population genetics is to understand how selection shapes the fate of new mutations. Here, we test the null hypothesis that insertion-deletion (indel) events in protein-coding regions occur randomly with respect to secondary structures. We identified indels across 11,444 sequence alignments in mouse, rat, human, chimp, and dog genomes and then quantified their overlap with four different types of secondary structure-alpha helices, beta strands, protein bends, and protein turns-predicted by deep-learning methods of AlphaFold2. Indels overlapped secondary structures 54% as much as expected and were especially underrepresented over beta strands, which tend to form internal, stable regions of proteins. In contrast, indels were enriched by 155% over regions without any predicted secondary structures. These skews were stronger in the rodent lineages compared to the primate lineages, consistent with population genetic theory predicting that natural selection will be more efficient in species with larger effective population sizes. Nonsynonymous substitutions were also less common in regions of protein secondary structure, although not as strongly reduced as in indels. In a complementary analysis of thousands of human genomes, we showed that indels overlapping secondary structure segregated at significantly lower frequency than indels outside of secondary structure. Taken together, our study shows that indels are selected against if they overlap secondary structure, presumably because they disrupt the tertiary structure and function of a protein.

DeepSS2GO: protein function prediction from secondary structure.

Predicting protein function is crucial for understanding biological life processes, preventing diseases and developing new drug targets. In recent years, methods based on sequence, structure and biological networks for protein function annotation have been extensively researched. Although obtaining a protein in three-dimensional structure through experimental or computational methods enhances the accuracy of function prediction, the sheer volume of proteins sequenced by high-throughput technologies presents a significant challenge. To address this issue, we introduce a deep neural network model DeepSS2GO (Secondary Structure to Gene Ontology). It is a predictor incorporating secondary structure features along with primary sequence and homology information. The algorithm expertly combines the speed of sequence-based information with the accuracy of structure-based features while streamlining the redundant data in primary sequences and bypassing the time-consuming challenges of tertiary structure analysis. The results show that the prediction performance surpasses state-of-the-art algorithms. It has the ability to predict key functions by effectively utilizing secondary structure information, rather than broadly predicting general Gene Ontology terms. Additionally, DeepSS2GO predicts five times faster than advanced algorithms, making it highly applicable to massive sequencing data. The source code and trained models are available at https://github.com/orca233/DeepSS2GO.

Scoring alignments by embedding vector similarity.

Sequence similarity is of paramount importance in biology, as similar sequences tend to have similar function and share common ancestry. Scoring matrices, such as PAM or BLOSUM, play a crucial role in all bioinformatics algorithms for identifying similarities, but have the drawback that they are fixed, independent of context. We propose a new scoring method for amino acid similarity that remedies this weakness, being contextually dependent. It relies on recent advances in deep learning architectures that employ self-supervised learning in order to leverage the power of enormous amounts of unlabelled data to generate contextual embeddings, which are vector representations for words. These ideas have been applied to protein sequences, producing embedding vectors for protein residues. We propose the E-score between two residues as the cosine similarity between their embedding vector representations. Thorough testing on a wide variety of reference multiple sequence alignments indicate that the alignments produced using the new $E$-score method, especially ProtT5-score, are significantly better than those obtained using BLOSUM matrices. The new method proposes to change the way alignments are computed, with far-reaching implications in all areas of textual data that use sequence similarity. The program to compute alignments based on various $E$-scores is available as a web server at e-score.csd.uwo.ca. The source code is freely available for download from github.com/lucian-ilie/E-score.

The physical and evolutionary energy landscapes of devolved protein sequences corresponding to pseudogenes.

Protein evolution is guided by structural, functional, and dynamical constraints ensuring organismal viability. Pseudogenes are genomic sequences identified in many eukaryotes that lack translational activity due to sequence degradation and thus over time have undergone "devolution." Previously pseudogenized genes sometimes regain their protein-coding function, suggesting they may still encode robust folding energy landscapes despite multiple mutations. We study both the physical folding landscapes of protein sequences corresponding to human pseudogenes using the Associative Memory, Water Mediated, Structure and Energy Model, and the evolutionary energy landscapes obtained using direct coupling analysis (DCA) on their parent protein families. We found that generally mutations that have occurred in pseudogene sequences have disrupted their native global network of stabilizing residue interactions, making it harder for them to fold if they were translated. In some cases, however, energetic frustration has apparently decreased when the functional constraints were removed. We analyzed this unexpected situation for Cyclophilin A, Profilin-1, and Small Ubiquitin-like Modifier 2 Protein. Our analysis reveals that when such mutations in the pseudogene ultimately stabilize folding, at the same time, they likely alter the pseudogenes' former biological activity, as estimated by DCA. We localize most of these stabilizing mutations generally to normally frustrated regions required for binding to other partners.

Assessing the functional impact of protein binding site definition.

Many biomedical applications, such as classification of binding specificities or bioengineering, depend on the accurate definition of protein binding interfaces. Depending on the choice of method used, substantially different sets of residues can be classified as belonging to the interface of a protein. A typical approach used to verify these definitions is to mutate residues and measure the impact of these changes on binding. Besides the lack of exhaustive data, this approach also suffers from the fundamental problem that a mutation introduces an unknown amount of alteration into an interface, which potentially alters the binding characteristics of the interface. In this study we explore the impact of alternative binding site definitions on the ability of a protein to recognize its cognate ligand using a pharmacophore approach, which does not affect the interface. The study also shows that methods for protein binding interface predictions should perform above approximately F-score = 0.7 accuracy level to capture the biological function of a protein.

The implications of APOBEC3-mediated C-to-U RNA editing for human disease.

Intra-organism biodiversity is thought to arise from epigenetic modification of constituent genes and post-translational modifications of translated proteins. Here, we show that post-transcriptional modifications, like RNA editing, may also contribute. RNA editing enzymes APOBEC3A and APOBEC3G catalyze the deamination of cytosine to uracil. RNAsee (RNA site editing evaluation) is a computational tool developed to predict the cytosines edited by these enzymes. We find that 4.5% of non-synonymous DNA single nucleotide polymorphisms that result in cytosine to uracil changes in RNA are probable sites for APOBEC3A/G RNA editing; the variant proteins created by such polymorphisms may also result from transient RNA editing. These polymorphisms are associated with over 20% of Medical Subject Headings across ten categories of disease, including nutritional and metabolic, neoplastic, cardiovascular, and nervous system diseases. Because RNA editing is transient and not organism-wide, future work is necessary to confirm the extent and effects of such editing in humans.

Genotype sampling for deep-learning assisted experimental mapping of a combinatorially complete fitness landscape.

MOTIVATION: Experimental characterization of fitness landscapes, which map genotypes onto fitness, is important for both evolutionary biology and protein engineering. It faces a fundamental obstacle in the astronomical number of genotypes whose fitness needs to be measured for any one protein. Deep learning may help to predict the fitness of many genotypes from a smaller neural network training sample of genotypes with experimentally measured fitness. Here I use a recently published experimentally mapped fitness landscape of more than 260 000 protein genotypes to ask how such sampling is best performed. RESULTS: I show that multilayer perceptrons, recurrent neural networks, convolutional networks, and transformers, can explain more than 90% of fitness variance in the data. In addition, 90% of this performance is reached with a training sample comprising merely ≈103 sequences. Generalization to unseen test data is best when training data is sampled randomly and uniformly, or sampled to minimize the number of synonymous sequences. In contrast, sampling to maximize sequence diversity or codon usage bias reduces performance substantially. These observations hold for more than one network architecture. Simple sampling strategies may perform best when training deep learning neural networks to map fitness landscapes from experimental data. AVAILABILITY AND IMPLEMENTATION: The fitness landscape data analyzed here is publicly available as described previously (Papkou et al. 2023). All code used to analyze this landscape is publicly available at https://github.com/andreas-wagner-uzh/fitness_landscape_sampling.

Polyphosphate attachment to lysine repeats is a non-covalent protein modification.

Polyphosphate (polyP) is a chain of inorganic phosphate that is present in all domains of life and affects diverse cellular phenomena, ranging from blood clotting to cancer. A study by Azevedo et al. described a protein modification whereby polyP is attached to lysine residues within polyacidic serine and lysine (PASK) motifs via what the authors claimed to be covalent phosphoramidate bonding. This was based largely on the remarkable ability of the modification to survive extreme denaturing conditions. Our study demonstrates that lysine polyphosphorylation is non-covalent, based on its sensitivity to ionic strength and lysine protonation and absence of phosphoramidate bond formation, as analyzed via 31P NMR. Ionic interaction with lysine residues alone is sufficient for polyP modification, and we present a new list of non-PASK lysine repeat proteins that undergo polyP modification. This work clarifies the biochemistry of polyP-lysine modification, with important implications for both studying and modulating this phenomenon. This Matters Arising paper is in response to Azevedo et al. (2015), published in Molecular Cell. See also the Matters Arising Response by Azevedo et al. (2024), published in this issue.

Analysis of AlphaMissense data in different protein groups and structural context.

Single amino acid substitutions can profoundly affect protein folding, dynamics, and function. The ability to discern between benign and pathogenic substitutions is pivotal for therapeutic interventions and research directions. Given the limitations in experimental examination of these variants, AlphaMissense has emerged as a promising predictor of the pathogenicity of missense variants. Since heterogenous performance on different types of proteins can be expected, we assessed the efficacy of AlphaMissense across several protein groups (e.g. soluble, transmembrane, and mitochondrial proteins) and regions (e.g. intramembrane, membrane interacting, and high confidence AlphaFold segments) using ClinVar data for validation. Our comprehensive evaluation showed that AlphaMissense delivers outstanding performance, with MCC scores predominantly between 0.6 and 0.74. We observed low performance on disordered datasets and ClinVar data related to the CFTR ABC protein. However, a superior performance was shown when benchmarked against the high quality CFTR2 database. Our results with CFTR emphasizes AlphaMissense's potential in pinpointing functional hot spots, with its performance likely surpassing benchmarks calculated from ClinVar and ProteinGym datasets.

DeepPPAPredMut: deep ensemble method for predicting the binding affinity change in protein-protein complexes upon mutation.

MOTIVATION: Protein-protein interactions underpin many cellular processes and their disruption due to mutations can lead to diseases. With the evolution of protein structure prediction methods like AlphaFold2 and the availability of extensive experimental affinity data, there is a pressing need for updated computational tools that can efficiently predict changes in binding affinity caused by mutations in protein-protein complexes. RESULTS: We developed a deep ensemble model that leverages protein sequences, predicted structure-based features, and protein functional classes to accurately predict the change in binding affinity due to mutations. The model achieved a correlation of 0.97 and a mean absolute error (MAE) of 0.35 kcal/mol on the training dataset, and maintained robust performance on the test set with a correlation of 0.72 and a MAE of 0.83 kcal/mol. Further validation using Leave-One-Out Complex (LOOC) cross-validation exhibited a correlation of 0.83 and a MAE of 0.51 kcal/mol, indicating consistent performance. AVAILABILITY AND IMPLEMENTATION: https://web.iitm.ac.in/bioinfo2/DeepPPAPredMut/index.html.

Expression of Insl3 Protein in Adult Danio rerio.

Insulin-like peptide 3 (INSL3) is a biomarker for Leydig cells in the testes of vertebrates, and it is principally involved in spermatogenesis through specific binding with the RXFP2 receptor. This study reports the insl3 gene transcript and the Insl3 prepropeptide expression in both non-reproductive and reproductive tissues of Danio rerio. An immunohistochemistry analysis shows that the hormone is present at a low level in the Leydig cells and germ cells at all stages of Danio rerio testis differentiation. Considering that the insl3 gene is transcribed in Leydig cells, our results highlight an autocrine and paracrine function of this hormone in the Danio rerio testis, adding new information on the Insl3 mode of action in reproduction. We also show that Insl3 and Rxfp2 belonging to Danio rerio and other vertebrate species share most of the amino acid residues involved in the ligand-receptor interaction and activation, suggesting a conserved mechanism of action during vertebrate evolution.

Associations of the obesity gene FTO variant with complications and comorbidities in patients with type 1 diabetes.

BACKGROUND AND AIMS: Because FTO gene is connected with the risk of obesity, cardiovascular disease and hypertension, as well as type 2 diabetes, we hypothesize that the rs9939609 FTO polymorphism may affect type 1 diabetes (T1D) complications and comorbidities. METHODS: We have investigated the associations of the FTO gene variant with the T1D and its complications and comorbidities, as well as the serum levels of pro- and anti-inflammatory markers and lipid profiles. RESULTS: The key results of our study are as follows: (1) the rs9939609 FTO polymorphism does not predispose individuals to T1D; (2) AA genotype is associated with an increased risk of overweight and obesity, retinopathy, hypertension, dyslipidemia and celiac disease; (3) AT genotype is associated with a decreased risk of retinopathy and celiac disease, whereas TT genotype is connected with decreased risk of dyslipidemia; (4) the FTO rs9939609 polymorphism affects the inflammatory status as well as lipid profile in T1D patients. CONCLUSIONS: Our results, for the first time, comprehensively indicate that the rs9939609 FTO polymorphism could be considered a genetic marker for increased susceptibility to T1D complications and comorbidities as well as suggests importance of FTO-mediated pathways in their etiology.