Stabilization of yeast hexokinase A by polyol osmolytes: correlation with the physicochemical properties of aqueous solutions

Osmolytes of the polyol series are known to accumulate in biological systems under stress and stabilize the structures of a wide variety of proteins. While increased surface tension of aqueous solutions has been considered an important factor in protein stabilization effect, glycerol is an exception, lowering the surface tension of water. To clarify this anomalous effect, the effect of a series of polyols on the thermal stability of a highly thermolabile two domain protein yeast hexokinase A has been investigated by differential scanning calorimetry and by monitoring loss in the biological activity of the enzyme as a function of time. A larger increase in the T(m) of domain 1 compared with that of domain 2, varying linearly with the number of hydroxyl groups in polyols, has been observed, sorbitol being the best stabilizer against both thermal as well as urea denaturation. Polyols help retain the activity of the enzyme considerably and a good correlation of the increase in T(m) (DeltaT(m)) and the retention of activity with the increase in the surface tension of polyol solutions, except glycerol, which breaks this trend, has been observed. However, the DeltaT(m) values show a linear correlation with apparent molal heat capacity and volume of aqueous polyol solutions including glycerol. These results suggest that while bulk solution properties contribute significantly to protein stabilization, interfacial properties are not always a good indicator of the stabilizing effect. A subtle balance of various weak binding and exclusion effects of the osmolytes mediated by water further regulates the stabilizing effect. Understanding these aspects is critical in the rational design of stable protein formulations.     

A single aquaporin gene encodes a water/glycerol/urea facilitator in Toxoplasma gondii with similarity to plant tonoplast intrinsic proteins

We describe a single aquaporin gene in Toxoplasma gondii which, surprisingly, has only 28% sequence similarity to the aquaglyceroporin of another apicomplexan parasite, Plasmodium falciparum. Sequence comparisons showed 47% similarity to water-specific plant aquaporins and the conservation of typical pore-forming residues. We established that the Toxoplasma aquaporin protein is a bifunctional membrane pore with intermediate water and high glycerol permeability. Furthermore, we identified hydroxyurea, an antineoplastic agent with inhibitory effects on parasite proliferation, as a permeant of this channel.     

A newly introduced salt bridge cluster improves structural and biophysical properties of de novo TIM barrels

Protein stability can be fine‐tuned by modifying different structural features such as hydrogen‐bond networks, salt bridges, hydrophobic cores, or disulfide bridges. Among these, stabilization by salt bridges is a major challenge in protein design and engineering since their stabilizing effects show a high dependence on the structural environment in the protein, and therefore are difficult to predict and model. In this work, we explore the effects on structure and stability of an introduced salt bridge cluster in the context of three different de novo TIM barrels. The salt bridge variants exhibit similar thermostability in comparison with their parental designs but important differences in the conformational stability at 25°C can be observed such as a highly stabilizing effect for two of the proteins but a destabilizing effect to the third. Analysis of the formed geometries of the salt bridge cluster in the crystal structures show either highly ordered salt bridge clusters or only single salt bridges. Rosetta modeling of the salt bridge clusters results in a good prediction of the tendency on stability changes but not the geometries observed in the three‐dimensional structures. The results show that despite the similarities in protein fold, the salt bridge clusters differently influence the structural and stability properties of the de novo TIM barrel variants depending on the structural background where they are introduced.     

Increased folding stability of TEM-1 beta-lactamase by in vitro selection

In vitro selections of stabilized proteins lead to more robust enzymes and, at the same time, yield novel insights into the principles of protein stability. We employed Proside, a method of in vitro selection, to find stabilized variants of TEM-1 β-lactamase from Escherichia coli. Proside links the increased protease resistance of stabilized proteins to the infectivity of a filamentous phage. Several libraries of TEM-1 β-lactamase variants were generated by error-prone PCR, and variants with increased protease resistance were obtained by raising temperature or guanidinium chloride concentration during proteolytic selections. Despite the small size of phage libraries, several strongly stabilizing mutations could be obtained, and a manual combination of the best shifted the profiles for thermal unfolding and temperature-dependent inactivation of β-lactamase by almost 20 °C to a higher temperature. The wild-type protein unfolds in two stages: from the native state via an intermediate of the molten-globule type to the unfolded form. In the course of the selections, the native protein was stabilized by 27 kJ mol^− 1 relative to the intermediate and the cooperativity of unfolding was strongly increased. Three of our stabilizing replacements (M182T, A224V, and R275L) had been identified independently in naturally occurring β-lactamase variants with extended substrate spectrum. In these variants, they acted as global suppressors of destabilizations caused by the mutations in the active site. The comparison between the crystal structure of our best variant and the crystal structure of the wild-type protein indicates that most of the selected mutations optimize helices and their packing. The stabilization by the E147G substitution is remarkable. It removes steric strain that originates from an overly tight packing of two helices in the wild-type protein. Such unfavorable van der Waals repulsions are not easily identified in crystal structures or by computational approaches, but they strongly reduce the conformational stability of a protein.     

Activation and stabilization of penicillin V acylase from streptomyces lavendulae in the presence of glycerol and glycols

Penicillin V acylase (EC 3.5.1.11) from Streptomyces lavendulae showed both enhanced activity and stability in mixed water/glycerol and water/glycols solvents. The catalytic activity was increased up to a critical concentration of these cosolvents, but further addition of the latter led to a gradual protein deactivation. The highest stabilizing effect was achieved in the presence of glycerol. Thermal stability was increased proportionally to the concentration of glycerol and glycols in the reaction mixture only if the amount added is below the threshold concentration. Reaction conditions that allow simultaneously enhanced activity and stability in the hydrolysis of penicillin V catalyzed by penicillin V acylase from S. lavendulae could be established.     

Effects of osmolytes on RNA secondary and tertiary structure stabilities and RNA-Mg2+ interactions

Osmolytes are small organic molecules accumulated by cells in response to osmotic stress. Although their effects on protein stability have been studied, there has been no systematic documentation of their influence on RNA. Here, the effects of nine osmolytes on the secondary and tertiary structure stabilities of six RNA structures of differing complexity and stability have been surveyed. Using thermal melting analysis, m-values (change in ΔG° of RNA folding per molal concentration of osmolyte) have been measured. All the osmolytes destabilize RNA secondary structure, although to different extents, probably because they favor solubilization of base surfaces. Osmolyte effects on tertiary structure, however, can be either stabilizing or destabilizing. We hypothesize that the stabilizing osmolytes have unfavorable interactions with the RNA backbone, which becomes less accessible to solvent in most tertiary structures. Finally, it was found that as a larger fraction of the negative charge of an RNA tertiary structure is neutralized by hydrated Mg²⁺, the RNA becomes less responsive to stabilizing osmolytes and may even be destabilized. The natural selection of osmolytes as protective agents must have been influenced by their effects on the stabilities of functional RNA structures, though in general, the effects of osmolytes on RNA and protein stabilities do not parallel each other. Our results also suggest that some osmolytes can be useful tools for studying intrinsically unstable RNA folds and assessing the mechanisms of Mg²⁺-induced RNA stabilization.     

Improving the stability of uricase from Aspergillus flavus by osmolytes: Use of response surface methodology for optimization of the enzyme stability

Urate oxidase (UOX) is used for the treatment of hyperuricemia and gout. The low stability of protein drugs such as UOX is a significant drawback in their liquid solutions. In this study, UOX from Aspergillus flavus was overexpressed and purified. Then, the effect of osmolytes (glycerol, sorbitol, and sucrose) on UOX thermal stability was studied. The enzyme stability was optimized by the aid of osmolytes using response surface methodology (RSM). The thermal stability results and Pareto analysis showed that sucrose has a stronger stabilization effect on UOX relative to the other osmolytes and increases the enzyme stability up to 2.4-fold, solely. RSM results indicated that the stability of UOX increased up to 4.9-fold in the presence of 12 % of glycerol and sucrose. The combinatorial effect results indicated that the osmolytes affect each other stabilization effects on UOX, probably through influencing the hydrogen bonds and mixing the osmolytes at different concentrations does not always lead to increased stability of UOX. The results revealed that the effect of sorbitol on UOX stability was negligible in the presence of glycerol and sucrose. Finally, this study results showed that RSM can be used to select a proper set of osmolytes for enzyme stabilization.     

Inferring Stabilizing Mutations from Protein Phylogenies: Application to Influenza Hemagglutinin

One selection pressure shaping sequence evolution is the requirement that a protein fold with sufficient stability to perform its biological functions. We present a conceptual framework that explains how this requirement causes the probability that a particular amino acid mutation is fixed during evolution to depend on its effect on protein stability. We mathematically formalize this framework to develop a Bayesian approach for inferring the stability effects of individual mutations from homologous protein sequences of known phylogeny. This approach is able to predict published experimentally measured mutational stability effects (ΔΔG values) with an accuracy that exceeds both a state-of-the-art physicochemical modeling program and the sequence-based consensus approach. As a further test, we use our phylogenetic inference approach to predict stabilizing mutations to influenza hemagglutinin. We introduce these mutations into a temperature-sensitive influenza virus with a defect in its hemagglutinin gene and experimentally demonstrate that some of the mutations allow the virus to grow at higher temperatures. Our work therefore describes a powerful new approach for predicting stabilizing mutations that can be successfully applied even to large, complex proteins such as hemagglutinin. This approach also makes a mathematical link between phylogenetics and experimentally measurable protein properties, potentially paving the way for more accurate analyses of molecular evolution.     

The inverse relationship between protein dynamics and thermal stability

Protein powders that are dehydrated or mixed with a glassy compound are known to have improved thermal stability. We present elastic and quasielastic neutron scattering measurements of the global dynamics of lysozyme and ribonuclease A powders. In the absence of solvation water, both protein powders exhibit largely harmonic motions on the timescale of the measurements. Upon partial hydration, quasielastic scattering indicative of relaxational processes appears at sufficiently high temperature. When the scattering spectrum are analyzed with the Kohlrausch-Williams-Watts formalism, the exponent beta decreases with increasing temperature, suggesting that multiple relaxation modes are emerging. When lysozyme was mixed with glycerol, its beta values were higher than the hydrated sample at comparable temperatures, reflecting the viscosity and stabilizing effects of glycerol.     

Production of 3-hydroxypropionic acid through propionaldehyde dehydrogenase PduP mediated biosynthetic pathway in Klebsiella pneumoniae

The pduP gene encodes a propionaldehyde dehydrogenase (PduP) was investigated for the role in 3-hydroxypropionic acid (3-HP) glycerol metabolism in Klebsiella pneumoniae. The enzyme assay showed that cell extracts from a pduP mutant strain lacked measurable dehydrogenase activity. Additionally, the mutant strain accumulated the cytotoxic intermediate metabolite 3-hydroxypropionaldehyde (3-HPA), causing both cell death and a lower final 3-HP titer. Ectopic expression of pduP restored normal cell growth to mutant. The enzymatic property of recombinant protein from Escherichia coli was examined, exhibiting a broad substrate specificity, being active on 3-HPA. The present work is thus the first to demonstrate the role of PduP in glycerol metabolism and biosynthesis of 3-HP.     

Biodiesel-derived crude glycerol bioconversion to animal feed: a sustainable option for a biodiesel refinery

This study examined the potential of producing an edible fungus, Rhizopus microsporus var. oligosporus, on biodiesel-derived crude glycerol. Prolific fungal growth was observed with a fungal biomass yield of 0.83 ± 0.02 (g biomass increase/g initial biomass) under optimal cultivation conditions (e.g. nonsterile crude glycerol at a concentration of 75% (w/v) with nutrient supplementation and without pH control). The potential of utilizing front-end processed banagrass (Pennisetum purpureum) juice as a source of nutrients for crude glycerol fermentation was evaluated with a 2.3-fold improvement in the fungal biomass yield. The glycerol-derived fungal biomass showed high amounts of threonine, one of the main limiting amino acids in non-ruminant feeds. An inexpensive fungal protein has the potential to reduce meat product prices by lowering the production costs of animal feeds. The application of fungal technology thus provides a unique sustainable option for biodiesel refineries by providing an additional source of revenue from fungal products.     

Engineering the Pattern of Protein Glycosylation Modulates the Thermostability of a GH11 Xylanase

Protein glycosylation is a common post-translational modification, the effect of which on protein conformational and stability is incompletely understood. Here we have investigated the effects of glycosylation on the thermostability of Bacillus subtilis xylanase A (XynA) expressed in Pichia pastoris. Intact mass analysis of the heterologous wild-type XynA revealed two, three, or four Hex_8–16GlcNAc₂ modifications involving asparagine residues at positions 20, 25, 141, and 181. Molecular dynamics (MD) simulations of the XynA modified with various combinations of branched Hex₉GlcNAc₂ at these positions indicated a significant contribution from protein-glycan interactions to the overall energy of the glycoproteins. The effect of glycan content and glycosylation position on protein stability was evaluated by combinatorial mutagenesis of all six potential N-glycosylation sites. The majority of glycosylated enzymes expressed in P. pastoris presented increased thermostability in comparison with their unglycosylated counterparts expressed in Escherichia coli. Steric effects of multiple glycosylation events were apparent, and glycosylation position rather than the number of glycosylation events determined increases in thermostability. The MD simulations also indicated that clustered glycan chains tended to favor less stabilizing glycan-glycan interactions, whereas more dispersed glycosylation patterns favored stabilizing protein-glycan interactions.     

Effects of osmolytes on arginine kinase from Euphausia superba: A study on thermal denaturation and aggregation

Investigations of energy-related enzymatic properties may provide valuable information about the mechanisms that are involved in the adaptation to extreme climatic environments. The protective effects of osmolytes on the thermal denaturation and aggregation of arginine kinase from E. superba (ESAK) was investigated. When the concentration of glycine, proline and glycerol increased, the relative activation was significantly enhanced, while the aggregation of ESAK during thermal denaturation was decreased. Spectrofluorometry results showed that the presence of these three osmolytes significantly decreased the tertiary structural changes of ESAK and that thermal denaturation directly induced ESAK aggregation. The results demonstrated that glycine, proline and glycerol not only prevented ESAK from inactivation and unfolding but also inhibited aggregation by stabilizing the ESAK conformation. We measured the ORF gene sequence of ESAK by RACE, and built the 3D structure of ESAK and osmolytes by homology models. The results showed that the docking energy was relatively low and that the clustering groups were spread to the surface of ESAK, indicating that osmolytes directly protect the surface of the protein. Our study provides important insight into the protective effects of osmolytes on ESAK folding.     

Thermal Stabilization of Dihydrofolate Reductase Using Monte Carlo Unfolding Simulations and Its Functional Consequences

Design of proteins with desired thermal properties is important for scientific and biotechnological applications. Here we developed a theoretical approach to predict the effect of mutations on protein stability from non-equilibrium unfolding simulations. We establish a relative measure based on apparent simulated melting temperatures that is independent of simulation length and, under certain assumptions, proportional to equilibrium stability, and we justify this theoretical development with extensive simulations and experimental data. Using our new method based on all-atom Monte-Carlo unfolding simulations, we carried out a saturating mutagenesis of Dihydrofolate Reductase (DHFR), a key target of antibiotics and chemotherapeutic drugs. The method predicted more than 500 stabilizing mutations, several of which were selected for detailed computational and experimental analysis. We find a highly significant correlation of r = 0.65–0.68 between predicted and experimentally determined melting temperatures and unfolding denaturant concentrations for WT DHFR and 42 mutants. The correlation between energy of the native state and experimental denaturation temperature was much weaker, indicating the important role of entropy in protein stability. The most stabilizing point mutation was D27F, which is located in the active site of the protein, rendering it inactive. However for the rest of mutations outside of the active site we observed a weak yet statistically significant positive correlation between thermal stability and catalytic activity indicating the lack of a stability-activity tradeoff for DHFR. By combining stabilizing mutations predicted by our method, we created a highly stable catalytically active E. coli DHFR mutant with measured denaturation temperature 7.2°C higher than WT. Prediction results for DHFR and several other proteins indicate that computational approaches based on unfolding simulations are useful as a general technique to discover stabilizing mutations.     

Chemical chaperone-mediated protein folding: stabilization of P22 tailspike folding intermediates by glycerol

Polyol co-solvents such as glycerol increase the thermal stability of proteins. This has been explained by preferential hydration favoring the more compact native over the denatured state. Although polyols are also expected to favor aggregation by the same mechanism, they have been found to increase the folding yields of some large, aggregation-prone proteins. We have used the homotrimeric phage P22 tailspike protein to investigate the origin of this effect. The folding of this protein is temperature-sensitive and limited by the stability of monomeric folding intermediates. At non-permissive temperature (>or=35 degrees C), tailspike refolding yields were increased significantly in the presence of 1-4 m glycerol. At low temperature, tailspike refolding is prevented when folding intermediates are destabilized by the addition of urea. Glycerol could offset the urea effect, suggesting that the polyol acts by stabilizing crucial folding intermediates and not by increasing solvent viscosity. The stabilization effect of glycerol on tailspike folding intermediates was confirmed in experiments using a temperature-sensitive folding mutant protein, by fluorescence measurements of subunit folding kinetics, and by temperature up-shift experiments. Our results suggest that the chemical chaperone effect of polyols observed in the folding of large proteins is due to preferential hydration favoring structure formation in folding intermediates.     

Cryopreservation of Steinernema carpocapsae and Heterorhabditis bacteriophora

A method for the cryopreservation of third-stage infective juveniles (IJ) of Steinernema carpocapsae and Heterorhabiditis bacteriophora was developed. Cryoprotection was achieved by incubating the nematodes in 22% glycerol (S. carpocapsae) or 14% glycerol (H. bacteriophora) for 24 hours, followed by 70% methanol at 0 C for 10 minutes. The viability of S. carpocapsae frozen in liquid nitrogen as 20 μl volumes spread over cover slip glass was > 80%. Survival of H. bacteriophora frozen on glass varied from 10 to 60% but was improved to > 80% by replacing the glass with filter paper. Cryopreservation and storage of 1-ml aliqots of S. carpocapsae IJ resulted in > 50% survival after 8 months; pathogenicity was retained and normal in vitro development took place. Trehalose and glycerol levels increased and glycogen levels decreased during incubation of S. carpocapsae IJ in glycerol. Normal levels of trehalose, glycerol and glycogen were restored during post freezing rehydration.     

Mannitol production from glycerol by resting cells of Candida magnoliae

Production of mannitol from glycerol by resting cells of Candida magnoliae under aerobic condition was investigated. The resting cells were suspended in aqueous solution of glycerol in Erlenmeyer flasks and incubated on rotary shaker. The samples were analyzed by ion exclusion–HPLC equipped with refractive index and UV detector. The resting cells of C. magnoliae produced mannitol from fructose, sucrose and glycerol but not from glucose. Addition of yeast extract and/or potassium phosphate to the glycerol solution adversely affected its conversion to mannitol. The conversion of glycerol to mannitol was dependent on oxygen availability. Using resting cells, the yield of mannitol was as high as 45%. This is probably the first report of conversion of glycerol to mannitol by osmophilic yeast.     

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.     

The effect of carbohydrate on the structure and stability of erythropoietin.

Erythropoietin is a glycoprotein hormone that stimulates the maturation of late erythroid progenitor cells. It has three N-linked and one O-linked carbohydrates which play an important role in the biosynthesis and biological activities of the protein. To determine the role the carbohydrate might have in maintaining the conformational stability of the protein, the protein expressed in mammalian cells (fully glycosylated), the asialo mammalian-expressed protein, and the protein expressed in Escherichia coli (no carbohydrate) were compared for their stability to guanidine HCl, pH, and temperature. Circular dichroism was used to follow protein unfolding. Both the intact and asialo mammalian-expressed proteins unfolded with a cooperative transition in guanidine HCl, with a midpoint at 1.75 M guanidine HCl. The E. coli-expressed material unfolded with a midpoint of 1.2 M guanidine HCl, and a delta G of unfolding which was 1.4 kcal/mol less than that of the two glycosylated molecules. The E. coli-derived protein was also significantly less stable to pH-induced conformational changes, showing a cooperative transition in 35% glycerol with a midpoint at pH 4.4, while both the intact and asialo mammalian-expressed molecules had a transition midpoint of pH 3.75 in the absence of glycerol, and approximately pH 3 in the presence of 35% glycerol. The E. coli-expressed molecule unfolded and precipitated upon heating to 44 degrees C, while the asialo and intact mammalian-expressed proteins remained soluble, with a Tm of 56 degrees C. From these experiments, the carbohydrate appears to play a critical role in stabilizing the erythropoietin molecule to denaturing conditions, and this increased stability does not depend on the presence of sialic acid.     

Comparison of the chaperon activity of glycerol and α-casein on amyloid formation of κ-casein in the presence of glycine and arginine

Amyloids are insoluble fibers which arise from inappropriately folded versions of proteins and have been associated with the pathology of many neurodegenerative diseases. α-Casein is one of the major components of the casein family which is known to show chaperone-like activity. Glycerol is a polyol compound which acts as a chemical chaperone to increase protein stability and inhibit protein aggregation. In this study, the effect of arginine and glycine on the chaperone ability of α-casein and glycerol against order aggregation of κ-casein was investigated and compared. We found that these additives reduced the chaperone ability of α-casein against the amyloid formation of κ-casein, especially in the presence of arginine. Importantly, our results show that the chaperone action of glycerol is enhanced in the presence of both arginine and glycine. Accordingly, our results suggest that these small molecules associated with glycerol, especially glycine, should be considered as a mechanism for the treatment of amyloid disease.