Soy protein hydrolysates induce menaquinone-7 biosynthesis by enhancing the biofilm formation of Bacillus subtilis natto.

Menaquinone-7 (MK-7) is a form of vitamin K2 with health-beneficial effects. A novel fermentation strategy based on combining soy protein hydrolysates (SPHs) with biofilm-based fermentation was investigated to enhance menaquinone-7 (MK-7) biosynthesis by Bacillus subtilis natto. Results showed the SPHs increased MK-7 yield by 199.4% in two-stage aeration fermentation as compared to the SP-based medium in submerged fermentation, which was related to the formation of robust biofilm with wrinkles and the enhancement of cell viability. Moreover, there was a significant correlation between key genes related to MK-7 and biofilm synthesis, and the quorum sensing (QS) related genes, Spo0A and SinR, were downregulated by 0.64-fold and 0.39-fold respectively, which promoted biofilm matrix synthesis. Meanwhile, SPHs also enhanced the MK-7 precursor, isoprene side chain, supply, and MK-7 assembly efficiency. Improved fermentation performances of bacterial cells during fermentation were attributed to abundant oligopeptides (Mw < 1 kDa) and moderate amino acids, particularly Arg, Asp, and Phe in SPHs. All these results revealed that SPHs were a potential and superior nitrogen source for MK-7 production by Bacillus subtilis natto.

Evaluating the dynamic effects of complex probiotics as cellulase replacements during fermentation of apple pomace.

Apple pomace, an abundant agricultural by-product with low utilization rates, often leads to environmental pollution if not properly managed. This study aimed to enhance the nutritional value of apple pomace by comparing the effects of solid-state fermentation using complex probiotics and cellulase preparation. Additionally, the study investigated the dynamic changes in various components throughout the fermentation process with complex probiotics. The results of single-strain solid-state fermentation tests indicated that Lactiplantibacillus plantarum DPH, Saccharomyces cerevisiae SC9, and Bacillus subtilis C9 were the optimal strains for fermenting the most effective substrate combination, comprising 73 % apple pomace and 20 % millet bran. The strains (complex probiotics) and a cellulase preparation were used for the solid-state fermentation of the apple pomace mixture for nine days, respectively. The contents of acid detergent fiber, neutral detergent fiber, hemicellulose, and insoluble dietary fiber decreased by up to 9.99 %, 9.59 %, 23.21 %, and 14.34 %, respectively. In contrast, the content of soluble dietary fiber significantly increased by up to 29.74 %. Both methods reduced cellulose crystallinity and modified the substrate's surface structure, resulting in a looser arrangement. Fermentation with complex probiotics for three or six days increased the abundance of lactic acid bacteria, which comprised >87 % of the total microbial population. Concurrently, the abundance of detrimental bacteria, such as Salmonella, Acetobacter, Escherichia, and Pantoea, significantly decreased. Furthermore, fermentation with complex probiotics for six or nine days enhanced antioxidant properties, leading to a significant increase in beneficial metabolites, including amino acids, organic acids, gamma-aminobutyric acid, serotonin. In conclusion, complex probiotics can effectively substitute for cellulase preparation in the solid-state fermentation of apple pomace, with a six-day fermentation period yielding optimal results. This study provides valuable insights into enhancing the value of apple pomace in the feed industry and the effective application of agro-industrial by-products.

How bacteria actively use passive physics to make biofilms.

Modern molecular microbiology elucidates the organizational principles of bacterial biofilms via detailed examination of the interplay between signaling and gene regulation. A complementary biophysical approach studies the mesoscopic dependencies at the cellular and multicellular levels with a distinct focus on intercellular forces and mechanical properties of whole biofilms. Here, motivated by recent advances in biofilm research and in other, seemingly unrelated fields of biology and physics, we propose a perspective that links the biofilm, a dynamic multicellular organism, with the physical processes occurring in the extracellular milieu. Using Bacillus subtilis as an illustrative model organism, we specifically demonstrate how such a rationale explains biofilm architecture, differentiation, communication, and stress responses such as desiccation tolerance, metabolism, and physiology across multiple scales-from matrix proteins and polysaccharides to macroscopic wrinkles and water-filled channels.

Biotransformation of Sumatriptan by Staphylococcus aureus, Bacillus subtilis, Pseudomonas aeruginosa and Salmonella enterica subsp. enterica.

This study aimed at the biotransformation of sumatriptan by Staphylococcus aureus, Bacillus subtilis, Pseudomonas aeruginosa and Salmonella enterica subsp. enterica and the identification of the drug metabolites by liquid chromatography-mass spectrometry. The drug was incubated with the organisms in tryptic soya broth at 37 °C. The broth was filtered and subjected to liquid chromatography-mass spectrometry. The metabolites identified by the use of mass spectral (+ve ion mode) fragmentation patterns were (3-methylphenyl)methanethiol (Bacillus subtilis), 1-(4-amino-3-ethylphenyl)-N-methylmethanesulfonamide (Salmonella enterica subsp. enterica) and 1-{4-amino-3-[(1E)-3-(dimethylamino)prop-1-en-1-yl]phenyl}methanesulfinamide (Salmonella enterica subsp. enterica, Bacillus subtilis, Pseudomonas aeruginosa, Staphylococcus aureus). These metabolites exhibit high gastrointestinal absorption, no blood-brain barrier permeability (except (3-methylphenyl)methanethiol), a bioavailability score of 0.55 and no inhibitory effect on CYP2C19, CYP2C9, CYP2D6, CYP3A4 or cytochrome P450 1A2 (except (3-methylphenyl)methanethiol), as determined by SwissADME software ver. 2024. The metabolites appear to be more toxic than the parent drug, as suggested by their calculated median lethal dose values. All four organisms under investigation transformed sumatriptan to different chemical substances that were more toxic than the parent drug.

[Efficient biosynthesis of guanidoacetic acid by a recombinant strain of Bacillus subtilis].

Guanidinoacetic acid, as an energetic substance, has a wide range of applications in the food, pharmaceutical, and feed industries. However, the biosynthesis of guanidinoacetic acid has not been applied in industrial production. In this study, we designed the synthetic route of guanidinoacetic acid in a food-grade strain of Bacillus subtilis. By regulating the expression of key enzymes, lifting feedback inhibition, and increasing membrane permeability, we achieved the efficient synthesis of guanidinoacetic acid by whole-cell catalysis. Firstly, the optimal L-arginine:glycine amidinotransferase was screened based on the phylogenetic tree, and the expression of the key enzyme was enhanced by a strategy combining strong promoter and genome integration. Secondly, the ornithine cycle for L-arginine synthesis in Corynebacterium glutamicum was introduced to alleviate the feedback inhibition of the enzyme by the byproduct L-ornithine, and the L-arginine degradation pathway was knocked down to enhance substrate regeneration. Thirdly, the expression of N-acetylmuramoyl-L-alanine amidase (LytC) was up-regulated to increase the cell membrane permeability. Finally, after optimization of whole-cell production conditions, strain Bs-13 achieved guanidinoacetic acid production at a titer of 13.1 g/L after 24 h, with a proudction rate of 0.54 g/(L·h) and a glycine conversion rate of 92.7%. The above strategy improved the production of guanidinoacetic acid and provided a reference for the biosynthesis of guanidinoacetic acid.

[Efficient whole-cell biosynthesis of D-mannose by recombinant Bacillus subtilis].

D-mannose is a natural hexose with great economic and application values in the food, medicine, and cosmetic fields. However, most biosynthesis methods of D-mannose rely on Escherichia coli as the host, which poses safety issues during the production process and imposes limitations on subsequent applications. This study compared the enzyme properties of mannose isomerases from multiple sources to select the most suitable source. B. subtilis 168/pMA5-EcMIaseA was constructed with "generally recognized as safe" (GRAS) Bacillus subtilis as the host and used as a whole-cell catalyst to synthesize D-mannose from d-fructose. Optimizing the conversion conditions such as culture temperature, pH, and substrate concentration increased the yield of D-mannose. The results showed that the conversion rates reached 27.75% and 27.22% and the yields of D-mannose were 138.74 g/L and 163.30 g/L after 6 h whole-cell transformation with d-fructose at the concentrations of 500 g/L and 600 g/L, respectively, in a 5 L fermentor. This study achieves the highest yield of D-mannose produced under the catalysis by recombinant B. subtilis that has ever been reported and provides a basis for the industrial production and application of D-mannose.

Dynamic in situ detection in iRhizo-Chip reveals diurnal fluctuations of Bacillus subtilis in the rhizosphere.

Effective colonization by microbe in the rhizosphere is critical for establishing a beneficial symbiotic relationship with the host plant. Bacillus subtilis, a soil-dwelling bacterium that is commonly found in association with plants and their rhizosphere, has garnered interest for its potential to enhance plant growth, suppress pathogens, and contribute to sustainable agricultural practices. However, research on the dynamic distribution of B. subtilis within the rhizosphere and its interaction mechanisms with plant roots remains insufficient due to limitations in existing in situ detection methodologies. To achieve dynamic in situ detection of the rhizosphere environment, we established iRhizo-Chip, a microfluidics-based platform. Using this device to investigate microbial behavior within the rhizosphere, we found obvious diurnal fluctuations in the growth of B. subtilis in the rhizosphere. Temporal dynamic analysis of rhizosphere dissolved oxygen (DO), pH, dissolved organic carbon, and reactive oxygen species showed that diurnal fluctuations in the growth of B. subtilis are potentially related to a variety of environmental factors. Spatial dynamic analysis also showed that the spatial distribution changes of B. subtilis and DO and pH were similar. Subsequently, through in vitro control experiments, we proved that rhizosphere DO and pH are the main driving forces for diurnal fluctuations in the growth of B. subtilis. Our results show that the growth of B. subtilis is driven by rhizosphere DO and pH, resulting in diurnal fluctuations, and iRhizo-Chip is a valuable tool for studying plant rhizosphere dynamics.

Double-Wing Motif Protein is a Novel Biofilm Regulatory Factor of the Plant Disease Biocontrol Agent, Bacillus subtilis.

Transposon mutagenesis screening of Bacillus subtilis YB-1471, a novel rhizosphere biocontrol agent of Fusarium crown rot (FCR) of wheat, resulted in the identification of orf04391, linked to reduced biofilm formation. The gene encodes a protein possessing a putative tertiary structure of a "double-wing" DNA-binding domain. Expression of orf04391 increased during biofilm development in stationary cultures and during rapid growth in shaking cultures. An orf04391 deletion strain showed reduced biofilm production related to lower levels of the extracellular matrix, and the mutant also had reduced sporulation, adhesion, root colonization, and FCR biocontrol efficiency. Transcriptome analysis of YB-1471 and Δorf04391 in stationary culture showed that the loss of orf04391 resulted in altered expression of numerous genes, including sinI, an initiator of biofilm formation. DNA binding was shown with his-tagged Orf04391 binding to the sinIR operon in vivo and in vitro. Orf04391 appears to be a transcriptional regulator of biofilm formation in B. subtilis through the Spo0A-SinI/SinR pathway.

De novo engineering of programmable and multi-functional biomolecular condensates for controlled biosynthesis.

There is a growing interest in the creation of engineered condensates formed via liquid-liquid phase separation (LLPS) to exert precise cellular control in prokaryotes. However, de novo design of cellular condensates to control metabolic flux or protein translation remains a challenge. Here, we present a synthetic condensate platform, generated through the incorporation of artificial, disordered proteins to realize specific functions in Bacillus subtilis. To achieve this, the "stacking blocks" strategy is developed to rationally design a series of LLPS-promoting proteins for programming condensates. Through the targeted recruitment of biomolecules, our investigation demonstrates that cellular condensates effectively sequester biosynthetic pathways. We successfully harness this capability to enhance the biosynthesis of 2'-fucosyllactose by 123.3%. Furthermore, we find that condensates can enhance the translation specificity of tailored enzyme fourfold, and can increase N-acetylmannosamine titer by 75.0%. Collectively, these results lay the foundation for the design of engineered condensates endowed with multifunctional capacities.

Integrating Deep Learning and Synthetic Biology: A Co-Design Approach for Enhancing Gene Expression via N-Terminal Coding Sequences.

N-terminal coding sequence (NCS) influences gene expression by impacting the translation initiation rate. The NCS optimization problem is to find an NCS that maximizes gene expression. The problem is important in genetic engineering. However, current methods for NCS optimization such as rational design and statistics-guided approaches are labor-intensive yield only relatively small improvements. This paper introduces a deep learning/synthetic biology codesigned few-shot training workflow for NCS optimization. Our method utilizes k-nearest encoding followed by word2vec to encode the NCS, then performs feature extraction using attention mechanisms, before constructing a time-series network for predicting gene expression intensity, and finally a direct search algorithm identifies the optimal NCS with limited training data. We took green fluorescent protein (GFP) expressed by Bacillus subtilis as a reporting protein of NCSs, and employed the fluorescence enhancement factor as the metric of NCS optimization. Within just six iterative experiments, our model generated an NCS (MLD62) that increased average GFP expression by 5.41-fold, outperforming the state-of-the-art NCS designs. Extending our findings beyond GFP, we showed that our engineered NCS (MLD62) can effectively boost the production of N-acetylneuraminic acid by enhancing the expression of the crucial rate-limiting GNA1 gene, demonstrating its practical utility. We have open-sourced our NCS expression database and experimental procedures for public use.

Impact of a Potent Strain of Plant Growth-Promoting Bacteria (PGPB), Bacillus subtilis S1 on Bacterial Community Composition, Enzymatic Activity, and Nitrogen Content in Cucumber Rhizosphere Soils.

Antagonistic bacterial strains from Bacillus spp. have been widely studied and utilized in the biocontrol of phytopathogens and the promotion of plant growth, but their impacts on the rhizosphere microecology when applied to crop plants are unclear. Herein, the effects of applying the antagonistic bacterium Bacillus subtilis S1 as a biofertilizer on the rhizosphere microecology of cucumbers were investigated. In a pot experiment on cucumber seedlings inoculated with S1, 3124 bacterial operational taxonomic units (OTUs) were obtained from the rhizosphere soils using high-throughput sequencing of 16S rRNA gene amplicons, and the most abundant phylum was Proteobacteria that accounted for 49.48% in the bacterial community. S1 treatment significantly reduced the abundances of soil bacterial taxa during a period of approximately 30 days but did not affect bacterial diversity in the rhizosphere soils of cucumbers. The enzymatic activities of soil nitrite reductase (S-Nir) and dehydrogenase (S-DHA) were significantly increased after S1 fertilization. However, the activities of soil urease (S-UE), cellulase (S-CL), and sucrase (S-SC) were significantly reduced compared to the control group. Additionally, the ammonium- and nitrate-nitrogen contents of S1-treated soil samples were significantly lower than those of the control group. S1 fertilization reshaped the rhizosphere soil bacterial community of cucumber plants. The S-CL activity and nitrate-nitrogen content in rhizosphere soil affected by S1 inoculation play important roles in altering the abundance of rhizosphere soil microbiota.

Control of bacterial cell wall autolysins by peptidoglycan crosslinking mode.

To withstand their internal turgor pressure and external threats, most bacteria have a protective peptidoglycan (PG) cell wall. The growth of this PG polymer relies on autolysins, enzymes that create space within the structure. Despite extensive research, the regulatory mechanisms governing these PG-degrading enzymes remain poorly understood. Here, we unveil a novel and widespread control mechanism of lytic transglycosylases (LTs), a type of autolysin responsible for breaking down PG glycan chains. Specifically, we show that LD-crosslinks within the PG sacculus act as an inhibitor of LT activity. Moreover, we demonstrate that this regulation controls the release of immunogenic PG fragments and provides resistance against predatory LTs of both bacterial and viral origin. Our findings address a critical gap in understanding the physiological role of the LD-crosslinking mode in PG homeostasis, highlighting how bacteria can enhance their resilience against environmental threats, including phage attacks, through a single structural PG modification.

Modelling protein complexes with crosslinking mass spectrometry and deep learning.

Scarcity of structural and evolutionary information on protein complexes poses a challenge to deep learning-based structure modelling. We integrate experimental distance restraints obtained by crosslinking mass spectrometry (MS) into AlphaFold-Multimer, by extending AlphaLink to protein complexes. Integrating crosslinking MS data substantially improves modelling performance on challenging targets, by helping to identify interfaces, focusing sampling, and improving model selection. This extends to single crosslinks from whole-cell crosslinking MS, opening the possibility of whole-cell structural investigations driven by experimental data. We demonstrate this by revealing the molecular basis of iron homoeostasis in Bacillus subtilis.

From dusty shelves toward the spotlight: growing evidence for Ap4A as an alarmone in maintaining RNA stability and proteostasis.

Bacteria thrive in diverse environments and must withstand various stresses. A key stress response mechanism is the reprogramming of macromolecular biosynthesis and metabolic processes through alarmones - signaling nucleotides that accumulate intracellularly in response to metabolic stress. Diadenosine tetraphosphate (Ap4A), a putative alarmone, is produced in a noncanonical reaction by universally conserved aminoacyl-tRNA synthetases. Ap4A is ubiquitous across all domains of life and accumulates during heat and oxidative stress. Despite its early discovery in 1966, Ap4A's alarmone status remained inconclusive. Recent discoveries identified Ap4A as a precursor to RNA 5' caps in Escherichia coli. Additionally, Ap4A was found to directly bind to and allosterically inhibit the purine biosynthesis enzyme inosine 5'-monophosphate dehydrogenase, regulating guanosine triphosphate levels and enabling heat resistance in Bacillus subtilis. These findings, along with previous research, strongly suggest that Ap4A plays a crucial role as an alarmone, warranting further investigation to fully elucidate its functions.

Bacillus subtilis JATP3 improved the immunity of weaned piglets by improving intestinal flora and producing citalopram.

The purpose of this study was to evaluate the ability of Bacillus subtilis JATP3 to stimulate immune response and improve intestinal health in piglets during the critical weaning period. Twelve 28-day-old weaned piglets were randomly divided into two groups. One group was fed a basal diet, while the other group was fed a basal diet supplemented with B. subtilis JATP3 (1 × 109 CFU/mL; 10 mL) for 28 days. The results revealed a significant increase in the intestinal villus gland ratio of weaned piglets following the inclusion of B. subtilis JATP3 (P < 0.05). Inclusion of a probiotic supplement improve the intestinal flora of jejunum and ileum of weaned piglets. Metabolomics analysis demonstrated a notable rise in citalopram levels in the jejunum and ileum, along with elevated levels of isobutyric acid and isocitric acid in the ileum. The results of correlation analysis show that indicated a positive correlation between citalopram and microbial changes. Furthermore, the probiotic-treated group exhibited a significant upregulation in the relative expression of Claudin, Zonula Occludens 1 (ZO-1), and Interleukin 10 (IL-10) in the jejunum and ileum, while displaying a noteworthy reduction in the relative expression of Interleukin 1ß (IL-1ß). Overall, these findings suggest that B. subtilis JATP3 can safeguard intestinal health by modulating the structure of the intestinal microbiota and their metabolites, wherein citalopram might be a key component contributing to the therapeutic effects of B. subtilis JATP3.

Membrane Localization of RNase Y Is Important for Global Gene Expression in Bacillus subtilis.

RNase Y is a key endoribonuclease that regulates global mRNA turnover and processing in Bacillus subtilis and likely many other bacteria. This enzyme is anchored to the cell membrane, creating a pseudo-compartmentalization that aligns with its role in initiating the decay of mRNAs primarily translated at the cell periphery. However, the reasons behind and the consequences of RNase Y's membrane attachment remain largely unknown. In our study, we examined a strain expressing wild-type levels of a cytoplasmic form of RNase Y from its chromosomal locus. This strain exhibits a slow-growth phenotype, similar to that of an RNase Y null mutant. Genome-wide data reveal a significant impact on the expression of hundreds of genes. While certain RNA substrates clearly depend on RNase Y's membrane attachment, others do not. We observed no correlation between mRNA stabilization in the mutant strains and the cellular location or function of the encoded proteins. Interestingly, the Y-complex, a specificity factor for RNase Y, also appears also recognize the cytoplasmic form of the enzyme, restoring wild-type levels of the corresponding transcripts. We propose that membrane attachment of RNase Y is crucial for its functional interaction with many coding and non-coding RNAs, limiting the cleavage of specific substrates, and potentially avoiding unfavorable competition with other ribonucleases like RNase J, which shares a similar evolutionarily conserved cleavage specificity.

High-level and -yield orotic acid production in Escherichia coli through systematic modular engineering and "Chaos to Order Cycles" fermentation.

Orotic acid is widely used in healthcare and cosmetic industries. However, orotic acid-producing microorganisms are auxotrophic, which results in inefficient microbial production. Herein, a plasmid-free, uninduced, non-auxotrophic orotic acid hyperproducer was constructed from Escherichia coli W3110. Initially, the orotic acid degradation pathway was blocked and the carbamoyl phosphate supply was enriched. Subsequently, pyr operon from Bacillus subtilis F126 was heterologously expressed and precursors' supply was optimized. Thereafter, pyrE was dynamically regulated to reconstruct the non-auxotrophic pathway. Employing fed-batch cultivation, orotic acid titer, yield, and productivity of strain Ora21 reached 182.5 g/L, 0.58 g/g, and 3.80 g/L/h, respectively, the highest levels reported so far. Finally, a novel "Chaos to Order Cycles (COC)" fermentation was developed, which effectively increased the yield to 0.63 g/g. This research is a remarkable achievement in orotic acid production by microbial fermentation and has vast potential for industrial applications.

Efficient Production of an Alginate Lyase in Bacillus subtilis with Combined Strategy: Vector and Host Selection, Promoter and Signal Peptide Screening, and Modification of a Translation Initiation Region.

Alginate lyases (ALys) whose degrading products, alginate oligosaccharides, exhibit various outstanding biochemical activities have aroused increasing interest of researchers in the marine bioresource field. However, their predominant sourcing from marine bacteria, with limited yields and unclear genetic backgrounds, presents a challenge for industrial production. In this study, ALys (Aly01) from Vibrio natriegens SK 42.001 was expressed in Bacillus subtilis (B. subtilis), a nonpathogenic microorganism recognized as generally safe (GRAS). This accomplishment was realized through a comprehensive strategy involving vector and host selection, promoter and signal peptide screening, and engineering of the ribosome binding site (RBS) and the N-terminal coding sequence (NCS). The optimal combination was identified as the pP43NMK and B. subtilis WB600. Among the 19 reported strong promoters, PnprE exhibited the best performance, showing intracellular enzyme activities of 4.47 U/mL. Despite expectations, dual promoter construction did not yield a significant increase. Further, SPydhT demonstrated the highest extracellular activity (1.33 U/mL), which was further improved by RBS/NCS engineering, reaching 4.58 U/mL. Finally, after fed-batch fermentation, the extracellular activity reached 18.01 U/mL, which was the highest of ALys with a high molecular weight expressed in B. subtilis. These findings are expected to offer valuable insights into the heterologous expression of ALys in B. subtilis.

Temperature-Dependent tRNA Modifications in Bacillales.

Transfer RNA (tRNA) modifications are essential for the temperature adaptation of thermophilic and psychrophilic organisms as they control the rigidity and flexibility of transcripts. To further understand how specific tRNA modifications are adjusted to maintain functionality in response to temperature fluctuations, we investigated whether tRNA modifications represent an adaptation of bacteria to different growth temperatures (minimal, optimal, and maximal), focusing on closely related psychrophilic (P. halocryophilus and E. sibiricum), mesophilic (B. subtilis), and thermophilic (G. stearothermophilus) Bacillales. Utilizing an RNA sequencing approach combined with chemical pre-treatment of tRNA samples, we systematically profiled dihydrouridine (D), 4-thiouridine (s4U), 7-methyl-guanosine (m7G), and pseudouridine (Ψ) modifications at single-nucleotide resolution. Despite their close relationship, each bacterium exhibited a unique tRNA modification profile. Our findings revealed increased tRNA modifications in the thermophilic bacterium at its optimal growth temperature, particularly showing elevated levels of s4U8 and Ψ55 modifications compared to non-thermophilic bacteria, indicating a temperature-dependent regulation that may contribute to thermotolerance. Furthermore, we observed higher levels of D modifications in psychrophilic and mesophilic bacteria, indicating an adaptive strategy for cold environments by enhancing local flexibility in tRNAs. Our method demonstrated high effectiveness in identifying tRNA modifications compared to an established tool, highlighting its potential for precise tRNA profiling studies.

Proteomic profiles revealed enzymatic activities associated with the flavor formation of salted shrimp paste influenced by Bacillus subtilis K-C3 inoculation.

Enzymatic proteomic profiles were examined to comprehend the predominant enzymes involved in the flavor development of salted shrimp paste influenced by Bacillus subtilis K-C3 inoculation (Inoc), compared to those without inoculation (CON). Inoc showed greater proteolytic, lipolytic, and chitinolytic activities than CON (P < 0.05) throughout 30 days of fermentation, indicating B. subtilis's ability to accelerate the fermentation rate and render distinctive flavor profiles to shrimp paste. Among 50 differential abundance proteins (DAPs), 24 DAPs were identified as potential key regulating enzymes, with a P-value < 0.05 and |FC| > 0.50, indicating their significance and regulating capacity within specific metabolic pathways. Notably, 27 and 23 DAPs were up-regulated in Inoc and CON, respectively. Moreover, gene ontology (GO) enrichment analysis revealed that hydrolases, involved in carbohydrate metabolic processes and proteolysis, were the most differentiating pathways between Inoc and CON. Both samples exhibited different flavor profiles. A greater abundance of N-containing volatile compounds with a lower total abundance of aldehydes, ketones, alcohols, and acids could suggest a more favorable flavor in Inoc, compared to CON. Principal component analysis (PCA) revealed a positive correlation between L-ascorbate peroxidase, carboxypeptidase, and tripeptidyl peptidase sed2, with proteolytic and lipolytic activities in Inoc (P < 0.05). Meanwhile, acids and alcohols were positively correlated with CON. Therefore, B. subtilis inoculation could produce a distinctive flavor with a desirable sensory perception of shrimp paste regarding its ability to release extracellular enzymes/proteins. B. subtilis K-C3 inoculation could be suggested in the production of shrimp paste to improve its flavor characteristics.