Identification of a rate-limiting step in a metabolic pathway using the kinetic model and in vitro experiment

OsTH1 is a key player in thiamin biosynthesis in rice

Thiamin is a vital nutrient that acts as a cofactor for several enzymes primarily localized in the mitochondria. These thiamin-dependent enzymes are involved in energy metabolism, nucleic acid biosynthesis, and antioxidant machinery. The enzyme HMP-P kinase/thiamin monophosphate synthase (TH1) holds a key position in thiamin biosynthesis, being responsible for the phosphorylation of HMP-P into HMP-PP and for the condensation of HMP-PP and HET-P to form TMP. Through mathematical kinetic model, we have identified TH1 as a critical player for thiamin biofortification in rice. We further focused on the functional characterization of OsTH1. Sequence and gene expression analysis, along with phylogenetic studies, provided insights into OsTH1 bifunctional features and evolution. The indispensable role of OsTH1 in thiamin biosynthesis was validated by heterologous expression of OsTH1 and successful complementation of yeast knock-out mutants impaired in thiamin production. We also proved that the sole OsTH1 overexpression in rice callus significantly improves B1 concentration, resulting in 50% increase in thiamin accumulation. Our study underscores the critical role of OsTH1 in thiamin biosynthesis, shedding light on its bifunctional nature and evolutionary significance. The significant enhancement of thiamin accumulation in rice callus upon OsTH1 overexpression constitutes evidence of its potential application in biofortification strategies.

Utilizing 5' UTR Engineering Enables Fine-Tuning of Multiple Genes within Operons to Balance Metabolic Flux in Bacillus subtilis

The application of synthetic biology tools to modulate gene expression to increase yield has been thoroughly demonstrated as an effective and convenient approach in industrial production. In this study, we employed a high-throughput screening strategy to identify a 5' UTR sequence from the genome of B. subtilis 168. This sequence resulted in a 5.8-fold increase in the expression level of EGFP. By utilizing the 5' UTR sequence to overexpress individual genes within the rib operon, it was determined that the genes ribD and ribAB serve as rate-limiting enzymes in the riboflavin synthesis pathway. Constructing a 5' UTR library to regulate EGFP expression resulted in a variation range in gene expression levels exceeding 100-fold. Employing the same 5' UTR library to regulate the expression of EGFP and mCherry within the operon led to a change in the expression ratio of these two genes by over 10,000-fold. So, employing a 5' UTR library to modulate the expression of the rib operon gene and construct a synthetic rib operon resulted in a 2.09-fold increase in riboflavin production. These results indicate that the 5' UTR sequence identified and characterized in this study can serve as a versatile synthetic biology toolkit for achieving complex metabolic network reconstruction. This toolkit can facilitate the fine-tuning of gene expression to produce target products.

Trace impurities in sodium phosphate influences the physiological activity of Escherichia coli in M9 minimal medium

In the field of applied microbiology, reproducibility and experimental variability are important factors that influence both basic research as well as process development for industrial applications. Experimental reproducibility and accuracy depend not only on culture conditions such as temperature and aeration but also on raw materials and procedures used for media preparation. The M9 minimal medium is one of the most common synthetic media for culturing Escherichia coli and other bacteria. This synthetic medium can be used to observe and evaluate the physiological activity of microbes under minimal nutritional requirements and determine the limiting factor for the desired phenotype. Although one of the advantages using the M9 medium is that its composition can be modulated, it is difficult to control presence of trace components and impurities from the reagents for preparing this medium. Herein, we showed that trace ingredients present in the reagents used for M9 media preparation affect the bacterial physiological activities (e.g., cell growth, substrate consumption, and byproduct formation). Additionally, we systematically identified the trace ingredient that influenced phenotypic differences. Our results showed that the selection of reagents and accuracy during reagent preparation is important for experimental reproducibility in the field of bio-engineering and systems biology focused on the systematic and continuous development of biomolecular systems (e.g., biorefinery, metabolic engineering, and synthetic biology).

Pichia pastoris growth-coupled heme biosynthesis analysis using metabolic modelling

Soy leghemoglobin is one of the most important and key ingredients in plant-based meat substitutes that can imitate the colour and flavour of the meat. To improve the high-yield production of leghemoglobin protein and its main component-heme in the yeast Pichia pastoris, glycerol and methanol cultivation conditions were studied. Additionally, in-silico metabolic modelling analysis of growth-coupled enzyme quantity, suggests metabolic gene up/down-regulation strategies for heme production. First, cultivations and metabolic modelling analysis of P. pastoris were performed on glycerol and methanol in different growth media. Glycerol cultivation uptake and production rates can be increased by 50% according to metabolic modelling results, but methanol cultivation-is near the theoretical maximum. Growth-coupled metabolic optimisation results revealed the best feasible upregulation (33 reactions) (1.47% of total reactions) and 66 downregulation/deletion (2.98% of total) reaction suggestions. Finally, we describe reaction regulation suggestions with the highest potential to increase heme production yields.

Fundamental behaviors emerge from simulations of a living minimal cell

We present a whole-cell fully dynamical kinetic model (WCM) of JCVI-syn3A, a minimal cell with a reduced genome of 493 genes that has retained few regulatory proteins or small RNAs. Cryo-electron tomograms provide the cell geometry and ribosome distributions. Time-dependent behaviors of concentrations and reaction fluxes from stochastic-deterministic simulations over a cell cycle reveal how the cell balances demands of its metabolism, genetic information processes, and growth, and offer insight into the principles of life for this minimal cell. The energy economy of each process including active transport of amino acids, nucleosides, and ions is analyzed. WCM reveals how emergent imbalances lead to slowdowns in the rates of transcription and translation. Integration of experimental data is critical in building a kinetic model from which emerges a genome-wide distribution of mRNA half-lives, multiple DNA replication events that can be compared to qPCR results, and the experimentally observed doubling behavior.

Genome-wide characterization, evolution, and expression profiling of FBA gene family in response to light treatments and abiotic stress in Nicotiana tabacum

Fructose 1,6-bisphosphate aldolase (FBA) as a key enzyme play crucial roles in glycolysis, gluconeogenesis and Calvin cycle processes in plants. However, limited information is known regarding FBA genes in Nicotiana tabacum. In this study, 16 FBAs were identified and characterized in Nicotiana tabacum. Phylogenetic analysis revealed that these genes can be categorized as type I (NtFBA1-10 located in chloroplast) and type II (NtFBA11-16 located in cytoplasm) subfamilies. According to the conserved motifs and gene structure analysis, NtFBA protein sequences had the highly homologous to FBAs in other species. Most members of the NtFBA gene family responded positively to NaHCO3 stress, especially the expression of NtFBA13/14 increased by 642%. In addition, the expression results of NtFBAs under five abiotic stress (light, NaCl, NaHCO3, drought, and cold) conditions were showed that NtFBA13/14 were highly up-regulated. qRT-PCR results showed that most of the NtFBAs expressed higher in leaves. NtFBA7/8 and NtFBA13/14 have important significance in photosynthesis and abiotic stress, respectively. This study provides a basis foundation for further elucidating the function of NtFBAs and the N. tabacum mechanism of resistance under abiotic stress.

Recent advances in metabolic engineering-integration of in silico design and experimental analysis of metabolic pathways

Microorganisms are widely used to produce valuable compounds. Because thousands of metabolic reactions occur simultaneously and many metabolic reactions are related to target production and cell growth, the development of a rational design method for metabolic pathway modification to optimize target production is needed. In this paper, recent advances in metabolic engineering are reviewed, specifically considering computational pathway modification design and experimental evaluation of metabolic fluxes by ¹³C-metabolic flux analysis. Computational tools for seeking effective gene deletion targets and recruiting heterologous genes are described in flux balance analysis approaches. A kinetic model and adaptive laboratory evolution were applied to identify and eliminate the rate-limiting step in metabolic pathways. Data science-based approaches for process monitoring and control are described to maximize the performance of engineered cells in bioreactors.

Engineering analysis of multienzyme cascade reactions for 3'-sialyllactose synthesis

Sialo‐oligosaccharides are important products of emerging biotechnology for complex carbohydrates as nutritional ingredients. Cascade bio‐catalysis is central to the development of sialo‐oligosaccharide production systems, based on isolated enzymes or whole cells. Multienzyme transformations have been established for sialo‐oligosaccharide synthesis from expedient substrates, but systematic engineering analysis for the optimization of such transformations is lacking. Here, we show a mathematical modeling‐guided approach to 3ʹ‐sialyllactose (3SL) synthesis from N‐acetyl‐ d‐neuraminic acid (Neu5Ac) and lactose in the presence of cytidine 5ʹ‐triphosphate, via the reactions of cytidine 5ʹ‐monophosphate‐Neu5Ac synthetase and α2,3‐sialyltransferase. The Neu5Ac was synthesized in situ from N‐acetyl‐ d‐mannosamine using the reversible reaction with pyruvate by Neu5Ac lyase or the effectively irreversible reaction with phosphoenolpyruvate by Neu5Ac synthase. We show through comprehensive time‐course study by experiment and modeling that, due to kinetic rather than thermodynamic advantages of the synthase reaction, the 3SL yield was increased (up to 75%; 10.4 g/L) and the initial productivity doubled (15 g/L/h), compared with synthesis based on the lyase reaction. We further show model‐based optimization to minimize the total loading of protein (saving: up to 43%) while maintaining a suitable ratio of the individual enzyme activities to achieve 3SL target yield (61%–75%; 7–10 g/L) and overall productivity (3–5 g/L/h). Collectively, our results reveal the principal factors of enzyme cascade efficiency for 3SL synthesis and highlight the important role of engineering analysis to make multienzyme‐catalyzed transformations fit for oligosaccharide production.