Search NASA⌕ Search

SEARCH · Search NASA

Results for “Protein”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Site specific N- and O-glycosylation mapping of the spike proteins of SARS-CoV-2 variants of concern

Abstract The glycosylation on the spike (S) protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, modulates the viral infection by altering conformational dynamics, receptor interaction and host immune responses. Several variants of concern (VOCs) of SARS-CoV-2 have evolved during the pandemic, and crucial mutations on the S protein of the virus have led to increased transmissibility and immune escape. In this study, we compare the site-specific glycosylation and overall glycomic profiles of the wild type Wuhan-Hu-1 strain (WT) S protein and five VOCs of SARS-CoV-2: Alpha, Beta, Gamma, Delta and Omicron. Interestingly, both N- and O-glycosylation sites on the S protein are highly conserved among the spike mutant variants, particularly at the sites on the receptor-binding domain (RBD). The conservation of glycosylation sites is noteworthy, as over 2 million SARS-CoV-2 S protein sequences have been reported with various amino acid mutations. Our detailed profiling of the glycosylation at each of the individual sites of the S protein across the variants revealed intriguing possible association of glycosylation pattern on the variants and their previously reported infectivity. While the sites are conserved, we observed changes in the N- and O-glycosylation profile across the variants. The newly emerged variants, which showed higher resistance to neutralizing antibodies and vaccines, displayed a decrease in the overall abundance of complex-type glycans with both fucosylation and sialylation and an increase in the oligomannose-type glycans across the sites. Among the variants, the glycosylation sites with significant changes in glycan profile were observed at both the N -terminal domain and RBD of S protein, with Omicron showing the highest deviation. The increase in oligomannose-type happens sequentially from Alpha through Delta. Interestingly, Omicron does not contain more oligomannose-type glycans compared to Delta but does contain more compared to the WT and other VOCs. O-glycosylation at the RBD showed lower occupancy in the VOCs in comparison to the WT. Our study on the sites and pattern of glycosylation on the SARS-CoV-2 S proteins across the VOCs may help to understand how the virus evolved to trick the host immune system. Our study also highlights how the SARS-CoV-2 virus has conserved both N - and O - glycosylation sites on the S protein of the most successful variants even after undergoing extensive mutations, suggesting a correlation between infectivity/ transmissibility and glycosylation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Top-down mass spectrometry and assigning internal fragments for determining disulfide bond positions in proteins

Disulfide bonds in proteins have a substantial impact on protein structure, stability, and biological activity. Localizing disulfide bonds is critical for understanding protein folding and higher-order structure. Conventional top-down mass spectrometry (TD-MS), where only terminal fragments are assigned for disulfide-intact proteins, can access disulfide information, but suffers from low fragmentation efficiency, thereby limiting sequence coverage. Here, we show that assigning internal fragments generated from TD-MS enhances the sequence coverage of disulfide-intact proteins by 20–60% by returning information from the interior of the protein sequence, which cannot be obtained by terminal fragments alone. Further, the inclusion of internal fragments can extend the sequence information of disulfide-intact proteins to near complete sequence coverage. Importantly, the enhanced sequence information that arise from the assignment of internal fragments can be used to determine the relative position of disulfide bonds and the exact disulfide connectivity between cysteines. The data presented here demonstrates the benefits of incorporating internal fragment analysis into the TD-MS workflow for analyzing disulfide-intact proteins, which would be valuable for characterizing biotherapeutic proteins such as monoclonal antibodies and antibody–drug conjugates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A gated graph transformer for protein complex structure quality assessment and its performance in CASP15

Abstract Motivation Proteins interact to form complexes to carry out essential biological functions. Computational methods such as AlphaFold-multimer have been developed to predict the quaternary structures of protein complexes. An important yet largely unsolved challenge in protein complex structure prediction is to accurately estimate the quality of predicted protein complex structures without any knowledge of the corresponding native structures. Such estimations can then be used to select high-quality predicted complex structures to facilitate biomedical research such as protein function analysis and drug discovery. Results In this work, we introduce a new gated neighborhood-modulating graph transformer to predict the quality of 3D protein complex structures. It incorporates node and edge gates within a graph transformer framework to control information flow during graph message passing. We trained, evaluated and tested the method (called DProQA) on newly-curated protein complex datasets before the 15th Critical Assessment of Techniques for Protein Structure Prediction (CASP15) and then blindly tested it in the 2022 CASP15 experiment. The method was ranked 3rd among the single-model quality assessment methods in CASP15 in terms of the ranking loss of TM-score on 36 complex targets. The rigorous internal and external experiments demonstrate that DProQA is effective in ranking protein complex structures. Availability and implementation The source code, data, and pre-trained models are available at https://github.com/jianlin-cheng/DProQA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

ForceGen: End-to-end de novo protein generation based on nonlinear mechanical unfolding responses using a language diffusion model

Through evolution, nature has presented a set of remarkable protein materials, including elastins, silks, keratins and collagens with superior mechanical performances that play crucial roles in mechanobiology. However, going beyond natural designs to discover proteins that meet specified mechanical properties remains challenging. Here, we report a generative model that predicts protein designs to meet complex nonlinear mechanical property-design objectives. Our model leverages deep knowledge on protein sequences from a pretrained protein language model and maps mechanical unfolding responses to create proteins. Via full-atom molecular simulations for direct validation, we demonstrate that the designed proteins are de novo, and fulfill the targeted mechanical properties, including unfolding energy and mechanical strength, as well as the detailed unfolding force-separation curves. Our model offers rapid pathways to explore the enormous mechanobiological protein sequence space unconstrained by biological synthesis, using mechanical features as the target to enable the discovery of protein materials with superior mechanical properties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Cell-free Scaled Production and Adjuvant Addition to a Recombinant Major Outer Membrane Protein from Chlamydia muridarum for Vaccine Development

Subunit vaccines offer advantages over more traditional inactivated or attenuated whole-cell-derived vaccines in safety, stability, and standard manufacturing. To achieve an effective protein-based subunit vaccine, the protein antigen often needs to adopt a native-like conformation. This is particularly important for pathogensurface antigens that are membrane-bound proteins. Cell-free methods have been successfully used to produce correctly folded functional membrane protein through the co-translation of nanolipoprotein particles (NLPs), commonly known as nanodiscs. This strategy can be used to produce subunit vaccines consisting of membrane proteins in a lipid-bound environment. However, cell-free protein production is often limited to small scale (<1 mL). The amount of protein produced in small-scale production runs is usually sufficient for biochemical and biophysical studies. However, the cell-free process needs to be scaled up, optimized, and carefully tested to obtain enough protein for vaccine studies in animal models. Other processes involved in vaccine production, such as purification, adjuvant addition, and lyophilization, need to be optimized in parallel. This paper reports the development of a scaled-up protocol to express, purify, and formulate a membrane-bound protein subunit vaccine

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1998

The Feasibility of Bulk Crystallization as an Industrial Purification and Production Technique for Proteins

Bulk crystallization in stirred vessels is used industrially for the recovery and purification of many inorganic and organic materials. Although much has been written on the crystallization of proteins for X-ray diffraction analysis, very little has been reported on the application of bulk crystallization in stirred vessels. In this study, a 1-liter, seeded, stirred, batch crystallizer was used with ovalbumin as a model protein to test the feasibility of this crystallization method as a recovery and purification process for proteins. Results were obtained for ovalbumin solubility, nucleation thresholds, crystal breakage and crystal growth kinetics in bulk solution under a range of operating conditions of pH and ammonium sulphate concentration (Judge et al., 1996). Experiments were also performed to determine the degree of purification that can be achieved by the crystallization of ovalbumin from a mixture of proteins. The effect of the presence of these proteins upon the ovalbumin crystal growth kinetics was also investigated (Judge et al., 1995). All of these aspects are essential for the design of bulk crystallization processes which have not previously been reported for proteins. Results from a second study that investigated the effect of structurally different proteins on the solubility, crystal growth rates and crystal purity of chicken egg white lysozyme are also presented (Judge et al., 1997). In this case face growth rates were measured using lysozyme purified by liquid chromatography and the effect of the addition of specific protein impurities were observed on the (110) and (101) crystal faces. In these two studies the results are presented to show the feasibility and purifying ability of crystallization as a production process for proteins.

Judge, Russell A.↗
NASA NTRS2001

Distribution and Spectroscopy of Green Fluorescent Protein and Acyl-CoA: Cholesterol Acytransferase in Sf21 Insect Cells

Acyl-CoA: cholesterol acyltransferase (ACAT) is thought to significantly participate in the pathway of cholesterol esterification that underlies the pathology of artherosclerosis. This enzyme is a membrane protein known to be preferentially bound within the endoplasmic reticulum of mammalian cells, from which location it esterifies cholesterol derived from low density lipoprotein. Cultures of insect cells were separately infected with baculovirus containing the gene for green fluroescent protein (GFP) and with baculovirus containing tandem genes for GFP and ACAT. These infected cultures expressed GFP and the fusion protein GCAT, respectively, with maximum expression occurring on the fourth day after infection. Extraction of GFP- and of GCAT-expressing cells with urea and detergent resulted in recovery of fluorescent protein in aqueous solution. Fluorescence spectra at neutral pH were identical for both GFP and GCAT extracts in aqueous solution, indicating unperturbed tertiary structure for the GFP moiety within GCAT. In a cholesterol esterification assay, GCAT demonstrated ACAT activity, but with less efficiency compared to native ACAT. It was hypothesized that the membrane protein ACAT would lead to differences in localization of GCAT compared to GFP within the respective expressing insect cells. The GFP marker directly and also within the fusion protein GCAT was accordingly used as the intracellular probe that was fluorescently analyzed by the new biophotonics technique of hyperspectral imaging. In that technique, fluorescence imaging was obtained from two dimensional arrays of cells, and regions of interest from within those images were then retrospectively analyzed for the emission spectra that comprises the image. Results of hyperspectral imaging of insect cells on day 4 postinfection showed that GCAT was preferentially localized to the cytoplasm of these cells compared to GFP. Furthermore, the emission spectra obtained for the localized GCAT displayed a peak blue shift from 518nm obtained in neutral aqueous solution to 505nm obtained in localized regions within the cells. This blue shift indicates change in the fluorescence coupling of the GFP moiety of GCAT. It is hypothesized that change in tertiary environment of GCAT, coincident with intracellular deposition of GCAT, follows from intracellular trafficking of GCAT leading to membrane interactions with the ACAT moiety, and/or self-assembly of GCAT, that alters the chromophore environment of the GFP moiety of GCAT. These findings introduce a new technique of biophotonic imaging to studies of intracellular protein trafficking and interactions. This technique of hyperspectral imaging could contribute to advancing the emergent field of proteomics. Because of the noninvasive nature of this technique, kinetic processes associated with intracellular protein trafficking, and interactions of proteins within cellular domains, can be considered for investigation within a single cell as well as a cell population.

Richmond, R. C.↗
NASA NTRS2003

Non-Genomic Origins of Proteins and Metabolism

It is proposed that evolution of inanimate matter to cells endowed with a nucleic acid- based coding of genetic information was preceded by an evolutionary phase, in which peptides not coded by nucleic acids were able to self-organize into networks capable of evolution towards increasing metabolic complexity. Recent findings that truly different, simple peptides (Keefe and Szostak, 2001) can perform the same function (such as ATP binding) provide experimental support for this mechanism of early protobiological evolution. The central concept underlying this mechanism is that the reproduction of cellular functions alone was sufficient for self-maintenance of protocells, and that self- replication of macromolecules was not required at this stage of evolution. The precise transfer of information between successive generations of the earliest protocells was unnecessary and, possibly, undesirable. The key requirement in the initial stage of protocellular evolution was an ability to rapidly explore a large number of protein sequences in order to discover a set of molecules capable of supporting self- maintenance and growth of protocells. Undoubtedly, the essential protocellular functions were carried out by molecules not nearly as efficient or as specific as contemporary proteins. Many, potentially unrelated sequences could have performed each of these functions at an evolutionarily acceptable level. As evolution progressed, however proteins must have performed their functions with increasing efficiency and specificity. This, in turn, put additional constraints on protein sequences and the fraction of proteins capable of performing their functions at the required level decreased. At some point, the likelihood of generating a sufficiently efficient set of proteins through a non-coded synthesis was so small that further evolution was not possible without storing information about the sequences of these proteins. Beyond this point, further evolution required coupling between proteins and informational polymers that is characteristic to all known forms of life. The emergence of such coupling must be postulated in any scenario of the origin of life, no matter whether it starts with RNA or proteins. To examine the evolutionary potential of non-genomic systems, a simple, computationally tractable model, which is still capable of capturing the essential features of the real system, has been studied computationally. Both constructive and destructive processes have been introduced into the model in a stochastic manner. Instead of assuming random reaction sets, only a suite of protobiologically plausible reactions has been considered. Peptides have been explicitly considered as protoenzymes and their catalytic efficiencies have been assigned on the basis of biochemical principles and experimental estimates. Simulations have been carried out using a novel approach (The Next Reaction Method) that is appropriate even for very low concentrations of reactants. Studies have focused on global autocatalytic processes and their diversity.

Pohorille, Andrew↗
NASA NTRS2000

A calmodulin binding protein from Arabidopsis is induced by ethylene and contains a DNA-binding motif

Calmodulin (CaM), a key calcium sensor in all eukaryotes, regulates diverse cellular processes by interacting with other proteins. To isolate CaM binding proteins involved in ethylene signal transduction, we screened an expression library prepared from ethylene-treated Arabidopsis seedlings with 35S-labeled CaM. A cDNA clone, EICBP (Ethylene-Induced CaM Binding Protein), encoding a protein that interacts with activated CaM was isolated in this screening. The CaM binding domain in EICBP was mapped to the C-terminus of the protein. These results indicate that calcium, through CaM, could regulate the activity of EICBP. The EICBP is expressed in different tissues and its expression in seedlings is induced by ethylene. The EICBP contains, in addition to a CaM binding domain, several features that are typical of transcription factors. These include a DNA-binding domain at the N terminus, an acidic region at the C terminus, and nuclear localization signals. In database searches a partial cDNA (CG-1) encoding a DNA-binding motif from parsley and an ethylene up-regulated partial cDNA from tomato (ER66) showed significant similarity to EICBP. In addition, five hypothetical proteins in the Arabidopsis genome also showed a very high sequence similarity with EICBP, indicating that there are several EICBP-related proteins in Arabidopsis. The structural features of EICBP are conserved in all EICBP-related proteins in Arabidopsis, suggesting that they may constitute a new family of DNA binding proteins and are likely to be involved in modulating gene expression in the presence of ethylene.

Non-NASA Center↗
NASA NTRS1998

The N-end rule pathway catalyzes a major fraction of the protein degradation in skeletal muscle

In skeletal muscle, overall protein degradation involves the ubiquitin-proteasome system. One property of a protein that leads to rapid ubiquitin-dependent degradation is the presence of a basic, acidic, or bulky hydrophobic residue at its N terminus. However, in normal cells, substrates for this N-end rule pathway, which involves ubiquitin carrier protein (E2) E214k and ubiquitin-protein ligase (E3) E3alpha, have remained unclear. Surprisingly, in soluble extracts of rabbit muscle, we found that competitive inhibitors of E3alpha markedly inhibited the 125I-ubiquitin conjugation and ATP-dependent degradation of endogenous proteins. These inhibitors appear to selectively inhibit E3alpha, since they blocked degradation of 125I-lysozyme, a model N-end rule substrate, but did not affect the degradation of proteins whose ubiquitination involved other E3s. The addition of several E2s or E3alpha to the muscle extracts stimulated overall proteolysis and ubiquitination, but only the stimulation by E3alpha or E214k was sensitive to these inhibitors. A similar general inhibition of ubiquitin conjugation to endogenous proteins was observed with a dominant negative inhibitor of E214k. Certain substrates of the N-end rule pathway are degraded after their tRNA-dependent arginylation. We found that adding RNase A to muscle extracts reduced the ATP-dependent proteolysis of endogenous proteins, and supplying tRNA partially restored this process. Finally, although in muscle extracts the N-end rule pathway catalyzes most ubiquitin conjugation, it makes only a minor contribution to overall protein ubiquitination in HeLa cell extracts.

Non-NASA Center↗
NASA NTRS1992

Cloning of the cDNA for U1 small nuclear ribonucleoprotein particle 70K protein from Arabidopsis thaliana

We cloned and sequenced a plant cDNA that encodes U1 small nuclear ribonucleoprotein (snRNP) 70K protein. The plant U1 snRNP 70K protein cDNA is not full length and lacks the coding region for 68 amino acids in the amino-terminal region as compared to human U1 snRNP 70K protein. Comparison of the deduced amino acid sequence of the plant U1 snRNP 70K protein with the amino acid sequence of animal and yeast U1 snRNP 70K protein showed a high degree of homology. The plant U1 snRNP 70K protein is more closely related to the human counter part than to the yeast 70K protein. The carboxy-terminal half is less well conserved but, like the vertebrate 70K proteins, is rich in charged amino acids. Northern analysis with the RNA isolated from different parts of the plant indicates that the snRNP 70K gene is expressed in all of the parts tested. Southern blotting of genomic DNA using the cDNA indicates that the U1 snRNP 70K protein is coded by a single gene.

NASA Discipline Plant Biology↗
NASA NTRS2012

High Protein Intake Improves Insulin Sensitivity but Exacerbates Bone Resorption in Immobility (WISE Study)

Inactivity, like bed rest (BR), causes insulin resistance (IR) and bone loss even in healthy subjects. High protein intake seems to mitigate this IR but might exacerbate bone loss. We hypothesized that high protein intake (animal:vegetable protein ratio: 60:40), isocaloric, compared to the control group plus high potassium intake would prevent IR without affecting bone turnover. After a 20-day ambulatory adaptation to controlled confinement and diet, 16 women participated in a 60-day, 6 deg head-down-tilt BR and were assigned randomly to one of the two groups. Control subjects (CON, n=8) received 1g/kg body mass/d dietary protein. Nutrition subjects (NUT, n=8) received 1.45g/kg body mass/d dietary protein plus 7.2g branched chain amino acids per day during BR. All subjects received 1670 kcal/d. Bed rest decreased glucose disposal by 35% (p<0.05) in CON. Isocaloric high protein intake prevented insulin resistance, but exacerbated bed rest induced increase in bone resorption markers C-telopeptide (> 30%) and Ntelopeptide (>20%) (both: p<0.001). Bone formation markers were unaffected by high protein intake. We conclude from these results that high protein intake might positively affect glucose tolerance, but might also foster bone loss. Further long-duration studies are mandatory before high protein intake for diabetic patients, who have an increased fracture risk, might be recommended.

Heer, Martina↗
DOE OSTI2024

Unveiling the nanoscale architectures and dynamics of protein assembly with in situ atomic force microscopy

Proteins play a vital role in different biological processes by forming complexes through precise folding with exclusive inter- and intra-molecular interactions. Understanding the structural and regulatory mechanisms underlying protein complex formation provides insights into biophysical processes. Furthermore, the principle of protein assembly gives guidelines for new biomimetic materials with potential applications in medicine, energy, and nanotechnology. Atomic force microscopy (AFM) is a powerful tool for investigating protein assembly and interactions across spatial scales (single molecules to cells) and temporal scales (milliseconds to days). It has significantly contributed to understanding nanoscale architectures, inter- and intra-molecular interactions, and regulatory elements that determine protein structures, assemblies, and functions. This review describes recent advancements in elucidating protein assemblies with in situ AFM. We discuss the structures, diffusions, interactions, and assembly dynamics of proteins captured by conventional and high-speed AFM in near-native environments and recent AFM developments in the multimodal high-resolution imaging, bimodal imaging, live cell imaging, and machine-learning-enhanced data analysis. These approaches show the significance of broadening the horizons of AFM and enable unprecedented explorations of protein assembly for biomaterial design and biomedical research.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structural evolution during gelation of pea and whey proteins envisaged by time-resolved ultra-small-angle x-ray scattering (USAXS)

Hydrogels from plant proteins commonly exhibit inferior gel strength compared to those from dairy proteins partially due to their distinct gel networks. How protein aggregates to form such networks in response to heat remains largely unknown. In this work, pea (PPI) and whey (WPI) protein isolate gels were produced at the same protein content and similar heating/cooling rate. The process was monitored using rheology, microscopy and in situ ultra-small-angle x-ray scattering (USAXS). Rheology showed an initial decrease in G' and G'' in PPI followed by a steady increase when the temperature surpassed ~60 °C whereas a much higher temperature (~80 °C) was required for WPI, both using 2 °C/min heating rate. Microscopy showed a coarse and heterogenous network in PPI, whereas for WPI, the network was finer and more continuous. In both gels, nano-sized spherical or ellipsoidal particles were present as the basic constituents. USAXS found individual protein was dominant in PPI or WPI solution at temperature below 57 °C. Their proportions decreased together with appearance of aggregates with an average R g of 9–10 nm in PPI and 6–7 nm in WPI at higher temperature. The size of the aggregates changed slightly during further heating and cooling, but their proportions increased. Power law exponents revealed the aggregates were mass fractals for WPI and PPI gels, and they became more compact during heating. Our findings suggested formation of primary aggregates in protein gel networks is a more organized process and provided theoretical guidance for production of high protein food gels with desirable texture.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Protein Assembly by Design

Proteins are nature’s primary building blocks for the construction of sophisticated molecular machines and dynamic materials, ranging from protein complexes such as photosystem II and nitrogenase that drive biogeochemical cycles to cytoskeletal assemblies and muscle fibers for motion. Such natural systems have inspired extensive efforts in the rational design of artificial protein assemblies in the last two decades. As molecular building blocks, proteins are highly complex, in terms of both their three-dimensional structures and chemical compositions. To enable control over the self-assembly of such complex molecules, scientists have devised many creative strategies by combining tools and principles of experimental and computational biophysics, supramolecular chemistry, inorganic chemistry, materials science, and polymer chemistry, among others. Owing to these innovative strategies, what started as a purely structure-building exercise two decades ago has, in short order, led to artificial protein assemblies with unprecedented structures and functions and protein-based materials with unusual properties. Furthermore, our goal in this review is to give an overview of this exciting and highly interdisciplinary area of research, first outlining the design strategies and tools that have been devised for controlling protein self-assembly, then describing the diverse structures of artificial protein assemblies, and finally highlighting the emergent properties and functions of these assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Hotspot Coevolution Is a Key Identifier of Near-Native Protein Complexes

Protein–protein interactions play a key role in mediating numerous biological functions, with more than half the proteins in living organisms existing as either homo- or hetero-oligomeric assemblies. Protein subunits that form oligomers minimize the free energy of the complex, but exhaustive computational search-based docking methods have not comprehensively addressed the challenge of distinguishing a natively bound complex from non-native forms. Current protein docking approaches address this problem by sampling multiple binding modes in proteins and scoring each mode, with the lowest-energy (or highest scoring) binding mode being regarded as a near-native complex. However, high-scoring modes often match poorly with the true bound form, suggesting a need for improvement of the scoring function. Here in this study, we propose a scoring function, KFC-E, that accounts for both conservation and coevolution of putative binding hotspot residues at protein–protein interfaces. We tested KFC-E on four benchmark sets of unbound examples and two benchmark sets of bound examples, with the results demonstrating a clear improvement over scores that examine conservation and coevolution across the entire interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The protein structurome of Orthornavirae and its dark matter

Metatranscriptomics is uncovering more and more diverse families of viruses with RNA genomes comprising the viral kingdom Orthornavirae in the realm Riboviria. Thorough protein annotation and comparison are essential to get insights into the functions of viral proteins and virus evolution. In addition to sequence- and hmm profile-based methods, protein structure comparison adds a powerful tool to uncover protein functions and relationships. We constructed an Orthornavirae “structurome” consisting of already annotated as well as unannotated (“dark matter”) proteins and domains encoded in viral genomes. We used protein structure modeling and similarity searches to illuminate the remaining dark matter in hundreds of thousands of orthornavirus genomes. The vast majority of the dark matter domains showed either “generic” folds, such as single α-helices, or no high confidence structure predictions. Nevertheless, a variety of lineage-specific globular domains that were new either to orthornaviruses in general or to particular virus families were identified within the proteomic dark matter of orthornaviruses, including several predicted nucleic acid-binding domains and nucleases. In addition, we identified a case of exaptation of a cellular nucleoside monophosphate kinase as an RNA-binding protein in several virus families. Notwithstanding the continuing discovery of numerous orthornaviruses, it appears that all the protein domains conserved in large groups of viruses have already been identified. The rest of the viral proteome seems to be dominated by poorly structured domains including intrinsically disordered ones that likely mediate specific virus-host interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression. Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to previous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding. Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

59 BASIC BIOLOGICAL SCIENCES↗