Protein folding intermediates and pathways studied by hydrogen exchange

Inside of the burst-phase intermediate of a protein folding. Hydration of hydrophobic groups

The thermal effect of the formation of the "burst-phase" folding intermediate has been studied using a titration calorimeter. It is shown that, unlike the total thermal effect of native structure formation, it can be both positive and negative depending on the temperature. The reasons for this paradoxical behavior are analyzed. A conclusion is drawn about the leading role of dehydration of non-polar groups in the first stage of folding.

The molecular determinants of classical pathway complement inhibition by OspEF-related proteins of Borrelia burgdorferi

The complement system serves as the first line of defense against invading pathogens by promoting opsonophagocytosis and bacteriolysis. Antibody-dependent activation of complement occurs through the classical pathway and relies on the activity of initiating complement proteases of the C1 complex, C1r and C1s. The causative agent of Lyme disease, Borrelia burgdorferi, expresses two paralogous outer surface lipoproteins of the OspEF-related protein family, ElpB and ElpQ, that act as specific inhibitors of classical pathway activation. We have previously shown that ElpB and ElpQ bind directly to C1r and C1s with high affinity and specifically inhibit C2 and C4 cleavage by C1s. To further understand how these novel protease inhibitors function, we carried out a series of hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments using ElpQ and full-length activated C1s as a model of Elp-protease interaction. Comparison of HDX-MS profiles between unbound ElpQ and the ElpQ/C1s complex revealed a putative C1s-binding site on ElpQ. HDX-MS-guided, site-directed ElpQ mutants were generated and tested for direct binding to C1r and C1s using surface plasmon resonance. Several residues within the C-terminal region of ElpQ were identified as important for protease binding, including a single conserved tyrosine residue that was required for ElpQ- and ElpB-mediated complement inhibition. Collectively, our study identifies key molecular determinants for classical pathway protease recognition by Elp proteins. This investigation improves our understanding of the unique complement inhibitory mechanism employed by Elp proteins which serve as part of a sophisticated complement evasion system present in Lyme disease spirochetes.

Protein dynamics govern the oxyferrous state lifetime of an artificial oxygen transport protein

It has long been known that the alteration of protein side chains that occlude or expose the heme cofactor to water can greatly affect the stability of the oxyferrous heme state. Here, we demonstrate that the rate of dynamically driven water penetration into the core of an artificial oxygen transport protein also correlates with oxyferrous state lifetime by reducing global dynamics, without altering the structure of the active site, via the simple linking of the two monomers in a homodimeric artificial oxygen transport protein using a glycine-rich loop. The tethering of these two helices does not significantly affect the active site structure, pentacoordinate heme-binding affinity, reduction potential, or gaseous ligand affinity. It does, however, significantly reduce the hydration of the protein core, as demonstrated by resonance Raman spectroscopy, backbone amide hydrogen exchange, and pKa shifts in buried histidine side chains. This further destabilizes the charge-buried entatic state and nearly triples the oxyferrous state lifetime. These data are the first direct evidence that dynamically driven water penetration is a rate-limiting step in the oxidation of these complexes. It furthermore demonstrates that structural rigidity that limits water penetration is a critical design feature in metalloenzyme construction and provides an explanation for both the failures and successes of earlier attempts to create oxygen-binding proteins.

Protein folding rate evolution upon mutations

Despite the spectacular success of cutting-edge protein fold prediction methods, many critical questions remain unanswered, including why proteins can reach their native state in a biologically reasonable time. A satisfactory answer to this simple question could shed light on the slowest folding rate of proteins as well as how mutations-amino-acid substitutions and/or post-translational modifications-might affect it. Preliminary results indicate that (i) Anfinsen's dogma validity ensures that proteins reach their native state on a reasonable timescale regardless of their sequence or length, and (ii) it is feasible to determine the evolution of protein folding rates without accounting for epistasis effects or the mutational trajectories between the starting and target sequences. These results have direct implications for evolutionary biology because they lay the groundwork for a better understanding of why, and to what extent, mutations-a crucial element of evolution and a factor influencing it-affect protein evolvability. Furthermore, they may spur significant progress in our efforts to solve crucial structural biology problems, such as how a sequence encodes its folding.

Deep mining of the protein energy landscape

For over half a century, it has been known that protein molecules naturally undergo extensive structural fluctuations, and that these internal motions are intimately related to their functional properties. The energy landscape view has provided a powerful framework for describing the various physical states that proteins visit during their lifetimes. This Perspective focuses on the commonly neglected and often disparaged axis of the protein energy landscape: entropy. Initially seen largely as a barrier to functionally relevant states of protein molecules, it has recently become clear that proteins retain considerable conformational entropy in the "native" state, and that this entropy can and often does contribute significantly to the free energy of fundamental protein properties, processes, and functions. NMR spectroscopy, molecular dynamics simulations, and emerging crystallographic views have matured in parallel to illuminate dynamic disorder of the "ground state" of proteins and their importance in not only transiting between biologically interesting structures but also greatly influencing their stability, cooperativity, and contribution to critical properties such as allostery.

Elucidating the mechanisms underlying protein conformational switching using NMR spectroscopy

How proteins switch between various ligand-free and ligand-bound structures has been a key biophysical question ever since the postulation of the Monod-Wyman-Changeux and Koshland-Nemethy-Filmer models over six decades ago. The ability of NMR spectroscopy to provide structural and kinetic information on biomolecular conformational exchange places it in a unique position as an analytical tool to interrogate the mechanisms of biological processes such as protein folding and biomolecular complex formation. In addition, recent methodological developments in the areas of saturation transfer and relaxation dispersion have expanded the scope of NMR for probing the mechanics of transitions in systems where one or more states constituting the exchange process are sparsely populated and 'invisible' in NMR spectra. In this review, we highlight some of the strategies available from NMR spectroscopy for examining the nature of multi-site conformational exchange, using five case studies that have employed NMR, either in isolation, or in conjunction with other biophysical tools.

Heterologous expression of novel SUMO proteases from Schizosaccharomyces pombe in E. coli: Catalytic domain identification and optimization of product yields

Small ubiquitin-related modifier (SUMO) proteins are efficiently used to target the soluble expression of various difficult-to-express proteins in E. coli. However, its utilization in large scale protein production is restricted by the higher cost of Ulp, which is required to cleave SUMO fusion tag from protein-of-interest to generate an authentic N-terminus. This study identified and characterized two novel SUMO proteases i.e., Ulp1 and Ulp2 from Schizosaccharomyces pombe. Codon-optimized gene sequences were cloned and expressed in E. coli. The sequence and structure of SpUlp1 and SpUlp2 catalytic domains were deduced using bioinformatics tools. Protein-protein interaction studies predicted the higher affinity of SpUlp1 towards SUMO compared to its counterpart from Saccharomyces cerevisiae (ScUlp1). The catalytic domain of SpUlp1 was purified using Ni-NTA chromatography with 83.33% recovery yield. Moreover, In vitro activity data further confirmed the fast-acting nature of SpUlp1 catalytic domain, where a 90% cleavage of fusion proteins was obtained within 1 h of incubation, indicating novelty and commercial relevance of S. pombe Ulp1. Biophysical characterization showed 8.8% α-helices, 36.7% β-sheets in SpUlp1SD. From thermal CD and fluorescence data, SpUlp1SD T_m was found to be 45 °C. Further, bioprocess optimization using fed-batch cultivation resulted in 3.5 g/L of SpUlp1SD production with Y_P/X of 77.26 mg/g DCW and volumetric productivity of 205.88 mg/L/h.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

Native State Hydrogen Exchange-Mass Spectrometry Methods to Probe Protein Folding and Unfolding

Native state hydrogen exchange (HX) methods provide high-resolution structural data on the rare and transient opening motions in proteins under native conditions. Mass spectrometry-based HX methods (HX-MS) have gained popularity because of their ability to delineate population distributions, which allow a direct determination of the mechanism of inter conversion of the partially folded states under native conditions. Various technological advancements have provided further impetus to the development of HX-MS-based experiments to study protein folding. Classical HX-MS studies use proteolytic digestion to produce fragments of the protein subsequent to HX in solution, in order to obtain structural data. New chemical fragmentation methods, which achieve the same result as proteolysis and cause minimal change to the HX pattern in the protein, provide an attractive alternative to proteolysis. Moreover, when used in conjunction with proteolysis, chemical fragmentation methods have significantly increased the structural resolution afforded by HX-MS studies, even bringing them at par with the single amino acid resolution observed in NMR-based measurements. Experiments based on one such chemical fragmentation method, electron transfer dissociation (ETD), are described in this chapter. The ETD HX-MS method is introduced using data from a protein which is inherently resistant to proteolytic digestion as example of how such an experiment can provide high-resolution structural data on the folding-unfolding transitions of the protein under native conditions.

Folding Intermediates, Heterogeneous Native Ensembles and Protein Function

Single domain proteins fold via diverse mechanisms emphasizing the intricate relationship between energetics and structure, which is a direct consequence of functional constraints and demands imposed at the level of sequence. On the other hand, elucidating the interplay between folding mechanisms and function is challenging in large proteins, given the inherent shortcomings in identifying metastable states experimentally and the sampling limitations associated with computational methods. Here, we show that free energy profiles and surfaces of large systems (>150 residues), as predicted by a statistical mechanical model, display a wide array of folding mechanisms with ubiquitous folding intermediates and heterogeneous native ensembles. Importantly, residues around the ligand binding or enzyme active site display a larger tendency to partially unfold and this manifests as intermediates or excited states along the folding coordinate in ligand binding domains, transcription repressors, and representative enzymes from all the six classes, including the SARS-CoV-2 receptor binding domain (RBD) of the spike protein and the protease M^pro. It thus appears that it is relatively easier to distill the imprints of function on the folding landscape of larger proteins as opposed to smaller systems. We discuss how an understanding of energetic-entropic features in ordered proteins can pinpoint specific avenues through which folding mechanisms, populations of partially structured states and function can be engineered.

Structural Characterization of the Cooperativity of Unfolding of a Heterodimeric Protein using Hydrogen Exchange-Mass Spectrometry

Little is known about how the sequence of structural changes in one chain of a heterodimeric protein is coupled to those in the other chain during protein folding and unfolding reactions, and whether individual secondary structural changes in the two chains occur in one or many coordinated steps. Here, the unfolding mechanism of a small heterodimeric protein, double chain monellin, has been characterized using hydrogen exchange-mass spectrometry. Transient structure opening, which enables HX, was found to be describable by a five state N ↔ I₁ ↔ I₂ ↔ I₃ ↔ U mechanism. Structural changes occur gradually in the first three steps, and cooperatively in the last step. β strands 2, 4 and 5, as well as the α-helix undergo transient unfolding during all three non-cooperative steps, while β1 and the two loops on both sides of the helix undergo transient unfolding during the first two steps. In the absence of GdnHCl, only β3 in chain A of the protein unfolds during the last cooperative step, while in the presence of 1 M GdnHCl, not only β3, but also β2 in chain B unfolds cooperatively. Hence, the extent of cooperative structural change and size of the cooperative unfolding unit increase when the protein is destabilized by denaturant. The naturally evolved two-chain variant of monellin folds and unfolds in a more cooperative manner than does a single chain variant created artificially, suggesting that increasing folding cooperativity, even at the cost of decreasing stability, may be a driving force in the evolution of proteins.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

Elucidating dynamic behavior of synthetic supramolecular polymers in water by hydrogen/deuterium exchange mass spectrometry

A comprehensive understanding of the structure, self-assembly mechanism, and dynamics of one-dimensional supramolecular polymers in water is essential for their application as biomaterials. Although a plethora of techniques are available to study the first two properties, there is a paucity in possibilities to study dynamic exchange of monomers between supramolecular polymers in solution. We recently introduced hydrogen/deuterium exchange mass spectrometry (HDX-MS) to characterize the dynamic nature of synthetic supramolecular polymers with only a minimal perturbation of the chemical structure. To further expand the application of this powerful technique some essential experimental aspects have been reaffirmed and the technique has been applied to a diverse library of assemblies. HDX-MS is widely applicable if there are exchangeable hydrogen atoms protected from direct contact with the solvent and if the monomer concentration is sufficiently high to ensure the presence of supramolecular polymers during dilution. In addition, we demonstrate that the kinetic behavior as probed by HDX-MS is influenced by the internal order within the supramolecular polymers and by the self-assembly mechanism.

2D NMR spectroscopy of refolding RNase Sa using polarization transfer from hyperpolarized water

Polarization transfer from hyperpolarized water through proton exchange is used to enhance the NMR signals of amide protons of the Ribonuclease Sa protein. Spectra of the refolding protein are measured within 6 s after dilution of the denaturant urea, at urea-dependent folding rates adjusted in the range of 0.3-0.8 s^-1. Peak patterns including a mixture of folded and unfolded protein at different ratios are observed. The changes in the observed signals indicate that each spectrum accesses a different point in the partial completion of the folding. A comparison to simulated 2D NMR spectra suggests a lower polarization transfer efficiency from water when the protein folds slowly, which may result from the molecular motions in the unfolded protein and the absence of long-range contacts. The ability to acquire 2D NMR spectra under different refolding conditions may open a new avenue for residue specific characterization of the folding process.

A Disordered Loop Mediates Heterogeneous Unfolding of an Ordered Protein by Altering the Native Ensemble

The high flexibility of long disordered or partially structured loops in folded proteins allows for entropic stabilization of native ensembles. Destabilization of such loops could alter the native ensemble or promote alternate conformations within the native ensemble if the ordered regions themselves are held together weakly. This is particularly true of downhill folding systems that exhibit weak unfolding cooperativity. Here, we combine experimental and computational methods to probe the response of the native ensemble of a helical, downhill folding domain PDD, which harbors an 11-residue partially structured loop, to perturbations. Statistical mechanical modeling points to continuous structural changes on both temperature and mutational perturbations driven by entropic stabilization of partially structured conformations within the native ensemble. Long time-scale simulations of the wild-type protein and two mutants showcase a remarkable conformational switching behavior wherein the parallel helices in the wild-type protein sample an antiparallel orientation in the mutants, with the C-terminal helix and the loop connecting the helices displaying high flexibility, disorder, and non-native interactions. We validate these computational predictions via the anomalous fluorescence of a native tyrosine located at the interface of the helices. Our observations highlight the role of long loops in determining the unfolding mechanisms, sensitivity of the native ensembles to mutational perturbations and provide experimentally testable predictions that can be explored in even two-state folding systems.

Salt Dependence Conformational Stability of the Dimeric SAM Domain of MAPKKK Ste11 from Budding Yeast: A Native-State H/D Exchange NMR Study

The sterile α motif, also called the SAM domain, is known to form homo or heterocomplexes that modulate diverse biological functions through the regulation of specific protein-protein interactions. The MAPK pathway of budding yeast Saccharomyces cerevisiae is comprised of a three-tier kinase system akin to mammals. The MAPKKK Ste11 protein of yeast contains a homodimer SAM domain, which is critical for transmitting cues to the downstream kinases. The structural stability of the dimeric Ste11 SAM is maintained by hydrophobic and ionic interactions at the interfacial amino acids. The urea-induced equilibrium-unfolding process of the Ste11 SAM domain is cooperative without evidence of any intermediate states. The native-state H/D exchange under subdenaturing conditions is a useful method for the detection of intermediate states of proteins. In the present study, we investigated the effect of ionic strength on the conformational stability of the dimer using the H/D exchange experiments. The hydrogen exchange behavior of the Ste11 dimer under physiological salt concentrations reveals two partially unfolded metastable intermediate states, which may be generated by a sequential and cooperative unfolding of the five helices present in the domain. These intermediates appear to be significant for the reversible unfolding kinetics via hydrophobic collapse. In contrast, higher ionic concentrations eliminate this cooperative interactions that stabilize the pairs of helices.

Assessing Site-Specific Enhancements Imparted by Hyperpolarized Water in Folded and Unfolded Proteins by 2D HMQC NMR

Hyperpolarized water can be a valuable aid in protein NMR, leading to amide group ¹H polarizations that are orders of magnitude larger than their thermal counterparts. Suitable procedures can exploit this to deliver 2D ¹H–¹⁵N correlations with good resolution and enhanced sensitivity. These enhancements depend on the exchange rates between the amides and the water, thereby yielding diagnostic information about solvent accessibility. This study applied this “HyperW” method to four proteins exhibiting a gamut of exchange behaviors: PhoA^(350–471), an unfolded 122-residue fragment; barstar, a fully folded ribonuclease inhibitor; R17, a 13.3 kDa system possessing folded and unfolded forms under slow interconversion; and drkN SH3, a protein domain whose folded and unfolded forms interchange rapidly and with temperature-dependent population ratios. For PhoA4^(350–471) HyperW sensitivity enhancements were ≥300×, as expected for an unfolded protein sequence. Though fully folded, barstar also exhibited substantial enhancements; these, however, were not uniform and, according to CLEANEX experiments, reflected the solvent-exposed residues. R17 showed the expected superposition of ≥100-fold enhancements for its unfolded form, coexisting with more modest enhancements for their folded counterparts. Unexpected, however, was the behavior of drkN SH3, for which HyperW enhanced the unfolded but, surprisingly, enhanced even more certain folded protein sites. These preferential enhancements were repeatedly and reproducibly observed. A number of explanations—including three-site exchange magnetization transfers between water and the unfolded and folded states; cross-correlated relaxation processes from hyperpolarized “structural” waters and labile side-chain protons; and the possibility that faster solvent exchange rates characterize certain folded sites over their unfolded counterparts—are considered to account for them.

The Pierced Lasso Topology Leptin has a Bolt on Dynamic Domain Composed by the Disordered Loops I and III

Leptin is an important signaling hormone, mostly known for its role in energy expenditure and satiety. Furthermore, leptin plays a major role in other proteinopathies, such as cancer, marked hyperphagia, impaired immune function, and inflammation. In spite of its biological relevance in human health, there are no NMR resonance assignments of the human protein available, obscuring high-resolution characterization of the soluble protein and/or its conformational dynamics, suggested as being important for receptor interaction and biological activity. Here, we report the nearly complete backbone resonance assignments of human leptin. Chemical shift-based secondary structure prediction confirms that in solution leptin forms a four-helix bundle including a pierced lasso topology. The conformational dynamics, determined on several timescales, show that leptin is monomeric, has a rigid four-helix scaffold, and a dynamic domain, including a transiently formed helix. The dynamic domain is anchored to the helical scaffold by a secondary hydrophobic core, pinning down the long loops of leptin to the protein body, inducing motional restriction without a well-defined secondary or tertiary hydrogen bond stabilized structure. This dynamic region is well suited for and may be involved in functional allosteric dynamics upon receptor binding.

Structural Characterization of an ACP from Thermotoga maritima: Insights into Hyperthermal Adaptation

Thermotoga maritima, a deep-branching hyperthermophilic bacterium, expresses an extraordinarily stable Thermotoga maritima acyl carrier protein (Tm-ACP) that functions as a carrier in the fatty acid synthesis system at near-boiling aqueous environments. Here, to understand the hyperthermal adaptation of Tm-ACP, we investigated the structure and dynamics of Tm-ACP by nuclear magnetic resonance (NMR) spectroscopy. The melting temperature of Tm-ACP (101.4 °C) far exceeds that of other ACPs, owing to extensive ionic interactions and tight hydrophobic packing. The D59 residue, which replaces Pro/Ser of other ACPs, mediates ionic clustering between helices III and IV. This creates a wide pocket entrance to facilitate the accommodation of long acyl chains required for hyperthermal adaptation of the T. maritima cell membrane. Tm-ACP is revealed to be the first ACP that harbor an amide proton hyperprotected against hydrogen/deuterium exchange for I15. The hydrophobic interactions mediated by I15 appear to be the key driving forces of the global folding process of Tm-ACP. Our findings provide insights into the structural basis of the hyperthermal adaptation of ACP, which might have allowed T. maritima to survive in hot ancient oceans.

Bridging protein structure, dynamics, and function using hydrogen/deuterium-exchange mass spectrometry

Much of our understanding of protein structure and mechanistic function has been derived from static high-resolution structures. As structural biology has continued to evolve it has become clear that high-resolution structures alone are unable to fully capture the mechanistic basis for protein structure and function in solution. Recently Hydrogen/Deuterium-exchange Mass Spectrometry (HDX-MS) has developed into a powerful and versatile tool for structural biologists that provides novel insights into protein structure and function. HDX-MS enables direct monitoring of a protein's structural fluctuations and conformational changes under native conditions in solution even as it is carrying out its functions. In this review, we focus on the use of HDX-MS to monitor these dynamic changes in proteins. We examine how HDX-MS has been applied to study protein structure and function in systems ranging from large, complex assemblies to intrinsically disordered proteins, and we discuss its use in probing conformational changes during protein folding and catalytic function. STATEMENT FOR A BROAD AUDIENCE: The biophysical and structural characterization of proteins provides novel insight into their functionalities. Protein motions, ranging from small scale local fluctuations to larger concerted structural rearrangements, often determine protein function. Hydrogen/Deuterium-exchange Mass Spectrometry (HDX-MS) has proven a powerful biophysical tool capable of probing changes in protein structure and dynamic protein motions that are often invisible to most other techniques.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Friction-Limited Folding of Disulfide-Reduced Monomeric SOD1

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

Structure of an Unfolding Intermediate of an RRM Domain of ETR-3 Reveals Its Native-like Fold

Prevalence of one or more partially folded intermediates during protein unfolding with different secondary and ternary conformations has been identified as an integral character of protein unfolding. These transition-state species need to be characterized structurally for elucidation of their folding pathways. We have determined the three-dimensional structure of an intermediate state with increased conformational space sampling under urea-denaturing condition. The protein unfolds completely at 10 M urea but retains residual secondary structural propensities with restricted motion. Here, we describe the native state, observable intermediate state, and unfolded state for ETR-3 RRM-3, which has canonical RRM fold. These observations can shed more light on unfolding events for RRM-containing proteins.

Theoretical and computational advances in protein misfolding

Misfolded proteins escape the cellular quality control mechanism and fail to fold properly or remain correctly folded leading to a loss in their functional specificity. Thus misfolding of proteins cause a large number of very different diseases ranging from errors in metabolism to various types of complex neurodegenerative diseases. A theoretical and computational perspective of protein misfolding is presented with a special emphasis on its salient features, mechanism and consequences. These insights quantitatively analyze different determinants of misfolding, that may be applied to design disease specific molecular targets.

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

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

Life in Phases: Intra- and Inter- Molecular Phase Transitions in Protein Solutions

Proteins, these evolutionarily-edited biological polymers, are able to undergo intramolecular and intermolecular phase transitions. Spontaneous intramolecular phase transitions define the folding of globular proteins, whereas binding-induced, intra- and inter- molecular phase transitions play a crucial role in the functionality of many intrinsically-disordered proteins. On the other hand, intermolecular phase transitions are the behind-the-scenes players in a diverse set of macrosystemic phenomena taking place in protein solutions, such as new phase nucleation in bulk, on the interface, and on the impurities, protein crystallization, protein aggregation, the formation of amyloid fibrils, and intermolecular liquid-liquid or liquid-gel phase transitions associated with the biogenesis of membraneless organelles in the cells. This review is dedicated to the systematic analysis of the phase behavior of protein molecules and their ensembles, and provides a description of the major physical principles governing intramolecular and intermolecular phase transitions in protein solutions.

Interlaboratory Comparison of Hydrogen-Deuterium Exchange Mass Spectrometry Measurements of the Fab Fragment of NISTmAb

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is an established, powerful tool for investigating protein-ligand interactions, protein folding, and protein dynamics. However, HDX-MS is still an emergent tool for quality control of biopharmaceuticals and for establishing dynamic similarity between a biosimilar and an innovator therapeutic. Because industry will conduct quality control and similarity measurements over a product lifetime and in multiple locations, an understanding of HDX-MS reproducibility is critical. To determine the reproducibility of continuous-labeling, bottom-up HDX-MS measurements, the present interlaboratory comparison project evaluated deuterium uptake data from the Fab fragment of NISTmAb reference material (PDB: 5K8A ) from 15 laboratories. Laboratories reported ∼89 800 centroid measurements for 430 proteolytic peptide sequences of the Fab fragment (∼78 900 centroids), giving ∼100% coverage, and ∼10 900 centroid measurements for 77 peptide sequences of the Fc fragment. Nearly half of peptide sequences are unique to the reporting laboratory, and only two sequences are reported by all laboratories. The majority of the laboratories (87%) exhibited centroid mass laboratory repeatability precisions of ⟨ sLab⟩ ≤ (0.15 ± 0.01) Da (1σx̅). All laboratories achieved ⟨sLab⟩ ≤ 0.4 Da. For immersions of protein at THDX = (3.6 to 25) °C and for D2O exchange times of tHDX = (30 s to 4 h) the reproducibility of back-exchange corrected, deuterium uptake measurements for the 15 laboratories is σreproducibility15 Laboratories( tHDX) = (9.0 ± 0.9) % (1σ). A nine laboratory cohort that immersed samples at THDX = 25 °C exhibited reproducibility of σreproducibility25C cohort( tHDX) = (6.5 ± 0.6) % for back-exchange corrected, deuterium uptake measurements.

Collapse Transitions of Proteins and the Interplay Among Backbone, Sidechain, and Solvent Interactions

Proteins can collapse into compact globules or form expanded, solvent-accessible, coil-like conformations. Additionally, they can fold into well-defined three-dimensional structures or remain partially or entirely disordered. Recent discoveries have shown that the tendency for proteins to collapse or remain expanded is not intrinsically coupled to their ability to fold. These observations suggest that proteins do not have to form compact globules in aqueous solutions. They can be intrinsically disordered, collapsed, or expanded, and even form well-folded, elongated structures. This ability to decouple collapse from folding is determined by the sequence details of proteins. In this review, we highlight insights gleaned from studies over the past decade. Using a polymer physics framework, we explain how the interplay among sidechains, backbone units, and solvent determines the driving forces for collapsed versus expanded states in aqueous solvents.

Multi-step conformational transitions in heat-treated protein therapeutics can be monitored in real time with temperature-controlled electrospray ionization mass spectrometry

Heat-induced conformational transitions are frequently used to probe the free energy landscapes of proteins. However, the extraction of information from thermal denaturation profiles pertaining to non-native protein conformations remains challenging due to their transient nature and significant conformational heterogeneity. Previously we developed a temperature-controlled electrospray ionization (ESI) source that allowed unfolding and association of biopolymers to be monitored by mass spectrometry (MS) in real time as a function of temperature. The scope of this technique is now extended to systems that undergo multi-step denaturation upon heat stress, as well as relatively small-scale conformational changes that are precursors to protein aggregation. The behavior of two therapeutic proteins (human antithrombin and an IgG1 monoclonal antibody) under heat-stress conditions is monitored in real time, providing evidence that relatively small-scale conformational changes in each system lead to protein oligomerization, followed by aggregation. Temperature-controlled ESI MS is particularly useful for the studies of heat-stressed multi-domain proteins such as IgG, where it allows distinct transitions to be observed. The ability of native temperature-controlled ESI MS to monitor both the conformational changes and oligomerization/degradation with high selectivity complements the classic calorimetric methods, lending itself as a powerful experimental tool for the thermostability studies of protein therapeutics.

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

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

Globular Protein Folding In Vitro and In Vivo

In vitro, computational, and theoretical studies of protein folding have converged to paint a rich and complex energy landscape. This landscape is sensitively modulated by environmental conditions and subject to evolutionary pressure on protein function. Of these environments, none is more complex than the cell itself, where proteins function in the cytosol, in membranes, and in different compartments. A wide variety of kinetic and thermodynamics experiments, ranging from single-molecule studies to jump kinetics and from nuclear magnetic resonance to imaging on the microscope, have elucidated how protein energy landscapes facilitate folding and how they are subject to evolutionary constraints and environmental perturbation. Here we review some recent developments in the field and refer the reader to some original work and additional reviews that cover this broad topic in protein science.

Dynamic diversity of synthetic supramolecular polymers in water as revealed by hydrogen/deuterium exchange

Numerous self-assembling molecules have been synthesized aiming at mimicking both the structural and dynamic properties found in living systems. Here we show the application of hydrogen/deuterium exchange (HDX) mass spectrometry (MS) to unravel the nanoscale organization and the structural dynamics of synthetic supramolecular polymers in water. We select benzene-1,3,5-tricarboxamide (BTA) derivatives that self-assemble in H₂O to illustrate the strength of this technique for supramolecular polymers. The BTA structure has six exchangeable hydrogen atoms and we follow their exchange as a function of time after diluting the H₂O solution with a 100-fold excess of D₂O. The kinetic H/D exchange profiles reveal that these supramolecular polymers in water are dynamically diverse; a notion that has previously not been observed using other techniques. In addition, we report that small changes in the molecular structure can be used to control the dynamics of synthetic supramolecular polymers in water.

In vitro characterization of six STUB1 variants in spinocerebellar ataxia 16 reveals altered structural properties for the encoded CHIP proteins

Spinocerebellar ataxia, autosomal recessive 16 (SCAR16) is caused by biallelic mutations in the STIP1 homology and U-box containing protein 1 (STUB1) gene encoding the ubiquitin E3 ligase and dimeric co-chaperone C-terminus of Hsc70-interacting protein (CHIP). It has been proposed that the disease mechanism is related to CHIP's impaired E3 ubiquitin ligase properties and/or interaction with its chaperones. However, there is limited knowledge on how these mutations affect the stability, folding, and protein structure of CHIP itself. To gain further insight, six previously reported pathogenic STUB1 variants (E28K, N65S, K145Q, M211I, S236T, and T246M) were expressed as recombinant proteins and studied using limited proteolysis, size-exclusion chromatography (SEC), and circular dichroism (CD). Our results reveal that N65S shows increased CHIP dimerization, higher levels of α-helical content, and decreased degradation rate compared with wild-type (WT) CHIP. By contrast, T246M demonstrates a strong tendency for aggregation, a more flexible protein structure, decreased levels of α-helical structures, and increased degradation rate compared with WT CHIP. E28K, K145Q, M211I, and S236T also show defects on structural properties compared with WT CHIP, although less profound than what observed for N65S and T246M. In conclusion, our results illustrate that some STUB1 mutations known to cause recessive SCAR16 have a profound impact on the protein structure, stability, and ability of CHIP to dimerize in vitro. These results add to the growing understanding on the mechanisms behind the disorder.

Impact of Turn Propensity on the Folding Rates of Z34C Protein: Implications for the Folding of Helix-Turn-Helix Motif

The rate-limiting step for the folding of the helix-turn-helix (HTH) protein, Z34C, involves β-turn region ²⁰DPNL²³. This reverse turn has been observed to be part of the transition state in the folding process for Z34C, influencing its folding rates. Molecular dynamics simulations were performed on this turn peptide and its two mutants, D20A and P21A, to study turn formation using GROMOS54A7 force field. We find that this region has a turn propensity of its own, and the highest turn propensity is observed for the wild-type, which correlates well with available experimental results. We also find that a slight unfavorable change in ΔG _{turn folding} causes a drastic change in the folding rates of HTH motif and a mechanistic interpretation is given. Implications of these observations for the folding of the HTH protein Z34C are discussed.

Toward a quantitative description of microscopic pathway heterogeneity in protein folding

Experimentally consistent statistical modeling of protein folding thermodynamics reveals unprecedented complexity with numerous parallel folding routes in five different proteins.

Ubiquitin utilizes an acidic surface patch to alter chromatin structure

Ubiquitylation of histone H2B, associated with gene activation, leads to chromatin decompaction through an unknown mechanism. We used a hydrogen-deuterium exchange strategy coupled with nuclear magnetic resonance spectroscopy to map the ubiquitin surface responsible for its structural effects on chromatin. Our studies revealed that a previously uncharacterized acidic patch on ubiquitin comprising residues Glu16 and Glu18 is essential for decompaction. These residues mediate promiscuous electrostatic interactions with the basic histone proteins, potentially positioning the ubiquitin moiety as a dynamic “wedge” that prevents the intimate association of neighboring nucleosomes. Using two independent cross-linking strategies and an oligomerization assay, we also showed that ubiquitin-ubiquitin contacts occur in the chromatin environment and are important for the solubilization of the chromatin polymers. Our work highlights a novel, chromatin-related aspect of the “ubiquitin code”, and sheds light on how the information rich ubiquitin modification can orchestrate different biochemical outcomes using different surface features.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.

NTL9 Folding at Constant pH: The Importance of Electrostatic Interaction and pH Dependence

The folding process of the N-terminal domain of ribosomal protein L9 (NTL9) was investigated at constant-pH computer simulations. Evaluation of the role of electrostatic interaction during folding was carried out by including a Debye-Hückel potential into a Cα structure-based model (SBM). In this study, the charges of the ionizable residues and the electrostatic potential are susceptible to the solution conditions, such as pH and ionic strength, as well as to the presence of charged groups. Simulations were performed under different pHs, and the results were validated by comparing them with experimental values of pKa and with denaturation experiment data. Also, the free energy profiles, Φ-values, and folding routes were calculated for each condition. It was shown how charges vary along the folding under different pH, which is subject to different scenarios. This study reveals how simplified models can capture essential physical features, reproducing experimental results, and presenting the role of electrostatic interactions before, during, and after the transition state.

Effect of circular permutations on transient partial unfolding in proteins

Under native conditions, proteins can undergo transient partial unfolding, which may cause proteins to misfold or aggregate. A change in sequence connectivity by circular permutation may affect the energetics of transient partial unfolding in proteins without altering the three-dimensional structures. Using Escherichia coli dihydrofolate reductase (DHFR) as a model system, we investigated how circular permutation affects transient partial unfolding in proteins. We constructed three circular permutants, CP18, CP37, and CP87, with the new N-termini at residue 18, 37, and 87, respectively, and probed transient partial unfolding by native-state proteolysis. The new termini in CP18, CP37, and CP87 are within, near, and distal to the Met20 loop, which is known to be dynamic and also part of the region that undergoes transient unfolding in wild-type DHFR. The stabilities of both native and partially unfolded forms of CP18 are similar to those of wild-type DHFR, suggesting that the influence of introducing new termini in a dynamic region to the protein is minimal. CP37 has a significantly more accessible partially unfolded form than wild-type DHFR, demonstrating that introducing new termini near a dynamic region may promote transient partial unfolding. CP87 has significantly destabilized native and partially unfolded forms, confirming that modification of the folded region in a partially unfolded form destabilizes the partially unfolded form similar to the native form. Our findings provide valuable guidelines to control transient partial unfolding in designing circular permutants in proteins.