QuanTest2: benchmarking multiple sequence alignments using secondary structure prediction

DiMA: sequence diversity dynamics analyser for viruses

Sequence diversity is one of the major challenges in the design of diagnostic, prophylactic, and therapeutic interventions against viruses. DiMA is a novel tool that is big data-ready and designed to facilitate the dissection of sequence diversity dynamics for viruses. DiMA stands out from other diversity analysis tools by offering various unique features. DiMA provides a quantitative overview of sequence (DNA/RNA/protein) diversity by use of Shannon's entropy corrected for size bias, applied via a user-defined k-mer sliding window to an input alignment file, and each k-mer position is dissected to various diversity motifs. The motifs are defined based on the probability of distinct sequences at a given k-mer alignment position, whereby an index is the predominant sequence, while all the others are (total) variants to the index. The total variants are sub-classified into the major (most common) variant, minor variants (occurring more than once and of incidence lower than the major), and the unique (singleton) variants. DiMA allows user-defined, sequence metadata enrichment for analyses of the motifs. The application of DiMA was demonstrated for the alignment data of the relatively conserved Spike protein (2,106,985 sequences) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the relatively highly diverse pol gene (2637) of the human immunodeficiency virus-1 (HIV-1). The tool is publicly available as a web server (https://dima.bezmialem.edu.tr), as a Python library (via PyPi) and as a command line client (via GitHub).

A computational strategy to improve the activity of tyrosine phenol-lyase for the synthesis of L-DOPA

Enzymes with high catalytic activity and stability are essential for industrial production, yet most natural enzymes do not meet these requirements. Therefore, efficient strategies for enzyme engineering are crucial. In this study, we developed a cost-effective computational design strategy to enhance the activity of tyrosine phenol-lyase (TPL) for the production of L-DOPA. By integrating structural analysis with computational design, and guided by our understanding of conformational flexibility of TPL, we identified a region where enhanced stability is most likely to facilitate enzyme activity. We screened stabilizing mutations by Cartesian_ddg in Rosetta. After identifying single stabilizing mutations, we grouped the nearby positions harboring multiple stabilizing mutations and calculated the energy of combinatorial variants. We found two promising groups where most variants exhibited lower calculated energy than the wild-type. Experimental validation showed five variants in these groups exhibit increased activity, with the two best variants showing catalytic activity enhancements of 1.8-fold and 1.6-fold compared to the wild-type enzyme.

Comparative Analysis of Dehydrins from Woody Plant Species

We conducted analyses on 253 protein sequences (Pfam00257) derived from 25 woody plant species, including trees, shrubs, and vines. Our goal was to gain insights into their architectural types, biochemical characteristics, and potential involvement in mitigating abiotic stresses, such as drought, cold, or salinity. The investigated protein sequences (253) comprised 221 angiosperms (85 trees/shrubs and 36 vines) and 32 gymnosperms. Our sequence analyses revealed the presence of seven architectural types: Kn, KnS, SKn, YnKn, YnSKn, FSKn, and FnKn. The FSKn type predominated in tree and shrub dehydrins of both gymnosperms and angiosperms, while the YnSKn type was more prevalent in vine dehydrins. The YnSKn and YnKn types were absent in gymnosperms. Gymnosperm dehydrins exhibited a shift towards more negative GRAVY scores and Fold Indexes. Additionally, they demonstrated a higher Lys content and lower His content. By analyzing promoter sequences in the angiosperm species, including trees, shrubs, and vines, we found that these dehydrins are induced by the ABA-dependent and light-responsive pathways. The presence of stress- and hormone-related cis-elements suggests a protective effect against dehydration, cold, or salinity. These findings could serve as a foundation for future studies on woody dehydrins, especially in the context of biotechnological applications.

Leveraging protein language models for accurate multiple sequence alignments

Multiple sequence alignment (MSA) is a critical step in the study of protein sequence and function. Typically, MSA algorithms progressively align pairs of sequences and combine these alignments with the aid of a guide tree. These alignment algorithms use scoring systems based on substitution matrices to measure amino acid similarities. Although successful, standard methods struggle on sets of proteins with low sequence identity: the so-called twilight zone of protein alignment. For these difficult cases, another source of information is needed. Protein language models are a powerful new approach that leverages massive sequence data sets to produce high-dimensional contextual embeddings for each amino acid in a sequence. These embeddings have been shown to reflect physicochemical and higher-order structural and functional attributes of amino acids within proteins. Here, we present a novel approach to MSA, based on clustering and ordering amino acid contextual embeddings. Our method for aligning semantically consistent groups of proteins circumvents the need for many standard components of MSA algorithms, avoiding initial guide tree construction, intermediate pairwise alignments, gap penalties, and substitution matrices. The added information from contextual embeddings leads to higher accuracy alignments for structurally similar proteins with low amino-acid similarity. We anticipate that protein language models will become a fundamental component of the next generation of algorithms for generating MSAs.

Towards the accurate alignment of over a million protein sequences: Current state of the art

Large-scale genomics requires highly scalable and accurate multiple sequence alignment methods. Results collected over this last decade suggest accuracy loss when scaling up over a few thousand sequences. This issue has been actively addressed with a number of innovative algorithmic solutions that combine low-level hardware optimization with novel higher-level heuristics. This review provides an extensive critical overview of these recent methods. Using established reference datasets we conclude that albeit significant progress has been achieved, a unified framework able to consistently and efficiently produce high-accuracy large-scale multiple alignments is still lacking.

Three critical regions of the erythromycin resistance methyltransferase, ErmE, are required for function supporting a model for the interaction of Erm family enzymes with substrate rRNA

6-Methyladenosine modification of DNA and RNA is widespread throughout the three domains of life and often accomplished by a Rossmann-fold methyltransferase domain which contains conserved sequence elements directing S-adenosylmethionine cofactor binding and placement of the target adenosine residue into the active site. Elaborations to the conserved Rossman-fold and appended domains direct methylation to diverse DNA and RNA sequences and structures. Recently, the first atomic-resolution structure of a ribosomal RNA adenine dimethylase (RRAD) family member bound to rRNA was solved, TFB1M bound to helix 45 of 12S rRNA. Since erythromycin resistance methyltransferases are also members of the RRAD family, and understanding how these enzymes recognize rRNA could be used to combat their role in antibiotic resistance, we constructed a model of ErmE bound to a 23S rRNA fragment based on the TFB1M-rRNA structure. We designed site-directed mutants of ErmE based on this model and assayed the mutants by in vivo phenotypic assays and in vitro assays with purified protein. Our results and additional bioinformatic analyses suggest our structural model captures key ErmE-rRNA interactions and indicate three regions of Erm proteins play a critical role in methylation: the target adenosine binding pocket, the basic ridge, and the α4-cleft.

Ancestral Sequence Reconstruction for Exploring Alkaloid Evolution

The complex and bioactive monoterpene indole alkaloids (MIAs) found in Catharanthus roseus and related species are the products of many millions of years of evolution through mutation and natural selection. Ancestral sequence reconstruction (ASR) is a method that combines phylogenetic analysis and experimental biochemistry to infer details about past events in protein evolution. Here, I propose that ASR could be leveraged to understand how enzymes catalyzing the formation of complex alkaloids arose over evolutionary time. I discuss the steps of ASR, including sequence selection, multiple sequence alignment, tree inference, and the generation and characterization of inferred ancestral enzymes.

Benchmarking Orthogroup Inference Accuracy: Revisiting Orthobench

Orthobench is the standard benchmark to assess the accuracy of orthogroup inference methods. It contains 70 expert-curated reference orthogroups (RefOGs) that span the Bilateria and cover a range of different challenges for orthogroup inference. Here, we leveraged improvements in tree inference algorithms and computational resources to reinterrogate these RefOGs and carry out an extensive phylogenetic delineation of their composition. This phylogenetic revision altered the membership of 31 of the 70 RefOGs, with 24 subject to extensive revision and 7 that required minor changes. We further used these revised and updated RefOGs to provide an assessment of the orthogroup inference accuracy of widely used orthogroup inference methods. Finally, we provide an open-source benchmarking suite to support the future development and use of the Orthobench benchmark.

RNAconTest: comparing tools for noncoding RNA multiple sequence alignment based on structural consistency

The importance of noncoding RNA sequences has become increasingly clear over the past decade. New RNA families are often detected and analyzed using comparative methods based on multiple sequence alignments. Accordingly, a number of programs have been developed for aligning and deriving secondary structures from sets of RNA sequences. Yet, the best tools for these tasks remain unclear because existing benchmarks contain too few sequences belonging to only a small number of RNA families. RNAconTest (RNA consistency test) is a new benchmarking approach relying on the observation that secondary structure is often conserved across highly divergent RNA sequences from the same family. RNAconTest scores multiple sequence alignments based on the level of consistency among known secondary structures belonging to reference sequences in their output alignment. Similarly, consensus secondary structure predictions are scored according to their agreement with one or more known structures in a family. Comparing the performance of 10 popular alignment programs using RNAconTest revealed that DAFS, DECIPHER, LocARNA, and MAFFT created the most structurally consistent alignments. The best consensus secondary structure predictions were generated by DAFS and LocARNA (via RNAalifold). Many of the methods specific to noncoding RNAs exhibited poor scalability as the number or length of input sequences increased, and several programs displayed substantial declines in score as more sequences were aligned. Overall, RNAconTest provides a means of testing and improving tools for comparative RNA analysis, as well as highlighting the best available approaches. RNAconTest is available from the DECIPHER website (http://DECIPHER.codes/Downloads.html).

Deep learning methods in protein structure prediction

Protein Structure Prediction is a central topic in Structural Bioinformatics. Since the '60s statistical methods, followed by increasingly complex Machine Learning and recently Deep Learning methods, have been employed to predict protein structural information at various levels of detail. In this review, we briefly introduce the problem of protein structure prediction and essential elements of Deep Learning (such as Convolutional Neural Networks, Recurrent Neural Networks and basic feed-forward Neural Networks they are founded on), after which we discuss the evolution of predictive methods for one-dimensional and two-dimensional Protein Structure Annotations, from the simple statistical methods of the early days, to the computationally intensive highly-sophisticated Deep Learning algorithms of the last decade. In the process, we review the growth of the databases these algorithms are based on, and how this has impacted our ability to leverage knowledge about evolution and co-evolution to achieve improved predictions. We conclude this review outlining the current role of Deep Learning techniques within the wider pipelines to predict protein structures and trying to anticipate what challenges and opportunities may arise next.

Kalign 3: multiple sequence alignment of large data sets

MOTIVATION

Kalign is an efficient multiple sequence alignment (MSA) program capable of aligning thousands of protein or nucleotide sequences. However, current alignment problems involving large numbers of sequences are exceeding Kalign's original design specifications. Here we present a completely re-written and updated version to meet current and future alignment challenges.

RESULTS

Kalign now uses a SIMD accelerated version of the bit-parallel Gene Myers algorithm to estimate pariwise distances, adopts a sequence embedding strategy and the bi-secting K-means algorithm to rapidly construct guide trees for thousands of sequences. The new version maintains high alignment accuracy on both protein and nucleotide alignments and scales better than other MSA tools.

AVAILABILITY

The source code of Kalign and code to reproduce the results are found here: https://github.com/timolassmann/kalign.