Evolutionary Information Research Articles

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows

BackgroundOver the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM.ResultsIn general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits’ genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets’ high effect.ConclusionOur results confirm the connection of milk production and health traits in dairy cattle via the animals’ metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits’ connection in more detail.

Prediction of protein-ATP binding residues using multi-view feature learning via contextual-based co-attention network

Accurately predicting protein-ATP binding residues is critical for protein function annotation and drug discovery. Computational methods dedicated to the prediction of binding residues based on protein sequence information have exhibited notable advancements in predictive accuracy. Nevertheless, these methods continue to grapple with several formidable challenges, including limited means of extracting more discriminative features and inadequate algorithms for integrating protein and residue information. To address the problems, we propose ATP-Deep, a novel protein-ATP binding residues predictor. ATP-Deep harnesses the capabilities of unsupervised pre-trained language models and incorporates domain-specific evolutionary context information from homologous sequences. It further refines the embedding at the residue level through integration with corresponding protein-level information and employs a contextual-based co-attention mechanism to adeptly fuse multiple sources of features. The performance evaluation results on the benchmark datasets reveal that ATP-Deep achieves an AUC of 0.954 and 0.951, respectively, surpassing the performance of the state-of-the-art model. These findings underscore the effectiveness of assimilating protein-level information and deploying a contextual-based co-attention mechanism grounded in context to bolster the prediction performance of protein-ATP binding residues.

DisPredict3.0: Prediction of intrinsically disordered regions/proteins using protein language model

Intrinsically disordered proteins (IDPs) or protein regions (IDRs) do not have a stable three-dimensional structure, even though they exhibit important biological functions. They are structurally and functionally very different from ordered proteins and can cause many critical diseases. Accurate identification of disordered proteins/regions significantly impacts fields such as drug design, protein engineering, protein design, and related research. However, experimental identification of IDRs is complex and time-consuming, necessitating the development of an accurate and efficient computational method. The recent development of deep learning methods for protein language models shows the ability to learn evolutionary information from billions of protein sequences. This motivates us to develop a computational method, named DisPredict3.0, to predict proteins’ disordered regions (IDRs) using evolutionary information from a protein language model. Compared to the state-of-the-art method in the CAID (2018) assessment, DisPredict3.0 has an improvement of 2.51 %, 16.13 %, 17.98 %, and 11.94 % in terms of AUC, F1-score, MCC, and kappa, respectively. In addition, in the CAID-2 assessment (2022), DisPredict3.0 shows promising results and is ranked first for disorder residue prediction on the Disorder-NOX dataset. The DisPredict3.0 webserver is available at https://bmll.cs.uno.edu.

PlantFUNCO: Integrative Functional Genomics Database Reveals Clues into Duplicates Divergence Evolution.

Evolutionary epigenomics and, more generally, evolutionary functional genomics, are emerging fields that study how non-DNA-encoded alterations in gene expression regulation are an important form of plasticity and adaptation. Previous evidence analyzing plants' comparative functional genomics has mostly focused on comparing same assay-matched experiments, missing the power of heterogeneous datasets for conservation inference. To fill this gap, we developed PlantFUN(ctional)CO(nservation) database, which is constituted by several tools and two main resources: interspecies chromatin states and functional genomics conservation scores, presented and analyzed in this work for three well-established plant models (Arabidopsis thaliana, Oryza sativa, and Zea mays). Overall, PlantFUNCO elucidated evolutionary information in terms of cross-species functional agreement. Therefore, providing a new complementary comparative-genomics source for assessing evolutionary studies. To illustrate the potential applications of this database, we replicated two previously published models predicting genetic redundancy in A. thaliana and found that chromatin states are a determinant of paralogs degree of functional divergence. These predictions were validated based on the phenotypes of mitochondrial alternative oxidase knockout mutants under two different stressors. Taking all the above into account, PlantFUNCO aim to leverage data diversity and extrapolate molecular mechanisms findings from different model organisms to determine the extent of functional conservation, thus, deepening our understanding of how plants epigenome and functional noncoding genome have evolved. PlantFUNCO is available at https://rocesv.github.io/PlantFUNCO.

SME-MFP: A novel spatiotemporal neural network with multiangle initialization embedding toward multifunctional peptides prediction

As a promising alternative to conventional antibiotic drugs in the biomedical field, functional peptide has been widely used in disease treatment owing to its low toxicity, high absorption rate, and biological activity. Recently, several machine learning methods have been developed for functional peptide prediction. However, the main research heavily relies on statistical features and few consider multifunctional peptide identification. So, we propose SME-MFP, a novel predictor in the imbalanced multi-label functional peptide datasets. First, we employ physicochemical and evolutionary information to represent the peptide sequence’s initialization features from multiple perspectives. Second, the features are fused and then put into spatial feature extractors, where the residual connection and multiscale convolutional neural network extract more discriminative features of different lengths’ peptide sequences. Besides, we also design AFT-based temporal feature extractors to fully capture the global interactions of the sequences. Finally, devising a new loss to replace the traditional cross entropy loss to settle the class imbalance problems. The results show that our framework not only enhances the model's ability to capture sequence features effectively, but also accuracy improves by 3.89% over existing methods on public peptide datasets.

Quantifying the effects of exceptional fossil preservation on the global availability of phylogenetic data in deep time.

Fossil deposits with exceptional preservation ("lagerstätten") provide important details not typically preserved in the fossil record, such that they hold an outsized influence on our understanding of biodiversity and evolution. In particular, the potential bias imparted by this so-called "lagerstätten effect" remains a critical, but underexplored aspect of reconstructing evolutionary relationships. Here, we quantify the amount of phylogenetic information available in the global fossil records of 1,327 species of non-avian theropod dinosaurs, Mesozoic birds, and fossil squamates (e.g., lizards, snakes, mosasaurs), and then compare the influence of lagerstätten deposits on phylogenetic information content and taxon selection in phylogenetic analyses to other fossil-bearing deposits. We find that groups that preserve a high amount of phylogenetic information in their global fossil record (e.g., non-avian theropods) are less vulnerable to a "lagerstätten effect" that leads to disproportionate representation of fossil taxa from one geologic unit in an evolutionary tree. Additionally, for each taxonomic group, we find comparable amounts of phylogenetic information in lagerstätten deposits, even though corresponding morphological character datasets vary greatly. Finally, we unexpectedly find that ancient sand dune deposits of the Late Cretaceous Gobi Desert of Mongolia and China exert an anomalously large influence on the phylogenetic information available in the squamate fossil record, suggesting a "lagerstätten effect" can be present in units not traditionally considered lagerstätten. These results offer a phylogenetics-based lens through which to examine the effects of exceptional fossil preservation on biological patterns through time and space, and invites further quantification of evolutionary information in the rock record.

Protein language models meet reduced amino acid alphabets.

Protein language models (PLMs), which borrowed ideas for modelling and inference from natural language processing, have demonstrated the ability to extract meaningful representations in an unsupervised way. This led to significant performance improvement in several downstream tasks. Clustering amino acids based on their physical-chemical properties to achieve reduced alphabets has been of interest in past research, but their application to PLMs or folding models is unexplored. Here, we investigate the efficacy of PLMs trained on reduced amino acid alphabets in capturing evolutionary information, and we explore how the loss of protein sequence information impacts learned representations and downstream task performance. Our empirical work shows that PLMs trained on the full alphabet and a large number of sequences capture fine details that are lost in alphabet reduction methods. We further show the ability of a structure prediction model(ESMFold) to fold CASP14 protein sequences translated using a reduced alphabet. For 10 proteins out of the 50 targets, reduced alphabets improve structural predictions with LDDT-Cα differences of up to 19%. Trained models and code are available at github.com/Ieremie/reduced-alph-PLM.

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

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

From interaction networks to interfaces, scanning intrinsically disordered regions using AlphaFold2

The revolution brought about by AlphaFold2 opens promising perspectives to unravel the complexity of protein-protein interaction networks. The analysis of interaction networks obtained from proteomics experiments does not systematically provide the delimitations of the interaction regions. This is of particular concern in the case of interactions mediated by intrinsically disordered regions, in which the interaction site is generally small. Using a dataset of protein-peptide complexes involving intrinsically disordered regions that are non-redundant with the structures used in AlphaFold2 training, we show that when using the full sequences of the proteins, AlphaFold2-Multimer only achieves 40% success rate in identifying the correct site and structure of the interface. By delineating the interaction region into fragments of decreasing size and combining different strategies for integrating evolutionary information, we manage to raise this success rate up to 90%. We obtain similar success rates using a much larger dataset of protein complexes taken from the ELM database. Beyond the correct identification of the interaction site, our study also explores specificity issues. We show the advantages and limitations of using the AlphaFold2 confidence score to discriminate between alternative binding partners, a task that can be particularly challenging in the case of small interaction motifs.

Concurrent prediction of RNA secondary structures with pseudoknots and local 3D motifs in an integer programming framework.

The prediction of RNA structure canonical base pairs from a single sequence, especially pseudoknotted ones, remains challenging in a thermodynamic models that approximates the energy of the local 3D motifs joining canonical stems. It has become more and more apparent in recent years that the structural motifs in the loops, composed of noncanonical interactions, are essential for the final shape of the molecule enabling its multiple functions. Our capacity to predict accurate 3D structures is also limited when it comes to the organization of the large intricate network of interactions that form inside those loops. We previously developed the integer programming framework RNA Motifs over Integer Programming (RNAMoIP) to reconcile RNA secondary structure and local 3D motif information available in databases. We further develop our model to now simultaneously predict the canonical base pairs (with pseudoknots) from base pair probability matrices with or without alignment. We benchmarked our new method over the all nonredundant RNAs below 150 nucleotides. We show that the joined prediction of canonical base pairs structure and local conserved motifs (i) improves the ratio of well-predicted interactions in the secondary structure, (ii) predicts well canonical and Wobble pairs at the location where motifs are inserted, (iii) is greatly improved with evolutionary information, and (iv) noncanonical motifs at kink-turn locations. The source code of the framework is available at https://gitlab.info.uqam.ca/cbe/RNAMoIP and an interactive web server at https://rnamoip.cbe.uqam.ca/.

Community-Aware Evolution Similarity for Link Prediction in Dynamic Social Networks

The link prediction problem is a time-evolving model in network science that has simultaneously abetted myriad applications and experienced extensive methodological improvement. Inferring the possibility of emerging links in dynamic social networks, also known as the dynamic link prediction task, is complex and challenging. In contrast to the link prediction in cross-sectional networks, dynamic link prediction methods need to cater to the actor-level temporal changes and associated evolutionary information regarding their micro- (i.e., link formation/deletion) and mesoscale (i.e., community formation) network structure. With the advent of abundant community detection algorithms, the research community has examined community-aware link prediction strategies in static networks. However, the same task in dynamic networks where, apart from the actors and links among them, their community pattern is also dynamic, is yet to be explored. Evolutionary community-aware information, including the associated link structure and temporal neighborhood changes, can effectively be mined to build dynamic similarity metrics for dynamic link prediction. This study aims to develop and integrate such dynamic features with machine learning algorithms for link prediction tasks in dynamic social networks. It also compares the performances of these features against well-known similarity metrics (i.e., ResourceAllocation) for static networks and a time series-based link prediction strategy in dynamic networks. These proposed features achieved high-performance scores, representing them as prospective candidates for both dynamic link prediction tasks and modeling the network growth.

Engineering Substrate Promiscuity of Nucleoside Phosphorylase Via an Insertions-Deletions Strategy.

Nucleoside phosphorylases (NPs) are the key enzymes in the nucleoside metabolism pathway and are widely employed for the synthesis of nucleoside analogs, which are difficult to access via conventional synthetic methods. NPs are generally classified as purine nucleoside phosphorylase (PNP) and pyrimidine or uridine nucleoside phosphorylase (PyNP/UP), based on their substrate preference. Here, based on the evolutionary information on the NP-I family, we adopted an insertions-deletions (InDels) strategy to engineer the substrate promiscuity of nucleoside phosphorylase AmPNPΔS2V102K, which exhibits both PNP and UP activities from a trimeric PNP (AmPNP) of Aneurinibacillus migulanus. Furthermore, the AmPNPΔS2V102K exerted phosphorylation activities toward arabinose nucleoside, fluorosyl nucleoside, and dideoxyribose, thereby broadening the unnatural-ribose nucleoside substrate spectrum of AmPNP. Finally, six purine nucleoside analogues were successfully synthesized, using the engineered AmPNPΔS2V102K instead of the traditional "two-enzymes PNP/UP" approach. These results provide deep insights into the catalytic mechanisms of the PNP and demonstrate the benefits of using the InDels strategy to achieve substrate promiscuity in an enzyme, as well as broadening the substrate spectrum of the enzyme.

Abstract IA010: Molecular archeology of cancer

Abstract Tumor development is driven by changes to the genome leading to fitness advantages underlying successive clonal expansions. As somatic changes occur across most or all cell cycles, the cancer genome carries an archeological record of its past. Over the past years, we have developed several approaches to mine that archeological record from the cancer genome, which we collectively call 'molecular archeology of cancer'. Using these approaches, we are able to infer the subclonal architecture of tumors and gain key insights into the order and timing of the genomic changes that occurred over their evolutionary history. We have applied these approaches in a large-scale pan-cancer setting, showing that intra-tumor heterogeneity is pervasive across cancers and that the timelines of tumor evolution span multiple years to decades. Key driver events in tumor evolution typically occur early, and copy number gains often accumulate as punctuated bursts, commonly after genome doubling. Late genome doubling is frequent in cancer evolution and is typically followed by an increase in the rate of copy number gains. Genome duplications affect the selection landscape of copy number losses, while only minimally impacting copy number gains. Extending our timing framework to time mutational signatures, we find that, as expected, environmentally associated mutations are associated with early tumor development. We also find evidence of episodic APOBEC mutagenesis across many tumors and clear timing patterns in many mutational processes that have unknown origins. Our approaches increase the evolutionary information that can be obtained from tumor genome sequences and, therefore, improve our understanding of the developmental history of cancer. Citation Format: Peter Van Loo. Molecular archeology of cancer [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Damage Repair: From Basic Science to Future Clinical Application; 2024 Jan 9-11; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2024;84(1 Suppl):Abstract nr IA010.

Evolution of protein domain repertoires of CALHM6.

Calcium (Ca2+) homeostasis is essential in conducting various cellular processes including nerve transmission, muscular movement, and immune response. Changes in Ca2+concentration in the cytoplasm are significant in bringing about various immune responses such as pathogen clearance and apoptosis. Various key players are involved in calcium homeostasis such as calcium binders, pumps, and channels. Sequence-based evolutionary information has recently been exploited to predict the biophysical behaviors of proteins, giving critical clues about their functionality. Ion channels are reportedly the first channels developed during evolution. Calcium homeostasis modulator protein 6 (CALHM6) is one such channel. Comprised of a single domain called Ca_hom_mod, CALHM6 is a stable protein interacting with various other proteins in calcium regulation. No previous attempt has been made to trace the exact evolutionary events in the domain of CALHM6, leaving plenty of room for exploring its evolution across a wide range of organisms. The current study aims to answer the questions by employing a computational-based strategy that used profile Hidden Markov Models (HMMs) to scan for the CALHM6 domain, integrated the data with a time-calibrated phylogenetic tree using BEAST and Mesquite, and visualized through iTOL. Around 4,000 domains were identified, and 14,000 domain gain, loss, and duplication events were observed at the end which also included various protein domains other than CALHM6. The data were analyzed concerning CALHM6 evolution as well as the domain gain, loss, and duplication of its interacting partners: Calpain, Vinculin, protein S100-A7, Thioredoxin, Peroxiredoxin, and Calmodulin-like protein 5. Duplication events of CALHM6 near higher eukaryotes showed its increasing complexity in structure and function. This in-silico phylogenetic approach applied to trace the evolution of CALHM6 was an effective approach to get a better understanding of the protein CALHM6.

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.

Population structure and connectivity among coastal and freshwater Kelp Gull (Larus dominicanus) populations from Patagonia.

The genetic identification of evolutionary significant units and information on their connectivity can be used to design effective management and conservation plans for species of concern. Despite having high dispersal capacity, several seabird species show population structure due to both abiotic and biotic barriers to gene flow. The Kelp Gull is the most abundant species of gull in the southern hemisphere. In Argentina it reproduces in both marine and freshwater environments, with more than 100,000 breeding pairs following a metapopulation dynamic across 140 colonies in the Atlantic coast of Patagonia. However, little is known about the demography and connectivity of inland populations. We aim to provide information on the connectivity of the largest freshwater colonies (those from Nahuel Huapi Lake) with the closest Pacific and Atlantic populations to evaluate if these freshwater colonies are receiving immigrants from the larger coastal populations. We sampled three geographic regions (Nahuel Huapi Lake and the Atlantic and Pacific coasts) and employed a reduced-representation genomic approach to genotype individuals for single-nucleotide polymorphisms (SNPs). Using clustering and phylogenetic analyses we found three genetic groups, each corresponding to one of our sampled regions. Individuals from marine environments are more closely related to each other than to those from Nahuel Huapi Lake, indicating that the latter population constitutes the first freshwater Kelp Gull colony to be identified as an evolutionary significant unit in Patagonia.

IProps: A comprehensive software tool for protein classification and analysis with automatic machine learning capabilities and model interpretation capabilities.

Protein classification is a crucial field in bioinformatics. The development of a comprehensive tool that can perform feature evaluation, visualization, automated machine learning, and model interpretation would significantly advance research in protein classification. However, there is a significant gap in the literature regarding tools that integrate all these essential functionalities. This paper presents iProps, a novel Python-based software package, meticulously crafted to fulfill these multifaceted requirements. iProps is distinguished by its proficiency in feature extraction, evaluation, automated machine learning, and interpretation of classification models. Firstly, iProps fully leverages evolutionary information and amino acid reduction information to propose or extend several numerical protein features that are independent of sequence length, including SC-PSSM, ORDip, TRC, CTDC-E, CKSAAGP-E, and so forth; at the same time, it also implements the calculation of 17 other numerical features within the software. iProps also provides feature combination operations for the aforementioned features to generate more hybrid features, and has added data balancing sampling processing as well as built-in classifier settings, among other functionalities. Thus, It can discern the most effective protein class recognition feature from a multitude of candidates, utilizing three automated machine learning algorithms to identify the most optimal classifiers and parameter settings. Furthermore, iProps generates a detailed explanatory report that includes 23 informative graphs derived from three interpretable models. To assess the performance of iProps, a series of numerical experiments were conducted using two well-established datasets. The results demonstrated that our software achieved superior recognition performance in every case. Beyond its contributions to bioinformatics, iProps broadens its applicability by offering robust data analysis tools that are beneficial across various disciplines, capitalizing on its automated machine learning and model interpretation capabilities. As an open-source platform, iProps is readily accessible and features an intuitive user interface, ensuring ease of use for individuals, even those without a background in programming. The source code of the software is available for download at the following website: https://github.com/LigosQ/iProps and https://gitee.com/LigosQ/iProps.

An Effective Computational Method for Predicting Self-Interacting Proteins Based on VGGNet Convolutional Neural Network and Gray-Level Co-occurrence Matrix.

Predicting Self-interacting proteins (SIPs) is a crucial area of research in predicting protein functions, as well as in understanding gene-disease and disease-drug associations. These interactions are integral to numerous cellular processes and play pivotal roles within cells. However, traditional methods for identifying SIPs through biological experiments are often expensive, time-consuming, and have long cycles. Therefore, the development of effective computational methods for accurately predicting SIPs is not only necessary but also presents a significant challenge. In this research, we introduce a novel computational prediction technique, VGGNGLCM, which leverages protein sequence data. This method integrates the VGGNet deep convolutional neural network (VGGN) with the Gray-Level Co-occurrence Matrix (GLCM) to detect Self-interacting proteins associations. Specifically, we initially utilized Position Specific Scoring Matrix (PSSM) to capture protein evolutionary information and integrated key features from PSSM using GLCM. We then employed VGGNet as a predictive classifier, leveraging its capabilities for powerful learning and classification prediction. Subsequently, the extracted features were input into the VGGNet deep convolutional neural network to identify Self-interacting proteins. To evaluate the performance of the VGGNGLCM model, we conducted experiments using yeast and human datasets, achieving average accuracies of 95.68% and 97.72% respectively. Additionally, we compared the prediction performance of the VGGNet classifier with that of the Convolutional Neural Network (CNN) and the state-of-the-art Support Vector Machine (SVM) using the same feature extraction method. We also compared the prediction ability of VGGNGLCM with other existing approaches. The comparison results further demonstrate the superior performance of VGGNGLCM over other prediction models in this domain. The experimental verification further strengthens the evidence that VGGNGLCM is effective and robust compared to existing methods. It also highlights the high accuracy and robustness of the VGGNGLCM model in predicting Self-interacting proteins (SIPs). Consequently, we believe that the VGGNGLCM method serves as a valuable computational tool and can catalyze extensive bioinformatics research related to SIPs prediction.

Evolutionary Structure Conservation and Covariance Scores.

Effective homology search for non-coding RNAs is frequently not possible via sequence similarity alone. Current methods leverage evolutionary information like structure conservation or covariance scores to identify homologs in organisms that are phylogenetically more distant. In this chapter, we introduce the theoretical background of evolutionary structure conservation and covariance score, and we show hands-on how current methods in the field are applied on example datasets.

Disentangling Stellar Age Estimates from Galactic Chemodynamical Evolution

Stellar ages are key for determining the formation history of the Milky Way, but are difficult to measure precisely. Furthermore, methods that use chemical abundances to infer ages may entangle the intrinsic evolution of stars with the chemodynamical evolution of the Galaxy. In this paper, we present a framework for making probabilistic predictions of stellar ages, and then quantify the contribution of both stellar evolution and Galactic chemical evolution to those predictions using SHapley Additive exPlanations. We apply this interpretable prediction framework to both a simulated Milky Way sample containing stars in a variety of evolutionary stages and an APOGEE-mocked sample of red clump stars. We find that in the former case, stellar evolution is the dominant driver for age estimates, while in the latter case, the more restricted evolutionary information causes the model to proxy ages through the chemical evolution model. We show that as a result of the use of nonintrinsic Galactic chemical information, trends estimated with the predicted ages, such as the age–metallicity relation, can deviate from the truth.