Improving the thermostability of methyl parathion hydrolase from Ochrobactrum sp. M231 using a computationally aided method

A Dual-Focus Workflow for Simultaneously Engineering High Activity and Thermal Stability in Methyl Parathion Hydrolase

Industrial fermentation applications typically require enzymes that exhibit high stability and activity at high temperatures. However, efforts to simultaneously improve these properties are usually limited by a trade-off between stability and activity. This report describes a computational strategy to enhance both activity and thermal stability of the mesophilic organophosphate-degrading enzyme, methyl parathion hydrolase (MPH). To predict hotspot mutation sites, we assembled a library of features associated with the target properties for each residue and then prioritized candidate sites by hierarchical clustering. Subsequent in silico screening with multiple algorithms to simulate selective pressures yielded a subset of 23 candidate mutations. Iterative parallel screening of mutations that improved thermal stability and activity yielded, MPHase-m5b, which exhibited 13.3 °C higher Tm and 4.2 times higher catalytic activity than wild-type (WT) MPH over a wide temperature range. Systematic analysis of crystal structures, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations revealed a wider entrance to the active site that increased substrate access with an extensive network of interactions outside the active site that reinforced αβ/βα sandwich architecture to improve thermal stability. This study thus provides an advanced, rational design framework to improve efficiency in engineering highly active, thermostable biocatalysts for industrial applications.

Effective remediation and decontamination of organophosphorus compounds using enzymes: From rational design to potential applications

Organophosphorus compounds (OPs) have been widely used in agriculture for decades because of their high insecticidal efficiency, which maintains and increases crop yields worldwide. More importantly, OPs, as typical chemical warfare agents, are a serious concern and significant danger for military and civilian personnel. The widespread use of OPs, superfluous and unreasonable use, has caused great harm to the environment and food chain. Developing efficient and environmentally friendly solutions for the decontamination of OPs is a long-term challenge. Microbial enzymes show potential application as natural and green biocatalysts. Thus, utilizing OP-degrading enzymes for environmental decontamination presents significant advantages, as these enzymes can rapidly hydrolyze OPs; are environmentally friendly, nonflammable, and noncorrosive; and can be discarded safely and easily. Here, the properties, structure and catalytic mechanism of various typical OP-degrading enzymes are reviewed. The methods and effects utilized to improve the expression level, catalytic performance and stability of OP-degrading enzymes were systematically summarized. In addition, the immobilization of OP-degrading enzymes was explicated emphatically, and the latest progress of cascade reactions based on immobilized enzymes was discussed. Finally, the latest applications of OP-degrading enzymes were summarized, including biosensors, nanozyme mimics and medical detoxification. This review provides guidance for the future development of OP-degrading enzymes and promotes their application in the field of environmental bioremediation and medicine.

Homologous Expression and Characterization of α-L-rhamnosidase from Aspergillus niger for the Transformation of Flavonoids

Aspergillus niger has been used for homologous and heterologous expressions of many protein products. In this study, the α-L-rhamnosidase from A. niger (Rha-N1, GenBank XP_001389086.1) was homologously expressed in A. niger 3.350 by Agrobacterium tumefaciens-mediated transformation. The enzyme activity of Rha-N1 was 0.658 U/mL, which was obtained by cultivation of engineered A. niger in a 5-L bioreactor. Rha-N1 was purified by affinity chromatography and characterized. The optimum temperature and optimum pH for Rha-N1 were 60 °C and 4.5, respectively. Enzyme activity was promoted by Al³⁺, Li⁺, Mg²⁺, and Ba²⁺ and was inhibited by Mn²⁺, Fe³⁺, Ca²⁺, Cu²⁺, and organic solvents. The result indicated that rutin was the most suitable substrate for Rha-N1 by comparison with the other two flavonoid substrates hesperidin and naringin. The transformed products of isoquercitrin, hesperetin-7-O-glucoside, and prunin were identified by LC-MS and ¹H-NMR.

Overview of a bioremediation tool: organophosphorus hydrolase and its significant application in the food, environmental, and therapy fields

In the past decades, the organophosphorus compounds had been widely used in the environment and food industries as pesticides. Owing to the life-threatening and long-lasting problems of organophosphorus insecticide (OPs), an effective detection and removal of OPs have garnered growing attention both in the scientific and practical fields in recent years. Bacterial organophosphorus hydrolases (OPHs) have been extensively studied due to their high specific activity against OPs. OPH could efficiently hydrolyze a broad range of substrates both including the OP pesticides and some nerve agents, suggesting a great potential for the remediation of OPs. In this review, the microbial identification, molecular modification, and practical application of OPHs were comprehensively discussed.Key points• Microbial OPH is a significant bioremediation tool against OPs.• Identification and molecular modification of OPH was discussed in detail.• The applications of OPH in food, environmental, and therapy fields are presented.

Molecular cloning and characterization of a methyl parathion hydrolase from an organophosphorus-degrading bacterium, Serratia marcescens MEW06

An organophosphorus-degrading bacterium MEW06, which exhibited excellent biodegradation capabilities towards 50 mg/L of methyl parathion (MP), paraoxon and dimethoate, was isolated from Sand Lake (Wuhan, China) and identified as Serratia marcescens subsp. marcescens based on physiological-biochemical characteristics and a 16S rDNA sequence-based phylogenetic tree. MEW06 genome contains a 31.09-kDa methyl parathion hydrolase (MPH) (MPHGM004539) that was 54.9% similar to Pseudomonas sp. WBC-3's MPH. RT-qPCR revealed that mphGM004539 gene expression was significant up-regulated when co-cultured with MP. mphGM004539 without signal peptide (mphGM004539Δsp) was successful cloned and expressed in Escherichia coli BL21 (DE3). Optimized specific enzyme activity of MPHGM004539ΔSP was 5.26 U/mg under 35°C and pH 11.0 conditions when MP as the substrate. Additionally, Co2+, Cd2+and Fe2+ increased the enzyme activity level. MP could be degraded by MPHGM004539ΔSP into p-nitrophenol probably by hydrolyzing the P-O ester bond. Virulence of MP towards Drosophila melanogaster W1118 was reduced by MEW06 or MPHGM004539ΔSP biodegradation. This is the first cloning and characterization of MPH from the organophosphorus-degrading bacterium S. marcescens. MEW06 and its MPH have potential roles in the bioremediation of organophosphorus pesticide-contaminated eco-systems.

Consensus design for improved thermostability of lipoxygenase from Anabaena sp. PCC 7120

Background

Lipoxygenase (LOX) from Anabaena sp. PCC 7120 (Ana-rLOX) offers important applications in the food industry, especially for improving aroma and dough rheological properties. However, industrial applications of LOXs have been limited by their poor thermostability. Herein, we report a bioinformatics-based consensus concept approach for the engineering of thermostable Ana-rLOX.

Results

A series of mutations (N130D, G260A, S437T, N130D/G260Q, N130D/S437Y) showed higher thermostability and activity than the wild-type enzyme. Thus, N130D/G260Q exhibited a 6.6-fold increase in half-life and 2.45 °C increase in unfolding temperature; N130D/S437Y showed a 10 °C increase in optimal temperature. The secondary structure did not change much that contributed to improved thermostability were investigated in detail using circular dichroism. Homology modeling suggested that enhanced thermostability and specific activity may result from favorable hydrophobic interactions.

Conclusions

A series of mutations were achieved, showing higher thermostability and activity than the wild-type enzyme by semi-rational mutagenesis with limited structure information. Our findings provide important new insights into molecular modifications aimed at improving Ana-rLOX thermostability and activity.

Electronic supplementary material

The online version of this article (10.1186/s12896-018-0468-4) contains supplementary material, which is available to authorized users.

Review: Engineering of thermostable enzymes for industrial applications

The catalytic properties of some selected enzymes have long been exploited to carry out efficient and cost-effective bioconversions in a multitude of research and industrial sectors, such as food, health, cosmetics, agriculture, chemistry, energy, and others. Nonetheless, for several applications, naturally occurring enzymes are not considered to be viable options owing to their limited stability in the required working conditions. Over the years, the quest for novel enzymes with actual potential for biotechnological applications has involved various complementary approaches such as mining enzyme variants from organisms living in extreme conditions (extremophiles), mimicking evolution in the laboratory to develop more stable enzyme variants, and more recently, using rational, computer-assisted enzyme engineering strategies. In this review, we provide an overview of the most relevant enzymes that are used for industrial applications and we discuss the strategies that are adopted to enhance enzyme stability and/or activity, along with some of the most relevant achievements. In all living species, many different enzymes catalyze fundamental chemical reactions with high substrate specificity and rate enhancements. Besides specificity, enzymes also possess many other favorable properties, such as, for instance, cost-effectiveness, good stability under mild pH and temperature conditions, generally low toxicity levels, and ease of termination of activity. As efficient natural biocatalysts, enzymes provide great opportunities to carry out important chemical reactions in several research and industrial settings, ranging from food to pharmaceutical, cosmetic, agricultural, and other crucial economic sectors.

Enhanced Thermostability of Glucose Oxidase through Computer-Aided Molecular Design

Glucose oxidase (GOD, EC.1.1.3.4) specifically catalyzes the reaction of β-d-glucose to gluconic acid and hydrogen peroxide in the presence of oxygen, which has become widely used in the food industry, gluconic acid production and the feed industry. However, the poor thermostability of the current commercial GOD is a key limiting factor preventing its widespread application. In the present study, amino acids closely related to the thermostability of glucose oxidase from Penicillium notatum were predicted with a computer-aided molecular simulation analysis, and mutant libraries were established following a saturation mutagenesis strategy. Two mutants with significantly improved thermostabilities, S100A and D408W, were subsequently obtained. Their protein denaturing temperatures were enhanced by about 4.4 °C and 1.2 °C, respectively, compared with the wild-type enzyme. Treated at 55 °C for 3 h, the residual activities of the mutants were greater than 72%, while that of the wild-type enzyme was only 20%. The half-lives of S100A and D408W were 5.13- and 4.41-fold greater, respectively, than that of the wild-type enzyme at the same temperature. This work provides novel and efficient approaches for enhancing the thermostability of GOD by reducing the protein free unfolding energy or increasing the interaction of amino acids with the coenzyme.

Improving the temperature characteristics and catalytic efficiency of a mesophilic xylanase from Aspergillus oryzae, AoXyn11A, by iterative mutagenesis based on in silico design

To improve the temperature characteristics and catalytic efficiency of a glycoside hydrolase family (GHF) 11 xylanase from Aspergillus oryzae (AoXyn11A), its variants were predicted based on in silico design. Firstly, Gly²¹ with the maximum B-factor value, which was confirmed by molecular dynamics (MD) simulation on the three-dimensional structure of AoXyn11A, was subjected to site-saturation mutagenesis. Thus, one variant with the highest thermostability, AoXyn11A^G21I, was selected from the mutagenesis library, E. coli/Aoxyn11A^G21X (X: any one of 20 amino acids). Secondly, based on the primary structure multiple alignment of AoXyn11A with seven thermophilic GHF11 xylanases, AoXyn11A^Y13F or AoXyn11A^G21I–Y13F, was designed by replacing Tyr¹³ in AoXyn11A or AoXyn11A^G21I with Phe. Finally, three variant-encoding genes, Aoxyn11A^G21I, Aoxyn11A^Y13F and Aoxyn11A^G21I–Y13F, were constructed by two-stage whole-plasmid PCR method, and expressed in Pichia pastoris GS115, respectively. The temperature optimum (T_opt) of recombinant (re) AoXyn11A^G21I–Y13F was 60 °C, being 5 °C higher than that of reAoXyn11A^G21I or reAoXyn11A^Y13F, and 10 °C higher than that of reAoXyn11A. The thermal inactivation half-life (t_1/2) of reAoXyn11A^G21I–Y13F at 50 °C was 240 min, being 40-, 3.4- and 2.5-fold longer than those of reAoXyn11A, reAoXyn11A^G21I and reAoXyn11A^Y13F. The melting temperature (T_m) values of reAoXyn11A, reAoXyn11A^G21I, reAoXyn11A^Y13F and reAoXyn11A^G21I–Y13F were 52.3, 56.5, 58.6 and 61.3 °C, respectively. These findings indicated that the iterative mutagenesis of both Gly21Ile and Tyr13Phe improved the temperature characteristics of AoXyn11A in a synergistic mode. Besides those, the catalytic efficiency (k_cat/K_m) of reAoXyn11A^G21I–Y13F was 473.1 mL mg⁻¹ s⁻¹, which was 1.65-fold higher than that of reAoXyn11A.

Engineering and introduction of de novo disulphide bridges in organophosphorus hydrolase enzyme for thermostability improvement

The organophosphorus hydrolase (OPH) has been used to degrade organophosphorus chemicals, as one of the most frequently used decontamination methods. Under chemical and thermal denaturing conditions, the enzyme has been shown to unfold. To utilize this enzyme in various applications, the thermal stability is of importance. The engineering of de novo disulphide bridges has been explored as a means to increase the thermal stability of enzymes in the rational method of protein engineering. In this study, Disulphide by Design software, homology modelling and molecular dynamics simulations were used to select appropriate amino acid pairs for the introduction of disulphide bridge to improve protein thermostability. The thermostability of the wild-type and three selected mutant enzymes were evaluated by half-life, delta G inactivation (ΔGi) and structural studies (fluorescence and far-UV CD analysis). Data analysis showed that half-life of A204C/T234C and T128C/E153C mutants were increased up to 4 and 24 min, respectively; however, for the G74C/A78C mutant, the half-life was decreased up to 9 min. For the T128C/E124C mutant, both thermal stability and Catalytic efficiency (kcat) were also increased. The half-life and ΔGi results were correlated to the obtained information from structural studies by circular dichroism (CD) spectrometry and extrinsic fluorescence experiments; as rigidity increased in A204C/T2234C and T128C/E153C mutants, half-life and ΔGi also increased. For G74C/A78C mutant, these parameters decreased due to its higher flexibility. The results were submitted a strong evidence for the possibility to improve the thermostability of OPH enzyme by introducing a disulphide bridge after bioinformatics design, even though this design would not be always successful.

Improving the Thermostability of Acidic Pullulanase from Bacillus naganoensis by Rational Design

Pullulanase (EC 3.2.1.41) plays an important role in the specific hydrolysis of branch points in amylopectin. Enhancing its thermostability is required for its industrial application. In this study, rational protein design was used to improve the thermostability of PulB from Bacillus naganoensis (AB231790.1), which has strong enzymatic properties. Three positive single-site mutants (PulB-D328H, PulB-N387D, and PulB-A414P) were selected from six mutants. After incubation at 65°C for 5 min, the residual activities of PulB-D328H, PulB-N387D, and PulB-A414P were 4.5-, 1.7-, and 1.47-fold higher than PulB-WT, and their Tm values (the temperature at which half protein molecule denature) were 1.8°C, 0.4°C, and 0.9°C higher than PulB-WT, respectively. Then the final combined mutant PulB-328/387/414 was constructed. The t1/2 of it was 12.9-fold longer than that of PulB-WT at 65°C and the total increase in Tm of it (5.0°C) was almost 60% greater than the sum of individual increases (3.1°C). In addition, kinetic studies revealed that the kcat and the kcat/Km of PulB-328/387/414 increased by 38.8% and 12.9%. The remarkable improvement in thermostability and the high catalytic efficiency of PulB-328/387/414 make it suitable for industrial applications.

Improvement in the thermostability of a type A feruloyl esterase, AuFaeA, from Aspergillus usamii by iterative saturation mutagenesis

Feruloyl or ferulic acid esterase (Fae, EC 3.1.1.73) catalyzes the hydrolysis of ester bonds between polysaccharides and phenolic acid compounds in xylan side chain. In this study, the thermostability of a type A feruloyl esterase (AuFaeA) from Aspergillus usamii was increased by iterative saturation mutagenesis (ISM). Two amino acids, Ser33 and Asn92, were selected for saturation mutagenesis according to the B-factors analyzed by B-FITTER software and ΔΔG values predicted by PoPMuSiC algorithm. After screening the saturation mutagenesis libraries constructed in Pichia pastoris, 15 promising variants were obtained. The best variant S33E/N92-4 (S33E/N92R) produced a T m value of 44.5 °C, the half-lives (t1/2) of 35 and 198 min at 55 and 50 °C, respectively, corresponding to a 4.7 °C, 2.33- and 3.96-fold improvement compared to the wild type. Additionally, the best S33 variant S33-6 (S33E) was thermostable at 50 °C with a t1/2 of 82 min, which was 32 min longer than that of the wild type. All the screened S33E/N92 variants were more thermostable than the best S33 variant S33-6 (S33E). This work would contribute to the further studies on higher thermostability modification of type A feruloyl esterases, especially those from fungi. The thermostable feruloyl esterase variants were expected to be potential candidates for industrial application in prompting the enzymic degradation of plant biomass materials at elevated temperatures.

Thermal stabilization of dihydrofolate reductase using monte carlo unfolding simulations and its functional consequences

Design of proteins with desired thermal properties is important for scientific and biotechnological applications. Here we developed a theoretical approach to predict the effect of mutations on protein stability from non-equilibrium unfolding simulations. We establish a relative measure based on apparent simulated melting temperatures that is independent of simulation length and, under certain assumptions, proportional to equilibrium stability, and we justify this theoretical development with extensive simulations and experimental data. Using our new method based on all-atom Monte-Carlo unfolding simulations, we carried out a saturating mutagenesis of Dihydrofolate Reductase (DHFR), a key target of antibiotics and chemotherapeutic drugs. The method predicted more than 500 stabilizing mutations, several of which were selected for detailed computational and experimental analysis. We find a highly significant correlation of r = 0.65–0.68 between predicted and experimentally determined melting temperatures and unfolding denaturant concentrations for WT DHFR and 42 mutants. The correlation between energy of the native state and experimental denaturation temperature was much weaker, indicating the important role of entropy in protein stability. The most stabilizing point mutation was D27F, which is located in the active site of the protein, rendering it inactive. However for the rest of mutations outside of the active site we observed a weak yet statistically significant positive correlation between thermal stability and catalytic activity indicating the lack of a stability-activity tradeoff for DHFR. By combining stabilizing mutations predicted by our method, we created a highly stable catalytically active E. coli DHFR mutant with measured denaturation temperature 7.2°C higher than WT. Prediction results for DHFR and several other proteins indicate that computational approaches based on unfolding simulations are useful as a general technique to discover stabilizing mutations.

All-atom molecular simulations have provided valuable insight into the workings of molecular machines and the folding and unfolding of proteins. However, commonly employed molecular dynamics simulations suffer from a limitation in accessible time scale, making it difficult to model large-scale unfolding events in a realistic amount of simulation time without employing unrealistically high temperatures. Here, we describe a rapid all-atom Monte Carlo simulation approach to simulate unfolding of the essential bacterial enzyme Dihydrofolate Reductase (DHFR) and all possible single point-mutants. We use these simulations to predict which mutants will be more thermodynamically stable (i.e., reside more often in the native folded state vs. the unfolded state) than the wild-type protein, and we confirm our predictions experimentally, creating several highly stable and catalytically active mutants. Thermally stable active engineered proteins can be used as a starting point in directed evolution experiments to evolve new functions on the background of this additional “reservoir of stability.” The stabilized enzyme may be able to accumulate a greater number of destabilizing yet functionally important mutations before unfolding, protease digestion, and aggregation abolish its activity.

Switching a newly discovered lactonase into an efficient and thermostable phosphotriesterase by simple double mutations His250Ile/Ile263Trp

OPHC2 is a thermostable organophosphate (OP) hydrolase in the β-lactamase superfamily. OPs are highly toxic synthetic chemicals with no natural analogs. How did OPHC2 acquire phosphotriesterase (PTE) activity remained unclear. In this study, an OPHC2 analogue, PoOPH was discovered from Pseudomonas oleovorans exhibiting high lactonase and esterase activities and latent PTE activity. Sequence analysis revealed conserved His250 and Ile263 and site-directed mutagenesis at these crucial residues enhanced PTE activity. The best variant PoOPHM2 carrying H250I/I263W mutations displayed 6,962- and 106-fold improvements in catalytic efficiency for methyl-parathion and ethyl-paraoxon degradation, whereas the original lactonase and esterase activities decreased dramatically. A 1.4 × 10(7) -fold of specificity inversion was achieved by only two residue substitutions. Significantly, thermostability of the variants was not compromised. Crystal structure of PoOPHM2 was determined at 2.25 Å resolution and docking studies suggested that the two residues in the binding pocket determine substrate recognition. Lastly, new organophosphorus hydrolases (OPHs) were discovered using simple double mutations. Among them, PpOPHM2 from Pseudomonas putida emerged as a new promising OPH with very high activity (41.0 U mg(-1) ) toward methyl-parathion. Our results offer a first scrutiny to PTE activity evolution of OPHs in β-lactamase superfamily and provide efficient and robust enzymes for OP detoxification.

Improving the secretion of a methyl parathion hydrolase in Pichia pastoris by modifying its N-terminal sequence

Pichia pastoris is commonly used to express and secrete target proteins, although not all recombinant proteins can be successfully produced. In this study, we used methyl parathion hydrolase (MPH) from Ochrobactrum sp. M231 as a model to study the importance of the N-terminus of the protein for its secretion. While MPH can be efficiently expressed intracellularly in P. pastoris, it is not secreted into the extracellular environment. Three MPH mutants (N₆₆-MPH, D₁₀-MPH, and N₉-MPH) were constructed through modification of its N-terminus, and the secretion of each by P. pastoris was improved when compared to wild-type MPH. The level of secreted D₁₀-MPH was increased to 0.21 U/mL, while that of N₉-MPH was enhanced to 0.16 U/mL. Although N₆₆-MPH was not enzymatically active, it was secreted efficiently, and was identified by SDS-PAGE. These results demonstrate that the secretion of heterologous proteins in P. pastoris may be improved by modifying their N-terminal structures.

Computational tools for designing and engineering enzymes

Protein engineering strategies aimed at constructing enzymes with novel or improved activities, specificities, and stabilities greatly benefit from in silico methods. Computational methods can be principally grouped into three main categories: bioinformatics; molecular modelling; and de novo design. Particularly de novo protein design is experiencing rapid development, resulting in more robust and reliable predictions. A recent trend in the field is to combine several computational approaches in an interactive manner and to complement them with structural analysis and directed evolution. A detailed investigation of designed catalysts provides valuable information on the structural basis of molecular recognition, biochemical catalysis, and natural protein evolution.

A comparison of organophosphate degradation genes and bioremediation applications

Organophosphates (OPs) form the bulk of pesticides that are currently in use around the world accounting for more than 30% of the world market. They also form the core for many nerve-based warfare agents including sarin and soman. The widespread use and the resultant build-up of OP pesticides and chemical nerve agents has led to the development of major health problems due to their extremely toxic interaction with any biological system that encounters them. Growing concern over the accumulation of OP compounds in our food products, in the soils from which they are harvested and in wastewater run-off has fuelled a growing interest in microbial biotechnology that provides cheap, efficient OP detoxification to supplement expensive chemical methods. In this article, we review the current state of knowledge of OP pesticide and chemical agent degradation and attempt to clarify confusion over identification and nomenclature of two major families of OP-degrading enzymes through a comparison of their structure and function. The isolation, characterization, utilization and manipulation of the major detoxifying enzymes and the molecular basis of degradation of OP pesticides and chemical nerve agents are discussed as well as the achievements and technological advancements made towards the bioremediation of such compounds.

Structure- and sequence-analysis inspired engineering of proteins for enhanced thermostability

Protein engineering strategies for increasing stability can be improved by replacing random mutagenesis and high-throughput screening by approaches that include bioinformatics and computational design. Mutations can be focused on regions in the structure that are most flexible and involved in the early steps of thermal unfolding. Sequence analysis can often predict the position and nature of stabilizing mutations, and may allow the reconstruction of thermostable ancestral sequences. Various computational tools make it possible to design stabilizing features, such as hydrophobic clusters and surface charges. Different methods for designing chimeric enzymes can also support the engineering of more stable proteins without the need of high-throughput screening.