Search NASA⌕ Search

SEARCH · Search NASA

Results for “Protein Biochemistry”

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
DOE OSTI2025

NGPINT V3: a containerized orchestration Python software for discovery of next-generation protein–protein interactions

Abstract Summary Batch yeast two-hybrid (Y2H) assays, leveraged with next-generation sequencing, have afforded successful innovations for the analysis of protein–protein interactions. NGPINT is a Conda-based software designed to process the millions of raw sequencing reads resulting from Y2H–next-generation interaction screens. Over time, increasing compatibility and dependency issues have prevented clean NGPINT installation and operation. A system-wide update was essential to continue effective use with its companion software, Y2H-SCORES. We present NGPINT V3, a containerized implementation built with both Singularity and Docker, allowing accessibility across virtually any operating system and computing environment. Availability and implementation This update includes streamlined dependencies and container images hosted on Sylabs (https://cloud.sylabs.io/library/schuyler/ngpint/ngpint) and Dockerhub (https://hub.docker.com/r/schuylerds/ngpint), facilitating easier adoption and integration into high-throughput and cloud-computing workflows. Full instructions and software can be also found in the GitHub repository https://github.com/Wiselab2/NGPINT_V3 and Zenodo https://doi.org/10.5281/zenodo.15256036.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Binding Free Energy Analysis of Colicin D, E3 and E8 to Their Respective Cognate Immunity Proteins Using Computational Simulations

Colicins are antimicrobial proteins produced by bacteria for the purpose of destroying neighboring bacteria. Colicin activity is neutralized by a specific cognate immunity protein in order to protect the host. This study investigates the structural and binding mechanisms underlying the interaction of colicin-D, -E3 and -E8 to their respective immunity proteins (ImD, Im3 and Im8) using structure prediction, molecular dynamics (MD) simulations and MM-PBSA approach of free energy calculations. High-confidence colicin-immunity (Col-Im) complex structures predicted using AlphaFold2 were subjected to MD simulations of 150 ns with GROMACS and were analyzed for the binding free energy calculation using gmx_MMPBSA. Results showed that the complex of Col_E3-Im3 exhibited the most favorable binding free energy, driven by strong van der Waals and electrostatic interactions. Col_D-ImD and Col_E8-Im8 also showed the favorable binding. Electrostatics and hydrogen bonding emerged as a key factor driving binding and stability, while polar solvation acted as a destabilizing factor across all systems. These outcomes provide an understanding of the molecular mechanisms of Col-Im systems, with potential applications for developing natural antimicrobials for food safety.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Functional Relevance of CASP16 Nucleic Acid Predictions as Evaluated by Structure Providers

ABSTRACT Accurate biomolecular structure prediction enables the prediction of mutational effects, the speculation of function based on predicted structural homology, the analysis of ligand binding modes, experimental model building, and many other applications. Such algorithms to predict essential functional and structural features remain out of reach for biomolecular complexes containing nucleic acids. Here, we report a quantitative and qualitative evaluation of nucleic acid structures for the CASP16 blind prediction challenge by 12 of the experimental groups who provided nucleic acid targets. Blind predictions accurately model secondary structure and some aspects of tertiary structure, including reasonable global folds for some complex RNAs; however, predictions often lack accuracy in the regions of highest functional importance. All models have inaccuracies in non‐canonical regions where, for example, the nucleic‐acid backbone bends, deviating from an A‐form helix geometry, or a base forms a non‐standard hydrogen bond (not a Watson‐Crick base pair). These bends and non‐canonical interactions are integral to forming functionally important regions such as RNA enzymatic active sites. Additionally, the modeling of conserved and functional interfaces between nucleic acids and ligands, proteins, or other nucleic acids remains poor. For some targets, the experimental structures may not represent the only structure the biomolecular complex occupies in solution or in its functional life cycle, posing a future challenge for the community.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Complete genomes of Asgard archaea reveal diverse integrated and mobile genetic elements

Asgard archaea are of great interest as the progenitors of Eukaryotes, but little is known about the mobile genetic elements (MGEs) that may shape their ongoing evolution. Here, we describe MGEs that replicate in Atabeyarchaeia, a wetland Asgard archaea lineage represented by two complete genomes. We used soil depth–resolved population metagenomic data sets to track 18 MGEs for which genome structures were defined and precise chromosome integration sites could be identified for confident host linkage. Additionally, we identified a complete 20.67 kbp circular plasmid and two family-level groups of viruses linked to Atabeyarchaeia, via CRISPR spacer targeting. Closely related 40 kbp viruses possess a hypervariable genomic region encoding combinations of specific genes for small cysteine-rich proteins structurally similar to restriction-homing endonucleases. One 10.9 kbp integrative conjugative element (ICE) integrates genomically into theAtabeyarchaeum deiterrae-1chromosome and has a 2.5 kbp circularizable element integrated within it. The 10.9 kbp ICE encodes an expressed Type IIG restriction-modification system with a sequence specificity matching an active methylation motif identified by Pacific Biosciences (PacBio) high-accuracy long-read (HiFi) metagenomic sequencing. Restriction-modification of Atabeyarchaeia differs from that of another coexisting Asgard archaea, Freyarchaeia, which has few identified MGEs but possesses diverse defense mechanisms, including DISARM and Hachiman, not found in Atabeyarchaeia. Overall, defense systems and methylation mechanisms of Asgard archaea likely modulate their interactions with MGEs, and integration/excision and copy number variation of MGEs in turn enable host genetic versatility.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Characterization of Humanized Mouse Model of Organophosphate Poisoning and Detection of Countermeasures via MALDI-MSI

Organophosphoate (OP) chemicals are known to inhibit the enzyme acetylcholinesterase (AChE). Studying OP poisoning is difficult because common small animal research models have serum carboxylesterase, which contributes to animals’ resistance to OP poisoning. Historically, guinea pigs have been used for this research; however, a novel genetically modified mouse strain (KIKO) was developed with nonfunctional serum carboxylase (Es1 KO) and an altered acetylcholinesterase (AChE) gene, which expresses the amino acid sequence of the human form of the same protein (AChE KI). KIKO mice were injected with 1xLD50 of an OP nerve agent or vehicle control with or without atropine. After one to three minutes, animals were injected with 35 mg/kg of the currently fielded Reactivator countermeasure for OP poisoning. Postmortem brains were imaged on a Bruker RapifleX ToF/ToF instrument. Data confirmed the presence of increased acetylcholine in OP-exposed animals, regardless of treatment or atropine status. More interestingly, we detected a small amount of Reactivator within the brain of both exposed and unexposed animals; it is currently debated if reactivators can cross the blood–brain barrier. Further, we were able to simultaneously image acetylcholine, the primary affected neurotransmitter, as well as determine the location of both Reactivator and acetylcholine in the brain. This study, which utilized sensitive MALDI-MSI methods, characterized KIKO mice as a functional model for OP countermeasure development.

2-PAM↗
DOE OSTI2025

Matrix Metalloproteinases as Candidate Antigenic Determinants for Anti‐Tumor Autoantibodies in Human Ovarian Cancer: A Post Hoc Analysis

Circulating antibodies in patients with cancer can facilitate the identification of accessible epitopes on autoantigens expressed by tumors. To identify previously unrecognized protein targets in ovarian cancer, we computationally assessed a heptapeptide consensus motif (VPELGHE, flanked by two cysteine residues yielding a cyclic nonapeptide under oxidizing conditions) previously discovered via phage display-based epitope mapping of autoantibodies in patients. Eight proteins associated with ovarian cancer encompass amino acid sequences similar to the consensus motif and were, therefore, considered as candidate native autoantigens. Among these candidate targets, however, matrix metalloproteinase 14 (MMP14) demonstrates gene expression that is both high and negatively correlated with survival in ovarian cancer patient cohorts. MMP14 protein levels are also stable in tumor versus non-tumor tissues. Moreover, the corresponding heptapeptide mimic in MMP14 occurs within an α-helical secondary structural element observed in its catalytic domain. These findings demonstrate that a subset of patient-derived autoantibodies may interact with a previously unknown antigenic epitope found in MMP14 and other MMPs, thereby providing opportunities for the development of new targeted agents.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Integrative mapping reveals molecular features underlying the mechanism of nucleocytoplasmic transport

Nuclear pore complexes (NPCs) enable rapid, selective, and robust nucleocytoplasmic transport. To explain how transport emerges from the system components and their interactions, we used experimental data and theoretical information to construct an integrative Brownian dynamics model of transport through an NPC, coupled to a kinetic model of transport in the cell. The model recapitulates key aspects of transport for a wide range of molecular cargoes, including preribosomes and viral capsids. Our model quantifies how flexible phenylalanine-glycine (FG) repeat proteins create an entropic barrier to passive diffusion and how this barrier is selectively lowered in facilitated diffusion by the many transient interactions of nuclear transport receptors with the FG repeats. Selective transport is enhanced by “fuzzy” multivalent interactions, redundant FG repeat mass, coupling to the energy-dependent RanGTP concentration gradient, and exponential dependence of transport kinetics on the transport barrier. Our model will facilitate rational modulation of the NPC and its artificial mimics.

Raveh, Barak [School of Computer Science and Engin↗
DOE OSTI2024

Investigating the mechanism of chloroplast singlet oxygen signaling in the Arabidopsis thaliana accelerated cell death 2 mutant

As sessile organisms, plants have evolved complex signaling mechanisms to sense stress and acclimate. This includes the use of reactive oxygen species (ROS) generated during dysfunctional photosynthesis to initiate signaling. One such ROS, singlet oxygen ( 1 O 2 ), can trigger retrograde signaling, chloroplast degradation, and programmed cell death. However, the signaling mechanisms are largely unknown. Several proteins (e.g. PUB4, OXI1, EX1) are proposed to play signaling roles across three Arabidopsis thaliana mutants that conditionally accumulate chloroplast 1 O 2 (fluorescent in blue light (flu), chlorina 1 (ch1), and plastid ferrochelatase 2 (fc2)). We previously demonstrated that these mutants reveal at least two chloroplast 1 O 2 signaling pathways (represented by flu and fc2/ch1). Here, we test if the 1 O 2 -accumulating lesion mimic mutant, accelerated cell death 2 (acd2), also utilizes these pathways. The pub4–6 allele delayed lesion formation in acd2 and restored photosynthetic efficiency and biomass. Conversely, an oxi1 mutation had no measurable effect on these phenotypes. acd2 mutants were not sensitive to excess light (EL) stress, yet pub4–6 and oxi1 both conferred EL tolerance within the acd2 background, suggesting that EL-induced 1 O 2 signaling pathways are independent from spontaneous lesion formation. Thus, 1 O 2 signaling in acd2 may represent a third (partially overlapping) pathway to control cellular degradation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Electronic Interactions Between the Receptor-Binding Domain of Omicron Variants and Angiotensin-Converting Enzyme 2: A Novel Amino Acid–Amino Acid Bond Pair Concept

SARS-CoV-2 remains a severe threat to worldwide public health, particularly as the virus continues to evolve and diversify into variants of concern (VOCs). Among these VOCs, Omicron variants exhibit unique phenotypic traits, such as immune evasion, transmissibility, and severity, due to numerous spike protein mutations and the rapid subvariant evolution. These Omicron subvariants have more than 15 mutations in the receptor-binding domain (RBD), a region of the SARS-CoV-2 spike protein that is important for recognition and binding with the angiotensin-converting enzyme 2 (ACE2) human receptor. To address the impact of these high numbers of Omicron mutations on the binding process, we have developed a novel method to precisely quantify amino acid interactions via the amino acid–amino acid bond pair (AABP). We applied this concept to investigate the interface interactions of the RBD–ACE2 complex in four Omicron Variants (BA.1, BA.2, BA.5, and XBB.1.16) with its Wild Type counterpart. Based on the AABP analysis, we have identified all the sites that are affected by mutation and have provided evidence that unmutated sites are also impacted by mutation. We have calculated that the binding between RBD and ACE2 is strongest in OV BA.1, followed by OV BA.2, WT, OV BA.5, and OV XBB.1.16. We also present the partial charge values for all 311 residues across these five models. Our analysis provides a detailed understanding of changes caused by mutation in each Omicron interface complex.

Biochemistry & Molecular Biology↗
DOE OSTI2025

LIPID DROPLET PROTEIN OF SEEDS is involved in the control of lipid droplet size in Arabidopsis seeds and seedlings

Abstract In oilseeds, energy-rich carbon is stored as triacylglycerols in organelles called lipid droplets (LDs). While several of the major biogenetic proteins involved in LD formation have been identified, the full repertoire of LD proteins and their functional roles remains incomplete. Here, we show that the low-abundance, seed-specific LD protein LIPID DROPLET PROTEIN OF SEEDS (LDPS) contains an amphipathic α-helix and proline hairpin motif that serves as an LD-targeting signal and a separate region that binds to the LD protein OLEOSIN 1 (OLEO1). Loss of LDPS function results in smaller LDs and less seed oil in comparison with wild type, while overexpression of LDPS results in an increase in LD size and seed oil content. Loss of LDPS function also results in an inability of LDs to undergo fusion during postgerminative seedling growth. Analysis of oleo1 and ldps single- and double-mutant seeds and freeze–thaw treatment of seeds revealed that OLEO1 suppresses the ability of LDPS to promote larger LDs. Collectively, our results identify LDPS as an important player in LD biology that functions together with OLEO1 to determine LD size in Arabidopsis (Arabidopsis thaliana) seeds and seedlings through a process that involves LD–LD fusion.

Biochemistry & Molecular Biology↗
DOE OSTI2025

In silico λ-dynamics predicts protein binding specificities to modified RNAs

Abstract RNA modifications shape gene expression through a variety of chemical changes to canonical RNA bases. Although numbering in the hundreds, only a few RNA modifications are well characterized, in part due to the absence of methods to identify modification sites. Antibodies remain a common tool to identify modified RNA and infer modification sites through straightforward applications. However, specificity issues can result in off-target binding and confound conclusions. This work utilizes in silico λ-dynamics to efficiently estimate binding free energy differences of modification-targeting antibodies between a variety of naturally occurring RNA modifications. Crystal structures of inosine and N6-methyladenosine (m6A) targeting antibodies bound to their modified ribonucleosides were determined and served as structural starting points. λ-Dynamics was utilized to predict RNA modifications that permit or inhibit binding to these antibodies. In vitro RNA-antibody binding assays supported the accuracy of these in silico results. High agreement between experimental and computed binding propensities demonstrated that λ-dynamics can serve as a predictive screen for antibody specificity against libraries of RNA modifications. More importantly, this strategy is an innovative way to elucidate how hundreds of known RNA modifications interact with biological molecules without the limitations imposed by in vitro or in vivo methodologies.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Q -score as a reliability measure for protein, nucleic acid and small-molecule atomic coordinate models derived from 3DEM maps

Atomic coordinate models are important for the interpretation of 3D maps produced with cryoEM and cryoET (3D electron microscopy; 3DEM). In addition to visual inspection of such maps and models, quantitative metrics can inform about the reliability of the atomic coordinates, in particular how well the model is supported by the experimentally determined 3DEM map. A recently introduced metric, Q-score, was shown to correlate well with the reported resolution of the map for well fitted models. Here, we present new statistical analyses of Q-score based on its application to ∼10 000 maps and models archived in the EMDB (Electron Microscopy Data Bank) and PDB (Protein Data Bank). Further, we introduce two new metrics based on Q-score to represent each map and model relative to all entries in the EMDB and those with similar resolution. We explore through illustrative examples of proteins, nucleic acids and small molecules how Q-scores can indicate whether the atomic coordinates are well fitted to 3DEM maps and also whether some parts of a map may be poorly resolved due to factors such as molecular flexibility, radiation damage and/or conformational heterogeneity. These examples and statistical analyses provide a basis for how Q-scores can be interpreted effectively in order to evaluate 3DEM maps and atomic coordinate models prior to publication and archiving.

B factors↗
DOE OSTI2024

Current State, Challenges, and Opportunities in Genome-Scale Resource Allocation Models: A Mathematical Perspective

Stoichiometric genome-scale metabolic models (generally abbreviated GSM, GSMM, or GEM) have had many applications in exploring phenotypes and guiding metabolic engineering interventions. Nevertheless, these models and predictions thereof can become limited as they do not directly account for protein cost, enzyme kinetics, and cell surface or volume proteome limitations. Lack of such mechanistic detail could lead to overly optimistic predictions and engineered strains. Initial efforts to correct these deficiencies were by the application of precursor tools for GSMs, such as flux balance analysis with molecular crowding. In the past decade, several frameworks have been introduced to incorporate proteome-related limitations using a genome-scale stoichiometric model as the reconstruction basis, which herein are called resource allocation models (RAMs). This review provides a broad overview of representative or commonly used existing RAM frameworks. This review discusses increasingly complex models, beginning with stoichiometric models to precursor to RAM frameworks to existing RAM frameworks. RAM frameworks are broadly divided into two categories: coarse-grained and fine-grained, with different strengths and challenges. Discussion includes pinpointing their utility, data needs, highlighting framework strengths and limitations, and appropriateness to various research endeavors, largely through contrasting their mathematical frameworks. Finally, promising future applications of RAMs are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1994

Three-dimensional structure of Schistosoma japonicum glutathione S-transferase fused with a six-amino acid conserved neutralizing epitope of gp41 from HIV

The 3-dimensional crystal structure of glutathione S-transferase (GST) of Schistosoma japonicum (Sj) fused with a conserved neutralizing epitope on gp41 (glycoprotein, 41 kDa) of human immunodeficiency virus type 1 (HIV-1) (Muster T et al., 1993, J Virol 67:6642-6647) was determined at 2.5 A resolution. The structure of the 3-3 isozyme rat GST of the mu gene class (Ji X, Zhang P, Armstrong RN, Gilliland GL, 1992, Biochemistry 31:10169-10184) was used as a molecular replacement model. The structure consists of a 4-stranded beta-sheet and 3 alpha-helices in domain 1 and 5 alpha-helices in domain 2. The space group of the Sj GST crystal is P4(3)2(1)2, with unit cell dimensions of a = b = 94.7 A, and c = 58.1 A. The crystal has 1 GST monomer per asymmetric unit, and 2 monomers that form an active dimer are related by crystallographic 2-fold symmetry. In the binding site, the ordered structure of reduced glutathione is observed. The gp41 peptide (Glu-Leu-Asp-Lys-Trp-Ala) fused to the C-terminus of Sj GST forms a loop stabilized by symmetry-related GSTs. The Sj GST structure is compared with previously determined GST structures of mammalian gene classes mu, alpha, and pi. Conserved amino acid residues among the 4 GSTs that are important for hydrophobic and hydrophilic interactions for dimer association and glutathione binding are discussed.

Schistosoma japonicum/enzymology/genetics/immunolo↗
DOE OSTI2025

Extracellular filaments revealed by affinity capture cryogenic-electron tomography

Cryogenic-electron tomography (cryo-ET) has provided an unprecedented glimpse into the nanoscale architecture of cells by combining cryogenic preservation of biological structures with electron tomography. Micropatterning of extracellular matrix proteins is increasingly used as a method to prepare adherent cell types for cryo-ET as it promotes optimal positioning of cells and subcellular regions of interest for vitrification, cryo-focused ion beam (cryo-FIB) milling, and data acquisition. Here we demonstrate a micropatterning workflow for capturing minimally adherent cell types, human T cells and Jurkat cells, for cryo-FIB and cryo-ET. Our affinity capture system facilitated the nanoscale imaging of Jurkat cells, revealing extracellular filamentous structures. It improved workflow efficiency by consistently producing grids with a sufficient number of well-positioned cells for an entire cryo-FIB session. Affinity capture can be extended to facilitate high-resolution imaging of other adherent and non-adherent cell types with cryo-ET.

Biochemistry↗
DOE OSTI2025

Supplementary table, figures and DNA sequences of sorghum gene models SbiRTx430.01G455400 and SbiRTx.02G006600 that feature primers, gRNAs and indels created

In-context promoter bashing via genome editing is a route to identify and characterize critical regulatory regions that govern expression of genes of interest. The outcomes of in-context promoter bashing can be used to inform editing strategies to modulate the expression of selected gene models in a desired fashion. Here we employed in-context promoter bashing to characterize the proximal upstream regulatory regions of sorghum genes encoding phosphoenolpyruvate carboxykinase (Sb.PEPCK.BS, SbiTx430.01G455400) and alanine aminotransferase (SbiTx430.02G006600, SbAlaAT.BS), two proteins involved in the PCK C4 pathway. Characterized germinal edits within the targeted regions upstream of these two genes ranged in size from 138 bp up to 1790 bp. A 138 bp within the Sb.PEPCK.BS upstream region and a 1643 bp element within the Sb.AlaAT.BS upstream region were determined to be important for maintenance of transcription levels. No change in development or various physiological parameters was observed in characterized lineages carrying promoter edits. However, significant changes in seed reserves and a reduction in 100 seed weight were consistently observed, under both greenhouse and field environments, in plants carrying an edit in the promoter of Sb.PEPCK.BS gene were significantly reduced in transcript accumulation for this gene.

Quach, Truyen [Center for Plant Science Innovation↗
DOE OSTI2024

Structural and functional insights into the interaction between the bacteriophage T4 DNA processing proteins gp32 and Dda

Abstract Bacteriophage T4 is a classic model system for studying the mechanisms of DNA processing. A key protein in T4 DNA processing is the gp32 single-stranded DNA-binding protein. gp32 has two key functions: it binds cooperatively to single-stranded DNA (ssDNA) to protect it from nucleases and remove regions of secondary structure, and it recruits proteins to initiate DNA processes including replication and repair. Dda is a T4 helicase recruited by gp32, and we purified and crystallized a gp32–Dda–ssDNA complex. The low-resolution structure revealed how the C-terminus of gp32 engages Dda. Analytical ultracentrifugation analyses were consistent with the crystal structure. An optimal Dda binding peptide from the gp32 C-terminus was identified using surface plasmon resonance. The crystal structure of the Dda–peptide complex was consistent with the corresponding interaction in the gp32–Dda–ssDNA structure. A Dda-dependent DNA unwinding assay supported the structural conclusions and confirmed that the bound gp32 sequesters the ssDNA generated by Dda. The structure of the gp32–Dda–ssDNA complex, together with the known structure of the gp32 body, reveals the entire ssDNA binding surface of gp32. gp32–Dda–ssDNA complexes in the crystal are connected by the N-terminal region of one gp32 binding to an adjacent gp32, and this provides key insights into this interaction.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Untargeted GC-MS Metabolic Profiling of Anaerobic Gut Fungi Reveals Putative Terpenoids and Strain-Specific Metabolites

Background/Objectives: Anaerobic gut fungi (Neocallimastigomycota) are biotechnologically relevant, lignocellulose-degrading microbes with under-explored biosynthetic potential for secondary metabolites. Untargeted metabolomic profiling with gas chromatography–mass spectrometry (GC-MS) was applied to two gut fungal strains, Anaeromyces robustus and Caecomyces churrovis, to establish a foundational metabolomic dataset to identify metabolites and provide insights into gut fungal metabolic capabilities. Methods: Gut fungi were cultured anaerobically in rumen-fluid-based media with a soluble substrate (cellobiose), and metabolites were extracted using the Metabolite, Protein, and Lipid Extraction (MPLEx) method, enabling metabolomic and proteomic analysis from the same cell samples. Samples were derivatized and analyzed via GC-MS, followed by compound identification by spectral matching to reference databases, molecular networking, and statistical analyses. Results: Distinct metabolites were identified between A. robustus and C. churrovis, including 2,3-dihydroxyisovaleric acid produced by A. robustus and maltotriitol, maltotriose, and melibiose produced by C. churrovis. C. churrovis may polymerize maltotriose to form an extracellular polysaccharide, like pullulan. GC-MS profiling potentially captured sufficiently volatile products of proteomically detected, putative non-ribosomal peptide synthetases and polyketide synthases of A. robustus and C. churrovis. The triterpene squalene and triterpenoid tetrahymanol were putatively identified in A. robustus and C. churrovis. Their conserved, predicted biosynthetic genes—squalene synthase and squalene tetrahymanol cyclase—were identified in A. robustus, C. churrovis, and other anaerobic gut fungal genera. Conclusions: This study provides a foundational, untargeted metabolomic dataset to unmask gut fungal metabolic pathways and biosynthetic potential and to prioritize future efforts for compound isolation and identification.

Biochemistry & Molecular Biology↗