Cosolvent effects on protein stability

Multicolor protein FRET with tryptophan, selective coumarin-cysteine labeling, and genetic acridonylalanine encoding

Site-specific fluorescence probes can be used to measure distances within proteins when used as part of a Förster resonance energy transfer (FRET) pair. Here we report the synthesis of a coumarin maleimide (Mcm-Mal) that is fluorogenic upon reaction with cysteine. We demonstrate that cysteine, acridonylalanine (Acd) double mutant proteins can be produced by unnatural amino acid mutagenesis and reacted with Mcm-Mal to generate Mcm/Acd labeled proteins for FRET studies. The Mcm/Acd FRET pair is minimally-perturbing, easy to install, and well-suited to studying protein distances in the 15-40 Å range. Furthermore, Mcm/Acd labeling can be combined with tryptophan fluorescence in three color FRET to monitor multiple interactions in one experiment.

Cosolvent Effects on Polymer Hydration Drive Hydrophobic Collapse

Water-mediated hydrophobic interactions play an important role in self-assembly processes, aqueous polymer solubility, and protein folding, to name a few. Cosolvents affect these interactions; however, the implications for hydrophobic polymer collapse and protein folding equilibria are not well-understood. This study examines cosolvent effects on the hydrophobic collapse equilibrium of a generic 32-mer hydrophobic polymer in urea, trimethylamine- N-oxide (TMAO), and acetone aqueous solutions using molecular dynamics simulations. Our results unveil a remarkable cosolvent-concentration-dependent behavior. Urea, TMAO, and acetone all shift the equilibrium toward collapsed structures below 2 M cosolvent concentration and, in turn, to unfolded structures at higher cosolvent concentrations, irrespective of the differences in cosolvent chemistry and the nature of cosolvent-water interactions. We find that weakly attractive polymer-water van der Waals interactions oppose polymer collapse in pure water, corroborating related observations reviewed by Ben-Amotz ( Annu. Rev. Phys. Chem. 2016, 67, 617-638). The cosolvents studied in the present work adsorb at the polymer/water interface and expel water molecules into the bulk, thereby effectively removing the dehydration energy penalty that opposes polymer collapse in pure water. At low cosolvent concentrations, this leads to cosolvent-induced stabilization of collapsed polymer structures. Only at sufficiently high cosolvent concentrations, polymer-cosolvent interactions favor polymer unfolding.

Aqueous ionic liquids and their influence on peptide conformations: denaturation and dehydration mechanisms

Low concentrated aqueous ionic liquids (ILs) and their influence on protein structures have attracted a lot of interest over the last few years. This can be mostly attributed to the fact that aqueous ILs, depending on the ion species involved, can be used as protein protectants or protein denaturants. Atomistic molecular dynamics (MD) simulations are performed in order to study the influence of different aprotic ILs on the properties of a short hairpin peptide. Our results reveal distinct binding and denaturation effects for 1-ethyl-3-methylimidazolium (EMIM) in combination with different anions, namely, chloride (CL), tetrafluoroborate (BF4) and acetate (ACE). The simulation outcomes demonstrate that the studied ILs with larger anions reveal a more pronounced accumulation behavior of the individual ion species around the peptide, which is accomplished by a stronger dehydration effect. We can relate these findings to the implications of the Kirkwood-Buff theory, which provides a thermodynamic explanation for the denaturation strength in terms of the IL accumulation behavior. The results for the spatial distribution functions, the binding energies and the local/bulk partition coefficients are in good agreement with metadynamics simulations in order to determine the energetically most stable peptide conformations. The free energy landscapes indicate a decrease of the denaturation strength in the order EMIM/ACE, EMIM/BF4 and EMIM/CL, which coincides with a decreasing size of the anion species. An analysis of the potential binding energies reveals that this effect is mainly of enthalpic nature.

Influence of cosolvents, self-crowding, temperature and pressure on the sub-nanosecond dynamics and folding stability of lysozyme

We studied the effects of temperature and hydrostatic pressure on the dynamical properties and folding stability of highly concentrated lysozyme solutions in the absence and presence of the osmolytes trimethylamine-N-oxide (TMAO) and urea. Elastic incoherent neutron scattering (EINS) was applied to determine the mean-squared displacement (MSD) of the protein's hydrogen atoms to yield insights into the effects of these cosolvents on the averaged sub-nanosecond dynamics in the pressure range from ambient up to 4000 bar. To evaluate the additional effect of self-crowding, two protein concentrations (80 and 160 mg mL^-1) were used. We observed a distinct effect of TMAO on the internal hydrogen dynamics, namely a reduced mobility. Urea, on the other hand, revealed no marked effect and consequently, no counteracting effect in an urea-TMAO mixture was observed. Different from the less concentrated protein solution, no significant effect of pressure on the MSD was observed for 160 mg mL^-1 lysozyme. The EINS experiments were complemented by Fourier-transform infrared (FTIR) spectroscopy measurements, which led to additional insights into the folding stability of lysozyme under the various environmental conditions. We observed a stabilization of the protein in the presence of the compatible osmolyte TMAO and a destabilization in the presence of urea against temperature and pressure for both protein concentrations. Additionally, we noticed a slight destabilizing effect upon self-crowding at very high protein concentration (160 mg mL^-1), which is attributable to transient destabilizing intermolecular interactions. Furthermore, a pressure-temperature diagram could be obtained for lysozyme at these high protein concentrations that mimics densely packed intracellular conditions.

Flory-type theories of polymer chains under different external stimuli

In this Review, we present a critical analysis of various applications of the Flory-type theories to a theoretical description of the conformational behavior of single polymer chains in dilute polymer solutions under a few external stimuli. Different theoretical models of flexible polymer chains in the supercritical fluid are discussed and analysed. Different points of view on the conformational behavior of the polymer chain near the liquid-gas transition critical point of the solvent are presented. A theoretical description of the co-solvent-induced coil-globule transitions within the implicit-solvent-explicit-co-solvent models is discussed. Several explicit-solvent-explicit-co-solvent theoretical models of the coil-to-globule-to-coil transition of the polymer chain in a mixture of good solvents (co-nonsolvency) are analysed and compared with each other. Finally, a new theoretical model of the conformational behavior of the dielectric polymer chain under the external constant electric field in the dilute polymer solution with an explicit account for the many-body dipole correlations is discussed. The polymer chain collapse induced by many-body dipole correlations of monomers in the context of statistical thermodynamics of dielectric polymers is analysed.

The temperature dependence of the Hofmeister series: thermodynamic fingerprints of cosolute-protein interactions

The Hofmeister series is a universal homologous series to rank ion-specific effects on biomolecular properties such as protein stability or aggregation propensity. Although this ranking is widely used, outliers and exceptions are discussed controversially and a molecular level understanding is still lacking. Studying the thermal unfolding equilibrium of RNase A, we here show that this ambiguity arises from the oversimplified approach to determine the ion rankings. Instead of measuring salt effects on a single point of the protein folding stability curve (e.g. the melting point T_m), we here consider the salt induced shifts of the entire protein 'stability curve' (the temperature dependence of the unfolding free energy change, ΔG_u(T)). We found multiple intersections of these curves, pinpointing a widely ignored fact: the Hofmeister cation and anion rankings are temperature dependent. We further developed a novel classification scheme of cosolute effects based on their thermodynamic fingerprints, reaching beyond salt effects to non-electrolytes.

Combined influence of ectoine and salt: spectroscopic and numerical evidence for compensating effects on aqueous solutions

Ectoine is an important osmolyte, which allows microorganisms to survive in extreme environmental salinity. The hygroscopic effects of ectoine in pure water can be explained by a strong water binding behavior whereas a study on the effects of ectoine in salty solution is yet missing. We provide Raman spectroscopic evidence that the influence of ectoine and NaCl are opposing and completely independent of each other. The effect can be explained by the formation of strongly hydrogen-bonded water molecules around ectoine which compensate the influence of the salt on the water dynamics. The mechanism is corroborated by first principles calculations and broadens our understanding of zwitterionic osmolytes in aqueous solution. Our findings allow us to provide a possible explanation for the relatively high osmolyte concentrations in halotolerant bacteria.

Biochemical Characterization of Recombinant Thermostable Cohnella sp. A01 β-Glucanase

Background

Typically, non-cellulytic glucanase, including fungi and yeast cell wall hydrolyzing enzymes, are released by some symbiotic fungi and plants during the mycoparasitic fungi attack on plants. These enzymes are known as the defense mechanisms of plants. This study intends to investigate the biochemical properties of β-1,6-glucanase (bg16M) from native thermophilic bacteria, Cohnella A01.

Methods

bg16M gene was cloned and expressed in E. coli BL21 (DE3). The enzyme was purified utilizing Ni-NTA nikcle sepharose column. Pustulan and laminarin were selected as substrates in enzyme assay. The purified bg16M enzyme was treated with different pH, temperature, metal ions, and detergents.

Results

The expressed protein, including 639 amino acids, showed a high similarity with the hydrolytic glycosylated family 30. The molecular weight of enzyme was 64 kDa, and purification yield was 46%. The bg16M demonstrated activity as 4.83 U/ml on laminarin and 2.88 U/ml on pustulan. The optimum pH and temperature of the enzyme were 8 and 50 °C, respectively. The enzyme had an appropriate stability at high temperatures and in the pH range of 7 to 9, showing acceptable stability, while it did not lose enzymatic activity completely at acidic or basic pH. None of the studied metal ions and chemical compounds was the activator of bg16M, and urea, SDS, and copper acted as enzyme inhibitors.

Conclusion

Biochemical characterization of this enzyme revealed that bg16M can be applied in beverage industries and medical sectors because of its high activity, as well as thermal and alkaline stability.

Effects of Cosolvents and Macromolecular Crowding on the Phase Transitions and Temperature-Pressure Stability of Chiral and Racemic Poly-Lysine

FTIR spectroscopy has been used to reveal the effects of different types of cosolvents (TMAO, urea) as well as macromolecular crowding (using the crowding agent Ficoll) on the temperature and pressure dependent structure of poly-L-lysine, poly-D-lysine and their racemic mixture. Compared to the effects of cosolvents on the unfolding transition of proteins, their effects on the α-helix to aggregated β-sheet transition of polylysine are quite small. High hydrostatic pressure has been found to favor the α-helical state over the aggregated β-sheet structure which is reflected in a volume decrease of ΔV=−32 mL mol−1, indicating that the packing mode is more efficient in the α-helical structure. Both, addition of urea and TMAO lead to a decrease in pressure stability of the aggregated β-sheet structure, which is accompanied by a three-fold decrease in ΔV, whereas the macromolecular crowder has little effect on the β-to-α transition. The more than 3 kbar higher β-to-α transition pressure of the racemic mixture compared with PLL confirms the drastic stabilization of β-sheet aggregates if the stereoisomers PLL and PDL are combined. Changes in hydration and packing of the polypeptide occurs upon interaction and fine packing of the polypeptide’s chains of opposed chirality, which are slightly modulated by the properties of cosolute and crowding, only. The underlying solvational and packing mechanisms observed here may be decisive factors responsible for the spontaneous protein aggregation in general and, as such, may shed additional light on the molecular basis of amyloid-associated diseases.

Aqueous ionic liquids influence the disulfide bond isoform equilibrium in conotoxin AuIB: a consequence of the Hofmeister effect?

The appearance of several disulfide bond isoforms in multiple cysteine containing venom peptides poses a significant challenge in their synthesis and purification under laboratory conditions. Recent experiments suggest that careful tuning of solvent and temperature conditions can propel the disulfide bond isoform equilibrium in favor of the most potent, native form. Certain aqueous ionic liquids (ILs) have proven significantly useful as solvents for this purpose, while exceptions have also been noted. To elucidate the molecular level origin behind such a preference, we report a detailed explicit solvent replica exchange molecular dynamics study of a conotoxin, AuIB, in pure water and four different aqueous IL solutions (~45-60% v/v). The ILs studied here are comprised of cations like 1-ethyl-3-methyl-imidazolium (Im₂₁⁺) or 1-butyl-3-methyl-imidazolium (Im₄₁⁺) coupled with either acetate (OAc^-) or chloride (Cl^-) as the counter anion. Our simulations unfold interesting features of the conformational spaces sampled by the peptide and its solvation in pure water and aqueous IL solutions. Detailed investigation into populations of the globular disulfide bond isoform of AuIB in aqueous IL solutions reveal distinct trends which might be related to the Hofmeister effect of the cation and anion of the IL and of specific interactions of the aqueous IL solutions with the peptide. In accordance with experimental observations, the aqueous [Im₂₁][OAc] solution is found to promote the highest globular isoform population in AuIB.

Changing relations between proteins and osmolytes: a choice of nature

The stabilization and destabilization of the protein in the presence of any additive is mainly attributed to its preferential exclusion from protein surface and its preferential binding to the protein surface, respectively.

Influence of compatible solute ectoine on distinct DNA structures: thermodynamic insights into molecular binding mechanisms and destabilization effects

We study ectoine-induced destabilization effects on DNA hairpins by a combination of atomistic molecular dynamics simulations, experiments, and theoretical approaches.

The effects of glycine, TMAO and osmolyte mixtures on the pressure dependent enzymatic activity of α-chymotrypsin

Different natural osmolytes modulate the pressure dependent enzyme kinetics in different ways.

Aqueous TMAO solutions as seen by theoretical THz spectroscopy: hydrophilic versus hydrophobic water

All THz resonances of aqueous TMAO solutions are computed and assigned based on ab initio molecular dynamics simulations.

Exploring the influence of natural cosolvents on the free energy and conformational landscape of filamentous actin and microtubules

Natural osmolytes have a significant influence on the temperature- and pressure-dependent stability of filamentous actin and microtubules.

Hydration in aqueous solutions of ectoine and hydroxyectoine

We explore the influence of the two osmolytes ectoine and hydroxyectoine on the structure of pure water and aqueous NaCl solutions using non-resonant X-ray Raman scattering spectroscopy at the oxygen K-edge.

The Effect of Natural Osmolyte Mixtures on the Temperature-Pressure Stability of the Protein RNase A

In biological cells, osmolytes appear as complex mixtures with variable compositions, depending on the particular environmental conditions of the organism. Based on various spectroscopic, thermodynamic and small-angle scattering data, we explored the effect of two different natural osmolyte mixtures, which are found in shallow-water and deep-sea shrimps, on the temperature and pressure stability of a typical monomeric protein, RNase A. Both natural osmolyte mixtures stabilize the protein against thermal and pressure denaturation. This effect seems to be mainly caused by the major osmolyte components of the osmolyte mixtures, i.e. by glycine and trimethylamine-N-oxide (TMAO), respectively. A minor compaction of the structure, in particular in the unfolded state, seems to be largely due to TMAO. Differences in thermodynamic properties observed for glycine and TMAO, and hence also for the two osmolyte mixtures, are most likely due to different solvation properties and interactions with the protein. Different from TMAO, glycine seems to interact with the amino acid side chains and/or the backbone of the protein, thus competing with hydration water and leading to a less hydrated protein surface.

Ab Initio Simulations of Water Dynamics in Aqueous TMAO Solutions: Temperature and Concentration Effects

We use ab initio molecular dynamics simulation to study the effect of hydrophobic groups on the dynamics of water molecules in aqueous solutions of trimethylamine N-oxide (TMAO). We show that hydrophobic groups induce a moderate (<2-fold) slowdown of water reorientation and hydrogen-bond dynamics in dilute solutions, but that this slowdown rapidly increases with solute concentration. In addition, the slowdown factor is found to vary very little with temperature, thus suggesting an entropic origin. All of these results are in quantitative agreement with prior classical molecular dynamics simulations and with the previously suggested excluded-volume model. The hydrophilic TMAO headgroup is found to affect water dynamics more strongly than the hydrophobic moiety, and the magnitude of this slowdown is very sensitive to the strength of the water-solute hydrogen-bond.

Conformational Ensembles of α-Synuclein Derived Peptide with Different Osmolytes from Temperature Replica Exchange Sampling

Intrinsically disordered proteins (IDP) are a class of proteins that do not have a stable three-dimensional structure and can adopt a range of conformations playing various vital functional role. Alpha-synuclein is one such IDP which can aggregate into toxic protofibrils and has been associated largely with Parkinson's disease (PD) along with other neurodegenerative diseases. Osmolytes are small organic compounds that can alter the environment around the proteins by acting as denaturants or protectants for the proteins. In the present study, we have conducted a series of replica exchange molecular dynamics simulations to explore the role of osmolytes, urea which is a denaturant and TMAO (trimethylamine N-oxide), a protecting osmolyte, in aggregation and conformations of the synuclein peptide. We observed that both the osmolytes have significantly distinct impacts on the peptide and led to transitions of the conformations of the peptide from one state to other. Our findings highlighted that urea attenuated peptide aggregation and resulted in the formation of extended peptide structures whereas TMAO led to compact and folded forms of the peptide.

Effects of Trimethylamine-N-oxide on the Conformation of Peptides and its Implications for Proteins

To provide insights into the stabilizing mechanisms of trimethylamine-N-oxide (TMAO) on protein structures, we perform all-atom molecular dynamics simulations of peptides and the Trp-cage miniprotein. The effects of TMAO on the backbone and charged residues of peptides are found to stabilize compact conformations, whereas effects of TMAO on nonpolar residues lead to peptide swelling. This suggests competing mechanisms of TMAO on proteins, which accounts for hydrophobic swelling, backbone collapse, and stabilization of charge-charge interactions. These mechanisms are observed in Trp cage.

Conditionally disordered proteins: bringing the environment back into the fold

For many proteins, biological function requires the folding of the polypeptide chain into a unique and persistent tertiary structure. This review concerns proteins that adopt a specific tertiary structure to function, but are otherwise partially or completely disordered. The biological cue for protein folding is environmental perturbation or minor post-translational modification. Hence, we term these proteins conditionally disordered. Many of these proteins recognize and bind other molecules, and conditional disorder has been hypothesized to allow for more nuanced control and regulation of binding processes. However, this remains largely unproven. The sequences of conditionally disordered proteins suggest their propensity to fold; yet, under the standard laboratory conditions, they do not do so, which may appear surprising. We argue that the surprise results from the failure to consider the role of the environment in protein structure formation and that conditional disorder arises as a natural consequence of the marginal stability of the folded state.

Putting the Piezolyte Hypothesis under Pressure

A group of small molecules that stabilize proteins against high hydrostatic pressure has been classified as piezolytes, a subset of stabilizing cosolutes. This distinction would imply that piezolytes counteract the effects of high hydrostatic pressure through effects on the volumetric properties of the protein. The purpose of this study was to determine if cosolutes proposed to be piezolytes have an effect on the volumetric properties of proteins through direct experimental measurements of volume changes upon unfolding of model proteins lysozyme and ribonuclease A, in solutions containing varying cosolute concentrations. Solutions containing the proposed piezolytes glutamate, sarcosine, and betaine were used, as well as solutions containing the denaturants guanidinium hydrochloride and urea. Changes in thermostability were monitored using differential scanning calorimetry whereas changes in volume were monitored using pressure perturbation calorimetry. Our findings indicate that increasing stabilizing cosolute concentration increases the stability and transition temperature of the protein, but does not change the temperature dependence of volume changes upon unfolding. The results suggest that the pressure stability of a protein in solution is not directly affected by the presence of these proposed piezolytes, and so they cannot be granted this distinction.

Effect of alcohol on the structure of cytochrome C: FCS and molecular dynamics simulations

Effect of ethanol on the size and structure of a protein cytochrome C (Cyt C) is investigated using fluorescence correlation spectroscopy (FCS) and molecular dynamics (MD) simulations. For FCS studies, Cyt C is covalently labeled with a fluorescent probe, alexa 488. FCS studies indicate that on addition of ethanol, the size of the protein varies non-monotonically. The size of Cyt C increases (i.e., the protein unfolds) on addition of alcohol (ethanol) up to a mole fraction of 0.2 (44.75% v/v) and decreases at higher alcohol concentration. In order to provide a molecular origin of this structural transition, we explore the conformational free energy landscape of Cyt C as a function of radius of gyration (R_g) at different compositions of water-ethanol binary mixture using MD simulations. Cyt C exhibits a minimum at R_g ∼ 13 Å in bulk water (0% alcohol). Upon increasing ethanol concentration, a second minimum appears in the free energy surface with gradually larger R_g up to χ_EtOH ∼ 0.2 (44.75% v/v). This suggests gradual unfolding of the protein. At a higher concentration of alcohol (χ_EtOH > 0.2), the minimum at large R_g vanishes, indicating compaction. Analysis of the contact map and the solvent organization around protein indicates a preferential solvation of the hydrophobic residues by ethanol up to χ_EtOH = 0.2 (44.75% v/v) and this causes the gradual unfolding of the protein. At high concentration (χ_EtOH = 0.3 (58% v/v)), due to structural organization in bulk water-ethanol binary mixture, the extent of preferential solvation by ethanol decreases. This causes a structural transition of Cyt C towards a more compact state.

Experimental Atom-by-Atom Dissection of Amide-Amide and Amide-Hydrocarbon Interactions in H2O

Quantitative information about amide interactions in water is needed to understand their contributions to protein folding and amide effects on aqueous processes and to compare with computer simulations. Here we quantify interactions of urea, alkylated ureas, and other amides by osmometry and amide-aromatic hydrocarbon interactions by solubility. Analysis of these data yields strengths of interaction of ureas and naphthalene with amide sp2O, amide sp2N, aliphatic sp3C, and amide and aromatic sp2C unified atoms in water. Interactions of amide sp2O with urea and naphthalene are favorable, while amide sp2O-alkylurea interactions are unfavorable, becoming more unfavorable with increasing alkylation. Hence, amide sp2O-amide sp2N interactions (proposed n-σ* hydrogen bond) and amide sp2O-aromatic sp2C (proposed n-π*) interactions are favorable in water, while amide sp2O-sp3C interactions are unfavorable. Interactions of all ureas with sp3C and amide sp2N are favorable and increase in strength with increasing alkylation, indicating favorable sp3C-amide sp2N and sp3C-sp3C interactions. Naphthalene results show that aromatic sp2C-amide sp2N interactions in water are unfavorable while sp2C-sp3C interactions are favorable. These results allow interactions of amide and hydrocarbon moieties and effects of urea and alkylureas on aqueous processes to be predicted or interpreted in terms of structural information. We predict strengths of favorable urea-benzene and N-methylacetamide interactions from experimental information to compare with simulations and indicate how amounts of hydrocarbon and amide surfaces buried in protein folding and other biopolymer processes and transition states can be determined from analysis of urea and diethylurea effects on equilibrium and rate constants.

Crowding in Cellular Environments at an Atomistic Level from Computer Simulations

The effects of crowding in biological environments on biomolecular structure, dynamics, and function remain not well understood. Computer simulations of atomistic models of concentrated peptide and protein systems at different levels of complexity are beginning to provide new insights. Crowding, weak interactions with other macromolecules and metabolites, and altered solvent properties within cellular environments appear to remodel the energy landscape of peptides and proteins in significant ways including the possibility of native state destabilization. Crowding is also seen to affect dynamic properties, both conformational dynamics and diffusional properties of macromolecules. Recent simulations that address these questions are reviewed here and discussed in the context of relevant experiments.

TMAO-Protein Preferential Interaction Profile Determines TMAO's Conditional In Vivo Compatibility

Trimethylamine N-oxide (TMAO) exemplifies how Nature uses the solute effect as a simple chemical strategy to cope with hydrodynamic pressure or urea stress to maintain proteostasis. It is a gut-microbe-generated metabolite that strongly promotes the development of atherosclerosis. It remains unclear how TMAO exerts its effects. In this study, we experimentally characterized the profile of the preferential interaction potential of TMAO with proteins, a thermodynamic key to understanding the effects of TMAO on protein processes and the distinction of TMAO among osmolytes. TMAO is thus found to be highly preferentially excluded from most types of protein surface, which explains why TMAO is a strong globular protein stabilizer and identifies the dominant stabilizing factor as the unfavorable interaction of TMAO with the hydrophobic surface exposed upon unfolding. We dissected the mechanism of the counteracting effects of TMAO and urea: the contrary feature of the interaction profiles of the two solutes maximizes the possibility for them to offset each other's perturbing effect on protein processes. The interaction profile also predicts that TMAO promotes aggregation of amyloidogenic intrinsically disordered peptide, as demonstrated here in Aβ42, and that TMAO has a strong potential to impact protein processes in the absence of stressors. Our data suggest that although TMAO is an evolutionally selected chemical chaperone for some organisms or organs, its compatibility in vivo is conditional and determined by its interaction profile with biopolymers and the nature of the essential biopolymer processes. Our thermodynamic framework plus the TMAO-protein interaction profile provides a basis for exploring the broad biological significance of TMAO, including its pathological impact in the absence of stressors. We argue for the general importance of controlling in vivo background solutes and the pathological significance of a control failure.

Aqueous ionic liquids and their effects on protein structures: an overview on recent theoretical and experimental results

Ionic liquids (ILs) are used in a variety of technological and biological applications. Recent experimental and simulation results reveal the influence of aqueous ionic liquids on the stability of protein and enzyme structures. Depending on different parameters like the concentration and the ion composition, one can observe distinct stabilization or denaturation mechanisms for various ILs. In this review, we summarize the main findings and discuss the implications with regard to molecular theories of solutions and specific ion effects. A preferential binding model is introduced in order to discuss protein-IL effects from a statistical mechanics perspective. The value of the preferential binding coefficient determines the strength of the ion influence and indicates a shift of the chemical equilibrium either to the native or the denatured state of the protein. We highlight the role of water in order to explain the self-association behavior of the IL species and discuss recent experimental and simulation results in the light of the observed binding effects.

Preferential binding effects on protein structure and dynamics revealed by coarse-grained Monte Carlo simulation

The effect of preferential binding of solute molecules within an aqueous solution on the structure and dynamics of the histone H3.1 protein is examined by a coarse-grained Monte Carlo simulation. The knowledge-based residue-residue and hydropathy-index-based residue-solvent interactions are used as input to analyze a number of local and global physical quantities as a function of the residue-solvent interaction strength (f). Results from simulations that treat the aqueous solution as a homogeneous effective solvent medium are compared to when positional fluctuations of the solute molecules are explicitly considered. While the radius of gyration (Rg) of the protein exhibits a non-monotonic dependence on solvent interaction over a wide range of f within an effective medium, an abrupt collapse in Rg occurs in a narrow range of f when solute molecules rapidly bind to a preferential set of sites on the protein. The structure factor S(q) of the protein with wave vector (q) becomes oscillatory in the collapsed state, which reflects segmental correlations caused by spatial fluctuations in solute-protein binding. Spatial fluctuations in solute binding also modify the effective dimension (D) of the protein in fibrous (D ∼ 1.3), random-coil (D ∼ 1.75), and globular (D ∼ 3) conformational ensembles as the interaction strength increases, which differ from an effective medium with respect to the magnitude of D and the length scale.

Simulated pressure denaturation thermodynamics of ubiquitin

Simulations of protein thermodynamics are generally difficult to perform and provide limited information. It is desirable to increase the degree of detail provided by simulation and thereby the potential insight into the thermodynamic properties of proteins. In this study, we outline how to analyze simulation trajectories to decompose conformation-specific, parameter free, thermodynamically defined protein volumes into residue-based contributions. The total volumes are obtained using established methods from Fluctuation Solution Theory, while the volume decomposition is new and is performed using a simple proximity method. Native and fully extended ubiquitin are used as the test conformations. Changes in the protein volumes are then followed as a function of pressure, allowing for conformation-specific protein compressibility values to also be obtained. Residue volume and compressibility values indicate significant contributions to protein denaturation thermodynamics from nonpolar and coil residues, together with a general negative compressibility exhibited by acidic residues.