Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 25 October 2024
Sec. Green and Sustainable Chemistry
This article is part of the Research Topic Green and Sustainable Chemistry Editor’s Pick 2024 View all 10 articles

Harnessing the potential of deep eutectic solvents in biocatalysis: design strategies using CO2 to formate reduction as a case study

Marijan Logaru&#x;i&#x;Marijan Logarušić1Karla &#x;ubarKarla Šubar1Maja Nikoli&#x;Maja Nikolić1Ana Jurinjak Tu&#x;ekAna Jurinjak Tušek1Anja Damjanovi&#x;Anja Damjanović1Mia Radovi&#x;Mia Radović1Ivana Radoj
i&#x; Redovnikovi&#x;Ivana Radojčić Redovniković1Polona 
nidar&#x;i
-PlazlPolona Žnidaršič-Plazl2Wolfgang KroutilWolfgang Kroutil3Marina Cvjetko Bubalo
Marina Cvjetko Bubalo1*
  • 1Faculty of Food Technology and Biotechnology, University of Zagreb, Zagreb, Croatia
  • 2Faculty of Chemistry and Chemical Technology, University of Ljubljana, Ljubljana, Slovenia
  • 3Institute of Chemistry, University of Graz, Field of Excellence BioHealth, BioTechMed Graz, Graz, Austria

Introduction: Deep eutectic solvents (DESs) have emerged as green solvents with versatile applications, demonstrating significant potential in biocatalysis. They often increase the solubility of poorly water-soluble substrates, serve as smart co-substrates, modulate enzyme stereoselectivity, and potentially improve enzyme activity and stability. Despite these advantages, screening for an optimal DES and determining the appropriate water content for a given biocatalytic reaction remains a complex and time-consuming process, posing a significant challenge.

Methods: This paper discusses the rational design of DES tailored to a given biocatalytic system through a combination of experimental screening and computational tools, guided by performance targets defined by solvent properties and process constraints. The efficacy of this approach is demonstrated by the reduction of CO2 to formate catalyzed by NADH-dependent formate dehydrogenase (FDH). By systematically analyzing FDH activity and stability, NADH stability (both long-term and short-term stability after solvent saturation with CO2), and CO2 solubility in initially selected glycerol-based DESs, we were able to skillfully guide the DES screening process.

Results and discussion: Considering trade-offs between experimentally determined performance metrics of DESs, 20% solution of choline chloride:glycerol in phosphate buffer (ChCl:Gly80%B) was identified as the most promising solvent system for a given reaction. Using ChCl:Gly as a co-solvent resulted in an almost 15-fold increase in FDH half-life compared to the reference buffer and stabilized the coenzyme after the addition of CO2. Moreover, the 20% addition of ChCl:Gly to the buffer improved the volumetric productivity of FDH-catalyzed CO2 reduction in a batch system compared to the reference buffer. The exceptional stability of the enzyme in this co-solvent system shows great potential for application in continuous operation, which can significantly improve process productivity. Additionally, based on easily measurable physicochemical solvent properties and molecular descriptors derived from COSMO-RS, QSAR models were developed, which successfully predicted enzyme activity and stability, as well as coenzyme stability in selected solvent systems with DESs.

1 Introduction

Medium engineering is one of the main components of biocatalysis engineering, which also includes substrate engineering, protein (enzyme) engineering, biocatalyst (formulation) engineering, biocatalytic cascade engineering, and reactor engineering (Sheldon and Pereira, 2017). Its purpose is to enhance the performance of biocatalysts, improve substrate solubility, and positively influence reaction equilibrium, so it has been a significant research focus for over the past 30 years (Castillo et al., 2016; Sheldon and Pereira, 2017; Sheldon et al., 2023). Accordingly, numerous enzymatic syntheses have been developed in both aqueous and non-aqueous media, including organic solvents, supercritical fluids, ionic liquids (ILs) and, more recently, deep eutectic solvents (DESs). The latter two, ILs and DESs, are of particular interest due to their high tunability as they can fulfil specific process requirements and meet “green solvent” criteria (Lozano et al., 2010; Villa et al., 2019; Panić et al., 2021; van Schie et al., 2021; Žnidaršič-Plazl, 2021a; Sheldon et al., 2023; Zhang et al., 2024). Eutectic systems have long been known as multi-component mixtures used in materials chemistry and engineering. However, in 2003 Abbott et al. coined the term “deep eutectic solvents” to describe mixtures that solidify at temperatures well below (deeper) than the crystallization points of their individual components, such that these compositions remain liquid even at room temperature (Abbott et al., 2003). In general, these solvents are formed by mixing two or more compounds that are normally solid at room temperature; when combined in certain molar ratios, they form a liquid solution due to molecular interactions (mainly hydrogen bonds) that lower the melting point of the mixture. In practice, DES can be easily prepared by combining hydrogen bond acceptors (HBAs) such as choline chloride and betaine with hydrogen bond donors (HBDs) such as polyols (e.g., glycerol, ethylene glycol, sorbitol), organic acids (e.g., citric, malic, oxalic acid), amino acids (e.g., alanine, proline), sugars (e.g., glucose, sucrose, trehalose) or amide urea with an atom economy of 100%. These combinations lead to different DESs, which almost always contain water to reduce the viscosity or to adjust certain properties of the DES (Hansen et al., 2021). The diversity of DES-forming components has led to many new and structurally different DESs, refining the initial definition by Abbott’s group. Recently, a (deep) eutectic solvent/system has been defined as a liquid system with eutectic properties that remains liquid at a given temperature, even if one component would normally be solid (Abranches and Coutinho, 2022). The wide range of their possible structural combinations, coupled with their sustainability and their distinctive physicochemical properties (non-volatility, non-flammability, easy and clean preparation, low to moderate toxicity), as well as the ability to fine-tune their solvent properties, make them ideal candidates for the development of efficient and sustainable processes or products (Cvjetko Bubalo et al., 2015).

In 2008, Gorke et al. published a pioneering work in which they demonstrated the potential of DESs as solvents or co-solvents for biocatalytic reactions (Gorke et al., 2008). It was shown that the media composed of urea, a strong HBD, with choline chloride, a HBA, does not denature hydrolases and even increases enzyme activity. Their study revealed that the components of DESs are significantly less denaturizing agents than expected, suggesting that the hydrogen bonding network in DESs lowers the chemical potential of the components. Since then, interest in the use of these solvents for various biocatalytic reactions has surged, as evidenced by the exponential growth in related publications over the past decade (according to Web of Science, more than 250 scientific papers have been published on this topic since 2009). DESs have shown considerable potential in enhancing biocatalysis through various mechanisms. These solvents can improve substrate solubility, act as smart co-substrates, and influence enzyme stereoselectivity. Furthermore, they are noted for their ability to improve enzyme activity and especially operational stability in reaction media (Zhang et al., 2024). In a recent study, a hydrophilic DES was used to tailor the properties of a copolymeric hydrogel utilized for enzyme immobilization. This addition not only exhibited superior mechanical properties of DES-infused hydrogel but also increased permeability to the specific substrate in a given biocatalytic reaction (Menegatti et al., 2024). Stabilization of enzymes, whether in homogeneous reactions or through enzyme immobilization, remains one of the major challenges in the development of long-term continuous biocatalytic processes. Overcoming this challenge is crucial as it can greatly increase the total turnover number and thus efficiency of biocatalysts, leading to the intensification of biocatalytic processes (Žnidaršič-Plazl, 2021b).

Despite the numerous advantages of DESs, the modification of the conventional aqueous medium with these solvents can significantly affect various aspects of a reaction, sometimes in undesirable ways. Therefore, finding the optimal DES for a given reaction remains a challenge. A promising method for rationally designing solvents for a specific biocatalytic system involves a complementary approach that integrates experimental screening with computational tools. If such an approach is successful, it would be possible to realize the full potential of these solvents while avoiding some problems associated with their use, such as high viscosity, product recovery issues and still quite unexplored methods of DES recovery and recycling. In this paper, we briefly review the current literature on DES-assisted biocatalysis and emphasize the need for strategic design of DESs for a given reaction. Furthermore, we propose and discuss complementary experimental and in silico methods to tackle the complexity of DES screening. Finally, this approach is illustrated by a case study of CO2 reduction catalyzed by formate dehydrogenase.

2 Optimizing the synergy between DES and biocatalysis: the quest for the ideal solvent

The synergy between DESs and biocatalysis is ideally suited for the efficient and sustainable production of commercially important products (Sheldon et al., 2023; Zhang et al., 2024). Biocatalysis, which is known to enable complex transformations with high regio-, chemo-, and enantioselectivity under mild and cost-effective conditions, could be greatly enhanced by DESs. These tunable solvents with their wide range of structural possibilities could be a particularly valuable aid in processes facing challenges such as enzyme stability/activity, cofactor stability, substrate solubility, product inhibition, and negative environmental impact. To date, a variety of DES-assisted biocatalytic reactions have been developed using versatile enzymes such as hydrolases (e.g., lipase, epoxide hydrolases, dehalogenases), oxidoreductases (e.g., alcohol dehydrogenases, laccases, peroxidases, monooxygenases), lyases (e.g., benzaldehyde lyase, phenolic acid decarboxylase), and transferases (e.g., amine transaminases) (Panić et al., 2021; Domingues et al., 2024; Zhang et al., 2024). In these reactions, DESs primarily act as (co-) solvents or additives. In some cases, they fulfil a dual function by acting as both solvents and (co-) substrates (Mourelle-Insua et al., 2019; Pätzold et al., 2019b). It should be noted that when DESs act as solvents, they are rarely used in their pure form. Instead, they are usually mixed with water, which is tightly incorporated into the DES cavities up to a proportion of about 50%. For example, Hammond et al. studied the effect of water on DES choline chloride:urea nanostructure, finding that it remains stable up to about 42% water (w/w) due to solvophobic sequestration of water into nanostructured domains around the cholinium cation. At 51% water (w/w), the structure disrupts, and water–water and DES–water interactions prevail, making the mixture resemble an aqueous solution of DES components (Hammond et al., 2017), which was later confirmed by Sapir and Harries (Sapir and Harries, 2020). Moreover, Nolasco et al. observed that in the same DES system containing <30% water (w/w), water molecules promote a strengthening of hydrogen bonds between choline chloride and urea (Nolasco et al., 2022). In biocatalysis, DES-water mixtures, whether in a water-in-DES or DES-in-water state, have been shown to be beneficial. These mixtures not only meet the enzyme’s hydration needs but also reduce the inherent viscosity of DESs to a more practical level (Sanchez-Fernandez et al., 2022).

One of the most prominent advantages of using DESs in biocatalysis is their versatility and tunability. By varying the component, an HBD and an HBA, and fine-tuning the water content, DESs can be tailored to the specific needs of biocatalytic processes. It used to be estimated that there are about 106 possible structural variations of DESs (Panić et al., 2021) but with the continuous reporting of new DES components, multicomponent DESs, and various possible molar ratios of components (including water), the number of possible structural combinations appears to be unlimited. This vast chemical space offers enormous potential for solvent design, but also poses a major challenge: How to identify the “ideal” DES that fulfils the various criteria of a given biocatalytic system? Obviously, this task is time-consuming and sometimes leads to contradictory results in DES performance with respect to specific targets that are not easy to reconcile. For example, when Wu et al. tested different cholinium chloride- and cholinium acetate-based DESs with four HBDs (urea, glycerol, acetamide, ethylene glycol) in three molar ratios and at different water contents for the activity and stability of horseradish peroxidase, it was found that increasing the DES concentration as a co-solvent in aqueous media made the enzyme much more stable but less active than in the reference buffer (Wu et al., 2014). This pattern has been confirmed in numerous cases: DESs that stabilize enzymes (usually co-solvents with low water content) are often poor media for catalytic reactions due to the reduced enzyme activity and high solvent viscosities, resulting in slow reaction rates (Taklimi et al., 2023). Additionally, our extensive experience with these solvents has taught us that hydrophobic substrates are readily soluble in acidic DESs containing organic acids as HBDs or in hydrophobic DESs based on terpenes and fatty acids, which are generally less suitable media for enzymes than hydrophilic DESs based on polyols, sugar alcohols and sugars (with the exception of lipases, which work well in hydrophobic DES (Elgharbawy et al., 2023; Ma et al., 2024). All of this implies that it is crucial to strike a perfect balance, which essentially means finding the DES candidate that satisfies all relevant factors for the application, such as stabilization and activation of enzymes or improved substrate solubility, in a balanced manner.

Traditionally, navigating this vast design space to determine an optimal DES for specific applications has mostly relied on trial-and-error methods (e.g., measuring rection rates and calculating reaction yields in different DESs), without systematic exploration of the influence of DES on specific reaction targets (e.g., enzyme behavior and substrate solubility in these solvents), and with structure-property relationships using available computational tools. This has hindered the strategic design of these inherently tailored solvents. In addition, the understanding of intermolecular interplay between DES components, water and reactants is still in its infancy although our comprehension of these solvents has evolved over the last decade.

To fully exploit the potential of DESs and accelerate the design of a DES with optimal properties, several important steps are required: (i) Obtaining a thorough understanding of the DES structure at the molecular level and its corresponding properties; (ii) Understanding the interactions of DES with reacting compounds (biocatalyst, substrate(s)/product(s), cofactor(s), co-substrate(s)) at the molecular level by combining experimental and computational approaches; (iii) Comprehensive experimental screening of promising DES candidates with respect to factors relevant to the process of interest (solubility and stability of reacting compounds, enzyme activity, and product/substrate inhibition); (iv) Utilizing computational tools to establish relationships between DES composition and desired properties, and possibly develop predictive Quantitative Structure-Property Relationships (QSPRs); (v) Development of sustainable methods for product recovery, and recycle/reuse of DES (Figure 1). In addition to the above steps, green chemistry principles should be considered when selecting DES, including consideration of the environmental footprint of the solvent, cost and ease of recycling (Abildskov et al., 2013). By favoring sustainable and non-toxic DESs, we could minimize environmental impact while optimizing the biocatalytic processes.

Figure 1
www.frontiersin.org

Figure 1. Rational design of DES for sustainable biocatalysis.

Molecular dynamics simulations have proven to be powerful in understanding thermodynamic and transport processes in DESs at the atomic level, providing insights into fundamental phenomena that may not be accessible through experiments (Bittner et al., 2024). Another computational tool, the Conductor-like Screening Model for Real Solvents (COSMO-RS), provides a computational approach to generate σ-profiles (molecular descriptors) of DESs, which provide essential information on hydrogen bonding, and electrostatic and dispersion interactions in solutions (Klamt, 2005). These descriptors enable the quantification of structural changes and are widely used for high-throughput screening of DES candidates regarding the solubility of organic and inorganic molecules in DESs, but are also very useful in development of QSPR models and machine learning to predict the physicochemical properties of DESs (Benguerba et al., 2019; Lemaoui et al., 2020; Lemaoui et al., 2022). Recently, we developed a robust and reliable QSPR model based on artificial neural networks to predict the ability of DESs to stabilize dehydrogenase (Radović et al., 2024). Thus, we have shown for the first time that such models with their high prediction accuracy provide a reliable means for in silico screening of DESs, obviating the need for labor-intensive experimental screening and paving the way for rational design of tailored solvents. In general, mathematical and simulation models allow for the exploration of a vast chemical space of DESs that would be impossible to capture experimentally, making them essential for optimizing these solvents for industrial applications. However, accurately modeling and simulating DESs at the molecular level is challenging due to the simultaneous occurrence of multiple interactions (Lemaoui et al., 2020). Concerning the ecological footprint of DESs, they generally exhibit low to moderate toxicity to vertebrates, invertebrates, and various animal cell lines (Lomba et al., 2021). Several authors reported that DESs produced from natural metabolites can be classified as “biodegradable” according to OECD guidelines (Radošević et al., 2015; Torregrosa-Crespo et al., 2020). This means that as long as natural metabolites such as choline, betaine, polyols, sugars, and amino acids are used in the preparation of DESs, these solvents are expected to have a low environmental impact. In terms of cost, DESs are generally easy and inexpensive to prepare from readily available and renewable materials with 100% atom economy. Depending on their primary constituents, the cost of DESs ranges from €7 to €100 kg−1, which is comparable to organic solvents (Panić et al., 2021). The cost flexibility allows solvent selection that matches the price of the product, making DESs an attractive option for large-scale applications (Rente et al., 2022).

It should be noted that the downstream processes of the biocatalytic reactions in DESs, which are often the most complex and costly part of the whole process, are not yet sufficiently addressed. One of the major challenges in DES-assisted biocatalytic reactions is the isolation of the product as well as the regeneration and reuse of the DES. In this context, an advantageous property of DESs — their non-volatility — becomes a problem unless the products are volatile and can be removed by evaporation, as has already been shown for the recovery of butyl butyrate from ILs (Pohar et al., 2012). Due to their low vapor pressure, it is practically impossible to remove DESs by evaporation. Therefore, various techniques have been proposed for the recovery of target compounds and recycling of DESs, including liquid-liquid extraction with a different solvent, solid-liquid extraction with macroporous resins, and the use of antisolvents (Panić et al., 2021; Zhang et al., 2024).

In this field, two studies stand out as they offer new perspectives for downstream processes with these solvents. First, Maugeri et al. showed the separation of alcohol and ester in DES after kinetic resolution, with the ester forming a separate phase, a viable method when the product or substrate is insoluble in DES (Maugeri et al., 2012). Secondly, in the study by Pätzold et al., DES compounds (menthol and dodecanoic acid) acted simultaneously as substrates and reaction solvent in lipase-catalyzed esterification for the synthesis of (-)-methyl dodecanoate, where the product was separated from the DES reaction mixture by a vacuum distillation step, and a second esterification reaction could be performed with the recovered (-)-menthol (Pätzold et al., 2019a). Both studies show that through innovative thinking in downstream processing, the unique properties of these solvents can be utilized, making them very attractive for sustainable biocatalysis.

3 Navigating DES screening complexity: enzymatic CO2 conversion to formate as a case study

Given the demonstrated potential of DESs as versatile media for biocatalysis, our study aimed to illustrate effective strategies for screening these solvents using the case of formate dehydrogenase (FDH)-catalyzed reduction of CO2.

The reduction of CO2 — whether by electrochemical or (bio) catalytic means — produces formate, the first stable intermediate product in the conversion of CO2. The resulting formate can be further converted to valuable chemicals such as formaldehyde and methanol through additional enzymatic reactions involving aldehyde dehydrogenase and alcohol dehydrogenase. However, formate is now increasingly recognized as an energy source in its own right (Reda et al., 2008; Villa et al., 2023). It is known that the enzyme FDH is able to reduce CO2 to formic acid using a coenzyme such as NADH or NADPH (Figure 2) (Villa et al., 2023). The complexity of this reaction arises from several limitations: (i) the low concentration of CO2 available for the enzyme, leading to a low reaction rate, (ii) the relative instability of FDH and the NADH coenzyme in aqueous environments, (iii) the acidification of the reaction medium, firstly by the dissolution of CO2 and secondly by formic acid formation, can cause further FDH inhibition and NADH degradation (Zhang et al., 2018; Calzadiaz-Ramirez and Meyer, 2022).

Figure 2
www.frontiersin.org

Figure 2. Formate dehydrogenase (FDH)-catalysed CO2 conversion to formate.

Removing these limitations is critical to making FDH-catalyzed reduction of CO2 to formate a practical and economically viable process for industrial and environmental applications. Among other approaches, switching to unconventional media could help solve the problems stated above. Zhang et al. were the first to show that ILs can act as cosolvents for enhanced conversion of CO2 to methanol catalyzed by NADH-dependent FDH: the 67.1% conversion achieved in 1-butyl-3-methylimidazolium tetrafluoroborate was more than twice that in phosphate buffer (24.3%) (Zhang et al., 2018). Later, the same group demonstrated the advantages of introducing DES, L-serine:glycerol at a molar ratio of 1:6, into the electro-enzymatic conversion of CO2 and achieved a 16-fold higher reaction yield compared to the control reaction in the buffer (Zhang et al., 2022).

Instead of the traditional approach of performing the reaction and then testing for suitability, which would be extremely complex given the vast number of possible DES candidates, here we set a performance targets defined by solvent properties and process constraints. By identifying DES that fulfil these predefined goals in a balanced way, we can significantly reduce the number of trials required and allow for a more strategic allocation of resources (Abildskov et al., 2013). Based on the above considerations, the design of the solvent system with DES was divided into several phases. First, we narrowed down the list of DES candidates based on the data available in the literature on dehydrogenase performance in these solvents as well as their sustainability attributes. Then, to gain a better understanding of DES potential in overcoming the limitations associated with the reaction, the influence of selected DESs on FDH activity and stability, NADH stability, and CO2 solubility was systematically investigated. A computational analysis was further performed to gain insight into the DES structure and the corresponding physicochemical properties affecting the measured parameters. For the most promising DES candidates, further experiments were performed to investigate NADH stability under acidification at CO2 saturation. Finally, considering the trade-offs between the results in different solvent systems with DESs in terms of performance objectives, the most suitable DES candidate was proposed. This solvent system was then tested for the reaction and compared with the results in the reference buffer. At this point, it should be mentioned that for the purpose of this study, we excluded in situ cofactor regeneration, e.g., based on chemical, electrochemical, photochemical, or enzymatic processes (Wichmann and Vasic-Racki, 2005) to maintain methodological clarity and to avoid potential confounding effects from new reactants.

3.1 DESs selection and characterization

Ensuring the stability of enzymes is crucial for their industrial implementation, as it increases operational stability, prolong activity and thus increases total turnover number, and improve cost-efficiency of the process (Žnidaršič-Plazl, 2021b; Woodley, 2022). To systematically evaluate DES candidates for performance goals, we first narrowed the list to those documented in the literature as beneficial for dehydrogenases, focusing specifically on those known to stabilize FDH. Studies suggest that polyol-based DESs, containing either choline chloride or betaine as HBA, are the most effective stabilizing media for various dehydrogenases, including FDH (Bittner et al., 2022; Gajardo-Parra et al., 2023). First, seven polyol-based DESs with either choline chloride or betaine as HBA were tested at three water contents (up to 50%, w/w) for their ability to stabilize FDH. It was confirmed that glycerol-based DES containing both tested HBAs are optimal candidates for stabilizing the enzyme upon prolonged incubation at room temperature (data not shown).

Glycerol-based DESs were recently demonstrated to stabilize NAD coenzymes (Radović et al., 2022). Besides, Leron and Li (2013), Leron et al. (2013) and Biswas et al. (2023) reported that among several choline chloride-based DES containing urea, ethylene glycol, and glycerol as HBDs, the one with glycerol had the highest CO2 solubility. The above studies also emphasize that water plays a crucial role in enzyme performance and CO2 solubility. Finally, concerning the toxicological footprint, glycerol-based DES are considered non-toxic and biodegradable (Radošević et al., 2015).

Based on the above considerations, we selected two glycerol (Gly)-based DESs, with choline chloride (ChCl) or betaine (B) as HBA in a molar ratio of 1:2, and prepared the corresponding solutions in water (10%–90%, w/w). In parallel, we also prepared solutions in 50 mM potassium phosphate buffer (pH 7.5) to keep the pH close to the enzyme’s optimal value and to possibly prevent a pH drop when CO2 is added to the reaction medium. As previously mentioned, DESs diluted with more than 50% water (w/w) can be considered aqueous solutions of DES components (Hammond et al., 2017). Nevertheless, these mixtures were included in the study, as a high water content within DES is often essential for enzymes to sustain their catalytic activity (Taklimi et al., 2023).

A total of 20 DES-based solvents were prepared and characterized for their physicochemical properties (pH, density, and viscosity) relevant to the reaction (Table 1). As expected, the densities and viscosities of the DES aqueous solutions were strongly influenced by the water/buffer content, peaking at mixtures with 10% water (up to 1.21 g cm-³ and 353.70 mPa s for B:Gly10%W, and 1.17 g cm-³ and 82.63 mPa s for ChCl:Gly10%B). In general, B:Gly-based mixtures were denser and more viscous than their ChCl:Gly-based counterparts. All mixtures tested had pH values ranging from 5.3 to 9.2, with betaine-based mixtures being more acidic than ChCl-based ones. Dissolving DES in buffer generally maintained the solutions at pH values between 7.5 and 8. All DES mixtures remained stable for 3 months under laboratory conditions and showed no signs of contamination or precipitation.

Table 1
www.frontiersin.org

Table 1. List of DES solutions in water/buffer (10%–90%, w/w) and buffer (50 mM potassium phosphate buffer pH 7.5) used for experimental screening, together with physicochemical properties and corresponding concentration of dissolved CO2 (cs), pH of CO2 saturated solutions (pH*), first-order degradation rate constant of NADH in CO2 saturated solutions (kNADH*), and residual FDH activity (ARes) (cFDH = 12.8 mg mL-1, t = 14 days, T = 30°C).

For statistical analysis and mathematical modelling (Section 3.4.), the identification of a molecular representation that converts the component structures into descriptive features for numerical evaluation is essential (Venkatraman et al., 2018). An advanced and accurate molecular representation is the σ-profile (sigma profile), an unnormalized histogram of the screened surface charge of a molecule (Klamt, 2005). σ-profiles are distinguished from other representations by the fact that they capture nuanced effects such as polarizability and electron density asymmetry (Abranches et al., 2022). The σ-profile can be divided into three key regions: (i) the HBD region with negative charge densities, (ii) the non-polar region with nearly neutral charge densities, and (iii) the HBA region with positive charge densities (Figure 3A). This division is based on the fact that each atom in an HBA or HBD molecule is identifiable by a distinct peak with a specific screening charge density (σ) value (Lemaoui et al., 2020).

Figure 3
www.frontiersin.org

Figure 3. σ-surface (A) and σ-profiles (B) of individual DES constituents (choline chloride, glycerol and water), and σ-profile of corresponding DES solutions in buffer (10%–90%, w/w).

Here, the σ-profile of each DES-based mixture was calculated using BOVIA COSMOtherm software: the σ-profile curves for each HBA and HBD were divided into 10 regions, and the area under each region was calculated considering the molar ratios of the components and the water content (Supplementary Table S1). For glycerol, the σ-profile reveals peaks at negative polar coordinates (left side) corresponding to the positively polar H atoms in the -OH group, while peaks at positive polar coordinates (right side) correspond to the O atoms in the -OH group (Figure 3A) (Cheng et al., 2018). These polar regions interact with opposite polar segments in a solution. The σ-profile’s extension into strongly polar regions (−0.022 e/Å2 < σ < −0.01 e/Å2 on the left and 0.01 e/Å2 < σ < 0.013 e/Å2 on the right) is asymmetric, indicating that glycerol has a larger positive polarity surface area. This electrostatic misfit suggests that glycerol tends to act as a HBD and thus shows greater affinity for HBA in a solution (Cui et al., 2021). In contrast, the σ-profile of water extends more symmetrically into strongly polar regions, indicating its balanced ability to act as both a HBD and HBA. Finally, for ChCl, as well as for betaine (Radović et al., 2024), the strongest peak appears in the non-polar region, corresponding to the cholinium cation, followed by a peak in the HBA region, associated with the Cl⁻ anion. Figure 3A illustrates the sigma surfaces of the DES components (choline chloride, glycerol and water) generated by TmoleX19. The colors represent a calculated charge gradient, ranging from charge-deficient to charge-dense regions: HBD regions are labelled as deep blue and HBA regions as deep red on the surface. Non-polar regions are marked in green (Quaid and Reza, 2023). Figure 3B shows that even small changes in the ChCl:Gly mixture, such as increasing the water content from 10% to 90%, result in solvents with different polarity distributions. This demonstrates that the software is capable of capturing nuanced phenomena, which is crucial for exploring the chemical landscape of these solvents and understanding their potential impact on enzyme and coenzyme behavior, as well as the solubility of reaction participants. We have recently demonstrated the same ability of the software for the system betaine:ethylene glycol with 3 water proportions (10, 30, 50% water, w/w) (Radović et al., 2024).

3.2 Activity and long-term stability of FDH in DESs

The activity of FDH in the reference buffer and in aqueous solutions of DESs at different concentrations were measured by monitoring the oxidation of sodium formate. For the stability test, enzyme solutions were incubated in selected solvent systems at 30°C for 14 days and residual enzyme activity (ARes) was measured at regular intervals (Table 1). Both in the reference buffer and in the mixtures with DESs, FDH inactivation followed first-order kinetics, allowing us to use this kinetic model to calculate the FDH half-life (t1/2, FDH) (Figure 4).

Figure 4
www.frontiersin.org

Figure 4. FDH relative activity (cFDH = 1.5 mg mL-1; cNADH = 0.1 mg mL-1; cformate = 10 mg mL-1, T = 25°C) and half-life (cFDH = 12.8 mg mL-1, T = 30°C), along with NADH half-life (cNADH = 0.03 mg mL-1, T = 25°C) in B:Gly and ChCl:Gly solutions in water/buffer (10%–90%, w/w) and the reference 50 mM potassium phosphate buffer (pH 7.5). FDH relative activities in DES mixtures are expressed as a percentage of the activity measured in the reference buffer.

As anticipated, FDH showed little or no activity in the DES mixtures containing ≤50% (w/w) water or buffer, while the activity increased with the addition of water and peaked in the highly diluted mixtures (90% water, w/w), although the values were still lower than those observed in the buffer (residual activities, AR, between 44% and 92%). Interestingly, it appears that dilutions with water generally resulted in better enzyme activity than their counterparts diluted with buffer, although buffered systems are closer to the enzyme’s optimal pH of 7.5 (determined experimentally, data not shown). For example, ChCl:Gly with 90% water content (ChCl:Gly90%W) with a pH of only 6.3 had an AR of about 92%, while its counterpart with the buffer (ChCl:Gly90%B) with a pH of 7.5, yielded the AR of only 44.5%.

DESs with a water/buffer content in the range of 30% - 80% (w/w) showed the best ability to stabilize the enzyme, with half-lives of up to 64 days, which is much higher than in the reference buffer, where t1/2, FDH was 2.1 days (Figure 4). On the other hand, highly diluted DESs (90% water/buffer, w/w) demonstrated a stabilizing effect on the enzyme comparable to that of the reference buffer.

The stability of FDH in buffered DES solutions (80% buffer, w/w) was also significantly increased compared to the reference buffer, with t1/2, FDH of 65.6 days and 29.9 days for B:Gly80%B and ChCl:Gly80%B, respectively. Again, DESs solutions with a water/buffer content of only 10% led to a faster destabilization of the enzyme. These results are consistent with previous findings that most dehydrogenases require more than 10% water content in DES to maintain their structural integrity. DESs absorb water in their hydrogen bonding network, reducing the availability of free water molecules required for enzyme hydration. This reduction in water activity can lead to dehydration and irreversible denaturation of the enzyme (Mourelle-Insua et al., 2019; Bittner et al., 2022; Radović et al., 2024). Overall, solvent systems with 80% buffer (w/w) exhibited an optimal balance between FDH activity and stability, which was particularly evident for B:Gly80%B, where the AR was 51.7% (Figure 4).

3.3 CO2 solubility and NADH stability in DESs

The saturated dissolved CO₂ concentrations (cs) in various DES solutions with water/buffer were evaluated after introduction of CO₂ at a flow rate of 100 mL min−1 until saturation. The results are presented in Table 1. All DESs mixtures with ≤50% water or buffer were poor media for dissolving CO2, with cs values between 282 and 913 mg L−1, which is lower than those observed in the reference buffer (1,029 mg L−1). Highly diluted DESs led to similar CO2 solubilities as in the buffer, with the highest improvements observed in ChCl:Gly90%B and B:Gly90%W, with values of 1,149 and 1,112 mg L−1, respectively.

To investigate the long-term stability of the NADH coenzyme, the changes in the UV-Vis absorption spectra of NADH were observed during a 14-day incubation at 25°C in the solvent systems described above. During incubation, the absorbance loss at 340 nm followed the first-order kinetics used to calculate the NADH half-lives shown in Figure 4. The results clearly show that DES composition plays an important role in the coenzyme degradation rate. Virtually all DES aqueous solutions, except ChCl:Gly90%B (t1/2, NADH = 2.9 days) and B:Gly90%W (t1/2, NADH = 0.8 days), stabilized the coenzyme compared to the reference buffer (t1/2, NADH = 4.6 days). In general, ChCl-based DESs were more suitable for coenzyme stabilization than betaine-based DESs. This was particularly pronounced for DESs with 10% buffer content (e.g., for ChCl:Gly10%B, t1/2, NADH was 44.2 days, while for B:Gly10%B, t1/2, NADH was 17.5 days). As evident from Figure 4, a higher DES content had a positive effect on the ability of solvent systems to stabilize the coenzyme for all DESs tested. This was most evidenced for B:Gly solutions in water, where the t1/2, NADH for B:Gly30%W and B:Gly90%W was 65.7 and 0.8 days, respectively.

To further navigate the DES screening, the short-term NADH stability in DES aqueous solutions saturated with CO2 was investigated. The presence of dissolved CO2 not only affects the intrinsic properties of the DESs but also the behavior of the system after substrate addition. It has been previously reported that the acidification of the reaction medium by dissolving CO2 leads to enhanced NADH degradation, which directly affects the conversion of CO2 to formic acid (Zhang et al., 2018). Therefore, NADH solutions in the solvent systems described above were monitored over a period of 90 min, and the corresponding degradation constants (kNADH*) were calculated using a first-order kinetic model (Table 1). In general, buffer-diluted DESs were equally or more successful in stabilization of the coenzyme over the tested period (kNADH* ≤ 0.003 min⁻1) than the reference buffer (kNADH* = 0.003 min⁻1), while water-based solvent systems with DESs were poor media in this regard, especially at high water contents (kNADH* up to 0.009 min⁻1). This effect is directly related to the inability of water-diluted DESs to maintain pH close to neutral. For example, after the introduction of CO2 into B:Gly80%W, pH decreased to 3.75, resulting in the highest observed kNADH* value of 0.008 min⁻1. Moreover, buffered solutions of ChCl-based DESs maintained a higher pH than betaine-based DES solutions after CO2 saturation, resulting in complete stabilization of NADH over the time tested, except for ChCl:Gly50%B (kNADH* = 0.003 min⁻1).

3.4 Statistical analysis of data and development of mathematical models

The correlations between the physicochemical properties (pH, viscosity and density) of the DESs used, the DES descriptors (σ-profiles), FDH performance, long-term NADH stability and CO2 solubility were analyzed using the Spearman correlation matrix, which was selected due to the non-normal data distribution. The analysis confirmed our above assumptions and findings from the experimental screening: the targeted properties of the reaction affected by the DES composition led to contradictory results (Table 2). First, both FDH activity and stability in tested solvent systems showed negative correlations with NADH stability. In addition, FDH activity values demonstrated negative correlations with all analyzed variables except CO2 solubility, favouring aqueous/buffered media over DES solutions with lower amounts of aqueous phase. Besides, negative correlations were found also for FDH and NADH stability with CO2 solubility (Table 2).

Table 2
www.frontiersin.org

Table 2. Spearman correlation matrix.

The Spearman correlation matrix shown in Table 2 highlights the complex interplay between the physicochemical properties of DES (pH, density, viscosity) and the DES descriptors with the performance characteristics of FDH and the stability of its cofactor NADH. As expected, significant negative correlations were observed between FDH activity and DES density/viscosity, likely due to the lower mobility and diffusion of substrates and enzymes. Conversely, NADH stability showed significant positive correlation with the density/viscosity of the solvent systems with DESs, which can be attributed to the slower overall dynamics of the solvents providing a less detrimental environment for dissolved NADH.

FDH activity and NADH stability are significantly impacted by pH value, while for FDH stability there is no correlation with this chemical property. It is well-established that enzyme activity (Bisswanger, 2017) and NADH stability (Zachos et al., 2019) are pH-dependent. However, the consistent reports on pH-independent ability of DESs to stabilize various enzymes remains puzzling. This intriguing observation hints at other mechanisms, such as direct interactions between DES components and proteins or with nearby water molecules, potentially altering the medium’s water activity (Damjanović et al., 2024).

Furthermore, the analysis revealed that nearly all DES descriptors significantly impacted the targeted properties. Specifically, S1mix and S2mix (HBD region, medium polarity), S3mix - S5mix (nonpolar region, positive charges), S6mixS8mix (nonpolar region, negative charges) and S9mix (HBA region) (Lemaoui et al., 2020), all demonstrated a significant influence on the properties being studied. FDH activity exhibited a negative correlation with all the descriptors, except for S2mix, which showed a positive correlation. Conversely, both FDH stability and NADH stability demonstrated positive correlations with all the descriptors, except for S2mix, which displayed a negative correlation. Interestingly, the analysis suggests an inverse relationship between enzyme activity and stability in DESs: solvents rich in HBA and non-polar domains stabilize the enzyme (and coenzyme), while HBD-rich solvents enhance enzyme activity but may lead to destabilization. These findings emphasize the delicate balance between enzyme activity and stability in DESs, driven by their specific compositional properties (as described by σ-descriptors). Similar interplay between the enzyme’s active and stable (but inactive) states, influenced by the water content in DESs, has been confirmed in our recent study on lysozyme behaviour in DESs based on various naturally occurring osmolytes (Damjanović et al., 2024). Insights into this relationship could guide the rational design of DESs for optimized biocatalysis applications: by performing similar statistical analyses on a larger set of DESs, it may be possible to predict an ideal σ-profile shape, and thereby identify or design the most suitable DES for a specific purpose.

According to discussed above, our next step was to see if it was possible to develop a simple QSAR model to summarize the relationship between the targeted properties (FDH activity, FDH stability, and NADH stability) and the DES descriptors (Figure 3; Supplementary Table S1) along with the physicochemical properties (Table 1) of the DESs using piecewise linear regression (PLR) (Figure 5). The latter is a powerful tool for modelling complex relationships in a simple and interpretable way, especially when the relationship between variables changes at a certain point. The input variables of the PLR models were selected based on the significant correlations in the Spearman correlation matrix. The relationship between the observed data and model predictions was also estimated using the coefficient of determination for prediction (Rpred2), the adjusted coefficient of determination for calibration (Rpred 2adj), the root mean square error of prediction (RMSEP), the ratio of prediction to deviation (RPD) and the ratio of the error range (RER).

Figure 5
www.frontiersin.org

Figure 5. Comparison between experimental data and piecewise linear model predicted data for (A) FDH activity (ARes), (B) FDH half-life (t1/2, FDH) and (C) NADH half-life (t1/2, NADH). (○) calibration data set (Δ) prediction data set.

As shown in Table 3 and Figure 5, the developed models describe the experimental data with high precision. The best agreement between the experimental data and the data predicted by the model was obtained for the FDH activity (Figure 5A) (Rcal2 = 0.914, Rcal2adj = 0.913, RMSEC = 7.203%, Rpred2 = 0.905, Rpred2adj = 903, RMSEP = 7.379%, RPD = 5.589, RER = 18.324). On the other hand, the largest scatter between the model and experimental data was found for the FDH half-life (Figure 5B) (Rcal2 = 0.829, Rcal2adj = 0.827, RMSEC = 4.775 days, Rpred2 = 0.738, Rpred2adj = 0.731, RMSEP = 4.944%, RPD 3.292, RER = 10.456).

Table 3
www.frontiersin.org

Table 3. Pricewise linear regression models for prediction of FDH stability, FDH activity and NADH stability based on the specific input variables. (Coefficient of determination for calibration (Rcal2), the adjusted coefficient of determination for calibration (Rcal2adj), the root mean square error for calibration (RMSEC), the coefficient of determination for prediction (Rpred2), the adjusted coefficient of determination for calibration (Rpred2adj), the root mean square error of prediction (RMSEP), the ratio of prediction to deviation (RPD),the ratio of the error range (RER)).

According to Hussain et al., an R2 value of 0.75 is considered significant, an R2 value of 0.50 is considered moderate, and an R2 value of 0.26 is considered weak (Hussain et al., 2018). Furthermore, models with RPD < 1.4 are considered non-reliable, those with RPD in the range from 1.4 to 2 are considered fair, while models with RPD > 2 are described as excellent models (Chang et al., 2001). Models with RER > 4 are acceptable for data screening, models with RER > 10 can be used for quality control, while models with RER > 15 can be used for quantification (Sim et al., 2023). Therefore, the PLR models developed for the prediction of FDH activity and NADH half-life based on Rpred2 can be considered substantial, while the model developed for the prediction of FDH half-life can be considered moderate. Based on the RPD values, all three models developed can be considered reliable. And based on the RER values, the model developed for the prediction of FDH activity can be used for quantification (RER = 18.324), while the other two modes can be used for quality control. Therefore, it can be concluded that the feasibility of mathematical models for predicting targeted properties or applications of DES using easily measurable physicochemical properties and chemical descriptors, as demonstrated here, could be valuable for both industrial applications and research efforts focusing on these solvents. Additionally, utilizing these QSPR models may not only assist in predicting the properties of interest but also provide valuable insights into the relationship between the structure of the DES and its measurable properties. By analyzing how various structural features influence the targeted properties, these models can help unravel the underlying mechanisms driving behavior of biomolecules in DESs. This understanding can inform the design and optimization of DESs, leading to more effective and tailored applications in various fields.

3.5 Trade-off between the performance of DESs with respect to target properties

We have demonstrated that the evaluation of different targets related to the tested reaction often leads to contradictory results regarding the optimal DES. For example, the enzyme dissolved in B:Gly30%B showed remarkable stability with a half-life of 40.1 days, while its relative activity was less than 2% of that in the buffer. This DES also showed average performance in the long-term stability of the coenzyme, with a half-life of 12.1 days. In contrast, the enzyme dissolved in B:Gly90%B maintained a high relative activity, which was 71.2% of that in the reference buffer, but was one of the worst candidates for enzyme stabilization with a half-life of only 6.0 days. In general, “concentrated” DESs (<50% water/buffer, w/w), with η > 18 mPa s showed high efficacy in stabilizing both the enzyme and the coenzyme, but also high viscosity, which poses significant challenges for scaling up processes with these solvents. Based on these observations, it is crucial to reconcile these results by finding a DES that optimally fulfils the desired properties and thus contributes to the overall efficiency of the process. A graphical representation shown in Figures 6A, B illustrates the trade-off between the performance of ChCl:Gly and B:Gly, respectively, with different buffer proportions with respect to the targeted property values (FDH activity/stability, CO2 solubility, NADH stability in saturated CO2 solutions, and compared to the reference buffer). The radar chart is bounded by the respective lower and upper limits of each target property, with ratings, ranging from 0 to 100, reflecting the performance of DES-based solvents relative to the best candidate for each target property. At this point, DESs diluted with water were omitted due to their poor ability to stabilise NADH after acidification of the medium due to CO2 introduction (Table 1).

Figure 6
www.frontiersin.org

Figure 6. Radar plot evaluating ChCl:Gly (A) and B:Gly (B) and their corresponding solutions in buffer (10%–90%, w/w) in terms of target properties (FDH activity and stability, CO2 solubility in DES systems, stability of NADH after solvent saturation with CO2, and DES systems’ viscosity). The radar chart is bounded by the specific lower and upper limits for each target property. Ratings, ranging from 0 to 100, reflect the performance of DES-based solvents relative to the best candidate for each target property.

As can be seen from Figure 6A, ChCl:Gly80%B has a balanced distribution across all target property values, which is crucial for the simultaneous optimization of the reaction where all design objectives are equally important. In particular, this DES, which in this case could be considered an additive rather than a solvent (Hammond et al., 2017), had the highest values for all target properties except for FDH activity, where the reference buffer resulted in the highest values. It should be emphasized that ChCl:Gly80%B had an almost 15-fold higher t1/2, FDH value compared to the buffer and stabilized the coenzyme more effectively when CO2 was dissolved in the solvent system. B:Gly80%B also showed similar balanced behavior to ChCl:Gly80%B, but was ineffective in stabilizing NADH after dissolving CO2.

3.6 FDH-catalyzed CO2 conversion to formate in the most promising solvent system

Since ChCl:Gly dissolved with 80% (w/w) buffer (ChCl:Gly80%B) was found as the most promising solvent system for the FDH-catalyzed conversion of CO2 to formate, the reaction was performed in a medium pre-saturated with CO2 and product formation was monitored over time. The amount of formate produced by the enzymatic CO2 reduction in the selected solvent system and in the reference buffer is shown in Figure 7. The reaction performed in the DES-supplemented medium yielded 26.5 μmol mL−1 of formate with a volumetric productivity of 6.6 μmol mL−1·h−1, while the reaction in the buffer yielded 22.7 μmol mL−1 with a volumetric productivity of 5.6 μmol mL−1·h−1. These differences can be attributed to the slightly lower solubility of CO2 in the buffer and the pronounced degradation of NADH in the buffer due to acidification by dissolved CO2 (Wu et al., 1986).

Figure 7
www.frontiersin.org

Figure 7. FDH-catalysed CO2 conversion to formate in ChCl:Gly80%B (20%, w/w) solution of ChCl:Gly in potassium phosphate buffer) and the reference 50 mM potassium phosphate buffer (pH 7.5) (cFDH = 28.8 mg mL-1, cNADH = 16.4 mg mL-1, T = 25°C).

NADH stability was monitored in a separate experiment without enzyme addition: after 4 h (corresponding to the reaction time), the NADH concentration in the buffer decreased to 78% of its initial value, while in the DES medium the NADH concentration remained above 98% (data not shown). The shape of the concentration vs. time curve reflects the interplay of several factors, including enzyme activity, CO2 solubility and availability, and NADH degradation dynamics. The variations in formate concentration between the two solvents observed during the first 120 min may be attributed, on one hand, to the higher enzyme activity in the buffer compared to ChCl:Gly80%B, and, on the other hand, to the enhanced CO2 solubility and NADH stability provided by the DES. It is presumed that this dynamic interplay results in the reaction rate being sometimes higher in the buffer and at other times in the DES.

In addition to improving volumetric productivity, the true potential of using ChCl:Gly80%B lies in its ability to stabilize FDH and thus extend the enzyme’s half-life by up to 15-fold compared to the buffer. Although this DES property is not fully utilized in a reaction lasting only 4 h, it is of great advantage in a continuous process (e.g., using an enzymatic membrane reactor). Under steady-state conditions, the volumetric productivity of the DES-assisted process could be significantly increased, as the productivity inversely correlates with the enzyme’s deactivation rate constant. Therefore, the productivity of the selected process under study was estimated for continuous operation mode comprising the enzyme inactivation rate. Considering the following assumptions: (i) that the process runs under steady-state conditions, (ii) that the enzyme is continuously deactivated over time according to first-order kinetics, and (iii) that the substrate concentration remains relatively constant, the reaction performed in buffer over 10-day period would yield 185.7 μmol mL−1 formate, while the reaction carried out in ChCl:Gly80%B over the same period would yield 639.0 μmol mL−1 formate, which is approximately a 3.5-fold improvement. The calculation of overall productivity included estimation of the Michaelis–Menten kinetic parameters by fitting the NADH concentration profiles to the differential equation describing the change in substrate concentration over time. The results showed that both the maximum reaction rate (vmax) and the NADH saturation constant (Ks) were higher for the solvent system with DES. In the buffer, the constant values were vmax = 0.091 μmol mL−1 min−1 and Ks = 6.746 μmol mL−1, while for the systems with DES, the constants were vmax = 0.133 μmol mL−1 min−1 and Ks = 12.133 μmol mL−1.

Finally, downstream processing was not part of this study. Nevertheless, based on the available literature, liquid–liquid extraction with the green solvent 2-methyltetrahydrofuran (Slater et al., 2016; Laitinen et al., 2021) appears to be a promising option for the separation of dilute aqueous formate solutions in the context of developing sustainable formic acid production.

4 Conclusion

DESs offer considerable potential for improving biocatalytic processes. Finding the ideal solvent requires a balanced approach that considers all relevant factors, from the behavior of all reaction participants in the medium to downstream processing and the overall economic and environmental sustainability of the process. This paper presents a complementary strategy that integrates experimental screening with computational tools. By using performance targets defined by solvent properties and process constraints, this approach aims to facilitate the rational design of DESs tailored to specific biocatalytic systems. The effectiveness of the proposed approach is demonstrated using a case study of the NADH-dependent formate dehydrogenase-catalyzed reduction of CO2 to formate.

By systematically analyzing the performance of 20 DES-based solvents in terms of enzyme stability, activity, co-enzyme stability, and product solubility, we successfully navigated the DES screening process. It has been shown that certain DESs are highly efficient in stabilizing FDH and the coenzyme NADH, while none of the DES candidates were able to improve enzyme activity or the solubility of CO2. Additionally, we have demonstrated that the evaluation of different targets related to the tested reaction often leads to contradictory results regarding the optimal DES. By analyzing the data using Spearman correlation coefficients and evaluating trade-offs between the performance of DES-based solvents with respect to target properties, we identified a candidate, ChCl:Gly80%B, that exhibited a balanced distribution across all target property values. Moreover, this DES displayed an almost 15-fold higher FDH half-life value compared to the buffer and was more effective in stabilizing the coenzyme after addition of CO2. In the FDH-catalyzed reduction of CO2, ChCl:Gly80%B outperformed the conventionally used buffered media in terms of volumetric productivity in a batch process. This solvent candidate proved to be suitable for use in continuous process due to the exceptional stability of FDH in this solvent. Furthermore, we developed a simple QSAR model to summarize the relationship between the targeted properties and the DES characteristics. These models demonstrate the feasibility of mathematical models for predicting the specific performances of DESs based on easily measurable physicochemical properties and molecular descriptors of the solvent.

However, to fully realize the benefits of DESs in biocatalysis, further research and development of computational tools and sustainable methods for product recovery is crucial. Finally, our combined experimental and computational approach increases reliability, optimizes resources, and accelerates DES-based solvent development for enzyme stabilization, as well as for other industrial applications of these green solvents. This approach improves the scalability and cost-effectiveness of DESs and represents a significant advance in DES-based industrial processes.

5 Materials and methods

5.1 Materials

For the enzymatic reactions performed in this study, the NADH-dependent formate dehydrogenase (FDH) from Pseudomonas sp. 101 was used (see Supplementary Information) (Tishkov et al., 1993). Carbon dioxide (CO₂) with a purity of 99.5% was acquired from Messer Croatia Plin (Zaprešić, Croatia), while all other chemicals were purchased from Sigma-Aldrich (St. Louis, Missouri, United States). All materials had a purity of at least 99% and were used as received without further purification.

5.2 DES preparation and physicochemical characterization

For the preparation of DESs based on betaine (B) and choline chloride (ChCl), the hydrogen bond acceptor (HBA) and the hydrogen bond donor (HBD) were mixed in a molar ratio of 1:2 with water or 50 mM potassium phosphate buffer (pH 7.5) in defined proportions (Table 1). Prior to use, ChCl was dried in a vacuum concentrator at 60°C for 24 h. The mixtures were stirred and heated to 60°C until a colourless and homogeneous liquid was formed. All prepared DESs were stored in sealed bottles until further use. The pH values of the prepared DESs were measured using a pH glass electrode (Mettler Toledo, Greifensee, Switzerland). The properties of the prepared DESs (pH, density and viscosity) were determined at 25°C. Density was determined using the pycnometric method and the viscosity using a rotary viscometer (Anton Paar ViscoQC 300, Ashland, Virginia, United States). All measurements were performed in triplicates. The σ-descriptors of the DESs were calculated according to Panić et al. (2022) (Supplementary Table S1).

5.3 FDH activity and stability assays

To determine FDH activity, the FDH enzyme, NAD+ coenzyme, and formate substrate were added sequentially to various solvent systems to achieve final concentrations of 1.5 mg mL−1, 0.1 mg mL−1, and 10 mg mL−1, respectively. The total working volume of the assay was 250 µL. The NADH formation rate was measured immediately in a 96-well plate for 5 min at 340 nm using a UV-Vis spectrophotometer (SpectraMax® ABS Plus, Molecular Devices, San Jose, CA, United States).

To evaluate FDH stability, stock solutions of the enzyme (12.8 mg mL−1) were prepared in the tested DES solutions and the reference buffer. These solutions were stored in sealed vials at 30°C in the dark. Aliquots were taken at regular intervals over a 14-day period and analyzed for FDH activity using the method described above.

The first-order deactivation rate constant (kFDH, day-1) was evaluated from the first-order kinetic model for the decrease of residual enzymatic activity over time (ARes (t)) (Equation 1):

kFDH=1tlnARes,0ARest(1)

where ARes is residual enzyme activity (either at time zero, AR,0, or at time t, AR(t)). The kinetic parameters were estimated by fitting the experimental data to the nonlinear equation using the Levenberg–Marquardt algorithm implemented in WR Mathematica 10.0 (Wolfram Research, Champaign, United States).

The half-life of the enzyme (t1/2, FDH, day) was calculated using the previously determined kFDH (Equation 1), according to Equation 2:

t1/2,FDH=ln2kFDH(2)

5.4 NADH stability assay

Measurements of the stability of the coenzyme NADH in various solvent systems (0.03 mg mL−1) were monitored for up to 14 days. The samples were stored in the dark at 25°C and the absorption spectra in the range from 230 to 400 nm were recorded regularly using a UV-Vis spectrophotometer (SpectraMax® ABS Plus, Molecular Devices, San Jose, CA, United States). Each measurement was carried out in triplicate. The absorption spectrum of NADH shows a characteristic peak with an absorption maximum at 260 nm, which is due to the adenosine monophosphate moiety, and another peak at 340 nm, which is due to the neutral nicotinamide moiety. The decrease in absorbance at 340 nm followed first-order kinetics (Equation 3):

kNADH=1tlnA0At(3)

where kNADH is the first-order degradation rate constant (day-1), A is the absorbance at 340 nm, either at time zero, A0, or at time t, At (Wu et al., 1986).

The NADH half-life (t1/2, NADH, day) was calculated using the previously determined kNADH (Equation 3), according to Equation 4:

t1/2,NADH=ln2kNADH(4)

Moreover, the stability of NADH (25 mg mL−1), added in tested solvent systems after their saturation with CO2, was assessed to evaluate the impact of dissolved CO2 on the coenzyme. Monitoring was accomplished by measuring the absorbance decrease at 340 nm using a UV-Vis spectrophotometer for 90 min. The absorbance decrease at 340 nm followed the first-order kinetics, and the first-order rate degradation constants (kNADH*, day-1) were calculated in the same manner as described above. For the FDH-catalyzed reduction of CO2 in the reference buffer and ChCl:Gly80%B, NADH stability was monitored throughout the entire course of the reaction (6 h).

5.5 CO2 solubility measurements

The CO2 solubility measurement, based on the setup described by Obert and Dave (Obert and Dave, 1999) was conducted in a glass tube where CO2 was bubbled through the tested solvent systems with DESs and the reference buffer (V = 5 mL) for 90 min (time sufficient to reach a saturation point) at a flow rate of 100 mL min−1. The gas was introduced using a small nozzle with an approximate diameter of 1 mm. To prevent CO2 loss, the glass tube was sealed with parafilm. The concentration of CO2 was continuously monitored using a Mettler Toledo (Greifensee, Switzerland) CO2 sensor InPro 5000i/120.

5.6 FDH-catalysed reduction of CO2

Prior to the reaction initiation, both 50 mM potassium phosphate buffer (pH 7.5) and ChCl:Gly80%B were saturated with CO₂ by bubbling both solvents for 90 min, resulting in CO₂ dissolved concentrations of 1,029 mg mL−1 and 1,057 mg mL−1, respectively. The enzymatic reaction was initiated by adding NADH and FDH to the CO2-saturated buffer or ChCl:Gly80%B to reach final concentrations of 16.4 mg mL−1 and 28.8 mg mL−1, respectively. The reactions occurred at a room temperature using a magnetic stirrer in small tubes with a working volume of 250 µL. Aliquots were sampled at specific intervals to monitor NADH consumption, which was measured spectrophotometrically at 340 nm, indicating enzyme activity in reducing CO₂. To validate the results and assess the stability of NADH under experimental conditions, control reactions were conducted without the FDH enzyme. These control mixtures allowed us to monitor NADH stability in both the buffer and ChCl:Gly80%B throughout the reaction period, as described in Section 5.4. To verify that the consumed NADH was utilized for formate synthesis, the concentration of formate was additionally determined using the method described by Lang and Lang (Singh et al., 2018). Briefly, samples (25 µL) containing formate were mixed with 50 µL of solution A, 2.5 µL of solution B, and 175 µL of 100% acetic anhydride. The mixture was incubated at 50°C for 2 h with occasional mixing. Formation of red color was subsequently measured spectrophotometrically at 515 nm using the SpectraMax® ABS Plus (Molecular Devices, San Jose, CA, United States). Solution A was prepared by dissolving 0.5 g of citric acid and 10 g of acetamide in 100 mL of isopropanol; solution B was prepared by dissolving 30 g of sodium acetate in 100 mL of water. For standard calibration, sodium formate dissolved in 50 mM potassium phosphate buffer (pH 7.5) was used.

5.7 Piecewise linear regression modelling

It was assumed that FDH activity, FDH and NADH stability can be described as a function of the DES physical properties and σ-profile of the mixture, expressed by a set of Simix descriptors:

AR or t1/2,FDH or t1/2,NADH=fρ,η,pH,Smix1,Smix2,Smix3,Smix4,Smix5,Smix6,Smix7,Smix8,Smix9,Smix10

Piecewise linear regression (PLR) models have been used to describe the relationship between input and output variables. Input variables were selected based on the Spearman correlation matrix (Equation 5).

AR or t1/2,FDH or t1/2,NADH=(b01+b11·ρ+b21·η+b31·pH+b41·Smix1+b51·Smix2+b61·Smix3+b71·Smix4+b81·Smix5+b91·Smix6+b101·Smix7+b111·Smix8+b121·Smix9+b131·Smix10)·AR or t1/2,FDH or t1/2,NADHbn+(b02+b12·ρ+b22·η+b32·pH+b42·Smix1+b52·Smix2+b62·Smix3+b72·Smix4+b82·Smix5+b92·Smix6+b102·Smix7+b112·Smix8+b122·Smix9+b132·Smix10)·AR or t1/2,FDH or t1/2,NADH>bn(5)

The PLR technique is based on estimating the parameters of two linear regression equations: one for dependent variable values (Y) less than or equal to the breakpoint (bn) and the other for dependent variable values (Y) higher than the breakpoint. The PLR parameters were estimated using the Levenberg-Marquardt algorithm implemented in the software Statistica 14.0 (Tibco Software Inc., Palo Alto, United States). The data set (63 data points for each output variable) was randomly split 70:30 into a calibration and a prediction data set. The applicability of the developed calibration models was estimated using the coefficient of determination for calibration (Rcal2), the adjusted coefficient of determination for calibration (Rcal2adj), and cosmo (RMSEC). Predictive performance of the models was estimated using the coefficient of determination for prediction (Rpred2), the adjusted coefficient of determination for calibration (Rpred2adj), the root mean square error of prediction (RMSEP), the ratio of prediction to deviation (RPD) and the ratio of the error range (RER) (Fearn, 2002).

5.8 Estimation of the productivity of the biocatalytic process

The productivity of the FDH-catalysed CO2 conversion to formate was estimated based on the enzyme reaction kinetics including the enzyme inactivation rate (Bisswanger, 2017). Michaelis–Menten kinetic parameters were estimated by fitting the NADH concentration profiles to the differential equation (Equation 6) using WR Mathematica 10.0 (Wolfram Research, Champaign, United States). The change in substrate concentration over time without enzyme inactivation reads:

dcNADHdt=vmax·cNADHKS+cNADH(6)

Considering enzyme inactivation, the following expression for the reaction rate v(t) (Equation 7) is obtained:

vt=vmax·cNADHKS+cNADH·ekFDH·t(7)

The productivity of the FDH-catalyzed CO2 conversion to formate can be calculated by integrating (Equation 7) from the beginning of the reaction (t = 0) to its end (t = tf) to obtain (Equation 8):

productivity=0tfvmax·cNADHKS+cNADH·ekFDH·tdt=vmax·cNADHKS+cNADH·1ekFDH·tfkFDH(8)

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

ML: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Writing–original draft, Writing–review and editing. KŠ: Data curation, Investigation, Writing–review and editing. MN: Data curation, Investigation, Writing–review and editing. AT: Writing–review and editing, Conceptualization, Formal Analysis, Methodology, Software. AD: Writing–review and editing, Investigation, Visualization. MR: Investigation, Writing–review and editing, Visualization. IR: Writing–review and editing, Data curation. PŽ-P: Writing–review and editing, Funding acquisition, Resources, Supervision. WK: Writing–review and editing, Resources. MC: Resources, Writing–original draft, Writing–review and editing, Conceptualization, Data curation, Formal Analysis, Funding acquisition, Methodology, Supervision, Visualization.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the Croatian Science Foundation (grant number IPS- 2022-02-3938), and the Slovenian Research and Innovation Agency through Grants J4-4562 and P2-0191. PŽ-P was also supported through the Interreg Central Europe project CE0200857 GreenChemForCE, co-funded by the European Union.

Acknowledgments

The authors are grateful for the continuous support of the Croatian Science Foundation for the research related to the application of deep eutectic solvents in biotechnology.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2024.1467810/full#supplementary-material

References

Abbott, A. P., Capper, G., Davies, D. L., Rasheed, R. K., and Tambyrajah, V. (2003). Novel solvent properties of choline chloride/urea mixtures. Chem. Commun., 70–71. doi:10.1039/B210714G

PubMed Abstract | CrossRef Full Text | Google Scholar

Abildskov, J., Van Leeuwen, M. B., Boeriu, C. G., and Van Den Broek, L. A. M. (2013). Computer-aided solvent screening for biocatalysis. J. Mol. Catal. B Enzym. 85 (86), 200–213. doi:10.1016/J.MOLCATB.2012.09.012

CrossRef Full Text | Google Scholar

Abranches, D. O., and Coutinho, J. A. P. (2022). Type V deep eutectic solvents: design and applications. Curr. Opin. Green Sustain. Chem. 35, 100612. doi:10.1016/J.COGSC.2022.100612

CrossRef Full Text | Google Scholar

Abranches, D. O., Zhang, Y., Maginn, E. J., and Colón, Y. J. (2022). Sigma profiles in deep learning: towards a universal molecular descriptor. Chem. Commun. 58, 5630–5633. doi:10.1039/D2CC01549H

CrossRef Full Text | Google Scholar

Benguerba, Y., Alnashef, I. M., Erto, A., Balsamo, M., and Ernst, B. (2019). A quantitative prediction of the viscosity of amine based DESs using Sσ-profile molecular descriptors. J. Mol. Struct. 1184, 357–363. doi:10.1016/J.MOLSTRUC.2019.02.052

CrossRef Full Text | Google Scholar

Bisswanger, H. (2017). Enzyme kinetics. doi:10.1002/9783527806461

CrossRef Full Text | Google Scholar

Biswas, R., Metya, A. K., Abebe, K. M., Gedf, S. A., and Melese, B. T. (2023). Carbon dioxide solubility in choline chloride-based deep eutectic solvents under diverse conditions. J. Mol. Model. 29, 236. doi:10.1007/S00894-023-05643-Z

PubMed Abstract | CrossRef Full Text | Google Scholar

Bittner, J. P., Smirnova, I., and Jakobtorweihen, S. (2024). Investigating biomolecules in deep eutectic solvents with molecular dynamics simulations: current state, challenges and future perspectives. Molecules 29, 703. doi:10.3390/MOLECULES29030703

PubMed Abstract | CrossRef Full Text | Google Scholar

Bittner, J. P., Zhang, N., Huang, L., Domínguez De María, P., Jakobtorweihen, S., and Kara, S. (2022). Impact of deep eutectic solvents (DESs) and individual des components on alcohol dehydrogenase catalysis: Connecting experimental data and molecular dynamics simulations. Green Chem. 24, 1120–1131. doi:10.1039/d1gc04059f

CrossRef Full Text | Google Scholar

Calzadiaz-Ramirez, L., and Meyer, A. S. (2022). Formate dehydrogenases for CO2 utilization. Curr. Opin. Biotechnol. 73, 95–100. doi:10.1016/J.COPBIO.2021.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Castillo, E., Casas-Godoy, L., and Sandoval, G. (2016). Medium-engineering: a useful tool for modulatinglipase activity and selectivity. Biocatalysis 1, 178–188. doi:10.1515/boca-2015-0013

CrossRef Full Text | Google Scholar

Chang, C.-W., Laird, D. A., Mausbach, M. J., and Hurburgh, C. R. (2001). Near-Infrared reflectance spectroscopy–principal components regression analyses of soil properties. Soil Sci. Soc. Am. J. 65, 480–490. doi:10.2136/SSSAJ2001.652480X

CrossRef Full Text | Google Scholar

Cheng, H., Liu, C., Zhang, J., Chen, L., Zhang, B., and Qi, Z.(2018). Screening deep eutectic solvents for extractive desulfurization of fuel based on COSMO-RS model.Chem. Eng. Process. - Process Intensif.125, 246–252. doi:10.1016/J.CEP.2018.02.006

CrossRef Full Text | Google Scholar

Cui, Z., Enjome Djocki, A. V., Yao, J., Wu, Q., Zhang, D., Nan, S., et al. (2021). COSMO-SAC-supported evaluation of natural deep eutectic solvents for the extraction of tea polyphenols and process optimization. J. Mol. Liq. 328, 115406. doi:10.1016/j.molliq.2021.115406

CrossRef Full Text | Google Scholar

Cvjetko Bubalo, M., Vidović, S., Radojčić Redovniković, I., and Jokić, S. (2015). Green solvents for green technologies. J. Chem. Technol. Biotechnol. 90, 1631–1639. doi:10.1002/jctb.4668

CrossRef Full Text | Google Scholar

Damjanović, A., Logarušić, M., Tumir, L.-M., Andreou, T., Bubalo, M. C., and Redovniković, I. R. (2024). Enhancing protein stability under stress: osmolyte-based deep eutectic solvents as a biocompatible and robust stabilizing medium for lysozyme under heat and cold shock. Phys. Chem. Chem. Phys. 26, 21040–21051. doi:10.1039/D4CP02275K

PubMed Abstract | CrossRef Full Text | Google Scholar

Domingues, L., Duarte, A. R. C., and Jesus, A. R. (2024). How can deep eutectic systems promote greener processes in medicinal chemistry and drug discovery? Pharmaceuticals 17, 221. doi:10.3390/PH17020221

PubMed Abstract | CrossRef Full Text | Google Scholar

Elgharbawy, A. A. M., Syed Putra, S. S., Khan, H. W., Azmi, N. A. N., Sani, M. S. A., Ab llah, N., et al. (2023). Menthol and fatty acid-based hydrophobic deep eutectic solvents as media for enzyme activation. Process 11, 547. doi:10.3390/PR11020547

CrossRef Full Text | Google Scholar

Fearn, T. (2002). Assessing calibrations: SEP, RPD, RER and R2. NIR News, 13, 12–13. doi:10.1255/NIRN.689

CrossRef Full Text | Google Scholar

Gajardo-Parra, N. F., Rodríguez, G., Arroyo-Avirama, A. F., Veliju, A., Happe, T., Canales, R. I., et al. (2023). Impact of deep eutectic solvents on kinetics and folding stability of formate dehydrogenase. Processes 11, 2815. doi:10.3390/pr11102815

CrossRef Full Text | Google Scholar

Gorke, J. T., Srienc, F., and Kazlauskas, R. J. (2008). Hydrolase-catalyzed biotransformations in deep eutectic solvents. Chem. Commun., 1235–1237. doi:10.1039/B716317G

PubMed Abstract | CrossRef Full Text | Google Scholar

Hammond, O. S., Bowron, D. T., Edler, K. J., Hammond, S., Edler, K. J., and B owron, D. T. (2017). The effect of water upon deep eutectic solvent nanostructure: an unusual transition from ionic mixture to aqueous solution. Angew. Chem. Int. Ed. 56, 9782–9785. doi:10.1002/ANIE.201702486

CrossRef Full Text | Google Scholar

Hansen, B. B., Spittle, S., Chen, B., Poe, D., Zhang, Y., Klein, J. M., et al. (2021). Deep eutectic solvents: a review of fundamentals and applications. Chem. Rev. 121, 1232–1285. doi:10.1021/acs.chemrev.0c00385

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussain, S., Fangwei, Z., Siddiqi, A. F., Ali, Z., and Shabbir, M. S. (2018). Structural equation model for evaluating factors affecting quality of social infrastructure projects. Sustain 10, 1415. doi:10.3390/SU10051415

CrossRef Full Text | Google Scholar

Klamt, A. (2005). COSMO-RS: from quantum chemistry to fluid phase thermodynamics and drug design. Chem. Eng. Data 51 (4)–1480. doi:10.1021/je0602317

CrossRef Full Text | Google Scholar

Laitinen, A. T., Parsana, V. M., Jauhiainen, O., Huotari, M., Van Den Broeke, L. J. P., De Jong, W., et al. (2021). Liquid-liquid extraction of formic acid with 2-methyltetrahydrofuran: experiments, process modeling, and economics. Ind. Eng. Chem. Res. 60, 5588–5599. doi:10.1021/acs.iecr.1c00159

PubMed Abstract | CrossRef Full Text | Google Scholar

Lemaoui, T., Boublia, A., Darwish, A. S., Alam, M., Park, S., Jeon, B. H., et al. (2022). Predicting the surface tension of deep eutectic solvents using artificial neural networks. ACS Omega 7, 32194–32207. doi:10.1021/acsomega.2c03458

PubMed Abstract | CrossRef Full Text | Google Scholar

Lemaoui, T., Darwish, A. S., Hammoudi, N. E. H., Abu Hatab, F., Attoui, A., Alnashef, I. M., et al. (2020). Prediction of electrical conductivity of deep eutectic solvents using COSMO-RS sigma profiles as molecular descriptors: a quantitative structure-property relationship study. Ind. Eng. Chem. Res. 59, 13343–13354. doi:10.1021/ACS.IECR.0C02542

CrossRef Full Text | Google Scholar

Leron, R. B., Caparanga, A., and Li, M. H. (2013). Carbon dioxide solubility in a deep eutectic solvent based on choline chloride and urea at T = 303.15–343.15 K and moderate pressures. J. Taiwan Inst. Chem. Eng. 44, 879–885. doi:10.1016/J.JTICE.2013.02.005

CrossRef Full Text | Google Scholar

Leron, R. B., and Li, M. H. (2013). Solubility of carbon dioxide in a choline chloride–ethylene glycol based deep eutectic solvent. Thermochim. Acta 551, 14–19. doi:10.1016/J.TCA.2012.09.041

CrossRef Full Text | Google Scholar

Lomba, L., Ribate, M. P., Zaragoza, E., Concha, J., Garralaga, M. P., Errazquin, D., et al. (2021). Deep eutectic solvents: are they safe? Appl. Sci. 11, 10061. doi:10.3390/APP112110061

CrossRef Full Text | Google Scholar

Lozano, P., Diego, T. D., and Iborra, J. L. (2010). Biocatalytic processes using ionic liquids and supercritical carbon dioxide. Handb. Green Chem., 51–73. doi:10.1002/9783527628698.HGC027

CrossRef Full Text | Google Scholar

Ma, G., Zhang, Z., Lu, Z., Qi, Y., Nian, B., and Hu, Y. (2024). Hydrophobic deep eutectic solvents-lipase synergistically catalyze the synthesis of Vitamin E succinate via hydrogen bonds. J. Mol. Liq. 394, 123711. doi:10.1016/J.MOLLIQ.2023.123711

CrossRef Full Text | Google Scholar

Maugeri, Z., Leitner, W., and Domínguez De María, P. (2012). Practical separation of alcohol–ester mixtures using Deep-Eutectic-Solvents. Tetrahedron Lett. 53, 6968–6971. doi:10.1016/J.TETLET.2012.10.044

CrossRef Full Text | Google Scholar

Menegatti, T., Kopač, T., and Žnidaršič-Plazl, P. (2024). Tuning mechanical characteristics and permeability of alginate hydrogel by polyvinyl alcohol and deep eutectic solvent addition. Bioengineering 11, 371. doi:10.3390/bioengineering11040371

PubMed Abstract | CrossRef Full Text | Google Scholar

Mourelle-Insua, Á., Lavandera, I., and Gotor-Fernández, V. (2019). A designer natural deep eutectic solvent to recycle the cofactor in alcohol dehydrogenase-catalysed processes. Green Chem. 21, 2946–2951. doi:10.1039/C9GC00318E

CrossRef Full Text | Google Scholar

Nolasco, M. M., Pedro, S. N., Vilela, C., Vaz, P. D., Ribeiro-Claro, P., Rudić, S., et al. (2022). Water in deep eutectic solvents: new insights from inelastic neutron scattering spectroscopy. Front. Phys. 10, 834571. doi:10.3389/fphy.2022.834571

CrossRef Full Text | Google Scholar

Obert, R., and Dave, B. C. (1999). Enzymatic conversion of carbon dioxide to methanol: enhanced methanol production in silica sol-gel matrices. J. Am. Chem. Soc. 121, 12192–12193. doi:10.1021/ja991899r

CrossRef Full Text | Google Scholar

Panić, M., Cvjetko Bubalo, M., and Radojčić Redovniković, I. (2021). Designing a biocatalytic process involving deep eutectic solvents. J. Chem. Technol. Biotechnol. 96, 14–30. doi:10.1002/JCTB.6545

CrossRef Full Text | Google Scholar

Panić, M., Radović, M., Cvjetko Bubalo, M., Radošević, K., Rogošić, M., Coutinho, J. A. P., et al. (2022). Prediction of pH value of aqueous acidic and basic deep eutectic solvent using COSMO-RS σ profiles’ molecular descriptors. Mol 27, 4489. doi:10.3390/MOLECULES27144489

CrossRef Full Text | Google Scholar

Pätzold, M., Burek, B. O., Liese, A., Bloh, J. Z., and Holtmann, D. (2019a). Product recovery of an enzymatically synthesized (−)-menthol ester in a deep eutectic solvent. Bioprocess Biosyst. Eng. 42, 1385–1389. doi:10.1007/s00449-019-02125-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Pätzold, M., Weimer, A., Liese, A., and Holtmann, D. (2019b). Optimization of solvent-free enzymatic esterification in eutectic substrate reaction mixture. Biotechnol. Rep. 22, e00333. doi:10.1016/J.BTRE.2019.E00333

CrossRef Full Text | Google Scholar

Pohar, A., Žnidaršič-Plazl, P., and Plazl, I. (2012). Integrated system of a microbioreactor and a miniaturized continuous separator for enzyme catalyzed reactions. Chem. Eng. J. 189–190, 376–382. doi:10.1016/J.CEJ.2012.02.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Quaid, T., and Reza, T. (2023). COSMO prediction of siloxane compounds absorption on type 3 and type 5 deep eutectic solvents. Chem. Eng. J. Adv. 14, 100489. doi:10.1016/J.CEJA.2023.100489

CrossRef Full Text | Google Scholar

Radošević, K., Cvjetko Bubalo, M., Gaurina Srček, V., Grgas, D., Landeka Dragičević, T., and Radojčić, R. I. (2015). Evaluation of toxicity and biodegradability of choline chloride based deep eutectic solvents. Ecotoxicol. Environ. Saf. 112, 46–53. doi:10.1016/J.ECOENV.2014.09.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Radović, M., Hok, L., Panić, M., Cvjetko Bubalo, M., Vianello, R., Vinković, M., et al. (2022). Deep eutectic solvents as a stabilising medium for NAD coenzyme: unravelling the mechanism behind coenzyme stabilisation effect. Green Chem. 24, 7661–7674. doi:10.1039/d2gc02656b

CrossRef Full Text | Google Scholar

Radović, M., Jurinjak Tušek, A., Kroutil, W., Reiter, T., Radojčić Redovniković, I., and Cvjetko Bubalo, M. (2024). Rational design of deep eutectic solvents for the stabilization of dehydrogenases: an artificial neural network prediction approach. Front. Chem. 12, 1436049. doi:10.3389/FCHEM.2024.1436049

PubMed Abstract | CrossRef Full Text | Google Scholar

Reda, T., Plugge, C. M., Abram, N. J., and Hirst, J. (2008). Reversible interconversion of carbon dioxide and formate by an electroactive enzyme. Proc. Natl. Acad. Sci. U. S. A. 105, 10654–10658. doi:10.1073/pnas.0801290105

PubMed Abstract | CrossRef Full Text | Google Scholar

Rente, D., Cvjetko Bubalo, M., Panić, M., Paiva, A., Caprin, B., Radojčić Redovniković, I., et al. (2022). Review of deep eutectic systems from laboratory to industry, taking the application in the cosmetics industry as an example. J. Clean. Prod. 380, 135147. doi:10.1016/J.JCLEPRO.2022.135147

CrossRef Full Text | Google Scholar

Sanchez-Fernandez, A., Basic, M., Xiang, J., Prevost, S., Jackson, A. J., and Dicko, C. (2022). Hydration in deep eutectic solvents induces non-monotonic changes in the conformation and stability of proteins. J. Am. Chem. Soc. 144, 23657–23667. doi:10.1021/jacs.2c11190

PubMed Abstract | CrossRef Full Text | Google Scholar

Sapir, L., and Harries, D. (2020). Restructuring a deep eutectic solvent by water: the nanostructure of hydrated choline chloride/urea. J. Chem. Theory Comput. 16, 3335–3342. doi:10.1021/acs.jctc.0c00120

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheldon, R. A., Bode, M. L., and Mathebula, N. (2023). Green and sustainable solvents for biocatalytic oxidations. Curr. Opin. Green Sustain. Chem. 39, 100741. doi:10.1016/J.COGSC.2022.100741

CrossRef Full Text | Google Scholar

Sheldon, R. A., and Pereira, P. C. (2017). Biocatalysis engineering: the big picture. Chem. Soc. Rev. 46, 2678–2691. doi:10.1039/C6CS00854B

PubMed Abstract | CrossRef Full Text | Google Scholar

Sim, J., McGoverin, C., Oey, I., Frew, R., and Kebede, B. (2023). Near-infrared reflectance spectroscopy accurately predicted isotope and elemental compositions for origin traceability of coffee. Food Chem. 427, 136695. doi:10.1016/J.FOODCHEM.2023.136695

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, R. K., Singh, R., Sivakumar, D., Kondaveeti, S., Kim, T., Li, J., et al. (2018). Insights into cell-free conversion of CO2 to chemicals by a multienzyme cascade reaction. ACS Catal. 8, 11085–11093. doi:10.1021/acscatal.8b02646

CrossRef Full Text | Google Scholar

Slater, C. S., Savelski, M. J., Hitchcock, D., and Cavanagh, E. J. (2016). Environmental analysis of the life cycle emissions of 2-methyl tetrahydrofuran solvent manufactured from renewable resources. J. Environ. Sci. Health. A. Tox. Hazard. Subst. Environ. Eng. 51, 487–494. doi:10.1080/10934529.2015.1128719

PubMed Abstract | CrossRef Full Text | Google Scholar

Taklimi, S. M., Divsalar, A., Ghalandari, B., Ding, X., Di Gioia, M. L., Omar, K. A., et al. (2023). Effects of deep eutectic solvents on the activity and stability of enzymes. J. Mol. Liq. 377, 121562. doi:10.1016/J.MOLLIQ.2023.121562

CrossRef Full Text | Google Scholar

Tishkov, V. I., Galkin, A. G., Marchenko, G. N., Egorova, O. A., Sheluho, D. V., Kulakova, L. B., et al. (1993). Catalytic properties and stability of a Pseudomonas sp.101 formate dehydrogenase mutants containing cys-255-ser and cys-255-met replacements. Biochem. Biophys. Res. Commun. 192, 976–981. doi:10.1006/BBRC.1993.1511

PubMed Abstract | CrossRef Full Text | Google Scholar

Torregrosa-Crespo, J., Marset, X., Guillena, G., Ramón, D. J., and María Martínez-Espinosa, R. (2020). New guidelines for testing “Deep eutectic solvents” toxicity and their effects on the environment and living beings. Sci. Total Environ. 704, 135382. doi:10.1016/J.SCITOTENV.2019.135382

PubMed Abstract | CrossRef Full Text | Google Scholar

van Schie, M. M. C. H., Spöring, J. D., Bocola, M., Domínguez de María, P., and Rother, D. (2021). Applied biocatalysis beyond just buffers – from aqueous to unconventional media. Options and guidelines. Green Chem. 23, 3191–3206. doi:10.1039/D1GC00561H

PubMed Abstract | CrossRef Full Text | Google Scholar

Venkatraman, V., Evjen, S., Knuutila, H. K., Fiksdahl, A., and Alsberg, B. K. (2018). Predicting ionic liquid melting points using machine learning. J. Mol. Liq. 264, 318–326. doi:10.1016/J.MOLLIQ.2018.03.090

CrossRef Full Text | Google Scholar

Villa, R., Alvarez, E., Porcar, R., Garcia-Verdugo, E., Luis, S. V., and Lozano, P. (2019). Ionic liquids as an enabling tool to integrate reaction and separation processes. Green Chem. 21, 6527–6544. doi:10.1039/C9GC02553G

CrossRef Full Text | Google Scholar

Villa, R., Nieto, S., Donaire, A., and Lozano, P. (2023). Direct biocatalytic processes for CO2 capture as a green tool to produce value-added chemicals. Mol 28, 5520. doi:10.3390/MOLECULES28145520

CrossRef Full Text | Google Scholar

Wichmann, R., and Vasic-Racki, D. (2005). Cofactor regeneration at the lab scale. Adv. Biochem. Eng. Biotechnol. 92, 225–260. doi:10.1007/B98911

PubMed Abstract | CrossRef Full Text | Google Scholar

Woodley, J. M. (2022). Ensuring the sustainability of biocatalysis. ChemSusChem 15, e202102683. doi:10.1002/CSSC.202102683

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, B. P., Wen, Q., Xu, H., and Yang, Z. (2014). Insights into the impact of deep eutectic solvents on horseradish peroxidase: activity, stability and structure. J. Mol. Catal. B. Enzym. Complete 101, 101–107. doi:10.1016/J.MOLCATB.2014.01.001

CrossRef Full Text | Google Scholar

Wu, J. T., Wu, L. H., and Knight, J. A. (1986). Stability of NADPH: effect of various factors on the kinetics of degradation. Clin. Chem. 32, 314–319. doi:10.1093/CLINCHEM/32.2.314

PubMed Abstract | CrossRef Full Text | Google Scholar

Zachos, I., Nowak, C., and Sieber, V. (2019). Biomimetic cofactors and methods for their recycling. Curr. Opin. Chem. Biol. 49, 59–66. doi:10.1016/J.CBPA.2018.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, N., Domínguez de María, P., and Kara, S. (2024). Biocatalysis for the synthesis of active pharmaceutical ingredients in deep eutectic solvents: state-of-the-art and prospects. Catal 14, 84. doi:10.3390/CATAL14010084

CrossRef Full Text | Google Scholar

Zhang, Z., Wang, H., Nie, Y., Zhang, X., and Ji, X. (2022). Natural deep eutectic solvents enhanced electro-enzymatic conversion of CO2 to methanol. Front. Chem. 10, 1–8. doi:10.3389/fchem.2022.894106

CrossRef Full Text | Google Scholar

Zhang, Z., Xu, B. H., Luo, J., Von Solms, N., He, H., Zhang, Y., et al. (2018). Ionic liquids as bifunctional cosolvents enhanced co2 conversion catalysed by nadh-dependent formate dehydrogenase. Catalysts 8, 304. doi:10.3390/catal8080304

CrossRef Full Text | Google Scholar

Žnidaršič-Plazl, P. (2021a). Let the biocatalyst flow. Acta Chim. Slov. 68, 1–16. doi:10.17344/ACSI.2020.6488

PubMed Abstract | CrossRef Full Text | Google Scholar

Žnidaršič-Plazl, P. (2021b). Biocatalytic process intensification via efficient biocatalyst immobilization, miniaturization, and process integration. Curr. Opin. Green Sustain. Chem. 32, 100546. doi:10.1016/J.COGSC.2021.100546

CrossRef Full Text | Google Scholar

Glossary

DES Deep eutectic solvent

FDH Formate dehydrogenase

NADH Nicotinamide adenine dinucleotide

HBA Hydrogen bond acceptor

HBD Hydrogen bond donor

ChCl Choline chloride

B Betaine

Gly Glycerol

%W % (w/w) of water

%B % (w/w) of buffer

η Dynamic viscosity

ρ Density

Simix σ-profile of the mixture

COSMO-RS Conductor-like Screening Model for Real Solvents

QSPR Quantitative structure property relationships

cs concentration of dissolved CO2

AR relative enzyme activity

ARes residual enzyme activity

t1/2, FDH FDH half-life

t1/2, NADH NADH half-life

kFDH first-order degradation rate constant of FDH

kNADH first-order degradation rate constant of NADH

kNADH* first-order degradation rate constant of NADH in CO2 saturated solutions

vmax maximum reaction rate

Ks NADH saturation constant

v(t) reaction rate

PLR piecewise linear regression

Rpred2 coefficient of determination for prediction

Rpred2adj adjusted coefficient of determination for prediction

RMSEP root mean square error of prediction

RPD ratio of prediction to deviation

RER ratio of the error range

Rcal2 coefficient of determination for calibration

Rcal2adj adjusted coefficient of determination for calibration

RMSEC root mean square error for calibration

Keywords: deep eutectic solvents, rational design, biocatalysis, mathematical modelling, QSAR, formate dehydrogenase, NADH, CO2 conversion

Citation: Logarušić M, Šubar K, Nikolić M, Jurinjak Tušek A, Damjanović A, Radović M, Radojčić Redovniković I, Žnidaršič-Plazl P, Kroutil W and Cvjetko Bubalo M (2024) Harnessing the potential of deep eutectic solvents in biocatalysis: design strategies using CO2 to formate reduction as a case study. Front. Chem. 12:1467810. doi: 10.3389/fchem.2024.1467810

Received: 20 July 2024; Accepted: 10 October 2024;
Published: 25 October 2024.

Edited by:

Richard Villanneau, Sorbonne Université, France

Reviewed by:

Andrea Mezzetta, University of Pisa, Italy
María Salomé Álvarez, University of Vigo, Spain

Copyright © 2024 Logarušić, Šubar, Nikolić, Jurinjak Tušek, Damjanović, Radović, Radojčić Redovniković, Žnidaršič-Plazl, Kroutil and Cvjetko Bubalo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Marina Cvjetko Bubalo, bWFyaW5hLmN2amV0a28uYnViYWxvQHBiZi51bml6Zy5ocg==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.