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

Biophysical description of Bromosulfophthalein interaction with the 28-kDa glutathione transferase from Schistosoma japonicum

Glutathione transferases (GSTs) are major detoxification enzymes vital for the survival and reproduction of schistosomes during infection in humans. Schistosoma encode two GST isoenzymes, the 26- and 28-kDa isoforms, that show different substrate specificities and cellular localisations. Bromosulfophthalein (BSP) has been identified and characterised as a potent 26-kDa Schistosoma japonicum GST (Sj26GST) inhibitor with an anthelmintic potential. This study describes the structure, function, and ligandin properties of the 28-kDa Schistosoma japonicum GST (Sj28GST) towards BSP. Enzyme kinetics show that BSP is a potent enzyme inhibitor, with a specific activity decreases from 60.4 µmol/min/mg to 0.0742 µmol/min/mg and an IC₅₀ in the micromolar range of 0.74 µM. Far-UV circular dichroism confirmed that purified Sj28GST follows a typical GST fold, which is predominantly alpha-helical. Fluorescence spectroscopy suggests that BSP binding occurs at a site distinct from the glutathione-binding site (G-site); however, the binding does not alter the local G-site environment. Isothermal titration calorimetry studies show that the binding of BSP to Sj28GST is exergonic (∆G°= -33 kJ/mol) and enthalpically-driven, with a stoichiometry of one BSP per dimer. The stability of Sj28GST (∆G_(H2O) = 4.7 kcal/mol) is notably lower than Sj26GST, owing to differences in the enzyme's dimeric interfaces. We conclude that Sj28GST shares similar biophysical characteristics with Sj26GST based on its kinetic properties and susceptibility to low concentrations of BSP. The study supports the potential benefits of re-purposing BSP as a potential drug or prodrug to mitigate the scourge of schistosomiasis.

Globin X: A highly stable intrinsically hexacoordinate globin

Several novel members of the vertebrate globin family were recently discovered with unique structural features that are not found in traditional penta-coordinate globins. Here we combine structural tools to better understand and recognize molecular determinants that contribute to the stability of hexacoordinate globin X (GbX) from Danio rerio (zebrafish). pH-induced unfolding data indicates increased stability of GbX with pH_mid of 1.9 ± 0.1 for met GbXWT, 2.4 ± 0.1 for met GbXC65A, and 3.4 ± 0.1 for GbXH90V. These results are in good agreement with GbX unfolding experiments using GuHCl, where a ΔG_unf 13.8 ± 2.5 kcal mol^-1 and 16.3 ± 2.6 kcal mol^-1 are observed for metGbXWT, and metGbXC65A constructs, respectively, and diminished stability is measured for GbXH90V, ΔG_unf = 9.5 ± 3.6 kcal mol^-1. The metGbXWT and metGbXC65A also exhibit high thermal stability (melting points of 118 °C and 107 °C, respectively). Native ion mobility - mass spectrometry (IM-MS) experiments showed a narrow charge state distribution (9-12+) characteristics of a native, structured protein; a single mobility band was observed for the native states. Collision induced unfolding IM-MS experiments showed a two-state transition, in good agreement with the solution studies. GbXWT retains the heme over a wide range of charge states, suggesting strong interactions between the prosthetic group and the apoprotein. The above results indicate that in addition to the disulfide bond and the heme iron hexa-coordination, other structural determinants enhance stability of this protein.

Conformational Stability and Denaturation Processes of Proteins Investigated by Electrophoresis under Extreme Conditions

The functional structure of proteins results from marginally stable folded conformations. Reversible unfolding, irreversible denaturation, and deterioration can be caused by chemical and physical agents due to changes in the physicochemical conditions of pH, ionic strength, temperature, pressure, and electric field or due to the presence of a cosolvent that perturbs the delicate balance between stabilizing and destabilizing interactions and eventually induces chemical modifications. For most proteins, denaturation is a complex process involving transient intermediates in several reversible and eventually irreversible steps. Knowledge of protein stability and denaturation processes is mandatory for the development of enzymes as industrial catalysts, biopharmaceuticals, analytical and medical bioreagents, and safe industrial food. Electrophoresis techniques operating under extreme conditions are convenient tools for analyzing unfolding transitions, trapping transient intermediates, and gaining insight into the mechanisms of denaturation processes. Moreover, quantitative analysis of electrophoretic mobility transition curves allows the estimation of the conformational stability of proteins. These approaches include polyacrylamide gel electrophoresis and capillary zone electrophoresis under cold, heat, and hydrostatic pressure and in the presence of non-ionic denaturing agents or stabilizers such as polyols and heavy water. Lastly, after exposure to extremes of physical conditions, electrophoresis under standard conditions provides information on irreversible processes, slow conformational drifts, and slow renaturation processes. The impressive developments of enzyme technology with multiple applications in fine chemistry, biopharmaceutics, and nanomedicine prompted us to revisit the potentialities of these electrophoretic approaches. This feature review is illustrated with published and unpublished results obtained by the authors on cholinesterases and paraoxonase, two physiologically and toxicologically important enzymes.

Loss of stability and unfolding cooperativity in hPGK1 upon gradual structural perturbation of its N-terminal domain hydrophobic core

Phosphoglycerate kinase has been a model for the stability, folding cooperativity and catalysis of a two-domain protein. The human isoform 1 (hPGK1) is associated with cancer development and rare genetic diseases that affect several of its features. To investigate how mutations affect hPGK1 folding landscape and interaction networks, we have introduced mutations at a buried site in the N-terminal domain (F25 mutants) that either created cavities (F25L, F25V, F25A), enhanced conformational entropy (F25G) or introduced structural strain (F25W) and evaluated their effects using biophysical experimental and theoretical methods. All F25 mutants folded well, but showed reduced unfolding cooperativity, kinetic stability and altered activation energetics according to the results from thermal and chemical denaturation analyses. These alterations correlated well with the structural perturbation caused by mutations in the N-terminal domain and the destabilization caused in the interdomain interface as revealed by H/D exchange under native conditions. Importantly, experimental and theoretical analyses showed that these effects are significant even when the perturbation is mild and local. Our approach will be useful to establish the molecular basis of hPGK1 genotype-phenotype correlations due to phosphorylation events and single amino acid substitutions associated with disease.

Evolution avoids a pathological stabilizing interaction in the immune protein S100A9

Stability constrains evolution. While much is known about constraints on destabilizing mutations, less is known about the constraints on stabilizing mutations. We recently identified a mutation in the innate immune protein S100A9 that provides insight into such constraints. When introduced into human S100A9, M63F simultaneously increases the stability of the protein and disrupts its natural ability to activate Toll-like receptor 4. Using chemical denaturation, we found that M63F stabilizes a calcium-bound conformation of hS100A9. We then used NMR to solve the structure of the mutant protein, revealing that the mutation distorts the hydrophobic binding surface of hS100A9, explaining its deleterious effect on function. Hydrogen-deuterium exchange (HDX) experiments revealed stabilization of the region around M63F in the structure, notably Phe37. In the structure of the M63F mutant, the Phe37 and Phe63 sidechains are in contact, plausibly forming an edge-face π-stack. Mutating Phe37 to Leu abolished the stabilizing effect of M63F as probed by both chemical denaturation and HDX. It also restored the biological activity of S100A9 disrupted by M63F. These findings reveal that Phe63 creates a molecular staple with Phe37 that stabilizes a nonfunctional conformation of the protein, thus disrupting function. Using a bioinformatic analysis, we found that S100A9 proteins from different organisms rarely have Phe at both positions 37 and 63, suggesting that avoiding a pathological stabilizing interaction indeed constrains S100A9 evolution. This work highlights an important evolutionary constraint on stabilizing mutations, namely, that they must avoid inappropriately stabilizing nonfunctional protein conformations.

Computationally Guided Redesign of a Heme-free Cytochrome with Native-like Structure and Stability

Metals can play key roles in stabilizing protein structures, but ensuring their proper incorporation is a challenge when a metalloprotein is overexpressed in a non-native cellular environment. Here, we have used computational protein design tools to redesign cytochrome b₅₆₂ (cyt b₅₆₂), which relies on the binding of its heme cofactor to achieve its proper fold, into a stable, heme-free protein. The resulting protein, ApoCyt, features only four mutations and no metal-ligand or covalent bonds, yet displays improved stability over cyt b₅₆₂. Mutagenesis studies and X-ray crystal structures reveal that the increase in stability is due to the computationally prescribed mutations, which stabilize the protein fold through a combination of hydrophobic packing interactions, hydrogen bonds, and cation-π interactions. Upon installation of the relevant mutations, ApoCyt is capable of assembling into previously reported, cytochrome-based trimeric and tetrameric assemblies, demonstrating that ApoCyt retains the structure and assembly properties of cyt b₅₆₂. The successful design of ApoCyt therefore enables further functional diversification of cytochrome-based assemblies and demonstrates that structural metal cofactors can be replaced by a small number of well-designed, non-covalent interactions.

Threading single proteins through pores to compare their energy landscapes

Protein function correlates with its structural dynamics. While theoretical approaches to studying protein energy landscapes are well developed, experimental methods that enable probing these landscapes of proteins remain challenging. We used solid-state nanopores to study the translocation behavior of three mutants of a helix bundle protein and quantified the number of energetically accessible conformational states for each mutant. We found that a slower-folding mutant with access to more conformational states translocates faster than a faster-folding mutant with a smaller number of accessible states, suggesting that ease of folding and ease of translocation are at odds in this case.

Translocation of proteins is correlated with structural fluctuations that access conformational states higher in free energy than the folded state. We use electric fields at the solid-state nanopore to control the relative free energy and occupancy of different protein conformational states at the single-molecule level. The change in occupancy of different protein conformations as a function of electric field gives rise to shifts in the measured distributions of ionic current blockades and residence times. We probe the statistics of the ionic current blockades and residence times for three mutants of the λ-repressor family in order to determine the number of accessible conformational states of each mutant and evaluate the ruggedness of their free energy landscapes. Translocation becomes faster at higher electric fields when additional flexible conformations are available for threading through the pore. At the same time, folding rates are not correlated with ease of translocation; a slow-folding mutant with a low-lying intermediate state translocates faster than a faster-folding two-state mutant. Such behavior allows us to distinguish among protein mutants by selecting for the degree of current blockade and residence time at the pore. Based on these findings, we present a simple free energy model that explains the complementary relationship between folding equilibrium constants and translocation rates.

Thermodynamics for the Self-Assembly of Alkylated Peptides

Self-assembling peptides form aggregates with various nanostructures such as spheres, sheets, and fibers and have potential applications in nanomedicine and drug delivery. The alkylation of peptides is a promising strategy for controlling the self-assembly of peptides. In this study, we investigated the thermodynamic properties associated with the aggregation of alkyl-chain-modified self-assembling peptides. The tripeptide sequence, KYF, which has been reported to form fibrous aggregates via self-assembly, was modified with various fatty acids at the N-terminus. The fibrous morphology of the aggregates was observed by transmission electron microscopy and atomic force microscopy. Thioflavin T fluorescence and circular dichroism spectroscopy revealed the formation of β-sheet structures. The critical micelle concentration and its temperature dependence were determined to obtain the thermodynamic parameters for aggregation. The results showed that the aggregation was an entropy-driven process at low temperatures, whereas it was enthalpy-driven at high temperatures. The negative heat capacity changes for aggregation suggested that hydrophobic interactions were the major driving force for self-assembly. Other entropic and enthalpic interactions were also contributed in part to the self-assembly. We individually identified the contributions of the peptide and alkyl chain moiety to the self-assembly. These contributions can be explained by the theoretical values for the self-assembly of each component. The results of this study provide fundamental insights into the design of self-associating peptides.

Conformational stability of ageritin, a metal binding ribotoxin-like protein of fungal origin

Ageritin is a ribotoxin-like protein of biotechnological interest, belonging to a family of ribonucleases from edible mushrooms. Its enzymatic activity is explicated through the hydrolysis of a single phosphodiester bond, located in the sarcin/ricin loop of ribosomes. Unlike other ribotoxins, ageritin activity requires divalent cations (Zn²⁺). Here we investigated the conformational stability of ageritin in the pH range 4.0-7.4, using calorimetric and spectroscopic techniques. We observed a high protein thermal stability at all pHs with a denaturation temperature of 78 °C. At pH 5.0 we calculated a value of 36 kJ mol^-1 for the unfolding Gibbs energy at 25 °C. We also analysed the thermodynamic and catalytic behaviour of S-pyridylethylated form, obtained by alkylating the single Cys18 residue, which is predicted to bind Zn²⁺. We show that this form possesses the same activity and structure of ageritin, but lower stability. In fact, the corresponding values of 52 °C and 14 kJ mol^-1 were found. Conservation of activity is consistent with the location of alkylation site on the opposite site of the catalytic site cleft. Inasmuch as Cys18 is part of a structurally stabilizing zinc-binding site, disrupted by cysteine alkylation, our results point to an important role of metal ions in ageritin stability.

Arsenal of nanobodies shows broad-spectrum neutralization against SARS-CoV-2 variants of concern in vitro and in vivo in hamster models

Nanobodies offer several potential advantages over mAbs for the control of SARS-CoV-2. Their ability to access cryptic epitopes conserved across SARS-CoV-2 variants of concern (VoCs) and feasibility to engineer modular, multimeric designs, make these antibody fragments ideal candidates for developing broad-spectrum therapeutics against current and continually emerging SARS-CoV-2 VoCs. Here we describe a diverse collection of 37 anti-SARS-CoV-2 spike glycoprotein nanobodies extensively characterized as both monovalent and IgG Fc-fused bivalent modalities. The nanobodies were collectively shown to have high intrinsic affinity; high thermal, thermodynamic and aerosolization stability; broad subunit/domain specificity and cross-reactivity across existing VoCs; wide-ranging epitopic and mechanistic diversity and high and broad in vitro neutralization potencies. A select set of Fc-fused nanobodies showed high neutralization efficacies in hamster models of SARS-CoV-2 infection, reducing viral burden by up to six orders of magnitude to below detectable levels. In vivo protection was demonstrated with anti-RBD and previously unreported anti-NTD and anti-S2 nanobodies. This collection of nanobodies provides a potential therapeutic toolbox from which various cocktails or multi-paratopic formats could be built to combat multiple SARS-CoV-2 variants.

Isolation and extensive characterization of a collection of 37 anti-SARS-CoV-2 spike glycoprotein nanobodies show broad neutralization efficacies in vitro and in vivo in a hamster model of SARS-CoV-2 infection.

Dynamics of integrin α5β1, fibronectin, and their complex reveal sites of interaction and conformational change

Integrin α5β1 mediates cell adhesion to the extracellular matrix by binding fibronectin (Fn). Selectivity for Fn by α5β1 is achieved through recognition of an RGD motif in the 10th type III Fn domain (Fn10) and the synergy site in the ninth type III Fn domain (Fn9). However, details of the interaction dynamics are unknown. Here, we compared synergy-site and Fn-truncation mutations for their α5β1-binding affinities and stabilities. We also interrogated binding of the α5β1 ectodomain headpiece fragment to Fn using hydrogen-deuterium exchange (HDX) mass spectrometry to probe binding sites and sites of integrin conformational change. Our results suggest the synergistic effect of Fn9 requires both specific residues and a folded domain. We found some residues considered important for synergy are required for stability. Additionally, we show decreases in fibronectin HDX are localized to a synergy peptide containing contacting residues in two β-strands, an intervening loop in Fn9, and the RGD-containing loop in Fn10, indicative of binding sites. We also identified binding sites in the α5-subunit β-propeller domain for the Fn9 synergy site and in the β1-subunit βI domain for Fn10 based on decreases in α5β1 HDX. Interestingly, the dominant effect of Fn binding was an increase in α5β1 deuterium exchange distributed over multiple sites that undergo changes in conformation or solvent accessibility and appear to be sites where energy is stored in the higher-energy, open-integrin conformation. Together, our results highlight regions important for α5β1 binding to Fn and dynamics associated with this interaction.

Tuning Protein Hydrogel Mechanics through Modulation of Nanoscale Unfolding and Entanglement in Postgelation Relaxation

Globular folded proteins are versatile nanoscale building blocks to create biomaterials with mechanical robustness and inherent biological functionality due to their specific and well-defined folded structures. Modulating the nanoscale unfolding of protein building blocks during network formation (in situ protein unfolding) provides potent opportunities to control the protein network structure and mechanics. Here, we control protein unfolding during the formation of hydrogels constructed from chemically cross-linked maltose binding protein using ligand binding and the addition of cosolutes to modulate protein kinetic and thermodynamic stability. Bulk shear rheology characterizes the storage moduli of the bound and unbound protein hydrogels and reveals a correlation between network rigidity, characterized as an increase in the storage modulus, and protein thermodynamic stability. Furthermore, analysis of the network relaxation behavior identifies a crossover from an unfolding dominated regime to an entanglement dominated regime. Control of in situ protein unfolding and entanglement provides an important route to finely tune the architecture, mechanics, and dynamic relaxation of protein hydrogels. Such predictive control will be advantageous for future smart biomaterials for applications which require responsive and dynamic modulation of mechanical properties and biological function.

Folding and Binding Mechanisms of the SH2 Domain from Crkl

SH2 domains are structural modules specialized in the recognition and binding of target sequences containing a phosphorylated tyrosine residue. They are mostly incorporated in the 3D structure of scaffolding proteins that represent fundamental regulators of several signaling pathways. Among those, Crkl plays key roles in cell physiology by mediating signals from a wide range of stimuli, and its overexpression is associated with several types of cancers. In myeloid cells expressing the oncogene BCR/ABL, one interactor of Crkl-SH2 is the focal adhesion protein Paxillin, and this interaction is crucial in leukemic transformation. In this work, we analyze both the folding pathway of Crkl-SH2 and its binding reaction with a peptide mimicking Paxillin, under different ionic strength and pH conditions, by using means of fluorescence spectroscopy. From a folding perspective, we demonstrate the presence of an intermediate along the reaction. Moreover, we underline the importance of the electrostatic interactions in the early event of recognition, occurring between the phosphorylated tyrosine of the Paxillin peptide and the charge residues of Crkl-SH2. Finally, we highlight a pivotal role of a highly conserved histidine residue in the stabilization of the binding complex. The experimental results are discussed in light of previous works on other SH2 domains.

The Denaturant- and Mutation-Induced Disassembly of Pseudomonas aeruginosa Hexameric Hfq Y55W Mutant

Although oligomeric proteins are predominant in cells, their folding is poorly studied at present. This work is focused on the denaturant- and mutation-induced disassembly of the hexameric mutant Y55W of the Qβ host factor (Hfq) from mesophilic Pseudomonas aeruginosa (Pae). Using intrinsic tryptophan fluorescence, dynamic light scattering (DLS), and high-performance liquid chromatography (HPLC), we show that the dissociation of Hfq Y55W occurs either under the effect of GuHCl or during the pre-denaturing transition, when the protein concentration is decreased, with both events proceeding through the accumulation of stable intermediate states. With an extremely low pH of 1.4, a low ionic strength, and decreasing protein concentration, the accumulated trimers and dimers turn into monomers. Also, we report on the structural features of monomeric Hfq resulting from a triple mutation (D9A/V43R/Y55W) within the inter-subunit surface of the protein. This globular and rigidly packed monomer displays a high thermostability and an oligomer-like content of the secondary structure, although its urea resistance is much lower.

αα-hub coregulator structure and flexibility determine transcription factor binding and selection in regulatory interactomes

Formation of transcription factor (TF)-coregulator complexes is a key step in transcriptional regulation, with coregulators having essential functions as hub nodes in molecular networks. How specificity and selectivity are maintained in these nodes remain open questions. In this work, we addressed specificity in transcriptional networks using complexes formed between TFs and αα-hubs, which are defined by a common αα-hairpin secondary structure motif, as a model. Using NMR spectroscopy and binding thermodynamics, we analyzed the structure, dynamics, stability, and ligand-binding properties of the Arabidopsis thaliana RST domains from TAF4 and known binding partner RCD1, and the TAFH domain from human TAF4, allowing comparison across species, functions, and architectural contexts. While these αα-hubs shared the αα-hairpin motif, they differed in length and orientation of accessory helices as well as in their thermodynamic profiles of ligand binding. Whereas biologically relevant RCD1-ligand pairs displayed high affinity driven by enthalpy, TAF4-ligand interactions were entropy driven and exhibited less binding-induced structuring. We in addition identified a thermal unfolding state with a structured core for all three domains, although the temperature sensitivity differed. Thermal stability studies suggested that initial unfolding of the RCD1-RST domain localized around helix 1, lending this region structural malleability, while effects in TAF4-RST were more stochastic, suggesting variability in structural adaptability upon binding. Collectively, our results support a model in which hub structure, flexibility, and binding thermodynamics contribute to αα-hub-TF binding specificity, a finding of general relevance to the understanding of coregulator-ligand interactions and interactome sizes.

Convergent Evolution of Multiple Mutations Improves the Viral Fitness of SARS-CoV-2 Variants by Balancing Positive and Negative Selection

Multiple mutations have been seen to undergo convergent evolution in SARS-CoV-2 variants of concern. One such evolution occurs in Beta, Gamma, and Omicron variants at three amino acid positions K417, E484, and N501 in the receptor binding domain of the spike protein. We examined the physical mechanisms underlying the convergent evolution of three mutations K417T/E484K/N501Y by delineating the individual and collective effects of mutations on binding to angiotensin converting enzyme 2 receptor, immune escape from neutralizing antibodies, protein stability, and expression. Our results show that each mutation serves a distinct function that improves virus fitness supporting its positive selection, even though individual mutations have deleterious effects that make them prone to negative selection. Compared to the wild-type, K417T escapes Class 1 antibodies and has increased stability and expression; however, it has decreased receptor binding. E484K escapes Class 2 antibodies; however, it has decreased receptor binding, stability, and expression. N501Y increases receptor binding; however, it has decreased stability and expression. When these mutations come together, the deleterious effects are mitigated due to the presence of compensatory effects. Triple mutant K417T/E484K/N501Y has increased receptor binding, escapes both Class 1 and Class 2 antibodies, and has similar stability and expression as that of the wild-type. These results show that the convergent evolution of multiple mutations enhances viral fitness on different fronts by balancing both positive and negative selection and improves the chances of selection of mutations together.

Conserved and divergent chaperoning effects of Hsp60/10 chaperonins on protein folding landscapes

The GroEL/ES chaperonin cavity surface charge properties, especially the negative charges, play an important role in its capacity to assist intracavity protein folding. Remarkably, the larger fraction of GroEL/ES negative charges are not conserved among different bacterial species, resulting in a large variation in negative-charge density in the GroEL/ES cavity across prokaryotes. Intriguingly, eukaryotic GroEL/ES homologs have the lowest negative-charge density in the chaperonin cavity. This prompted us to investigate if GroEL’s chaperoning mechanism changed during evolution. Using a model in vivo GroEL/ES substrate, we show that the ability of GroEL/ES to buffer entropic traps in the folding pathway of its substrate was partially dependent upon the negative-charge density inside its cavity. While this activity of GroEL/ES was found to be essential for Escherichia coli, it has been perfected in some organisms and diminished in others. However, irrespective of their charges, all the tested homologs retained their ability to regulate polypeptide chain collapse and remove enthalpic traps from folding pathways. The ability of these GroEL/ES homologs to buffer mutational variations in a model substrate correlated with their negative-charge density. Thus, Hsp60/10 chaperonins in different organisms may have changed to accommodate a different spectrum of mutations on their substrates.

Mechanism of Osmolyte Stabilization-Destabilization of Proteins: Experimental Evidence

In this work, we investigated the influence of stabilizing (N,N,N-trimethylglycine) and destabilizing (urea) osmolytes on the hydration spheres of biomacromolecules in folded forms (trpzip-1 peptide and hen egg white lysozyme—hewl) and unfolded protein models (glycine—GLY and N-methylglycine—NMG) by means of infrared spectroscopy. GLY and NMG were clearly limited as minimal models for unfolded proteins and should be treated with caution. We isolated the spectral share of water changed simultaneously by the biomacromolecule/model molecule and the osmolyte, which allowed us to provide unambiguous experimental arguments for the mechanism of stabilization/destabilization of proteins by osmolytes. In the case of both types of osmolytes, the decisive factor determining the equilibrium folded/unfolded state of protein was the enthalpy effect exerted on the hydration spheres of proteins in both forms. In the case of stabilizing osmolytes, enthalpy was also favored by entropy, as the unfolded state of a protein was more entropically destabilized than the folded state.

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔG_u°'(H₂O), of UBA(1) is 2.4 kcal mol^-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, k_f, in the absence of a denaturant of 13,000 s^-1 and a Tanford β-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone ¹⁵NH, ¹³CO, and ¹³Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

Recipes for Inducing Cold Denaturation in an Otherwise Stable Protein

Although cold denaturation is a fundamental phenomenon common to all proteins, it can only be observed in a handful of cases where it occurs at temperatures above the freezing point of water. Understanding the mechanisms that determine cold denaturation and the rules that permit its observation is an important challenge. A way to approach them is to be able to induce cold denaturation in an otherwise stable protein by means of mutations. Here, we studied CyaY, a relatively stable bacterial protein with no detectable cold denaturation and a high melting temperature of 54 °C. We have characterized for years the yeast orthologue of CyaY, Yfh1, a protein that undergoes cold and heat denaturation at 5 and 35 °C, respectively. We demonstrate that, by transferring to CyaY the lessons learnt from Yfh1, we can induce cold denaturation by introducing a restricted number of carefully designed mutations aimed at destabilizing the overall fold and inducing electrostatic frustration. We used molecular dynamics simulations to rationalize our findings and demonstrate the individual effects observed experimentally with the various mutants. Our results constitute the first example of rationally designed cold denaturation and demonstrate the importance of electrostatic frustration on the mechanism of cold denaturation.

Functionalization of a symmetric protein scaffold: Redundant folding nuclei and alternative oligomeric folding pathways

Successful de novo protein design ideally targets specific folding kinetics, stability thermodynamics, and biochemical functionality, and the simultaneous achievement of all these criteria in a single step design is challenging. Protein design is potentially simplified by separating the problem into two steps: (a) an initial design of a protein "scaffold" having appropriate folding kinetics and stability thermodynamics, followed by (b) appropriate functional mutation-possibly involving insertion of a peptide functional "cassette." This stepwise approach can also separate the orthogonal effects of the "stability/function" and "foldability/function" tradeoffs commonly observed in protein design. If the scaffold is a protein architecture having an exact rotational symmetry, then there is the potential for redundant folding nuclei and multiple equivalent sites of functionalization; thereby enabling broader functional adaptation. We describe such a "scaffold" and functional "cassette" design strategy applied to a β-trefoil threefold symmetric architecture and a heparin ligand functionality. The results support the availability of redundant folding nuclei within this symmetric architecture, and also identify a minimal peptide cassette conferring heparin affinity. The results also identify an energy barrier of destabilization that switches the protein folding pathway from monomeric to trimeric, thereby identifying another potential advantage of symmetric protein architecture in de novo design.

In-Cell Structural Biology by NMR: The Benefits of the Atomic Scale

In-cell structural biology aims at extracting structural information about proteins or nucleic acids in their native, cellular environment. This emerging field holds great promise and is already providing new facts and outlooks of interest at both fundamental and applied levels. NMR spectroscopy has important contributions on this stage: It brings information on a broad variety of nuclei at the atomic scale, which ensures its great versatility and uniqueness. Here, we detail the methods, the fundamental knowledge, and the applications in biomedical engineering related to in-cell structural biology by NMR. We finally propose a brief overview of the main other techniques in the field (EPR, smFRET, cryo-ET, etc.) to draw some advisable developments for in-cell NMR. In the era of large-scale screenings and deep learning, both accurate and qualitative experimental evidence are as essential as ever to understand the interior life of cells. In-cell structural biology by NMR spectroscopy can generate such a knowledge, and it does so at the atomic scale. This review is meant to deliver comprehensive but accessible information, with advanced technical details and reflections on the methods, the nature of the results, and the future of the field.

Molten globule-like transition state of protein barnase measured with calorimetric force spectroscopy

Understanding the molecular forces driving the unfolded polypeptide chain to self-assemble into a functional native structure remains an open question. However, identifying the states visited during protein folding (e.g., the transition state between the unfolded and native states) is tricky due to their transient nature. Here, we introduce calorimetric force spectroscopy in a temperature jump optical trap to determine the enthalpy, entropy, and heat capacity of the transition state of protein barnase. We find that the transition state has the properties of a dry molten globule, that is, high free energy and low configurational entropy, being structurally similar to the native state. This experimental single-molecule study characterizes the thermodynamic properties of the transition state in funneled energy landscapes.

Understanding how proteins fold into their native structure is a fundamental problem in biophysics, crucial for protein design. It has been hypothesized that the formation of a molten globule intermediate precedes folding to the native conformation of globular proteins; however, its thermodynamic properties are poorly known. We perform single-molecule pulling experiments of protein barnase in the range of 7 ^∘C to 37 ^∘C using a temperature-jump optical trap. We derive the folding free energy, entropy and enthalpy, and heat capacity change (ΔCp = 1,050 ± 50 cal/mol·K) at low ionic strength conditions. From the measured unfolding and folding kinetic rates, we also determine the thermodynamic properties of the transition state, finding a significant change in ΔCp (∼90%) between the unfolded and the transition states. In contrast, the major change in enthalpy (∼80%) occurs between the transition and native states. These results highlight a transition state of high energy and low configurational entropy structurally similar to the native state, in agreement with the molten globule hypothesis.

Molecular understanding of calorimetric protein unfolding experiments

Testing and predicting protein stability gained importance because proteins, including antibodies, became pharmacologically relevant in viral and cancer therapies. Isothermal scanning calorimetry is the principle method to study protein stability. Here, we use the excellent experimental heat capacity C_p(T) data from the literature for a critical inspection of protein unfolding as well as for the test of a new cooperative model. In the relevant literature, experimental temperature profiles of enthalpy, H_cal(T), entropy, S_cal(T), and free energy, G_cal(T) are missing. First, we therefore calculate the experimental H_cal(T), S_cal(T), and G_cal(T) from published C_p(T) thermograms. Considering only the unfolding transition proper, the heat capacity and all thermodynamic functions are zero in the region of the native protein. In particular, the free energy of the folded proteins is also zero and G_cal(T) displays a trapezoidal temperature profile when cold denaturation is included. Second, we simulate the DSC-measured thermodynamic properties with a new molecular model based on statistical-mechanical thermodynamics. The model quantifies the protein cooperativity and predicts the aggregate thermodynamic variables of the system with molecular parameters only. The new model provides a perfect simulation of all thermodynamic properties, including the observed trapezoidal G_cal(T) temperature profile. Importantly, the new cooperative model can be applied to a broad range of protein sizes, including antibodies. It predicts not only heat and cold denaturation but also provides estimates of the unfolding kinetics and allows a comparison with molecular dynamics calculations and quasielastic neutron scattering experiments.

Effects of ionic strength on the folding and stability of SAMP1, a ubiquitin-like halophilic protein

Our knowledge of the folding behavior of proteins from extremophiles is limited at this time. These proteins may more closely resemble the primordial proteins selected in early evolution under extreme conditions. The small archaeal modifier protein 1 (SAMP1) studied in this report is an 87-residue protein with a β-grasp fold found in the halophile Haloferax volcanii from the Dead Sea. To gain insight into the effects of salt on the stability and folding mechanism of SAMP1, we conducted equilibrium and kinetic folding experiments as a function of sodium chloride concentration. The results revealed that increasing ionic strength accelerates refolding and slows down unfolding of SAMP1, giving rise to a pronounced salt-induced stabilization. With increasing NaCl concentration, the rate of folding observed via a combination of continuous-flow (0.1-2 ms time range) and stopped-flow measurements (>2 ms) exhibited a >100-fold increase between 0.1 and 1.5 M NaCl and leveled off at higher concentrations. Using the Linderström-Lang smeared charge formalism to model electrostatic interactions in ground and transition states encountered during folding, we showed that the observed salt dependence is dominated by Debye-Hückel screening of electrostatic repulsion among numerous negatively charged residues. Comparisons are also drawn with three well-studied mesophilic members of the β-grasp superfamily: protein G, protein L, and ubiquitin. Interestingly, the folding rate of SAMP1 in 3 M sodium chloride is comparable to that of protein G, ubiquitin, and protein L at lower ionic strength. The results indicate the important role of electrostatic interactions in protein folding and imply that proteins have evolved to minimize unfavorable charge-charge interactions under their specific native conditions.

Metal Ion Binding Induces Local Protein Unfolding and Destabilizes Human Carbonic Anhydrase II

Human carbonic anhydrase II (HCA) is a robust metalloprotein and an excellent biological model system to study the thermodynamics of metal ion coordination. Apo-HCA binds one zinc ion or two copper ions. We studied these binding processes at five temperatures (15-35 °C) using isothermal titration calorimetry, yielding thermodynamic parameters corrected for pH and buffer effects. We then sought to identify binding-induced structural changes. Our data suggest that binding at the active site organizes 6-8 residues; however, copper binding near the N-terminus results in a net unfolding of 6-7 residues. This surprising destabilization was confirmed using circular dichroism and protein stability measurements. Metal binding induced unfolding may represent an important regulatory mechanism, but it may be easily missed by NMR and X-ray crystallography. Thus, in addition to highlighting a highly novel binding-induced unfolding event, we demonstrate the value of calorimetry for studying the structural implications of metal binding.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

Biophysical and Structural Characterization of Ribulose-5-phosphate Epimerase from Leishmania donovani

Pentose phosphate pathway (PPP) plays a crucial role in the maintenance of NADPH/NADP⁺ homeostasis and provides protection against oxidative stress through detoxification of the reactive oxygen species. Ribulose-5-phosphate epimerase (RPE) participates in catalysis of the interconversion of ribulose-5-phosphate (Ru5P) to xylulose-5-phosphate (Xu5P) during PPP, however the structural attributes of this enzyme are still underexplored in many human pathogens including leishmanial parasites. The present study focuses upon cloning, purification and characterization of RPE of Leishmania donovani (LdRPE) using various biophysical and structural approaches. Sequence analysis has shown the presence of trypanosomatid-specific insertions at the N-terminus that are absent in humans and other eukaryotes. Gel filtration chromatography indicated recombinant LdRPE to exist as a dimer in the solution. Circular dichroism studies revealed a higher alpha helical content at physiological pH and temperature that comparatively varies with changing these parameters. Additionally, intrinsic fluorescence and quenching studies of LdRPE have depicted that tryptophan residues are mainly buried in the hydrophobic regions, and the recombinant enzyme is moderately tolerant to urea. Moreover, homology modeling was employed to generate the three-dimensional structure of LdRPE followed by molecular docking with the substrate, product, and substrate analogues. The modeled structure of LdRPE unravelled the presence of conserved active site residues as well as a single binding pocket for the substrate and product, while an in silico study suggested binding of substrate analogues into a similar pocket with more affinity than the substrate. Additionally, molecular dynamics simulation analysis has deciphered complexes of LdRPE with most of the ligands exhibiting more stability than its apo form and lesser fluctuations in active site residues in the presence of ligands. Altogether, our study presents structural insights into leishmanial RPE that could provide the basis for its implication to develop potent antileishmanials.

Denaturant- or ligand-induced change in protein volume by pressure shift assay

A complete thermodynamic description of protein-ligand binding includes parameters related to pressure and temperature. The changes in protein volume and compressibility upon binding a ligand are pressure-related parameters that are...

Engineered Metal-Binding Sites to Probe Protein Folding Transition States: Psi Analysis

The formation of the transition state ensemble (TSE) represents the rate-limiting step in protein folding. The TSE is the least populated state on the pathway, and its characterization remains a challenge. Properties of the TSE can be inferred from the effects on folding and unfolding rates for various perturbations. A difficulty remains on how to translate these kinetic effects to structural properties of the TSE. Several factors can obscure the translation of point mutations in the frequently used method, "mutational Phi analysis." We take a complementary approach in "Psi analysis," employing rationally inserted metal binding sites designed to probe pairwise contacts in the TSE. These contacts can be confidently identified and used to construct structural models of the TSE. The method has been applied to multiple proteins and consistently produces a considerably more structured and native-like TSE than Phi analysis. This difference has significant implications to our understanding of protein folding mechanisms. Here we describe the application of the method and discuss how it can be used to study other conformational transitions such as binding.

Prediction and Validation of a Protein's Free Energy Surface Using Hydrogen Exchange and (Importantly) Its Denaturant Dependence

The denaturant dependence of hydrogen-deuterium exchange (HDX) is a powerful measurement to identify the breaking of individual H-bonds and map the free energy surface (FES) of a protein including the very rare states. Molecular dynamics (MD) can identify each partial unfolding event with atomic-level resolution. Hence, their combination provides a great opportunity to test the accuracy of simulations and to verify the interpretation of HDX data. For this comparison, we use Upside, our new and extremely fast MD package that is capable of folding proteins with an accuracy comparable to that of all-atom methods. The FESs of two naturally occurring and two designed proteins are so generated and compared to our NMR/HDX data. We find that Upside's accuracy is considerably improved upon modifying the energy function using a new machine-learning procedure that trains for proper protein behavior including realistic denatured states in addition to stable native states. The resulting increase in cooperativity is critical for replicating the HDX data and protein stability, indicating that we have properly encoded the underlying physiochemical interactions into an MD package. We did observe some mismatch, however, underscoring the ongoing challenges faced by simulations in calculating accurate FESs. Nevertheless, our ensembles can identify the properties of the fluctuations that lead to HDX, whether they be small-, medium-, or large-scale openings, and can speak to the breadth of the native ensemble that has been a matter of debate.

Experimental Characterization of the Interaction between the N-Terminal SH3 Domain of Crkl and C3G

Crkl is a protein involved in the onset of several cancer pathologies that exerts its function only through its protein-protein interaction domains, a SH2 domain and two SH3 domains. SH3 domains are small protein interaction modules that mediate the binding and recognition of proline-rich sequences. One of the main physiological interactors of Crkl is C3G (also known as RAPGEF1), an interaction with key implications in regulating cellular growth and differentiation, cell morphogenesis and adhesion processes. Thus, understanding the interaction between Crkl and C3G is fundamental to gaining information about the molecular determinants of the several cancer pathologies in which these proteins are involved. In this paper, through a combination of fast kinetics at different experimental conditions and site-directed mutagenesis, we characterize the binding reaction between the N-SH3 domain of Crkl and a peptide mimicking a specific portion of C3G. Our results show a clear effect of pH on the stability of the complex, due to the protonation of negatively charged residues in the binding pocket of N-SH3. Our results are discussed under the light of previous work on SH3 domains.

Protein nanofibrils and their use as building blocks of sustainable materials

The development towards a sustainable society requires a radical change of many of the materials we currently use. Besides the replacement of plastics, derived from petrochemical sources, with renewable alternatives, we will also need functional materials for applications in areas ranging from green energy and environmental remediation to smart foods. Proteins could, with their intriguing ability of self-assembly into various forms, play important roles in all these fields. To achieve that, the code for how to assemble hierarchically ordered structures similar to the protein materials found in nature must be cracked. During the last decade it has been demonstrated that amyloid-like protein nanofibrils (PNFs) could be a steppingstone for this task. PNFs are formed by self-assembly in water from a range of proteins, including plant resources and industrial side streams. The nanofibrils display distinct functional features and can be further assembled into larger structures. PNFs thus provide a framework for creating ordered, functional structures from the atomic level up to the macroscale. This review address how industrial scale protein resources could be transformed into PNFs and further assembled into materials with specific mechanical and functional properties. We describe what is required from a protein to form PNFs and how the structural properties at different length scales determine the material properties. We also discuss potential chemical routes to modify the properties of the fibrils and to assemble them into macroscopic structures.

A newly introduced salt bridge cluster improves structural and biophysical properties of de novo TIM barrels

Protein stability can be fine‐tuned by modifying different structural features such as hydrogen‐bond networks, salt bridges, hydrophobic cores, or disulfide bridges. Among these, stabilization by salt bridges is a major challenge in protein design and engineering since their stabilizing effects show a high dependence on the structural environment in the protein, and therefore are difficult to predict and model. In this work, we explore the effects on structure and stability of an introduced salt bridge cluster in the context of three different de novo TIM barrels. The salt bridge variants exhibit similar thermostability in comparison with their parental designs but important differences in the conformational stability at 25°C can be observed such as a highly stabilizing effect for two of the proteins but a destabilizing effect to the third. Analysis of the formed geometries of the salt bridge cluster in the crystal structures show either highly ordered salt bridge clusters or only single salt bridges. Rosetta modeling of the salt bridge clusters results in a good prediction of the tendency on stability changes but not the geometries observed in the three‐dimensional structures. The results show that despite the similarities in protein fold, the salt bridge clusters differently influence the structural and stability properties of the de novo TIM barrel variants depending on the structural background where they are introduced.

PDB Code(s): 7OSU, 7OT7, 7OSV, 7OT8 and 7P12;

The Nuts and Bolts of SARS-CoV-2 Spike Receptor-Binding Domain Heterologous Expression

COVID-19 is a highly infectious disease caused by a newly emerged coronavirus (SARS-CoV-2) that has rapidly progressed into a pandemic. This unprecedent emergency has stressed the significance of developing effective therapeutics to fight the current and future outbreaks. The receptor-binding domain (RBD) of the SARS-CoV-2 surface Spike protein is the main target for vaccines and represents a helpful "tool" to produce neutralizing antibodies or diagnostic kits. In this work, we provide a detailed characterization of the native RBD produced in three major model systems: Escherichia coli, insect and HEK-293 cells. Circular dichroism, gel filtration chromatography and thermal denaturation experiments indicated that recombinant SARS-CoV-2 RBD proteins are stable and correctly folded. In addition, their functionality and receptor-binding ability were further evaluated through ELISA, flow cytometry assays and bio-layer interferometry.

SARS-CoV-2 beta variant substitutions alter spike glycoprotein receptor binding domain structure and stability

The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) and the subsequent COVID-19 pandemic have visited a terrible cost on the world in the forms of disease, death, and economic turmoil. The rapid development and deployment of extremely effective vaccines against SARS-CoV-2 have seemingly brought within reach the end of the pandemic. However, the virus has acquired mutations. and emerging variants of concern are more infectious and reduce the efficacy of existing vaccines. Although promising efforts to combat these variants are underway, the evolutionary pressures leading to these variants are poorly understood. To that end, here we have studied the effects on the structure and function of the SARS-CoV-2 spike glycoprotein receptor-binding domain of three amino-acid substitutions found in several variants of concern, including alpha (B.1.1.7), beta (B.1.351), and gamma (P.1). We found that these substitutions alter the receptor-binding domain structure, stability, and ability to bind to angiotensin converting enzyme 2, in such a way as to possibly have opposing and compensatory effects. These findings provide new insights into how these variants of concern may have been selected for infectivity while maintaining the structure and stability of the receptor binding domain.

Hofmeister effects on protein stability are dependent on the nature of the unfolded state

The interpretation of a salt's effect on protein stability traditionally discriminates low concentration regimes (<0.3 M), dominated by electrostatic forces, and high concentration regimes, generally described by ion-specific Hofmeister effects. However, increased theoretical and experimental studies have highlighted observations of the Hofmeister phenomena at concentration ranges as low as 0.001 M. Reasonable quantitative predictions of such observations have been successfully achieved throughout the inclusion of ion dispersion forces in classical electrostatic theories. This molecular description is also on the basis of quantitative estimates obtained resorting to surface/bulk solvent partition models developed for ion-specific Hofmeister effects. However, the latter are limited by the availability of reliable structures representative of the unfolded state. Here, we use myoglobin as a model to explore how ion-dependency on the nature of the unfolded state affects protein stability, combining spectroscopic techniques with molecular dynamic simulations. To this end, the thermal and chemical stability of myoglobin was assessed in the presence of three different salts (NaCl, (NH₄)₂SO₄ and Na₂SO₄), at physiologically relevant concentrations (0-0.3 M). We observed mild destabilization of the native state induced by each ion, attributed to unfavorable neutralization and hydrogen-bonding with the protein side-chains. Both effects, combined with binding of Na⁺, Cl^- and SO₄^2- to the thermally unfolded state, resulted in an overall destabilization of the protein. Contrastingly, ion binding was hindered in the chemically unfolded conformation, due to occupation of the binding sites by urea molecules. Such mechanistic action led to a lower degree of destabilization, promoting surface tension effects that stabilized myoglobin according to the Hofmeister series. Therefore, we demonstrate that Hofmeister effects on protein stability are modulated by the heterogeneous physico-chemical nature of the unfolded state. Altogether, our findings evidence the need to characterize the structure of the unfolded state when attempting to dissect the molecular mechanisms underlying the effects of salts on protein stability.

The A39G FF domain folds on a volcano-shaped free energy surface via separate pathways

Conformational dynamics play critical roles in protein folding, misfolding, function, misfunction, and aggregation. While detecting and studying the different conformational states populated by protein molecules on their free energy surfaces (FESs) remain a challenge, NMR spectroscopy has emerged as an invaluable experimental tool to explore the FES of a protein, as conformational dynamics can be probed at atomic resolution over a wide range of timescales. Here, we use chemical exchange saturation transfer (CEST) to detect "invisible" minor states on the energy landscape of the A39G mutant FF domain that exhibited "two-state" folding kinetics in traditional experiments. Although CEST has mostly been limited to studies of processes with rates between ∼5 to 300 s^-1 involving sparse states with populations as low as ∼1%, we show that the line broadening that is often associated with minor state dips in CEST profiles can be exploited to inform on additional conformers, with lifetimes an order of magnitude shorter and populations close to 10-fold smaller than what typically is characterized. Our analysis of CEST profiles that exploits the minor state linewidths of the 71-residue A39G FF domain establishes a folding mechanism that can be described in terms of a four-state exchange process between interconverting states spanning over two orders of magnitude in timescale from ∼100 to ∼15,000 μs. A similar folding scheme is established for the wild-type domain as well. The study shows that the folding of this small domain proceeds through a pair of sparse, partially structured intermediates via two discrete pathways on a volcano-shaped FES.

Folding Steps in the Fibrillation of Functional Amyloid: Denaturant Sensitivity Reveals Common Features in Nucleation and Elongation

Functional bacterial amyloids (FuBA) are intrinsically disordered proteins (IDPs) which rapidly and efficiently aggregate, forming extremely stable fibrils. The conversion from IDP to amyloid is evolutionarily optimized and likely couples folding to association. Many FuBA contain several imperfect repeat sequences which contribute to the stability of mature FuBA fibrils. Aggregation can be considered an intermolecular extension of the process of intramolecular protein folding which has traditionally been studied using chemical denaturants. Here we employ denaturants to investigate folding steps during fibrillation of CsgA and FapC. We quantify protein compactification (i.e. the extent of burial of otherwise exposed surface area upon association of proteins) during different stages of fibrillation based on the dependence of fibrillation rate constants on the denaturant concentration (m-values) determined from fibrillation curves. For both proteins, urea mainly affects nucleation and elongation (not fragmentation), consistent with the fact that these steps involve both intra- and intermolecular association. The two steps have similar m-values, indicating that activation steps in nucleation and elongation involve the same level of folding. Surprisingly, deletion of two or three repeats from FapC leads to larger m-values (i.e. higher compactification) during the activation step of fibril growth. This observation is extended by SAXS analysis of the fibrils which indicates that weakening of the amyloidogenic core caused by repeat deletions causes a larger portion of normally unstructured regions of the protein to be included into the amyloid backbone. We conclude that the sensitivity of fibrillation to denaturants can provide useful insight into molecular mechanisms of aggregation.

Surface residues and non-additive interactions stabilize a consensus homeodomain protein

Despite the widely reported success of consensus design in producing highly stabilized proteins, little is known about the physical mechanisms underlying this stabilization. Here we explore the potential sources of stabilization by performing a systematic analysis of the 29 substitutions that we previously found to collectively stabilize a consensus homeodomain compared to an extant homeodomain. By separately introducing groups of consensus substitutions that alter or preserve charge state, occur at varying degrees of residue burial, and occur at positions of varying degrees of conservation, we determine the extent to which these three features contribute to the consensus stability enhancement. Surprisingly, we find that the largest total contribution to stability comes from consensus substitutions on the protein surface and that the largest per-substitution contributions come from substitutions that maintain charge state. This finding suggests that although consensus proteins are often enriched in charged residues, consensus stabilization does not result primarily from interactions involving charged residues. Although consensus substitutions at strongly conserved positions also contribute disproportionately to stabilization, significant stabilization is also contributed from substitutions at weakly conserved positions. Furthermore, we find that identical consensus substitutions show larger stabilizing effects when introduced into the consensus background than when introduced into an extant homeodomain, indicating that synergistic, stabilizing interactions among the consensus residues contribute to consensus stability enhancement of the homeodomain. By measuring DNA binding affinity for the same set of variants, we find that although consensus design of the homeodomain increases both affinity and folding stability, it does so using a largely non-overlapping set of substitutions.

Single-Molecule Force Spectroscopy Studies of Missense Titin Mutations That Are Likely Causing Cardiomyopathy

The giant muscle protein titin plays important roles in heart function. Mutations in titin have emerged as a major cause of familial cardiomyopathy. Missense mutations have been identified in cardiomyopathy patients; however, it is challenging to distinguish disease-causing mutations from benign ones. Given the importance of titin mechanics in heart function, it is critically important to elucidate the mechano-phenotypes of cardiomyopathy-causing mutations found in the elastic I-band part of cardiac titin. Using single-molecule atomic force microscopy (AFM) and equilibrium chemical denaturation, we investigated the mechanical and thermodynamic effects of two missense mutations, R57C-I94 and S22P-I84, found in the elastic I-band part of cardiac titin that were predicted to be likely causing cardiomyopathy by bioinformatics analysis. Our AFM results showed that mutation R57C had a significant destabilization effect on the I94 module. R57C reduced the mechanical unfolding force of I94 by ∼30-40 pN, accelerated the unfolding kinetics, and decelerated the folding. These effects collectively increased the unfolding propensity of I94, likely resulting in altered titin elasticity. In comparison, S22P led to only modest destabilization of I84, with a decrease in unfolding force by ∼10 pN. It is unlikely that such a modest destabilization would lead to a change in titin elasticity. These results will serve as the first step toward elucidating mechano-phenotypes of cardiomyopathy-causing mutations in the elastic I-band.

Functional-State Dependence of Picosecond Protein Dynamics

We examine temperature-dependent picosecond dynamics of two benchmarking proteins lysozyme and cytochrome c using temperature-dependent terahertz permittivity measurements. We find that a double Arrhenius temperature dependence with activation energies E₁ ∼ 0.1 kJ/mol and E₂ ∼ 10 kJ/mol fits the folded and ligand-free state response. The higher activation energy is consistent with the so-called protein dynamical transition associated with beta relaxations at the solvent-protein interface. The lower activation energy is consistent with correlated structural motions. When the structure is removed by denaturing, the lower-activation-energy process is no longer present. Additionally, the lower-activation-energy process is diminished with ligand binding but not for changes in the internal oxidation state. We suggest that the lower-energy activation process is associated with collective structural motions that are no longer accessible with denaturing or binding.

Effect of urea on the hydration and aggregation of hydrophobic and amphiphilic solute models: Implications to protein aggregation

Despite the emergence of a molecular picture of urea's protein unfolding mechanism in the past few decades, less is known about its action mechanism on protein aggregation. This is especially relevant for understanding the aggregation of amyloid proteins and peptides, implicated in several neurodegenerative diseases. While urea is believed to weaken the hydrophobic effect, a picture consistent with the decrease in the excess chemical potential of sufficiently large alkanes, interactions with protein polar side chains and backbone atoms are also important. Here, we study, through molecular dynamics, the hydration and aggregation of several alkanes and amphiphilic "mutants" of n-dodecane, in an 8M aqueous urea solution, aiming at getting insight into urea's mode of action. A size-dependent crossover temperature is found, above which the hydration of the alkanes is favored in the aqueous urea solution. The hydration of the alkanes is enhanced via entropy, with the enthalpy opposing hydration, consistent with experiments. The reason is that although solute-solvent interactions are favorable, these are overwhelmed by urea-water and urea-urea interactions. In contrast, water-water interactions and entropy are favored by a water depletion around the solute and a reduced water depletion around methane explains its exceptional solubility decrease. Furthermore, we show that while urea favors the hydration of n-dodecane and the amphiphilic mutants, it slightly enhances and reduces, respectively, the aggregation of the alkanes and the amphiphilic mutants. Thus, opposite to the common view, our results show that urea does not necessarily weaken hydrophobic interactions despite solvation being favored.

Determining folding and binding properties of the C-terminal SH2 domain of SHP2

SH2 domains are a class of protein–protein interaction modules with the function to recognize and bind sequences characterized by the presence of a phosphorylated tyrosine. SHP2 is a protein phosphatase involved in the Ras‐ERK1/2 signaling pathway that possess two SH2 domains, namely, N‐SH2 and C‐SH2, that mediate the interaction of SHP2 with various partners and determine the regulation of its catalytic activity. One of the main interactors of the SH2 domains of SHP2 is Gab2, a scaffolding protein with critical role in determining cell differentiation. Despite their key biological role and the importance of a correct native fold to ensure it, the mechanism of binding of SH2 domains with their ligands and the determinants of their stability have been poorly characterized. In this article, we present a comprehensive kinetic study of the folding of the C‐SH2 domain and the binding mechanism with a peptide mimicking a region of Gab2. Our data, obtained at different pH and ionic strength conditions and supported by site‐directed mutagenesis, highlight the role of electrostatic interactions in the early events of recognition. Interestingly, our results suggest a key role of a highly conserved histidine residue among SH2 family in the interaction with negative charges carried by the phosphotyrosine of Gab2. Moreover, the analysis of the equilibrium and kinetic folding data of C‐SH2 describes a complex mechanism implying a change in rate‐limiting step at high denaturant concentrations. Our data are discussed under the light of previous works on N‐SH2 domain of SHP2 and other SH2 domains.

Frustrated peptide chains at the fibril tip control the kinetics of growth of amyloid-β fibrils

Amyloid fibrillization is an exceedingly complex process in which incoming peptide chains bind to the fibril while concertedly folding. The coupling between folding and binding is not fully understood. We explore the molecular pathways of association of Aβ40 monomers to fibril tips by combining time-resolved in situ scanning probe microscopy with molecular modeling. The comparison between experimental and simulation results shows that a complex supported by nonnative contacts is present in the equilibrium structure of the fibril tip and impedes fibril growth in a supersaturated solution. The unraveling of this frustrated state determines the rate of fibril growth. The kinetics of growth of freshly cut fibrils, in which the bulk fibril structure persists at the tip, complemented by molecular simulations, indicate that this frustrated complex comprises three or four monomers in nonnative conformations and likely is contained on the top of a single stack of peptide chains in the fibril structure. This pathway of fibril growth strongly deviates from the common view that the conformational transformation of each captured peptide chain is templated by the previously arrived peptide. The insights into the ensemble structure of the frustrated complex may guide the search for suppressors of Aβ fibrillization. The uncovered dynamics of coupled structuring and assembly during fibril growth are more complex than during the folding of most globular proteins, as they involve the collective motions of several peptide chains that are not guided by a funneled energy landscape.

Quantification of Conformational Entropy Unravels Effect of Disordered Flanking Region in Coupled Folding and Binding

Intrinsic disorder (ID) constitutes a new dimension to the protein structure-function relationship. The ability to undergo conformational changes upon binding is a key property of intrinsically disordered proteins and remains challenging to study using conventional methods. A 1994 paper by R. S. Spolar and M. T. Record presented a thermodynamic approach for estimating changes in conformational entropy based on heat capacity changes, allowing quantification of residues folding upon binding. Here, we adapt the method for studies of intrinsically disordered proteins. We integrate additional data to provide a broader experimental foundation for the underlying relations and, based on >500 protein-protein complexes involving disordered proteins, reassess a key relation between polar and nonpolar surface area changes, previously determined using globular protein folding. We demonstrate the improved suitability of the adapted method to studies of the folded αα-hub domain RST from radical-induced cell death 1, whose interactome is characterized by ID. From extensive thermodynamic data, quantifying the conformational entropy changes upon binding, and comparison to the NMR structure, the adapted method improves accuracy for ID-based studies. Furthermore, we apply the method, in conjunction with NMR, to reveal hitherto undetected effects of interaction-motif context. Thus, inclusion of the disordered context of the DREB2A RST-binding motif induces structuring of the binding motif, resulting in major enthalpy-entropy compensation in the interaction interface. This study, also evaluating additional interactions, demonstrates the strength of the ID-adapted Spolar-Record thermodynamic approach for dissection of structural features of ID-based interactions, easily overlooked in traditional studies, and for translation of these into mechanistic knowledge.