Rational Design of Thermostable Carbonic Anhydrase Mutants Using Molecular Dynamics Simulations

Computational design of α-amylase from Bacillus licheniformis to increase its activity and stability at high temperatures

The thermostable α-amylase derived from Bacillus licheniformis (BLA) has multiple advantages, including enhancing the mass transfer rate and by reducing microbial contamination in starch hydrolysis. Nonetheless, the application of BLA is constrained by the accessibility and stability of enzymes capable of achieving high conversion rates at elevated temperatures. Moreover, the thermotolerance of BLA requires further enhancement. Here, we developed a computational strategy for constructing small and smart mutant libraries to identify variants with enhanced thermostability. Initially, molecular dynamics (MD) simulations were employed to identify the regions with high flexibility. Subsequently, FoldX, a computational design predictor, was used to design mutants by rigidifying highly flexible residues, whereas the simultaneous decrease in folding free energy assisted in improving thermostability. Through the utilization of MD and FoldX, residues K251, T277, N278, K319, and E336, situated at a distance of 5 Å from the catalytic triad, were chosen for mutation. Seventeen mutants were identified and characterized by evaluating enzymatic characteristics and kinetic parameters. The catalytic efficiency of the E271L/N278K mutant reached 184.1 g L-1 s-1, which is 1.88-fold larger than the corresponding value determined for the WT. Furthermore, the most thermostable mutant, E336S, exhibited a 1.43-fold improvement in half-life at 95 ℃, compared with that of the WT. This study, by combining computational simulation with experimental verification, establishes that potential sites can be computationally predicted to increase the activity and stability of BLA and thus provide a possible strategy by which to guide protein design.

Advancing sustainable biotechnology through protein engineering

The push for industrial sustainability benefits from the use of enzymes as a replacement for traditional chemistry. Biological catalysts, especially those that have been engineered for increased activity, stability, or novel function, and are often greener than alternative chemical approaches. This Review highlights the role of engineered enzymes (and identifies directions for further engineering efforts) in the application areas of greenhouse gas sequestration, fuel production, bioremediation, and degradation of plastic wastes.

Enantioselectivity in Vanadium-Dependent Haloperoxidases of Different Marine Sources for Sulfide Oxidation to Sulfoxides

This study explores the reasons behind the variations in the enantioselectivity of the sulfoxidation of methyl phenyl sulfide by marine-derived vanadium-dependent haloperoxidases (VHPOs). Twelve new VHPOs of marine organisms were overexpressed, purified, and tested for their ability to oxidize sulfide. Most of these marine enzymes exhibited nonenantioselective behavior, underscoring the uniqueness of AnVBPO from the brown seaweed Ascophyllum nodosum and CpVBPO from the red seaweed Corallina pilulifera, which produce (R)- and (S)-sulfoxides, respectively. The enantioselective sulfoxidation pathway is likely due to direct oxygen transfer within the VHPO active site. This was demonstrated through molecular docking and molecular dynamics simulations, which revealed differences in the positioning of sulfide within AnVBPO and CpVBPO, thus explaining their distinct enantioselectivities. Nonenantioselective VHPOs probably follow a different oxidation pathway, initiating with sulfide oxidation to form a positively charged radical. Further insights were gained from studying the catalytic effect of VO43- on H2O2-driven sulfoxidation. This research improves the understanding of VHPO-mediated sulfoxidation and aids in developing biocatalysts for sulfoxide synthesis.

Bacterial α-CAs: a biochemical and structural overview

Carbonic anhydrases belonging to the α-class are widely distributed in bacterial species. These enzymes have been isolated from bacteria with completely different characteristics including both Gram-negative and Gram-positive strains. α-CAs show a considerable similarity when comparing the biochemical, kinetic and structural features, with only small differences which reflect the diverse role these enzymes play in Nature. In this chapter, we provide a comprehensive overview on bacterial α-CA data, with a highlight to their potential biomedical and biotechnological applications.

Alpha Carbonic Anhydrase from Nitratiruptor tergarcus Engineered for Increased Activity and Thermostability

The development of carbon capture and storage technologies has resulted in a rising interest in the use of carbonic anhydrases (CAs) for CO2 fixation at elevated temperatures. In this study, we chose to rationally engineer the α-CA (NtCA) from the thermophilic bacterium Nitratiruptor tergarcus, which has been previously suggested to be thermostable by in silico studies. Using a combination of analyses with the DEEPDDG software and available structural knowledge, we selected residues in three regions, namely, the catalytic pocket, the dimeric interface and the surface, in order to increase thermostability and CO2 hydration activity. A total of 13 specific mutations, affecting seven amino acids, were assessed. Single, double and quadruple mutants were produced in Escherichia coli and analyzed. The best-performing mutations that led to improvements in both activity and stability were D168K, a surface mutation, and R210L, a mutation in the dimeric interface. Apart from these, most mutants showed improved thermostability, with mutants R210K and N88K_R210L showing substantial improvements in activity, up to 11-fold. Molecular dynamics simulations, focusing particularly on residue fluctuations, conformational changes and hydrogen bond analysis, elucidated the structural changes imposed by the mutations. Successful engineering of NtCA provided valuable lessons for further engineering of α-CAs.

Computer-directed rational design enhanced the thermostability of carbonyl reductase LsCR for the synthesis of ticagrelor precursor

Carbonyl reductases are useful for producing optically active alcohols from their corresponding prochiral ketones. Herein, we applied a computer-assisted strategy to increase the thermostability of a previously constructed carbonyl reductase, LsCRM4 (N101D/A117G/F147L/E145A), which showed an outstanding activity in the synthesis of the ticagrelor precursor (1S)-2-chloro-1-(3,4-difluorophenyl)ethanol. The stability changes introduced by mutations at the flexible sites were predicted using the computational tools FoldX, I-Mutant 3.0, and DeepDDG, which demonstrated that 12 virtually screened mutants could be thermally stable; 11 of these mutants exhibited increased thermostability. Then a superior mutant LsCRM4-V99L/D150F was screened out from the library that was constructed by iteratively combining the beneficial sites, which showed a 78% increase in activity and a 17.4°C increase in melting temperature compared to LsCRM4. Our computer-assisted design and combinatorial strategy dramatically increased the efficiency of thermostable enzyme production.

A new horizon in the phosphorylated sites of AGA: the structural impact of C163S mutation in aspartylglucosaminuria through molecular dynamics simulation

Aspartylglucosaminuria (AGU) is a lysosomal storage disorder caused by insufficient aspartylglucosaminidase (AGA) activity leading to chronic neurodegeneration. We utilized the PhosphoSitePlus tool to identify the AGA protein's phosphorylation sites. The phosphorylation was induced on the specific residue of the three-dimensional AGA protein, and the structural changes upon phosphorylation were studied via molecular dynamics simulation. Furthermore, the structural behaviour of C163S mutation and C163S mutation with adjacent phosphorylation was investigated. We have examined the structural impact of phosphorylated forms and C163S mutation in AGA. Molecular dynamics simulations (200 ns) exposed patterns of deviation, fluctuation, and change in compactness of Y178 phosphorylated AGA protein (Y178-p), T215 phosphorylated AGA protein (T215-p), T324 phosphorylated AGA protein (T324-p), C163S mutant AGA protein (C163S), and C163S mutation with Y178 phosphorylated AGA protein (C163S-Y178-p). Y178-p, T215-p, and C163S mutation demonstrated an increase in intramolecular hydrogen bonds, leading to greater compactness of the AGA forms. Principle component analysis (PCA) and Gibbs free energy of the phosphorylated/C163S mutation structures exhibit transition in motion/orientation than Wild type (WT). T215-p may be more dominant among these than the other studied phosphorylated forms. It might contribute to hydrolyzing L-asparagine functioning as an asparaginase, thereby regulating neurotransmitter activity. This study revealed structural insights into the phosphorylation of Y178, T215, and T324 in AGA protein. Additionally, it exposed the structural changes of the C163S mutation and C163S-Y178-p of AGA protein. This research will shed light on a better understanding of AGA's phosphorylated mechanism.Communicated by Ramaswamy H. Sarma.

Biocatalysis of CO2 and CH4: Key enzymes and challenges

Mitigating greenhouse gas emissions is a critical challenge for promoting global sustainability. The utilization of CO2 and CH4 as substrates for the production of valuable products offers a promising avenue for establishing an eco-friendly economy. Biocatalysis, a sustainable process utilizing enzymes to facilitate biochemical reactions, plays a significant role in upcycling greenhouse gases. This review provides a comprehensive overview of the enzymes and associated reactions involved in the biocatalytic conversion of CO2 and CH4. Furthermore, the challenges facing the field are discussed, paving the way for future research directions focused on developing robust enzymes and systems for the efficient fixation of CO2 and CH4.

Synthesis, In Vivo Anticonvulsant Activity Evaluation and In Silico Studies of Some Quinazolin-4(3H)-One Derivatives

Two series, "a" and "b", each consisting of nine chemical compounds, with 2,3-disubstituted quinazolin-4(3H)-one scaffold, were synthesized and evaluated for their anticonvulsant activity. They were investigated as dual potential positive allosteric modulators of the GABAA receptor at the benzodiazepine binding site and inhibitors of carbonic anhydrase II. Quinazolin-4(3H)-one derivatives were evaluated in vivo (D1-3 = 50, 100, 150 mg/kg, administered intraperitoneally) using the pentylenetetrazole (PTZ)-induced seizure model in mice, with phenobarbital and diazepam, as reference anticonvulsant agents. The in silico studies suggested the compounds act as anticonvulsants by binding on the allosteric site of GABAA receptor and not by inhibiting the carbonic anhydrase II, because the ligands-carbonic anhydrase II predicted complexes were unstable in the molecular dynamics simulations. The mechanism targeting GABAA receptor was confirmed through the in vivo flumazenil antagonism assay. The pentylenetetrazole experimental anticonvulsant model indicated that the tested compounds, 1a-9a and 1b-9b, present a potential anticonvulsant activity. The evaluation, considering the percentage of protection against PTZ, latency until the onset of the first seizure, and reduction in the number of seizures, revealed more favorable results for the "b" series, particularly for compound 8b.

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

A Bibliometric Analysis: Current Perspectives and Potential Trends of Enzyme Thermostability from 1991-2022

Thermostability is considered a crucial parameter to evaluate the viability of enzymes in industrial applications. Over the past 31 years, many studies have been reported on the thermostability of enzymes. However, there is no systematic bibliometric analysis of publications on the thermostability of enzymes. In this study, 16,035 publications related to the thermostability of enzymes were searched and collected, showing an increasing annual trend. China contributed the most publications, while the United States had the highest citation count. International Journal of Biological Macromolecules is the most productive journal in the research field. Moreover, Chinese acad sci and Khosro Khajeh are the most active institutions and prolific authors in the field, respectively. Analysis of references with the strongest citation bursts and keyword co-occurrences, magnetic nanoparticles, metal-organic frameworks, molecular dynamics, and rational design are current hot spots and significant future research directions. This study is the first comprehensive bibliometric analysis summarizing trends and developments in enzyme thermostability research. Our findings could provide scholars with an understanding of the fundamental knowledge framework of the field and identify recent potential hotspots and research trends that could facilitate the discovery of collaboration opportunities.

Understanding the basis of thermostability for enzyme "Nanoluc" towards designing industry-competent engineered variants

As a leading contender in the study of luminescence, nanoluciferase has recently attracted attention and proven effective in a wide variety of research areas. Although numerous attempts have been made to improve activity, there has yet to be a thorough exploration of further possibilities to improve thermostability. In this study, protein engineering in tandem with molecular dynamics simulation at various temperatures (300 K, 400 K, 450 K and 500 K) was used to improve our understanding of nanoluciferase dynamics and identification of factors that could significantly enhance the thermostability. Based on these, three novel mutations have been narrowed down, which were hypothesised to improve thermostability. Root mean square deviation and root mean square fluctuation studies confirmed higher stability of mutant at high temperature. Solvent-accessible surface area and protein unfolding studies revealed a decreased tendency of mutant to unfold at higher temperatures. Further free energy landscape and principal component analysis was adapted to get deeper insights into the thermodynamic and structural behavior of these proteins at elevated temperature. Thus, this study provides a deeper insight into the dynamic factors for thermostability and introduces a novel, enhanced nanoluciferase candidate with potential use in industry.Communicated by Ramaswamy H. Sarma.

Improved Solubility and Stability of a Thermostable Carbonic Anhydrase via Fusion with Marine-Derived Intrinsically Disordered Solubility Enhancers

Carbonic anhydrase (CA), an enzyme catalyzing the reversible hydration reaction of carbon dioxide (CO2), is considered a promising biocatalyst for CO2 reduction. The α-CA of Thermovibrio ammonificans (taCA) has emerged as a compelling candidate due to its high thermostability, a critical factor for industrial applications. However, the low-level expression and poor in vitro solubility have hampered further utilization of taCA. Recently, these limitations have been addressed through the fusion of the NEXT tag, a marine-derived, intrinsically disordered small peptide that enhances protein expression and solubility. In this study, the solubility and stability of NEXT-taCA were further investigated. When the linker length between the NEXT tag and the taCA was shortened, the expression level decreased without compromising solubility-enhancing performance. A comparison between the NEXT tag and the NT11 tag demonstrated the NEXT tag's superiority in improving both the expression and solubility of taCA. While the thermostability of taCA was lower than that of the extensively engineered DvCA10, the NEXT-tagged taCA exhibited a 30% improvement in long-term thermostability compared to the untagged taCA, suggesting that enhanced solubility can contribute to enzyme thermostability. Furthermore, the bioprospecting of two intrinsically disordered peptides (Hcr and Hku tags) as novel solubility-enhancing fusion tags was explored, demonstrating their performance in improving the expression and solubility of taCA. These efforts will advance the practical application of taCA and provide tools and insights for enzyme biochemistry and bioengineering.

Rational design engineering of a more thermostable Sulfurihydrogenibium yellowstonense carbonic anhydrase for potential application in carbon dioxide capture technologies

Implementing hyperthermostable carbonic anhydrases into CO2 capture and storage technologies in order to increase the rate of CO2 absorption from the industrial flue gases is of great importance from technical and economical points of view. The present study employed a combination of in silico tools to further improve thermostability of a known thermostable carbonic anhydrase from Sulfurihydrogenibium yellowstonense. Experimental results showed that our rationally engineered K100G mutant not only retained the overall structure and catalytic efficiency but also showed a 3 °C increase in the melting temperature and a two-fold improvement in the enzyme half-life at 85 °C. Based on the molecular dynamics simulation results, rearrangement of salt bridges and hydrogen interactions network causes a reduction in local flexibility of the K100G variant. In conclusion, our study demonstrated that thermostability can be improved through imposing local structural rigidity by engineering a single-point mutation on the surface of the enzyme.

Enhanced Activity and Stability of an Acetyl Xylan Esterase in Hydrophilic Alcohols through Site-Directed Mutagenesis

Current demands for the development of suitable biocatalysts showing high process performance is stimulated by the need to replace current chemical synthesis with cleaner alternatives. A drawback to the use of biocatalysts for unique applications is their low performance in industrial conditions. Hence, enzymes with improved performance are needed to achieve innovative and sustainable biocatalysis. In this study, we report the improved performance of an engineered acetyl xylan esterase (BaAXE) in a hydrophilic organic solvent. The structure of BaAXE was partitioned into a substrate-binding region and a solvent-affecting region. Using a rational design approach, charged residues were introduced at protein surfaces in the solvent-affecting region. Two sites present in the solvent-affecting region, A12D and Q143E, were selected for site-directed mutagenesis, which generated the mutants MUT12, MUT143 and MUT12-143. The mutants MUT12 and MUT143 reported lower Km (0.29 mM and 0.27 mM, respectively) compared to the wildtype (0.41 mM). The performance of the mutants in organic solvents was assessed after enzyme incubation in various strengths of alcohols. The mutants showed improved activity and stability compared to the wild type in low strengths of ethanol and methanol. However, the activity of MUT143 was lost in 40% methanol while MUT12 and MUT12-143 retained over 70% residual activity in this environment. Computational analysis links the improved performance of MUT12 and MUT12-143 to novel intermolecular interactions that are absent in MUT143. This work supports the rationale for protein engineering to augment the characteristics of wild-type proteins and provides more insight into the role of charged residues in conferring stability.

Computational analysis of the effect of Gly100Ala mutation on the thermostability of SazCA

Maintaining the protein stability upon mutation is a challenging task in protein engineering. In the present computational study, we induced a single point Gly100Ala mutation in SazCA and examined the factors governing the stability and flexibility of the mutated form, and compared it to that of the wildtype using molecular dynamics simulations. We observed higher structural stability and lesser residual mobility in the mutated SazCA. Improved H-bonding due to Gly100Ala was observed. Ala100 was responsible for the increased helical contents in the mutated SazCA while Gly100 compromised the secondary structure contents in the wildtype. A strong network of salt bridges and high local ordering of the solvent molecules at the protein surface contributed to the enhanced stability of the mutated protein. Our simulations conclusively highlight Gly100Ala mutation as a step towards designing a more robust and thermostable SazCA.Communicated by Ramaswamy H. Sarma.

Microbial carbonic anhydrase mediated carbon capture, sequestration & utilization: A sustainable approach to delivering bio-renewables

In the recent scenario, anthropogenic interventions have alarmingly disrupted climatic conditions. The persistent change in the climate necessitates carbon neutrality. Efficient ways of carbon capture and sequestration could be employed for sustainable product generation. Carbonic anhydrase (CA) is an enzyme that reversibly catalyzes the conversion of carbon dioxide to bicarbonate ions, further utilized by cells for metabolic processes. Hence, utilizing CA from microbial sources for carbon sequestration and the corresponding delivery of bio-renewables could be the eco-friendly approach. Consequently, the microbial CA and amine-based carbon capture chemicals are synergistically applied to enhance carbon capture efficiency and eventual utilization. This review comprehends recent developments coupled with engineering techniques, especially in microbial CA, to create integrated systems for CO₂ sequestration. It envisages developing sustainable approaches towards mitigating environmental CO₂ from industries and fossil fuels to generate bio-renewables and other value-added chemicals.

Enhanced Thermal Stability of Polyphosphate-Dependent Glucomannokinase by Directed Evolution

Polyphosphate-dependent glucomannokinase (PPGMK) is able to utilize inorganic polyphosphate to synthesize mannose-6-phosphate (M6P) instead of highly costly ATP. This enzyme was modified and designed by combining error-prone PCR (EP-PCR) and site-directed saturation mutagenesis. Two mutants, H92L/A138V and E119V, were screened out from the random mutation library, and we used site-specific saturation mutations to find the optimal amino acid at each site. Finally, we found the optimal combination mutant, H92K/E119R. The thermal stability of H92K/E119R increased by 5.4 times at 50 °C, and the half-life at 50 °C increased to 243 min. Moreover, the enzyme activity of H92K/E119R increased to 16.6 U/mg, and its enzyme activity is twice that of WT. We analyzed the structure of the mutant using molecular dynamics simulation. We found that the shortening of the hydrogen bond distance and the formation of salt bridges can firmly connect the α-helix and β-sheet and improve the stability of the PPGMK structure.

Carbonic anhydrase to boost CO2 sequestration: Improving carbon capture utilization and storage (CCUS)

CO₂ Capture Utilization and Storage (CCUS) is a fundamental strategy to mitigate climate change, and carbon sequestration, through absorption, can be one of the solutions to achieving this goal. In nature, carbonic anhydrase (CA) catalyzes the CO₂ hydration to bicarbonates. Targeting the development of novel biotechnological routes which can compete with traditional CO₂ absorption methods, CA utilization has presented a potential to expand as a promising catalyst for CCUS applications. Driven by this feature, the search for novel CAs as biocatalysts and the utilization of enzyme improvement techniques, such as protein engineering and immobilization methods, has resulted in suitable variants able to catalyze CO₂ absorption at relevant industrial conditions. Limitations related to enzyme recovery and recyclability are still a concern in the field, affecting cost efficiency. Under different absorption approaches, CA enhances both kinetics and CO₂ absorption yields, besides reduced energy consumption. However, efforts directed to process optimization and demonstrative plants are still limited. A recent topic with great potential for development is the CA utilization in accelerated weathering, where industrial residues could be re-purposed towards becoming carbon sequestrating agents. Furthermore, research of new solvents has identified potential candidates for integration with CA in CO₂ capture, and through techno-economic assessments, CA can be a path to increase the competitiveness of alternative CO₂ absorption systems, offering lower environmental costs. This review provides a favorable scenario combining the enzyme and CO₂ capture, with possibilities in reaching an industrial-like stage in the future.

Computational design of Lactobacillus Acidophilus α-L-rhamnosidase to increase its structural stability

α-L-rhamnosidase catalyzes hydrolysis of the terminal α-L-rhamnose from various natural rhamnoglycosides, including naringin and hesperidin, and has various applications such as debittering of citrus juices in the food industry and flavonoid derhamnosylation in the pharmaceutical industry. However, its activity is lost at high temperatures, limiting its usage. To improve Lactobacillus acidophilus α-L-rhamnosidase stability, we employed molecular dynamics (MD) to identify a highly flexible region, as evaluated by its root mean square fluctuation (RMSF) value, and computational protein design (Rosetta) to increase rigidity and favorable interactions of residues in highly flexible regions. MD results show that five regions have the highest flexibilities and were selected for design by Rosetta. Twenty-one designed mutants with the best ΔΔG at each position and ΔΔG < 0 REU were simulated at high temperature. Eight designed mutants with ΔRMSF of highly flexible regions lower than -10.0% were further simulated at the optimum temperature of the wild type. N88Q, N202V, G207D, Q209M, N211T and Y213K mutants were predicted to be more stable and could maintain their native structures better than the wild type due to increased hydrogen bond interactions of designed residues and their neighboring residues. These designed mutants are promising enzymes with high potential for stability improvement.

Characterization of a Thermostable and Surfactant-Tolerant Chondroitinase B from a Marine Bacterium Microbulbifer sp. ALW1

Chondroitinase plays an important role in structural and functional studies of chondroitin sulfate (CS). In this study, a new member of chondroitinase B of PL6 family, namely ChSase B6, was cloned from marine bacterium Microbulbifer sp. ALW1 and subjected to enzymatic and structural characterization. The recombinant ChSase B6 showed optimum activity at 40 °C and pH 8.0, with enzyme kinetic parameters of Km and Vmax against chondroitin sulfate B (CSB) to be 7.85 µg/mL and 1.21 U/mg, respectively. ChSase B6 demonstrated thermostability under 60 °C for 2 h with about 50% residual activity and good pH stability under 4.0-10.0 for 1 h with above 60% residual activity. In addition, ChSase B6 displayed excellent stability against the surfactants including Tween-20, Tween-80, Trion X-100, and CTAB. The degradation products of ChSase B6-treated CSB exhibited improved antioxidant ability as a hydroxyl radical scavenger. Structural analysis and site-directed mutagenesis suggested that the conserved residues Lys248 and Arg269 were important for the activity of ChSase B6. Characterization, structure, and molecular dynamics simulation of ChSase B6 provided a guide for further tailoring for its industrial application for chondroitin sulfate bioresource development.

Integrating dynamics into enzyme engineering

Enzyme engineering has become a widely adopted practice in research labs and industry. In parallel, the past decades have seen tremendous strides in characterizing the dynamics of proteins, using a growing array of methodologies. Importantly, links have been established between the dynamics of proteins and their function. Characterizing the dynamics of an enzyme prior to, and following, its engineering is beginning to inform on the potential of 'dynamic engineering', i.e. the rational modification of protein dynamics to alter enzyme function. Here we examine the state of knowledge at the intersection of enzyme engineering and protein dynamics, describe current challenges and highlight pioneering work in the nascent area of dynamic engineering.

Thermostabilization mechanisms in thermophilic versus mesophilic three-helix bundle proteins

The engineered three‐helix bundle, UVF, is thermostabilized entropically due to heightened, native‐state dynamics. However, it is unclear whether this thermostabilization strategy is observed in natural proteins from thermophiles. We performed all‐atom, explicit solvent molecular dynamics simulations of two three‐helix bundles from thermophilic H. butylicus (2lvsN and 2lvsC) and compared their dynamics to a mesophilic three‐helix bundle, the Engrailed homeodomain (EnHD). Like UVF, 2lvsC had heightened native dynamics, which it maintained without unfolding at 100°C. Shortening and rigidification of loops in 2lvsN and 2lvsC and increased surface hydrogen bonds in 2lvsN were observed, as is common in thermophilic proteins. A buried disulfide and salt bridge in 2lvsN and 2lvsC, respectively, provided some stabilization, and addition of a homologous disulfide bond in EnHD slowed unfolding. The transferability and commonality of stabilization strategies among members of the three‐helix bundle fold suggest that these strategies may be general and deployable in designing thermostable proteins.

Carbonic anhydrase for CO2 capture, conversion and utilization

Carbonic anhydrase (CA) enzymes, catalyzing the CO₂ hydration at a high turnover number, can be employed in expediting CO₂ capture, conversion and utilization to aid in carbon neutrality. Despite extensive research over the last decade, there remain challenges in CA-related technologies due to poor stability and suboptimal use of CAs. Herein, we discuss recent advances in CA stabilization by protein engineering and enzyme immobilization, and shed light on state-of-the-art of in vitro and in vivo CA-mediated CO₂ conversion for improved production of value-added chemicals using CO₂ as a feedstock.

In Silico Prediction Methods for Site-Saturation Mutagenesis

Directed enzyme evolution has proven to be a powerful means to endow biocatalysts with novel catalytic repertoires. Apart from completely random gene mutagenesis, site-directed or site-saturation mutagenesis requires a semi-rational selection of the amino acid positions or the substituted residues, which can dramatically reduce the screening efforts in protein engineering. To this end, in silico prediction methods play a pivotal role in targeting site-saturation mutagenesis. In this chapter, we provide two distinct computational methods, (a) conformational dynamics-guided design and (b) protein-ligand interaction fingerprinting analysis, to identify specific positions for site-saturation mutagenesis toward manipulating substrate specificity/stereoselectivity of an alcohol dehydrogenase, and improving activity of a carboxylic acid reductase, respectively.

Unlocking the Stereoselectivity and Substrate Acceptance of Enzymes: Proline-Induced Loop Engineering Test

Protein stability and evolvability influence each other. Although protein dynamics play essential roles in various catalytically important properties, their high flexibility and diversity makes it difficult to incorporate such properties into rational engineering. Therefore, how to unlock the potential evolvability in a user-friendly rational design process remains a challenge. In this endeavor, we describe a method for engineering an enantioselective alcohol dehydrogenase. It enables synthetically important substrate acceptance for 4-chlorophenyl pyridine-2-yl ketone, and perfect stereocontrol of both (S)- and (R)-configured products. Thermodynamic analysis unveiled the subtle interaction between enzyme stability and evolvability, while computational studies provided insights into the origin of selectivity and substrate recognition. Preparative-scale synthesis of the (S)-product (73 % yield; >99 % ee) was performed on a gram-scale. This proof-of-principle study demonstrates that interfaced proline residues can be rationally engineered to unlock evolvability and thus provide access to new biocatalysts with highly improved catalytic performance.

Semirational engineering of an aldo-keto reductase KmAKR for overcoming trade-offs between catalytic activity and thermostability

Enzyme engineering usually generates trade-offs between activity, stability, and selectivity. Herein, we report semirational engineering of an aldo-keto reductase (AKR) KmAKR for simultaneously enhancing its thermostability and catalytic activity. Previously, we constructed KmAKR_M9 (W297H/Y296W/K29H/Y28A/T63M/A30P/T302S/N109K/S196C), which showed outstanding activity towards t-butyl 6-chloro-(3R,5S)-dihydroxyhexanoate ((3R,5S)-CDHH), and t-butyl 6-cyano-(3R,5R)-dihydroxyhexanoate, the key chiral building blocks of rosuvastatin and atorvastatin. Under the guidance of computer-aided design including consensus residues analysis and molecular dynamics (MD) simulations, K164, S182, S232, and Q266 were dug out for their thermostability conferring roles, generating the "best" mutant KmAKR_M13 (W297H/Y296W/K29H/Y28A/T63M/A30P/T302S/N109K/S196C/K164E/S232A/S182H/Q266D). The T_m and T₅₀ ¹⁵ values of KmAKR_M13 were 10.4 and 6.1°C higher than that of KmAKR_M9 , respectively. Moreover, it displayed a significantly elevated organic solvent tolerance over KmAKR_M9 . Structural analysis indicated that stabilization of the α-helixes mainly contributed to thermostability enhancement. Under the optimized conditions, KmAKR_M13 completely asymmetrically reduced 400 g/l t-butyl 6-chloro-(5S)-hydroxy-3-oxohexanoate ((5S)-CHOH) in 8.0 h at a high substrate to catalyst ratio (S/C) of 106.7 g/g, giving diastereomerically pure (3R,5S)-CDHH (>99.5% d.e._P ) with a space-time yield (STY) of 449.2 g/l·d.

In Silico Investigation of Potential Applications of Gamma Carbonic Anhydrases as Catalysts of CO2 Biomineralization Processes: A Visit to the Thermophilic Bacteria Persephonella hydrogeniphila, Persephonella marina, Thermosulfidibacter takaii, and Thermus thermophilus

Carbonic anhydrases (CAs) have been identified as ideal catalysts for CO₂ sequestration. Here, we report the sequence and structural analyses as well as the molecular dynamics (MD) simulations of four γ-CAs from thermophilic bacteria. Three of these, Persephonella marina, Persephonella hydrogeniphila, and Thermosulfidibacter takaii originate from hydrothermal vents and one, Thermus thermophilus HB8, from hot springs. Protein sequences were retrieved and aligned with previously characterized γ-CAs, revealing differences in the catalytic pocket residues. Further analysis of the structures following homology modeling revealed a hydrophobic patch in the catalytic pocket, presumed important for CO₂ binding. Monitoring of proton shuttling residue His69 (P. marina γ-CA numbering) during MD simulations of P. hydrogeniphila and P. marina’s γ-CAs (γ-PhCA and γ-PmCA), showed a different behavior to that observed in the γ-CA of Escherichia coli, which periodically coordinates Zn²⁺. This work also involved the search for hotspot residues that contribute to interface stability. Some of these residues were further identified as key in protein communication via betweenness centrality metric of dynamic residue network analysis. T. takaii’s γ-CA showed marginally lower thermostability compared to the other three γ-CA proteins with an increase in conformations visited at high temperatures being observed. Hydrogen bond analysis revealed important interactions, some unique and others common in all γ-CAs, which contribute to interface formation and thermostability. The seemingly thermostable γ-CA from T. thermophilus strangely showed increased unsynchronized residue motions at 423 K. γ-PhCA and γ-PmCA were, however, preliminarily considered suitable as prospective thermostable CO₂ sequestration agents.

Engineering stable carbonic anhydrases for CO2 capture: a critical review

Targeted inhibition of misregulated protein-protein interactions (PPIs) has been a promising area of investigation in drug discovery and development for human diseases. However, many constraints remain, including shallow binding surfaces and dynamic conformation changes upon interaction. A particularly challenging aspect is the undesirable off-target effects caused by inherent structural similarity among the protein families. To tackle this problem, phage display has been used to engineer PPIs for high-specificity binders with improved binding affinity and greatly reduced undesirable interactions with closely related proteins. Although general steps of phage display are standardized, library design is highly variable depending on experimental contexts. Here in this review, we examined recent advances in the structure-based combinatorial library design and the advantages and limitations of different approaches. The strategies described here can be explored for other protein-protein interactions and aid in designing new libraries or improving on previous libraries.

Molecular dynamics simulations identify the regions of compromised thermostability in SazCA

The present study examined the structure and dynamics of the most active and thermostable carbonic anhydrase, SazCA, probed using molecular dynamics simulations. The molecular system was described by widely used biological force-fields (AMBER, CHARMM22, CHARMM36, and OPLS-AA) in conjunction with TIP3P water model. The comparison of molecular dynamics simulation results suggested AMBER to be a suitable choice to describe the structure and dynamics of SazCA. In addition to this, we also addressed the effect of temperature on the stability of SazCA. We performed molecular dynamics simulations at 313, 333, 353, 373, and 393 K to study the relationship between thermostability and flexibility in SazCA. The amino acid residues VAL98, ASN99, GLY100, LYS101, GLU145, and HIS207 were identified as the most flexible residues from root-mean-square fluctuations. The salt bridge analysis showed that ion-pairs ASP113-LYS81, ASP115-LYS81, ASP115-LYS114, GLU144-LYS143, and GLU144-LYS206, were responsible for the compromised thermal stability of SazCA.

High-level production in a plant system of a thermostable carbonic anhydrase and its immobilization on microcrystalline cellulose beads for CO2 capture

KEY MESSAGE

Plant-produced SazCA and its application to CO2 capture. Technologies that rely on chemical absorption or physical adsorption have been developed to capture CO2 from industrial flue gases and sequester it at storage sites. Carbonic anhydrases (CAs), metalloenzymes, that catalyze the reversible hydration of CO2 have recently received attention as biocatalysts in the capture of CO2 from flue gases, but their cost presents a major obstacle for use at an industrial scale. This cost, however, can be reduced either by producing a long-lasting enzyme suitable for CO2 capture or by lowering production costs. High-level expression, easy purification, and immobilization of CAs from Sulfurihydrogenibium azorense (SazCA) were investigated in a plant system. Fusion of the 60-amino acid-long ectodomain (M-domain) of the human receptor-type tyrosine-protein phosphatase C increased the levels of SazCA accumulation. Fusion of the cellulose-binding module (CBM3) from Clostridium thermocellum resulted in tight binding of recombinant protein to microcrystalline cellulose beads, enabling easy purification. The chimeric fusion protein, BMC-SazCA, which consisted of SazCA with the M and CBM3 domains, was expressed in tobacco (Nicotiana benthamiana), giving a recombinant protein yield in leaf extracts of 350 mg/kg fresh weight. BMC-SazCA produced in planta was active in the presence of various chemicals used in CO2 capture. Immobilization of BMC-SazCA on the surface of microcrystalline cellulose beads extended its heat stability, allowing its reuse in multiple rounds of the CO2 hydration reaction. These results suggest that production of SazCA in plants has great potential for CA-based CO2 sequestration and mineralization.

Enhanced Thermostability and Enzymatic Activity of Cel6A Variants from Thermobifida fusca by Empirical Domain Engineering (Short Title: Domain Engineering of Cel6A)

Cellulases are a set of lignocellulolytic enzymes, capable of producing eco-friendly low-cost renewable bioethanol. However, low stability and hydrolytic activity limit their wide-scale applicability at the industrial scale. In this work, we report the domain engineering of endoglucanase (Cel6A) of Thermobifida fusca to improve their catalytic activity and thermal stability. Later, enzymatic activity and thermostability of the most efficient variant named as Cel6A.CBC was analyzed by molecular dynamics simulations. This variant demonstrated profound activity against soluble and insoluble cellulosic substrates like filter paper, alkali-treated bagasse, regenerated amorphous cellulose (RAC), and bacterial microcrystalline cellulose. The variant Cel6A.CBC showed the highest catalysis of carboxymethyl cellulose (CMC) and other related insoluble substrates at a pH of 6.0 and a temperature of 60 °C. Furthermore, a sound rationale was observed between experimental findings and molecular modeling of Cel6A.CBC which revealed thermostability of Cel6A.CBC at 26.85, 60.85, and 74.85 °C as well as structural flexibility at 126.85 °C. Therefore, a thermostable derivative of Cel6A engineered in the present work has enhanced biological performance and can be a useful construct for the mass production of bioethanol from plant biomass.

Molecular Dynamics Analysis of a Rationally Designed Aldehyde Dehydrogenase Gives Insights into Improved Activity for the Non-Native Cofactor NAD

The aldehyde dehydrogenase from Thermoplasma acidophilum was previously implemented as a key enzyme in a synthetic cell-free reaction cascade for the production of alcohols. In order to engineer the enzyme's cofactor specificity from NADP⁺ to NAD⁺, we identified selectivity-determining residues with the CSR-SALAD tool and investigated further positions based on the crystal structure. Stepwise combination of the initially discovered six point mutations allowed us to monitor the cross effects of each mutation, resulting in a final variant with reduced K_M for the non-native cofactor NAD⁺ (from 18 to 0.6 mM) and an increased activity for the desired substrate d-glyceraldehyde (from 0.4 to 1.5 U/mg). Saturation mutagenesis of the residues at the entrance of the substrate pocket could eliminate substrate inhibition. Molecular dynamics simulations showed a significant gain of flexibility at the cofactor binding site for the final variant. The concomitant increase in stability against isobutanol and only a minor reduction in its temperature stability render the final variant a promising candidate for future optimization of our synthetic cell-free enzymatic cascade.

Rational molecular design for improving digestive enzyme resistance of beta-glucosidase from Trichoderma viride based on inhibition of bound state formation

Beta-glucosidase (BGL1) is widely used in animal feed industries. However, degradation caused by digestive enzymes in the intestine hampers its application. Improving the resistance of feed enzymes against proteases is crucial in livestock farming. To improve the resistance of beta-glucosidase against pepsin and trypsin, a rational molecular design based on the inhibition of bound-state formation and secondary design was developed. The strategy includes: (1) prediction of the interaction surface of the pepsin-BGL1 complex structure, (2) prediction of key amino acids affecting the formation of the complex, (3) optimization of pepsin-resistant mutants by structural evaluation, (4) secondary molecular design based on pepsin-resistant mutants, and optimization of pepsin and trypsin-resistant mutants. Two BGL1 protein mutants (BGL1^Q627C and BGL1^{Q627C/R543H/R646W}) were constructed, and then mutated and wild-type BGL1s were expressed in Pichia pastoris. The half-life of BGL1^Q627C and BGL1^{Q627C/R543H/R646W} were 1.36 and 1.51 times that of the wild type upon pepsin exposure, respectively. For trypsin resistance, the half-life were 0.93 and 1.53 times that of the wild type, respectively. Compare to those of the wild type, most of the basic enzymatic properties of both mutants were not significantly changed except for increased Michaelis constants. The rational design method can be used as a guide for modifying other feed enzymes.

Characterization and High-Level Periplasmic Expression of Thermostable α-Carbonic Anhydrase from Thermosulfurimonas Dismutans in Escherichia Coli for CO2 Capture and Utilization

Carbonic anhydrase (CA) is a diffusion-controlled enzyme that rapidly catalyzes carbon dioxide (CO₂) hydration. CA has been considered as a powerful and green catalyst for bioinspired CO₂ capture and utilization (CCU). For successful industrial applications, it is necessary to expand the pool of thermostable CAs to meet the stability requirement under various operational conditions. In addition, high-level expression of thermostable CA is desirable for the economical production of the enzyme. In this study, a thermostable CA (tdCA) of Thermosulfurimonas dismutans isolated from a deep-sea hydrothermal vent was expressed in Escherichia coli and characterized in terms of expression level, solubility, activity and stability. tdCA showed higher solubility, activity, and stability compared to those of CA from Thermovibrio ammonificans, one of the most thermostable CAs, under low-salt aqueous conditions. tdCA was engineered for high-level expression by the introduction of a point mutation and periplasmic expression via the Sec-dependent pathway. The combined strategy resulted in a variant showing at least an 8.3-fold higher expression level compared to that of wild-type tdCA. The E. coli cells with the periplasmic tdCA variant were also investigated as an ultra-efficient whole-cell biocatalyst. The engineered bacterium displayed an 11.9-fold higher activity compared to that of the recently reported system with a halophilic CA. Collectively these results demonstrate that the highly expressed periplasmic tdCA variant, either in an isolated form or within a whole-cell platform, is a promising biocatalyst with high activity and stability for CCU applications.