Secondary Structure Alignment Research Articles

REDalign: accurate RNA structural alignment using residual encoder-decoder network

BackgroundRNA secondary structural alignment serves as a foundational procedure in identifying conserved structural motifs among RNA sequences, crucially advancing our understanding of novel RNAs via comparative genomic analysis. While various computational strategies for RNA structural alignment exist, they often come with high computational complexity. Specifically, when addressing a set of RNAs with unknown structures, the task of simultaneously predicting their consensus secondary structure and determining the optimal sequence alignment requires an overwhelming computational effort of O(L6)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$O(L^6)$$\\end{document} for each RNA pair. Such an extremely high computational complexity makes these methods impractical for large-scale analysis despite their accurate alignment capabilities.ResultsIn this paper, we introduce REDalign, an innovative approach based on deep learning for RNA secondary structural alignment. By utilizing a residual encoder-decoder network, REDalign can efficiently capture consensus structures and optimize structural alignments. In this learning model, the encoder network leverages a hierarchical pyramid to assimilate high-level structural features. Concurrently, the decoder network, enhanced with residual skip connections, integrates multi-level encoded features to learn detailed feature hierarchies with fewer parameter sets. REDalign significantly reduces computational complexity compared to Sankoff-style algorithms and effectively handles non-nested structures, including pseudoknots, which are challenging for traditional alignment methods. Extensive evaluations demonstrate that REDalign provides superior accuracy and substantial computational efficiency.ConclusionREDalign presents a significant advancement in RNA secondary structural alignment, balancing high alignment accuracy with lower computational demands. Its ability to handle complex RNA structures, including pseudoknots, makes it an effective tool for large-scale RNA analysis, with potential implications for accelerating discoveries in RNA research and comparative genomics.

RNA threading with secondary structure and sequence profile.

RNA threading aims to identify remote homologies for template-based modeling of RNA 3D structure. Existing RNA alignment methods primarily rely on secondary structure alignment. They are often time- and memory-consuming, limiting large-scale applications. In addition, the accuracy is far from satisfactory. Using RNA secondary structure and sequence profile, we developed a novel RNA threading algorithm, named RNAthreader. To enhance the alignment process and minimize memory usage, a novel approach has been introduced to simplify RNA secondary structures into compact diagrams. RNAthreader employs a two-step methodology. Initially, integer programming and dynamic programming are combined to create an initial alignment for the simplified diagram. Subsequently, the final alignment is obtained using dynamic programming, taking into account the initial alignment derived from the previous step. The benchmark test on 80 RNAs illustrates that RNAthreader generates more accurate alignments than other methods, especially for RNAs with pseudoknots. Another benchmark, involving 30 RNAs from the RNA-Puzzles experiments, exhibits that the models constructed using RNAthreader templates have a lower average RMSD than those created by alternative methods. Remarkably, RNAthreader takes less than two hours to complete alignments with ∼5000 RNAs, which is 3-40 times faster than other methods. These compelling results suggest that RNAthreader is a promising algorithm for RNA template detection. https://yanglab.qd.sdu.edu.cn/RNAthreader.

A Multiplexing Activity-Based Protein-Profiling Platform for Dissection of a Native Bacterial Xyloglucan-Degrading System.

Bacteria and yeasts grow on biomass polysaccharides by expressing and excreting a complex array of glycoside hydrolase (GH) enzymes. Identification and annotation of such GH pools, which are valuable commodities for sustainable energy and chemistries, by conventional means (genomics, proteomics) are complicated, as primary sequence or secondary structure alignment with known active enzymes is not always predictive for new ones. Here we report a "low-tech", easy-to-use, and sensitive multiplexing activity-based protein-profiling platform to characterize the xyloglucan-degrading GH system excreted by the soil saprophyte, Cellvibrio japonicus, when grown on xyloglucan. A suite of activity-based probes bearing orthogonal fluorophores allows for the visualization of accessory exo-acting glycosidases, which are then identified using biotin-bearing probes. Substrate specificity of xyloglucanases is directly revealed by imbuing xyloglucan structural elements into bespoke activity-based probes. Our ABPP platform provides a highly useful tool to dissect xyloglucan-degrading systems from various sources and to rapidly select potentially useful ones. The observed specificity of the probes moreover bodes well for the study of other biomass polysaccharide-degrading systems, by modeling probe structures to those of desired substrates.

Compensatory Base Changes in ITS2 Secondary Structure Alignment, Modelling, and Molecular Phylogeny: An Integrated Approach to Improve Species Delimitation in Tulasnella (Basidiomycota).

The delimitation of species of Tulasnella has been extensively studied, mainly at the morphological (sexual and asexual states) and molecular levels-showing ambiguity between them. An integrative species concept that includes characteristics such as molecular, ecology, morphology, and other information is crucial for species delimitation in complex groups such as Tulasnella. The aim of this study is to test evolutionary relationships using a combination of alignment-based and alignment-free distance matrices as an alternative molecular tool to traditional methods, and to consider the secondary structures and CBCs from ITS2 (internal transcribed spacer) sequences for species delimitation in Tulasnella. Three phylogenetic approaches were plotted: (i) alignment-based, (ii) alignment-free, and (iii) a combination of both distance matrices using the DISTATIS and pvclust libraries from an R package. Finally, the secondary structure consensus was modeled by Mfold, and a CBC analysis was obtained to complement the species delimitation using 4Sale. The phylogenetic tree results showed delimited monophyletic clades in Tulasnella spp., where all 142 Tulasnella sequences were divided into two main clades A and B and assigned to seven species (T. asymmetrica, T. andina, T. eichleriana ECU6, T. eichleriana ECU4 T. pinicola, T. violea), supported by bootstrap values from 72% to 100%. From the 2D secondary structure alignment, three types of consensus models with helices and loops were obtained. Thus, T. albida belongs to type I; T. eichleriana, T. tomaculum, and T. violea belong to type II; and T. asymmetrica, T. andina, T. pinicola, and T. spp. (GER) belong to type III; each type contains four to six domains, with nine CBCs among these that corroborate different species.

BpRNA-align: improved RNA secondary structure global alignment for comparing and clustering RNA structures.

Ribonucleic acid (RNA) is a polymeric molecule that is fundamental to biological processes, with structure being more highly conserved than primary sequence and often key to its function. Advances in RNA structure characterization have resulted in an increase in the number of accurate secondary structures. The task of uncovering common RNA structural motifs with a collective function through structural comparison, providing a level of similarity, remains challenging and could be used to improve RNA secondary structure databases and discover new RNA families. In this work, we present a novel secondary structure alignment method, bpRNA-align. bpRNA-align is a customized global structural alignment method, utilizing an inverted (gap extend costs more than gap open) and context-specific affine gap penalty along with a structural, feature-specific substitution matrix to provide similarity scores. We evaluate our similarity scores in comparison to other methods, using affinity propagation clustering, applied to a benchmarking data set of known structure types. bpRNA-align shows improvement in clustering performance over a broad range of structure types.

Whole genome analyses for c-type cytochromes associated with respiratory chains in the extreme thermophile, Thermus thermophilus.

In thermophilic microorganisms, c-type cytochrome (cyt) proteins mainly function in the respiratory chain as electron carriers. Genome analyses at the beginning of this century revealed a variety of genes harboring the heme c motif. Here, we describe the results of surveying genes with the heme c motif, CxxCH, in a genome database comprising four strains of Thermus thermophilus, including strain HB8, and the confirmation of 19 c-type cytochromes among 27 selected genes. We analyzed the 19 genes, including the expression of four, by a bioinformatics approach to elucidate their individual attributes. One of the approaches included an analysis based on the secondary structure alignment pattern between the heme c motif and the 6th ligand. The predicted structures revealed many cyt c domains with fewer β-strands, such as mitochondrial cyt c, in addition to the β-strand unique to Thermus inserted in cyt c domains, as in T. thermophilus cyt c552 and caa3 cyt c oxidase subunit IIc. The surveyed thermophiles harbor potential proteins with a variety of cyt c folds. The gene analyses led to the development of an index for the classification of cyt c domains. Based on these results, we propose names for T. thermophilus genes harboring the cyt c fold.

In silico analysis of missense variants of the C1qA gene related to infection and autoimmune diseases

ObjectivesC1q is a key activator of the classical pathway of the complement system and exerts consequences relating to opsonization and phagocytosis. The C1qA gene is one of three genes encoding the C1q molecule. Defects in C1q, and especially in C1qA, have been linked to an increased susceptibility to infection, sepsis, and systemic lupus erythematosus. These defects could arise from missense single nucleotide polymorphisms (SNPs) and their deleterious impacts on protein structure and function. Thus, identifying high-risk missense SNPs in C1qA has become a necessity if we are to identify appropriate measures for prevention and management of affected patients. MethodsA comprehensive in silico study was conducted to screen the 184 missense SNPs in the C1qA gene using different tools with different algorithms and approaches. We investigated the impact of SNPs on protein function, stability, and structure. In addition, we identified the location of the SNPs on protein domains, secondary structure alignment, and the phylogenetic conservation of their positions. ResultsOf the 184 missense SNPs, 10 SNPs were predicted to be the most damaging to protein function and structure. ConclusionTen missense SNPs were predicted to have the highest risk of damaging protein function and structure, thus leading to infection, sepsis, and systemic lupus erythematosus. These 10 SNPs constitute the best candidates for further experimental investigations.

Biophysical studies on latarcins: antimicrobial peptides from spider venom

The ant spider Lachesena tarabaevi produces one of the most potent venom that is largely composed of membrane active peptides. A group of seven linear antimicrobial peptides isolated from the venom of this spider which belong to the latarcin family are reported to be cytolytic and highly membrane active[1]. However, there is no single study where the secondary structure, membrane alignment of all peptides have been systematically addressed. Here, using circular dichroism(CD) spectroscopy, oriented CD, solid-state 15N-NMR and fluorescence spectroscopy, we have characterized the secondary structure and helix alignment of all seven peptides in their membrane-bound state [2].

Bioinformatics and 3D homology modelling of AsAlaDH from Amycolatopsis sulphurea

Alanine dehydrogenase (AlaDH) (E.C.1.4.1.1) is an enzyme that catalyzes the interconversion of pyruvate and alanine. This enzyme has the key catalytic role for the sporulation of microorganism and synthesis of the many amino acids, proteins, and peptidoglycan layers in the microorganisms. Amycolatopsis sulphurea one of the strains of Amycolatopsis genus within the family Pseudonocardiaceae has capable to produce different antibiotics such as Ristocetin, Vancomycin, and Epoxyquinomicin as well as to biodegrade the bioplastic (poly-lactic acid (PLA) films). The 3D homology model of Alanine dehydrogenase from Amycolatopsis sulphurea was carried out through I-TASSER. The interaction of L-alanine and active site amino acids of the enzyme was determined by docking in silico via AutoDock Vina program. Protein secondary structures were predicted with EMBOSS tool garnier. Structural and functional analysis and determination of Physico-chemical properties of AsAlaDH were performed by using different bioinformatics tools. The secondary structure and multiple alignment analysis of alanine dehydrogenase displayed that there are conserved amino acid residues of AlaDH's from different microorganisms.

A systematic approach to inserting split inteins for Boolean logic gate engineering and basal activity reduction

Split inteins are powerful tools for seamless ligation of synthetic split proteins. Yet, their use remains limited because the already intricate split site identification problem is often complicated by the requirement of extein junction sequences. To address this, we augment a mini-Mu transposon-based screening approach and devise the intein-assisted bisection mapping (IBM) method. IBM robustly reveals clusters of split sites on five proteins, converting them into AND or NAND logic gates. We further show that the use of inteins expands functional sequence space for splitting a protein. We also demonstrate the utility of our approach over rational inference of split sites from secondary structure alignment of homologous proteins, and that basal activities of highly active proteins can be mitigated by splitting them. Our work offers a generalizable and systematic route towards creating split protein-intein fusions for synthetic biology.

First phylogenetic analysis of Dryophthorinae (Coleoptera, Curculionidae) based on structural alignment of ribosomal DNA reveals Cenozoic diversification.

Dryophthorinae is an economically important, ecologically distinct, and ubiquitous monophyletic group of pantropical weevils with more than 1,200 species in 153 genera. This study provides the first comprehensive phylogeny of the group with the aim to provide insights into the process and timing of diversification of phytophagous insects, inform classification and facilitate predictions. The taxon sampling is the most extensive to date and includes representatives of all five dryophthorine tribes and all but one subtribe. The phylogeny is based on secondary structural alignment of 18S and 28S rRNA totaling 3,764 nucleotides analyzed under Bayesian and maximum likelihood inference. We used a fossil‐calibrated relaxed clock model with two approaches, node‐dating and fossilized birth‐death models, to estimate divergence times for the subfamily. All tribes except the species‐rich Rhynchophorini were found to be monophyletic, but higher support is required to ascertain the paraphyly of Rhynchophorini with more confidence. Nephius is closely related to Dryophthorini and Stromboscerini, and there is strong evidence for paraphyly of Sphenophorina. We find a large gap between the divergence of Dryophthorinae from their sister group Platypodinae in the Jurassic‐Cretaceous boundary and the diversification of extant species in the Cenozoic, highlighting the role of coevolution with angiosperms in this group.

DSARna: RNA Secondary Structure Alignment Based on Digital Sequence Representation

With increasing applications and development of high-throughput sequencing, knowledge of the primary structure of RNA has expanded exponentially. Moreover, the function of RNA is determined by the secondary or higher RNA structure, and similar structures are related to similar functions, such as the secondary clover structure of tRNA. Therefore, RNA structure alignment is an important subject in computational biology and bioinformatics to predict function accurately. However, the traditional RNA structure alignment algorithms have some drawbacks such as high complexity and easy loss of secondary structure information. To study R,,NA secondary structure alignment according to the shortcomings of existing secondary structure alignment algorithms and the characteristics of RNA secondary structure. We propose a new digital sequence RNA structure representation algorithm named "DSARna". Then based on a dynamic programming algorithm, the scoring matrix and binary path matrix are simultaneously constructed. The backtracking path is identified in the path matrix, and the optimal result is predicted according to the path length. Upon comparison with the existing SimTree algorithm through experimental analysis, the proposed method showed higher accuracy and could ensure that the structural information is not easily lost in terms of improved specificity, sensitivity, and the Matthews correlation coefficient.

SCOT: Rethinking the classification of secondary structure elements.

Secondary structure classification is one of the most important issues in structure-based analyses due to its impact on secondary structure prediction, structural alignment and protein visualization. There are still open challenges concerning helix and sheet assignments which are currently not addressed by a single multi-purpose software. We introduce SCOT (Secondary structure Classification On Turns) as a novel secondary structure element assignment software which supports the assignment of turns, right-handed α-, 310- and π-helices, left-handed α- and 310-helices, 2.27- and polyproline II helices, β-sheets and kinks. We demonstrate that the introduction of helix Purity values enables a clear differentiation between helix classes. SCOT's unique strengths are highlighted by comparing it to six state-of-the-art methods (DSSP, STRIDE, ASSP, SEGNO, DISICL and SHAFT). The assignment approaches were compared concerning geometric consistency, protein structure quality and flexibility dependency and their impact on secondary structure element-based structural alignments. We show that only SCOT's combination of hydrogen bonds, geometric criteria and dihedral angles enables robust assignments independent of the structure quality and flexibility. We demonstrate that this combination and the elaborate kink detection lead to SCOT's clear superiority for protein alignments. As the resulting helices and strands are provided in a PDB conform output format, they can immediately be used for structure alignment algorithms. Taken together, the application of our new method and the straight-forward visualization using the accompanying PyMOL scripts enable the comprehensive analysis of regular backbone geometries in proteins. https://this-group.rocks. Supplementary data are available at Bioinformatics online.

Pathogen spread and globalization: The case of Pestivirus heterogeneity in southern Italy

Pestiviruses are responsible for widespread diseases affecting cattle, pigs and other ruminants, presenting a wide range of clinical manifestations, with significant impact on animal production. Given the recent various reports of a relatively high number of new strains and atypical genomic variants, in the present study, ninety-seven genomic sequences from southern Italy have been evaluated applying the palindromic nucleotide substitutions method, based on 5’-UTR secondary structure alignment and computing genetic distance among strains in the internal ribosome entry site. Sequence analysis revealed a highly heterogeneous virus population, indicating the introduction of virus variants of Bovine viral diarrhea virus and Border disease virus species from foreign countries. The application of different analytical procedures was useful to avoid interpretation difficulties. Circulation of heterogeneous virus populations showed the need for more accurate epidemiological investigations and stringent veterinary controls to protect animal health and welfare.

Computational Approach for Structural Feature Determination of Grapevine NHX Antiporters.

Plant NHX antiporters are responsible for monovalent cation/H+ exchange across cellular membranes and play therefore a critical role for cellular pH regulation, Na+ and K+ homeostasis, and salt tolerance. Six members of grapevine NHX family (VvNHX1-6) have been structurally characterized. Phylogenetic analysis revealed their organization in two groups: VvNHX1-5 belonging to group I (vacuolar) and VvNHX6 belonging to group II (endosomal). Conserved domain analysis of these VvNHXs indicates the presence of different kinds of domains. Out of these, two domains function as monovalent cation-proton antiporters and one as the aspartate-alanine exchange; the remaining are not yet with defined function. Overall, VvNHXs proteins are typically made of 11-13 putative transmembrane regions at their N-terminus which contain the consensus amiloride-binding domain in the 3rd TM domain and a cation-binding site in between the 5th and 6th TM domain, followed by a hydrophilic C-terminus that is the target of several and diverse regulatory posttranslational modifications. Using a combination of primary structure analysis, secondary structure alignments, and the tertiary structural models, the VvNHXs revealed mainly 18 α helices although without β sheets. Homology modeling of the 3D structure showed that VvNHX antiporters are similar to the bacterial sodium proton antiporters MjNhaP1 (Methanocaldococcus jannaschii) and PaNhaP (Pyrococcus abyssi).

Large-Scale Analyses of Site-Specific Evolutionary Rates across Eukaryote Proteomes Reveal Confounding Interactions between Intrinsic Disorder, Secondary Structure, and Functional Domains.

Various structural and functional constraints govern the evolution of protein sequences. As a result, the relative rates of amino acid replacement among sites within a protein can vary significantly. Previous large-scale work on Metazoan (Animal) protein sequence alignments indicated that amino acid replacement rates are partially driven by a complex interaction among three factors: intrinsic disorder propensity; secondary structure; and functional domain involvement. Here, we use sequence-based predictors to evaluate the effects of these factors on site-specific sequence evolutionary rates within four eukaryotic lineages: Metazoans; Plants; Saccharomycete Fungi; and Alveolate Protists. Our results show broad, consistent trends across all four Eukaryote groups. In all four lineages, there is a significant increase in amino acid replacement rates when comparing: (i) disordered vs. ordered sites; (ii) random coil sites vs. sites in secondary structures; and (iii) inter-domain linker sites vs. sites in functional domains. Additionally, within Metazoans, Plants, and Saccharomycetes, there is a strong confounding interaction between intrinsic disorder and secondary structure—alignment sites exhibiting both high disorder propensity and involvement in secondary structures have very low average rates of sequence evolution. Analysis of gene ontology (GO) terms revealed that in all four lineages, a high fraction of sequences containing these conserved, disordered-structured sites are involved in nucleic acid binding. We also observe notable differences in the statistical trends of Alveolates, where intrinsically disordered sites are more variable than in other Eukaryotes and the statistical interactions between disorder and other factors are less pronounced.

Mutations and functional analysis of 14-3-3 stress response protein from Triticum aestivum: An evolutionary analysis through in silico structural biochemistry approach

Wheat (Triticum aestivum), having high nutritional values is one of the staple food of most of the countries in the world. The productivity of the crop decreases drastically when it encounters various abiotic stresses, most common of which are heat, drought, flood and salinity. There is a crucial role of stress response proteins for the survival of the crops in stress conditions. So the study of wheat stress response proteins is of great importance to raise wheat production in different stress conditions. In this study, we analysed 14-3-3 protein, a stress response protein that is expressed in three major stresses, for example heat, drought and salinity and helps the plants to survive in those conditions. Effect of mutations in the 14-3-3 sequence was predicted using its domain, secondary structure and multiple sequence alignment of amino acid sequences from wheat and its related species. The functional diversity of the protein in different species was correlated with mutations, change in secondary structure and the evolutionary relatedness of the protein in different species. This is the first novel work for analysing the mutational effect on the structure and function of a stress response protein (14-3-3) from Triticum aestivum and its related species.

Bovine viral diarrhea virus type 1 current taxonomy according to palindromic nucleotide substitutions method

Pestivirus bovine viral diarrhea virus type 1 species is responsible for cosmopolitan diseases affecting cattle and other ruminants, presenting a wide range of clinical manifestations, with relevant impact on zootechnic production. Understanding genomic characteristic and virus taxonomy is fundamental in order to sustain control and prophylactic programs. Given the recent various studies reporting a relatively high number of new strains, in particular from Asian countries, in the present study, four hundred-eighty-two genomic sequences have been evaluated applying the palindromic nucleotide substitutions method for genotyping. Based on the secondary structure alignment and computing genetic distance among strains in the 5′ untranslated region of Pestivirus RNA, the current taxonomy of the species was reviewed. Twenty-two genotypes have been identified, applying a nomenclature based on divergence in the genus.

Biochemical and genetic characterization of a novel metallo-\u03b2-lactamase from marine bacterium Erythrobacter litoralis HTCC 2594

Metallo-β-lactamases (MBLs) are a group of enzymes that can inactivate most commonly used β-lactam-based antibiotics. Among MBLs, New Delhi metallo-β-lactamase-1 (NDM-1) constitutes an urgent threat to public health as evidenced by its success in rapidly disseminating worldwide since its first discovery. Here we report the biochemical and genetic characteristics of a novel MBL, ElBla2, from the marine bacterium Erythrobacter litoralis HTCC 2594. This enzyme has a higher amino acid sequence similarity to NDM-1 (56%) than any previously reported MBL. Enzymatic assays and secondary structure alignment also confirmed the high similarity between these two enzymes. Whole genome comparison of four Erythrobacter species showed that genes located upstream and downstream of elbla2 were highly conserved, which may indicate that elbla2 was lost during evolution. Furthermore, we predicted two prophages, 13 genomic islands and 25 open reading frames related to insertion sequences in the genome of E. litoralis HTCC 2594. However, unlike NDM-1, the chromosome encoded ElBla2 did not locate in or near these mobile genetic elements, indicating that it cannot transfer between strains. Finally, following our phylogenetic analysis, we suggest a reclassification of E. litoralis HTCC 2594 as a novel species: Erythrobacter sp. HTCC 2594.

CompPSA: A Component-Based Pairwise RNA Secondary Structure Alignment Algorithm

CompPSA: A Component-Based Pairwise RNA Secondary Structure Alignment Algorithm