Search NASA⌕ Search

SEARCH · Search NASA

Results for “Network motifs”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

The Conservation of Structure and Mechanism of Catalytic Action in a Family of Thiamin Pyrophosphate (TPP)-dependent Enzymes

Thiamin pyrophosphate (TPP)-dependent enzymes are a divergent family of TPP and metal ion binding proteins that perform a wide range of functions with the common decarboxylation steps of a -(O=)C-C(OH)- fragment of alpha-ketoacids and alpha- hydroxyaldehydes. To determine how structure and catalytic action are conserved in the context of large sequence differences existing within this family of enzymes, we have carried out an analysis of TPP-dependent enzymes of known structures. The common structure of TPP-dependent enzymes is formed at the interface of four alpha/beta domains from at least two subunits, which provide for two metal and TPP-binding sites. Residues around these catalytic sites are conserved for functional purpose, while those further away from TPP are conserved for structural reasons. Together they provide a network of contacts required for flip-flop catalytic action within TPP-dependent enzymes. Thus our analysis defines a TPP-action motif that is proposed for annotating TPP-dependent enzymes for advancing functional proteomics.

Dominiak, P.↗
DOE OSTI2024

HSV ‐1 ICP0 dimer domain adopts a novel β‐barrel fold

Abstract Infected cell protein 0 (ICP0) is an immediate‐early regulatory protein of herpes simplex virus 1 (HSV‐1) that possesses E3 ubiquitin ligase activity. ICP0 transactivates viral genes, in part, through its C‐terminal dimer domain (residues 555–767). Deletion of this dimer domain results in reduced viral gene expression, lytic infection, and reactivation from latency. Since ICP0's dimer domain is associated with its transactivation activity and efficient viral replication, we wanted to determine the structure of this specific domain. The C‐terminus of ICP0 was purified from bacteria and analyzed by X‐ray crystallography to solve its structure. Each subunit or monomer in the ICP0 dimer is composed of nine β‐strands and two α‐helices. Interestingly, two adjacent β‐strands from one monomer “reach” into the adjacent subunit during dimer formation, generating two β‐barrel‐like structures. Additionally, crystallographic analyses indicate a tetramer structure is formed from two β‐strands of each dimer, creating a “stacking” of the β‐barrels. The structural protein database searches indicate the fold or structure adopted by the ICP0 dimer is novel. The dimer is held together by an extensive network of hydrogen bonds. Computational analyses reveal that ICP0 can either form a dimer or bind to SUMO1 via its C‐terminal SUMO‐interacting motifs but not both. Understanding the structure of the dimer domain will provide insights into the activities of ICP0 and, ultimately, the HSV‐1 life cycle.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Lipid sponge droplets as programmable synthetic organelles

Living cells segregate molecules and reactions in various subcellular compartments known as organelles. Spatial organization is likely essential for expanding the biochemical functions of synthetic reaction systems, including artificial cells. Many studies have attempted to mimic organelle functions using lamellar membrane-bound vesicles. However, vesicles typically suffer from highly limited transport across the membranes and an inability to mimic the dense membrane networks typically found in organelles such as the endoplasmic reticulum. Here, we describe programmable synthetic organelles based on highly stable nonlamellar sponge phase droplets that spontaneously assemble from a single-chain galactolipid and nonionic detergents. Due to their nanoporous structure, lipid sponge droplets readily exchange materials with the surrounding environment. In addition, the sponge phase contains a dense network of lipid bilayers and nanometric aqueous channels, which allows different classes of molecules to partition based on their size, polarity, and specific binding motifs. The sequestration of biologically relevant macromolecules can be programmed by the addition of suitably functionalized amphiphiles to the droplets. We demonstrate that droplets can harbor functional soluble and transmembrane proteins, allowing for the colocalization and concentration of enzymes and substrates to enhance reaction rates. Droplets protect bound proteins from proteases, and these interactions can be engineered to be reversible and optically controlled. Our results show that lipid sponge droplets permit the facile integration of membrane-rich environments and self-assembling spatial organization with biochemical reaction systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Bonding similarities and differences between Y–Sb–Te and Sc–Sb–Te phase-change memory materials

The scandium (Sc) – alloyed Sb 2 Te 3 phase-change alloy has recently been found to enable ultrafast crystal nucleation due to the formation of Sc-stabilized octahedral motifs in the amorphous phase, rendering cache-type phase-change memory feasible. When yttrium (Y) is added, however, non-octahedral bonding patterns form in the amorphous Sb 2 Te 3 -based network even though Y has a valence electron configuration similar to that of Sc and also forms perfect octahedral bonding environments with tellurium in the YTe crystal. Here we elucidate the origin of this difference between Sc and Y, by carrying out thorough ab initio simulations and orbital-based bonding analyses on amorphous Y–Sb–Te and Sc–Sb–Te compounds. We also demonstrate how the smooth overlap of atomic positions (SOAP) similarity kernel can be used to quantify the structural similarity of local motifs in the amorphous phase with respect to various crystalline yttrium and scandium tellurides, both in the nearest-neighbor shell and beyond. We discover that the bonding contrast of Y- and Sc-centered structural motifs in amorphous Sb 2 Te 3 stems from their parent crystals at high Te concentrations. The larger atomic radius of Y and the weaker charge transfer when bonded with Te is found to allow more Te neighbors and cause a more open bonding environment, leading to higher coordination numbers and non-octahedral environments. We discuss the implications of the different local environments for practical applications in memory devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nickel(ii) complexes based on dithiolate–polyamine binary ligand systems: crystal structures, hirshfeld surface analysis, theoretical study, and catalytic activity study on photocatalytic hydrogen generation

To ascertain the influence of binary ligand systems [1,1-dicyanoethylene-2,2-dithiolate (i-mnt -2 ) and polyamine {tetraen = tris(2-aminoethyl)amine, tren = diethylene triamine and opda = o-phenylenediamine}] on the coordination modes of the Ni(II) metal center and resulting supramolecular architectures, a series of nickel(II) thiolate complexes [Ni(tetraen)(i-mnt)](DMSO) (1), [Ni 2 (tren) 2 (i-mnt) 2 ] (2), and [Ni 2 (i-mnt) 2 (opda) 2 ] n (3) have been synthesized in high yield in one step in water and structurally characterized by single crystal X-ray crystallography and spectroscopic techniques. X-ray diffraction studies disclose the diverse i-mnt -2 coordination to the Ni +2 center in the presence of active polyamine ligands, forming a slightly distorted octahedral geometry (NiN 4 S 2 ) in 1, square planar (NiS 4 ) and distorted octahedral geometries (NiN 6 ) in the bimetallic co-crystallized aggregate of cationic [Ni(tren) 2 ] +2 and anionic [Ni(i-mnt) 2 ] -2 in 2, and a one dimensional (1D) polymeric chain along the [100] axis in 3, having consecutive square planar (NiS 4 ) and octahedral (NiN 6 ) coordination kernels. The N–H···O, N–H···S, N–H···N, N–H···S, N–H···N, and N–H···O type hydrogen bonds stabilize the supramolecular assemblies in 1, 2, and 3 respectively imparting interesting graph-set-motifs. The molecular Hirshfeld surface analyses (HS) and 2D fingerprint plots were utilized for decoding all types of non-covalent contacts in the crystal networks. Atomic HS analysis of the Ni +2 centers reveals significant Ni–N metal–ligand interactions compared to Ni–S interactions. We have also studied the unorthodox interactions observed in the solid state structures of 1–3 by QTAIM and NBO analyses. Moreover, all the complexes proved to be highly active water reduction co-catalysts (WRC) in a photo-catalytic hydrogen evolution process involving iridium photosensitizers, wherein 2 and 3 having a square planar arrangement around the nickel center(s) – were found to be the most active ones, achieving 1000 and 1119 turnover numbers (TON), respectively.

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

BCL10 Mutations Define Distinct Dependencies Guiding Precision Therapy for DLBCL

Activated B cell–like diffuse large B-cell lymphomas (ABC-DLBCL) have unfavorable outcomes and chronic activation of CARD11–BCL10–MALT1 (CBM) signal amplification complexes that form due to polymerization of BCL10 subunits, which is affected by recurrent somatic mutations in ABC-DLBCLs. Herein, we show that BCL10 mutants fall into at least two functionally distinct classes: missense mutations of the BCL10 CARD domain and truncation of its C-terminal tail. Truncating mutations abrogated a motif through which MALT1 inhibits BCL10 polymerization, trapping MALT1 in its activated filament-bound state. CARD missense mutations enhanced BCL10 filament formation, forming glutamine network structures that stabilize BCL10 filaments. Mutant forms of BCL10 were less dependent on upstream CARD11 activation and thus manifested resistance to BTK inhibitors, whereas BCL10 truncating but not CARD mutants were hypersensitive to MALT1 inhibitors. Therefore, BCL10 mutations are potential biomarkers for BTK inhibitor resistance in ABC-DLBCL, and further precision can be achieved by selecting therapy based on specific biochemical effects of distinct mutation classes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

De novo design of protein structure and function with RFdiffusion

Abstract There has been considerable recent progress in designing new proteins using deep-learning methods 1–9 . Despite this progress, a general deep-learning framework for protein design that enables solution of a wide range of design challenges, including de novo binder design and design of higher-order symmetric architectures, has yet to be described. Diffusion models 10,11 have had considerable success in image and language generative modelling but limited success when applied to protein modelling, probably due to the complexity of protein backbone geometry and sequence–structure relationships. Here we show that by fine-tuning the RoseTTAFold structure prediction network on protein structure denoising tasks, we obtain a generative model of protein backbones that achieves outstanding performance on unconditional and topology-constrained protein monomer design, protein binder design, symmetric oligomer design, enzyme active site scaffolding and symmetric motif scaffolding for therapeutic and metal-binding protein design. We demonstrate the power and generality of the method, called RoseTTAFold diffusion (RFdiffusion), by experimentally characterizing the structures and functions of hundreds of designed symmetric assemblies, metal-binding proteins and protein binders. The accuracy of RFdiffusion is confirmed by the cryogenic electron microscopy structure of a designed binder in complex with influenza haemagglutinin that is nearly identical to the design model. In a manner analogous to networks that produce images from user-specified inputs, RFdiffusion enables the design of diverse functional proteins from simple molecular specifications.

Science & Technology - Other Topics↗
DOE OSTI2020

Reinterpreting π-stacking

The nature of π–π interactions has long been debated. The term “π-stacking” is considered by some to be a misnomer, in part because overlapping π-electron densities are thought to incur steric repulsion, and the physical origins of the widely-encountered “slip-stacked” motif have variously been attributed to either sterics or electrostatics, in competition with dispersion. Here, we use quantum-mechanical energy decomposition analysis to investigate π–π interactions in supramolecular complexes of polycyclic aromatic hydrocarbons, ranging in size up to realistic models of graphene, and for comparison we perform the same analysis on stacked complexes of polycyclic saturated hydrocarbons, which are cyclohexane-based analogues of graphane. Our results help to explain the short-range structure of liquid hydrocarbons that is inferred from neutron scattering, trends in melting-point data, the interlayer separation of graphene sheets, and finally band gaps and observation of molecular plasmons in graphene nanoribbons. Analysis of intermolecular forces demonstrates that aromatic π–π interactions constitute a unique and fundamentally quantum-mechanical form of non-bonded interaction. Not only do stacked π–π architectures enhance dispersion, but quadrupolar electrostatic interactions that may be repulsive at long range are rendered attractive at the intermolecular distances that characterize π-stacking, as a result of charge penetration effects. The planar geometries of aromatic sp 2 carbon networks lead to attractive interactions that are “served up on a molecular pizza peel”, and adoption of slip-stacked geometries minimizes steric (rather than electrostatic) repulsion. The slip-stacked motif therefore emerges not as a defect induced by electrostatic repulsion but rather as a natural outcome of a conformational landscape that is dominated by van der Waals interactions (dispersion plus Pauli repulsion), and is therefore fundamentally quantum-mechanical in its origins. This reinterpretation of the forces responsible for π-stacking has important implications for the manner in which non-bonded interactions are modeled using classical force fields, and for rationalizing the prevalence of the slip-stacked π–π motif in protein crystal structures.

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

Direct observation of 3D atomic packing in monatomic amorphous materials

Liquids and solids are two fundamental states of matter. Although the structure of crystalline solids has long been solved by crystallography, our understanding of the 3D atomic structure of liquids and amorphous materials remained speculative. In 1952, Frank hypothesized that icosahedral order is the prevalent atomic motif in monatomic liquids. Over the past six decades, there have been a great deal of experimental, computational and theoretical studies to understand the structure of liquids and amorphous materials. A polytetrahedral packing model was proposed to explain the 3D atomic structure of monatomic liquids and amorphous materials, in which icosahedral order is a key feature. The icosahedral order has also been found to play an important role in the structure of metallic glasses and quasicrystals. Here we advance atomic electron tomography to determine for the first time the 3D atomic positions in monatomic amorphous materials, including a Ta thin film and two Pd nanoparticles. Despite the different chemical composition and synthesis methods, we observe that pentagonal bipyramids are the most abundant atomic motifs in these amorphous materials. Contrary to the traditional understanding, the pentagonal bipyramids do not assemble full icosahedra, but are closely connected with each other to form networks that extend to medium-range scale. Our molecular dynamic simulations, coupled with experimental results, further reveal that the 3D atomic structure of monatomic liquids is similar to the experimental amorphous material. Moreover, we find that the pentagonal bipyramid networks, prevalent in monatomic liquids, rapidly grow in size and form a fraction of full icosahedra during the quench from a liquid to a metallic glass state, providing a possible new picture of structural evolution through the glass transition.

Yuan, Yakun↗
DOE OSTI2021

Size-Dependent Onset of Nitric Acid Dissociation in Cs + ·(HNO 3 )(H 2 O) n =0–11 Clusters at 20 K

We report the water-mediated charge separation of nitric acid upon incorporation into size-selected Cs + ∙(HNO 3 )(H 2 O) n=0-11 clusters at 20 K. Dramatic spectral changes are observed in the range n=7-9 that are traced to the formation of many isomeric structures associated with intermediate transfer of the acidic proton to the water network. This transfer is complete by n=10, which exhibits much simpler vibrational band patterns consistent with those expected for a tri-coordinated hydronium ion (the Eigen motif) along with the NO stretching bands predicted for a hydrated NO 3 – anion that is directly complexed to the Cs + cation. Theoretical analysis of the n=10 spectrum indicates that the dissociated ions adopt a solvent-separated ion-pair configuration such that the Cs + and H 3 O + cations flank the NO 3 – anion in a microhydrated salt bridge. In conclusion, this charge separation motif is evidently assisted by the electrostatic stabilization of the product NO 3 – /H 3 O + ion pair by the proximal metal ion.

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

Constructing Regulatory Networks to Compare Axenic and Interspecies Microbial Gene Transcription

In this preliminary study, we constructed gene regulatory networks (GRNs) from transcriptional expression data of axenic and interspecies microbial cultures with the goal of predicting how cocultivation affected greenhouse gas respiration by these species. The specific strains of Methylotuvimicrobium alkaliphilum 20Z, a methylotroph, and Cyanobacterium stanieri HL-69, a phototroph, were chosen for their viability in industrial bioprocessing. We ranked directed interactions between gene pairs based on the ability of the input gene to predict the expression of a target gene relative to their transcriptomes. While we were able to identify topological differences between conditions, our initial findings require validation through experimental analysis and further modeling. We aimed to develop a systematic thresholding approach to optimize the accuracy of our networks. We filtered out trial networks separately from top gene interactions of the scored rankings. Parameters of unfiltered and filtered networks were used to test and develop thresholding approaches. Knee point detection of edge weight distributions was explored as an approach for separating significant interactions from insignificant interactions in unfiltered networks. While knee detection failed to produce analogous networks for broad cross-condition comparisons, the results informed us about the proportions of significant edges present in unfiltered networks. We also calculated the average mean degree for nodes in a selection of trial networks to find a thresholding value characteristic to all groups. While we did not reach a definitive conclusion, we gained insight into the coregulatory structures of our groups and made critical evaluations of systematic methods for filtering networks. We recommend an iterative process for the inference of GRNs, where the most significant results from preliminary explorations are used to improve the efficiency with which regulatory motifs are chosen for experimental characterization. Experimental results can then inform the framework of adjusted models to improve broad interpretations of GRNs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Rapid Cis–Trans Coevolution Driven by a Novel Gene Retroposed from a Eukaryotic Conserved CCR4–NOT Component in Drosophila

Young, or newly evolved, genes arise ubiquitously across the tree of life, and they can rapidly acquire novel functions that influence a diverse array of biological processes. Previous work identified a young regulatory duplicate gene in Drosophila, Zeus that unexpectedly diverged rapidly from its parent, Caf40, an extremely conserved component in the CCR4–NOT machinery in post-transcriptional and post-translational regulation of eukaryotic cells, and took on roles in the male reproductive system. This neofunctionalization was accompanied by differential binding of the Zeus protein to loci throughout the Drosophila melanogaster genome. However, the way in which new DNA-binding proteins acquire and coevolve with their targets in the genome is not understood. Here, by comparing Zeus ChIP-Seq data from D. melanogaster and D. simulans to the ancestral Caf40 binding events from D. yakuba, a species that diverged before the duplication event, we found a dynamic pattern in which Zeus binding rapidly coevolved with a previously unknown DNA motif, which we term Caf40 and Zeus-Associated Motif (CAZAM), under the influence of positive selection. Interestingly, while both copies of Zeus acquired targets at male-biased and testis-specific genes, D. melanogaster and D. simulans proteins have specialized binding on different chromosomes, a pattern echoed in the evolution of the associated motif. Using CRISPR-Cas9-mediated gene knockout of Zeus and RNA-Seq, we found that Zeus regulated the expression of 661 differentially expressed genes (DEGs). Our results suggest that the evolution of young regulatory genes can be coupled to substantial rewiring of the transcriptional networks into which they integrate, even over short evolutionary timescales. Our results thus uncover dynamic genome-wide evolutionary processes associated with new genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Expanding the Scope of Bacterial CRISPR Activation with PAM-Flexible dCas9 Variants

CRISPR-Cas transcriptional tools have been widely applied for programmable regulation of complex biological networks. In comparison to eukaryotic systems, bacterial CRISPR activation (CRISPRa) has stringent target site requirements for effective gene activation. While genes may not always have an NGG protospacer adjacent motif (PAM) at the appropriate position, PAM-flexible dCas9 variants can expand the range of targetable sites. Here we systematically evaluate a panel of PAM-flexible dCas9 variants for their ability to activate bacterial genes. We observe that dxCas9-NG provides a high dynamic range of gene activation for sites with NGN PAMs while dSpRY permits modest activity across almost any PAM. Similar trends were observed for heterologous and endogenous promoters. For all variants tested, improved PAM-flexibility comes with the trade-off that CRISPRi-mediated gene repression becomes less effective. Weaker CRISPR interference (CRISPRi) gene repression can be partially rescued by expressing multiple sgRNAs to target many sites in the gene of interest. Furthermore, our work provides a framework to choose the most effective dCas9 variant for a given set of gene targets, which will further expand the utility of CRISPRa/i gene regulation in bacterial systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The SEEL motif and members of the MYB-related REVEILLE transcription factor family are important for the expression of LORELEI in the synergid cells of the Arabidopsis female gametophyte

Synergid cells in the micropylar end of the female gametophyte are required for critical cell–cell signaling interactions between the pollen tube and the ovule that precede double fertilization and seed formation in flowering plants. LORELEI (LRE) encodes a putative GPI-anchored protein that is expressed primarily in the synergid cells, and together with FERONIA, a receptor-like kinase, it controls pollen tube reception by the receptive synergid cell. Still, how LRE expression is controlled in synergid cells remains poorly characterized. We identified candidate cis-regulatory elements enriched in LRE and other synergid cell-expressed genes. One of the candidate motifs (‘TAATATCT’) in the LRE promoter was an uncharacterized variant of the Evening Element motif that we named as the Short Evening Element-like (SEEL) motif. Deletion or point mutations in the SEEL motif of the LRE promoter resulted in decreased reporter expression in synergid cells, demonstrating that the SEEL motif is important for expression of LRE in synergid cells. Additionally, we found that LRE expression is decreased in the loss of function mutants of REVEILLE (RVE) transcription factors, which are clock genes known to bind the SEEL and other closely related motifs. Here, we propose that RVE transcription factors regulate LRE expression in synergid cells by binding to the SEEL motif in the LRE promoter. Identification of cis-regulatory elements and transcription factors involved in the expression of LRE will serve as a foundation to characterize the gene regulatory networks in synergid cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

AI-Powered Knowledge Graphs for Neuromorphic and Energy-Efficient Computing

The surge in scientific literature obscures breakthroughs and hinders the discovery of new research paths. We propose an artificial intelligence (AI) powered framework using large language models (LLMs) and knowledge graphs (KGs) to automate parts of scientific discovery, focusing on energy-efficient AI circuits. Our hybrid approach combines LLMs, structured data, and ontology-based reasoning to construct a comprehensive knowledge graph that integrates insights across computational neuroscience, spiking neuron models, learning rules, architectural motifs, and neuromorphic device technologies. This multi-domain representation enables the generation of hypotheses that connect biological function with implementable, energy-efficient hardware architectures. Using KG embeddings and graph neural networks, the framework generates hypotheses for novel circuits, validates them through optimization on exascale HPC systems, and with tools like SuperNeuro and Fugu, the most promising designs will be prototyped in hardware. This open-source system aims to accelerate discoveries and bridging neuroscience with hardware innovation, drive collaboration, and unlock new opportunities in low-power AI computing.

Gautam, Ashish [ORNL]↗
DOE OSTI2022

METABOLIC: high-throughput profiling of microbial genomes for functional traits, metabolism, biogeochemistry, and community-scale functional networks

Background Advances in microbiome science are being driven in large part due to our ability to study and infer microbial ecology from genomes reconstructed from mixed microbial communities using metagenomics and single-cell genomics. Such omics-based techniques allow us to read genomic blueprints of microorganisms, decipher their functional capacities and activities, and reconstruct their roles in biogeochemical processes. Currently available tools for analyses of genomic data can annotate and depict metabolic functions to some extent; however, no standardized approaches are currently available for the comprehensive characterization of metabolic predictions, metabolite exchanges, microbial interactions, and microbial contributions to biogeochemical cycling. Results We present METABOLIC (METabolic And BiogeOchemistry anaLyses In miCrobes), a scalable software to advance microbial ecology and biogeochemistry studies using genomes at the resolution of individual organisms and/or microbial communities. The genome-scale workflow includes annotation of microbial genomes, motif validation of biochemically validated conserved protein residues, metabolic pathway analyses, and calculation of contributions to individual biogeochemical transformations and cycles. The community-scale workflow supplements genome-scale analyses with determination of genome abundance in the microbiome, potential microbial metabolic handoffs and metabolite exchange, reconstruction of functional networks, and determination of microbial contributions to biogeochemical cycles. METABOLIC can take input genomes from isolates, metagenome-assembled genomes, or single-cell genomes. Results are presented in the form of tables for metabolism and a variety of visualizations including biogeochemical cycling potential, representation of sequential metabolic transformations, community-scale microbial functional networks using a newly defined metric “MW-score” (metabolic weight score), and metabolic Sankey diagrams. METABOLIC takes ~ 3 h with 40 CPU threads to process ~ 100 genomes and corresponding metagenomic reads within which the most compute-demanding part of hmmsearch takes ~ 45 min, while it takes ~ 5 h to complete hmmsearch for ~ 3600 genomes. Tests of accuracy, robustness, and consistency suggest METABOLIC provides better performance compared to other software and online servers. To highlight the utility and versatility of METABOLIC, we demonstrate its capabilities on diverse metagenomic datasets from the marine subsurface, terrestrial subsurface, meadow soil, deep sea, freshwater lakes, wastewater, and the human gut. Conclusion METABOLIC enables the consistent and reproducible study of microbial community ecology and biogeochemistry using a foundation of genome-informed microbial metabolism, and will advance the integration of uncultivated organisms into metabolic and biogeochemical models. METABOLIC is written in Perl and R and is freely available under GPLv3 at https://github.com/AnantharamanLab/METABOLIC.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Hypernetwork science via high-order hypergraph walks

We propose high-order hypergraph walks as a framework to generalize graph-based network science techniques to hypergraphs. Edge incidence in hypergraphs is quantitative, yielding hypergraph walks with both length and width. Graph methods which then generalize to hypergraphs include connected component analyses, graph distance-based metrics such as closeness centrality, and motif-based measures such as clustering coefficients. We apply high-order analogs of these methods to real world hypernetworks, and show they reveal nuanced and interpretable structure that cannot be detected by graph-based methods. Lastly, we apply three generative models to the data and find that basic hypergraph properties, such as density and degree distributions, do not necessarily control these new structural measurements. Our work demonstrates how analyses of hypergraph-structured data are richer when utilizing tools tailored to capture hypergraph-native phenomena, and suggests one possible avenue towards that end.

97 MATHEMATICS AND COMPUTING↗