Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “evolutionary 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Inferring Type II-P Supernova Progenitor Masses from Plateau Luminosities

Abstract Connecting observations of core-collapse supernova explosions to the properties of their massive star progenitors is a long-sought, and challenging, goal of supernova science. Recently, Barker et al. presented bolometric light curves for a landscape of progenitors from spherically symmetric neutrino-driven core-collapse supernova (CCSN) simulations using an effective model. They find a tight relationship between the plateau luminosity of the Type II-P CCSN light curve and the terminal iron-core mass of the progenitor. Remarkably, this allows us to constrain progenitor properties with photometry alone. We analyze a large observational sample of Type II-P CCSN light curves and estimate a distribution of iron-core masses using the relationship of Barker et al. The inferred distribution matches extremely well with the distribution of iron-core masses from stellar evolutionary models and namely, contains high-mass iron cores that suggest contributions from very massive progenitors in the observational data. We use this distribution of iron-core masses to infer minimum and maximum masses of progenitors in the observational data. Using Bayesian inference methods to locate optimal initial mass function parameters, we find M min = 9.8 − 0.27 + 0.37 and M max = 24.0 − 1.9 + 3.9 solar masses for the observational data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Parallel derivative-free optimization for simulation-based design of behind-the-meter energy systems

In this work, the integrated design and dispatch of behind-the-meter or distributed resources (e.g. stationary battery storage and solar PV generation) is considered. A simulation-based framework is employed, generating high-fidelity results with closed-loop predictive control at a fine resolution, at the expense of high computational cost (several minutes to a few hours per design point). To address this challenge, parallel derivative-free design methods are considered. Four methods are compared, including state-of-the-art surrogate-based methods (Radial-Basis Functions and Gaussian processes) and sampling strategies, an evolutionary-based method, and a simple sequential grid refinement method. As a case study, two types of design problem with increasing complexity are considered, namely, the design of behind-the-meter resources (three design variables) and the inclusion of grid capacity (four design variables). The second yields a constrained design problem for which violations can only be determined after solving the computationally expensive simulation. For the three-dimensional case, all methods present a good performance, achieving a solution within 1% of the optimum after the first iteration, with the sequential grid refinement exhibiting the fastest convergence and achieving the best final objective value. This indicates that the parallel evaluation of multiple sampling points may be more important than the choice of method for small decision spaces. For the four-dimensional constrained case, the Genetic Algorithm presents the best tradeoff between performance and computational effort, while the rough objective function terrain generated by constraint violation penalties reduces the performance of surrogate-based methods. Contour plots with flat regions indicate flexibility in the optimal design and highlight the importance of characterizing the solution space.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Evolutionary Divergence of Marinobacter Strains in Cryopeg Brines as Revealed by Pangenomics

Marinobacter spp. are cosmopolitan in saline environments, displaying a diverse set of metabolisms that allow them to competitively occupy these environments, some of which can be extreme in both salinity and temperature. Here, we introduce a distinct cluster of Marinobacter genomes, composed of novel isolates and in silico assembled genomes obtained from subzero, hypersaline cryopeg brines, relic seawater-derived liquid habitats within permafrost sampled near Utqiaġvik, Alaska. Using these new genomes and 45 representative publicly available genomes of Marinobacter spp. from other settings, we assembled a pangenome to examine how the new extremophile members fit evolutionarily and ecologically, based on genetic potential and environmental source. This first genus-wide genomic analysis revealed that Marinobacter spp. in general encode metabolic pathways that are thermodynamically favored at low temperature, cover a broad range of organic compounds, and optimize protein usage, e.g., the Entner–Doudoroff pathway, the glyoxylate shunt, and amino acid metabolism. The new isolates contributed to a distinct clade of subzero brine-dwelling Marinobacter spp. that diverged genotypically and phylogenetically from all other Marinobacter members. The subzero brine clade displays genomic characteristics that may explain competitive adaptations to the extreme environments they inhabit, including more abundant membrane transport systems (e.g., for organic substrates, compatible solutes, and ions) and stress-induced transcriptional regulatory mechanisms (e.g., for cold and salt stress) than in the other Marinobacter clades. We also identified more abundant signatures of potential horizontal transfer of genes involved in transcription, the mobilome, and a variety of metabolite exchange systems, which led to considering the importance of this evolutionary mechanism in an extreme environment where adaptation via vertical evolution is physiologically rate limited. Assessing these new extremophile genomes in a pangenomic context has provided a unique view into the ecological and evolutionary history of the genus Marinobacter , particularly with regard to its remarkable diversity and its opportunism in extremely cold and saline environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

An evolutionary algorithm for designing microbial communities via environmental modification

Despite a growing understanding of how environmental composition affects microbial communities, it remains difficult to apply this knowledge to the rational design of synthetic multispecies consortia. This is because natural microbial communities can harbour thousands of different organisms and environmental substrates, making up a vast combinatorial space that precludes exhaustive experimental testing and computational prediction. Here, we present a method based on the combination of machine learning and metabolic modelling that selects optimal environmental compositions to produce target community phenotypes. In this framework, dynamic flux balance analysis is used to model the growth of a community in candidate environments. A genetic algorithm is then used to evaluate the behaviour of the community relative to a target phenotype, and subsequently adjust the environment to allow the organisms to approach this target. We apply this iterative process to thousands of in silico communities of varying sizes, showing how it can rapidly identify environments that yield desired taxonomic compositions and patterns of metabolic exchange. Moreover, this combination of approaches produces testable predictions for the assembly of experimental microbial communities with specific properties and can facilitate rational environmental design processes for complex microbiomes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Kingdom-Wide Analysis of Fungal Protein-Coding and tRNA Genes Reveals Conserved Patterns of Adaptive Evolution

Protein-coding genes evolved codon usage bias due to the combined but uneven effects of adaptive and nonadaptive influences. Studies in model fungi agree on codon usage bias as an adaptation for fine-tuning gene expression levels; however, such knowledge is lacking for most other fungi. Our comparative genomics analysis of over 450 species supports codon usage and transfer RNAs (tRNAs) as coadapted for translation speed and this is most likely a realization of convergent evolution. Rather than drift, phylogenetic reconstruction inferred adaptive radiation as the best explanation for the variation of interspecific codon usage bias. Although the phylogenetic signals for individual codon and tRNAs frequencies are lower than expected by genetic drift, we found remarkable conservation of highly expressed genes being codon optimized for translation by the most abundant tRNAs, especially by inosine-modified tRNAs. As an application, we present a sequence-to-expression neural network that uses codons to reliably predict highly expressed transcripts. The kingdom Fungi, with over a million species, includes many key players in various ecosystems and good targets for biotechnology. Collectively, our results have implications for better understanding the evolutionary success of fungi, as well as informing the biosynthetic manipulation of fungal genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Kinetic modeling of microbial growth, enzyme activity, and gene deletions: An integrated model of β-glucosidase function in Cellvibrio japonicus

Understanding the complex growth and metabolic dynamics in microorganisms requires advanced kinetic models containing both metabolic reactions and enzymatic regulation to predict phenotypic behaviors under different conditions and perturbations. Most current kinetic models lack gene expression dynamics and are separately calibrated to distinct media, which consequently makes them unable to account for genetic perturbations or multiple substrates. This challenge limits our ability to gain a comprehensive understanding of microbial processes towards advanced metabolic optimizations that are desired for many biotechnology applications. Here, we present an integrated computational and experimental approach for the development and optimization of mechanistic kinetic models for microbial growth and metabolic and enzymatic dynamics. Our approach integrates growth dynamics, gene expression, protein secretion, and gene-deletion phenotypes. We applied this methodology to build a dynamic model of the growth kinetics in batch culture of the bacterium Cellvibrio japonicus grown using either cellobiose or glucose media. The model parameters were inferred from an experimental data set using an evolutionary computation method. The resulting model was able to explain the growth dynamics of C. japonicus using either cellobiose or glucose media and was also able to accurately predict the metabolite concentrations in the wild-type strain as well as in β-glucosidase gene deletion mutant strains. We validated the model by correctly predicting the non-diauxic growth and metabolite consumptions of the wild-type strain in a mixed medium containing both cellobiose and glucose, made further predictions of mutant strains growth phenotypes when using cellobiose and glucose media, and demonstrated the utility of the model for designing industrially-useful strains. Importantly, the model is able to explain the role of the different β-glucosidases and their behavior under genetic perturbations. Here, this integrated approach can be extended to other metabolic pathways to produce mechanistic models for the comprehensive understanding of enzymatic functions in multiple substrates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Deep Reinforcement Learning Based Control of Wind Turbines for Fast Frequency Response

In order to fulfill vital auxiliary grid services, such as load regulation, spin and non-spin reserve provision, and frequency support during emergencies, there is often a requirement for certain wind farms to operate in de-loaded modes. Leveraging the swift response capabilities of wind farms, this study demonstrates that reserving power in de-loaded modes can significantly enhance power grid stability and reliability during system contingencies. Controlling wind farms optimally for frequency support is intricate due to the nonlinearity of models and controllers and the complexity of wind farm interactions with power systems. Here, to address this challenge, this paper introduces a novel approach that integrates wind turbines into reinforcement learning-based solutions for frequency response. This innovative methodology utilizes the state-of-the-art reinforcement learning algorithm known as the surrogate-gradient-based evolutionary strategy. The proposed learning-based algorithm provides continuous control of wind farm output to rapidly stabilize system frequency and prevent unnecessary trips of under-frequency load shedding relays. To facilitate efficient training, parallel computing techniques are employed. The proposed methodology is evaluated on a modified IEEE-39 bus system, and simulation results reveal its efficacy in reliably supporting power system frequency and preventing the need for unnecessary load shedding.

Gao, Wei [Argonne National Laboratory (ANL), Argon↗
DOE OSTI2024

Nitrogen acquisition and retention pathways in sustainable perennial bioenergy grass cropping systems

Abstract Perennial tall grasses show promise as bioenergy crops due to high productivity and efficient nutrient use. Ongoing research on bioenergy grasses seeks to reduce their reliance on synthetic nitrogen (N) fertilizer, the manufacture of which relies on fossil fuel combustion. Excessive use of fertilizers also causes adverse environmental consequences and leads to the evolutionary loss of plant traits beneficial to sustainable N cycle. Notably, perennial tall grasses have exhibited the potential to maintain high biomass yield without the need for N fertilizer or causing soil N depletion. Perennial grasses can be adept at interacting with their microbial partners to facilitate N acquisition and retention via mechanisms such as biological N fixation and nitrification inhibition. These inherent N management traits should be preserved and optimized at the this early stage of bioenergy grass breeding programs. This review examines the impact of external N on bioenergy grass production and explores the potential of leveraging advantageous N‐cycling attributes of perennial tall grasses, laying groundwork for future management and research efforts. With minimized dependency on external N input, the cultivation of perennial energy grasses will pave the way toward more resilient agricultural systems and play a significant role in addressing key global energy and environmental challenges.

Agriculture↗
DOE OSTI2021

VirION2: a short- and long-read sequencing and informatics workflow to study the genomic diversity of viruses in nature

Microbes play fundamental roles in shaping natural ecosystem properties and functions, but do so under constraints imposed by their viral predators. However, studying viruses in nature can be challenging due to low biomass and the lack of universal gene markers. Though metagenomic short-read sequencing has greatly improved our virus ecology toolkit—and revealed many critical ecosystem roles for viruses—microdiverse populations and fine-scale genomic traits are missed. Some of these microdiverse populations are abundant and the missed regions may be of interest for identifying selection pressures that underpin evolutionary constraints associated with hosts and environments. Though long-read sequencing promises complete virus genomes on single reads, it currently suffers from high DNA requirements and sequencing errors that limit accurate gene prediction. Here we introduce VirION2, an integrated short- and long-read metagenomic wet-lab and informatics pipeline that updates our previous method (VirION) to further enhance the utility of long-read viral metagenomics. Using a viral mock community, we first optimized laboratory protocols (polymerase choice, DNA shearing size, PCR cycling) to enable 76% longer reads (now median length of 6,965 bp) from 100-fold less input DNA (now 1 nanogram). Using a virome from a natural seawater sample, we compared viromes generated with VirION2 against other library preparation options (unamplified, original VirION, and short-read), and optimized downstream informatics for improved long-read error correction and assembly. VirION2 assemblies combined with short-read based data (‘enhanced’ viromes), provided significant improvements over VirION libraries in the recovery of longer and more complete viral genomes, and our optimized error-correction strategy using long- and short-read data achieved 99.97% accuracy. In the seawater virome, VirION2 assemblies captured 5,161 viral populations (including all of the virus populations observed in the other assemblies), 30% of which were uniquely assembled through inclusion of long-reads, and 22% of the top 10% most abundant virus populations derived from assembly of long-reads. Viral populations unique to VirION2 assemblies had significantly higher microdiversity means, which may explain why short-read virome approaches failed to capture them. These findings suggest the VirION2 sample prep and workflow can help researchers better investigate the virosphere, even from challenging low-biomass samples. Our new protocols are available to the research community on protocols.io as a ‘living document’ to facilitate dissemination of updates to keep pace with the rapid evolution of long-read sequencing technology.

Long-reads↗
DOE OSTI2025

Phylogenomic discovery and engineering of nitrogen fixation into the bioenergy woody crop poplar

Biological nitrogen fixation (BNF) is a key process enabling plants in specific lineages to convert atmospheric dinitrogen (N₂) into bioavailable ammonia through symbioses with diazotrophic microbes. Expanding this capability beyond native nitrogen-fixing clades into non-nodulating crops would reduce synthetic fertilizer use, lowering energy inputs and environmental impacts in agriculture. Supported by DOE Funding Award DE-SC0018247, the NitFix project advanced foundational knowledge required to engineer root-nodule symbioses in new host species. The team generated the most comprehensive phylogenomic analysis to date of all known nodulating lineages, resolving the evolutionary history of nitrogen-fixing symbiosis and identifying core gene suites retained across nodulating taxa. Through multimodal genomics, transcriptomics, and functional analyses in Medicago truncatula and related species, the project mapped regulatory networks underlying nodule organogenesis, bacterial infection, and nitrogen-fixation efficiency. Key discoveries include the identification of conserved signaling modules for rhizobial recognition, transcription factors controlling nodule differentiation, and metabolic pathways integrating fixed nitrogen into plant growth. The project also developed enabling tools—including optimized transformation pipelines, gene-editing workflows, and imaging-based phenotyping—to accelerate engineering efforts in emerging models. Together, these results refine the mechanistic framework of symbiotic nitrogen fixation and highlight transferable components essential for rewiring these traits into non-nodulating crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Comparative Analysis of TRGBs (CATs) from Unsupervised, Multi-halo-field Measurements: Contrast is Key

The tip of the red giant branch (TRGB) is an apparent discontinuity of the luminosity function (LF) due to the end of the red giant evolutionary phase and is used to measure distances in the local universe. In practice, tip localization via edge detection response (EDR) relies on several methods applied on a case-by-case basis. It is hard to evaluate how individual choices affect a distance estimation using only a single host field while also avoiding confirmation bias. To devise a standardized approach, we compare unsupervised, algorithmic analyses of the TRGB in multiple halo fields per galaxy. We first optimize methods for the lowest field-to-field dispersion, including spatial filtering, smoothing, and weighting of LF, color band selection, and tip selection based on the number of likely RGB stars and the ratio of stars below versus above the tip (R). We find R, which we call the tip contrast, to be the most important indicator of the quality of EDR measurements; higher R selection can decrease field-to-field dispersion. Further, since R is found to correlate with the age or metallicity of the stellar population based on theoretical modeling, it might result in a displacement of the detected tip magnitude. We find a tip-contrast relation with a slope of -0.023 ± 0.0046 mag/ratio, an ~5σ result that can be used to correct these variations in the detections. When using TRGB to establish a distance ladder, consistent TRGB standardization using tip-contrast relation across rungs is vital to make robust cosmological measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

An Evolutionary Algorithm Approach to Multi-Pass ERL Optics Design

An Energy Recovery Experiment at CEBAF (ER@CEBAF) is aimed at demonstrating high energy, low current, multi-pass energy recovery at the existing 12 GeV CEBAF accelerator. The beam break-up instability, limiting the maximum beam current, can be controlled through minimizing beta functions for the lowest energy pass, which gives a preference to strongly focusing optics, e.g. a semi-periodic FODO lattice. On the other hand, one needs to limit beta function excursions, caused by under focusing, at the higher energy passes, which in turn favors weakly focusing linac optics. Balancing both effects is the main objective of proposed multi-pass linac optics optimization. Here, we discuss an optics design process for ER@CEBAF transverse optics using a genetic algorithm.

Neththikumara, I.↗
DOE OSTI2024

Heterologous expression of Arabidopsis laccase2, laccase4 and peroxidase52 driven under developing xylem specific promoter DX15 improves saccharification in populus

Secondary cell wall holds considerable potential as it has gained immense momentum to replace the lignocellulosic feedstock into fuels. Lignin one of the components of secondary cell wall tightly holds the polysaccharides thereby enhancing the recalcitrance and complexity in the biomass. Laccases (LAC) and peroxidases (PRX) are the major phenyl-oxidases playing key functions during the polymerization of monolignols into lignin. Yet, the functions of laccase and peroxidases gene families remained largely unknown. Hence, the objective of this conducted study is to understand the role of specific LAC and PRX in Populus wood formation and to further investigate how the altered Lac and Prx expression affects biomass recalcitrance and plant growth. This study of heterologous expression of Arabidopsis Lac and Prx genes was conducted in poplar to avoid any otherwise occurring co-suppression mechanism during the homologous overexpression of highly expressed native genes. In the pursuit of optimizing lignocellulosic biomass for biofuel production, the present study focuses on harnessing the enzymatic potential of Arabidopsis thaliana Laccase2, Laccase4, and Peroxidase52 through heterologous expression. We overexpressed selected Arabidopsis laccase2 (AtLac2), laccase4 (AtLac4), and peroxidase52 (AtPrx52) genes, based on their high transcript expression respective to the differentiating xylem tissues in the stem, in hybrid poplar (cv. 717) expressed under the developing xylem tissue-specific promoter, DX15 characterized the transgenic populus for the investigation of growth phenotypes and recalcitrance efficiency. Bioinformatics analyses conducted on AtLac2 and AtLac4 and AtPrx52, revealed the evolutionary relationship between the laccase gene and peroxidase gene homologs, respectively. Transgenic poplar plant lines overexpressing the AtLac2 gene (AtLac2-OE) showed an increase in plant height without a change in biomass yield as compared to the controls; whereas AtLac4-OE and AtPrx52-OE transgenic lines did not show any such observable growth phenotypes compared to their respective controls. The changes in the levels of lignin content and S/G ratios in the transgenic poplar resulted in a significant increase in the saccharification efficiency as compared to the control plants. Overall, saccharification efficiency was increased by 35–50%, 21–42%, and 8–39% in AtLac2-OE, AtLac4-OE, and AtPrx52-OE transgenic poplar lines, respectively, as compared to their controls. Moreover, the bioengineered plants maintained normal growth and development, underscoring the feasibility of this approach for biomass improvement without compromising overall plant fitness. This study also sheds light on the potential of exploiting regulatory elements of DX15 to drive targeted expression of lignin-modifying enzymes, thereby providing a promising avenue for tailoring biomass for improved biofuel production. These findings contribute to the growing body of knowledge in synthetic biology and plant biotechnology, offering a sustainable solution to address the challenges associated with lignocellulosic biomass recalcitrance.

09 BIOMASS FUELS↗
DOE OSTI2025

Warming is Associated With More Encoded Antimicrobial Resistance Genes and Transcriptions Within Five Drug Classes in Soil Bacteria: A Case Study and Synthesis

ABSTRACT The effect of warming on anti‐microbial resistance (AMR) genes in the environment has critical implications for public health but is little studied. We collected published soil bacterial genomes from the BV‐BRC database and tested the correlation between reported optimal growth temperature and the number of encoded AMR genes. Furthermore, we tested the relationship between temperature and AMR gene transcription in a natural ecosystem by analysing soil transcriptomes from a warming manipulation experiment in an Alaskan boreal forest. We hypothesised that there is a positive relationship between warming and AMR prevalence in gene content in bacterial genomes and transcriptomic sequences, and that this effect would vary by drug class. Regarding the bacterial genomes, we found a positive relationship between the fraction of encoded AMR genes and the reported optimal temperature of soil bacteria. The drug classes tetracycline and lincosamide/macrolide/streptogramin had the strongest positive relationship with reported optimal temperature. For the case study in a natural ecosystem, we found 61 significantly upregulated AMR gene‐associated transcripts spanning eight drug classes in warmed plots. In the Alaskan soil samples, we found that warming elicited the strongest positive effect on transcripts targeting lincosamide/streptogramin, beta‐lactam and phenicol/quinolone antibiotics. Overall, higher temperatures were linked to AMR gene prevalence.

Hacopian, Melanie T. [Department of Ecology and Ev↗
DOE OSTI2025

GenomeFace v1.0

GenomeFace is meta-genome binning software. Metagenomic binning, the process of grouping DNA sequences into taxonomic units, is critical for understanding the functions, interactions, and evolutionary dynamics of microbial communities. We propose a deep learning approach to binning using two neural networks, one based on composition and another on environmental abundance, dynamically weighting the contribution of each based on characteristics of the input data. Trained on over 43,000 prokaryotic genomes, our network for composition-based binning is inspired by metric learning techniques used for facial recognition. Using a task-specific, multi-GPU accelerated algorithm to cluster the embeddings produced by our network, our binner leverages marker genes observed to be universally present in nearly all taxa to grade and select optimal clusters of sequences from a hierarchy of candidates. We evaluate our approach on four simulated datasets with known ground truth. Our linear time integration of marker genes recovers more near complete genomes than state of the art but computationally infeasible solutions using them, while being over an order of magnitude faster. Finally, we demonstrate the scalability and acuity of our approach by testing it on three of the largest metagenome assemblies ever performed. Compared to other binners, we produced 47%-183% more near complete genomes. From these datasets, we find over the genomes of over 3000 new candidate species which have never been previously cataloged, representing a potential 4% expansion of the known bacterial tree of life.

Lettich, Richard [Lawrence Berkeley National Labor↗
DOE OSTI2021

The PshX subunit of the photochemical reaction center from Heliobacterium modesticaldum acts as a low-energy antenna

The anoxygenic phototrophic bacterium Heliobacterium modesticaldum contains a photochemical reaction center protein complex (called the HbRC) consisting of a homodimer of the PshA polypeptide and two copies of a newly discovered polypeptide called PshX, which is a single transmembrane helix that binds two bacteriochlorophyll g molecules. Here, to assess the function of PshX, we produced a ΔpshX strain of Hbt. modesticaldum by leveraging the endogenous Hbt. modesticaldum Type I-A CRISPR-Cas system to aid in mutant selection. We optimized this system by separating the homologous recombination and CRISPR-based selection steps into two plasmid transformations, allowing for markerless gene replacement. Fluorescence and low-temperature absorbance of the purified HbRC from the wild-type and ΔpshX strains showed that the bacteriochlorophylls bound by PshX have the lowest site energies in the entire HbRC. This indicates that PshX acts as a low-energy antenna subunit, participating in entropy-assisted uphill energy transfer towards the P 800 special bacteriochlorophyll g pair. We further discuss the role that PshX may play in stability of the HbRC, its conservation in other heliobacterial species, and the evolutionary pressure to produce and maintain single-TMH subunits in similar locations in other reaction centers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Prolonging genetic circuit stability through adaptive evolution of overlapping genes

Abstract The development of synthetic biological circuits that maintain functionality over application-relevant time scales remains a significant challenge. Here, we employed synthetic overlapping sequences in which one gene is encoded or ‘entangled’ entirely within an alternative reading frame of another gene. In this design, the toxin-encoding relE was entangled within ilvA, which encodes threonine deaminase, an enzyme essential for isoleucine biosynthesis. A functional entanglement construct was obtained upon modification of the ribosome-binding site of the internal relE gene. Using this optimized design, we found that the selection pressure to maintain functional IlvA stabilized the production of burdensome RelE for >130 generations, which compares favorably with the most stable kill-switch circuits developed to date. This stabilizing effect was achieved through a complete alteration of the allowable landscape of mutations such that mutations inactivating the entangled genes were disfavored. Instead, the majority of lineages accumulated mutations within the regulatory region of ilvA. By reducing baseline relE expression, these more ‘benign’ mutations lowered circuit burden, which suppressed the accumulation of relE-inactivating mutations, thereby prolonging kill-switch function. Overall, this work demonstrates the utility of sequence entanglement paired with an adaptive laboratory evolution campaign to increase the evolutionary stability of burdensome synthetic circuits.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression. Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to previous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding. Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

59 BASIC BIOLOGICAL SCIENCES↗