Protein stability in the presence of cosolutes

Stability of Protein Pharmaceuticals: Recent Advances

There have been significant advances in the formulation and stabilization of proteins in the liquid state over the past years since our previous review. Our mechanistic understanding of protein-excipient interactions has increased, allowing one to develop formulations in a more rational fashion. The field has moved towards more complex and challenging formulations, such as high concentration formulations to allow for subcutaneous administration and co-formulation. While much of the published work has focused on mAbs, the principles appear to apply to any therapeutic protein, although mAbs clearly have some distinctive features. In this review, we first discuss chemical degradation reactions. This is followed by a section on physical instability issues. Then, more specific topics are addressed: instability induced by interactions with interfaces, predictive methods for physical stability and interplay between chemical and physical instability. The final parts are devoted to discussions how all the above impacts (co-)formulation strategies, in particular for high protein concentration solutions.'

"Cooperative collapse" of the denatured state revealed through Clausius-Clapeyron analysis of protein denaturation phase diagrams

Protein phase diagrams have a unique potential to identify the presence of additional thermodynamic states even when non-2-state character is not readily apparent from the experimental observables used to follow protein unfolding transitions. Two-state analysis of the von Willebrand factor A3 domain has previously revealed a discrepancy in the calorimetric enthalpy obtained from thermal unfolding transitions as compared with Gibbs-Helmholtz analysis of free energies obtained from the Linear Extrapolation Method (Tischer and Auton, Prot Sci 2013; 22(9):1147-60). We resolve this thermodynamic conundrum using a Clausius-Clapeyron analysis of the urea-temperature phase diagram that defines how Δ H and the urea m-value interconvert through the slope of c_m versus T, ( ∂ c m / ∂ T ) = Δ H / ( m T ) . This relationship permits the calculation of Δ H at low temperature from m-values obtained through iso-thermal urea denaturation and high temperature m-values from Δ H obtained through iso-urea thermal denaturation. Application of this equation uncovers sigmoid transitions in both cooperativity parameters as temperature is increased. Such residual thermal cooperativity of Δ H and the m-value confirms the presence of an additional state which is verified to result from a cooperative phase transition between urea-expanded and thermally-compact denatured states. Comparison of the equilibria between expanded and compact denatured ensembles of disulfide-intact and carboxyamidated A3 domains reveals that introducing a single disulfide crosslink does not affect the presence of the additional denatured state. It does, however, make a small thermodynamically favorable free energy (∼-13 ± 1 kJ/mol) contribution to the cooperative denatured state collapse transition as temperature is raised and urea concentration is lowered. The thermodynamics of this "cooperative collapse" of the denatured state retain significant compensations between the enthalpy and entropy contributions to the overall free energy.

Osmolyte depletion viewed in terms of the dividing membrane and its work of expansion against osmotic pressure

How osmolytes enhance the folding, binding, and self-assembly of biological macromolecules at a microscopic scale has long been a matter of debate. Ambiguities persist on the key interpretive concepts, such as the "effective membrane" (which marks the boundary of the volume from which osmolytes are excluded) and the "free energy of exclusion" of osmolytes from biomolecular surfaces. In this paper, we formulate these elusive concepts based upon chemical thermodynamics and rigorous statistical thermodynamics (the Kirkwood-Buff theory). Positioning of the membrane at the osmotic dividing surface is crucial in order not to affect the thermodynamics of solvation. The notion of the free energy (work) of excluding osmolytes is refined to the expansion work against the osmotic pressure, which indeed describes the change of solvation free energy at dilute osmolyte concentrations.

Protein Stability: Enhancement and Measurement

This article defines protein stability, emphasizes its importance and surveys the field of protein stabilization, with summary reference to a selection of 2009-2015 publications. One can enhance stability by, in particular, protein engineering strategies and by chemical modification (including conjugation) in solution. General protocols are set out on how to measure a given protein's (1) kinetic thermal stability, and (2) oxidative stability, and (3) how to undertake chemical modification of a protein in solution.

Salinity Tolerance Mechanism of Economic Halophytes From Physiological to Molecular Hierarchy for Improving Food Quality

Soil salinity is becoming the key constraints factor to agricultural production. Therefore, the plant especially the crops possessing capacities of salt tolerance will be of great economic significance. The adaptation or tolerance of plant to salinity stress involves a series of physiological, metabolic and molecular mechanisms. Halophytes are the kind of organisms which acquire special salt tolerance mechanisms to respond to the salt tress and ensure normal growth and development under saline conditions in their lengthy evolutionary adaptation, so understanding how halophytes respond to salinity stress will provide us with methods and tactics to foster and develop salt resistant varieties of crops. The strategies in physiological and molecular level adopted by halophytes are various including the changes in photosynthetic and transpiration rate, the sequestration of Na+ to extracellular or vacuole, the regulation of stomata aperture and stomatal density, the accumulation and synthesis of the phytohormones as well as the relevant gene expression underlying these physiological traits, such as the stress signal transduction, the regulation of the transcription factors, the activation and expression of the transporter genes, the activation or inhibition of the synthetases and so on. This review focuses on the research advances of the regulating mechanisms in halophytes from physiological to molecular, which render the halophytes tolerance and adaption to salinity stress.

Synergy in protein-osmolyte mixtures

Virtually all taxa use osmolytes to protect cells against biochemical stress. Osmolytes often occur in mixtures, such as the classical combination of urea with TMAO (trimethylamine N-oxide) in cartilaginous fish or the cocktail of at least six different osmolytes in the kidney. The concentration patterns of osmolyte mixtures found in vivo make it likely that synergy between them plays an important role. Using statistical mechanical n-component Kirkwood-Buff theory, we show from first principles that synergy in protein-osmolyte systems can arise from two separable sources: (1) mutual alteration of protein surface solvation and (2) effects mediated through bulk osmolyte chemical activities. We illustrate both effects in a four-component system with the experimental example of the unfolding of a notch ankyrin domain in urea-TMAO mixtures, which make urea a less effective denaturant and TMAO a more effective stabilizer. Protein surface effects are primarily responsible for this synergy. The specific patterns of surface solvation point to denatured state expansion as the main factor, as opposed to direct competition.

Infinitely dilute partial molar properties of proteins from computer simulation

A detailed understanding of temperature and pressure effects on an infinitely dilute protein's conformational equilibrium requires knowledge of the corresponding infinitely dilute partial molar properties. Established molecular dynamics methodologies generally have not provided a way to calculate these properties without either a loss of thermodynamic rigor, the introduction of nonunique parameters, or a loss of information about which solute conformations specifically contributed to the output values. Here we implement a simple method that is thermodynamically rigorous and possesses none of the above disadvantages, and we report on the method's feasibility and computational demands. We calculate infinitely dilute partial molar properties for two proteins and attempt to distinguish the thermodynamic differences between a native and a denatured conformation of a designed miniprotein. We conclude that simple ensemble average properties can be calculated with very reasonable amounts of computational power. In contrast, properties corresponding to fluctuating quantities are computationally demanding to calculate precisely, although they can be obtained more easily by following the temperature and/or pressure dependence of the corresponding ensemble averages.

Protein stabilization and counteraction of denaturing effect of urea by glycine betaine

The counteraction of the denaturing effect of urea by osmolytes has been one of the most studied problems of osmolyte action. However, the possibility of synergy in osmolyte mixtures has often been neglected. Here, we report synergy in the glycine betaine (GB)-urea mixture by using a model peptide. The results show that in the GB-urea mixture, GB acts as a stronger osmolyte and urea becomes a weaker denaturing agent. This is reflected by an increase in the exclusion of GB from the peptide surface and a decrease in interactions between the peptide and urea. The cause of this synergistic behaviour includes direct interactions between GB and urea through hydrogen bonding, van der Waals interactions between them and strengthening of hydrogen bonding network of water in the GB-urea mixture. The results obtained in this study provide insights into osmolyte induced counteraction of the denaturing effect of urea.

Protein domain definition should allow for conditional disorder

Proteins are often classified in a binary fashion as either structured or disordered. However this approach has several deficits. Firstly, protein folding is always conditional on the physiochemical environment. A protein which is structured in some circumstances will be disordered in others. Secondly, it hides a fundamental asymmetry in behavior. While all structured proteins can be unfolded through a change in environment, not all disordered proteins have the capacity for folding. Failure to accommodate these complexities confuses the definition of both protein structural domains and intrinsically disordered regions. We illustrate these points with an experimental study of a family of small binding domains, drawn from the RNA polymerase of mumps virus and its closest relatives. Assessed at face value the domains fall on a structural continuum, with folded, partially folded, and near unstructured members. Yet the disorder present in the family is conditional, and these closely related polypeptides can access the same folded state under appropriate conditions. Any heuristic definition of the protein domain emphasizing conformational stability divides this domain family in two, in a way that makes no biological sense. Structural domains would be better defined by their ability to adopt a specific tertiary structure: a structure that may or may not be realized, dependent on the circumstances. This explicitly allows for the conditional nature of protein folding, and more clearly demarcates structural domains from intrinsically disordered regions that may function without folding.

Urea-temperature phase diagrams capture the thermodynamics of denatured state expansion that accompany protein unfolding

We have analyzed the thermodynamic properties of the von Willebrand factor (VWF) A3 domain using urea-induced unfolding at variable temperature and thermal unfolding at variable urea concentrations to generate a phase diagram that quantitatively describes the equilibrium between native and denatured states. From this analysis, we were able to determine consistent thermodynamic parameters with various spectroscopic and calorimetric methods that define the urea-temperature parameter plane from cold denaturation to heat denaturation. Urea and thermal denaturation are experimentally reversible and independent of the thermal scan rate indicating that all transitions are at equilibrium and the van't Hoff and calorimetric enthalpies obtained from analysis of individual thermal transitions are equivalent demonstrating two-state character. Global analysis of the urea-temperature phase diagram results in a significantly higher enthalpy of unfolding than obtained from analysis of individual thermal transitions and significant cross correlations describing the urea dependence of ΔH0 and ΔCP0 that define a complex temperature dependence of the m-value. Circular dichroism (CD) spectroscopy illustrates a large increase in secondary structure content of the urea-denatured state as temperature increases and a loss of secondary structure in the thermally denatured state upon addition of urea. These structural changes in the denatured ensemble make up ∼40% of the total ellipticity change indicating a highly compact thermally denatured state. The difference between the thermodynamic parameters obtained from phase diagram analysis and those obtained from analysis of individual thermal transitions illustrates that phase diagrams capture both contributions to unfolding and denatured state expansion and by comparison are able to decipher these contributions.

The disordered C-terminal domain of human DNA glycosylase NEIL1 contributes to its stability via intramolecular interactions

NEIL1 [Nei (endonuclease VIII)-like protein 1], one of the five mammalian DNA glycosylases that excise oxidized DNA base lesions in the human genome to initiate base excision repair, contains an intrinsically disordered C-terminal domain (CTD; ~100 residues), not conserved in its Escherichia coli prototype Nei. Although dispensable for NEIL1's lesion excision and AP lyase activities, this segment is required for efficient in vivo enzymatic activity and may provide an interaction interface for many of NEIL1's interactions with other base excision repair proteins. Here, we show that the CTD interacts with the folded domain in native NEIL1 containing 389 residues. The CTD is poised for local folding in an ordered structure that is induced in the purified fragment by osmolytes. Furthermore, deletion of the disordered tail lacking both Tyr and Trp residues causes a red shift in NEIL1's intrinsic Trp-specific fluorescence, indicating a more solvent-exposed environment for the Trp residues in the truncated protein, which also exhibits reduced stability compared to the native enzyme. These observations are consistent with stabilization of the native NEIL1 structure via intramolecular, mostly electrostatic, interactions that were disrupted by mutating a positively charged (Lys-rich) cluster of residues (amino acids 355-360) near the C-terminus. Small-angle X-ray scattering (SAXS) analysis confirms the flexibility and dynamic nature of NEIL1's CTD, a feature that may be critical to providing specificity for NEIL1's multiple, functional interactions.

Identifying core features of adaptive metabolic mechanisms for chronic heat stress attenuation contributing to systems robustness

The contribution of metabolism to heat stress may play a significant role in defining robustness and recovery of systems; either by providing the energy and metabolites required for cellular homeostasis, or through the generation of protective osmolytes. However, the mechanisms by which heat stress attenuation could be adapted through metabolic processes as a stabilizing strategy against thermal stress are still largely unclear. We address this issue through metabolomic and transcriptomic profiles for populations along a thermal cline where two seagrass species, Zostera marina and Zostera noltii, were found in close proximity. Significant changes captured by these profile comparisons could be detected, with a larger response magnitude observed in northern populations to heat stress. Sucrose, fructose, and myo-inositol were identified to be the most responsive of the 29 analyzed organic metabolites. Many key enzymes in the Calvin cycle, glycolysis and pentose phosphate pathways also showed significant differential expression. The reported comparison suggests that adaptive mechanisms are involved through metabolic pathways to dampen the impacts of heat stress, and interactions between the metabolome and proteome should be further investigated in systems biology to understand robust design features against abiotic stress.

Comparison of chemical and thermal protein denaturation by combination of computational and experimental approaches. II

Chemical and thermal denaturation methods have been widely used to investigate folding processes of proteins in vitro. However, a molecular understanding of the relationship between these two perturbation methods is lacking. Here, we combined computational and experimental approaches to investigate denaturing effects on three structurally different proteins. We derived a linear relationship between thermal denaturation at temperature T(b) and chemical denaturation at another temperature T(u) using the stability change of a protein (ΔG). For this, we related the dependence of ΔG on temperature, in the Gibbs-Helmholtz equation, to that of ΔG on urea concentration in the linear extrapolation method, assuming that there is a temperature pair from the urea (T(u)) and the aqueous (T(b)) ensembles that produces the same protein structures. We tested this relationship on apoazurin, cytochrome c, and apoflavodoxin using coarse-grained molecular simulations. We found a linear correlation between the temperature for a particular structural ensemble in the absence of urea, T(b), and the temperature of the same structural ensemble at a specific urea concentration, T(u). The in silico results agreed with in vitro far-UV circular dichroism data on apoazurin and cytochrome c. We conclude that chemical and thermal unfolding processes correlate in terms of thermodynamics and structural ensembles at most conditions; however, deviations were found at high concentrations of denaturant.

Volume exclusion and H-bonding dominate the thermodynamics and solvation of trimethylamine-N-oxide in aqueous urea

Trimethylamine-N-oxide (TMAO) and urea represent the extremes among the naturally occurring organic osmolytes in terms of their ability to stabilize/destabilize proteins. Their mixtures are found in nature and have generated interest in terms of both their physiological role and their potential use as additives in various applications (crystallography, drug formulation, etc.). Here we report experimental density and activity coefficient data for aqueous mixtures of TMAO with urea. From these data we derive the thermodynamics and solvation properties of the osmolytes, using Kirkwood–Buff theory. Strong hydrogen-bonding at the TMAO oxygen, combined with volume exclusion, accounts for the thermodynamics and solvation of TMAO in aqueous urea. As a result, TMAO behaves in a manner that is surprisingly similar to that of hard-spheres. There are two mandatory solvation sites. In plain water, these sites are occupied with water molecules, which are seamlessly replaced by urea, in proportion to its volume fraction. We discuss how this result gives an explanation both for the exceptionally strong exclusion of TMAO from peptide groups and for the experimentally observed synergy between urea and TMAO.

Osmolyte effects on protein stability and solubility: a balancing act between backbone and side-chains

In adaptation biology the discovery of intracellular osmolyte molecules that in some cases reach molar levels, raises questions of how they influence protein thermodynamics. We've addressed such questions using the premise that from atomic coordinates, the transfer free energy of a native protein (ΔG(tr,N)) can be predicted by summing measured water-to-osmolyte transfer free energies of the protein's solvent exposed side chain and backbone component parts. ΔG(tr,D) is predicted using a self avoiding random coil model for the protein, and ΔG(tr,D)-ΔG(tr,N), predicts the m-value, a quantity that measures the osmolyte effect on the N⇌D transition. Using literature and newly measured m-values we show 1:1 correspondence between predicted and measured m-values covering a range of 12 kcal/mol/M in protein stability for 46 proteins and 9 different osmolytes. Osmolytes present a range of side chain and backbone effects on N and D solubility and protein stability key to their biological roles.