Protein stabilization by the rational design of surface charge-charge interactions

Hyperthermostability of prawn ferritin nanocage facilitates its application as a robust nanovehicle for nutraceuticals

The favorable physicochemical properties are essential for the application of protein-based nanovehicles in the field of biomaterials. Herein, we found that the thermal stability of Marsupenaeus japonicus ferritin (MjFer) (T_m = 109.1 ± 0.4 °C) is markedly higher than human H-chain ferritin (HuHF) (T_m = 87.7 ± 0.3 °C), although they share a high structural similarity. Multiple results indicated that the promoted thermal stability of MjFer is mainly derived from the salt bridges located at the C₃ interface. Consequently, MjFer exhibits strong protective effects on encapsulated curcumin upon exposure at high temperatures. In contrast, most of the curcumin loaded HuHF composites precipitated rapidly under the same conditions. These findings elucidated the molecular mechanism of the hyperthermostability of MjFer and illustrated that MjFer could act as a robust insulation nanocarrier for bioactive compounds against various thermal treatments.

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.

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.

A critical review of five machine learning-based algorithms for predicting protein stability changes upon mutation

A number of machine learning (ML)-based algorithms have been proposed for predicting mutation-induced stability changes in proteins. In this critical review, we used hypothetical reverse mutations to evaluate the performance of five representative algorithms and found all of them suffer from the problem of overfitting. This approach is based on the fact that if a wild-type protein is more stable than a mutant protein, then the same mutant is less stable than the wild-type protein. We analyzed the underlying issues and suggest that the main causes of the overfitting problem include that the numbers of training cases were too small, and the features used in the models were not sufficiently informative for the task. We make recommendations on how to avoid overfitting in this important research area and improve the reliability and robustness of ML-based algorithms in general.

Coimmobilization of β-Agarase and α-Neoagarobiose Hydrolase for Enhancing the Production of 3,6-Anhydro-l-galactose

Here we report a simple and efficient method to produce 3,6-anhydro-l-galactose (l-AHG) and agarotriose (AO3) in one step by a multienzyme system with the coimmobilized β-agarase AgWH50B and α-neoagarobiose hydrolase K134D. K134D was obtained by AgaWH117 mutagenesis and showed improved thermal stability when immobilized via covalent bonds on functionalized magnetic nanoparticles. The obtained multienzyme biocatalyst was characterized by Fourier transform infrared spectroscopy (FTIR). Compared with free agarases, the coimmobilized agarases exhibited a relatively higher agarose-to-l-AHG conversion efficiency. The yield of l-AHG obtained with the coimmobilized agarases was 40.6%, which was 6.5% higher than that obtained with free agarases. After eight cycles, the multienzyme biocatalyst still preserved 46.4% of the initial activity. To the best of our knowledge, this is the first report where two different agarases were coimmobilized. These results demonstrated the feasibility of the new method to fabricate a new multienzyme system onto magnetic nanoparticles via covalent bonds to produce l-AHG.

Characterization of an arylsulfatase from a mutant library of Pseudoalteromonas carrageenovora arylsulfatase

A library of Pseudoalteromonas carrageenovora arylsulfatase mutants was constructed by introducing random mutagenesis using error-prone PCR. After screening, one mutant strain was obtained whose arylsulfatase had improved thermal stability. Protein sequence analysis revealed one amino acid substitution of H260L. The mutant arylsulfatase (named H260L) retained higher residual activity than wild-type enzyme (named WT) after incubation at 45, 50, 55 and 60°C for 60min. Thermal inactivation analysis showed that the half-life (t_1/2) value at 55°C for H260L was 40.6min, while that of WT was 9.1min. When p-nitrophenyl sulfate was used as a substrate, the optimal reaction temperature and pH for the mutant enzyme were 55°C and pH 8.0, respectively. H260L was stable over the pH range of 6.0-9.0. Inhibition assay with EDTA indicated that metal ions play an important role during the catalytic process of the mutant enzyme. The desulfation ratio against agar of Gracilaria lemaneiformis was 82%.

Protein-stabilizing and cell-penetrating properties of α-helix domain of 30Kc19 protein

The protein-stabilizing and cell-penetrating activities of Bombyx mori 30Kc19 α-helix domain (30Kc19α) are investigated. Recently, 30Kc19 protein has been studied extensively as it has both protein-stabilizing and cell-penetrating properties. However, it is unknown which part of 30Kc19 is responsible for those properties. 30Kc19 protein is composed of two distinct domains, an α-helix N-terminal domain (30Kc19α) and a β-trefoil C-terminal domain (30Kc19β). The authors construct and produce truncated forms of 30Kc19 to demonstrate their biological functions. Interestingly, 30Kc19α was shown to be responsible for both the protein-stabilizing and cell-penetrating properties of 30Kc19 protein. 30Kc19α shows even higher protein delivery activity than did whole 30Kc19 protein and has low cytotoxicity when added to cell culture medium. Therefore, based on its multifunctional properties, 30Kc19α can be developed as a novel candidate for a therapeutic protein carrier into various cells and tissues.

Adding an appropriate amino acid during crosslinking results in more stable crosslinked enzyme aggregates

Carrier free immobilization, especially crosslinked enzyme aggregates (CLEAs), has become an important design for biocatalysis in several areas. Adding amino acids during formation of CLEAs was found to give biocatalysts more stable at 55 °C and in the presence of 60% acetonitrile. The half-lives of CLEAs prepared with and without Arg addition were 21 and 15 h (subtilisin) and 4 and 1.6 h (α-chymotrypsin) at 55 °C, respectively. The corresponding half-lives during acetonitrile presence were 4.1 and 3.0 h (subtilisin) and 39 and 22 min (α-chymotrypsin), respectively. CLEAs made with Arg had higher percentages of alpha helix. CLEAs made by adding Lys, Ala, or Asp also were more stable. In the case of Thermomyces lanuginosus lipase (TLL), CLEA with Ala was even more stable than CLEA with Arg. The addition of a suitable amino acid, thus, enhances CLEA stabilities. The results are discussed in the light of earlier results on chemical modification of proteins and the observation that the Arg/Lys ratio is invariably high in the case of enzymes from thermophiles.

Thermodynamics of Ion Pair Formations Between Charged Poly(Amino Acid)s

Electrostatic interactions between the positively and negatively charged amino acids in proteins play an important role in macromolecular stability, binding, and recognition. Numerous amino acids in proteins are ionizable and may exist in negatively (e.g., Glu, Asp, Cys, Tyr) or positively (e.g., Arg, Lys, His, Orn) charged form dependent on pH and their pKas. In this work, isothermal titration calorimetry was used to determine the average standard values of thermodynamic parameters (the Gibbs free energy, enthalpy, entropy, and the heat capacity) of interaction between the positively charged amino acid homopolymers (polyarginine, polylysine, and polyornithine) and the negatively charged homopolymers (polyaspartic and polyglutamic acids). These values are of potential use in the computational models of interacting proteins and other biological macromolecules. The study showed that oppositely charged poly(amino acid)s bound each other with the stoichiometry of one positive to one negative charge. Arginine bound to the negatively charged amino acids with exothermic enthalpy and higher affinity than lysine. This result also suggests that positive charges in proteins should not be considered entirely equivalent if carried by lysine or arginine. The difference in binding energy of arginine and lysine association with the negatively charged amino acids was attributed to the enthalpy of the second ionic hydrogen bond formation between the guanidine and carboxylic groups. Despite the favorable enthalpic contribution, all such ion pair formation reactions were largely entropy-driven. Consistent with previously observed ionic interactions, the positive heat capacity was always observed during the amino acid ion pair formation.

The role of surface electrostatics on the stability, function and regulation of human cystathionine β-synthase, a complex multidomain and oligomeric protein

Human cystathionine β-synthase (hCBS) is a key enzyme of sulfur amino acid metabolism, controlling the commitment of homocysteine to the transsulfuration pathway and antioxidant defense. Mutations in hCBS cause inherited homocystinuria (HCU), a rare inborn error of metabolism characterized by accumulation of toxic homocysteine in blood and urine. hCBS is a complex multidomain and oligomeric protein whose activity and stability are independently regulated by the binding of S-adenosyl-methionine (SAM) to two different types of sites at its C-terminal regulatory domain. Here we study the role of surface electrostatics on the complex regulation and stability of hCBS using biophysical and biochemical procedures. We show that the kinetic stability of the catalytic and regulatory domains is significantly affected by the modulation of surface electrostatics through noticeable structural and energetic changes along their denaturation pathways. We also show that surface electrostatics strongly affect SAM binding properties to those sites responsible for either enzyme activation or kinetic stabilization. Our results provide new insight into the regulation of hCBS activity and stability in vivo with implications for understanding HCU as a conformational disease. We also lend experimental support to the role of electrostatic interactions in the recently proposed binding modes of SAM leading to hCBS activation and kinetic stabilization.

A method to rationally increase protein stability based on the charge-charge interaction, with application to lipase LipK107

We report a suite of enzyme redesign protocol based on the surface charge-charge interaction calculation, which is potentially applied to improve the stability of an enzyme without compromising its catalytic activity. Together with the experimental validation, we have released a suite of enzyme redesign algorithm Enzyme Thermal Stability System, written based on our model, for open access to meet the needs in wet labs. Lipk107, a lipase of a versatile industrial use, was chosen to test our software. Our calculation determined that four residues, D113, D149, D213, and D253, located on the surface of LipK107 were critical to the stability of the enzyme. The model was validated with mutagenesis at these four residues followed by stability and activity tests. LipK107 mutants D113A and D149K were more resistant to thermal inactivation with ∼10°C higher half-inactivation temperature than wild-type LipK107. Moreover, mutant D149K exhibited significant retention in residual activity under constant heat, showing a 14-fold increase in the half-inactivation time at 50°C. Activity tests showed that these mutants retained the equal or higher specific activity, among which noteworthy was the mutant D253A with as much as 20% higher activity. We suggest that our protocol could be used as a general guideline to redesign protein enzymes with increased stabilities and enhanced activities.

Engineering proteins for thermostability through rigidifying flexible sites

Engineering proteins for thermostability is an exciting and challenging field since it is critical for broadening the industrial use of recombinant proteins. Thermostability of proteins arises from the simultaneous effect of several forces such as hydrophobic interactions, disulfide bonds, salt bridges and hydrogen bonds. All of these interactions lead to decreased flexibility of polypeptide chain. Structural studies of mesophilic and thermophilic proteins showed that the latter need more rigid structures to compensate for increased thermal fluctuations. Hence flexibility can be an indicator to pinpoint weak spots for enhancing thermostability of enzymes. A strategy has been proven effective in enhancing proteins' thermostability with two steps: predict flexible sites of proteins firstly and then rigidify these sites. We refer to this approach as rigidify flexible sites (RFS) and give an overview of such a method through summarizing the methods to predict flexibility of a protein, the methods to rigidify residues with high flexibility and successful cases regarding enhancing thermostability of proteins using RFS.

Endoglucanases: insights into thermostability for biofuel applications

Obtaining bioethanol from cellulosic biomass involves numerous steps, among which the enzymatic conversion of the polymer to individual sugar units has been a main focus of the biotechnology industry. Among the cellulases that break down the polymeric cellulose are endoglucanases that act synergistically for subsequent hydrolytic reactions. The endoglucanases that have garnered relatively more attention are those that can withstand high temperatures, i.e., are thermostable. Although our understanding of thermostability in endoglucanases is incomplete, some molecular features that are responsible for increased thermostability have been recently identified. This review focuses on the investigations of endoglucanases and their implications for biofuel applications.

PROTS-RF: a robust model for predicting mutation-induced protein stability changes

The ability to improve protein thermostability via protein engineering is of great scientific interest and also has significant practical value. In this report we present PROTS-RF, a robust model based on the Random Forest algorithm capable of predicting thermostability changes induced by not only single-, but also double- or multiple-point mutations. The model is built using 41 features including evolutionary information, secondary structure, solvent accessibility and a set of fragment-based features. It achieves accuracies of 0.799,0.782, 0.787, and areas under receiver operating characteristic (ROC) curves of 0.873, 0.868 and 0.862 for single-, double- and multiple- point mutation datasets, respectively. Contrary to previous suggestions, our results clearly demonstrate that a robust predictive model trained for predicting single point mutation induced thermostability changes can be capable of predicting double and multiple point mutations. It also shows high levels of robustness in the tests using hypothetical reverse mutations. We demonstrate that testing datasets created based on physical principles can be highly useful for testing the robustness of predictive models.

In silico classification of proteins from acidic and neutral cytoplasms

Protein acidostability is a common problem in biopharmaceutical and other industries. However, it remains a great challenge to engineer proteins for enhanced acidostability because our knowledge of protein acidostabilization is still very limited. In this paper, we present a comparative study of proteins from bacteria with acidic (AP) and neutral cytoplasms (NP) using an integrated statistical and machine learning approach. We construct a set of 393 non-redundant AP-NP ortholog pairs and calculate a total of 889 sequence based features for these proteins. The pairwise alignments of these ortholog pairs are used to build a residue substitution propensity matrix between APs and NPs. We use Gini importance provided by the Random Forest algorithm to rank the relative importance of these features. A scoring function using the 10 most significant features is developed and optimized using a hill climbing algorithm. The accuracy of the score function is 86.01% in predicting AP-NP ortholog pairs and is 76.65% in predicting non-ortholog AP-NP pairs, suggesting that there are significant differences between APs and NPs which can be used to predict relative acidostability of proteins. The overall trends uncovered in the study can be used as general guidelines for designing acidostable proteins. To best of our knowledge, this work represents the first systematic comparative study of the acidostable proteins and their non-acidostable orthologs.

Increasing sequence diversity with flexible backbone protein design: the complete redesign of a protein hydrophobic core

Protein design tests our understanding of protein stability and structure. Successful design methods should allow the exploration of sequence space not found in nature. However, when redesigning naturally occurring protein structures, most fixed backbone design algorithms return amino acid sequences that share strong sequence identity with wild-type sequences, especially in the protein core. This behavior places a restriction on functional space that can be explored and is not consistent with observations from nature, where sequences of low identity have similar structures. Here, we allow backbone flexibility during design to mutate every position in the core (38 residues) of a four-helix bundle protein. Only small perturbations to the backbone, 1-2 Å, were needed to entirely mutate the core. The redesigned protein, DRNN, is exceptionally stable (melting point >140°C). An NMR and X-ray crystal structure show that the side chains and backbone were accurately modeled (all-atom RMSD = 1.3 Å).

Electrostatic contribution of surface charge residues to the stability of a thermophilic protein: benchmarking experimental and predicted pKa values

Optimization of the surface charges is a promising strategy for increasing thermostability of proteins. Electrostatic contribution of ionizable groups to the protein stability can be estimated from the differences between the pKa values in the folded and unfolded states of a protein. Using this pKa-shift approach, we experimentally measured the electrostatic contribution of all aspartate and glutamate residues to the stability of a thermophilic ribosomal protein L30e from Thermococcus celer. The pKa values in the unfolded state were found to be similar to model compound pKas. The pKa values in both the folded and unfolded states obtained at 298 and 333 K were similar, suggesting that electrostatic contribution of ionizable groups to the protein stability were insensitive to temperature changes. The experimental pKa values for the L30e protein in the folded state were used as a benchmark to test the robustness of pKa prediction by various computational methods such as H++, MCCE, MEAD, pKD, PropKa, and UHBD. Although the predicted pKa values were affected by crystal contacts that may alter the side-chain conformation of surface charged residues, most computational methods performed well, with correlation coefficients between experimental and calculated pKa values ranging from 0.49 to 0.91 (p<0.01). The changes in protein stability derived from the experimental pKa-shift approach correlate well (r = 0.81) with those obtained from stability measurements of charge-to-alanine substituted variants of the L30e protein. Our results demonstrate that the knowledge of the pKa values in the folded state provides sufficient rationale for the redesign of protein surface charges leading to improved protein stability.

Thermophilic proteins: insight and perspective from in silico experiments

Proteins from thermophilic and hyperthermophilic organisms are stable and function at high temperatures (50-100 °C). The importance of understanding the microscopic mechanisms underlying this thermal resistance is twofold: it is key for acquiring general clues on how proteins maintain their fold stable and for targeting those medical and industrial applications that aim at designing enzymes that can work under harsh conditions. In this tutorial review we first provide the general background of protein thermostability by specifically focusing on the structural and thermodynamic peculiarities; next, we discuss how computational studies based on Molecular Dynamics simulations can broaden and refine our knowledge on such special class of proteins.

PROTS: a fragment based protein thermo-stability potential

Designing proteins with enhanced thermo-stability has been a main focus of protein engineering because of its theoretical and practical significance. Despite extensive studies in the past years, a general strategy for stabilizing proteins still remains elusive. Thus effective and robust computational algorithms for designing thermo-stable proteins are in critical demand. Here we report PROTS, a sequential and structural four-residue fragment based protein thermo-stability potential. PROTS is derived from a nonredundant representative collection of thousands of thermophilic and mesophilic protein structures and a large set of point mutations with experimentally determined changes of melting temperatures. To the best of our knowledge, PROTS is the first protein stability predictor based on integrated analysis and mining of these two types of data. Besides conventional cross validation and blind testing, we introduce hypothetical reverse mutations as a means of testing the robustness of protein thermo-stability predictors. In all tests, PROTS demonstrates the ability to reliably predict mutation induced thermo-stability changes as well as classify thermophilic and mesophilic proteins. In addition, this white-box predictor allows easy interpretation of the factors that influence mutation induced protein stability changes at the residue level.

Increasing protein stability: importance of DeltaC(p) and the denatured state

Increasing the conformational stability of proteins is an important goal for both basic research and industrial applications. In vitro selection has been used successfully to increase protein stability, but more often site-directed mutagenesis is used to optimize the various forces that contribute to protein stability. In previous studies, we showed that improving electrostatic interactions on the protein surface and improving the beta-turn sequences were good general strategies for increasing protein stability, and used them to increase the stability of RNase Sa. By incorporating seven of these mutations in RNase Sa, we increased the stability by 5.3 kcal/mol. Adding one more mutation, D79F, gave a total increase in stability of 7.7 kcal/mol, and a melting temperature 28 degrees C higher than the wild-type enzyme. Surprisingly, the D79F mutation lowers the change in heat capacity for folding, DeltaC(p), by 0.6 kcal/mol/K. This suggests that this mutation stabilizes structure in the denatured state ensemble. We made other mutants that give some insight into the structure present in the denatured state. Finally, the thermodynamics of folding of these stabilized variants of RNase Sa are compared with those observed for proteins from thermophiles.

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

Here, we report the application of a computational approach that allows the rational design of enzymes with enhanced thermostability while retaining full enzymatic activity. The approach is based on the optimization of the energy of charge-charge interactions on the protein surface. We experimentally tested the validity of the approach on 2 human enzymes, acylphosphatase (AcPh) and Cdc42 GTPase, that differ in size (98 vs. 198-aa residues, respectively) and tertiary structure. We show that the designed proteins are significantly more stable than the corresponding WT proteins. The increase in stability is not accompanied by significant changes in structure, oligomerization state, or, most importantly, activity of the designed AcPh or Cdc42. This success of the design methodology suggests that it can be universally applied to other enzymes, on its own or in combination with the other strategies based on redesign of the interactions in the protein core.