Search NASA⌕ Search

SEARCH · Search NASA

Results for “Amino acid sequence analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Molecular modeling of calmodulin: a comparison with crystallographic data

Two methods of side-chain placement on a modeled protein have been examined. Two molecular models of calmodulin were constructed that differ in the treatment of side chains prior to optimization of the molecule. A virtual bond analysis program developed by Purisima and Scheraga was used to determine the backbone conformation based on 2.2 angstroms resolution C alpha coordinates for the molecules. In the first model, side chains were initially constructed in an extended conformation. In the second model, a conformational grid search technique was employed. Calcium ions were treated explicitly during energy optimization using CHARMM. The models are compared to a recently published refined crystal structure of calmodulin. The results indicate that the initial choices for side-chains, but also significant effects on the main-chain conformation and supersecondary structure. The conformational differences are discussed. Analysis of these and other methods makes possible the formulation of a methodology for more appropriate side-chain placement in modeled proteins.

NASA Discipline Exobiology↗
DOE OSTI2025

Cholesterol-dependent enzyme activity of human TSPO1

The amino acid sequence of the tryptophan-rich sensory proteins (TSPO) is substantially conserved throughout all kingdoms of life. Human mitochondrial TSPO1 (HsTSPO1) binds to porphyrins and steroids, although its interactions with these molecules remains unknown.HsTSPO1 is associated with numerous physiological and pathological disorders, but the underlying molecular mechanisms are unknown. Here, we disclose the finding of human mitochondrial TSPO as a cholesterol-dependent protoporphyrin IX oxygenase. The results of our biochemical characterization are consistent with structural data and evolutionary analysis. The dependence ofHsTSPO1 activity on cholesterol may be the result of the coevolution of this membrane protein with the membrane system. Our study provides a molecular foundation for comprehending the various roles played by mitochondrial TSPO in normal physiological and pathological situations.

Science & Technology - Other Topics↗
NASA NTRS2000

An early ethylene up-regulated gene encoding a calmodulin-binding protein involved in plant senescence and death

35S-Labeled calmodulin (CaM) was used to screen a tobacco anther cDNA library. A positive clone (NtER1) with high homology to an early ethylene-up-regulated gene (ER66) in tomato, and an Arabidopsis homolog was isolated and characterized. Based on the helical wheel projection, a 25-mer peptide corresponding to the predicted CaM-binding region of NtER1 (amino acids 796-820) was synthesized. The gel-mobility shift assay showed that the peptide formed a stable complex with CaM only in the presence of Ca(2+). CaM binds to NtER1 with high affinity (K(d) approximately 12 nm) in a calcium-dependent manner. Tobacco flowers at different stages of development were treated with ethylene or with 1-methylcyclopropene for 2 h before treating with ethylene. Northern analysis showed that the NtER1 was rapidly induced after 15 min of exposure to ethylene. However, the 2-h 1-methylcyclopropene treatment totally blocked NtER1 expression in flowers at all stages of development, suggesting that NtER1 is an early ethylene-up-regulated gene. The senescing leaves and petals had significantly increased NtER1 induction as compared with young leaves and petals, implying that NtER1 is developmentally regulated and acts as a trigger for senescence and death. This is the first documented evidence for the involvement of Ca(2+)/CaM-mediated signaling in ethylene action.

Non-NASA Center↗
DOE OSTI2020

Uncovering DNA-PKcs ancient phylogeny, unique sequence motifs and insights for human disease

DNA-dependent protein kinase catalytic subunit (DNA-PKcs) is a key member of the phosphatidylinositol-3 kinase-like (PIKK) family of protein kinases with critical roles in DNA-double strand break repair, transcription, metastasis, mitosis, RNA processing, and innate and adaptive immunity. The absence of DNA-PKcs from many model organisms has led to the assumption that DNA-PKcs is a vertebrate-specific PIKK. Here, we find that DNA-PKcs is widely distributed in invertebrates, fungi, plants, and protists, and that threonines 2609, 2638, and 2647 of the ABCDE cluster of phosphorylation sites are highly conserved amongst most Eukaryotes. Furthermore, we identify highly conserved amino acid sequence motifs and domains that are characteristic of DNA-PKcs relative to other PIKKs. These include residues in the Forehead domain and a novel motif we have termed YRPD, located in an α helix C-terminal to the ABCDE phosphorylation site loop. Combining sequence with biochemistry plus structural data on human DNA-PKcs unveils conserved sequence and conformational features with functional insights and implications. The defined generally progressive DNA-PKcs sequence diversification uncovers conserved functionality supported by Evolutionary Trace analysis, suggesting that for many organisms both functional sites and evolutionary pressures remain identical due to fundamental cell biology. The mining of cancer genomic data and germline mutations causing human inherited disease reveal that robust DNA-PKcs activity in tumors is detrimental to patient survival, whereas germline mutations compromising function are linked to severe immunodeficiency and neuronal degeneration. We anticipate that these collective results will enable ongoing DNA-PKcs functional analyses with biological and medical implications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Metagenome-assembled genomes from topsoils collected during NEON campaign in East River, CO (06/14/2018-06/28/2018)

The Watershed Function Science Focus Area (WF SFA) at Lawrence Berkeley National Lab is working to build a mechanistic understanding of the distribution and dynamics of biogeochemical processes in mountainous watersheds and their response to perturbation. In June 2018, the NEON (National Ecological Observatory Network) Airborne Observatory Platform (AOP) performed a taskable airborne imaging campaign to collect visible to shortwave infrared (VSWIR) imaging spectroscopy and LiDAR data across 330 km2 in the Upper East River at Crested Butte, CO. We conducted a parallel ground sampling campaign to sample vegetation traits, as well as soil physical, chemical, and microbiological characteristics. We collected these samples from 438 sites across 12 locations spanning much of the elevation, topographic, and geologic variability across the study area. A subset of 250 samples were used for soil metagenomics which is presented here. In addition, at each site, vegetation samples were collected to measure species-specific leaf water content and leaf mass area, foliar elemental composition and foliar CN stable isotope ratios. Soil samples were collected to measure soil physical properties which include bulk density and soil texture analysis. A suite of soil chemical properties was measured from the samples collected at each site, including pH, organic matter, concentrations exchangeable cations, total elemental composition, and the concentrations of extractable N pools (e.g. total free amino acids, ammonium, nitrate, dissolved organic N, and total dissolved N). Additionally, we have measured soil microbial biomass CN stoichiometry. Here, we present 1982 metagenome-assembled genomes (MAGs) for the bacterial and archaeal community from topsoil collected from during NEON 2018 campaign. All metagenomes were sequenced at JGI (Joint Genome Institute) (GOLD Study ID: Gs0149986). Metagenomes were assembled using JGI Metagenome Workflow (10.1128/mSystems.00804-20). The dataset includes (1) zip files for 1982 MAG fasta files (neon_genomes1-5.tar.gz, split into 5 tarballs to keep tarballs under 0.5 GB), (2) neon_Gs0149986_samples_soilproperties_metagenomes.csv: the sample information together with the accession numbers for the underlying metagenomes and the associated soil physical and chemical measurements in NMDC (National Microbiome Data Collaborative) compliant format, (3) neon_Gs0149986.kml: location bounding box file for the sampled locations, (4) samples.csv: sample metadata file used to register Internationall Generic Sample Numbers (IGSNs), (5) flmd.csv: file level metadata file, and (6) dd.csv: data dictionary file. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2022

A mutagenic analysis of NahE, a hydratase-aldolase in the naphthalene degradative pathway

NahE is a hydratase-aldolase that converts o-substituted trans-benzylidenepyruvates (H, OH, or CO 2 - ) to benzaldehyde, salicylaldehyde, or 2-carboxybenzaldehyde, respectively, and pyruvate. The enzyme is in a bacterial degradative pathway for naphthalene, which is a toxic and persistent environmental contaminant. Sequence, crystallographic, and mutagenic analysis identified the enzyme as a member of the N-acetylneuraminate lyase (NAL) subgroup in the aldolase superfamily. As such, it has a conserved lysine (Lys183) and tyrosine (Tyr155), for Schiff base formation, as well as a GXXGE motif for binding of the pyruvoyl carboxylate group. A crystal structure of the selenomethionine derivative of NahE shows these active site elements along with nearby residues that might be involved in the mechanism and/or specificity. Mutations of five active site amino acids (Thr65, Trp128, Tyr155, Asn157, and Asn281) were constructed and kinetic parameters measured in order to assess the effect(s) on catalysis. Here, the results show that the two Trp128 mutants (Phe and Tyr) have the least effect on catalysis, whereas amino acids with bulky side chains at Thr65 (Val) and Asn281 (Leu) have the greatest effect. Changing Tyr155 to Phe and Asn157 to Ala also hinders catalysis, and the effects fall in between these extremes. These observations are put into a structural context using a crystal structure of the Schiff base of the reaction intermediate. Trapping experiments with substrate, Na(CN)BH 3 , and wild type enzyme and selected mutants mostly paralleled the kinetic analysis, and identified two salicylaldehyde-modified lysines: the active site lysine (Lys183) and one outside the active site (Lys279). The latter could be responsible for the observed inhibition of NahE by salicylaldehyde. Together, the results provide new insights into the NahE-catalyzed reaction.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1999

Mechanosensitive channels in bacteria as membrane tension reporters

The purpose of this short review is to discuss recent data on the molecular structure and mechanism of gating of MscL, a mechanosensitive channel of large conductance from Escherichia coli. MscL is the first isolated molecule shown to convert mechanical stress of the membrane into a simple response, the opening of a large aqueous pore. The functional complex appears to be a stable homo-pentamer of 15-kDa subunits, the gating transitions in which are driven by stretch forces conveyed through the lipid bilayer. We have measured the open probability of MscL and the kinetics of transitions as a function of membrane tension. The parameters extracted from the single-channel current recordings and dose-response curves such as the energy difference between the closed, open, and intermediate conducting states, and the transition-related changes in protein dimensions suggest a large conformational rearrangement of the channel complex. The estimations show that in native conditions MscL openings could be driven primarily by forces of osmotic nature. The thermodynamic and spatial parameters reasonably correlate with the available data on the structure of a single MscL subunit and multimeric organization of the complex. Combined with the functional analysis of mutations, these data give grounds to hypotheses on the nature of the channel mechanosensitivity.

Review↗
DOE OSTI2023

National Prevalence of Caprine Prion Protein Genetic Variability at Codons 146, 211, and 222 in Goat Herds in the United States

Scrapie is a neurodegenerative disease that impacts sheep and goats, characterized by gradual and progressive changes in neurological function. Recent research shows that the scrapie incubation period is significantly influenced by specific variations in amino acids within the prion protein gene (PRNP). The objective of this study was to estimate the national prevalence of caprine PRNP genetic variability at codons 146, 211, and 222 in goat populations across the United States. A total of 3052 blood, ear tissue, and brain tissue samples were collected from goats from 50 states. The participating states were categorized into four Veterinary Service (VS) district regions. The samples underwent DNA extraction, and the PRNP variants corresponding to codons 146, 211, and 222 were amplified and sequenced. The analysis of PRNP variants, when compared to the PRNP reference sequence, revealed seven alleles in twelve genotypes. The homozygous 146NN, 211RR, and 222QQ alleles, which have been linked to an increased risk of scrapie, were found to be the most prevalent among all the goats. The heterozygous 222QK, 211RQ, 146SD, 146ND, and 146NS alleles and the homozygous 222KK, 146SS, and 146DD alleles, known to be associated with reduced scrapie susceptibility and a prolonged incubation period after experimental challenge, were found in 1.098% (222QK), 2.33% (211RQ), 0.58% (146SD), 3.13% (146ND), 20.68% (146NS), 0.005% (222KK), 3.31% (146SS), and 0.67% (146DD) of goats, respectively. The 222QK allele was found most frequently in goats tested from the east (VS District 1, 1.59%) and southwest (VS District 4, 1.08%) regions, whereas the 211RQ allele was found most often in goats tested from the Midwest (VS District 2, 8.03%) and east (VS District 1, 6.53%) regions. The 146NS allele was found most frequently in goats tested from the northwest (VS District 3, 29.02%) and southwest (VS District 4, 20.69%) regions. Our results showed that the prevalence of less susceptible genotypes at PRNP codon 146 may be sufficient to use genetic susceptibility testing in some herds. This may reduce the number of goats removed as part of a herd clean-up plan and may promote the selective breeding goats for less susceptible alleles in high-risk herds at the national level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Learning protein fitness models from evolutionary and assay-labeled data

Machine learning-based models of protein fitness typically learn from either unlabeled, evolutionarily related sequences or variant sequences with experimentally measured labels. For regimes where only limited experimental data are available, recent work has suggested methods for combining both sources of information. Toward that goal, we propose a simple combination approach that is competitive with, and on average outperforms more sophisticated methods. Our approach uses ridge regression on site-specific amino acid features combined with one probability density feature from modeling the evolutionary data. Within this approach, we find that a variational autoencoder-based probability density model showed the best overall performance, although any evolutionary density model can be used. Moreover, our analysis highlights the importance of systematic evaluations and sufficient baselines.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Diploid-dominant life cycles characterize the early evolution of Fungi

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane– and cell-wall–associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Simultaneous carbon catabolite repression governs sugar and aromatic co-utilization in Pseudomonas putida M2

ABSTRACT Pseudomonas putida have emerged as promising biocatalysts for the conversion of sugars and aromatic compounds obtained from lignocellulosic biomass. Understanding the role of carbon catabolite repression (CCR) in these strains is critical to optimize biomass conversion to fuels and chemicals. The CCR functioning in P. putida M2, a strain capable of consuming both hexose and pentose sugars as well as aromatic compounds, was investigated by cultivation experiments, proteomics, and CRISPRi-based gene repression. Strain M2 co-utilized sugars and aromatic compounds simultaneously; however, during cultivation with glucose and aromatic compounds ( p- coumarate and ferulate) mixture, intermediates (4-hydroxybenzoate and vanillate) accumulated, and substrate consumption was incomplete. In contrast, xylose-aromatic consumption resulted in transient intermediate accumulation and complete aromatic consumption, while xylose was incompletely consumed. Proteomics analysis revealed that glucose exerted stronger repression than xylose on the aromatic catabolic proteins. Key glucose (Eda) and xylose (XylX) catabolic proteins were also identified at lower abundance during cultivation with aromatic compounds implying simultaneous catabolite repression by sugars and aromatic compounds. Reduction of crc expression via CRISPRi led to faster growth and glucose and p -coumarate uptake in the CRISPRi strains compared to the control, while no difference was observed on xylose+ p -coumarate. The increased abundances of Eda and amino acid biosynthesis proteins in the CRISPRi strain further supported these observations. Lastly, small RNAs (sRNAs) sequencing results showed that CrcY and CrcZ homologues levels in M2, previously identified in P. putida strains, were lower under strong CCR (glucose+ p -coumarate) condition compared to when repression was absent ( p -coumarate or glucose only). IMPORTANCE A newly isolated Pseudomonas putida strain, P. putida M2, can utilize both hexose and pentose sugars as well as aromatic compounds making it a promising host for the valorization of lignocellulosic biomass. Pseudomonads have developed a regulatory strategy, carbon catabolite repression, to control the assimilation of carbon sources in the environment. Carbon catabolite repression may impede the simultaneous and complete metabolism of sugars and aromatic compounds present in lignocellulosic biomass and hinder the development of an efficient industrial biocatalyst. This study provides insight into the cellular physiology and proteome during mixed-substrate utilization in P. putida M2. The phenotypic and proteomics results demonstrated simultaneous catabolite repression in the sugar-aromatic mixtures, while the CRISPRi and sRNA sequencing demonstrated the potential role of the crc gene and small RNAs in carbon catabolite repression.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Cross-reactive immunogenicity of group A streptococcal vaccines designed using a recurrent neural network to identify conserved M protein linear epitopes

The M protein of group A streptococci (Strep A) is a major virulence determinant and protective antigen. The N-terminal sequence of the protein defines the more than 200 M types of Strep A and also contains epitopes that elicit opsonic antibodies, some of which cross-react with heterologous M types. Current efforts to develop broadly protective M protein-based vaccines are directed at identifying potential cross-protective epitopes located in the N-terminal regions of cluster-related M proteins for use as vaccine antigens. In this study, we have used a comprehensive approach using the recurrent neural network ABCpred and IEDB epitope conservancy analysis tools to predict 16 residue linear B-cell epitopes from 117 clinically relevant M types of Strep A (~88% of global Strep A infections). Furthermore, to examine the immunogenicity of these epitope-based vaccines, nine peptides that together shared ≥60% sequence identity with 37 heterologous M proteins were incorporated into two recombinant hybrid protein vaccines, in which the epitopes were repeated 2 or 3 times, respectively. The combined immune responses of immunized rabbits showed that the vaccines elicited significant levels of antibodies against all nine vaccine epitopes present in homologous N-terminal 1–50 amino acid synthetic M peptides, as well as cross-reactive antibodies against 16 of 37 heterologous M peptides predicted to contain similar epitopes. The epitope-specificity of the cross-reactive antibodies was confirmed by ELISA inhibition assays and functional opsonic activity was assayed in HL-60-based bactericidal assays. The results provide important information for the future design of broadly protective M protein-based Strep A vaccines.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Amino acid substrate specificities and tissue expression profiles of the nine CYP79A encoding genes in Sorghum bicolor

Cytochrome P450s of the CYP79 family catalyze two N-hydroxylation reactions, converting a selected number of amino acids into the corresponding oximes. The sorghum genome (Sorghum bicolor) harbours nine CYP79A encoding genes, and here sequence comparisons of the CYP79As along with their substrate recognition sites (SRSs) are provided. The substrate specificity of previously uncharacterized CYP79As was investigated by transient expression in Nicotiana benthamiana and subsequent transformation of the oximes formed into the corresponding stable oxime glucosides catalyzed by endogenous UDPG-glucosyltransferases (UGTs). CYP79A61 uses phenylalanine as a substrate, whereas CYP79A91, CYP79A93, and CYP79A95 use valine and isoleucine as substrates, with CYP79A93 showing the ability also to use phenylalanine. CYP79A94 uses isoleucine as a substrate. Analysis of 249 sorghum transcriptomes from two different sorghum cultivars showed the expression levels and tissue-specific expression of the CYP79As. CYP79A1 is the committed gene in dhurrin formation and was the highest expressed gene in most tissues/organs. CYP79A61 was primarily expressed in fully developed leaf blades and leaf sheaths. CYP79A91 and CYP79A92 were expressed mainly in roots >200 cm below ground, while CYP79A93 and CYP79A94 were most highly expressed in the leaf collar and leaf sheath, respectively. Here, the possible signalling effects of the oximes and their metabolites produced in different sorghum tissues are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1992

Development and use of domain-specific antibodies in a characterization of the large subunits of soybean photosystem 1

The molecular architecture of the soybean photosystem 1 reaction center complex was examined using a combination of surface labeling and immunological methodology on isolated thylakoid membranes. Synthetic peptides (12 to 14 amino acids in length) were prepared which correspond to the N-terminal regions of the 83 and 82.4 kDa subunits of photosystem 1 (the PsaA and PsaB proteins, respectively). Similarly, a synthetic peptide was prepared corresponding to the C-terminal region of the PsaB subunit. These peptides were conjugated to a carrier protein, and were used for the production of polyclonal antibodies in rabbits. The resulting sera could distinguish between the PsaA and PsaB photosystem 1 subunits by Western blot analysis, and could identify appropriate size classes of cyanogen bromide cleavage fragments as predicted from the primary sequences of these two subunits. When soybean thylakoid membranes were surface-labeled with N-hydroxysuccinimidobiotin, several subunits of the complete photosystem 1 lipid/protein complex incorporated label. These included the light harvesting chlorophyll proteins of photosystem 1, and peptides thought to aid in the docking of ferredoxin to the complex during photosynthetic electron transport. However, the PsaA and PsaB subunits showed very little biotinylation. When these subunits were examined for the domains to which biotin did attach, most of the observed label was associated with the N-terminal domain of the PsaA subunit, as identified using a domain-specific polyclonal antisera.

NASA Program NSCORT↗
DOE OSTI2023

Crystal structure of the Arabidopsis SPIRAL2 C-terminal domain reveals a p80-Katanin-like domain

Epidermal cells of dark-grown plant seedlings reorient their cortical microtubule arrays in response to blue light from a net lateral orientation to a net longitudinal orientation with respect to the long axis of cells. The molecular mechanism underlying this microtubule array reorientation involves katanin, a microtubule severing enzyme, and a plant-specific microtubule associated protein called SPIRAL2. Katanin preferentially severs longitudinal microtubules, generating seeds that amplify the longitudinal array. Upon severing, SPIRAL2 binds nascent microtubule minus ends and limits their dynamics, thereby stabilizing the longitudinal array while the lateral array undergoes net depolymerization. To date, no experimental structural information is available for SPIRAL2 to help inform its mechanism. To gain insight into SPIRAL2 structure and function, we determined a 1.8 Å resolution crystal structure of the Arabidopsis thaliana SPIRAL2 C-terminal domain. The domain is composed of seven core α-helices, arranged in an α-solenoid. Amino-acid sequence conservation maps primarily to one face of the domain involving helices α1, α3, α5, and an extended loop, the α6-α7 loop. The domain fold is similar to, yet structurally distinct from the C-terminal domain of Ge-1 (an mRNA decapping complex factor involved in P-body localization) and, surprisingly, the C-terminal domain of the katanin p80 regulatory subunit. The katanin p80 C-terminal domain heterodimerizes with the MIT domain of the katanin p60 catalytic subunit, and in metazoans, binds the microtubule minus-end factors CAMSAP3 and ASPM. Structural analysis predicts that SPIRAL2 does not engage katanin p60 in a mode homologous to katanin p80. The SPIRAL2 structure highlights an interesting evolutionary convergence of domain architecture and microtubule minus-end localization between SPIRAL2 and katanin complexes, and establishes a foundation upon which structure-function analysis can be conducted to elucidate the role of this domain in the regulation of plant microtubule arrays.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Distributed Many-to-Many Protein Sequence Alignment using Sparse Matrices

Identifying similar protein sequences is a core step in many computational biology pipelines such as detection of homologous protein sequences, generation of similarity protein graphs for downstream analysis, functional annotation, and gene location. Performance and scalability of protein similarity search have proven to be a bottleneck in many bioinformatics pipelines due to increase in cheap and abundant sequencing data. This work presents a new distributed-memory software PASTIS. PASTIS relies on sparse matrix computations for efficient identification of possibly similar proteins. We use distributed sparse matrices for scalability and show that the sparse matrix infrastructure is a great fit for protein similarity search when coupled with a fully-distributed dictionary of sequences that allow remote sequence requests to be fulfilled. Our algorithm incorporates the unique bias in amino acid sequence substitution in search without altering basic sparse matrix model, and in turn, achieves ideal scaling up to millions of protein sequences.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Molecular Architecture of Early Dissemination and Massive Second Wave of the SARS-CoV-2 Virus in a Major Metropolitan Area

We sequenced the genomes of 5,085 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains causing two coronavirus disease 2019 (COVID-19) disease waves in metropolitan Houston, TX, an ethnically diverse region with 7 million residents. The genomes were from viruses recovered in the earliest recognized phase of the pandemic in Houston and from viruses recovered in an ongoing massive second wave of infections. The virus was originally introduced into Houston many times independently. Virtually all strains in the second wave have a Gly614 amino acid replacement in the spike protein, a polymorphism that has been linked to increased transmission and infectivity. Patients infected with the Gly614 variant strains had significantly higher virus loads in the nasopharynx on initial diagnosis. We found little evidence of a significant relationship between virus genotype and altered virulence, stressing the linkage between disease severity, underlying medical conditions, and host genetics. Some regions of the spike protein—the primary target of global vaccine efforts—are replete with amino acid replacements, perhaps indicating the action of selection. We exploited the genomic data to generate defined single amino acid replacements in the receptor binding domain of spike protein that, importantly, produced decreased recognition by the neutralizing monoclonal antibody CR3022. Our report represents the first analysis of the molecular architecture of SARS-CoV-2 in two infection waves in a major metropolitan region. The findings will help us to understand the origin, composition, and trajectory of future infection waves and the potential effect of the host immune response and therapeutic maneuvers on SARS-CoV-2 evolution.

59 BASIC BIOLOGICAL SCIENCES↗