Search NASA⌕ Search

SEARCH · Search NASA

Results for “regulatory network”

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 91 records · Page 5

Coding and noncoding variants in EBF3 are involved in HADDS and simplex autism

Abstract Background Previous research in autism and other neurodevelopmental disorders (NDDs) has indicated an important contribution of protein-coding (coding) de novo variants (DNVs) within specific genes. The role of de novo noncoding variation has been observable as a general increase in genetic burden but has yet to be resolved to individual functional elements. In this study, we assessed whole-genome sequencing data in 2671 families with autism (discovery cohort of 516 families, replication cohort of 2155 families). We focused on DNVs in enhancers with characterized in vivo activity in the brain and identified an excess of DNVs in an enhancer named hs737. Results We adapted the fitDNM statistical model to work in noncoding regions and tested enhancers for excess of DNVs in families with autism. We found only one enhancer (hs737) with nominal significance in the discovery (p = 0.0172), replication (p = 2.5 × 10 −3 ), and combined dataset (p = 1.1 × 10 −4 ). Each individual with a DNV in hs737 had shared phenotypes including being male, intact cognitive function, and hypotonia or motor delay. Our in vitro assessment of the DNVs showed they all reduce enhancer activity in a neuronal cell line. By epigenomic analyses, we found that hs737 is brain-specific and targets the transcription factor gene EBF3 in human fetal brain. EBF3 is genome-wide significant for coding DNVs in NDDs (missense p = 8.12 × 10 −35 , loss-of-function p = 2.26 × 10 −13 ) and is widely expressed in the body. Through characterization of promoters bound by EBF3 in neuronal cells, we saw enrichment for binding to NDD genes (p = 7.43 × 10 −6 , OR = 1.87) involved in gene regulation. Individuals with coding DNVs have greater phenotypic severity (hypotonia, ataxia, and delayed development syndrome [HADDS]) in comparison to individuals with noncoding DNVs that have autism and hypotonia. Conclusions In this study, we identify DNVs in the hs737 enhancer in individuals with autism. Through multiple approaches, we find hs737 targets the gene EBF3 that is genome-wide significant in NDDs. By assessment of noncoding variation and the genes they affect, we are beginning to understand their impact on gene regulatory networks in NDDs.

59 BASIC BIOLOGICAL SCIENCES↗

Differential Etv2 threshold requirement for endothelial and erythropoietic development

Endothelial and erythropoietic lineages arise from a common developmental progenitor. Etv2 is a master transcriptional regulator required for the development of both lineages. However, the mechanisms through which Etv2 initiates the gene-regulatory networks (GRNs) for endothelial and erythropoietic specification and how the two GRNs diverge downstream of Etv2 remain incompletely understood. Here, by analyzing a hypomorphic Etv2 mutant, we demonstrate different threshold requirements for initiation of the downstream GRNs for endothelial and erythropoietic development. We show that Etv2 functions directly in a coherent feedforward transcriptional network for vascular endothelial development, and a low level of Etv2 expression is sufficient to induce and sustain the endothelial GRN. In contrast, Etv2 induces the erythropoietic GRN indirectly via activation of Tal1, which requires a significantly higher threshold of Etv2 to initiate and sustain erythropoietic development. These results provide important mechanistic insight into the divergence of the endothelial and erythropoietic lineages.

59 BASIC BIOLOGICAL SCIENCES↗

Engineering activatable promoters for scalable and multi-input CRISPRa/i circuits

Dynamic, multi-input gene regulatory networks (GRNs) are ubiquitous in nature. Multilayer CRISPR-based genetic circuits hold great promise for building GRNs akin to those found in naturally occurring biological systems. We develop an approach for creating high-performing activatable promoters that can be assembled into deep, wide, and multi-input CRISPR-activation and -interference (CRISPRa/i) GRNs. By integrating sequence-based design and in vivo screening, we engineer activatable promoters that achieve up to 1,000-fold dynamic range in an Escherichia coli-based cell-free system. These components enable CRISPRa GRNs that are six layers deep and four branches wide. We show the generalizability of the promoter engineering workflow by improving the dynamic range of the light-dependent EL222 optogenetic system from 6-fold to 34-fold. Additionally, high dynamic range promoters enable CRISPRa systems mediated by small molecules and protein–protein interactions. We apply these tools to build input-responsive CRISPRa/i GRNs, including feedback loops, logic gates, multilayer cascades, and dynamic pulse modulators. Our work provides a generalizable approach for the design of high dynamic range activatable promoters and enables classes of gene regulatory functions in cell-free systems.

42 ENGINEERING↗

Sorghum bicolor BTx623 Nitrogen Grown Conditions Set2 Gene Expression Profiling

Dhurrin, a cyanogenic glucoside, plays an important role in Sorghum bicolor physiology and defense. The concentration of dhurrin in sorghum is influenced by both nitrogen status and stage of plant organ development. While nitrogen resupply activates the expression of genes for dhurrin biosynthesis, the molecular mechanisms underlying this regulation remain unclear. In this study, we investigated the transcriptional response of sorghum to nitrogen resupply following growth under nitrogen-limiting conditions. Using a time-course design, we measured hydrogen cyanide potential (HCNp), growth, and nitrate content at 0-, 2-, 6-, 12-, 24-, 36-, 48-, and 60-h after resupply and collected tissue for RNAseq analysis in parallel for analysis of gene expression and construction of gene regulatory networks (GRNs). HCNp (mg g−1 DW) increased significantly in leaf and stem tissues following nitrogen resupply, with increases in the leaf partially driven by continued declines in controls under ongoing nitrogen stress. Expression of the dhurrin pathway genes was upregulated in leaves from 24 h after nitrogen resupply, with diel expression patterns observable over the remaining time points. No upregulation was observed in roots or stems, suggesting that developmental context overrides environmental cues. GRN analysis identified candidate transcription factors regulating dhurrin biosynthesis genes, including members of the MYB, bZIP, and GARP-type transcription factor families. Some of these candidate transcription factors may be involved in relieving senescence-associated suppression of dhurrin biosynthesis and link nitrogen signaling to pathway activation. These findings provide new insight into the nitrogen-responsive regulation of dhurrin in sorghum, highlighting candidate regulators for future functional characterization.

cyanogenic glucoside↗

Sorghum bicolor BTx623 Nitrogen Grown Conditions Gene Expression Profiling

Dhurrin, a cyanogenic glucoside, plays an important role in Sorghum bicolor physiology and defense. The concentration of dhurrin in sorghum is influenced by both nitrogen status and stage of plant organ development. While nitrogen resupply activates the expression of genes for dhurrin biosynthesis, the molecular mechanisms underlying this regulation remain unclear. In this study, we investigated the transcriptional response of sorghum to nitrogen resupply following growth under nitrogen-limiting conditions. Using a time-course design, we measured hydrogen cyanide potential (HCNp), growth, and nitrate content at 0-, 2-, 6-, 12-, 24-, 36-, 48-, and 60-h after resupply and collected tissue for RNAseq analysis in parallel for analysis of gene expression and construction of gene regulatory networks (GRNs). HCNp (mg g−1 DW) increased significantly in leaf and stem tissues following nitrogen resupply, with increases in the leaf partially driven by continued declines in controls under ongoing nitrogen stress. Expression of the dhurrin pathway genes was upregulated in leaves from 24 h after nitrogen resupply, with diel expression patterns observable over the remaining time points. No upregulation was observed in roots or stems, suggesting that developmental context overrides environmental cues. GRN analysis identified candidate transcription factors regulating dhurrin biosynthesis genes, including members of the MYB, bZIP, and GARP-type transcription factor families. Some of these candidate transcription factors may be involved in relieving senescence-associated suppression of dhurrin biosynthesis and link nitrogen signaling to pathway activation. These findings provide new insight into the nitrogen-responsive regulation of dhurrin in sorghum, highlighting candidate regulators for future functional characterization.

cyanogenic glucoside↗

Sorghum SbGhd7 is a major regulator of floral transition and directly represses genes crucial for flowering activation

Molecular genetic understanding of flowering time regulation is crucial for sorghum development. GRAIN NUMBER, PLANT HEIGHT AND HEADING DATE 7 (SbGhd7) is one of the six classical loci conferring photoperiod sensitivity of sorghum flowering. However, its functions remain poorly studied. The molecular functions of SbGhd7 were characterized. The gene regulatory network controlled by SbGhd7 was constructed and validated. The biological roles of SbGhd7 and its major targets were studied. SbGhd7overexpression (OE) completely prevented sorghum flowering. Additionally, we show thatSbGhd7is a major negative regulator of flowering, binding to the promoter motif TGAATG(A/T)(A/T/C) and repressing transcription of the major florigen FLOWERING LOCUS T 10(SbFT10) and floral activators EARLY HEADING DATE(SbEhd1),FLAVIN-BINDING,KELCH REPEAT,F-BOX1(SbFKF1) and EARLY FLOWERING 3(SbELF3). Reinforcing the direct effect ofSbGhd7,SbEhd1OE activated the promoters of three functional florigens (SbFT1,SbFT8andSbFT10), dramatically accelerating flowering. Here our studies demonstrate that SbGhd7 is a major repressor of sorghum flowering by directly and indirectly targeting genes for flowering activation. The mechanism appears ancient. Our study extends the current model of floral transition regulation in sorghum and provides a framework for a comprehensive understanding of sorghum photoperiod response.

59 BASIC BIOLOGICAL SCIENCES↗

A novel regulator of the fungal phosphate starvation response revealed by transcriptional profiling and DNA affinity purification sequencing

Cells must accurately sense and respond to nutrients to compete for resources and establish growth. Phosphate is a critical nutrient source necessary for signaling, energy metabolism, and synthesis of nucleic acids, phospholipids, and cellular metabolites. During phosphate limitation, fungi import phosphate from the environment and liberate phosphate from phosphate-containing molecules in the cell. In the model filamentous fungus Neurospora crassa, the phosphate starvation response is regulated by the conserved transcription factor NUC-1. The activity of NUC-1 is repressed by a complex of the cyclin-dependent kinase MDK-1 and the cyclin PREG when phosphate is plentiful. When phosphate is limiting, NUC-1 repression by MDK-1/PREG is relieved by the cyclin-dependent kinase inhibitor NUC-2. We investigated the global response of N. crassa to phosphate starvation. During phosphate starvation, NUC-1 directly activated the expression of genes encoding phosphatases, nucleases, and a phosphate transporter and directly repressed genes associated with the ribosome. Additionally, NUC-1 indirectly activated the expression of an uncharacterized transcription factor, which we named nuc-3. NUC-3 directly repressed the expression of genes involved in phosphate acquisition and liberation after an extended period of phosphate starvation. Additionally, NUC-3 directly repressed the expression of the cyclin-dependent kinase inhibitor nuc-2. Thus, through the combination of NUC-3 direct repression of genes in the phosphate starvation response and nuc-2, an activator of the phosphate starvation response, NUC-3 serves to act as a brake on the phosphate starvation response after an extended period of phosphate starvation. This braking mechanism could reduce transcription, a phosphate-intensive process, under conditions of extended phosphate limitation.IMPORTANCEFungi have evolved regulatory networks to respond to available nutrients. Phosphate is often a limiting nutrient for fungi that is critical for many cellular functions, including nucleic acid and phospholipid biosynthesis, cell signaling, and energy metabolism. The fungal response to phosphate limitation is important in interactions with plants and animals. We investigated the global transcriptional response to phosphate starvation and the role of a major transcriptional regulator, NUC-1, in the model filamentous fungus Neurospora crassa. Our data show that NUC-1 is a bifunctional transcription factor that directly activates phosphate acquisition genes, while directly repressing genes associated with phosphate-intensive processes. NUC-1 indirectly regulates an uncharacterized transcription factor, which we named nuc-3. NUC-3 directly represses phosphate acquisition genes and nuc-2, an activator of the phosphate starvation response, during extended periods of phosphate starvation. Thus, NUC-3 acts as a brake on the phosphate starvation response to reduce phosphate-intensive activities, like transcriptional activation, when phosphate starvation persists.

DNA affinity purification sequencing↗

Nitrogen Status Rewires Transcriptional Regulation of Dhurrin, a Dual‐Purpose Defense Metabolite in Sorghum bicolor

Dhurrin, a cyanogenic glucoside, plays an important role in Sorghum bicolor physiology and defense. The concentration of dhurrin in sorghum is influenced by both nitrogen status and stage of plant organ development. While nitrogen resupply activates the expression of genes for dhurrin biosynthesis, the molecular mechanisms underlying this regulation remain unclear. In this study, we investigated the transcriptional response of sorghum to nitrogen resupply following growth under nitrogen-limiting conditions. Using a time-course design, we measured hydrogen cyanide potential (HCNp), growth, and nitrate content at 0-, 2-, 6-, 12-, 24-, 36-, 48-, and 60-h after resupply and collected tissue for RNAseq analysis in parallel for analysis of gene expression and construction of gene regulatory networks (GRNs). HCNp (mg g −1 DW) increased significantly in leaf and stem tissues following nitrogen resupply, with increases in the leaf partially driven by continued declines in controls under ongoing nitrogen stress. Expression of the dhurrin pathway genes was upregulated in leaves from 24 h after nitrogen resupply, with diel expression patterns observable over the remaining time points. No upregulation was observed in roots or stems, suggesting that developmental context overrides environmental cues. GRN analysis identified candidate transcription factors regulating dhurrin biosynthesis genes, including members of the MYB, bZIP, and GARP-type transcription factor families. Some of these candidate transcription factors may be involved in relieving senescence-associated suppression of dhurrin biosynthesis and link nitrogen signaling to pathway activation. These findings provide new insight into the nitrogen-responsive regulation of dhurrin in sorghum, highlighting candidate regulators for future functional characterization.

S. bicolor↗

The trans-regulatory landscape of gene networks in plants

The transcriptional effector domains of transcription factors play a key role in controlling gene expression; however, their functional nature is poorly understood, hampering our ability to explore this fundamental dimension of gene regulatory networks. To map the trans-regulatory landscape in a complex eukaryote, we systematically characterized the putative transcriptional effector domains of over 400 Arabidopsis thaliana transcription factors for their capacity to modulate transcription. We demonstrate that transcriptional effector activity can be integrated into gene regulatory networks capable of elucidating the functional dynamics underlying gene expression patterns. We further show how our characterized domains can enhance genome engineering efforts and reveal how plant transcriptional activators share regulatory features conserved across distantly related eukaryotes. Our results provide a framework to systematically characterize the regulatory role of transcription factors at a genome-scale in order to understand the transcriptional wiring of biological systems.

59 BASIC BIOLOGICAL SCIENCES↗

Through the lens of bioenergy crops: advances, bottlenecks, and promises of plant engineering

Advances in engineering of bioenergy crops were driven over the past years by adapting technological breakthroughs and accelerating conventional applications but also exposed intriguing challenges. New tools revealed rich interconnectivity in the exponentially growing and dynamic 'big' omics data' of metabolomes, transcriptomes, and genomes at previously inaccessible magnitude (global, cross-species, meta-) and resolution (single cell). Insights enabled fresh hypotheses and stimulated disciplines such as functional genomics with discovery of broad regulatory networks and their determinants, that is, DNA parts, including promoters, regulatory elements, and transcription factors. Their rational design, assembly into increasingly complex blueprints, and installation into diverse chassis is an existing frontier that may benefit from emerging technologies to address bottlenecks. Interweaving nature-inspired to fully synthetic parts has already allowed building of fine-tuned regulatory circuits, or new-to-nature metabolic routes insulated from the biological context of the chassis species. Similarly, developments and the evolving need for unifying principles in plant transformation and species-agnostic technologies highlight future opportunities for engineering the next generation of bioenergy plants.

60 APPLIED LIFE SCIENCES↗

Gene network centrality analysis identifies key regulators coordinating day-night metabolic transitions in Synechococcus elongatus PCC 7942 despite limited accuracy in predicting direct regulator-gene interactions

Synechococcus elongatus PCC 7942 is a model organism for studying circadian regulation and bioproduction, where precise temporal control of metabolism significantly impacts photosynthetic efficiency and CO 2 -to-bioproduct conversion. Despite extensive research on core clock components, our understanding of the broader regulatory network orchestrating genome-wide metabolic transitions remains incomplete. We address this gap by applying machine learning tools and network analysis to investigate the transcriptional architecture governing circadian-controlled gene expression. While our approach showed moderate accuracy in predicting individual transcription factor-gene interactions - a common challenge with real expression data - network-level topological analysis successfully revealed the organizational principles of circadian regulation. Our analysis identified distinct regulatory modules coordinating day-night metabolic transitions, with photosynthesis and carbon/nitrogen metabolism controlled by day-phase regulators, while nighttime modules orchestrate glycogen mobilization and redox metabolism. Through network centrality analysis, we identified potentially significant but previously understudied transcriptional regulators: HimA as a putative DNA architecture regulator, and TetR and SrrB as potential coordinators of nighttime metabolism, working alongside established global regulators RpaA and RpaB. This work demonstrates how network-level analysis can extract biologically meaningful insights despite limitations in predicting direct regulatory interactions. The regulatory principles uncovered here advance our understanding of how cyanobacteria coordinate complex metabolic transitions and may inform metabolic engineering strategies for enhanced photosynthetic bioproduction from CO 2 .

59 BASIC BIOLOGICAL SCIENCES↗

Spurious regulatory connections dictate the expression-fitness landscape of translation factors

During steady-state cell growth, individual enzymatic fluxes can be directly inferred from growth rate by mass conservation, but the inverse problem remains unsolved. Perturbing the flux and expression of a single enzyme could have pleiotropic effects that may or may not dominate the impact on cell fitness. Here, we quantitatively dissect the molecular and global responses to varied expression of translation termination factors (peptide release factors, RFs) in the bacterium Bacillus subtilis. While endogenous RF expression maximizes proliferation, deviations in expression lead to unexpected distal regulatory responses that dictate fitness reduction. Molecularly, RF depletion causes expression imbalance at specific operons, which activates master regulators and detrimentally overrides the transcriptome. Through these spurious connections, RF abundances are thus entrenched by focal points within the regulatory network, in one case located at a single stop codon. Such regulatory entrenchment suggests that predictive bottom-up models of expression-fitness landscapes will require near-exhaustive characterization of parts.

59 BASIC BIOLOGICAL SCIENCES↗

Global transcriptome profile of the developmental principles of in vitro iPSC-to-motor neuron differentiation

Abstract Background Human induced pluripotent stem cells (iPSC) have opened new avenues for regenerative medicine. Consequently, iPSC-derived motor neurons have emerged as potentially viable therapies for spinal cord injuries and neurodegenerative disorders including Amyotrophic Lateral Sclerosis. However, direct clinical application of iPSC bears in itself the risk of tumorigenesis and other unforeseeable genetic or epigenetic abnormalities. Results Employing RNA-seq technology, we identified and characterized gene regulatory networks triggered by in vitro chemical reprogramming of iPSC into cells with the molecular features of motor neurons (MNs) whose function in vivo is to innervate effector organs. We present meta-transcriptome signatures of 5 cell types: iPSCs, neural stem cells, motor neuron progenitors, early motor neurons, and mature motor neurons. In strict response to the chemical stimuli, along the MN differentiation axis we observed temporal downregulation of tumor growth factor-β signaling pathway and consistent activation of sonic hedgehog, Wnt/β-catenin, and Notch signaling. Together with gene networks defining neuronal differentiation (neurogenin 2, microtubule-associated protein 2, Pax6, and neuropilin-1), we observed steady accumulation of motor neuron-specific regulatory genes, including Islet-1 and homeobox protein HB9. Interestingly, transcriptome profiling of the differentiation process showed that Ca 2+ signaling through cAMP and LPC was downregulated during the conversion of the iPSC to neural stem cells and key regulatory gene activity of the pathway remained inhibited until later stages of motor neuron formation. Pathways shaping the neuronal development and function were well-represented in the early motor neuron cells including, neuroactive ligand-receptor interactions, axon guidance, and the cholinergic synapse formation. A notable hallmark of our in vitro motor neuron maturation in monoculture was the activation of genes encoding G-coupled muscarinic acetylcholine receptors and downregulation of the ionotropic nicotinic acetylcholine receptors expression. We observed the formation of functional neuronal networks as spontaneous oscillations in the extracellular action potentials recorded on multi-electrode array chip after 20 days of differentiation. Conclusions Detailed transcriptome profile of each developmental step from iPSC to motor neuron driven by chemical induction provides the guidelines to novel therapeutic approaches in the re-construction efforts of muscle innervation.

59 BASIC BIOLOGICAL SCIENCES↗

Machine learning uncovers independently regulated modules in the Bacillus subtilis transcriptome

The transcriptional regulatory network (TRN) of Bacillus subtilis coordinates cellular functions of fundamental interest, including metabolism, biofilm formation, and sporulation. Here, we use unsupervised machine learning to modularize the transcriptome and quantitatively describe regulatory activity under diverse conditions, creating an unbiased summary of gene expression. We obtain 83 independently modulated gene sets that explain most of the variance in expression and demonstrate that 76% of them represent the effects of known regulators. The TRN structure and its condition-dependent activity uncover putative or recently discovered roles for at least five regulons, such as a relationship between histidine utilization and quorum sensing. The TRN also facilitates quantification of population-level sporulation states. As this TRN covers the majority of the transcriptome and concisely characterizes the global expression state, it could inform research on nearly every aspect of transcriptional regulation in B. subtilis.

59 BASIC BIOLOGICAL SCIENCES↗

Eliater: a Python package for estimating outcomes of perturbations in biomolecular networks

We introduce Eliater, a Python package for estimating the effect of perturbation of an upstream molecule on a downstream molecule in a biomolecular network. The estimation takes as input a biomolecular network, observational biomolecular data, and a perturbation of interest, and outputs an estimated quantitative effect of the perturbation. We showcase the functionalities of Eliater in a case study of Escherichia coli transcriptional regulatory network.

59 BASIC BIOLOGICAL SCIENCES↗

TroR is the primary regulator of the iron homeostasis transcription network in the halophilic archaeon Haloferax volcanii

Abstract Maintaining the intracellular iron concentration within the homeostatic range is vital to meet cellular metabolic needs and reduce oxidative stress. Previous research revealed that the haloarchaeon Halobacterium salinarum encodes four diphtheria toxin repressor (DtxR) family transcription factors (TFs) that together regulate the iron response through an interconnected transcriptional regulatory network (TRN). However, the conservation of the TRN and the metal specificity of DtxR TFs remained poorly understood. Here we identified and characterized the TRN of Haloferax volcanii for comparison. Genetic analysis demonstrated that Hfx. volcanii relies on three DtxR transcriptional regulators (Idr, SirR, and TroR), with TroR as the primary regulator of iron homeostasis. Bioinformatics and molecular approaches revealed that TroR binds a conserved cis-regulatory motif located ∼100 nt upstream of the start codon of iron-related target genes. Transcriptomics analysis demonstrated that, under conditions of iron sufficiency, TroR repressed iron uptake and induced iron storage mechanisms. TroR repressed the expression of one other DtxR TF, Idr. This reduced DtxR TRN complexity relative to that of Hbt. salinarum appeared correlated with natural variations in iron availability. Based on these data, we hypothesize that variable environmental conditions such as iron availability appear to select for increasing TRN complexity.

59 BASIC BIOLOGICAL SCIENCES↗

SOX2-driven enhancer landscape defines the transcriptional architecture of retinogenesis

Retinal neurogenesis is mediated by the coordinated activities of a complex gene regulatory network (GRN) of transcription factors (TFs) in multipotent retinal progenitor cells (RPCs). How this GRN mechanistically guides neural competence remains poorly understood. In this study, we present integrated transcriptional, genetic and genomic analyses to uncover the regulatory mechanisms of SOX2, a key factor in establishing neural identity in RPCs. We show that SOX2 is preferentially enriched in the RPC-specific enhancer landscape associated with essential regulators of retinogenesis. Disruption of SOX2 expression impairs retinogenesis, marked by a selective loss of enhancer activity near genes essential for RPC proliferation and lineage specification. We identified the RPC transcription factor VSX2 as a binding partner for SOX2 and, together, SOX2 and VSX2 co-target a core, retina-specific chromatin repertoire characterized by enhanced TF binding and robust chromatin accessibility. This cooperative binding establishes a shared SOX2-VSX2 transcriptional code that promotes the expression of crucial regulators of neurogenesis while repressing the acquisition of alternative lineage cell fate. Our data illuminate fundamental biological insights on how transcription factors act in concert to drive chromatin-based genetic programs underlying retinal neural identity.

Chromatin↗

Sensitivity Analysis of Genome-Scale Metabolic Flux Prediction

TRIMER, Transcription Regulation Integrated with MEtabolic Regulation, is a genome-scale modeling pipeline targeting at metabolic engineering applications. Using TRIMER, regulated metabolic reactions can be effectively predicted by integrative modeling of metabolic reactions with a Transcription Factor (TF)-gene regulatory network (TRN), which is modeled via a Bayesian network (BN). In this paper, we focus on sensitivity analysis of metabolic flux prediction for uncertainty quantification of BN structures for TRN modeling in TRIMER. We propose a computational strategy to construct the uncertainty class of TRN models based on the inferred regulatory order uncertainty given transcriptomic expression data. With that, we analyze the prediction sensitivity of the TRIMER pipeline for the metabolite yields of interest. The obtained sensitivity analyses can guide Optimal Experimental Design (OED) to help acquire new data that can enhance TRN modeling and achieve specific metabolic engineering objectives, including metabolite yield alterations. Here we have performed small- and large-scale simulated experiments, demonstrating the effectiveness of our developed sensitivity analysis strategy for BN structure learning to quantify the edge importance in terms of metabolic flux prediction uncertainty reduction and its potential to effectively guide OED.

59 BASIC BIOLOGICAL SCIENCES↗