Cosolvent effects on protein stability

Effect of TMAO and Urea on Dimers and Tetramers of Amyloidogenic Heptapeptides (23FGAILSS29)

Human islet amyloid polypeptide (hIAPP) (1-37) is an intrinsically disordered protein that is released with insulin by β-cells found in the pancreas. Under certain environmental conditions, hIAPP can aggregate, which leads to β-cell death. FGAILSS (23-29) residues of the hIAPP protein form β sheets, which may be toxic species in type 2 diabetes (T2D) patients. All-atom molecular dynamics (MD) simulations have been used to analyze the effect of two distinct types of osmolytes trimethylamine N-oxide (TMAO) and urea on two and four FGAILSS heptapeptides. TMAO leads the individual peptide toward an extended conformation with a higher radius of gyration and favors the formation of antiparallel β-sheets with an increase in its concentration. However, urea mostly shows compaction of individual peptides except at 4.0 M in the case of a tetramer but does not show aggregation behavior as a whole. TMAO leads both the dimer and tetramer toward the native state with an increase in its concentration. Moreover, both the dimer and tetramer show irregular behavior in urea. The tetramer in 4.0 M urea shows the maximum fraction of native contacts due to the formation of antiparallel β-sheets. This formation of antiparallel β-sheets favors the aggregation of peptides. TMAO forms a smaller number of hydrogen bonds with peptides as compared to urea as the exclusion of TMAO and accumulation of urea around the peptides have occurred in the first solvation shell (FSS). Principal component analysis (PCA) results suggest that the minima in the free energy landscape (FEL) plot are homogeneous for a particular conformation in TMAO with smaller basins, while in urea, the dimer shows minima mostly for extended conformations. For a 4.0 M urea concentration, the tetramer shows the minimum for antiparallel β-sheets, which indicates the aggregation behavior of the tetramer, and for a higher concentration, it shows minima with wider basins of extended conformations.

The effects of cosolutes and crowding on the kinetics of protein condensate formation based on liquid-liquid phase separation: a pressure-jump relaxation study

Biomolecular assembly processes based on liquid-liquid phase separation (LLPS) are ubiquitous in the biological cell. To fully understand the role of LLPS in biological self-assembly, it is necessary to characterize also their kinetics of formation and dissolution. Here, we introduce the pressure-jump relaxation technique in concert with UV/Vis and FTIR spectroscopy as well as light microscopy to characterize the evolution of LLPS formation and dissolution in a time-dependent manner. As a model system undergoing LLPS we used the globular eye-lens protein γD-crystallin. As cosolutes and macromolecular crowding are known to affect the stability and dynamics of biomolecular condensates in cellulo, we extended our kinetic study by addressing also the impact of urea, the deep-sea osmolyte trimethylamine-N-oxide (TMAO) and a crowding agent on the transformation kinetics of the LLPS system. As a prerequisite for the kinetic studies, the phase diagram of γD-crystallin at the different solution conditions also had to be determined. The formation of the droplet phase was found to be a very rapid process and can be switched on and off on the 1-4 s timescale. Theoretical treatment using the Johnson-Mehl-Avrami-Kolmogorov model indicates that the LLPS proceeds via a diffusion-limited nucleation and growth mechanism at subcritical protein concentrations, a scenario which is also expected to prevail within biologically relevant crowded systems. Compared to the marked effect the cosolutes take on the stability of the LLPS region, their effect at biologically relevant concentrations on the phase transformation kinetics is very small, which might be a particular advantage in the cellular context, as a fast switching capability of the transition should not be compromised by the presence of cellular cosolutes.

Trehalose Inhibits the Heat-Induced Formation of the Amyloid-Like Structure of Soluble Proteins Isolated from Human Cataract Lens

The age-dependent loss of solubility and aggregation of crystallins constitute the pathological hallmarks of cataract. Several biochemical and biophysical factors are responsible for the reduction of crystallins' solubility and formation of irreversible protein aggregates, which display amyloid-like characteristics. The present study reports the heat-induced aggregation of soluble proteins isolated from human cataract lenses and the formation of amyloid-like structures. Exposure of protein at 55 °C for 4 h resulted in extensive (≈ 60%) protein aggregation. The heat-induced protein aggregates displayed substantial (≈ 20 nm) redshift in the wavelength of maximum absorption (λ_max) of Congo red (CR) and increase in Thioflavin T (ThT) fluorescence emission intensity, indicating the presence of amyloid-like structures in the heat-induced protein aggregates. Subsequently, the addition of trehalose resulted in substantial inhibition of heat-induced aggregation and the formation of amyloid-like structure. The ability of trehalose to inhibit the heat-induced aggregation was found to be linearly dependent upon its concentration used. The optimum effect was observed in the presence of 30-40% (w/v) trehalose where the aggregated was found to be reduced from 60 to 30%. The present study demonstrated the ability to trehalose to inhibit the protein aggregation and interfere with the formation of amyloid-like structures.

SDS-induced multi-stage unfolding of a small globular protein through different denatured states revealed by single-molecule fluorescence

Ionic surfactants such as sodium dodecyl sulfate (SDS) unfold proteins in a much more diverse yet effective way than chemical denaturants such as guanidium chloride (GdmCl). But how these unfolding processes compare on a molecular level is poorly understood. Here, we address this question by scrutinising the unfolding pathway of the globular protein S6 in SDS and GdmCl with single-molecule Förster resonance energy transfer (smFRET) spectroscopy. We show that the unfolding mechanism in SDS is strikingly different and convoluted in comparison to denaturation in GdmCl. In contrast to the reversible two-state unfolding behaviour in GdmCl characterised by kinetics on the timescale of seconds, SDS demonstrated not one, but four distinct regimes of interactions with S6, dependent on the surfactant concentration. At ≤1 mM SDS, S6 and surfactant molecules form quasi-micelles on a minute timescale; at millimolar [SDS], the protein denatures through an unfolded/denatured ensemble of highly heterogeneous states on a multi-second timescale; at tens of millimolar of SDS, the protein unfolds into a micelle-packed conformation on the second timescale; and >50 mM SDS, the protein unfolds with millisecond timescale dynamics. We propose a detailed model for multi-stage unfolding of S6 in SDS, which involves at least three different types of denatured states with different level of compactness and dynamics and a continually changing landscape of interactions between protein and surfactant. Our results highlight the great potential of single-molecule fluorescence as a direct probe of nanoscale protein structure and dynamics in chemically complex surfactant environments.

Study on effects of co-solvents on the structure of DhaA by molecular dynamics simulation

With the increasing application of enzymes in various research fields, the choices of co-solvents in enzymatic preparations which directly related to the catalytic activity have been attracted attention. Thus, researching on the stabilization or destabilization behaviors of enzymes in different solvents is extremely essential. In this study, the structural changes of DhaA in two typical aprotic co-solvents (acetonitrile and tetrahydrofuran) were firstly investigated by molecular dynamics (MD) simulation. The simulation results revealed the strong van der Waals force between co-solvents and DhaA which could induce the structural change of enzyme. Interestingly, the differences of molecular size and the electrostatic force with enzyme of two co-solvents led to quite different influences on DhaA. As for acetonitrile, solvent molecules could penetrate into the catalytic site of DhaA which promoted by the electrostatic interaction. On the contrary, tetrahydrofuran molecules were mainly distributed around the catalytic site due to the relative weak electrostatic interaction and steric resistance effect. It can be concluded that different co-solvent can affect the key domains, substrate pathway and catalytic pocket of DhaA.Communicated by Ramaswamy H. Sarma.

Heterologous expression, purification, and refolding of SRY protein: role of L-arginine as analyzed by simulation and practical study

Escherichia coli is a widely-used cell factory for recombinant protein production, nevertheless, high amount of produced protein is seen in aggregated form. The purpose of this study was to improve the solubility of recombinant bovine sex-determining region Y protein (rbSRY) by exploring the effect of temperature, inducer, and water-arginine mixed solvent. Codon-optimized rbSRY expressed in Rosetta-gami B (DE3) pLysS and purified by NI-NTA His-select affinity chromatography in the native and denaturing conditions. A three-dimensional model of SRY was built and studied through molecular dynamics simulations in water and in the presence of L-arginine as co-solvent. Results indicated the significant effects of temperature and IPTG concentration (P < 0.001) on the solubility of rbSRY. The binding activity of native, inclusion bodies and refolded fractions to anti-rbSRY monoclonal antibody were concentration-dependent (P < 0.001). Based on molecular modeling results, the propensity of fragments in the N-terminal domain to form β-sheet and the relative instability of α-helices in terminal domains are the probable reasons for the high aggregation potential of SRY, which are mitigated in the presence of L-arginine. Altogether, our rbSRY protein was properly produced and applying appropriate culture conditions could help enhance its solubility, refold inclusion bodies, and improve its activity upon refolding.

Conformational dynamics of superoxide dismutase (SOD1) in osmolytes: a molecular dynamics simulation study

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease caused by the misfolding of Cu, Zn superoxide dismutase (SOD1). Several earlier studies have shown that monomeric apo SOD1 undergoes significant local unfolding dynamics and is the predecessor for aggregation. Here, we have employed atomistic molecular dynamics (MD) simulations to study the structure and dynamics of monomeric apo and holo SOD1 in water, aqueous urea and aqueous urea-TMAO (trimethylamine oxide) solutions. Loop IV (zinc-binding loop) and loop VII (electrostatic loop) of holo SOD1 are considered as functionally important loops as they are responsible for the structural stability of holo SOD1. We found larger local unfolding of loop IV and VII of apo SOD1 as compared to holo SOD1 in water. Urea induced more unfolding in holo SOD1 than apo SOD1, whereas the stabilization of both the form of SOD1 was observed in ternary solution (i.e. water/urea/TMAO solution) but the extent of stabilization was higher in holo SOD1 than apo SOD1. The partially unfolded structures of apo SOD1 in water, urea and holo SOD1 in urea were identified by the exposure of the hydrophobic cores, which are highly dynamic and these may be the initial events of aggregation in SOD1. Our simulation studies support the formation of aggregates by means of the local unfolding of monomeric apo SOD1 as compared to holo SOD1 in water.

Discriminating changes in protein structure using tyrosine conjugation

Chemical modification of proteins has been crucial in engineering protein-based therapies, targeted biopharmaceutics, molecular probes, and biomaterials. Here, we explore the use of a conjugation-based approach to sense alternative conformational states in proteins. Tyrosine has both hydrophobic and hydrophilic qualities, thus allowing it to be positioned at protein surfaces, or binding interfaces, or to be buried within a protein. Tyrosine can be conjugated with 4-phenyl-3H-1,2,4-triazole-3,5(4H)-dione (PTAD). We hypothesized that individual protein conformations could be distinguished by labeling tyrosine residues in the protein with PTAD. We conjugated tyrosine residues in a well-folded protein, bovine serum albumin (BSA), and quantified labeled tyrosine with liquid chromatography with tandem mass spectrometry. We applied this approach to alternative conformations of BSA produced in the presence of urea. The amount of PTAD labeling was found to relate to the depth of each tyrosine relative to the protein surface. This study demonstrates a new use of tyrosine conjugation using PTAD as an analytic tool able to distinguish the conformational states of a protein.

Influence of TMAO as co-solvent on the gelation of silica-PNIPAm core-shell nanogels at intermediate volume fractions

We study the structure and dynamics of poly(N-isopropylacrylamide) (PNIPAm) core-shell nanogels dispersed in aqueous trimethylamine N-oxide (TMAO) solutions by means of small-angle X-ray scattering and X-ray photon correlation spectroscopy (XPCS). Upon increasing the temperature above the lower critical solution temperature of PNIPAm at 33 °C, a colloidal gel is formed as identified by an increase of I(q) at small q as well as a slowing down of sample dynamics by various orders of magnitude. With increasing TMAO concentration the gelation transition shifts linearly to lower temperatures. Above a TMAO concentration of approximately 0.40 mol/L corresponding to a 1 : 1 ratio of TMAO and NIPAm groups, collapsed PNIPAm states are found for all temperatures without any gelation transition. This suggests that reduction of PNIPAm-water hydrogen bonds due to the presence of TMAO results in a stabilisation of the collapsed PNIPAm state and suppresses gelation of the nanogel.

Insights into the Sensitivity of Arginine Concentration to Preserve the Folded Form of Insulin Monomer under Thermal Stress

Arginine, although popularly known as aggregation suppressor additive, has been found to quench proteins' structure and function by destabilizing their conformations. Driven by such controversial evidence, in this work we performed a series of atomistic molecular dynamics simulations of insulin monomer, a biologically active hormone protein, in arginine solution of varying concentrations (0.5, 1, and 2 M) at ambient and elevated temperature (400 K) to explore the arginine concentration driven structure-based stability of the protein. Our study reveals that the flexibility of the protein's structure is dependent on the arginine concentration, and among all the used solutions, 2 M arginine, a "neutral crowder" that mimics the cellular environment, can preserve the native folded form of the protein at ambient temperature in an excellent manner. Further, while the protein unfolds at 400 K in pure water, this solution worked satisfactorily to preserve the protein's folded conformation more firmly than the other solutions. The replica-exchange MD of insulin in 2 M arginine solution further supports the fact. In this aspect an important issue in molecular pharmacology is to identify and recognize the physical origin of the stability of a protein, i.e, in this case, how arginine directs the conformational flexibility of the protein and preserves its native folded form. We identified that the exclusion of arginine from the protein surface increases the local structuration of water around the protein, thereby preserving its "biological water" layer, and makes the protein more hydrated at 2 M concentration as compared to the other arginine solutions. Additionally, our microscopic investigation on the interactions of the protein-solvation layer revealed that the structural heterogeneity of the protein surface, arising from the differential physicochemical nature of the amino acid residues, controls the favorable formation of sluggish water-arginine mixed solvation layer at higher arginine concentration that helps the protein to maintain its structural rigidity. Importantly, apart from the protein-solvent hydrogen-bonding interactions, the anion-pi interactions, established between the carboxyl group of arginine and the aromatic amino acid residues of insulin, were recognized to facilitate the protein to maintain its native folded form at the experimental temperatures.

Submillisecond Freezing Permits Cryoprotectant-Free EPR Double Electron-Electron Resonance Spectroscopy

Double electron-electron resonance (DEER) EPR spectroscopy is a powerful method for obtaining distance distributions between pairs of engineered nitroxide spin-labels in proteins and other biological macromolecules. These measurements require the use of cryogenic temperatures (77 K or less) to prolong the phase memory relaxation time (T_m ) sufficiently to enable detection of a DEER echo curve. Generally, a cryoprotectant such as glycerol is added to protein samples to facilitate glass formation and avoid protein clustering (which can result in a large decrease in T_m ) during relatively slow flash freezing in liquid N₂ . However, cryoprotectants are osmolytes and can influence protein folding/unfolding equilibria, as well as species populations in weak multimeric systems. Here we show that submillisecond rapid freezing, achieved by high velocity spraying of the sample onto a rapidly spinning, liquid nitrogen cooled copper disc obviates the requirement for cryoprotectants and permits high quality DEER data to be obtained in absence of glycerol. We demonstrate this approach on five different protein systems: protein A, the metastable drkN SH3 domain, urea-unfolded drkN SH3, HIV-1 reverse transcriptase, and the transmembrane domain of HIV-1 gp41 in lipid bicelles.

Mitigating effects of osmolytes on the interactions between nanoparticles and supported lipid bilayer

To maintain osmotic balance, cells usually produce neutral solutes (i.e., osmolytes), together with charged species to cope with salinity stress. Osmolytes are known to be important in stabilizing/destabilizing macromolecules (e.g., proteins) via depletion /accumulation around their surfaces. To better understand the physiological fate of nanoparticles (NPs), we investigated the effect of osmolytes [(urea and trimethylamine N-oxide (TMAO)] and specific anions (NO₃^- and F^-) on the interactions between NPs and supported lipid bilayers (SLBs). Carboxylated polystyrene NPs (60 nm) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) were chosen as model NPs and lipid. Quartz crystal microbalance with dissipation monitoring (QCM-D) was used to quantify NP deposition dynamics. Microscale thermophoresis (MST) was used to characterize the affinity between DOPC vesicles (or NPs) and osmolytes. Our results show that osmolytes are capable of protecting SLBs from NP-induced disruption. Upon NP deposition onto supported vesicle layers (SVLs), the leakage of encapsulated dyes decreased with the addition of osmolytes. The combination of kosmotropes (TMAO and F^-) are more efficient than that of chaotropes (urea and NO₃^-) in weakening the hydrophobic interaction between NPs and SLBs by preferential binding to NPs and/or SLBs.

Quantitative Interpretation of Solvent Paramagnetic Relaxation for Probing Protein-Cosolute Interactions

Protein-small cosolute molecule interactions are ubiquitous and known to modulate the solubility, stability, and function of many proteins. Characterization of such transient weak interactions at atomic resolution remains challenging. In this work, we develop a simple and practical NMR method for extracting both energetic and dynamic information on protein-cosolute interactions from solvent paramagnetic relaxation enhancement (sPRE) measurements. Our procedure is based on an approximate (non-Lorentzian) spectral density that behaves exactly at both high and low frequencies. This spectral density contains two parameters, one global related to the translational diffusion coefficient of the paramagnetic cosolute, and the other residue specific. These parameters can be readily determined from sPRE data, and then used to calculate analytically a concentration normalized equilibrium average of the interspin distance, ⟨r^-6⟩_norm, and an effective correlation time, τ_C, that provide measures of the energetics and dynamics of the interaction at atomic resolution. We compare our approach with existing ones, and demonstrate the utility of our method using experimental ¹H longitudinal and transverse sPRE data recorded on the protein ubiquitin in the presence of two different nitroxide radical cosolutes, at multiple static magnetic fields. The approach for analyzing sPRE data outlined here provides a powerful tool for deepening our understanding of extremely weak protein-cosolute interactions.

Protein adaptation to high hydrostatic pressure: Computational analysis of the structural proteome

Hydrostatic pressure has a vital role in the biological adaptation of the piezophiles, organisms that live under high hydrostatic pressure. However, the mechanisms by which piezophiles are able to adapt their proteins to high hydrostatic pressure is not well understood. One proposed hypothesis is that the volume changes of unfolding (ΔVTot ) for proteins from piezophiles is distinct from those of nonpiezophilic organisms. Since ΔVTot defines pressure dependence of stability, we performed a comprehensive computational analysis of this property for proteins from piezophilic and nonpiezophilic organisms. In addition, we experimentally measured the ΔVTot of acylphosphatases and thioredoxins belonging to piezophilic and nonpiezophilic organisms. Based on this analysis we concluded that there is no difference in ΔVTot for proteins from piezophilic and nonpiezophilic organisms. Finally, we put forward the hypothesis that increased concentrations of osmolytes can provide a systemic increase in pressure stability of proteins from piezophilic organisms and provide experimental thermodynamic evidence in support of this hypothesis.

Dynamical Model for the Counteracting Effects of Trimethylamine N-Oxide on Urea in Aqueous Solutions under Pressure

Of cosolutes found in living cells, urea denatures and trimethylamine N-oxide (TMAO) stabilizes proteins; furthermore, these effects cancel at a 2:1 ratio of urea to TMAO. Interestingly, cartilaginous fish use urea and TMAO as osmolytes at similar ratios at the ocean surface but with increasing fractions of TMAO at increasing depths. Here, molecular dynamics simulations of aqueous solutions with different urea:TMAO ratios show that the diffusion coefficients of water in the solutions vary with pressure if the urea:TMAO ratio is constant, but strikingly, they are almost pressure independent at the ratio found in these fish as a function of depth. This suggests that this ratio may be maintaining a homeostasis of water dynamics. In addition, diffusion is determined by hydrogen-bond lifetimes of the different species in the solution. Based on these observations, a dynamical model in terms of hydrogen-bond lifetimes is developed for the hydrogen bonding propensities of cosolutes and water in an aqueous solution to proteins. This model provides an explanation for both the counteracting effects of TMAO on urea denaturation and the depth-dependent urea:TMAO ratio found in cartilaginous fish.

Protein Preferential Solvation in Water:Glycerol Mixtures

For proteins in solvent mixtures, the relative abundances of each solvent in their solvation shell have a critical impact on their properties. Preferential solvation of a series of proteins in water-glycerol mixtures is studied here over a broad range of solvent compositions via classical molecular dynamics simulations. Our simulation results reveal that the differences between shell and bulk compositions exhibit dramatic changes with solvent composition, temperature, and protein nature. In contrast with the simple and widely used picture where glycerol is completely excluded from the protein interface, we show that for aqueous solutions with less than 50% glycerol in volume, protein solvation shells have approximately the same composition as the bulk solvent and proteins are in direct contact with glycerol. We further demonstrate that at high glycerol concentration, glycerol depletion from the solvation shell is due to an entropic factor arising from the reduced accessibility of bulky glycerol molecules in protein cavities. The resulting molecular picture is important to understand protein activity and cryopreservation in mixed aqueous solvents.

Stability of the chaperonin system GroEL-GroES under extreme environmental conditions

The chaperonin system GroEL-GroES is present in all kingdoms of life and rescues proteins from improper folding and aggregation upon internal and external stress conditions, including high temperatures and pressures. Here, we set out to explore the thermo- and piezostability of GroEL, GroES and the GroEL-GroES complex in the presence of cosolvents, nucleotides and salts employing quantitative FTIR spectroscopy and small-angle X-ray scattering. Owing to its high biological relevance and lack of data, our focus was especially on the effect of pressure on the chaperonin system. The experimental results reveal that the GroEL-GroES complex is remarkably temperature stable with an unfolding temperature beyond 70 °C, which can still be slightly increased by compatible cosolutes like TMAO. Conversely, the pressure stability of GroEL and hence the GroEL-GroES complex is rather limited and much less than that of monomeric proteins. Whereas GroES is pressure stable up to ∼5 kbar, GroEl and the GroEl-GroES complex undergo minor structural changes already beyond 1 kbar, which can be attributed to a dissociation-induced conformational drift. Quite unexpectedly, no significant unfolding of GroEL is observed even up to 10 kbar, however, i.e., the subunits themselves are very pressure stable. As for the physiological relevance, the structural integrity of the chaperonin system is retained in a relatively narrow pressure range, from about 1 to 1000 bar, which is just the pressure range encountered by life on Earth.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

Correction: High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation

Correction for ‘High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation’ by Hsuan-Lei Sung et al., Phys. Chem. Chem. Phys., 2020, DOI: 10.1039/d0cp01921f.

Distinctive behavior and two-dimensional vibrational dynamics of water molecules inside glycine solvation shell

We present a first principles molecular dynamics study of a deuterated aqueous solution of a single glycine moiety to explore the structure, dynamics, and two-dimensional infrared spectra of water molecules found in the solvation shell of glycine.

The multifaceted effects of DMSO and high hydrostatic pressure on the kinetic constants of hydrolysis reactions catalyzed by α-chymotrypsin

The cosolvent DMSO and high pressure have antagonistic effects on the kinetic constants of α-chymotrypsin-catalyzed hydrolysis reactions.

Hydration in aqueous osmolyte solutions: the case of TMAO and urea

X-ray Raman scattering spectroscopy and first principles simulations reveal details of the hydration and hydrogen-bond topology of trimethylamine N-oxide (TMAO) and urea in aqueous solutions.

Statistical field theory of ion–molecular solutions

Schematic representation of the multipolar molecule surrounded by salt ions in a dielectric solvent medium.

The shift in urea orientation at protein surfaces at low pH is compatible with a direct mechanism of protein denaturation

The protonation of acidic side-chains promotes a orientational shift of urea molecules, but only locally, with the interactions with other protein moieties being preserved.

Trimethylamine N-oxide (TMAO) resists the compression of water structure by magnesium perchlorate: terrestrial kosmotrope vs. Martian chaotrope

Neutron diffraction and computational modelling provide insight into water structure.

Ameliorating amyloid aggregation through osmolytes as a probable therapeutic molecule against Alzheimer's disease and type 2 diabetes

Large numbers of neurological and metabolic disorders occurring in humans are induced by the aberrant growth of aggregated or misfolded proteins.

Selection and screening strategies in directed evolution to improve protein stability

Protein stability is not only fundamental for experimental, industrial, and therapeutic applications, but is also the baseline for evolving novel protein functions. For decades, stability engineering armed with directed evolution has continued its rapid development and inevitably poses challenges. Generally, in directed evolution, establishing a reliable link between a genotype and any interpretable phenotype is more challenging than diversifying genetic libraries. Consequently, we set forth in a small picture to emphasize the screening or selection techniques in protein stability-directed evolution to secure the link. For a more systematic review, two main branches of these techniques, namely cellular or cell-free display and stability biosensors, are expounded with informative examples.

Efficacy of several additives to modulate the phase behavior of biomedical polymers: A comprehensive and comparative outlook

Several new classes of polymeric materials are being introduced with unique properties. Thermoresponsive polymers (TRPs) are one of the most fascinating and emerging class of biomaterials in biomedical research. The design of TRPs with good response to temperature and its ability to exhibit coil to globular transition behavior near to physiological temperature made them more promising materials in the field of biomaterials and biomedicines. Instead of numerous studies on TRPs, the mechanistic interplay among several additives and TRPs is still not understood clearly and completely. The lack of complete understanding of biomolecular interactions of various additives with TRPs is limiting their applications in interdisciplinary science as well as pharmaceutical industry. There is a great need to provide a collective and comprehensive information of various additives and their behavior on widely accepted biopolymers, TRPs such as poly(N-isopropylacrylamide) (PNIPAM), poly(vinyl methyl ether) (PVME), poly(N-vinylcaprolactum) (PVCL) and poly(ethylene glycol)-poly(propylene glycol)-poly(ethylene glycol) (PEG-PPG-PEG) in aqueous solution. Obviously, as the literature on the influence of various additives on TRPs is very vast, therefore we focus our review only on these four selected TRPs. Additives such as polyols, methylamines, surfactants and denaturants basically made the significant changes in water structure associated to polymer via their entropy variation which is the direct influence of their directly or indirectly binding abilities. Eventually, this review addresses a brief overview of the most recent literature of applications based phase behavior of four selected TRPs in response to external stimuli. The work enhances the knowledge for use of TRPs in the advanced development of drug delivery system and in many more pharmaceutical applications. These kinds of studies provide powerful impact in exploring the utility range of polymeric materials in various field of science.

Effect of Osmolytes on Conformational Behavior of Intrinsically Disordered Protein α-Synuclein

α-Synuclein is an intrinsically disordered protein whose function in a healthy brain is poorly understood. It is genetically and neuropathologically linked to Parkinson's disease (PD). PD is manifested after the accumulation of plaques of α-synuclein aggregates in the brain cells. Aggregates of α-synuclein are very toxic and lead to the disruption of cellular homeostasis and neuronal death. α-Synuclein can also contribute to disease propagation as it may exert noxious effects on neighboring cells. Understanding the mechanism of α-synuclein aggregation will facilitate the problem of dealing with neurodegenerative diseases in general and that of PD in particular. Here, we have used molecular dynamics simulations to investigate the behavior of α-synuclein at various temperatures and in different concentrations of urea and trimethyl amine oxide. The residue region from 61 to 95 of α-synuclein is experimentally known as amyloidogenic. In our study, we have identified some other regions, which also have the propensity to form an aggregate besides this known sequence. Urea being a denaturant interacts more with these regions of α-synuclein through hydrogen bond formation and inhibits the β-sheet formation, whereas trimethyl amine oxide itself does not interact much with the protein and stabilizes the protein by preferentially distributing water molecules on the surface of the protein.

Hofmeister Order of Anions on Protein Stability Originates from Lifshitz-van der Waals Dispersion Interaction with the Protein Phase

The mechanism underpinning the Hofmeister order of anions on protein stability and other physical and biological processes has been a mystery since its discovery in 1888. In the present study, we investigated electrostatic and Lifshitz-van der Waals (L-vdW) dispersion (electrodynamic) interactions between Hofmeister salts and four monomeric globular proteins. It is shown that structure-stabilizing salts exerted positive L-vdW pressure, whereas structure-destabilizing salts exerted negative L-vdW pressure on proteins. The relative order of the L-vdW pressure followed the Hofmeister series and it overshadowed the electrostatic pressure at high salt concentrations. The net change in the thermal denaturation temperature (ΔT_d) of proteins in 0.8 M Hofmeister salt solutions followed a linear relationship (r² > 0.8) with the net electrodynamic pressure regardless of the physicochemical differences between proteins. This study also revealed that segregation of anions into structure stabilizers and destabilizers depended on the dielectric susceptibility of the anion in the ultraviolet region: ions having absorbtion spectrum in the ultraviolet region (e.g., Cl^-, Br^-, I^-, and SCN^-) exerted a negative electrodynamic pressure, whereas those with absorbtion spectrum only in the infrared region, for example, SO₄^2-, exerted a positive electrodynamic pressure. The lack of ultraviolet absorption of SO₄^2- ions was because of quenching of ultraviolet radiation by water at below 170 nm.

Dynamics of TMAO and urea in the hydration shell of the protein SNase

Using all-atom molecular dynamics simulations of aqueous solutions of the globular protein SNase, the dynamic behavior of water molecules and cosolvents (trimethylamine-N-oxide (TMAO) and urea) in the hydration shell of the protein was studied for different solvent compositions. TMAO is a potent protein-stabilizing osmolyte, whereas urea is known to destabilize proteins. For molecules that are initially located in successive narrow layers at a given distance from the protein, the mean displacements and the distribution of displacements for short time intervals are calculated. For molecules that are initially located in solvation shells of a given thickness around the protein, the characteristic residence times in these shells are determined to characterize the dynamic behavior of the solvent molecules as a function of the distance to the protein. A combined consideration of these characteristics allows to reveal additional features of the dynamics of the cosolvents. It is shown that TMAO molecules leave the nearest vicinity of the protein faster than urea molecules, despite the fact that the mobility of TMAO molecules, measured by their mean displacements, is lower than that of urea. Moreover, we show that the rate of release of TMAO molecules from the hydration shell is lower in ternary (TMAO + urea + H₂O) solvent mixtures than in the binary ones. This is consistent with a recent observation that the fraction of TMAO near the protein decreases in the presence of urea. From the analysis of the decay of the number of particles initially located in the region of the first peak of the distribution function of solvent molecules around the protein, we estimated that about 20 water molecules and 6-7 urea molecules stay near the protein for more than 1000 ps.

Temperature, Hydrostatic Pressure, and Osmolyte Effects on Liquid-Liquid Phase Separation in Protein Condensates: Physical Chemistry and Biological Implications

Liquid-liquid phase separation (LLPS) of proteins and other biomolecules play a critical role in the organization of extracellular materials and membrane-less compartmentalization of intra-organismal spaces through the formation of condensates. Structural properties of such mesoscopic droplet-like states were studied by spectroscopy, microscopy, and other biophysical techniques. The temperature dependence of biomolecular LLPS has been studied extensively, indicating that phase-separated condensed states of proteins can be stabilized or destabilized by increasing temperature. In contrast, the physical and biological significance of hydrostatic pressure on LLPS is less appreciated. Summarized here are recent investigations of protein LLPS under pressures up to the kbar-regime. Strikingly, for the cases studied thus far, LLPSs of both globular proteins and intrinsically disordered proteins/regions are typically more sensitive to pressure than the folding of proteins, suggesting that organisms inhabiting the deep sea and sub-seafloor sediments, under pressures up to 1 kbar and beyond, have to mitigate this pressure-sensitivity to avoid unwanted destabilization of their functional biomolecular condensates. Interestingly, we found that trimethylamine-N-oxide (TMAO), an osmolyte upregulated in deep-sea fish, can significantly stabilize protein droplets under pressure, pointing to another adaptive advantage for increased TMAO concentrations in deep-sea organisms besides the osmolyte's stabilizing effect against protein unfolding. As life on Earth might have originated in the deep sea, pressure-dependent LLPS is pertinent to questions regarding prebiotic proto-cells. Herein, we offer a conceptual framework for rationalizing the recent experimental findings and present an outline of the basic thermodynamics of temperature-, pressure-, and osmolyte-dependent LLPS as well as a molecular-level statistical mechanics picture in terms of solvent-mediated interactions and void volumes.

Conformational behavior of a semiflexible dipolar chain with a variable relative size of charged groups via molecular dynamics simulations

The conformational behavior of an isolated semiflexible dipolar chain has been studied by molecular dynamics simulations. The dipolar chain was modeled as a backbone chain of charged beads, each containing an oppositely charged unit connected to it by a rigid spring. The main focus was on the effect of the backbone chain rigidity and the size of the charged groups on the morphology of the collapsed states of the chain formed in low-polar media where the electrostatic interactions are essential. It has been found that the stable globular conformations of the long chain of N = 256 backbone beads are a toroid and an elliptical globule. The macroscopic parameters (such as the radius of gyration and shape factors) as well as the local characteristics of these conformations (radial density distributions of ions, orientational correlations of chain segments, dipoles etc.) are studied depending on the chain stiffness. The regions of stability of a torus and an elliptical globule are found for the dipolar chains with variable dipole length and stiffness, which depend on the strength of electrostatic interactions. It has been shown that a size asymmetry of oppositely charged beads destabilizes globular states favoring elongated chain conformations. A coexistence of various metastable states was demonstrated for shorter chains of N = 128, 64, and 32.

Effects of in vivo conditions on amyloid aggregation

One of the grand challenges of biophysical chemistry is to understand the principles that govern protein misfolding and aggregation, which is a highly complex process that is sensitive to initial conditions, operates on a huge range of length- and timescales, and has products that range from protein dimers to macroscopic amyloid fibrils. Aberrant aggregation is associated with more than 25 diseases, which include Alzheimer's, Parkinson's, Huntington's, and type II diabetes. Amyloid aggregation has been extensively studied in the test tube, therefore under conditions that are far from physiological relevance. Hence, there is dire need to extend these investigations to in vivo conditions where amyloid formation is affected by a myriad of biochemical interactions. As a hallmark of neurodegenerative diseases, these interactions need to be understood in detail to develop novel therapeutic interventions, as millions of people globally suffer from neurodegenerative disorders and type II diabetes. The aim of this review is to document the progress in the research on amyloid formation from a physicochemical perspective with a special focus on the physiological factors influencing the aggregation of the amyloid-β peptide, the islet amyloid polypeptide, α-synuclein, and the hungingtin protein.

Using a FRET Library with Multiple Probe Pairs To Drive Monte Carlo Simulations of α-Synuclein

We describe a strategy for experimentally-constraining computational simulations of intrinsically disordered proteins (IDPs), using α-synuclein, an IDP with a central role in Parkinson's disease pathology, as an example. Previously, data from single-molecule Förster Resonance Energy Transfer (FRET) experiments have been effectively utilized to generate experimentally constrained computational models of IDPs. However, the fluorophores required for single-molecule FRET experiments are not amenable to the study of short-range (<30 Å) interactions. Using ensemble FRET measurements allows one to acquire data from probes with multiple distance ranges, which can be used to constrain Monte Carlo simulations in PyRosetta. To appropriately employ ensemble FRET data as constraints, we optimized the shape and weight of constraining potentials to afford ensembles of structures that are consistent with experimental data. We also used this approach to examine the structure of α-synuclein in the presence of the compacting osmolyte trimethylamine-N-oxide. Despite significant compaction imparted by 2 M trimethylamine-N-oxide, the underlying ensemble of α-synuclein remains largely disordered and capable of aggregation, also in agreement with experimental data. These proof-of-concept experiments demonstrate that our modeling protocol enables one to efficiently generate experimentally constrained models of IDPs that incorporate atomic-scale detail, allowing one to study an IDP under a variety of conditions.