Rational stabilization of enzymes by computational redesign of surface charge-charge interactions

Enhanced catalytic performance and pH stability of Streptomyces Laccase Y230R and its degradation of malachite green

The escalating threat of malachite green (MG) pollution poses significant risks to ecosystems. Saturation mutation targeting Tyr230 of small laccase (SLAC) from Streptomyces coelicolor yielded Y230R, exhibiting a remarkable 104 % increase in specific activity. Notably, this mutation achieved dual enhancements in both activity and pH stability. Molecular dynamics simulation revealed higher structural stability of Y230R compared to wild-type (WT) across varying pH levels. The increased count of hydrogen bonds in Y230R compared to WT may be contribute to its stability. Y230R demonstrated superior catalytic efficiency (67.0 %) in MG decolorization, maintaining over 90 % activity after 30 min incubation in MG solution (500 mg/L), highlighting enhanced tolerance compared to WT. Molecular docking analysis attributed the differential catalytic effects on MG and ABTS to structural disparities and hydrogen bonding. Y230R stands as a promising composite mutant for future laccase engineering and industrial applications.

Computer-aided rational design strategy based on protein surface charge to improve the thermal stability of a novel esterase from Geobacillus jurassicus

OBJECTIVES

Although Geobacillus are significant thermophilic bacteria source, there are no reports of thermostable esterase gene in Geobacillus jurassicus or rational design strategies to increase the thermal stability of esterases.

RESULTS

Gene gju768 showed a highest similarity of 15.20% to esterases from Geobacillus sp. with detail enzymatic properties. Using a combination of Gibbs Unfolding Free Energy (∆∆G) calculator and the distance from the mutation site to the catalytic site (DsCα-Cα) to screen suitable mutation sites with elimination of negative surface charge, the mutants (D24N, E221Q, and E253Q) displayed stable mutants with higher thermal stability than the wild-type (WT). Mutant E253Q exhibited the best thermal stability, with a half-life (T1/2) at 65 °C of 32.4 min, which was 1.8-fold of the WT (17.9 min).

CONCLUSION

Cloning of gene gju768 and rational design based on surface charge engineering contributed to the identification of thermostable esterase from Geobacillus sp. and the exploration of evolutionary strategies for thermal stability.

Functional expression, purification, biochemical and biophysical characterizations, and molecular dynamics simulation of a histidine acid phosphatase from Saccharomyces cerevisiae

A histidine acid phosphatase (HAP) (PhySc) with 99.50% protein sequence similarity with PHO5 from Saccharomyces cerevisiae was expressed functionally with the molecular mass of ∼110 kDa through co-expression along with the set of molecular chaperones dnaK, dnaJ, GroESL. The purified HAP illustrated the optimum activity of 28.75 ± 0.39 U/mg at pH 5.5 and 40 ˚C. The Km and Kcat values towards calcium phytate were 0.608 ± 0.09 mM and 650.89 ± 3.6 s- 1. The half-lives (T1/2) at 55 and 60 ˚C were 2.75 min and 55 s, respectively. The circular dichroism (CD) demonstrated that PhySc includes 30.5, 28.1, 21.3, and 20.1% of random coils, α-Helix, β-Turns, and β-Sheet, respectively. The Tm recorded by CD for PhySc was 56.5 ± 0.34˚C. The molecular docking illustrated that His59 and Asp322 act as catalytic residues in the PhySc. MD simulation showed that PhySc at 40 ˚C has higher structural stability over those of the temperatures 60 and 80 ˚C that support the thermodynamic in vitro investigations. Secondary structure content results obtained from MD simulation indicated that PhySc consists of 34.03, 33.09, 17.5, 12.31, and 3.05% of coil, helix, turn, sheet, and helix310, respectively, which is almost consistent with the experimental results.

Computation-Based Design of Salt Bridges in PETase for Enhanced Thermostability and Performance for PET Degradation

Polyethylene terephthalate (PET) is one of the most widely used plastics, and the accumulation of PET poses a great threat to the environment. IsPETase can degrade PET rapidly at moderate temperatures, but its application is greatly limited by the low stability. Herein, molecular dynamics (MD) simulations combined with a sequence alignment strategy were adopted to introduce salt bridges into the flexible region of IsPETase to improve its thermal stability. In the designed variants, the Tm values of IsPETaseI168R/S188D and IsPETaseI168R/S188E were 7.4 and 8.7 °C higher than that of the wild type, respectively. The release of products degraded by IsPETaseI168R/S188E was 4.3 times that of the wild type. Tertiary structure characterization demonstrated that the structure of the variants IsPETaseI168R/S188D and IsPETaseI168R/S188E became more compact. Extensive MD simulations verified that a stable salt bridge was formed between the residue R168 and D186 in IsPETaseI168R/S188D , while in IsPETaseI168R/S188E an R168-D186-E188 salt bridge network was observed. These results confirmed that the proposed computation-based salt bridge design strategy could efficiently generate variants with enhanced thermal stability for the long-term degradation of PET, which would be helpful for the design of enzymes with improved stability.

Computational aided design of a halotolerant CMP kinase for enzymatic synthesis of cytidine triphosphate

The current biocatalytic method of industrial Cytidine triphosphate (CTP) production suffers from reaction rate loss. It is caused by gradually increasing acetate salt concentration, which inhibits enzyme activities and decreases the final yield. This work gave a possible solution to this problem through computational aided design of CMP kinase (CMPK), an enzyme in the CTP production system, to increase its stability in solution with high acetate salt concentration. Enlightened by the features of natural halophilic enzymes, the basic and neutral surface residues were replaced with acidic amino acids. This protein design strategy effectively increased the activity of CMPK in the working condition (acetate concentration over 1200 mM). The halotolerant CMPK was applied in fed-batch production of CTP. The maximum titer was 201.4 ± 1.6 mM, and the productivity was 12.6 mM L^-1 h^-1, increased 26.4% and 27.8% from the process using wild-type CMPK, respectively.

Charge reversal mutations in mesophilic-thermophilic orthologous protein pairs and their role in enhancing coulombic interaction energy

Proteins from thermophilic organisms are a matter of immense interest for decades because of its application in fields like de-novo protein design, thermostable variants of biocatalysts etc. Previous studies have found several sequence and structural adaptations related to thermal stability, while charge reversal study remains ignored. Here we address whether charge reversal mutations naturally occur in mesophilic-thermophilic/hyperthermophilic orthologous proteins. Do they contribute to thermal stability? Our systematic study on 1550 mesophilic-thermophilic/hyperthermophilic orthologous protein pairs with remarkable structural and topological similarity, shows gain in coulombic interaction energy in thermophilic/hyperthermophilic proteins at short range associated with partially exposed and buried charge reversal mutations, which may enhance thermostability. Our findings call forth its application in future protein engineering studies. Communicated by Ramaswamy H. Sarma.

Protein engineering via Bayesian optimization-guided evolutionary algorithm and robotic experiments

Directed protein evolution applies repeated rounds of genetic mutagenesis and phenotypic screening and is often limited by experimental throughput. Through in silico prioritization of mutant sequences, machine learning has been applied to reduce wet lab burden to a level practical for human researchers. On the other hand, robotics permits large batches and rapid iterations for protein engineering cycles, but such capacities have not been well exploited in existing machine learning-assisted directed evolution approaches. Here, we report a scalable and batched method, Bayesian Optimization-guided EVOlutionary (BO-EVO) algorithm, to guide multiple rounds of robotic experiments to explore protein fitness landscapes of combinatorial mutagenesis libraries. We first examined various design specifications based on an empirical landscape of protein G domain B1. Then, BO-EVO was successfully generalized to another empirical landscape of an Escherichia coli kinase PhoQ, as well as simulated NK landscapes with up to moderate epistasis. This approach was then applied to guide robotic library creation and screening to engineer enzyme specificity of RhlA, a key biosynthetic enzyme for rhamnolipid biosurfactants. A 4.8-fold improvement in producing a target rhamnolipid congener was achieved after examining less than 1% of all possible mutants after four iterations. Overall, BO-EVO proves to be an efficient and general approach to guide combinatorial protein engineering without prior knowledge.

Thermal stabilization of diverse biologics using reversible hydrogels

Improving the thermal stability of biologics, including vaccines, is critical to reduce the economic costs and health risks associated with the cold chain. Here, we designed a versatile, safe, and easy-to-use reversible PEG-based hydrogel platform formed via dynamic covalent boronic ester cross-linking for the encapsulation, stabilization, and on-demand release of biologics. Using these reversible hydrogels, we thermally stabilized a wide range of biologics up to 65°C, including model enzymes, heat-sensitive clinical diagnostic enzymes (DNA gyrase and topoisomerase I), protein-based vaccines (H5N1 hemagglutinin), and whole viruses (adenovirus type 5). Our data support a generalized protection mechanism for the thermal stabilization of diverse biologics using direct encapsulation in reversible hydrogels. Furthermore, preliminary toxicology data suggest that the components of our hydrogel are safe for in vivo use. Our reversible hydrogel platform offers a simple material solution to mitigate the costs and risks associated with reliance on a continuous cold chain for biologic transport and storage.

A Novel NADP(H)-Dependent 7alpha-HSDH: Discovery and Construction of Substrate Selectivity Mutant by C-Terminal Truncation

7α-Hydroxysteroid dehydrogenase (7α-HSDH) plays an important role in the biosynthesis of tauroursodeoxycholic acid (TUDCA) using complex substrate chicken bile powder as raw material. However, chicken bile powder contains 4.74% taurocholic acid (TCA), and a new by-product tauroursocholic acid (TUCA) will be produced, having the risk of causing colorectal cancer. Here, we obtained a novel NADP(H)-dependent 7α-HSDH with good thermostability from Ursus thibetanus gut microbiota (named St-2-2). St-2-2 could catalyze taurochenodeoxycholic acid (TCDCA) and TCA with the catalytic activity of 128.13 and 269.39 U/mg, respectively. Interestingly, by a structure-based C-terminal truncation strategy, St-2-2△C10 only remained catalytic activity on TCDCA (14.19 U/mg) and had no activity on TCA. As a result, it can selectively catalyze TCDCA in waste chicken bile powder. MD simulation and structural analysis indicated that enhanced surface hydrophilicity and improved C-terminal rigidity affected the entry and exit of substrates. Hydrogen bond interactions between different subunits and interaction changes in Phe249 of the C-terminal loop inverted the substrate catalytic activity. This is the first report on substrate selectivity of 7α-HSDH by C-terminal truncation strategy and it can be extended to other 7α-HSDHs (J-1-1, S1-a-1).

Thermostability engineering of industrial enzymes through structure modification

Thermostability is an essential requirement of enzymes in the industrial processes to catalyze the reactions at high temperatures; thus, enzyme engineering through directed evolution, semi-rational design and rational design are commonly employed to construct desired thermostable mutants. Several strategies are implemented to fulfill enzymes' thermostability demand including decreasing the entropy of the unfolded state through substitutions Gly → Xxx or Xxx → Pro, hydrogen bond, salt bridge, introducing two different simultaneous interactions through single mutant, hydrophobic interaction, filling the hydrophobic cavity core, decreasing surface hydrophobicity, truncating loop, aromatic-aromatic interaction and introducing positively charged residues to enzyme surface. In the current review, horizons about compatibility between secondary structures and substitutions at preferable structural positions to generate the most desirable thermostability in industrial enzymes are broadened. KEY POINTS: • Protein engineering is a powerful tool for generating thermostable industrial enzymes. • Directed evolution and rational design are practical approaches in enzyme engineering. • Substitutions in preferable structural positions can increase thermostability.

Tuning Enzyme Thermostability via Computationally Guided Covalent Stapling and Structural Basis of Enhanced Stabilization

Enhancing the thermostability of enzymes without impacting their catalytic function represents an important yet challenging goal in protein engineering and biocatalysis. We recently introduced a novel method for enzyme thermostabilization that relies on the computationally guided installation of genetically encoded thioether "staples" into a protein via cysteine alkylation with the noncanonical amino acid O-2-bromoethyl tyrosine (O2beY). Here, we demonstrate the functionality of an expanded set of electrophilic amino acids featuring chloroacetamido, acrylamido, and vinylsulfonamido side-chain groups for protein stapling using this strategy. Using a myoglobin-based cyclopropanase as a model enzyme, our studies show that covalent stapling with p-chloroacetamido-phenylalanine (pCaaF) provides higher stapling efficiency and enhanced stability (thermodynamic and kinetic) compared to the other stapled variants and the parent protein. Interestingly, molecular simulations of conformational flexibility of the cross-links show that the pCaaF staple allows fewer energetically feasible conformers than the other staples, and this property may be a broader indicator of stability enhancement. Using this strategy, pCaaF-stapled variants with significantly enhanced stability against thermal denaturation (ΔTm' = +27 °C) and temperature-induced heme loss (ΔT50 = +30 °C) were obtained while maintaining high levels of catalytic activity and stereoselectivity. Crystallographic analyses of singly and doubly stapled variants provide key insights into the structural basis for stabilization, which includes both direct interactions of the staples with protein residues and indirect interactions through adjacent residues involved in heme binding. This work expands the toolbox of protein stapling strategies available for protein stabilization.

Metagenomic mining and structure-function studies of a hyper-thermostable cellobiohydrolase from hot spring sediment

Enzymatic breakdown is an attractive cellulose utilisation method with a low environmental load. Its high temperature operation could promote saccharification and lower contamination risk. Here we report a hyper-thermostable cellobiohydrolase (CBH), named HmCel6A and its variant HmCel6A-3SNP that were isolated metagenomically from hot spring sediments and expressed in Escherichia coli. They are classified into glycoside hydrolases family 6 (GH6). HmCel6A-3SNP had three amino acid replacements to HmCel6A (P88S/L230F/F414S) and the optimum temperature at 95 °C, while HmCel6A did it at 75 °C. Crystal structure showed conserved features among GH6, a (β/α)₈-barrel core and catalytic residues, and resembles TfCel6B, a bacterial CBH II of Thermobifida fusca, that had optimum temperature at 60 °C. From structure-function studies, we discuss unique structural features that allow the enzyme to reach its high thermostability level, such as abundance of hydrophobic and charge-charge interactions, characteristic metal bindings and disulphide bonds. Moreover, structure and surface plasmon resonance analysis with oligosaccharides suggested that the contribution of an additional tryptophan located at the tunnel entrance could aid in substrate recognition and thermostability. These results may help to design efficient enzymes and saccharification methods for cellulose working at high temperatures.

Bacteria from hot springs are known for highly thermostable enzymes, which may have industrial potential. Here, a unique thermostable cellobiohydrolase is reported that can breakdown cellulose at temperature up to 95 degrees Celsius.

Enhanced thermal and alkaline stability of L-lysine decarboxylase CadA by combining directed evolution and computation-guided virtual screening

As an important monomer for bio-based nylons PA5X, cadaverine is mainly produced by enzymatic decarboxylation of L-lysine. A key issue with this process is the instability of L-lysine decarboxylase (CadA) during the reaction due to the dissociation of CadA subunits with the accumulation of alkaline cadaverine. In this work, we attempted to improve the thermal and alkaline stability of CadA by combining directed evolution and computation-guided virtual screening. Interestingly, site 477 residue located at the protein surface and not the decamer interface was found as a hotspot in directed evolution. By combinatorial mutagenesis of the positive mutations obtained by directed evolution and virtual screening with the previously reported T88S mutation, K477R/E445Q/T88S/F102V was generated as the best mutant, delivering 37% improvement of cadaverine yield at 50 ºC and pH 8.0. Molecular dynamics simulations suggested the improved rigidity of regional structures, increased number of salt bridges, and enhancement of hydrogen bonds at the multimeric interface as possible origins of the improved stability of the mutant. Using this four-point mutant, 160.7 g/L of cadaverine was produced from 2.0 M Lysine hydrochloride at 50 °C without pH regulation, with a conversion of 78.5%, whereas the wild type produced 143.7 g/L cadaverine, corresponding to 70% conversion. This work shows the combination of directed evolution and virtual screening as an efficient protein engineering strategy.

Enhancing the Catalytic Performance of Candida antarctica Lipase B by Chemical Modification With Alkylated Betaine Ionic Liquids

Various betaine ionic liquids composed of different chain lengths and different anions were designed and synthesized to modify Candida antarctica lipase B (CALB). The results showed that the catalytic activity of all modified lipases improved under different temperature and pH conditions, while also exhibiting enhanced thermostability and tolerance to organic solvents. With an increase in ionic liquid chain length, the modification effect was greater. Overall, CALB modified by [BetaineC₁₆][H₂PO₄] performed best, with the modified CALB enzyme activity increased 3-fold, thermal stability increased 1.5-fold when stored at 70°C for 30 min, with tolerance increased 2.9-fold in 50% DMSO and 2.3-fold in 30% mercaptoethanol. Fluorescence and circular dichroism (CD) spectroscopic analysis showed that the introduction of an ionic liquid caused changes in the microenvironment surrounding some fluorescent groups and the secondary structure of the CALB enzyme protein. In order to establish the enzyme activity and stability change mechanisms of the modified CALB, the structures of CALB modified with [BetaineC₄][Cl] and [BetaineC₁₆][Cl] were constructed, while the reaction mechanisms were studied by molecular dynamics simulations. Results showed that the root mean square deviation (RMSD) and total energy of modified CALB were less than those of native CALB, indicating that modified CALB has a more stable structure. Root mean square fluctuation (RMSF) calculations showed that the rigidity of modified CALB was enhanced. Solvent accessibility area (SASA) calculations exhibited that both the hydrophilicity and hydrophobicity of the modified enzyme-proteins were improved. The increase in radial distribution function (RDF) of water molecules confirmed that the number of water molecules around the active sites also increased. Therefore, modified CALB has enhanced structural stability and higher hydrolytic activity.

Computer-Aided Targeted Mutagenesis of Thermoclostridium caenicola d-Allulose 3-Epimerase for Improved Thermostability

d-Allulose 3-epimerase (DAEase) is a key enzyme in d-allulose bioproduction. DAEase from Thermoclostridium caenicola suffers from poor thermostability, hampering its large-scale applications in industry. In this study, mutants A70P, G107P, F155Y, and D162T with increased melting point temperature (T_m) were obtained by targeted mutagenesis based on the calculation of the free energy of folding. The optimal single-point mutant G107P showed 11.08 h, 5, and 5.70 °C increases in the values of half-life (t_1/2) at 60 °C, the optimum temperature (T_opt), and T_m, respectively. Beneficial mutations were combined by ordered recombination mutagenesis, and the combinational mutant Var3 (G107P/F155Y/D162T/A70P) was generated with ΔT_opt of 10 °C and ΔT_m of 12.25 °C. Its t_1/2 value at 65 °C was more than 140 times higher than that of the wild-type enzyme. Molecular dynamics simulations and homology modeling analysis indicated that the enhanced overall rigidity, increased hydrogen bonds between subunits, and redistributed surface electrostatic charges might be responsible for the improved thermostability of the mutant Var3.

Enhancing thermostabilization of a newly discovered α-amylase from Bacillus cereus GL96 by combining computer-aided directed evolution and site-directed mutagenesis

This study has described the characterization of a new a-amylase from the recently isolated Bacillus cereus GL96. Subsequently, an in-silico approach was taken into account to redesign the enzyme to meet higher thermal stability. Finally, the engineered enzyme was constructed experimentally using side-directed mutagenesis (SDM) and characterized accordingly. The enzyme was stable over pH 4-11, with the highest activity at 9.5. The temperature profile of the wild-type enzyme showed optimum activity at 50 °C plus 40% of stability at temperatures up to 70 °C. The in-silico result was indicated D162W, D162R, and D162K as the three stabilizing mutations. Among them, D162K showed better results, especially in the molecular dynamics simulation, and therefore, it was constructed by SDM. This variant was shown 5 °C higher optimum temperature (55 °C) with increasing activity than the native enzyme. In addition, it was significantly more stable than the native form. For example, while the latter almost wholly lost its function at a temperature above 70 °C, the D162K can retain more than 40% of its initial activity up to 80 °C. Considering the promising properties that the mutant enzyme showed, it can be considered for further investigation to meet the industrial requirement completely.

Coarse-Grained Simulations of Protein Folding: Bridging Theory and Experiments

Computational coarse-grained models play a fundamental role as a research tool in protein folding, and they are important in bridging theory and experiments. Folding mechanisms are generally discussed using the energy landscape framework, which is well mapped within a class of simplified structure-based models. In this chapter, simplified computer models are discussed with special focus on structure-based ones.

Quantifying Protein Electrostatic Interactions in Cells by Nuclear Magnetic Resonance Spectroscopy

Most proteins perform their functions in cells. How the cellular environment modulates protein interactions is an important question. In this work, electrostatic interactions between protein charges were studied using in-cell nuclear magnetic resonance (NMR) spectroscopy. A total of eight charge pairs were introduced in protein GB3. Compared to the charge pair electrostatic interactions in a buffer, five charge pairs in cells displayed no apparent changes whereas three pairs had the interactions weakened by more than 70%. Further investigation suggests that the transfer free energy is responsible for the electrostatic interaction modulation. Both the transfer free energy of the folded state and that of the unfolded state can contribute to the cellular environmental effect on protein electrostatics, although the latter is generally larger (more negative) than the former. Our work highlights the importance of direct in-cell studies of protein interactions and thus protein function.

High surfactant-tolerant β-mannanase isolated from Dynastes hercules larvae excrement, and identification of its hotspot using site-directed mutagenesis and molecular dynamics simulations

The β-mannanase from Bacillus subtilis HM7 (Man26HM7) isolated from Dynastes hercules larvae excrement was cloned and expressed in Escherichia coli. Biochemical characterization shows that optimal pH and temperature for catalysis are 6.0 and 50 °C, respectively. Man26HM7 displayed excellent surfactant stability by retaining 70% of initial activity in 1%(w/v) SDS, and more than 90% of initial activity in 1%(w/v) Triton X-100 and Tween 80. Results from amino acid sequence alignment and molecular modeling suggest residue 238 of β-mannanase as a hotspot of SDS-tolerance. Mutagenesis at the equivalent residue of another homolog, β-mannanase from Bacillus subtilis CAe24 (Man26CAe24), significantly enhanced the SDS stability of this enzyme. Comparative computational analysis, including molecular docking and molecular dynamics simulation, were then performed to compute the binding free energy of SDS to Man26HM7, Man26CAe24, and variant enzymes. The results suggest that residue 238 of Man26HM7 is involved in SDS binding to the hydrophobic surface of β-mannanase. This study provides not only the promising application of Man26HM7 in detergent and cleaning products but also valuable information for enhancing the surfactant stability of β-mannanase by enzyme engineering.

Surface residues and non-additive interactions stabilize a consensus homeodomain protein

Despite the widely reported success of consensus design in producing highly stabilized proteins, little is known about the physical mechanisms underlying this stabilization. Here we explore the potential sources of stabilization by performing a systematic analysis of the 29 substitutions that we previously found to collectively stabilize a consensus homeodomain compared to an extant homeodomain. By separately introducing groups of consensus substitutions that alter or preserve charge state, occur at varying degrees of residue burial, and occur at positions of varying degrees of conservation, we determine the extent to which these three features contribute to the consensus stability enhancement. Surprisingly, we find that the largest total contribution to stability comes from consensus substitutions on the protein surface and that the largest per-substitution contributions come from substitutions that maintain charge state. This finding suggests that although consensus proteins are often enriched in charged residues, consensus stabilization does not result primarily from interactions involving charged residues. Although consensus substitutions at strongly conserved positions also contribute disproportionately to stabilization, significant stabilization is also contributed from substitutions at weakly conserved positions. Furthermore, we find that identical consensus substitutions show larger stabilizing effects when introduced into the consensus background than when introduced into an extant homeodomain, indicating that synergistic, stabilizing interactions among the consensus residues contribute to consensus stability enhancement of the homeodomain. By measuring DNA binding affinity for the same set of variants, we find that although consensus design of the homeodomain increases both affinity and folding stability, it does so using a largely non-overlapping set of substitutions.

Improving the alkaline stability of pepsin through rational protein design using renin, an alkaline-stable aspartic protease, as a structural and functional reference

The present study sought to identify the structural determinants of aspartic protease structural stability and activity at elevated pH. Various hypotheses have been published regarding the features responsible for the unusual alkaline structural stability of renin, however, few structure-function studies have verified these claims. Using pepsin as a model system, and renin as a template for functional and structural alkaline stability, a rational re-design of pepsin was undertaken to identify residues contributing to the alkaline instability of pepsin-like aspartic proteases in regards to both structure and function. We constructed 13 mutants based on this strategy. Among them, mutants D159 L and D60A led to an increase in activity at elevated pH levels (p ≤ 0.05) and E4V and H53F were shown to retain native-like structure at elevated pH (p ≤ 0.05). Previously suggested carboxyl groups Asp¹¹, Asp¹¹⁸, and Glu¹³ were individually shown not to be responsible for the structural instability or lack of activity at neutral pH in pepsin. The importance of the β-barrel to structural stability was highlighted as the majority of the stabilizing residues identified, and 39% of the weakly conserved residues in the N-terminal lobe, were located in β-sheet strands of the barrel. The results of the present study indicate that alkaline stabilization of pepsin will require reduction of electrostatic repulsions and an improved understanding of the role of the hydrogen bonding network of the characteristic β-barrel.

Ionic Atmosphere Effect on the Absorption Spectrum of a Flavoprotein: A Reminder to Consider Solution Ions

This study utilizes the FMN-dependent NADH:quinone oxidoreductase from Pseudomonas aeruginosa PAO1 to investigate the effect of introducing an active site negative charge on the flavin absorption spectrum both in the absence and presence of a long-range electrostatic potential coming from solution ions. There were no observed changes in the flavin UV-visible spectrum when an active site tyrosine (Y277) becomes deprotonated in vitro. These results could only be reproduced computationally using average solvent electrostatic configuration (ASEC) QM/MM simulations that include both positive and negative solution ions. The same calculations performed with minimal ions to neutralize the total protein charge predicted that deprotonating Y277 would significantly alter the flavin absorption spectrum. Analyzing the distribution of solution ions indicated that the ions reorganize around the protein surface upon Y277 deprotonation to cancel the effect of the tyrosinate on the flavin absorption spectrum. Additional biochemical experiments were performed to test this hypothesis.

Surface salt bridges contribute to the extreme thermal stability of an FN3-like domain from a thermophilic bacterium

This study uses differential scanning calorimetry, X-ray crystallography, and molecular dynamics simulations to investigate the structural basis for the high thermal stability (melting temperature 97.5°C) of a FN3-like protein domain from thermophilic bacteria Thermoanaerobacter tengcongensis (FN3tt). FN3tt adopts a typical FN3 fold with a three-stranded beta sheet packing against a four-stranded beta sheet. We identified three solvent exposed arginine residues (R23, R25, and R72), which stabilize the protein through salt bridge interactions with glutamic acid residues on adjacent strands. Alanine mutation of the three arginine residues reduced melting temperature by up to 22°C. Crystal structures of the wild type (WT) and a thermally destabilized (∆Tm -19.7°C) triple mutant (R23L/R25T/R72I) were found to be nearly identical, suggesting that the destabilization is due to interactions of the arginine residues. Molecular dynamics simulations showed that the salt bridge interactions in the WT were stable and provided a dynamical explanation for the cooperativity observed between R23 and R25 based on calorimetry measurements. In addition, folding free energy changes computed using free energy perturbation molecular dynamics simulations showed high correlation with melting temperature changes. This work is another example of surface salt bridges contributing to the enhanced thermal stability of thermophilic proteins. The molecular dynamics simulation methods employed in this study may be broadly useful for in silico surface charge engineering of proteins.

Enhancement of protein thermostability by three consecutive mutations using loop-walking method and machine learning

We developed a method to improve protein thermostability, “loop-walking method”. Three consecutive positions in 12 loops of Burkholderia cepacia lipase were subjected to random mutagenesis to make 12 libraries. Screening allowed us to identify L7 as a hot-spot loop having an impact on thermostability, and the P233G/L234E/V235M mutant was found from 214 variants in the L7 library. Although a more excellent mutant might be discovered by screening all the 8000 P233X/L234X/V235X mutants, it was difficult to assay all of them. We therefore employed machine learning. Using thermostability data of the 214 mutants, a computational discrimination model was constructed to predict thermostability potentials. Among 7786 combinations ranked in silico, 20 promising candidates were selected and assayed. The P233D/L234P/V235S mutant retained 66% activity after heat treatment at 60 °C for 30 min, which was higher than those of the wild-type enzyme (5%) and the P233G/L234E/V235M mutant (35%).

Improving the Thermostability of Xylanase A from Bacillus subtilis by Combining Bioinformatics and Electrostatic Interactions Optimization

The rational improvement of the enzyme catalytic activity is one of the most significant challenges in biotechnology. Most conventional strategies used to engineer enzymes involve selecting mutations to increase their thermostability. Determining good criteria for choosing these substitutions continues to be a challenge. In this work, we combine bioinformatics, electrostatic analysis, and molecular dynamics to predict beneficial mutations that may improve the thermostability of XynA from Bacillus subtilis. First, the Tanford-Kirkwood surface accessibility method is used to characterize each ionizable residue contribution to the protein native state stability. Residues identified to be destabilizing were mutated with the corresponding residues determined by the consensus or ancestral sequences at the same locations. Five mutants (K99T/N151D, K99T, S31R, N151D, and K154A) were investigated and compared with 12 control mutants derived from experimental approaches from the literature. Molecular dynamics results show that the mutants exhibited folding temperatures in the order K99T > K99T/N151D > S31R > N151D > WT > K154A. The combined approaches employed provide an effective strategy for low-cost enzyme optimization needed for large-scale biotechnological and medical applications.

Improving the Thermostability and Activity of Transaminase From Aspergillus terreus by Charge-Charge Interaction

Transaminases that promote the amination of ketones into amines are an emerging class of biocatalysts for preparing a series of drugs and their intermediates. One of the main limitations of (R)-selective amine transaminase from Aspergillus terreus (At-ATA) is its weak thermostability, with a half-life (t _1/2) of only 6.9 min at 40°C. To improve its thermostability, four important residue sites (E133, D224, E253, and E262) located on the surface of At-ATA were identified using the enzyme thermal stability system (ETSS). Subsequently, 13 mutants (E133A, E133H, E133K, E133R, E133Q, D224A, D224H, D224K, D224R, E253A, E253H, E253K, and E262A) were constructed by site-directed mutagenesis according to the principle of turning the residues into opposite charged ones. Among them, three substitutions, E133Q, D224K, and E253A, displayed higher thermal stability than the wild-type enzyme. Molecular dynamics simulations indicated that these three mutations limited the random vibration amplitude in the two α-helix regions of 130-135 and 148-158, thereby increasing the rigidity of the protein. Compared to the wild-type, the best mutant, D224K, showed improved thermostability with a 4.23-fold increase in t _1/2 at 40°C, and 6.08°C increase in T 50 10 . Exploring the three-dimensional structure of D224K at the atomic level, three strong hydrogen bonds were added to form a special "claw structure" of the α-helix 8, and the residues located at 151-156 also stabilized the α-helix 9 by interacting with each other alternately.

Surface charged amino acid-based strategy for rational engineering of kinetic stability and specific activity of enzymes: Linking experiments with computational modeling

A rational workflow for engineering kinetically stable enzymes with good specific activity by surface charged amino acids engineering was proposed based on systematically analyzing the results of mutating 44 negatively charged surface amino acids of a thermophilic β-mannanase (ManAK). Computational data, combined with experimental results indicated that percentage side-chain solvent accessibility (PSSA), changes in Gibbs free energy of unfolding (∆∆G_mut) and root-mean-square fluctuations (RMSF) could be suitable for screening kinetically stable mutants. A combinational standard (∆∆G_mut < -0.5 kJ/mol and RMSF >0.68 Å) resulted a decrease in the proportion of destabilizing mutants to 12.5%. The perturbations of substrate affinity and specific activity caused by mutation were weakened as the shortest distance from C_α of mutated site to C_α of catalytic sites (Ds_Cα-Cα) increased. Results indicated that hotspot zones contributing to the local stability and integrity of catalytic motif at elevated temperatures might be widely distributed across spatial structure of the protein, while the mutation perturbation on enzyme specific activity demonstrated a gradually weakening trend from the catalytic core to the protein surface. These findings further our understanding of the structural-functional relationships of protein and highlight a deduced workflow to engineering industrially useful enzymes.

Electrostatic interaction optimization improves catalytic rates and thermotolerance on xylanases

Understanding the aspects that contribute to improving proteins' biochemical properties is of high relevance for protein engineering. Properties such as the catalytic rate, thermal stability, and thermal resistance are crucial for applying enzymes in the industry. Different interactions can influence those biochemical properties of an enzyme. Among them, the surface charge-charge interactions have been a target of particular attention. In this study, we employ the Tanford-Kirkwood solvent accessibility model using the Monte Carlo algorithm (TKSA-MC) to predict possible interactions that could improve stability and catalytic rate of a WT xylanase (XynA^WT) and its M6 xylanase (XynA^M6) mutant. The modeling prediction indicates that mutating from a lysine in position 99 to a glutamic acid (K99E) favors the native state stabilization in both xylanases. Our lab results showed that mutated xylanases had their thermotolerance and catalytic rate increased, which conferred higher processivity of delignified sugarcane bagasse. The TKSA-MC approach employed here is presented as an efficient computational-based design strategy that can be applied to improve the thermal resistance of enzymes with industrial and biotechnological applications.

Computational design of noncanonical amino acid-based thioether staples at N/C-terminal domains of multi-modular pullulanase for thermostabilization in enzyme catalysis

Enzyme thermostabilization is considered a critical and often obligatory step in biosynthesis, because thermostability is a significant property of enzymes that can be used to evaluate their feasibility for industrial applications. However, conventional strategies for thermostabilizing enzymes generally introduce non-covalent interactions and/or natural covalent bonds caused by natural amino acid substitutions, and the trade-off between the activity and stability of enzymes remains a challenge. Here, we developed a computationally guided strategy for constructing thioether staples by incorporating noncanonical amino acid (ncAA) into the more flexible N/C-terminal domains of the multi-modular pullulanase from Bacillus thermoleovorans (BtPul) to enhance its thermostability. First, potential thioether staples located in the N/C-terminal domains of BtPul were predicted using RosettaMatch. Next, eight variants involving stable thioether staples were precisely predicted using FoldX and Rosetta ddg_monomer. Six positive variants were obtained, of which T73(O2beY)-171C had a 157% longer half-life at 70 °C and an increase of 7.0 °C in T m, when compared with the wild-type (WT). T73(O2beY)-171C/T126F/A72R exhibited an even more improved thermostability, with a 211% increase in half-life at 70 °C and a 44% enhancement in enzyme activity compared with the WT, which was attributed to further optimization of the local interaction network. This work introduces and validates an efficient strategy for enhancing the thermostability and activity of multi-modular enzymes.

Improving the kinetic stability of a hyperthermostable β-mannanase by a rationally combined strategy

Feasible and easily accessible methods for the rational design of enzyme engineering strategies remain to be established. Thus, a new rationally combined strategy based on disulfide bond engineering and HotSpot Wizard 3.0 was proposed and experimentally demonstrated to be effective using a hyperthermostable β-mannanase. Ten of 42 mutants showed prominent enhancement of kinetic stability with 26.4%-39.9% increases in t_1/2 (75 °C) compared with the parent enzyme ManAKH. The best mutant, D273-V308, showed apparent increases in both optimal temperature (5 °C) and T₅₀ (6.8 °C), as well as advanced catalytic efficiency. The low rate of inactive mutants and the high rate of positive mutants indicated that newly introduced screening factors (distance from catalytic residues, Gibbs free energy term, molecular simulation, and visual inspections) greatly enhance the design of thermostable β-mannanase. Moreover, these findings further advance the industrial application of β-mannanase (ManAK) in food and food-related applications.

Improved thermostability of creatinase from Alcaligenes Faecalis through non-biased phylogenetic consensus-guided mutagenesis

Background

Enzymatic quantification of creatinine has become an essential method for clinical evaluation of renal function. Although creatinase (CR) is frequently used for this purpose, its poor thermostability severely limits industrial applications. Herein, we report a novel creatinase from Alcaligenes faecalis (afCR) with higher catalytic activity and lower K_M value, than currently used creatinases. Furthermore, we developed a non-biased phylogenetic consensus method to improve the thermostability of afCR.

Results

We applied a non-biased phylogenetic consensus method to identify 59 candidate consensus residues from 24 creatinase family homologs for screening afCR mutants with improved thermostability. Twenty-one amino acids of afCR were selected to mutagenesis and 11 of them exhibited improved thermostability compared to the parent enzyme (afCR-M0). Combination of single-site mutations in sequential screens resulted in a quadruple mutant D17V/T199S/L6P/T251C (M4-2) which showed ~ 1700-fold enhanced half-life at 57 °C and a 4.2 °C higher T₅₀¹⁵ than that of afCR-M0. The mutant retained catalytic activity equivalent to afCR-M0, and thus showed strong promise for application in creatinine detection. Structural homology modeling revealed a wide range of potential molecular interactions associated with individual mutations that contributed to improving afCR thermostability.

Conclusions

Results of this study clearly demonstrated that the non-biased-phylogenetic consensus design for improvement of thermostability in afCR is effective and promising in improving the thermostability of more enzymes.

Computation-aided engineering of starch-debranching pullulanase from Bacillus thermoleovorans for enhanced thermostability

Pullulanases are widely used in food, medicine, and other industries because they specifically hydrolyze α-1,6-glycosidic linkages in starch and oligosaccharides. In addition, high-temperature thermostable pullulanase has multiple advantages, including decreasing saccharification solution viscosity accompanied with enhanced mass transfer and reducing microbial contamination in starch hydrolysis. However, thermophilic pullulanase availability remains limited. Additionally, most do not meet starch-manufacturing requirements due to weak thermostability. Here, we developed a computation-aided strategy to engineer the thermophilic pullulanase from Bacillus thermoleovorans. First, three computational design predictors (FoldX, I-Mutant 3.0, and dDFIRE) were combined to predict stability changes introduced by mutations. After excluding conserved and catalytic sites, 17 mutants were identified. After further experimental verification, we confirmed six positive mutants. Among them, the G692M mutant had the highest thermostability improvement, with 3.8 °C increased T_m and 2.1-fold longer half-life than the wild type at 70 °C. We then characterized the mechanism underlying increased thermostability, such as rigidity enhancement, closer conformation, and strengthened motion correlation using root mean square fluctuation (RMSF), principal component analysis (PCA), dynamic cross-correlation map (DCCM), and free energy landscape (FEL) analysis. KEY POINTS: • A computation-aided strategy was developed to engineer pullulanase thermostability. • Seventeen mutants were identified by combining three computational design predictors. • The G692M mutant was obtained with increased T_mand half-life at 70 °C.

Identification of a chitosanase from the marine metagenome and its molecular improvement based on evolution data

Chitooligosaccharides have important application value in the fields of food and agriculture. Chitosanase can degrade chitosan to obtain chitooligosaccharides. The marine metagenome contains many genes related to the degradation of chitosan. However, it is difficult to mine valuable genes from large gene resources. This study proposes a method to screen chitosanases directly from the marine metagenome. Chitosanase gene chis1754 was identified from the metagenome and heterologously expressed in Escherichia coli. The optimal temperature and pH of CHIS1754 were 55 °C and 5.5, respectively. A mutant, CHIS1754T, with 15 single point mutations designed based on molecular evolution data was also expressed in E. coli. The results indicated that the thermal stability of CHIS1754T was significantly improved, as the T_m showed an increase of ~ 7.63 °C. Additionally, the k_cat/K_m of CHIS1754T was 4.8-fold higher than that of the wild type. This research provides new theories and foundations for the excavation, modification, and industrial application of chitosanases. KEY POINTS: A chitosanase gene, chis1754, was firstly identified from marine metagenome. A multi-site mutant was designed to improve enzyme stability and activity. The k_cat/K_mof the designed mutant was 4.8-fold higher than that of the wild type.

Approaching boiling point stability of an alcohol dehydrogenase through computationally-guided enzyme engineering

Enzyme instability is an important limitation for the investigation and application of enzymes. Therefore, methods to rapidly and effectively improve enzyme stability are highly appealing. In this study we applied a computational method (FRESCO) to guide the engineering of an alcohol dehydrogenase. Of the 177 selected mutations, 25 mutations brought about a significant increase in apparent melting temperature (ΔTm ≥ +3 °C). By combining mutations, a 10-fold mutant was generated with a Tm of 94 °C (+51 °C relative to wild type), almost reaching water's boiling point, and the highest increase with FRESCO to date. The 10-fold mutant's structure was elucidated, which enabled the identification of an activity-impairing mutation. After reverting this mutation, the enzyme showed no loss in activity compared to wild type, while displaying a Tm of 88 °C (+45 °C relative to wild type). This work demonstrates the value of enzyme stabilization through computational library design.

Robust ω-Transaminases by Computational Stabilization of the Subunit Interface

Transaminases are attractive catalysts for the production of enantiopure amines. However, the poor stability of these enzymes often limits their application in biocatalysis. Here, we used a framework for enzyme stability engineering by computational library design (FRESCO) to stabilize the homodimeric PLP fold type I ω-transaminase from Pseudomonas jessenii. A large number of surface-located point mutations and mutations predicted to stabilize the subunit interface were examined. Experimental screening revealed that 10 surface mutations out of 172 tested were indeed stabilizing (6% success), whereas testing 34 interface mutations gave 19 hits (56% success). Both the extent of stabilization and the spatial distribution of stabilizing mutations showed that the subunit interface was critical for stability. After mutations were combined, 2 very stable variants with 4 and 6 mutations were obtained, which in comparison to wild type (T _m ^app = 62 °C) displayed T _m ^app values of 80 and 85 °C, respectively. These two variants were also 5-fold more active at their optimum temperatures and tolerated high concentrations of isopropylamine and cosolvents. This allowed conversion of 100 mM acetophenone to (S)-1-phenylethylamine (>99% enantiomeric excess) with high yield (92%, in comparison to 24% with the wild-type transaminase). Crystal structures mostly confirmed the expected structural changes and revealed that the most stabilizing mutation, I154V, featured a rarely described stabilization mechanism: namely, removal of steric strain. The results show that computational interface redesign can be a rapid and powerful strategy for transaminase stabilization.

Rozhkova A, N. Zorov I, Korotkova O, P. Sinitsyn A, Schwaneberg U, D. Davari M. Engineering Robust Cellulases for Tailored Lignocellulosic Degradation Cocktails

Lignocellulosic biomass is a most promising feedstock in the production of second-generation biofuels. Efficient degradation of lignocellulosic biomass requires a synergistic action of several cellulases and hemicellulases. Cellulases depolymerize cellulose, the main polymer of the lignocellulosic biomass, to its building blocks. The production of cellulase cocktails has been widely explored, however, there are still some main challenges that enzymes need to overcome in order to develop a sustainable production of bioethanol. The main challenges include low activity, product inhibition, and the need to perform fine-tuning of a cellulase cocktail for each type of biomass. Protein engineering and directed evolution are powerful technologies to improve enzyme properties such as increased activity, decreased product inhibition, increased thermal stability, improved performance in non-conventional media, and pH stability, which will lead to a production of more efficient cocktails. In this review, we focus on recent advances in cellulase cocktail production, its current challenges, protein engineering as an efficient strategy to engineer cellulases, and our view on future prospects in the generation of tailored cellulases for biofuel production.

Creating Flavin Reductase Variants with Thermostable and Solvent-Tolerant Properties by Rational-Design Engineering

We have employed computational approaches-FireProt and FRESCO-to predict thermostable variants of the reductase component (C₁ ) of (4-hydroxyphenyl)acetate 3-hydroxylase. With the additional aid of experimental results, two C₁ variants, A166L and A58P, were identified as thermotolerant enzymes, with thermostability improvements of 2.6-5.6 °C and increased catalytic efficiency of 2- to 3.5-fold. After heat treatment at 45 °C, both of the thermostable C₁ variants remain active and generate reduced flavin mononucleotide (FMNH^- ) for reactions catalyzed by bacterial luciferase and by the monooxygenase C₂ more efficiently than the wild type (WT). In addition to thermotolerance, the A166L and A58P variants also exhibited solvent tolerance. Molecular dynamics (MD) simulations (6 ns) at 300-500 K indicated that mutation of A166 to L and of A58 to P resulted in structural changes with increased stabilization of hydrophobic interactions, and thus in improved thermostability. Our findings demonstrated that improvements in the thermostability of C₁ enzyme can lead to broad-spectrum uses of C₁ as a redox biocatalyst for future industrial applications.

Protein thermal stability

Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.

NADPH biosensor-based identification of an alcohol dehydrogenase variant with improved catalytic properties caused by a single charge reversal at the protein surface

Alcohol dehydrogenases (ADHs) are used in reductive biotransformations for the production of valuable chiral alcohols. In this study, we used a high-throughput screening approach based on the NADPH biosensor pSenSox and fluorescence-activated cell sorting (FACS) to search for variants of the NADPH-dependent ADH of Lactobacillus brevis (LbADH) with improved activity for the reduction of 2,5-hexanedione to (2R,5R)-hexanediol. In a library of approx. 1.4 × 10⁶ clones created by random mutagenesis we identified the variant LbADH^K71E. Kinetic analysis of the purified enzyme revealed that LbADH^K71E had a ~ 16% lowered K_M value and a 17% higher V_max for 2,5-hexanedione compared to the wild-type LbADH. Higher activities were also observed for the alternative substrates acetophenone, acetylpyridine, 2-hexanone, 4-hydroxy-2-butanone, and methyl acetoacetate. K71 is solvent-exposed on the surface of LbADH and not located within or close to the active site. Therefore, K71 is not an obvious target for rational protein engineering. The study demonstrates that high-throughput screening using the NADPH biosensor pSenSox represents a powerful method to find unexpected beneficial mutations in NADPH-dependent alcohol dehydrogenases that can be favorable in industrial biotransformations.

The Role of Electrostatics and Folding Kinetics on the Thermostability of Homologous Cold Shock Proteins

Understanding which aspects contribute to the thermostability of proteins is a challenge that has persisted for decades, and it is of great relevance for protein engineering. Several types of interactions can influence the thermostability of a protein. Among them, the electrostatic interactions have been a target of particular attention. Aiming to explore how this type of interaction can affect protein thermostability, this paper investigated four homologous cold shock proteins from psychrophilic, mesophilic, thermophilic, and hyperthermophilic organisms using a set of theoretical methodologies. It is well-known that electrostatics as well as hydrophobicity are key-elements for the stabilization of these proteins. Therefore, both interactions were initially analyzed in the native structure of each protein. Electrostatic interactions present in the native structures were calculated with the Tanford-Kirkwood model with solvent accessibility, and the amount of hydrophobic surface area buried upon folding was estimated by measuring both folded and extended structures. On the basis of Energy Landscape Theory, the local frustration and the simplified alpha-carbon structure-based model were modeled with a Debye-Hückel potential to take into account the electrostatics and the effects of an implicit solvent. Thermodynamic data for the structure-based model simulations were collected and analyzed using the Weighted Histogram Analysis and Stochastic Diffusion methods. Kinetic quantities including folding times, transition path times, folding routes, and Φ values were also obtained. As a result, we found that the methods are able to qualitatively infer that electrostatic interactions play an important role on the stabilization of the most stable thermophilic cold shock proteins, showing agreement with the experimental data.

pH and Charged Mutations Modulate Cold Shock Protein Folding and Stability: A Constant pH Monte Carlo Study

The folding and stability of proteins is a fundamental problem in several research fields. In the present paper, we have used different computational approaches to study the effects caused by changes in pH and for charged mutations in cold shock proteins from Bacillus subtilis (Bs-CspB). First, we have investigated the contribution of each ionizable residue for these proteins to their thermal stability using the TKSA-MC, a Web server for rational mutation via optimizing the protein charge interactions. Based on these results, we have proposed a new mutation in an already optimized Bs-CspB variant. We have evaluated the effects of this new mutation in the folding energy landscape using structure-based models in Monte Carlo simulation at constant pH, SBM-CpHMC. Our results using this approach have indicated that the charge rearrangements already in the unfolded state are critical to the thermal stability of Bs-CspB. Furthermore, the conjunction of these simplified methods was able not only to predict stabilizing mutations in different pHs but also to provide essential information about their effects in each stage of protein folding.

Rational Design of Chymotrypsin Inhibitor 2 by Optimizing Non-Native Interactions

Rational design of proteins via mutagenesis is crucial for several biotechnological applications. A significant challenge of the computational strategies used to predict optimized mutations is to understand the influence of each amino acid during the folding process. In the present work, chymotrypsin inhibitor 2 (CI2) and several of its designed mutants have been simulated using a non-native hydrophobic and electrostatic potential as a structure-based Cα model. Through these simulations, we could identify the most critical folding stage to accelerate CI2 and also the charged residues responsible for providing its thermostability. The replacement of ionizable residues for hydrophobic ones tended to promote the formation of the CI2 secondary structure in the early transition state, which speeds up folding. However, this same replacement destabilized the native structure, and there was a decrease in the protein thermostability. Such a simple method proved to be capable of providing valuable information about thermodynamics and kinetics of CI2 and its mutations, thus being a fast alternative to the study of rational protein design.