Search NASA⌕ Search

SEARCH · Search NASA

Results for “Amino Acid Sequence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Matrix Metalloproteinases as Candidate Antigenic Determinants for Anti‐Tumor Autoantibodies in Human Ovarian Cancer: A Post Hoc Analysis

Circulating antibodies in patients with cancer can facilitate the identification of accessible epitopes on autoantigens expressed by tumors. To identify previously unrecognized protein targets in ovarian cancer, we computationally assessed a heptapeptide consensus motif (VPELGHE, flanked by two cysteine residues yielding a cyclic nonapeptide under oxidizing conditions) previously discovered via phage display-based epitope mapping of autoantibodies in patients. Eight proteins associated with ovarian cancer encompass amino acid sequences similar to the consensus motif and were, therefore, considered as candidate native autoantigens. Among these candidate targets, however, matrix metalloproteinase 14 (MMP14) demonstrates gene expression that is both high and negatively correlated with survival in ovarian cancer patient cohorts. MMP14 protein levels are also stable in tumor versus non-tumor tissues. Moreover, the corresponding heptapeptide mimic in MMP14 occurs within an α-helical secondary structural element observed in its catalytic domain. These findings demonstrate that a subset of patient-derived autoantibodies may interact with a previously unknown antigenic epitope found in MMP14 and other MMPs, thereby providing opportunities for the development of new targeted agents.

Biochemistry & Molecular Biology↗
DOE OSTI2022

Zwartia hollandica gen. nov., sp. nov., Jezberella montanilacus gen. nov., sp. nov. and Sheuella amnicola gen. nov., comb. nov., representing the environmental GKS98 (betIII) cluster

We present two strains affiliated with the GKS98 cluster. This phylogenetically defined cluster is representing abundant, mainly uncultured freshwater bacteria, which were observed by many cultivation-independent studies on the diversity of bacteria in various freshwater lakes and streams. Bacteria affiliated with the GKS98 cluster were detected by cultivation-independent methods in freshwater systems located in Europe, Asia, Africa and the Americas. The two strains, LF4-65 T (=CCUG 56422 T =DSM 107630 T ) and MWH-P2sevCIIIb T (=CCUG 56420 T =DSM 107629 T ), are aerobic chemoorganotrophs, both with genome sizes of 3.2 Mbp and G+C values of 52.4 and 51.0 mol%, respectively. Phylogenomic analyses based on concatenated amino acid sequences of 120 proteins suggest an affiliation of the two strains with the family Alcaligenaceae and revealed Orrella amnicola and Orrella marina (= Algicoccus marinus ) as being the closest related, previously described species. However, the calculated phylogenomic trees clearly suggest that the current genus Orrella represents a polyphyletic taxon. Based on the branching order in the phylogenomic trees, as well as the revealed phylogenetic distances and chemotaxonomic traits, we propose to establish the new genus Zwartia gen. nov. and the new species Z. hollandica sp. nov. to harbour strain LF4-65 T and the new genus Jezberella gen. nov. and the new species J. montanilacus sp. nov. to harbour strain MWH-P2sevCIIIb T . Furthermore, we propose the reclassification of the species Orrella amnicola in the new genus Sheuella gen. nov. The new genera Zwartia, Jezberella and Sheuella together represent taxonomically the GKS98 cluster.

Microbiology↗
DOE OSTI2021

Solid-Binding Proteins: Bridging Synthesis, Assembly, and Function in Hybrid and Hierarchical Materials Fabrication

There is considerable interest in the development of hybrid organic–inorganic materials because of the potential for harvesting the unique capabilities that each system has to offer. Proteins are an especially attractive organic component owing to the high amount of chemical information encoded in their amino acid sequence, their amenability to molecular and computational (re)design, and the many structures and functions they specify. Genetic installation of solid-binding peptides (SBPs) within protein frameworks affords control over the position and orientation of adhesive and morphogenetic segments, and a path toward predictive synthesis and assembly of functional materials and devices, all while harnessing the built-in properties of the host scaffold. Furthermore, we review the current understanding of the mechanisms through which SBPs bind to technologically relevant interfaces, with an emphasis on the variables that influence the process, and highlight the last decade of progress in the use of solid-binding proteins for hybrid and hierarchical materials synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Spatial top-down proteomics for the functional characterization of human kidney

Background: The Human Proteome Project has credibly detected nearly 93% of the roughly 20,000 proteins which are predicted by the human genome. However, the proteome is enigmatic, where alterations in amino acid sequences from polymorphisms and alternative splicing, errors in translation, and post-translational modifications result in a proteome depth estimated at several million unique proteoforms. Recently mass spectrometry has been demonstrated in several landmark efforts mapping the human proteoform landscape in bulk analyses. Herein, we developed an integrated workflow for characterizing proteoforms from human tissue in a spatially resolved manner by coupling laser capture microdissection, nanoliter-scale sample preparation, and mass spectrometry imaging. Results: Using healthy human kidney sections as the case study, we focused our analyses on the major functional tissue units including glomeruli, tubules, and medullary rays. After laser capture microdissection, these isolated functional tissue units were processed with microPOTS (microdroplet processing in one-pot for trace samples) for sensitive top-down proteomics measurement. This provided a quantitative database of 616 proteoforms that was further leveraged as a library for mass spectrometry imaging with near-cellular spatial resolution over the entire section. Notably, several mitochondrial proteoforms were found to be differentially abundant between glomeruli and convoluted tubules, and further spatial contextualization was provided by mass spectrometry imaging confirming unique differences identified by microPOTS, and further expanding the field-of-view for unique distributions such as enhanced abundance of a truncated form (1-74) of ubiquitin within cortical regions. Conclusions: We developed an integrated workflow to directly identify proteoforms and reveal their spatial distributions. Where of the 20 differentially abundant proteoforms identified as discriminate between tubules and glomeruli by microPOTS, the vast majority of tubular proteoforms were of mitochondrial origin (8 of 10) where discriminate proteoforms in glomeruli were primarily hemoglobin subunits (9 of 10). These trends were also identified within ion images demonstrating spatially resolved characterization of proteoforms that has the potential to reshape discovery-based proteomics because the proteoforms are the ultimate effector of cellular functions. Applications of this technology have the potential to unravel etiology and pathophysiology of disease states, informing on biologically active proteoforms, which remodel the proteomic landscape in chronic and acute disorders.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Real-time structural motif searching in proteins using an inverted index strategy

Biochemical and biological functions of proteins are the product of both the overall fold of the polypeptide chain, and, typically, structural motifs made up of smaller numbers of amino acids constituting a catalytic center or a binding site that may be remote from one another in amino acid sequence. Detection of such structural motifs can provide valuable insights into the function(s) of previously uncharacterized proteins. Technically, this remains an extremely challenging problem because of the size of the Protein Data Bank (PDB) archive. Existing methods depend on a clustering by sequence similarity and can be computationally slow. We have developed a new approach that uses an inverted index strategy capable of analyzing >170,000 PDB structures with unmatched speed. The efficiency of the inverted index method depends critically on identifying the small number of structures containing the query motif and ignoring most of the structures that are irrelevant. Our approach (implemented at motif.rcsb.org ) enables real-time retrieval and superposition of structural motifs, either extracted from a reference structure or uploaded by the user. Herein, we describe the method and present five case studies that exemplify its efficacy and speed for analyzing 3D structures of both proteins and nucleic acids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Crystal structure of the Arabidopsis SPIRAL2 C-terminal domain reveals a p80-Katanin-like domain

Epidermal cells of dark-grown plant seedlings reorient their cortical microtubule arrays in response to blue light from a net lateral orientation to a net longitudinal orientation with respect to the long axis of cells. The molecular mechanism underlying this microtubule array reorientation involves katanin, a microtubule severing enzyme, and a plant-specific microtubule associated protein called SPIRAL2. Katanin preferentially severs longitudinal microtubules, generating seeds that amplify the longitudinal array. Upon severing, SPIRAL2 binds nascent microtubule minus ends and limits their dynamics, thereby stabilizing the longitudinal array while the lateral array undergoes net depolymerization. To date, no experimental structural information is available for SPIRAL2 to help inform its mechanism. To gain insight into SPIRAL2 structure and function, we determined a 1.8 Å resolution crystal structure of the Arabidopsis thaliana SPIRAL2 C-terminal domain. The domain is composed of seven core α-helices, arranged in an α-solenoid. Amino-acid sequence conservation maps primarily to one face of the domain involving helices α1, α3, α5, and an extended loop, the α6-α7 loop. The domain fold is similar to, yet structurally distinct from the C-terminal domain of Ge-1 (an mRNA decapping complex factor involved in P-body localization) and, surprisingly, the C-terminal domain of the katanin p80 regulatory subunit. The katanin p80 C-terminal domain heterodimerizes with the MIT domain of the katanin p60 catalytic subunit, and in metazoans, binds the microtubule minus-end factors CAMSAP3 and ASPM. Structural analysis predicts that SPIRAL2 does not engage katanin p60 in a mode homologous to katanin p80. The SPIRAL2 structure highlights an interesting evolutionary convergence of domain architecture and microtubule minus-end localization between SPIRAL2 and katanin complexes, and establishes a foundation upon which structure-function analysis can be conducted to elucidate the role of this domain in the regulation of plant microtubule arrays.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Biophysical and biochemical evidence for the role of acetate kinases (AckAs) in an acetogenic pathway in pathogenic spirochetes

Unraveling the metabolism of Treponema pallidum is a key component to understanding the pathogenesis of the human disease that it causes, syphilis. For decades, it was assumed that glucose was the sole carbon/energy source for this parasitic spirochete. But the lack of citric-acid-cycle enzymes suggested that alternative sources could be utilized, especially in microaerophilic host environments where glycolysis should not be robust. Recent bioinformatic, biophysical, and biochemical evidence supports the existence of an acetogenic energy-conservation pathway in T . pallidum and related treponemal species. In this hypothetical pathway, exogenous D-lactate can be utilized by the bacterium as an alternative energy source. Herein, we examined the final enzyme in this pathway, acetate kinase (named TP0476), which ostensibly catalyzes the generation of ATP from ADP and acetyl-phosphate. We found that TP0476 was able to carry out this reaction, but the protein was not suitable for biophysical and structural characterization. We thus performed additional studies on the homologous enzyme (75% amino-acid sequence identity) from the oral pathogen Treponema vincentii , TV0924. This protein also exhibited acetate kinase activity, and it was amenable to structural and biophysical studies. We established that the enzyme exists as a dimer in solution, and then determined its crystal structure at a resolution of 1.36 Å, showing that the protein has a similar fold to other known acetate kinases. Mutation of residues in the putative active site drastically altered its enzymatic activity. A second crystal structure of TV0924 in the presence of AMP (at 1.3 Å resolution) provided insight into the binding of one of the enzyme’s substrates. On balance, this evidence strongly supported the roles of TP0476 and TV0924 as acetate kinases, reinforcing the hypothesis of an acetogenic pathway in pathogenic treponemes.

Deka, Ranjit K.↗
DOE OSTI2024

Chickpea NCR13 disulfide cross-linking variants exhibit profound differences in antifungal activity and modes of action

Small cysteine-rich antifungal peptides with multi-site modes of action (MoA) have potential for development as biofungicides. In particular, legumes of the inverted repeat-lacking clade express a large family of nodule-specific cysteine-rich (NCR) peptides that orchestrate differentiation of nitrogen-fixing bacteria into bacteroids. These NCRs can form two or three intramolecular disulfide bonds and a subset of these peptides with high cationicity exhibits antifungal activity. However, the importance of intramolecular disulfide pairing and MoA against fungal pathogens for most of these plant peptides remains to be elucidated. Our study focused on a highly cationic chickpea NCR13, which has a net charge of +8 and contains six cysteines capable of forming three disulfide bonds. NCR13 expression in Pichia pastoris resulted in formation of two peptide folding variants, NCR13_PFV1 and NCR13_PFV2, that differed in the pairing of two out of three disulfide bonds despite having an identical amino acid sequence. The NMR structure of each PFV revealed a unique three-dimensional fold with the PFV1 structure being more compact but less dynamic. Surprisingly, PFV1 and PFV2 differed profoundly in the potency of antifungal activity against several fungal plant pathogens and their multi-faceted MoA. PFV1 showed significantly faster fungal cell-permeabilizing and cell entry capabilities as well as greater stability once inside the fungal cells. Additionally, PFV1 was more effective in binding fungal ribosomal RNA and inhibiting protein translation in vitro. Furthermore, when sprayed on pepper and tomato plants, PFV1 was more effective in reducing disease symptoms caused by Botrytis cinerea, causal agent of gray mold disease in fruits, vegetables, and flowers. In conclusion, our work highlights the significant impact of disulfide pairing on the antifungal activity and MoA of NCR13 and provides a structural framework for design of novel, potent antifungal peptides for agricultural use.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

RCSB Protein Data Bank: visualizing groups of experimentally determined PDB structures alongside computed structure models of proteins

Recent advances in Artificial Intelligence and Machine Learning (e.g., AlphaFold, RosettaFold, and ESMFold) enable prediction of three-dimensional (3D) protein structures from amino acid sequences alone at accuracies comparable to lower-resolution experimental methods. These tools have been employed to predict structures across entire proteomes and the results of large-scale metagenomic sequence studies, yielding an exponential increase in available biomolecular 3D structural information. Given the enormous volume of this newly computed biostructure data, there is an urgent need for robust tools to manage, search, cluster, and visualize large collections of structures. Equally important is the capability to efficiently summarize and visualize metadata, biological/biochemical annotations, and structural features, particularly when working with vast numbers of protein structures of both experimental origin from the Protein Data Bank (PDB) and computationally-predicted models. Moreover, researchers require advanced visualization techniques that support interactive exploration of multiple sequences and structural alignments. This paper introduces a suite of tools provided on the RCSB PDB research-focused web portal RCSB. org, tailor-made for efficient management, search, organization, and visualization of this burgeoning corpus of 3D macromolecular structure data.

3D visualization↗
DOE OSTI2025

AlgaeOrtho, a bioinformatics tool for processing ortholog inference results in algae

Introduction: Microalgae constitute a prominent feedstock for producing biofuels and biochemicals by virtue of their prolific reproduction, high bioproduct accumulation, and the ability to grow in brackish and saline water. However, naturally occurring wild type algal strains are rarely optimal for industrial use; therefore, bioengineering of algae is necessary to generate superior performing strains that can address production challenges in industrial settings, particularly the bioenergy and bioproduct sectors. One of the crucial steps in this process is deciding on a bioengineering target: namely, which gene/protein to differentially express. These targets are often orthologs which are defined as genes/proteins originating from a common ancestor in divergent species. Although bioinformatics tools for the identification of protein orthologs already exist, processing the output from such tools is nontrivial, especially for a researcher with little or no bioinformatics experience. Methods: The present study introduces AlgaeOrtho, a user-friendly tool that builds upon the SonicParanoid orthology inference tool (based on an algorithm that identifies potential protein orthologs based on amino acid sequences) and the PhycoCosm database from JGI (Joint Genome Institute) to help researchers identify orthologs of their proteins of interest in multiple diverse algal species. Results: The output of this application includes a table of the putative orthologs of their protein of interest, a heatmap showing sequence similarity (%), and an unrooted tree of the putative protein orthologs. Notably, the tool would be instrumental in identifying novel bioengineering targets in different algal strains, including targets in not-fully annotated algal species, since it does not depend on existing protein annotations. We tested AlgaeOrtho using three case studies, for which orthologs of proteins relevant to bioengineering targets, were identified from diverse algal species, demonstrating its ease of use and utility for bioengineering researchers. Discussion: This tool is unique in the protein ortholog identification space as it can visualize putative orthologs, as desired by the user, across several algal species.

09 BIOMASS FUELS↗
DOE OSTI2022

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Two Alternative Splicing Variants of AtERF73/HRE1, HRE1α and HRE1β, Have Differential Transactivation Activities in Arabidopsis

AtERF73/HRE1 is an AP2/ERF transcription factor in Arabidopsis and has two distinct alternative splicing variants, HRE1α and HRE1β. In this study, we examined the differences between the molecular functions of HRE1α and HRE1β. We found that HRE1α and HRE1β are both involved in hypoxia response and root development and have transactivation activity. Two conserved motifs in the C-terminal region of HRE1α and HRE1β, EELL and LWSY-like, contributed to their transactivation activity, specifically the four E residues in the EELL motif and the MGLWS amino acid sequence at the end of the LWSY-like motif. The N-terminal region of HRE1β also showed transactivation activity, mediated by the VDDG motif, whereas that of HRE1α did not. The transactivation activity of HRE1β was stronger than that of HRE1α in Arabidopsis protoplasts. Both transcription factors transactivated downstream genes via the GCC box. RNA-sequencing analysis further supported that both HRE1α and HRE1β might regulate gene expression associated with the hypoxia stress response, although they may transactivate different subsets of genes in downstream pathways. Our results, together with previous studies, suggested that HRE1α and HRE1β differentially transactivate downstream genes in hypoxia response and root development in Arabidopsis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Characterization of Humanized Mouse Model of Organophosphate Poisoning and Detection of Countermeasures via MALDI-MSI

Organophosphoate (OP) chemicals are known to inhibit the enzyme acetylcholinesterase (AChE). Studying OP poisoning is difficult because common small animal research models have serum carboxylesterase, which contributes to animals’ resistance to OP poisoning. Historically, guinea pigs have been used for this research; however, a novel genetically modified mouse strain (KIKO) was developed with nonfunctional serum carboxylase (Es1 KO) and an altered acetylcholinesterase (AChE) gene, which expresses the amino acid sequence of the human form of the same protein (AChE KI). KIKO mice were injected with 1xLD50 of an OP nerve agent or vehicle control with or without atropine. After one to three minutes, animals were injected with 35 mg/kg of the currently fielded Reactivator countermeasure for OP poisoning. Postmortem brains were imaged on a Bruker RapifleX ToF/ToF instrument. Data confirmed the presence of increased acetylcholine in OP-exposed animals, regardless of treatment or atropine status. More interestingly, we detected a small amount of Reactivator within the brain of both exposed and unexposed animals; it is currently debated if reactivators can cross the blood–brain barrier. Further, we were able to simultaneously image acetylcholine, the primary affected neurotransmitter, as well as determine the location of both Reactivator and acetylcholine in the brain. This study, which utilized sensitive MALDI-MSI methods, characterized KIKO mice as a functional model for OP countermeasure development.

2-PAM↗
DOE OSTI2021

deeprob/ThioesteraseEnzymeSpecificity: EnZymClass-first-release

Characterizing the functional properties of plant acyl-ACP thioesterases (TEs), a key enzyme class used in the production of renewable oleochemicals in microbial hosts, experimentally, can be an expensive and time consuming process since it requires manual screening of thousands of candidates in a database. Using amino acid sequence to computationally predict an enzyme’s function might accelerate this process; however obtaining the necessary amount of information on previously characterized enzymes and their respective sequences required by standard Machine Learning (ML) based approaches to accurately infer sequence-function relationships can be prohibitive, especially with a low-throughput testing cycle. Experimental noise, unbalanced dataset where high sequence similarity does not always imply identical functional properties will further prevent robust prediction performance. Herein we present a ML method, Ensemble method for enZyme Classification (EnZymClass), that is specifically designed to address these issues. We used EnZymClass to classify TEs into short, long and mixed free fatty acid substrate specificity categories. While general guidelines for inferring substrate specificity have been proposed before, prediction of chain-length preference from primary sequence has remained elusive for plant acyl-ACP TEs. By applying EnZymClass to a subset of TEs in the ThYme database, we identified two medium chain TEs, ClFatB3 and CwFatB2, with previously uncharacterized activity in E. coli fatty acid production hosts.

Banerjee, Deepro↗
DOE OSTI2021

Designed proteins for pH switchable antibody purification

Modified Fc-binding domain polypeptides or polypeptides including one or more modified Fc-binding domains are provided. In some examples, a modified Fc-binding domain includes or consists of the amino acid sequence of any one of SEQ ID NOs: 4, 6, or 8. Nucleic acids encoding the modified Fc-binding domains or polypeptides including one or more of the modified Fc-binding domains are also provided. In some embodiments, the nucleic acids are included in a vector and may be operably linked to a promoter. Methods for purifying a polypeptide including one or more Fc regions are also provided. Kits and reagents including the modified Fc-binding domain polypeptides or polypeptides including one or more modified Fc-binding domains are provided. In some examples, the modified Fc-binding domain polypeptides or polypeptides including one or more modified Fc-binding domains are linked to a solid support.

Jha, Roshan Kumar↗
DOE OSTI2022

Evolution of tRNA synthetases

The disclosure provides amino acid sequence variants of orthogonal aminoacyl-tRNA synthetases (AARSs) having increased activity and selectivity compared to previous AARSs, and methods of producing the same.

Liu, David R.↗
DOE OSTI2020

Short-wavelength-sensitive 2 (Sws2) visual photopigment models combined with atomistic molecular simulations to predict spectral peaks of absorbance

For many species, vision is one of the most important sensory modalities for mediating essential tasks that include navigation, predation and foraging, predator avoidance, and numerous social behaviors. The vertebrate visual process begins when photons of the light interact with rod and cone photoreceptors that are present in the neural retina. Vertebrate visual photopigments are housed within these photoreceptor cells and are sensitive to a wide range of wavelengths that peak within the light spectrum, the latter of which is a function of the type of chromophore used and how it interacts with specific amino acid residues found within the opsin protein sequence. Minor differences in the amino acid sequences of the opsins are known to lead to large differences in the spectral peak of absorbance (i.e. the λ max value). In our prior studies, we developed a new approach that combined homology modeling and molecular dynamics simulations to gather structural information associated with chromophore conformation, then used it to generate statistical models for the accurate prediction of λ max values for photopigments derived from Rh1 and Rh2 amino acid sequences. In the present study, we test our novel approach to predict the λ max of phylogenetically distant Sws2 cone opsins. To build a model that can predict the λ max using our approach presented in our prior studies, we selected a spectrally-diverse set of 11 teleost Sws2 photopigments for which both amino acid sequence information and experimentally measured λ max values are known. The final first-order regression model, consisting of three terms associated with chromophore conformation, was sufficient to predict the λ max of Sws2 photopigments with high accuracy. This study further highlights the breadth of our approach in reliably predicting λ max values of Sws2 cone photopigments, evolutionary-more distant from template bovine RH1, and provided mechanistic insights into the role of known spectral tuning sites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Snekmer: a scalable pipeline for protein sequence fingerprinting based on amino acid recoding

Abstract Motivation The vast expansion of sequence data generated from single organisms and microbiomes has precipitated the need for faster and more sensitive methods to assess evolutionary and functional relationships between proteins. Representing proteins as sets of short peptide sequences (kmers) has been used for rapid, accurate classification of proteins into functional categories; however, this approach employs an exact-match methodology and thus may be limited in terms of sensitivity and coverage. We have previously used similarity groupings, based on the chemical properties of amino acids, to form reduced character sets and recode proteins. This amino acid recoding (AAR) approach simplifies the construction of protein representations in the form of kmer vectors, which can link sequences with distant sequence similarity and provide accurate classification of problematic protein families. Results Here, we describe Snekmer, a software tool for recoding proteins into AAR kmer vectors and performing either (i) construction of supervised classification models trained on input protein families or (ii) clustering for de novo determination of protein families. We provide examples of the operation of the tool against a set of nitrogen cycling families originally collected using both standard hidden Markov models and a larger set of proteins from Uniprot and demonstrate that our method accurately differentiates these sequences in both operation modes. Availability and implementation Snekmer is written in Python using Snakemake. Code and data used in this article, along with tutorial notebooks, are available at http://github.com/PNNL-CompBio/Snekmer under an open-source BSD-3 license. Supplementary information Supplementary data are available at Bioinformatics Advances online.

59 BASIC BIOLOGICAL SCIENCES↗