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

Harnessing Environmental Ca2+ for Extracellular Protein Thermostabilization

Ca²⁺ is the third-most prevalent metal ion in the environment. EF hands are common Ca²⁺-binding motifs found in both extracellular and intracellular proteins of eukaryotes and prokaryotes. Cytoplasmic EF hand proteins often mediate allosteric control of signal transduction pathway components in response to intracellular Ca²⁺ concentration fluctuations by coupling Ca²⁺ binding to changes in protein structure. We show that an extracellular structural Ca²⁺-binding site mediates protein thermostabilization by such conformational coupling as well. Binding Ca²⁺ to the EF hand of the extracellular (periplasmic) Escherichia coli glucose-galactose binding protein thermostabilizes this protein by ∼17 K relative to its Ca²⁺-free form. Using statistical thermodynamic analysis of a fluorescent conjugate of ecGGBP that reports simultaneously on ligand binding and multiple conformational states, we found that its Ca²⁺-mediated stabilization is determined by conformational coupling mechanisms in two independent conformational exchange reactions. Binding to folded and unfolded states determines the maximum Ca²⁺-mediated stability. A disorder → order transition accompanies the formation of the Ca²⁺ complex in the folded state and dictates the minimum Ca²⁺ concentration at which the Ca²⁺-bound state becomes dominant. Similar transitions also encode the structural changes necessary for Ca²⁺-mediated control elements in signal transduction pathways. Ca²⁺-mediated thermostabilization and allosteric control, therefore, share a fundamental conformational coupling mechanism, which may have implications for the evolution of EF hands.

Target identification and validation of natural products with label-free methodology: A critical review from 2005 to 2020

Natural products (NPs) have been an important source of therapeutic drugs in clinic use and contributed many chemical probes for research. The usefulness of NPs is however often marred by the incomplete understanding of their direct cellular targets. A number of experimental methods for drug target identification have been developed over the years. One class of methods, termed "label-free" methodology, exploits the energetic and biophysical features accompanying the association of macromolecules with drugs and other compounds in their native forms. Herein we review the working principles, assay implementations, and key applications of the most important approaches, and also give examples where they have been applied to NPs. We also assess the key advantages and limitations of each method. Furthermore, we address when and how the label-free methodology can be particularly useful considering some of the unique features of NP chemistry and bioactivation.

A Highly Conserved Iron-Sulfur Cluster Assembly Machinery between Humans and Amoeba Dictyostelium discoideum: The Characterization of Frataxin

Several biological activities depend on iron-sulfur clusters ([Fe-S]). Even though they are well-known in several organisms their function and metabolic pathway were poorly understood in the majority of the organisms. We propose to use the amoeba Dictyostelium discoideum, as a biological model to study the biosynthesis of [Fe-S] at the molecular, cellular and organism levels. First, we have explored the D. discoideum genome looking for genes corresponding to the subunits that constitute the molecular machinery for Fe-S cluster assembly and, based on the structure of the mammalian supercomplex and amino acid conservation profiles, we inferred the full functionality of the amoeba machinery. After that, we expressed the recombinant mature form of D. discoideum frataxin protein (DdFXN), the kinetic activator of this pathway. We characterized the protein and its conformational stability. DdFXN is monomeric and compact. The analysis of the secondary structure content, calculated using the far-UV CD spectra, was compatible with the data expected for the FXN fold, and near-UV CD spectra were compatible with the data corresponding to a folded protein. In addition, Tryptophan fluorescence indicated that the emission occurs from an apolar environment. However, the conformation of DdFXN is significantly less stable than that of the human FXN, (4.0 vs. 9.0 kcal mol-1, respectively). Based on a sequence analysis and structural models of DdFXN, we investigated key residues involved in the interaction of DdFXN with the supercomplex and the effect of point mutations on the energetics of the DdFXN tertiary structure. More than 10 residues involved in Friedreich's Ataxia are conserved between the human and DdFXN forms, and a good correlation between mutational effect on the energetics of both proteins were found, suggesting the existence of similar sequence/function/stability relationships. Finally, we integrated this information in an evolutionary context which highlights particular variation patterns between amoeba and humans that may reflect a functional importance of specific protein positions. Moreover, the complete pathway obtained forms a piece of evidence in favor of the hypothesis of a shared and highly conserved [Fe-S] assembly machinery between Human and D. discoideum.

Quantitating denaturation by formic acid: imperfect repeats are essential to the stability of the functional amyloid protein FapC

Bacterial functional amyloids are evolutionarily optimized to aggregate, so much so that the extreme robustness of functional amyloid makes it very difficult to examine their structure-function relationships in a detailed manner. Previous work has shown that functional amyloids are resistant to conventional chemical denaturants, but they dissolve in formic acid (FA) at high concentrations. However, systematic investigation requires a quantitative analysis of FA's ability to denature proteins. Amyloid formed by Pseudomonas sp. protein FapC provides an excellent model to investigate FA denaturation. It contains three imperfect repeats, and stepwise removal of these repeats slows fibrillation and increases fragmentation during aggregation. However, the link to stability is unclear. We first calibrated FA denaturation using three small, globular, and acid-resistant proteins. This revealed a linear relationship between the concentration of FA and the free energy of unfolding with a slope of mFA+pH (the combined contribution of FA and FA-induced lowering of pH), as well as a robust correlation between protein size and mFA+pH We then measured the solubilization of fibrils formed from different FapC variants with varying numbers of repeats as a function of the concentration of FA. This revealed a decline in the number of residues driving amyloid formation upon deleting at least two repeats. The midpoint of denaturation declined with the removal of repeats. Complete removal of all repeats led to fibrils that were solubilized at FA concentrations 2-3 orders of magnitude lower than the repeat-containing variants, showing that at least one repeat is required for the stability of functional amyloid.

Apolipoprotein E4 exhibits intermediates with domain interaction

ApoE4(C112R) is the strongest risk factor for Alzheimer's disease, while apoE3(C112) is considered normal. The C112R substitution is believed to alter the interactions between the N-terminal (NTD) and the C-terminal domain (CTD) leading to major functional differences. Here we investigate how the molecular property of the residue at position 112 affects domain interaction using an array of C112X substitutions with arginine, alanine, threonine, valine, leucine and isoleucine as 'X'. We attempt to determine the free energy of domain interaction (∆G_INT) from stabilities of the NTD (∆G_NTD) and CTD (∆G_CTD) in the full-length apoE, and the stabilities of fragments of the NTD (∆G_NTF) and CTD (∆G_CTF), using the relationship, ∆G_INT = ∆G_NTD + ∆G_CTD - ∆G_NTF - ∆G_CTF. We find that although ∆G_NTD is strongly dependent on the C112X substitutions, ∆G_NTD - ∆G_NTF is small. Furthermore, ∆G_CTD remains nearly the same as ∆G_CTF. Therefore, ∆G_INT is estimated to be small and similar for the apoE isoforms. However, stability of domain interaction monitored by urea dependent changes in interdomain Forster Resonance Energy Transfer (FRET) is found to be strongly dependent on C112X substitutions. ApoE4 exhibits the highest mid-point of denaturation of interdomain FRET. To resolve the apparently contradictory observations, we hypothesize that higher interdomain FRET in apoE4 in urea may involve 'intermediate' states. Enhanced fluorescence of bis-ANS and susceptibility to proteolytic cleavage support that apoE4, specifically, the NTD of apoE4 harbor 'intermediates' in both native and mildly denaturing conditions. The intermediates could hold key to the pathological functions of apoE4.

Hidden kinetic traps in multidomain folding highlight the presence of a misfolded but functionally competent intermediate

Much of our current knowledge on protein folding is based on work focused on isolated domains. In this study, using a combination of NMR and kinetic experiments, we depict the folding pathway of a multidomain construct comprising two PDZ domains in tandem, belonging to the protein Whirlin. We demonstrate the presence of a misfolded intermediate that competes with productive folding. Interestingly, we show that, unexpectedly, this misfolded state retains the native-like functional ability to bind its physiological ligand, representing a clear example of a functionally competent misfolded state. On the basis of these results and a comparative analysis of the amino acidic sequences of Whirlin from different species, we propose a possible physiological role of the misfolded intermediate.

Although more than 75% of the proteome is composed of multidomain proteins, current knowledge of protein folding is based primarily on studies of isolated domains. In this work, we describe the folding mechanism of a multidomain tandem construct comprising two distinct covalently bound PDZ domains belonging to a protein called Whirlin, a scaffolding protein of the hearing apparatus. In particular, via a synergy between NMR and kinetic experiments, we demonstrate the presence of a misfolded intermediate that competes with productive folding. In agreement with the view that tandem domain swapping is a potential source of transient misfolding, we demonstrate that such a kinetic trap retains native-like functional activity, as shown by the preserved ability to bind its physiological ligand. Thus, despite the general knowledge that protein misfolding is intimately associated with dysfunction and diseases, we provide a direct example of a functionally competent misfolded state. Remarkably, a bioinformatics analysis of the amino acidic sequence of Whirlin from different species suggests that the tendency to perform tandem domain swapping between PDZ1 and PDZ2 is highly conserved, as demonstrated by their unexpectedly high sequence identity. On the basis of these observations, we discuss on a possible physiological role of such misfolded intermediate.

Global analysis of protein stability by temperature and chemical denaturation

The stability of a protein is a fundamental property that determines under which conditions, the protein is functional. Equilibrium unfolding with denaturants requires preparation of several samples and only provides the free energy of folding when performed at a single temperature. The typical sample requirement is around 0.5-1 mg of protein. If the stability of many proteins or protein variants needs to be determined, substantial protein production may be needed. Here we have determined the stability of acyl-coenzyme A binding protein at pH 5.3 and chymotrypsin inhibitor 2 at pH 3 and pH 6.25 by combined temperature and denaturant unfolding. We used a setup where tryptophan fluorescence is measured in quartz capillaries where only 10 μl is needed. Temperature unfolding of a series of 15 samples at increasing denaturant concentrations provided accurate and precise thermodynamic parameters. We find that the number of samples may be further reduced and less than 10 μg of protein in total are needed for reliable stability measurements. For assessment of stability of protein purified in small scale e.g. in micro plate format, our method will be highly applicable. The routine for fitting the experimental data is made available as a python notebook.

Following Structural Changes by Thermal Denaturation Using Trapped Ion Mobility Spectrometry-Mass Spectrometry

The behavior of biomolecules as a function of the solution temperature is often crucial to assessing their biological activity and function. While heat-induced changes of biomolecules are traditionally monitored using optical spectroscopy methods, their conformational changes and unfolding transitions remain challenging to interpret. In the present work, the structural transitions of bovine serum albumin (BSA) in native conditions (100 mM aqueous ammonium acetate) were investigated as a function of the starting solution temperature (T ∼ 23-70 °C) using a temperature-controlled nanoelectrospray ionization source (nESI) coupled to a trapped ion mobility spectrometry-mass spectrometry (TIMS-MS) instrument. The charge state distribution of the monomeric BSA changed from a native-like, narrow charge state ([M + 12H]12+ to [M + 16H]16+ at ∼23 °C) and narrow mobility distribution toward an unfolded-like, broad charge state (up to [M + 46H]46+ at ∼70 °C) and broad mobility distribution. Inspection of the average charge state and collision cross section (CCS) distribution suggested a two-state unfolding transition with a melting temperature Tm ∼ 56 ± 1 °C; however, the inspection of the CCS profiles at the charge state level as a function of the solution temperature showcases at least six structural transitions (T1-T7). If the starting solution concentration is slightly increased (from 2 to 25 μM), this method can detect nonspecific BSA dimers and trimers which dissociate early (Td ∼ 34 ± 1 °C) and may disturb the melting curve of the BSA monomer. In a single experiment, this technology provides a detailed view of the solution, protein structural landscape (mobility vs solution temperature vs relative intensity for each charge state).

Salt-induced LCST-type thermal gelation of methylcellulose: quantifying non-specific interactions via fluctuation theory

What drives the phase separation of water-soluble polymers in the presence of electrolytes was quantified on a molecular scale via statistical thermodynamic fluctuation theory. Quantifying polymer-water and polymer-salt interactions enabled us to identify the dominant interaction for phase separation. As a model system, the lower critical solution temperature (LCST) type thermal gelation of methylcellulose (MC) in aqueous salt solutions was chosen. The Kirkwood-Buff integrals for intermolecular interactions, calculated from the published calorimetric and volumetric data, showed that (1) the accumulation of salts around MC molecules (favourable interaction between salts and MC) inhibits thermal gelation and the depletion of salts from MC (unfavourable interaction between salts and MC) promotes the gelation, and (2) this salt-MC interaction is the dominant factor (50-100 times stronger than the water-MC interaction). This insight from the fluctuation theory is at odds with the age-old consensus regarding the driving force of thermal gelation: water structure change in the presence of salts induces the promotion or inhibition of thermal gelation. However, our conclusion is founded upon the ability of the fluctuation theory to quantify water-MC and salt-MC interaction independently via the Kirkwood-Buff integrals. Flory-Huggins (FH) theory, on the contrary, could not separate these two interactions owing to the lack of a thermodynamic degree of freedom because the lattice solution is assumed to be fully packed. In addition, the dominant contribution from salt depletion poses difficulty for the χ parameter, which is essentially the difference of contact energies. Our approach, requiring calorimetric and volumetric data alone as input, provides a simple and versatile method towards elucidating the effect of cosolvents on biopolymer phase separation of physiological importance.

Larger active site in an ancestral hydroxynitrile lyase increases catalytically promiscuous esterase activity

Hydroxynitrile lyases (HNL's) belonging to the α/β-hydrolase-fold superfamily evolved from esterases approximately 100 million years ago. Reconstruction of an ancestral hydroxynitrile lyase in the α/β-hydrolase fold superfamily yielded a catalytically active hydroxynitrile lyase, HNL1. Several properties of HNL1 differ from the modern HNL from rubber tree (HbHNL). HNL1 favors larger substrates as compared to HbHNL, is two-fold more catalytically promiscuous for ester hydrolysis (p-nitrophenyl acetate) as compared to mandelonitrile cleavage, and resists irreversible heat inactivation to 35 °C higher than for HbHNL. We hypothesized that the x-ray crystal structure of HNL1 may reveal the molecular basis for the differences in these properties. The x-ray crystal structure solved to 1.96-Å resolution shows the expected α/β-hydrolase fold, but a 60% larger active site as compared to HbHNL. This larger active site echoes its evolution from esterases since related esterase SABP2 from tobacco also has a 38% larger active site than HbHNL. The larger active site in HNL1 likely accounts for its ability to accept larger hydroxynitrile substrates. Site-directed mutagenesis of HbHNL to expand the active site increased its promiscuous esterase activity 50-fold, consistent with the larger active site in HNL1 being the primary cause of its promiscuous esterase activity. Urea-induced unfolding of HNL1 indicates that it unfolds less completely than HbHNL (m-value = 0.63 for HNL1 vs 0.93 kcal/mol·M for HbHNL), which may account for the ability of HNL1 to better resist irreversible inactivation upon heating. The structure of HNL1 shows changes in hydrogen bond networks that may stabilize regions of the folded structure.

Chaperone action of a cofactor in protein folding

Previously, we have reported that ATP accelerates the folding and unfolding of Escherichia coli glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which is a glycolytic enzyme utilizing NAD⁺ as a cofactor. Because ATP and NAD⁺ share the ADP moiety, we hypothesized that NAD⁺ also accelerates the folding of GAPDH and that the common structural motif between ATP and NAD⁺ is responsible for the chaperone activity. After confirming that NAD⁺ indeed accelerates the folding of GAPDH, we examined the chaperone activity of the structural fragments of NAD⁺ (ADP, AMP, adenosine, and nicotinamide monophosphate). Our finding showed that ADP and AMP significantly speed up the folding of GAPDH, while adenosine and nicotinamide monophosphate do not. ADP and AMP also dramatically speed up the unfolding of GAPDH by selectively stabilizing a transition state in which GAPDH has a partially unfolded conformation. Similar to the previously reported effect of ATP on the equilibrium unfolding of GAPDH, a partially unfolded intermediate also accumulates in the presence of ADP and AMP. Based on the effect of the structural fragments of NAD⁺ on the folding of GAPDH, we identified that AMP is the structural determinant of the chaperone activity of ATP and NAD⁺ . Also, we propose a plausible model to explain how NAD⁺ accelerates the folding of GAPDH through a stepwise development of molecular interactions with the protein.

Site-specific ubiquitination affects protein energetics and proteasomal degradation

Changes in the cellular environment modulate protein energy landscapes to drive important biology, with consequences for signaling, allostery and other vital processes. The effects of ubiquitination are particularly important because of their potential influence on degradation by the 26S proteasome. Moreover, proteasomal engagement requires unstructured initiation regions that many known proteasome substrates lack. To assess the energetic effects of ubiquitination and how these manifest at the proteasome, we developed a generalizable strategy to produce isopeptide-linked ubiquitin within structured regions of a protein. The effects on the energy landscape vary from negligible to dramatic, depending on the protein and site of ubiquitination. Ubiquitination at sensitive sites destabilizes the native structure and increases the rate of proteasomal degradation. In well-folded proteins, ubiquitination can even induce the requisite unstructured regions needed for proteasomal engagement. Our results indicate a biophysical role of site-specific ubiquitination as a potential regulatory mechanism for energy-dependent substrate degradation.

Cooperative unfolding of a single-layer β-sheet protein, CPAP G-box

CPAP is a centriolar protein and its C-terminal domain, G-box or TCP, has a very unique structure that comprises a single-layer β-sheet without hydrophobic core packing. Here we characterized its biophysical properties, including its stability against chemical denaturation. Interestingly, upon urea-induced equilibrium unfolding, the CPAP G-box showed cooperative unfolding behavior that is the hallmark of globular proteins. We analyzed the m-value, a measure of the cooperative transition, from the urea-induced unfolding and found that the estimated m-value from surface burial upon folding is consistent with the experimental value, supporting the two-state unfolding. Next, we constructed deletion mutants of the terminal β-strands and found that the mutants showed reduced stability. The unique structure and characteristics of CPAP G-box provides an interesting opportunity to observe how the core-less flat β-sheet protein can be folded in solution.

Dissecting the statistical properties of the linear extrapolation method of determining protein stability

The linear extrapolation method to determine protein stability from denaturant-induced unfolding experiments is based on the observation that the free energy of unfolding is often a linear function of the denaturant concentration. The value in the absence of denaturant is then estimated by extrapolation from this linear relationship. Parameters and their confidence intervals are typically estimated by nonlinear least-squares regression. We have compared different methods for calculating confidence intervals and found that a simple method based on linear theory gives accurate results. We have also compared three different parameterizations of the linear extrapolation method and show that the most commonly used form is problematic since the stability and m-value are correlated in the nonlinear least-squares analysis. Parameter correlation can in some cases causes problems in the estimation of confidence intervals and regions and should be avoided when possible. Two alternative parameterizations show much less correlation between parameters.

Searching for Antipneumococcal Targets: Choline-Binding Modules as Phagocytosis Enhancers

Choline-binding proteins (CBPs) from Streptococcus pneumoniae comprise a family of modular polypeptides involved in essential events of this pathogen. They recognize the choline residues present in the teichoic and lipoteichoic acids of the cell wall using the so-called choline-binding modules (CBMs). The importance of CBPs in pneumococcal physiology points to them as novel targets to combat antimicrobial resistances shown by this organism. In this work we have tested the ability of exogenously added CBMs to act as CBP inhibitors by competing with the latter for the binding to the choline molecules in the bacterial surface. First, we carried out a thorough physicochemical characterization of three native CBMs, namely C-LytA, C-Cpl1, and C-CbpD, and assessed their affinity for choline and macromolecular, pneumococcal cell-wall mimics. The interaction with these substrates was evaluated by molecular modeling, analytical ultracentrifugation, surface plasmon resonance, and fluorescence and circular dichroism spectroscopies. Van't Hoff thermal analyses unveiled the existence of one noncanonical choline binding site in each of the C-Cpl1 and C-CbpD proteins, leading in total to 5 ligand-binding sites per dimer and 4 sites per monomer, respectively. Remarkably, the binding affinities of the CBMs do not directly correlate with their native oligomeric state or with the number of choline-binding sites, suggesting that choline recognition by these modules is a complex phenomenon. On the other hand, the exogenous addition of CBMs to pneumococcal planktonic cultures caused extensive cell-chaining probably as a consequence of the inhibition of CBP attachment to the cell wall. This was accompanied by bacterial aggregation and sedimentation, causing an enhancement of bacterial phagocytosis by peritoneal macrophages. In addition, the rational design of an oligomeric variant of a native CBM led to a substantial increase in its antibacterial activity by multivalency effects. These results suggest that CBMs might constitute promising nonlytic antimicrobial candidates based on the natural induction of the host defense system.

Oligomerization of a symmetric β-trefoil protein in response to folding nucleus perturbation

Gene duplication and fusion events in protein evolution are postulated to be responsible for the common protein folds exhibiting internal rotational symmetry. Such evolutionary processes can also potentially yield regions of repetitive primary structure. Repetitive primary structure offers the potential for alternative definitions of critical regions, such as the folding nucleus (FN). In principle, more than one instance of the FN potentially enables an alternative folding pathway in the face of a subsequent deleterious mutation. We describe the targeted mutation of the carboxyl-terminal region of the (internally located) FN of the de novo designed purely-symmetric β-trefoil protein Symfoil-4P. This mutation involves wholesale replacement of a repeating trefoil-fold motif with a "blade" motif from a β-propeller protein, and postulated to trap that region of the Symfoil-4P FN in a nonproductive folding intermediate. The resulting protein (termed "Bladefoil") is shown to be cooperatively folding, but as a trimeric oligomer. The results illustrate how symmetric protein architectures have potentially diverse folding alternatives available to them, including oligomerization, when preferred pathways are perturbed.

De novo design of protein logic gates

The design of modular protein logic for regulating protein function at the posttranscriptional level is a challenge for synthetic biology. Here, we describe the design of two-input AND, OR, NAND, NOR, XNOR, and NOT gates built from de novo-designed proteins. These gates regulate the association of arbitrary protein units ranging from split enzymes to transcriptional machinery in vitro, in yeast and in primary human T cells, where they control the expression of the TIM3 gene related to T cell exhaustion. Designed binding interaction cooperativity, confirmed by native mass spectrometry, makes the gates largely insensitive to stoichiometric imbalances in the inputs, and the modularity of the approach enables ready extension to three-input OR, AND, and disjunctive normal form gates. The modularity and cooperativity of the control elements, coupled with the ability to de novo design an essentially unlimited number of protein components, should enable the design of sophisticated posttranslational control logic over a wide range of biological functions.

Alanine to serine substitutions drive thermal adaptation in a psychrophilic diatom cytochrome c6

In this study, we investigate the thermodynamic mechanisms by which electron transfer proteins adapt to environmental temperature by directly comparing the redox properties and folding stability of a psychrophilic cytochrome c and a mesophilic homolog. Our model system consists of two cytochrome c₆ proteins from diatoms: one adapted specifically to polar environments, the other adapted generally to surface ocean environments. Direct electrochemistry shows that the midpoint potential for the mesophilic homolog is slightly higher at all temperatures measured. Cytochrome c₆ from the psychrophilic diatom unfolds with a melting temperature 10.4 °C lower than the homologous mesophilic cytochrome c₆. Changes in free energy upon unfolding are identical, within error, for the psychrophilic and mesophilic protein; however, the chemical unfolding transition of the psychrophilic cytochrome c₆ is more cooperative than for the mesophilic cytochrome c₆. Substituting alanine residues found in the mesophile with serine found in corresponding positions of the psychrophile demonstrates that burial of the polar serine both decreases the thermal stability and decreases the midpoint potential. The mutagenesis data, combined with differences in the m-value of chemical denaturation, suggest that differences in solvent accessibility of the hydrophobic core underlie the adaptation of cytochrome c₆ to differing environmental temperature.

Friction-Limited Folding of Disulfide-Reduced Monomeric SOD1

The folding reaction of a stable monomeric variant of Cu/Zn superoxide dismutase (mSOD1), an enzyme responsible for the conversion of superoxide free radicals into hydrogen peroxide and oxygen, is known to be among the slowest folding processes that adhere to two-state behavior. The long lifetime, ∼10 s, of the unfolded state presents ample opportunities for the polypeptide chain to transiently sample nonnative structures before the formation of the productive folding transition state. We recently observed the formation of a nonnative structure in a peptide model of the C-terminus of SOD1, a sequence that might serve as a potential source of internal chain friction-limited folding. To test for friction-limited folding, we performed a comprehensive thermodynamic and kinetic analysis of the folding mechanism of mSOD1 in the presence of the viscogens glycerol and glucose. Using a, to our knowledge, novel analysis of the folding reactions, we found the disulfide-reduced form of the protein that exposes the C-terminal sequence, but not its disulfide-oxidized counterpart that protects it, experiences internal chain friction during folding. The sensitivity of the internal friction to the disulfide bond status suggests that one or both of the cross-linked regions play a critical role in driving the friction-limited folding. We speculate that the molecular mechanisms giving rise to the internal friction of disulfide-reduced mSOD1 might play a role in the amyotrophic lateral sclerosis-linked aggregation of SOD1.

Thermal and Chemical Unfolding of a Monoclonal IgG1 Antibody: Application of the Multistate Zimm-Bragg Theory

The thermal unfolding of a recombinant monoclonal antibody IgG1 (mAb) was measured with differential scanning calorimetry (DSC). The DSC thermograms reveal a pretransition at 72°C with an unfolding enthalpy of ΔH_cal ∼200-300 kcal/mol and a main transition at 85°C with an enthalpy of ∼900-1000 kcal/mol. In contrast to small single-domain proteins, mAb unfolding is a complex reaction that is analyzed with the multistate Zimm-Bragg theory. For the investigated mAb, unfolding is characterized by a cooperativity parameter σ ∼6 × 10^-5 and a Gibbs free energy of unfolding of g_nu ∼100 cal/mol per amino acid. The enthalpy of unfolding provides the number of amino acid residues ν participating in the unfolding reaction. On average, ν∼220 ± 50 amino acids are involved in the pretransition and ν∼850 ± 30 in the main transition, accounting for ∼90% of all amino acids. Thermal unfolding was further studied in the presence of guanidineHCl. The chemical denaturant reduces the unfolding enthalpy ΔH_cal and lowers the midpoint temperature T_m. Both parameters depend linearly on the concentration of denaturant. The guanidineHCl concentrations needed to unfold mAb at 25°C are predicted to be 2-3 M for the pretransition and 5-7 M for the main transition, varying with pH. GuanidineHCl binds to mAb with an exothermic binding enthalpy, which partially compensates the endothermic mAb unfolding enthalpy. The number of guanidineHCl molecules bound upon unfolding is deduced from the DSC thermograms. The bound guanidineHCl-to-unfolded amino acid ratio is 0.79 for the pretransition and 0.55 for the main transition. The pretransition binds more denaturant molecules and is more sensitive to unfolding than the main transition. The current study shows the strength of the Zimm-Bragg theory for the quantitative description of unfolding events of large, therapeutic proteins, such as a monoclonal antibody.

Regulation of Photocycle Kinetics of Photoactive Yellow Protein by Modulating Flexibility of the β-Turn

The role of the significant flexibility of the β-turn in photoactive yellow protein (PYP) due to Gly115 was studied. G115A and G115P mutations were observed to accelerate the photocycle and shift the equilibrium between the late photocycle intermediate (pB) and its precursor (pR) toward pR. Thermodynamic analysis of dark-state recovery from pB demonstrated that the transition state (pB^⧧) has a negative change in transition heat capacity, suggesting that an exposed hydrophobic surface of pB is buried in pB^⧧. Fourier transform infrared spectroscopy showed that the structural ensemble of pB is populated by the compact structure in G115P. Taken together, the rigid structure induced by mutation of Gly115 facilitates its transition to pB^⧧, which adopts a substantially more compact structure as opposed to the ensemble-averaged structure of pB. The photocycle kinetics of PYP may be fine-tuned by modulating the flexibility of the 115 loop to activate an appropriate number of transducer proteins.

Tissue-Specificity of Dystrophin-Actin Interactions: Isoform-Specific Thermodynamic Stability and Actin-Binding Function of Tandem Calponin-Homology Domains

Genetic mutations in Duchenne muscular dystrophy (DMD) gene affecting the expression of dystrophin protein lead to a number of muscle disorders collectively called dystrophinopathies. In addition to muscle dystrophin, mutations in brain-specific dystrophin isoforms, in particular those that are expressed in the brain cortex and Purkinje neurons, result in cognitive impairment associated with DMD. These isoforms carry minor variations in the flanking region of the N-terminal actin-binding domain (ABD1) of dystrophin, which is composed of two calponin-homology (CH) domains in tandem. Determining the effect of these sequence variations is critical for understanding the mechanisms that govern varied symptoms of the disease. We studied the impact of differences in the N-terminal flanking region on the structure and function of dystrophin tandem CH domain isoforms. The amino acid changes did not affect the global structure of the protein but drastically affected the thermodynamic stability, with the muscle isoform more stable than the brain and Purkinje isoforms. Actin binding investigated with actin from different sources (skeletal muscle, smooth muscle, cardiac muscle, and platelets) revealed that the muscle isoform binds to filamentous actin (F-actin) with a lower affinity compared to the brain and Purkinje isoforms, and a similar trend was observed with actin from different sources. In addition, all isoforms showed a higher affinity to smooth muscle actin in comparison to actin from other sources. In conclusion, tandem CH domain isoforms might be using minor sequence variations in the N-terminal flanking regions to modulate their thermodynamic stability and actin-binding function, thus leading to specificity in dystrophin-actin interactions in various tissues.

Application of machine learning to predict monomer retention of therapeutic proteins after long term storage

An important aspect of initial developability assessments as well formulation development and selection of therapeutic proteins is the evaluation of data obtained under accelerated stress condition, i.e. at elevated temperatures. We propose the application of artificial neural networks (ANNs) to predict long term stability in real storage condition from accelerated stability studies and other high-throughput biophysical properties e.g. the first apparent temperature of unfolding (T_m). Our models have been trained on therapeutic relevant proteins, including monoclonal antibodies, in various pharmaceutically relevant formulations. Further, we developed network architectures with good prediction power using the least amount of input features, i.e. experimental effort to train the network. This provides an empiric means to highlight the most important parameters in the prediction of real-time protein stability. Further, several models were developed by a different validation means (i.e. leave-one-protein-out cross-validation) to test the robustness and the limitations of our approach. Finally, we apply surrogate machine learning algorithms (e.g. linear regression) to build trust in the ANNs decision making procedure and to highlight the connection between the leading inputs and the outputs.

Describing Complex Structure-Function Relationships in Biomolecules at Equilibrium

One of the great ambitions of structural biology is to describe structure-function relationships quantitatively. Statistical thermodynamics is a powerful, general tool for computing the behavior of biological macromolecules at equilibrium because it establishes a direct link between structure and function. Complex behavior emerges as equilibria of multiple reactions are coupled. Analytical treatment of linked equilibria scales poorly with increasing numbers of reactions and states as the algebraic constructs rapidly become unwieldy. We therefore developed a generalizable, but straightforward computational method to handle arbitrarily complex systems. To demonstrate this approach, we collected a multidimensional fluorescence landscape of an engineered fluorescent glucose biosensor and showed that its features could be modeled with ten intricately linked ligand-binding and conformational exchange reactions. This protein represents a minimalist model of sufficient complexity to encompass fundamental biomolecular structure-function relationships: two-state and multistate conformational ensembles, conformational hierarchies, osmolytes, coupling between different binding sites and coupling between ligand binding and conformational change. The successful fit of this complex, multifaceted system demonstrates generality of the method.

Investigation of guanidinium chloride-induced unfolding pathway of sphingosine kinase 1

Sphingosine kinase 1 (SphK1) is a lipid kinase which plays vital role in the regulation of varieties of biological processes including, cell growth, apoptosis and mitogenesis. In the present study, we investigated the guanidinium chloride (GdmCl)-induced denaturation of SphK1 at pH 8.0 and 25 °C using two different spectroscopic probes, i.e., mean residue ellipticity at 222 nm ([θ]₂₂₂) and fluorescence emission maxima (λ_max). A significant overlap between the transition curves obtained from both the spectral properties indicate that GdmCl-induced unfolding of SphK1 follows two-state process i.e., Native (N) ⇌ Denatured (D) state. Interestingly, a visible protein aggregation was observed at low concentrations of GdmCl ([GdmCl] ≤ 1.5 M). The analysis of transition curves was done to estimate the thermodynamic parameters associated with the stability of SphK1. To complement our experimental findings, 100 ns molecular dynamics (MD) simulations were performed. Spectroscopic studies together with MD simulations provided mechanistic insights of unfolding pathway of SphK1 along with its stability parameters.

Advancing Therapeutic Protein Discovery and Development through Comprehensive Computational and Biophysical Characterization

Therapeutic protein candidates should exhibit favorable properties that render them suitable to become drugs. Nevertheless, there are no well-established guidelines for the efficient selection of proteinaceous molecules with desired features during early stage development. Such guidelines can emerge only from a large body of published research that employs orthogonal techniques to characterize therapeutic proteins in different formulations. In this work, we share a study on a diverse group of proteins, including their primary sequences, purity data, and computational and biophysical characterization at different pH and ionic strength. We report weak linear correlations between many of the biophysical parameters. We suggest that a stability comparison of diverse therapeutic protein candidates should be based on a computational and biophysical characterization in multiple formulation conditions, as the latter can largely determine whether a protein is above or below a certain stability threshold. We use the presented data set to calculate several stability risk scores obtained with an increasing level of analytical effort and show how they correlate with protein aggregation during storage. Our work highlights the importance of developing combined risk scores that can be used for early stage developability assessment. We suggest that such scores can have high prediction accuracy only when they are based on protein stability characterization in different solution conditions.

Determining the Free Energies of Outer Membrane Proteins in Lipid Bilayers

The thermodynamic stabilities of membrane proteins are of fundamental interest to provide a biophysical description of their structure-function relationships because energy determines conformational populations. In addition, structure-energy relationships can be exploited in membrane protein design and in synthetic biology. To determine the thermodynamic stability of a membrane protein, it is not sufficient to be able to unfold and refold the molecule: establishing path independence of this reaction is essential. Here we describe the procedures required to measure and verify path independence for the folding of outer membrane proteins in large unilamellar vesicles.

The Effect of Proline cis-trans Isomerization on the Folding of the C-Terminal SH2 Domain from p85

SH2 domains are protein domains that modulate protein-protein interactions through a specific interaction with sequences containing phosphorylated tyrosines. In this work, we analyze the folding pathway of the C-terminal SH2 domain of the p85 regulatory subunit of the protein PI3K, which presents a proline residue in a cis configuration in the loop between the βE and βF strands. By employing single and double jump folding and unfolding experiments, we demonstrate the presence of an on-pathway intermediate that transiently accumulates during (un)folding. By comparing the kinetics of folding of the wild-type protein to that of a site-directed variant of C-SH2 in which the proline was replaced with an alanine, we demonstrate that this intermediate is dictated by the peptidyl prolyl cis-trans isomerization. The results are discussed in the light of previous work on the effect of peptidyl prolyl cis-trans isomerization on folding events.

Thermodynamics of amyloid fibril formation from chemical depolymerization

Amyloid fibrils are homo-molecular protein polymers that play an important role in disease and biological function. While much is known about their kinetics and mechanisms of formation, the origin and magnitude of their thermodynamic stability has received significantly less attention. This is despite the fact that the thermodynamic stability of amyloid fibrils is an important determinant of their lifetimes and processing in vivo. Here we use depolymerization by chemical denaturants of amyloid fibrils of two different proteins (PI3K-SH3 and glucagon) at different concentrations and show that the previously applied isodesmic linear polymerization model is an oversimplification that does not capture the concentration dependence of chemical depolymerization of amyloid fibrils. We show that cooperative polymerization, which is compatible with the picture of amyloid formation as a nucleated polymerization process, is able to quantitatively describe the thermodynamic data. We use this combined experimental and conceptual framework in order to probe the ionic strength dependence of amyloid fibril stability. In combination with previously published data on the ionic strength dependence of amyloid fibril growth kinetics, our results provide strong evidence for the product-like nature of the transition state of amyloid fibril growth.

Thermal and Chemical Unfolding of Lysozyme. Multistate Zimm-Bragg Theory Versus Two-State Model

Thermal and chemical unfolding of lysozyme in the presence of the guanidine HCl denaturant is a model system to compare the conventional two-state model of protein unfolding with the multistate Zimm-Bragg theory. The two-state model is shown to be the noncooperative limit of the Zimm-Bragg theory. In particular, the Zimm-Bragg theory provides a molecular interpretation of the empirical linear extrapolation method (LEM) of the two-state model. Differential scanning calorimetry (DSC) experiments reported in the literature are analyzed with both methods. Lysozyme unfolding is associated with a large endothermic enthalpy that decreases significantly upon addition of guanidine HCl. In contrast, the Gibbs free energy of unfolding is small, negative, and independent of the guanidine HCl concentration, contradicting, in part, the conclusions of the LEM. The unfolding enthalpy is compensated by an even larger entropy term. The multistate Zimm-Bragg theory predicts a larger conformational enthalpy and a smaller Gibbs free energy than the two-state model. The Zimm-Bragg theory provides the protein cooperativity parameter, the average length of independently folding protein domains, and the Gibbs free energy of unfolding of individual amino acid residues. Guanidine HCl binding to lysozyme is exothermic and counteracts the endothermic unfolding enthalpy. The number of bound denaturant molecules is determined from the decrease in enthalpy and is extrapolated to the guanidine HCl-to-amino acid stoichiometry at complete lysozyme unfolding. Chemical unfolding isotherms measured with circular dichroism (CD) spectroscopy are analyzed with both models. The chemical Zimm-Bragg theory is a cooperative molecular model, yielding the guanidine HCl binding constant and the protein cooperativity parameter. It allows a quantitative comparison between thermal and chemical protein unfolding. The two reactions have almost identical changes in Gibbs free energy. However, thermal unfolding is significantly more cooperative than chemical unfolding. Finally, distinct differences are observed in thermal unfolding between DSC and CD spectroscopy.

Human aldose reductase unfolds through an intermediate

Background: Human aldose reductase (hAR) is the first and rate-limiting enzyme of the polyol pathway. For the development of secondary complications of diabetes in chronic hyperglycemic conditions, one of the critical factors is the increased flux of glucose through the polyol pathway.  Due to this clinical implication, hAR attracted considerable attention from the drug discovery perspective. In spite of extensive characterization in the context of biochemical and structural aspects, we know very little about the unfolding behavior of hAR. This study reports equilibrium unfolding studies of hAR.

Methods: We carried out thermal denaturation and chemical-induced equilibrium unfolding studies of hAR monitored by circular dichroism and fluorescence spectroscopy. 

Results: Thermal denaturation studies presented a classical picture of two-state unfolding from native to the denatured state. The data was used to derive thermodynamic parameters and study the thermostability of hAR. Chemical induced equilibrium unfolding studies led us to discover an intermediate state, which gets populated at 3.5-4.0 M and 0.7-2.0 M of urea and GuHCl, respectively. Thermodynamic parameters derived from chemical-induced unfolding are in agreement with those obtained from thermal denaturation of hAR.

Conclusion: This study revealed that aldose reductase unfolds from native to the unfolded state via an intermediate. Assessment of the thermodynamic stability of native, intermediate, and unfolded states shows that significant energy barriers separate these states, which ensures the cooperativity of unfolding. As hAR functions in cells that are under osmotic and oxidative stress, these in vitro findings may have implications for its native conformation under the physiological state.

Stabilization of backbone-circularized protein is attained by synergistic gains in enthalpy of folded structure and entropy of unfolded structure

Backbone circularization is an effective technique for protein stabilization. Here, we investigated the effect of a connector, an engineered segment that connects two protein termini, on the conformational stability of previously designed circularized variants of granulocyte colony-stimulating factor (G-CSF). Heat tolerance and chemical denaturation analyses revealed that aggregation resistance and thermodynamic stability of the circularized variants were superior to those of linear G-CSF. Crystal structure and molecular dynamics (MD) simulation of the most thermodynamically stable variant (C166) revealed a high number of intramolecular hydrogen bonds in both the connector region and Helix D adjacent to the connector region in the folded structure. MD simulations and theoretical calculations involving different force fields indicated a reduction in the main chain entropy of C166 in the unfolded state and increase in the intramolecular hydrogen bond energy of C166 in the folded structure. Although backbone circularization is usually considered to alter chain entropy of the unfolded state, the data indicated that it could also improve the conformational enthalpy of the folded state. Further structural examination of the connector region confirmed that protein design based on a statistical analysis of local structures is an effective approach for predicting an optimum connector length to improve the conformational stability of backbone-circularized proteins. Protein design using backbone circularization with an optimum connector length will be useful for the development of effective and safe protein therapeutics. DATABASE: Structural data are available in Protein Data Bank under the accession number 5ZO6.

Methionine in a protein hydrophobic core drives tight interactions required for assembly of spider silk

Web spiders connect silk proteins, so-called spidroins, into fibers of extraordinary toughness. The spidroin N-terminal domain (NTD) plays a pivotal role in this process: it polymerizes spidroins through a complex mechanism of dimerization. Here we analyze sequences of spidroin NTDs and find an unusually high content of the amino acid methionine. We simultaneously mutate all methionines present in the hydrophobic core of a spidroin NTD from a nursery web spider’s dragline silk to leucine. The mutated NTD is strongly stabilized and folds at the theoretical speed limit. The structure of the mutant is preserved, yet its ability to dimerize is substantially impaired. We find that side chains of core methionines serve to mobilize the fold, which can thereby access various conformations and adapt the association interface for tight binding. Methionine in a hydrophobic core equips a protein with the capacity to dynamically change shape and thus to optimize its function.

Spider silk is of interest in material science research. Here the authors show that the tight binding of a spider silk protein domain relies on the amino acid methionine, which is abundant in the domain core where it facilitates dynamic shape adaption of the binding interface.

Folding and Dynamics Are Strongly pH-Dependent in a Psychrophile Frataxin

Protein dynamics, folding, and thermodynamics represent a central aspect of biophysical chemistry. pH, temperature, and denaturant perturbations inform our understanding of diverse contributors to stability and rates. In this work, we performed a thermodynamic analysis using a combined experimental and computational approach to gain insights into the role of electrostatics in the folding reaction of a psychrophile frataxin variant from Psychromonas ingrahamii. This folding reaction is strongly modulated by pH with a single, narrow, and well-defined transition state with ∼80% compactness, ∼70% electrostatic interactions, and ∼60% hydration shell compared to the native state (αD = 0.82, αH = 0.67, and αΔCp = 0.59). Our results are best explained by a two-proton/two-state model with very different pKa values of the native and denatured states (∼5.5 and ∼8.0, respectively). As a consequence, the stability strongly increases from pH 8.0 to 6.0 (|ΔΔG°| = 5.2 kcal mol-1), mainly because of a decrease in the TΔS°. Variation of ΔH° and ΔS° at pH below 7.0 is dominated by a change in ΔHf⧧ and ΔSf⧧, while at pH above 7.0, it is governed by ΔHu⧧ and ΔSu⧧. Molecular dynamics simulations showed that these pH modulations could be explained by the fluctuations of two regions, rich in electrostatic contacts, whose dynamics are pH-dependent and motions are strongly correlated. Results presented herein contribute to the understanding of the stability and dynamics of this frataxin variant, pointing to an intrinsic feature of the family topology to support different folding mechanisms.

The Cytoplasm-Entry Domain of Antibacterial CdiA Is a Dynamic α-Helical Bundle with Disulfide-Dependent Structural Features

Many Gram-negative bacterial species use contact-dependent growth inhibition (CDI) systems to compete with neighboring cells. CDI+ strains express cell-surface CdiA effector proteins, which carry a toxic C-terminal region (CdiA-CT) that is cleaved from the effector upon transfer into the periplasm of target bacteria. The released CdiA-CT consists of two domains. The C-terminal domain is typically a nuclease that inhibits cell growth, and the N-terminal "cytoplasm-entry" domain mediates toxin translocation into the target-cell cytosol. Here, we use NMR and circular dichroism spectroscopic approaches to probe the structure, stability, and dynamics of the cytoplasm-entry domain from Escherichia coli STEC_MHI813. Chemical shift analysis reveals that the CdiA-CTMHI813 entry domain is composed of a C-terminal helical bundle and a dynamic N-terminal region containing two disulfide linkages. Disruption of the disulfides by mutagenesis or chemical reduction destabilizes secondary structure over the N-terminus, but has no effect on the C-terminal helices. Although critical for N-terminal structure, the disulfides have only modest effects on global thermodynamic stability, and the entry domain exhibits characteristics of a molten globule. We find that the disulfides form in vivo as the entry domain dwells in the periplasm of inhibitor cells prior to target-cell recognition. CdiA-CTMHI813 variants lacking either disulfide still kill target bacteria, but disruption of both bonds abrogates growth inhibition activity. We propose that the entry domain's dynamic structural features are critical for function. In its molten globule-like state, the domain resists degradation after delivery, yet remains pliable enough to unfold for membrane translocation.

Protein stability engineering insights revealed by domain-wide comprehensive mutagenesis

The accurate prediction of protein stability upon sequence mutation is an important but unsolved challenge in protein engineering. Large mutational datasets are required to train computational predictors, but traditional methods for collecting stability data are either low-throughput or measure protein stability indirectly. Here, we develop an automated method to generate thermodynamic stability data for nearly every single mutant in a small 56-residue protein. Analysis reveals that most single mutants have a neutral effect on stability, mutational sensitivity is largely governed by residue burial, and unexpectedly, hydrophobics are the best tolerated amino acid type. Correlating the output of various stability-prediction algorithms against our data shows that nearly all perform better on boundary and surface positions than for those in the core and are better at predicting large-to-small mutations than small-to-large ones. We show that the most stable variants in the single-mutant landscape are better identified using combinations of 2 prediction algorithms and including more algorithms can provide diminishing returns. In most cases, poor in silico predictions were tied to compositional differences between the data being analyzed and the datasets used to train the algorithm. Finally, we find that strategies to extract stabilities from high-throughput fitness data such as deep mutational scanning are promising and that data produced by these methods may be applicable toward training future stability-prediction tools.

Isotonic concentrations of excipients control the dimerization rate of a therapeutic immunoglobulin G1 antibody during refrigerated storage based on their rank order of native-state interaction

Inert co-solutes, or excipients, are often included in protein biologic formulations to adjust the tonicity of liquid dosage forms intended for subcutaneous delivery. Despite the low concentration of their use, many of these excipients alter protein-protein interactions such as dimerization and aggregation rates of high concentration monoclonal antibody (mAb) therapeutics to varying extents during long-term refrigerated clinical storage, challenging the formulation scientist to make informed excipient selections at the earliest stages of development when protein supply and time are often limited. The objectives of this study were to better understand how isotonic concentrations of excipients influence the dimerization rates of a model mAb stored at refrigerated and room temperatures and explore protein sparing biophysical methods capable of predicting this dependence. Despite their prevalence of use in the biopharmaceutical industry, methods for assessing conformational stability such differential scanning calorimetry and isothermal equilibrium unfolding showed little predictive power and we highlight some of the assumptions and technical challenges of their use with mAbs. Conversely, measures of colloidal stability of the native-state such as preferential interaction coefficients measured by vapor pressure osmometry and solubility assessed by polyethylene-glycol induced precipitation correlated reasonably well with the mAb dimerization data and are most consistent with the excipients tested minimizing dimerization by interacting favorably with the residues comprising the protein-protein association interface.

Context-Dependent Energetics of Loop Extensions in a Family of Tandem-Repeat Proteins

Consensus-designed tetratricopeptide repeat proteins are highly stable, modular proteins that are strikingly amenable to rational engineering. They therefore have tremendous potential as building blocks for biomaterials and biomedicine. Here, we explore the possibility of extending the loops between repeats to enable further diversification, and we investigate how this modification affects stability and folding cooperativity. We find that extending a single loop by up to 25 residues does not disrupt the overall protein structure, but, strikingly, the effect on stability is highly context-dependent: in a two-repeat array, destabilization is relatively small and can be accounted for purely in entropic terms, whereas extending a loop in the middle of a large array is much more costly because of weakening of the interaction between the repeats. Our findings provide important and, to our knowledge, new insights that increase our understanding of the structure, folding, and function of natural repeat proteins and the design of artificial repeat proteins in biotechnology.

Consensus sequence design as a general strategy to create hyperstable, biologically active proteins

Consensus sequence design offers a promising strategy for designing proteins of high stability while retaining biological activity since it draws upon an evolutionary history in which residues important for both stability and function are likely to be conserved. Although there have been several reports of successful consensus design of individual targets, it is unclear from these anecdotal studies how often this approach succeeds and how often it fails. Here, we attempt to assess generality by designing consensus sequences for a set of six protein families with a range of chain lengths, structures, and activities. We characterize the resulting consensus proteins for stability, structure, and biological activities in an unbiased way. We find that all six consensus proteins adopt cooperatively folded structures in solution. Strikingly, four of six of these consensus proteins show increased thermodynamic stability over naturally occurring homologs. Each consensus protein tested for function maintained at least partial biological activity. Although peptide binding affinity by a consensus-designed SH3 is rather low, K _m values for consensus enzymes are similar to values from extant homologs. Although consensus enzymes are slower than extant homologs at low temperature, they are faster than some thermophilic enzymes at high temperature. An analysis of sequence properties shows consensus proteins to be enriched in charged residues, and rarified in uncharged polar residues. Sequence differences between consensus and extant homologs are predominantly located at weakly conserved surface residues, highlighting the importance of these residues in the success of the consensus strategy.

Using Single-Molecule Chemo-Mechanical Unfolding to Simultaneously Probe Multiple Structural Parameters in Protein Folding

While single-molecule force spectroscopy has greatly advanced the study of protein folding, there are limitations to what can be learned from studying the effect of force alone. We developed a novel technique, chemo-mechanical unfolding, that combines multiple perturbants—force and chemical denaturant—to more fully characterize the folding process by simultaneously probing multiple structural parameters—the change in end-to-end distance, and solvent accessible surface area. Here, we describe the theoretical background, experimental design, and data analysis for chemo-mechanical unfolding experiments probing protein folding thermodynamics and kinetics. This technique has been applied to characterize parallel protein folding pathways, the protein denatured state, protein folding on the ribosome, and protein folding intermediates.