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 325 records · Page 18

Rotational structures and the wobbling mode in {sup 167}Ta

Excited states in the neutron-deficient nucleus {sup 167}Ta were studied through the {sup 120}Sn({sup 51}V,4n) reaction. Twelve rotational bands have been observed and the relative excitation energy of each sequence is now known owing to the multiple interband connections. Several quasineutron alignments were observed that aided in the quasiparticle assignments of these bands. The resulting interpretation is in line with observations in neighboring nuclei. Trends in the wobbling phonon energy seen in {sup 161,163,165,167}Lu and {sup 167}Ta are also discussed and particle-rotor model calculations (assuming constant moments of inertia) are found to be inconsistent with the experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

VITAL NMR: Using Chemical Shift Derived Secondary Structure Information for a Limited Set of Amino Acids to Assess Homology Model Accuracy

Homology modeling is a powerful tool for predicting protein structures, whose success depends on obtaining a reasonable alignment between a given structural template and the protein sequence being analyzed. In order to leverage greater predictive power for proteins with few structural templates, we have developed a method to rank homology models based upon their compliance to secondary structure derived from experimental solid-state NMR (SSNMR) data. Such data is obtainable in a rapid manner by simple SSNMR experiments (e.g., (13)C-(13)C 2D correlation spectra). To test our homology model scoring procedure for various amino acid labeling schemes, we generated a library of 7,474 homology models for 22 protein targets culled from the TALOS+/SPARTA+ training set of protein structures. Using subsets of amino acids that are plausibly assigned by SSNMR, we discovered that pairs of the residues Val, Ile, Thr, Ala and Leu (VITAL) emulate an ideal dataset where all residues are site specifically assigned. Scoring the models with a predicted VITAL site-specific dataset and calculating secondary structure with the Chemical Shift Index resulted in a Pearson correlation coefficient (-0.75) commensurate to the control (-0.77), where secondary structure was scored site specifically for all amino acids (ALL 20) using STRIDE. This method promises to accelerate structure procurement by SSNMR for proteins with unknown folds through guiding the selection of remotely homologous protein templates and assessing model quality.

59 BASIC BIOLOGICAL SCIENCES↗

Physically motivated global alignment method for electron tomography

Electron tomography is widely used for nanoscale determination of 3-D structures in many areas of science. Determining the 3-D structure of a sample from electron tomography involves three major steps: acquisition of sequence of 2-D projection images of the sample with the electron microscope, alignment of the images to a common coordinate system, and 3-D reconstruction and segmentation of the sample from the aligned image data. The resolution of the 3-D reconstruction is directly influenced by the accuracy of the alignment, and therefore, it is crucial to have a robust and dependable alignment method. In this paper, we develop a new alignment method which avoids the use of markers and instead traces the computed paths of many identifiable ‘local’ center-of-mass points as the sample is rotated. Compared with traditional correlation schemes, the alignment method presented here is resistant to cumulative error observed from correlation techniques, has very rigorous mathematical justification, and is very robust since many points and paths are used, all of which inevitably improves the quality of the reconstruction and confidence in the scientific results.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

De novo transcriptome assembly databases for the central nervous system of the medicinal leech

The study of non-model organisms stands to benefit greatly from genetic and genomic data. For a better understanding of the molecular mechanisms driving neuronal development, and to characterize the entire leech Hirudo medicinalis central nervous system (CNS) transcriptome we combined Trinity for de-novo assembly and Illumina HiSeq2000 for RNA-Seq. We present a set of 73,493 de-novo assembled transcripts for the leech, reconstructed from RNA collected, at a single ganglion resolution, from the CNS. This set of transcripts greatly enriches the available data for the leech. Here, we share two databases, such that each dataset allows a different type of search for candidate homologues. The first is the raw set of assembled transcripts. This set allows a sequence-based search. A comprehensive analysis of which revealed 22,604 contigs with high e-values, aligned versus the Swiss-Prot database. This analysis enabled the production of the second database, which includes correlated sequences to annotated transcript names, with the confidence of BLAST best hit.

59 BASIC BIOLOGICAL SCIENCES↗

Alignment of the ATLAS Inner Detector in Run 2

The performance of the ATLAS Inner Detector alignment has been studied using pp collision data at s√=13TeV collected by the ATLAS experiment during Run 2 (2015–2018) of the Large Hadron Collider (LHC). The goal of the detector alignment is to determine the detector geometry as accurately as possible and correct for time-dependent movements. The Inner Detector alignment is based on the minimization of track-hit residuals in a sequence of hierarchical levels, from global mechanical assembly structures to local sensors. Subsequent levels have increasing numbers of degrees of freedom; in total there are almost 750,000. The alignment determines detector geometry on both short and long timescales, where short timescales describe movements within an LHC fill. The performance and possible track parameter biases originating from systematic detector deformations are evaluated. Momentum biases are studied using resonances decaying to muons or to electrons. The residual sagitta bias and momentum scale bias after alignment are reduced to less than ~0.1 TeV -1 and 0.9×10 -3 , respectively. Impact parameter biases are also evaluated using tracks within jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparing Top-Down Proteoform Identification: Deconvolution, PrSM Overlap, and PTM Detection

Generating top-down tandem mass spectra (MS/MS) for complex mixtures of proteoforms has become possible through improvements in fractionation, on-line separation, dissociation, and mass analysis. The algorithms to match tandem mass spectra to sequences have undergone a parallel evolution, with both spectral alignment and peak matching being paired with diverse methods for scoring proteoform-spectral matches (PrSMs). This study assesses state-of-the-art algorithms for top-down identification through three distinct challenges. The first is identifying a large yield of PrSMs while controlling false discovery rate (FDR) in identifying thousands of proteoforms from complex cell lysates via four software workflows: ProSight Proteome Discoverer, TopPIC, Informed Proteomics, and pTop. The second is the deconvolution of data from both Thermo Orbitrap-class and Bruker maXis Q-TOF instruments to produce consistent precursor charge and mass determinations while generating fragment mass lists to optimize identification. The third attempts to detect diverse post-translational modifications (PTMs) in proteoforms from cow milk and human ovarian tissue. The data demonstrate that existing software suites produce admirable sensitivity, in some cases identifying a third of collected tandem mass spectra with FDR controlled below 2%; the overlap in these PrSMs, however, illustrates real value in searching data with multiple search engines. Differences among identification workflows seem to result from each search algorithm incorporating its own deconvolution algorithm. By transmitting deconvolution data from multiple deconvolution routes (Thermo Xtract, Bruker Auto MSn, Mascot Distiller, TopFD, and FLASHDeconv) to the downstream TopPIC search algorithm, we were able to detect common causes of deconvolution disagreement. The detection of PTMs was very inconsistent among search algorithms, with some workflows suggesting as little as 1% of PrSMs from cow’s milk were singly-phosphorylated while other workflows found that 18% of PrSMs were singly-phosphorylated. Taken together, these results make a strong argument for top-down researchers to adopt a standard practice of analyzing each MS/MS experiment with at least two different search engines.

59 BASIC BIOLOGICAL SCIENCES↗

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based translated sequence annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long-read sequencing data and in the context of pseudogenes.

59 BASIC BIOLOGICAL SCIENCES↗

Crystallization and preliminary crystallographic analysis of the type IIL restriction enzyme MmeI in complex with DNA

Type IIL restriction enzymes have rejuvenated the search for user-specified DNA binding and cutting. By aligning and contrasting the highly comparable amino-acid sequences yet diverse recognition specificities across the family of enzymes, amino acids involved in DNA binding have been identified and mutated to produce alternative binding specificities. To date, the specificity of MmeI (a type IIL restriction enzyme) has successfully been altered at positions 3, 4 and 6 of the asymmetric TCCRAC (where R is a purine) DNA-recognition sequence. To further understand the structural basis of MmeI DNA-binding specificity, the enzyme has been crystallized in complex with its DNA substrate. The crystal belonged to space group P1, with unit-cell parameters a = 61.73, b = 94.96, c = 161.24 {angstrom}, {alpha} = 72.79, {beta} = 89.12, {gamma} = 71.68{sup o}, and diffracted to 2.6 {angstrom} resolution when exposed to synchrotron radiation. The structure promises to reveal the basis of MmeI DNA-binding specificity and will complement efforts to create enzymes with novel specificities.

59 BASIC BIOLOGICAL SCIENCES↗

A Microfluidic DNA Library Preparation Platform for Next-Generation Sequencing

Next-generation sequencing (NGS) is emerging as a powerful tool for elucidating genetic information for a wide range of applications. Unfortunately, the surging popularity of NGS has not yet been accompanied by an improvement in automated techniques for preparing formatted sequencing libraries. To address this challenge, we have developed a prototype microfluidic system for preparing sequencer-ready DNA libraries for analysis by Illumina sequencing. Our system combines droplet-based digital microfluidic (DMF) sample handling with peripheral modules to create a fully-integrated, sample-in library-out platform. In this report, we use our automated system to prepare NGS libraries from samples of human and bacterial genomic DNA. E. coli libraries prepared on-device from 5 ng of total DNA yielded excellent sequence coverage over the entire bacterial genome, with .99% alignment to the reference genome, even genome coverage, and good quality scores. Furthermore, we produced a de novo assembly on a previously unsequenced multi-drug resistant Klebsiella pneumoniae strain BAA-2146 (KpnNDM). The new method described here is fast, robust, scalable, and automated. Our device for library preparation will assist in the integration of NGS technology into a wide variety of laboratories, including small research laboratories and clinical laboratories.

59 BASIC BIOLOGICAL SCIENCES↗

Design and calculations for HWR Strongback

The PIP-II Superconductive Linac is a linear accelerator that provides the first stage of beam acceleration. The velocity increase is achieved through a sequence of different cryomodules, where all the cavities are aligned to form a continuous string. The cryomodule team is currently developing a spare for the HWR cryomodule. This spare will not be an exact replica of the existing one, but rather will be based on the SSR f650 cryomodule design created for PIP-II . Each cavity , magnet , and coupler form a unit, with a total of eight units per string. These components, together with additional equipment (such as piping , the thermal shield , and beamline interconnections ), are supported by the Strongback . The aim of this project is to develop a brand-new Strongback design for HWR cavities, and to perform the related analyses and calculations in order to meet the required specifications.

43 PARTICLE ACCELERATORS↗

GenomeDepot: data management system for microbial comparative genomics

Summary GenomeDepot is an open-source 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 websites for user-defined genome collections that include bioinformatic tools for interactive genome browsing, Basic Local Alignment Search Tool (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. Availability and implementation GenomeDepot is open source and distributed under the GNU General Public License via GitHub (https://github.com/aekazakov/genome-depot). GenomeDepot is implemented in Python and was tested in Ubuntu Linux. Full installation instructions and documentation are available at https://aekazakov.github.io/genome-depot/. GenomeDepot demo server is freely accessible at https://iseq.lbl.gov/demogd/.

Kazakov, Alexey [Lawrence Berkeley National Labora↗

Spatially Aligned Binary Single-Site Catalyst on Defective SiO 2 for Cascading Reactions

Capitalizing on the success of single-atom catalysts (SACs), dual-atom catalysts (DACs) have emerged as a new frontier in heterogeneous catalysis. However, most SACs and DACs studies seek to uniformly distribute the catalytic sites on the support material, which can hinder their effectiveness in intricate multistep cascading reactions. Particularly, it is a grand challenge to precisely control the spatial distribution of two different single sites forming binary sites so that reactants and intermediates contact the catalytic sites in the exact sequence required by the reaction steps. Here, in this work, we report a new type of binary single-site catalyst, Cu 1 –Zr 1 @SiO 2 , with Cu 1 and Zr 1 sites spatially aligned with the reaction sequence of the cascade reactions. The catalyst is synthesized by a modified reverse microemulsion approach, with single Cu sites anchored by nonbridging oxygen hole centers, which were induced by doping single Zr sites into SiO 2 . Low-energy ion scattering spectroscopy (LEIS) reveals that the outermost surface of the catalyst contains only Cu single sites, while the Zr sites are dispersed in the bulk. The catalytic performance is demonstrated in ethanol conversion to butenes, a model cascade reaction which includes ethanol dehydrogenation and aldol condensation steps. The precisely spatially controlled binary sites enable ethanol to first undergo dehydrogenation to acetaldehyde on Cu sites, followed by aldol condensation of acetaldehyde on Zr sites. As a result, C 3+ olefins selectivity as high as 77.0% (56.0% selectivity of butenes) is achieved by suppressing ethylene formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of the Impact of Concentration and Extraction Methods on the Targeted Sequencing of Human Viruses from Wastewater

Sequencing human viruses in wastewater is challenging due to their low abundance compared to the total microbial background. This study compared the impact of four virus concentration/extraction methods (Innovaprep, Nanotrap, Promega, and Solids extraction) on probe-capture enrichment for human viruses followed by sequencing. Different concentration/extraction methods yielded distinct virus profiles. Innovaprep ultrafiltration (following solids removal) had the highest sequencing sensitivity and richness, resulting in the successful assembly of several near-complete human virus genomes. However, it was less sensitive in detecting SARS-CoV-2 by digital polymerase chain reaction (dPCR) compared to Promega and Nanotrap. Across all preparation methods, astroviruses and polyomaviruses were the most highly abundant human viruses, and SARS-CoV-2 was rare. These findings suggest that sequencing success can be increased using methods that reduce nontarget nucleic acids in the extract, though the absolute concentration of total extracted nucleic acid, as indicated by Qubit, and targeted viruses, as indicated by dPCR, may not be directly related to targeted sequencing performance. Further, using broadly targeted sequencing panels may capture viral diversity but risks losing signals for specific low-abundance viruses. Overall, this study highlights the importance of aligning wet lab and bioinformatic methods with specific goals when employing probe-capture enrichment for human virus sequencing from wastewater.

59 BASIC BIOLOGICAL SCIENCES↗

PNNL-CompBio/perseq

Per sequence functional classification and Functional annotation in addition to linking those functions to taxonomic groups presents a challenge due to many established tools performing one task and not the other creating disparate downstream data. PerSeq aims to address these challenges by performing local alignments against a reference database annotated for both aspects, function and taxonomy, to annotate sequences individually. Downstream results are merged into count tables facilitating complex analyses based on observed function or functional potential in addition to being able to separate out function by taxonomic level to determine metabolic contributions by present organisms.

Brown, Joe↗

Profiling expression strategies for a type III polyketide synthase in a lysate-based, cell-free system

Abstract Some of the most metabolically diverse species of bacteria (e.g., Actinobacteria) have higher GC content in their DNA, differ substantially in codon usage, and have distinct protein folding environments compared to tractable expression hosts like Escherichia coli . Consequentially, expressing biosynthetic gene clusters (BGCs) from these bacteria in E. coli often results in a myriad of unpredictable issues with regard to protein expression and folding, delaying the biochemical characterization of new natural products. Current strategies to achieve soluble, active expression of these enzymes in tractable hosts can be a lengthy trial-and-error process. Cell-free expression (CFE) has emerged as a valuable expression platform as a testbed for rapid prototyping expression parameters. Here, we use a type III polyketide synthase from Streptomyces griseus , RppA, which catalyzes the formation of the red pigment flaviolin, as a reporter to investigate BGC refactoring techniques. We applied a library of constructs with different combinations of promoters and rppA coding sequences to investigate the synergies between promoter and codon usage. Subsequently, we assess the utility of cell-free systems for prototyping these refactoring tactics prior to their implementation in cells. Overall, codon harmonization improves natural product synthesis more than traditional codon optimization across cell-free and cellular environments. More importantly, the choice of coding sequences and promoters impact protein expression synergistically, which should be considered for future efforts to use CFE for high-yield protein expression. The promoter strategy when applied to RppA was not completely correlated with that observed with GFP, indicating that different promoter strategies should be applied for different proteins. In vivo experiments suggest that there is correlation, but not complete alignment between expressing in cell free and in vivo. Refactoring promoters and/or coding sequences via CFE can be a valuable strategy to rapidly screen for catalytically functional production of enzymes from BCGs, which advances CFE as a tool for natural product research.

59 BASIC BIOLOGICAL SCIENCES↗

Intrinsic alignments in redMaPPer clusters – I. Central galaxy alignments and angular segregation of satellites

The shapes of cluster central galaxies are not randomly oriented, but rather exhibit coherent alignments with the shapes of their parent clusters as well as with the surrounding large-scale structures. In this work, we aim to identify the galaxy and cluster quantities that most strongly predict the central galaxy alignment phenomenon among a large parameter space with a sample of 8237 clusters and 94 817 members within 0.1 < z < 0.35, based on the red-sequence Matched-filter Probabilistic Percolation cluster catalogue constructed from the Sloan Digital Sky Survey. We first quantify the alignment between the projected central galaxy shapes and the distribution of member satellites, to understand what central galaxy and cluster properties most strongly correlate with these alignments. Next, we investigate the angular segregation of satellites with respect to their central galaxy major axis directions, to identify the satellite properties that most strongly predict their angular segregation. We find that central galaxies are more aligned with their member galaxy distributions in clusters that are more elongated and have higher richness, and for central galaxies with larger physical size, higher luminosity and centring probability, and redder colour. Satellites with redder colour, higher luminosity, located closer to the central galaxy, and with smaller ellipticity show a stronger angular segregation towards their central galaxy major axes. Lastly, we provide physical explanations for some of the identified correlations, and discuss the connection to theories of central galaxy alignments, the impact of primordial alignments with tidal fields, and the importance of anisotropic accretion.

79 ASTRONOMY AND ASTROPHYSICS↗

Identification and removal of ribosomal RNA sequences from metatranscriptomes

Here, we present riboPicker, a robust framework for the rapid, automated identification and removal of ribosomal RNA sequences from metatranscriptomic datasets. The results can be exported for subsequent analysis, and the databases used for the web-based version are updated on a regular basis. riboPicker categorizes rRNA-like sequences and provides graphical visualizations and tabular outputs of ribosomal coverage, alignment results and taxonomic classifications.

59 BASIC BIOLOGICAL SCIENCES↗

Decoupled and semi-decoupled bands in Hg 197 and Hg 199

Here, the structure of decoupled positive-parity and semidecoupled negative-parity bands up to high spin in 197,199 Hg has been studied using multinucleon transfer reactions with the Gammasphere array. A positive-parity band in 199 Hg with the (vi 13/2 ) 3 decoupled configuration that deexcites to levels in the negative-parity sequences is newly established. The semidecoupled sequences in both isotopes are extended to higher spin. Rotation-aligned angular momenta and crossing frequencies determined from experiment are in fair agreement with those from cranking calculations. Finally, a transition from oblate collective rotation at low spin to triaxial noncollective shapes at high spin is indicated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗