Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Protein 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Characterizing proteases in sorghum drought response using activity-based protein profiling (ABPP)

Sorghum (Sorghum bicolor) is a C4 grass species that is a major food and bioenergy crop grown worldwide. Understanding the response of sorghum to climate change is therefore critical to future sustainable food and bioenergy production. Major advancements in genome sequencing have provided a wealth of genomic data for DOE-relevant plant species, but the functions of many of these genes remain unknown. Genotype to phenotype characterization of sorghum and other bioenergy crops is needed to identify key genes responsible for plant resilience to drought and other environmental stressors. Proteases, which modulate protein degradation and cellular signaling processes by catalyzing the hydrolysis of proteins into smaller proteins, peptides, or amino acids, are highly implicated in plant drought stress responses. We developed a method for using activity-based protein profiling (ABPP) of selected protease activities (serine and cysteine proteases) to characterize drought stress responses in sorghum roots and leaves. Analysis of these ABPP data identified several known serine and cysteine proteases as well as uncharacterized sorghum proteins which may have protease activities that can be explored in future validation studies.

09 BIOMASS FUELS↗
DOE OSTI2020

Single-Cell Transcriptomic Analysis of Tumor-Derived Fibroblasts and Normal Tissue-Resident Fibroblasts Reveals Fibroblast Heterogeneity in Breast Cancer

Cancer-associated fibroblasts (CAFs) are a prominent stromal cell type in solid tumors and molecules secreted by CAFs play an important role in tumor progression and metastasis. CAFs coexist as heterogeneous populations with potentially different biological functions. Although CAFs are a major component of the breast cancer stroma, molecular and phenotypic heterogeneity of CAFs in breast cancer is poorly understood. In this study, we investigated CAF heterogeneity in triple-negative breast cancer (TNBC) using a syngeneic mouse model, BALB/c-derived 4T1 mammary tumors. Using single-cell RNA sequencing (scRNA-seq), we identified six CAF subpopulations in 4T1 tumors including: 1) myofibroblastic CAFs, enriched for α-smooth muscle actin and several other contractile proteins; 2) ‘inflammatory’ CAFs with elevated expression of inflammatory cytokines; and 3) a CAF subpopulation expressing major histocompatibility complex (MHC) class II proteins that are generally expressed in antigen-presenting cells. Comparison of 4T1-derived CAFs to CAFs from pancreatic cancer revealed that these three CAF subpopulations exist in both tumor types. Interestingly, cells with inflammatory and MHC class II-expressing CAF profiles were also detected in normal breast/pancreas tissue, suggesting that these phenotypes are not tumor microenvironment-induced. This work enhances our understanding of CAF heterogeneity, and specifically targeting these CAF subpopulations could be an effective therapeutic approach for treating highly aggressive TNBCs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Unique Structural Features Relate to Evolutionary Adaptation of Cytochrome P450 in the Abyssal Zone

Cytochromes P450 (CYPs) form one of the largest enzyme superfamilies, with similar structural folds yet biological functions varying from synthesis of physiologically essential compounds to metabolism of myriad xenobiotics. Sterol 14α-demethylases (CYP51s) represent a very special P450 family, regarded as a possible evolutionary progenitor for all currently existing P450s. In metazoans CYP51 is critical for the biosynthesis of sterols including cholesterol. Here we determined the crystal structures of ligand-free CYP51s from the abyssal fish Coryphaenoides armatus and human-. Comparative sequence–structure–function analysis revealed specific structural elements that imply elevated conformational flexibility, uncovering a molecular basis for faster catalytic rates, lower substrate selectivity, and intrinsic resistance to inhibition. In addition, the C. armatus structure displayed a large-scale repositioning of structural segments that, in vivo, are immersed in the endoplasmic reticulum membrane and border the substrate entrance (the FG arm, >20 Å, and the β4 hairpin, >15 Å). The structural distinction of C. armatus CYP51, which is the first structurally characterized deep sea P450, suggests stronger involvement of the membrane environment in regulation of the enzyme function. We interpret this as a co-adaptation of the membrane protein structure with membrane lipid composition during evolutionary incursion to life in the deep sea.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Adipose-specific knockout of Protein Kinase D1 suppresses de novo lipogenesis in mice via SREBP1c-dependent signaling

Highlights: • Protein Kinase D1 is a functional regulator in adipocyte and white adipose tissue. • PKD1 drives de novo lipogenesis in adipocyte. • SREBP1c signaling regulates lipogenesis in adipocyte. • PKD1-mediated de novo lipogenesis relies on SREBP1c. Having healthy adipose tissue is essential for metabolic health, as excessive adipose tissue in the body can cause its dysregulation and driving chronic metabolic diseases. Protein kinase D1 (PKD1) is considered to be a key kinase in signal transduction, which regulates multiple cellular functions, but its physiological functions in adipose are still not fully understood. This study aimed at elucidating the function of adipocyte PKD1 on lipogenesis. From RNA-Sequencing data, we found that the fatty acid biosynthesis pathway in white adipose tissue lacking PKD1 was significantly affected. Critical rate-limiting enzymes for de novo lipogenesis in adipocytes, such as FASN, ACCα, and SCD1, were significantly repressed after deleting PKD1 in vivo and in vitro. Further studies revealed that blockade of PKD1 significantly increased phosphorylation of SREBP1c at serine 372 site. Co-immunoprecipitation analysis showed that PKD1 interacts with SREBP1c in vitro and in vivo. Importantly, overexpression of SREBP1c reversed the inhibition of FASN and ACCα expression caused by PKD1 silencing. Together, adipocyte PKD1 promotes de novo lipogenesis via SREBP1c-dependent manner in visceral white adipose tissue and might provide a new target for the development of anti-obesity therapies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Changes in DNA Methylation in Response to 6-Benzylaminopurine Affect Allele-Specific Gene Expression in Populus Tomentosa

Cytokinins play important roles in the growth and development of plants. Physiological and photosynthetic characteristics are common indicators to measure the growth and development in plants. However, few reports have described the molecular mechanisms of physiological and photosynthetic changes in response to cytokinin, particularly in woody plants. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression in response to the external environment. In this study, we examined genome-wide DNA methylation variation and transcriptional variation in poplar (Populus tomentosa) after short-term treatment with the synthetic cytokinin 6-benzylaminopurine (6-BA). We identified 460 significantly differentially methylated regions (DMRs) in response to 6-BA treatment. Transcriptome analysis showed that 339 protein-coding genes, 262 long non-coding RNAs (lncRNAs), and 15,793 24-nt small interfering RNAs (siRNAs) were differentially expressed under 6-BA treatment. Among these, 79% were differentially expressed between alleles in P. tomentosa, and 102,819 allele-specific expression (ASE) loci in 19,200 genes were detected showing differences in ASE levels after 6-BA treatment. Combined DNA methylation and gene expression analysis demonstrated that DNA methylation plays an important role in regulating allele-specific gene expression. To further investigate the relationship between these 6-BA-responsive genes and phenotypic variation, we performed SNP analysis of 460 6-BA-responsive DMRs via re-sequencing using a natural population of P. tomentosa, and we identified 206 SNPs that were significantly associated with growth and wood properties. Association analysis indicated that 53% of loci with allele-specific expression had primarily dominant effects on poplar traits. Our comprehensive analyses of P. tomentosa DNA methylation and the regulation of allele-specific gene expression suggest that DNA methylation is an important regulator of imbalanced expression between allelic loci.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Integration of high-density genetic mapping with transcriptome analysis uncovers numerous agronomic QTL and reveals candidate genes for the control of tillering in sorghum

Phenotypes such as branching, photoperiod sensitivity, and height were modified during plant domestication and crop improvement. Here, we perform quantitative trait locus (QTL) mapping of these and other agronomic traits in a recombinant inbred line (RIL) population derived from an interspecific cross between Sorghum propinquum and Sorghum bicolor inbred Tx7000. Using low-coverage Illumina sequencing and a bin-mapping approach, we generated ~1920 bin markers spanning ~875 cM. Phenotyping data were collected and analyzed from two field locations and one greenhouse experiment for six agronomic traits, thereby identifying a total of 30 QTL. Many of these QTL were penetrant across environments and co-mapped with major QTL identified in other studies. Other QTL uncovered new genomic regions associated with these traits, and some of these were environment-specific in their action. To further dissect the genetic underpinnings of tillering, we complemented QTL analysis with transcriptomics, identifying 6189 genes that were differentially expressed during tiller bud elongation. We identified genes such as Dormancy Associated Protein 1 (DRM1) in addition to various transcription factors that are differentially expressed in comparisons of dormant to elongating tiller buds and lie within tillering QTL, suggesting that these genes are key regulators of tiller elongation in sorghum. Our study demonstrates the usefulness of this RIL population in detecting domestication and improvement-associated genes in sorghum, thus providing a valuable resource for genetic investigation and improvement to the sorghum community.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii : Exploring Constitutive and Plastic Responses Using Reciprocal Transplants

Abstract The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST–QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis–environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

3′ TagSeq↗
DOE OSTI2025

Different chemical scaffolds bind to L-phe site in Mycobacterium tuberculosis Phe-tRNA synthetase

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mt), is one of the deadliest infectious diseases. The rise of multidrug-resistant strains represents a major public health threat, requiring new therapeutic options. Bacterial aminoacyl-tRNA synthetases (aaRS) have been shown to be highly promising drug targets, including for TB treatment. These enzymes play an essential role in translating the DNA gene code into protein sequence by attaching specific amino acid to their cognate tRNAs. They have multiple binding sites that can be targeted for inhibitor discovery: amino acid binding pocket, ATP binding pocket, tRNA binding site and an editing domain. Recently we reported several high-resolution structures of M. tuberculosis phenylalanyl-tRNA synthetase (MtPheRS) complexed with tRNA Phe and either L-Phe or a nonhydrolyzable phenylalanine adenylate analog. Here, in this study, using Nucleic Magnetic Resonance (NMR) and Surface Plasmon Resonance (SPR) we identified fragments that bind to MtPheRS and we determined crystal structures of their complexes with MtPheRS/tRNA Phe . All the binders interact with the L-Phe amino acid binding site. The analysis of interactions of the new compounds combined with adenylate analog structure provides insights for the rational design of antituberculosis drugs. The 3 ' arm of the tRNA Phe in all the structures was disordered with exception of one complex with D-735 compound. In this structure the 3' CCA end of the acceptor stem is observed in the editing domain of MtPheRS providing insights regarding the post-transfer editing activity of class II aaRS.

Gade, Priyanka [Univ. of Chicago, IL (United State↗
DOE OSTI2025

Intracellular Biomacromolecule Delivery by Stimuli-Responsive Protein Vesicles Loaded by Hydrophobic Ion Pairing

Proteins can perform ideal therapeutic functions. However, their large size and significant surface hydrophilicity and charge prohibit them from reaching intracellular targets. These chemical features also render them poorly encapsulated by nanoparticles used for intracellular delivery. In this work, a novel combination of protein vesicles and hydrophobic ion pairing (HIP) was used to load protein cargo and achieve cytosolic delivery to overcome the limitations of previous protein vesicle properties. Protein vesicles are thermally self-assembling nanoparticles made from elastin-like polypeptide (ELP) fused to an arginine-rich leucine zipper and a globular protein fused to a glutamate-rich leucine zipper. To impart stimuli-responsive disassembly, physiological stability, and small size, the ELP sequence was modified to include histidine and tyrosine residues. HIP was used to load and release protein cargo requiring endosomal escape for cytosolic function. HIP vesicles enabled delivery of cytochrome c, a cytosolically active protein, and a significant reduction in viability in both a traditional two-dimensional (2D) human cancer cell line culture and a biomimetic three-dimensional (3D) organoid model of acute myeloid leukemia. By examining the uptake of positively and negatively charged fluorescent protein cargos loaded by HIP, this work revealed the necessity of HIP for cytosolic cargo delivery and how HIP loading influences protein vesicle self-assembly and disassembly using microscopy, small-angle X-ray scattering, and nanoparticle tracking analysis. HIP protein vesicles have the potential to broaden the use of intracellular proteins as therapeutics for various diseases and extend protein vesicles to deliver other biomacromolecules, as the strategy developed here resulted in the first cytosolic protein cargo delivery using protein vesicles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Data for Phylogenetic diversity of light-dependent phosphorylation of Thr78 in Rubisco activase

Rubisco activase is an ATP-dependent chaperone that facilitates dissociation of inhibitory sugar phosphates from the catalytic sites of Rubisco during photosynthesis. In Arabidopsis, Rubisco activase is negatively regulated by dark-dependent phosphorylation of Thr78. The prevalence of Thr78 in Rubisco activase was investigated across sequences from 91 plant species, finding that 29 (∼32%) species shared a threonine in the same position. Analysis of seven C3 species with an antibody raised against a Thr78 phospho-peptide demonstrated that this position is phosphorylated in multiple genera. However, light-dependent dephosphorylation of Thr78 was observed only in Arabidopsis. Further, phosphorylation of Thr78 could not be detected in any of the four C4 grass species examined. The results suggest that despite conservation of Thr78 in Rubisco activase from a wide range of species, a regulatory role for phosphorylation at this site is more limited. This provides a case study for how variation in post-translational regulation can amplify functional divergence across the phylogeny of plants beyond what is explained by sequence variation in a metabolically important protein.

photosynthesis↗
DOE OSTI2025

Phylogenetic diversity of light-dependent phosphorylation of Thr78 in Rubisco activase

Rubisco activase is an ATP-dependent chaperone that facilitates dissociation of inhibitory sugar phosphates from the catalytic sites of Rubisco during photosynthesis. In Arabidopsis, Rubisco activase is negatively regulated by dark-dependent phosphorylation of Thr78. The prevalence of Thr78 in Rubisco activase was investigated across sequences from 91 plant species, finding that 29 (∼32%) species shared a threonine in the same position. Analysis of seven C3 species with an antibody raised against a Thr78 phospho-peptide demonstrated that this position is phosphorylated in multiple genera. However, light-dependent dephosphorylation of Thr78 was observed only in Arabidopsis. Further, phosphorylation of Thr78 could not be detected in any of the four C4 grass species examined. The results suggest that despite conservation of Thr78 in Rubisco activase from a wide range of species, a regulatory role for phosphorylation at this site is more limited. Furthermore, this provides a case study for how variation in post-translational regulation can amplify functional divergence across the phylogeny of plants beyond what is explained by sequence variation in a metabolically important protein.

Arabidopsis↗
DOE OSTI2020

Single-Cell Genomics of Novel Actinobacteria With the Wood–Ljungdahl Pathway Discovered in a Serpentinizing System

Serpentinite-hosted systems represent modern-day analogs of early Earth environments. In these systems, water-rock interactions generate highly alkaline and reducing fluids that can contain hydrogen, methane, and low-molecular-weight hydrocarbons-potent reductants capable of fueling microbial metabolism. In this study, we investigated the microbiota of Hakuba Happo hot springs (~50°C; pH~10.5–11), located in Nagano (Japan), which are impacted by the serpentinization process. Analysis of the 16S rRNA gene amplicon sequences revealed that the bacterial community comprises Nitrospirae (47%), “Parcubacteria” (19%), Deinococcus-Thermus (16%), and Actinobacteria (9%), among others. Notably, only 57 amplicon sequence variants (ASV) were detected, and fifteen of these accounted for 90% of the amplicons. Among the abundant ASVs, an early-branching, uncultivated actinobacterial clade identified as RBG-16-55-12 in the SILVA database was detected. Ten single-cell genomes (average pairwise nucleotide identity: 0.98–1.00; estimated completeness: 33–93%; estimated genome size: ~2.3 Mb) that affiliated with this clade were obtained. Taxonomic classification using single copy genes indicates that the genomes belong to the actinobacterial class-level clade UBA1414 in the Genome Taxonomy Database. Based on metabolic pathway predictions, these actinobacteria are anaerobes, capable of glycolysis, dissimilatory nitrate reduction and CO 2 fixation via the Wood–Ljungdahl (WL) pathway. Several other genomes within UBA1414 and two related class-level clades also encode the WL pathway, which has not yet been reported for the Actinobacteria phylum. For the Hakuba actinobacterium, the energy metabolism related to the WL pathway is likely supported by a combination of the Rnf complex, group 3b and 3d [NiFe]-hydrogenases, [FeFe]-hydrogenases, and V-type (H + /Na + pump) ATPase. The genomes also harbor a form IV ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) complex, also known as a RubisCO-like protein, and contain signatures of interactions with viruses, including clustered regularly interspaced short palindromic repeat (CRISPR) regions and several phage integrases. This is the first report and detailed genome analysis of a bacterium within the Actinobacteria phylum capable of utilizing the WL pathway. The Hakuba actinobacterium is a member of the clade UBA1414/RBG-16-55-12, formerly within the group “OPB41.” We propose to name this bacterium ‘Candidatus Hakubanella thermoalkaliphilus.'

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Birth of protein folds and functions in the virome

The rapid evolution of viruses generates proteins that are essential for infectivity and replication but with unknown functions, due to extreme sequence divergence. Here, using a database of 67,715 newly predicted protein structures from 4,463 eukaryotic viral species, we found that 62% of viral proteins are structurally distinct and lack homologues in the AlphaFold database. Among the remaining 38% of viral proteins, many have non-viral structural analogues that revealed surprising similarities between human pathogens and their eukaryotic hosts. Structural comparisons suggested putative functions for up to 25% of unannotated viral proteins, including those with roles in the evasion of innate immunity. In particular, RNA ligase T-like phosphodiesterases were found to resemble phage-encoded proteins that hydrolyse the host immune-activating cyclic dinucleotides 3',3'- and 2',3'-cyclic GMP-AMP (cGAMP). Experimental analysis showed that RNA ligase T homologues encoded by avian poxviruses similarly hydrolyse cGAMP, showing that RNA ligase T-mediated targeting of cGAMP is an evolutionarily conserved mechanism of immune evasion that is present in both bacteriophage and eukaryotic viruses. Together, the viral protein structural database and analyses presented here afford new opportunities to identify mechanisms of virus–host interactions that are common across the virome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

MAL33 drives natural variation in maltose metabolism in Saccharomyces eubayanus

Maltose is one of the most abundant sugars in brewer’s wort, and its efficient utilization is critical for successful fermentation. However, maltose consumption varies naturally among Saccharomyces eubayanus strains isolated from different host trees, such as Quercus and Nothofagus. To identify the genetic determinants underlying these phenotypic differences, we performed bulk segregant analysis (BSA) and quantitative trait loci (QTL) mapping using an F 2 offspring derived from QC18 (Quercus-associated) and CL467.1 (Nothofagus-associated) strains. QTL mapping identified two significant genomic regions on subtelomeric loci of chromosomes V-R and XVI-L, each containing complete MAL loci composed of MAL32 (encoding maltase), MAL31 (transporter), and MAL33 (transcriptional activator) genes. Comparative polymorphism analyses identified mutations in MAL32 and MAL33 of QC18, including frameshift mutations resulting in premature stop codons. Functional validation demonstrated that the heterologous expression of MAL33 ChrV from CL467.1 fully restored maltose utilization in QC18, indicating the functional presence of MAL33 cis-regulatory sequences and MAL32 and MAL31 genes in QC18. While structural protein predictions identified truncation and impaired functionality in the maltose-responsive activation domain of Mal33p from QC18, overexpression of QC18’s own MAL33 ChrV allele also improved maltose metabolism, suggesting dosage-dependent transcriptional limitations rather than complete functional loss. These results indicate that allelic variations in the maltose-responsive activation domain of Mal33p result in differences in maltose consumption between strains. Here, we hypothesized that reduced maltose metabolism in QC18 is an adaptive response to the distinct sugar composition in Quercus robur bark, contrasting with the starch-rich environment of Nothofagus pumilio. These findings highlight subtelomeric MAL gene diversity as a reservoir of genetic variation, representing a key evolutionary mechanism that influences maltose adaptation among natural Saccharomyces isolates.

evolutionary plasticity↗
DOE OSTI2020

Fatty Acid and Alcohol Metabolism in Pseudomonas putida: Functional Analysis Using Random Barcode Transposon Sequencing.

With its ability to catabolize a wide variety of carbon sources and a growing engineering toolkit, Pseudomonas putida KT2440 is emerging as an important chassis organism for metabolic engineering. Despite advances in our understanding of the organism, many gaps remain in our knowledge of the genetic basis of its metabolic capabilities. The gaps are particularly noticeable in our understanding of both fatty acid and alcohol catabolism, where many paralogs putatively coding for similar enzymes coexist, making biochemical assignment via sequence homology difficult. To rapidly assign function to the enzymes responsible for these metabolisms, we leveraged random barcode transposon sequencing (RB-Tn-Seq). Global fitness analyses of transposon libraries grown on 13 fatty acids and 10 alcohols produced strong phenotypes for hundreds of genes. Fitness data from mutant pools grown on fatty acids of varying chain lengths indicated specific enzyme substrate preferences and enabled us to hypothesize that DUF1302/DUF1329 family proteins potentially function as esterases. From the data, we also postulate catabolic routes for the two biogasoline molecules isoprenol and isopentanol, which are catabolized via leucine metabolism after initial oxidation and activation with coenzyme A (CoA). Because fatty acids and alcohols may serve as both feedstocks and final products of metabolic-engineering efforts, the fitness data presented here will help guide future genomic modifications toward higher titers, rates, and yields. IMPORTANCE To engineer novel metabolic pathways into P. putida , a comprehensive understanding of the genetic basis of its versatile metabolism is essential. Here, we provide functional evidence for the putative roles of hundreds of genes involved in the fatty acid and alcohol metabolism of the bacterium. These data provide a framework facilitating precise genetic changes to prevent product degradation and to channel the flux of specific pathway intermediates as desired.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Quantifying the Long‐Range Coupling of Electronic Properties in Proteins with ab initio Molecular Dynamics**

Abstract The delicate interplay of covalent and non‐covalent interactions in proteins is inherently quantum mechanical and highly dynamic in nature. To directly interrogate the evolving nature of the electronic structure of proteins, we carry out 100‐ps‐scale ab initio molecular dynamics simulations of three representative small proteins with range‐separated hybrid density functional theory. We quantify the nature and length‐scale of the coupling of residue‐specific charge probability distributions in these proteins. While some nonpolar residues exhibit expectedly narrow charge distributions, most polar and charged residues exhibit broad, multimodal distributions. Even for nonpolar residues, we observe sequence‐specific deviations corresponding to charge accumulation or depletion that would be challenging to capture in a fixed charge force field. We quantify the effect of residue‐residue interactions on charge distributions first with linear cross‐correlations. We then show how additional insight can be gained from evaluating the mutual information of charge distributions. We show that a significant number of residues couple most strongly with residues that are distant in both sequence and space over a range of secondary structures including α‐helical, β‐sheet, disulfide bridging, and lasso motifs. The mutual information analysis is necessary to capture coupling between some polar and charged residues that would be otherwise missed.

Yang, Zhongyue↗
DOE OSTI2023

A cellular selection identifies elongated flavodoxins that support electron transfer to sulfite reductase

Flavodoxins (Flds) mediate the flux of electrons between oxidoreductases in diverse metabolic pathways. To investigate whether Flds can support electron transfer to a sulfite reductase (SIR) that evolved to couple with a ferredoxin, we evaluated the ability of Flds to transfer electrons from a ferredoxin-NADP reductase (FNR) to a ferredoxin-dependent SIR using growth complementation of an Escherichia coli strain with a sulfur metabolism defect. We show that Flds from cyanobacteria complement this growth defect when coexpressed with an FNR and an SIR that evolved to couple with a plant ferredoxin. When we evaluated the effect of peptide insertion on Fld-mediated electron transfer, we observed a sensitivity to insertions within regions predicted to be proximal to the cofactor and partner binding sites, while a high insertion tolerance was detected within loops distal from the cofactor and within regions of helices and sheets that are proximal to those loops. Bioinformatic analysis showed that natural Fld sequence variability predicts a large fraction of the motifs that tolerate insertion of the octapeptide SGRPGSLS. In conclusion, these results represent the first evidence that Flds can support electron transfer to assimilatory SIRs, and they suggest that the pattern of insertion tolerance is influenced by interactions with oxidoreductase partners.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Spatial turnover of soil viral populations and genotypes overlain by cohesive responses to moisture in grasslands

Viruses shape microbial communities, food web dynamics, and carbon and nutrient cycling in diverse ecosystems. However, little is known about the patterns and drivers of viral community composition, particularly in soil, precluding a predictive understanding of viral impacts on terrestrial habitats. To investigate soil viral community assembly processes, here we analyzed 43 soil viromes from a rainfall manipulation experiment in a Mediterranean grassland in California. We identified 5,315 viral populations (viral operational taxonomic units [vOTUs] with a representative sequence ≥10 kbp) and found that viral community composition exhibited a highly significant distance–decay relationship within the 200-m 2 field site. This pattern was recapitulated by the intrapopulation microheterogeneity trends of prevalent vOTUs (detected in ≥90% of the viromes), which tended to exhibit negative correlations between spatial distance and the genomic similarity of their predominant allelic variants. Although significant spatial structuring was also observed in the bacterial and archaeal communities, the signal was dampened relative to the viromes, suggesting differences in local assembly drivers for viruses and prokaryotes and/or differences in the temporal scales captured by viromes and total DNA. Despite the overwhelming spatial signal, evidence for environmental filtering was revealed in a protein-sharing network analysis, wherein a group of related vOTUs predicted to infect actinobacteria was shown to be significantly enriched in low-moisture samples distributed throughout the field. Overall, our results indicate a highly diverse, dynamic, active, and spatially structured soil virosphere capable of rapid responses to changing environmental conditions.

59 BASIC BIOLOGICAL SCIENCES↗