SEARCH · Engineering Papers

Results for “biological system modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2020

One Electron Multiple Proton Transfer in Model Organic Donor–Acceptor Systems: Implications for High-Frequency EPR

EPR spectroscopy is an important spectroscopic method for identification and characterization of radical species involved in many biological reactions. The tyrosyl radical is one of the most studied amino acid radical intermediates in biology. Often in conjunction with histidine residues, it is involved in many fundamental biological electron and proton transfer processes, such as in the water oxidation in photosystem II. As biological processes are typically extremely complicated and hard to control, molecular bio-mimetic model complexes are often used to clarify the mechanisms of the biological reactions. In this work, we present theoretical calculations to investigate the sensitivity of magnetic resonance parameters to proton-coupled electron transfer events, as well as conformational substates of the molecular constructs which mimic the tyrosine-histidine (Tyr-His) pairs found in a large variety of proteins. Upon oxidation of the phenol, the Tyr analog, these complexes can perform not only one-electron one-proton transfer (EPT), but also one-electron two-proton transfers (E2PT). It is shown that in aprotic environment the g X -components of the electronicg-tensor are extremely sensitive to the first proton transfer from the phenoxyl oxygen to the imidazole nitrogen (EPT product), leading to a significant increase of the g X -value of up to 0.003, but are not sensitive to the second proton transfer (E2PT). In the latter case, the change of the g X -value is much smaller (ca. 0.0001), which is too small to be distinguished even by high-frequency EPR. The 14 N hyperfine values are also too similar to allow differentiation between the different protonation states in EPT and E2PT. The magnetic resonance parameters were also calculated as a function of the rotation angles around single bonds. It was demonstrated that rotation of the phenoxyl group results in large positive changes (> 0.001) in the g X -values. Analysis of the data reveals that the main source of these changes is related to the strength of the H-bond between phenoxyl oxygen and the proton(s) on N 1 and N 2 positions of the imidazole.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

2025 Workshop on Envisioning Frontiers in AI and Computing for Biological Research: Position Papers

This workshop aims to identify key research directions for transforming biology using artificial intelligence (AI), machine learning (ML) and computational methods to facilitate the discovery of new behaviors, mechanisms, and designs of biological processes relevant to DOE missions, underpinning a broader U.S. bioeconomy. By developing novel AI/ML technologies to analyze and interpret complex biological data, researchers can organize and simulate biological processes at various scales as well as advance predictive understanding and manipulation of biological systems. This integration of computation, experimentation, and next-generation experimental technologies can lead to discoveries in new biological behaviors and mechanisms relevant to DOE missions. The focus is on how advanced computational and mathematical methods can impact this mission by exploring digital twins, foundation models, automated laboratory experiments, modeling of complex living systems, and data-driven approaches for the biodesign of plants and microbial systems. While data management is important, it is not the primary focus of this workshop, which will assess the current state, trends, and AI/ML challenges at the interface between biology and computational science to identify opportunities for high-impact research at their intersection. The goal is to define research needs and opportunities that align with biological sciences, computational sciences, and applied mathematics research.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

The impact of physiologically relevant temperatures on physical properties of thylakoid membranes: a molecular dynamics study

Thylakoid membranes are energy-converting membranes with a unique lipid composition. Though the membranes are primarily composed of proteins, their photosynthetic function is strongly influenced by the lipid constituents. Here we characterize, with molecular dynamics (MD) simulations, lipid bilayers with compositions representative of plant thylakoid membranes. We determine, in a wide range of temperatures, the physical parameters of the model membranes which are relevant for the photosynthetic function. We found a marked impact of temperature on membrane permeability due to a combination of increased compressibility and curvature of the membrane at elevated temperatures. With increasing temperatures, we observed increasingly smeared transmembrane density profiles of the membrane forming lipid headgroups predicting increased membrane flexibility. The diffusion coefficient of the lipids increased with temperature without apparent specificity for lipid species. Instead of a comprehensive experimental dataset in the relevant temperature range, we quantitatively compared and validated our MD results with MD simulations on a dipalmitoylphosphatidylcholine model system.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Genomes to ecosystem function: Targeting critical knowledge gaps in methanogenesis and translation to updated global biogeochemical models

Natural freshwater temperate wetland systems currently represent the largest natural source of atmospheric methane, but are relatively understudied using systems biology tools (e.g. meta-omics) compared to other high producing methane systems (e.g. peat, tropical, or reconstructed wetlands). Using field investigations at the NOAA operated sentinel site on Lake Erie, methane producing activities and responses to geochemical conditions will be determined along seasonal and spatial gradients (Objective 1). Here, a combination of high-throughput activity and gas measurements, combined with high-resolution systems biology and analytical methods, will provide in depth knowledge of the microbiological, chemical, and physical constraints on methane production in wetlands. Using laboratory microcosms, the formation of anoxic microsites and their capacity to facilitate methane production in wetland soils will be simulated (Objective 2). This objective will validate the findings from the field investigations, offering a more controlled environment for teasing out the role of different yet interrelated variables. Lastly, these field and laboratory data will be used for multi-scale, process-level evaluation of an ecosystem biogeochemical model that accommodates these newly identified processes and parameterizes representation of these processes along relevant environmental gradients (Objective 3).This research will identify multiple interacting geochemical, ecological, and metabolic constraints that are poorly understood, oversimplified, or missing in global biogeochemical methane models. This proposal targets the role of oxygen limitation on methane processes in soil domains, to improve reactive transport models of microbial carbon cycling across terrestrial-aquatic soils and generate data on nutrient cycling activities in Great Lake wetlands. This information could provide new insights into the microbial controllers of Lake Erie eutrophication.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Deriving spatially explicit direct and indirect interaction networks from animal movement data

Abstract Quantifying spatiotemporally explicit interactions within animal populations facilitates the understanding of social structure and its relationship with ecological processes. Data from animal tracking technologies (Global Positioning Systems [“GPS”]) can circumvent longstanding challenges in the estimation of spatiotemporally explicit interactions, but the discrete nature and coarse temporal resolution of data mean that ephemeral interactions that occur between consecutive GPS locations go undetected. Here, we developed a method to quantify individual and spatial patterns of interaction using continuous‐time movement models (CTMMs) fit to GPS tracking data. We first applied CTMMs to infer the full movement trajectories at an arbitrarily fine temporal scale before estimating interactions, thus allowing inference of interactions occurring between observed GPS locations. Our framework then infers indirect interactions—individuals occurring at the same location, but at different times—while allowing the identification of indirect interactions to vary with ecological context based on CTMM outputs. We assessed the performance of our new method using simulations and illustrated its implementation by deriving disease‐relevant interaction networks for two behaviorally differentiated species, wild pigs ( Sus scrofa ) that can host African Swine Fever and mule deer ( Odocoileus hemionus ) that can host chronic wasting disease. Simulations showed that interactions derived from observed GPS data can be substantially underestimated when temporal resolution of movement data exceeds 30‐min intervals. Empirical application suggested that underestimation occurred in both interaction rates and their spatial distributions. CTMM‐Interaction method, which can introduce uncertainties, recovered majority of true interactions. Our method leverages advances in movement ecology to quantify fine‐scale spatiotemporal interactions between individuals from lower temporal resolution GPS data. It can be leveraged to infer dynamic social networks, transmission potential in disease systems, consumer–resource interactions, information sharing, and beyond. The method also sets the stage for future predictive models linking observed spatiotemporal interaction patterns to environmental drivers.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

The matrix product approximation for the dynamic cavity method

Stochastic dynamics of classical degrees of freedom, defined on vertices of locally tree-like graphs, can be studied in the framework of the dynamic cavity method which is exact for tree graphs. Such models correspond for example to spin-glass systems, Boolean networks, neural networks, and other technical, biological, and social networks. The central objects in the cavity method are edge messages—conditional probabilities of two vertex variable trajectories. In this paper, we discuss a rather pedagogical derivation for the dynamic cavity method, give a detailed account of the novel matrix product edge message (MPEM) algorithm for the solution of the dynamic cavity equation as introduced, and present optimizations and extensions. Matrix product approximations of the edge messages are constructed recursively in an iteration over time. Computation costs and precision can be tuned by controlling the matrix dimensions of the MPEM in truncations. Without truncations, the dynamics is exact. Data for Glauber–Ising dynamics shows a linear growth of computation costs in time. In contrast to Monte Carlo simulations, the approach has a much better error scaling. Hence, it gives for example access to low probability events and decaying observables like temporal correlations. Here, we discuss optimized truncation schemes and an extension that allows to capture models which have a continuum time limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

NGEE Arctic Authorship Guidelines

Authorship Guidelines were developed to help facilitate trust among team members as we span multiple institutions, scientific disciplines, and career stages. NGEE Arctic was built on a foundation of open science, data sharing, and collaboration. In Phase 4 of the project, it was particularly important to keep this foundation in mind as we develop new collaborations across the Arctic. Included in this package is one *.pdf. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

Iversen, Colleen [ORNL] (ORCID:0000000182933450)↗

DOE OSTI2026

Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

Modeling↗

DOE OSTI2025

Soil oxygen dynamics: a key mediator of tile drainage impacts on coupled hydrological, biogeochemical, and crop systems

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Dominant controls of cold-season precipitation variability over the high mountains of Asia

A robust understanding of the sub-seasonal cold season (November–March) precipitation variability over the High Mountains of Asia (HMA) is lacking. Here, we identify dynamic and thermodynamic pathways through which natural modes of climate variability establish their teleconnections over the HMA. First, we identify evaporative sources that contribute to the cold season precipitation over the HMA and surrounding areas. The predominant moisture contribution comes from the mid-latitude regions, including the Mediterranean/Caspian Seas and Mediterranean land. Second, we establish that several tropical and extratropical forcings display a sub-seasonally fluctuating influence on precipitation distribution over the region during the cold season. Many of them varyingly interact, so their impacts cannot be explained independently or at seasonal timescales. Lastly, a single set of evaporative sources is not identifiable as the key determinant in propagating a remote teleconnection because the sources of moisture anomalies depend on the pattern of sub-seasonally varying dynamical forcing in the atmosphere.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Proteomic and metabolic disturbances in lignin-modified Brachypodium distachyon

Abstract Lignin biosynthesis begins with the deamination of phenylalanine and tyrosine (Tyr) as a key branch point between primary and secondary metabolism in land plants. Here, we used a systems biology approach to investigate the global metabolic responses to lignin pathway perturbations in the model grass Brachypodium distachyon. We identified the lignin biosynthetic protein families and found that ammonia-lyases (ALs) are among the most abundant proteins in lignifying tissues in grasses. Integrated metabolomic and proteomic data support a link between lignin biosynthesis and primary metabolism mediated by the ammonia released from ALs that is recycled for the synthesis of amino acids via glutamine. RNA interference knockdown of lignin genes confirmed that the route of the canonical pathway using shikimate ester intermediates is not essential for lignin formation in Brachypodium, and there is an alternative pathway from Tyr via sinapic acid for the synthesis of syringyl lignin involving yet uncharacterized enzymatic steps. Our findings support a model in which plant ALs play a central role in coordinating the allocation of carbon for lignin synthesis and the nitrogen available for plant growth. Collectively, these data also emphasize the value of integrative multiomic analyses to advance our understanding of plant metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Semantic Property Graph for Scalable Knowledge Graph Analytics

Graphs are a natural and fundamental representation to describe entities, relationships, activities, and evolution of complex systems. Many domains such as communication, citation, procurement, biology, social media, and transportation can be modeled as a set of entities and their relationships. Resource Description Framework (RDF) and Labeled Property Graph (LPG) are two of the most used data models to encode information in a graph. Both models are similar in terms of using basic graph elements such as nodes and edges but differ in terms of the modeling approach, expressibility, serialization, and target applications. RDF is a flexible data exchange model for expressing information about entities but it tends to a have high memory footprint and inefficient storage, which does not make it a natural choice to perform scalable graph analytics. In contrast, LPG has gained traction as a reliable model to perform scalable graph analytic tasks such as sub-graph matching, network alignment, and real-time knowledge graph query. It provides efficient storage, fast traversal, and flexibility to model various real-world domains. At the same time, the LPG lacks the support of a formal knowledge representation such as an ontology to provide automated knowledge inference. We propose Semantic Property Graph (SPG) as a logical projection of reified RDF into the LPG model. SPG continues to use RDF ontology to define the type hierarchy of the projected graph and validate it against a given ontology. We present a framework to convert reified RDF graphs into SPG using two different computing environments. We also present cloud-based graph migration capabilities using Amazon Web Services.

Purohit, Sumit↗

DOE OSTI2022

ITeM: Independent temporal motifs to summarize and compare temporal networks

We report networks are a fundamental and flexible way of representing various complex systems. Many domains such as communication, citation, procurement, biology, social media, and transportation can be modeled as a set of entities and their relationships. Temporal networks are a specialization of general networks where every relationship occurs at a discrete time. The temporal evolution of such networks is as important to understand as the structure of the entities and relationships. We present the Independent Temporal Motif (ITeM) to characterize temporal graphs from different domains. ITeMs can be used to model the structure and the evolution of the graph. In contrast to existing work, ITeMs are edge-disjoint directed motifs that measure the temporal evolution of ordered edges within the motif. For a given temporal graph, we produce a feature vector of ITeM frequencies and the time it takes to form the ITeM instances. We apply this distribution to measure the similarity of temporal graphs. We show that ITeM has higher accuracy than other motif frequency-based approaches. We define various ITeM-based metrics that reveal salient properties of a temporal network. We also present importance sampling as a method to efficiently estimate the ITeM counts. We present a distributed implementation of the ITeM discovery algorithm using Apache Spark and GraphFrame. We evaluate our approach on both synthetic and real temporal networks.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge or artificial intelligence driven, and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary- structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models, which contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly, while refining them against the density maps. Here, we introduce such an adaptive decision making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, and the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two protein systems, Adenylate Kinase and Carbon Monoxide Dehydrogenase. For the test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2 - 3 Å, map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that is inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from Chimpanzee adenovirus. The overhead for decision making remaining low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install the R-MDFF on different computing environments, from local Linux based workstations to High Performance Computing (HPC) environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Generative design of de novo proteins based on secondary-structure constraints using an attention-based diffusion model

We report two generative deep-learning models that predict amino acid sequences and 3D protein structures on the basis of secondary-structure design objectives via either the overall content or the per-residue structure. Both models are robust regarding imperfect inputs and offer de novo design capacity because they can discover new protein sequences not yet discovered from natural mechanisms or systems. The residue-level secondary-structure design model generally yields higher accuracy and more diverse sequences. These findings suggest unexplored opportunities for protein designs and functional outcomes within the vast amino acid sequences beyond known proteins. Our models, based on an attention-based diffusion model and trained on a dataset extracted from experimentally known 3D protein structures, offer numerous downstream applications in the conditional generative design of various biological or engineering systems. Future work could include additional conditioning and an exploration of other functional properties of the generated proteins for various properties beyond structural objectives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Cadmium alters whole animal ionome and promotes the re-distribution of iron in intestinal cells of Caenorhabditis elegans

The chronic exposure of humans to the toxic metal cadmium (Cd), either occupational or from food and air, causes various diseases, including neurodegenerative conditions, dysfunction of vital organs, and cancer. While the toxicology of Cd and its effect on the homeostasis of biologically relevant elements is increasingly recognized, the spatial distribution of Cd and other elements in Cd toxicity-caused diseases is still poorly understood. Here, we use Caenorhabditis elegans as a non-mammalian multicellular model system to determine the distribution of Cd at the tissue and cellular resolution and its effect on the internal levels and the distribution of biologically relevant elements. Using inductively coupled plasma-mass spectrophotometry (ICP-MS), we show that exposure of worms to Cd not only led to its internal accumulation but also significantly altered the C. elegans ionome. Specifically, Cd treatment was associated with increased levels of toxic elements such as arsenic (As) and rubidium (Rb) and a decreased accumulation of essential elements such as zinc (Zn), copper (Cu), manganese (Mn), calcium (Ca), cobalt (Co) and, depending on the Cd-concentration used in the assay, iron (Fe). We regarded these changes as an ionomic signature of Cd toxicity in C. elegans. We also show that supplementing nematode growth medium with Zn but not Cu, rescues Cd toxicity and that mutant worms lacking Zn transporters CDF-1 or SUR-7, or both are more sensitive to Cd toxicity. Finally, using synchrotron X-Ray fluorescence Microscopy (XRF), we showed that Cd significantly alters the spatial distribution of mineral elements. The effect of Cd on the distribution of Fe was particularly striking: while Fe was evenly distributed in intestinal cells of worms grown without Cd, in the presence of Cd, Fe, and Cd co-localized in punctum-like structures in the intestinal cells. Together, this study advances our understanding of the effect of Cd on the accumulation and distribution of biologically relevant elements. Considering that C. elegans possesses the principal tissues and cell types as humans, our data may have important implications for future therapeutic developments aiming to alleviate Cd-related pathologies in humans.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Quantifying microbial roles in environmental iron oxidation via an integrated kinetics, `omics and metabolic modeling study (Final Report)

Iron oxyhydroxides are extremely reactive components of environmental systems, and therefore exert a strong influence on biogeochemical cycles. These oxyhydroxides strongly adsorb many biologically-relevant elements, including organic carbon and phosphate, as well as a wide range of metals including uranium and actinide species. Thus, the formation mechanism of iron oxyhydroxides are key to understanding both nutrient and contaminant cycling. Microorganisms can catalyze iron oxidation and promote the formation of Fe biominerals and thus are increasingly recognized as important players in biogeochemical cycling. However, it is completely unknown how much of environmental iron oxidation is biologically mediated versus abiotic, and various challenges in studying microbial iron oxidation have hindered accurate incorporation into hydrobiogeochemical models. The overarching goal of our work was to quantify and constrain microbial iron oxidation rates and use ‘omics to gain insight into the controls on this process, while developing tools to enable integration of biotic iron oxidation into hydrobiogeochemical models. Our work focused on the Savannah River Site (SRS) in South Carolina, where extensive microbial iron oxidation has been observed. At Tims Branch, part of the Argonne National Laboratory Wetland Hydrobiogeochemistry Science Focus Area (Argonne SFA), where groundwater discharges into a stream, iron-oxidizing microbial mats form and appear to be a major sink of uranium. In the wetlands that surround Tims Branch, there are wide swaths of iron microbial mats and flocs (mobilized mat). We measured biotic and abiotic iron oxidation rates using mats and water sampled from these sites and found that iron oxidation is primarily carried out by chemolithotrophic microorganisms. The resulting rate constants can be incorporated into models. These mats were characterized by metagenomics and metatranscriptomics, which showed that aerobic chemolithotrophs were the dominant iron-oxidizing bacteria (FeOB), and these included Gallionellaceae and Leptothrix, and possibly Rhodoferax, which is known as an Fe-reducer but may also oxidize Fe(II). This demonstrated that diverse FeOB can coexist and suggests that there are a range of niches and therefore drivers of chemolithotrophic iron oxidation. Analysis of reconstructed genomes strongly suggests that a major factor in diversity is carbon source, as genomes contained varied pathways for autotrophy and heterotrophy. We performed an in-depth analysis of Leptothrix ochracea genomes, since this sheath-former is one of the primary mat builders, yet its physiology remained unresolved. A combination of genomics, transcriptomics, and metabolic modeling suggest that L. ochracea grows mixotrophically using a combination of Fe(II) and organics for energy and both inorganic and organic carbon to create biomass. This contrasts with the largely autotrophic Gallionellaceae (Gallionella, Sideroxydans, and Ferriphaselus) also present in the mats and flocs. Remarkably, multiple FeOB, both Leptothrix and Gallionellaceae, showed activity in response to Fe(II) in live mat incubations. We tracked the gene expression of individual MAGs to Fe(II) and found that various autotrophic and heterotrophic FeOB responded to Fe(II), increasing expression of both carbon fixation and organic utilization genes. The results of the integrated field, kinetics, and omics studies give detailed insight into 1) the taxa that oxidize Fe, and 2) how they connect Fe, C, and N cycles. Towards the goal of connecting omics data to hydrobiogeochemical models, we worked with the KBase team to create a template metabolic model for chemolithotrophic iron oxidation. We initially modeled the well-characterized isolate Gallionellaceae Sideroxydans lithotrophicus, and also applied the model to the mixotroph L. ochracea. In all, we have characterized diverse FeOB in a representative wetland system and solved key problems that enable better incorporation of iron-oxidizing microbes into hydrobiogeochemical models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Covalent Drug Binding in Live Cells Monitored by Mid-Infrared Quantum Cascade Laser Spectroscopy: Photoactive Yellow Protein as a Model System

The detection of drug-target interactions in live cells enables analysis of therapeutic compounds in a native cellular environment. Recent advances in spectroscopy and molecular biology have facilitated the development of genetically encoded vibrational probes like nitriles that can sensitively report on molecular interactions. Nitriles are powerful tools for measuring electrostatic environments within condensed media like proteins, but such measurements in live cells have been hindered by low signal-to-noise ratios. In this study, we design a spectrometer based on a double-beam quantum cascade laser (QCL)-based transmission infrared (IR) source with balanced detection that can significantly enhance sensitivity to nitrile vibrational probes embedded in proteins within cells compared to a conventional FTIR spectrometer. Here, using this approach, we detect small-molecule binding in Escherichia coli, with particular focus on the interaction between para-Coumaric acid (pCA) and nitrile-incorporated photoactive yellow protein (PYP). This system effectively serves as a model for investigating covalent drug binding in a cellular environment. Notably, we observe large spectral shifts of up to 15 cm –1 for nitriles embedded in PYP between the unbound and drug-bound states directly within bacteria, in agreement with observations for purified proteins. Such large spectral shifts are ascribed to the changes in the hydrogen-bonding environment around the local environment of nitriles, accurately modeled through high-level molecular dynamics simulations using the AMOEBA force field. Our findings underscore the QCL spectrometer’s ability to enhance sensitivity for monitoring drug–protein interactions, offering new opportunities for advanced methodologies in drug development and biochemical research.

chromophores↗