Search NASA⌕ Search

SEARCH · Search NASA

Results for “Biotechnology & Applied Microbiology”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Microfluidic droplets with amended culture media cultivate a greater diversity of soil microorganisms

ABSTRACT Uncultivated but abundant soil microorganisms have untapped potential for producing broad ranges of natural products, as well as for bioremediation. However, cultivating soil microorganisms while maintaining a broad microorganism diversity to enable phenotyping and functional analysis of as diverse individual isolates as possible remains challenging. In this study, we developed and tested the ability of several culture media formulations that contain defined soil metabolites or soil extracts to maintain microorganism diversity during culture. We also assessed their performance in microfluidic droplet cultivation where single-soil microorganism isolates were encapsulated and cultivated in picoliter-volume water-in-oil emulsion droplets to enable clonal growth needed for downstream functional analyses. Our results show that droplet cultivation with media supplemented by soil extract or soil metabolites enables the recovery of soil microorganisms with higher diversity (up to 1.5-fold higher richness) compared to bulk cultivation methods. Importantly, 1.7-fold more of less abundant (<1%) phyla and 11-fold more of unique genera were recovered, demonstrating the utility of this method for interrogating highly diverse soil microorganisms for broad ranges of applications. IMPORTANCE Although soil microorganisms hold a significant value in bioproduction and bioremediation, only a small fraction—less than 1%—can be cultured under specific media and cultivation conditions. This indicates that there are ample opportunities in harvesting the diverse environmental microorganisms if isolating and recovering these uncultured microorganisms are possible. This paper presents a new cultivation technique composed of isolating single-soil microorganism cell from anin situsoil microorganism community in microfluidic droplets and conducting in-droplet cultivation in media supplemented by soil extract or soil metabolites. This method enables the recovery of a broader diversity of the original microorganism community, laying the groundwork for a high-throughput phenotyping of these diverse microorganisms from their natural habitats.

Biotechnology & Applied Microbiology↗
DOE OSTI2024

Metabolic engineering of Caldicellulosiruptor bescii for 2,3-butanediol production from unpretreated lignocellulosic biomass and metabolic strategies for improving yields and titers

ABSTRACT The platform chemical 2,3-butanediol (2,3-BDO) is used to derive products, such as 1,3-butadiene and methyl ethyl ketone, for the chemical and fuel production industries. Efficient microbial 2,3-BDO production at industrial scales has not been achieved yet for various reasons, including product inhibition to host organisms, mixed stereospecificity in product formation, and dependence on expensive substrates (i.e., glucose). In this study, we explore engineering of a 2,3-BDO pathway inCaldicellulosiruptor bescii, an extremely thermophilic (optimal growth temperature = 78°C) and anaerobic bacterium that can break down crystalline cellulose and hemicellulose into fermentable C 5 and C 6 sugars. In addition, C. besciigrows on unpretreated plant biomass, such as switchgrass. Biosynthesis of 2,3-BDO involves three steps: two molecules of pyruvate are condensed into acetolactate; acetolactate is decarboxylated to acetoin, and finally, acetoin is reduced to 2,3-BDO.C. besciinatively produces acetoin; therefore, in order to complete the 2,3-BDO biosynthetic pathway,C. besciiwas engineered to produce a secondary alcohol dehydrogenase (sADH) to catalyze the final step. Two previously characterized, thermostable sADH enzymes with high affinity for acetoin, one from a bacterium and one from an archaeon, were tested independently. When either sADH was present inC. bescii,the recombinant strains were able to produce up to 2.5-mM 2,3-BDO from crystalline cellulose and xylan and 0.2-mM 2,3-BDO directly from unpretreated switchgrass. This serves as the basis for higher yields and productivities, and to this end, limiting factors and potential genetic targets for further optimization were assessed using the genome-scale metabolic model ofC. bescii. IMPORTANCE Lignocellulosic plant biomass as the substrate for microbial synthesis of 2,3-butanediol is one of the major keys toward cost-effective bio-based production of this chemical at an industrial scale. However, deconstruction of biomass to release the sugars for microbial growth currently requires expensive thermochemical and enzymatic pretreatments. In this study, the thermo-cellulolytic bacteriumCaldicellulosiruptor besciiwas successfully engineered to produce 2,3-butanediol from cellulose, xylan, and directly from unpretreated switchgrass. Genome-scale metabolic modeling ofC. besciiwas applied to adjust carbon and redox fluxes to maximize productivity of 2,3-butanediol, thereby revealing bottlenecks that require genetic modifications.

Biotechnology & Applied Microbiology↗
DOE OSTI2019

Comparative Genomics of Rumen Butyrivibrio spp. Uncovers a Continuum of Polysaccharide-Degrading Capabilities

Plant polysaccharide breakdown by microbes in the rumen is fundamental to digestion in ruminant livestock. Bacterial species belonging to the rumen genera Butyrivibrio and Pseudobutyrivibrio are important degraders and utilizers of lignocellulosic plant material. These bacteria degrade polysaccharides and ferment the released monosaccharides to yield short-chain fatty acids that are used by the ruminant for growth and the production of meat, milk, and fiber products. Although rumen Butyrivibrio and Pseudobutyrivibrio species are regarded as common rumen inhabitants, their polysaccharide-degrading and carbohydrate-utilizing enzymes are not well understood. In this study, we analyzed the genomes of 40 Butyrivibrio and 6 Pseudobutyrivibrio strains isolated from the plant-adherent fraction of New Zealand dairy cows to explore the polysaccharide-degrading potential of these important rumen bacteria. Comparative genome analyses combined with phylogenetic analysis of their 16S rRNA genes and short-chain fatty acid production patterns provide insight into the genomic diversity and physiology of these bacteria and divide Butyrivibrio into 3 species clusters. Rumen Butyrivibrio bacteria were found to encode a large and diverse spectrum of degradative carbohydrate-active enzymes (CAZymes) and binding proteins. In total, 4,421 glycoside hydrolases (GHs), 1,283 carbohydrate esterases (CEs), 110 polysaccharide lyases (PLs), 3,605 glycosyltransferases (GTs), and 1,706 carbohydrate-binding protein modules (CBM) with predicted activities involved in the depolymerization and transport of the insoluble plant polysaccharides were identified. Butyrivibrio genomes had similar patterns of CAZyme families but varied greatly in the number of genes within each category in the Carbohydrate-Active Enzymes database (CAZy), suggesting some level of functional redundancy. These results suggest that rumen Butyrivibrio species occupy similar niches but apply different degradation strategies to be able to coexist in the rumen.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Enhanced transport of bacteria along root systems by protists can impact plant health

ABSTRACT Soil protists have been shown to contribute to the structure and function of the rhizosphere in a variety of ways. Protists are key contributors to nutrient cycling through the microbial loop, where biomass is digested by protists and otherwise stored nutrients are returned to the environment. Protists have also been shown to feed on plant pathogenic bacteria and alter root microbiomes in ways that may benefit plants. Recently, a mechanism involving bacterial transport, facilitated by protists, has been hypothesized to contribute to the spatial distribution of bacteria in the rhizosphere. Here, we observe the differential abilities of three soil protists: a ciliate (Colpodasp.), a flagellate (Cercomonassp.), and a naked amoeba (Acanthamoeba castellanii) to transport nitrogen-fixingSinorhizobium melilotito infectible root tips. Co-inoculation of protists plusS. melilotiresulted in the movement of bacteria, as measured by the presence of nitrogen-fixing nodules, up to 15 cm farther down the root systems when compared to plants inoculated withS. melilotialone. Co-inoculation of the ciliate, Colpodasp., withS. meliloti, resulted in shoot weights that were similar to plants that grew in nitrogen-replete potting mix.Colpodasp.-feeding style and motility likely contributed to their success at transporting bacteria through the rhizosphere. We observed that the addition of protists alone without the co-inoculum ofS. melilotiresulted in plants with larger shoot weights than control plants. Follow-up experiments showed that protists plus their associated microbiomes were aiding in plant health, likely through means of nutrient cycling. IMPORTANCE Protists represent a significant portion of the rhizosphere microbiome and have been shown to contribute to plant health, yet they are understudied compared to their bacterial and fungal counterparts. This study elucidates their role in the rhizosphere community and suggests a mechanism by which protists can be used to move bacteria along plant roots. We found that the co-inoculation of protists with nitrogen-fixing beneficial bacteria,Sinorhizobium meliloti, resulted in nodules farther down the roots when compared to plants inoculated with S.melilotialone, and shoot weights similar to plants that received nitrogen fertilizer. These data illustrate the ability of protists to transport viable bacteria to uninhabited regions of the root system.

Biotechnology & Applied Microbiology↗
DOE OSTI2025

Polyyne production is regulated by the transcriptional regulators PgnC and GacA in Pseudomonas protegens Pf-5

ABSTRACT Polyynes produced by bacteria have promising applications in agriculture and medicine due to their potent antimicrobial activities. Polyyne biosynthetic genes have been identified inPseudomonasandBurkholderia. However, the molecular mechanisms underlying the regulation of polyyne biosynthesis remain largely unknown. In this study, we used a soil bacteriumPseudomonas protegensPf-5, which was recently reported to produce polyyne called protegenin, as a model to investigate the regulation of bacterial polyyne production. Our results show that Pf-5 controls polyyne production at both the pathway-specific level and a higher global level. Mutation ofpgnC, a transcriptional regulatory gene located in the polyyne biosynthetic gene cluster, abolished polyyne production. Gene expression analysis revealed that PgnC directly activates the promoter of polyyne biosynthetic genes. The production of polyyne also requires a global regulator GacA. Mutation ofgacAdecreased the translation of PgnC, which is consistent with the result thatpgnCleader mRNA bound directly to RsmE, an RNA-binding protein negatively regulated by GacA. These results suggest that GacA induces the expression of the PgnC regulator, which in turn activates polyyne biosynthesis. Additionally, the polyyne-producing strain of Pf-5, but not the polyyne-nonproducing strain, could inhibit a broad spectrum of bacteria including both Gram-negative and Gram-positive bacteria. IMPORTANCE Antimicrobial metabolites produced by bacteria are widely used in agriculture and medicine to control plant, animal, and human pathogens. Although bacteria-derived polyynes have been identified as potent antimicrobials for decades, the molecular mechanisms by which bacteria regulate polyyne biosynthesis remain understudied. In this study, we found that polyyne biosynthesis is directly activated by a pathway-specific regulator PgnC, which is induced by a global regulator GacA through the RNA-binding protein RsmE inPseudomonas protegens. To our knowledge, this work is the first comprehensive study of the regulatory mechanisms of bacterial polyyne biosynthesis at both pathway-specific level and global level. The discovered molecular mechanisms can help us optimize polyyne production for agricultural or medical applications.

Biotechnology & Applied Microbiology↗
DOE OSTI2022

Differential Response of Bacterial Microdiversity to Simulated Global Change

Microbial communities play an integral role in global biogeochemical cycling, but our understanding of how global change will affect microbial community structure and functioning remains limited. Microbiome analyses typically aggregate large amounts of genetic diversity which may obscure finer variation in traits.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Multi-technique characterization of iron reduction by an Antarctic Shewanella : an analog system for putative Martian biosignature identification

ABSTRACT Microbes from terrestrial extreme environments enable testing of biosignature production in conditions relevant to astrobiological targets. Mars, which was likely more conducive to life during early warmer and wetter epochs, has inspired missions that search for signs of early life in the surficial rock record, including mineral or organic biosignatures. Microbial iron reduction is a common and ancient metabolism that may have also operated on other rocky celestial bodies. To investigate biosignature production during iron reduction, aShewanellasp. (strain BF02_Schw) isolated from a subglacial discharge known as Blood Falls, Antarctica, was incubated with the electron acceptor ferrihydrite (Fh). Biosignatures associated with Fh reduction were identified using a suite of techniques currently utilized or proposed for Mars missions, including X-ray diffraction and infrared, Mössbauer, and Raman spectroscopy. The biotic origin of features was validated by transcriptional changes observed between treatments with and without Fh and comparison to killed controls. In live treatments, Fh was reduced to magnetite and goethite, both detected in Martian lacustrine basins. Several soluble and volatile metabolites were also detected, including riboflavin and dimethyl sulfide (DMS), which could be astrobiological indicators of active microbial processes. While none of the identified biosignatures individually would serve as definitive proof of life (past or present), detecting concomitant features associated with known terrestrial biotic processes would provide compelling rationale for more targeted life detection missions. Terrestrial extremophiles can support the exploration of astrobiologically relevant microbial processes, validation of life detection instrumentation, and potentially the discovery of new biomarkers. IMPORTANCE Culture-based experiments with terrestrial extremophiles can elucidate biosignatures that may be analogous to those produced under extraterrestrial conditions, and thus inform sampling and technology strategies for future missions. Here, we demonstrate the production of several biosignatures under iron-reducing conditions byShewanellasp. BF02_Schw, originally isolated from an Antarctic analog feature. These biosignatures could be detectable using flight-ready instrumentation. Growth experiments with terrestrial extremophiles can identify biosignatures measurable by current methodologies and inform the development and optimization of techniques for detecting extant or extinct life on other worlds.

Biotechnology & Applied Microbiology↗
DOE OSTI2021

Redundancy in Aromatic O -Demethylation and Ring-Opening Reactions in Novosphingobium aromaticivorans and Their Impact in the Metabolism of Plant-Derived Phenolics

ABSTRACT Lignin is a plant heteropolymer composed of phenolic subunits. Because of its heterogeneity and recalcitrance, the development of efficient methods for its valorization still remains an open challenge. One approach to utilize lignin is its chemical deconstruction into mixtures of monomeric phenolic compounds, followed by biological funneling into a single product. Novosphingobium aromaticivorans DSM 12444 has been previously engineered to produce 2-pyrone-4,6-dicarboxylic acid (PDC) from depolymerized lignin by simultaneously metabolizing multiple aromatics through convergent routes involving the intermediates 3-methoxygallic acid (3-MGA) and protocatechuic acid (PCA). We investigated enzymes predicted to be responsible for O -demethylation and oxidative aromatic ring opening, two critical reactions involved in the metabolism of phenolic compounds by N. aromaticivorans . The results showed the involvement of DesA in O -demethylation of syringic and vanillic acids, LigM in O- demethylation of vanillic acid and 3-MGA, and a new O- demethylase, DmtS, in the conversion of 3-MGA into gallic acid (GA). In addition, we found that LigAB was the main aromatic ring-opening dioxygenase involved in 3-MGA, PCA, and GA metabolism and that a previously uncharacterized dioxygenase, LigAB2, had high activity with GA. Our results indicate a metabolic route not previously identified in N. aromaticivorans that involves O -demethylation of 3-MGA to GA. We predict that this pathway channels ∼15% of the carbon flow from syringic acid, with the rest following ring opening of 3-MGA. The new knowledge obtained in this study allowed for the creation of an improved engineered strain for the funneling of aromatic compounds that exhibits stoichiometric conversion of syringic acid into PDC. IMPORTANCE For lignocellulosic biorefineries to effectively contribute to reduction of fossil fuel use, they need to become efficient at producing chemicals from all major components of plant biomass. Making products from lignin will require engineering microorganisms to funnel multiple phenolic compounds to the chemicals of interest, and N. aromaticivorans is a promising chassis for this technology. The ability of N. aromaticivorans to efficiently and simultaneously degrade many phenolic compounds may be linked to having functionally redundant aromatic degradation pathways and enzymes with broad substrate specificity. A detailed knowledge of aromatic degradation pathways is thus essential to identify genetic engineering targets to maximize product yields. Furthermore, knowledge of enzyme substrate specificity is critical to redirect flow of carbon to desired pathways. This study described an uncharacterized pathway in N. aromaticivorans and the enzymes that participate in this pathway, allowing the engineering of an improved strain for production of PDC from lignin.

Biotechnology & Applied Microbiology↗
DOE OSTI2020

Combining Random Mutagenesis and Metabolic Engineering for Enhanced Tryptophan Production in Synechocystis sp. Strain PCC 6803

Tryptophan (Trp) is an essential aromatic amino acid that has value as an animal feed supplement, as the amount found in plant-based sources is insufficient. An alternative to production by engineered microbial fermentation is to have tryptophan biosynthesized by a photosynthetic microorganism that could replace or supplement both the plant and industrially used microbes. We selected Synechocystis sp. strain PCC 6803, a model cyanobacterium, as the host and studied metabolic engineering and random mutagenesis approaches. Previous work on engineering heterotrophic microbes for improved Trp titers has targeted allosteric feedback regulation in enzymes 3-deoxy-d-arabinoheptulosonate 7-phosphate synthase (DAHPS) and anthranilate synthase (AS) as major bottlenecks in the shikimate pathway. In this work, the genes encoding feedback-resistant enzymes from Escherichia coli, aroG fbr and trpE fbr , were overexpressed in the host wild-type (WT) strain. Separately, the WT strain was subjected to random mutagenesis and selection using an amino acid analog to isolate tryptophan-overproducing strains. The randomly mutagenized strains were sequenced in order to identify the mutations that resulted in the desirable phenotypes. Interestingly, the tryptophan overproducers had mutations in the gene encoding chorismate mutase (CM), which catalyzes the conversion of chorismate to prephenate. The best tryptophan overproducer from random mutagenesis was selected as a host for metabolic engineering where aroG fbr and trpE fbr were overexpressed. The best strain developed produced 212 ± 23 mg/liter of tryptophan after 10 days of photoautotrophic growth under 3% (vol/vol) CO2. We demonstrated that a combination of random mutagenesis and metabolic engineering was superior to either individual approach.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Quantile-Dependent Expressivity of Serum Uric Acid Concentrations

“Quantile-dependent expressivity” occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g., serum uric acid) is high or low relative to its distribution. Analyses were performed to test whether serum uric acid heritability is quantile-specific and whether this could explain some reported gene-environment interactions. Methods. Serum uric acid concentrations were analyzed from 2151 sibships and 12,068 offspring-parent pairs from the Framingham Heart Study. Quantile-specific heritability from offspring-parent regression slopes (β OP , h 2 = 2β OP /(1 + r spouse )) and full-sib regression slopes (β FS ,h 2 = {(1 + 8r spouse β FS ) 0.5 - 1}/(2r spouse ) was robustly estimated by quantile regression with nonparametric significance assigned from 1000 bootstrap samples. Results. Quantile-specific h 2 (±SE) increased with increasing percentiles of the offspring’s sex- and age-adjusted uric acid distribution when estimated from β OP (P trend = 0.001): 0.34 ± 0.03 at the 10th, 0.36 ± 0.03 at the 25th, 0.41 ± 0.03 at the 50th, 0.46 ± 0.04 at the 75th, and 0.49 ± 0.05 at the 90th percentile and when estimated from β FS (P trend = 0.006). This is consistent with the larger genetic effect size of (1) the SLC2A9 rs11722228 polymorphism in gout patients vs. controls, (2) the ABCG2 rs2231142 polymorphism in men vs. women, (3) the SLC2A9 rs13113918 polymorphism in obese patients prior to bariatric surgery vs. two-year postsurgery following 29 kg weight loss, (4) the ABCG2 rs6855911 polymorphism in obese vs. nonobese women, and (5) the LRP2 rs2544390 polymorphism in heavier drinkers vs. abstainers. Quantile-dependent expressivity may also explain the larger genetic effect size of an SLC2A9/PKD2/ABCG2 haplotype for high vs. low intakes of alcohol, chicken, or processed meats. Conclusions. Heritability of serum uric acid concentrations is quantile-specific.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

DeepGRN: prediction of transcription factor binding site across cell-types using attention-based deep neural networks

Abstract Background Due to the complexity of the biological systems, the prediction of the potential DNA binding sites for transcription factors remains a difficult problem in computational biology. Genomic DNA sequences and experimental results from parallel sequencing provide available information about the affinity and accessibility of genome and are commonly used features in binding sites prediction. The attention mechanism in deep learning has shown its capability to learn long-range dependencies from sequential data, such as sentences and voices. Until now, no study has applied this approach in binding site inference from massively parallel sequencing data. The successful applications of attention mechanism in similar input contexts motivate us to build and test new methods that can accurately determine the binding sites of transcription factors. Results In this study, we propose a novel tool (named DeepGRN) for transcription factors binding site prediction based on the combination of two components: single attention module and pairwise attention module. The performance of our methods is evaluated on the ENCODE-DREAM in vivo Transcription Factor Binding Site Prediction Challenge datasets. The results show that DeepGRN achieves higher unified scores in 6 of 13 targets than any of the top four methods in the DREAM challenge. We also demonstrate that the attention weights learned by the model are correlated with potential informative inputs, such as DNase-Seq coverage and motifs, which provide possible explanations for the predictive improvements in DeepGRN. Conclusions DeepGRN can automatically and effectively predict transcription factor binding sites from DNA sequences and DNase-Seq coverage. Furthermore, the visualization techniques we developed for the attention modules help to interpret how critical patterns from different types of input features are recognized by our model.

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

Optimal dimensionality selection for independent component analysis of transcriptomic data

Independent component analysis is an unsupervised machine learning algorithm that separates a set of mixed signals into a set of statistically independent source signals. Applied to high-quality gene expression datasets, independent component analysis effectively reveals both the source signals of the transcriptome as co-regulated gene sets, and the activity levels of the underlying regulators across diverse experimental conditions. Two major variables that affect the final gene sets are the diversity of the expression profiles contained in the underlying data, and the user-defined number of independent components, or dimensionality, to compute. Availability of high-quality transcriptomic datasets has grown exponentially as high-throughput technologies have advanced; however, optimal dimensionality selection remains an open question. We computed independent components across a range of dimensionalities for four gene expression datasets with varying dimensions (both in terms of number of genes and number of samples). We computed the correlation between independent components across different dimensionalities to understand how the overall structure evolves as the number of user-defined components increases. We then measured how well the resulting gene clusters reflected known regulatory mechanisms, and developed a set of metrics to assess the accuracy of the decomposition at a given dimension. We found that over-decomposition results in many independent components dominated by a single gene, whereas under-decomposition results in independent components that poorly capture the known regulatory structure. From these results, we developed a new method, called OptICA, for finding the optimal dimensionality that controls for both over- and under-decomposition. Specifically, OptICA selects the highest dimension that produces a low number of components that are dominated by a single gene. We show that OptICA outperforms two previously proposed methods for selecting the number of independent components across four transcriptomic databases of varying sizes. OptICA avoids both over-decomposition and under-decomposition of transcriptomic datasets resulting in the best representation of the organism’s underlying transcriptional regulatory network.

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

G4Boost: a machine learning-based tool for quadruplex identification and stability prediction

Background: G-quadruplexes (G4s), formed within guanine-rich nucleic acids, are secondary structures involved in important biological processes. Although every G4 motif has the potential to form a stable G4 structure, not every G4 motif would, and accurate energy-based methods are needed to assess their structural stability. Here, we present a decision tree-based prediction tool, G4Boost, to identify G4 motifs and predict their secondary structure folding probability and thermodynamic stability based on their sequences, nucleotide compositions, and estimated structural topologies. Results: G4Boost predicted the quadruplex folding state with an accuracy greater then 93% and an F1-score of 0.96, and the folding energy with an RMSE of 4.28 and R 2 of 0.95 only by the means of sequence intrinsic feature. G4Boost was successfully applied and validated to predict the stability of experimentally-determined G4 structures, including for plants and humans. Conclusion: G4Boost outperformed the three machine-learning based prediction tools, DeepG4, Quadron, and G4RNA Screener, in terms of both accuracy and F1-score, and can be highly useful for G4 prediction to understand gene regulation across species including plants and humans.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Profiling the BLAST bioinformatics application for load balancing on high-performance computing clusters

Abstract Background The Basic Local Alignment Search Tool (BLAST) is a suite of commonly used algorithms for identifying matches between biological sequences. The user supplies a database file and query file of sequences for BLAST to find identical sequences between the two. The typical millions of database and query sequences make BLAST computationally challenging but also well suited for parallelization on high-performance computing clusters. The efficacy of parallelization depends on the data partitioning, where the optimal data partitioning relies on an accurate performance model. In previous studies, a BLAST job was sped up by 27 times by partitioning the database and query among thousands of processor nodes. However, the optimality of the partitioning method was not studied. Unlike BLAST performance models proposed in the literature that usually have problem size and hardware configuration as the only variables, the execution time of a BLAST job is a function of database size, query size, and hardware capability. In this work, the nucleotide BLAST application BLASTN was profiled using three methods: shell-level profiling with the Unix “time” command, code-level profiling with the built-in “profiler” module, and system-level profiling with the Unix “gprof” program. The runtimes were measured for six node types, using six different database files and 15 query files, on a heterogeneous HPC cluster with 500+ nodes. The empirical measurement data were fitted with quadratic functions to develop performance models that were used to guide the data parallelization for BLASTN jobs. Results Profiling results showed that BLASTN contains more than 34,500 different functions, but a single function, RunMTBySplitDB, takes 99.12% of the total runtime. Among its 53 child functions, five core functions were identified to make up 92.12% of the overall BLASTN runtime. Based on the performance models, static load balancing algorithms can be applied to the BLASTN input data to minimize the runtime of the longest job on an HPC cluster. Four test cases being run on homogeneous and heterogeneous clusters were tested. Experiment results showed that the runtime can be reduced by 81% on a homogeneous cluster and by 20% on a heterogeneous cluster by re-distributing the workload. Discussion Optimal data partitioning can improve BLASTN’s overall runtime 5.4-fold in comparison with dividing the database and query into the same number of fragments. The proposed methodology can be used in the other applications in the BLAST+ suite or any other application as long as source code is available.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Pickaxe: a Python library for the prediction of novel metabolic reactions

Abstract Background Biochemical reaction prediction tools leverage enzymatic promiscuity rules to generate reaction networks containing novel compounds and reactions. The resulting reaction networks can be used for multiple applications such as designing novel biosynthetic pathways and annotating untargeted metabolomics data. It is vital for these tools to provide a robust, user-friendly method to generate networks for a given application. However, existing tools lack the flexibility to easily generate networks that are tailor-fit for a user’s application due to lack of exhaustive reaction rules, restriction to pre-computed networks, and difficulty in using the software due to lack of documentation. Results Here we present Pickaxe, an open-source, flexible software that provides a user-friendly method to generate novel reaction networks. This software iteratively applies reaction rules to a set of metabolites to generate novel reactions. Users can select rules from the prepackaged JN1224min ruleset, derived from MetaCyc, or define their own custom rules. Additionally, filters are provided which allow for the pruning of a network on-the-fly based on compound and reaction properties. The filters include chemical similarity to target molecules, metabolomics, thermodynamics, and reaction feasibility filters. Example applications are given to highlight the capabilities of Pickaxe: the expansion of common biological databases with novel reactions, the generation of industrially useful chemicals from a yeast metabolome database, and the annotation of untargeted metabolomics peaks from an E. coli dataset. Conclusion Pickaxe predicts novel metabolic reactions and compounds, which can be used for a variety of applications. This software is open-source and available as part of the MINE Database python package ( https://pypi.org/project/minedatabase/ ) or on GitHub ( https://github.com/tyo-nu/MINE-Database ). Documentation and examples can be found on Read the Docs ( https://mine-database.readthedocs.io/en/latest/ ). Through its documentation, pre-packaged features, and customizable nature, Pickaxe allows users to generate novel reaction networks tailored to their application.

59 BASIC BIOLOGICAL SCIENCES↗