Denaturant m values and heat capacity changes: relation to changes in accessible surface areas of protein unfolding

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

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

Thermodynamic origin of α-helix stabilization by side-chain cross-links in a small protein

Peptide cross-linking has been widely explored as a means of constraining short sequences into stable folded conformations, most commonly α-helices. The prevailing hypothesis for the origin of helix stabilization is an entropic effect resulting from backbone pre-organization; however, obtaining direct evidence bearing on this hypothesis is challenging. Here, we compare the folding thermodynamics of a small helix-rich protein domain and analogues containing one of three common cross-linking motifs. Analysis of the folding free energy landscapes of linear vs. cyclized species reveal consistent trends in the effect of cyclization on folding energetics, as well as subtle differences based on the chemistry of the cross link. Stabilization in all three systems arises entirely from a reduction in the entropic penalty of folding that more than compensates for an enthalpic destabilization of the folded state.

Stabilizing Leaf and Branch Compost Cutinase (LCC) with Glycosylation: Mechanism and Effect on PET Hydrolysis

Cutinases are polyester hydrolases that show a remarkable capability to hydrolyze polyethylene terephthalate (PET) to its monomeric units. This revelation has stimulated research aimed at developing sustainable and green cutinase-catalyzed PET recycling methods. Leaf and branch compost cutinase (LCC) is particularly suited toward these ends given its relatively high PET hydrolysis activity and thermostability. Any practical enzymatic PET recycling application will require that the protein have kinetic stability at or above the PET glass transition temperature (T_g, i.e., 70 °C). This paper elucidates the thermodynamics and kinetics of LCC conformational and colloidal stability. Aggregation emerged as a major contributor that reduces LCC kinetic stability. In its native state, LCC is prone to aggregation owing to electrostatic interactions. Further, with increasing temperature, perturbation of LCC's tertiary structure and corresponding exposure of hydrophobic domains leads to rapid aggregation. Glycosylation was employed in an attempt to impede LCC aggregation. Owing to the presence of three putative N-glycosylation sites, expression of native LCC in Pichia pastoris resulted in the production of glycosylated LCC (LCC-G). LCC-G showed improved stability to native state aggregation while increasing the temperature for thermal induced aggregation by 10 °C. Furthermore, stabilization against thermal aggregation resulted in improved catalytic PET hydrolysis both at its optimum temperature and concentration.

Unravelling the Carbohydrate-Binding Preferences of the Carbohydrate-Binding Modules of AMP-Activated Protein Kinase

The β subunit of adenosine monophosphate (AMP)-activated protein kinase (AMPK), which exists as two isoforms (β1 and β2) in humans, has a carbohydrate-binding module (CBM) that interacts with glycogen. Although the β1- and β2-CBMs are structurally similar, with strictly conserved ligand-contact residues, they show different carbohydrate affinities. β2-CBM shows the strongest affinity for both branched and unbranched oligosaccharides and it has recently been shown that a Thr insertion into β2-CBM (Thr101) forms a pocket to accommodate branches. This insertion does not explain why β2-CBM binds all carbohydrates with stronger affinity. Herein, it is shown that residue 134 (Val for β2 and Thr for β1), which does not come into contact with a carbohydrate, appears to account for the affinity difference. Characterisation by NMR spectroscopy, however, suggests that mutant β2-Thr101Δ/Val134Thr differs from that of β1-CBM, and mutant β1-Thr101ins/Thr134Val differs from that of β2-CBM. Furthermore, these mutants are less stable to chemical denaturation, relative to that of wild-type β-CBMs, which confounds the affinity analyses. To support the importance of Thr101 and Val134, the ancestral CBM has been constructed. This CBM retains Thr101 and Val134, which suggests that the extant β1-CBM has a modest loss of function in carbohydrate binding. Because the ancestor bound carbohydrate with equal affinity to that of β2-CBM, it is concluded that residue 134 plays an indirect role in carbohydrate binding.

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.

Precise estimation of transfer free energies for ionic species between similar media

Two-dimensional umbrella sampling is combined with molecular dynamics to calculate correction-free estimates of transfer properties for individual ions.

Conformational sampling of the intrinsically disordered dsRBD-1 domain from Arabidopsis thaliana DCL1

Partial folding and stability of DCL1-dsRBD1.

Structure and stability of the Human respiratory syncytial virus M2-1 RNA-binding core domain reveals a compact and cooperative folding unit

Human syncytial respiratory virus is a nonsegmented negative-strand RNA virus with serious implications for respiratory disease in infants, and has recently been reclassified into a new family, Pneumoviridae. One of the main reasons for this classification is the unique presence of a transcriptional antiterminator, called M2-1. The puzzling mechanism of action of M2-1, which is a rarity among antiterminators in viruses and is part of the RNA polymerase complex, relies on dissecting the structure and function of this multidomain tetramer. The RNA-binding activity is located in a monomeric globular `core' domain, a high-resolution crystal structure of which is now presented. The structure reveals a compact domain which is superimposable on the full-length M2-1 tetramer, with additional electron density for the C-terminal tail that was not observed in the previous models. Moreover, its folding stability was determined through chemical denaturation, which shows that the secondary and tertiary structure unfold concomitantly, which is indicative of a two-state equilibrium. These results constitute a further step in the understanding of this unique RNA-binding domain, for which there is no sequence or structural counterpart outside this virus family, in addition to its implications in transcription regulation and its likeliness as an antiviral target.

The histidine phosphocarrier protein, HPr, binds to the highly thermostable regulator of sigma D protein, Rsd, and its isolated helical fragments

The phosphotransferase system (PTS) controls the preferential use of sugars in bacteria and it is also involved in other processes, such as chemotaxis. It is formed by a protein cascade in which the first two proteins are general (namely, EI and HPr) and the others are sugar-specific permeases. The Rsd protein binds specifically to the RNA polymerase (RNAP) σ⁷⁰ factor. We first characterized the conformational stability of Escherichia coli Rsd. And second, we delineated the binding regions of Streptomyces coelicolor, HPr^sc, and E. coli Rsd, by using fragments derived from each protein. To that end, we used several biophysical probes, namely, fluorescence, CD, NMR, ITC and BLI. Rsd had a free energy of unfolding of 15 kcal mol^-1 at 25 °C, and a thermal denaturation midpoint of 103 °C at pH 6.5. The affinity between Rsd and HPr^sc was 2 μM. Interestingly enough, the isolated helical-peptides, comprising the third (RsdH3) and fourth (RsdH4) Rsd helices, also interacted with HPr^sc in a specific manner, and with affinities similar to that of the whole Rsd. Moreover, the isolated peptide of HPr^sc, HPr^9-30, comprising the active site, His15, also was bound to intact Rsd with similar affinity. Therefore, binding between Rsd and HPr^sc was modulated by the two helices H3 and H4 of Rsd, and the regions around the active site of HPr^sc. This implies that specific fragments of Rsd and HPr^sc can be used to interfere with other protein-protein interactions (PPIs) of each other protein.

Effects of Y361-auto-phosphorylation on structural plasticity of the HIPK2 kinase domain

The dual-specificity activity of the homeodomain interacting protein kinase 2 (HIPK2) is regulated by cis-auto-phosphorylation of tyrosine 361 (Y361) on the activation loop. Inhibition of this process or substitution of Y361 with nonphosphorylatable amino acid residues result in aberrant HIPK2 forms that show altered functionalities, pathological-like cellular relocalization, and accumulation into cytoplasmic aggresomes. Here, we report an in vitro characterization of wild type HIPK2 kinase domain and of two mutants, one at the regulating Y361 (Y361F, mimicking a form of HIPK2 lacking Y361 phosphorylation) and another at the catalytic lysine 228 (K228A, inactivating the enzyme). Gel filtration and thermal denaturation analyzes along with equilibrium binding experiments and kinase assays performed in the presence or absence of ATP-competitors were performed. The effects induced by mutations on overall stability, oligomerization and activity support the existence of different conformations of the kinase domain linked to Y361 phosphorylation. In addition, our in vitro data are consistent with both the cross-talk between the catalytic site and the activation loop of HIPK2 and the aberrant activities and accumulation previously reported for the Y361 nonphosphorylated HIPK2 in mammalian cells.

Prediction of New Stabilizing Mutations Based on Mechanistic Insights from Markov State Models

Protein stabilization is fundamental to enzyme function and evolution, yet understanding the determinants of a protein's stability remains a challenge. This is largely due to a shortage of atomically detailed models for the ensemble of relevant protein conformations and their relative populations. For example, the M182T substitution in TEM β-lactamase, an enzyme that confers antibiotic resistance to bacteria, is stabilizing but the precise mechanism remains unclear. Here, we employ Markov state models (MSMs) to uncover how M182T shifts the distribution of different structures that TEM adopts. We find that M182T stabilizes a helix that is a key component of a domain interface. We then predict the effects of other mutations, including a novel stabilizing mutation, and experimentally test our predictions using a combination of stability measurements, crystallography, NMR, and in vivo measurements of bacterial fitness. We expect our insights and methodology to provide a valuable foundation for protein design.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign ¹H and ¹⁵N backbone amide and ¹³C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

A Cooperative Folding Unit as the Structural Link for Energetic Coupling within a Protein

Previously, we demonstrated that binding of a ligand to Escherichia coli cofactor-dependent phosphoglycerate mutase (dPGM), a homodimeric protein, is energetically coupled with dimerization. The equilibrium unfolding of dPGM occurs with a stable, monomeric intermediate. Binding of several nonsubstrate metabolites stabilizes the dimeric native form over the monomeric intermediate, reducing the population of the intermediate. Both the active site and the dimer interface appear to be unfolded in the intermediate. We hypothesized that a loop containing residues 118-152 was responsible for the energetic coupling between the dimer interface and the distal active site and was unfolded in the intermediate. Here, we investigated the structure of the dPGM intermediate by probing side-chain interactions and solvent accessibility of the peptide backbone. By comparing the effect of a mutation on the global stability and the stability of the intermediate, we determine an equilibrium φ value (φ_eq value), which provides information about whether side-chain interactions are retained or lost in the intermediate. Hydrogen/deuterium exchange coupled with mass spectrometry (HDX-MS) was used to investigate differences in the solvent accessibility of the peptide backbone in the intermediate and native forms of dPGM. The results of φ_eq value analysis and HDX-MS reveal the least stable folding unit of dPGM, which is unfolded in the intermediate and links the active site to the dimer interface. The structure of the intermediate reveals how the cooperative network of residues in dPGM gives rise to the observed energetic coupling between dimerization and ligand binding.

Cooperative protein unfolding. A statistical-mechanical model for the action of denaturants

Knowledge of protein stability is of utmost importance in various fields of biotechnology. Protein stability can be assessed in solution by increasing the concentration of denaturant and recording the structural changes with spectroscopic or thermodynamic methods. The standard interpretation of the experimental data is to assume a 2-state equilibrium between completely folded and completely unfolded protein molecules. Here we propose a cooperative model based on the statistical-mechanical Zimm-Bragg theory. In this model protein unfolding is driven by the weak binding of a rather small number of denaturant molecules, inducing the cooperative unfolding with multiple dynamic intermediates. The modified Zimm-Bragg theory is applied to published thermodynamic and spectroscopic data leading to the following conclusions. (i) The binding constant K_D is correlated with the midpoint concentration, c₀, of the unfolding reaction according to c₀≅1/K_D. The better the binding of denaturant the lower is the concentration to achieve unfolding. (ii) The binding constant K_D agrees with direct thermodynamic measurements. A rather small number of bound denaturants suffices to induce the cooperative unfolding of the whole protein. (iii) Chemical unfolding occurs in the concentration range Δc_D=c_end-c_ini. The theory predicts the unfolding energy per amino acid residue as g_nu=RTK_D(c_end-c_ini). The Gibbs free energy of an osmotic gradient of the same size is ΔG_Diff=-RTln(c_end/c_ini). In all examples investigated ΔG_Diff exactly balances the unfolding energy g_nu. The total unfolding energy is thus close to zero. (iv) Protein cooperativity in chemical unfolding is rather low with cooperativity parameters σ≥3x10^-3. As a consequence, the theory predicts a dynamic mixture of conformations during the unfolding reaction. The probabilities of individual conformations are easily accessible via the partition function Z(c_D,σ).

Conservation of folding and association within a family of spidroin N-terminal domains

Web spiders synthesize silk fibres, nature’s toughest biomaterial, through the controlled assembly of fibroin proteins, so-called spidroins. The highly conserved spidroin N-terminal domain (NTD) is a pH-driven self-assembly device that connects spidroins to super-molecules in fibres. The degree to which forces of self-assembly is conserved across spider glands and species is currently unknown because quantitative measures are missing. Here, we report the comparative investigation of spidroin NTDs originating from the major ampullate glands of the spider species Euprosthenops australis, Nephila clavipes, Latrodectus hesperus, and Latrodectus geometricus. We characterized equilibrium thermodynamics and kinetics of folding and self-association using dynamic light scattering, stopped-flow fluorescence and circular dichroism spectroscopy in combination with thermal and chemical denaturation experiments. We found cooperative two-state folding on a sub-millisecond time scale through a late transition state of all four domains. Stability was compromised by repulsive electrostatic forces originating from clustering of point charges on the NTD surface required for function. pH-driven dimerization proceeded with characteristic fast kinetics yielding high affinities. Results showed that energetics and kinetics of NTD self-assembly are highly conserved across spider species despite the different silk mechanical properties and web geometries they produce.

Long-term protein packaging in bio-ionic liquids: Improved catalytic activity and enhanced stability of cytochrome C against multiple stresses

There is a considerable interest in the use of structurally stable and catalytically active enzymes, such as cytochrome C (Cyt C), in the pharmaceutical and fine chemical industries. However, harsh process conditions, such as temperature, pH, and presence of organic solvents, are the major barriers to the effective use of enzymes in biocatalysis. Herein, we demonstrate the suitability of bio-based ionic liquids (ILs) formed by the cholinium cation and dicarboxylate-based anions as potential media for enzymes, in which remarkable enhanced activity and improved stability of Cyt C against multiple stresses were obtained. Among the several bio-ILs studied, an exceptionally high catalytic activity (> 50-fold) of Cyt C was observed in aqueous solutions of cholinium glutarate ([Ch][Glu]; 1g/mL) as compared to the commonly used phosphate buffer solutions (pH 7.2), and > 25-fold as compared to aqueous solutions of cholinium dihydrogen phosphate ([Ch][Dhp]; 0.5g/mL) -the best known IL for long term stability of Cyt C. The catalytic activity of the enzyme in presence of bio-ILs was retained against several external stimulus, such as chemical denaturants (H2O2 and GuHCl), and temperatures up to 120 °C. The observed enzyme activity is in agreement with its structural stability, as confirmed by UV-Vis, circular dichroism (CD), and Fourier transform infrared (FT-IR) spectroscopies. Taking advantage of the multi-ionization states of di/tri-carboxylic acids, the pH was switched from acidic to basic by the addition of the corresponding carboxylic acid and choline hydroxide, respectively. The activity was found to be maximum at a 1:1 ratio of [Ch][carboxylate], with a pH in the range from 3 to 5.5. Moreover, it was found that the bio-ILs studied herein protect the enzyme against protease digestion and allow long-term storage (at least for 21 weeks) at room temperature. An attempt by molecular docking was also made to better understand the efficacy of the investigated bio-ILs towards the enhanced activity and long term stability of Cyt C. The results showed that dicarboxylates anions interact with the active site's amino acids of the enzyme through H-bonding and electrostatic interactions, which are responsible for the observed enhancement of the catalytic activity. Finally, it is demonstrated that Cyt C can be successfully recovered from the aqueous solution of bio-ILs and reused without compromising its yield, structural integrity and catalytic activity, thereby overcoming the major limitations in the use of IL-protein systems in biocatalysis.

Foldability of a Natural De Novo Evolved Protein

The de novo evolution of protein-coding genes from noncoding DNA is emerging as a source of molecular innovation in biology. Studies of random sequence libraries, however, suggest that young de novo proteins will not fold into compact, specific structures typical of native globular proteins. Here we show that Bsc4, a functional, natural de novo protein encoded by a gene that evolved recently from noncoding DNA in the yeast S. cerevisiae, folds to a partially specific three-dimensional structure. Bsc4 forms soluble, compact oligomers with high β sheet content and a hydrophobic core, and undergoes cooperative, reversible denaturation. Bsc4 lacks a specific quaternary state, however, existing instead as a continuous distribution of oligomer sizes, and binds dyes indicative of amyloid oligomers or molten globules. The combination of native-like and non-native-like properties suggests a rudimentary fold that could potentially act as a functional intermediate in the emergence of new folded proteins de novo.

Stability and aggregation propensity do not fully account for the association of various germline variable domain gene segments with light chain amyloidosis

Variable domain (VL) gene segments exhibit variable tendencies to be associated with light chain amyloidosis (AL). While few of them are very frequent in AL and give rise to most of the amyloidogenic light chains compiled at the sequence databases, other are rarely found among the AL cases. To analyze to which extent these tendencies depend on folding stability and aggregation propensity of the germline VL protein, we characterized VL proteins encoded by four AL-associated germline gene segments and one not associated to AL. We found that the AL-associated germline rVL proteins differ widely in conformational stability and propensity to in vitro amyloid aggregation. While in vitro the amyloid formation kinetics of these proteins correlate well with their folding stabilities, the folding stability does not clearly correlate with their germline's frequencies in AL. We conclude that the association of the VL genes segments to amyloidosis is not determined solely by the folding stability and aggregation propensity of the germline VL protein. Other factors, such as the frequencies of destabilizing mutations and susceptibility to proteolysis, must play a role in determining the light chain amyloidogenicity.

Addressing the role of the α-helical extension in the folding of the third PDZ domain from PSD-95

PDZ domains are one of the most important protein-protein interaction domains in human. While presenting a conserved three dimensional structure, a substantial number of PDZ domains display structural extensions suggested to be involved in their folding and binding mechanisms. The C-terminal α-helix extension (α3) of the third PDZ domain from PSD-95 (PDZ3) has been reported to have a role in function of the domain as well as in the stabilization of the native fold. Here we report an evaluation of the effect of the truncation of this additional helix on the folding and unfolding kinetics of PDZ3. Fluorescent variants of full length and truncated PDZ3 were produced and stopped-flow fluorescence measurements were made under different experimental conditions (pH, ionic strength and temperature) to investigate the folding kinetics of the respective variant. The results show that folding of PDZ3 is robust and that the mechanism is only marginally affected by the truncation, which contributes to a destabilization of the native state, but otherwise do not change the overall observed kinetics. Furthermore, the increase in the unfolding rate constants, but not the folding rate constant upon deletion of α3 suggests that the α-helical extension is largely unstructured in the folding transition state.

Comparative analyses of the thermodynamic RNA binding signatures of different types of RNA recognition motifs

RNA recognition motifs (RRMs) are structurally versatile domains important in regulation of alternative splicing. Structural mechanisms of sequence-specific recognition of single-stranded RNAs (ssRNAs) by RRMs are well understood. The thermodynamic strategies are however unclear. Therefore, we utilized microcalorimetry and semi-empirical analyses to comparatively analyze the cognate ssRNA binding thermodynamics of four different RRM domains, each with a different RNA binding mode. The different binding modes are: canonical binding to the β-sheet surface; canonical binding with involvement of N- and C-termini; binding to conserved loops; and binding to an α-helix. Our results identify enthalpy as the sole and general force driving association at physiological temperatures. Also, networks of weak interactions are a general feature regulating stability of the different RRM–ssRNA complexes. In agreement, non-polyelectrolyte effects contributed between ∼75 and 90% of the overall free energy of binding in the considered complexes. The various RNA binding modes also displayed enormous heat capacity differences, that upon dissection revealed large differential changes in hydration, conformations and dynamics upon binding RNA. Altogether, different modes employed by RRMs to bind cognate ssRNAs utilize various thermodynamics strategies during the association process.

Kinetic stability of membrane proteins

Although membrane proteins constitute an important class of biomolecules involved in key cellular processes, study of the thermodynamic and kinetic stability of their structures is far behind that of soluble proteins. It is known that many membrane proteins become unstable when removed by detergent extraction from the lipid environment. In addition, most of them undergo irreversible denaturation, even under mild experimental conditions. This process was found to be associated with partial unfolding of the polypeptide chain exposing hydrophobic regions to water, and it was proposed that the formation of kinetically trapped conformations could be involved. In this review, we will describe some of the efforts toward understanding the irreversible inactivation of membrane proteins. Furthermore, its modulation by phospholipids, ligands, and temperature will be herein discussed.

Apparent activation energies of protein-protein complex dissociation in the gas-phase determined by electrospray mass spectrometry

We have developed a method to determine apparent activation energies of dissociation for ionized protein-protein complexes in the gas phase using electrospray ionization mass spectrometry following the Rice-Ramsperger-Kassel-Marcus quasi-equilibrium theory. Protein-protein complexes were formed in solution, transferred into the gas phase, and separated from excess free protein by ion mobility filtering. Afterwards, complex disassembly was initiated by collision-induced dissociation with step-wise increasing energies. Relative intensities of ion signals were used to calculate apparent activation energies of dissociation in the gas phase by applying linear free energy relations. The method was developed using streptavidin tetramers. Experimentally determined apparent gas-phase activation energies for dissociation ([Formula: see text]) of complexes consisting of Fc parts from immunoglobulins (IgG-Fc) and three closely related protein G' variants (IgG-Fc•protein G'e, IgG-Fc•protein G'f, and IgG-Fc•protein G'g) show the same order of stabilities as can be inferred from their in-solution binding constants. Differences in stabilities between the protein-protein complexes correspond to single amino acid residue exchanges in the IgG-binding regions of the protein G' variants. Graphical abstract Electrospray mass spectrometry and collision-induced dissociation delivers apparent activation energies and supramolecular bond force constants of protein-protein complexes in the gas phase.

Biophysical analysis of the effect of chemical modification by 4-oxononenal on the structure, stability, and function of binding immunoglobulin protein (BiP)

Binding immunoglobulin protein (BiP) is a molecular chaperone important for the folding of numerous proteins, which include millions of immunoglobulins in human body. It also plays a key role in the unfolded protein response (UPR) in the endoplasmic reticulum. Free radical generation is a common phenomenon that occurs in cells under healthy as well as under stress conditions such as ageing, inflammation, alcohol consumption, and smoking. These free radicals attack the cell membranes and generate highly reactive lipid peroxidation products such as 4-oxononenal (4-ONE). BiP is a key protein that is modified by 4-ONE. In this study, we probed how such chemical modification affects the biophysical properties of BiP. Upon modification, BiP shows significant tertiary structural changes with no changes in its secondary structure. The protein loses its thermodynamic stability, particularly, that of the nucleotide binding domain (NBD) where ATP binds. In terms of function, the modified BiP completely loses its ATPase activity with decreased ATP binding affinity. However, modified BiP retains its immunoglobulin binding function and its chaperone activity of suppressing non-specific protein aggregation. These results indicate that 4-ONE modification can significantly affect the structure-function of key proteins such as BiP involved in cellular pathways, and provide a molecular basis for how chemical modifications can result in the failure of quality control mechanisms inside the cell.

The effect of microwave on the interaction of flavour compounds with G-actin from grass carp (Catenopharyngodon idella)

BACKGROUND

In order to investigate the influence of non-thermal effects of microwaves on the flavour of fish and meat products, the G-actin of grass carp in ice baths was exposed to different microwave powers (0, 100, 300 or 500 W); the surface hydrophobicity, sulfhydryl contents, secondary structures and adsorption capacity of G-actin to ketones were determined.

RESULTS

As microwave power increased from 0 to 300 W, the surface hydrophobicity, total and reactive sulfhydryls increased; α-helix, β-sheet and random coil fractions turned into β-turn fractions. As microwave power increased from 300 to 500 W, however, hydrophobicity and sulfhydryl contents decreased; β-turn and random coil fractions turned into α-helix and β-sheet fractions. The tendencies of adsorbed capacity of ketones were similar to hydrophobicity and sulfhydryl contents.

CONCLUSION

The increased adsorbing of ketones could be attributed to the unfolding of secondary structures by revealing new binding sites, including thiol groups and hydrophobic binding sites. The decreased binding capacity was related to the refolding and aggregation of protein. The results suggested that microwave powers had obvious effects on the flavour retention and proteins structures in muscle foods. © 2017 Society of Chemical Industry.

Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force, Temperature, and Urea

While it is widely appreciated that the denatured state of a protein is a heterogeneous conformational ensemble, there is still debate over how this ensemble changes with environmental conditions. Here, we use single-molecule chemo-mechanical unfolding, which combines force and urea using the optical tweezers, together with traditional protein unfolding studies to explore how perturbants commonly used to unfold proteins (urea, force, and temperature) affect the denatured-state ensemble. We compare the urea m-values, which report on the change in solvent accessible surface area for unfolding, to probe the denatured state as a function of force, temperature, and urea. We find that while the urea- and force-induced denatured states expose similar amounts of surface area, the denatured state at high temperature and low urea concentration is more compact. To disentangle these two effects, we use destabilizing mutations that shift the T_m and C_m. We find that the compaction of the denatured state is related to changing temperature as the different variants of acyl-coenzyme A binding protein have similar m-values when they are at the same temperature but different urea concentration. These results have important implications for protein folding and stability under different environmental conditions.

l-Proline and RNA Duplex m-Value Temperature Dependence

The temperature dependence of l-proline interactions with the RNA dodecamer duplex surface exposed after unfolding was quantified using thermal and isothermal titration denaturation monitored by uv-absorbance. The m-value quantifying proline interactions with the RNA duplex surface area exposed after unfolding was measured using RNA duplexes with GC content ranging between 17 and 83%. The m-values from thermal denaturation decreased with increasing GC content signifying increasingly favorable proline interactions with the exposed RNA surface area. However, m-values from isothermal titration denaturation at 25.0 °C were independent of GC content and less negative than those from thermal denaturation. The m-value from isothermal titration denaturation for a 50% GC RNA duplex decreased (became more negative) as the temperature increased and was in nearly exact agreement with the m-value from thermal denaturation. Since RNA duplex transition temperatures increased with GC content, the more favorable proline interactions with the high GC content duplex surface area observed from thermal denaturation resulted from the temperature dependence of proline interactions rather than the RNA surface chemical composition. The enthalpy contribution to the m-value was positive and small (indicating a slight increase in duplex unfolding enthalpy with proline) while the entropic contribution to the m-value was positive and increased with temperature. Our results will facilitate proline's use as a probe of solvent accessible surface area changes during biochemical reactions at different reaction temperatures.

Heat-induced native dimerization prevents amyloid formation by variable domain from immunoglobulin light-chain REI

Amyloid light-chain (AL) amyloidosis is a protein-misfolding disease characterized by accumulation of immunoglobulin light chains (LCs) into amyloid fibrils. Dimerization of a full length or variable domain (V_L ) of LC serves to stabilize the native state and prevent the formation of amyloid fibrils. We here analyzed the thermodynamic properties of dimerization and unfolding reactions by nonamyloidogenic V_L from REI LC or its monomeric Y96K mutant using sedimentation velocity and circular dichroism. The data indicate that the equilibrium shifts to native dimerization for wild-type REI V_L by elevating temperature due to the negative enthalpy change for dimer dissociation (-81.2 kJ·mol^-1 ). The Y96K mutation did not affect the stability of the monomeric native state but increased amyloidogenicity. These results suggest that the heat-induced native homodimerization is the major factor preventing amyloid formation by wild-type REI V_L . Heat-induced native oligomerization may be an efficient strategy to avoid the formation of misfolded aggregates particularly for thermostable proteins that are used at elevated temperatures under conditions where other proteins tend to misfold.

DATABASE

Structural data are available in the Protein Data Bank under the accession numbers 5XP1 and 5XQY.

Three-Dimensional Domain Swapping Changes the Folding Mechanism of the Forkhead Domain of FoxP1

The forkhead family of transcription factors (Fox) controls gene transcription during key processes such as regulation of metabolism, embryogenesis, and immunity. Structurally, Fox proteins feature a conserved DNA-binding domain known as forkhead. Interestingly, solved forkhead structures of members from the P subfamily (FoxP) show that they can oligomerize by three-dimensional domain swapping, whereby structural elements are exchanged between adjacent subunits, leading to an intertwined dimer. Recent evidence has largely stressed the biological relevance of domain swapping in FoxP, as several disease-causing mutations have been related to impairment of this process. Here, we explore the equilibrium folding and binding mechanism of the forkhead domain of wild-type FoxP1, and of two mutants that hinder DNA-binding (R53H) and domain swapping (A39P), using size-exclusion chromatography, circular dichroism, and hydrogen-deuterium exchange mass spectrometry. Our results show that domain swapping of FoxP1 occurs at micromolar protein concentrations within hours of incubation and is energetically favored, in contrast to classical domain-swapping proteins. Also, DNA-binding mutations do not significantly affect domain swapping. Remarkably, equilibrium unfolding of dimeric FoxP1 follows a three-state N2 ↔ 2I ↔ 2U folding mechanism in which dimer dissociation into a monomeric intermediate precedes protein unfolding, in contrast to the typical two-state model described for most domain-swapping proteins, whereas the A39P mutant follows a two-state N ↔ U folding mechanism consistent with the second transition observed for dimeric FoxP1. Also, the free-energy change of the N ↔ U in A39P FoxP1 is ∼2 kcal⋅mol(-1) larger than the I ↔ U transition of both wild-type and R53H FoxP1. Finally, hydrogen-deuterium exchange mass spectrometry reveals that the intermediate strongly resembles the native state. Our results suggest that domain swapping in FoxP1 is at least partially linked to monomer folding stability and follows an unusual three-state folding mechanism, which might proceed via transient structural changes rather than requiring complete protein unfolding as do most domain-swapping proteins.

Pseudomonas aeruginosa cytochrome c551 denaturation by five systematic urea derivatives that differ in the alkyl chain length

Reversible denaturation of Pseudomonas aeruginosa cytochrome c551 (PAc551) could be followed using five systematic urea derivatives that differ in the alkyl chain length, i.e. urea, N-methylurea (MU), N-ethylurea (EU), N-propylurea (PU), and N-butylurea (BU). The BU concentration was the lowest required for the PAc551 denaturation, those of PU, EU, MU, and urea being gradually higher. Furthermore, the accessible surface area difference upon PAc551 denaturation caused by BU was found to be the highest, those by PU, EU, MU, and urea being gradually lower. These findings indicate that urea derivatives with longer alkyl chains are stronger denaturants. In this study, as many as five systematic urea derivatives could be applied for the reversible denaturation of a single protein, PAc551, for the first time, and the effects of the alkyl chain length on protein denaturation were systematically verified by means of thermodynamic parameters.

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.

Urea-Dependent Adenylate Kinase Activation following Redistribution of Structural States

Proteins are often functionally dependent on conformational changes that allow them to sample structural states that are sparsely populated in the absence of a substrate or binding partner. The distribution of such structural microstates is governed by their relative stability, and the kinetics of their interconversion is governed by the magnitude of associated activation barriers. Here, we have explored the interplay among structure, stability, and function of a selected enzyme, adenylate kinase (Adk), by monitoring changes in its enzymatic activity in response to additions of urea. For this purpose we used a ³¹P NMR assay that was found useful for heterogeneous sample compositions such as presence of urea. It was found that Adk is activated at low urea concentrations whereas higher urea concentrations unfolds and thereby deactivates the enzyme. From a quantitative analysis of chemical shifts, it was found that urea redistributes preexisting structural microstates, stabilizing a substrate-bound open state at the expense of a substrate-bound closed state. Adk is rate-limited by slow opening of substrate binding domains and the urea-dependent redistribution of structural states is consistent with a model where the increased activity results from an increased rate-constant for domain opening. In addition, we also detected a strong correlation between the catalytic free energy and free energy of substrate (ATP) binding, which is also consistent with the catalytic model for Adk. From a general perspective, it appears that urea can be used to modulate conformational equilibria of folded proteins toward more expanded states for cases where a sizeable difference in solvent-accessible surface area exists between the states involved. This effect complements the action of osmolytes, such as trimethylamine N-oxide, that favor more compact protein states.

The Folding Pathway of the KIX Domain

The KIX domain is an 89-residues globular domain with an important role in mediating protein-protein interactions. The presence of two distinct binding sites in such a small domain makes KIX a suitable candidate to investigate the effect of the potentially divergent demands between folding and function. Here, we report an extensive mutational analysis of the folding pathway of the KIX domain, based on 30 site-directed mutants, which allow us to assess the structures of both the transition and denatured states. Data reveal that, while the transition state presents mostly native-like interactions, the denatured state is somewhat misfolded. We mapped some of the non-native contacts in the denatured state using a second round of mutagenesis, based on double mutant cycles on 15 double mutants. Interestingly, such a misfolding arises from non-native interactions involving the residues critical for the function of the protein. The results described in this work appear to highlight the diverging demands between folding and function that may lead to misfolding, which may be observed in the early stages of folding.

Energetics of side-chain partitioning of β-signal residues in unassisted folding of a transmembrane β-barrel protein

The free energy of water-to-interface amino acid partitioning is a major contributing factor in membrane protein folding and stability. The interface residues at the C terminus of transmembrane β-barrels form the β-signal motif required for assisted β-barrel assembly in vivo but are believed to be less important for β-barrel assembly in vitro. Here, we experimentally measured the thermodynamic contribution of all 20 amino acids at the β-signal motif to the unassisted folding of the model β-barrel protein PagP. We obtained the partitioning free energy for all 20 amino acids at the lipid-facing interface (ΔΔG⁰_w,i(φ)) and the protein-facing interface (ΔΔG⁰_w,i(π)) residues and found that hydrophobic amino acids are most favorably transferred to the lipid-facing interface, whereas charged and polar groups display the highest partitioning energy. Furthermore, the change in non-polar surface area correlated directly with the partitioning free energy for the lipid-facing residue and inversely with the protein-facing residue. We also demonstrate that the interface residues of the β-signal motif are vital for in vitro barrel assembly, because they exhibit a side chain–specific energetic contribution determined by the change in nonpolar accessible surface. We further establish that folding cooperativity and hydrophobic collapse are balanced at the membrane interface for optimal stability of the PagP β-barrel scaffold. We conclude that the PagP C-terminal β-signal motif influences the folding cooperativity and stability of the folded β-barrel and that the thermodynamic contributions of the lipid- and protein-facing residues in the transmembrane protein β-signal motif depend on the nature of the amino acid side chain.

Fluorescence correlation spectroscopy reveals a cooperative unfolding of monomeric amyloid-β 42 with a low Gibbs free energy

The amyloid-beta peptide (Aβ) plays a major role in the progression of Alzheimer's disease. Due to its high toxicity, the 42 amino acid long isoform Aβ42 has become of considerable interest. The Aβ42 monomer is prone to aggregation down to the nanomolar range which makes conventional structural methods such as NMR or X-ray crystallography infeasible. Conformational information, however, will be helpful to understand the different aggregation pathways reported in the literature and will allow to identify potential conditions that favour aggregation-incompetent conformations. In this study, we applied fluorescence correlation spectroscopy (FCS) to investigate the unfolding of Alexa Fluor 488 labelled monomeric Aβ42 using guanidine hydrochloride as a denaturant. We show that our Aβ42 pre-treatment and the low-nanomolar concentrations, typically used for FCS measurements, strongly favour the presence of monomers. Our results reveal that there is an unfolding/folding behaviour of monomeric Aβ42. The existence of a cooperative unfolding curve suggests the presence of structural elements with a Gibbs free energy of unfolding of about 2.8 kcal/mol.

Simulative and Experimental Characterization of a pH-Dependent Clamp-like DNA Triple-Helix Nanoswitch

Here we couple experimental and simulative techniques to characterize the structural/dynamical behavior of a pH-triggered switching mechanism based on the formation of a parallel DNA triple helix. Fluorescent data demonstrate the ability of this structure to reversibly switch between two states upon pH changes. Two accelerated, half microsecond, MD simulations of the system having protonated or unprotonated cytosines, mimicking the pH 5.0 and 8.0 conditions, highlight the importance of the Hoogsteen interactions in stabilizing the system, finely depicting the time-dependent disruption of the hydrogen bond network. Urea-unfolding experiments and MM/GBSA calculations converge in indicating a stabilization energy at pH 5.0, 2-fold higher than that observed at pH 8.0. These results validate the pH-controlled behavior of the designed structure and suggest that simulative approaches can be successfully coupled with experimental data to characterize responsive DNA-based nanodevices.

Tuning the Continuum of Structural States in the Native Ensemble of a Regulatory Protein

The mesoscale nature of proteins allows for an efficient coupling between environmental cues and conformational changes, enabling their function as molecular transducers. Delineating the precise structural origins of such a connection and the expected spectroscopic response has, however, been challenging. In this work, we perform a combination of urea-temperature double perturbation experiments and theoretical modeling to probe the conformational landscape of Cnu, a natural thermosensor protein. We observe unique ensemble signatures that point to a continuum of conformational substates in the native ensemble and that respond intricately to perturbations upon monitoring secondary and tertiary structures, distances between an intrinsic FRET pair, and hydrodynamic volumes. Binding assays further reveal a weakening of the Cnu functional complex with temperature, highlighting the molecular origins of signal transduction critical for pathogenic response in enterobacteriaceae.

Creating a Homeodomain with High Stability and DNA Binding Affinity by Sequence Averaging

There is considerable interest in generating proteins with both high stability and high activity for biomedical and industrial purposes. One approach that has been used successfully to increase the stability of linear repeat proteins is consensus design. It is unclear the extent over which the consensus design approach can be used to produce folded and hyperstable proteins, and importantly, whether such stabilized proteins would retain function. Here we extend the consensus strategy to design a globular protein. We show that a consensus-designed homeodomain (HD) sequence adopts a cooperatively folded homeodomain structure. The unfolding free energy of the consensus-HD is 5 kcal·mol^-1 higher than that of the naturally occurring engrailed-HD from Drosophila melanogaster. Remarkably, the consensus-HD binds the engrailed-HD cognate DNA in a similar mode as the engrailed-HD with approximately 100-fold higher affinity. ¹⁵N relaxation studies show a decrease in ps-ns backbone dynamics in the free state of consensus-HD, suggesting that increased affinity is not a result of increased plasticity. In addition to demonstrating the potential for consensus design of globular proteins with increased stability, these results demonstrate that greatly stabilized proteins can bind cognate substrates with increased affinities, showing that high stability is compatible with function.