Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex system biology”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Single-cell and spatial omics in plants: from cellular atlases to regulatory mechanisms

Single-cell RNA sequencing (scRNA-seq) has transformed transcriptomic studies by enabling gene expression profiling at the resolution of individual cells within and across a broad range of tissue types, revealing cellular heterogeneity that is obscured in bulk tissue transcriptomes. Over the past decade, improvements in microfluidics and library preparation have drastically increased throughput, allowing tens of thousands of cells to be assayed in a single experiment. Although initially developed in animal systems, scRNA-seq has rapidly emerged as a powerful and widely adopted approach in plant biology. Beyond transcriptomics, the integration of single-cell data with chromatin accessibility, proteomics, metabolomics, and spatial omics is enabling a system-level understanding of plant gene regulation and cellular organization. Network-based analytical frameworks further support the reconstruction of gene regulatory networks and the interpretation of complex single-cell data. In this review, we summarize the current technological landscape of plant single-cell studies, discuss key experimental and analytical challenges, and review emerging strategies for validating single-cell discoveries. We also discuss future directions in applying single-cell technologies to woody perennials plants and bioenergy-relevant crops, emphasizing their potential to accelerate the discovery of cell type-specific regulatory mechanisms underlying growth, stress resilience, and biomass production.

Li, Miaomiao [ORNL] (ORCID:0000000321326168)↗
DOE OSTI2026

Function, Structure, and Regulation of Nitrogen Fixation-like Metalloproteins for Nitrogen, Energy, Carbon, and Sulfur Metabolism

Nitrogenases (N 2 ases) and nitrogen fixation-like (NFL) systems play distinct roles in nitrogen, carbon, sulfur, and energy metabolism based on their fundamental differences in structure and metallocofactor identity. As new NFL systems have recently been identified and characterized, striking parallels and differences compared to N 2 ase structure, catalysis, and regulation have emerged. NFL systems use metallocofactors that span from simple [4Fe-4S] clusters to complex clusters akin to FeMo-co, previously only thought to occur in N 2 ase. This review describes the present state of knowledge on the function, structure, catalytic mechanisms, and regulation of NFL systems that perform distinct biological roles across all three domains of life. Recent advancements in N 2 ase spectroscopic techniques for probing metallocofactor structure and electronic states guide current and future work on how each NFL system catalyzes its specific biological reaction(s). Key knowledge gaps and needed areas of research for uncovering the specific metallocofactors and structural motifs that are at the heart of NFL system reaction specificity, along with how these systems are regulated, are discussed.

Bacteria↗
DOE OSTI2023

Optimizing the electron transport chain to sustainably improve photosynthesis

Abstract Genetically improving photosynthesis is a key strategy to boosting crop production to meet the rising demand for food and fuel by a rapidly growing global population in a warming climate. Many components of the photosynthetic apparatus have been targeted for genetic modification for improving photosynthesis. Successful translation of these modifications into increased plant productivity in fluctuating environments will depend on whether the electron transport chain (ETC) can support the increased electron transport rate without risking overreduction and photodamage. At present atmospheric conditions, the ETC appears suboptimal and will likely need to be modified to support proposed photosynthetic improvements and to maintain energy balance. Here, I derive photochemical equations to quantify the transport capacity and the corresponding reduction level based on the kinetics of redox reactions along the ETC. Using these theoretical equations and measurements from diverse C3/C4 species across environments, I identify several strategies that can simultaneously increase the transport capacity and decrease the reduction level of the ETC. These strategies include increasing the abundances of reaction centers, cytochrome b6f complexes, and mobile electron carriers, improving their redox kinetics, and decreasing the fraction of secondary quinone–nonreducing photosystem II reaction centers. I also shed light on several previously unexplained experimental findings regarding the physiological impacts of the abundances of the cytochrome b6f complex and plastoquinone. The model developed, and the insights generated from it facilitate the development of sustainable photosynthetic systems for greater crop yields.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape: Modeling Archive

This modeling archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025). This archive contains model input files and outputs from landscape-scale simulations conducted using ELM, the land model component of the Department of Energy’s Energy Exascale Earth System Model (E3SM), at the Council NGEE Arctic field site (Council Road mile marker 71) on Alaska’s Seward Peninsula. Input data and model output from two sets of ELM simulations are provided. The first set of simulations were conducted with the two default ELM Arctic plant functional types (PFTs; broadleaf deciduous boreal shrub and a C3 grass) and the second set of simulations were conducted with a set of nine Arctic-specific PFTs including nonvascular mosses and lichens, graminoids, forbs, evergreen dwarf shrubs, three height classes of deciduous shrubs (dwarf, low, and low to tall), and deciduous alder shrubs (Sulman et al., 2021). Parameter names and major parameter changes in the Arctic-specific PFT configuration are described in Sulman et al. (2021) and archived in the Sulman et al. (2021) dataset (see below). Simulations were spatially explicit, covering an approximately 6.4X3.3 km domain at the Council site with a spatial resolution of 100 m for a total of 2,112 simulated grid cells under each ELM PFT configuration. The modeling archive contains meteorological forcing (seven *.nc files and one *.txt file), a domain definition file (one *.nc files), land surface configuration files (two *.nc files), parameter files (two *.nc files), annual ELM output files spanning 1980-2014 (68 *.nc files), and a User’s Guide (*pdf file). Additional information on the provided files is in the “Modeling Archive Contents” section of the User’s Guide. Model outputs are aggregated to the column scale (i.e. PFT-specific outputs are not provided here).

Murphy, Bailey [ORNL] (ORCID:0000000203995221)↗
DOE OSTI2025

Maps of growing season gross primary production and net ecosystem exchange for Council Road Mile Marker 71, Seward Peninsula, Alaska, [2017-2023]

This data archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025a). Murphy et al. (2025a) evaluated whether incorporating observed Arctic vegetation heterogeneity into ELM, the land model of the Department of Energy’s Energy Exascale Earth System Model (E3SM), improved simulations of tundra carbon cycling. The associated model archive can be found at Murphy et al. (2025b). The study focused on the spatial patterns and net landscape-level growing season productivity and carbon uptake. As part of this evaluation, observationally derived maps of average growing season (June–August) net ecosystem exchange (NEE) and gross primary production (GPP) were developed for the same domain. These maps, which form the dataset described here, integrate eddy covariance flux tower, remote sensing, and vegetation community data to provide spatially explicit benchmarks for model evaluation. The maps provide spatially explicit estimates of average growing season NEE and GPP across 13 tundra vegetation communities within the study domain. By combining flux tower observations with Airborne Visible-Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) hyperspectral imagery and drone-based normalized difference vegetation index (NDVI), these maps capture the heterogeneity of carbon fluxes associated with different Arctic vegetation types. While they represent average seasonal conditions rather than interannual variability, the maps provide a unique dataset for evaluating model performance, comparing vegetation community contributions to landscape-scale carbon cycling, and supporting regional analyses of Arctic carbon dynamics. This data archive contains 5 m resolution maps of vegetation communities, vegetation community average growing season GPP, and vegetation community average growing season NEE (three *.tif files), a User’s Guide (*pdf file), and Table 1 of the User’s Guide displaying vegetation community coverage and average growing season NEE and GPP values (*.csv file).

Murphy, Bailey [ORNL] (ORCID:0000000203995221)↗
DOE OSTI2025

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). 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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Introduction: Neuromorphic Materials

The explosive growth in data collection and the need to process it efficiently, as well as the desire to automate increasingly complex tasks in transportation, medical care, manufacturing, security and many other fields have motivated a growing interest in neuromorphic computing. Unlike the binary, transistorbased ON/OFF logic gates and separate logic and memory functionalities employed in digital computing, neuromorphic computing is inspired by animal brains that use interconnected synapses and neurons to perform processing, storage and transmission of information at the same location, while only consuming ~20 W or less of power. Motivated by the brain’s efficiency, adaptability, self-learning and resiliency qualities, neuromorphic computing can be broadly defined as an approach to processing and storing information using hardware and algorithms inspired by models of biological neural systems. Present research in neuromorphic computing encompasses approaches that vary significantly in their degree of neuro-inspiration, from systems that only incorporate features such as asynchronous, event-driven operation or use crossbar arrays of non-volatile memory (NVM) elements to accelerate deep neural networks (DNNs), to designs that embrace the extreme parallelism, sparsity, reconfigurability, adaptability, complexity and stochasticity observed in nervous systems. The term ‘neuromorphic’ computing is often credited to Carver Mead, who in the 1980s investigated Si-based analog electronics to replicate functions of the animal retina. Earlier important advances in this field include the work of Frank Rosenblatt, who proposed the concept of the perceptron, Bernard Widrow, who used this concept to build one of the first analog neural networks, the Adaline and many other researchers (see ref. 6 for an historical perspective on neuromorphic computing). With the recent increase in the use of artificial intelligence and large language models, and rising concerns over the associated energy costs, interest in neuromorphic hardware has expanded rapidly. According to some estimates, driven largely by the drastic growth in the training use of artificial intelligence (AI) models using the current computing architectures, the energy cost of computing is projected to reach the energy supply worldwide by 2045. Furthermore, while this is not a realistic outcome, it means that, if more efficient computing technologies are not developed -- soon -- the world will soon become one where demand for energy and market constraints limit the continued increase of societal access to AI and cloud services from data centers. Data centers used for training and use of these models consume hundreds of terawatt hours of electricity, already past 4% of the US electricity demand.

Circuits↗
DOE OSTI2026

Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra: Modeling Archive

This dataset is a model archive of the paper Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra (in prep) to support a modeling study investigating how projected increases in Arctic temperature and precipitation will jointly influence hydrologic conditions in ice-rich tundra landscapes. With this dataset, this study is to address the research question: Will Arctic tundra landscapes become wetter or drier with increasing precipitation and temperature in the future when thaw-induced ground subsidence and associated microtopographic evolution are represented? The simulations focus on ice-wedge polygon tundra, a widespread form of ice-rich permafrost terrain that is highly sensitive to thaw-driven landscape change. This dataset contains model input and output data for four study watersheds in Alaska: Anaktuvuk, Utqiagvik (formerly Barrow), Brooks Foothills, and Prudhoe Bay. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.5), a physics-rich integrated surface–subsurface hydrologic model. For each watershed, ten modeling cases were performed representing two landscape evolution conditions (with subsidence and without subsidence) combined with five climate forcing scenarios derived from Shared Socioeconomic Pathways (SSP5, SSP5 with precipitation trend, SSP2, SSP2 with precipitation trend, and SSP2 with double precipitation trend). Particularly, for each watershed under the forcing SSP2 with precipitation trend, there are two additional simulations considering spatially heterogeneous subsidence distributions: one assumes randomly distributed scaling and the other includes elevation dependent distribution scaling. These simulations span 1980–2099 and include spin-up runs (1980–2009) followed by transient projections (2010–2099). To facilitate reproducibility of simulations, all datasets are organized by watershed. For each study watershed, the dataset contains: (1) Pre-partitioned mesh files for 32-core modeling (.par.32.XX), located in EACH_WATERSHED/mesh/basin; and also a non-partitioned mesh file (.exo) located in EACH_WATERSHED/mesh; (2) Climate forcings corresponding to the five SSP scenarios (.h5), located in EACH_WATERSHED/data; (3) Final states (.h5) from column spin-up modeling used to initialize historical watershed-scale spin-up runs from 1980 to 2009, located in EACH_WATERSHED/PreSpinupHistorical; (4) Final states (.h5) of historical watershed-scale spin-up runs from 1980 to 2009 used to initialize projection runs, located in EACH_WATERSHED/Spinup_daymetERA5; (5) ATS modeling input files (.xml), located in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/inputfiles; (6) ATS modeling output files (.dat), located in in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/combined_obs; (7) For the Brooks Foothills watershed, additional spatial model outputs are provided (.h5) for selected years (2033 and 2093) used to generate spatial figures in this study, located in Brooksfoothills/EACH_SIMULATION_SCENARIO/results-WITH/WITHOUT_SUBSIDENCE-year2033/2093. All data files with suffix .h5 can be accessible through Python h5py, and all data files with suffix of .dat can be imported by Python pandas. Mesh file with .exo can be visualized through Paraview or read by Python netCDF. 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.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION↗
DOE OSTI2021

Systematic investigations on iron cycling in phosphorus/siderophore systems: Synergism or antagonism?

Synergisms between microbial exudates on Fe (hydr)oxide dissolution as an effective Fe acquisition pathway have been recently addressed and vigorously debated. However, Fe liberation mechanisms and where siderophores and phosphorus (P) coexist received little attentions. Current study systematically investigated ferrihydrite dissolution in the presence of desferrioxamine B (DFOB) (a kind of fungally-derived siderophores) and inorganic/organic phosphorus (orthophosphate, Pi; myo-inositol hexaphosphate, IHP), as a function of solute pH, reaction time and reagent content. Reacted solids were characterized by N 2 -BET adsorption, zeta (ζ) potential analysis, sequential extraction analysis (SEDEX), field emission scanning electron microscopy (FESEM), high-resolution transmission electron microscopy (HRTEM), X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), micro Raman spectroscopy and attenuated total reflectance-Fourier transform infrared spectroscopy (ATR-FTIR). Our results indicate that upon reaction with P-only or (DFOB + P) systems, interfacial complexation partially switched from monolayer bidentate-binuclear surface complexes to ternary complexes, or laterally transformed into amorphous Fe–P precipitates. The Fe-Pi complex precipitated more readily under acidic conditions, and Fe-IHP complex preferentially nucleated in neutral-alkaline environments. Phosphorus slightly promoted Fe release from minerals and fixation to the leached layer or interfacial liquid zone initially, but subsequently prevented further attacks from protons and DFOB. The co-effects of P and DFOB likely correspond to two successive scenarios: 1) DFOB is preferentially attracted toward ferrihydrite surfaces by negative electrical fields induced by adsorbed phosphorus and can act synergistically with labile P–Fe complexes, resulting in intensive temporal dissolution of Fh; 2) subsequent Fe shuttling to DFOB can be prohibited by stabilized, passive P/Fe–P layers. Finally, our results emphasize the antagonism between P compounds and siderophores (i.e., DFOB here) on ferrihydrite dissolution to improve understanding of the biologically-mediated Fe cycling in natural systems.

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

Reciprocity relations in a biologically inspired nanomagnonic system with dipolar coupling

In this study, magnetosome chains in magnetotactic bacteria present ideal nanomagnonic model systems for studying collective resonance modes of dipolar-coupled single domain particles in relation to their spatial arrangement. Using microresonator-based ferromagnetic resonance (FMR) spectroscopy, electron microscopy, and micromagnetic modeling, we here provide insights into the complex magnonic activity within a single magnetosome chain. While the angular dependence of its FMR spectrum is dominated by twofold symmetry features due to the uniaxial anisotropy of linear chain segments, we also observed an unexpected behavior such as interrupted lines and flat bands due to the intricate geometrical details of this particular chain, such as a cross-like structural anomaly where a pair of particles is oriented perpendicular to the main axis of the chain and thus breaks the prevailing axial dipolar coupling symmetry. Such a cross junction formed by four particles exhibits interesting magnonic network properties. Notably, we observe reciprocity in the sense that the spectral response of one particle to an excitation of another one is identical to the response of the latter given an excitation of the former. Furthermore, we have identified that magnonic coupling between A and B can be facilitated via a dark state, as in magnonic stimulated Raman adiabatic passage, and that this dark-state coupling can be made non-reciprocal between A and B by breaking the symmetry of the spatial arrangement of the four particles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Quantitative trait locus mapping combined with variant and transcriptome analyses identifies a cluster of gene candidates underlying the variation in leaf wax between upland and lowland switchgrass ecotypes

Switchgrass (Panicum virgatum L.) is a promising warm-season candidate energy crop. It occurs in two ecotypes, upland and lowland, which vary in a number of phenotypic traits, including leaf glaucousness. To initiate trait mapping, two F 2 mapping populations were developed by crossing two different F 1 sibs derived from a cross between the tetraploid lowland genotype AP13 and the tetraploid upland genotype VS16, and high-density linkage maps were generated. Quantitative trait locus (QTL) analyses of visually scored leaf glaucousness and of hydrophobicity of the abaxial leaf surface measured using a drop shape analyzer identified highly significant colocalizing QTL on chromosome 7K (Chr07K). Using a multipronged approach, we identified a cluster of genes including Pavir.7KG077009, which encodes a Type III polyketide synthase-like protein, and Pavir.7KG013754 and Pavir.7KG030500, two highly similar genes that encode putative acyl-acyl carrier protein (ACP) thioesterases, as strong candidates underlying the QTL. The lack of homoeologs for any of the three genes on Chr07N, the relatively low level of identity with other switchgrass KCS proteins and thioesterases, as well as the organization of the surrounding region suggest that Pavir.7KG077009 and Pavir.7KG013754/Pavir.7KG030500 were duplicated into a fast-evolving chromosome region, which led to their neofunctionalization. Furthermore, sequence analyses showed all three genes to be absent in the two upland compared to the two lowland accessions analyzed. This study provides an example of and practical guide for trait mapping and candidate gene identification in a complex genetic system by combining QTL mapping, transcriptomics and variant analysis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A comparative genomics approach identifies contact-dependent growth inhibition as a virulence determinant

Emerging evidence suggests the Pseudomonas aeruginosa accessory genome is enriched with uncharacterized virulence genes. Identification and characterization of such genes may reveal novel pathogenic mechanisms used by particularly virulent isolates. In this work, we utilized a mouse bacteremia model to quantify the virulence of 100 individual P. aeruginosa bloodstream isolates and performed whole-genome sequencing to identify accessory genomic elements correlated with increased bacterial virulence. From this work, we identified a specific contact-dependent growth inhibition (CDI) system enriched among highly virulent P. aeruginosa isolates. CDI systems contain a large exoprotein (CdiA) with a C-terminal toxin (CT) domain that can vary between different isolates within a species. Prior work has revealed that delivery of a CdiA-CT domain upon direct cell-to-cell contact can inhibit replication of a susceptible target bacterium. Aside from mediating interbacterial competition, we observed our virulence-associated CdiA-CT domain to promote toxicity against mammalian cells in culture and lethality during mouse bacteremia. Structural and functional studies revealed this CdiA-CT domain to have in vitro tRNase activity, and mutations that abrogated this tRNAse activity in vitro also attenuated virulence. Furthermore, CdiA contributed to virulence in mice even in the absence of contact-dependent signaling. Overall, our findings indicate that this P. aeruginosa CDI system functions as both an interbacterial inhibition system and a bacterial virulence factor against a mammalian host. These findings provide an impetus for continued studies into the complex role of CDI systems in P. aeruginosa pathogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Integrative Circuit-Host Modeling of a Genetic Switch in Varying Environments

Synthetic biology is advancing into a new phase where real-world applications are emphasized. There is hence an urgent need for mathematical modeling that can quantitatively describe the behaviors of genetic devices in natural, fluctuating environments. We utilize an integrative circuit-host modeling framework to examine the dynamics of a genetic switch and its host cell in varying environments. For both steady-state and transient cases, we find increasing nutrient reduces the bistability region of the phase space and eventually drives the switch from bistability to monostability. In response, cellular growth and proteome partitioning experience the same transition. Antibiotic perturbations cause the similar circuit and host responses as nutrient variations. However, one difference is the trend of growth rate, which augments with nutrient but declines with antibiotic levels. The framework provides a mechanistic scheme to account for both the dynamic and static characteristics of the circuit-host system upon environmental perturbations, underscoring the intimacy of gene circuits and their hosts and elucidating the complexity of circuit behaviors arising from environmental variations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Medial packing and elastic asymmetry stabilize the double-gyroid in block copolymers

Triply-periodic networks are among the most complex and functionally valuable self-assembled morphologies, yet they form in nearly every class of biological and synthetic soft matter building blocks. In contrast to simpler assembly motifs – spheres, cylinders, layers – networks require molecules to occupy variable local environments, confounding attempts to understand their formation. Here, we examine the double-gyroid network phase by using a geometric formulation of the strong stretching theory of block copolymer melts, a prototypical soft self-assembly system. The theory establishes the direct link between molecular packing, assembly thermodynamics and the medial map, a generic measure of the geometric center of complex shapes. We show that “medial packing” is essential for stability of double-gyroid in strongly-segregated melts, reconciling a long-standing contradiction between infinite- and finite-segregation theories. Additionally, we find a previously unrecognized non-monotonic dependence of network stability on the relative entropic elastic stiffness of matrix-forming to tubular-network forming blocks. The composition window of stable double-gyroid widens for both large and small elastic asymmetry, contradicting intuitive notions that packing frustration is localized to the tubular domains. This study demonstrates the utility of optimized medial tessellations for understanding soft-molecular assembly and packing frustration via an approach that is readily generalizable far beyond gyroids in neat block copolymers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Synthetic Biology of Plants and Microbes for Agriculture, Environment, and Future Applications

Agriculture is under pressure to provide food for a growing population and the feedstock required to drive the bioeconomy. Methods to breed and genetically modify plants are inadequate to keep pace. When engineering crops, traits are painstakingly introduced into plants one-at-a-time, combine unpredictably, and are continuously expressed. Synthetic biology is changing these paradigms with new genome construction tools, computer aided design (CAD), and artificial intelligence (AI). “Smart plants” contain circuits that respond to environmental change, alter morphology, or respond to threats. Further, the plant and associated microbes (fungi, bacteria, archaea) are now being viewed by genetic engineers as a holistic system. Historically, plant health has been enhanced by many natural and laboratory-evolved soil microbes marketed to enhance growth, provide nutrients, or confer pest/stress resistance. Synthetic biology has expanded the number of species that can be engineered, increased the complexity of engineered functions, controlled environmental release, and assembled stable consortia. New CAD tools will manage genetic engineering projects spanning multiple plant genomes (nucleus, chloroplast, mitochondrion) and the thousands of genomes of associated bacteria/fungi. Here, this review covers advanced genetic engineering techniques to drive the next agricultural revolution, as well as push plant engineering into new realms for manufacturing, infrastructure, sensing, and remediation.

Clauer, Phillip [Massachusetts Inst. of Technology↗
DOE OSTI2022

Automating Genetic Algorithm Mutations for Molecules Using a Masked Language Model

Inspired by the evolution of biological systems, genetic algorithms have been applied to generate solutions for optimization problems in a variety of scientific and engineering disciplines. For a given problem, a suitable genome representation must be defined along with a mutation operator to generate subsequent generations. Unlike natural systems which display a variety of complex rearrangements (e.g. mobile genetic elements), mutation for genetic algorithms commonly utilizes only random point-wise changes. Furthermore, generalizing beyond point-wise mutations poses a key difficulty as useful genome rearrangements depend on the representation and problem domain. To move beyond the limitations of manually defined point-wise changes, here we propose the use of techniques from masked language models to automatically generate mutations. As a first step, common subsequences within a given population are used to generate a vocabulary. The vocabulary is then used to tokenize each genome. A masked language model is trained on the tokenized data in order to generate possible rearrangements (i.e. mutations). In order to illustrate the proposed strategy, we use string representations of molecules and use a genetic algorithm to optimize for drug-likeness and synthesizability. Finally, our results show that moving beyond random point-wise mutations accelerates genetic algorithm optimization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Quantitative dissection of Agrobacterium T-DNA expression in single plant cells reveals density-dependent synergy and antagonism

Agrobacterium pathogenesis, which involves transferring T-DNA into plant cells, is the cornerstone of plant genetic engineering. As the applications that rely on Agrobacterium increase in sophistication, it becomes critical to achieve a quantitative and predictive understanding of T-DNA expression at the level of single plant cells. Here we examine if a classic Poisson model of interactions between pathogens and host cells holds true for Agrobacterium infecting Nicotiana benthamiana. Systematically challenging this model revealed antagonistic and synergistic density-dependent interactions between bacteria that do not require quorum sensing. Using various approaches, we studied the molecular basis of these interactions. To overcome the engineering constraints imposed by antagonism, we created a dual binary vector system termed ‘BiBi’, which can improve the efficiency of a reconstituted complex metabolic pathway in a predictive fashion. Our findings illustrate how combining theoretical models with quantitative experiments can reveal new principles of bacterial pathogenesis, impacting both fundamental and applied plant biology.

Alamos, Simon↗