Search PubMedSearch

PubMed · 41959221

Fast, accurate construction of multiple sequence alignments from protein language embeddings.

Abstract

Multiple sequence alignment (MSA) is a foundational task in computational biology, underpinning protein structure prediction, evolutionary analysis, and domain annotation. Traditional MSA algorithms rely on pairwise amino acid substitution matrices derived from conserved protein families. While effective for aligning closely related sequences, these scoring schemes struggle in the low-identity "twilight zone." Here, we present a new approach for constructing MSAs leveraging amino acid embeddings generated by protein language models (PLMs), which capture rich evolutionary and contextual information from massive and diverse sequence datasets. We introduce a windowed reciprocal-weighted embedding similarity metric that is surprisingly effective in identifying corresponding amino acids across sequences. Building on this metric, we develop ARIES (Alignment via RecIprocal Embedding Similarity), an algorithm that constructs a PLM-generated template embedding and aligns each sequence to this template via dynamic time warping in order to build a global MSA. Across diverse benchmark datasets, ARIES achieves higher accuracies than existing state-of-the-art approaches, especially in low-identity regimes where traditional methods degrade, while scaling almost linearly with the number of sequences to be aligned. Together, these results provide the first large-scale demonstration of the power of PLMs for accurate and scalable MSA construction across protein families of varying sizes and levels of similarity, highlighting the potential of PLMs to transform comparative sequence analysis.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Minh Hoang, Isabel Armour-Garb, Mona Singh. 2026-03-13. Fast, accurate construction of multiple sequence alignments from protein language embeddings.. https://doi.org/10.64898/2026.01.02.697423

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Decoding TnsC Filament Assembly in CRISPR-Associated Transposons Using Interpretable Deep Learning and Molecular Simulations.

CRISPR-associated transposons (CASTs) enable programmable DNA integration, yet how the TnsC regulator forms processive filaments on DNA to coordinate RNA-guided transposition in type V-K CAST systems remains unknown. Here, we integrate large-scale molecular simulations, interpretable deep learning using graph attention networks (GATs), and causal inference analyses to define the molecular determinants of TnsC filament nucleation and elongation. We show that TnsC nucleates by inducing localized DNA deformation that propagates along extended filaments, with Granger causality revealing that TnsC motions precede and predict DNA deformation. Interpretable GAT models demonstrate that elongation is determined during early recognition between incoming and DNA-bound subunits, followed by structural reorganization that regenerates the recruitment interface and enables processive assembly. These results elucidate the molecular mechanism of processive TnsC filament assembly and explain why isolated TnsC filaments preferentially elongate in the 5' → 3' direction, while accessory transposition factors can reshape the interaction landscape and alter filament growth polarity. Together, these findings advance our understanding of CAST function and inform the engineering of programmable DNA integration platforms. Beyond CAST systems, this work introduces an interpretable GAT approach as a general and transferable deep learning strategy for uncovering molecular mechanisms in biological systems, while demonstrating the power of causal inference for dissecting directional relationships in molecular dynamics.

Deep Learning

Using deep learning models as a genetic architecture for the simulation of breeding schemes.

In several simulation studies, long-term selection led to the rapid depletion of genetic variance. These outcomes differ from real-life observations that we aim to replicate, thereby highlighting a fundamental limitation of current classical quantitative genetic simulation models. Deep learning (DL) models have demonstrated promising results in capturing complex interactions essential for maintaining genetic variance; thus, we hypothesize that DL-based genetic simulation models may preserve more genetic variance than classical models, because the biological pathways underlying complex traits exhibit interactions that classical models ignore. The primary objective of this study was to introduce alternative DL-based genetic simulation models and compare them with classical genetic simulation models in terms of their retention of additive genetic variance under truncation selection in a simulated full-sib pig breeding scheme using real haplotypes as founders. After 20 generations of directional truncation selection, the classical models (A, ADAA, and ADAAADDD) retained between 55% and 64% of their initial additive genetic variance. In contrast, while the DL_simple model lost all its additive variance, the DL medium retained 92% to 98% of its additive variance, and the DL_complex model's initial additive variance increased by 296% to 314%. This paper introduces DL-based genetic simulation models and concludes that their ability to retain additive genetic variance depends on the models' architectural complexity. When sufficiently complex, DL-based models exhibit greater retention of additive genetic variance because they intrinsically capture epistatic interactions that are converted into additive variance, as selection progresses, thus, affirming the role of non-additive genetic effects in maintaining long-term genetic variation.

Deep Learning

High-accuracy SNV calling for bacterial isolates using deep learning with AccuSNV.

Accurate detection of mutations within bacterial species is critical for fundamental studies of microbial evolution, reconstruction of transmission events, and identification of antimicrobial resistance mutations. Although many tools have been developed to identify single-nucleotide variants (SNVs) from whole-genome sequencing, they often suffer from high false-positive rates owing to the complexity of bacterial genomes and the need for different filtering cutoffs across sample types and sequencing depths. As data sets increase in size, the manual filtering required for high accuracy presents a significant obstacle. Here, we present AccuSNV, a novel deep learning-based tool for high-precision and automated bacterial SNV calling. Unlike traditional methods that process one sample at a time, AccuSNV leverages a convolutional neural network (CNN) that integrates alignment information across multiple samples, enhancing precision through learned across-sample patterns. We evaluate AccuSNV against seven popular SNV-calling tools using simulated data from six bacterial species with varied sequencing depths, numbers of isolates, mutations, and divergence levels. To further validate its real-world utility, we test AccuSNV on multiple curated bacterial data sets containing reported SNVs. In both simulated and real-world scenarios, AccuSNV consistently achieves the best performance. Moreover, AccuSNV provides comprehensive user-friendly downstream analysis modules and outputs, including mutation annotation information, phylogenetic inference, d N/d S calculations, and optional manual filtering. Together with the automated deep learning-based calling, these features make AccuSNV broadly accessible to users with different levels of computational expertise.

Deep Learning