Improving the Thermostability and Activity of a Thermophilic Subtilase by Incorporating Structural Elements of Its Psychrophilic Counterpart

Enhancing Machine-Learning Prediction of Enzyme Catalytic Temperature Optima through Amino Acid Conservation Analysis

Enzymes play a crucial role in various industrial production and pharmaceutical developments, serving as catalysts for numerous biochemical reactions. Determining the optimal catalytic temperature (Topt) of enzymes is crucial for optimizing reaction conditions, enhancing catalytic efficiency, and accelerating the industrial processes. However, due to the limited availability of experimentally determined Topt data and the insufficient accuracy of existing computational methods in predicting Topt, there is an urgent need for a computational approach to predict the Topt values of enzymes accurately. In this study, using phosphatase (EC 3.1.3.X) as an example, we constructed a machine learning model utilizing amino acid frequency and protein molecular weight information as features and employing the K-nearest neighbors regression algorithm to predict the Topt of enzymes. Usually, when conducting engineering for enzyme thermostability, researchers tend not to modify conserved amino acids. Therefore, we utilized this machine learning model to predict the Topt of phosphatase sequences after removing conserved amino acids. We found that the predictive model's mean coefficient of determination (R2) value increased from 0.599 to 0.755 compared to the model based on the complete sequences. Subsequently, experimental validation on 10 phosphatase enzymes with undetermined optimal catalytic temperatures shows that the predicted values of most phosphatase enzymes based on the sequence without conservative amino acids are closer to the experimental optimal catalytic temperature values. This study lays the foundation for the rapid selection of enzymes suitable for industrial conditions.

Genomic analyses reveal a low-temperature adapted clade in Halorubrum, a widespread haloarchaeon across global hypersaline environments

BACKGROUND

Cold-adapted archaea have diverse ecological roles in a wide range of low-temperature environments. Improving our knowledge of the genomic features that enable psychrophiles to grow in cold environments helps us to understand their adaptive responses. However, samples from typical cold regions such as the remote Arctic and Antarctic are rare, and the limited number of high-quality genomes available leaves us with little data on genomic traits that are statistically associated with cold environmental conditions.

RESULTS

In this study, we examined the haloarchaeal genus Halorubrum and defined a new clade that represents six isolates from polar and deep earth environments ('PD group' hereafter). The genomic G + C content and amino acid composition of this group distinguishes it from other Halorubrum and the trends are consistent with the established genomic optimization of psychrophiles. The cold adaptation of the PD group was further supported by observations of increased flexibility of proteins encoded across the genome and the findings of a growth test.

CONCLUSIONS

The PD group Halorubrum exhibited denser genome packing, which confers higher metabolic potential with constant genome size, relative to the reference group, resulting in significant differences in carbon, nitrogen and sulfur metabolic patterns. The most marked feature was the enrichment of genes involved in sulfur cycling, especially the production of sulfite from organic sulfur-containing compounds. Our study provides an updated view of the genomic traits and metabolic potential of Halorubrum and expands the range of sources of cold-adapted haloarchaea.

Improving the catalytic activity of a detergent-compatible serine protease by rational design

Serine proteases are among the most important biological additives in various industries such as detergents, leather, animal feed and food. A serine protease gene, Fgapt4, from Fusarium graminearum 2697 was identified, cloned and expressed in Pichia pastoris. The optimal pH and temperature of FgAPT4 were 8.5 and 40°C, respectively. The relative activity was >30% even at 10°C. It had a wide range of pH stability (4.0-12.0) and detergent compatibility. To improve the catalytic activity, a strategy combining molecular docking and evolutionary analysis was adopted. Twelve amino acid residue sites and three loops (A, B and C) were selected as potential hot spots that might play critical roles in the enzyme's functional properties. Twenty-eight mutants targeting changes in individual sites or loops were designed, and mutations with good performance were combined. The best mutant was FgAPT4-M3 (Q70N/D142S/A143S/loop C). The specific activity and catalytic efficiency of FgAPT4-M3 increased by 1.6 (1008.5 vs. 385.9 U/mg) and 2.2-fold (3565.1 vs. 1106.3/s/mM), respectively. Computational analyses showed that the greater flexibility of the substrate pocket may be responsible for the increased catalytic activity. In addition, its application in detergents indicated that FgAPT4-M3 has great potential in washing.

Auto- and Hetero-Catalytic Processing of the N-Terminal Propeptide Promotes the C-Terminal Fibronectin Type III Domain-Mediated Dimerization of a Thermostable Vpr-like Protease

Microbial Vpr-like proteases are extracellular multidomain subtilases with diverse functions and can form oligomers, but their maturation and oligomerization mechanisms remain to be elucidated. Here, we report a novel Vpr-like protease (BTV) from thermophilic bacterium Brevibacillus sp. WF146. The BTV precursor comprises a signal peptide, an N-terminal propeptide, a subtilisin-like catalytic domain with an inserted protease-associated (PA) domain, two tandem fibronectin type III domains (Fn1 and Fn2), and a C-terminal propeptide. The BTV proform (pro-BTV) could be autoprocessed into the mature form (mBTV) via two intermediates lacking the N- or C-terminal propeptide, respectively, and the C-terminal propeptide delays the autocatalytic maturation of the enzyme. By comparison, pro-BTV is more efficiently processed into mBTV by protease TSS from strain WF146. Purified mBTV is a Ca²⁺-dependent thermostable protease, showing optimal activity at 60°C and retaining more than 60% of activity after incubation at 60°C for 8 h. The PA domain is important for enzyme stability and contributes to the substrate specificity of BTV by restricting the access of protein substrates to the active site. The proform and mature form of BTV exist as a monomer and a homodimer, respectively, and the dimerization is mediated by the Fn1 and Fn2 domains. The N-terminal propeptide of BTV not only acts as intramolecular chaperone and enzymatic inhibitor but also inhibits the homodimerization of the enzyme. The removal of the N-terminal propeptide leads to a structural adjustment of the enzyme and thus promotes enzyme dimerization. IMPORTANCE Vpr-like proteases are widely distributed in bacteria and fungi and are involved in processing lantibiotics, degrading collagen, keratin, and fibrin, and pathogenesis of microbes. The dissection of the roles of individual domains in enzyme maturation and oligomerization is crucial for understanding the action mechanisms of these multidomain proteases. Our results demonstrate that hetero-catalytic maturation of the extracellular Vpr-like protease BTV of Brevibacillus sp. WF146 is more efficient than autocatalytic maturation of the enzyme. Moreover, we found that the C-terminal tandem fibronectin type III domains rather than the PA domain mediate the dimerization of mature BTV, while the N-terminal propeptide inhibits the dimerization of the BTV proform. This study provides new insight into the activation and oligomerization mechanisms of Vpr-like proteases.

Maturation process and characterization of a novel thermostable and halotolerant subtilisin-like protease with high collagenolytic but low gelatinolytic activity

Enzymatic degradation of collagen is of great industrial and environmental significance; however, little is known about thermophile-derived collagenolytic proteases. Here, we report a novel collagenolytic protease (TSS) from thermophilic Brevibacillus sp. WF146. The TSS precursor comprises a signal peptide, an N-terminal propeptide, a subtilisin-like catalytic domain, a β-jelly roll (βJR) domain, and a prepeptidase C-terminal (PPC) domain. The maturation of TSS involves a stepwise autoprocessing of the N-terminal propeptide and the PPC domain, and the βJR rather than the PPC domain is necessary for correct folding of the enzyme. Purified mature TSS displayed optimal activity at 70°C and pH 9.0, a half-life of 1.5 h at 75°C, and an increased thermostability with rising salinity up to 4 M. TSS possesses an increased number of surface acidic residues and ion pairs, as well as four Ca 2+ -binding sites, which contribute to its high thermostability and halotolerance. At high temperatures, TSS exhibited high activity toward insoluble type I collagen and azocoll, but showed a low gelatinolytic activity, with a strong preference for Arg and Gly at the P1 and P1’ positions, respectively. Both the βJR and PPC domains could bind but not swell collagen, and thus facilitate TSS-mediated collagenolysis via improving the accessibility of the enzyme to the substrate. Additionally, TSS has the ability to efficiently degrade fish scale collagen at high temperatures.

IMPORTANCE

Proteolytic degradation of collagen at high temperatures has the advantages of increasing degradation efficiency and minimizing the risk of microbial contamination. Reports on thermostable collagenolytic proteases are limited, and their maturation and catalytic mechanisms remain to be elucidated. Our results demonstrate that the thermophile-derived TSS matures in an autocatalytic manner, and represents one of the most thermostable collagenolytic proteases reported so far. At elevated temperatures, TSS prefers hydrolyzing insoluble heat-denatured collagen rather than gelatin, providing new insight into the mechanism of collagen degradation by thermostable collagenolytic proteases. Moreover, TSS has the potential to be used in recycling collagen-rich wastes such as fish scales.

Enhancing extracellular production of recombinant proteins in Escherichia coli by co-expressing with a haloarchaeal protein containing a putative LolA-like domain

Escherichia coli represents one of the most widely used hosts for recombinant protein production, but its limited capacity for producing extracellular proteins is often cited as a drawback. NJ7G_0991 is an extracellular protein of the haloarchaeon Natrinema sp. J7-2 and comprises a signal peptide, a putative LolA-like domain, and a C-terminal domain of unknown function. Here, we found that the full-length (0991) and the C-terminal domain-deletion variant (0991ΔC) of NJ7G_0991, but not its signal peptide-deletion variant (0991ΔS), were efficiently released into the culture supernatant of E. coli without extensive cell lysis as determined by β-galactosidase activity assay. After lysozyme treatment, E. coli cells producing 0991 or 0991ΔC, but not 0991ΔS, were converted from rod-shaped forms to spheres, suggesting that the secretion of 0991 or 0991ΔC into the periplasm leads to an increase of outer membrane permeability of E. coli. A pelB signal peptide was fused to the N-terminus of the LolA-like domain, and the resulting variant PelB-0991ΔC could be released into the culture supernatant of E. coli more efficiently than 0991ΔC. By using PelB-0991ΔC as a co-expression partner, the extracellular production level of a recombinant thermostable subtilase WF146 could be enhanced by up to 14-fold, and the extracellular concentration of an active site variant of WF146 (WF146-SA) reached up to 129 mg/l. To the best of our knowledge, this is the first report on archaeal protein-based co-expression system for extracellular production of recombinant proteins in E. coli. KEY POINTS: • The haloarchaeal protein NJ7G_0991 can be efficiently released into the culture supernatant of E. coli. • The recombinant NJ7G_0991 increases the outer membrane permeability of E. coli. • The LolA-like domain of NJ7G_0991 can be used as a co-expression partner to improve extracellular production of recombinant proteins in E. coli.

Sortase mutants with improved protein thermostability and enzymatic activity obtained by consensus design

Staphylococcus aureus sortase A (SaSrtA) is an enzyme that anchors proteins to the cell surface of Gram-positive bacteria. During the transpeptidation reaction performed by SaSrtA, proteins containing an N-terminal glycine can be covalently linked to another protein with a C-terminal LPXTG motif (X being any amino acid). Since the sortase reaction can be performed in vitro as well, it has found many applications in biotechnology. Although sortase-mediated ligation has many advantages, SaSrtA is limited by its low enzymatic activity and dependence on Ca2+. In our study, we evaluated the thermodynamic stability of the SaSrtA wild type and found the enzyme to be stable. We applied consensus analysis to further improve the enzyme's stability while at the same time enhancing the enzyme's activity. As a result, we found thermodynamically improved, more active and Ca2+-independent mutants. We envision that these new variants can be applied in conjugation reactions in low Ca2+ environments.

Droplet-based optofluidic systems for measuring enzyme kinetics

The study of enzyme kinetics is of high significance in understanding metabolic networks in living cells and using enzymes in industrial applications. To gain insight into the catalytic mechanisms of enzymes, it is necessary to screen an enormous number of reaction conditions, a process that is typically laborious, time-consuming, and costly when using conventional measurement techniques. In recent times, droplet-based microfluidic systems have proved themselves to be of great utility in large-scale biological experimentation, since they consume a minimal sample, operate at high analytical throughput, are characterized by efficient mass and heat transfer, and offer high levels of integration and automation. The primary goal of this review is the introduction of novel microfluidic tools and detection methods for use in high-throughput and sensitive analysis of enzyme kinetics. The first part of this review focuses on introducing basic concepts of enzyme kinetics and describing most common microfluidic approaches, with a particular focus on segmented flow. Herein, the key advantages include accurate control over the flow behavior, efficient mass and heat transfer, multiplexing, and high-level integration with detection modalities. The second part describes the current state-of-the-art platforms for high-throughput and sensitive analysis of enzyme kinetics. In addition to our categorization of recent advances in measuring enzyme kinetics, we have endeavored to critically assess the limitations of each of these detection approaches and propose strategies to improve measurements in droplet-based microfluidics. Graphical abstract.

Engineering enzyme catalysis: an inverse approach

Enzymes' inherent chirality confers their exquisite enantiomeric specificity and makes their use as green alternatives to chiral metal complexes or chiral organocatalysts invaluable to the fine chemical industry. The most prevalent way to alter enzyme activity in terms of regioselectivity and stereoselectivity for both industry and fundamental research is to engineer the enzyme. In a recent article by Keinänen et al., published in Bioscience Reports 2018, 'Controlling the regioselectivity and stereoselectivity of FAD-dependent polyamine oxidases with the use of amine-attached guide molecules as conformational modulators', an inverse approach was presented that focuses on the manipulation of the enzyme substrate rather than the enzyme. This approach not only uncovered dormant enantioselectivity in related enzymes but allowed for its control by the use of guide molecules simply added to the reaction solution or covalently linked to an achiral scaffold molecule.

Computation-Aided Rational Deletion of C-Terminal Region Improved the Stability, Activity, and Expression Level of GH2 β-Glucuronidase

In this study, computation-aided design on the basis of structural analysis was employed to rationally identify a highly dynamic C-terminal region that regulates the stability, expression level, and activity of a GH2 fungal glucuronidase from Aspergillus oryzae Li-3 (PGUS). Then, four mutants with a precisely truncated C-terminal region in different lengths were constructed; among them, mutant D591-604 with a 3.8-fold increase in half-life at 65 °C and a 6.8 kJ/mol increase in Gibbs free energy showed obviously improved kinetic and thermodynamic stability in comparison to PGUS. Mutants D590-604 and D591-604 both showed approximately 2.4-fold increases in the catalytic efficiency k_cat/ K_m and 1.8-fold increases in the expression level. Additionally, the expression level of PGUS was doubled through a C-terminal region swap with bacterial GUS from E. coli (EGUS). Finally, the robust PGUS mutants D590-604 and D591-604 were applied in the preparation of glycyrrhetinic acid with 4.0- and 4.4-fold increases in concentration through glycyrrhizin hydrolysis by a fed-batch process.

Trade-offs with stability modulate innate and mutationally acquired drug resistance in bacterial dihydrofolate reductase enzymes

Structural stability is a major constraint on the evolution of protein sequences. However, under strong directional selection, mutations that confer novel phenotypes but compromise structural stability of proteins may be permissible. During the evolution of antibiotic resistance, mutations that confer drug resistance often have pleiotropic effects on the structure and function of antibiotic-target proteins, usually essential metabolic enzymes. In the present study, we show that trimethoprim (TMP)-resistant alleles of dihydrofolate reductase from Escherichia coli (EcDHFR) harboring the Trp30Gly, Trp30Arg or Trp30Cys mutations are significantly less stable than the wild-type, making them prone to aggregation and proteolysis. This destabilization is associated with a lower expression level, resulting in a fitness cost and negative epistasis with other TMP-resistant mutations in EcDHFR. Using structure-based mutational analysis, we show that perturbation of critical stabilizing hydrophobic interactions in wild-type EcDHFR enzyme explains the phenotypes of Trp30 mutants. Surprisingly, though crucial for the stability of EcDHFR, significant sequence variation is found at this site among bacterial dihydrofolate reductases (DHFRs). Mutational and computational analyses in EcDHFR and in DHFR enzymes from Staphylococcus aureus and Mycobacterium tuberculosis demonstrate that natural variation at this site and its interacting hydrophobic residues modulates TMP resistance in other bacterial DHFRs as well, and may explain the different susceptibilities of bacterial pathogens to TMP. Our study demonstrates that trade-offs between structural stability and function can influence innate drug resistance as well as the potential for mutationally acquired drug resistance of an enzyme.

Release of an HtrA-Like Protease from the Cell Surface of Thermophilic Brevibacillus sp. WF146 via Substrate-Induced Autoprocessing of the N-terminal Membrane Anchor

High-temperature requirement A (HtrA)-like proteases participate in protein quality control in prokaryotes and eukaryotes by degrading damaged proteins; however, little is known about HtrAs produced by thermophiles. HtrAw is an HtrA-like protease of thermophilic Brevibacillus sp. WF146. The intact form of HtrAw (iHtrAw) consisting of a transmembrane segment-containing N-terminal domain, a trypsin-like protease domain, and a C-terminal PDZ domain was produced in Escherichia coli. Purified iHtrAw itself is unable to cleave the N-terminal domain, but requires protein substrates to autoprocess the N-terminal domain intermolecularly, yielding a short form (sHtrAw). Mutation at the substrate-binding site in the PDZ domain affects the conversion of iHtrAw to sHtrAw. Deletion analysis revealed that the N-terminal domain is not necessary for enzyme folding, activity, and thermostability. Compared with other known HtrAs, HtrAw contains an additional Ca²⁺-binding Dx[DN]xDG motif important for enzyme stability and/or activity. When produced in an htrA/htrB double deletion mutant of Bacillus subtilis, iHtrAw localized predominantly to the cell pellet, and the amount of sHtrAw in the culture supernatant increased at elevated temperatures. Moreover, HtrAw increased the heat resistance of the B. subtilis mutant. In strain WF146, HtrAw exists in both a cell-associated intact form and a cell-free short form; an increase in growth temperature enhanced HtrAw production and the amount of cell-free short form. Release of the short form of HtrAw from the membrane may have the advantage of allowing the enzyme to freely access and degrade damaged proteins surrounding the bacterium living at high temperatures.

Four Inserts within the Catalytic Domain Confer Extra Stability and Activity to Hyperthermostable Pyrolysin from Pyrococcus furiosus

Pyrolysin from the hyperthermophilic archaeon Pyrococcus furiosus is the prototype of the pyrolysin family of the subtilisin-like serine protease superfamily (subtilases). It contains four inserts (IS147, IS29, IS27, and IS8) of unknown function in the catalytic domain. We performed domain deletions and showed that three inserts are either essential (IS147 and IS27) or important (IS8) for efficient maturation of pyrolysin at high temperatures, whereas IS29 is dispensable. The large insert IS147 contains Ca3 and Ca4, two calcium-binding Dx[DN]xDG motifs that are conserved in many pyrolysin-like proteases. Mutagenesis revealed that the Ca3 site contributes to enzyme thermostability and the Ca4 site is necessary for pyrolysin to fold into a maturation-competent conformation. Mature insert-deletion variants were characterized and showed that IS29 and IS8 contribute to enzyme activity and stability, respectively. In the presence of NaCl, pyrolysin undergoes autocleavage at two sites: one within IS29 and the other in IS27 Disrupting the ion pairs in IS27 and IS8 induces autocleavage in the absence of salts. Interestingly, autocleavage products combine noncovalently to form an active, nicked enzyme that is resistant to SDS and urea denaturation. Additionally, a single mutation in IS29 increases resistance to salt-induced autocleavage and further increases enzyme thermostability. Our results suggest that these extra structural elements play a crucial role in adapting pyrolysin to hyperthermal environments.IMPORTANCE Pyrolysin-like proteases belong to the subtilase superfamily and are characterized by large inserts and long C-terminal extensions; however, the role of the inserts in enzyme function is unclear. Our results demonstrate that four inserts in the catalytic domain of hyperthermostable pyrolysin contribute to the folding, maturation, stability, and activity of the enzyme at high temperatures. The modification of extra structural elements in pyrolysin-like proteases is a promising strategy for modulating global structure stability and enzymatic activity of this class of protease.