Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence analysis”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Impact-Driven Sampling Strategies for Hybrid Attack Graphs

Cyber-Physical Systems (CPSs) have a large input space, with discrete and continuous elements across multiple layers. Hybrid Attack Graph (HAG) provide a flexible and efficient approach to generate attack sequences for a CPS. Analysis and testing of large-scale HAGs are prohibitively costly. We propose a dimension reduction via property-preserving multi-layer graph sampling algorithms. Existing property-preserving graph sampling approaches generate a representative subgraph of an original large-sized graph while preserving the key properties, such as node and edge distribution, clustering coefficients, and betweenness. On the other hand, we propose impact-driven sampling strategies to transform the input data to a lower-dimensional representation while retaining key properties of the data.

Subasi, Omer↗
DOE OSTI2017

Identification of chalcone isomerase in the basal land plants reveals an ancient evolution of enzymatic cyclization activity for synthesis of flavonoids

Flavonoids ubiquitously distribute to the terrestrial plants and chalcone isomerase (CHI)- catalyzed intramolecular and stereospecific cyclization of chalcones is a committed step in the production of flavonoids. However, so far the bona fide CHIs are found only in vascular plants, and their origin and evolution remains elusive. We conducted transcriptomic and/or genomic sequence search, subsequent phylogenetic analysis, and detailed biochemical and genetic characterization to explore the potential existence of CHI proteins in the basal bryophyte liverwort species and the lycophyte Selaginella moellendorffii. We found that both liverwort and Selaginella species possess canonical CHI-fold proteins that cluster with their corresponding higher plant counterparts. Among them, some members exhibited bona fide CHI activity, which catalyze stereospecific cyclization of both 60- hydroxychalcone and 60-deoxychalcone, yielding corresponding 5-hydroxy and 5- deoxyflavanones, resembling the typical type II CHIs currently known to be ‘specific’ for legume plants. Expressing those primitive bona fide CHIs in the Arabidopsis chi mutant restores the seed coat transparent testa phenotype and the accumulation of flavonoids. These findings, in contrast to our current understanding of the evolution of enzymatic CHIs, suggest that emergence of the bona fide type II CHIs is an ancient evolution event that occurred before the divergence of liverwort lineages.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

diBELLA: Distributed Long Read to Long Read Alignment

© 2019 ACM. We present a parallel algorithm and scalable implementation for genome analysis, specifically the problem of finding overlaps and alignments for data from "third generation" long read sequencers [29]. While long sequences of DNA offer enormous advantages for biological analysis and insight, current long read sequencing instruments have high error rates and therefore require different approaches to analysis than their short read counterparts. Our work focuses on an efficient distributed-memory parallelization of an accurate single-node algorithm for overlapping and aligning long reads. We achieve scalability of this irregular algorithm by addressing the competing issues of increasing parallelism, minimizing communication, constraining the memory footprint, and ensuring good load balance. The resulting application, diBELLA, is the first distributed memory overlapper and aligner specifically designed for long reads and parallel scalability. We describe and present analyses for high level design trade-offs and conduct an extensive empirical analysis that compares performance characteristics across state-of-the-art HPC systems as well as a commercial cloud architectures, highlighting the advantages of state-of-the-art network technologies.

Ellis, Marquita↗
DOE OSTI2026

CAHS: Context-Aware Homology Search

Protein homology search is foundational to bioinformatics: it supports annotation transfer, structure/function inference, and evolutionary analysis over rapidly expanding sequence repositories (e.g., UniProtKB). Profile hidden Markov models (pHMMs), as implemented in HMMER, remain the most widely trusted approach because they provide statistically calibrated E-values; however, their gap behavior is fixed once a profile is trained, despite biological evidence that insertion/deletion tolerance varies across flexible loops and intrinsically disordered regions. We present CAHS (Context-Aware Homology Search), a lightweight query-time adapter for pHMM search that incorporates learned and biologically motivated signals without changing HMMER's downstream search pipeline or its calibrated E-value reporting. Given a query sequence, CAHS computes per-residue representations from a protein language model and a disorder predictor, maps these to profile coordinates, and modulates only match-state transition rows (gap-open and gap-extension probabilities) while preserving Plan7 constraints. We comprehensively evaluate CAHS across six structurally diverse protein families and multi-domain architectures against a 570k-sequence target corpus. CAHS expands detection capability, retrieving thousands of additional remote homologs at relaxed thresholds by maintaining alignment quality through flexible regions. For multi-domain proteins, context-aware modulation resolves 94% of fragmented alignments. Crucially, CAHS preserves hit-set invariance at stringent operating points (E<10-10), demonstrating increased statistical confidence without inflating false positives. Furthermore, sharper statistical distinction between homologs and background noise during early filter stages yields up to a 3.87× acceleration in end-to-end wall-clock time on high-performance computing clusters. Overall, CAHS illustrates a practical AI-for-science design pattern: augmenting a trusted probabilistic model with query-specific learned signals to improve interpretable, reproducible inference in data-rich biology.

Bhattaram, Swethasree [Georgia Institute of Techno↗
DOE OSTI2020

Distributed Berkeley Efficient Long-Read to Long-Read Aligner and Overlapper (DiBELLA) v1.0.0

We present a parallel algorithm and scalable implementation for genome analysis, specifically the problem of finding overlaps and alignments for data from "third generation" long read sequencers. While long sequences of DNA offer enormous advantages for biological analysis and insight, current long read sequencing instruments have high error rates and therefore require different approaches to analysis than their short read counterparts. Our work focuses on an efficient distributed-memory parallelization of an accurate single-node algorithm for overlapping and aligning long reads. We achieve scalability of this irregular algorithm by addressing the competing issues of increasing parallelism, minimizing communication, constraining the memory footprint, and ensuring good load balance. The resulting application, DiBELLA, is the first distributed memory overlapper and aligner specifically designed for long reads and parallel scalability.

Ellis, Marquita↗
DOE OSTI2025

GenomeDepot v1.0

GenomeDepot is a web-based platform for annotation, management, and comparative analysis of microbial genomic sequences and associated data including ortholog families, protein domains, operons, regulatory interactions, strain taxonomy, and sample metadata. GenomeDepot supports rapid creation of web-sites for user-defined genome collections that include bioinformatic tools for interactive genome browsing, BLAST search, annotation search, comparative genomic neighborhood visualization, and sequence download. Gene function annotations are generated by a customizable annotation pipeline. The pipeline runs annotation tools in Conda environments and can be easily extended with additional user-specified tools.

Kazakov, Alexey [Lawrence Berkeley National Labora↗
DOE OSTI2022

Structural variant analysis of a cancer reference cell line sample using multiple sequencing technologies

The cancer genome is commonly altered with thousands of structural rearrangements including insertions, deletions, translocation, inversions, duplications, and copy number variations. Thus, structural variant (SV) characterization plays a paramount role in cancer target identification, oncology diagnostics, and personalized medicine. As part of the SEQC2 Consortium effort, the present study established and evaluated a consensus SV call set using a breast cancer reference cell line and matched normal control derived from the same donor, which were used in our companion benchmarking studies as reference samples. We systematically investigated somatic SVs in the reference cancer cell line by comparing to a matched normal cell line using multiple NGS platforms including Illumina short-read, 10X Genomics linked reads, PacBio long reads, Oxford Nanopore long reads, and high-throughput chromosome conformation capture (Hi-C). We established a consensus SV call set of a total of 1788 SVs including 717 deletions, 230 duplications, 551 insertions, 133 inversions, 146 translocations, and 11 breakends for the reference cancer cell line. To independently evaluate and cross-validate the accuracy of our consensus SV call set, we used orthogonal methods including PCR-based validation, Affymetrix arrays, Bionano optical mapping, and identification of fusion genes detected from RNA-seq. We evaluated the strengths and weaknesses of each NGS technology for SV determination, and our findings provide an actionable guide to improve cancer genome SV detection sensitivity and accuracy. A high-confidence consensus SV call set was established for the reference cancer cell line. A large subset of the variants identified was validated by multiple orthogonal methods.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Small and mighty: adaptation of superphylum Patescibacteria to groundwater environment drives their genome simplicity

The newly defined superphylum Patescibacteria such as Parcubacteria (OD1) and Microgenomates (OP11) has been found to be prevalent in groundwater, sediment, lake, and other aquifer environments. Recently increasing attention has been paid to this diverse superphylum including > 20 candidate phyla (a large part of the candidate phylum radiation, CPR) because it refreshed our view of the tree of life. However, adaptive traits contributing to its prevalence are still not well known. Results Here, we investigated the genomic features and metabolic pathways of Patescibacteria in groundwater through genome-resolved metagenomics analysis of >600 Gbp sequence data. We observed that, while the members of Patescibacteria have reduced genomes (~ 1 Mbp) exclusively, functions essential to growth and reproduction such as genetic information processing were retained. Surprisingly, they have sharply reduced redundant and nonessential functions, including specific metabolic activities and stress response systems. The Patescibacteria have ultra-small cells and simplified membrane structures, including flagellar assembly, transporters, and two-component systems. Despite the lack of CRISPR viral defense, the bacteria may evade predation through deletion of common membrane phage receptors and other alternative strategies, which may explain the low representation of prophage proteins in their genomes and lack of CRISPR. By establishing the linkages between bacterial features and the groundwater environmental conditions, our results provide important insights into the functions and evolution of this CPR group.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Effect of fiber reinforcement on the process parameters of injection molding manufacturing process

Use of fillers with the resin is common in injection molding manufacturing technique to enhance the stiffness and strength of the manufactured components. However, the addition of fillers changes the rheological and mechanical properties of the material, required adjustments in the process parameters used for the manufacturing. The present study aims to investigate this impact of fillers on the process parameters of the injection molding manufacturing process. Neat Acrylonitrile Butadiene Styrene (ABS) and ABS reinforced with short glass fibers are used as the materials to perform the comprehensive numerical analysis including complete process sequence, namely filling, packing, cooling and warpage. Parametric studies are conducted to optimize the manufacturing process for minimal warpage (or shrinkage) in the produced part. The selected process parameters for the parametric study include cooling time, packing pressure, and mold temperature, all of which directly impact the quality of the manufactured parts. The optimal values obtained for these parameters are compared to examine the effect fiber reinforcement has on the manufacturing process.

Garg, Nikhil↗
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 OSTI2016

Comparing residue clusters from thermophilic and mesophilic enzymes reveals adaptive mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. As a result, the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.

09 BIOMASS FUELS↗
DOE OSTI2016

Rapid discovery and functional characterization of terpene synthases from four endophytic xylariaceae

Endophytic fungi are ubiquitous plant endosymbionts that establish complex and poorly understood relationships with their host organisms. Many endophytic fungi are known to produce a wide spectrum of volatile organic compounds (VOCs) with potential energy applications, which have been described as "mycodiesel". Many of these mycodiesel hydrocarbons are terpenes, a chemically diverse class of compounds produced by many plants, fungi, and bacteria. Due to their high energy densities, terpenes, such as pinene and bisabolene, are actively being investigated as potential "drop-in" biofuels for replacing diesel and aviation fuel. In this study, we rapidly discovered and characterized 26 terpene synthases (TPSs) derived from four endophytic fungi known to produce mycodiesel hydrocarbons. The TPS genes were expressed in an E. coli strain harboring a heterologous mevalonate pathway designed to enhance terpene production, and their product profiles were determined using Solid Phase Micro-Extraction (SPME) and GC-MS. Lastly, out of the 26 TPS’s profiled, 12 TPS’s were functional, with the majority of them exhibiting both monoterpene and sesquiterpene synthase activity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Differential CLE peptide perception by plant receptors implicated from structural and functional analyses of TDIF-TDR interactions

Tracheary Element Differentiation Inhibitory Factor (TDIF) belongs to the family of post-translationally modified CLE (CLAVATA3/embryo surrounding region (ESR)-related) peptide hormones that control root growth and define the delicate balance between stem cell proliferation and differentiation in SAM (shoot apical meristem) or RAM (root apical meristem). In Arabidopsis, Tracheary Element Differentiation Inhibitory Factor Receptor (TDR) and its ligand TDIF signaling pathway is involved in the regulation of procambial cell proliferation and inhibiting its differentiation into xylem cells. Here we present the crystal structures of the extracellular domains (ECD) of TDR alone and in complex with its ligand TDIF resolved at 2.65 Åand 2.75 Å respectively. These structures provide insights about the ligand perception and specific interactions between the CLE peptides and their cognate receptors. Our in vitro biochemical studies indicate that the interactions between the ligands and the receptors at the C-terminal anchoring site provide conserved binding. While the binding interactions occurring at the N-terminal anchoring site dictate differential binding specificities between different ligands and receptors. Our studies will open different unknown avenues of TDR-TDIF signaling pathways that will enhance our knowledge in this field highlighting the receptor ligand interaction, receptor activation, signaling network, modes of action and will serve as a structure function relationship model between the ligand and the receptor for various similar leucine-rich repeat receptor-like kinases (LRR-RLKs).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Molecular basis of surface anchored protein A deficiency in the Staphylococcus aureus strain Wood 46

Protein A in Staphylococcus aureus is encoded by the spa (staphylococcal protein A) gene and binds to immunoglobulin (Ig). The S. aureus strain Wood 46 has been variously reported as protein A-deficient and/or spa negative and used as a control in animal models of staphylococcal infections. The results of this study indicate that Wood 46 has normal spa expression but transcribes very low levels of the srtA gene which encodes the sortase A (SrtA) enzyme. This is consistent with unique mutations in the srtA promoter. In this study, a low level of sortase A explains deficient anchoring of proteins with an LPXTG motif, such as protein A, fibrinogen-binding protein and fibronectin-binding proteins A and B on to the peptidoglycan cell wall. The activity of secreted protein A is an important consideration for use of Wood 46 in functional experiments and animal models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Crystal structure of a Thermus aquaticus diversity-generating retroelement variable protein

Diversity-generating retroelements (DGRs) are widely distributed in bacteria, archaea, and microbial viruses, and bring about unparalleled levels of sequence variation in target proteins. While DGR variable proteins share low sequence identity, the structures of several such proteins have revealed the C-type lectin (CLec)-fold as a conserved scaffold for accommodating massive sequence variation. This conservation has led to the suggestion that the CLec-fold may be useful in molecular surface display applications. Thermostability is an attractive feature in such applications, and thus we studied the variable protein of a DGR encoded by a prophage of the thermophile Thermus aquaticus. We report here the 2.8 Å resolution crystal structure of the variable protein from the T. aquaticus DGR, called TaqVP, and confirm that it has a CLec-fold. Remarkably, its variable region is nearly identical in structure to those of several other CLec-fold DGR variable proteins despite low sequence identity among these. TaqVP was found to be thermostable, which appears to be a property shared by several CLec-fold DGR variable proteins. These results provide impetus for the pursuit of the DGR variable protein CLec-fold in molecular display applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

X-ray crystal structure of the Escherichia coli RadD DNA repair protein bound to ADP reveals a novel zinc ribbon domain

Genome maintenance is an essential process in all cells. In prokaryotes, the RadD protein is important for survival under conditions that include DNA-damaging radiation. Precisely how RadD participates in genome maintenance remains unclear. Here we present a high-resolution X-ray crystal structure of ADP-bound Escherichia coli RadD, revealing a zinc-ribbon element that was not modelled in a previous RadD crystal structure. Insights into the mode of nucleotide binding and additional structure refinement afforded by the new RadD model will help to drive investigations into the activity of RadD as a genome stability and repair factor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Initiation of RNA Polymerization and Polymerase Encapsidation by a Small dsRNA Virus

During the replication cycle of double-stranded (ds) RNA viruses, the viral RNA-dependent RNA polymerase (RdRP) replicates and transcribes the viral genome from within the viral capsid. How the RdRP molecules are packaged within the virion and how they function within the confines of an intact capsid are intriguing questions with answers that most likely vary across the different dsRNA virus families. In this study, we have determined a 2.4 Å resolution structure of an RdRP from the human picobirnavirus (hPBV). In addition to the conserved polymerase fold, the hPBV RdRP possesses a highly flexible 24 amino acid loop structure located near the C-terminus of the protein that is inserted into its active site. In vitro RNA polymerization assays and site-directed mutagenesis showed that: (1) the hPBV RdRP is fully active using both ssRNA and dsRNA templates; (2) the insertion loop likely functions as an assembly platform for the priming nucleotide to allow de novo initiation; (3) RNA transcription by the hPBV RdRP proceeds in a semi-conservative manner; and (4) the preference of virus-specific RNA during transcription is dictated by the lower melting temperature associated with the terminal sequences. Co-expression of the hPBV RdRP and the capsid protein (CP) indicated that, under the conditions used, the RdRP could not be incorporated into the recombinant capsids in the absence of the viral genome. Additionally, the hPBV RdRP exhibited higher affinity towards the conserved 5’-terminal sequence of the viral RNA, suggesting that the RdRP molecules may be encapsidated through their specific binding to the viral RNAs during assembly.

59 BASIC BIOLOGICAL SCIENCES↗