Comparison Of Protein Sequences Research Articles

Multi-omics analysis provides new insights into the changes of important nutrients and fructose metabolism in loquat bud sport mutant.

Bud sport is a common and stable somatic variation in perennial fruit trees, and often leads to significant modification of fruit traits and affects the breeding value. To investigate the impact of bud sport on the main metabolites in the fruit of white-fleshed loquat, we conducted a multi-omics analysis of loquat fruits at different developmental stages of a white-fleshed bud sport mutant of Dongting loquat (TBW) and its wild type (TBY). The findings from the detection of main fruit quality indices and metabolites suggested that bud sport resulted in a reduction in the accumulation of carotenoids, fructose, titratable acid and terpenoids at the mature stage of TBW, while leading to the accumulation of flavonoids, phenolic acids, amino acids and lipids. The comparably low content of titratable acid further enhances the balanced and pleasent taste profile of TBW. Expression patterns of differentially expressed genes involved in fructose metabolism exhibited a significant increase in the expression level of S6PDH (EVM0006243, EVM0044405) prior to fruit maturation. The comparison of protein sequences and promoter region of S6PDH between TBY and TBW revealed no structural variations that would impact gene function or expression, indicating that transcription factors may be responsible for the rapid up-regulation of S6PDH before maturation. Furthermore, correlation analysis helped to construct a comprehensive regulatory network of fructose metabolism in loquat, including 23 transcription factors, six structural genes, and nine saccharides. Based on the regulatory network and existing studies, it could be inferred that transcription factors such as ERF, NAC, MYB, GRAS, and bZIP may promote fructose accumulation in loquat flesh by positively regulating S6PDH. These findings improve our understanding of the nutritional value and breeding potential of white-fleshed loquat bud sport mutant, as well as serve as a foundation for exploring the genes and transcription factors that regulate fructose metabolism in loquat.

Diversity, origin, and evolution of the ESCRT systems.

Endosomal sorting complexes required for transport (ESCRT) play key roles in protein sorting between membrane-bounded compartments of eukaryotic cells. Homologs of many ESCRT components are identifiable in various groups of archaea, especially in Asgardarchaeota, the archaeal phylum that is currently considered to include the closest relatives of eukaryotes, but not in bacteria. We performed a comprehensive search for ESCRT protein homologs in archaea and reconstructed ESCRT evolution using the phylogenetic tree of Vps4 ATPase (ESCRT IV) as a scaffold and using sensitive protein sequence analysis and comparison of structural models to identify previously unknown ESCRT proteins. Several distinct groups of ESCRT systems in archaea outside of Asgard were identified, including proteins structurally similar to ESCRT-I and ESCRT-II, and several other domains involved in protein sorting in eukaryotes, suggesting an early origin of these components. Additionally, distant homologs of CdvA proteins were identified in Thermoproteales which are likely components of the uncharacterized cell division system in these archaea. We propose an evolutionary scenario for the origin of eukaryotic and Asgard ESCRT complexes from ancestral building blocks, namely, the Vps4 ATPase, ESCRT-III components, wH (winged helix-turn-helix fold) and possibly also coiled-coil, and Vps28-like domains. The last archaeal common ancestor likely encompassed a complex ESCRT system that was involved in protein sorting. Subsequent evolution involved either simplification, as in the TACK superphylum, where ESCRT was co-opted for cell division, or complexification as in Asgardarchaeota. In Asgardarchaeota, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was already established.IMPORTANCEAll eukaryotic cells possess complex intracellular membrane organization. Endosomal sorting complexes required for transport (ESCRT) play a central role in membrane remodeling which is essential for cellular functionality in eukaryotes. Recently, it has been shown that Asgard archaea, the archaeal phylum that includes the closest known relatives of eukaryotes, encode homologs of many components of the ESCRT systems. We employed protein sequence and structure comparisons to reconstruct the evolution of ESCRT systems in archaea and identified several previously unknown homologs of ESCRT subunits, some of which can be predicted to participate in cell division. The results of this reconstruction indicate that the last archaeal common ancestor already encoded a complex ESCRT system that was involved in protein sorting. In Asgard archaea, ESCRT systems evolved toward greater complexity, and in particular, the connection between ESCRT and the ubiquitin system that was previously considered a eukaryotic signature was established.

MMV method: a new approach to compare protein sequences under binary representation

In the present work, a new form of descriptor using minimal moment vector (MMV) is introduced to compare protein sequences in the frequency domain under their component wise binary representations. From every sequence, 20 different binary component sequences are formed, each corresponding to 20 amino acids. Each such vector is now shifted from the time domain to the frequency domain by applying the Fast Fourier Transform (FFT). Next, the power spectrum calculated from the FFT values for each component sequence is so normalized that the sum of the components equals 1. The descriptor is defined as a 20-component vector composed of the 20 second-order minimal moments calculated from the normalized spectrum of the 20 component sequences. Once the descriptor is known, the distance matrix is created by applying the Euclidean Distance measure. The phylogenetic tree is generated by applying the unweighted pair group method with the arithmetic mean (UPGMA) algorithm using Molecular Evolutionary Genetics Analysis11 (MEGA11) software. In this work, the datasets used for similarity studies are 9 NADH dehydrogenase 5 (ND5), 12 Baculoviruses, 24 Transferrins (TF) proteins, and 50 Spike Protein of coronavirus. A qualitative measure using rationalized perception is used to compare the effectiveness of the proposed method. Quantitative measure based on symmetric distance (SD) is used to compare the phylogenetic trees of the present method with those obtained by other methods. It is observed that the phylogenetic trees generated by the proposed technique are at par with their known biological references, and they produce results better than those of the earlier methods. Communicated by Ramaswamy H. Sarma

Comprehensive amino acid composition analysis of seed storage proteins of cereals and legumes: identification and understanding of intrinsically disordered and allergenic peptides

The seed storage proteins of cereal and legumes are the primary source of amino acids which are required for sustaining the nitrogen and carbon demands during germination and growth. Humans derive most of their dietary proteins from storage proteins in form of a wide variety of foods, for consumption. The amino acid content of most of these proteins is biased and the need for this biasness is not understood. The high abundance of proline, glutamine, and cysteine in cereals makes the gluten fraction viscoelastic. The cereal proteins have less charge and legume proteins have more charge on them. Their non-polar amino acid distribution has large variations. These characteristics are strongly responsible for the partial and complete unfolding of several domains of the storage proteins. Many of the storage proteins share a highly conserved structural feature within the cupin superfamily spread across all kingdoms of life. The intrinsically disordered viscoelastic proteins help in making dough which is vital for the quality of bread. Unfolded regions harbor more immunogenic sequences and cause food-related allergies and intolerance. We have discussed these properties in terms of comparison of cereal and legume storage protein sequences and allergy. Our study supports the findings that large disordered regions contain allergen-representative peptides. Interestingly, a high number of allergen-representative peptides were cleavable by digestive enzymes. Furthermore, unfolded storage proteins mimic microbial immunogens to induce a memory immune response. Results findings can be used to guide the understanding of immunological characteristics of storage proteins and may assist in treatment decisions for food allergy. Communicated by Ramaswamy H. Sarma

Cloning, Identification, and Functional Analysis of the Foxl2 Gene in Procambarus clarkii.

Procambarus clarkii is the most widely distributed freshwater shrimp in China, with important economic value and great potential for development. The forkheadboxL2 (Foxl2) gene has been found to be involved in the reproductive development of many crustaceans. To understand the role of the Foxl2 gene in the gonad development of P. clarkii, we designed CDS-specific primers for the P. clarkii Foxl2 (PcFoxl2) gene and cloned its CDS sequence using RT-PCR. The nucleotide and protein sequence information was then analyzed through bioinformatics analysis. The expression and subcellular localization of PcFoxl2 in various tissues were detected using qRT-PCR and in situ hybridization. The effects of PcFoxl2 knockdown on gonad development were investigated using RNA interference. The results showed that the CDS length of the PcFoxl2 gene was 1614 bp and encoded 537 amino acids. Protein sequence comparison and phylogenetic analysis showed that PcFoxl2 was the closest relative to Crayfish. qRT-PCR analysis indicated that the expression level of PcFoxl2 in the testis was significantly higher (>40 fold) than that in the ovary (p < 0.01). The in situ hybridization results showed that PcFoxl2 was expressed in both the cytoplasm and the nucleus of egg cells, and that the expression was strongest in egg cells at the early stage of yolk synthesis, while weak in the secondary oocytes. The positive signal was strongest in the spermatocyte nucleolus, while only a trace signal was observed in the cytoplasm. After interfering with the PcFoxl2 gene using dsRNA, the expression of PcFoxl2 in the RNA interference group was significantly lower than that in the control group, and this interference effect lasted for one week. Moreover, the gonad index of the experimental group was significantly lower than that of the control group (p < 0.05) after 10 days of P. clarkii cultivation following PcFoxl2 knockdown. The expression levels of the nanos and S3a genes, which are related to gonad development, decreased significantly after PcFoxl2 gene interference. The results suggest that the Foxl2 gene is involved in the growth and development of gonads, particularly in the development of testis, and is related to the early development of oocytes. This study provides a theoretical basis for the artificial breeding of P. clarkii.

Molecular characterization and allergenicity assessment of different samples of Mung Bean

Legumes represent a promising nutritional alternative source of proteins to meat and dairy products. Additionally, Novel Foods (Regulation EU 2015/2283) can help meet the rising protein demand. However, despite their benefits, emerging allergenicity risks must be considered. The aim of this work is the molecular characterization of the Novel Food Mung bean protein isolate for allergenicity prediction with High Resolution Mass Spectrometry analysis. The assessment of the allergenicity was evaluated in silico by comparing protein sequences of the Novel Food with other known legume allergens, using bioinformatic databases. The results highlighted similarity higher than 60 % of the protein structure of Mung bean with two known allergens of soybean and pea. Furthermore, enzymatic hydrolysis effects on allergenic potential was evaluated by immunoblotting analysis using sera of patients allergic to legumes. The protein hydrolysates obtained showed a high nutritional quality and a reduced allergenic potential, making them suitable for hypoallergenic food formulations.

Use of 2D FFT and DTW in Protein Sequence Comparison.

Protein sequence comparison remains a challenging work for the researchers owing to the computational complexity due to the presence of 20 amino acids compared with only four nucleotides in Genome sequences. Further, protein sequences of different species are of different lengths; it throws additional changes to the researchers to develop methods, specially alignment-free methods, to compare protein sequences. In this work, an efficient technique to compare protein sequences is developed by a graphical representation. First, the classified grouping of 20 amino acids with a cardinality of 4 based on polar class is considered to narrow down the representational range from 20 to 4. Then a unit vector technique based on a two-quadrant Cartesian system is proposed to provide a new two-dimensional graphical representation of the protein sequence. Now, two approaches are proposed to cope with the varying lengths of protein sequences from various species: one uses Dynamic Time Warping (DTW), while the other one uses a two-dimensional Fast Fourier Transform (2D FFT). Next, the effectiveness of these two techniques is analyzed using two evaluation criteria-quantitative measures based on symmetric distance (SD) and computational speed. An analysis is performed on five data sets of 9 ND4, 9 ND5, 9 ND6, 12 Baculovirus, and 24 TF proteins under the two methods. It is found that the FFT-based method produces the same results as DTW but in less computational time. It is found that the result of the proposed method agrees with the known biological reference. Further, the present method produces better clustering than the existing ones.

Protein sequence comparison based on representation on a finite dimensional unit hypercube

Numerous techniques are used to compare protein sequences based on the values of the physiochemical properties of amino acids. In this work, a single physical/chemical property value based non-binary representation of protein sequences is obtained on a 20 × 20-dimensional unit hypercube. The represented vector expressed in the matrix form is taken as the descriptor. The generalized NTV metric, which is an extension of the NTV metric used for polynucleotide space is taken as a distance measure. Based on this distance measure, a distance matrix is obtained for protein sequence comparison. Using this distance matrix, phylogenetic trees are drawn by using Molecular Evolutionary Genetics Analysis 11 (MEGA11) software applying the neighbor-joining method. Data sets used in this current work are 9-ND4, 9-ND5, 9-ND6, 24 TF-LF proteins, 27 different viruses and 127 proteins from the protein kinase C (PKC) family. Two sets of phylogenetic trees are obtained – one based on property value of polarity and the other based on property value of molecular weight. They are found to be exactly the same. Similar results also hold for other single property value based representation. The present trees are individually tested for efficiency based on the criterion of rationalized perception and computational time. The results of the present method are compared with those obtained earlier by other methods on the same protein sequences using assessment criteria of Symmetric distance (SD), Correlation coefficient, and Rationalized perception. In all the cases, the present results are found to be better than the results of other methods under comparison. Communicated by Ramaswamy H. Sarma

ATP binding cassette transporters and uridine diphosphate glycosyltransferases are ancient protein families that evolved roles in herbicide resistance through exaptation.

ATP-binding cassette (ABC) transporters actively transport various substances across membranes, while uridine diphosphate (UDP) glycosyltransferases (UGTs) are proteins that catalyse the chemical modification of various organic compounds. Both of these protein superfamilies have been associated with conferring herbicide resistance in weeds. Little is known about the evolutionary history of these protein families in the Archaeplastida. To infer the evolutionary histories of these protein superfamilies, we compared protein sequences collected from 10 species which represent distinct lineages of the Archaeplastida-the lineage including glaucophyte algae, rhodophyte algae, chlorophyte algae and the streptophytes-and generated phylogenetic trees. We show that ABC transporters were present in the last common ancestor of the Archaeplastida which lived 1.6 billion years ago, and the major clades identified in extant plants were already present then. Conversely, we only identified UGTs in members of the streptophyte lineage, which suggests a loss of these proteins in earlier diverging Archaeplastida lineages or arrival of UGTs into a common ancestor of the streptophyte lineage through horizontal gene transfer from a non-Archaeplastida eukaryote lineage. We found that within the streptophyte lineage, most diversification of the UGT protein family occurred in the vascular lineage, with 17 of the 20 clades identified in extant plants present only in vascular plants. Based on our findings, we conclude that ABC transporters and UGTs are ancient protein families which diversified during Archaeplastida evolution, which may have evolved for developmental functions as plants began to occupy new environmental niches and are now being selected to confer resistance to a diverse range of herbicides in weeds.

A wild boar cathelicidin peptide derivative inhibits severe acute respiratory syndrome coronavirus-2 and its drifted variants

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a clear threat to humanity. It has infected over 200 million and killed 4 million people worldwide, and infections continue with no end in sight. To control the pandemic, multiple effective vaccines have been developed, and global vaccinations are in progress. However, the virus continues to mutate. Even when full vaccine coverage is achieved, vaccine-resistant mutants will likely emerge, thus requiring new annual vaccines against drifted variants analogous to influenza. A complimentary solution to this problem could be developing antiviral drugs that inhibit SARS CoV-2 and its drifted variants. Host defense peptides represent a potential source for such an antiviral as they possess broad antimicrobial activity and significant diversity across species. We screened the cathelicidin family of peptides from 16 different species for antiviral activity and identified a wild boar peptide derivative that inhibits SARS CoV-2. This peptide, which we named Yongshi and means warrior in Mandarin, acts as a viral entry inhibitor. Following the binding of SARS-CoV-2 to its receptor, the spike protein is cleaved, and heptad repeats 1 and 2 multimerize to form the fusion complex that enables the virion to enter the cell. A deep learning-based protein sequence comparison algorithm and molecular modeling suggest that Yongshi acts as a mimetic to the heptad repeats of the virus, thereby disrupting the fusion process. Experimental data confirm the binding of Yongshi to the heptad repeat 1 with a fourfold higher affinity than heptad repeat 2 of SARS-CoV-2. Yongshi also binds to the heptad repeat 1 of SARS-CoV-1 and MERS-CoV. Interestingly, it inhibits all drifted variants of SARS CoV-2 that we tested, including the alpha, beta, gamma, delta, kappa and omicron variants.

Differences in amino acid composition of carrot α-tubulin potentially confer the resistance to dinitroaniline herbicides

Aim. To reveal the features of amino acid composition of carrot α-tubulin isotypes that potentially determine natural tolerance to dinitroaniline herbicides. Methods. Literature and database search. Comparison of protein sequences and structures: multiple sequence alignment, phylogenetic profiling, protein and ligand structure modeling, etc. Results. Genomic and proteomic analysis of Daucus carota has revealed at least eight unique isotypes of α-tubulin that differ in amino acid sequences and gene loci. Remarkable differences in amino acid composition of the dinitroanilinebinding-like (DBL) region of analyzed α-tubulin have been revealed, which may be the reason of its natural resistance to these compounds. Сonclusions. Differences in amino acids at positions of canonical mutations – Cys4 (TBA1, 2, 3, 6, 7 and 8), Thr53 (TBA6), Ile202 (TBA1 and 7) and Met202 (TBA5), as well as previously undescribed non-canonical substitutions – Ile4 (TBA4 and 5), Cys52 (TBA6), Ser201 (TBA1, 2, 3 and 8) and Val194 (TBA4 and 5), were noted as potentially associated with natural tolerance of the carrot to dinitroaniline herbicides.

Identification of a Novel ABC-2 Family Transporter Protein in Bacillus megaterium DSM319 Genome using Bioinformatics

Bacillus megaterium is a gram positive rod shaped bacterium that is greatly responsible for the formation of spores aerobically and found throughout the world. The size of respective bacterium genome is 5.48 Mb and contains about 5,858 exons. It has two strains named as QMB1551 and DSM319 having GC amount 38%. Regardless, a large number of the genes do not have any anticipated function, being depicted as hypothetical proteins. It is one of the most significant industrial bacterium that is employed for the generation of various useful products. It has played a promising role in the field of biotechnology as a production of genetically engineered proteins. Bioinformatics is one of the most advanced field that employed various productive tools in order to analyze the sequence, structure and function of proteins. In order to predict the function of a hypothetical protein of B. megaterium DSM319, BLAST was utilized for sequence comparison of hypothetical proteins of the bacterium dataset against the subject sequence of dataset of biotechnological important proteins. The best hit obtained from BLASTp results for sequence similarity was ABC-2 family transporter proteins having known functions. With that about 80% of similarity a daunorubicin ABC transporter permease (Anaerobacillus alkaliacustris) having query cover 100% and e-value 1 e-153 is obtained. The permease enzyme facilitates the ABC transporter to carry out export of daunorubicin that is produced by Bacillus spp. as a secondary metabolite. From the current study, it is predicted that B. megaterium having a protein known as daunorubicin ABC transporter permease. Moreover, this bacterium might be able to produce daunorubicin as a secondary metabolite. Daunorubicin has a great significance in medical sciences especially in anticancer therapy.

Proteome-wide comparison of tertiary protein structures reveals molecular mimicry in Plasmodium-human interactions

IntroductionMolecular mimicry is a strategy used by parasites to evade the host’s immune system and facilitate transmission to a new host. To date, high-throughput examples of molecular mimicry have been limited to comparing protein sequences. However, recent advances in the prediction of tertiary structural models, led by Deepmind’s AlphaFold, enable the comparison of thousands of proteins from parasites and their hosts at the structural level, allowing for the identification of more mimics. Here, we present the first proteome-level search for tertiary structure similarity between proteins from Plasmodium falciparum, a malaria-causing parasite, and humans.MethodsWe assembled a database of experimentally-characterized protein tertiary structures (from the Protein Data Bank) and AlphaFold-generated protein tertiary structures from P. falciparum, human, and 15 negative control species, i.e., species not infected by P. falciparum. We aligned human and control structures to the parasite structures using Foldseek. ResultsWe identified molecular mimicry in three proteins that have been previously proposed as mediators of Plasmodium-human interactions. By extending this approach to all P. falciparum proteins, we identified an additional 41 potential mimics that are supported by additional experimental data. DiscussionOur findings demonstrate a valuable application of AlphaFold-derived tertiary structural models, and we discuss key considerations for its effective use in other host-parasite systems.

Understanding Snail Mucus Biosynthesis and Shell Biomineralisation through Genomic Data Mining of the Reconstructed Carbohydrate and Glycan Metabolic Pathways of the Giant African Snail (Achatina fulica).

The giant African snail (Order Stylommatophora: Family Achatinidae), Achatina fulica (Bowdich, 1822), is the most significant and invasive land snail pest. The ecological adaptability of this snail involves high growth rate, reproductive capacity, and shell and mucus production, driven by several biochemical processes and metabolism. The available genomic information for A. fulica provides excellent opportunities to hinder the underlying processes of adaptation, mainly carbohydrate and glycan metabolic pathways toward the shell and mucus formation. The authors analysed the 1.78 Gb draft genomic contigs of A. fulica to identify enzyme-coding genes and reconstruct biochemical pathways related to the carbohydrate and glycan metabolism using a designed bioinformatic workflow. Three hundred and seventy-seven enzymes involved in the carbohydrate and glycan metabolic pathways were identified based on the KEGG pathway reference in combination with protein sequence comparison, structural analysis, and manual curation. Fourteen complete pathways of carbohydrate metabolism and seven complete pathways of glycan metabolism supported the nutrient acquisition and production of the mucus proteoglycans. Increased copy numbers of amylases, cellulases, and chitinases highlighted the snail advantage in food consumption and fast growth rate. The ascorbate biosynthesis pathway identified from the carbohydrate metabolic pathways of A. fulica was involved in the shell biomineralisation process in association with the collagen protein network, carbonic anhydrases, tyrosinases, and several ion transporters. Thus, our bioinformatic workflow was able to reconstruct carbohydrate metabolism, mucus biosynthesis, and shell biomineralisation pathways from the A. fulica genome and transcriptome data. These findings could reveal several evolutionary advantages of the A. fulica snail, and will benefit the discovery of valuable enzymes for industrial and medical applications.

E067 Four SNPs associated with monocyte/macrophage cell lineage uniquely associated with CRPS-1 in discovery and replication cohorts and suggest predisposition to regional osteopaenia and digit misperception

Abstract Background/Aims Do patients with CRPS1 have overexpressed SNPs in their exome sequencing? Do these SNP alleles correlate with phenotype? Methods We performed exome sequencing in a discovery cohort of well-characterised patients with CRPS using the Agilent 51M kit, sequencing to an average of 50-fold coverage. We used the vcf, bam and bam.bai files for genome wide SNP allele frequency assessment using the fSNPd programme. Deviations from expected frequency against 3 normative data sources were assessed for significance using a Chi-squared test with two tails and Yates correction. The resulting P values were subjected to a false discovery rate (FDR) correction as approximately 100,000 SNPs were assessed. Significant P-value SNPs were then filtered using bioinformatics, seeking SNPs where the alternative allele was predicted to have deleterious effects on the translated protein. We used the Human Genome Browser, NCBI BLASTP for protein sequence comparisons, Conserved Domains (CD search) for detecting if the amino acid change occurred within a known protein domain, and the SIFT and REVEL programs to predict if mis-sense mutations were deleterious. SNPs identified in the Discovery cohort were individually Sanger sequenced using genomic DNA in a Replication cohort. Patients with fibromyalgia were recruited as controls. Samples were collected from all the cohorts between 2014-19 after NRES approval (12/EE/40486; Cambridge Central). A subgroup of Cambridge patients were assessed for pain phenotype, including novel clinical signs, and radiographic osteopaenia (2 independent, blinded assessors (MSK Radiologist and Rheumatologist Consultants)). Results In total 5 SNPs were found to be overexpressed in the Discovery cohort (n = 37). All 5 SNPs were to be found to be functionally active and expressed in cells from the macrophage-monocyte lineage. 4 of these 5 SNPs were then found to be overrepresented in the Replication cohort (n = 48). The Fibromyalgia cohort did not overexpress these SNPs.24/84 (28.6%) patients expressed these SNPs. More males (8/15 (53.3%)) than females (16/69 (23.2%)) expressed SNPs. 7/27 patients from the phenotyped Cambridge cohort expressed these SNPs. 2/15 had osteopaenia as judged by 2 independent blinded assessors., both having the same SNP. 2/13 without osteopaenia expressed a SNP. There was a 1 in 16 chance of such clustering. Patients with SNPs had more depression and anxiety; less activity; and more likely to have digit misperception. There were no differences in functioning, pain severity, or examination findings. Conclusion Patients with CRPS-1 overexpress the rare allele of 4 SNPs. Remarkably, all 4 SNPs are expressed and functionally active in the same pathway within Monocyte-Macrophage cell lineages. Patients who carry these SNPs are more likely to be male and demonstrate radiographic osteopaenia, anxiety, depression, less activity and digit misperception. These new targets direct further basic and clinical research into pathophysiological mechanisms underpinning CRPS-type 1 and possibly other diseases. Disclosure N. Shenker: Honoraria; Lilly, Roche, Pfizer.

FAS: assessing the similarity between proteins using multi-layered feature architectures.

Protein sequence comparison is a fundamental element in the bioinformatics toolkit. When sequences are annotated with features such as functional domains, transmembrane domains, low complexity regions or secondary structure elements, the resulting feature architectures allow better informed comparisons. However, many existing schemes for scoring architecture similarities cannot cope with features arising from multiple annotation sources. Those that do fall short in the resolution of overlapping and redundant feature annotations. Here, we introduce FAS, a scoring method that integrates features from multiple annotation sources in a directed acyclic architecture graph. Redundancies are resolved as part of the architecture comparison by finding the paths through the graphs that maximize the pair-wise architecture similarity. In a large-scale evaluation on more than 10000 human-yeast ortholog pairs, architecture similarities assessed with FAS are consistently more plausible than those obtained using e-values to resolve overlaps or leaving overlaps unresolved. Three case studies demonstrate the utility of FAS on architecture comparison tasks: benchmarking of orthology assignment software, identification of functionally diverged orthologs, and diagnosing protein architecture changes stemming from faulty gene predictions. With the help of FAS, feature architecture comparisons can now be routinely integrated into these and many other applications. FAS is available as python package: https://pypi.org/project/greedyFAS/.

Protein Sequence Comparison Method Based on 3-ary Huffman Coding

Based on 3-ary Huffman coding algorithm, we propose a digital mapping method of protein sequence. Firstly, a 3-ary Huffman tree is defined by the frequency characteristic of 20 amino acids in given protein sequences. The 0-2 codes of 20 amino acids constructed by the 3-ary Huffman tree can convert long protein sequences into one-to-one 0-2 digital sequences. According to the frequency characteristic and the distribution information of 0-2 codes of 20 amino acids in the 0-2 digital sequences, we design the 40-dimensional vectors to characterize the protein sequences. Next, the proposed digital mapping method is used to perform three separate applications, similarity comparison of nine ND6 proteins, evolutionary trend analysis of the 2009 pandemic Human influenza A (H1N1) viruses from January 2020 to June 2022, and the evolution analysis of 95 coronavirus genes. The results illustrate the utility of the proposed method.

Comparative analysis of SilA-laccase mediated degradation of ciprofloxacin, norfloxacin and ofloxacin and interpretation of the possible catalytic mechanism

Fluoroquinolones (FQs) are the most commonly used antimicrobial drugs and regardless of their advantages in the healthcare sector, the pollution of these antimicrobial drugs in the environment has big concerns about human and environmental health. The presence of these antibiotic drugs even at the lowest concentrations in the environment has resulted in the emergence and spread of antibiotic resistance. Hence, it is necessary to remediate these pollutants from the environment. Previously alkaline laccase (SilA) from Streptomyces ipomoeae has been demonstrated to show degrading potentials against two of the FQs, Ciprofloxacin (CIP) and Norfloxacin (NOR); however, the molecular mechanism was not elucidated in detail. In this study, we have analyzed the possible molecular catalytic mechanism of FQ degrading SilA-laccase for the degradation of the FQs, CIP, NOR and Ofloxacin (OFL) using three-dimensional protein structure modeling, molecular docking and molecular dynamic (MD) studies. The comparative protein sequence analysis revealed the presence of tetrapeptide conserved catalytic motif, His102-X-His104-Gly105. After evaluating the active site of the enzyme in depth using CDD, COACH and S-site tools, we have identified the catalytic triad composed of three conserved amino acid residues, His102, Val103 and Tyr108 with which ligands interacted during the catalysis process. By analyzing the MD trajectories, it is revealed that the highest degradation potential of SilA is for CIP followed by NOR and OFL. Ultimately, this study provides the possible comparative catalytic mechanism for the degradation of CIP, NOR and OFL by the SilA enzyme. Communicated by Ramaswamy H. Sarma

Identification and genetic characterization of a new totivirus from Bursera graveolens in western Ecuador.

The complete genomic sequence of a previously uncharacterized virus provisionally named "Bursera graveolens associated totivirus 1" (BgTV-1) was obtained from Bursera graveolens (Kunth) Triana & Planch., a tree known as "palo santo" in Ecuador. The BgTV-1 genome is a monopartite double-stranded RNA (dsRNA) that is 4794 nucleotides (nt) long (GenBank accession number ON988291). Phylogenetic analysis of the capsid protein (CP) and RNA-dependent RNA polymerase (RdRp) placed BgTV-1 in a clade with other plant-associated totiviruses. Amino acid (aa) sequence comparisons of putative BgTV-1 proteins showed the highest sequence similarity to those of taro-associated totivirus L (QFS21890.1-QFS21891.1) and Panax notoginseng virus A (YP_009225664.1- YP_009225665.1), with 51.4% and 49.8% identity, respectively, in the CP and 56.4% and 55.2% identity, respectively, in the RdRp. BgTV-1 was not detected in total RNA from either of the two endophytic fungi cultured from BgTV-1-positive B. graveolens leaves, suggesting that BgTV-1 may be a plant-infecting totivirus. Based on its distinct host and the low aa sequence similarity between the CP of BgTV-1 and its counterparts from the closest relatives, the virus described in this study should be assigned as a new member of the genus Totivirus.

Demonstration of the Sequence Alignment to Predict Across Species Susceptibility Tool for Rapid Assessment of Protein Conservation.

The US Environmental Protection Agency Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool is a fast, freely available, online screening application that allows researchers and regulators to extrapolate toxicity information across species. For biological targets in model systems such as human cells, mice, rats, and zebrafish, toxicity data are available for a variety of chemicals. Through the evaluation of protein target conservation, this tool can be used to extrapolate data generated from such model systems to thousands of other species lacking toxicity data, yielding predictions of relative intrinsic chemical susceptibility. The latest releases of the tool (versions 2.0-6.1) have incorporated new features that allow for the rapid synthesis, interpretation, and use of the data for publication plus presentation-quality graphics. Among these features are customizable data visualizations and a comprehensive summary report designed to summarize SeqAPASS data for ease of interpretation. This paper describes the protocol to guide users through submitting jobs, navigating the various levels of protein sequence comparisons, and interpreting and displaying the resulting data. New features of SeqAPASS v2.0-6.0 are highlighted. Furthermore, two use-cases focused on transthyretin and opioid receptor protein conservation using this tool are described. Finally, SeqAPASS' strengths and limitations are discussed to define the domain of applicability for the tool and highlight different applications for cross-species extrapolation.