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

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.

Insights into the Conformation and Self-Association of a Concentrated Monoclonal Antibody using Isothermal Chemical Denaturation and Nuclear Magnetic Resonance

The purpose of this investigation was to highlight the utility of nuclear magnetic resonance (NMR) as a multi-attribute method for the characterization of therapeutic antibodies. In this case study, we compared results from isothermal chemical denaturation (ICD) and NMR with standard methods to relate conformational states of a model monoclonal antibody (mAb1) with protein-protein interactions (PPI) that lead to self - association in concentrated solutions. The increase in aggregation rate and relative viscosity for mAb1 was found to be both concentration and pH dependent. The free energy of unfolding (∆G⁰) from ICD and thermal analysis in dilute solutions indicated that although the native state predominated between pH 4 - pH 7, it was disrupted at the CH₂ and unfolded noncooperatively under acidic conditions. One-dimensional (1D) ¹H NMR and two-dimensional (2D) ¹³C-¹H NMR performed, in concentrated solutions, confirmed that PPI between pH 4-7 occurred while mAb1 was in the native state. NMR corroborated that mAb1 maintained a dominant native state at formulation-relevant conditions at the tested pH range, had increased global molecular tumbling dynamics at lower pH and confirmed increased PPI at higher pH conditions. This report aligns and compares typical characterization of an IgG1 with assessment of structure by NMR and provided a more precise assessment and deeper insight into the conformation of an IgG1 in concentrated solutions.

Improving the alkaline stability of pepsin through rational protein design using renin, an alkaline-stable aspartic protease, as a structural and functional reference

The present study sought to identify the structural determinants of aspartic protease structural stability and activity at elevated pH. Various hypotheses have been published regarding the features responsible for the unusual alkaline structural stability of renin, however, few structure-function studies have verified these claims. Using pepsin as a model system, and renin as a template for functional and structural alkaline stability, a rational re-design of pepsin was undertaken to identify residues contributing to the alkaline instability of pepsin-like aspartic proteases in regards to both structure and function. We constructed 13 mutants based on this strategy. Among them, mutants D159 L and D60A led to an increase in activity at elevated pH levels (p ≤ 0.05) and E4V and H53F were shown to retain native-like structure at elevated pH (p ≤ 0.05). Previously suggested carboxyl groups Asp¹¹, Asp¹¹⁸, and Glu¹³ were individually shown not to be responsible for the structural instability or lack of activity at neutral pH in pepsin. The importance of the β-barrel to structural stability was highlighted as the majority of the stabilizing residues identified, and 39% of the weakly conserved residues in the N-terminal lobe, were located in β-sheet strands of the barrel. The results of the present study indicate that alkaline stabilization of pepsin will require reduction of electrostatic repulsions and an improved understanding of the role of the hydrogen bonding network of the characteristic β-barrel.

Synergistic stabilization of a double mutant in chymotrypsin inhibitor 2 from a library screen in E. coli

Most single point mutations destabilize folded proteins. Mutations that stabilize a protein typically only have a small effect and multiple mutations are often needed to substantially increase the stability. Multiple point mutations may act synergistically on the stability, and it is often not straightforward to predict their combined effect from the individual contributions. Here, we have applied an efficient in-cell assay in E. coli to select variants of the barley chymotrypsin inhibitor 2 with increased stability. We find two variants that are more than 3.8 kJ mol-1 more stable than the wild-type. In one case, the increased stability is the effect of the single substitution D55G. The other case is a double mutant, L49I/I57V, which is 5.1 kJ mol-1 more stable than the sum of the effects of the individual mutations. In addition to demonstrating the strength of our selection system for finding stabilizing mutations, our work also demonstrate how subtle conformational effects may modulate stability.

Interconversion of Unexpected Thiol States Affects the Stability, Structure, and Dynamics of Antibody Engineered for Site-Specific Conjugation

Antibody-drug conjugates have become one of the most actively developed classes of drugs in recent years. Their great potential comes from combining the strengths of large and small molecule therapeutics: the exquisite specificity of antibodies and the highly potent nature of cytotoxic compounds. More recently, the approach of engineering antibody-drug conjugate scaffolds to achieve highly controlled drug to antibody ratios has focused on substituting or inserting cysteines to facilitate site-specific conjugation. Herein, we characterize an antibody scaffold engineered with an inserted cysteine that formed an unexpected disulfide bridge during manufacture. A combination of mass spectrometry and biophysical techniques have been used to understand how the additional disulfide bridge forms, interconverts, and changes the stability and structural dynamics of the antibody intermediate. This quantitative and structurally resolved model of the local and global changes in structure and dynamics associated with the engineering and subsequent disulfide-bonded variant can assist future engineering strategies.

The Right-Handed Parallel β-Helix Topology of Erwinia chrysanthemi Pectin Methylesterase Is Intimately Associated with Both Sequential Folding and Resistance to High Pressure

The complex topologies of large multi-domain globular proteins make the study of their folding and assembly particularly demanding. It is often characterized by complex kinetics and undesired side reactions, such as aggregation. The structural simplicity of tandem-repeat proteins, which are characterized by the repetition of a basic structural motif and are stabilized exclusively by sequentially localized contacts, has provided opportunities for dissecting their folding landscapes. In this study, we focus on the Erwinia chrysanthemi pectin methylesterase (342 residues), an all-β pectinolytic enzyme with a right-handed parallel β-helix structure. Chemicals and pressure were chosen as denaturants and a variety of optical techniques were used in conjunction with stopped-flow equipment to investigate the folding mechanism of the enzyme at 25 °C. Under equilibrium conditions, both chemical- and pressure-induced unfolding show two-state transitions, with average conformational stability (ΔG° = 35 ± 5 kJ·mol⁻¹) but exceptionally high resistance to pressure (P_m = 800 ± 7 MPa). Stopped-flow kinetic experiments revealed a very rapid (τ < 1 ms) hydrophobic collapse accompanied by the formation of an extended secondary structure but did not reveal stable tertiary contacts. This is followed by three distinct cooperative phases and the significant population of two intermediate species. The kinetics followed by intrinsic fluorescence shows a lag phase, strongly indicating that these intermediates are productive species on a sequential folding pathway, for which we propose a plausible model. These combined data demonstrate that even a large repeat protein can fold in a highly cooperative manner.

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

•

The TIM barrel is a versatile fold to understand structure-stability relationships.

•

A collection of de novo TIM barrels with improved hydrophobic cores was designed.

•

DeNovoTIMs are reversible in chemical and thermal unfolding, which is uncommon in TIM barrels.

•

Epistatic effects play a central role in DeNovoTIMs stabilization.

•

DeNovoTIMs navigate a previously uncharted region of the stability landscape.

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.

Dissecting Inter-domain Cooperativity in the Folding of a Multi Domain Protein

Correct protein folding underlies all cellular functions. While there are detailed descriptions and a good understanding of protein folding pathways for single globular domains there is a paucity of quantitative data regarding folding of multidomain proteins. We have here investigated the folding of a three-domain supramodule from the protein PSD-95, consisting of one PDZ domain, one SH3 domain and one guanylate kinase-like (GK) domain. This supramodule has previously been shown to work as one functional unit with regard to ligand binding. We used equilibrium and kinetic folding experiments to demonstrate that the PDZ domain folds faster and independently from the SH3-GK tandem, which folds as one cooperative unit. However, concurrent folding of the PDZ domain slows down folding of SH3-GK by non-native interactions, resulting in an off-pathway folding intermediate. Our data contribute to an emerging description of multidomain protein folding in which individual domains cannot a priori be viewed as separate folding units.

Protons Are Fast and Smart; Proteins Are Slow and Dumb: On the Relationship of Electrospray Ionization Charge States and Conformations

We present simple considerations of how differences in time scales of motions of protons, the lightest and fastest chemical moiety, and the much longer time scales associated with the dynamics of proteins, among the heaviest and slowest analytes, may allow many protein conformations from solution to be kinetically trapped during the process of electrospraying protein solutions into the gas phase. In solution, the quantum nature of protons leads them to change locations by tunneling, an instantaneous process; moreover, the Grotthuss mechanism suggests that these small particles can respond nearly instantaneously to the dynamic motions of proteins that occur on much longer time scales. A conformational change is accompanied by favorable or unfavorable variations in the free energy of the system, providing the impetus for solvent ↔ protein proton exchange. Thus, as thermal distributions of protein conformations interconvert, protonation states rapidly respond, as specific acidic and basic sites are exposed or protected. In the vacuum of the mass spectrometer, protons become immobilized in locations that are specific to the protein conformations from which they were incorporated. In this way, conformational states from solution are preserved upon electrospraying them into the gas phase. These ideas are consistent with the exquisite sensitivity of electrospray mass spectra to small changes of the local environment that alter protein structure in solution. We might remember this approximation for the protonation of proteins in solution with the colloquial expression-protons are fast and smart; proteins are slow and dumb.

Peptide-MHC I complex stability measured by nanoscale differential scanning fluorimetry reveals molecular mechanism of thermal denaturation

Recombinant major histocompatibility complex class I molecules are used in diagnostic and therapeutic approaches in cancer immunotherapy, with many studies exploring their binding to antigenic peptides. Current techniques for kinetic peptide binding studies are hampered by high sample consumption, low throughput, interference with protein stability, and/or high background signal. Here, we validate nanoscale differential scanning fluorimetry (nanoDSF), a method using the tryptophan fluorescence of class I molecules, for class I/peptide binding, and we use it to determine the molecular mechanism of the thermal denaturation of HLA-A*02:01.

Improving the Thermostability of Xylanase A from Bacillus subtilis by Combining Bioinformatics and Electrostatic Interactions Optimization

The rational improvement of the enzyme catalytic activity is one of the most significant challenges in biotechnology. Most conventional strategies used to engineer enzymes involve selecting mutations to increase their thermostability. Determining good criteria for choosing these substitutions continues to be a challenge. In this work, we combine bioinformatics, electrostatic analysis, and molecular dynamics to predict beneficial mutations that may improve the thermostability of XynA from Bacillus subtilis. First, the Tanford-Kirkwood surface accessibility method is used to characterize each ionizable residue contribution to the protein native state stability. Residues identified to be destabilizing were mutated with the corresponding residues determined by the consensus or ancestral sequences at the same locations. Five mutants (K99T/N151D, K99T, S31R, N151D, and K154A) were investigated and compared with 12 control mutants derived from experimental approaches from the literature. Molecular dynamics results show that the mutants exhibited folding temperatures in the order K99T > K99T/N151D > S31R > N151D > WT > K154A. The combined approaches employed provide an effective strategy for low-cost enzyme optimization needed for large-scale biotechnological and medical applications.

Transferrin receptor targeting by de novo sheet extension

The de novo design of polar protein-protein interactions is challenging because of the thermodynamic cost of stripping water away from the polar groups. Here, we describe a general approach for designing proteins which complement exposed polar backbone groups at the edge of beta sheets with geometrically matched beta strands. We used this approach to computationally design small proteins that bind to an exposed beta sheet on the human transferrin receptor (hTfR), which shuttles interacting proteins across the blood-brain barrier (BBB), opening up avenues for drug delivery into the brain. We describe a design which binds hTfR with a 20 nM Kd, is hyperstable, and crosses an in vitro microfluidic organ-on-a-chip model of the human BBB. Our design approach provides a general strategy for creating binders to protein targets with exposed surface beta edge strands.

A nanobody toolbox targeting dimeric coiled-coil modules for functionalization of designed protein origami structures

Coiled-coil (CC) dimers are widely used in protein design because of their modularity and well-understood sequence-structure relationship. In CC protein origami design, a polypeptide chain is assembled from a defined sequence of CC building segments that determine the self-assembly of protein cages into polyhedral shapes, such as the tetrahedron, triangular prism, or four-sided pyramid. However, a targeted functionalization of the CC modules could significantly expand the versatility of protein origami scaffolds. Here, we describe a panel of single-chain camelid antibodies (nanobodies) directed against different CC modules of a de novo designed protein origami tetrahedron. We show that these nanobodies are able to recognize the same CC modules in different polyhedral contexts, such as isolated CC dimers, tetrahedra, triangular prisms, or trigonal bipyramids, thereby extending the ability to functionalize polyhedra with nanobodies in a desired stoichiometry. Crystal structures of five nanobody-CC complexes in combination with small-angle X-ray scattering show binding interactions between nanobodies and CC dimers forming the edges of a tetrahedron with the nanobody entering the tetrahedral cavity. Furthermore, we identified a pair of allosteric nanobodies in which the binding to the distant epitopes on the antiparallel homodimeric APH CC is coupled via a strong positive cooperativity. A toolbox of well-characterized nanobodies specific for CC modules provides a unique tool to target defined sites in the designed protein structures, thus opening numerous opportunities for the functionalization of CC protein origami polyhedra or CC-based bionanomaterials.

Self-reconfigurable high-weight-per-volume-gelatin films for all-solution-processed on-skin electronics with ultra-conformal contact

Although conventional skin-attachable electronics exhibit good functionalities, their direct attachment (without any adhesive) to human skin with sufficient conformal contact is challenging. Herein, all-solution-processed on-skin electronics based on self-reconfigurable high-weight-per- volume-gelatin (HWVG) film constructed using an effective, biocompatible water absorption-evaporation technique are demonstrated. Completely conformal contact of self-reconfigurable HWVG films is realized by rapidly inducing anisotropic swelling in the perpendicular direction and covering any curvature on the skin without spatial gap or void after shrinking. A sufficiently thin HWVG film (~2 um) exhibited higher adhesion owing to van der Waals force and the carboxylic acid and amine groups in HWVG film form cross-linkages through intermolecular bonds with human skin. Self-reconfigurable HWVG films with high biocompatibility are optimized to afford a superior efficiency of 87.83 % at a concentration of 20 % (w/v) and a storage modulus of 1822 MPa at 36.5 °C. Furthermore, functional nanoelectrodes consisting of self-reconfigurable silver nanowires/HWVG films for high-performance on-skin sensors allowing the detection of sensitive motion and electrophysiological signals, as well as an armband-type sensor system incorporated with a smartphone for health-care monitoring are demonstrated. Outstanding performances, including stability, reliability, flexibility, re-usability, biocompatibility, and permeability of on-skin electronics based on HWVG films can open-up a prospective route to realizing breathable human-machine interfaces based on biocompatible materials and processes.

Anion-Specific Effects on the Alkaline State of Cytochrome c

Specific effects of anions on the structure, thermal stability, and peroxidase activity of native (state III) and alkaline (state IV) cytochrome c (cyt c) have been studied by the UV-VIS absorbance spectroscopy, intrinsic tryptophan fluorescence, and circular dichroism. Thermal and isothermal denaturation monitored by the tryptophan fluorescence and circular dichroism, respectively, implied lower stability of cyt c state IV in comparison with the state III. The pK_a value of alkaline isomerization of cyt c depended on the present salts, i.e., kosmotropic anions increased and chaotropic anions decreased pK_a (Hofmeister effect on protein stability). The peroxidase activity of cyt c in the state III, measured by oxidation of guaiacol, showed clear dependence on the salt position in the Hofmeister series, while cyt c in the alkaline state lacked the peroxidase activity regardless of the type of anions present in the solution. The alkaline isomerization of cyt c in the presence of 8 M urea, measured by Trp59 fluorescence, implied an existence of a high-affinity non-native ligand for the heme iron even in a partially denatured protein conformation. The conformation of the cyt c alkaline state in 8 M urea was considerably modulated by the specific effect of anions. Based on the Trp59 fluorescence quenching upon titration to alkaline pH in 8 M urea and molecular dynamics simulation, we hypothesize that the Lys79 conformer is most likely the predominant alkaline conformer of cyt c. The high affinity of the sixth ligand for the heme iron is likely a reason of the lack of peroxidase activity of cyt c in the alkaline state.

Halogen Atoms in the Protein-Ligand System. Structural and Thermodynamic Studies of the Binding of Bromobenzotriazoles by the Catalytic Subunit of Human Protein Kinase CK2

Binding of a family of brominated benzotriazoles to the catalytic subunit of human protein kinase CK2 (hCK2α) was used as a model system to assess the contribution of halogen bonding to protein–ligand interaction. CK2 is a constitutively active pleiotropic serine/threonine protein kinase that belongs to the CMGC group of eukaryotic protein kinases (EPKs). Due to the addiction of some cancer cells, CK2 is an attractive and well-characterized drug target. Halogenated benzotriazoles act as ATP-competitive inhibitors with unexpectedly good selectivity for CK2 over other EPKs. We have characterized the interaction of bromobenzotriazoles with hCK2α by X-ray crystallography, low-volume differential scanning fluorimetry, and isothermal titration calorimetry. Properties of free ligands in solution were additionally characterized by volumetric and RT-HPLC measurements. Thermodynamic data indicate that the affinity increases with bromo substitution, with greater contributions from 5- and 6-substituents than 4- and 7-substituents. Except for 4,7-disubstituted compounds, the bromobenzotriazoles adopt a canonical pose with the triazole close to lysine 68, which precludes halogen bonding. More highly substituted benzotriazoles adopt many additional noncanonical poses, presumably driven by a large hydrophobic contribution to binding. Some noncanonical ligand orientations allow the formation of halogen bonds with the hinge region. Consistent with a predominantly hydrophobic interaction, the isobaric heat capacity decreases upon ligand binding, the more so the higher the substitution.

Interdomain linkers tailor the stability of immunoglobulin repeats in polyproteins

Linkers in polyproteins are considered as mere spacers between two adjacent domains. However, a series of studies using single-molecule force spectroscopy have recently reported distinct thermodynamic stability of I27 in polyproteins with varying linkers and indicated the vital role of linkers in domain stability. A flexible glycine rich linker (-(GGG)_n, n ≥ 3) featured unfolding at lower forces than the regularly used arg-ser (RS) based linker. Interdomain interactions among I27 domains in Gly-rich linkers were suggested to lead to reduced domain stability. However, the negative impact of inter domain interactions on domain stability is thermodynamically counter-intuitive and demanded thorough investigations. Here, using an array of ensemble equilibrium experiments and in-silico measurements with I27 singlet and doublets with two aforementioned linkers, we delineate that the inter-domain interactions in fact raise the stability of the polyprotein with RS linker. More surprisingly, a highly flexible Gly-rich linker has no interference on the stability of polyprotein. Overall, we conclude that flexible linkers are preferred in a polyprotein for maintaining domain's independence.

The Folding Pathway of 6aJL2

One-third of the reported cases of light chain amyloidosis are related to the germ line λ6 family; remarkably, healthy individuals express this type of protein in just 2% of the peripheral blood and bone marrow B-cells. The appearance of the disease has been related to the inherent properties of this protein family. A recombinant representative model for λ6 proteins called 6aJL2 containing the amino acid sequence encoded by the 6a and JL2 germ line genes was previously designed and synthesized to study the properties of this family. Previous work on 6aJL2 suggested a simple two-state folding model at 25 °C; no intermediate could be identified either by kinetics or by fluorescence and circular dichroism equilibrium studies, although the presence of an intermediate that is populated at ∼2.4 M urea was suggested by size exclusion chromatography. In this study we employed classic equilibrium and kinetic experiments and analysis to elucidate the detailed folding mechanism of this protein. We identify species that are kinetically accessible and/or are populated at equilibrium. We describe the presence of intermediate and native-like species and propose a five-species folding mechanism at 25 °C at short incubation times, similar to and consistent with those observed in other proteins of this fold. The formation of intermediates in the mechanism of 6aJL2 is faster than that proposed for a V_κ light chain, which could be an important distinction in the amyloidogenic potential of both germ lines.

Linking thermodynamics and measurements of protein stability

We review the background, theory and general equations for the analysis of equilibrium protein unfolding experiments, focusing on denaturant and heat-induced unfolding. The primary focus is on the thermodynamics of reversible folding/unfolding transitions and the experimental methods that are available for extracting thermodynamic parameters. We highlight the importance of modelling both how the folding equilibrium depends on a perturbing variable such as temperature or denaturant concentration, and the importance of modelling the baselines in the experimental observables.

Suppression of Electrostatic Mediated Antibody Liquid-Liquid Phase Separation by Charged and Noncharged Preferentially Excluded Excipients

Isotonic concentrations of inert cosolutes or excipients are routinely used in protein therapeutic formulations to minimize physical instabilities including aggregation, particulation, and precipitation that are often manifested during drug substance/product manufacture and long-term storage. Despite their prevalent use within the biopharmaceutical industry, a more detailed understanding for how excipients modulate the specific protein-protein interactions responsible for these instabilities is still needed so that informed formulation decisions can be made at the earliest stages of development when protein supply and time are limited. In the present report, subisotonic concentrations of the five common formulation excipients, sucrose, proline, sorbitol, glycerol, arginine hydrochloride, and the denaturant urea, were studied for their effect on the room temperature liquid-liquid phase separation of a model monoclonal antibody (mAb-B). Although each excipient lowered the onset temperatures of mAb-B liquid-liquid phase separation to different extents, all six were found to be preferentially excluded from the native state monomer by vapor pressure osmometry, and no apparent correlations to the excipient dependence of mAb-B melting temperatures were observed. These results and those of the effects of solution pH, addition of salt, and impact of a small number of charge mutations were most consistent with a mechanism of local excipient accumulation, to an extent dependent on their type, with the specific residues that mediate mAb-B electrostatic protein-protein interactions. These findings suggest that selection of excipients on the basis of their interaction with the solvent exposed residues of the native state may at times be a more effective strategy for limiting protein-protein interactions at pharmaceutically relevant storage conditions than choosing those that are excluded from the residues of the native state interior.

Urea titration of a lipase from Pseudomonas sp. reveals four different conformational states, with a stable partially folded state explaining its high aggregation propensity

The conversion of soluble proteins into amyloid fibrils has importance in protein chemistry, biology, biotechnology and medicine. A novel lipase from Pseudomonas sp. was previously shown to have an extremely high aggregation propensity. It was therefore herein studied to elucidate the physicochemical and structural determinants of this extreme behaviour. Amyloid-like structures were found to form in samples up to 2.5-3.0 M using Thioflavin T fluorescence and Congo red binding assays. However, dynamic light scattering (DLS), static light scattering and turbidimetry revealed the existence of aggregates up to 4.0 M urea, without amyloid-like structure. Two monomeric conformational states were detected with intrinsic fluorescence, 8-anilinonaphthalene-1-sulfonate (ANS) binding and circular dichroism. These were further characterized in 7.5 M and 4.5 M urea using enzymatic activity measurements, tryptophan fluorescence quenching, DLS and nuclear magnetic resonance (NMR) and were found to consist of a largely disordered and a partially folded state, respectively, with the latter appearing stable, cooperative, fairly compact, non-active, α-helical, with largely buried hydrophobic residues. The persistence of a stable structure up to high concentrations of urea, in the absence of sequence characteristics typical of a high intrinsic aggregation propensity, explains the high tendency of this enzyme to form amyloid-like structures.

Double Mutant Cycles as a Tool to Address Folding, Binding, and Allostery

Quantitative measurement of intramolecular and intermolecular interactions in protein structure is an elusive task, not easy to address experimentally. The phenomenon denoted 'energetic coupling' describes short- and long-range interactions between two residues in a protein system. A powerful method to identify and quantitatively characterize long-range interactions and allosteric networks in proteins or protein-ligand complexes is called double-mutant cycles analysis. In this review we describe the thermodynamic principles and basic equations that underlie the double mutant cycle methodology, its fields of application and latest employments, and caveats and pitfalls that the experimentalists must consider. In particular, we show how double mutant cycles can be a powerful tool to investigate allosteric mechanisms in protein binding reactions as well as elusive states in protein folding pathways.

Exploring the Evolutionary History of Kinetic Stability in the α-Lytic Protease Family

In addition to encoding the tertiary fold and stability, the primary sequence of a protein encodes the folding trajectory and kinetic barriers that determine the speed of folding. How these kinetic barriers are encoded is not well understood. Here, we use evolutionary sequence variation in the α-lytic protease (αLP) protein family to probe the relationship between sequence and energy landscape. αLP has an unusual energy landscape: the native state of αLP is not the most thermodynamically favored conformation and, instead, remains folded due to a large kinetic barrier preventing unfolding. To fold, αLP utilizes an N-terminal pro region similar in size to the protease itself that functions as a folding catalyst. Once folded, the pro region is removed, and the native state does not unfold on a biologically relevant time scale. Without the pro region, αLP folds on the order of millennia. A phylogenetic search uncovers αLP homologs with a wide range of pro region sizes, including some with no pro region at all. In the resulting phylogenetic tree, these homologs cluster by pro region size. By studying homologs naturally lacking a pro region, we demonstrate they can be thermodynamically stable, fold much faster than αLP, yet retain the same fold as αLP. Key amino acids thought to contribute to αLP's extreme kinetic stability are lost in these homologs, supporting their role in kinetic stability. This study highlights how the entire energy landscape plays an important role in determining the evolutionary pressures on the protein sequence.

Conformational Plasticity and DNA-Binding Specificity of the Eukaryotic Transcription Factor Pax5

The eukaryotic transcription factor Pax5 has a DNA-binding Paired domain composed of two independent helical bundle subdomains joined by a flexible linker. Previously, we showed distinct biophysical properties of the N-terminal (NTD) and C-terminal (CTD) subdomains, with implications for how these two regions cooperate to distinguish nonspecific and cognate DNA sites [Perez-Borrajero, C., et al. (2016) J. Mol. Biol. 428, 2372-2391]. In this study, we combined experimental methods and molecular dynamics (MD) simulations to dissect the mechanisms underlying the functional differences between the Pax5 subdomains. Both subdomains showed a similar dependence of DNA-binding affinity on ionic strength. However, due to a greater contribution of non-ionic interactions, the NTD bound its cognate DNA half-site with an affinity approximately 10-fold higher than that of the CTD with its half-site. These interactions involve base-mediated contacts as evidenced by nuclear magnetic resonance spectroscopy-monitored chemical shift perturbations. Isothermal titration calorimetry revealed that favorable enthalpic and compensating unfavorable entropic changes were substantially larger for DNA binding by the NTD than by the CTD. Complementary MD simulations indicated that the DNA recognition helix H3 of the NTD is particularly flexible in the absence of DNA and undergoes the largest changes in conformational dynamics upon binding. Overall, these data suggest that the differences observed for the subdomains of Pax5 are due to the coupling of DNA binding with dampening of motions in the NTD required for specific base contacts. Thus, the conformational plasticity of the Pax5 Paired domain underpins the differing roles of its subdomains in association with nonspecific versus cognate DNA sites.

Unveiling the Folding Mechanism of PDZ Domains

Understanding the mechanism of folding of single domain proteins demands a complete characterization of their equilibrium and kinetic properties. By using a well-studied class of protein domain, the PDZ domain, here we exemplify the typical procedure to address this problem.

A lowly populated, transient β-sheet structure in monomeric Aβ1-42 identified by multinuclear NMR of chemical denaturation

Chemical denaturation is a well-established approach for probing the equilibrium between folded and unfolded states of proteins. We demonstrate applicability of this method to the detection of a small population of a transiently folded structural element in a system that is often considered to be intrinsically fully disordered. The ¹H^N, ¹⁵N, ¹³C^α, and ¹³C' chemical shifts of Aβ^1-40 and Aβ^1-42 peptides and their M35-oxidized variants were monitored as a function of urea concentration and compared to analogous urea titrations of synthetic pentapeptides of homologous sequence. Fitting of the chemical shift titrations yields a 10 ± 1% population for a structured element at the C-terminus of Aβ^1-42 that folds with a cooperativity of m = 0.06 kcal/mol·M. The fit also yields the chemical shifts of the folded state and, using a database search, for Aβ^1-42 these shifts identified an antiparallel intramolecular β-sheet for residues I32-A42, linked by a type I' β-turn at G37 and G38. The structure is destabilized by oxidation of M35. Paramagnetic relaxation rates and two previously reported weak, medium-range NOE interactions are consistent with this transient β-sheet. Introduction of the requisite A42C mutation and tagging with MTSL resulted in a small stabilization of this β-sheet. Chemical shift analysis suggests a C-terminal β-sheet may be present in Aβ^1-40 too, but the turn type at G37 is not type I'. The approach to derive Transient Structure from chemical Denaturation by NMR (TSD-NMR), demonstrated here for Aβ peptides, provides a sensitive tool for identifying the presence of lowly populated, transiently ordered elements in proteins that are considered to be intrinsically disordered, and permits extraction of structural data for such elements.

Leishmanial aspartyl-tRNA synthetase: Biochemical, biophysical and structural insights

Aminoacyl tRNA synthetases (aaRSs) are integral components of protein biosynthesis along with several non-canonical cellular processes. Inhibition studies of aaRSs presented these enzymes as promising drug targets in many pathogens, however aspartyl tRNA synthetase has not been studied in trypanosomatids despite its essentiality. Hence, full-length ORF of Leishmania donovani aspartyl tRNA synthetase (LdaspRS) was cloned and purified to homogeneity followed by molecular mass determination. The aminoacylation assay established that the purified protein performs its function optimally at physiological pH and temperature. The kinetic parameters of LdaspRS revealed the affinity of l-aspartate towards the enzyme to be very much lower than the cofactor. Our study also highlights the moonlighting function of LdaspRS to stimulate the pro-inflammatory cytokines and nitric oxide generation by host macrophage. Furthermore, CD and intrinsic tryptophan fluorescence measurements showed the changes in structural conformation at varying pH, denaturants and ligands. The modelled LdaspRS structure presented all the specific characteristics of class II aaRSs, while in silico study suggested binding of pyrimidine-derived inhibitors in its cofactor binding site with high affinity followed by validation using MD simulation. Altogether, this study could provide a platform for exploring LdaspRS to develop potential therapeutics against leishmaniasis.

Zinc shapes the folding landscape of p53 and establishes a pathway for reactivating structurally diverse cancer mutants

Missense mutations in the p53 DNA-binding domain (DBD) contribute to half of new cancer cases annually. Here we present a thermodynamic model that quantifies and links the major pathways by which mutations inactivate p53. We find that DBD possesses two unusual properties—one of the highest zinc affinities of any eukaryotic protein and extreme instability in the absence of zinc—which are predicted to poise p53 on the cusp of folding/unfolding in the cell, with a major determinant being available zinc concentration. We analyze the 20 most common tumorigenic p53 mutations and find that 80% impair zinc affinity, thermodynamic stability, or both. Biophysical, cell-based, and murine xenograft experiments demonstrate that a synthetic zinc metallochaperone rescues not only mutations that decrease zinc affinity, but also mutations that destabilize DBD without impairing zinc binding. The results suggest that zinc metallochaperones have the capability to treat 120,500 patients annually in the U.S.

Solute-Solvent Interactions in Modern Physical Organic Chemistry: Supramolecular Polymers as a Muse

Interactions between solvents and solutes are a cornerstone of physical organic chemistry and have been the subject of investigations over the last century. In recent years, a renewed interest in fundamental aspects of solute-solvent interactions has been sparked in the field of supramolecular chemistry in general and that of supramolecular polymers in particular. Although solvent effects in supramolecular chemistry have been recognized for a long time, the unique opportunities that supramolecular polymers offer to gain insight into solute-solvent interactions have become clear relatively recently. The multiple interactions that hold the supramolecular polymeric structure together are similar in strength to those between solute and solvent. The cooperativity found in ordered supramolecular polymers leads to the possibility of amplifying these solute-solvent effects and will shed light on extremely subtle solvation phenomena. As a result, many exciting effects of solute-solvent interactions in modern physical organic chemistry can be studied using supramolecular polymers. Our aim is to put the recent progress into a historical context and provide avenues toward a more comprehensive understanding of solvents in multicomponent supramolecular systems.

Potassium Glutamate and Glycine Betaine Induce Self-Assembly of the PCNA and β-Sliding Clamps

Sliding clamps are oligomeric ring-shaped proteins that increase the efficiency of DNA replication. The stability of the Escherichia coli β-clamp, a homodimer, is particularly remarkable. The dissociation equilibrium constant of the β-clamp is of the order of 10 pM in buffers of moderate ionic strength. Coulombic electrostatic interactions have been shown to contribute to this remarkable stability. Increasing NaCl concentration in the assay buffer results in decreased dimer stability and faster subunit dissociation kinetics in a way consistent with simple charge-screening models. Here, we examine non-Coulombic ionic effects on the oligomerization properties of sliding clamps. We determined relative diffusion coefficients of two sliding clamps using fluorescence correlation spectroscopy. Replacing NaCl by KGlu, the primary cytoplasmic salt in E. coli, results in a decrease of the diffusion coefficient of these proteins consistent with the formation of protein assemblies. The UV-vis spectrum of the β-clamp labeled with tetramethylrhodamine shows the characteristic absorption band of dimers of rhodamine when KGlu is present in the buffer. This suggests that KGlu induces the formation of assemblies that involve two or more rings stacked face-to-face. Results can be quantitatively explained on the basis of unfavorable interactions between KGlu and the functional groups on the protein surface, which drive biomolecular processes that bury exposed surface. Similar results were obtained with the Saccharomyces cerevisiae PCNA sliding clamp, suggesting that KGlu effects are not specific to the β-clamp. Clamp association is also promoted by glycine betaine, a zwitterionic compound that accumulates intracellularly when E. coli is exposed to high concentrations of extracellular solute. Possible biological implications are discussed.

Benefits and constrains of covalency: the role of loop length in protein stability and ligand binding

Protein folding is governed by non-covalent interactions under the benefits and constraints of the covalent linkage of the backbone chain. In the current work we investigate the influence of loop length variation on the free energies of folding and ligand binding in a small globular single-domain protein containing two EF-hand subdomains—calbindin D_9k. We introduce a linker extension between the subdomains and vary its length between 1 to 16 glycine residues. We find a close to linear relationship between the linker length and the free energy of folding of the Ca²⁺-free protein. In contrast, the linker length has only a marginal effect on the Ca²⁺ affinity and cooperativity. The variant with a single-glycine extension displays slightly increased Ca²⁺ affinity, suggesting that the slightly extended linker allows optimized packing of the Ca²⁺-bound state. For the extreme case of disconnected subdomains, Ca²⁺ binding becomes coupled to folding and assembly. Still, a high affinity between the EF-hands causes the non-covalent pair to retain a relatively high apparent Ca²⁺ affinity. Our results imply that loop length variation could be an evolutionary option for modulating properties such as protein stability and turnover without compromising the energetics of the specific function of the protein.

Cryo vs Thermo: Duality of Ethylene Glycol on the Stability of Proteins

Osmolytes are known to stabilize proteins under stress conditions. Thermal denaturation studies on globular proteins (β-lactoglobulin, cytochrome c, myoglobin, α-chymotrypsin) in the presence of ethylene glycol (EG), a polyol class of osmolyte, demonstrate a unique property of EG. EG stabilizes proteins against cold denaturation and destabilizes them during heat-induced denaturation. Further, chemical denaturation experiments performed at room temperature and at a sub-zero temperature (-10 °C) show that EG stabilizes the proteins at subzero temperature but destabilizes them at room temperature. The experiments carried out in the presence of glycerol, however, showed that glycerol stabilizes proteins against all of the denaturing conditions. This differential effect has not been reported for any other polyol class of osmolyte and might be specific to EG. Moreover, molecular dynamics simulations of all of the four proteins were carried out at three different temperatures, 240, 300, and 340 K, in the absence and presence of EG (20 and 40%). The results suggest that EG preferably accumulates around the hydrophobic residues and reduces the hydrophobic hydration of the proteins at a low temperature leading to stabilization of the proteins. At 340 K, the preferential hydration of the proteins is significantly reduced and the preferential binding of EG destabilizes the proteins like common denaturants.

HDX and In Silico Docking Reveal that Excipients Stabilize G-CSF via a Combination of Preferential Exclusion and Specific Hotspot Interactions

Assuring the stability of therapeutic proteins is a major challenge in the biopharmaceutical industry, and a better molecular understanding of the mechanisms through which formulations influence their stability is an ongoing priority. While the preferential exclusion effects of excipients are well known, the additional presence and impact of specific protein-excipient interactions have proven to be more elusive to identify and characterize. We have taken a combined approach of in silico molecular docking and hydrogen deuterium exchange-mass spectrometry (HDX-MS) to characterize the interactions between granulocyte colony-stimulating factor (G-CSF), and some common excipients. These interactions were related to their influence on the thermal-melting temperatures (T_m) for the nonreversible unfolding of G-CSF in liquid formulations. The residue-level interaction sites predicted in silico correlated well with those identified experimentally and highlighted the potential impact of specific excipient interactions on the T_m of G-CSF.

Origin of heat capacity increment in DNA folding: The hydration effect

BACKGROUND

Understanding DNA folding thermodynamics is crucial for prediction of DNA thermal stability. It is now well established that DNA folding is accompanied by a decrease of the heat capacity ∆c_{p, F}, however its molecular origin is not understood. In analogy to protein folding it has been assumed that this is due to dehydration of DNA constituents, however no evidence exists to support this conclusion.

METHODS

Here we analyze partial molar heat capacity of nucleic bases and nucleosides in aqueous solutions obtained from calorimetric experiments and calculate the hydration heat capacity contribution ∆c_p^hyd.

RESULTS

We present hydration heat capacity contributions of DNA constituents and show that they correlate with the solvent accessible surface area. The average contribution for nucleic base dehydration is +0.56 J mol^-1 K^-1 Å^-2 and can be used to estimate the ∆c_{p, F} contribution for DNA folding.

CONCLUSIONS

We show that dehydration is one of the major sources contributing to the observed ∆c_{p, F} increment in DNA folding. Other possible sources contributing to the overall ∆c_{p, F} should be significant but appear to compensate each other to high degree. The calculated ∆c_p^hyd for duplexes and noncanonical DNA structures agree excellently with the overall experimental ∆c_{p, F} values. By contrast, empirical parametrizations developed for proteins result in poor ∆c_p^hyd predictions and should not be applied to DNA folding.

GENERAL SIGNIFICANCE

Heat capacity is one of the main thermodynamic quantities that strongly affects thermal stability of macromolecules. At the molecular level the heat capacity in DNA folding stems from removal of water from nucleobases.

Ensemble-Based Thermodynamics of the Fuzzy Binding between Intrinsically Disordered Proteins and Small-Molecule Ligands

In contrast to the "lock-and-key" model underlying the long-term success of structural biology and rational drug design, intrinsically disordered proteins (IDPs) exist in an ensemble of highly heterogeneous conformations even after binding with small-molecule ligands. It remains controversial how to characterize the thermodynamics of such fuzzy interactions. Here, we derive an ensemble-based thermodynamic framework to analyze the apparent affinity between IDPs and ligands. It is shown that the apparent affinity is related to the interaction free energy between the individual conformation and ligand in a way similar to Jarzynski's equality in nonequilibrium statistics. The oncoprotein c-Myc is adopted as an example to demonstrate the related properties, for example, the distribution of conformation-ligand interaction free energy, the entropic contribution from the ensemble, the conformation shift under ligand binding, and how to control the error under a limited number of sampled conformations.

Thermodynamic and Mechanistic Insights into Coupled Binding and Unwinding of Collagen by Matrix Metalloproteinase 1

Local unwinding of the collagen triple helix is a necessary step for initiating the collagen degradation cascade in extracellular matrices. A few matrix metalloproteinases (MMPs) are known to support this key process, but its energetic aspects remain unknown. Here, we captured the thermodynamics of the triple helix unwinding by monitoring interactions between a collagen peptide and MMP-1(E200A) - an active-site mutant of an archetypal vertebrate collagenase - at increasing temperatures, using isothermal titration calorimetry (ITC). Coupled binding and unwinding manifests as a curved relationship between the total enthalpy change and temperature of the reaction, producing increasingly negative heat capacity change (ΔΔC_p ≈ -36.3 kcal/molK²). A specially designed solid-phase binding and cleavage assay (SPBCA) reported strain in the catalytically relevant unwound state, suggesting that this state is distinct from the horizon of sampled conformations of the collagenase-susceptible site. MMP-1 appears to blend selected fit with induced fit mechanisms to catalyse collagen unwinding prior to cleavage of individual collagen chains.