NetSurfP-2.0: Improved prediction of protein structural features by integrated deep learning

Prediction of deleterious non-synonymous SNPs of human MDC1 gene: an in silico approach

MDC1 (Mediator of DNA damage Checkpoint protein 1) functions to facilitate the localization of numerous DNA damage response (DDR) components to DNA double-strand break sites. MDC1 is an integral component in preserving genomic stability and appropriate DDR regulation. There haven't been systematic investigations of MDC1 mutations that induce cancer and genomic instability. Variations in nsSNPs have the potential to modify the protein chemistry and their function. Describing functional SNPs in disease-associated genes presents a significant conundrum for investigators, it is possible to assess potential functional SNPs before conducting larger population examinations. Multiple sequences and structure-based bioinformatics strategies were implemented in the current in-silico investigation to discern potential nsSNPs of the MDC1 genes. The nsSNPs were identified with SIFT, SNAP2, Align GVGD, PolyPhen-2, and PANTHER, and their stability was determined with MUpro. The conservation, solvent accessibility, and structural effects of the mutations were identified with ConSurf, NetSurfP-2.0, and SAAFEC-SEQ respectively. Cancer-related analysis of the nsSNPs was conducted using cBioPortal and TCGA web servers. The present study appraised five nsSNPs (P1426T, P69S, P194R, P203L, and H131Y) as probably mutilating due to their existence in highly conserved regions and propensity to deplete protein stability. The nsSNPs P194R, P203L, and H131Y were concluded as deleterious and possibly damaging from the 5 prediction tools. The functional nsSNP P194R mutation is associated with skin cutaneous melanoma while no significant records were found for other nsSNPs. The present study concludes that the highly deleterious P194R mutations can potentially induce genomic instability and contribute to various cancers' pathogenesis. Developing drugs targeting these mutations can undoubtedly be advantageous in large population-based studies, particularly in the development of precision medicine.

The rate of adaptive molecular evolution in wild and domesticated Saccharomyces cerevisiae populations

Through its fermentative capacities, Saccharomyces cerevisiae was central in the development of civilisation during the Neolithic period, and the yeast remains of importance in industry and biotechnology, giving rise to bona fide domesticated populations. Here, we conduct a population genomic study of domesticated and wild populations of S. cerevisiae. Using coalescent analyses, we report that the effective population size of yeast populations decreased since the divergence with S. paradoxus. We fitted models of distributions of fitness effects to infer the rate of adaptive (ω a ) and non-adaptive (ω na ) non-synonymous substitutions in protein-coding genes. We report an overall limited contribution of positive selection to S. cerevisiae protein evolution, albeit with higher rates of adaptive evolution in wild compared to domesticated populations. Our analyses revealed the signature of background selection and possibly Hill-Robertson interference, as recombination was found to be negatively correlated withω na and positively correlated withω a . However, the effect of recombination onω a was found to be labile, as it is only apparent after removing the impact of codon usage bias on the synonymous site frequency spectrum and disappears if we control for the correlation withω na , suggesting that it could be an artefact of the decreasing population size. Furthermore, the rate of adaptive non-synonymous substitutions is significantly correlated with the residue solvent exposure, a relation that cannot be explained by the population's demography. Together, our results provide a detailed characterisation of adaptive mutations in protein-coding genes across S. cerevisiae populations.

Caffeine improves mitochondrial dysfunction in the white matter of neonatal rats with hypoxia-ischemia through deacetylation: a proteomic analysis of lysine acetylation

Aims

White matter damage (WMD) is linked to both cerebral palsy and cognitive deficits in infants born prematurely. The focus of this study was to examine how caffeine influences the acetylation of proteins within the neonatal white matter and to evaluate its effectiveness in treating white matter damage caused by hypoxia-ischemia.

Main methods

We employed a method combining affinity enrichment with advanced liquid chromatography and mass spectrometry to profile acetylation in proteins from the white matter of neonatal rats grouped into control (Sham), hypoxic-ischemic (HI), and caffeine-treated (Caffeine) groups.

Key findings

Our findings included 1,999 sites of lysine acetylation across 1,123 proteins, with quantifiable changes noted in 1,342 sites within 689 proteins. Analysis of these patterns identified recurring sequences adjacent to the acetylation sites, notably YKacN, FkacN, and G * * * GkacS. Investigation into the biological roles of these proteins through Gene Ontology analysis indicated their involvement in a variety of cellular processes, predominantly within mitochondrial locations. Further analysis indicated that the acetylation of tau (Mapt), a protein associated with microtubules, was elevated in the HI condition; however, caffeine treatment appeared to mitigate this over-modification, thus potentially aiding in reducing oxidative stress, inflammation in the nervous system, and improving mitochondrial health. Caffeine inhibited acetylated Mapt through sirtuin 2 (SITR2), promoted Mapt nuclear translocation, and improved mitochondrial dysfunction, which was subsequently weakened by the SIRT2 inhibitor, AK-7.

Significance

Caffeine-induced changes in lysine acetylation may play a key role in improving mitochondrial dysfunction and inhibiting oxidative stress and neuroinflammation.

Enhancing subcellular protein localization mapping analysis using Sc2promap utilizing attention mechanisms

BACKGROUND

Aberrant protein localization is a prominent feature in many human diseases and can have detrimental effects on the function of specific tissues and organs. High-throughput technologies, which continue to advance with iterations of automated equipment and the development of bioinformatics, enable the acquisition of large-scale data that are more pattern-rich, allowing for the use of a wider range of methods to extract useful patterns and knowledge from them.

METHODS

The proposed sc2promap (Spatial and Channel for SubCellular Protein Localization Mapping) model, designed to proficiently extract meaningful features from a vast repository of single-channel grayscale protein images for the purposes of protein localization analysis and clustering. Sc2promap incorporates a prediction head component enriched with supplementary protein annotations, along with the integration of a spatial-channel attention mechanism within the encoder to enables the generation of high-resolution protein localization maps that encapsulate the fundamental characteristics of cells, including elemental cellular localizations such as nuclear and non-nuclear domains.

RESULTS

Qualitative and quantitative comparisons were conducted across internal and external clustering evaluation metrics, as well as various facets of the clustering results. The study also explored different components of the model. The research outcomes conclusively indicate that, in comparison to previous methods, Sc2promap exhibits superior performance.

CONCLUSIONS

The amalgamation of the attention mechanism and prediction head components has led the model to excel in protein localization clustering and analysis tasks.

GENERAL SIGNIFICANCE

The model effectively enhances the capability to extract features and knowledge from protein fluorescence images.

Salp-J Colony Optimization-based advanced hybrid ensemble deep predictor with LSTM for protein structure prediction

Protein structure prediction (PSP) is a key concern in computational biology, which is considered a challenging task that is vital to determine the structure and the protein function since each protein possesses a definite shape, whereas the protein secondary structure prediction (PSSP) is the foundation for three-dimensional PSP. An Advanced hybrid ensemble deep predictor is utilized for predicting the structure of a protein using Long-Short Term Memory (LSTM), in which the performance of the predictor is improved for obtaining the features through the Salp-J Colony Optimization, which is developed by integrating the features of three optimizations the exploration behavior of Ulmaris, the immune system of virus colony and the teamwork of salp for solution update that helps to predict the accurate protein structure. The proposed method achieved the value of 99.1% accuracy, 99.5% sensitivity, 98.85% specificity, and 0.9% error at the 80% of training percentage 90 using CullPDB. Similarly, in Protein Net, the attained value of accuracy is 97.27%, sensitivity is 98.13%, specificity is 97%, and error is 2.7% concerning training percentage 90%.Communicated by Ramaswamy H. Sarma.

Evaluate the potential use of TonB-dependent receptor protein as a subunit vaccine against Aeromonas veronii infection in Nile tilapia (Oreochromis niloticus)

Aeromonas veronii is an emerging bacterial pathogen that causes serious systemic infections in cultured Nile tilapia (Oreochromis niloticus), leading to massive deaths. Therefore, there is an urgent need to identify effective vaccine candidates to control the spread of this emerging disease. TonB-dependent receptor (Tdr) of A. veronii, which plays a role in the virulence factor of the organism, could be useful in terms of protective antigens for vaccine development. This study aims to evaluate the potential use of Tdr protein as a novel subunit vaccine against A. veronii infection in Nile tilapia. The Tdr gene from A. veronii was cloned into the pET28b expression vector, and the recombinant protein was subsequently produced in Escherichia coli strain BL21 (DE3). Tdr was expressed as an insoluble protein and purified by affinity chromatography. Antigenicity test indicated that this protein was recognized by serum from A. veronii infected fish. When Nile tilapia were immunized with the Tdr protein, specific antibody levels increased significantly (p-value <0.05) at 7 days post-immunization (dpi), and peaked at 21 dpi compared to antibody levels at 0 dpi. Furthermore, bacterial agglutination activity was observed in the fish serum immunized with the Tdr protein, indicating that specific antibodies in the serum can detect Tdr on the bacterial cell surface. These results suggest that Tdr protein has potential as a vaccine candidate. However, challenging tests with A.veronii in Nile tilapia needs to be investigated to thoroughly evaluate its protective efficacy for future applications.

Evaluating generalizability of artificial intelligence models for molecular datasets

Deep learning has made rapid advances in modeling molecular sequencing data. Despite achieving high performance on benchmarks, it remains unclear to what extent deep learning models learn general principles and generalize to previously unseen sequences. Benchmarks traditionally interrogate model generalizability by generating metadata based (MB) or sequence-similarity based (SB) train and test splits of input data before assessing model performance. Here, we show that this approach mischaracterizes model generalizability by failing to consider the full spectrum of cross-split overlap, i.e., similarity between train and test splits. We introduce Spectra, a spectral framework for comprehensive model evaluation. For a given model and input data, Spectra plots model performance as a function of decreasing cross-split overlap and reports the area under this curve as a measure of generalizability. We apply Spectra to 18 sequencing datasets with associated phenotypes ranging from antibiotic resistance in tuberculosis to protein-ligand binding to evaluate the generalizability of 19 state-of-the-art deep learning models, including large language models, graph neural networks, diffusion models, and convolutional neural networks. We show that SB and MB splits provide an incomplete assessment of model generalizability. With Spectra, we find as cross-split overlap decreases, deep learning models consistently exhibit a reduction in performance in a task- and model-dependent manner. Although no model consistently achieved the highest performance across all tasks, we show that deep learning models can generalize to previously unseen sequences on specific tasks. Spectra paves the way toward a better understanding of how foundation models generalize in biology.

Structural effect of the H992D/H418D mutation of angiotensin-converting enzyme in the Indian population: implications for health and disease

The Non synonymous SNPs (nsSNPs) of the renin-angiotensin-system (RAS) pathway, unique to the Indian population were investigated in view of its importance as an endocrine system. nsSNPs of the RAS pathway genes were mined from the IndiGenome database. Damaging nsSNPs were predicted using SIFT, PredictSNP, SNP and GO, Snap2 and Protein Variation Effect Analyzer. Loss of function was predicted based on protein stability change using I mutant, PremPS and CONSURF. The structural impact of the nsSNPs was predicted using HOPE and Missense3d followed by modeling, refinement, and energy minimization. Molecular Dynamics studies were carried out using Gromacsv2021.1. 23 Indian nsSNPs of the RAS pathway genes were selected for structural analysis and 8 were predicted to be damaging. Further sequence analysis showed that HEMGH zinc binding motif changes to HEMGD in somatic ACE-C domain (sACE-C) H992D and Testis ACE (tACE) H418D resulted in loss of zinc coordination, which is essential for enzymatic activity in this metalloprotease. There was a loss of internal interactions around the zinc coordination residues in the protein structural network. This was also confirmed by Principal Component Analysis, Free Energy Landscape and residue contact maps. Both mutations lead to broadening of the AngI binding cavity. The H992D mutation in sACE-C is likely to be favorable for cardiovascular health, but may lead to renal abnormalities with secondary impact on the heart. H418D in tACE is potentially associated with male infertility.Communicated by Ramaswamy H. Sarma.

Pathological mutations reveal the key role of the cytosolic iRhom2 N-terminus for phosphorylation-independent 14-3-3 interaction and ADAM17 binding, stability, and activity

The protease ADAM17 plays an important role in inflammation and cancer and is regulated by iRhom2. Mutations in the cytosolic N-terminus of human iRhom2 cause tylosis with oesophageal cancer (TOC). In mice, partial deletion of the N-terminus results in a curly hair phenotype (cub). These pathological consequences are consistent with our findings that iRhom2 is highly expressed in keratinocytes and in oesophageal cancer. Cub and TOC are associated with hyperactivation of ADAM17-dependent EGFR signalling. However, the underlying molecular mechanisms are not understood. We have identified a non-canonical, phosphorylation-independent 14-3-3 interaction site that encompasses all known TOC mutations. Disruption of this site dysregulates ADAM17 activity. The larger cub deletion also includes the TOC site and thus also dysregulated ADAM17 activity. The cub deletion, but not the TOC mutation, also causes severe reductions in stimulated shedding, binding, and stability of ADAM17, demonstrating the presence of additional regulatory sites in the N-terminus of iRhom2. Overall, this study contrasts the TOC and cub mutations, illustrates their different molecular consequences, and reveals important key functions of the iRhom2 N-terminus in regulating ADAM17.

ACPScanner: Prediction of Anticancer Peptides by Integrated Machine Learning Methodologies

Novel therapeutic alternatives for cancer treatment are increasingly attracting global research attention. Although chemotherapy remains a primary clinical solution, it often results in significant side effects for patients. In recent years, anticancer peptides (ACPs) have emerged as promising candidates for highly specific anticancer drugs, and a number of computational approaches have been developed to identify ACPs. However, existing methods do not recognize specific types of anticancer function. In this article, we propose ACPScanner, an integrated approach to predict ACPs and non-ACPs at first and then predict several specific activity types for potential ACPs. We incorporate sequential, physicochemical properties, secondary structural information, and deep representation learning embeddings which are generated from artificial intelligence methods to build feature space. Customized deep learning and statistical learning methods are combined to form an integral architecture for the comprehensive two-level prediction task. To the best of our knowledge, ACPScanner is the first approach for specific ACP activity prediction. The comparative evaluation illustrates that ACPScanner achieves competitive prediction performance in both prediction phases in independent testings. We establish a web server at http://acpscanner.denglab.org to provide convenient usage of ACPScanner and make the predictive framework, source code, and data sets publicly available.

DiscoTope-3.0: improved B-cell epitope prediction using inverse folding latent representations

Accurate computational identification of B-cell epitopes is crucial for the development of vaccines, therapies, and diagnostic tools. However, current structure-based prediction methods face limitations due to the dependency on experimentally solved structures. Here, we introduce DiscoTope-3.0, a markedly improved B-cell epitope prediction tool that innovatively employs inverse folding structure representations and a positive-unlabelled learning strategy, and is adapted for both solved and predicted structures. Our tool demonstrates a considerable improvement in performance over existing methods, accurately predicting linear and conformational epitopes across multiple independent datasets. Most notably, DiscoTope-3.0 maintains high predictive performance across solved, relaxed and predicted structures, alleviating the need for experimental structures and extending the general applicability of accurate B-cell epitope prediction by 3 orders of magnitude. DiscoTope-3.0 is made widely accessible on two web servers, processing over 100 structures per submission, and as a downloadable package. In addition, the servers interface with RCSB and AlphaFoldDB, facilitating large-scale prediction across over 200 million cataloged proteins. DiscoTope-3.0 is available at: https://services.healthtech.dtu.dk/service.php?DiscoTope-3.0.

Most Monogenic Disorders Are Caused by Mutations Altering Protein Folding Free Energy

Revealing the molecular effect that pathogenic missense mutations have on the corresponding protein is crucial for developing therapeutic solutions. This is especially important for monogenic diseases since, for most of them, there is no treatment available, while typically, the treatment should be provided in the early development stages. This requires fast targeted drug development at a low cost. Here, we report an updated database of monogenic disorders (MOGEDO), which includes 768 proteins and the corresponding 2559 pathogenic and 1763 benign mutations, along with the functional classification of the corresponding proteins. Using the database and various computational tools that predict folding free energy change (ΔΔG), we demonstrate that, on average, 70% of pathogenic cases result in decreased protein stability. Such a large fraction indicates that one should aim at in silico screening for small molecules stabilizing the structure of the mutant protein. We emphasize that knowledge of ΔΔG is essential because one wants to develop stabilizers that compensate for ΔΔG, but do not make protein over-stable, since over-stable protein may be dysfunctional. We demonstrate that, by using ΔΔG and predicted solvent exposure of the mutation site, one can develop a predictive method that distinguishes pathogenic from benign mutations with a success rate even better than some of the leading pathogenicity predictors. Furthermore, hydrophobic-hydrophobic mutations have stronger correlations between folding free energy change and pathogenicity compared with others. Also, mutations involving Cys, Gly, Arg, Trp, and Tyr amino acids being replaced by any other amino acid are more likely to be pathogenic. To facilitate further detection of pathogenic mutations, the wild type of amino acids in the 768 proteins mentioned above was mutated to other 19 residues (14,847,817 mutations), the ΔΔG was calculated with SAAFEC-SEQ, and 5,506,051 mutations were predicted to be pathogenic.

Structure-based evaluation of the envelope domain III-nonstructural protein 1 (EDIII-NS1) fusion as a dengue virus vaccine candidate

The lack of effective medicines or vaccines, combined with climate change and other environmental factors, annually subjects a significant proportion of the world's inhabitants to the risk of dengue virus (DENV) infection. These conditions increase the likelihood of exposure to mosquito-borne diseases such as dengue fever. Hence, many research approaches tend to develop efficient vaccine candidates against the dengue virus. Therefore, we used immunoinformatics and bioinformatics to design a construction for developing a candidate vaccine against dengue virus serotypes. In this study, the in silico structure, containing the non-structural protein 1 region (NS1) (consensus and epitope), the envelope domain III protein (EDIII) as the structural part of the virus construction, and the bc-loop of envelope domain II (EDII) as the neutralizing and protected epitope, were employed. We utilized in silico tools to enhance the immunogenicity and effectiveness of dengue virus vaccine candidates. Evaluations included refining and validating physicochemical characteristics, B and T-cell epitopes, homology modeling, and the three-dimensional structure to assess the designed vaccine's quality. In silico results for tertiary structure prediction and validation revealed high-quality modeling for all vaccine constructs. Additionally, the instructed model demonstrated stability throughout molecular dynamics simulation. The results of the immune simulation suggested that the titers of IgG and IgM could be raised to desirable values following injection into in vivo models. It can be concluded that the designed construct effectively induce humoral and cellular immunity and can be proposed as effective vaccine candidate against four dengue serotypes.Communicated by Ramaswamy H. Sarma.

Maternal oxytocin treatment at birth increases epigenetic age in male offspring

Exogenous oxytocin (OT) is widely used to induce or augment labor with little understanding of the impact on offspring development. In rodent models, including the prairie vole (Microtus ochrogaster), it has been shown that oxytocin administered to mothers can affect the nervous system of the offspring with long lasting behavioral effects especially on sociality. Here, we examined the hypothesis that perinatal oxytocin exposure could have epigenetic and transcriptomic consequences. Prairie voles were exposed to exogenous oxytocin, through injections given to the mother just prior to birth, and were studied at the time of weaning. The outcome of this study revealed increased epigenetic age in oxytocin-exposed animals compared to the saline-exposed group. Oxytocin exposure led to 900 differentially methylated CpG sites (annotated to 589 genes), and 2 CpG sites (2 genes) remained significantly different after correction for multiple comparisons. Differentially methylated CpG sites were enriched in genes known to be involved in regulation of gene expression and neurodevelopment. Using RNA-sequencing we also found 217 nominally differentially expressed genes (p<0.05) in nucleus accumbens, a brain region involved in reward circuitry and social behavior; after corrections for multiple comparisons 6 genes remained significantly differentially expressed. Finally, we found that maternal oxytocin administration led to widespread alternative splicing in the nucleus accumbens. These results indicate that oxytocin exposure during birth may have long lasting epigenetic consequences. A need for further investigation of how oxytocin administration impacts development and behavior throughout the lifespan is supported by these outcomes.

PS4: a next-generation dataset for protein single-sequence secondary structure prediction

Protein secondary structure prediction is a subproblem of protein folding. A light-weight algorithm capable of accurately predicting secondary structure from only the protein residue sequence could provide useful input for tertiary structure prediction, alleviating the reliance on multiple sequence alignments typically seen in today's best-performing models. Unfortunately, existing datasets for secondary structure prediction are small, creating a bottleneck. We present PS4, a dataset of 18,731 nonredundant protein chains and their respective secondary structure labels. Each chain is identified, and the dataset is nonredundant against other secondary structure datasets commonly seen in the literature. We perform ablation studies by training secondary structure prediction algorithms on the PS4 training set and obtains state-of-the-art accuracy on the CB513 test set in zero shots.

Global profiling of protein lactylation in Caenorhabditis elegans

Lactylation, as a novel posttranslational modification, is essential for studying the functions and regulation of proteins in physiological and pathological processes, as well as for gaining in-depth knowledge on the occurrence and development of many diseases, including tumors. However, few studies have examined the protein lactylation of one whole organism. Thus, we studied the lactylation of global proteins in Caenorhabditis elegans to obtain an in vivo lactylome. Using an MS-based platform, we identified 1836 Class I (localization probabilities > 0.75) lactylated sites in 487 proteins. Bioinformatics analysis showed that lactylated proteins were mainly located in the cytoplasm and involved in the tricarboxylic acid cycle (TCA cycle) and other metabolic pathways. Then, we evaluated the conservation of lactylation in different organisms. In total, 41 C. elegans proteins were lactylated and homologous to lactylated proteins in humans and rats. Moreover, lactylation on H4K80 was conserved in three species. An additional 238 lactylated proteins were identified in C. elegans for the first time. This study establishes the first lactylome database in C. elegans and provides a basis for studying the role of lactylation.

A novel homozygous PIGO mutation associated with severe infantile epileptic encephalopathy, profound developmental delay and psychomotor retardation: structural and 3D modelling investigations and genotype-phenotype correlation

The PIGO gene encodes the GPI-ethanolamine phosphate transferase 3, which is crucial for the final synthetic step of the glycosylphosphatidylinositol-anchor serving to attach various proteins to their cell surface. These proteins are intrinsic for normal neuronal and embryonic development. In the current research work, a clinical investigation was conducted on a patient from a consanguineous family suffering from epileptic encephalopathy, characterized by severe seizures, developmental delay, hypotonia, ataxia and hyperphosphatasia. Molecular analysis was performed using Whole Exome Sequencing (WES). The molecular investigation revealed a novel homozygous variant c.1132C > T in the PIGO gene, in which a highly conserved Leucine was changed to a Phenylalanine (p.L378F). To investigate the impact of the non-synonymous mutation, a 3D structural model of the PIGO protein was generated using the AlphaFold protein structure database as a resource for template-based tertiary structure modeling. A structural analysis by applying some bioinformatic tools on both variants 378L and 378F models predicted the pathogenicity of the non-synonymous mutation and its potential functional and structural effects on PIGO protein. We also discussed the phenotypic and genotypic variability associated with the PIGO deficiency. To our best knowledge, this is the first report of a patient diagnosed with infantile epileptic encephalopathy showing a high elevation of serum alkaline phosphatase level. Our findings, therefore, widen the genotype and phenotype spectrum of GPI-anchor deficiencies and broaden the cohort of patients with PIGO associated epileptic encephalopathy with an elevated serum alkaline phosphatase level.

Comprehensive analysis of non-synonymous missense SNPs of human galactose mutarotase (GALM) gene: an integrated computational approach

Missense Non-synonymous single nucleotide polymorphisms (nsSNPs) of Galactose Mutarotase (GALM) are associated with the Novel type of Galactosemia (Galactosemia type 4) together with symptoms such as high blood galactose levels and eye cataracts. The objective of the present study was to identify deleterious nsSNPs of GALM recorded on the dbSNP database through comprehensive insilico analysis. Among the 319 missense nsSNPs reported, various insilco tools predicted R78S, R82G, A163E, P210S, Y281C, E307G and F339C as the most deleterious mutations. Structural analysis, PTM analysis and molecular dynamics simulations (MDS) were carried out to understand the effect of these mutations on the structural and physicochemical properties of the GALM protein. The residues R82G and E307G were found to be part of the binding site that resulted in decreased surface accessibility. Replacing the charged wild-type residue with a neutral mutant type affected its substrate binding. All 7 mutations were found to increase the rigidity of the protein structure, which is unfavorable during ligand binding. The mutation F339E made the protein structure more rigid than all the other mutations. Y281 is a phosphorylated site, and therefore, less significant structural changes were observed when compared to other mutations; however, it may have significant differences in the usual functioning of the protein. In summary, the structural and functional analysis of missense SNPs of GALM is important to reduce the number of potential mutations to be evaluated in vitro to understand the association with some genetic diseases.Communicated by Ramaswamy H. Sarma.

GraphPart: homology partitioning for biological sequence analysis

When splitting biological sequence data for the development and testing of predictive models, it is necessary to avoid too-closely related pairs of sequences ending up in different partitions. If this is ignored, performance of prediction methods will tend to be overestimated. Several algorithms have been proposed for homology reduction, where sequences are removed until no too-closely related pairs remain. We present GraphPart, an algorithm for homology partitioning that divides the data such that closely related sequences always end up in the same partition, while keeping as many sequences as possible in the dataset. Evaluation of GraphPart on Protein, DNA and RNA datasets shows that it is capable of retaining a larger number of sequences per dataset, while providing homology separation on a par with reduction approaches.

RNA sequencing of peripheral blood in amyotrophic lateral sclerosis reveals distinct molecular subtypes: Considerations for biomarker discovery

AIM

Amyotrophic lateral sclerosis (ALS) is a heterogeneous neurodegenerative disease with limited therapeutic options. A key factor limiting the development of effective therapeutics is the lack of disease biomarkers. We sought to assess whether biomarkers for diagnosis, prognosis or cohort stratification could be identified by RNA sequencing (RNA-seq) of ALS patient peripheral blood.

METHODS

Whole blood RNA-seq data were generated for 96 Australian sporadic ALS (sALS) cases and 48 healthy controls (NCBI GEO accession GSE234297). Differences in sALS-control gene expression, transcript usage and predicted leukocyte proportions were assessed, with pathway analysis used to predict the activity state of biological processes. Weighted Gene Co-expression Network Analysis (WGCNA) and machine learning algorithms were applied to search for diagnostic and prognostic gene expression patterns. Unsupervised clustering analysis was employed to determine whether sALS patient subgroups could be detected.

RESULTS

Two hundred and forty-five differentially expressed genes were identified in sALS patients relative to controls, with enrichment of immune, metabolic and stress-related pathways. sALS patients also demonstrated switches in transcript usage across a small set of genes. We established a classification model that distinguished sALS patients from controls with an accuracy of 78% (sensitivity: 79%, specificity: 75%) using the expression of 20 genes. Clustering analysis identified four patient subgroups with gene expression signatures and immune cell proportions reflective of distinct peripheral effects.

CONCLUSIONS

Our findings suggest that peripheral blood RNA-seq can identify diagnostic biomarkers and distinguish molecular subtypes of sALS patients however, its prognostic value requires further investigation.

Understanding the impact of structural modifications at the NNAT gene's post-translational acetylation site: in silico approach for predicting its drug-interaction role in anorexia nervosa

PURPOSE

Anorexia nervosa (AN) is a neuropsychological public health concern with a socially disabling routine and affects a person's healthy relationship with food. The role of the NNAT (Neuronatin) gene in AN is well established. The impact of mutation at the protein's post-translational modification (PTM) site has been exclusively associated with the worsening of the protein's biochemical dynamics.

METHODS

To understand the relationship between genotype and phenotype, it is essential to investigate the appropriate molecular stability of protein required for proper biological functioning. In this regard, we investigated the PTM-acetylation site of the NNAT gene in terms of 19 other specific amino acid probabilities in place of wild type (WT) through various in silico algorithms. Based on the highest pathogenic impact computed through the consensus classifier tool, we generated 3 residue-specific (K59D, P, W) structurally modified 3D models of NNAT. These models were further tested through the AutoDock Vina tool to compute the molecular drug binding affinities and inhibition constant (Ki) of structural variants and WT 3D models.

RESULTS

With trained in silico machine learning algorithms and consensus classifier; the three structural modifications (K59D, P, W), which were also the most deleterious substitution at the acetylation site of the NNAT gene, showed the highest structural destabilization and decreased molecular flexibility. The validation and quality assessment of the 3D model of these structural modifications and WT were performed. They were further docked with drugs used to manage AN, it was found that the ΔGbind (kcal/mol) values and the inhibition constants (Ki) were relatively lower in structurally modified models as compared to WT.

CONCLUSION

We concluded that any future structural variation(s) at the PTM-acetylation site of the NNAT gene due to possible mutational consequences, will serve as a basis to explore its relationship with the propensity of developing AN.

LEVEL OF EVIDENCE

No level of evidence-open access bioinformatics research.

Genetic insights into obesity: in silico identification of pathogenic SNPs in MBOAT4 gene and their structural molecular dynamics consequences

Membrane Bound O-Acyltransferase Domain-Containing 4 (MBOAT4) protein catalyzes ghrelin acylation, leading to prominent ghrelin activity, hence characterizing its role as an anti-obesity target. We extracted 625 exonic SNPs from the ENSEMBL database and one phenotype-based missense mutation associated with obesity (A46T) from the HGMD (Human Gene Mutation Database). These were differentiated on deleterious missense SNPs of the MBOAT4 gene through MAF (minor allele frequency: <0.01) cut-off criteria in relation to some bioinformatics-based supervised machine learning tools. We found 8 rare-coding and harmful missense SNPs. The consensus classifier (PredictSNP) tool predicted that the SNP (G57S, C: rs561065025) was the most pathogenic. Several trained in silico algorithms have predicted decreased protein stability [ΔΔG (kcal/mol)] function in the presence of these rare-coding pathogenic mutations in the MBOAT4 gene. Then, a stereochemical quality check (i.e. validation and assessment) of the 3D model was performed, followed by a blind cavity docking approach, used to search for druggable cavities and molecular interactions with citrus flavonoids of the Rutaceae family, ranked with energetic estimations. Significant interactions with Phloretin 3',5'-Di-C-Glucoside were also observed at R304, W306, N307, A311, L314 and H338 with (iGEMDOCK: -95.82 kcal/mol and AutoDock: -7.80 kcal/mol). The RMSD values and other variables of MD simulation analyses on this protein further validated its significant interactions with the above flavonoids. The MBOAT4 gene and its molecular interactions could serve as an interventional future anti-obesity target. The current study's findings will benefit future prospects for large population-based studies and drug development, particularly for generating personalized medicine.Communicated by Ramaswamy H. Sarma.

Transcript-Level In Silico Analysis of Alzheimer's Disease-Related Gene Biomarkers and Their Evaluation with Bioactive Flavonoids to Explore Therapeutic Interactions

Alzheimer's disease (AD) is a progressive brain disorder that can significantly affect the quality of life. We used a variety of in silico tools to investigate the transcript-level mutational impact of exonic missense rare variations (single nucleotide polymorphisms, SNPs) on protein function and to identify potential druggable protein cavities that correspond to potential therapeutic targets for the management of AD. According to the NIA-AA (National Institute on Aging-Alzheimer's Association) framework, we selected three AD biomarker genes (APP, NEFL, and MAPT). We systematically screened transcript-level exonic rare SNPs from these genes with a minor allele frequency of 1% in 1KGD (1000 Genomes Project Database) and gnomAD (Genome Aggregation Database). With downstream functional effect predictions, a single variation (rs182024939: K > N) of the MAPT gene with nine transcript SNPs was identified as the most pathogenic variation from the large dataset of mutations. The machine learning consensus classifier predictor categorized these transcript-level SNPs as the most deleterious variations, resulting in a large decrease in protein structural stability (ΔΔG kcal/mol). The bioactive flavonoid library was screened for drug-likeness and toxicity risk. Virtual screening of eligible flavonoids was performed using the MAPT protein. Identification of druggable protein-binding cavities showed VAL305, GLU655, and LYS657 as consensus-interacting residues present in the MAPT-docked top-ranked flavonoid compounds. The MM/PB(GB)SA analysis indicated hesperetin (-5.64 kcal/mol), eriodictyol (-5.63 kcal/mol), and sakuranetin (-5.60 kcal/mol) as the best docked flavonoids with the near-native binding pose. The findings of this study provide important insights into the potential of hesperetin as a promising flavonoid that can be utilized for further rational drug design and lead optimization to open new gateways in the field of AD therapeutics.

Most monogenic disorders are caused by mutations altering protein folding free energy

Revealing the molecular effect that pathogenic missense mutations cause on the corresponding protein is crucial for developing therapeutic solutions. This is especially important for monogenic diseases since, for most of them, there is no treatment available, while typically, the treatment should be provided in the early development stages. This requires fast, targeted drug development at a low cost. Here, we report a database of monogenic disorders (MOGEDO), which includes 768 proteins, the corresponding 2559 pathogenic and 1763 benign mutations, along with the functional classification of the corresponding proteins. Using the database and various computational tools that predict folding free energy change (ΔΔG), we demonstrate that, on average, 70% of pathogenic cases result in decreased protein stability. Such a large fraction indicates that one should aim at in-silico screening for small molecules stabilizing the structure of the mutant protein. We emphasize that knowledge of ΔΔG is essential because one wants to develop stabilizers that compensate for ΔΔG but not to make protein over-stable since over-stable protein may be dysfunctional. We demonstrate that using ΔΔG and predicted solvent exposure of the mutation site; one can develop a predictive method that distinguishes pathogenic from benign mutation with a success rate even better than some of the leading pathogenicity predictors. Furthermore, hydrophobic-hydrophobic mutations have stronger correlations between folding free energy change and pathogenicity compared with others. Also, mutations involving Cys, Gly, Arg, Trp and Tyr amino acids being replaced by any other amino acid are more likely to be pathogenic. To facilitate further detection of pathogenic mutations, the wild type of amino acids in the 768 proteins mentioned above was mutated to other 19 residues (14,847,817 mutations), and the ΔΔG was calculated with SAAFEC-SEQ, and 5,506,051 mutations were predicted to be pathogenic.

Naive Prediction of Protein Backbone Phi and Psi Dihedral Angles Using Deep Learning

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

A novel vaccine candidate against A. baumannii based on a new OmpW family protein (OmpW2); structural characterization, antigenicity and epitope investigation, and in-vivo analysis

A. baumannii is an MDR pathogen whose SARS-CoV-2 has recently increased its mortality rate in hospitalized patients. So, the virulence factors investigation and design of a vaccine against this bacterium seem to be critical. In this regard, the OmpW2 protein was structurally characterized by this study, and its B-T cell epitopes were mapped by bioinformatic tools. In-vivo analyses were employed to verify the immunogenicity of this protein and its selected epitopes. The results indicated that OmpW2 is a conserved virulent antigen, not toxic for the host, and not similar to the human or mouse proteome. A putative interaction between OmpW2 and a Fe-S-cluster redox enzyme was detected. Based on the results, OmpW2 belongs to the OmpW superfamily and eight beta sheets have been predicted in its tight beta-barrel structure. Several exposed epitopes were detected in the OmpW2 sequence and structure, and a sub-unit potential vaccine was generated based on the epitopes. The ELISA results indicated that after the second booster vaccination of BALB/c mice with the whole OmpW2 protein or its sub-unit fragment, the IgG titer significantly raised (p < 0.05). The mortality rate and the bacterial burden in the lung, liver, kidney, and spleen in both passive and active immunized mice were significantly decreased (p ≤ 0.001). In-vivo experiments confirmed that the OmpW2 whole protein and its sub-unit fragment induce the host immune system and can be applied to design a commercial vaccine or diagnostic kit.

Identification of the effects of pathogenic genetic variations of human CYP2C9 and CYP2D6: an in silico approach

Genetic variations in the human cytochrome P450 family 2 subfamily C member 9 (CYP2C9) and cytochrome P450 family 2 subfamily D member 6 (CYP2D6) genes may affect drug metabolism and lead to alterations in phenotypes. Genetic variations are associated with toxicity, adverse drug reactions, inefficient treatment. Various in silico tools were combined to investigate the deleterious effects of missense non-synonymous single nucleotide polymorphisms (nsSNPs) of the human CYP2C9 and CYP2D6. The structural and functional effects of the high-risk non-synonymous SNPs in the human CYP2C9 and CYP2D6 were predicted by numerous computational mutation analysis methods. Out of 24 pathogenic missense SNPs in the CYP2C9, 22 nsSNPs had a decreasing effect on protein stability and 13 SNPs were showed to be located at conserved positions. Out of 27 high-risk deleterious non-synonymous SNPs in the human CYP2D6, 21 SNPs decreased protein stability and 16 nsSNPs were predicted to be positioned at conserved regions. Our present study suggests that the identified functional SNPs may affect drug metabolism associated with CYP2C9 and CYP2D6 enzymes.

Deciphering Protein Secretion from the Brain to Cerebrospinal Fluid for Biomarker Discovery

Cerebrospinal fluid (CSF) is an essential matrix for the discovery of neurological disease biomarkers. However, the high dynamic range of protein concentrations in CSF hinders the detection of the least abundant protein biomarkers by untargeted mass spectrometry. It is thus beneficial to gain a deeper understanding of the secretion processes within the brain. Here, we aim to explore if and how the secretion of brain proteins to the CSF can be predicted. By combining a curated CSF proteome and the brain elevated proteome of the Human Protein Atlas, brain proteins were classified as CSF or non-CSF secreted. A machine learning model was trained on a range of sequence-based features to differentiate between CSF and non-CSF groups and effectively predict the brain origin of proteins. The classification model achieves an area under the curve of 0.89 if using high confidence CSF proteins. The most important prediction features include the subcellular localization, signal peptides, and transmembrane regions. The classifier generalized well to the larger brain detected proteome and is able to correctly predict novel CSF proteins identified by affinity proteomics. In addition to elucidating the underlying mechanisms of protein secretion, the trained classification model can support biomarker candidate selection.

Comprehensive analysis predicting effects of deleterious SNPs of human progesterone receptor gene on its structure and functions: a computational approach

Progesterone receptor plays a crucial role in the development of the mammary gland and breast cancer. Single nucleotide polymorphisms (SNPs) within its gene, PGR, are associated with the risk of miscarriages and preterm birth as well as many cancers across different populations. The main aim of this work is to investigate the most deleterious SNPs in the PGR gene to identify potential biomarkers for various disease susceptibility and treatments. Both sequence and structure-based computational approaches were adopted and in total 11 nsSNPs have been filtered out of 674 nsSNPs along with seven non-coding SNPs. R740Q, I744T and D746E belonged to a mutation cluster. R740Q, D746E along with S865L altered H-bond interactions within the receptor. The same mutations have been found to be associated with several cancers including uterine and breast cancer among others. It is, therefore, possible that the high-risk SNPs associated with cancers may exert their effect by causing changes in the protein structure, particularly in its bonding patterns, and thus affecting its function. In addition, seven non-coding SNPs that were located in the UTR region created a new miRNA site while three SNPs disrupted a conserved miRNA site. These high-risk SNPs can play an instrumental role in generating a dataset of the PGR gene's SNPs. Thus, the present study may pave the way to design and develop novel therapeutics for overcoming the challenges associated with certain cancers and pregnancy that result from a change in the protein structure and function due to the SNP mutations in the PGR gene.Communicated by Ramaswamy H. Sarma.

DichroIDP: a method for analyses of intrinsically disordered proteins using circular dichroism spectroscopy

Intrinsically disordered proteins (IDPs) are comprised of significant numbers of residues that form neither helix, sheet, nor any other canonical type of secondary structure. They play important roles in a broad range of biological processes, such as molecular recognition and signalling, largely due to their chameleon-like ability to change structure from unordered when free in solution to ordered when bound to partner molecules. Circular dichroism (CD) spectroscopy is a widely-used method for characterising protein secondary structures, but analyses of IDPs using CD spectroscopy have suffered because the methods and reference datasets used for the empirical determination of secondary structures do not contain adequate representations of unordered structures. This work describes the creation, validation and testing of a standalone Windows-based application, DichroIDP, and a new reference dataset, IDP175, which is suitable for analyses of proteins containing significant amounts of disordered structure. DichroIDP enables secondary structure determinations of IDPs and proteins containing intrinsically disordered regions.

Systematic analysis of the impact of phosphorylation and O-GlcNAcylation on protein subcellular localization

The subcellular localization of proteins is critical for their functions in eukaryotic cells and is tightly correlated with protein modifications. Here, we comprehensively investigate the nuclear-cytoplasmic distributions of the phosphorylated, O-GlcNAcylated, and non-modified forms of proteins to dissect the correlation between protein distribution and modifications. Phosphorylated and O-GlcNAcylated proteins have overall higher nuclear distributions than non-modified ones. Different distributions among the phosphorylated, O-GlcNAcylated, and non-modified forms of proteins are associated with protein size, structure, and function, as well as local environment and adjacent residues around modification sites. Moreover, we perform site-mutagenesis experiments using phosphomimetic and phospho-null mutants of two proteins to validate the proteomic results. Additionally, the effects of the OGT/OGA inhibition on glycoprotein distribution are systematically investigated, and the distribution changes of glycoproteins are related to their abundance changes under the inhibitions. Systematic investigation of the relationship between protein modification and localization advances our understanding of protein functions.

In silico analysis of sporozoite surface antigen 1 of Theileria annulata (TaSPAG1) for multi-epitope vaccine design against theileriosis

Tropical theileriosis is a protozoan infection caused by Theileria annulata, which significantly affects cattle worldwide. This study was aimed to analyze the TaSPAG1 protein and design a novel multi-epitope vaccine candidate. Online tools were employed for the prediction of Physico-chemical properties, antigenicity, allergenicity, solubility, transmembrane domains and signal peptide, posttranslational modification (PTM) sites, secondary and tertiary structures as well as intrinsically disordered regions, followed by identification and screening of potential linear and conformational B-cell epitopes and those peptides having affinity to bind bovine major histocompatibility complex class I (MHC-I) molecules. Next, a multi-epitope vaccine construct was designed and analyzed. This 907-residue protein was hydrophilic (GRAVY: -0.399) and acidic (pI: 5.04) in nature, with high thermotolerance (aliphatic: 71.27). Also, 5 linear and 12 conformational B-cell epitopes along with 8 CTL epitopes were predicted for TaSPAG1. The 355-residue vaccine candidate had a MW of about 35 kDa and it was antigenic, non-allergenic, soluble and stable, which was successfully interacted with cattle MHC-I molecule and finally cloned into the pET28a(+) vector. Further wet studies are required to assess the vaccine efficacy in cattle.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-023-00153-5.

Pooled tagging and hydrophobic targeting of endogenous proteins for unbiased mapping of unfolded protein responses

System-level understanding of proteome organization and function requires methods for direct visualization and manipulation of proteins at scale. We developed an approach enabled by high-throughput gene tagging for the generation and analysis of complex cell pools with endogenously tagged proteins. Proteins are tagged with HaloTag to enable visualization or direct perturbation. Fluorescent labeling followed by in situ sequencing and deep learning-based image analysis identifies the localization pattern of each tag, providing a bird's-eye-view of cellular organization. Next, we use a hydrophobic HaloTag ligand to misfold tagged proteins, inducing spatially restricted proteotoxic stress that is read out by single cell RNA sequencing. By integrating optical and perturbation data, we map compartment-specific responses to protein misfolding, revealing inter-compartment organization and direct crosstalk, and assigning proteostasis functions to uncharacterized genes. Altogether, we present a powerful and efficient method for large-scale studies of proteome dynamics, function, and homeostasis.

Deep structured learning for variant prioritization in Mendelian diseases

Effective computer-aided or automated variant evaluations for monogenic diseases will expedite clinical diagnostic and research efforts of known and novel disease-causing genes. Here we introduce MAVERICK: a Mendelian Approach to Variant Effect pRedICtion built in Keras. MAVERICK is an ensemble of transformer-based neural networks that can classify a wide range of protein-altering single nucleotide variants (SNVs) and indels and assesses whether a variant would be pathogenic in the context of dominant or recessive inheritance. We demonstrate that MAVERICK outperforms all other major programs that assess pathogenicity in a Mendelian context. In a cohort of 644 previously solved patients with Mendelian diseases, MAVERICK ranks the causative pathogenic variant within the top five variants in over 95% of cases. Seventy-six percent of cases were solved by the top-ranked variant. MAVERICK ranks the causative pathogenic variant in hitherto novel disease genes within the first five candidate variants in 70% of cases. MAVERICK has already facilitated the identification of a novel disease gene causing a degenerative motor neuron disease. These results represent a significant step towards automated identification of causal variants in patients with Mendelian diseases.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

ATCUN-like Copper Site in βB2-Crystallin Plays a Protective Role in Cataract-Associated Aggregation

Cataract is the leading cause of blindness worldwide, and it is caused by crystallin damage and aggregation. Senile cataractous lenses have relatively high levels of metals, while some metal ions can directly induce the aggregation of human γ-crystallins. Here, we evaluated the impact of divalent metal ions in the aggregation of human βB2-crystallin, one of the most abundant crystallins in the lens. Turbidity assays showed that Pb2+, Hg2+, Cu2+, and Zn2+ ions induce the aggregation of βB2-crystallin. Metal-induced aggregation is partially reverted by a chelating agent, indicating the formation of metal-bridged species. Our study focused on the mechanism of copper-induced aggregation of βB2-crystallin, finding that it involves metal-bridging, disulfide-bridging, and loss of protein stability. Circular dichroism and electron paramagnetic resonance (EPR) revealed the presence of at least three Cu2+ binding sites in βB2-crystallin, one of them with spectroscopic features typical for Cu2+ bound to an amino-terminal copper and nickel (ATCUN) binding motif, which is found in Cu transport proteins. The ATCUN-like Cu binding site is located at the unstructured N-terminus of βB2-crystallin, and it could be modeled by a peptide with the first six residues in the protein sequence (NH2-ASDHQF-). Isothermal titration calorimetry indicates a nanomolar Cu2+ binding affinity for the ATCUN-like site. An N-truncated form of βB2-crystallin is more susceptible to Cu-induced aggregation and is less thermally stable, indicating a protective role for the ATCUN-like site. EPR and X-ray absorption spectroscopy studies reveal the presence of a copper redox active site in βB2-crystallin that is associated with metal-induced aggregation and formation of disulfide-bridged oligomers. Our study demonstrates metal-induced aggregation of βB2-crystallin and the presence of putative copper binding sites in the protein. Whether the copper-transport ATCUN-like site in βB2-crystallin plays a functional/protective role or constitutes a vestige from its evolution as a lens structural protein remains to be elucidated.

Impact of highly deleterious non-synonymous polymorphisms on GRIN2A protein's structure and function

GRIN2A is a gene that encodes NMDA receptors found in the central nervous system and plays a pivotal role in excitatory synaptic transmission, plasticity and excitotoxicity in the mammalian central nervous system. Changes in this gene have been associated with a spectrum of neurodevelopmental disorders such as epilepsy. Previous studies on GRIN2A suggest that non-synonymous single nucleotide polymorphisms (nsSNPs) can alter the protein's structure and function. To gain a better understanding of the impact of potentially deleterious variants of GRIN2A, a range of bioinformatics tools were employed in this study. Out of 1320 nsSNPs retrieved from the NCBI database, initially 16 were predicted as deleterious by 9 tools. Further assessment of their domain association, conservation profile, homology models, interatomic interaction, and Molecular Dynamic Simulation revealed that the variant I463S is likely to be the most deleterious for the structure and function of the protein. Despite the limitations of computational algorithms, our analyses have provided insights that can be a valuable resource for further in vitro and in vivo research on GRIN2A-associated diseases.

Improved structure-related prediction for insufficient homologous proteins using MSA enhancement and pre-trained language model

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

CONTACT

guofei@csu.edu.cn, jj.tang@siat.ac.cn.

Deep Proteogenomics of a Photosynthetic Cyanobacterium

Cyanobacteria, the evolutionary ancestors of plant chloroplasts, contribute substantially to the Earth's biogeochemical cycles and are of great interest for a sustainable economy. Knowledge of protein expression is the key to understanding cyanobacterial metabolism; however, proteome studies in cyanobacteria are limited and cover only a fraction of the theoretical proteome. Here, we performed a comprehensive proteogenomic analysis of the model cyanobacterium Synechocystis sp. PCC 6803 to characterize the expressed (phospho)proteome, re-annotate known and discover novel open reading frames (ORFs). By mapping extensive shotgun mass spectrometry proteomics data onto a six-frame translation of the Synechocystis genome, we refined the genomic annotation of 64 ORFs, including eight completely novel ORFs. Our study presents the largest reported (phospho)proteome dataset for a unicellular cyanobacterium, covering the expression of about 80% of the theoretical proteome under various cultivation conditions, such as nitrogen or carbon limitation. We report 568 phosphorylated S/T/Y sites that are present on numerous regulatory proteins, including the transcriptional regulators cyAbrB1 and cyAbrB2. We also catalogue the proteins that have never been detected under laboratory conditions and found that a large portion of them is plasmid-encoded. This dataset will serve as a resource, providing dedicated information on growth condition-dependent protein expression and phosphorylation.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

Prediction of human protein interactome of dengue virus non-structural protein 5 (NS5) and its downstream immunological implications

The non-structural protein 5 (NS5) is the most conserved protein among flaviviruses, a family that includes the dengue virus. It functions both as an RNA-dependent RNA polymerase and an RNA-methyltransferase and is therefore essential for the replication of viral RNA. The discovery that dengue virus NS5 protein (DENV-NS5) can also localize to the nucleus has resulted in renewed interest in its potential roles at the host-virus interface. In this study, we have used two complementary computational approaches in parallel - one based on linear motifs (ELM) and another based on tertiary structure of the protein (DALI) - to predict the host proteins that DENV-NS5 might interact with. Of the 42 human proteins predicted by both these methods, 34 are novel. Pathway analysis of these 42 human proteins shows that they are involved in key host cellular processes related to cell cycle regulation, proliferation, protein degradation, apoptosis, and immune responses. A focused analysis of transcription factors that directly interact with the predicted DENV-NS5 interacting proteins was performed, followed by the identification of downstream genes that are differentially expressed after dengue infection using previously published RNA-seq data. Our study provides unique insights into the DENV-NS5 interaction network and delineates mechanisms whereby DENV-NS5 could impact the host-virus interface. The novel interactors identified in this study could be potentially targeted by NS5 to modulate the host cellular environment in general, and the immune response in particular, thereby extending the role of DENV-NS5 beyond its known enzymatic functions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03569-0.

AttSec: protein secondary structure prediction by capturing local patterns from attention map

BACKGROUND

Protein secondary structures that link simple 1D sequences to complex 3D structures can be used as good features for describing the local properties of protein, but also can serve as key features for predicting the complex 3D structures of protein. Thus, it is very important to accurately predict the secondary structure of the protein, which contains a local structural property assigned by the pattern of hydrogen bonds formed between amino acids. In this study, we accurately predict protein secondary structure by capturing the local patterns of protein. For this objective, we present a novel prediction model, AttSec, based on transformer architecture. In particular, AttSec extracts self-attention maps corresponding to pairwise features between amino acid embeddings and passes them through 2D convolution blocks to capture local patterns. In addition, instead of using additional evolutionary information, it uses protein embedding as an input, which is generated by a language model.

RESULTS

For the ProteinNet DSSP8 dataset, our model showed 11.8% better performance on the entire evaluation datasets compared with other no-evolutionary-information-based models. For the NetSurfP-2.0 DSSP8 dataset, it showed 1.2% better performance on average. There was an average performance improvement of 9.0% for the ProteinNet DSSP3 dataset and an average of 0.7% for the NetSurfP-2.0 DSSP3 dataset.

CONCLUSION

We accurately predict protein secondary structure by capturing the local patterns of protein. For this objective, we present a novel prediction model, AttSec, based on transformer architecture. Although there was no dramatic accuracy improvement compared with other models, the improvement on DSSP8 was greater than that on DSSP3. This result implies that using our proposed pairwise feature could have a remarkable effect for several challenging tasks that require finely subdivided classification. Github package URL is https://github.com/youjin-DDAI/AttSec .

Improvement in salt-tolerance of Aspergillus oryzae γ-glutamyl transpeptidase via protein chimerization with Aspergillus sydowii homolog

γ-Glutamyl transpeptidase is one of the key enzymes involved in glutamate production during high-salt fermentation of soy sauce and miso by koji mold, Aspergillus oryzae. However, the activity of γ-glutamyl transpeptidase from A. oryzae (AOggtA) is markedly reduced in the presence of NaCl, thus classifying it as a non-salt-tolerant enzyme. In contrast, the homologous protein from the xerophilic mold, A. sydowii (ASggtA) maintains its activity under high-salt conditions. Therefore, in this study, a chimeric enzyme, ASAOggtA, was designed and engineered to improve salt-tolerance in AOggtA by swapping the N-terminal region, based on sequence and structure comparisons between salt-tolerant ASggtA and non-salt-tolerant AOggtA. The parental AOggtA and ASggtA and their chimera, ASAOggtA, were heterologously expressed in A. oryzae and purified. The chimeric enzyme inherited the superior activity and stability from each of the two parent enzymes. ASAOggtA showed > 2-fold greater tolerance than AOggtA in the presence of 18% NaCl. In addition, the chimera showed a broader range of pH stability and greater thermostability than ASggtA. AOggtA and ASAOggtA were sy over the range pH 3.0 to pH 10.5. Thermal stability was found to be in the order AOggtA (57.5 °C, t_1/2 = 32.5 min) > ASAOggtA (55 °C, t_1/2 = 20.5 min) > ASggtA (50 °C, t_1/2 = 12.5 min). The catalytic and structural characteristics indicated that non-salt-tolerant AOggtA would not undergo irreversible structural changes in the presence of NaCl, but rather a temporary conformational change, which might result in reducing the substrate binding and catalytic activity, on the basis of kinetic properties. In addition, the chimeric enzyme showed hydrolytic activity toward L-glutamine that was as high as that of AOggtA. The newly-designed chimeric ASAOggtA might have potential applications in high-salt fermentation, such as miso and shoyu, to increase the content of the umami-flavor amino acid, L-glutamate.

Zasp52 strengthens whole embryo tissue integrity through supracellular actomyosin networks

During morphogenesis, large-scale changes of tissue primordia are coordinated across an embryo. In Drosophila, several tissue primordia and embryonic regions are bordered or encircled by supracellular actomyosin cables, junctional actomyosin enrichments networked between many neighbouring cells. We show that the single Drosophila Alp/Enigma-family protein Zasp52, which is most prominently found in Z-discs of muscles, is a component of many supracellular actomyosin structures during embryogenesis, including the ventral midline and the boundary of the salivary gland placode. We reveal that Zasp52 contains within its central coiled-coil region a type of actin-binding motif usually found in CapZbeta proteins, and this domain displays actin-binding activity. Using endogenously-tagged lines, we identify that Zasp52 interacts with junctional components, including APC2, Polychaetoid and Sidekick, and actomyosin regulators. Analysis of zasp52 mutant embryos reveals that the severity of the embryonic defects observed scales inversely with the amount of functional protein left. Large tissue deformations occur where actomyosin cables are found during embryogenesis, and in vivo and in silico analyses suggest a model whereby supracellular Zasp52-containing cables aid to insulate morphogenetic changes from one another.

Align-gram: Rethinking the Skip-gram Model for Protein Sequence Analysis

The inception of next generations sequencing technologies have exponentially increased the volume of biological sequence data. Protein sequences, being quoted as the 'language of life', has been analyzed for a multitude of applications and inferences. Owing to the rapid development of deep learning, in recent years there have been a number of breakthroughs in the domain of Natural Language Processing. Since these methods are capable of performing different tasks when trained with a sufficient amount of data, off-the-shelf models are used to perform various biological applications. In this study, we investigated the applicability of the popular Skip-gram model for protein sequence analysis and made an attempt to incorporate some biological insights into it. We propose a novel k-mer embedding scheme, Align-gram, which is capable of mapping the similar k-mers close to each other in a vector space. Furthermore, we experiment with other sequence-based protein representations and observe that the embeddings derived from Align-gram aids modeling and training deep learning models better. Our experiments with a simple baseline LSTM model and a much complex CNN model of DeepGoPlus shows the potential of Align-gram in performing different types of deep learning applications for protein sequence analysis.

Dps-dependent in vivo mutation enhances long-term host adaptation in Vibrio cholerae

As one of the most successful pathogenic organisms, Vibrio cholerae (V. cholerae) has evolved sophisticated regulatory mechanisms to overcome host stress. During long-term colonization by V. cholerae in adult mice, many spontaneous nonmotile mutants (approximately 10% at the fifth day post-infection) were identified. These mutations occurred primarily in conserved regions of the flagellar regulator genes flrA, flrC, and rpoN, as shown by Sanger and next-generation sequencing, and significantly increased fitness during colonization in adult mice. Intriguingly, instead of key genes in DNA repair systems (mutS, nfo, xthA, uvrA) or ROS and RNS scavenging systems (katG, prxA, hmpA), which are generally thought to be associated with bacterial mutagenesis, we found that deletion of the cyclin gene dps significantly increased the mutation rate (up to 53% at the fifth day post-infection) in V. cholerae. We further determined that the dpsD65A and dpsF46E point mutants showed a similar mutagenesis profile as the Δdps mutant during long-term colonization in mice, which strongly indicated that the antioxidative function of Dps directly contributes to the development of V. cholerae nonmotile mutants. Methionine metabolism pathway may be one of the mechanism for ΔflrA, ΔflrC and ΔrpoN mutant increased colonization in adult mice. Our results revealed a new phenotype in which V. cholerae fitness increases in the host gut via spontaneous production nonmotile mutants regulated by cyclin Dps, which may represent a novel adaptation strategy for directed evolution of pathogens in the host.

High risk genetic variants of human insulin receptor substrate 1(IRS1) infer structural instability and functional interference

Insulin receptor substrate 1(IRS1) is a signaling adapter protein encoded by the IRS1 gene. This protein delivers signals from insulin and insulin-like growth factor-1(IGF-1) receptors to the phosphatidylinositol 3-kinases (P13K)/protein kinase B (Akt) and Extracellular signal-regulated kinases (Erk) - Mitogen-activated protein (MAP) kinase pathways, which regulate particular cellular processes. Mutations in this gene have been linked to type 2 diabetes mellitus, a heightened risk of insulin resistance, and an increased likelihood of developing a number of different malignancies. The structure and function of IRS1 could be severely compromised as a result of single nucleotide polymorphism (SNP) type genetic variants. In this study, we focused on identification of the most harmful non-synonymous SNPs (nsSNPs) of the IRS1 gene as well as prediction of their structural and functional consequences. Six different algorithms made the initial prediction that 59 of the 1142 IRS1 nsSNPs would have a negative impact on the protein structure. In-depth evaluations detected 26 nsSNPs located inside the functional domains of IRS1. Following that, 16 nsSNPs were identified as more harmful based on conservation profile, hydrophobic interaction, surface accessibility, homology modelling, and inter-atomic interactions. Following an in-depth analysis of protein stability, M249T (rs373826433), I223T (rs1939785175) and V204G (rs1574667052) were identified as three most deleterious SNPs and were subjected to molecular dynamics simulation for further insights. These findings will help us understand the implications for disease susceptibility, cancer progression, and the efficacy of therapeutic development against IRS1 gene mutants.Communicated by Ramaswamy H. Sarma.

Unsupervised cross-domain translation via deep learning and adversarial attention neural networks and application to music-inspired protein designs

Taking inspiration from nature about how to design materials has been a fruitful approach, used by humans for millennia. In this paper we report a method that allows us to discover how patterns in disparate domains can be reversibly related using a computationally rigorous approach, the AttentionCrossTranslation model. The algorithm discovers cycle- and self-consistent relationships and offers a bidirectional translation of information across disparate knowledge domains. The approach is validated with a set of known translation problems, and then used to discover a mapping between musical data-based on the corpus of note sequences in J.S. Bach's Goldberg Variations created in 1741-and protein sequence data-information sampled more recently. Using protein folding algorithms, 3D structures of the predicted protein sequences are generated, and their stability is validated using explicit solvent molecular dynamics. Musical scores generated from protein sequences are sonified and rendered into audible sound.

IGPRED-MultiTask: A Deep Learning Model to Predict Protein Secondary Structure, Torsion Angles and Solvent Accessibility

Protein secondary structure, solvent accessibility and torsion angle predictions are preliminary steps to predict 3D structure of a protein. Deep learning approaches have achieved significant improvements in predicting various features of protein structure. In this study, IGPRED-Multitask, a deep learning model with multi task learning architecture based on deep inception network, graph convolutional network and a bidirectional long short-term memory is proposed. Moreover, hyper-parameters of the model are fine-tuned using Bayesian optimization, which is faster and more effective than grid search. The same benchmark test data sets as in the OPUS-TASS paper including TEST2016, TEST2018, CASP12, CASP13, CASPFM, HARD68, CAMEO93, CAMEO93_HARD, as well as the train and validation sets, are used for fair comparison with the literature. Statistically significant improvements are observed in secondary structure prediction on 4 datasets, in phi angle prediction on 2 datasets and in psi angel prediction on 3 datasets compared to the state-of-the-art methods. For solvent accessibility prediction, TEST2016 and TEST2018 datasets are used only to assess the performance of the proposed model.

Rapid adaptations of Legionella pneumophila to the human host

Legionella pneumophila are host-adapted bacteria that infect and reproduce primarily in amoeboid protists. Using similar infection mechanisms, they infect human macrophages, and cause Legionnaires' disease, an atypical pneumonia, and the milder Pontiac fever. We hypothesized that, despite the similarities in infection mechanisms, the hosts are different enough that there exist high-selective value mutations that would dramatically increase the fitness of Legionella inside the human host. By comparing a large number of isolates from independent infections, we identified two genes, mutated in three unrelated patients, despite the short duration of the incubation period (2-14 days). One is a gene coding for an outer membrane protein (OMP) belonging to the OmpP1/FadL family. The other is a gene coding for an EAL-domain-containing protein involved in cyclic-di-GMP regulation, which in turn modulates flagellar activity. The clinical strain, carrying the mutated EAL-domain-containing homologue, grows faster in macrophages than the wild-type strain, and thus appears to be better adapted to the human host. As human-to-human transmission is very rare, fixation of these mutations into the population and spread into the environment is unlikely. Therefore, parallel evolution - here mutations in the same genes observed in independent human infections - could point to adaptations to the accidental human host. These results suggest that despite the ability of L. pneumophila to infect, replicate in and exit from macrophages, its human-specific adaptations are unlikely to be fixed in the population.