Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homology modelling”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

The Transitional Protoplanetary Disk Frequency as a Function of Age: Disk Evolution in the Coronet Cluster, Taurus, and Other 1--8 Myr-old Regions

We present Spitzer 3.6-24 micron photometry and spectroscopy for stars in the 1-3 Myr-old Coronet Cluster, expanding upon the survey of Sicilia-Aguilar et al. (2008). Using sophisticated radiative transfer models, we analyze these new data and those from Sicilia-Aguilar et al. (2008) to identify disks with evidence for substantial dust evolution consistent with disk clearing: transitional disks. We then analyze data in Taurus and others young clusters - IC 348, NGC 2362, and eta Cha -- to constrain the transitional disk frequency as a function of time. Our analysis confirms previous results finding evidence for two types of transitional disks -- those with inner holes and those that are homologously depleted. The percentage of disks in the transitional phase increases from approx.15-20% at 1-2 Myr to > 50% at 5-8 Myr; the mean transitional disk lifetime is closer to approx. 1 Myr than 0.1-0.5 Myr, consistent with previous studies by Currie et al. (2009) and Sicilia-Aguilar et al. (2009). In the Coronet Cluster and IC 348, transitional disks are more numerous for very low-mass M3--M6 stars than for more massive K5-M2 stars, while Taurus lacks a strong spectral type-dependent frequency. Assuming standard values for the gas-to-dust ratio and other disk properties, the lower limit for the masses of optically-thick primordial disks is Mdisk approx. 0.001-0.003 M*. We find that single color-color diagrams do not by themselves uniquely identify transitional disks or primordial disks. Full SED modeling is required to accurately assess disk evolution for individual sources and inform statistical estimates of the transitional disk population in large samples using mid-IR colors.

Currie, Thayne↗
DOE OSTI2021

DNA asymmetry promotes SUMO modification of the single-stranded DNA-binding protein RPA

Repair of DNA double-stranded breaks by homologous recombination (HR) is dependent on DNA end resection and on post-translational modification of repair factors. In budding yeast, single-stranded DNA is coated by replication protein A (RPA) following DNA end resection, and DNA–RPA complexes are then SUMO-modified by the E3 ligase Siz2 to promote repair. Here, we show using enzymatic assays that DNA duplexes containing 3' single-stranded DNA overhangs increase the rate of RPA SUMO modification by Siz2. The SAP domain of Siz2 binds DNA duplexes and makes a key contribution to this process as highlighted by models and a crystal structure of Siz2 and by assays performed using protein mutants. Enzymatic assays performed using DNA that can accommodate multiple RPA proteins suggest a model in which the SUMO-RPA signal is amplified by successive rounds of Siz2-dependent SUMO modification of RPA and dissociation of SUMO-RPA at the junction between single- and double-stranded DNA. Our results provide insights on how DNA architecture scaffolds a substrate and E3 ligase to promote SUMO modification in the context of DNA repair.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Natural transformation as a tool in Acinetobacter baylyi : Streamlined engineering and mutational analysis

Natural transformation and homologous recombination in a soil bacterium, Acinetobacter baylyi ADP1, occur with exceptionally high efficiency. These genetic features can be harnessed to address a wide variety of fundamental and applied scientific topics. Recent advances in synthetic biology and laboratory evolution have led to renewed appreciation for the use of A. baylyi as a model organism. To complement several review articles that highlight new tool sets, this chapter focuses on simple protocols and examples of transformation assays that facilitate genetic analysis and engineering. Whole genome sequencing often reveals extensive genetic variation among closely related isolates that can confound the association of genotypic and phenotypic changes. In A. baylyi, such associations can be deciphered in unique ways by directly transforming cells with linear DNA fragments. The resulting allelic replacement, which occurs at high frequency, rapidly generates desired mutants via targeted chromosomal editing. Diverse screening and selection methods can be used to test hypotheses and streamline experimental strategies to reveal the significance of specific DNA sequences. Moreover, large procedural variations are well tolerated, and techniques can be readily adapted for new purposes. Furthermore, one goal of highlighting natural transformation methodology in A. baylyi is to expand the community of researchers using this versatile bacterial host.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1997

Identification of the prooxidant site of human ceruloplasmin: a model for oxidative damage by copper bound to protein surfaces

Free transition metal ions oxidize lipids and lipoproteins in vitro; however, recent evidence suggests that free metal ion-independent mechanisms are more likely in vivo. We have shown previously that human ceruloplasmin (Cp), a serum protein containing seven Cu atoms, induces low density lipoprotein oxidation in vitro and that the activity depends on the presence of a single, chelatable Cu atom. We here use biochemical and molecular approaches to determine the site responsible for Cp prooxidant activity. Experiments with the His-specific reagent diethylpyrocarbonate (DEPC) showed that one or more His residues was specifically required. Quantitative [14C]DEPC binding studies indicated the importance of a single His residue because only one was exposed upon removal of the prooxidant Cu. Plasmin digestion of [14C]DEPC-treated Cp (and N-terminal sequence analysis of the fragments) showed that the critical His was in a 17-kDa region containing four His residues in the second major sequence homology domain of Cp. A full length human Cp cDNA was modified by site-directed mutagenesis to give His-to-Ala substitutions at each of the four positions and was transfected into COS-7 cells, and low density lipoprotein oxidation was measured. The prooxidant site was localized to a region containing His426 because CpH426A almost completely lacked prooxidant activity whereas the other mutants expressed normal activity. These observations support the hypothesis that Cu bound at specific sites on protein surfaces can cause oxidative damage to macromolecules in their environment. Cp may serve as a model protein for understanding mechanisms of oxidant damage by copper-containing (or -binding) proteins such as Cu, Zn superoxide dismutase, and amyloid precursor protein.

NASA Discipline Regulatory Physiology↗
DOE OSTI2024

Characterization of lignin-degrading enzyme PmdC, which catalyzes a key step in the synthesis of polymer precursor 2-pyrone-4,6-dicarboxylic acid

Pyrone-2,4-dicarboxylic acid (PDC) is a valuable polymer precursor that can be derived from the microbial degradation of lignin. The key enzyme in the microbial production of PDC is 4-carboxy-2-hydroxymuconate-6-semialdehyde (CHMS) dehydrogenase, which acts on the substrate CHMS. We present the crystal structure of CHMS dehydrogenase (PmdC from Comamonas testosteroni) bound to the cofactor NADP, shedding light on its three-dimensional architecture, and revealing residues responsible for binding NADP. Using a combination of structural homology, molecular docking, and quantum chemistry calculations, we have predicted the binding site of CHMS. Key histidine residues in a conserved sequence are identified as crucial for binding the hydroxyl group of CHMS and facilitating dehydrogenation with NADP. Mutating these histidine residues results in a loss of enzyme activity, leading to a proposed model for the enzyme's mechanism. These findings are expected to help guide efforts in protein and metabolic engineering to enhance PDC yields in biological routes to polymer feedstock synthesis.

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

Computational Modeling To Adapt Neutralizing Antibody

Monoclonal antibodies (mAbs) is the leading therapy for viral infections because it provides immediate protection and can be administered at higher levels than in a natural immune response. Finding mAbs that neutralize a broad spectrum of viral targets has proven difficult because many species and strains exist and blanket targeting is a slow and laborious process to experimentally screen 108 variants. A new method is needed to rapidly redesign mAbs for homologous targets. This project speeds up redesign using structure-based computational design to reduce the mAbs search space to a manageable level and screen mutants at a much higher rate than in experiments. Computation will also provide critical knowledge about the fundamental interactions. The project will adapt S230, a human antibody that neutralizes SARS-CoV, to neutralize SARS-COV-2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Cell-Free Protein Synthesis for High-Throughput Biosynthetic Pathway Prototyping

Biological systems provide a sustainable and complimentary approach to synthesizing useful chemical products. Metabolic engineers seeking to establish economically viable biosynthesis platforms strive to increase product titers, rates, and yields. Despite continued advances in genetic tools and metabolic engineering techniques, cellular workflows remain limited in throughput. It may take months to test dozens of unique pathway designs even in a robust model organism, such as Escherichia coli. In contrast, cell-free protein synthesis enables the rapid generation of enzyme libraries that can be combined to reconstitute metabolic pathways in vitro for biochemical synthesis in days rather than weeks. Cell-free reactions thereby enable comparison of hundreds to thousands of unique combinations of enzyme homologs and concentrations, which can quickly identify the most productive pathway variants to test in vivo or further characterize in vitro. Furthermore, this cell-free pathway prototyping strategy provides a complementary approach to accelerate cellular metabolic engineering efforts toward highly productive strains for metabolite production.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A positive feedback loop involving the Spa2 SHD domain contributes to focal polarization

Focal polarization is necessary for finely arranged cell-cell interactions. The yeast mating projection, with its punctate polarisome, is a good model system for this process. We explored the critical role of the polarisome scaffold protein Spa2 during yeast mating with a hypothesis motivated by mathematical modeling and tested by in vivo experiments. Our simulations predicted that two positive feedback loops generate focal polarization, including a novel feedback pathway involving the N-terminal domain of Spa2. We characterized the latter using loss-of-function and gain-of-function mutants. The N-terminal region contains a Spa2 Homology Domain (SHD) which is conserved from yeast to humans, and when mutated largely reproduced the spa2Δ phenotype. Our work together with published data show that the SHD domain recruits Msb3/4 that stimulates Sec4-mediated transport of Bud6 to the polarisome. There, Bud6 activates Bni1-catalyzed actin cable formation, recruiting more Spa2 and completing the positive feedback loop. We demonstrate that disrupting this loop at any point results in morphological defects. Gain-of-function perturbations partially restored focal polarization in a spa2 loss-of-function mutant without restoring localization of upstream components, thus supporting the pathway order. Thus, we have collected data consistent with a novel positive feedback loop that contributes to focal polarization during pheromone-induced polarization in yeast.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2003

Identification and Characterization of Extremophile Microorganisms with Significance to Astrobiology

It is now well recognized that microorganisms thrive in extreme ecological conditions such as geothermal vents, polar region, acid and alkaline lakes, and the cold pressurized depth of the ocean floor of this planet. Morphological, physiological, biochemical and genetic adaptations to extreme environments by these extremophile microorganisms have generated immense interest amongst astrobiologists who increasingly believe in the existence of extraterrestrial life. The evidence collected by NASA's space probe Galileo suggested the presence of liquid water and volcanic activity on Mars and Jupiter's satellite Europa. Volcanic activity provides some of the heat necessary to keep the water on Europa from freezing that could provide important dissolved chemicals needed by living organisms. The possibility of the existence of hypersaline alkaline lakes and evaporites confined within closed volcanic basins and impact craters on Mars, and a layer of liquid water under the ice on Europa provide sufficient 'raison d'etre' to study microorganisms in similar extreme environments on Earth, which could provide us with a model that would help establish the existence of extraterrestrial life on other planetary bodies. The objectives of the summer research project were as follows: (1) application of molecular approaches to help establish new species of extremophile microorganisms isolated from a hypersaline alkaline lake; and (2) identification of a major cold-shock gene (cspA) homolog from a psychrotolerant microorganism, PmagG1.

Bej, Asim K.↗
DOE OSTI2025

Nonthermal Ionization of Kilonova Ejecta: Observable Impacts

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass (≲0.05 M ⊙ ) and high ejecta velocity (≳0.1c). We show that, even at very early times (≲2 days), regions of ejecta fall below critical density and temperature thresholds at which nonlocal thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized r-process elements. We find that incorporating ionization from high-energy electrons produces an “inverted” and “blended” ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta, and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light-curve properties of AT 2017gfo, such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity (∼0.3 c) ejecta components, our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

High energy astrophysics↗
DOE OSTI2021

BETO 2021 Peer Review - Bio-Insecticides from Thermochemical Biomass Conversion

This work is developing a sustainable and inexpensive bioinsecticide that leverages the inherent chemical functionality created from thermochemical biomass conversion. The impact of this project lies at the intersection of energy and agricultural sustainability. Thermochemical conversion of biomass to fuels remains an attractive pathway, and bioinsecticides, isolated from a fraction of upgraded bio-oils, are a high-valued coproduct that can increase biorefinery profitability. Existing insecticide products are facing significant pressure from regulators and consumers due to negative health and environmental impacts. Bioinsecticides, produced through the deconstruction of biomass, can offer a safer more environmentally benign alternative due to the chemical homology with biologically degraded lignin. To enable thermochemical coproduction of bioinsecticides insecticidal activity must be competitive with existing products, separations of the bioinsecticide must be technically viable, and the overall process economics must be improved. Using vacuum distillation, we have separated a bioinsecticide active ingredient that has similar activity to current commercial products, and technoeconomic modeling has shown it can be produced at market competitive prices. This work has established a bioinsecticide fraction that is ready for further, more resource intensive development by addressing early technical, economic, and toxicological risks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Three-dimensional structure of Schistosoma japonicum glutathione S-transferase fused with a six-amino acid conserved neutralizing epitope of gp41 from HIV

The 3-dimensional crystal structure of glutathione S-transferase (GST) of Schistosoma japonicum (Sj) fused with a conserved neutralizing epitope on gp41 (glycoprotein, 41 kDa) of human immunodeficiency virus type 1 (HIV-1) (Muster T et al., 1993, J Virol 67:6642-6647) was determined at 2.5 A resolution. The structure of the 3-3 isozyme rat GST of the mu gene class (Ji X, Zhang P, Armstrong RN, Gilliland GL, 1992, Biochemistry 31:10169-10184) was used as a molecular replacement model. The structure consists of a 4-stranded beta-sheet and 3 alpha-helices in domain 1 and 5 alpha-helices in domain 2. The space group of the Sj GST crystal is P4(3)2(1)2, with unit cell dimensions of a = b = 94.7 A, and c = 58.1 A. The crystal has 1 GST monomer per asymmetric unit, and 2 monomers that form an active dimer are related by crystallographic 2-fold symmetry. In the binding site, the ordered structure of reduced glutathione is observed. The gp41 peptide (Glu-Leu-Asp-Lys-Trp-Ala) fused to the C-terminus of Sj GST forms a loop stabilized by symmetry-related GSTs. The Sj GST structure is compared with previously determined GST structures of mammalian gene classes mu, alpha, and pi. Conserved amino acid residues among the 4 GSTs that are important for hydrophobic and hydrophilic interactions for dimer association and glutathione binding are discussed.

Schistosoma japonicum/enzymology/genetics/immunolo↗
DOE OSTI2022

Neutralizing epitopes on Clostridioides difficile toxin A revealed by the structures of two camelid VHH antibodies

Toxin A (TcdA) and toxin B (TcdB) are two key virulence factors secreted by Clostridioides difficile, which is listed as an urgent threat by the CDC. These two large homologous exotoxins are mainly responsible for diseases associated with C. difficile infection (CDI) with symptoms ranging from diarrhea to life threatening pseudomembranous colitis. Single-domain camelid antibodies (VHHs) AH3 and AA6 are two potent antitoxins against TcdA, which when combined with two TcdB-targeting VHHs showed effective protection against both primary and recurrent CDI in animal models. Here, we report the co-crystal structures of AH3 and AA6 when they form complexes with the glucosyltransferase domain (GTD) and a fragment of the delivery and receptor-binding domain (DRBD) of TcdA, respectively. Based on these structures, we find that AH3 binding enhances the overall stability of the GTD and interferes with its unfolding at acidic pH, and AA6 may inhibit the pH-dependent conformational changes in the DRBD that is necessary for pore formation of TcdA. These studies reveal two functionally critical epitopes on TcdA and shed new insights into neutralizing mechanisms and potential development of epitope-focused vaccines against TcdA.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2010

Interactions of Ku70/80 with Double-Strand DNA: Energetic, Dynamics, and Functional Implications

Space radiation is a proficient inducer of DNA damage leading to mutation, aberrant cell signaling, and cancer formation. Ku is among the first responding proteins in nucleus to recognize and bind the DNA double strand breaks (DSBs) whenever they are introduced. Once loaded Ku works as a scaffold to recruit other repair factors of non-homologous end joining and facilitates the following repair processes. The crystallographic study of the Ku70/80 heterodimer indicate the core structure of this protein shows virtually no conformational change after binding with DNA. To investigate the dynamical features as well as the energetic characteristics of Ku-DNA binding, we conduct multi-nanosecond molecular dynamics simulations of a modeled Ku70/80 structure and several complexes with two 24-bp DNA duplexes. Free energy calculations show significant energy differences between the complexes with Ku bound at DSBs and those with Ku associated at an internal site of a chromosome. The results also reveal detailed interactions between different nucleotides and the amino acids along the DNA-binding cradle of Ku, indicating subtle binding preference of Ku at specific DNA sequences. The covariance matrix analyses along the trajectories demonstrate the protein is stimulated to undergo correlated motions of different domains once bound to DNA ends. Additionally, principle component analyses identify these low frequency collective motions suitable for binding with and translocation along duplex DNA. It is proposed that the modification of dynamical properties of Ku upon binding with DSBs may provide a signal for the further recruitment of other repair factors such as DNA-PKcs, XLF, and XRCC4.

Hu, Shaowen↗
DOE OSTI2025

Accessible, uniform protein property prediction with a scikit-learn based toolset AIDE

Summary Protein property prediction via machine learning with and without labeled data is becoming increasingly powerful, yet methods are disparate and capabilities vary widely over applications. The software presented here, “Artificial Intelligence Driven protein Estimation (AIDE)”, enables instantiating, optimizing, and testing many zero-shot and supervised property prediction methods for variants and variable length homologs in a single, reproducible notebook or script by defining a modular, standardized application programming interface (API), i.e. drop-in compatible with scikit-learn transformers and pipelines. Availability and implementation AIDE is an installable, importable python package inheriting from scikit-learn classes and API and is installable on Windows, Mac, and Linux. Many of the wrapped models internal to AIDE will be effectively inaccessible without a GPU, and some assume CUDA. The newest stable, tested version can be found at https://github.com/beckham-lab/aide_predict and a full user guide and API reference can be found at https://beckham-lab.github.io/aide_predict/. Static versions of both at the time of writing can be found on Zenodo.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Growth-associated polyhydroxybutyrate accumulation in Azospira suillum PS during aerobic and perchlorate respiration

Polyhydroxyalkanoates (PHAs) are widespread microbial storage polymers increasingly recognized for roles beyond carbon and energy storage, including redox homeostasis and stress physiology. While PHA accumulation is classically associated with stationary-phase metabolism under severe nutrient imbalance, comparatively little is known about growth-associated PHA synthesis in facultative anaerobes with unusual respiratory strategies. Here, we investigated polyhydroxybutyrate (PHB) metabolism in Azospira suillum PS, a genetically tractable perchlorate-reducing bacterium capable of both aerobic and anaerobic respiration. Using physiological growth experiments, targeted gene deletions, PHB extractions, and intracellular redox measurements, we examined PHB accumulation under varying respiratory and nutrient conditions. We demonstrate that PHB accumulates during exponential growth under moderately nitrogen-limited conditions, both aerobically and during perchlorate respiration, representing a rare example of growth-associated PHB synthesis under anaerobic conditions. Genomic analysis revealed four phaC homologs, one of which could not be deleted under the experimental conditions tested and co-localized with phaB and phaR. Redox profiling further revealed a strong positive correlation between PHB accumulation and intracellular NADPH/NADP+ ratios. Together, these findings expand the physiological contexts in which PHB synthesis is known to occur and highlight perchlorate-respiring bacteria as underexplored model systems for studying growth-integrated carbon and redox storage strategies.

Meier, David A O↗
NASA NTRS2006

Geochemistry of Pallasite Olivines and the Origin of Main-Group Pallasites

Main-group pallasites (PMG) are mixtures of iron-nickel metal and magnesian olivine thought to have been formed at the core-mantle boundary of an asteroid [1]. Some have anomalous metal compositions (PMG-am) and a few have atypically ferroan olivines (PMG-as) [2]. PMG metal is consistent with an origin as a late fractionate of the IIIAB iron core [2]. Most PMG olivines have very similar Fe/Mg ratios, likely due to subsolidus redox reaction with the metal [3]. In contrast, minor and trace elements show substantial variation, which may be explained by either: (i) PMG were formed at a range of depths in the parent asteroid; the element variations reflect variations in igneous evolution with depth, (ii) the pallasite parent asteroid was chemically heterogeneous; the heterogeneity partially survived igneous processing, or (iii) PMG represent the core-mantle boundaries of several distinct parent asteroids [4, 5]. We have continued doing major, minor and trace elements by EMPA and INAA on a wider suite of PMG olivines, and have begun doing precise oxygen isotope analyses to test these hypotheses. Manganese is homologous with Fe(2+), and can be used to distinguish between magmatic and redox processes as causes for Fe/Mg variations. PMG olivines have a range in molar 1000*Mn/Mg of 2.3-4.6 indicating substantial igneous fractionation in olivines with very similar Fe/Mg (0.138-0.148). The Mg-Mn-Fe distributions can be explained by a fractional crystallization-reduction model; higher Mn/Mg ratios reflect more evolved olivines while Fe/Mg is buffered by redox reactions with the metal. There is a positive association between Mn/Mg and Sc content that is consistent with igneous fractionation. However, most PMG olivines fall within a narrow Mn/Mg range (3.0-3.6), but these show a substantial range in Sc (1.00-2.29 micro-g/g). Assuming fractional crystallization, this Sc range could have resulted from approx.65% crystallization of an ultramafic magma. This is inconsistent with formation at the core-mantle boundary of a single asteroid [4]. One alternative is that the PMG are fragments of several asteroids, and these could have had different initial Sc contents, Mn/Mg and differences in igneous history. Our preliminary O isotope data and those of [6, 7] do not support this, although the coverage of PMG olivines is incomplete. The PMG-as Springwater is not easily fit in any scenario. Its olivine has among the highest Mn/Mg suggesting it is one of the most evolved, but the lowest Sc content suggesting it is the least evolved. The O isotopic composition of Springwater olivine is the same as that of other PMG. Thus there is no indication that it represents a distinct parent asteroid. Our preliminary O isotopic data favor a single PMG parent asteroid. In this case, the olivines are more likely melt-residues, and that the parent asteroid was initially heterogeneous in chemical, but not isotopic, composition.

Mittlefehldt, D. W.↗
DOE OSTI2022

Gene flow and introgression are pervasive forces shaping the evolution of bacterial species

Abstract Background Although originally thought to evolve clonally, studies have revealed that most bacteria exchange DNA. However, it remains unclear to what extent gene flow shapes the evolution of bacterial genomes and maintains the cohesion of species. Results Here, we analyze the patterns of gene flow within and between >2600 bacterial species. Our results show that fewer than 10% of bacterial species are truly clonal, indicating that purely asexual species are rare in nature. We further demonstrate that the taxonomic criterion of ~95% genome sequence identity routinely used to define bacterial species does not accurately represent a level of divergence that imposes an effective barrier to gene flow across bacterial species. Interruption of gene flow can occur at various sequence identities across lineages, generally from 90 to 98% genome identity. This likely explains why a ~95% genome sequence identity threshold has empirically been judged as a good approximation to define bacterial species. Our results support a universal mechanism where the availability of identical genomic DNA segments required to initiate homologous recombination is the primary determinant of gene flow and species boundaries in bacteria. We show that these barriers of gene flow remain porous since many distinct species maintain some level of gene flow, similar to introgression in sexual organisms. Conclusions Overall, bacterial evolution and speciation are likely shaped by similar forces driving the evolution of sexual organisms. Our findings support a model where the interruption of gene flow—although not necessarily the initial cause of speciation—leads to the establishment of permanent and irreversible species borders.

59 BASIC BIOLOGICAL SCIENCES↗