Thermally denatured state determines refolding in lipase: mutational analysis

Evolving Robust and Interpretable Enzymes for the Bioethanol Industry

Obtaining a robust and applicable enzyme for bioethanol production is a dream for biorefinery engineers. Herein, we describe a general method to evolve an all-round and interpretable enzyme that can be directly employed in the bioethanol industry. By integrating the transferable protein evolution strategy InSiReP 2.0 (In Silico guided Recombination Process), enzymatic characterization for actual production, and computational molecular understanding, the model cellulase PvCel5A (endoglucanase II Cel5A from Penicillium verruculosum) was successfully evolved to overcome the remaining challenges of low ethanol and temperature tolerance, which primarily limited biomass transformation and bioethanol yield. Remarkably, application of the PvCel5A variants in both first- and second-generation bioethanol production processes (i. Conventional corn ethanol fermentation combined with the in situ pretreatment process; ii. cellulosic ethanol fermentation process) resulted in a 5.7-10.1 % increase in the ethanol yield, which was unlikely to be achieved by other optimization techniques.

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Machine Learning Enables Selection of Epistatic Enzyme Mutants for Stability Against Unfolding and Detrimental Aggregation

Machine learning (ML) has pervaded most areas of protein engineering, including stability and stereoselectivity. Using limonene epoxide hydrolase as the model enzyme and innov'SAR as the ML platform, comprising a digital signal process, we achieved high protein robustness that can resist unfolding with concomitant detrimental aggregation. Fourier transform (FT) allows us to take into account the order of the protein sequence and the nonlinear interactions between positions, and thus to grasp epistatic phenomena. The innov'SAR approach is interpolative, extrapolative and makes outside-the-box, predictions not found in other state-of-the-art ML or deep learning approaches. Equally significant is the finding that our approach to ML in the present context, flanked by advanced molecular dynamics simulations, uncovers the connection between epistatic mutational interactions and protein robustness.

Esterification of Polymeric Carbohydrate Through Congener Cutinase-Like Biocatalyst

Cutinase-like enzymes (CLEs) are bi-functional hydrolases, which share the conserved catalytic site of lipase and consensus pentapeptide sequence of cutinase. Here, we have genetically replaced the canonical amino acids (CAA) by their non-canonical fluorinated surrogates to biosynthesize a novel class of congener biocatalyst for esterification of polymeric carbohydrate with long-chain fatty acid. It is a new enzyme-engineering approach used to manipulate industrially relevant biocatalyst through genetic incorporation of new functionally encoded non-canonical amino acids (NCAA). Global fluorination of CLE improved its catalytic, functional, and structural stability. Molecular docking studies confirmed that the fluorinated CLE (FCLE) had developed a binding affinity towards different fatty acids compared with the parent CLE. Importantly, FCLE could catalyze starch oleate synthesis in 24 h with a degree of substitution of 0.3 ± 0.001. Biophysical and microscopic analysis substantiated the efficient synthesis of the ester by FCLE. Our data represent the first step in the generation of an industrially relevant fluorous multifunctional enzyme for facile synthesis of high fatty acid starch esters.

Thermostability of Lipase A and Dynamic Communication Based on Residue Interaction Network

OBJECTIVE

Dynamic communication caused by mutation affects protein stability. The main objective of this study is to explore how mutations affect communication and to provide further insight into the relationship between heat resistance and signal propagation of Bacillus subtilis lipase (Lip A).

METHODS

The relationship between dynamic communication and Lip A thermostability is studied by long-time MD simulation and residue interaction network. The Dijkstra algorithm is used to get the shortest path of each residue pair. Subsequently, time-series frequent paths and spatio-temporal frequent paths are mined through an Apriori-like algorithm.

RESULTS

Time-series frequent paths show that the communication between residue pairs, both in wild-type lipase (WTL) and mutant 6B, becomes chaotic with an increase in temperature; however, more residues in 6B can maintain stable communication at high temperature, which may be associated with the structural rigidity. Furthermore, spatio-temporal frequent paths reflect the interactions among secondary structures. For WTL at 300K, β7, αC, αB, the longest loop, αA and αF contact frequently. The 310-helix between β3 and αA is penetrated by spatio-temporal frequent paths. At 400K, only αC can be frequently transmitted. For 6B, when at 300K, αA and αF are in more tight contact by spatio-temporal frequent paths though I157M and N166Y. Moreover, the rigidity of the active site His156 and the C-terminal of Lip A are increased, as reflected by the spatio-temporal frequent paths. At 400K, αA and αF, 310-helix between β3 and αA, the longest loop, and the loop where the active site Asp133 is located can still maintain stable communication.

CONCLUSION

From the perspective of residue dynamic communication, it is obviously found that mutations cause changes in interactions between secondary structures and enhance the rigidity of the structure, contributing to the thermal stability and functional activity of 6B.

Insight into the Structure and Activity of Surface-Engineered Lipase Biofluids

Despite a successful application of solvent-free liquid protein (biofluids) concept to a number of commercial enzymes, the technical advantages of enzyme biofluids as hyperthermal stable biocatalysts cannot be fully utilized as up to 90-99% of native activities are lost when enzymes were made into biofluids. With a two-step strategy (site-directed mutagenesis and synthesis of variant biofluids) on Bacillus subtilis lipase A (BsLA), we elucidated a strong dependency of structure and activity on the number and distribution of polymer surfactant binding sites on BsLA surface. Here, it is demonstrated that improved BsLA variants can be engineered via site-mutagenesis by a rational design, either with enhanced activity in aqueous solution in native form, or with improved physical property and increased activity in solvent-free system in the form of a protein liquid. This work answered some fundamental questions about the surface characteristics for construction of biofluids, useful for identifying new strategies for developing advantageous biocatalysts.

Structural and functional insights into thermostable and organic solvent stable variant Pro247-Ser of Bacillus lipase

Thermostability of enzymes is an important issue in protein engineering and has been studied in detail. Still there is no hard and fast rule to define the conditions which will provide thermal stability. Understanding the various factors and mechanism responsible for thermal stability will add on new insights into our present knowledge in this area. Pro247-Ser variant was constructed based on homology modelling and rational design. It exhibited 60 fold increase in thermal stability at 60°C and+0.7M shift in C_1/2 value for urea denaturation as compared to WT. Variant displayed noticeable tolerance to organic solvents. With decrease in K_m, catalytic efficiency of Pro247-Ser variant was increased by 12 fold. The activity and stability assay including circular dichroism and fluorescence spectroscopy favoured increased thermal performance of variant. Hydrolytic activity of variant was found to be high in comparison to control for all p-nitrophenol esters investigated. The immobilized variant enzyme demonstrated nearly two fold enhanced conversion of methyl oleate than WT enzyme. The additional molecular interactions of variant residue might contribute to increased thermostability of lipase. The homology modeling predicted formation of additional hydrogen bonds between Ser247/O-Thr251/OG1 as well as Ser247/O-Glu250/N.

Synthesis of paramagnetic dendritic silica nanomaterials with fibrous pore structure (Fe3O4@KCC-1) and their application in immobilization of lipase from Candida rugosa with enhanced catalytic activity and stability

Paramagnetic mesoporous fibrous silica (Fe₃O₄@KCC-1) was prepared and its surface was functionalized with 3-aminopropyltriethoxysilane (APTES).

Destabilizing an interacting motif strengthens the association of a designed ankyrin repeat protein with tubulin

Affinity maturation by random mutagenesis and selection is an established technique to make binding molecules more suitable for applications in biomedical research, diagnostics and therapy. Here we identified an unexpected novel mechanism of affinity increase upon in vitro evolution of a tubulin-specific designed ankyrin repeat protein (DARPin). Structural analysis indicated that in the progenitor DARPin the C-terminal capping repeat (C-cap) undergoes a 25° rotation to avoid a clash with tubulin upon binding. Additionally, the C-cap appears to be involved in electrostatic repulsion with tubulin. Biochemical and structural characterizations demonstrated that the evolved mutants achieved a gain in affinity through destabilization of the C-cap, which relieves the need of a DARPin conformational change upon tubulin binding and removes unfavorable interactions in the complex. Therefore, this specific case of an order-to-disorder transition led to a 100-fold tighter complex with a subnanomolar equilibrium dissociation constant, remarkably associated with a 30% decrease of the binding surface.

Application of Rigidity Theory to the Thermostabilization of Lipase A from Bacillus subtilis

Protein thermostability is a crucial factor for biotechnological enzyme applications. Protein engineering studies aimed at improving thermostability have successfully applied both directed evolution and rational design. However, for rational approaches, the major challenge remains the prediction of mutation sites and optimal amino acid substitutions. Recently, we showed that such mutation sites can be identified as structural weak spots by rigidity theory-based thermal unfolding simulations of proteins. Here, we describe and validate a unique, ensemble-based, yet highly efficient strategy to predict optimal amino acid substitutions at structural weak spots for improving a protein’s thermostability. For this, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state has been found as the cause for thermostability. When applied prospectively to lipase A from Bacillus subtilis, we achieved both a high success rate (25% over all experimentally tested mutations, which raises to 60% if small-to-large residue mutations and mutations in the active site are excluded) in predicting significantly thermostabilized lipase variants and a remarkably large increase in those variants’ thermostability (up to 6.6°C) based on single amino acid mutations. When considering negative controls in addition and evaluating the performance of our approach as a binary classifier, the accuracy is 63% and increases to 83% if small-to-large residue mutations and mutations in the active site are excluded. The gain in precision (predictive value for increased thermostability) over random classification is 1.6-fold (2.4-fold). Furthermore, an increase in thermostability predicted by our approach significantly points to increased experimental thermostability (p < 0.05). These results suggest that our strategy is a valuable complement to existing methods for rational protein design aimed at improving thermostability.

Protein thermostability is a crucial factor for biotechnological enzyme applications. However, performance studies of computational approaches for predicting effects of mutations on protein (thermo)stability have suggested that there is still room for improvement. We describe and validate a novel and unique strategy to predict optimal amino acid substitutions at structural weak spots. At variance with other rational approaches, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state is the underlying cause for thermostability. When applied prospectively on lipase LipA from Bacillus subtilis, a high success rate in predicting thermostabilized lipase variants and a remarkably large increase in their thermostability is achieved. This demonstrates the value of the novel strategy, which extends the existing portfolio of methods for rational protein design.

Strategies of preserving lipids in microorganism after fermentation

Microbial oil is accumulated by microorganisms as stored energy during cultivation, and will be degraded to generate energy when they are not able to obtain external energy source. As the lipid is the desired production, the in situ degradation of oil by microorganisms after fermentation or during downstream processing has to be prevented. This study investigated the effect of pH, thermal, salinity, bead milling, ultrasonication, and microwave treatment on the viability of Trichosporon oleaginosus when it was used for lipid production. The cells in broth were completely inactivated with the treatment of pH 1 and 2 for 1h, temperature 80 °C for 10 min, ultrasonication for 8 min, and microwave for 6 min, respectively. It was observed that these treatments had no impact on final product (biodiesel) composition and were considered as safe and efficient methods to preserve lipid in cells.

Structural Rigidity and Protein Thermostability in Variants of Lipase A from Bacillus subtilis

Understanding the origin of thermostability is of fundamental importance in protein biochemistry. Opposing views on increased or decreased structural rigidity of the folded state have been put forward in this context. They have been related to differences in the temporal resolution of experiments and computations that probe atomic mobility. Here, we find a significant (p = 0.004) and fair (R² = 0.46) correlation between the structural rigidity of a well-characterized set of 16 mutants of lipase A from Bacillus subtilis (BsLipA) and their thermodynamic thermostability. We apply the rigidity theory-based Constraint Network Analysis (CNA) approach, analyzing directly and in a time-independent manner the statics of the BsLipA mutants. We carefully validate the CNA results on macroscopic and microscopic experimental observables and probe for their sensitivity with respect to input structures. Furthermore, we introduce a robust, local stability measure for predicting thermodynamic thermostability. Our results complement work that showed for pairs of homologous proteins that raising the structural stability is the most common way to obtain a higher thermostability. Furthermore, they demonstrate that related series of mutants with only a small number of mutations can be successfully analyzed by CNA, which suggests that CNA can be applied prospectively in rational protein design aimed at higher thermodynamic thermostability.

Understanding the thermostability and activity of Bacillus subtilis lipase mutants: insights from molecular dynamics simulations

Improving the thermostability of industrial enzymes is an important protein engineering challenge. Point mutations, induced to increase thermostability, affect the structure and dynamics of the target protein in several ways and thus can also affect its activity. There appears to be no general rules for improving the thermostabilty of enzymes without adversely affecting their enzymatic activity. We report MD simulations, of wild type Bacillus subtilis lipase (WT) and its six progressively thermostable mutants (2M, 3M, 4M, 6M, 9M, and 12M), performed at different temperatures, to address this issue. Less thermostable mutants (LTMs), 2M to 6M, show WT-like dynamics at all simulation temperatures. However, the two more thermostable mutants (MTMs) show the required flexibility at appropriate temperature ranges and maintain conformational stability at high temperature. They show a deep and rugged free-energy landscape, confining them within a near-native conformational space by conserving noncovalent interactions, and thus protecting them from possible aggregation. In contrast, the LTMs having marginally higher thermostabilities than WT show greater probabilities of accessing non-native conformations, which, due to aggregation, have reduced possibilities of reverting to their respective native states under refolding conditions. Our analysis indicates the possibility of nonadditive effects of point mutations on the conformational stability of LTMs.

Thermostabilization of Bacillus subtilis lipase A by minimizing the structural deformation caused by packing enhancement

Enzyme thermostabilization is a critical research topic due to potential industrial benefits. Among the various reasons to increase enzyme thermostability, enhancement of residual packing at the core of the enzyme structure has been commonly accepted as a successful strategy. However, structural changes that occur with residual packing enhancement may decrease enzyme activity. In this study, a strategy to minimize structural deformation by calculating the overlapping packing volume of a single-point mutation followed by applying a double-point mutation was suggested. Four double mutants, A38V_K23A, A75V_T83A, G80A_N106A, and G172A_V100A, were selected for the in vitro experiment; three of the four showed enhancements in both thermostability and catalytic activity. In particular, G80A_N106A showed 2.78 times higher catalytic activity compared with wild type.

Efficient and robust analysis of biomacromolecular flexibility using ensembles of network topologies based on fuzzy noncovalent constraints

We describe an approach (ENT(FNC)) for performing rigidity analyses of biomacromolecules on ensembles of network topologies (ENT) generated from a single input structure. The ENT is based on fuzzy noncovalent constraints, which considers thermal fluctuations of biomacromolecules without actually sampling conformations. Definitions for fuzzy noncovalent constraints were derived from persistency data from molecular dynamics (MD) simulations. A very good agreement between local flexibility and rigidity characteristics from ENT(FNC) and MD simulations-generated ensembles is found. Regarding global characteristics, convincing results were obtained when relative thermostabilities of citrate synthase and lipase A structures were computed. The ENT(FNC) approach significantly improves the robustness of rigidity analyses, is highly efficient, and does not require a protein-specific parameterization. Its low computational demand makes it especially valuable for the analysis of large data sets, e.g., for data-driven protein engineering.

Improved thermal stability of lipase in W/O microemulsion by temperature-sensitive polymers

Lipase is active at the water-oil interface and thus very useful for many applications in non-aqueous media. However, the use of lipase is often limited due to the heat inactivation which is mainly caused by the irreversible aggregation among lipase molecules. The temperature-sensitive polymers can spontaneously form complexes with lipases at higher temperature in the confined spaces of the water in oil microemulsion. With cooling, lipases are released from the complexes and refold into the native state. In this way, the thermal stability of lipase in a microemulsion is effectively improved, and so is the stability of lipase at ambient temperature. Apart from proving the effectiveness and generality of this method, the temperature-sensitive polymers/lipase microemulsion represents a simple and efficient system which could be used in practical applications, since lipase retains the interfacial activity in this system. Moreover, the influences of some factors on the improvement are discussed and the mechanism of this method is suggested after exploring the process by dynamic light scattering and fluorescence measurements.

Role of active site rigidity in activity: MD simulation and fluorescence study on a lipase mutant

Relationship between stability and activity of enzymes is maintained by underlying conformational flexibility. In thermophilic enzymes, a decrease in flexibility causes low enzyme activity while in less stable proteins such as mesophiles and psychrophiles, an increase in flexibility is associated with enhanced enzyme activity. Recently, we identified a mutant of a lipase whose stability and activity were enhanced simultaneously. In this work, we probed the conformational dynamics of the mutant and the wild type lipase, particularly flexibility of their active site using molecular dynamic simulations and time-resolved fluorescence techniques. In contrast to the earlier observations, our data show that active site of the mutant is more rigid than wild type enzyme. Further investigation suggests that this lipase needs minimal reorganization/flexibility of active site residues during its catalytic cycle. Molecular dynamic simulations suggest that catalytically competent active site geometry of the mutant is relatively more preserved than wild type lipase, which might have led to its higher enzyme activity. Our study implies that widely accepted positive correlation between conformation flexibility and enzyme activity need not be stringent and draws attention to the possibility that high enzyme activity can still be accomplished in a rigid active site and stable protein structures. This finding has a significant implication towards better understanding of involvement of dynamic motions in enzyme catalysis and enzyme engineering through mutations in active site.

Biophysical characterization of mutants of Bacillus subtilis lipase evolved for thermostability: factors contributing to increased activity retention

Previously, Lipase A from Bacillus subtilis was subjected to in vitro directed evolution using iterative saturation mutagenesis, with randomization sites chosen on the basis of the highest B-factors available from the crystal structure of the wild-type (WT) enzyme. This provided mutants that, unlike WT enzyme, retained a large part of their activity after heating above 65 °C and cooling down. Here, we subjected the three best mutants along with the WT enzyme to biophysical and biochemical characterization. Combining thermal inactivation profiles, circular dichroism, X-ray structure analyses and NMR experiments revealed that mutations of surface amino acid residues counteract the tendency of Lipase A to undergo precipitation under thermal stress. Reduced precipitation of the unfolding intermediates rather than increased conformational stability of the evolved mutants seems to be responsible for the activity retention.

Probing protein stability and proteolytic resistance by loop scanning: a comprehensive mutational analysis

Improvement in protein thermostability was often found to be associated with increase in its proteolytic resistance as revealed by comparative studies of homologous proteins from extremophiles or mutational studies. Structural elements of protein responsible for this association are not firmly established although loops are implicated indirectly due to their structural role in protein stability. To get a better insight, a detailed study of protein wide mutants and their influence on stability and proteolytic resistance would be helpful. To generate such a data set, a model protein, Bacillus subtilis lipase was subjected to loop scanning site-saturation mutagenesis on 86 positions spanning all loops including termini. Upon screening of ~16,000 clones, 17 single mutants with improved thermostability were identified with increment in apparent melting temperature (Tm(app) ) by 1-6°C resulting in an increase in free energy of unfolding (ΔG(unf) ) by 0.04-1.16 kcal/mol. Proteolytic resistance of all single mutants upon incubation with nonspecific protease, Subtilisin A, was determined. Upon comparison, post-proteolysis residual activities as well as kinetics of proteolysis of mutants showed excellent correlation with ΔG(unf) , (r > 0.9), suggesting that proteolysis was strongly correlated with the global stability of this protein. This significant correlation in this set, with least possible sequence changes (single aa substitution), while covering >60% of protein surface strongly argues for the covariance of these two variables. Compared to studies from extremophiles, with large sequence heterogeneity, the observed correlation in such a narrow sequence space (ΔΔG(unf) = 1.57 kcal⁻¹) justifies the robustness of this relation.

In vitro evolved non-aggregating and thermostable lipase: structural and thermodynamic investigation

Rational and in vitro evolutionary approaches to improve either protein stability or aggregation resistance were successful, but empirical rules for simultaneous improvement of both stability and aggregation resistance under denaturing conditions are still to be ascertained. We have created a robust variant of a lipase from Bacillus subtilis named "6B" using multiple rounds of in vitro evolution. T(m) and optimum activity temperature of 6B is 78 °C and 65 °C, respectively, which is ~22 °C and 30 °C higher than that of wild-type lipase. Most significantly, 6B does not aggregate upon heating. Physical basis of remarkable thermostability and non-aggregating behavior of 6B was explored using X-ray crystallography, NMR and differential scanning calorimetry. Our structural investigations highlight the importance of tightening of mobile regions of the molecule such as loops and helix termini to attain higher thermostability. Accordingly, NMR studies suggest a very rigid structure of 6B lipase. Further investigation suggested that reduction/perturbation of the large hydrophobic patches present in the wild-type protein structure, decreased propensity of amino acid sequence for aggregation and absence of aggregation-prone intermediate during thermal unfolding of 6B can account for its resistance to aggregation. Overall, our study suggest that better anchoring of the loops with the rest of the protein molecule through mutations particularly on the sites that perturb/disturb the exposed hydrophobic patches can simultaneously increase protein stability and aggregation resistance.

Stability curves of laboratory evolved thermostable mutants of a Bacillus subtilis lipase

Shape of the protein stability curves changes to achieve higher melting temperature. Broadly, these changes have been classified as upward shift (increased G(s)), rightward shift (increase in T(s)) and flattening of the stability curves (decrease in C(p)). Comparative studies on homologous mesophilic-thermophilic protein pairs highlighted the differential contribution of these three strategies amongst proteins. But unambiguous way of identification of the strategies, which will be preferred for a protein, is still not achieved. We have performed comparative thermodynamic studies using differential scanning calorimeter (DSC) on thermostable variants of a lipase from Bacillus subtilis. These variants are products of 1, 2, 3 and 4 rounds of directed evolution and harbor mutations having definite contribution in thermostability unlike natural thermophilic proteins. We have shown that upward and rightward shift in stability curves are prime strategies in this lipase. Our results along with that from the other study on laboratory evolved xylanase A suggest that optimization of suboptimal thermodynamic parameters is having a dominant influence in selection of thermodynamic strategies for higher thermostability.