Application of molecular dynamics simulation in the field of food enzymes: improving the thermal-stability and catalytic ability

Enzymes can produce high-quality food with low pollution, high function, high acceptability, and medical aid. However, most enzymes, in their native form, do not meet the industrial requirements. Sequence-based and structure-based methods are the two main strategies used for enzyme modification. Molecular Dynamics (MD) simulation is a sufficiently comprehensive technology, from a molecular perspective, which has been widely used for structure information analysis and enzyme modification. In this review, we summarize the progress and development of MD simulation, particularly for software, force fields, and a standard procedure. Subsequently, we review the application of MD simulation in various food enzymes for thermostability and catalytic improvement was reviewed in depth. Finally, the limitations and prospects of MD simulation in food enzyme modification research are discussed. This review highlights the significance of MD simulation and its prospects in food enzyme modification.

Brave new surfactant world revisited by thermoalkalophilic lipases: computational insights into the role of SDS as a substrate analog

Non-ionic surfactants were shown to stabilize the active conformation of thermoalkalophilic lipases by mimicking the lipid substrate while the catalytic interactions formed by anionic surfactants have not been well characterized. In this study, we combined μs-scale molecular dynamics (MD) simulations and lipase activity assays to analyze the effect of ionic surfactant, sodium dodecyl sulfate (SDS), on the structure and activity of thermoalkalophilic lipases. Both the open and closed lipase conformations that differ in geometry were recruited to the MD analysis to provide a broader understanding of the molecular effect of SDS on the lipase structure. Simulations at 298 K showed the potential of SDS for maintaining the active lipase through binding to the sn-1 acyl-chain binding pocket in the open conformation or transforming the closed conformation to an open-like state. Consistent with MD findings, experimental analysis showed increased lipase activity upon SDS incubation at ambient temperature. Notably, the lipase cores stayed intact throughout 2 μs regardless of an increase in the simulation temperature or SDS concentration. However, the surface structures were unfolded in the presence of SDS and at elevated temperature for both conformations. Simulations of the dimeric lipase were also carried out and showed reduced flexibility of the surface structures which were unfolded in the monomer, indicating the insulating role of dimer interactions against SDS. Taken together, this study provides insights into the possible substrate mimicry by the ionic surfactant SDS for the thermoalkalophilic lipases without temperature elevation, underscoring SDS's potential for interfacial activation at ambient temperatures.

Thermostable lipases and their dynamics of improved enzymatic properties

Thermal stability is one of the most desirable characteristics in the search for novel lipases. The search for thermophilic microorganisms for synthesising functional enzyme biocatalysts with the ability to withstand high temperature, and capacity to maintain their native state in extreme conditions opens up new opportunities for their biotechnological applications. Thermophilic organisms are one of the most favoured organisms, whose distinctive characteristics are extremely related to their cellular constituent particularly biologically active proteins. Modifications on the enzyme structure are critical in optimizing the stability of enzyme to thermophilic conditions. Thermostable lipases are one of the most favourable enzymes used in food industries, pharmaceutical field, and actively been studied as potential biocatalyst in biodiesel production and other biotechnology application. Particularly, there is a trade-off between the use of enzymes in high concentration of organic solvents and product generation. Enhancement of the enzyme stability needs to be achieved for them to maintain their enzymatic activity regardless the environment. Various approaches on protein modification applied since decades ago conveyed a better understanding on how to improve the enzymatic properties in thermophilic bacteria. In fact, preliminary approach using advanced computational analysis is practically conducted before any modification is being performed experimentally. Apart from that, isolation of novel extremozymes from various microorganisms are offering great frontier in explaining the crucial native interaction within the molecules which could help in protein engineering. In this review, the thermostability prospect of lipases and the utility of protein engineering insights into achieving functional industrial usefulness at their high temperature habitat are highlighted. Similarly, the underlying thermodynamic and structural basis that defines the forces that stabilize these thermostable lipase is discussed. KEY POINTS: • The dynamics of lipases contributes to their non-covalent interactions and structural stability. • Thermostability can be enhanced by well-established genetic tools for improved kinetic efficiency. • Molecular dynamics greatly provides structure-function insights on thermodynamics of lipase.

Modifying the catalytic preference of alpha-amylase toward n-alkanes for bioremediation purposes using in silico strategies

Since the beginning of oil exploration, whole ecosystems have been affected by accidents and bad practices involving petroleum compounds. In this sense, bioremediation stands out as the cheapest and most eco-friendly alternatives to reverse the damage done in oil-impacted areas. However, more efforts must be made to engineer enzymes that could be used in the bioremediation process. Interestingly, a recent work described that α-amylase, one of the most evolutionary conserved enzymes, was able to promiscuously degrade n-alkanes, a class of molecules abundant in the petroleum admixture. Considering that α-amylase is expressed in almost all known organisms, and employed in numerous biotechnological processes, using it can be a great leap toward more efficient applications of enzyme or microorganism-consortia bioremediation approaches. In this work, we employed a strict computational approach to design new α-amylase mutants with potentially enhanced catalytic efficiency toward n-alkanes. Using in silico techniques, such as molecular docking, molecular dynamics, metadynamics, and residue-residue interaction networks, we generated mutants potentially more efficient for degrading n-alkanes, L183Y, and N314A. Our results indicate that the new mutants have an increased binding rate for tetradecane, the longest n-alkane previously tested, which can reside in the catalytic center for more extended periods. Additionally, molecular dynamics and network analysis showed that the new mutations have no negative impact on protein structure than the WT. Our results aid in solidifying this enzyme as one more tool in the petroleum bioremediation toolbox.

Effect of mutation at oxyanion hole residu (H110F) on activity of Lk4 lipase

Mutant of lipase at oxyanion hole (H110 F) was constructed. The gene was highly expressed in Eschericia coli BL21 (DE3) and the recombinant protein was purified using Ni-NTA affinity chromatography. The activity of mutant enzyme was significantly increased compared to that the wild type. Further comparison showed that both of the enzymes exhibited same optimum pH and temperature, and showed highest lipolytic activity on pNP-decanoate (C10). The wild type appeared lost of activity on C14 and C16 substrates meanwhile the mutant still showed activity up to 20 %. In the presence of non polar organic solvent such as n-hexane, the wild type became inactive enzyme meanwhile the mutant still remained 50 % of its activity. The results suggested that mutation at oxyanion hole (H110 F) caused enzyme-substrate interaction change resulting on elevation of activity, better activity toward longer carbon chain substrate and improving the activity in the present of non polar organic solvent.

Main Structural Targets for Engineering Lipase Substrate Specificity

Microbial lipases represent one of the most important groups of biotechnological biocatalysts. However, the high-level production of lipases requires an understanding of the molecular mechanisms of gene expression, folding, and secretion processes. Stable, selective, and productive lipase is essential for modern chemical industries, as most lipases cannot work in different process conditions. However, the screening and isolation of a new lipase with desired and specific properties would be time consuming, and costly, so researchers typically modify an available lipase with a certain potential for minimizing cost. Improving enzyme properties is associated with altering the enzymatic structure by changing one or several amino acids in the protein sequence. This review detailed the main sources, classification, structural properties, and mutagenic approaches, such as rational design (site direct mutagenesis, iterative saturation mutagenesis) and direct evolution (error prone PCR, DNA shuffling), for achieving modification goals. Here, both techniques were reviewed, with different results for lipase engineering, with a particular focus on improving or changing lipase specificity. Changing the amino acid sequences of the binding pocket or lid region of the lipase led to remarkable enzyme substrate specificity and enantioselectivity improvement. Site-directed mutagenesis is one of the appropriate methods to alter the enzyme sequence, as compared to random mutagenesis, such as error-prone PCR. This contribution has summarized and evaluated several experimental studies on modifying the substrate specificity of lipases.

Functional expression of a novel methanol-stable esterase from Geobacillus subterraneus DSM13552 for biocatalytic synthesis of cinnamyl acetate in a solvent-free system

BACKGROUND

Esterases are widely distributed in nature and have important applications in medical, industrial and physiological. Recently, the increased demand for flavor esters has prompted the search of catalysts like lipases and esterases. Esterases from thermophiles also show thermal stability at elevated temperatures and have become enzymes of special interest in biotechnological applications. Although most of esterases catalyzed reactions are carried out in toxic and inflammable organic solvents, the solvent-free system owning many advantages such as low cost and easy downstream processing.

RESULTS

The gene estGSU753 from Geobacillus subterraneus DSM13552 was cloned, sequenced and overexpressed into Escherichia coli BL21 (DE3). The novel gene has an open reading frame of 753 bp and encodes 250-amino-acid esterase (EstGSU753). The sequence analysis showed that the protein contains a catalytic triad formed by Ser97, Asp196 and His226, and the Ser of the active site is located in the conserved motif Gly95-X-Ser97-X-Gly99 included in most esterases and lipases. The protein catalyzed the hydrolysis of pNP-esters of different acyl chain lengths, and the enzyme specific activity was 70 U/mg with the optimum substrate pNP-caprylate. The optimum pH and temperature of the recombinant enzyme were 8.0 and 60 °C respectively. The resulting EstGSU753 showed remarkable stability against methanol. After the incubation at 50% methanol for 9 days, EstGSU753 retained 50% of its original activity. Even incubation at 90% methanol for 35 h, EstGSU753 retained 50% of its original activity. Also, the preliminary study of the transesterification shows the potential value in synthesis of short-chain flavor esters in a solvent-free system, and more than 99% conversion was obtained in 6 h (substrate: cinnamyl alcohol, 1.0 M).

CONCLUSIONS

This is the first report of esterase gene cloning from Geobacillus subterraneus with detailed enzymatic properties. This methanol-stable esterase showed potential value in industrial applications especially in the perfume industry.

Understanding thermal and organic solvent stability of thermoalkalophilic lipases: insights from computational predictions and experiments

Bacillus thermocatenulatus lipase (BTL2), a member of the isolated lipase family known as thermoalkalophilic lipases, carries potential for industrial applications owing to its ability to catalyze versatile reactions under extreme conditions. This study investigates the molecular effects of distinct solvents on the stability of BTL2 at different temperatures, aiming to contribute to lipase use in industrial applications. Initially, molecular dynamic (MD) simulations were carried out to address for the molecular impacts of distinct solvents on the structural stability of BTL2 at different temperatures. Two lipase conformations representing the active and inactive forms were simulated in 5 solvents including water, ethanol, methanol, cyclohexane, and toluene. Low temperature simulations showed that polar solvents led to enhanced lid fluctuations compared with non-polar solvents reflecting a more dynamic equilibrium between active and inactive lipase conformations in polar solvents including water, while the overall structure of the lipase in both forms became more rigid in non-polar solvents than they were in polar solvent. Notably, the native lipase fold was maintained in non-polar solvents even at high temperatures, indicating an enhancement of lipase's thermostability in non-polar organic solvents. Next, we conducted experiments for which BTL2 was expressed in a heterologous host and purified to homogeneity, and its thermostability in different solvents was assessed. Parallel to the computational findings, experimental results suggested that non-polar organic solvents contributed to BTL2's thermostability at concentrations as high as 70% (v/v). Altogether, this study provides beneficial insights to the lipase use under extreme conditions. Graphical Abstract.

Computational tools for enzyme improvement: why everyone can - and should - use them

This review presents computational methods that experimentalists can readily use to create smart libraries for enzyme engineering and to obtain insights into protein-substrate complexes. Computational tools have the reputation of being hard to use and inaccurate compared to experimental methods in enzyme engineering, yet they are essential to probe datasets of ever-increasing size and complexity. In recent years, bioinformatics groups have made a huge leap forward in providing user-friendly interfaces and accurate algorithms for experimentalists. These methods guide efficient experimental planning and allow the enzyme engineer to rationalize time and resources. Computational tools nevertheless face challenges in the realm of transient modern technology.

Predicting the impact of mutations on the specific activity of Bacillus thermocatenulatus lipase using a combined approach of docking and molecular dynamics

Lipases are important biocatalysts owing to their ability to catalyze diverse reactions with exceptional substrate specificities. A combined docking and molecular dynamics (MD) approach was applied to study the chain-length selectivity of Bacillus thermocatenulatus lipase (BTL2) towards its natural substrates (triacylglycerols). A scoring function including electrostatic, van der Waals (vdW) and desolvation energies along with conformational entropy was developed to predict the impact of mutation. The native BTL2 and its 6 mutants (F17A, V175A, V175F, D176F, T178V and I320F) were experimentally analyzed to determine their specific activities towards tributyrin (C4) or tricaprylin (C8), which were used to test our approach. Our scoring methodology predicted the chain-length selectivity of BTL2 with 85.7% (6/7) accuracy with a positive correlation between the calculated scores and the experimental activity values (r = 0.82, p = 0.0004). Additionally, the impact of mutation on activity was predicted with 75% (9/12) accuracy. The described study represents a fast and reliable approach to accurately predict the effect of mutations on the activity and selectivity of lipases and also of other enzymes. Copyright © 2016 John Wiley & Sons, Ltd.

Zinc Modulates Self-Assembly of Bacillus thermocatenulatus Lipase

Thermoalkalophilic lipases are prone to aggregation from their dimer interface to which structural zinc is very closely located. Structural zinc sites have been shown to induce protein aggregation, but the interaction between zinc and aggregation tendency in thermoalkalophilic lipases remains elusive. Here we delineate the interplay between zinc and aggregation of the lipase from Bacillus thermocatenulatus (BTL2), which is taken to be a representative of thermoalkalophilic lipase. Results showed that zinc removal disrupted the BTL2 dimer, leading to monomer formation and reduced thermostability manifesting as a link between zinc and dimerization that leads to thermostability, while zinc addition induced aggregation. Biochemical and kinetic characterizations of zinc-induced aggregates showed that the aggregates obtained from the early and late stages of aggregation had differential characteristics. In the early stages, the aggregates were soluble and possessed native-like structures, while in the late stages, the aggregates became insoluble and showed fibrillar characteristics with binding affinities for Congo red and thioflavin T. The impact of temperature on zinc-induced aggregation was further investigated, and it was found that the native-like early aggregates could completely dissociate into functional lipase forms at high temperatures while dissociation of the late aggregates was limited. To this end, we report that the zinc-induced aggregation of BTL2 can be reversed by temperature switches and initiated by ordered aggregates in the early stages that gain fibrillar-like features over time. Insights revealed by this work contributes to the knowledge of aggregation mechanisms that exist in thermophilic proteins, reflecting the potential use of metal addition and/or removal to fine-tune aggregation tendency.

Adaptational properties and applications of cold-active lipases from psychrophilic bacteria

Psychrophilic microorganisms are cold-adapted with distinct properties from other thermal classes thriving in cold conditions in large areas of the earth's cold environment. Maintenance of functional membranes, evolving cold-adapted enzymes and synthesizing a range of structural features are basic adaptive strategies of psychrophiles. Among the cold-evolved enzymes are the cold-active lipases, a group of microbial lipases with inherent stability-activity-flexibility property that have engaged the interest of researchers over the years. Current knowledge regarding these cold-evolved enzymes in psychrophilic bacteria proves a display of high catalytic efficiency with low thermal stability, which is a differentiating feature with that of their mesophilic and thermophilic counterparts. Improvement strategies of their adaptive structural features have significantly benefited the enzyme industry. Based on their homogeneity and purity, molecular characterizations of these enzymes have been successful and their properties make them unique biocatalysts for various industrial and biotechnological applications. Although, strong association of lipopolysaccharides from Antarctic microorganisms with lipid hydrolases pose a challenge in their purification, heterologous expression of the cold-adapted lipases with affinity tags simplifies purification with higher yield. The review discusses these cold-evolved lipases from bacteria and their peculiar properties, in addition to their potential biotechnological and industrial applications.

A review of metabolic and enzymatic engineering strategies for designing and optimizing performance of microbial cell factories

Microbial cell factories (MCFs) are of considerable interest to convert low value renewable substrates to biofuels and high value chemicals. This review highlights the progress of computational models for the rational design of an MCF to produce a target bio-commodity. In particular, the rational design of an MCF involves: (i) product selection, (ii) de novo biosynthetic pathway identification (i.e., rational, heterologous, or artificial), (iii) MCF chassis selection, (iv) enzyme engineering of promiscuity to enable the formation of new products, and (v) metabolic engineering to ensure optimal use of the pathway by the MCF host. Computational tools such as (i) de novo biosynthetic pathway builders, (ii) docking, (iii) molecular dynamics (MD) and steered MD (SMD), and (iv) genome-scale metabolic flux modeling all play critical roles in the rational design of an MCF. Genome-scale metabolic flux models are of considerable use to the design process since they can reveal metabolic capabilities of MCF hosts. These can be used for host selection as well as optimizing precursors and cofactors of artificial de novo biosynthetic pathways. In addition, recent advances in genome-scale modeling have enabled the derivation of metabolic engineering strategies, which can be implemented using the genomic tools reviewed here as well.

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.

The conserved lid tryptophan, W211, potentiates thermostability and thermoactivity in bacterial thermoalkalophilic lipases

We hypothesize that aggregation of thermoalkalophilic lipases could be a thermostability mechanism. The conserved tryptophans (W211, W234) in the lid are of particular interest owing to their previous involvements in aggregation and thermostability mechanisms in many other proteins. The thermoalkalophilic lipase from Bacillus thermocatenulatus (BTL2) and its mutants (W211A, W234A) were expressed and purified to homogeneity. We found that, when aggregated, BTL2 is more thermostable than its non-aggregating form, showing that aggregation potentiates thermostability in the thermoalkalophilic lipase. Among the two lid mutants, the W211A lowered aggregation tendency drastically and resulted in a much less thermostable variant of BTL2, which indicated that W211 stabilizes the intermolecular interactions in BTL2 aggregates. Further thermoactivity and CD spectroscopy analyses showed that W211A also led to a strong decrease in the optimal and the melting temperature of BTL2, implying stabilization by W211 also to the intramolecular interactions. The other lid mutant W234A had no effects on these properties. Finally, we analyzed the molecular basis of these experimental findings in-silico using the dimer (PDB ID: 1KU0) and the monomer (PDB ID: 2W22) lipase structures. The computational analyses confirmed that W211 stabilized the intermolecular interactions in the dimer lipase and it is critical to the stability of the monomer lipase. Explicitly W211 confers stability to the dimer and the monomer lipase through distinct aromatic interactions with Y273-Y282 and H87-P232 respectively. The insights revealed by this work shed light not only on the mechanism of thermostability and its relation to aggregation but also on the particular role of the conserved lid tryptophan in the thermoalkalophilic lipases.