Search NASA⌕ Search

SEARCH · Search NASA

Results for “genetic optimization”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Bioexcitons by Design: How Do We Get There?

Biological pigment-protein complexes (PPCs) exhibit a remarkable ability to tune the optical properties of biological excitons (bioexcitons) through specific pigment-protein interactions. While such fine-tuning allows natural systems (e.g., photosynthetic proteins) to carry out their native functions with near-optimal performance, native function itself is often suboptimal for applications such as biofuel production or quantum technology development. Furthermore, this perspective offers a look at near-term prospects for the rational reoptimization of PPC bioexcitons for new functions using site-directed mutagenesis. The primary focus is on the "structure-spectrum" challenge of understanding the relationships between structural features and spectroscopic properties. While recent examples demonstrate that site-directed mutagenesis can be used to tune nearly all key bioexciton parameters (e.g., site energies, inter-pigment couplings, and electronic-vibrational interactions), critical challenges remain before we achieve truly rational design of bioexciton properties.

09 BIOMASS FUELS↗
DOE OSTI2023

Monomeric crystal structure of the vaccine carrier protein CRM 197 and implications for vaccine development

CRM 197 is a genetically detoxified mutant of diphtheria toxin (DT) that is widely used as a carrier protein in conjugate vaccines. Protective immune responses to several bacterial diseases are obtained by coupling CRM 197 to glycans from these pathogens. Wild-type DT has been described in two oligomeric forms: a monomer and a domain-swapped dimer. Their proportions depend on the chemical conditions and especially the pH, with a large kinetic barrier to interconversion. A similar situation occurs in CRM 197 , where the monomer is preferred for vaccine synthesis. Despite 30 years of research and the increasing application of CRM 197 in conjugate vaccines, until now all of its available crystal structures have been dimeric. Here, CRM 197 was expressed as a soluble, intracellular protein in an Escherichia coli strain engineered to have an oxidative cytoplasm. The purified product, called EcoCRM, remained monomeric throughout crystallization. The structure of monomeric EcoCRM is reported at 2.0 Å resolution with the domain-swapping hinge loop (residues 379–387) in an extended, exposed conformation, similar to monomeric wild-type DT. The structure enables comparisons across expression systems and across oligomeric states, with implications for monomer–dimer interconversion and for the optimization of conjugation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Dynamics-Based Peptide–MHC Binding Optimization by a Convolutional Variational Autoencoder: A Use-Case Model for CASTELO

An unsolved challenge in the development of antigen-specific immunotherapies is determining the optimal antigens to target. Comprehension of antigen–major histocompatibility complex (MHC) binding is paramount toward achieving this goal. Here, we apply CASTELO, a combined machine learning-molecular dynamics (ML-MD) approach, to identify per-residue antigen binding contributions and then design novel antigens of increased MHC-II binding affinity for a type 1 diabetes-implicated system. We build upon a small-molecule lead optimization algorithm by training a convolutional variational autoencoder (CVAE) on MD trajectories of 48 different systems across four antigens and four HLA serotypes. We develop several new machine learning metrics including a structure-based anchor residue classification model as well as cluster comparison scores. ML-MD predictions agree well with experimental binding results and free energy perturbation-predicted binding affinities. Moreover, ML-MD metrics are independent of traditional MD stability metrics such as contact area and root-mean-square fluctuations (RMSF), which do not reflect binding affinity data. Finally, our work supports the role of structure-based deep learning techniques in antigen-specific immunotherapy design.

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

Vaccine Efficacy of a Replication-Competent Interferon-Expressing Porcine Reproductive and Respiratory Syndrome (PRRS) Virus Against NADC-34 Challenge

Background/Objectives: Porcine reproductive and respiratory syndrome virus (PRRSV) significantly impedes swine production due to rapid genetic variation and suppression of antiviral interferon (IFN) responses, leading to ineffective immunity. To address this, we developed IFNmix, a replication-competent PRRSV modified live vaccine (MLV) candidate co-expressing three Type I IFN subclasses (IFNα, IFNβ, IFNδ) to enhance antiviral immunity. Methods: In two independent in vivo experiments, we compared the protection of IFNmix and a commercial PRRSV MLV vaccine during challenge with a virulent PRRSV strain. Clinical signs, antibody and cytokine production, viral replication, and lung pathology in IFNmix-vaccinated pigs were compared to those of commercial PRRSV vaccines and controls. Results: Pigs vaccinated with IFNmix exhibited similar anti-PRRSV antibody development, serum viral loads, lung lesions, and cytokine responses post-challenge with the virulent NADC34 strain, with comparable or lower body temperatures and weight gain, to pigs vaccinated with the commercial vaccines. While IFNmix showed early viral load reduction compared to the commercial vaccine (Days 7–14 post-challenge), it demonstrated similar efficacy in controlling PRRSV replication and lung pathology. Conclusions: These findings suggest that IFNmix, by expressing multiple IFNs, can potentially enhance innate and adaptive immune responses, offering a promising approach to improving PRRSV vaccine efficacy. Further studies are needed to evaluate IFNmix against a broader range of PRRSV strains and to optimize its attenuation and immunogenicity.

Miller, Laura C. (ORCID:0000000289469416)↗
DOE OSTI2022

“Multiagent” Screening Improves Directed Enzyme Evolution by Identifying Epistatic Mutations

Enzyme evolution has enabled numerous advances in biotechnology and synthetic biology, yet still requires many iterative rounds of screening to identify optimal mutant sequences. This is due to the sparsity of the fitness landscape, which is caused by epistatic mutations that only offer improvements when combined with other mutations. We report an approach that incorporates diverse substrate analogues in the screening process, where multiple substrates act like multiple agents navigating the fitness landscape, identifying epistatic mutant residues without a need for testing the entire combinatorial search space. We initially validate this approach by engineering a malonyl-CoA synthetase and identify numerous epistatic mutations improving activity for several diverse substrates. The majority of these mutations would have been missed upon screening for a single substrate alone. We expect that this approach can accelerate a wide array of enzyme engineering programs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

CO2 Valorization via Rewiring Carbon Metabolic Network

Energizing the world sustainably requires revolutionizing the way we harness natural resources. New renewable technologies include the microbiological upgrading of CO2, the primary greenhouse gas that causes climate change. This project aims to develop a new bioeconomy by taking up CO2 and converting it into high-value products using an autotrophic microbe. Specifically, the team constructed an acetogenic bacterium, Clostridium ljungdahlii, and developed gas fermentation processes that convert syngas (CO2/H2/CO) to the platform compound 3-hydroxybutyrate (3HB). The project has achieved several key successes, including the design and construction of 3HB genes in the gas-fermenting bacterium, the identification of a native 3HB pathway, the conversion of syngas to 3HB through strain engineering, and the improvement of the conversion through rational and iterative strain development. The optimized strain has achieved industry-relevant production rates (0.1 g/L/h) and titers (9.2 g/L), and sustained continuous gas fermentation for 3HB production for > 50 days in a bench-top gas bioreactor. The project has both industrial and scientific implications, including the exploitation of host microbes and processes to improve industrial-scale production of value-added compounds from CO2, and the promotion of gas fermentation technology. The team is also strengthening R&D synergies with other teams and industrial partners in the CO2 consortium, developing cutting-edge approaches for rapid and rational engineering of gas-fermenting hosts, and dissecting the thermodynamics and kinetics of processes for biological syngas valorization. So far, one patent application and one research article have resulted from this work, and the team anticipates additional significant outcomes from this project.

3-hydroxybutyrate↗
DOE OSTI2024

Systems-Level Modeling for CRISPR-Based Metabolic Engineering

The CRISPR-Cas system has enabled the development of sophisticated, multigene metabolic engineering programs through the use of guide RNA-directed activation or repression of target genes. To optimize biosynthetic pathways in microbial systems, we need improved models to inform design and implementation of transcriptional programs. Recent progress has resulted in new modeling approaches for identifying gene targets and predicting the efficacy of guide RNA targeting. Genome-scale and flux balance models have successfully been applied to identify targets for improving biosynthetic production yields using combinatorial CRISPR-interference (CRISPRi) programs. Here, the advent of new approaches for tunable and dynamic CRISPR activation (CRISPRa) promises to further advance these engineering capabilities. Once appropriate targets are identified, guide RNA prediction models can lead to increased efficacy in gene targeting. Developing improved models and incorporating approaches from machine learning may be able to overcome current limitations and greatly expand the capabilities of CRISPR-Cas9 tools for metabolic engineering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Data for "Viral-mediated Delivery of Morphogenic Regulators Enables Leaf Transformation in Sorghum bicolor (L.)"

Recent advancements in monocot transformation, using leaf tissue as explant material, have expanded the number of grass species capable of transgenesis. However, the complexity of vectors and reliance on inducible excision of essential morphogenic regulators have so far limited widespread application. Plant RNA viruses, such as Foxtail Mosaic Virus (FoMV), present a unique opportunity to express morphogenic regulator genes, such as Babyboom ( Bbm ), Wuschel2 ( Wus2 ), Wuschel-like homeobox protein 2a ( Wox2a ) and the GROWTH-REGULATING FACTOR 4 (GRF4) GRF-INTERACTING FACTOR 1 (GIF1) fusion protein transiently in leaf explant tissues. Furthermore, altruistic delivery of conventional and viral vectors could provide opportunities to simplify vectors used for leaf transformation—facilitating vector optimization and reducing reliance on morphogenic regulator gene integration. In this study, both viral and conventional T-DNA vectors were tested for their ability to promote the formation of embryonic calli, a critical step in leaf transformation protocols, using Sorghum bicolor leaf explants. Although conventional leaf transformation vectors yielded viable embryonic calli (43.2 ± 2.9%: GRF4-GIF1, 50.2 ± 3%: Bbm/Wus2), altruistic conventional vectors employing the GRF4-GIF1 morphogenic regulator resulted in improved efficiencies (61.3 ± 4.7%). Altruistic delivery was further enhanced with the use of viral vectors employing both GRF4-GIF1 and Bbm/Wus2 regulators, resulting in 75.1 ± 2.3% and 79.2 ± 2.5% embryonic calli formation, respectively. Embryonic calli generated from both conventional and viral vectors produced shoots expressing fluorescent reporters, which were confirmed using molecular analysis. This work provides an important proof-of-concept for the use of both altruistic vectors and viral-expressed morphogenic regulators for improving plant transformation.

gene editing↗
DOE OSTI2025

Systematic Mapping of Bacterial CRISPRa Systems for Synergistic Gene Activation Reveals Antagonistic Effects

CRISPR gene activation (CRISPRa) tools have shown great promise for bacterial strain engineering but often require customization for each intended application. Our goal is to create generalizable CRISPRa tools that can overcome previous limitations of gene activation in bacteria. In eukaryotic cells, multiple activators can be combined for synergistic gene activation. To identify potential effectors for synergistic activation in bacteria, we systematically characterized bacterial activator proteins with a set of engineered synthetic promoters. We found that optimal target sites for different activators could vary by up to 200 bases in the region upstream of the transcription start site (TSS). These optimal target sites qualitatively matched previous reports for each activator, but the precise targeting rules varied between different promoters. By characterizing targeting rules in the same promoter context, we were able to test activator combinations with each effector positioned at its optimal target site. We did not find any activator combinations that produced synergistic activation, and we found that many combinations were antagonistic. Furthermore, this systematic investigation highlights fundamental mechanistic differences between bacterial and eukaryotic transcriptional activation systems and suggests that alternative strategies will be necessary for strong bacterial gene activation at arbitrary endogenous targets.

CRISPR activation↗
DOE OSTI2026

Data for Comparison of Genotyping Assays for Detection of Targeted CRISPR/Cas Mutagenesis in Highly Polyploid Sugarcane

Sugarcane ( Saccharum spp.) is an important biofuel feedstock and a leading source of global table sugar. Saccharum hybrid cultivars are highly polyploid (2n = 100–130), containing large numbers of functionally redundant hom(e)ologs in their genomes. Genome editing with sequence-specific nucleases holds tremendous promise for sugarcane breeding. However, identification of plants with the desired level of co-editing within a pool of primary transformants can be difficult. While DNA sequencing provides direct evidence of targeted mutagenesis, it is cost-prohibitive as a primary screening method in sugarcane and most other methods of identifying mutant lines have not been optimized for use in highly polyploid species. In this study, non-sequencing methods of mutant screening, including capillary electrophoresis (CE), Cas9 RNP assay, and high-resolution melt analysis (HRMA), were compared to assess their potential for CRISPR/Cas9-mediated mutant screening in sugarcane. These assays were used to analyze sugarcane lines containing mutations at one or more of six sgRNA target sites. All three methods distinguished edited lines from wild type, with co-mutation frequencies ranging from 2% to 100%. Cas9 RNP assays were able to identify mutant sugarcane lines with as low as 3.2% co-mutation frequency, and samples could be scored based on undigested band intensity. CE was highlighted as the most comprehensive assay, delivering precise information on both mutagenesis frequency and indel size to a 1 bp resolution across all six targets. This represents an economical and comprehensive alternative to sequencing-based genotyping methods which could be applied in other polyploid species.

Genomics↗
DOE OSTI2025

Plant genotype and rhizobia strain combinations strongly influence the transcriptome under heavy metal stress conditions in Medicago truncatula

Heavy metals such as cadmium (Cd) and mercury (Hg) pose significant threats to plant health and food safety as they are absorbed from the environment. Legumes are generally considered sensitive to heavy metals but possess standing genetic variation for accumulation and tolerance to toxic ions. We conducted a transcriptomic analysis on hydroponically and soil grown Medicago truncatula plants to investigate gene expression responses to Cd and Hg exposure in roots, leaves, and nodules. By using plant genotypes with varying metal tolerance or accumulation levels, we observed distinct clustering of gene ontologies, indicating tissue-specific, genotype-specific, and metal-specific gene expression patterns. Considering the symbiotic relationship between legumes and nitrogen-fixing bacteria, we further examined plant phenotypes and transcriptomes of plant genotypes with contrasting Hg accumulation levels and inoculated them with high or low Hg-tolerant Sinorhizobium medicae strains that have presence-absence variation for a mercury reductase (Mer) operon. Host plants inoculated with the Hg-tolerant rhizobia strain possessing a Mer operon exhibited less reduction in nodule number and plant biomass. A smaller reduction in iron (Fe) distribution in nodules after Hg stress was measured using X-ray Fluorescence (XRF) imaging. Dual transcriptome (host plant and bacteria) analysis of nodules revealed a remarkable decrease in the number of differentially expressed genes (DEGs) and clustering of gene ontologies in plants inoculated with the Hg-tolerant rhizobia strain, including symbiosis related genes. This finding suggests that the Hg-tolerant rhizobia strain has the potential to mitigate Hg stress in host plants. Furthermore, we observed genotype by-genotype interactions between the high Hg accumulating plant genotype and the Hg-tolerant rhizobia strain. These findings provide insights into enhancing plant resilience in contaminated environments through optimizing legume-rhizobia interactions for heavy metal tolerance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Ecosystems and Networks Integrated with Genes and Molecular Assemblies (ENIGMA). Final report

The goal of the Stahl lab is to use a combination of isotopic, genomic, and activity-based analyses in the field, of selected field isolates, and defined synthetic communities to identify sources and sinks of different nitrogen species in the field, focusing on the FRC at the ORNL. The most general research objective is to characterize variants in the pathway for denitrification among field isolates and synthetic microbial assemblies as a foundation for predicting the biotic and abiotic controls of nitrogen and carbon transformations in the field. Specifically, we have focused on characterizing environmental controls of nitrous oxide emissions, an important greenhouse gas and ozone depleting molecule. Studies have mapped the flux of nitrous oxide and controlling genetic elements within the FRC subsurface, indicating the importance of different enzyme variants in nitrous oxide production and consumption in controlling surface emissions. We have developed methods to quantify and catalog the isotopic signatures imparted on nitrous oxide by different nitrous oxide reductase variants in isolated organisms and model synthetic communities. We have developed methods to quantify and catalog the nitrous oxide affinity of field relevant isolates and their biomass yield when respiring nitrous oxide on different carbon substrates. We continue to expand the collection of field-relevant microbial isolates within ENIGMA by focused enrichment and isolation from field material. These objectives have been facilitated by the development of new technologies, including: 1) the development of a multiplexed automated optical density measurement system for quantifying growth kinetics of microorganisms dependent on gaseous substrates in monoculture or gas exchange in synthetic communities and 2) the development of a microcalorimetric analysis format to quantify the impact of environmental stress on microbial maintenance energy. Since field relevant processes are mediated by assemblies of organisms interacting with each other and their local environment, we have also collaborated with other ENIGMA investigators on the assembly and characterization of synthetic communities. These studies have evaluated the role of organisms lacking a complete pathway for denitrification in the control of field processes, measuring process rates of assemblies of organisms that individually lack the full denitrification pathway but mediate the complete denitrification when growing together. These efforts have so far resulted in the publication of 18 peer reviewed papers ranging from ecological and physiological drivers of pathway optimization to field-scale process characterization.

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

A Conjugated Oligoelectrolyte Exhibiting Room Temperature Spin-Correlated Radical Pair Character for Biological Sensing

We report a water-soluble conjugated oligoelectrolyte (COE) composed of carbazole-benzophenone, COE-CbzBP, that exhibits photogenerated spin-correlated radical pair (SCRP) behavior sensitive to static electric fields from DNA but not from lipid bilayers. The SCRP forms from a thermally activated, spin-polarized state enabled by partial π-conjugation disruption at the donor–acceptor (carbazole-benzophenone) nitrogen–carbon (N–C) junction, which facilitates a twisted intramolecular charge-transfer (TICT) geometry. This state minimizes the singlet–triplet energy gap (ΔE ST = 0.12 eV), radical–pair exchange coupling (J RP ∼ ΔE ST /2), and charge separation free energy (ΔG CS ) in both DNA (−0.19 eV) and lipid bilayers (−0.55 eV). Room-temperature continuous-wave electron paramagnetic resonance (CW-EPR) reveals a photogenerated spin-polarized singlet for COE-CbzBP that splits upon DNA association, consistent with modulation of J RP and hyperfine coupling (A x ), presumably via electric field-spin coupling. No spin-polarized signal was observed under dark, cryogenic conditions, or in liposomes, but was quenched by the spin trap 4-POBN. Transient absorption and spectroelectrochemistry confirmed magnetic-field sensitive long-lived excited-state absorption features attributed to charge-separated states 3 [Cbz •+ -BP •– ]*, which were lengthened by DNA, and quenched in lipid bilayers and 4-POBN. Quantum chemical simulations show that planar geometries (lipid-like) increase ΔE ST by 0.31 eV compared to TICT-optimized structures. This geometry-dependent modulation explains the absence of SCRP signatures in rigid environments, underscoring the importance of TICT states, minimized ΔE ST , and favorable ΔG CS for achieving room-temperature SCRP generation. These findings establish design principles for TICT-enabled molecules exhibiting qubit-like behavior that operate under ambient and biologically relevant conditions, with direct implications for quantum information science (QIS).

Aromatic compounds↗
DOE OSTI2022

TULIP: An RNA-seq-based Primary Tumor Type Prediction Tool Using Convolutional Neural Networks

Background: With cancer as one of the leading causes of death worldwide, accurate primary tumor type prediction is critical in identifying genetic factors that can inhibit or slow tumor progression. There have been efforts to categorize primary tumor types with gene expression data using machine learning, and more recently with deep learning, in the last several years. Methods In this paper, we developed four 1-dimensional (1D) Convolutional Neural Network (CNN) models to classify RNA-seq count data as one of 17 highly represented primary tumor types or 32 primary tumor types regardless of imbalanced representation. Additionally, we adapted the models to take as input either all Ensembl genes (60,483) or protein coding genes only (19,758). Unlike previous work, we avoided selection bias by not filtering genes based on expression values. RNA-seq count data expressed as FPKM-UQ of 9,025 and 10,940 samples from The Cancer Genome Atlas (TCGA) were downloaded from the Genomic Data Commons (GDC) corresponding to 17 and 32 primary tumor types respectively for training and validating the models. Results: All 4 1D-CNN models had an overall accuracy of 94.7% to 97.6% on the test dataset. Further evaluation indicates that the models with protein coding genes only as features performed with better accuracy compared to the models with all Ensembl genes for both 17 and 32 primary tumor types. For all models, the accuracy by primary tumor type was above 80% for most primary tumor types. Conclusions: We packaged all 4 models as a Python-based deep learning classification tool called TULIP (TUmor CLassIfication Predictor) for performing quality control on primary tumor samples and characterizing cancer samples of unknown tumor type. Further optimization of the models is needed to improve the accuracy of certain primary tumor types.

Jones, Sara↗
DOE OSTI2021

Quantifying the effect of hydride microstructure on zirconium alloys embrittlement using image analysis

Since hydride microstructure can significantly influence Zr alloy nuclear fuel cladding’s ductility, a new metric has been developed to quantify hydride microstructure in 2D micrographs and relate it to crack propagation. As cladding failure usually results from a hoop stress, this new metric, called the Radial Hydride Continuous Path (RHCP), is based on quantifying the continuity of brittle hydride particles along the radial direction of the cladding tube. Compared to previous metrics, this approach more closely relates to the propensity of a crack to propagate radially through the cladding tube thickness. The RHCP takes into account hydride length, orientation, and connectivity to choose the optimal path for crack propagation through the cladding thickness. The RHCP can therefore be more closely linked to hydride embrittlement of the Zr alloy material, thus creating a relationship between material structure, properties, and performance. The new definition, along with previously proposed metrics such as the Radial Hydride Fraction (RHF), the Hydride Continuity Coefficient (HCC), and the Radial Hydride Continuity Factor (RHCF), have been implemented and automated in MATLAB. These codes are made available with this publication as supplementary materials. These metrics were verified by comparing their predictions of hydride morphology against expected values in simple cases, and the implementation of the new metric was validated by comparing its predictions with manual measurements of hydride microstructure performed on ImageJ. The RHCP was also validated against experimental measurements of fracture behavior and it was shown to correlate with cladding failure better than previous metrics. The information provided by these metrics will help accurately assess cladding integrity during operation, transportation, and storage.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Li + Diffusion in Amorphous and Crystalline Al 2 O 3 for Battery Electrode Coatings

Al 2 O 3 is often applied protectively to lithium-ion battery anode and cathode materials to inhibit surface degradation, suppress dendrite formation, and relieve mechanical stresses. Given the very high intrinsic band gap and diffusion barrier of the material, the mechanism that allows Li diffusion through these coatings is not well understood, and widely varying laboratory results indicate that there may be dependencies on morphology and stoichiometry. Using nudged elastic band calculations and ab initio molecular dynamics, we perform a systematic investigation across Al 2 O 3 structures, both crystalline and amorphous, and at various concentrations of Li + to uncover the optimal parameters for maximally diffusive coatings. We find a correlation between the low proximity of Li + to Al 3+ and the low Li + migration barrier. Although barriers are the lowest in the highly diffusive one-dimensional channels of crystalline θ-Al 2 O 3 , the system is structurally delicate and subject to detrimental distortion as the Li + content is increased. The α-Al 2 O 3 lattice is, conversely, highly stable against distortion at all Li + concentrations but disadvantageous for Li + migration. In amorphous systems, unscreened Li + –Li + Coulomb repulsion and pre-emptive occupation of “trapping sites” combine to lower the energy barriers as a function of increasing concentration. One of our most important findings is that Al-deficient materials can sharply increase Li + movement, and we predict that an amorphous material with a combination of high Li + concentration and Al deficiency would enable highly Li + -conductive protective coatings for electrodes.

25 ENERGY STORAGE↗
DOE OSTI2021

Computational design and analysis of modular cells for large libraries of exchangeable product synthesis modules

Microbial metabolism can be harnessed to produce a large library of useful chemicals from renewable resources such as plant biomass. However, it is laborious and expensive to create microbial biocatalysts to produce each new product. To tackle this challenge, we have recently developed modular cell (ModCell) design principles that enable rapid generation of production strains by assembling a modular (chassis) cell with exchangeable production modules to achieve overproduction of target molecules. Previous computational ModCell design methods are limited to analyze small libraries of around 20 products. In this study, we developed a new computational method, named ModCell-HPC, that can design modular cells for large libraries with hundreds of products with a highly-parallel and multi-objective evolutionary algorithm and enable us to elucidate modular design properties. We demonstrated ModCell-HPC to design Escherichia coli modular cells towards a library of 161 endogenous production modules. From these simulations, we identified E. coli modular cells with few genetic manipulations that can produce dozens of molecules in a growth-coupled manner with different types of fermentable sugars. These designs revealed key genetic manipulations at the chassis and module levels to accomplish versatile modular cells, involving not only in the removal of major by-products but also modification of branch points in the central metabolism. We further found that the effect of various sugar degradation on redox metabolism results in lower compatibility between a modular cell and production modules for growth on pentoses than hexoses. To better characterize the degree of compatibility, we developed a method to calculate the minimal set cover, identifying that only three modular cells are all needed to couple with up 85 compatible production modules. By determining the unknown compatibility contribution metric, we further elucidated the design features that allow an existing modular cell to be re-purposed towards production of new molecules. Altogether, ModCell-HPC is a useful tool for understanding modularity of biological systems and guiding more efficient and generalizable design of modular cells that help reduce research and development cost in biocatalysis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Complementing a continuum thermodynamic approach to constitutive modeling with symbolic regression

Here, an interpretable machine learning method, physics-informed genetic programming-based symbolic regression (P-GPSR), is integrated into a continuum thermodynamic approach to developing constitutive models. The proposed strategy for combining a thermodynamic analysis with P-GPSR is demonstrated by generating a yield function for an idealized material with voids, i.e., the Gurson yield function. First, a thermodynamic-based analysis is used to derive model requirements that are exploited in a custom P-GPSR implementation as fitness criteria or are strongly enforced in the solution. The P-GPSR implementation improved accuracy, generalizability, and training time compared to the same GPSR code without physics-informed fitness criteria. The yield function generated through the P-GPSR framework is in the form of a composite function that describes a class of materials and is characteristically more interpretable than GPSR-derived equations. The physical significance of the input functions learned by P-GPSR within the composite function is acquired from the thermodynamic analysis. Fundamental explanations of why the implemented P-GPSR capabilities improve results over a conventional GPSR algorithm are provided.

42 ENGINEERING↗