Investigating the Counteracting Effect of Trehalose on Urea-Induced Protein Denaturation Using Molecular Dynamics Simulation

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.'

Temperature dependent molecular dynamics simulation study to understand the stabilizing effect of NADES on the protein β-Lactoglobulin

The thermal stability of a protein is an important concern for its practical application in food processing industries. In this study, we have carried out classical molecular dynamics simulations to systematically investigate the effect of NADES (natural deep eutectic solvent) on the stabilization of the protein β-Lactoglobulin (BLG) at different temperatures. This study sheds light on the very aspects of NADES composed of betaine and sorbitol on the stability of the protein. NADES provides better stability to the protein up to a temperature of 400 K than in water. It is observed that the protein starts to unfold above temperature 400 K in spite of the presence of NADES which is quiet evident from the root mean square deviation (RMSD) and radius of gyration (Rg) plots. The decreasing average solvent accessible surface area (SASA) values and increasing intra-protein hydrogen bonds indicate better stability of the protein in NADES medium than in water at temperatures 300 K and 400 K. At high temperatures viz. 450 K and 500 K the number and distribution of solvent species (betaine and sorbitol) around the protein surface show an increment that are evident from the calculations of solvation shell, radial and spatial distribution functions. Increased number of betaine molecules that interact with the protein through electrostatic interaction may lead to destabilization of the protein at these temperatures. This study suggests that NADES could be used as an ideal medium for thermal stability of the protein BLG up to a temperature of 400 K. Beyond this temperature, NADES used for this study fails to exert stabilization effect on the protein.

Counteraction Effects of Ammonium-Based Ionic Liquids on Urea-Induced Denaturation of α-Lactalbumin: A Comprehensive Molecular Simulation Study

Ionic liquids (ILs) are known to stabilize protein conformations in aqueous medium. Importantly, ILs can also act as refolding additives in urea-driven denaturation of proteins. However, despite the importance of the problem, detailed microscopic understanding of the counteraction effects of ILs on urea-induced protein denaturation remains elusive. In this work, atomistic molecular dynamics (MD) simulations of the protein α-lactalbumin have been carried out in pure aqueous medium, in 8 M binary urea-water solution and in ternary urea-IL-water solutions containing ammonium-based ethyl ammonium acetate (EAA) as the IL at different concentrations (1-4 M). Attempts have been made to quantify detailed molecular-level understanding of the origin behind the counteraction effects of the IL on urea-induced partial unfolding of the protein. The calculations revealed significant conformational changes of the protein with multiple free energy minima due to its partial unfolding in binary urea-water solution. The counteraction effect of the IL was evident from the enhanced structural rigidity of the protein with propensity to transform into a single native free energy minimum state in ternary urea-IL-water solutions. Such an effect has been found to be associated with preferential direct binding of the IL components with the protein and simultaneous expulsion of urea from the interface, thereby providing additional stabilization of the protein in ternary solutions. Most importantly, modified rearrangement of the hydrogen bond network at the interface due to the formation of stronger protein-cation (PC) and protein-anion (PA) hydrogen bonds by breaking relatively weaker protein-urea (PU) and protein-water (PW) hydrogen bonds has been recognized as the microscopic origin behind the counteraction effects of EAA on urea-induced partial unfolding of the protein.

Stability of the Aβ42 Peptide in Mixed Solutions of Denaturants and Proline

Amyloid β-peptide (Aβ) is responsible for the neuronal damage and death of a patient with Alzheimer's disease (AD). Aβ42 oligomeric forms are dominant neurotoxins and are related to neurodegeneration. Their different forms are related to various pathological conditions in the brain. We investigated Aβ42 peptides in different environments of proline, urea, and GdmCl solutions (in pure and mixed binary forms) through atomistic molecular dynamics simulations. Preferential exclusion from the protein surface and facile formation of a large number of weak molecular interactions are the driving forces for the osmolyte's action. We have focused on these interactions between peptide monomers and pure/mixed osmolytes and denaturants. Urea, as usual, denatures the peptide strongly compared to the GdmCl by accumulation around the peptide. GdmCl shows lesser build-up around protein in contrast to urea but is involved in destabilizing the salt bridge formation of Asp23 and Lys28. Proline as an osmolyte protects the peptide from aggregation when mixed with urea and GdmCl solutions. In mixed solutions of two denaturants and osmolyte plus denaturant, the peptide shows enhanced stability as compared to pure denaturant urea solution. The enhanced stability of peptides in proline may be attributed to its exclusion from the peptide surface and favoring salt bridge formation.

Conformational Features and Hydration Dynamics of Proteins in Cosolvents: A Perspective from Computational Approaches

The importance of solvent in stabilizing protein structures has long been recognized. Water is the common solvent for proteins, and hydration is elemental in governing protein stability, flexibility, and function through various interactions. The addition of small organic molecules known as cosolvents may deploy stabilization (folding) or destabilization (unfolding) effects on native protein conformations. Despite exhaustive literature, the molecular mechanism by which cosolvents regulate protein conformations and dynamics is controversial. Specifically, the cosolvent behavior has been unpredictable with the nature and concentrations that lead to protein stabilizing/destabilizing effects as it changes in water content near the vicinity of proteins. With the massive development of computational resources, advancement of computational methods, and the availability of numerous experimental techniques, various theoretical and computational studies of proteins in a mixture of solvents have been instigated. The growing interest in such studies has been to unravel the underlying mechanism of protein folding and cosolvent/solvent-protein interactions that have significant implications in biomedical and biotechnological applications. In this mini-review, apart from the brief overview of important theories and force-field model-based cosolvent effects on proteins, we present the current state of knowledge and recent advances in the field to describe cosolvent-guided conformational features of proteins and hydration dynamics from computational approaches. The mini-review further explains the mechanistic details of protein stability in various popularly used cosolvents, including limitations of present studies and future outlooks. The counteracting effects of cosolvent on the proteins in the mixture of stabilizing and destabilizing cosolvents are also presented and discussed.

Preparation of a novel regenerated silk fibroin-based hydrogel for extrusion bioprinting

Three-dimensional (3D) bioprinting technology, allowing rapid prototyping and personalized customization, has received much attention in recent years, while regenerated silk fibroin (RSF) has also been widely investigated for its excellent biocompatibility, processibility, and comprehensive mechanical properties. However, due to the difficulty in curing RSF aqueous solution and the tendency of conformational transition of RSF chains under shearing, it is rather complicated to fabricate RSF-based materials with high mechanical strength through extrusion bioprinting. To solve this problem, a printable hydrogel with thixotropy was prepared from regenerated silk fibroin with high-molecular-weight (HMWRSF) combined with a small amount of hydroxypropyl methylcellulose (HPMC) in urea containing aqueous solution. It was found that the introduction of urea could not only vary the solid content of the hydrogel to benefit the mechanical properties of the 3D-bioprinted pre-cured hydrogels or 3D-bioprinted sponges, but also expand the "printable window" of this system. Indeed, the printability and rheological properties could be modulated by varying the solid content, the heating time, the urea/HMWRSF weight ratio, etc. Moreover, the microstructure of nanospheres stacked in these lyophilized 3D-bioprinted sponges was interesting to observe, which indicated the existence of microhydrogels and both "the reversible network" and "the irreversible network" in this HMWRSF-based pre-cured hydrogel. Like other HMWRSF materials fabricated in other ways, these 3D-bioprinted HMWRSF-based sponges exhibited good cytocompatibility for dental pulp mesenchymal stem cells. This work may inspire the design of functional HMWRSF-based materials by regulating the relationship between structure and properties.

In silico studies of the human IAPP in the presence of osmolytes

The human islet amyloid polypeptide or amylin is secreted along with insulin by pancreatic islets. Under the drastic environmental conditions, amylin can aggregate to form amyloid fibrils. This amyloid plaque of hIAPP in the pancreatic cells is the cause of type II diabetes. Early stages of amylin aggregates are more cytotoxic than the matured fibrils. Here, we have used the all-atom molecular dynamic simulation to see the effect of water, TMAO, urea and urea/TMAO having ratio 2:1 of different concentrations on the amylin protein. Our study suggest that the amylin protein forms β-sheets in its monomeric form and may cause the aggregation of protein through the residue 13-17 and the C-terminal region. α-Helical content of protein increases with an increase in TMAO concentration by decreasing the SASA value of protein, increase in intramolecular hydrogen bonds and on making the short-range hydrophobic interactions. Electrostatic potential surfaces show that hydrophobic groups are buried and normalised configurational entropy of backbone, and side-chain atoms is lesser in the presence of TMAO, whereas opposite behaviour is obtained in the case of urea. Counteraction effect of TMAO using Kast model towards urea is also observed in ternary solution of urea/TMAO.

Effects of Ionic Liquids on Aqueous Urea Solutions: Insights into the Ionic Liquid-Assisted Protein Renaturation

Ionic liquids (ILs) are designer solvents that find wide applications in various areas. Recently, ILs have been shown to induce the refolding of certain proteins that were previously denatured under the treatment of urea. A molecular-level understanding of the counteracting mechanism of ILs on urea-induced protein denaturation remains elusive. In this study, we employ atomistic molecular dynamics simulations to investigate the ternary urea-water-IL solution in comparison to the aqueous urea solution to understand how the presence of ILs can modulate the structure, energetics, and dynamics of urea-water solutions. Our results show that the ions of the IL used, ethylammonium nitrate (EAN), interact strongly with urea and disrupt the urea aggregates that were known to stabilize the unfolded state of the proteins. Results also suggest a disruption in urea-water interaction that releases more free water molecules in solution. We subsequently strengthened these findings by simulating a model peptide in the absence and presence of EAN, which showed broken versus intact secondary structure in urea solution. Analyses show that these changes were accomplished by the added IL, which enforced a gradual displacement of urea from the peptide surface by water. We propose that the ILs facilitate protein renaturation by breaking down the urea aggregates and increasing the amount of free water molecules around the protein.

New insight into the molecular mechanism of the trehalose effect on urate oxidase stability

Urate oxidase (EC 1.7.3.3) is a key enzyme in the purine metabolism which is applied in the treatment of gout and also, as a diagnostic reagent for uric acid detection. In the current study, the trehalose (TRE) effects as an additive on the structural stability and function of uricase were investigated. For recombinant expression of UOX in E. coli BL21 cells, firstly the coding sequence was subcloned into the pET-28a vector and after induction with IPTG, the recombinant UOX was purified by affinity chromatography using a Ni-NTA agarose column. To specify the trehalose effects on the urate oxidase (UOX) structure, optimum pH, optimum temperature, kinetic and thermodynamic parameters and also, the intrinsic fluorescence of UOX in the absence and presence of trehalose were examined. The UOX half-life is 24.32 min at 40 °C, whereas the UOX-TRE has a higher half-life (32.09 min) at this temperature. Generally, our findings confirm that trehalose has a protective effect on the enzyme structure. Optimum pH and temperature were 9 and 25 °C, respectively for both the naked and treated enzymes and their activity retained 42.18 and 64.80%, respectively after 48 h of incubation at room temperature. Also, theoretical results indicate that the random coil of the enzyme was converted to α-helix and β-sheet in the presence of trehalose which may preserve the integrity of the active site of the enzyme and increased the enzymatic activity. The MD simulation results indicated greater stability of the uricase structure in the presence of trehalose.Communicated by Ramaswamy H. Sarma.

Surface hydration and preferential interaction directs the charged amino acids-induced changes in protein stability

In the present study, we investigate the interaction of amino acid osmolytes, Arg, Lys, Asp and Glu, and a denaturant, guanidinium chloride (Gdm) with proteins. To achieve this, molecular dynamics (MD) simulation of RNase A and α-lactalbumin was performed in the presence of three charged amino acids Arg, Lys, and Asp and the molecular mechanism of amino acid-induced (de)stabilization of the proteins was examined by combining with our earlier report on Glu. As Arg has the side chain similar to that of Gdm and destabilizes the proteins, MD simulation was carried out in the presence of Gdm as well. Radial distribution function and hydration fraction around the protein surface reveals that preferential hydration increases upon the addition of any of the cosolvent; however, the extent of increase is more in the presence of stabilizing cosolvents (stAAs: Lys, Asp and Glu) compared to destabilizing cosolvents (Arg and Gdm). Moreover, the preferential interaction of Arg and Gdm with the proteins is higher than that of stAAs. Residue-level interaction analysis suggests that stAAs preferably interacts with charged amino acids of the proteins whereas Arg and Gdm interactions could be found on almost all the surface exposed residues which might provide higher preferential interaction for these residues. From the results, we propose that the net outcome of preferential hydration versus preferential interaction of the amino acids might determine their effect on the stability of proteins.

Synergistic Catalytic Synthesis of Gemini Lipoamino Acids Based on Multiple Hydrogen-Bonding Interactions in Natural Deep Eutectic Solvents-Enzyme System

In the previous studies, gemini lipoamino acids (GLAA) were always synthesized by complex multistep organic synthesis, which involved a large number of byproducts and organic solvents. To develop a straightforward, efficient, and renewable synthesis strategy for GLAA, in this study, a type of novel green solvents, natural deep eutectic solvents (NADESs), were adopted as the solvents for these reactions. Five commercial enzymes were involved in the enzyme screening section, and Candida antarctica lipase B (CALB) tended to have the best performance in NADESs systems. The optimization procedure was performed using the Taguchi crossed array method and the highest yield of GLAA (59.14 ± 0.51%) was obtained in choline chloride-glycerol (C-Gly). The purification procedure was carried out with ethyl acetate and water, and the isolate yield ranged from 86.31 ± 2.36 to 91.34 ± 2.26%. With 10 times recycling, the yield of GLAA in C-Gly decreased from 59.14 ± 0.51 to 51.31 ± 0.68%. Interestingly, a synergistic effect of CALB and NADESs was found in the enzymatic synthesis of GLAA, which can be attributed to fatty acids being activated by chloride ions via hydrogen-bonding interactions and resulting in an enhancement in its electron-attracting ability.

ATP Controls the Aggregation of Aβ16-22 Peptides

The oligomerization of Aβ_16-22 peptide, which is the hydrophobic core region of full-length Aβ_1-42, causes Alzheimer's disease (AD). This progressive neurodegenerative disease affects over 44 million people worldwide. However, very few synthesized drug molecules are available to inhibit the aggregation of Aβ. Recently, experimental studies have shown that the biological ATP molecule prevents Aβ fibrillation at the millimolar scale; however, the significance of ATP molecules on Aβ fibrillation and the mechanism behind it remain elusive. We have carried out a total of 7.5 μs extensive all-atom molecular dynamics and 8.82 μs of umbrella sampling in explicit water using AMBER14SB, AMBER99SB-ILDN, and AMBER-FB15 force fields for Aβ_16-22 peptide, to investigate the role of ATP on the disruption of Aβ_16-22 prefibrils. From various analyses, such as secondary structure analysis, residue-wise contact map, SASA, and interaction energies, we have observed that, in the presence of ATP, the aggregation of Aβ_16-22 peptide is very unfavorable. Moreover, the biological molecule ATP interacts with the Aβ_16-22 peptide via hydrogen bonding, π-π stacking, and NH-π interactions which, ultimately, prevent the aggregation of Aβ_16-22 peptide. Hence, we assume that the deficiency of ATP may cause Alzheimer's disease (AD).

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.

Trehalose mediated stabilisation of cellobiase aggregates from the filamentous fungus Penicillium chrysogenum

Extracellular fungal cellobiases develop large stable aggregates by reversible concentration driven interaction. In-vitro addition of trehalose resulted in bigger cellobiase assemblies with increased stability against heat and dilution induced dissociation. In presence of 0.1 M trehalose, the size of aggregates increased from 344 nm to 494 nm. The increase in size was also observed in zymography of cellobiase. Activation energy of the trehalose stabilised enzyme (Ea = 220.9 kJ/mol) as compared to control (Ea = 257.734 kJ/mol), suggested enhanced thermostability and also showed increased resistance to chaotropes. Purified cellobiase was found to contain 196.27 μg of sugar/μg of protein. It was proposed that presence of glycan on protein's surface impedes and delays trehalose docking. Consequently, self-association of cellobiase preceded coating by trehalose leading to stabilisation of bigger cellobiase aggregates. In unison with the hypothesis, ribosylated BSA failed to get compacted by trehalose and developed into bigger aggregates with average size increasing from 210 nm to 328 nm. Wheat Germ Lectin, in presence of trehalose, showed higher molecular weight assemblies in DLS, native-PAGE and fluorescence anisotropy. This is the first report of cross-linking independent stabilisation of purified fungal glycosidases providing important insights towards understanding the aggregation and stability of glycated proteins.

Trehalose Mediated Inhibition of Lactate Dehydrogenase from Rabbit Muscle. The Application of Kramers' Theory in Enzyme Catalysis

Lactate dehydrogenase (LDH) catalyzes the reduction of pyruvate to lactate by using NADH. LDH kinetics has been proposed to be dependent on the dynamics of a loop over the active site. Kramers' theory has been useful in the study of enzyme catalysis dependent on large structural dynamics. In this work, LDH kinetics was studied in the presence of trehalose and at different temperatures. In the absence of trehalose, temperature increase raised exponentially the LDH V_max and revealed a sigmoid transition of K_m toward a low-affinity state similar to protein unfolding. Notably, LDH V_max diminished when in the presence of trehalose, while pyruvate affinity increased and the temperature-mediated binding site transition was hindered. The effect of trehalose on k_cat was viscosity dependent as described by Kramers' theory since V_max correlated inversely with the viscosity of the medium. As a result, activation energy ( E_a) for pyruvate reduction was dramatically increased by trehalose presence. This work provides experimental evidence that the dynamics of a structural component in LDH is essential for catalysis, i.e., the closing of the loop on the active site. While the trehalose mediated-increased of pyruvate affinity is proposed to be due to the compaction and/or increase of structural order at the binding site.

Glutamate Induced Thermal Equilibrium Intermediate and Counteracting Effect on Chemical Denaturation of Proteins

When organisms are subjected to stress conditions, one of their adaptive responses is accumulation of small organic molecules called osmolytes. These osmolytes affect the structure and stability of the biological macromolecules including proteins. The present study examines the effect of a negatively charged amino acid osmolyte, glutamate (Glu), on two model proteins, ribonuclease A (RNase A) and α-lactalbumin (α-LA), which have positive and negative surface charges at pH 7, respectively. These proteins follow two-state unfolding transitions during both heat and chemical induced denaturation processes. The addition of Glu stabilizes the proteins against temperature and induces an early equilibrium intermediate during unfolding. The stability is found to be enthalpy-driven, and the free energy of stabilization is more for α-LA compared to RNase A. The decrease in the partial molar volume and compressibility of both of the proteins in the presence of Glu suggests that the proteins attain a more compact state through surface hydration which could provide a more stable conformation. This is also supported by molecule dynamic simulation studies which demonstrate that the water density around the proteins is increased upon the addition of Glu. Further, the intermediates could be completely destabilized by lower concentrations (∼0.5 M) of guanidinium chloride and salt. However, urea subverts the Glu-induced intermediate formed by α-LA, whereas it only slightly destabilizes in the case of RNase A which has a positive surface charge and could possess charge-charge interactions with Glu. This suggests that, apart from hydration, columbic interactions might also contribute to the stability of the intermediate. Gdm-induced denaturation of RNase A and α-LA in the absence and the presence of Glu at different temperatures was carried out. These results also show the Glu-induced stabilization of both of the proteins; however, all of the unfolding transitions followed two-state transitions during chemical denaturation. The extent of stability exerted by Glu is higher for RNase A at higher temperature, whereas it provides more stability for α-LA at lower temperature. Thus, the experiments indicate that Glu induces a thermal equilibrium intermediate and increases the thermodynamic stability of proteins irrespective of their surface charges. The extent of stability varies between the proteins in a temperature-dependent manner.

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.

Can an ammonium-based room temperature ionic liquid counteract the urea-induced denaturation of a small peptide?

The folding/unfolding equilibrium of proteins in aqueous medium can be altered by adding small organic molecules generally termed as co-solvents.

Molecular dynamics approach to understand the denaturing effect of a millimolar concentration of dodine on a λ-repressor and counteraction by trehalose

Here, we use a molecular dynamics approach to calculate the spatial distribution function of the ternary water–dodine–trehalose (1.0 M) system.

Ammonium based stabilizers effectively counteract urea-induced denaturation in a small protein: insights from molecular dynamics simulations

Room temperature ionic liquids (IL) and deep eutectic solvents (DES) are known to aid the conformational stability and activity of proteins and enzymes in aqueous solutions.

DMSO enhanced conformational switch of an interfacial enzyme

Interfacial proteins function in unique heterogeneous solvent environments, such as water-oil interfaces. One important example is microbial lipase, which is activated in an oil-water emulsion phase and has many important enzymatic functions. A unique aprotic dipolar organic solvent, dimethyl sulfoxide (DMSO), has been shown to increase the activity of lipases, but the mechanism behind this enhancement is still unknown. Here, all-atom molecular dynamics simulations of lipase in a binary solution were performed to examine the effects of DMSO on the dynamics of the gating mechanism. The amphiphilic α5 region of the lipase was a focal point for the analysis, since the structural ordering of α5 has been shown to be important for gating under other perturbations. Compared to the closed-gorge ensemble in an aqueous environment, the conformational ensemble shifts towards open-gorge structures in the presence of DMSO solvents. Increased width of the access channel is particularly prevalent in 45% and 60% DMSO concentrations (w/w). As the amount of DMSO increases, the α5 region of the lipase becomes more α-helical, as we previously observed in studies that address water-oil interfacial and high pressure activation. We believe that the structural ordering of α5 plays an essential role on gating and lipase activity.

Proteins in binary solvents

Proteins in living organisms exist in complex aqueous solutions or embedded in membranes. In solution, proteins are surrounded by a tightly bound hydration layer, which is more ordered and less mobile than bulk water. As a consequence, water plays a major role in controlling protein structure stability, conformational flexibility, dynamics, and functionality, but it also appears that protein surface regulates the structuring of the surrounding water. The presence of cosolvents can modify the hydration layer characteristics and then the whole protein structural and dynamical properties. Because cytoplasm or biological liquids are complex solutions, the knowledge of the solvation shell characteristics in mixed solvents should be considered as a crucial step in describing biological processes at molecular level. This review reports on recent studies on the structural and thermodynamic properties of model proteins dissolved in binary solvent mixtures by small-angle neutron scattering (SANS) and differential scanning microcalorimetry (DSC) techniques. We will show that contrast variation SANS experiments allow to acquire a direct knowledge of both protein structure and protein solvation shell (in terms of low-resolution shape and solvent/cosolvent composition), while DSC experiments provide information on all the relevant thermodynamic properties. We will focus on two main points. First, an extended description of the thermodynamic model used to define the equilibria between water and cosolvent molecules in the protein solvation shell will be presented. Second, the determination of the peculiar characteristics of the protein solvation layer, which will be illustrated by considering different systems. As a conclusion, we will show that the investigation of structure and thermodynamics of proteins in binary aqueous mixtures is an important way to understand the role of hydration in protein stability and activity.

Counteraction of Trehalose on N, N-Dimethylformamide-Induced Candida rugosa Lipase Denaturation: Spectroscopic Insight and Molecular Dynamic Simulation

Candida rugosa lipase (CRL) has been widely used as a biocatalyst for non-aqueous synthesis in biotechnological applications, which, however, often suffers significant loss of activity in organic solvent. Experimental results show that trehalose could actively counteract the organic-solvent-induced protein denaturation, while the molecular mechanisms still don’t unclear. Herein, CRL was used as a model enzyme to explore the effects of trehalose on the retention of enzymatic activity upon incubation in N,N-dimethylformamide (DMF). Results showed that both catalytic activity and conformation changes of CRL influenced by DMF solvent were inhibited by trehalose in a dose-dependent fashion. The simulations further indicated that the CRL protein unfolded in binary DMF solution, but retained the native state in the ternary DMF/trehalose system. Trehalose as the second osmolyte added into binary DMF solution decreased DMF-CRL hydrogen bonds efficiently, whereas increased the intermolecular hydrogen bondings between DMF and trehalose. Thus, the origin of its denaturing effects of DMF on protein is thought to be due to the preferential exclusion of trehalose as well as the intermolecular hydrogen bondings between trehalose and DMF. These findings suggest that trehalose protect the CRL protein from DMF-induced unfolding via both indirect and direct interactions.

Effects of monohydric alcohols and polyols on the thermal stability of a protein

The thermal stability of a protein is lowered by the addition of a monohydric alcohol, and this effect becomes larger as the size of hydrophobic group in an alcohol molecule increases. By contrast, it is enhanced by the addition of a polyol possessing two or more hydroxyl groups per molecule, and this effect becomes larger as the number of hydroxyl groups increases. Here, we show that all of these experimental observations can be reproduced even in a quantitative sense by rigid-body models focused on the entropic effect originating from the translational displacement of solvent molecules. The solvent is either pure water or water-cosolvent solution. Three monohydric alcohols and five polyols are considered as cosolvents. In the rigid-body models, a protein is a fused hard spheres accounting for the polyatomic structure in the atomic detail, and the solvent is formed by hard spheres or a binary mixture of hard spheres with different diameters. The effective diameter of cosolvent molecules and the packing fractions of water and cosolvent, which are crucially important parameters, are carefully estimated using the experimental data of properties such as the density of solid crystal of cosolvent, parameters in the pertinent cosolvent-cosolvent interaction potential, and density of water-cosolvent solution. We employ the morphometric approach combined with the integral equation theory, which is best suited to the physical interpretation of the calculation result. It is argued that the degree of solvent crowding in the bulk is the key factor. When it is made more serious by the cosolvent addition, the solvent-entropy gain upon protein folding is magnified, leading to the enhanced thermal stability. When it is made less serious, the opposite is true. The mechanism of the effects of monohydric alcohols and polyols is physically the same as that of sugars. However, when the rigid-body models are employed for the effect of urea, its addition is predicted to enhance the thermal stability, which conflicts with the experimental fact. We then propose, as two essential factors, not only the solvent-entropy gain but also the loss of protein-solvent interaction energy upon protein folding. The competition of changes in these two factors induced by the cosolvent addition determines the thermal-stability change.