Search NASA⌕ Search

SEARCH · Search NASA

Results for “cell biology 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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Redox Poise during Rhodospirillum rubrum Phototrophic Growth Drives Large-scale Changes in Macromolecular Synthesis Pathways

During photoheterotrophic growth on organic substrates, purple nonsulfur photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by multiple means, including oxidation of organic substrates, oxidation of inorganic electron donors (e.g., H2), and by reverse electron flow from the photosynthetic electron transport chain. These electrons are stored as reduced electron-carrying cofactors (e.g., NAD(P)H and ferredoxin). The overall ratio of oxidized to reduced cofactors (e.g., NAD(P)+:NAD(P)H), or ’redox poise’, is difficult to understand or predict, as are the cellular processes for dissipating these reducing equivalents. Using physics-based models that capture mass action kinetics consistent with the thermodynamics of reactions and pathways, a range of redox conditions for heterophototrophic growth are evaluated, from conditions in which the NADP+/NADPH levels approach thermodynamic equilibrium to conditions in which the NADP+/NADPH ratio is far above the typical physiological values. Modeling predictions together with experimental measurements indicate that the redox poise of the cell results in large-scale changes in the activity of biosynthetic pathways and, thus, changes in cell macromolecule levels (DNA, RNA, proteins, and fatty acids). Furthermore, modeling predictions indicate that during phototrophic growth, reverse electron flow from the quinone pool is a minor contributor to the production of reduced cofactors (e.g., NAD(P)H) compared to other oxidative processes (H2 and carbon substrate oxidation). Instead, the quinone pool primarily operates to aid ATP production. The high level of ATP, in turn, drives reduction processes even when NADPH levels are relatively low compared to NADP+ by coupling ATP hydrolysis to the reductive processes. The model, in agreement with experimental measurements of macromolecule ratios of cells growing on different carbon substrates, indicates that the dynamics of nucleotide versus lipid and protein production is likely a significant mechanism of balancing oxidation and reduction in the cell.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A bromodomain-independent mechanism of gene regulation by the BET inhibitor JQ1: direct activation of nuclear receptor PXR

Abstract Bromodomain and extraterminal (BET) proteins are extensively studied in multiple pathologies, including cancer. BET proteins modulate transcription of various genes, including those synonymous with cancer, such as MYC. Thus, BET inhibitors are a major area of drug development efforts. (+)-JQ1 (JQ1) is the prototype inhibitor and is a common tool to probe BET functions. While showing therapeutic promise, JQ1 is not clinically usable, partly due to metabolic instability. Here, we show that JQ1 and the BET-inactive (−)-JQ1 are agonists of pregnane X receptor (PXR), a nuclear receptor that transcriptionally regulates genes encoding drug-metabolizing enzymes such as CYP3A4, which was previously shown to oxidize JQ1. A PXR-JQ1 co-crystal structure identified JQ1′s tert-butyl moiety as a PXR anchor and explains binding by (−)-JQ1. Analogs differing at the tert-butyl lost PXR binding, validating our structural findings. Evaluation in liver cell models revealed both PXR-dependent and PXR-independent modulation of CYP3A4 expression by BET inhibitors. We have characterized a non-BET JQ1 target, a mechanism of physiological JQ1 instability, a biological function of (−)-JQ1, and BET-dependent transcriptional regulation of drug metabolism genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Improved recovery of cell-cycle gene expression in Saccharomyces cerevisiae from regulatory interactions in multiple omics data

Gene expression is regulated by DNA-binding transcription factors (TFs). Together with their target genes, these factors and their interactions collectively form a gene regulatory network (GRN), which is responsible for producing patterns of transcription, including cyclical processes such as genome replication and cell division. However, identifying how this network regulates the timing of these patterns, including important interactions and regulatory motifs, remains a challenging task. Results We employed four in vivo and in vitro regulatory data sets to investigate the regulatory basis of expression timing and phase-specific patterns cell-cycle expression in Saccharomyces cerevisiae . Specifically, we considered interactions based on direct binding between TF and target gene, indirect effects of TF deletion on gene expression, and computational inference. We found that the source of regulatory information significantly impacts the accuracy and completeness of recovering known cell-cycle expressed genes. The best approach involved combining TF-target and TF-TF interactions features from multiple datasets in a single model. In addition, TFs important to multiple phases of cell-cycle expression also have the greatest impact on individual phases. Important TFs regulating a cell-cycle phase also tend to form modules in the GRN, including two sub-modules composed entirely of unannotated cell-cycle regulators ( STE12-TEC1 and RAP1-HAP1-MSN4 ). Conclusion Our findings illustrate the importance of integrating both multiple omics data and regulatory motifs in order to understand the significance regulatory interactions involved in timing gene expression. This integrated approached allowed us to recover both known cell-cycles interactions and the overall pattern of phase-specific expression across the cell-cycle better than any single data set. Likewise, by looking at regulatory motifs in the form of TF-TF interactions, we identified sets of TFs whose co-regulation of target genes was important for cell-cycle expression, even when regulation by individual TFs was not. Overall, this demonstrates the power of integrating multiple data sets and models of interaction in order to understand the regulatory basis of established biological processes and their associated gene regulatory networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Cryo-electron tomography provides topological insights into mutant huntingtin exon 1 and polyQ aggregates

Huntington disease (HD) is a neurodegenerative trinucleotide repeat disorder caused by an expanded poly-glutamine (polyQ) tract in the mutant huntingtin (mHTT) protein. The formation and topology of filamentous mHTT inclusions in the brain (hallmarks of HD implicated in neurotoxicity) remain elusive. Using cryo-electron tomography and subtomogram averaging, here we show that mHTT exon 1 and polyQ-only aggregates in vitro are structurally heterogenous and filamentous, similar to prior observations with other methods. Yet, we find filaments in both types of aggregates under ~2 nm in width, thinner than previously reported, and regions forming large sheets. In addition, our data show a prevalent subpopulation of filaments exhibiting a lumpy slab morphology in both aggregates, supportive of the polyQ core model. This provides a basis for future cryoET studies of various aggregated mHTT and polyQ constructs to improve their structure-based modeling as well as their identification in cells without fusion tags.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A network-enabled pipeline for gene discovery and validation in non-model plant species

Identifying key regulators of important genes in non-model crop species is challenging due to limited multi-omics resources. To address this, we introduce the network-enabled gene discovery pipeline NEEDLE, a user-friendly tool that systematically generates coexpression gene network modules, measures gene connectivity, and establishes network hierarchy to pinpoint key transcriptional regulators from dynamic transcriptome datasets. After validating its accuracy with two independent datasets, we applied NEEDLE to identify transcription factors (TFs) regulating the expression of cellulose synthase-like F6 ( CSLF6 ), a crucial cell wall biosynthetic gene, in Brachypodium and sorghum. Our analyses uncover regulators of CSLF6 and also shed light on the evolutionary conservation or divergence of gene regulatory elements among grass species. These results highlight NEEDLE’s capability to provide biologically relevant TF predictions and demonstrate its value for non-model plant species with dynamic transcriptome datasets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

KOGUT — Knowledge Oriented Graph Unified Transformer KOGUT implements the Relational Graph Transformer (RelGT) architecture for knowledge graph link prediction in biological domains, with a primary focus on microbial growth media prediction. While the original RelGT (arXiv:2505.10960) targets relational tables, time series, and multi-table databases, KOGUT adapts this architecture for heterogeneous biological knowledge graphs, providing first-in-class AI predictive models for microbial cultivation. Key Adaptations Beyond Original RelGT: - Knowledge Graph Focus: Applied to biological KGs with semantic node types (taxa, chemicals, media, phenotypes, environments) versus generic relational database tables, trained on the KG-Microbe knowledge graph (1.3M entities, 2.9M edges, 24 relation types). - Multimodal Node Encoding: Integrates node labels, categories, descriptions, and synonyms from KG metadata through learned embedding layers—adapting relational column features to graph node attributes with textual semantics. - Extended K-Hop Subgraph Strategy: Optimized neighborhood sampling (3-hop default, configurable up to 200 nodes) tuned for sparse biological networks, building on the original local-global attention framework with biological relation preservation. - Biolink Predicate Preservation: Type-specific transformations for 24 biological edge semantics (occurs_in, consumes, produces, has_phenotype, subclass_of) beyond standard relational foreign keys, enabling multi-relation link prediction. - Inductive Learning Support: Enables zero-shot predictions for novel taxa through feature-based embeddings (temperature, oxygen requirements, gram stain, cell shape), extending the original transductive relational benchmark scope to uncultured microorganisms. CheapSOTA Performance Optimizations (This Distribution): - VQ-EMA Centroid Attention: Vector quantization with exponential moving average for improved global context modeling (+5-10% MRR improvement). - HDF5 Precomputed Data Loading: One-time preprocessing of k-hop subgraphs to eliminate redundant graph traversals (2-5× training speedup). - Distributed Data Parallel Training: Multi-GPU support for scaling to larger knowledge graphs (tested on 4× NVIDIA A100 GPUs at NERSC Perlmutter). - Mixed Precision Training: Automatic mixed precision (AMP) for memory efficiency and faster training. Advantages Over Standard Knowledge Graph Embedding Models: Combines RelGT's proven multi-element tokenization (features, type, hop, structure) with graph-native biological representations, enabling interpretable link prediction across heterogeneous entities that standard embedding models (TransE, RotatE, ComplEx) and table-based transformers cannot directly model. Achieves near-perfect performance on microbial growth media prediction (MRR: 0.9966, Precision@1: 0.9932, Hit@10: 1.0000) while maintaining explainability through attention-based reasoning over biological pathways. Training Data: - KG-Microbe merged knowledge graph: 1,379,337 nodes, 2,960,472 edges - 24 biological relation types including taxonomic hierarchies, metabolic interactions, phenotype associations, and environmental relationships - Primary prediction task: Growth media suitability for microbial taxa (biolink:occurs_in, 50K edges) - Multi-relation capability: Predicts links for any of the 24 relation types, including chemical consumption/production, phenotype associations, and taxonomic classification Citation: Original RelGT Architecture: Dwivedi et al., "Relational Graph Transformer", arXiv:2505.10960, 2025 KOGUT Implementation: Knowledge Oriented Graph Unified Transformer for Microbial Growth Media Prediction Developed at Lawrence Berkeley National Laboratory (LBNL) Trained on NERSC Perlmutter supercomputer

Joachimiak, Marcin [Lawrence Berkeley National Lab↗
DOE OSTI2021

Dual arginine recognition of LRRK2 phosphorylated Rab GTPases

Parkinson’s-disease-associated LRRK2 is a multidomain Ser/Thr kinase that phosphorylates a subset of Rab GTPases to control their effector functions. Rab GTPases are the prime regulators of membrane trafficking in eukaryotic cells. Rabs exert their biological effects by recruitment of effector proteins to subcellular compartments via their Rab-binding domain (RBD). Effectors are modular and typically contain additional domains that regulate various aspects of vesicle formation, trafficking, fusion, and organelle dynamics. The RBD of effectors is typically an α-helical coiled coil that recognizes the GTP conformation of the switch 1 and switch 2 motifs of Rabs. LRRK2 phosphorylates Rab8a at T72 (pT72) of its switch 2 α-helix. This post-translational modification enables recruitment of RILPL2, an effector that regulates ciliogenesis in model cell lines. A newly identified RBD motif of RILPL2, termed the X-cap, has been shown to recognize the phosphate via direct interactions between an arginine residue (R132) and pT72 of Rab8a. Here, we show that a second distal arginine (R130) is also essential for phospho-Rab binding by RILPL2. Through structural, biophysical, and cellular studies, we find that R130 stabilizes the primary R132:pT72 salt bridge through favorable enthalpic contributions to the binding affinity. These findings may have implications for the mechanism by which LRRK2 activation leads to assembly of phospho-Rab complexes and subsequent control of their membrane trafficking functions in cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Evaluating Limits of Machine Learning-Assisted Raman Spectroscopy in Classification of Biological Samples

Machine learning (ML)-assisted Raman spectroscopy has become a powerful analytical tool for the classification and identification of analytes; however, technical challenges impacting its detection accuracy have not been thoroughly investigated. This study explores experimental factors affecting classification performance. Among the evaluated ML models, ML algorithms show minimal impact on classification accuracy. Instead, experimental factors, including spectral similarity between tested samples and data quality, dominate detection performance. Increases in spectral noise and spectral similarity significantly reduce classification accuracy. In well-controlled samples with low experimental noise, ML-assisted Raman spectroscopy can discriminate lipid mixtures with a composition difference of 1.85 mol %. To assess the effect of biological heterogeneity, we analyzed single-cell Raman spectra from Saccharomyces cerevisiae strains carrying single, double, or triple gene mutations. Intrinsic cell-to-cell variability introduced substantial spectral differences, severely reducing the accuracy of multiclass classification of these genetically similar strains at the single-cell level. Averaging Raman spectra across multiple cells improved classification accuracy by reducing this spectral variability. We also assess the effectiveness of transfer learning across different Raman spectrometers, specifically by applying an ML model trained on one instrument to another Raman spectrometer. Transfer learning can be improved with proper instrument calibration, highlighting the importance of instrument standardization. Overall, our results demonstrate that data quality and spectral similarity are the primary bottlenecks in ML-assisted Raman spectroscopy. Careful attention to sample preparation, data acquisition, measurement conditions, and instrument calibration is critical to achieving robust and reliable classification performance.

Fungi↗
DOE OSTI2025

A framework for challenges and solutions in biodesign research

The bioeconomy represents an advanced economic paradigm that builds upon previous agricultural, industrial, and digital economic models. It seeks to tackle critical global challenges such as resource scarcity, escalating healthcare demands, and environmental degradation. At the heart of the bioeconomy is biomanufacturing, which uses natural or engineered enzymes or cell factories built from ​biological components like promoters, terminators, regulatory sequences, reporters, and functional genes into various chassis hosts (including animal, microbial, plant, and de novo systems) to create products such as food, energy, medicine, materials, chemicals, and engineered tissue/organs. An enabler of biomanufacturing is biodesign – also known as biosystems design and closely related to synthetic biology or engineering biology. This interdisciplinary field aims to understand and predictably modify existing life forms or create entirely new biological entities/systems using rational engineering strategies and automated design tools. Through these capabilities, biodesign supports the discovery, optimization, and creation of efficient platforms for biomanufacturing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Streamlined spatial and environmental expression signatures characterize the minimalist duckweed Wolffia australiana

Single-cell genomics permits a new resolution in the examination of molecular and cellular dynamics, allowing global, parallel assessments of cell types and cellular behaviors through development and in response to environmental circumstances, such as interaction with water and the light–dark cycle of the Earth. Here, we leverage the smallest, and possibly most structurally reduced, plant, the semiaquaticWolffia australiana, to understand dynamics of cell expression in these contexts at the whole-plant level. We examined single-cell-resolution RNA-sequencing data and foundWolffiacells divide into four principal clusters representing the above- and below-water-situated parenchyma and epidermis. Although these tissues share transcriptomic similarity with model plants, they display distinct adaptations thatWolffiahas made for the aquatic environment. Within this broad classification, discrete subspecializations are evident, with select cells showing unique transcriptomic signatures associated with developmental maturation and specialized physiologies. Assessing this simplified biological system temporally at two key time-of-day (TOD) transitions, we identify additional TOD-responsive genes previously overlooked in whole-plant transcriptomic approaches and demonstrate that the core circadian clock machinery and its downstream responses can vary in cell-specific manners, even in this simplified system. Distinctions between cell types and their responses to submergence and/or TOD are driven by expression changes of unexpectedly few genes, characterizingWolffiaas a highly streamlined organism with the majority of genes dedicated to fundamental cellular processes.Wolffiaprovides a unique opportunity to apply reductionist biology to elucidate signaling functions at the organismal level, for which this work provides a powerful resource.

Biochemistry & Molecular Biology↗
DOE OSTI2022

Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification

To identify and stop active HIV transmission chains new epidemiological techniques are needed. Here, we describe the development of a multi-biomarker augmentation to phylogenetic inference of the underlying transmission history in a local population. HIV biomarkers are measurable biological quantities that have some relationship to the amount of time someone has been infected with HIV. To train our model, we used five biomarkers based on real data from serological assays, HIV sequence data, and target cell counts in longitudinally followed, untreated patients with known infection times. The biomarkers were modeled with a mixed effects framework to allow for patient specific variation and general trends, and fit to patient data using Markov Chain Monte Carlo (MCMC) methods. Subsequently, the density of the unobserved infection time conditional on observed biomarkers were obtained by integrating out the random effects from the model fit. This probabilistic information about infection times was incorporated into the likelihood function for the transmission history and phylogenetic tree reconstruction, informed by the HIV sequence data. To critically test our methodology, we developed a coalescent-based simulation framework that generates phylogenies and biomarkers given a specific or general transmission history. Testing on many epidemiological scenarios showed that biomarker augmented phylogenetics can reach 90% accuracy under idealized situations. Under realistic within-host HIV-1 evolution, involving substantial within-host diversification and frequent transmission of multiple lineages, the average accuracy was at about 50% in transmission clusters involving 5–50 hosts. Realistic biomarker data added on average 16 percentage points over using the phylogeny alone. Using more biomarkers improved the performance. Shorter temporal spacing between transmission events and increased transmission heterogeneity reduced reconstruction accuracy, but larger clusters were not harder to get right. More sequence data per infected host also improved accuracy. We show that the method is robust to incomplete sampling and that adding biomarkers improves reconstructions of real HIV-1 transmission histories. The technology presented here could allow for better prevention programs by providing data for locally informed and tailored strategies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Vitamin E Acetate Causes Softening of Pulmonary Surfactant Membrane Models

The popularity of electronic cigarettes and vaping products has launched the outbreak of a condition affecting the respiratory system of users, known as electronic-cigarette/vaping-associated lung injury (EVALI). The build-up of vitamin E acetate (VEA), a diluent of some illicit vaping oils, in the bronchoalveolar lavage of patients with EVALI provided circumstantial evidence as a target for investigation. In this work, we provide a fundamental characterization of the interaction of VEA with lung cells and pulmonary surfactant (PS) models to explore the mechanisms by which vaping-related lung injuries may be present. We first confirm the localization and uptake of VEA in pulmonary epithelial cells. Further, as PS is vitally responsible for the biophysical functions of the lungs, we explore the effect of added VEA on three increasingly complex models of PS: dipalmitoylphosphatidylcholine (DPPC), a lipid-only synthetic PS, and the biologically derived extract Curosurf. Using high-resolution techniques of small-angle X-ray scattering, small-angle neutron scattering, neutron spin–echo spectroscopy, and neutron reflectometry, we compare the molecular-scale behaviors of these membranes to the bulk viscoelastic properties of surfactant monolayer films as studied by Langmuir monolayer techniques. While VEA does not obviously alter the structure or organization of PS membranes, a consistent softening of membrane systems—regardless of compositional complexity—provides a biophysical explanation for the respiratory distress associated with EVALI and yields a new perspective on the behavior of the PS system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Nanoengineering of Complex Materials

We are primarily interested in two important biomimetic principles, namely, self-assembly based on supramolecular chemistry and the use of these structures to template in organic phases as observed in biomineralization. The projects proposed are inspired by some of our achievements in the previous funding period described above, and also by preliminary data gathered recently which we describe below. In addition to our broad objective of acquiring knowledge in supramolecular self-assembly and templating of organic–inorganic hybrid materials, we also propose here work on biomolecular structures based on DNA, peptide nucleic acids, and peptides which can serve as models for artificial materials based on biomimetic principles. The proposed work is divided into three sections, supramolecular systems, hybrid systems, and biomolecular and biological systems. In terms of functions, we are using systems that have relevance to solar cells, catalysis, dynamic/active materials, and ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Saccharide analysis of onion outer epidermal walls

Abstract Background Epidermal cell walls have special structural and biological roles in the life of the plant. Typically they are multi-ply structures encrusted with waxes and cutin which protect the plant from dehydration and pathogen attack. These characteristics may also reduce chemical and enzymatic deconstruction of the wall for sugar analysis and conversion to biofuels. We have assessed the saccharide composition of the outer epidermal wall of onion scales with different analytical methods. This wall is a particularly useful model for cell wall imaging and mechanics. Results Epidermal walls were depolymerized by acidic methanolysis combined with 2M trifluoracetic acid hydrolysis and the resultant sugars were analyzed by high-performance anion-exchange chromatography with pulsed amperometric detection (HPAEC-PAD). Total sugar yields based on wall dry weight were low (53%). Removal of waxes with chloroform increased the sugar yields to 73% and enzymatic digestion did not improve these yields. Analysis by gas chromatography/mass spectrometry (GC/MS) of per- O -trimethylsilyl (TMS) derivatives of the sugar methyl glycosides produced by acidic methanolysis gave a high yield for galacturonic acid (GalA) but glucose (Glc) was severely reduced. In a complementary fashion, GC/MS analysis of methyl alditols produced by permethylation gave substantial yields for glucose and other neutral sugars, but GalA was severely reduced. Analysis of the walls by 13 C solid-state NMR confirmed and extended these results and revealed 15% lipid content after chloroform extraction (potentially cutin and unextractable waxes). Conclusions Although exact values vary with the analytical method, our best estimate is that polysaccharide in the outer epidermal wall of onion scales is comprised of homogalacturonan (~ 50%), cellulose (~ 20%), galactan (~ 10%), xyloglucan (~ 10%) and smaller amounts of other polysaccharides. Low yields of specific monosaccharides by some methods may be exaggerated in epidermal walls impregnated with waxes and cutin and call for cautious interpretation of the results.

09 BIOMASS FUELS↗
DOE OSTI2024

Rapid, antibiotic incubation-free determination of tuberculosis drug resistance using machine learning and Raman spectroscopy

Tuberculosis (TB) is the world’s deadliest infectious disease, with over 1.5 million deaths and 10 million new cases reported anually. The causative organism Mycobacterium tuberculosis (Mtb) can take nearly 40 d to culture, a required step to determine the pathogen’s antibiotic susceptibility. Both rapid identification and rapid antibiotic susceptibility testing of Mtb are essential for effective patient treatment and combating antimicrobial resistance. Here, we demonstrate a rapid, culture-free, and antibiotic incubation-free drug susceptibility test for TB using Raman spectroscopy and machine learning. We collect few-to-single-cell Raman spectra from over 25,000 cells of the Mtb complex strain Bacillus Calmette-Guérin (BCG) resistant to one of the four mainstay anti-TB drugs, isoniazid, rifampicin, moxifloxacin, and amikacin, as well as a pan-susceptible wildtype strain. By training a neural network on this data, we classify the antibiotic resistance profile of each strain, both on dried samples and on patient sputum samples. On dried samples, we achieve >98% resistant versus susceptible classification accuracy across all five BCG strains. In patient sputum samples, we achieve ~79% average classification accuracy. We develop a feature recognition algorithm in order to verify that our machine learning model is using biologically relevant spectral features to assess the resistance profiles of our mycobacterial strains. Finally, we demonstrate how this approach can be deployed in resource-limited settings by developing a low-cost, portable Raman microscope that costs <$5,000. We show how this instrument and our machine learning model enable combined microscopy and spectroscopy for accurate few-to-single-cell drug susceptibility testing of BCG.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Study of wrap mode impact on Pseudomonas aeruginosa motion in the chemotactic field of fungi

An agent-based discrete computational model biologically calibrated to Pseudomonas aeruginosa migration is used to explore the impacts of bacterial reversals and wrap mode on the efficiency of motion in different environments, both with and without chemotaxis. It is first shown that wrap mode increases the exploration of continuous multimodal chemotactic profiles such as those produced by biologically relevant fungal networks. For cells undergoing a run-reverse pattern, it is shown that the bacteria are likely to remain at the first local chemoattractant maximal production site on a hypha they find. However, with wrap mode, the bacteria can more easily escape these local sites to further explore their neighboring environment along the fungi, suggesting that wrap mode may be beneficial for migration along the fungi in liquid. In a different set of simulations of bacterial motion close to an isolated chemotactic source, wrap mode is shown to increase the ability of a bacterium to reorient toward the source while reducing the overall motion required for similar chemotactic efficiency as a run-reverse strategy, suggesting a potential metabolic benefit. In contrast, model simulations show that wrap mode can increase the rate of dispersal of P. aeruginosa in a nonchemotactic environment.

Hansen, Austin [University of California, Riversid↗
DOE OSTI2023

An approach to learn regulation to maximize growth and entropy production rates in metabolism

Elucidating cell regulation remains a challenging task due to the complexity of metabolism and the difficulty of experimental measurements. Here we present a method for prediction of cell regulation to maximize cell growth rate while maintaining the solvent capacity of the cell. Prediction is formulated as an optimization problem using a thermodynamic framework that can leverage experimental data. We develop a formulation and variable initialization procedure that allows for computing solutions of the optimization with an interior point method. The approach is applied to photoheterotrophic growth of Rhodospirilium rubrum using ethanol as a carbon source, which has applications to biosynthesis of ethylene production. Growth is captured as the rate of synthesis of amino acids into proteins, and synthesis of nucleotide triphoshaptes into RNA and DNA. The method predicts regulation that produces a high rate of protein and RNA synthesis while DNA synthesis is reduced close to zero in agreement with production of DNA being turned off for much of the cell cycle.

59 BASIC BIOLOGICAL SCIENCES↗