Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Base 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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

An adaptable pentaloop defines a robust neomycin-B RNA aptamer with conditional ligand-bound structures

Aptamers can be highly specific for their targets, which implies precise molecular recognition between aptamer and target. However, as small polymers, their structures are more subject to environmental conditions than the more constrained longer RNAs such as those that constitute the ribosome. To understand the balance between structural and environmental factors in establishing ligand specificity of aptamers, we examined the RNA aptamer (NEO1A) previously reported as specific for neomycin-B. We show that NEO1A can recognize other aminoglycosides with similar affinities as for neomycin-B and its aminoglycoside specificity is strongly influenced by ionic strength and buffer composition. NMR and 2-aminopurine (2AP) fluorescence studies of the aptamer identified a flexible pentaloop and a stable binding pocket. Consistent with a well structured binding pocket, docking analysis results correlated with experimental measures of the binding energy for most ligands. Steady state fluorescence studies of 2AP-substituted aptamers confirmed that A16 moves to a more solvent accessible position upon ligand binding while A14 moves to a less solvent accessible position, which is most likely a base stack. Analysis of binding affinities of NEO1A sequence variants showed that the base in position 16 interacts differently with each ligand and the interaction is a function of the buffer constituents. Our results show that the pentaloop provides NEO1A with the ability to adapt to external influences on its structure, with the critical base at position 16 adjusting to incorporate each ligand into a stable pocket by hydrophobic interactions and/or hydrogen bonds depending on the ligand and the ionic environment.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Combining de novo and reference-guided assembly with scaffold_builder

Genome sequencing has become routine, however genome assembly still remains a challenge despite the computational advances in the last decade. In particular, the abundance of repeat elements in genomes makes it difficult to assemble them into a single complete sequence. Identical repeats shorter than the average read length can generally be assembled without issue. However, longer repeats such as ribosomal RNA operons cannot be accurately assembled using existing tools. The application Scaffold_builder was designed to generate scaffolds – super contigs of sequences joined by N-bases – based on the similarity to a closely related reference sequence. This is independent of mate-pair information and can be used complementarily for genome assembly, e.g. when matepairs are not available or have already been exploited. Scaffold_builder was evaluated using simulated pyrosequencing reads of the bacterial genomes Escherichia coli 042, Lactobacillus salivarius UCC118 and Salmonella enterica subsp. enterica serovar Typhi str. P-stx-12. Moreover, we sequenced two genomes from Salmonella enterica serovar Typhimurium LT2 G455 and Salmonella enterica serovar Typhimurium SDT1291 and show that Scaffold_builder decreases the number of contig sequences by 53% while more than doubling their average length. Scaffold_builder is written in Python and is available at http://edwards.sdsu.edu/scaffold_builder. A web-based implementation is additionally provided to allow users to submit a reference genome and a set of contigs to be scaffolded.

Salmonella enterica serovar typhimurium↗
DOE OSTI2014

k-merSNP discovery: Software for alignment-and reference-free scalable SNP discovery, phylogenetics, and annotation for hundreds of microbial genomes

With the flood of whole genome finished and draft microbial sequences, we need faster, more scalable bioinformatics tools for sequence comparison. An algorithm is described to find single nucleotide polymorphisms (SNPs) in whole genome data. It scales to hundreds of bacterial or viral genomes, and can be used for finished and/or draft genomes available as unassembled contigs or raw, unassembled reads. The method is fast to compute, finding SNPs and building a SNP phylogeny in minutes to hours, depending on the size and diversity of the input sequences. The SNP-based trees that result are consistent with known taxonomy and trees determined in other studies. The approach we describe can handle many gigabases of sequence in a single run. The algorithm is based on k-mer analysis.

Source record↗
DOE OSTI2026

Data for FUN-PROSE: A Deep Learning Approach to Predict Condition-Specific Gene Expression in Fungi

mRNA levels of all genes in a genome is a critical piece of information defining the overall state of the cell in a given environmental condition. Being able to reconstruct such condition-specific expression in fungal genomes is particularly important to metabolically engineer these organisms to produce desired chemicals in industrially scalable conditions. Most previous deep learning approaches focused on predicting the average expression levels of a gene based on its promoter sequence, ignoring its variation across different conditions. Here we present FUN-PROSE—a deep learning model trained to predict differential expression of individual genes across various conditions using their promoter sequences and expression levels of all transcription factors. We train and test our model on three fungal species and get the correlation between predicted and observed condition-specific gene expression as high as 0.85. We then interpret our model to extract promoter sequence motifs responsible for variable expression of individual genes. We also carried out input feature importance analysis to connect individual transcription factors to their gene targets. A sizeable fraction of both sequence motifs and TF-gene interactions learned by our model agree with previously known biological information, while the rest corresponds to either novel biological facts or indirect correlations.

Genomics↗
NASA NTRSDate not supplied

Autonomous Control for Arbitrary Thruster Configurations and Mass Properties in Special Euclidean Group SE(3)

Most current methods for determining maneuvers and thrust firing sequences depend on explicit and predetermined commands generated by a combination of on-board systems and ground-based human-in-the-loop methods. For spacecraft and space structures with changing mass properties and thruster configurations, such as the Deep Space Gateway as it changes configurations throughout its lifetime, determining these commands can be time-consuming and computationally intensive. However, recent work within the Lie GroupSE (3) has offered ways of autonomously determining the location, power, precision, and capabilities of thrusters in any arbitrary position. Furthermore, a method for determining thruster firing sequences based on an arbitrary control input (both translational and rotational in a coupled, 6-element vector) and arbitrary thruster configurations has also recently been developed. When combining these methods, any spacecraft with any mass properties and thruster configurations can be understood in terms of controllability limits and thruster firing sequences can be generated quickly and with low computational load, thus extending the autonomous capabilities of deep space missions. In this work, this method is presented and explored in terms of computational load, robustness in the presence of uncertainty, and overall performance. The capabilities of this method are also examined in the case of the Deep Space Gateway both in fully controllable configurations and uncontrollable configurations.

Control↗
NASA NTRSDate not supplied

Autonomous Control for Arbitrary Thruster Configurations and Mass Properties in Special Euclidean Group SE(3)

Most current methods for determining maneuvers and thrust firing sequences depend on explicit and predetermined commands generated by a combination of on-board systems and ground-based human-in-the-loop methods. For spacecraft and space structures with changing mass properties and thruster configurations, such as the Deep Space Gateway as it changes configurations throughout its lifetime, determining these commands can be time-consuming and computationally intensive. However, recent work within the Lie group SE(3) has offered ways of autonomously determining the location, power, precision, and capabilities of thrusters in any arbitrary position. Furthermore, a method for determining thruster firing sequences based on an arbitrary control input (both translational and rotational in a coupled, 6-element vector) and arbitrary thruster configurations has also recently been developed. When combining these methods, any spacecraft with any mass properties and thruster configurations can be understood in terms of controllability limits and thruster firing sequences can be generated quickly and with low computational load, thus extending the autonomous capabilities of deep space missions. In this work, this method is presented and explored in terms of convergence time to the desired pose. The capabilities of this method are also examined in the case of the Deep Space Gateway both in fully controllable configurations and uncontrollable configurations.

SE(3)↗
DOE OSTI2014

Inference of Expanded Lrp-Like Feast/Famine Transcription Factor Targets in a Non-Model Organism Using Protein Structure-Based Prediction

Widespread microbial genome sequencing presents an opportunity to understand the gene regulatory networks of non-model organisms. This requires knowledge of the binding sites for transcription factors whose DNA-binding properties are unknown or difficult to infer. We adapted a protein structure-based method to predict the specificities and putative regulons of homologous transcription factors across diverse species. As a proof-of-concept we predicted the specificities and transcriptional target genes of divergent archaeal feast/famine regulatory proteins, several of which are encoded in the genome of Halobacterium salinarum. This was validated by comparison to experimentally determined specificities for transcription factors in distantly related extremophiles, chromatin immunoprecipitation experiments, and cis-regulatory sequence conservation across eighteen related species of halobacteria. Through this analysis we were able to infer that Halobacterium salinarum employs a divergent local trans-regulatory strategy to regulate genes (carA and carB) involved in arginine and pyrimidine metabolism, whereas Escherichia coli employs an operon. The prediction of gene regulatory binding sites using structure-based methods is useful for the inference of gene regulatory relationships in new species that are otherwise difficult to infer.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Draft Genome Sequence of Comamonas thiooxydans Strain S23 T (DSM 17888 T ), a Thiosulfate-Oxidizing Bacterium Isolated from a Sulfur Spring in India

The genus Comamonas contains species isolated from various environments, such as termite guts, wetlands, activated sludge, soil, humans, and fresh water. Here, we report the draft genome sequence of Comamonas thiooxydans strain S23 T capable of oxidizing thiosulfate under mixotrophic growth conditions. Based upon draft genome sequencing, the genome is 5.3 Mb and encodes 4,767 proteins. The Comamonas thiooxydans whole-genome sequence will help understand the metabolic diversity in sulfur oxidation pathways.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Deletions at short direct repeats and base substitutions are characteristic mutations for bleomycin-induced double- and single-strand breaks, respectively, in a human shuttle vector system

Using the radiomimetic drug, bleomycin, we have determined the mutagenic potential of DNA strand breaks in the shuttle vector pZ189 in human fibroblasts. The bleomycin treatment conditions used produce strand breaks with 3'-phosphoglycolate termini as > 95% of the detectable dose-dependent lesions. Breaks with this end group represent 50% of the strand break damage produced by ionizing radiation. We report that such strand breaks are mutagenic lesions. The type of mutation produced is largely determined by the type of strand break on the plasmid (i.e. single versus double). Mutagenesis studies with purified DNA forms showed that nicked plasmids (i.e. those containing single-strand breaks) predominantly produce base substitutions, the majority of which are multiples, which presumably originate from error-prone polymerase activity at strand break sites. In contrast, repair of linear plasmids (i.e. those containing double-strand breaks) mainly results in deletions at short direct repeat sequences, indicating the involvement of illegitimate recombination. The data characterize the nature of mutations produced by single- and double-strand breaks in human cells, and suggests that deletions at direct repeats may be a 'signature' mutation for the processing of DNA double-strand breaks.

NASA Discipline Radiation Health↗
DOE OSTI2022

Multi-head attention-based U-Nets for predicting protein domain boundaries using 1D sequence features and 2D distance maps

Abstract The information about the domain architecture of proteins is useful for studying protein structure and function. However, accurate prediction of protein domain boundaries (i.e., sequence regions separating two domains) from sequence remains a significant challenge. In this work, we develop a deep learning method based on multi-head U-Nets (called DistDom) to predict protein domain boundaries utilizing 1D sequence features and predicted 2D inter-residue distance map as input. The 1D features contain the evolutionary and physicochemical information of protein sequences, whereas the 2D distance map includes the structural information of proteins that was rarely used in domain boundary prediction before. The 1D and 2D features are processed by the 1D and 2D U-Nets respectively to generate hidden features. The hidden features are then used by the multi-head attention to predict the probability of each residue of a protein being in a domain boundary, leveraging both local and global information in the features. The residue-level domain boundary predictions can be used to classify proteins as single-domain or multi-domain proteins. It classifies the CASP14 single-domain and multi-domain targets at the accuracy of 75.9%, 13.28% more accurate than the state-of-the-art method. Tested on the CASP14 multi-domain protein targets with expert annotated domain boundaries, the average per-target F1 measure score of the domain boundary prediction by DistDom is 0.263, 29.56% higher than the state-of-the-art method.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Sequence Impedance Modeling of Grid-Forming Converters: Preprint

Grid-forming control of inverter-based resources (IBRs) has been identified as a critical technology for operating power systems with high levels of inverter-based resources. This paper presents sequence impedance modeling of grid-forming inverters to evaluate its small-signal stability properties. Droop based grid-forming control with full consideration of the multiloop structure is studied, and the impact of individual controller on the inverter impedance characteristics is discussed. The developed sequence impedance model is compared with that of the grid-following inverter. It is found from the developed sequence impedance models that grid-forming inverters are less prone to harmonic resonance problems during operation with weaker grids. The developed models are validated using PSCAD simulations.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

GMFOLD: Subgraph matching for high-throughput DNA-aptamer secondary structure classification and machine learning interpretability

Aptamers are oligonucleotide receptors that bind to their targets with high affinity. Here, we consider aptamers comprised of single-stranded DNA that undergo target-binding-induced conformational changes, giving rise to unique secondary and tertiary structures. Given a specific aptamer primary sequence, there are well-established computational tools (notably mfold) to predict the secondary structure via free energy minimization algorithms. While mfold generates secondary structures for individual sequences, there is a need for a high-throughput process whereby thousands of DNA structures can be predicted in real-time for use in an interactive setting, when combined with aptamer selections that generate candidate pools that are too large to be experimentally interrogated. We developed a new Python code for high-throughput aptamer secondary structure determination (GMfold). GMfold uses subgraph matching methods to group aptamer candidates by secondary structure similarities. We also improve an open-source code, SeqFold, to incorporate subgraph matching concepts. We represent each secondary structure as a lowest-energy bipartite subgraph matching of the DNA graph to itself. These new tools enable thousands of DNA sequences to be compared based on their secondary structures, using machine-learning algorithms. This process is advantageous when analyzing sequences that arise from aptamer selections via systematic evolution of ligands by exponential enrichment (SELEX). This work is a building block for future machine-learning-informed DNA-aptamer selection processes to identify aptamers with improved target affinity and selectivity and advance aptamer biosensors and therapeutics.

Aptamer↗
DOE OSTI2019

High spatiotemporal variability of bacterial diversity over short time scales with unique hydrochemical associations within a shallow aquifer

Understanding microbial community structure and function within the subsurface is critical to assessing overall quality and maintenance of groundwater; however, the factors that determine microbial community assembly, structure, and function in groundwater systems and their impact on water quality remains poorly understood. In this study, three shallow wells (FW301, FW303, FW305) in a non-contaminated shallow aquifer in the ENIGMA-Oak Ridge Field Research Center (Oak Ridge, TN) were sampled approximately 3 times a week over a period of three months to measure changes in groundwater geochemistry and microbial diversity. It was expected that the sampled microbial diversity from two historic field wells (FW301, FW303) would be relatively stable, while diversity from a newer well (FW305) would be less stable over time. The wells displayed some degree of hydrochemical variability over time unique to each well, with FW303 being overall the most stable well and FW301 being the most dynamic based upon dissolved oxygen, conductivity, and nitrate. Community analysis via ss-rRNA paired-end sequencing and distribution-based clustering revealed higher OTU richness, diversity, and variability in groundwater communities of FW301 than the other two wells for diversity binned over all time points. Microbial community composition of a given well was on average > 50% dissimilar to any other well at a given time (days), yet, functional gene diversity as measured with GeoChip remained relatively constant. Similarities in community structure across wells were observed with respect to the presence of 20 shared bacterial groups in all samples in all wells, although at varying levels over the tested time period. Similarity percentage (SIMPER) analysis revealed that variability in FW301 was largely attributed to low abundance, highly-transient populations, while variability in the most hydrochemically stable well (FW303) was due to fluctuations in more highly abundant and frequently present taxa. Additionally, the youngest well FW305 showed a dramatic shift in community composition towards the end of the sampling period that was not observed in the other wells, suggesting possible succession events over time. Time-series analysis using vector auto-regressive models and Granger causality showed unique relationships between richness and geochemistry over time in each well. Finally, these results indicate temporally dynamic microbial communities over short time scales, with day-to-day population shifts in local community structure influenced by available source community diversity and local groundwater hydrochemistry.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

EvoProtGrad (Directed Evolution for Proteins with Gradients) [SWR-23-48]

A Python package for directed evolution on a protein sequence with gradient-based discrete Markov chain monte carlo (MCMC). Users are able to compose custom models that map sequence to function with pretrained models, including protein language models (PLMs), to guide and constrain search. Our package natively integrates with the HuggingFace platform and supports PLMs from transformers. Our MCMC sampler identifies promising amino acids to mutate via model gradients taken with respect to the input (i.e., sensitivity analysis). We allow users to compose their own custom target function for MCMC by leveraging the Product of Experts MCMC paradigm. Each model is an "expert" that contributes its own knowledge about the protein's fitness landscape to the overall target function. The sampler is designed to be more efficient and effective than brute force and random search while maintaining most of the generality and flexibility. Additional information can be found in the related publication: https://iopscience.iop.org/article/10.1088/2632-2153/accacd

Emami, Patrick↗
DOE OSTI2020

Geochemistry, RNA radiocarbon signature, and metatranscriptome data from Rifle, Colorado groundwater well

Understanding the different sources of carbon subsurface microbes use is highly informative yet remains challenging. The goal of this work was to develop and test a cultivation-independent method that analyzes both the radiocarbon and ribosomal RNA (rRNA) signature of unamplified RNA from the same environmental sample, and that collectively characterizes community carbon utilization under in situ conditions. The method was applied to a filtered groundwater sample from a geochemically reduced zone of a shallow alluvial aquifer near Rife, Colorado (CO), where multiple potential carbon sources co-occur. The radiocarbon and RNA-Seq results, while addressing diverse sample properties, indicated a strikingly similar conclusion about the quantitative role of autotrophy in the subsurface microbial community.This data package provides (1) ribosomal RNA (rRNA) sequences for groundwater microbial populations, (2) the taxonomic identity and relative abundance of groundwater microbial populations based on these rRNA sequences, (3) radiocarbon signatures of the rRNA indicating the source of carbon used by the groundwater microbial populations (provided in both Word docx file and pdf file), (4) geochemical measurements for the monitoring well corresponding to these measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

The Dark Side of the Mushroom Spring Microbial Mat: Life in the Shadow of Chlorophototrophs. I. Microbial Diversity Based on 16S rRNA Gene Amplicons and Metagenomic Sequencing

Microbial-mat communities in the effluent channels of Octopus and Mushroom Springs within the Lower Geyser Basin at Yellowstone National Park have been studied for nearly 50 years. The emphasis has mostly focused on the chlorophototrophic bacterial organisms of the phyla Cyanobacteria and Chloroflexi . In contrast, the diversity and metabolic functions of the heterotrophic community in the microoxic/anoxic region of the mat are not well understood. In this study we analyzed the orange-colored undermat of the microbial community of Mushroom Spring using metagenomic and rRNA-amplicon (iTag) analyses. Our analyses disclosed a highly diverse community exhibiting a high degree of unevenness, strongly dominated by a single taxon, the filamentous anoxygenic phototroph, Roseiflexus spp. The second most abundant organisms belonged to the Thermotogae , which have been hypothesized to be a major source of H-2 from fermentation that could enable photomixotrophic metabolism by Chloroflexus and Roseiflexus spp. Other abundant organisms include two members of the Armatimonadetes (OP10); Thermocrinis sp.; and phototrophic and heterotrophic members of the Chloroflexi . Further, an Atribacteria (OP9/JS1) member; a sulfate-reducing Therrnodesulfovibrio sp.; a Planctomycetes member; a member of the EM3 group tentatively affiliated with the Thermotogae , as well as a putative member of the Arrninicenantes (OP8) represented ≥ 1% of the reads. Archaea were not abundant in the iTag analysis, and no metagenomic bin representing an archaeon was identified. A high microdiversity of 16S rRNA gene sequences was identified for the dominant taxon, Roseiflexus spp. Previous studies demonstrated that highly similar Synechococcus variants in the upper layer of the mats represent ecological species populations with specific ecological adaptations. In conclusion, this study suggests that similar putative ecotypes specifically adapted to different niches occur within the undermat community, particularly for Roseiflexus spp.

59 BASIC BIOLOGICAL SCIENCES↗