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

Human apo-metallothionein 1a is not a random coil: Evidence from guanidinium chloride, high temperature, and acidic pH unfolding studies

The structures of apo-metallothioneins (apo-MTs) have been relatively elusive due to their fluxional, disordered state which has been difficult to characterize. However, intrinsically disordered protein (IDP) structures are rather diverse, which raises questions about where the structure of apo-MTs fit into the protein structural spectrum. In this paper, the unfolding transitions of apo-MT1a are discussed with respect to the effect of the chemical denaturant GdmCl, temperature conditions, and pH environment. Cysteine modification in combination with electrospray ionization mass spectrometry was used to probe the unfolding transition of apo-MT1a in terms of cysteine exposure. Circular dichroism spectroscopy was also used to monitor the change in secondary structure as a function of GdmCl concentration. For both of these techniques, cooperative unfolding was observed, suggesting that apo-MT1a is not a random coil. More GdmCl was required to unfold the protein backbone than to expose the cysteines, indicating that cysteine exposure is likely an early step in the unfolding of apo-MT1a. MD simulations complement the experimental results, suggesting that apo-MT1a adopts a more compact structure than expected for a random coil. Overall, these results provide further insight into the intrinsically disordered structure of apo-MT.

The Streptococcus phage protein paratox is an intrinsically disordered protein

The bacteriophage protein paratox (Prx) blocks quorum sensing in its streptococcal host by directly binding the signal receptor and transcription factor ComR. This reduces the ability of Streptococcus to uptake environmental DNA and protects phage DNA from damage by recombination. Past work characterizing the Prx:ComR molecular interaction revealed that paratox adopts a well-ordered globular fold when bound to ComR. However, solution-state biophysical measurements suggested that Prx may be conformationally dynamic. To address this discrepancy, we investigated the stability and dynamic properties of Prx in solution using circular dichroism, nuclear magnetic resonance, and several fluorescence-based protein folding assays. Our work shows that under dilute buffer conditions Prx is intrinsically disordered. We also show that the addition of kosmotropic salts or protein stabilizing osmolytes induces Prx folding. However, the solute stabilized fold is different from the conformation Prx adopts when it is bound to ComR. Furthermore, we have characterized Prx folding thermodynamics and folding kinetics through steady-state fluorescence and stopped flow kinetic measurements. Our results show that Prx is a highly dynamic protein in dilute solution, folding and refolding within the 10 ms timescale. Overall, our results demonstrate that the streptococcal phage protein Prx is an intrinsically disordered protein in a two-state equilibrium with a solute-stabilized folded form. Furthermore, the solute-stabilized fold is likely the predominant form of Prx in a solute-crowded bacterial cell. Finally, our work suggests that Prx binds and inhibits ComR, and thus quorum sensing in Streptococcus, by a combination of conformational selection and induced-fit binding mechanisms.

Protein Binding Leads to Reduced Stability and Solvated Disorder in the Polystyrene Nanoparticle Corona

Understanding the conformation of proteins in the nanoparticle corona has important implications in how organisms respond to nanoparticle-based drugs. These proteins coat the nanoparticle surface, and their properties will influence the nanoparticle's interaction with cell targets and the immune system. While some coronas are thought to be disordered, two key unanswered questions are the degree of disorder and solvent accessibility. Here, a model is developed for protein corona disorder in polystyrene nanoparticles of varying size. For two different proteins, it is found that binding affinity decreases as nanoparticle size increases. The stoichiometry of binding, along with changes in the hydrodynamic size, supports a highly solvated, disordered protein corona anchored at a small number of attachment sites. The scaling of the stoichiometry versus nanoparticle size is consistent with disordered polymer dimensions. Moreover, it is found that proteins are destabilized less in the presence of larger nanoparticles, and hydrophobic exposure decreases at lower curvatures. The observations hold for proteins on flat polystyrene surfaces, which have the lowest hydrophobic exposure. The model provides an explanation for previous observations of increased amyloid fibrillation rates in the presence of larger nanoparticles, and it may rationalize how cell receptors can recognize protein disorder in therapeutic nanoparticles.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape-all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs-questioning the differences between them and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the consensus protein derived from our full Ribonuclease H sequence alignment is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted set of sequences is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order sequence correlations using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Backbone Modification in a Protein Hydrophobic Core

Targeted protein backbone modification can recreate tertiary structures reminiscent of folds found in nature on artificial scaffolds with improved biostability. Incorporation of altered monomers in such entities is typically limited to sites distant from the hydrophobic core to avoid potential disruptions to folding. This is limiting, as it is advantageous in some applications to incorporate artificial connectivity at buried sites. Here, we report an examination of protein backbone modification targeted specifically to hydrophobic core positions and its impacts on tertiary folded structure and fold stability. Different artificial monomer types are placed at core, core-flanking, or solvent-exposed positions in a compact three-helix protein. Effects on structure and folding energetics are assessed by NMR spectroscopy and biophysical methods. Results show that artificial residues can be well accommodated in the hydrophobic core of a defined tertiary fold, with effects on stability only modestly larger than identical changes at solvent-exposed sites. Collectively, these results provide new insights into folding behavior of protein-like artificial chains as well as strategies for the design of such molecules.

Protein-adaptive differential scanning fluorimetry using conformationally responsive dyes

Differential scanning fluorimetry (DSF) is a technique that reports protein thermal stability via the selective recognition of unfolded states by fluorogenic dyes. However, DSF applications remain limited by protein incompatibilities with existing DSF dyes. Here we overcome this obstacle with the development of a protein-adaptive DSF platform (paDSF) that combines a dye library 'Aurora' with a streamlined procedure to identify protein-dye pairs on demand. paDSF was successfully applied to 94% (66 of 70) of proteins, tripling the previous compatibility and delivering assays for 66 functionally and biochemically diverse proteins, including 10 from severe acute respiratory syndrome coronavirus 2. We find that paDSF can be used to monitor biological processes that were previously inaccessible, demonstrated for the interdomain allostery of O-GlcNAc transferase. The chemical diversity and varied selectivities of Aurora dyes suggest that paDSF functionality may be readily extended. paDSF is a generalizable tool to interrogate protein stability, dynamics and ligand binding.

How hydrophobicity, side chains, and salt affect the dimensions of disordered proteins

Despite the generally accepted role of the hydrophobic effect as the driving force for folding, many intrinsically disordered proteins (IDPs), including those with hydrophobic content typical of foldable proteins, behave nearly as self-avoiding random walks (SARWs) under physiological conditions. Here, we tested how temperature and ionic conditions influence the dimensions of the N-terminal domain of pertactin (PNt), an IDP with an amino acid composition typical of folded proteins. While PNt contracts somewhat with temperature, it nevertheless remains expanded over 10-58°C, with a Flory exponent, ν, >0.50. Both low and high ionic strength also produce contraction in PNt, but this contraction is mitigated by reducing charge segregation. With 46% glycine and low hydrophobicity, the reduced form of snow flea anti-freeze protein (red-sfAFP) is unaffected by temperature and ionic strength and persists as a near-SARW, ν ~ 0.54, arguing that the thermal contraction of PNt is due to stronger interactions between hydrophobic side chains. Additionally, red-sfAFP is a proxy for the polypeptide backbone, which has been thought to collapse in water. Increasing the glycine segregation in red-sfAFP had minimal effect on ν. Water remained a good solvent even with 21 consecutive glycine residues (ν > 0.5), and red-sfAFP variants lacked stable backbone hydrogen bonds according to hydrogen exchange. Similarly, changing glycine segregation has little impact on ν in other glycine-rich proteins. These findings underscore the generality that many disordered states can be expanded and unstructured, and that the hydrophobic effect alone is insufficient to drive significant chain collapse for typical protein sequences.

Design and Thermodynamics Principles to Program the Cooperativity of Molecular Assemblies

Most functional nanosystems in living organisms are constructed using multimeric assemblies that provide multiple advantages over their monomeric counterparts such as cooperative or anti-cooperative responses, integration of multiple signals and self-regulation. Inspired by these natural nanosystems, chemists have been synthesizing self-assembled supramolecular systems over the last 50 years with increasing complexity with applications ranging from biosensing, drug delivery, synthetic biology, and system chemistry. Although many advances have been made concerning the design principles of novel molecular architectures and chemistries, little is still known, however, about how to program their dynamic of assembly so that they can assemble at the required concentration and with the right sensitivity. Here, we used synthetic DNA assemblies and double-mutant cycle analysis to explore the thermodynamic basis to program the cooperativity of molecular assemblies. The results presented here exemplify how programmable molecular assemblies can be efficiently built by fusing interacting domains and optimizing their compaction. They may also provide the rational basis for understanding the thermodynamic and mechanistic principles driving the evolution of multimeric biological complexes.

Stepwise introduction of stabilizing mutations reveals nonlinear additive effects in de novo TIM barrels

Over the past decades, the TIM-barrel fold has served as a model system for the exploration of how changes in protein sequences affect their structural, stability, and functional characteristics, and moreover, how this information can be leveraged to design proteins from the ground up. After numerous attempts to design de novo proteins with this specific fold, sTIM11 was the first validated de novo design of an idealized four-fold symmetric TIM barrel. Subsequent efforts to enhance the stability of this initial design resulted in the development of DeNovoTIMs, a family of de novo TIM barrels with various stabilizing mutations. In this study, we present an investigation into the biophysical and thermodynamic effects upon introducing a varying number of stabilizing mutations per quarter along the sequence of a four-fold symmetric TIM barrel. We compared the base design DeNovoTIM0 without any stabilizing mutations with variants containing mutations in one, two, three, and all four quarters-designated TIM1q, TIM2q, TIM3q, and DeNovoTIM6, respectively. This analysis revealed a stepwise and nonlinear change in the thermodynamic properties that correlated with the number of mutated quarters, suggesting positive nonadditive effects. To shed light on the significance of the location of stabilized quarters, we engineered two variants of TIM2q which contain the same number of mutations but positioned in different quarter locations. Characterization of these TIM2q variants revealed that the mutations exhibit varying effects on the overall protein stability, contingent upon the specific region in which they are introduced. These findings emphasize that the amount and location of stabilized interfaces among the four quarters play a crucial role in shaping the conformational stability of these four-fold symmetric TIM barrels. Analysis of de novo proteins, as described in this study, enhances our understanding of how sequence variations can finely modulate stability in both naturally occurring and computationally designed proteins.

Design of complicated all-α protein structures

A wide range of de novo protein structure designs have been achieved, but the complexity of naturally occurring protein structures is still far beyond these designs. Here, to expand the diversity and complexity of de novo designed protein structures, we sought to develop a method for designing 'difficult-to-describe' α-helical protein structures composed of irregularly aligned α-helices like globins. Backbone structure libraries consisting of a myriad of α-helical structures with five or six helices were generated by combining 18 helix-loop-helix motifs and canonical α-helices, and five distinct topologies were selected for de novo design. The designs were found to be monomeric with high thermal stability in solution and fold into the target topologies with atomic accuracy. This study demonstrated that complicated α-helical proteins are created using typical building blocks. The method we developed will enable us to explore the universe of protein structures for designing novel functional proteins.

Differential effects of disulfide bond formation in TEM-1 versus CTX-M-9 β-lactamase

To investigate how disulfide bonds can impact protein energy landscapes, we surveyed the effects of adding or removing a disulfide in two β-lactamase enzymes, TEM-1 and CTX-M-9. The homologs share a structure and 38% sequence identity, but only TEM-1 contains a native disulfide bond. They also differ in thermodynamic stability and in the number of states populated at equilibrium: CTX-M-9 is two-state whereas TEM-1 has an additional intermediate state. We hypothesized that the disulfide bond is the major underlying determinant for these observed differences in their energy landscapes. To test this, we removed the disulfide bridge from TEM-1 and introduced a disulfide bridge at the same location in CTX-M-9. This modest change to sequence modulates the stabilities-and therefore populations-of TEM-1's equilibrium states and, more surprisingly, creates a novel third state in CTX-M-9. Unlike TEM-1's partially folded intermediate, this third state is a higher-order oligomer with reduced cysteines that retains the native fold and is fully active. Sub-denaturing concentrations of urea shifts the equilibrium to the monomeric form, allowing the disulfide bond to form. Interestingly, comparing the stability of the oxidized monomer with a variant lacking cysteines reveals the disulfide is neither stabilizing nor destabilizing in CTX-M-9, in contrast with the observed stabilization in TEM-1. Thus, we can conclude that engineering disulfide bonds is not always an effective stabilization strategy even when analogous disulfides exist in more stable structural homologs. This study also illustrates how homo-oligomerization can result from a small number of mutations, suggesting complex formation might be easily accessed during a protein family's evolution.

Inhibition and structural insights of leishmanial glutamyl-tRNA synthetase for designing potent therapeutics

Aminoacyl-tRNA synthetases (aaRSs), essential components of the protein synthesizing machinery, have been often chosen for devising therapeutics against parasitic diseases. Due to their relevance in drug development, the current study was designed to explore functional and structural aspects of Leishmania donovani glutamyl-tRNA synthetase (LdGluRS). Hence, LdGluRS was cloned into an expression vector and purified to homogeneity using chromatographic techniques. Purified protein showed maximum enzymatic activity at physiological pH, with more binding capacity towards its cofactor (Adenosine triphosphate, 0.06 ± 0.01 mM) than the cognate substrate (L-glutamate, 9.5 ± 0.5 mM). Remarkably, salicylate inhibited LdGluRS competitively with respect to L-glutamate and exhibited druglikeness with negligible effect on human macrophages. The protein possessed more α-helices (43 %) than β-sheets (12 %), whereas reductions in thermal stability and cofactor-binding affinity, along with variation in mode of inhibition after mutation signified the role of histidine (H60) as a catalytic residue. LdGluRS could also generate a pro-inflammatory milieu in human macrophages by upregulating cytokines. The docking study demonstrated the placement of salicylate into LdGluRS substrate-binding site, and the complex was found to be stable during molecular dynamics (MD) simulation. Altogether, our study highlights the understanding of molecular inhibition and structural features of glutamyl-tRNA synthetase from kinetoplastid parasites.

Salt Effects on Caffeine across Concentration Regimes

Salts affect the solvation thermodynamics of molecules of all sizes; the Hofmeister series is a prime example in which different ions lead to salting-in or salting-out of aqueous proteins. Early work of Tanford led to the discovery that the solvation of molecular surface motifs is proportional to the solvent accessible surface area (SASA), and later studies have shown that the proportionality constant varies with the salt concentration and type. Using multiscale computer simulations combined with vapor-pressure osmometry on caffeine-salt solutions, we reveal that this SASA description captures a rich set of molecular driving forces in tertiary solutions at changing solute and osmolyte concentrations. Central to the theoretical work is a new potential energy function that depends on the instantaneous surface area, salt type, and concentration. Used in, e.g., Monte Carlo simulations, this allows for a highly efficient exploration of many-body interactions and the resulting thermodynamics at elevated solute and salt concentrations.

From haemadin to haemanorm: Synthesis and characterization of full-length haemadin from the leech Haemadipsa sylvestris and of a novel bivalent, highly potent thrombin inhibitor (haemanorm)

Hirudin from Hirudo medicinalis is a bivalent α-Thrombin (αT) inhibitor, targeting the enzyme active site and exosite-I, and is currently used in anticoagulant therapy along with its simplified analogue hirulog. Haemadin, a small protein (57 amino acids) isolated from the land-living leech Haemadipsa sylvestris, selectively inhibits αT with a potency identical to that of recombinant hirudin (KI  = 0.2 pM), with which it shares a common disulfide topology and overall fold. At variance with hirudin, haemadin targets exosite-II and therefore (besides the free protease) it also blocks thrombomodulin-bound αT without inhibiting the active intermediate meizothrombin, thus offering potential advantages over hirudin. Here, we produced in reasonably high yields and pharmaceutical purity (>98%) wild-type haemadin and the oxidation resistant Met5 → nor-Leucine analogue, both inhibiting αT with a KI of 0.2 pM. Thereafter, we used site-directed mutagenesis, spectroscopic, ligand-displacement, and Hydrogen/Deuterium Exchange-Mass Spectrometry techniques to map the αT regions relevant for the interaction with full-length haemadin and with the synthetic N- and C-terminal peptides Haem(1-10) and Haem(45-57). Haem(1-10) competitively binds to/inhibits αT active site (KI  = 1.9 μM) and its potency was enhanced by 10-fold after Phe3 → β-Naphthylalanine exchange. Conversely to full-length haemadin, haem(45-57) displays intrinsic affinity for exosite-I (KD  = 1.6 μM). Hence, we synthesized a peptide in which the sequences 1-9 and 45-57 were joined together through a 3-Glycine spacer to yield haemanorm, a highly potent (KI  = 0.8 nM) inhibitor targeting αT active site and exosite-I. Haemanorm can be regarded as a novel class of hirulog-like αT inhibitors with potential pharmacological applications.

Circular permutation at azurin's active site slows down its folding

Circular permutation (CP) is a technique by which the primary sequence of a protein is rearranged to create new termini. The connectivity of the protein is altered but the overall protein structure generally remains unperturbed. Understanding the effect of CP can help design robust proteins for numerous applications such as in genetic engineering, optoelectronics, and improving catalytic activity. Studies on different protein topologies showed that CP usually affects protein stability as well as unfolding rates. Though a significant number of proteins contain metals or other cofactors, reports of metalloprotein CPs are rare. Thus, we chose a bacterial metalloprotein, azurin, and its CP within the metal-binding site (cpF114). We studied the stabilities, folding, and unfolding rates of apo- and Zn2+-bound CP azurin using fluorescence and circular dichroism. The introduced CP had destabilizing effects on the protein. Also, the folding of the Zn2+-CP protein was much slower than that of the Zn2+-WT or apo-protein. We compared this study to our previously reported azurin-cpN42, where we had observed an equilibrium and kinetic intermediate. cpF114 exhibits an apparent two-state equilibrium unfolding but has an off-pathway kinetic intermediate. Our study hinted at CP as a method to modify the energy landscape of proteins to alter their folding pathways. WT azurin, being a faster folder, may have evolved to optimize the folding rate of metal-bound protein compared to its CPs, albeit all of them have the same structure and function. Our study underscores that protein sequence and protein termini positions are crucial for metalloproteins. TOC Figure. (Top) Zn2+-azurin WT structure (PDB code: 1E67) and 2-D topology diagram of Zn2+-cpF114 azurin. (Bottom) Cartoon diagram representing folding (red arrows) and unfolding (blue arrows) of apo- and Zn2+- WT and cpF114 azurins. The width of the arrows represents the rate of the corresponding processes.

FireProt 2.0: web-based platform for the fully automated design of thermostable proteins

Thermostable proteins find their use in numerous biomedical and biotechnological applications. However, the computational design of stable proteins often results in single-point mutations with a limited effect on protein stability. However, the construction of stable multiple-point mutants can prove difficult due to the possibility of antagonistic effects between individual mutations. FireProt protocol enables the automated computational design of highly stable multiple-point mutants. FireProt 2.0 builds on top of the previously published FireProt web, retaining the original functionality and expanding it with several new stabilization strategies. FireProt 2.0 integrates the AlphaFold database and the homology modeling for structure prediction, enabling calculations starting from a sequence. Multiple-point designs are constructed using the Bron-Kerbosch algorithm minimizing the antagonistic effect between the individual mutations. Users can newly limit the FireProt calculation to a set of user-defined mutations, run a saturation mutagenesis of the whole protein or select rigidifying mutations based on B-factors. Evolution-based back-to-consensus strategy is complemented by ancestral sequence reconstruction. FireProt 2.0 is significantly faster and a reworked graphical user interface broadens the tool's availability even to users with older hardware. FireProt 2.0 is freely available at http://loschmidt.chemi.muni.cz/fireprotweb.

Domain swap facilitates structural transitions of spider silk protein C-terminal domains

Domain swap is a mechanism of protein dimerization where the two interacting domains exchange parts of their structure. Web spiders make use of the process in the connection of C-terminal domains (CTDs) of spidroins, the soluble protein building blocks that form tough silk fibers. Besides providing connectivity and solubility, spidroin CTDs are responsible for inducing structural transitions during passage through an acidified assembly zone within spinning ducts. The underlying molecular mechanisms are elusive. Here, we studied the folding of five homologous spidroin CTDs from different spider species or glands. Four of these are domain-swapped dimers formed by five-helix bundles from spidroins of major and minor ampullate glands. The fifth is a dimer that lacks domain swap, formed by four-helix bundles from a spidroin of a flagelliform gland. Spidroins from this gland do not undergo structural transitions whereas the others do. We found a three-state mechanism of folding and dimerization that was conserved across homologues. In chemical denaturation experiments the native CTD dimer unfolded to a dimeric, partially structured intermediate, followed by full unfolding to denatured monomers. The energetics of the individual folding steps varied between homologues. Contrary to the common belief that domain swap stabilizes protein assemblies, the non-swapped homologue was most stable and folded four orders of magnitude faster than a swapped variant. Domain swap of spidroin CTDs induces an entropic penalty to the folding of peripheral helices, thus unfastening them for acid-induced unfolding within a spinning duct, which primes them for refolding into alternative structures during silk formation.

The variable domain from dynamin-related protein 1 promotes liquid-liquid phase separation that enhances its interaction with cardiolipin-containing membranes

Dynamins are an essential superfamily of mechanoenzymes that remodel membranes and often contain a "variable domain" important for regulation. For the mitochondrial fission dynamin, dynamin-related protein 1, a regulatory role for the variable domain (VD) is demonstrated by gain- and loss-of-function mutations, yet the basis for this is unclear. Here, the isolated VD is shown to be intrinsically disordered and undergo a cooperative transition in the stabilizing osmolyte trimethylamine N-oxide. However, the osmolyte-induced state is not folded and surprisingly appears as a condensed state. Other co-solutes including known molecular crowder Ficoll PM 70, also induce a condensed state. Fluorescence recovery after photobleaching experiments reveal this state to be liquid-like indicating the VD undergoes a liquid-liquid phase separation under crowding conditions. These crowding conditions also enhance binding to cardiolipin, a mitochondrial lipid, which appears to promote phase separation. Since dynamin-related protein 1 is found assembled into discrete punctate structures on the mitochondrial surface, the inference from the present work is that these structures might arise from a condensed state involving the VD that may enable rapid tuning of mechanoenzyme assembly necessary for fission.

Characterization of the folding and binding properties of the PTB domain of FRS2 with phosphorylated and unphosphorylated ligands

PTB (PhosphoTyrosine Binding) domains are protein domains that exert their function by binding phosphotyrosine residues on other proteins. They are commonly found in a variety of signaling proteins and are important for mediating protein-protein interactions in numerous cellular processes. PTB domains can also exhibit binding to unphosphorylated ligands, suggesting that they have additional binding specificities beyond phosphotyrosine recognition. Structural studies have reported that the PTB domain from FRS2 possesses this peculiar feature, allowing it to interact with both phosphorylated and unphosphorylated ligands, such as TrkB and FGFR1, through different topologies and orientations. In an effort to elucidate the dynamic and functional properties of these protein-protein interactions, we provide a complete characterization of the folding mechanism of the PTB domain of FRS2 and the binding process to peptides mimicking specific regions of TrkB and FGFR1. By analyzing the equilibrium and kinetics of PTB folding, we propose a mechanism implying the presence of an intermediate along the folding pathway. Kinetic binding experiments performed at different ionic strengths highlighted the electrostatic nature of the interaction with both peptides. The specific role of single amino acids in early and late events of binding was pinpointed by site-directed mutagenesis. These results are discussed in light of previous experimental works on these protein systems.

Folding stabilities of ribosome-bound nascent polypeptides probed by mass spectrometry

The folding of most proteins occurs during the course of their translation while their tRNA-bound C termini are embedded in the ribosome. How the close proximity of nascent proteins to the ribosome influences their folding thermodynamics remains poorly understood. Here, we have developed a mass spectrometry-based approach for determining the stabilities of nascent polypeptide chains using methionine oxidation as a folding probe. This approach enables quantitative measurement subglobal folding stabilities of ribosome nascent chains within complex protein mixtures and extracts. To validate the methodology, we analyzed the folding thermodynamics of three model proteins (dihydrofolate reductase, chemotaxis protein Y, and DNA polymerase IV) in soluble and ribosome-bound states. The data indicate that the ribosome can significantly alter the stability of nascent polypeptides. Ribosome-induced stability modulations were highly variable among different folding domains and were dependent on localized charge distributions within nascent polypeptides. The results implicated electrostatic interactions between the ribosome surface and nascent polypeptides as the cause of ribosome-induced stability modulations. The study establishes a robust proteomic methodology for analyzing localized stabilities within ribosome-bound nascent polypeptides and sheds light on how the ribosome influences the thermodynamics of protein folding.

Architectural digest: Thermodynamic stability and domain structure of a consensus monomeric globin

Artificial proteins representing the consensus of a set of homologous sequences have attracted attention for their increased thermodynamic stability and conserved activity. Here, we applied the consensus approach to a b-type heme-binding protein to inspect the contribution of a dissociable cofactor to enhanced stability and the chemical consequences of creating a generic heme environment. We targeted the group 1 truncated hemoglobin (TrHb1) subfamily of proteins for their small size (∼120 residues) and ease of characterization. The primary structure, derived from a curated set of ∼300 representative sequences, yielded a highly soluble consensus globin (cGlbN) enriched in acidic residues. Optical and NMR spectroscopies revealed high-affinity heme binding in the expected site and in two orientations. At neutral pH, proximal and distal iron coordination was achieved with a pair of histidine residues, as observed in some natural TrHb1s, and with labile ligation on the distal side. As opposed to studied TrHb1s, which undergo additional folding upon heme binding, cGlbN displayed the same extent of secondary structure whether the heme was associated with the protein or not. Denaturation required guanidine hydrochloride and showed that apo- and holoprotein unfolded in two transitions-the first (occurring with a midpoint of ∼2 M) was shifted to higher denaturant concentration in the holoprotein (∼3.7 M) and reflected stabilization due to heme binding, while the second transition (∼6.2 M) was common to both forms. Thus, the consensus sequence stabilized the protein but exposed the existence of two separately cooperative subdomains within the globin architecture, masked as one single domain in TrHb1s with typical stabilities. The results suggested ways in which specific chemical or thermodynamic features may be controlled in artificial heme proteins.

Protein Binding Leads to Reduced Stability and Solvated Disorder in the Polystyrene Nanoparticle Corona

Understanding the conformation of proteins in the nanoparticle corona has important implications in how organisms respond to nanoparticle-based drugs. These proteins coat the nanoparticle surface, and their properties will influence the nanoparticle's interaction with cell targets and the immune system. While some coronas are thought to be disordered, two key unanswered questions are the degree of disorder and solvent accessibility. Here, using a comprehensive thermodynamic approach, along with supporting spectroscopic experiments, we develop a model for protein corona disorder in polystyrene nanoparticles of varying size. For two different proteins, we find that binding affinity decreases as nanoparticle size increases. The stoichiometry of binding, along with changes in the hydrodynamic size, support a highly solvated, disordered protein corona anchored at a small number of enthalpically-driven attachment sites. The scaling of the stoichiometry vs. nanoparticle size is consistent disordered polymer dimensions. Moreover, we find that proteins are destabilized less severely in the presence of larger nanoparticles, and this is supported by measurements of hydrophobic exposure, which becomes less pronounced at lower curvatures. Our observations hold for flat polystyrene surfaces, which, when controlled for total surface area, have the lowest hydrophobic exposure of all systems. Our model provides an explanation for previous observations of increased amyloid fibrillation rates in the presence of larger nanoparticles, and it may rationalize how cell receptors can recognize protein disorder in therapeutic nanoparticles.

Unraveling the folding and dimerization properties of the human FoxP subfamily of transcription factors

Human FoxP proteins share a highly conserved DNA-binding domain that dimerizes via three-dimensional domain swapping, although showing varying oligomerization propensities among its members. Here, we present an experimental and computational characterization of all human FoxP proteins to unravel how their amino acid substitutions impact their folding and dimerization mechanism. We solved the crystal structure of the forkhead domain of FoxP4 to then perform a comparison across all members, finding that their sequence changes impact not only the structural heterogeneity of their forkhead domains but also the protein-protein association energy barrier. Lastly, we demonstrate that the accumulation of a monomeric intermediate is an oligomerization-dependent feature rather than a common aspect of monomers and dimers in this protein subfamily.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins

A protein sequence encodes its energy landscape - all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs - questioning the differences and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the overall consensus protein is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted region is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order couplings using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.

Ordered-domain unfolding of thermophilic isolated β subunit ATP synthase

The flexibility of the ATP synthase's β subunit promotes its role in the ATP synthase rotational mechanism, but its domains stability remains unknown. A reversible thermal unfolding of the isolated β subunit (Tβ) of the ATP synthase from Bacillus thermophilus PS3, tracked through circular dichroism and molecular dynamics, indicated that Tβ shape transits from an ellipsoid to a molten globule through an ordered unfolding of its domains, preserving the β-sheet residual structure at high temperature. We determined that part of the stability origin of Tβ is due to a transversal hydrophobic array that crosses the β-barrel formed at the N-terminal domain and the Rossman fold of the nucleotide-binding domain (NBD), while the helix bundle of the C-terminal domain is the less stable due to the lack of hydrophobic residues, and thus the more flexible to trigger the rotational mechanism of the ATP synthase.

Protein salvage and repurposing in evolution: Phospholipase D toxins are stabilized by a remodeled scrap of a membrane association domain

The glycerophosphodiester phosphodiesterase (GDPD)-like SMaseD/PLD domain family, which includes phospholipase D (PLD) toxins in recluse spiders and actinobacteria, evolved anciently in bacteria from the GDPD. The PLD enzymes retained the core (β/α)8 barrel fold of GDPD, while gaining a signature C-terminal expansion motif and losing a small insertion domain. Using sequence alignments and phylogenetic analysis, we infer that the C-terminal motif derives from a segment of an ancient bacterial PLAT domain. Formally, part of a protein containing a PLAT domain repeat underwent fusion to the C terminus of a GDPD barrel, leading to attachment of a segment of a PLAT domain, followed by a second complete PLAT domain. The complete domain was retained only in some basal homologs, but the PLAT segment was conserved and repurposed as the expansion motif. The PLAT segment corresponds to strands β7-β8 of a β-sandwich, while the expansion motif as represented in spider PLD toxins has been remodeled as an α-helix, a β-strand, and an ordered loop. The GDPD-PLAT fusion led to two acquisitions in founding the GDPD-like SMaseD/PLD family: (1) a PLAT domain that presumably supported early lipase activity by mediating membrane association, and (2) an expansion motif that putatively stabilized the catalytic domain, possibly compensating for, or permitting, loss of the insertion domain. Of wider significance, messy domain shuffling events can leave behind scraps of domains that can be salvaged, remodeled, and repurposed.

The variable domain from the mitochondrial fission mechanoenzyme Drp1 promotes liquid-liquid phase separation

Dynamins are an essential superfamily of mechanoenzymes that remodel membranes and often contain a "variable domain" (VD) important for regulation. For the mitochondrial fission dynamin, Drp1, a regulatory role for the VD is demonstrated by mutations that can elongate, or fragment, mitochondria. How the VD encodes inhibitory and stimulatory activity is unclear. Here, isolated VD is shown to be intrinsically disordered (ID) yet undergoes a cooperative transition in the stabilizing osmolyte TMAO. However, the TMAO stabilized state is not folded and surprisingly appears as a condensed state. Other co-solutes including known molecular crowder Ficoll PM 70, also induce a condensed state. Fluorescence recovery after photobleaching experiments reveal this state to be liquid-like indicating the VD undergoes a liquid-liquid phase separation under crowding conditions. These crowding conditions also enhance binding to cardiolipin, a mitochondrial lipid, raising the possibility that phase separation may enable rapid tuning of Drp1 assembly necessary for fission.

Monomeric and dimeric states of human ZO1-PDZ2 are functional partners of the SARS-CoV-2 E protein

The Envelope (E) protein of SARS-CoV-2 plays a key role in virus maturation, assembly, and virulence mechanisms. The E protein is characterized by the presence of a PDZ-binding motif (PBM) at its C-terminus that allows it to interact with several PDZ-containing proteins in the intracellular environment. One of the main binding partners of the SARS-CoV-2 E protein is the PDZ2 domain of ZO1, a protein with a crucial role in the formation of epithelial and endothelial tight junctions (TJs). In this work, through a combination of analytical ultracentrifugation analysis and equilibrium and kinetic folding experiments, we show that ZO1-PDZ2 domain is able to fold in a monomeric state, an alternative form to the dimeric conformation that is reported to be functional in the cell for TJs assembly. Importantly, surface plasmon resonance (SPR) data indicate that the PDZ2 monomer is fully functional and capable of binding the C-terminal portion of the E protein of SARS-CoV-2, with a measured affinity in the micromolar range. Moreover, we present a detailed computational analysis of the complex between the C-terminal portion of E protein with ZO1-PDZ2, both in its monomeric conformation (computed as a high confidence AlphaFold2 model) and dimeric conformation (obtained from the Protein Data Bank), by using both polarizable and nonpolarizable simulations. Together, our results indicate both the monomeric and dimeric states of PDZ2 to be functional partners of the E protein, with similar binding mechanisms, and provide mechanistic and structural information about a fundamental interaction required for the replication of SARS-CoV-2.

Critical Beginnings: Selective Tuning of Solubility and Structural Accuracy of Newly Synthesized Proteins by the Hsp70 Chaperone System

Proteins are particularly prone to aggregation immediately after release from the ribosome, and it is therefore important to elucidate the role of chaperones during these key steps of protein life. The Hsp70 and trigger factor (TF) chaperone systems interact with nascent proteins during biogenesis and immediately post-translationally. It is unclear, however, whether these chaperones can prevent formation of soluble and insoluble aggregates. Here, we address this question by monitoring the solubility and structural accuracy of globin proteins biosynthesized in an Escherichia coli cell-free system containing different concentrations of the bacterial Hsp70 and TF chaperones. We find that Hsp70 concentrations required to grant solubility to newly synthesized proteins are extremely sensitive to client-protein sequence. Importantly, Hsp70 concentrations yielding soluble client proteins are insufficient to prevent formation of soluble aggregates. In fact, for some aggregation-prone protein variants, avoidance of soluble-aggregate formation demands Hsp70 concentrations that exceed cellular levels in E. coli. In all, our data highlight the prominent role of soluble aggregates upon nascent-protein release from the ribosome and show the limitations of the Hsp70 chaperone system in the case of highly aggregation-prone proteins. These results demonstrate the need to devise better strategies to prevent soluble-aggregate formation upon release from the ribosome.

Metal-binding and circular permutation-dependent thermodynamic and kinetic stability of azurin

Native topology is known to determine the folding kinetics and the energy landscape of proteins. Furthermore, the circular permutation (CP) of proteins alters the order of the secondary structure connectivity while retaining the three-dimensional structure, making it an elegant and powerful approach to altering native topology. Previous studies elucidated the influence of CP in proteins with different folds such as Greek key β-barrel, β-sandwich, β-α-β, and all α-Greek key. CP mainly affects the protein stability and unfolding kinetics, while folding kinetics remains mostly unaltered. However, the effect of CP on metalloproteins is yet to be elaborately studied. The active site of metalloproteins poses an additional complexity in studying protein folding. Here, we investigate a CP variant (cpN42) of azurin-in both metal-free and metal-bound (holo) forms. As observed earlier in other proteins, apo-forms of wild-type (WT) and cpN42 fold with similar rates. In contrast, zinc-binding accelerates the folding of WT but decelerates the folding of cpN42. On zinc-binding, the spontaneous folding rate of WT increases by >250 times that of cpN42, which is unprecedented and the highest for any CP to date. On the other hand, zinc-binding reduces the spontaneous unfolding rate of cpN42 by ~100 times, making the WT and CP azurins unfold at similar rates. Our study demonstrates metal binding as a novel way to modulate the unfolding and folding rates of CPs compared to their WT counterparts. We hope our study increases the understanding of the effect of CP on the folding mechanism and energy landscape of metalloproteins.

Multiple Modes of Zinc Binding to Histatin 5 Revealed by Buffer-Independent Thermodynamics

Histatin 5 (Hist5) is an antimicrobial peptide found in human saliva as part of the innate immune system. Hist5 can bind several metal ions in vitro, and Zn²⁺ has been shown to function as an inhibitory switch to regulate the peptide's biological activity against the opportunistic fungal pathogen Candida albicans in cell culture. Here, we studied Zn²⁺ binding to Hist5 at four temperatures from 15 to 37 °C using isothermal titration calorimetry to obtain thermodynamic parameters that were corrected for competing buffer effects. Hist5 bound Zn²⁺ with a buffer-dependent association constant of ∼10⁵ M^-1 and a buffer-independent association constant of ∼6 × 10⁶ M^-1 at pH 7.4 and at all temperatures tested. Zn²⁺ binding was both enthalpically and entropically favorable, with larger entropic contributions at 15 °C and larger enthalpic contributions at 37 °C. Additionally, the Zn:Hist5 binding stoichiometry increased from 1:1 to 2:1 as temperature increased. The enthalpy-entropy compensation and the variable stoichiometry lead us to propose a model in which the Zn-Hist5 complex exists in an equilibrium between two distinct binding modes with different Zn:Hist5 stoichiometries. The in-depth thermodynamic analysis presented herein may help illuminate the biophysical basis for Zn-dependent changes in the antifungal activity of Hist5.

Protein Stability─Analysis of Heat and Cold Denaturation without and with Unfolding Models

Protein stability is important in many areas of life sciences. Thermal protein unfolding is investigated extensively with various spectroscopic techniques. The extraction of thermodynamic properties from these measurements requires the application of models. Differential scanning calorimetry (DSC) is less common, but is unique as it measures directly a thermodynamic property, that is, the heat capacity C_p(T). The analysis of C_p(T) is usually performed with the chemical equilibrium two-state model. This is not necessary and leads to incorrect thermodynamic consequences. Here we demonstrate a straightforward model-independent evaluation of heat capacity experiments in terms of protein unfolding enthalpy ΔH(T), entropy ΔS(T), and free energy ΔG(T)). This now allows the comparison of the experimental thermodynamic data with the predictions of different models. We critically examined the standard chemical equilibrium two-state model, which predicts a positive free energy for the native protein, and diverges distinctly from the experimental temperature profiles. We propose two new models which are equally applicable to spectroscopy and calorimetry. The Θ_U(T)-weighted chemical equilibrium model and the statistical-mechanical two-state model provide excellent fits of the experimental data. They predict sigmoidal temperature profiles for enthalpy and entropy, and a trapezoidal temperature profile for the free energy. This is illustrated with experimental examples for heat and cold denaturation of lysozyme and β-lactoglobulin. We then show that the free energy is not a good criterion to judge protein stability. More useful parameters are discussed, including protein cooperativity. The new parameters are embedded in a well-defined thermodynamic context and are amenable to molecular dynamics calculations.

Modeling the role of charged residues in thermophilic proteins by rotamer and dynamic cross correlation analysis

Discerning the determinants of protein thermostability is very important both from the theoretical and applied perspective. Different lines of evidence seem to indicate that a dynamical network of salt bridges/charged residues plays a fundamental role in the thermostability of enzymes. In this work, we applied measures of dynamic variance, like the Gini coefficients, Kullback-Leibler (KL) divergence and dynamic cross correlation (DCC) coefficients to compare the behavior of 3 pairs of homologous proteins from the thermophilic bacterium Thermus thermophilus and mesophilic Escherichia coli. Molecular dynamic (MD) simulations of these proteins were performed at 303 K and 363 K. In the characterization of their side chain rotamer distributions, the corresponding Gini coefficients and KL-divergence both revealed significant correlations with temperature. Similarly, a DCC analysis revealed a higher trend to de-correlate the movement of charged residues at higher temperatures in the thermophilic proteins, when compared with their mesophilic homologues. These results highlight the importance of dynamic electrostatic network interactions for the thermostability of enzymes.

Salt-Specific Suppression of the Cold Denaturation of Thermophilic Multidomain Initiation Factor 2

Thermophilic proteins and enzymes are attractive for use in industrial applications due to their resistance against heat and denaturants. Here, we report on a thermophilic protein that is stable at high temperatures (Ttrs, hot 67 °C) but undergoes significant unfolding at room temperature due to cold denaturation. Little is known about the cold denaturation of thermophilic proteins, although it can significantly limit their applications. We investigated the cold denaturation of thermophilic multidomain protein translation initiation factor 2 (IF2) from Thermus thermophilus. IF2 is a GTPase that binds to ribosomal subunits and initiator fMet-tRNAfMet during the initiation of protein biosynthesis. In the presence of 9 M urea, measurements in the far-UV region by circular dichroism were used to capture details about the secondary structure of full-length IF2 protein and its domains during cold and hot denaturation. Cold denaturation can be suppressed by salt, depending on the type, due to the decreased heat capacity. Thermodynamic analysis and mathematical modeling of the denaturation process showed that salts reduce the cooperativity of denaturation of the IF2 domains, which might be associated with the high frustration between domains. This characteristic of high interdomain frustration may be the key to satisfying numerous diverse contacts with ribosomal subunits, translation factors, and tRNA.

FkpA enhances membrane protein folding using an extensive interaction surface

Outer membrane protein (OMP) biogenesis in gram-negative bacteria is managed by a network of periplasmic chaperones that includes SurA, Skp, and FkpA. These chaperones bind unfolded OMPs (uOMPs) in dynamic conformational ensembles to suppress aggregation, facilitate diffusion across the periplasm, and enhance folding. FkpA primarily responds to heat-shock stress, but its mechanism is comparatively understudied. To determine FkpA chaperone function in the context of OMP folding, we monitored the folding of three OMPs and found that FkpA, unlike other periplasmic chaperones, increases the folded yield but decreases the folding rate of OMPs. The results indicate that FkpA behaves as a chaperone and not as a folding catalyst to influence the OMP folding trajectory. Consistent with the folding assay results, FkpA binds all three uOMPs as determined by sedimentation velocity (SV) and photo-crosslinking experiments. We determine the binding affinity between FkpA and uOmpA171 by globally fitting SV titrations and find it to be intermediate between the known affinities of Skp and SurA for uOMP clients. Notably, complex formation steeply depends on the urea concentration, suggesting an extensive binding interface. Initial characterizations of the complex using photo-crosslinking indicate that the binding interface spans the entire FkpA molecule. In contrast to prior findings, folding and binding experiments performed using subdomain constructs of FkpA demonstrate that the full-length chaperone is required for full activity. Together these results support that FkpA has a distinct and direct effect on OMP folding that it achieves by utilizing an extensive chaperone-client interface to tightly bind clients.

Exploration of seryl tRNA synthetase to identify potent inhibitors against leishmanial parasites

Aminoacyl-tRNA synthetases are crucial enzymes for cellular protein metabolism and have been considered as an attractive target for development of new antimicrobials. In the current study, seryl tRNA synthetase of Leishmania donovani (LdSerRS) and its mutants were purified and characterized through biochemical and structural methods. Purified LdSerRS was found to be enzymatically active and exhibited more alpha helices in secondary structure. The enzymatic activity of purified protein was observed as highest near physiological temperature and pH. Mutation in ATP binding residues (R295 and E297) demonstrated reduction in the affinity for cofactor with no significant deviation in secondary structure. In vitro inhibition studies with ureidosulfocoumarin derivatives helped to identify Comp 5l as a specific inhibitor for leishmanial SerRS that showed lesser potency towards purified HsSerRS. The identified compound presented competitive mode of inhibition for LdSerRS and also revealed druglikeness along with very low toxicity for human macrophages. Structural analysis of protein and ligand complex depicted the binding of Comp 5l into the cofactor binding site of LdSerRS with high affinity succeeded by validation employing molecular dynamics simulations. Altogether, our study presents a promising scaffold to explore small molecules to target the enzymatic activity of leishmanial SerRS to develop the specific therapeutics.

Clamp Interactions with +3/+6 Duplex and Upstream-to-Downstream Allosteric Effects in Late Steps of Forming a Stable RNA Polymerase-Promoter Open Complex

Stable 37 °C open complexes (OC) of E. coli RNA polymerase (RNAP) at λP_R and T7A1 promoters form at similar rates but have very different lifetimes. To understand the downstream interactions responsible for OC lifetime, how promoter sequence directs them and when they form, we report lifetimes of stable OC and unstable late (I₂) intermediates for promoters with different combinations of λP_R (L) and T7A1 (T) discriminators, core promoters and UP elements. I₂ lifetimes are similarly short, while stable OC lifetimes differ greatly, determined largely by the discriminator and modulated by core-promoter and UP elements. The free energy change ΔG₃^o for I₂ → stable OC is approximately -4 kcal more favorable for L-discriminator than for T-discriminator promoters. Downstream-truncation at +6 (DT+6) greatly destabilizes OC at L-discriminator but not T-discriminator promoters, making all ΔG₃^o values similar (approximately -4 kcal). Urea reduces OC lifetime greatly by affecting ΔG₃^o. We deduce that urea acts by disfavoring coupled folding of key elements of the β'-clamp, that I₂ is an open-clamp OC, and that clamp-closing in I₂ → stable OC involves coupled folding. Differences in ΔG₃^o between downstream-truncated and full-length promoters yield contributions to ΔG₃^o from interactions with downstream mobile elements (DME) including β-lobe and β'-jaw, more favorable for L-discriminator than for T-discriminator promoters. We deduce how competition between far-downstream DNA and σ⁷⁰ region 1.1 affects ΔG₃^o values. We discuss variant-specific ΔG₃^o contributions in terms of the allosteric network by which differences in discriminator and -10 sequence are sensed and transmitted downstream to affect DME-duplex interactions in I₂ → stable OC.

Collagen-like Motifs of SasG: A Novel Fold for Protein Mechanical Strength

The Staphylococcus aureus surface protein G (SasG) is associated with host colonisation and biofilm formation. As colonisation occurs at the liquid-substrate interface bacteria are subject to a myriad of external forces and, presumably as a consequence, SasG displays extreme mechanical strength. This mechanical phenotype arises from the B-domain; a repetitive region composed of alternating E and G5 subdomains. These subdomains have an unusual structure comprising collagen-like regions capped by triple-stranded β-sheets. To identify the determinants of SasG mechanical strength, we characterised the mechanical phenotype and thermodynamic stability of 18 single substitution variants of a pseudo-wildtype protein. Visualising the mechanically-induced transition state at a residue-level by ϕ-value analysis reveals that the main force-bearing regions are the N- and C-terminal 'Mechanical Clamps' and their side-chain interactions. This is tailored by contacts at the pseudo-hydrophobic core interface. We also describe a novel mechanical motif - the collagen-like region and show that glycine to alanine substitutions, analogous to those found in Osteogenesis Imperfecta (brittle bone disease), result in a significantly reduced mechanical strength.

Protein Unfolding—Thermodynamic Perspectives and Unfolding Models

We review the key steps leading to an improved analysis of thermal protein unfolding. Thermal unfolding is a dynamic cooperative process with many short-lived intermediates. Protein unfolding has been measured by various spectroscopic techniques that reveal structural changes, and by differential scanning calorimetry (DSC) that provides the heat capacity change Cp(T). The corresponding temperature profiles of enthalpy ΔH(T), entropy ΔS(T), and free energy ΔG(T) have thus far been evaluated using a chemical equilibrium two-state model. Taking a different approach, we demonstrated that the temperature profiles of enthalpy ΔH(T), entropy ΔS(T), and free energy ΔG(T) can be obtained directly by a numerical integration of the heat capacity profile Cp(T). DSC thus offers the unique possibility to assess these parameters without resorting to a model. These experimental parameters now allow us to examine the predictions of different unfolding models. The standard two-state model fits the experimental heat capacity peak quite well. However, neither the enthalpy nor entropy profiles (predicted to be almost linear) are congruent with the measured sigmoidal temperature profiles, nor is the parabolic free energy profile congruent with the experimentally observed trapezoidal temperature profile. We introduce three new models, an empirical two-state model, a statistical–mechanical two-state model and a cooperative statistical-mechanical multistate model. The empirical model partially corrects for the deficits of the standard model. However, only the two statistical–mechanical models are thermodynamically consistent. The two-state models yield good fits for the enthalpy, entropy and free energy of unfolding of small proteins. The cooperative statistical–mechanical multistate model yields perfect fits, even for the unfolding of large proteins such as antibodies.

Determining the Conformational Stability of a Protein from Urea and Thermal Unfolding Curves

This article contains protocols for determining the conformational stability of a globular protein from either urea or thermal unfolding curves. Circular dichroism is the optical spectroscopic technique most commonly used to monitor protein unfolding. These protocols describe how to analyze data from an unfolding curve to obtain the thermodynamic parameters necessary to calculate conformational stability, and how to determine differences in stability between protein variants. Curr. Protoc. Protein Sci. 71:28.4.1-28.4.14. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Determining protein conformational stability from urea-induced unfolding curves Support Protocol 1: Preparing a urea stock solution Support Protocol 2: Analyzing urea unfolding curves Basic Protocol 2: Determining the conformational stability of a protein from thermal unfolding curves Support Protocol 3: Analyzing thermal unfolding curves Support Protocol 4: Determining differences in conformational stability for protein variants.

Ocular protein optineurin shows reversibility from unfolded states and exhibits chaperone-like activity

Optineurin (OPTN) is a multifunctional, ubiquitously expressed cytoplasmic protein, mutants of which are associated with primary open-angle glaucoma (POAG) and amyotrophic lateral sclerosis (ALS). The most abundant heat shock protein crystallin, known for its remarkable thermodynamic stability and chaperoning activity, allows ocular tissues to withstand stress. The presence of OPTN in ocular tissues is intriguing. Interestingly, OPTN also harbors heat shock elements in its promoter region. Sequence analysis of OPTN exhibits intrinsically disordered regions and nucleic acid binding domains. These properties hinted that OPTN might be endowed with sufficient thermodynamic stability and chaperoning activity. However, these attributes of OPTN have not yet been explored. Here, we studied these properties through thermal and chemical denaturation experiments and monitored the processes using CD, fluorimetry, differential scanning calorimetry, and dynamic light scattering. We found that upon heating, OPTN reversibly forms higher-order multimers. OPTN also displayed a chaperone-like function by reducing the thermal aggregation of bovine carbonic anhydrase. It regains its native secondary structure, RNA-binding property, and melting temperature (T _m) after refolding from a thermally as well as chemically denatured state. From our data, we conclude that OPTN, with its unique ability to revert from the stress-mediated unfolded state and its unique chaperoning function, is a valuable protein of the ocular tissues.

Temperature dependent mechanical unfolding and refolding of a protein studied by thermo-regulated optical tweezers

Temperature is a useful system variable to gather kinetic and thermodynamic information from proteins. Usually, free energy and the associated entropic and enthalpic contributions are obtained by quantifying the conformational equilibrium based on melting experiments performed in bulk conditions. Such experiments are suitable only for those small single-domain proteins whose side reactions of irreversible aggregation are unlikely to occur. Here, we avoid aggregation by pulling single-protein molecules in a thermo-regulated optical tweezers. Thus, we are able to explore the temperature dependence of the thermodynamic and kinetic parameters of MJ0366 from Methanocaldococcus jannaschii at the single-molecule level. By performing force-ramp experiments between 2°C and 40°C, we found that MJ0366 has a nonlinear dependence of free energy with temperature and a specific heat change of 2.3 ± 1.2 kcal/mol∗K. These thermodynamic parameters are compatible with a two-state unfolding/refolding mechanism for MJ0366. However, the kinetics measured as a function of the temperature show a complex behavior, suggesting a three-state folding mechanism comprising a high-energy intermediate state. The combination of two perturbations, temperature and force, reveals a high-energy species in the folding mechanism of MJ0366 not detected in force-ramp experiments at constant temperature.

Understanding the molecular basis of folding cooperativity through a comparative analysis of a multidomain protein and its isolated domains

Although cooperativity is a well-established and general property of folding, our current understanding of this feature in multidomain folding is still relatively limited. In fact, there are contrasting results indicating that the constituent domains of a multidomain protein may either fold independently on each other or exhibit interdependent supradomain phenomena. To address this issue, here we present the comparative analysis of the folding of a tandem repeat protein, comprising two contiguous PDZ domains, in comparison to that of its isolated constituent domains. By analyzing in detail the equilibrium and kinetics of folding at different experimental conditions, we demonstrate that despite each of the PDZ domains in isolation being capable of independent folding, at variance with previously characterized PDZ tandem repeats, the full-length construct folds and unfolds as a single cooperative unit. By exploiting quantitatively, the comparison of the folding of the tandem repeat to those observed for its constituent domains, as well as by characterizing a truncated variant lacking a short autoinhibitory segment, we successfully rationalize the molecular basis of the observed cooperativity and attempt to infer some general conclusions for multidomain systems.

Conformationally responsive dyes enable protein-adaptive differential scanning fluorimetry

Flexible in vitro methods alter the course of biological discoveries. Differential Scanning Fluorimetry (DSF) is a particularly versatile technique which reports protein thermal unfolding via fluorogenic dye. However, applications of DSF are limited by widespread protein incompatibilities with the available DSF dyes. Here, we enable DSF applications for 66 of 70 tested proteins (94%) including 10 from the SARS-CoV2 virus using a chemically diverse dye library, Aurora, to identify compatible dye-protein pairs in high throughput. We find that this protein-adaptive DSF platform (paDSF) not only triples the previous protein compatibility, but also fundamentally extends the processes observable by DSF, including interdomain allostery in O-GlcNAc Transferase (OGT). paDSF enables routine measurement of protein stability, dynamics, and ligand binding.

Molecular recognition requires dimerization of a VHH antibody

Camelid heavy-chain-only antibodies are a unique class of antibody that possesses only a single variable domain (termed VHH) for antigen recognition. Despite their apparent canonical mechanism of target recognition, where a single VHH domain binds a single target, an anti-caffeine VHH has been observed to possess 2:1 stoichiometry. Here, the structure of the anti-caffeine VHH/caffeine complex enabled the generation and biophysical analysis of variants that were used to better understand the role of VHH homodimerization in caffeine recognition. VHH interface mutants and caffeine analogs, which were examined to probe the mechanism of caffeine binding, suggested caffeine recognition is only possible with the VHH dimer species. Correspondingly, in the absence of caffeine, the anti-caffeine VHH was found to form a dimer with a dimerization constant comparable to that observed with VH:VL domains in conventional antibody systems, which was most stable near physiological temperature. While the VHH:VHH dimer structure (at 1.13 Å resolution) is reminiscent of conventional VH:VL heterodimers, the homodimeric VHH possesses a smaller angle of domain interaction, as well as a larger amount of apolar surface area burial. To test the general hypothesis that the short complementarity-determining region-3 (CDR3) may help drive VHH:VHH homodimerization, an anti-picloram VHH domain containing a short CDR3 was generated and characterized, which revealed it also existed as dimer species in solution. These results suggest homodimer-driven recognition may represent a more common method of VHH ligand recognition, opening opportunities for novel VHH homodimer affinity reagents and helping to guide their use in chemically induced dimerization applications.

Fused in sarcoma undergoes cold denaturation: Implications for phase separation

The mediation of liquid-liquid phase separation (LLPS) for fused in sarcoma (FUS) protein is generally attributed to the low-complexity, disordered domains and is enhanced at low temperature. The role of FUS folded domains on the LLPS process remains relatively unknown since most studies are mainly based on fragmented FUS domains. Here, we investigate the effect of metabolites on full-length (FL) FUS LLPS using turbidity assays and differential interference contrast (DIC) microscopy, and explore the behavior of the folded domains by nuclear magnetic resonance (NMR) spectroscopy. FL FUS LLPS is maximal at low concentrations of glucose and glutamate, moderate concentrations of NaCl, Zn²⁺ , and Ca²⁺ and at the isoelectric pH. The FUS RNA recognition motif (RRM) and zinc-finger (ZnF) domains are found to undergo cold denaturation above 0°C at a temperature that is determined by the conformational stability of the ZnF domain. Cold unfolding exposes buried nonpolar residues that can participate in LLPS-promoting hydrophobic interactions. Therefore, these findings constitute the first evidence that FUS globular domains may have an active role in LLPS under cold stress conditions and in the assembly of stress granules, providing further insight into the environmental regulation of LLPS.

Guanidinium–amino acid hydrogen-bonding interactions in protein crystal structures: implications for guanidinium-induced protein denaturation

Structural analysis of guanidinium–amino acid interaction pairs in protein crystal structures is coupled with an effective scheme for classifying the optimized pairs, to gain understanding of the guanidinium:protein hydrogen bonding modes.

The cooperative folding of annexin A2 relies on a transient nonnative intermediate

Annexins (Anxs) are a family of highly homologous proteins that bind and aggregate lipid vesicles in the presence of calcium. All members of the family contain a variable N-terminus determining specific functions, followed by a conserved core region responsible for the general calcium-dependent lipid-binding property. The core structure consists of four homologous domains (DI-DIV), each consisting of a right-handed super-helix of five α-helices. We present data from a combination of site-directed mutagenesis, NMR, and circular dichroism showing that the G25-D34 region of the N-terminus as well as the contacts between residues D38A, R63A, and Q67A of AnxA2-DI are crucial for the autonomous folding and stability of DI of AnxA2. However, we also show that the folding of the full-length protein is very robust in that mutations and truncations that disrupted the folding of AnxA2-DI did not abolish the folding of full-length AnxA2, only lowering its thermal stability. This robustness of the folding of full-length AnxA2 is likely to be mediated by the existence of at least one transient nonnative intermediate as suggested by our kinetic data using stopped-flow fluorescence experiments. We also show that hydrophobic amino acids in AnxA2-DI involved in interfacial contacts with AnxA2-DIV are important for the cooperative folding and stability of the full-length protein. Mutating all of the V57E-V98R-G101Y residues in AnxA2-DI did not affect the folding of the domain, only its stability, but prevented the cooperative folding of the full-length protein. Our collective results favor a highly cooperative and robust folding process mediated by alternative intermediate steps. Since AnxA2 is a multifunctional protein involved in several steps of the progression of cell transformation, these data on structure and folding pathways are therefore crucial to designing anticancer drugs targeting AnxA2.

Recombinant Globular Domain of TcpA Pilin from Vibrio cholerae El Tor: Recovery from Inclusion Bodies and Structural Characterization

The production of recombinant proteins in Escherichia coli cells is often hampered by aggregation of newly synthesized proteins and formation of inclusion bodies. Here we propose the use of transverse urea gradient electrophoresis (TUGE) in testing the capability of folding of a recombinant protein from inclusion bodies dissolved in urea. A plasmid encoding the amino acid sequence 55–224 of TcpA pilin (C-terminal globular domain: TcpA-C) from Vibrio cholerae El Tor enlarged by a His-tag on its N-terminus was expressed in E. coli cells. The major fraction (about 90%) of the target polypeptide was detected in cell debris. The polypeptide was isolated from the soluble fraction and recovered from inclusion bodies after their urea treatment. Some structural properties of the polypeptide from each sample proved identical. The refolding protocol was developed on the basis of TUGE data and successfully used for the protein large-scale recovery from inclusion bodies. Spectral, hydrodynamic, and thermodynamic characteristics of the recombinant TcpA recovered from inclusion bodies indicate the presence of a globular conformation with a pronounced secondary structure and a rigid tertiary structure, which is promising for the design of immunodiagnostics preparations aimed to assess the pilin level in different strains of V. cholerae and to develop cholera vaccines.