Modeling hydration mechanisms of enzymes in nonpolar and polar organic solvents

Water-mediated active conformational transitions of lipase on organic solvent interfaces

When it comes to enzyme stability and their application in organic solvents, enzyme biocatalysis has emerged as a popular substitute for conventional chemical processes. However, the demand for enzymes exhibiting improved stability remains a persistent challenge. Organic solvents can significantly impacts enzyme properties, thereby limiting their practical application. This study focuses on Lipase Thermomyces lanuginose, through molecular dynamics simulations and experiments, we quantified the effect of different solvent-lipase interfaces on the interfacial activation of lipase. Revealed molecular views of the complex solvation processes through the minimum distance distribution function. Solvent-protein interactions were used to interpret the factors influencing changes in lipase conformation and enzyme activity. We found that water content is crucial for enzyme stability, and the optimum water content for lipase activity was 35 % in the presence of benzene-water interface, which is closely related to the increase of its interfacial activation angle from 78° to 102°. Methanol induces interfacial activation in addition to significant competitive inhibition and denaturation at low water content. Our findings shed light on the importance of understanding solvent effects on enzyme function and provide practical insights for enzyme engineering and optimization in various solvent-lipase interfaces.

A Study on the Regioselective Acetylation of Flavonoid Aglycons Catalyzed by Immobilized Lipases

This study aimed to explore the capacity of immobilized lipases on the acetylation of six aglycon flavonoids, namely myricetin, quercetin, luteolin, naringenin, fisetin and morin. For this purpose, lipase B from Candida antarctica (CaLB) and lipase from Thermomyces lanuginosus (TLL) were immobilized onto the surface of ZnOFe nanoparticles derived from an aqueous olive leaf extract. Various factors affecting the conversion of substrates and the formation of monoesterified and diesterified products, such as the amount of biocatalyst and the molar ratio of the substrates and reaction solvents were investigated. Both CaLB and TLL-ZnOFe achieved 100% conversion yield of naringenin to naringenin acetate after 72 h of reaction time, while TLL-ZnOFe achieved higher conversion yields of quercetin, morin and fisetin (73, 85 and 72% respectively). Notably, CaLB-ZnOFe displayed significantly lower conversion yields for morin compared with TLL-ZnOFe. Molecular docking analysis was used to elucidate this discrepancy, and it was revealed that the position of the hydroxyl groups of the B ring on morin introduced hindrances on the active site of CaLB. Finally, selected flavonoid esters showed significantly higher antimicrobial activity compared with the original compound. This work indicated that these lipase-based nanobiocatalysts can be successfully applied to produce lipophilic derivatives of aglycon flavonoids with improved antimicrobial activity.

Investigating Biomolecules in Deep Eutectic Solvents with Molecular Dynamics Simulations: Current State, Challenges and Future Perspectives

Deep eutectic solvents (DESs) have recently gained increased attention for their potential in biotechnological applications. DESs are binary mixtures often consisting of a hydrogen bond acceptor and a hydrogen bond donor, which allows for tailoring their properties for particular applications. If produced from sustainable resources, they can provide a greener alternative to many traditional organic solvents for usage in various applications (e.g., as reaction environment, crystallization agent, or storage medium). To navigate this large design space, it is crucial to comprehend the behavior of biomolecules (e.g., enzymes, proteins, cofactors, and DNA) in DESs and the impact of their individual components. Molecular dynamics (MD) simulations offer a powerful tool for understanding thermodynamic and transport processes at the atomic level and offer insights into their fundamental phenomena, which may not be accessible through experiments. While the experimental investigation of DESs for various biotechnological applications is well progressed, a thorough investigation of biomolecules in DESs via MD simulations has only gained popularity in recent years. Within this work, we aim to provide an overview of the current state of modeling biomolecules with MD simulations in DESs and discuss future directions with a focus for optimizing the molecular simulations and increasing our fundamental knowledge.

Influence of elevated pressure and pressurized fluids on microenvironment and activity of enzymes

Enzymes have great potential in bioprocess engineering due to their green and mild reaction conditions. However, there are challenges to their application, such as enzyme extraction and purification costs, enzyme recovery, and long reaction time. Enzymatic reaction rate enhancement and enzyme immobilization have the potential to overcome some of these challenges. Application of high pressure (e.g., hydrostatic pressure, supercritical carbon dioxide) has been shown to increase the activity of some enzymes, such as lipases and cellulases. Under high pressure, enzymes undergo multiple alterations simultaneously. High pressure reduces the bond lengths of molecules of reaction components and causes a reduction in the activation volume of enzyme-substrate complex. Supercritical CO2 interacts with enzyme molecules, catalyzes structural changes, and removes some water molecules from the enzyme's hydration layer. Interaction of scCO2 with the enzyme also leads to an overall change in secondary structure content. In the extreme, such changes may lead to enzyme denaturation, but enzyme activation and stabilization have also been observed. Immobilization of enzymes onto silica and zeolite-based supports has been shown to further stabilize the enzyme and provide resistance towards perturbation under subjection to high pressure and scCO2.

Solvent induced conformational changes for the altered activity of laccase: A molecular dynamics study

The growing demands of solvent-based industries like paint, pharmaceutical, petrochemical, paper and pulp, etc., have directly increased the release of effluents that are rich in hazardous aromatic compounds in the environment. A sustainable biotechnological approach utilizing laccases as biocatalyst enable in biodegradation of these aromatic toxin-rich effluents. However, this enzymatic process is ineffective as laccases lose their stability and catalytic activity at high organic solvent concentrations. In this study, molecular dynamic simulations of a novel solvent tolerant laccase, DLac from Cerrena sp. RSD1 was performed to explore the molecular-level understanding of DLac in 30%(v/v) acetone and acetonitrile. Solvent-induced conformational changes were analyzed via protein structure network, which was illustrated with respect to cliques and communities. In the presence of acetonitrile, the cliques around the active site and substrate-binding site were disjoined, thus the communities lost their network integrity. Whereas with acetone, the community near the substrate-binding site gained new residues and formed a rigidified network that corresponded to enhanced DLac's activity. Moreover, prominent solvent binding sites were speculated, which can be probable mutation targets to further improve solvent tolerance and catalytic activity. The molecular basis behind solvent induced catalytic activity will further aid in engineering laccase for its industrial application.

Using Molecular Simulation to Guide Protein Engineering for Biocatalysis in Organic Solvents

Biocatalysis in organic solvents (OSs) is very appealing for the industry in producing bulk and/or fine chemicals, such as pharmaceuticals, biodiesel, and fragrances. The poor performance of enzymes in OSs (e.g., reduced activity, insufficient stability, and deactivation) negates OSs' excellent solvent properties. Molecular dynamics (MD) simulations provide a complementary method to study the relationship between enzymes dynamics and the stability in OSs. Here we describe computational procedure for MD simulation of enzymes in OSs with an example of Bacillus subtilis lipase A (BSLA) in dimethyl sulfoxide (DMSO) cosolvent with software GROMACS. We discuss main essential practical issues considered (such as choice of force field, parameterization, simulation setup, and trajectory analysis). The core part of this protocol (enzyme-OS system setup, analysis of structural-based and solvation-based observables) is transferable to other enzymes and any OS systems. Combining with experimental studies, the obtained molecular knowledge is most likely to guide researchers to access rational protein engineering approaches to tailor OS resistant enzymes and expand the scope of biocatalysis in OS media. Finally, we discuss potential solutions to overcome the remaining challenges of computational biocatalysis in OSs and briefly draw future directions for further improvement in this field.

Efficient Enzymatic Synthesis of Lipophilic Phenolic Glycoside Azelaic Acid Esters and Their Depigmenting Activity

In this paper, an enzymatic route for synthesizing phenolic glycoside azelaic acid esters was successfully set up via lipase-catalyzed esterification and transesterification. Among the lipases tested, Candida antarctica lipase B (Novozyme 435) showed the highest activity in catalyzing esterification and Thermomyces lanuginosus (Lipozyme TLIM) gave the highest substrate conversion in catalyzing transesterification for the synthesis of ester. The addition of 4A molecular sieves into the reaction system is found to be an effective method for in situ absorption of the byproduct water and methanol, with which the substrate conversions of the enzymatic esterification and transesterification were 98.7 and 95.1%, respectively. Also, the main product ratios in transesterification were above 99.0% with lipozyme TLIM as a catalyst because the hydrolysis reaction was hindered. The results of the physical and biological properties indicate that all esters had higher Clog p values than their parent compounds. Also, the esters showed higher intracellular tyrosinase inhibitory and depigmentating activities than phenolic glycosides, azelaic acid (AA), and their physical mixtures due to their higher membrane penetration and tyrosinase inhibitory effects. In particular, piceid 6″-O-azelaic acid ester (PIA) showed the strongest inhibitory effect against melanin production. Its inhibitory rate was 77.4% at a concentration of 0.25 mM, about 4.2 times higher than that of arbutin (18.5%).

Green Synthesized Magnetic Nanoparticles as Effective Nanosupport for the Immobilization of Lipase: Application for the Synthesis of Lipophenols

In this work, hybrid zinc oxide-iron oxide (ZnOFe) magnetic nanoparticles were synthesized employing Olea europaea leaf aqueous extract as a reducing/chelating and capping medium. The resulting magnetic nanoparticles were characterized by basic spectroscopic and microscopic techniques, namely, X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), fourier-transform infrared (FTIR) and atomic force microscopy (AFM), exhibiting a spherical shape, average size of 15-17 nm, and a functionalized surface. Lipase from Thermomyces lanuginosus (TLL) was efficiently immobilized on the surface of ZnOFe nanoparticles through physical absorption. The activity of immobilized lipase was found to directly depend on the enzyme to support the mass ratio, and also demonstrated improved pH and temperature activity range compared to free lipase. Furthermore, the novel magnetic nanobiocatalyst (ZnOFe-TLL) was applied to the preparation of hydroxytyrosyl fatty acid esters, including derivatives with omega-3 fatty acids, in non-aqueous media. Conversion yields up to 90% were observed in non-polar solvents, including hydrophobic ionic liquids. Different factors affecting the biocatalyst performance were studied. ZnOFe-TLL was reutilized for eight subsequent cycles, exhibiting 90% remaining esterification activity (720 h of total operation at 50 °C). The green synthesized magnetic nanoparticles, reported here for the first time, are excellent candidates as nanosupports for the immobilization of enzymes with industrial interest, giving rise to nanobiocatalysts with elevated features.

Advances in deep eutectic solvents and water: applications in metal- and biocatalyzed processes, in the synthesis of APIs, and other biologically active compounds

This review highlights recent advances in metal- and biocatalyzed transformations, in the synthesis of APIs and other biologically active compounds, when employing deep eutectic solvents and water as environmentally responsible solvents.

What do we learn from enzyme behaviors in organic solvents? - Structural functionalization of ionic liquids for enzyme activation and stabilization

Enzyme activity in nonaqueous media (e.g. conventional organic solvents) is typically lower than in water by several orders of magnitude. There is a rising interest of developing new nonaqueous solvent systems that are more "water-like" and more biocompatible. Therefore, we need to learn from the current state of nonaqueous biocatalysis to overcome its bottleneck and provide guidance for new solvent design. This review firstly focuses on the discussion of how organic solvent properties (such as polarity and hydrophobicity) influence the enzyme activity and stability, and how these properties impact the enzyme's conformation and dynamics. While hydrophobic organic solvents usually lead to the maintenance of enzyme activity, solvents carrying functional groups like hydroxys and ethers (including crown ethers and cyclodextrins) can lead to enzyme activation. Ionic liquids (ILs) are designable solvents that can conveniently incorporate these functional groups. Therefore, we systematically survey these ether- and/or hydroxy-functionalized ILs, and find most of them are highly compatible with enzymes leading to high activity and stability. In particular, ILs carrying both ether and tert-alcohol groups are among the most enzyme-activating solvents. Future direction is to learn from enzyme behaviors in both water and nonaqueous media to design biocompatible "water-like" solvents.

Modeling Alcohol Dehydrogenase Catalysis in Deep Eutectic Solvent/Water Mixtures

The use of oxidoreductases (EC1) in non-conventional reaction media has been increasingly explored. In particular, deep eutectic solvents (DESs) have emerged as a novel class of solvents. Herein, an in-depth study of bioreduction with an alcohol dehydrogenase (ADH) in the DES glyceline is presented. The activity and stability of ADH in mixtures of glyceline/water with varying water contents were measured. Furthermore, the thermodynamic water activity and viscosity of mixtures of glyceline/water have been determined. For a better understanding of the observations, molecular dynamics simulations were performed to quantify the molecular flexibility, hydration layer, and intraprotein hydrogen bonds of ADH. The behavior of the enzyme in DESs follows the classic dependence of water activity (aW ) in non-conventional media. At low aW values (<0.2), ADH does not show any activity; at higher aW values, the activity was still lower than that in pure water due to the high viscosities of the DES. These findings could be further explained by increased enzyme flexibility with increasing water content.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

Lid opening and conformational stability of T1 Lipase is mediated by increasing chain length polar solvents

The dynamics and conformational landscape of proteins in organic solvents are events of potential interest in nonaqueous process catalysis. Conformational changes, folding transitions, and stability often correspond to structural rearrangements that alter contacts between solvent molecules and amino acid residues. However, in nonaqueous enzymology, organic solvents limit stability and further application of proteins. In the present study, molecular dynamics (MD) of a thermostable Geobacillus zalihae T1 lipase was performed in different chain length polar organic solvents (methanol, ethanol, propanol, butanol, and pentanol) and water mixture systems to a concentration of 50%. On the basis of the MD results, the structural deviations of the backbone atoms elucidated the dynamic effects of water/organic solvent mixtures on the equilibrium state of the protein simulations in decreasing solvent polarity. The results show that the solvent mixture gives rise to deviations in enzyme structure from the native one simulated in water. The drop in the flexibility in H₂O, MtOH, EtOH and PrOH simulation mixtures shows that greater motions of residues were influenced in BtOH and PtOH simulation mixtures. Comparing the root mean square fluctuations value with the accessible solvent area (SASA) for every residue showed an almost correspondingly high SASA value of residues to high flexibility and low SASA value to low flexibility. The study further revealed that the organic solvents influenced the formation of more hydrogen bonds in MtOH, EtOH and PrOH and thus, it is assumed that increased intraprotein hydrogen bonding is ultimately correlated to the stability of the protein. However, the solvent accessibility analysis showed that in all solvent systems, hydrophobic residues were exposed and polar residues tended to be buried away from the solvent. Distance variation of the tetrahedral intermediate packing of the active pocket was not conserved in organic solvent systems, which could lead to weaknesses in the catalytic H-bond network and most likely a drop in catalytic activity. The conformational variation of the lid domain caused by the solvent molecules influenced its gradual opening. Formation of additional hydrogen bonds and hydrophobic interactions indicates that the contribution of the cooperative network of interactions could retain the stability of the protein in some solvent systems. Time-correlated atomic motions were used to characterize the correlations between the motions of the atoms from atomic coordinates. The resulting cross-correlation map revealed that the organic solvent mixtures performed functional, concerted, correlated motions in regions of residues of the lid domain to other residues. These observations suggest that varying lengths of polar organic solvents play a significant role in introducing dynamic conformational diversity in proteins in a decreasing order of polarity.

To Keep or Not to Keep? The Question of Crystallographic Waters for Enzyme Simulations in Organic Solvent

The use of enzymes in non-aqueous solvents expands the use of biocatalysts to hydrophobic substrates, with the ability to tune selectivity of reactions through solvent selection. Non-aqueous enzymology also allows for fundamental studies on the role of water and other solvents in enzyme structure, dynamics, and function. Molecular dynamics simulations serve as a powerful tool in this area, providing detailed atomic information about the effect of solvents on enzyme properties. However, a common protocol for non-aqueous enzyme simulations does not exist. If you want to simulate enzymes in non-aqueous solutions, how many and which crystallographic waters do you keep? In the present work, this question is addressed by determining which crystallographic water molecules lead most quickly to an equilibrated protein structure. Five different methods of selecting and keeping crystallographic waters are used in order to discover which crystallographic waters lead the protein structure to reach an equilibrated structure more rapidly in organic solutions. It is found that buried waters contribute most to rapid equilibration in organic solvent, with slow-diffusing waters giving similar results.

Organic co-solvents affect activity, stability and enantioselectivity of haloalkane dehalogenases

Haloalkane dehalogenases are microbial enzymes with a wide range of biotechnological applications, including biocatalysis. The use of organic co-solvents to solubilize their hydrophobic substrates is often necessary. In order to choose the most compatible co-solvent, the effects of 14 co-solvents on activity, stability and enantioselectivity of three model enzymes, DbjA, DhaA, and LinB, were evaluated. All co-solvents caused at high concentration loss of activity and conformational changes. The highest inactivation was induced by tetrahydrofuran, while more hydrophilic co-solvents, such as ethylene glycol and dimethyl sulfoxide, were better tolerated. The effects of co-solvents at low concentration were different for each enzyme-solvent pair. An increase in DbjA activity was induced by the majority of organic co-solvents tested, while activities of DhaA and LinB decreased at comparable concentrations of the same co-solvent. Moreover, a high increase of DbjA enantioselectivity was observed. Ethylene glycol and 1,4-dioxane were shown to have the most positive impact on the enantioselectivity. The favorable influence of these co-solvents on both activity and enantioselectivity makes DbjA suitable for biocatalytic applications. This study represents the first investigation of the effects of organic co-solvents on the biocatalytic performance of haloalkane dehalogenases and will pave the way for their broader use in industrial processes.

Effects of water content on the tetrahedral intermediate of chymotrypsin - trifluoromethylketone in polar and non-polar media: observations from molecular dynamics simulation

The work uses MD simulation to study effects of five water contents (3 %, 10 %, 20 %, 50 %, 100 % v/v) on the tetrahedral intermediate of chymotrypsin - trifluoromethyl ketone in polar acetonitrile and non-polar hexane media. The water content induced changes in the structure of the intermediate, solvent distribution and H-bonding are analyzed in the two organic media. Our results show that the changes in overall structure of the protein almost display a clear correlation with the water content in hexane media while to some extent U-shaped/bell-shaped dependence on the water content is observed in acetonitrile media with a minimum/maximum at 10-20 % water content. In contrast, the water content change in the two organic solvents does not play an observable role in the stability of catalytic hydrogen-bond network, which still exhibits high stability in all hydration levels, different from observations on the free enzyme system [Zhu L, Yang W, Meng YY, Xiao X, Guo Y, Pu X, Li M (2012) J Phys Chem B 116(10):3292-3304]. In low hydration levels, most water molecules mainly distribute near the protein surface and an increase in the water content could not fully exclude the organic solvent from the protein surface. However, the acetonitrile solvent displays a stronger ability to strip off water molecules from the protein than the hexane. In a summary, the difference in the calculated properties between the two organic solvents is almost significant in low water content (<10 %) and become to be small with increasing water content. In addition, some structural properties at 10 ~ 20 % v/v hydration zone, to large extent, approach to those in aqueous solution.

Calcium-ion-induced stabilization of the protease from Bacillus cereus WQ9-2 in aqueous hydrophilic solvents: effect of calcium ion binding on the hydration shell and intramolecular interactions

The neutral protease WQ from Bacillus cereus is stable in various aqueous organic mixtures, with the exception of those containing acetonitrile (ACN) and dimethylformamide (DMF). The stability of the enzyme in aqueous hydrophilic solvents was dramatically enhanced with the addition of calcium ions, with the degree of improvement in the half-life relative to different solutions ranging from fourfold to more than 70-fold. Studies of the kinetic constants showed that calcium ions induced slight conformational changes in the active site of the enzyme in aqueous ACN. We investigated the molecular mechanisms underlying this stabilizing effect by employing a combination of biophysical techniques and molecular dynamics simulation. In aqueous ACN, the intrinsic fluorescence and circular dichroism analysis demonstrated that the addition of calcium ions induced a relatively compact conformation and maintained both the native-like microenvironment near the tryptophan residues and the secondary structure. Alternatively, homology modeling confirmed the location of four calcium-ion-binding sites in the enzyme, and molecular dynamics simulation revealed that three other calcium ions were bound to the surface of the enzyme. Calcium ions, known as a type of kosmotrope, can strongly bond with water molecules, thus aiding in the formation of the regional hydration shell required for the maintenance of enzyme activity. In addition, the introduction of calcium ions resulted in the formation of additional ionic interactions, providing propitious means for protein stabilization. Thus, the stronger intramolecular interactions were also expected to contribute partially to the enhanced stability of the enzyme in an aqueous organic solvent.

[Influence of dioxane on the hydration shell of polypeptides]

Our interest in modifying influence of the low-molecular organic compounds on the hydration shell of biopolymers is due to a well-known fact that the former can act as regulators of the enzymatic activity changed the hydration shell of the latter. Dioxane was chosen because of its wide application in non-aqueous biocatalysis. In the present study we investigate the mechanisms of the dioxane influence on water of the first hydration layer of the model polypeptides by FTIR-spectroscopy during simultaneous sorption of water-dioxane vapors into the polypeptide films at low water activity. It was found that modification of the hydration shell of the studied polypeptides was mainly due to the indirect ordering ofpolypeptide secondary structure at penetration of dioxane molecules at these conditions.