High-yield tannase strain screening method based on genome engineering

Abstract

The invention provides a high-yield tannase strain screening method based on genome engineering, which comprises the following steps: acquiring a signal peptide sequence set of tannase, and performing sequence analysis on amino acid composition of a hydrophobic core region to form a signal peptide conformation feature data set; if the identification efficiency score of the transport channel is higher than a preset threshold value, performing fusion expression on the corresponding hydrophobic core region sequence and a tannase gene to obtain a recombinant expression vector; optimizing the amino acid composition of the hydrophobic core region according to the signal peptide sequence of which the exposure degree of the cleavage site is higher than a threshold value, and generating an optimized signal peptide sequence set; an optimized signal peptide sequence set is adopted, a recombinant expression vector is reconstructed, and a fermentation process formula is improved; and determining a final tannase strain production scheme according to the improved fermentation process formula.

Description

Technical Field

[0001] The present invention relates to the field of information technology, and in particular to a method for screening high-yield tannase strains based on genome engineering. Background Art

[0002] The application of genome engineering in industrial enzyme production is crucial, particularly in optimizing microbial strains for high-yield secretory expression of tannases. This is crucial for cost control and efficiency improvements in industries such as food, feed, and pharmaceuticals. Genetically engineering the secretory pathway of a strain can significantly increase the extracellular production of the target enzyme and reduce fermentation and purification costs. However, existing signal peptide optimization methods typically rely on single-sequence design or low-throughput screening, making it difficult to systematically address the complex regulatory mechanisms of secretion efficiency. Furthermore, a lack of in-depth understanding of the dynamic effects of signal peptide conformational changes has resulted in limited improvements in secretion efficiency, making it difficult to meet the demands of industrial production. As key components in regulating protein secretion, the amino acid composition of the hydrophobic core region of signal peptides directly influences secretion efficiency. However, current research faces three key challenges: First, the mechanism by which signal peptide conformational changes influence the recognition efficiency of transport channels on the endoplasmic reticulum membrane remains unclear, limiting the design of efficient signal peptides. Second, the accumulation of folding intermediates in the secretory pathway varies with signal peptide sequence, potentially leading to protein misfolding or degradation. Finally, the ease of exposure of the signal peptide cleavage site influences enzyme maturation and release, but current methods struggle to precisely control this process. These challenges have resulted in low efficiency in the construction and screening of high-throughput signal peptide libraries in practical applications, making it difficult to increase the extracellular secretion of tannase and reduce purification costs. Therefore, optimizing the amino acid composition of the hydrophobic core region of the signal peptide to regulate the effects of conformational changes on transport channel recognition, folding intermediate accumulation, and cleavage site exposure, combined with the screening technology of high-throughput signal peptide libraries and tannase gene fusion expression, has become a key issue in improving the extracellular secretion of tannase and reducing purification costs in fermentation processes. Summary of the Invention

[0003] The present invention provides a method for screening high-yield tannase strains based on genome engineering, which mainly comprises:

[0004] Obtain a collection of tannase signal peptide sequences and perform sequence analysis on the amino acid composition of the hydrophobic core region to form a signal peptide conformational feature dataset;

[0005] Based on the signal peptide conformational feature dataset, we analyzed the effects of changes in the amino acid composition of the hydrophobic core region on the conformational stability of the signal peptide and the changes in the binding affinity of the signal peptide to the transport channel on the endoplasmic reticulum membrane. We then determined the transport channel recognition efficiency score by combining the analysis results of conformational stability and binding affinity changes.

[0006] If the transport channel recognition efficiency score is higher than the preset threshold, the corresponding hydrophobic core region sequence is fused with the tannase gene to obtain a recombinant expression vector;

[0007] Fermentation experiments were conducted based on a high-throughput signal peptide library containing recombinant expression vectors to detect the amount of tannase secreted outside the cells, and the accumulation of folding intermediates was assessed by analyzing the intracellular content of tannase precursor protein.

[0008] The relationship between the accumulation level of intracellular folding intermediates and the secretion efficiency of the extracellular mature enzyme was analyzed. If the accumulation level of intracellular intermediates was lower than the threshold, the cleavage site exposure was calculated based on the solvent-accessible surface area of the signal peptide, and the cleavage efficiency score was generated based on the enzyme cleavage time in the in vitro cleavage experiment.

[0009] Optimize the amino acid composition of the hydrophobic core region based on the signal peptide sequence whose cleavage site exposure is higher than the threshold, and generate an optimized signal peptide sequence set;

[0010] Using an optimized signal peptide sequence set, reconstructing the recombinant expression vector, and improving the fermentation process formula;

[0011] The final tannase strain production plan was determined based on the improved fermentation process formula.

[0012] The technical solution provided by the embodiment of the present invention may have the following beneficial effects:

[0013] The present invention discloses a method for screening high-yield tannase strains based on genome engineering. By analyzing the amino acid composition and conformational characteristics of the hydrophobic core region of the signal peptide sequence and its effect on the recognition efficiency of the transport channel, a high-efficiency signal peptide sequence is screened. The screened signal peptide sequence is fused with the tannase gene for expression, a recombinant expression vector is constructed, and a high-throughput fermentation experiment is carried out. By analyzing the relationship between the accumulation level of intracellular folding intermediates and the secretion efficiency of the extracellular mature enzyme, the cleavage efficiency of the signal peptide sequence is optimized. Finally, the present invention uses the optimized signal peptide sequence to reconstruct the expression vector, optimizes the fermentation conditions through the response surface analysis method, and achieves efficient secretory expression of tannase. This method can significantly increase the yield of tannase, reduce purification costs, and provide reliable technical support for industrial production. BRIEF DESCRIPTION OF THE DRAWINGS

[0014] Figure 1 The present invention is a flow chart of a method for screening high-yielding tannase strains based on genome engineering. DETAILED DESCRIPTION

[0015] In order to make the objectives, technical solutions and advantages of the present invention more clear, the present invention is described in detail below with reference to the accompanying drawings and specific embodiments.

[0016] like Figure 1In this embodiment, a method for screening high-yielding tannase strains based on genome engineering may specifically include:

[0017] S101. Obtain a signal peptide sequence set of tannase, and perform sequence analysis on the amino acid composition of the hydrophobic core region to form a signal peptide conformation feature dataset.

[0018] A set of tannase signal peptide sequences was obtained according to a preset genome database, and the signal peptide sequence set was subjected to amino acid residue coding standardization to obtain a standardized signal peptide sequence; for the standardized signal peptide sequence, the amino acid residues were scored for hydrophobicity, the local average hydrophobicity value was calculated, and an amino acid hydrophobicity score matrix was obtained; a hydrophobicity threshold was set according to the amino acid hydrophobicity score matrix, the boundary of the hydrophobic core region was located, the types and occurrence frequencies of amino acid residues in the hydrophobic core region were counted, and a hydrophobic core region amino acid composition feature vector was generated; for the standardized signal peptide sequence, the three conformational probability values of α-helix, β-fold, and random coil of the amino acid residues were calculated to obtain a secondary structure prediction probability matrix; based on the secondary structure prediction probability matrix, the φ-ψ dihedral angle values of the amino acid residues in each conformational segment were calculated, and the stability score of each conformational segment was scored in combination with the amino acid residue tendency parameter; the hydrophobic core region amino acid composition feature vector, the secondary structure prediction probability matrix, and the conformational segment stability score were integrated to construct a signal peptide sequence feature data set.

[0019] For example, based on the tannase sequence in a preset genomic database, a protein sequence alignment algorithm was used to obtain a set of tannase signal peptide sequences. The sequence set was then standardized for amino acid residue encoding to generate a standardized signal peptide sequence dataset. For the standardized signal peptide sequence dataset, the hydrophobicity of the amino acid residues in the entire sequence was scored using the Cait-Doolittle hydrophobicity calculation method. The local average hydrophobicity values were calculated using a heptapeptide sliding window to obtain an amino acid hydrophobicity score matrix. Based on the amino acid hydrophobicity score matrix, a hydrophobicity threshold was set to locate the boundaries of the hydrophobic core region. The types and frequencies of amino acid residues in the hydrophobic core region were counted to generate a hydrophobic core region amino acid composition feature vector. For the standardized signal peptide sequence, the Zhu-Deng secondary structure prediction algorithm was used to calculate the probability values of the three conformations of α-helix, β-sheet, and random coil for each amino acid residue in the sequence to obtain a secondary structure prediction probability matrix. Based on the secondary structure prediction probability matrix, the φ-ψ dihedral angle values of the amino acid residues in each conformational segment were calculated, and the stability score of each conformational segment was calculated in combination with the amino acid residue propensity parameters. A signal peptide sequence feature dataset was constructed by integrating the hydrophobic core region amino acid composition feature vector, secondary structure prediction probability matrix, and conformational segment stability score. Sequence analysis of tannase signal peptides began with data preprocessing. The amino acid residues were standardized, using single-letter codes instead of three-letter codes, such as ALA for alanine and LEU for leucine, to ensure consistent sequence representation. Based on the standardized sequences, the Kayt-Doolittle hydrophobicity calculation method was applied, assigning different hydrophobicity parameters to different amino acid residues. For example, isoleucine was assigned a value of 4.5, valine was assigned a value of 4.2, leucine was assigned a value of 3.8, and phenylalanine was assigned a value of 2.8. When the local average hydrophobicity values were calculated using a heptapeptide sliding window, the sequence ILVFAML had a hydrophobicity score of 3.6, while the sequence RKDENG had a hydrophobicity score of -2.8, reflecting the differences in hydrophobic properties between different segments. To locate the hydrophobic core region, a hydrophobicity threshold of 2.5 was set, and the sequence was scanned. A hydrophobic core region was identified when the average hydrophobicity score of seven consecutive amino acid residues exceeded the threshold. Amino acid composition analysis of the hydrophobic core region revealed that hydrophobic amino acids such as isoleucine, leucine, valine, and phenylalanine were significantly more frequent than other amino acids, with isoleucine accounting for 28%, leucine for 24%, and valine for 18%. For secondary structure prediction, the Zhu-Deng algorithm calculates the probability of each site forming an α-helix, β-sheet, or random coil based on the physicochemical properties of the amino acid residues and the influence of neighboring residues. For example, the probability of forming an α-helix at position 4 in the sequence AIAALLV is 0.82, the probability of forming a β-sheet is 0.12, and the probability of forming a random coil is 0.06, indicating that this site favors an α-helical conformation. Conformational stability is scored based on the φ-ψ dihedral angles of the amino acid residues. In an α-helical conformation, the φ angle is approximately -57 degrees, and the ψ angle is approximately -47 degrees.The stability scores of the conformational segments were calculated by combining conformational propensity parameters of amino acid residues, such as the α-helical propensity of alanine (1.42) and leucine (1.21). The hydrophobic core region characteristics, secondary structure predictions, and conformational stability scores were integrated to construct a multidimensional dataset for subsequent functional analysis of signal peptides. In practical applications, these characteristic data can reflect the structural characteristics and functional properties of signal peptide sequences, providing data support for studies of protein secretion and targeted trafficking mechanisms.

[0020] S102. Based on the signal peptide conformational feature dataset, analyze the effects of changes in the amino acid composition of the hydrophobic core region on the conformational stability of the signal peptide and the changes in the binding affinity of the signal peptide to the transport channel on the endoplasmic reticulum membrane. Determine the transport channel recognition efficiency score by combining the results of the conformational stability and binding affinity changes.

[0021] The hydrophobic core region amino acid sequence is obtained from the signal peptide conformational feature data set, the variation sites in the hydrophobic core region amino acid sequence are identified, and the difference values of the amino acid physicochemical properties are obtained; the protein conformational parameters are calculated using the molecular dynamics method based on the difference values of the amino acid physicochemical properties, and a conformational stability feature matrix is generated; based on the conformational stability feature matrix, the spatial matching degree of the binding interface between the hydrophobic core region of the signal peptide and the endoplasmic reticulum membrane transport channel is calculated to obtain a binding site coordinate set; based on the binding site coordinate set, the interaction energy value at the binding site is obtained by the molecular force field calculation method to generate a binding affinity feature vector; based on the conformational stability feature matrix and the binding affinity feature vector, a multi-dimensional data fusion method is used to calculate the transport channel recognition efficiency score.

[0022] For example, for a signal peptide conformational feature dataset, the amino acid sequence of the hydrophobic core region is extracted, and the amino acid variation sites are identified using a sequence alignment method, and the difference in amino acid physicochemical properties before and after the variation is calculated. Based on the difference in amino acid physicochemical properties, the rate of change of the protein main chain dihedral angle, the rate of change of the side chain conformation, and the rate of change of the number of intramolecular hydrogen bonds are calculated using molecular dynamics methods to generate a conformational stability feature matrix. Based on the conformational stability feature matrix, a structural superposition algorithm is used to calculate the spatial matching degree of the binding interface between the hydrophobic core region of the signal peptide and the endoplasmic reticulum membrane transport channel to obtain a set of binding site coordinates. Based on the binding site coordinate set, the hydrogen bond energy value, hydrophobic interaction energy value, and van der Waals interaction energy value at each binding site are obtained using a molecular force field calculation method to generate a binding affinity feature vector. Based on the conformational stability feature matrix and the binding affinity feature vector, a multi-dimensional data fusion method is used to calculate the interaction area, site exposure, and steric hindrance parameters between the signal peptide and the transport channel to generate a transport channel recognition efficiency score. Analysis of the amino acid composition of the hydrophobic core region of the signal peptide began with sequence alignment. Multiple sequence alignment algorithms were used to identify amino acid variants. For example, when leucine was replaced with alanine at position 15, the calculated differences in physicochemical properties revealed a 2.6-unit decrease in hydrophobicity and a 0.8-unit decrease in β-sheet tendency. This change in amino acid composition directly impacts the structural stability of the signal peptide. Molecular dynamics simulations, using 300-nanosecond trajectories, revealed that the fluctuation range of the main chain dihedral angle φ in the hydrophobic core region increased from ±10 degrees to ±25 degrees, indicating a significant increase in structural flexibility. Conversely, the number of intramolecular hydrogen bonds decreased from 12 to 8, with the presence of four key stabilizing hydrogen bonds decreasing from 85% to 45% of the time, indicating a significant decrease in conformational stability. Structural superposition analysis revealed a significant change in the spatial fit between the signal peptide's hydrophobic core region and the binding interface of the transport channel, with the average interatomic distance increasing from 3.8 angstroms to 4.5 angstroms. The spatial coordinates of the key binding site revealed a 2.2-angstrom gap in the previously tightly contacted hydrophobic pocket, resulting in a 28% reduction in the hydrophobic interaction area. Molecular force field calculations revealed a decrease in the average energy of individual hydrogen bonds from -5.8 kcal / mol to -3.2 kcal / mol, and a decrease in the total hydrophobic interaction energy from -42.5 kcal / mol to -28.6 kcal / mol. Regarding van der Waals interactions, the interaction energy decreased from -15.4 kcal / mol to -8.8 kcal / mol due to increased steric hindrance. Data fusion analysis revealed a correlation between conformational changes and weakened binding. The interaction interface area decreased from 850 square angstroms to 620 square angstroms, and the relative solvent-accessible area of the key recognition site increased from 15% to 35%, indicating a significant increase in binding site exposure.At the same time, steric hindrance parameters revealed that the complementarity at the interface between the signal peptide and the transport channel decreased, while van der Waals repulsion increased 2.6-fold, resulting in a decrease in the transport channel recognition efficiency score from 0.82 to 0.45. This multi-dimensional data analysis reveals how changes in amino acid composition alter the interaction pattern between the signal peptide and the transport channel by affecting conformational stability. In practical applications, these parameter changes can predict the transport efficiency of signal peptides, providing data support for the optimized design of protein secretion signal peptides.

[0023] S103. If the transport channel recognition efficiency score is higher than a preset threshold, the corresponding hydrophobic core region sequence is fused with the tannase gene to obtain a recombinant expression vector.

[0024] Obtain a transport channel recognition efficiency score, and determine whether it exceeds a preset threshold based on the score. If it exceeds the preset threshold, derive a hydrophobic core region sequence from a signal peptide conformational feature dataset and generate a target DNA sequence; based on the target DNA sequence, obtain a tannase gene sequence from a genome database using a sequence homology alignment tool, and determine a sequence splicing site for the tannase gene sequence; based on the splicing site, perform a directionally connected connection between the target DNA sequence and the tannase gene sequence using a gene splicing tool to obtain a fusion gene sequence; based on the fusion gene sequence, perform sequence optimization on the codon composition, select a matching expression regulatory element for the optimized fusion gene sequence, and construct a recombinant expression vector using a vector design tool.

[0025] For example, a numerical comparator is used to determine whether the transport channel recognition efficiency score exceeds a preset threshold. The corresponding hydrophobic core region sequence is derived from a signal peptide conformational feature dataset, and a nucleotide sequence encoding conversion tool is used to generate a target DNA sequence. Based on the target DNA sequence, a sequence homology alignment tool is used to retrieve the tannase gene sequence from a genomic database, and a sequence analysis tool is used to determine suitable sequence splicing sites. Based on the identified splicing sites, a gene splicing tool is used to perform a directionally linked ligation of the target DNA sequence with the tannase gene sequence to generate a fusion gene sequence. Based on the fusion gene sequence, a gene expression prediction tool is used to analyze its open reading frame integrity, codon usage frequency, and expression level, and a sequence optimization algorithm is used to adjust the codon composition. Based on the optimized fusion gene sequence, a promoter screening tool is used to select matching expression regulatory elements, and a recombinant expression vector containing the complete expression unit is constructed using a vector design tool. Based on the recombinant expression vector sequence, a nucleic acid sequence analysis tool is used to verify the insertion direction, reading frame integrity, and correctness of the expression regulatory elements in the fusion gene. The transport pathway recognition efficiency score was set at a preset threshold of 0.75; a score exceeding this value indicates a good match between the signal peptide and the transport pathway. A typical hydrophobic core region sequence, LVFAMLIL, was subjected to nucleotide code conversion to generate the target DNA sequence, CTGGTGTTTGCTATGCTGATCCTG. This conversion maintains codon translation efficiency. Sequence homology comparison revealed 35% similarity between the target DNA sequence and the tannase gene sequence in the N-terminal leader region, with high sequence conservation from nucleotides 15 to 45. Analysis of the secondary structure of this region revealed a natural breakpoint after nucleotide 42, with a GC content of 48%, making it a suitable splicing site. At the identified splicing site, the target DNA sequence was fused to the tannase gene sequence using a directed ligation strategy. The ligation reaction introduced a GAATTC recognition site at the 5' end and a GGATCC recognition site at the 3' end, forming the directed ligation fusion gene sequence. Sequence analysis revealed uniform nucleotide composition before and after the fusion site, with no significant GC content mutations. Gene expression prediction indicated an open reading frame integrity score of 0.92 and a codon adaptation index of 0.85 for the fusion gene sequence. An optimization algorithm was used to adjust rare codons, such as replacing AGG with CGT, to align codon usage frequency with host preference. The predicted expression level of the optimized sequence increased by 25%. During promoter screening, a strong promoter sequence with a transcriptional activity of 2.6 units was selected, with sequence integrity scores of 0.88 and 0.92 for the -35 and -10 regions, respectively. The spacer between this promoter and the ribosome binding site is 8 nucleotides, falling within the optimal range.Through vector design, the fusion gene sequence was inserted 15 nucleotides downstream of the promoter to construct a complete expression unit. Vector verification results showed that the fusion gene was inserted into the vector in the forward orientation and in the correct orientation relative to the promoter. Open reading frame analysis confirmed that there were no frameshift mutations between the start codon ATG and the stop codon TAG, and the full sequence length was 1236 nucleotides. The transcription terminator has a standard stem-loop structure, and its predicted free energy value is -12.5 kcal / mol, indicating a strong transcription termination signal. These parameters verified the correct construction of the recombinant expression vector.

[0026] S104. Fermentation experiments were conducted based on a high-throughput signal peptide library containing a recombinant expression vector to detect the amount of tannase secreted outside the cell, and the accumulation level of folding intermediates was evaluated by analyzing the intracellular tannase precursor protein content.

[0027] A site-directed mutagenesis tool is used to obtain amino acid replacement data at conserved sites in the signal peptide sequence, and a mutation site combination matrix is generated based on the amino acid replacement data; a recombinant expression vector is constructed based on the mutation site combination matrix using a gene synthesis tool, and the recombinant expression vector is transformed to obtain an expression strain library; for the expression strain library, a fermentation monitoring device is used to record dissolved oxygen concentration and cell density data, and a fermentation culture supernatant is obtained through a bioreactor control device; for the fermentation culture supernatant, secretory tannase content data is measured to obtain protein secretion efficiency; based on the protein secretion efficiency, proteins suitable for predetermined requirements are screened for cell disruption, intracellular protein extracts are obtained, and the accumulation level of folding intermediates is evaluated.

[0028] Exemplarily, based on the recombinant expression vector sequence, a site-directed mutagenesis tool is used to replace amino acids at conserved sites in the signal peptide sequence to generate a mutation site combination matrix. Based on the mutation site combination matrix, multiple recombinant expression vectors are constructed using gene synthesis tools, and an expression strain library is obtained using an efficient transformation method. Based on the expression strain library, a fermentation monitoring device is used to record the dissolved oxygen concentration, cell density, and culture medium pH during the culture process, and the specified fermentation conditions are maintained by a bioreactor control device. For the fermentation culture supernatant, an activity detection method is used to determine the secretory tannase content, and the protein secretion efficiency is calculated by a fluorescence quantitative method. Based on the cell density data, a cell separation method is used to collect the fermentation bacteria, and an intracellular protein extract is obtained by an ultrasonic disruption method. For the intracellular protein extract, a protein separation method is used to obtain the tannase precursor protein, and the folding intermediate content is determined by a structure detection method. Based on the secretory tannase content and the precursor protein content data, the signal peptide-mediated protein secretion and transport efficiency is evaluated by a correlation calculation method. Site-directed mutagenesis of the signal peptide sequence began with the selection of conserved sites. Sequence alignment identified key amino acid positions. For example, leucine at position 8 in the hydrophobic core region is 95% conserved, and mutations in this position significantly affect signal peptide function. Amino acid substitutions were made at this site, replacing leucine with alanine, valine, or isoleucine, generating a comprehensive dataset of mutation sites. During gene synthesis, codon selection at each mutation site took into account host preference. For example, in Escherichia coli, leucine prefers the CTG codon, with a usage frequency of 50%. By optimizing codon combinations, an expression vector library containing 128 mutants was constructed, and the resulting transformed strains achieved an average transformation efficiency of 5×10⁶ transformants per microgram of DNA. During fermentation monitoring, dissolved oxygen concentration was maintained at 40% saturation, cell density, measured by optical density, increased from an initial 0.2 to 4.8, and culture medium pH fluctuated between 7.2 and 6.8. These parameters reflect cell growth and metabolic activity, directly influencing protein expression and secretion. Tannase activity in the culture supernatant was measured using the substrate conversion method, with the unit of enzyme activity defined as the amount of enzyme required to convert 1 micromole of substrate per minute. Fluorescence quantitative analysis revealed that the proportion of secreted tannase in the total protein increased from 12% to 28%, indicating a significant improvement in signal peptide-mediated protein secretion efficiency. Proteins with significantly improved secretion efficiency were selected for cell disruption, with the specific increase set at 16%. Cell disruption was performed using ultrasonication, with 20 cycles of 15 seconds of sonication at 25% power and 45 seconds of rest at 4°C. The resulting intracellular protein extract had a total protein concentration of 2.8 mg / mL. The tannase precursor protein, isolated by immunoaffinity chromatography, was over 90% pure. Structural analysis revealed the presence of a partial folding intermediate within the precursor protein, characterized by an increased exposure of hydrophobic regions, leading to an increased propensity for protein aggregation.Circular dichroism analysis revealed a decrease in α-helix content from 25% to 18% and an increase in β-sheet content from 35% to 42%, indicating significant changes in protein conformation. Correlation analysis revealed that as secreted tannase levels increased, intracellular precursor protein accumulation decreased, exhibiting a significant negative correlation with a correlation coefficient of -0.85. This relationship reflects the efficiency of signal peptide-mediated protein transmembrane transport and provides a quantitative basis for evaluating signal peptide function.

[0029] S105. Analyze the relationship between the accumulation level of intracellular folding intermediates and the secretion efficiency of the extracellular mature enzyme. If the accumulation level of the intracellular intermediates is below the threshold, calculate the cleavage site exposure based on the solvent-accessible surface area of the signal peptide, and combine it with the enzyme cleavage time in the in vitro cleavage experiment to generate a cleavage efficiency score.

[0030] Based on the extracellular secretion amount of tannase and the intracellular folding intermediate content, standardized data is obtained through a data normalization method, and a quantitative prediction equation is obtained using a nonlinear regression algorithm; for the intracellular folding intermediate content, a numerical comparator is used to determine whether it is lower than a preset cumulative threshold. If it is lower than the preset cumulative threshold, the corresponding sample sequence number is obtained; based on the sample sequence number, the three-dimensional structure model coordinates are obtained through a molecular modeling method; for the three-dimensional structure model coordinates, the solvent accessibility values of the amino acid residues are obtained through a surface area calculation tool to determine the degree of exposure of the cleavage site; based on the in vitro enzymatic cleavage experimental data, the substrate consumption curve is calculated through a reaction kinetics method, the enzymatic cleavage reaction rate is obtained using a numerical fitting method, and the cleavage efficiency score is calculated based on the degree of exposure of the cleavage site and the enzymatic cleavage reaction rate.

[0031] For example, based on the extracellular secretion amount of tannase and the intracellular folding intermediate content, data normalization methods are used to eliminate dimensional differences, and a nonlinear regression algorithm is used to generate a quantitative prediction equation. A numerical aligner is used to determine whether the intracellular folding intermediate content is below a preset cumulative threshold, and a data screening method is used to obtain the sequence numbers of samples that meet the criteria. Based on the sample sequence numbers, secondary structure data for the corresponding sequences are extracted from a signal peptide database, and a three-dimensional structural model is constructed using molecular modeling methods. For the three-dimensional structural model, atomic coordinates are adjusted using a protein structure optimization tool, and a stable conformation is obtained using an energy minimization method. Based on the coordinates of the stable conformation, the solvent accessibility values of the amino acid residues are obtained using a surface area calculation tool, and the exposure of the cleavage site is calculated using a spatial analysis method. Based on the in vitro enzymatic cleavage experimental data, the substrate consumption curve is calculated using a reaction kinetic method, and the enzymatic cleavage reaction rate is obtained using a numerical fitting method. Based on the exposure of the cleavage site and the enzymatic cleavage reaction rate, a cleavage efficiency score is calculated using a multi-parameter weighted method, and a standardized score is obtained using normalization. During data normalization, the raw data for extracellular secretion of tannase ranged from 0.5 to 4.8 mg / mL, while the intracellular folding intermediate content ranged from 0.2 to 2.5 mg / mL. These data were converted to a standard distribution between 0 and 1 using maximum and minimum normalization. Nonlinear regression analysis revealed a negative exponential correlation between the two, with a correlation coefficient of -0.92 and a coefficient of determination of 0.85 for the prediction equation. A threshold for folding intermediate accumulation was set at 0.8 mg / mL; values above this threshold indicate abnormal intracellular accumulation of the protein. Data screening revealed that 82 of the 128 samples had folding intermediate content below the threshold, and the signal peptide sequences corresponding to these samples were subsequently analyzed. Molecular modeling revealed that the secondary structure of the signal peptide was predicted to contain 65% α-helices, 15% β-sheets, and 20% random coils. The root mean square deviation of the initial structure obtained by homology modeling compared to the known signal peptide crystal structure was 1.2 angstroms, demonstrating high structural accuracy. Structural optimization using molecular force field calculations reduced the interatomic interaction energy from an initial 125 kcal / mol to 42 kcal / mol. The protein backbone dihedral angle distribution fell within the allowed region of the Ramachandran plot, with the root mean square fluctuation of the main chain atoms less than 0.8 angstroms. Solvent accessibility analysis revealed an average exposure of 68% for hydrophilic residues and 22% for hydrophobic residues on the signal peptide surface. Within the five amino acid residues surrounding the cleavage site, the serine protease recognition site was exposed to 85%, significantly higher than the average for other regions. Kinetic experiments revealed that the substrate concentration decreased linearly from an initial 100 μmol / L to 50 μmol / L after 15 minutes and to 25 μmol / L after 30 minutes, demonstrating typical first-order kinetics. The reaction rate constant was 0.062 μmin, and the half-life was 11.2 minutes.Based on the exposure of the cleavage site and the enzyme digestion reaction rate, weighted coefficients of 0.6 and 0.4 were used for the calculation. The cleavage efficiency scores ranged from 0.35 to 0.92, with 35% of the sequences scoring greater than 0.75. These high-efficiency sequences possessed more optimal structural features around the cleavage site, demonstrating faster enzyme digestion rates and higher processing efficiency.

[0032] S106. Optimize the amino acid composition of the hydrophobic core region based on the signal peptide sequence whose cleavage site exposure degree is higher than a threshold value, and generate an optimized signal peptide sequence set.

[0033] The exposure degree of the cleavage site of the signal peptide sequence is numerically compared with the preset threshold to obtain a set of efficient cleavage sequences; sequence structure analysis is performed on the efficient cleavage sequence set to obtain the hydrophobic core region amino acid composition data, hydrophobicity index and side chain volume parameters; the amino acid residues are grouped using a physicochemical feature clustering method, and the spatial arrangement characteristic values are obtained through a residue distribution statistical method; the sequence conformational stability is predicted based on the spatial arrangement characteristic values, and a stable conformation sequence is obtained through a structure evaluation method; based on the stable conformation sequence and tannase secretion data, the hydrophobic core region amino acid composition is optimized through a genetic algorithm to obtain an optimized signal peptide sequence set.

[0034] Exemplarily, based on the cleavage efficiency score, a numerical aligner is used to screen signal peptide sequences with a cleavage site exposure degree higher than a preset threshold to obtain a set of efficient cleavage sequences. For the set of efficient cleavage sequences, the hydrophobic core region amino acid composition data is extracted using a sequence structure analysis tool, and the hydrophobicity index and side chain volume parameters of each site are obtained using an amino acid feature calculation method. Based on the hydrophobicity index and side chain volume parameters, the amino acid residues are grouped using a physicochemical feature clustering method, and the spatial arrangement characteristic values are calculated using a residue distribution statistics method. Based on the spatial arrangement characteristic values, the molecular dynamics method is used to predict the sequence conformational stability, and the stable conformation sequence is screened using a structure evaluation method. Based on the stable conformation sequence and the tannase secretion data, the hydrophobic core region amino acid composition is optimized using a genetic algorithm to generate candidate amino acid combinations. For the candidate amino acid combinations, a sequence splicing tool is used to construct a complete signal peptide sequence, and the secondary structure characteristics of the sequence are verified using a structure prediction method. Based on the secondary structure characteristics, the optimized sequences are scored and ranked using a structure scoring method to obtain an optimized signal peptide sequence set. The exposure threshold of the cleavage site was set at 0.8, and 128 signal peptide sequences were screened. Among them, the exposure of 45 sequences exceeded the threshold, and the cleavage efficiency scores of these sequences were all above 0.75. In the efficient cleavage sequences, the five amino acid residues around the cleavage site are mostly polar or charged residues, which increases the binding affinity of the sequence to the signal peptidase. In the analysis of the amino acid composition of the hydrophobic core region, the hydrophobicity index of each amino acid residue was calculated using the Caltech-Doolittle scale. The hydrophobicity index of isoleucine was the highest at 4.5, and the lowest was aspartic acid at -3.5. The side chain volume parameters showed that tryptophan had the largest volume of 163 cubic angstroms and glycine had the smallest volume of 60 cubic angstroms. These parameters directly affect the spatial arrangement characteristics of the sequence. Amino acid cluster analysis classified the 20 amino acids into four groups: a strongly hydrophobic group consisting of isoleucine, leucine, and valine; a weakly hydrophobic group consisting of alanine, methionine, and phenylalanine; a polar group consisting of serine, threonine, and tyrosine; and a charged group consisting of lysine, arginine, and aspartic acid. Spatial distribution characteristics revealed a periodic distribution of hydrophobic residues in the core region, with a strongly hydrophobic amino acid appearing every three to four residues. Molecular dynamics simulations revealed a root mean square fluctuation of less than 1.2 angstroms within 300 nanoseconds, indicating good structural stability. The backbone hydrogen bond network is intact, and the hydrophobic core region forms a stable α-helical structure, which facilitates conformational maintenance during transmembrane transport. Genetic algorithm optimization was performed using 100 iterations, with the 20% of sequences with the highest fitness selected as parents in each generation. Optimization objectives included a summed hydrophobicity index between 35 and 45, a volume difference of less than 50 cubic angstroms between adjacent residues, and a separation of at least four charged residues. Optimization yielded 28 candidate amino acid combinations. Sequence splicing adopts flexible connection mode, and two glycine residues are added at both ends of the hydrophobic core region as buffer zones.Secondary structure prediction revealed that the optimized sequence had an alpha helix content of 75%, exceeding the 65% of the original sequence. Structural scoring was performed using a multi-parameter weighted approach, including helical integrity (weighted 0.4), hydrophobicity distribution (weighted 0.3), and steric hindrance (weighted 0.3). Twelve high-scoring sequences were identified, all with scores exceeding 0.85.

[0035] S107. Use the optimized signal peptide sequence set to reconstruct the recombinant expression vector and improve the fermentation process formula.

[0036] An optimized signal peptide sequence combination is received, where the signal peptide sequence combination is generated by a sequence optimization module; a vector construction tool is used to perform gene assembly according to the signal peptide sequence combination to obtain a recombinant expression vector library; a parallel fermentation device is used to record temperature, dissolved oxygen, pH value, stirring rate, ventilation rate and nutrient source concentration parameters for the recombinant expression vector library, and fermentation process data is obtained through an online monitoring tool; a principal component analysis method is used to process the fermentation process data, and three main parameters, namely temperature, dissolved oxygen and nutrient source concentration, are extracted; the influence weights of the main parameters on tannase production are calculated through a multivariate statistical method; a response surface experimental plan is constructed based on the influence weights using a central composite experimental design method, and a fermentation process formula is improved through a multi-objective optimization method.

[0037] For example, based on a set of optimized signal peptide sequences, each sequence was assembled using a vector construction tool, and a recombinant expression vector library was constructed using a plasmid transformation method to generate a fermentation strain combination. For each fermentation strain combination, six process parameters, namely temperature, dissolved oxygen, pH, agitation rate, aeration rate, and nutrient source concentration, were recorded using a parallel fermentation apparatus. Fermentation process data was then acquired using an online monitoring tool. Principal component analysis was used to extract the three key parameters, temperature, dissolved oxygen, and nutrient source concentration, from the fermentation process data. Multivariate statistical methods were used to calculate the weights of the parameters' influence on tannase production. Based on the weights of the key parameters, a response surface experiment was constructed using a central composite experimental design, and the effects of parameter combinations were predicted using a fermentation process simulation tool. Based on the predicted parameter combination effects, multiple process recipes were generated using an orthogonal experimental method, and the recipes were verified using a fermentation simulation tool. For the verified process recipes, a cost accounting tool was used to calculate the consumption of consumables during the purification process, and a multi-objective optimization method was used to balance yield and cost. Based on the yield-cost balance results, the process recipes were ranked using a comprehensive scoring method to obtain the optimal fermentation process recipe. During the construction of the recombinant expression vector, the optimized signal peptide sequence was inserted into the vector's multiple cloning site via directional cloning. Each recombinant vector insert ranged in length from 750 to 850 base pairs. A total of 96 strains were generated, establishing a complete fermentation strain library. A parallel fermentation system monitored all 96 fermentation units simultaneously. The temperature was controlled between 25 and 37 degrees Celsius, the dissolved oxygen level was maintained between 20% and 60% saturation, the pH was adjusted between 6.0 and 7.5, the agitation rate varied from 200 to 800 revolutions per minute, the aeration rate was adjusted between 0.5 and 2.0 volumes per minute, and the carbon source concentration varied from 5 to 20 grams per liter. Fermentation data revealed that tannase activity exhibited a typical S-shaped growth curve over time. Principal component analysis revealed that the cumulative contribution of the first three principal components reached 85.6%, with temperature contributing 42.3%, dissolved oxygen 24.8%, and carbon source concentration 18.5%. The influence weights of these parameters on tannase yield were 0.45, 0.32, and 0.23, respectively, indicating that temperature is the most critical process parameter. A central composite experimental design generated 27 parameter combinations, with five temperature levels ranging from 28 to 34 degrees Celsius, three dissolved oxygen levels ranging from 30% to 50%, and three carbon source concentration levels ranging from 8 to 15 grams per liter. Fermentation simulation results showed that the highest tannase yield occurred at 32 degrees Celsius, 40% dissolved oxygen, and 12 grams per liter of carbon source. Orthogonal experiments further optimized the process to produce eight recipes, four of which achieved tannase yields exceeding 3.5 grams per liter. Purification cost calculations revealed that chromatography medium consumption per gram of tannase ranged from 15 to 25 ml, and buffer usage ranged from 200 to 350 ml, for a total consumables cost ranging from 85 to 140 yuan.A multi-objective optimization approach weighed yield against cost, weighting yield at 0.6 and cost at 0.4. Comprehensive scoring revealed that the optimal process recipe was achieved at 32°C, 40% dissolved oxygen, 12 grams per liter of carbon source, a stirring rate of 600 rpm, an aeration rate of 1.2 volumes per minute, and a pH of 6.8. This recipe achieved a tannase yield of 3.8 grams per liter, a purification cost of 95 yuan per gram, and a comprehensive score of 0.92, significantly outperforming other recipes.

[0038] S108. Determine the final tannase strain production plan based on the improved fermentation process formula.

[0039] A bioreactor control device is used to control fermentation parameters, which include temperature, dissolved oxygen, pH value, stirring rate, ventilation rate and nutrient concentration; based on the fermentation parameters, a parameter change curve is recorded using a process monitoring tool to obtain a process stability score; based on the process stability score, an online monitoring tool is used to measure the amount of tannase secretion in the culture medium, and a tannase concentration value is obtained using a fluorescence quantitative method; based on the tannase concentration value, a separation and purification tool is used to perform a purification process, and tannase purification yield data is obtained using a chromatographic separation method; based on the tannase purification yield data, a cost accounting tool is used to calculate the production cost to obtain a final tannase strain production plan.

[0040] For example, based on the optimal fermentation process recipe, process parameters are set using a fermentation process control tool. The bioreactor control device maintains temperature, dissolved oxygen, pH, agitation rate, aeration rate, and nutrient concentration within the set ranges. Based on the bioreactor operating data, a process monitoring tool is used to record the fermentation process parameter change curves. Statistical methods are used to calculate the parameter fluctuation ranges and generate a process stability score. Based on the process stability score, an online monitoring tool is used to measure tannase secretion in real time. Fluorescence quantification is used to calculate tannase concentration and obtain protein yield data. Based on the protein yield data, mass spectrometry is used to identify the target protein components in the culture medium, and protein quantification methods are used to calculate tannase purity. Based on the tannase purity data, a separation and purification tool is used to set purification process parameters, and process validation methods are used to obtain purification yield data. Based on the purification yield data, a cost accounting tool is used to calculate consumables consumption and equipment operating costs to generate tannase scaled production cost data. Based on the scaled production cost data, the production process is optimized using a comprehensive evaluation method to obtain the final tannase strain production plan. Stable control of process parameters is crucial in large-scale fermentation production. The temperature was maintained at 32 ± 0.2°C, the dissolved oxygen level was controlled at 40% ± 2% saturation, the pH was stabilized at 6.8 ± 0.1, the agitation rate was set at 600 ± 20 rpm, the aeration rate was maintained at 1.2 ± 0.05 volume per minute, and the carbon source concentration was controlled at 12 ± 0.5 g / L. Precise control of these parameters directly impacted the stability of fermentation yield. Process monitoring data revealed that the cell density curve exhibited a typical hyperbolic growth pattern during the fermentation cycle, reaching a maximum of 8.5 absorbance units at 48 hours. Parameter fluctuation analysis revealed that the coefficient of variation for temperature was 0.6%, the coefficient of variation for dissolved oxygen was 2.8%, and the coefficient of variation for pH was 1.2%. All parameters were within control ranges, resulting in a process stability score of 0.92. Online monitoring revealed that tannase secretion exhibited a phased change over fermentation time, beginning a rapid increase at 36 hours and reaching a maximum of 4.2 g / L at 72 hours. Real-time tracking of target protein expression using fluorescent labeling revealed a good linear relationship between fluorescence intensity and protein concentration, with a correlation coefficient of 0.98. Mass spectrometry analysis identified the protein components in the culture medium, with tannase accounting for 85% of the total protein, and the main impurity proteins being proteases and lipases secreted by the cells. The purification process uses a two-step chromatography method, with a recovery rate of 92% in the first affinity chromatography step and 88% in the second ion exchange chromatography step, resulting in a final product purity exceeding 98%. Consumables consumed in the purification process mainly include chromatography media, buffers, and filter membranes. The chromatography media has a service life of 50 cycles, 280 ml of buffer is consumed per gram of product, and the deep filter membrane uses an area of 0.5 square meters per kilogram of product. Equipment operating costs include electricity consumption, cooling water consumption, and compressed air consumption, with a total cost of 75 yuan per gram of product.The comprehensive evaluation was based on multiple dimensions, including yield stability (weighted 0.3), product purity (weighted 0.3), process reliability (weighted 0.2), and production cost (weighted 0.2). The finalized production plan achieved reasonable control of production costs while ensuring product quality, with an overall score of 0.89, meeting the requirements for large-scale production.

[0041] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for screening high-yield tannase strains based on genome engineering, characterized in that: The method comprises: A signal peptide conformational feature dataset was constructed based on a collection of tannase signal peptide sequences. The effects of changes in the amino acid composition of the hydrophobic core region on the conformational stability of the signal peptide and changes in the binding affinity of the signal peptide to the transport channel on the endoplasmic reticulum membrane were analyzed. The transport channel recognition efficiency score was determined based on the analysis results of conformational stability and binding affinity changes. If the transport channel recognition efficiency score was higher than a preset threshold, the corresponding hydrophobic core region sequence was fused with the tannase gene to generate a recombinant expression vector. Fermentation experiments were carried out based on the recombinant expression vector. The cleavage site exposure was calculated according to the solvent-accessible surface area of the signal peptide to optimize the amino acid composition of the hydrophobic core region and determine the final tannase strain production plan.

2. The method according to claim 1, characterized in that The signal peptide conformation feature data set is constructed based on the tannase signal peptide sequence set, including: Obtaining a tannase signal peptide sequence set according to a preset genome database, and performing amino acid residue coding standardization on the signal peptide sequence set to obtain a standardized signal peptide sequence; For the standardized signal peptide sequence, the hydrophobicity of the amino acid residues is scored, the hydrophobicity threshold is set, the boundary of the hydrophobic core region is located, the types and occurrence frequencies of the amino acid residues in the hydrophobic core region are counted, and a hydrophobic core region amino acid composition feature vector is generated; For the standardized signal peptide sequence, the probability values of three conformations of amino acid residues, namely, α-helix, β-sheet and random coil, are calculated to obtain a secondary structure prediction probability matrix; Each conformational segment was scored for stability based on the secondary structure prediction probability matrix; The signal peptide sequence feature dataset was constructed by integrating the hydrophobic core region amino acid composition feature vector, secondary structure prediction probability matrix and conformational segment stability score.

3. The method according to claim 2, characterized in that The analysis of the effect of changes in the amino acid composition of the hydrophobic core region on the conformational stability of the signal peptide and the changes in the binding affinity of the signal peptide to the transport channel on the endoplasmic reticulum membrane; and the determination of the transport channel recognition efficiency score based on the analysis results of the conformational stability and the changes in the binding affinity include: Obtain the hydrophobic core region amino acid sequence based on the signal peptide conformational feature dataset, identify the variation sites in the hydrophobic core region amino acid sequence, and obtain the difference values of the amino acid physicochemical properties; Calculating protein conformational parameters based on the difference values of the physicochemical properties of the amino acids to generate a conformational stability characteristic matrix; Based on the conformational stability characteristic matrix, the spatial matching degree between the hydrophobic core region of the signal peptide and the binding interface of the endoplasmic reticulum membrane transport channel is calculated to obtain a set of binding site coordinates; Calculating the interaction energy value at the binding site according to the binding site coordinate set to generate a binding affinity feature vector; A transport channel recognition efficiency score is calculated based on the conformational stability feature matrix and the binding affinity feature vector.

4. The method according to claim 1, wherein If the transport channel recognition efficiency score is higher than a preset threshold, the corresponding hydrophobic core region sequence is fused with the tannase gene to obtain a recombinant expression vector, comprising: Obtaining a transport channel recognition efficiency score, and determining whether the score exceeds a preset threshold based on the score; if the score exceeds the preset threshold, deriving a hydrophobic core region sequence from the signal peptide conformational feature dataset and generating a target DNA sequence; Determine the tannase gene sequence splicing site according to the target DNA sequence, and perform a directionally connected connection between the target DNA sequence and the tannase gene sequence using a gene splicing tool to obtain a fusion gene sequence; According to the fusion gene sequence, the codon composition is sequence optimized, matching expression regulatory elements are selected for the optimized fusion gene sequence, and a recombinant expression vector is constructed using a vector design tool.

5. The method according to claim 1 or 4, characterized in that The fermentation experiment based on the recombinant expression vector is carried out, and the cleavage site exposure is calculated according to the solvent accessible surface area of the signal peptide to optimize the amino acid composition of the hydrophobic core region, and the final tannase strain production plan is determined, including: Fermentation experiments were conducted based on a high-throughput signal peptide library containing recombinant expression vectors to detect the amount of tannase secreted outside the cells, and the accumulation of folding intermediates was assessed by analyzing the intracellular content of tannase precursor protein. The relationship between the accumulation level of intracellular folding intermediates and the secretion efficiency of the extracellular mature enzyme was analyzed. If the accumulation level of intracellular intermediates was lower than the threshold, the cleavage site exposure was calculated based on the solvent-accessible surface area of the signal peptide, and the cleavage efficiency score was generated based on the enzyme cleavage time in the in vitro cleavage experiment. Optimize the amino acid composition of the hydrophobic core region based on the signal peptide sequence whose cleavage site exposure is higher than the threshold, and generate an optimized signal peptide sequence set; Using an optimized signal peptide sequence set, reconstructing the recombinant expression vector, and improving the fermentation process formula; The final tannase strain production plan was determined based on the improved fermentation process formula.

6. The method according to claim 5, characterized in that The fermentation experiment is conducted based on a high-throughput signal peptide library containing a recombinant expression vector to detect the amount of tannase secreted outside the cell, and the accumulation level of folding intermediates is evaluated by analyzing the content of intracellular tannase precursor protein, including: Obtaining amino acid substitution data for conserved sites in the signal peptide sequence, and generating a mutation site combination matrix based on the amino acid substitution data; According to the mutation site combination matrix, a recombinant expression vector is constructed using a gene synthesis tool. The recombinant expression vector is transformed to obtain an expression strain library. A fermentation monitoring device is used to record dissolved oxygen concentration and cell density data. The fermentation culture supernatant is obtained using a bioreactor control device, and the secretory tannase content data is measured to obtain protein secretion efficiency. Proteins suitable for predetermined requirements are screened based on protein secretion efficiency, and the cells are broken to obtain intracellular protein extracts to evaluate the accumulation level of folding intermediates.

7. The method according to claim 5, characterized in that The analysis is performed to determine the relationship between the accumulation level of intracellular folding intermediates and the secretion efficiency of the extracellular mature enzyme. If the accumulation level of the intracellular intermediates is lower than the threshold, the cleavage site exposure is calculated based on the solvent-accessible surface area of the signal peptide, and the cleavage efficiency score is generated in combination with the enzyme cleavage time in the in vitro cleavage experiment, including: According to the extracellular secretion amount of tannase and the intracellular folding intermediate content, the standardized data were obtained by data normalization method, and the quantitative prediction equation was obtained by nonlinear regression algorithm. Determine whether the intracellular folding intermediate content value is lower than a preset cumulative threshold. If it is lower than the preset cumulative threshold, obtain the corresponding sample sequence number and obtain the three-dimensional structure model coordinates through molecular modeling methods; Based on the coordinates of the three-dimensional structure model, the solvent accessibility values of the amino acid residues are obtained using a surface area calculation tool to determine the exposure degree of the cleavage site; Based on the in vitro enzymatic cleavage experimental data, the substrate consumption curve was calculated by the reaction kinetics method, the enzymatic cleavage reaction rate was obtained by the numerical fitting method, and the cleavage efficiency score was calculated based on the degree of exposure of the cleavage site and the enzymatic cleavage reaction rate.

8. The method according to claim 5, characterized in that The method optimizes the amino acid composition of the hydrophobic core region according to the signal peptide sequence whose cleavage site exposure degree is higher than the threshold value to generate an optimized signal peptide sequence set, including: The exposure degree of the cleavage site of the signal peptide sequence is numerically compared with the preset threshold to obtain a set of efficient cleavage sequences; Performing sequence structure analysis on the efficient cleavage sequence set to obtain hydrophobic core region amino acid composition data and hydrophobicity index and side chain volume parameters; The physicochemical feature clustering method is used to group amino acid residues, obtain spatial arrangement feature values, predict sequence conformation stability, and obtain stable conformation sequences; According to the stable conformation sequence and tannase secretion data, the amino acid composition of the hydrophobic core region was optimized by genetic algorithm to obtain the optimized signal peptide sequence set.

9. The method according to claim 5, characterized in that The method of using the optimized signal peptide sequence set to reconstruct the recombinant expression vector and improve the fermentation process formula includes: receiving an optimized signal peptide sequence combination generated by a sequence optimization module; Gene assembly is performed using a vector construction tool according to the signal peptide sequence combination to obtain a recombinant expression vector library, and temperature, dissolved oxygen, pH value, stirring rate, aeration rate and nutrient source concentration parameters are recorded to obtain fermentation process data; The principal component analysis method was used to process the fermentation process data, and the three main parameters of temperature, dissolved oxygen and nutrient source concentration were extracted. The influence weights of the main parameters on tannase production were calculated using multivariate statistical methods. A response surface experimental plan was constructed, and the fermentation process formula was improved through a multi-objective optimization method.

10. The method according to claim 5, characterized in that The final tannase strain production plan is determined according to the improved fermentation process formula, including: The fermentation parameters are controlled by a bioreactor control device, including temperature, dissolved oxygen, pH value, stirring rate, aeration rate and nutrient concentration; According to the fermentation parameters, the parameter change curve is recorded by a process monitoring tool to obtain a process stability score value, the tannase secretion amount in the culture medium is measured by an online monitoring tool, and the tannase concentration value is obtained by a fluorescence quantitative method; According to the tannase concentration value, a separation and purification tool is used to perform a purification process, the tannase purification yield data is obtained through a chromatographic separation method, the production cost is calculated using a cost accounting tool, and the final tannase strain production plan is obtained.

Images