Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence features”

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 37 records · Page 2

Ultracompact hybrid stars consistent with multimessenger astrophysics

In this work, we consider the consequences of phase transition in dense QCD on the properties of compact stars and implications for the observational program in gravitational wave and x-ray astrophysics. The key underlying assumption of our modeling is a strong first-order phase transition past the point where the hadronic branch of compact stars reaches the two-solar mass limit. Furthermore, our analysis predicts ultracompact stars with very small radii—in the range of 6–9 km—living on compact star sequences that are entirely consistent with the current multimessenger data. We show that sequences featuring two-solar mass hadronic stars consistent with radio-pulsar observations are also consistent with the inferences of large radii for massive neutron stars by NICER x-ray observations of neutron stars and the small radii predicted by gravitational waves analysis of the binary neutron star inspiral event GW170817 for our models that feature a strong first-order QCD phase transition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Amplification and adaptation of centromeric repeats in polyploid switchgrass species

Centromeres in most higher eukaryotes are composed of long arrays of satellite repeats from a single satellite repeat family. Why centromeres are dominated by a single satellite repeat and how the satellite repeats originate and evolve are among the most intriguing and long-standing questions in centromere biology. We identified eight satellite repeats in the centromeres of tetraploid switchgrass (Panicum virgatum). Seven repeats showed characteristics associated with classical centromeric repeats with monomeric lengths ranging from 166 to 187 bp. Interestingly, these repeats share an 80-bp DNA motif. We demonstrate that this 80-bp motif may dictate translational and rotational phasing of the centromeric repeats with the cenH3 nucleosomes. The sequence of the last centromeric repeat, Pv156, is identical to the 5S ribosomal RNA genes. We demonstrate that a 5S ribosomal RNA gene array was recruited to be the functional centromere for one of the switchgrass chromosomes. Our findings reveal that certain types of satellite repeats, which are associated with unique sequence features and are composed of monomers in mono-nucleosomal length, are favorable for centromeres. Centromeric repeats may undergo dynamic amplification and adaptation before the centromeres in the same species become dominated by the best adapted satellite repeat.

59 BASIC BIOLOGICAL SCIENCES↗

Design of intrinsically disordered protein variants with diverse structural properties

Intrinsically disordered proteins (IDPs) perform a broad range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Computational design of IDPs with specific conformational properties has, however, been difficult because of their substantial dynamics and structural complexity. We describe a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs that differ in compaction, long-range contacts, and propensity to phase separate. We experimentally tested and validated our designs and analyzed the sequence features that determine conformations. We show how our results are captured by a machine learning model, enabling us to speed up the algorithm. Our work expands the toolbox for computational protein design and will facilitate the design of proteins whose functions exploit the many properties afforded by protein disorder.

Science & Technology - Other Topics↗

Soil Metabolomics Predict Microbial Taxa as Biomarkers of Moisture Status in Soils from a Tidal Wetland

We present observations from a laboratory-controlled study on the impacts of extreme wetting and drying on a wetland soil microbiome. Our approach was to experimentally challenge the soil microbiome to understand impacts on anaerobic carbon cycling processes as the system transitions from dryness to saturation and vice-versa. Specifically, we tested for impacts on stress responses related to shifts from wet to drought conditions. We used a combination of high-resolution data for small organic chemical compounds (metabolites) and biological (community structure based on 16S rRNA gene sequencing) features. Using a robust correlation-independent data approach, we further tested the predictive power of soil metabolites for the presence or absence of taxa. Here, we demonstrate that taking an untargeted, multidimensional data approach to the interpretation of metabolomics has the potential to indicate the causative pathways selecting for the observed bacterial community structure in soils.

16S rRNA gene↗

An electrostatic selection mechanism controls sequential kinase signaling downstream of the T cell receptor

The sequence of events that initiates T cell signaling is dictated by the specificities and order of activation of the tyrosine kinases that signal downstream of the T cell receptor. Using a platform that combines exhaustive point-mutagenesis of peptide substrates, bacterial surface-display, cell sorting, and deep sequencing, we have defined the specificities of the first two kinases in this pathway, Lck and ZAP-70, for the T cell receptor ζ chain and the scaffold proteins LAT and SLP-76. We find that ZAP-70 selects its substrates by utilizing an electrostatic mechanism that excludes substrates with positively-charged residues and favors LAT and SLP-76 phosphosites that are surrounded by negatively-charged residues. This mechanism prevents ZAP-70 from phosphorylating its own activation loop, thereby enforcing its strict dependence on Lck for activation. The sequence features in ZAP-70, LAT, and SLP-76 that underlie electrostatic selectivity likely contribute to the specific response of T cells to foreign antigens.

59 BASIC BIOLOGICAL SCIENCES↗

Improving deep learning-based protein distance prediction in CASP14

Abstract Motivation Accurate prediction of residue–residue distances is important for protein structure prediction. We developed several protein distance predictors based on a deep learning distance prediction method and blindly tested them in the 14th Critical Assessment of Protein Structure Prediction (CASP14). The prediction method uses deep residual neural networks with the channel-wise attention mechanism to classify the distance between every two residues into multiple distance intervals. The input features for the deep learning method include co-evolutionary features as well as other sequence-based features derived from multiple sequence alignments (MSAs). Three alignment methods are used with multiple protein sequence/profile databases to generate MSAs for input feature generation. Based on different configurations and training strategies of the deep learning method, five MULTICOM distance predictors were created to participate in the CASP14 experiment. Results Benchmarked on 37 hard CASP14 domains, the best performing MULTICOM predictor is ranked 5th out of 30 automated CASP14 distance prediction servers in terms of precision of top L/5 long-range contact predictions [i.e. classifying distances between two residues into two categories: in contact (<8 Angstrom) and not in contact otherwise] and performs better than the best CASP13 distance prediction method. The best performing MULTICOM predictor is also ranked 6th among automated server predictors in classifying inter-residue distances into 10 distance intervals defined by CASP14 according to the precision of distance classification. The results show that the quality and depth of MSAs depend on alignment methods and sequence databases and have a significant impact on the accuracy of distance prediction. Using larger training datasets and multiple complementary features improves prediction accuracy. However, the number of effective sequences in MSAs is only a weak indicator of the quality of MSAs and the accuracy of predicted distance maps. In contrast, there is a strong correlation between the accuracy of contact/distance predictions and the average probability of the predicted contacts, which can therefore be more effectively used to estimate the confidence of distance predictions and select predicted distance maps. Availability and implementation The software package, source code and data of DeepDist2 are freely available at https://github.com/multicom-toolbox/deepdist and https://zenodo.org/record/4712084#.YIIM13VKhQM. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

Parallel Embeddings: a Visualization Technique for Contrasting Learned Representations

We introduce ``Parallel Embeddings'', a new technique that generalizes the classical Parallel Coordinates visualization technique from a sequence of features within a single matrix to a sequence of multiple feature matrices. This is accomplished by separately embedding each matrix in the sequence into a 2-D frame, juxtaposing these 2-D representations, then connecting points between frames. We outline how this visualization can be used for several ``model comparison'' tasks in machine learning. We compare user performance with Parallel Embeddings to a baseline using the TensorFlow Embedding Projector for estimating model accuracy and understanding conceptual model differences. We found that users were more accurate with Parallel Embeddings for model comparison and that users learned how to use Parallel Embeddings more quickly than TensorFlow Embedding Projector. Furthermore, users' analytical process using Embedding Projector was positively affected by using Parallel Embeddings beforehand.

Model comparison, machine learing, machine learnin↗

Mutations in the Second Alternative Oxidase Gene: A New Approach to Group Aspergillus niger Strains

Alternative oxidase is a terminal oxidase in the branched mitochondrial electron transport chain of most fungi including Aspergillus niger (subgenus Circumdati, section Nigri). A second, paralogous aox gene (aoxB) is extant in some A. niger isolates but also present in two divergent species of the subgenus Nidulantes—A. calidoustus and A. implicatus—as well as in Penicillium swiecickii. Black aspergilli are cosmopolitan opportunistic fungi that can cause diverse mycoses and acute aspergillosis in immunocompromised individuals. Amongst the approximately 75 genome-sequenced A. niger strains, aoxB features considerable sequence variation. Five mutations were identified that rationally affect transcription or function or terminally modify the gene product. One mutant allele that occurs in CBS 513.88 and A. niger neotype strain CBS 554.65 involves a chromosomal deletion that removes exon 1 and intron 1 from aoxB. Another aoxB allele results from retrotransposon integration. Three other alleles result from point mutations: a missense mutation of the start codon, a frameshift, and a nonsense mutation. A. niger strain ATCC 1015 has a full-length aoxB gene. The A. niger sensu stricto complex can thus be subdivided into six taxa according to extant aoxB allele, which may facilitate rapid and accurate identification of individual species

59 BASIC BIOLOGICAL SCIENCES↗

India National Gas Hydrate Program Expedition 02 Summary of Scientific Results: Numerical Simulation of Reservoir Response to Depressurization

The India National Gas Hydrate Program Expedition 02 (NGHP-02) discovered gas hydrate at high saturation in sand reservoirs at several sites in the deepwater Bay of Bengal. To assess the potential response of those deposits to scientific depressurization experiments, comprehensive geologic models were constructed to enable numerical simulation for two sites. Both sites (NGHP-02-09 and NGHP-02-16) feature thick sequences of thinly-interbedded reservoir and non-reservoir facies at sub-seafloor depths less than 300m and sub-sea depths of 2400m or more. These settings pose significant challenges to current modeling capabilities. First, the thinly-interbedded reservoir architecture complicates the determination of basic reservoir parameters from both log and core data due to measurement resolution issues. Secondly, the fine scale variation in sediment properties imparts great contrasts in key parameters over very short distances, creating high gradients at multiple scales and varying orientations that necessitate careful design of high-definition simulation grids. Thirdly, the deposits include internal sources of water, as well as a range of complex boundary conditions, including variable permeability within the overlying mud-rich “seals,” that complicate reservoir depressurization. Lastly, because of the unique combination of great water depth and relatively shallow sub-seafloor depth: models designed to maximize the dissociation rate impose large pressure drawdowns on relatively low-strength sediments. This condition renders the proper evaluation and integration of the geomechanical response to hydrate dissociation critical. In this report, we review the history of gas hydrate reservoir simulation, discuss methods for creating geologic input models, and summarize the key findings and implications of the collaborative NGHP-02 numerical simulation effort. Together, the studies confirm the viability of the modeled accumulations for scientific testing and identify key challenges related to the selection of specific test sites and the design of test wells.

Boswell, Ray↗

3D-equivariant graph neural networks for protein model quality assessment

Quality assessment (QA) of predicted protein tertiary structure models plays an important role in ranking and using them. With the recent development of deep learning end-to-end protein structure prediction techniques for generating highly confident tertiary structures for most proteins, it is important to explore corresponding QA strategies to evaluate and select the structural models predicted by them since these models have better quality and different properties than the models predicted by traditional tertiary structure prediction methods. We develop EnQA, a novel graph-based 3D-equivariant neural network method that is equivariant to rotation and translation of 3D objects to estimate the accuracy of protein structural models by leveraging the structural features acquired from the state-of-the-art tertiary structure prediction method—AlphaFold2. We train and test the method on both traditional model datasets (e.g. the datasets of the Critical Assessment of Techniques for Protein Structure Prediction) and a new dataset of high-quality structural models predicted only by AlphaFold2 for the proteins whose experimental structures were released recently. Our approach achieves state-of-the-art performance on protein structural models predicted by both traditional protein structure prediction methods and the latest end-to-end deep learning method—AlphaFold2. It performs even better than the model QA scores provided by AlphaFold2 itself. The results illustrate that the 3D-equivariant graph neural network is a promising approach to the evaluation of protein structural models. Integrating AlphaFold2 features with other complementary sequence and structural features is important for improving protein model QA.

59 BASIC BIOLOGICAL SCIENCES↗

Per-Phase Control for CHB Converters with Negative-Sequence Current and DC Voltage Balancing Control

This paper presents a per-phase controller for a grid-connected Cascaded H-Bridge (CHB) converter featuring negative sequence current regulation and DC bus voltage balancing capabilities. Unlike conventional three-phase controllers, each phase of the CHB can be separately and controlled to regulate the DC voltage and reactive power flow. per-phase negative sequence current control strategy is incorporated into the proposed controller to ensure compliance with IEEE Standard 2800-2022, which requires the absorption of negative-sequence reactive current during grid voltage disturbances. Furthermore, a feed-forward control-based DC bus voltage balancing method is employed, which eliminates the need for the integrator commonly used in traditional approaches. The performance of the proposed controller is validated under various unbalanced grid scenarios and DC bus voltage imbalances through electromagnetic transient (EMT) simulation using a three-phase, three-level grid-connected CHB–Dual Active Bridge (DAB) converter testbed developed on MATLAB/Simulink. Additionally, the performance of the negative-sequence current regulation is investigated through simulation under different unbalanced grid voltage scenarios.

negative sequence current↗

Better understanding and prediction of antiviral peptides through primary and secondary structure feature importance

The emergence of viral epidemics throughout the world is of concern due to the scarcity of available effective antiviral therapeutics. The discovery of new antiviral therapies is imperative to address this challenge, and antiviral peptides (AVPs) represent a valuable resource for the development of novel therapies to combat viral infection. We present a new machine learning model to distinguish AVPs from non-AVPs using the most informative features derived from the physicochemical and structural properties of their amino acid sequences. To focus on those features that are most likely to contribute to antiviral performance, we filter potential features based on their importance for classification. These feature selection analyses suggest that secondary structure is the most important peptide sequence feature for predicting AVPs. Our Feature-Informed Reduced Machine Learning for Antiviral Peptide Prediction (FIRM-AVP) approach achieves a higher accuracy than either the model with all features or current state-of-the-art single classifiers. Understanding the features that are associated with AVP activity is a core need to identify and design new AVPs in novel systems. The FIRM-AVP code and standalone software package are available at https://github.com/pmartR/FIRM-AVP with an accompanying web application at https://msc-viz.emsl.pnnl.gov/AVPR.

59 BASIC BIOLOGICAL SCIENCES↗

IMG/VR v4: an expanded database of uncultivated virus genomes within a framework of extensive functional, taxonomic, and ecological metadata

Viruses are widely recognized as critical members of all microbiomes. Metagenomics enables large-scale exploration of the global virosphere, progressively revealing the extensive genomic diversity of viruses on Earth and highlighting the myriad of ways by which viruses impact biological processes. IMG/VR provides access to the largest collection of viral sequences obtained from (meta)genomes, along with functional annotation and rich metadata. A web interface enables users to efficiently browse and search viruses based on genome features and/or sequence similarity. Here, for this work, we present the fourth version of IMG/VR, composed of >15 million virus genomes and genome fragments, a ≈6-fold increase in size compared to the previous version. These clustered into 8.7 million viral operational taxonomic units, including 231 408 with at least one high-quality representative. Viral sequences in IMG/VR are now systematically identified from genomes, metagenomes, and metatranscriptomes using a new detection approach (geNomad), and IMG standard annotation are complemented with genome quality estimation using CheckV, taxonomic classification reflecting the latest taxonomic standards, and microbial host taxonomy prediction. IMG/VR v4 is available at https://img.jgi.doe.gov/vr, and the underlying data are available to download at https://genome.jgi.doe.gov/portal/IMG_VR.

59 BASIC BIOLOGICAL SCIENCES↗

Prediction of Specificity of α-Conotoxins to Subtypes of Human Nicotinic Acetylcholine Receptors with Semi-supervised Machine Learning

Conotoxins are a family of highly toxic neurotoxins composed of cysteine-rich peptides produced by marine cone snails. The most lethal cone snail species to humans is Conus geographus, with fatality rates of up to ∼65% from a single sting, which is caused mostly by the activity of α-conotoxins against human nicotinic acetylcholine receptors (nAChRs). While sequence-based machine learning (ML) classifiers have been trained to identify targets of conotoxins binding voltage-gated ion channels, no ML model has been built to predict the subtype-specific nAChR targets of α-conotoxins. Here, we trained an ML model in a semi-supervised manner to predict the specificity of α-conotoxin binding toward different human nAChR subtypes to overcome the challenge of limited data in subtype-specific nAChR targets of α-conotoxins and the issue that one α-conotoxin can bind multiple nAChR subtypes with high selectivity. We considered additional features of sequences of α-conotoxins in training our ML model, including the secondary structure propensities and electrostatic properties, which resulted in better prediction capability for the ML model. Notably, we identify that most α-conotoxins bind to α3β2, α1γδ, and α7 subtypes of human nAChRs. Our findings from this study provide a framework for predicting targets of various kinds of toxins.

59 BASIC BIOLOGICAL SCIENCES↗

India National Gas Hydrate Program Expedition 02 Summary of Scientific Results: Gas Hydrate Systems Along the Eastern Continental Margin of India

The India National Gas Hydrate Program Expedition 02 (NGHP-02) discovered gas hydrate at high saturation in sand reservoirs at several sites in the deepwater Bay of Bengal. To assess the potential response of those deposits to scientific depressurization experiments, comprehensive geologic models were constructed to enable numerical simulation for two sites. Both sites (NGHP-02-09 and NGHP-02-16) feature thick sequences of thinly-interbedded reservoir and non-reservoir facies at sub-seafloor depths less than 300m and sub-sea depths of 2400m or more. These settings pose significant challenges to current modeling capabilities. First, the thinly-interbedded reservoir architecture complicates the determination of basic reservoir parameters from both log and core data due to measurement resolution issues. Secondly, the fine scale variation in sediment properties imparts great contrasts in key parameters over very short distances, creating high gradients at multiple scales and varying orientations that necessitate careful design of high-definition simulation grids. Thirdly, the deposits include internal sources of water, as well as a range of complex boundary conditions, including variable permeability within the overlying mud-rich “seals,” that complicate reservoir depressurization. Lastly, because of the unique combination of great water depth and relatively shallow sub-seafloor depth: models designed to maximize the dissociation rate impose large pressure drawdowns on relatively low-strength sediments. This condition renders the proper evaluation and integration of the geomechanical response to hydrate dissociation critical. In this report, we review the history of gas hydrate reservoir simulation, discuss methods for creating geologic input models, and summarize the key findings and implications of the collaborative NGHP-02 numerical simulation effort. Together, the studies confirm the viability of the modeled accumulations for scientific testing and identify key challenges related to the selection of specific test sites and the design of test wells.

Collett, Tim↗

G4Boost: a machine learning-based tool for quadruplex identification and stability prediction

Background: G-quadruplexes (G4s), formed within guanine-rich nucleic acids, are secondary structures involved in important biological processes. Although every G4 motif has the potential to form a stable G4 structure, not every G4 motif would, and accurate energy-based methods are needed to assess their structural stability. Here, we present a decision tree-based prediction tool, G4Boost, to identify G4 motifs and predict their secondary structure folding probability and thermodynamic stability based on their sequences, nucleotide compositions, and estimated structural topologies. Results: G4Boost predicted the quadruplex folding state with an accuracy greater then 93% and an F1-score of 0.96, and the folding energy with an RMSE of 4.28 and R 2 of 0.95 only by the means of sequence intrinsic feature. G4Boost was successfully applied and validated to predict the stability of experimentally-determined G4 structures, including for plants and humans. Conclusion: G4Boost outperformed the three machine-learning based prediction tools, DeepG4, Quadron, and G4RNA Screener, in terms of both accuracy and F1-score, and can be highly useful for G4 prediction to understand gene regulation across species including plants and humans.

59 BASIC BIOLOGICAL SCIENCES↗

A tale of two tools: mcres.py, a stochastic resonance generator, and grokres.py, a resonance quality assurance tool

We detail two software tools, now integrated into the fudge code system. The first tool, mcres.py, can be used to generate stochastic ensembles of resonances which are both consistent with the expectations of the Gaussian Orthogonal Ensemble of Random Matrix Theory and with the level densities and widths encoded in ENDF formatted files. The second tool, grokres.py, can be used to assess global and local features of sequences of resonances found in ENDF files and make comparisons to known results from Random Matrix Theory. We apply these tools to 54 Fe and other nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗