Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence Alignment”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Genome-wide Transcription Factor DNA Binding Sites and Gene Regulatory Networks in Clostridium thermocellum

Clostridium thermocellum is a thermophilic bacterium recognized for its natural ability to effectively deconstruct cellulosic biomass. While there is a large body of studies on the genetic engineering of this bacterium and its physiology to-date, there is limited knowledge in the transcriptional regulation in this organism and thermophilic bacteria in general. The study herein is the first report of a high-throughput application of DNA-affinity purification sequencing (DAP-seq) to transcription factors (TFs) from a thermophile. We applied DAP-seq to >90 TFs in C. thermocellum and detected genome-wide binding sites for 11 of them. We then compiled and aligned DNA binding sequences from these TFs to deduce the primary DNA-binding sequence motifs for each TF. These binding motifs are further validated with electrophoretic mobility shift assay (EMSA) and are used to identify individual TFs’ regulatory targets in C. thermocellum. Our results led to the discovery of novel, uncharacterized TFs as well as homologues of previously studied TFs including RexA-, LexA- and LacI-type TFs. We then used these data to reconstruct gene regulatory networks for the 11 TFs individually, which resulted in a global network encompassing the TFs with some interconnections. As gene regulation governs and constrains how bacteria behave, our findings shed light on the roles of TFs delineated by their regulons, and potentially provides a means to enable rational, advanced genetic engineering of C. thermocellum and other organisms alike towards a desired phenotype.

09 BIOMASS FUELS↗
DOE OSTI2023

Sample Type Adaptable RNA depletion performance data [Slides]

Sequencing data was aligned to rRNA data sets QIIME_16S_MiDAS_4.8.1 and SILVA_138.1_LSURef_NR99 using BWA. RNA used was a composite (pool) of wastewater RNA samples collected at LANL between April 2022 and December 2022. RNA sample was split into four aliquots (STAR depletion and STAR depletion no-probe-control as well as RiboZero and RiboZero no-probe-control). All sequencing libraries were prepared from the depleted RNA samples using the same library prep method and sequenced on Illumina platforms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

CoverM: read alignment statistics for metagenomics

SUMMARY: Genome-centric analysis of metagenomic samples is a powerful method for understanding the function of microbial communities. Calculating read coverage is a central part of analysis, enabling differential coverage binning for recovery of genomes and estimation of microbial community composition. Coverage is determined by processing read alignments to reference sequences of either contigs or genomes. Per-reference coverage is typically calculated in an ad-hoc manner, with each software package providing its own implementation and specific definition of coverage. Here we present a unified software package CoverM which calculates several coverage statistics for contigs and genomes in an ergonomic and flexible manner. It uses "Mosdepth arrays" for computational efficiency and avoids unnecessary I/O overhead by calculating coverage statistics from streamed read alignment results. AVAILABILITY AND IMPLEMENTATION: CoverM is free software available at https://github.com/wwood/coverm. CoverM is implemented in Rust, with Python (https://github.com/apcamargo/pycoverm) and Julia (https://github.com/JuliaBinaryWrappers/CoverM_jll.jl) interfaces.

Aroney, Samuel T N↗
DOE OSTI2012

The Crystal Structure of the Extracellular 11-heme Cytochrome UndA Reveals a Conserved 10-heme Motif and Defined Binding Site for Soluble Iron Chelates.

Members of the genus Shewanella translocate deca- or undeca-heme cytochromes to the external cell surface thus enabling respiration using extracellular minerals and polynuclear Fe(III) chelates. The high resolution structure of the first undeca-heme outer membrane cytochrome, UndA, reveals a crossed heme chain with four potential electron ingress/egress sites arranged within four domains. Sequence and structural alignment of UndA and the deca-heme MtrF reveals the extra heme of UndA is inserted between MtrF hemes 6 and 7. The remaining UndA hemes can be superposed over the heme chain of the decaheme MtrF, suggesting that a ten heme core is conserved between outer membrane cytochromes. The UndA structure is the first outer membrane cytochrome to be crystallographically resolved in complex with substrates, an Fe(III)-nitrilotriacetate dimer or an Fe(III)-citrate trimer. The structural resolution of these UndA-Fe(III)-chelate complexes provides a rationale for previous kinetic measurements on UndA and other outer membrane cytochromes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Identification of 2-Hydroxyacyl-CoA Synthases with High Acyloin Condensation Activity for Orthogonal One-Carbon Bioconversion

One-carbon (C1) compounds are emerging as cost-effective and potentially carbon-negative feedstocks for biomanufacturing, which require efficient, versatile metabolic platforms for the synthesis of value-added products. Synthetic formyl-CoA elongation (FORCE) pathways allow diverse product synthesis from C1 compounds via iterative C1 elongation, operating independently from the host metabolism with reduced engineering complexity and improved theoretical yields. However, a major bottleneck was identified as the suboptimal kinetics of the core C1–C1 condensation enzyme, 2-hydroxyacyl-CoA synthase (HACS), catalyzing the acyloin condensation reaction between formaldehyde and formyl-CoA. Furthermore, we used a combinatorial approach of bioprospecting and rational protein engineering to identify multiple HACS variants with significantly improved activities toward C1 substrates. Sequence and structure alignment of the active variants elucidated the key regions for the catalytic function, which were targeted for mutagenesis, leading to improved catalytic efficiency. In parallel, a consecutive round of bioprospecting for homologs with high similarity with active variants revealed a highly active HACS variant exhibiting up to 7-fold improvement in catalytic efficiency (k cat /K M ) and 14-fold improvement in the FORCE pathway flux in vivo compared to the previous reports. Upon further optimization of the downstream pathway, the orthogonal C1-to-product bioconversion system showed a metabolic flux of up to 700 μM glycolate OD –1 h –1 (2.1 mmol gDCW –1 h –1 ) and an industrially relevant glycolate titer, rate, and yield of 5.2 g L –1 (67.8 mM), 0.22 g L –1 h –1 , and 94% carbon yield, respectively.

2-hydroxyacyl-CoA synthase↗
DOE OSTI2018

Navigating the maze: the impact of configurability in bioinformatics software

The bioinformatics software domain contains thousands of applications for automating tasks such as the pairwise alignment of DNA sequences, building and reasoning about metabolic models or simulating growth of an organism. Its end users range from sophisticated developers to those with little computational experience. In response to their needs, developers provide many options to customize the way their algorithms are tuned. Yet there is little or no automated help for the user in determining the consequences or impact of the options they choose. In this paper we describe our experience working with configurable bioinformatics tools. We find limited documentation and help for combining and selecting options along with variation in both functionality and performance. We also find previously undetected faults. We summarize our findings with a set of lessons learned, and present a roadmap for creating automated techniques to interact with bioinformatics software. We believe these will generalize to other types of scientific software.

McDevitt, Mikaela Cashman↗
DOE OSTI2018

Classification of Complete Proteomes of Different Organisms and Protein Sets Based on Their Protein Distributions in Terms of Some Key Attributes of Proteins

The existence of complete genome sequences makes it important to develop different approaches for classification of large-scale data sets and to make extraction of biological insights easier. Here, we propose an approach for classification of complete proteomes/protein sets based on protein distributions on some basic attributes. We demonstrate the usefulness of this approach by determining protein distributions in terms of two attributes: protein lengths and protein intrinsic disorder contents (ID). The protein distributions based on L and ID are surveyed for representative proteome organisms and protein sets from the three domains of life. The two-dimensional maps (designated as fingerprints here) from the protein distribution densities in the LD space defined by ln( L ) and ID are then constructed. The fingerprints for different organisms and protein sets are found to be distinct with each other, and they can therefore be used for comparative studies. As a test case, phylogenetic trees have been constructed based on the protein distribution densities in the fingerprints of proteomes of organisms without performing any protein sequence comparison and alignments. The phylogenetic trees generated are biologically meaningful, demonstrating that the protein distributions in the LD space may serve as unique phylogenetic signals of the organisms at the proteome level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

PASTIS v0.1

PASTIS is a distributed-memory library for aligning large protein sequence sets and constructing protein similarity networks.

Buluc, Aydin↗
DOE OSTI2023

Structural models and functional annotations for the Sphagnum divinum proteome

This dataset contains the structural models for the primary transcripts of the Sphagnum divinum proteome. Additionally, for a subset of these proteins, sequence and structural alignment results are provided. This dataset represents the most thorough structural study of a Sphagnum species, also known as peat mosses, by providing three-dimensional atomic resolution structures of the majority of the encoded proteins as well as structural alignment results used in the application of annotating the proteome. References (DOI) AlphaFold v2 Monomer: https://doi.org/10.1038/s41586-021-03819-2. References (DOI) US-align2: https://doi.org/10.1038/s41592-022-01585-1

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Metatransciptomic Analysis Data for Interactive effects of depth and differential irrigation on soil microbiome composition and functioning

RNA was collected from soil at different depths and after three different levels of irrigation T1 100% of normal field irrigation, T4: 18.75% or normal irrigation and T5: unirrigated controls. Total RNA was isolated using the Zymo Quick-RNA fecal/soil microbe miniprep (catalog no. R2040), incorporating the DNase I treatment using Zymo’s DNase I kit (catalog no. E1010). To increase the yield of RNA, we modified the manufacturer’s instructions by first doubling the amount of soil per extraction (from 0.25 g to 0.5 g) and by performing extractions in triplicate before pooling separate extractions together. Certain soil samples (largely those from deeper soil layers) had low yield (< 100 ng per extraction) so additional rounds of extraction were performed to obtain sufficient RNA. RNA concentration was assessed using a Qubit RNA HS assay kit (Thermo Fisher) and RNA quality was determined using an Agilent 2100 BioAnalyzer (Agilent; Santa Clara, CA). The resultant RNA samples were then sequenced by GENEWIZ using Illumina technology (GENEWIZ; South Plainfield, NJ). Sequences were then aligned to a soil metagenome previously obtained from the same site using the Burrows-Wheeler aligner (BWA). SAM files were then converted to raw counts using HTSeq.

Soil microbiome, metatranscriptomics↗
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 OSTI2015

Characterization of MtoD from Sideroxydans lithotrophicus : a cytochrome c electron shuttle used in lithoautotrophic growth

The autotrophic Sideroxydans lithotrophicus ES-1 can grow by coupling the oxidation of ferrous iron to the reduction of oxygen. Soluble ferrous iron is oxidized at the surface of the cell by an MtoAB porin-cytochrome complex that functions as an electron conduit through the outer membrane. Electrons are then transported to the cytoplasmic membrane where they are used to generate proton motive force (PMF) (for ATP synthesis) and NADH for autotrophic processes such as carbon fixation. As part of the mtoAB gene cluster, S. lithotrophicus also contains the gene mtoD that is proposed to encode a cytochrome c protein. We isolated mtoD from a Shewanella oneidensis expression system where the mtoD gene was expressed on a pBAD plasmid vector. Biochemical, biophysical, and crystallographic characterization of the purified MtoD revealed it as an 11 kDa monomeric protein containing a single heme. Sequence and structural alignment indicated that MtoD belonged to the class-1 cytochrome c family and had a similar fold to ferricytochrome c552 family, however the MtoD heme is bis-histidine coordinated and is substantially more exposed than the hemes of other family members. The reduction potential of the MtoD heme at pH 7 was +155 mV vs. Standard Hydrogen Electrode, which is approximately 100 mV lower than that of mitochondrial cytochrome c. Consideration of the properties of MtoD in the context of the potential respiratory partners identified from the genome suggests that MtoD could associate to multiple electron transfer partners as the primary periplasmic electron shuttle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2011

Fast Identification and Removal of Sequence Contamination from Genomic and Metagenomic Datasets

High-throughput sequencing technologies have strongly impacted microbiology, providing a rapid and cost-effective way of generating draft genomes and exploring microbial diversity. However, sequences obtained from impure nucleic acid preparations may contain DNA from sources other than the sample. Those sequence contaminations are a serious concern to the quality of the data used for downstream analysis, causing misassembly of sequence contigs and erroneous conclusions. Therefore, the removal of sequence contaminants is a necessary and required step for all sequencing projects. We developed DeconSeq, a robust framework for the rapid, automated identification and removal of sequence contamination in longer-read datasets (150 bp mean read length). DeconSeq is publicly available as standalone and web-based versions. The results can be exported for subsequent analysis, and the databases used for the web-based version are automatically updated on a regular basis. DeconSeq categorizes possible contamination sequences, eliminates redundant hits with higher similarity to non-contaminant genomes, and provides graphical visualizations of the alignment results and classifications. Using DeconSeq, we conducted an analysis of possible human DNA contamination in 202 previously published microbial and viral metagenomes and found possible contamination in 145 (72%) metagenomes with as high as 64% contaminating sequences. This new framework allows scientists to automatically detect and efficiently remove unwanted sequence contamination from their datasets while eliminating critical limitations of current methods. DeconSeq’s web interface is simple and user-friendly. The standalone version allows offline analysis and integration into existing data processing pipelines. DeconSeq’s results reveal whether the sequencing experiment has succeeded, whether the correct sample was sequenced, and whether the sample contains any sequence contamination from DNA preparation or host. In addition, the analysis of 202 metagenomes demonstrated significant contamination of the non-human associated metagenomes, suggesting that this method is appropriate for screening all metagenomes. DeconSeq is available at http://deconseq.sourceforge.net/.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Towards Scalable Optimal Sequence Homology Detection

Abstract—The field of bioinformatics and computational biol- ogy is experiencing a data revolution — experimental techniques to procure data have increased in throughput, improved in accuracy and reduced in costs. This has spurred an array of high profile sequencing and data generation projects. While the data repositories represent untapped reservoirs of rich information critical for scientific breakthroughs, the analytical software tools that are needed to analyze large volumes of such sequence data have significantly lagged behind in their capacity to scale. In this paper, we address homology detection, which is a funda- mental problem in large-scale sequence analysis with numerous applications. We present a scalable framework to conduct large- scale optimal homology detection on massively parallel super- computing platforms. Our approach employs distributed memory work stealing to effectively parallelize optimal pairwise alignment computation tasks. Results on 120,000 cores of the Hopper Cray XE6 supercomputer demonstrate strong scaling and up to 2.42 × 107 optimal pairwise sequence alignments computed per second (PSAPS), the highest reported in the literature.

Daily, Jeffrey A.↗
DOE OSTI2015

Structural Insights into Separase Architecture and Substrate Recognition through Computational Modelling of Caspase-Like and Death Domains

Separases are large proteins that mediate sister chromatid disjunction in all eukaryotes. They belong to clan CD of cysteine peptidases and contain a well-conserved C-terminal catalytic protease domain similar to caspases and gingipains. However, unlike other well-characterized groups of clan CD peptidases, there are no high-resolution structures of separases and the details of their regulation and substrate recognition are poorly understood. Here we undertook an in-depth bioinformatical analysis of separases from different species with respect to their similarity in amino acid sequence and protein fold in comparison to caspases, MALT-1 proteins (mucosa-associated lymphoidtissue lymphoma translocation protein 1) and gingipain-R. A comparative model of the single C-terminal caspase-like domain in separase from C. elegans suggests similar binding modes of substrate peptides between these protein subfamilies, and enables differences in substrate specificity of separase proteins to be rationalised. We also modelled a newly identified putative death domain, located N-terminal to the caspase-like domain. The surface features of this domain identify potential sites of protein-protein interactions. Notably, we identified a novel conserved region with the consensus sequence WWxxRxxLD predicted to be exposed on the surface of the death domain, which we termed the WR motif. We envisage that findings from our study will guide structural and functional studies of this important protein family.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Use of 16S rRNA gene for identification of a broad range of clinically relevant bacterial pathogens

According to World Health Organization statistics of 2011, infectious diseases remain in the top five causes of mortality worldwide. However, despite sophisticated research tools for microbial detection, rapid and accurate molecular diagnostics for identification of infection in humans have not been extensively adopted. Time-consuming culture-based methods remain to the forefront of clinical microbial detection. The 16S rRNA gene, a molecular marker for identification of bacterial species, is ubiquitous to members of this domain and, thanks to ever-expanding databases of sequence information, a useful tool for bacterial identification. In this study, we assembled an extensive repository of clinical isolates (n = 617), representing 30 medically important pathogenic species and originally identified using traditional culture-based or non-16S molecular methods. This strain repository was used to systematically evaluate the ability of 16S rRNA for species level identification. To enable the most accurate species level classification based on the paucity of sequence data accumulated in public databases, we built a Naïve Bayes classifier representing a diverse set of high-quality sequences from medically important bacterial organisms. We show that for species identification, a model-based approach is superior to an alignment based method. Overall, between 16S gene based and clinical identities, our study shows a genus-level concordance rate of 96% and a species-level concordance rate of 87.5%. We point to multiple cases of probable clinical misidentification with traditional culture based identification across a wide range of gram-negative rods and gram-positive cocci as well as common gram-negative cocci.

60 APPLIED LIFE SCIENCES↗